The lunar Askaryan technique is a method to study the highest-energy cosmic rays, and their predicted counterparts, the ultra-high-energy neutrinos. By observing the Moon with a radio telescope, and searching for the characteristic nanosecond-scale Askaryan pulses emitted when a high-energy particle interacts in the outer layers of the Moon, the visible lunar surface can be used as a detection area. Several previous experiments, at Parkes, Goldstone, Kalyazin, Westerbork, the ATCA, Lovell, LOFAR, and the VLA, have developed the necessary techniques to search for these pulses, but existing instruments have lacked the necessary sensitivity to detect the known flux of cosmic rays from such a distance. This will change with the advent of the SKA. The Square Kilometre Array (SKA) will be the world's most powerful radio telescope. To be built in southern Africa, Australia and New Zealand during the next decade, it will have an unsurpassed sensitivity over the key 100 MHz to few-GHZ band. We introduce a planned experiment to use the SKA to observe the highest-energy cosmic rays and, potentially, neutrinos. The estimated event rate will be presented, along with the predicted energy and directional resolution. Prospects for directional studies with phase 1 of the SKA will be discussed, as will the major technical challenges to be overcome to make full use of this powerful instrument. Finally, we show how phase 2 of the SKA could provide a vast increase in the number of detected cosmic rays at the highest energies, and thus to provide new insight into their spectrum and origin.

The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely large detector of ultra-high-energy particles. The high time resolution required to detect these pulses, which have a duration of around a nanosecond, puts this technique in a regime quite different from other forms of radio astronomy, with a unique set of associated technical challenges which have been addressed in a series of experiments by various groups. Implementing the methods and techniques developed by these groups for detecting lunar Askaryan pulses will be important for a future experiment with the Square Kilometre Array (SKA), which is expected to have sufficient sensitivity to allow the first positive detection using this technique. Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and the exclusion of spurious events from radio-frequency interference. We review the progress in each of these areas, and consider the further progress expected for future application with the SKA.

We report a limit on the ultra-high-energy neutrino flux based on a non-detection of radio pulses from neutrino-initiated particle cascades in the Moon, in observations with the Parkes radio telescope undertaken as part of the LUNASKA project. Due to the improved sensitivity of these observations, which had an effective duration of 127 hours and a frequency range of 1.2-1.5 GHz, this limit extends to lower neutrino energies than those from previous lunar radio experiments, with a detection threshold below 10^20 eV. The calculation of our limit allows for the possibility of lunar-origin pulses being misidentified as local radio interference, and includes the effect of small-scale lunar surface roughness. The targeting strategy of the observations also allows us to place a directional limit on the neutrino flux from the nearby radio galaxy Centaurus A.

Observations of the properties of dense molecular clouds are critical in understanding the process of star-formation. One of the most important, but least understood, is the role of the magnetic fields. We discuss the possibility of using high-resolution, high-sensitivity radio observations with the SKA to measure for the first time the in-situ synchrotron radiation from these molecular clouds. If the cosmic-ray (CR) particles penetrate clouds as expected, then we can measure the B-field strength directly using radio data. So far, this signature has never been detected from the collapsing clouds themselves and would be a unique probe of the magnetic field. Dense cores are typically ~0.05 pc in size, corresponding to ~arcsec at ~kpc distances, and flux density estimates are ~mJy at 1 GHz. The SKA should be able to readily detect directly, for the first time, along lines-of-sight that are not contaminated by thermal emission or complex foreground/background synchrotron emission. Polarised synchrotron may also be detectable providing additional information about the regular/turbulent fields.

The origin of the most energetic particles in nature, the ultra-high-energy (UHE) cosmic rays, is still a mystery. Only the most energetic of these have sufficiently small angular deflections to be used for directional studies, and their flux is so low that even the 3,000 km^2 Pierre Auger detector registers only about 30 cosmic rays per year of these energies. A method to provide an even larger aperture is to use the lunar Askaryan technique, in which ground-based radio telescopes search for the nanosecond radio flashes produced when a cosmic ray interacts with the Moon's surface. The technique is also sensitive to UHE neutrinos, which may be produced in the decays of topological defects from the early universe. Observations with existing radio telescopes have shown that this technique is technically feasible, and established the required procedure: the radio signal should be searched for pulses in real time, compensating for ionospheric dispersion and filtering out local radio interference, and candidate events stored for later analysis. For the Square Kilometre Array (SKA), this requires the formation of multiple tied-array beams, with high time resolution, covering the Moon, with either SKA1-LOW or SKA1-MID. With its large collecting area and broad bandwidth, the SKA will be able to detect the known flux of UHE cosmic rays using the visible lunar surface - millions of square km - as the detector, providing sufficient detections of these extremely rare particles to address the mystery of their origin.

We describe an experiment using the Parkes radio telescope in the 1.2-1.5 GHz frequency range as part of the LUNASKA project, to search for nanosecond-scale pulses from particle cascades in the Moon, which may be triggered by ultra-high-energy astroparticles. Through the combination of a highly sensitive multi-beam radio receiver, a purpose-built backend and sophisticated signal-processing techniques, we achieve sensitivity to radio pulses with a threshold electric field strength of 0.0053 $\mu$V/m/MHz, lower than previous experiments by a factor of three. We observe no pulses in excess of this threshold in observations with an effective duration of 127 hours. The techniques we employ, including compensating for the phase, dispersion and spectrum of the expected pulse, are relevant for future lunar radio experiments.

Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultra-high-energy cosmic-ray source. We apply several methods of determining jet power and affirm the consistency of various power estimates of ~ 1 x 10^43 erg s^-1. Employing scaling relations based on previous results for 3C 31, we estimate particle number densities in the jets, encompassing available radio through X-ray observations. Our model is compatible with the jets ingesting ~ 3 x 10^21 g s^-1 of matter via external entrainment from hot gas and ~ 7 x 10^22 g s^-1 via internal entrainment from jet-contained stars. This leads to an imbalance between the internal lobe pressure available from radiating particles and magnetic field, and our derived external pressure. Based on knowledge of the external environments of other FR I sources, we estimate the thermal pressure in the giant lobes as 1.5 x 10^-12 dyn cm^-2, from which we deduce a lower limit to the temperature of ~ 1.6 x 10^8 K. Using dynamical and buoyancy arguments, we infer ~ 440-645 Myr and ~ 560 Myr as the sound-crossing and buoyancy ages of the giant lobes respectively, inconsistent with their spectral ages. We re-investigate the feasibility of particle acceleration via stochastic processes in the lobes, placing new constraints on the energetics and on turbulent input to the lobes. The same 'very hot' temperatures that allow self-consistency between the entrainment calculations and the missing pressure also allow stochastic UHECR acceleration models to work.

The Moon is used as a target volume for ultra-high energy neutrino searches with terrestrial radio telescopes. The LUNASKA project has conducted observations with the Parkes and ATCA telescopes; and, most recently, with both of them in combination. We present an analysis of the data obtained from these searches, including validation and calibration results for the Parkes-ATCA experiment, as well as a summary of prospects for future observations.

We present a time-dependent approach to the one-zone hadronic model in the case where the photon spectrum is produced by ultrarelativistic protons interacting with soft photons that are produced from protons and low magnetic fields. Assuming that protons are injected at a certain rate in a homogeneous spherical volume containing a magnetic field, the evolution of the system can be described by five coupled kinetic equations, for protons, electrons, photons, neutrons, and neutrinos. Photopair and photopion interactions are modelled using the results of Monte-Carlo simulations and, in particular from the SOPHIA code for the latter. The coupling of energy losses and injection introduces a self-consistency in our approach and allows the study of the comparative relevancy of processes at various conditions, the efficiency of the conversion of proton luminosity to radiation, the resulting neutrino spectra, and the effects of time variability on proton injection, among other topics. We present some characteristic examples of the temporal behaviour of the system and show that this can be very different from the one exhibited by leptonic models. Furthermore, we argue that, contrary to the wide-held belief, there are parameter regimes where the hadronic models can become quite efficient. However, to keep the free parameters at a minimum and facilitate an in-depth study of the system, we have only concentrated on the case where protons are injected; i.e., we did not consider the effects of a co-accelerated leptonic component.

The most sensitive method for detecting neutrinos at the very highest energies is the lunar Cherenkov technique, which employs the Moon as a target volume, using conventional radio telescopes to monitor it for nanosecond-scale pulses of Cherenkov radiation from particle cascades in its regolith. Multiple-antenna radio telescopes are difficult to effectively combine into a single detector for this purpose, while single antennas are more susceptible to false events from radio interference, which must be reliably excluded for a credible detection to be made. We describe our progress in excluding such interference in our observations with the single-antenna Parkes radio telescope, and our most recent experiment (taking place the week before the ICRC) using it in conjunction with the Australia Telescope Compact Array, exploiting the advantages of both types of telescope.

We present images of the Sagittarius (Sgr) B giant molecular cloud at 2368 and 1384 MHz obtained using new, multi-configuration Australia Telescope Compact Array (ATCA) observations. We have combined these observations with archival single-dish observations yielding images at resolutions of 47" by 14" and 27" by 8" at 1384 and 2368 MHz respectively. These observations were motivated by our theoretical work (Protheroe et al. 2008) indicating the possibility that synchrotron emission from secondary electrons and positrons created in hadronic cosmic ray (CR) collisions with the ambient matter of the Sgr B2 cloud could provide a detectable (and possibly linearly polarized) non-thermal radio signal. We find that the only detectable non-thermal emission from the Sgr B region is from a strong source to the south of Sgr B2, which we label Sgr B2 Southern Complex (SC). We find Sgr B2(SC) integrated flux densities of 1.2+/-0.2 Jy at 1384 MHz and 0.7+/-0.1 Jy at 2368 MHz for a source of FWHM size at 1384 MHz of ~54". Despite its non-thermal nature, the synchrotron emission from this source is unlikely to be dominantly due to secondary electrons and positrons. We use polarization data to place 5-sigma upper limits on the level of polarized intensity from the Sgr B2 cloud of 3.5 and 3 mJy/beam at 1384 and 2368 MHz respectively. We also use the angular distribution of the total intensity of archival 330 MHz VLA and the total intensity and polarized emission of our new 1384 MHz and 2368 MHz data to constrain the diffusion coefficient for transport of the parent hadronic CRs into the dense core of Sgr B2 to be no larger than about 1% of that in the Galactic disk. Finally, we have also used the data to perform a spectral and morphological study of the features of the Sgr B cloud and compare and contrast these to previous studies.

LUNASKA (Lunar UHE Neutrino Astrophysics with the Square Kilometre Array) is an ongoing project conducting lunar Cherenkov observations in order to develop techniques for detecting neutrinos with the next generation of radio telescopes. Our current observing campaign is with the 64-metre Parkes radio telescope, using a multibeam receiver with 300 MHz of bandwidth from 1.2-1.5 GHz. Here we provide an overview of the various factors that must be considered in the signal processing for such an experiment. We also briefly describe the flux limits which we expect to set with our current observations, including a directional limit for Centaurus A.

We present the first observational limits to the ultra-high energy (UHE) neutrino flux from the Galactic Center, and from Centaurus A which is the nearest active galactic nucleus (AGN). These results are based on our "Lunar UHE Neutrino Astrophysics using the Square Kilometer Array" (LUNASKA) project experiments at the Australia Telescope Compact Array (ATCA). We also derive limits for the previous experiments and compare these limits with expectations for acceleration and super-heavy dark matter models of the origin of UHE cosmic rays.

We describe the design, performance, sensitivity and results of our recent experiments using the Australia Telescope Compact Array (ATCA) for lunar Cherenkov observations with a very wide (600 MHz) bandwidth and nanosecond timing, including a limit on an isotropic neutrino flux. We also make a first estimate of the effects of small-scale surface roughness on the effective experimental aperture, finding that contrary to expectations, such roughness will act to increase the detectability of near-surface events over the neutrino energy-range at which our experiment is most sensitive (though distortions to the time-domain pulse profile may make identification more difficult). The aim of our "Lunar UHE Neutrino Astrophysics using the Square Kilometer Array" (LUNASKA) project is to develop the lunar Cherenkov technique of using terrestrial radio telescope arrays for ultra-high energy (UHE) cosmic ray (CR) and neutrino detection, and in particular to prepare for using the Square Kilometer Array (SKA) and its path-finders such as the Australian SKA Pathfinder (ASKAP) and the Low Frequency Array (LOFAR) for lunar Cherenkov experiments.

The amplitude of the magnetic field near the Galactic Centre has been uncertain by two orders of magnitude for several decades. On a scale of approximately 100 pc fields of approximately 1000 microG have been reported, implying a magnetic energy density more than 10,000 times stronger than typical for the Galaxy. Alternatively, the assumption of pressure equilibrium between the various phases of the Galactic Centre interstellar medium (including turbulent molecular gas; the contested "very hot" plasma; and the magnetic field) suggests fields of approximately 100 microG over approximately 400 pc size scales. Finally, assuming equipartition, fields of only approximately 6 microG have been inferred from radio observations for 400 pc scales. Here we report a compilation of previous data that reveals a down-break in the region's non-thermal radio spectrum (attributable to a transition from bremsstrahlung to synchrotron cooling of the in situ cosmic-ray electron population). We show that the spectral break requires that the Galactic Centre field be at least 50 microG on 400 pc scales, lest the synchrotron-emitting electrons produce too much gamma-ray emission given existing constraints. Other considerations support a field of 100 microG, implying that > 10% of the Galaxy's magnetic energy is contained in only < 0.05% of its volume.

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.

This contribution describes the experimental set-up implemented by the LUNASKA project at the Australia Telescope Compact Array (ATCA) to enable the radio-telescope to be used to search for pulses of coherent Cherenkov radiation from UHE particle interactions in the Moon with an unprecedented bandwidth, and hence sensitivity. Our specialised hardware included analogue de-dispersion filters to coherently correct for the dispersion expected of a ~nanosecond pulse in the Earth's ionosphere over our wide (600 MHz) bandwidth, and FPGA-based digitising boards running at 2.048 GHz for pulse detection. The trigger algorithm is described, as are the methods used discriminate between terrestrial RFI and true lunar pulses. We also outline the next stage of hardware development expected to be used in our 2010 observations.

The lunar Cherenkov technique is a method to use radio-telescopes to detect ultra-high energy cosmic rays (CR) and neutrinos ($\nu$). By observing the short-duration ($\sim$few nanosecond) pulses of coherent Cherenkov radiation emitted from particle cascades via the Askaryan Effect in the Moon's outer layers (nominally the regolith), the primary particles initiating the cascades may be identified. Our collaboration (LUNASKA) aims to develop the technique to be used with the next generation of giant radio-arrays. Here, we present the results of our two preliminary UHE particle searches using this technique with three antennas at the Australia Telescope Compact Array (ATCA) during February and May 2008.

In this paper, we study the diffusive shock acceleration of cosmic-ray protons and nuclei, taking into account all the relevant interaction processes with photon backgrounds. We investigate how the competition between protons and nuclei is modified by the acceleration parameters such as the acceleration rate, its rigidity dependence, the photon density and the confinement capability of the sources. We find that in the case of interaction-limited acceleration processes protons are likely to be accelerated to higher energies than nuclei, whereas for confinement-limited acceleration nuclei are accelerated to higher energies than protons. Finally, we discuss our results in the context of possible astrophysical accelerators, and in the light of recent cosmic-ray data.

LUNASKA (Lunar UHE Neutrino Astrophysics with the Square Kilometre Array) is a theoretical and experimental project developing the lunar Cherenkov technique for the next generation of giant radio-telescope arrays. Here we report on a series of observations with ATCA (the Australia Telescope Compact Array). Our current observations use three of the six 22m ATCA antennas with a 600 MHz bandwidth at 1.2-1.8 GHz, analogue dedispersion filters to correct for the typical night-time ionospheric dispersion, and state-of-the-art 2 GHz FPGA-based digital pulse detection hardware. We have observed so as to maximise the UHE neutrino sensitivity in the region surrounding the galactic centre and to Centaurus A, to which current limits on the highest-energy neutrinos are relatively weak.

The LUNASKA (Lunar UHE Neutrino Astrophysics with the Square Kilometre Array) project is a theoretical and experimental project developing the lunar Cherenkov technique for the next generation of giant radio-telescope arrays. This contribution presents our simulation results on the directional dependence of the technique for UHE neutrino detection. In particular, these indicate that both the instantaneous sensitivities and time-integrated limits from lunar Cherenkov experiments such as those at Parkes, Goldstone, Kalyazin and ATCA are highly anisotropic. We study the regions of the sky which have not been probed by either these or other experiments, and present the expected sky coverage of future experiments with the SKA. Our results show how the sensitivity of Lunar Cherenkov observations to potential astrophysical sources of UHE particles may be maximised by choosing appropriate observations dates and antenna-beam pointing positions.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

Recent ammonia (1,1) inversion line data on the Galactic star forming region Sgr B2 show that the column density is consistent with a radial Gaussian density profile with a standard deviation of 2.75 pc. Deriving a formula for the virial mass of spherical Gaussian clouds, we obtain a virial mass of 1.9 million solar masses for Sgr B2. For this matter distribution, a reasonable magnetic field and an impinging flux of cosmic rays of solar neighbourhood intensity, we predict the expected synchrotron emission from the Sgr B2 giant molecular cloud due to secondary electrons and positrons resulting from cosmic ray interactions, including effects of losses due to pion production collisions during diffusive propagation into the cloud complex. We assemble radio continuum data at frequencies between 330 MHz and 230 GHz. From the spectral energy distribution the emission appears to be thermal at all frequencies. Before using these data to constrain the predicted synchrotron flux, we first model the spectrum as free-free emission from the known ultra compact HII regions plus emission from an envelope or wind with a radial density gradient. This severely constrains the possible synchrotron emission by secondary electrons to quite low flux levels. The absence of a significant contribution by secondary electrons is almost certainly due to multi-GeV energy cosmic rays being unable to penetrate far into giant molecular clouds. This would also explain why 100 MeV--GeV gamma-rays (from neutral pion decay or bremsstrahlung by secondary electrons) were not observed from Sgr B2 by EGRET, while TeV energy gamma-rays were observed, being produced by higher energy cosmic rays which more readily penetrate giant molecular clouds.

We report radio continuum observations with the Australia Telescope Compact Array of two molecular clouds. The impetus for these observations is a search for synchrotron radiation by cosmic ray secondary electrons/positrons in a region of enhanced density and possibly high magnetic field. We present modelling which shows that there should be an appreciable flux of synchrotron above the more diffuse, galactic synchrotron background. The starless core G333.125-0.562 and infrared source IRAS 15596-5301 were observed at 1384 and 2368 MHz. For G333.125-0.562, we find no significant levels of radio emission from this source at either frequency, nor any appreciable polarisation: we place an upper limit on the radio continuum flux from this source of 0.5 mJy per beam at both 1384 and 2368 MHz. Due to the higher than expected flux density limits, we also obtained archival ATCA data at 8640 MHz for this cloud and place an upper limit on the flux density of 50 micro-Jy per beam. Assuming the cosmic ray spectrum is similar to that near the Sun, and given the cloud's molecular density and mass, we place an upper limit on the magnetic field of 500 micro-G. IRAS 15596-5301, with an RMS of 50 micro-Jy per beam at 1384 MHz, shows an HII region consistent with optically thin free-free emission already detected at 4800 MHz. We use the same prescription as G333 to constrain the magnetic field from this cloud to be less than 500 micro-G. We find that these values are not inconsistent with the view that magnetic field values scale with the average density of the molecular cloud.

We use computer simulations to obtain the directional-dependence of the lunar Cherenkov technique for ultra-high energy (UHE) neutrino detection. We calculate the instantaneous effective area of past lunar Cherenkov experiments as a function of neutrino arrival direction, and hence the directional-dependence of the combined limit imposed by GLUE and the experiment at Parkes. We also determine the directional dependence of the aperture of future planned experiments with ATCA, ASKAP and the SKA to a UHE neutrino flux, and calculate the potential annual exposure to astronomical objects as a function of angular distance from the lunar trajectory through celestial coordinates.

We present simulation results for the detection of ultra-high energy (UHE) cosmic ray (CR) and neutrino interactions in the Moon by radio-telescopes. We simulate the expected radio signal at Earth from such interactions, expanding on previous work to include interactions in the sub-regolith layer for single dish and multiple telescope systems. For previous experiments at Parkes, Goldstone, and Kalyazin we recalculate the sensitivity to an isotropic flux of UHE neutrinos. Our predicted sensitivity for future experiments using the Australia Telescope Compact Array (ATCA) and the Australian SKA Pathfinder (ASKAP) indicate these instruments will be able to detect the more optimistic UHE neutrino flux predictions, while the Square Kilometre Array (SKA) will also be sensitive to all bar one prediction of a diffuse `cosmogenic', or `GZK', neutrino flux. Current uncertainties concerning the structure and roughness of the lunar surface prevents an accurate calculation of the sensitivity of the lunar Cherenkov technique for UHE cosmic ray astronomy at high frequencies. However, below 200 MHz we find that the proposed SKA low-frequency aperture array should be able to detect events above 56 EeV at a rate about 30 times that of the current Pierre Auger Observatory. This would allow directional analysis of UHE cosmic rays, and investigation of correlations with putative cosmic ray source populations, to be conducted with very high statistics.

A promising method for the detection of UHE neutrinos is the Lunar Cherenkov technique, which utilises Earth-based radio telescopes to detect the coherent Cherenkov radiation emitted when a UHE neutrino interacts in the outer layers of the Moon. The LUNASKA project aims to overcome the technological limitations of past experiments to utilise the next generation of radio telescopes in the search for these elusive particles. To take advantage of broad-bandwidth data from potentially thousands of antennas requires advances in signal processing technology. Here we describe recent developments in this field and their application in the search for UHE neutrinos, from a preliminary experiment using the first stage of an upgrade to the Australia Telescope Compact Array, to possibilities for fully utilising the completed Square Kilometre Array. We also explore a new real time technique for characterising ionospheric pulse dispersion which specifically measures ionospheric electron content that is line of sight to the moon.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

TenTen is a proposed array of Imaging Atmospheric Cherenkov Telescopes (IACT) optimized for the gamma ray energy regime of 10 TeV to 100 TeV, but with a threshold of ~1 to a few TeV. It will offer a collecting area of 10 km2 above energies of 10 TeV. In the initial phase, a cell of 3 to 5 modest-sized telescopes, each with 10-30 m2 mirror area, is suggested for an Australian site. A possible expansion of the array could comprise many such cells. Here we present work on configuration and technical issues from our simulation studies of the array. Working topics include array layout, telescope size and optics, camera field of view, telescope trigger system, electronics, and site surveys.

Remnants of AGN jets and their surrounding cocoons leave colossal magnetohydrodynamic (MHD) fossil structures storing total energies ~10^{60} erg. The original active galacic nucleus (AGN) may be dead but the fossil will retain its stable configuration resembling the reversed-field pinch (RFP) encountered in laboratory MHD experiments. In an RFP the longitudinal magnetic field changes direction at a critical distance from the axis, leading to magnetic re-connection there, and to slow decay of the large-scale RFP field. We show that this field decay induces large-scale electric fields which can accelerate cosmic rays with an E^{-2} power-law up to ultra-high energies with a cut-off depending on the fossil parameters. The cut-off is expected to be rigidity dependent, implying the observed composition would change from light to heavy close to the cut-off if one or two nearby AGN fossils dominate. Given that several percent of the universe's volume may house such slowly decaying structures, these fossils may even re-energize ultra-high energy cosmic rays from distant/old sources, offsetting the ``GZK-losses'' due to interactions with photons of the cosmic microwave background radiation and giving evidence of otherwise undetectable fossils. In this case the composition would remain light to the highest energies if distant sources or fossils dominated, but otherwise would be mixed. It is hoped the new generation of cosmic ray experiments such as the Pierre Auger Observatory and ultra-high energy neutrino telescopes such as ANITA and lunar Cherenkov experiments will clarify this.

The exciting results from H.E.S.S. point to a new population of gamma-ray sources at energies E > 10 TeV, paving the way for future studies and new discoveries in the multi-TeV energy range. Connected with these energies is the search for sources of PeV cosmic-rays (CRs) and the study of multi-TeV gamma-ray production in a growing number of astrophysical environments. TenTen is a proposed stereoscopic array (with a suggested site in Australia) of modest-sized (10 to 30m^2) Cherenkov imaging telescopes with a wide field of view (8 to 10deg diameter) optimised for the E~10 to 100 TeV range. TenTen will achieve an effective area of ~10 km^2 at energies above 10 TeV. We outline here the motivation for TenTen and summarise key performance parameters.

The lunar Cherenkov technique, which aims to detect the coherent Cherenkov radiation produced when UHE particles interact in the lunar regolith, was first attempted with the Parkes radio-telescope in 1995, though the theory was not sufficiently developed at this time to calculate a limit on the UHE neutrino flux from the non-observation. Since then, the technique has evolved to include experiments utilising lower frequencies, wider bandwidths, and entire arrays of antenna. We develop a simulation to analyse the full range of experiments, and calculate the UHE neutrino flux limit from the Parkes experiment, including the directional dependence. Our results suggest a methodology for planning future observations, and demonstrate how to utilise all available information on the nature of radio pulses from the Moon for the detection of UHE particles.

The first search for ultra-high energy (UHE) neutrinos using a radio telescope was conducted by Hankins, Ekers and O'Sullivan (1996). This was a search for nanosecond duration radio Cherenkov pulses from electromagnetic cascades initiated by ultra-high energy (UHE) neutrino interactions in the lunar regolith, and was made using a broad-bandwidth receiver fitted to the Parkes radio telescope, Australia. At the time, no simulations were available to convert the null result into a neutrino flux limit. Since then, similar experiments at Goldstone, USA, and Kalyazin, Russia, have also recorded null results, and computer simulations have been used to model the experimental sensitivities of these two experiments and put useful limits on the UHE neutrino flux. Proposed future experiments include the use of broad-bandwidth receivers, making the sensitivity achieved by the Parkes experiment highly relevant to the future prospects of this field. We have therefore calculated the effective aperture for the Parkes experiment and found that when pointing at the lunar limb, the effective aperture at all neutrino energies was superior to single-antenna, narrow-bandwidth experiments, and that the detection threshold was comparable to that of the double-antenna experiment at Goldstone. However, because only a small fraction of the observing time was spent pointing the limb, the Parkes experiment places only comparatively weak limits on the UHE neutrino flux. Future efforts should use multiple telescopes and broad-bandwidth receivers.

The HESS instrument has observed a diffuse flux of ~ TeV gamma rays from a large solid angle around the Galactic center (GC). This emission is correlated with the distribution of gas in the region suggesting that the gamma rays originate in collisions between cosmic ray hadrons (CRHs) and ambient matter. Of particular interest, HESS has detected gamma rays from the Sagittarius (Sgr) B Molecular Cloud Complex. Prompted by the suggestion of a hadronic origin for the gamma rays, we have examined archival 330 and 74 MHz Very Large Array radio data and 843 MHz Sydney University Molonglo Sky Survey data covering Sgr B, looking for synchrotron emission from secondary electrons and positrons (expected to be created in the same interactions that supply the observed gamma rays). Intriguingly, we have uncovered non-thermal emission, but at a level exceeding expectation. Adding to the overall picture, recent observations by the Atacama Pathfinder Experiment telescope show that the cosmic ray ionization rate is ten times greater in the Sgr B2 region of Sgr B than the local value. Lastly, Sgr B2 is also a very bright X-ray source. We examine scenarios for the spectra of CRHs and/or primary electrons that would reconcile all these different data. We determine that (i) a hard (~ E^-2.2), high-energy (> TeV) population CRHs is unavoidably required by the HESS gamma ray data and (ii) the remaining broad-band, non-thermal phenomenology is explained either by a rather steep (~ E^-2.9) spectrum of primary electrons or a (~ E^-2.7) population of CRHs. No single, power-law population of either leptons or hadrons can explain the totality of broadband, non-thermal Sgr B phenomenology.

We present calculations of the spectral and temporal radiative signatures expected from ultrarelativistic protons in compact sources. The coupling between the protons and the leptonic component is assumed to occur via Bethe-Heitler pair production. This process is treated by modeling the results of Monte-Carlo simulations and incorporating them in a time-dependent kinetic equation, that we subsequently solve numerically. Thus, the present work is, in many respects, an extension of the leptonic `one-zone' models to include hadrons. Several examples of astrophysical importance are presented, such as the signature resulting from the cooling of relativistic protons on an external black-body field and that of their cooling in the presence of radiation from injected electrons. We also investigate and refine the threshold conditions for the 'Pair Production/Synchrotron' feedback loop which operates when relativistic protons cool efficiently on the synchrotron radiation of the internally produced Bethe-Heitler pairs. We demonstrate that an additional component of injected electrons lowers the threshold for this instability.

The origin and nature of the highest energy cosmic ray events is currently the subject of intense investigation by giant air shower arrays and fluorescent detectors. These particles reach energies well beyond what can be achieved in ground-based particle accelerators and hence they are fundamental probes for particle physics as well as astrophysics. Because of the scarcity of these high-energy particles, larger and larger ground-based detectors have been built. The new generation of digital radio telescopes may play an important role in this, if properly designed. Radio detection of cosmic ray showers has a long history but was abandoned in the 1970's. Recent experimental developments together with sophisticated air shower simulations incorporating radio emission give a clearer understanding of the relationship between the air shower parameters and the radio signal, and have led to resurgence in its use. Observations of air showers by the SKA could, because of its large collecting area, contribute significantly to measuring the cosmic ray spectrum at the highest energies. Because of the large surface area of the moon, and the expected excellent angular resolution of the SKA, using the SKA to detect radio Cherenkov emission from neutrino-induced cascades in lunar regolith will be potentially the most important technique for investigating cosmic ray origin at energies above the photoproduction cut-off. (abridged)

The giant radio galaxy M87 is usually classified as a Fanaroff-Riley class I source, suggesting that M87 is a mis-aligned BL Lac object. Its unresolved nuclear region emits strong non-thermal emission from radio to X-rays which has been interpreted as synchrotron radiation. In an earlier paper we predicted M87 as a source of detectable gamma ray emission in the context of the hadronic Synchrotron-Proton Blazar (SPB) model. The subsequent tentative detection of TeV energy photons by the HEGRA-telescope array would, if confirmed, make it the first radio galaxy to be detected at TeV-energies. We discuss the emission from the unresolved nuclear region of M87 in the context of the SPB model, and give examples of possible model representations of its non-simultaneous spectral energy distribution. The low-energy component can be explained as synchrotron radiation by a primary relativistic electron population that is injected together with energetic protons into a highly magnetized emission region. We find that the gamma-ray power output is dominated either by mu/pi-synchrotron or proton synchrotron radiation depending on whether the primary electron synchrotron component peaks at low or high energies, respectively. The predicted gamma-ray luminosity peaks at around 100 GeV at a level comparable to that of the low-energy hump, and this makes M87 a promising candidate source for the newly-commissioned high-sensitivity low-threshold Cherenkov telescopes H.E.S.S., VERITAS, MAGIC and CANGAROO III. (abridged)

I discuss the shape of the high energy end of the spectrum of particles arising from diffusive shock acceleration in the presence of (i) additional diffusive escape from the accelerator, (ii) continuous energy losses, (iii) energy changes arising from interactions. The form of the spectrum near cut-off is sensitive to these processes as well as to the momentum-dependence of the diffusion coefficients and the compression ratio, and so the spectrum of any radiation emitted by the accelerated particles may reflect the physical conditions of the acceleration region. Results presented in this paper have applications in interpreting the spectral energy distributions of many types of astrophysical object including supernova remnants (SNR), active galactic nuclei (AGN) and acceleration sources of ultra high energy cosmic rays (UHE CR). Except for extremely nearby sources, spectral features imprinted on the spectrum of UHE CR during the acceleration process will be largely eroded during propagation, but the spectrum of UHE neutrinos produced in interactions of UHE CR with radiation, both during cosmic ray acceleration and subsequent propagation through the cosmic microwave background radiation, contains sufficient information to determine the cut-off momentum of the UHE CR just after acceleration for reasonable assumptions. Observation of these UHE neutrinos by the Pierre Auger Observatory may help in identifying the sources of the highest energy cosmic rays.

Cosmic rays with energies above $10^{18}$ eV are currently of considerable interest in astrophysics and are to be further studied in a number of projects which are either currently under construction or the subject of well-developed proposals. This paper aims to discuss some of the physics of such particles in terms of current knowledge and information from particle astrophysics at other energies.

The existence of intrinsic obscuration of Fanaroff-Riley I objects is a controversial topic. M87, the nearest such object, is puzzling in that although it has very massive central black hole it has a relatively low luminosity, suggesting it is in a dormant state. Despite of its proximity to us (16 Mpc) it is not known with certainty whether or not M87 has a dusty torus. Infrared observations indicate that if a torus exists in M87 it must have a rather low infrared luminosity. Using arguments from unification theory of active galactic nuclei, we have earlier suggested that the inner parsec-scale region of M87 could still harbour a small torus sufficiently cold such that its infrared emission is dwarfed by the jet emission. The infrared emission from even a small cold torus could affect through photon-photon pair production interactions the escape of 100 GeV to TeV energy gamma rays from the central parsec of M87. The TeV gamma-ray flux from the inner jet of M87 has recently been predicted in the context of the Synchrotron Proton Blazar (SPB) model to extend up to at least 100 GeV (Protheroe, Donea, Reimer, 2002). Subsequently, the detection of gamma-rays above 730 GeV by the HEGRA Cherenkov telescopes has been reported. We discuss the interactions of gamma-rays produced in the inner jet of M87 with the weak infrared radiation expected from a possible dusty small-scale torus, and show that the HEGRA detection shows that the temperature of any torus surrounding the gamma-ray emission region must be cooler than about 250 K. We suggest that if no gamma-rays are in future detected during extreme flaring activity in M87 at other wavelength, this may be expected because of torus heating.

Motivated by unification schemes of active galactic nuclei, we review evidence for the existence of a small-scale dust torus in M87, a Fanaroff-Riley Class I radio galaxy. Since there is no direct evidence of any thermal emission from its torus we consider indirect evidence, such as BLR activity and ageing arguments to model the cold dust structure of M87. In the context of the jet -- accretion disk -- torus symbiosis we discuss the interactions of GeV and TeV gamma-rays produced in the jet of M87 with the infrared radiation fields external to the jet, produced by a less active torus. A thin and cold torus with less defined outer boundaries could still posses problems to some of the TeV emission from the jet.

I re-examine the brightness temperature problem in PKS 0405-385 which is an extreme intra-day variable radio quasar with an inferred brightness temperature of $\sim 5 \times 10^{14}$ K at 5 GHz, well above the Compton catastrophe limit of $\sim 10^{11}$ K reached when the synchrotron photon energy density exceeds the energy density of the magnetic field. If one takes into account the uncertainty in the distance to the ionized clouds responsible for interstellar scintillation causing rapid intra-day variability in PKS 0405-385 it is possible that the brightness temperature could be as low as $\sim 10^{13}$ K at 5 GHz, or even lower. The radio spectrum can be fitted by optically thin emission from mono-energetic electrons, or an electron spectrum with a low-energy cut-off such that the critical frequency of the lowest energy electrons is above the radio frequencies of interest. If one observes optically thin emission along a long narrow emission region, the average energy density in the emission region can be many orders of magnitude lower than calculated from the observed intensity if one assumed a spherical emission region. I discuss the physical conditions in the emission region and find that the Compton catastrophe can then be avoided using a reasonable Doppler factor. I also show that MeV to 100 GeV gamma-ray emission at observable flux levels should be expected from extreme intra-day variable sources such as PKS 0405-385.

The unresolved nuclear region of M87 emits strong non-thermal emission from radio to X-rays. Assuming this emission to originate in the pc scale jet aligned at $\theta \sim 30^\circ$ to the line of sight, we interpret this emission in the context of the Synchrotron Proton Blazar (SPB) model. We find the observed nuclear jet emission to be consistent with M87 being a mis-aligned BL Lac Object and predict gamma-ray emission extending up to at least 100 GeV at a level easily detectable by GLAST and MAGIC, and possibly by VERITAS depending on whether it is high-frequency or low-frequency peaked. Predicted neutrino emission is below the sensitivity of existing and planned neutrino telescopes. Ultra-high energy neutrons produced in pion photoproduction interactions decay into protons after escaping from the host galaxy. Because energetic protons are deflected by the intergalactic magnetic field, the protons from the decay of neutrons emitted in all directions, including along the jet axis where the Doppler factor and hence emitted neutron energies are higher, can contribute to the observed ultra-high energy cosmic rays. We consider the propagation of these cosmic ray protons to Earth and conclude that M87 could account for the observed flux if the extragalactic magnetic field topology were favourable.

We investigate ways of accurately simulating the propagation of energetic charged particles over small times where the standard Monte Carlo approximation to diffusive transport breaks down. We find that a small-angle scattering procedure with appropriately chosen step-lengths and scattering angles gives accurate results, and we apply this to the simulation of propagation upstream in relativistic shock acceleration.

The relationship between observed variability time and emission region geometry is explored for the case of emission by relativistic jets. The approximate formula for the jet-frame size of the emission region, $R'=Dc\Delta t_{\rm obs}$ is shown to lead to large systematic errors when used together with observed luminosity and assumed or estimated Doppler factor $D$ to estimate the jet-frame photon energy density. These results have implications for AGN models in which low-energy photons are targets for interaction of high energy particles and photons, e.g. synchrotron-self Compton models and hadronic blazar models, as well as models of intra-day variable sources in which the photon energy density imposes a brightness temperature limit through Compton scattering. The actual relationship between emission region geometry and observed variability is discussed for a variety of geometries including cylinders, spheroids, bent, helical and conical jet structures, and intrinsic variability models including shock excitation. The effects of time delays due to finite particle acceleration and radiation time scales are also discussed.

We calculate the spectral energy distribution (SED) of electromagnetic radiation and the spectrum of high energy neutrinos from BL Lac objects in the context of the Synchrotron Proton Blazar Model. In this model, the high energy hump of the SED is due to accelerated protons, while most of the low energy hump is due to synchrotron radiation by co-accelerated electrons. To accelerate protons to sufficiently high energies to produce the high energy hump, rather high magnetic fields are required. Assuming reasonable emission region volumes and Doppler factors, we then find that in low-frequency peaked BL Lacs (LBLs), which have higher luminosities than high-frequency peaked BL Lacs (HBLs), there is a significant contribution to the high frequency hump of the SED from pion photoproduction and subsequent cascading, including synchrotron radiation by muons. In contrast, in HBLs we find that the high frequency hump of the SED is dominated by proton synchrotron radiation. We are able to model the SED of typical LBLs and HBLs, and to model the famous 1997 flare of Markarian 501. We also calculate the expected neutrino output of typical BL Lac objects, and estimate the diffuse neutrino intensity due to all BL Lacs. Because pion photoproduction is inefficient in HBLs, as protons lose energy predominantly by synchrotron radiation, the contribution of LBLs dominates the diffuse neutrino intensity. We suggest that nearby LBLs may well be observable with future high-sensitivity TeV gamma-ray telescopes.

The radiation fields external to the jets and originating from within a few parsecs from the black hole, are discussed in this paper. They are the direct radiation from an accretion disk in symbiosis with jets, the radiation field from the broad line region (BLR) surrounding the accretion disk, and the infrared radiation from a dusty torus. The jet/disk symbiosis modifies the energetics in the central parsec of AGN such that for a given accretion rate, a powerful jet would occur at the expense of the disk luminosity, and consequently the disk would less efficiently ionize the BRL clouds or heat the dust in the torus, thereby affecting potentially important target photon fields for interactions of gamma-rays, accelerated electrons and protons along the jet. Motivated by unification schemes of active galactic nuclei, we briefly review the evidence for the existence of broad line regions and small-scale dust tori in BL Lacs and Fanaroff-Riley Class I (FR-I) radio galaxies. We propose that an existing jet-accretion disk symbiosis can be extrapolated to provide a large scale-symbiosis between other important dusty constituents of the blazar/FR-I family. In the present paper, we discuss in the context of this symbiosis interactions of GeV and TeV gamma-rays produced in the jet with the various radiation fields external to the jet in quasars and blazars, taking account the anisotropy of the radiation.

We report on the status of the construction of an array of four 10 m atmospheric Cherenkov telescopes for gamma-ray astronomy, near Woomera, in South Australia -- the CANGAROO-III project. The first telescope of this array is the upgraded version of the CANGAROO-II 7 m telescope and has been in operation since March 2000. The second telescope, an improved version of the first, is being constructed for installation in late 2001. Stereoscopic observation of sub TeV gamma-rays with the two 10 m telescopes will begin in 2002 and the full array will be operational in 2004.

The Synchrotron Proton Blazar model is a promising model to explain high energy emission from gamma-ray loud BL Lac objects like Mkn 421. In contrast to leptonic models, the hadronic explanation of gamma-ray emission predicts ultrahigh energy neutrinos. The predicted neutrino spectra from a typical High-energy cutoff BL Lac Object (HBL) and a Low-energy cutoff BL Lac Object (LBL) are presented. We find that cooling due to muon synchrotron radiation causes a cutoff of the neutrino spectrum at $\sim 10^{18}$ eV, with the exception of $\nu_\mu$ from kaon decay which may extend to higher energies if meson production takes place in the secondary resonance region of the cross section. The impact of the neutrino output from both source populations to the diffuse neutrino background is discussed.

We consider the galactic population of gamma-ray pulsars as possible sources of cosmic rays at and just above the ``knee'' in the observed cosmic ray spectrum at $10^{15}$--$10^{16}$ eV. We suggest that iron nuclei may be accelerated in the outer gaps of pulsars, and then suffer partial photo-disintegration in the non-thermal radiation fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei are injected into the expanding supernova remnant. We compute the spectra of nuclei escaping from supernova remnants into the interstellar medium, taking into account the observed population of radio pulsars. Our calculations, which include a realistic model for acceleration and propagation of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the observed cosmic rays in the knee region. Such a contribution of heavy nuclei to the cosmic ray spectrum at the knee can significantly increase the average value of $<\ln A>$ with increasing energy as is suggested by recent observations.

We apply the synchrotron proton blazar (SPB) model to the April 1997 flare of Markarian 501 and find we are able to fit the observed spectral energy distribution. We explore the effect of target photon density on the high energy part of the spectral energy distribution (SED) for fixed assumed magnetic field, emission region size and Doppler factor and find that the luminosity and peak frequency of the high energy part of the SED may depend on the luminosity of the low energy part of the SED in high-frequency peaked BL Lac objects (HBL).

I give a brief review of high energy gamma-ray signatures of dark matter. The decay of massive $X$-particles and subsequent hadronization have been suggested as the origin of the highest energy cosmic rays. Propagation over cosmological distances to Earth (as would be the case in some topological defect origin models for the $X$-particles) results in potentially observable gamma-ray fluxes at GeV energies. Massive relic particles on the other hand, would cluster in galaxy halos, including that of our Galaxy, and may give rise to anisotropic gamma ray and cosmic ray signals at ultra high energies. Future observations above 100 Gev of gamma rays due to WIMP annihilation in the halo of the Galaxy may be used to place constraints on supersymmetry parameter space.

We discuss the various contributions to the jet luminosity in proton blazar models of active galactic nuclei and describe a method of estimating the jet luminosity from the observed spectral energy distribution (SED) and the fitted model parameters. We apply this to a synchrotron proton blazar (SPB) model for Markarian 501.

We consider the implications of the recent determination of the universal infrared background for the propagation of photons up to 20 TeV from the active galaxy Markarian 501 as observed by HEGRA. At 20 TeV the mean free path for photon-photon collisions on the infrared background would be much shorter than the distance to Markarian 501, implying absorption factors of the order of exp(-10), or greater, and consequently an excessive power output for this active galaxy. Possible solutions of this problem are discussed.

We present a new calculation of the propagation of protons with energies above $10^{19}$ eV over distances of up to several hundred Mpc. The calculation is based on a Monte Carlo approach using the event generator SOPHIA for the simulation of hadronic nucleon-photon interactions and a realistic integration of the particle trajectories in a random extragalactic magnetic field. Accounting for the proton scattering in the magnetic field affects noticeably the nucleon energy as a function of the distance to their source and allows us to give realistic predictions on arrival energy, time delay, and arrival angle distributions and correlations as well as secondary particle production spectra.

(abr.) The spectral energy distribution (SED) of blazars typically has a double-humped appearance usually interpreted in terms of synchrotron self-Compton models. In proton blazar models, the SED is instead explained in terms of acceleration of protons and subsequent cascading. We discuss a variation of the Synchrotron Proton Blazar model, first proposed by M\"ucke & Protheroe (1999), in which the low energy part of the SED is mainly synchrotron radiation by electrons co-accelerated with protons which produce the high energy part of the SED mainly asproton synchrotron radiation. Using a Monte Carlo/numerical technique to simulate the interactions and subsequent cascading of the accelerated protons, we are able to fit the observed SED of Markarian 501 during the April 1997 flare. We find that the emerging cascade spectra initiated by gamma-rays from $\pi^0$ decay and by $e^\pm$ from $\mu^\pm$ decay turn out to be relatively featureless. Synchrotron radiation produced by $\mu^\pm$ from $\pi^\pm$ decay, and even more importantly by protons, and subsequent synchrotron-pair cascading, is able to reproduce well the high energy part of the SED. For this fit we find that synchrotron radiation by protons dominates the TeV emission, pion photoproduction being less important with the consequence that we predict a lower neutrino flux than in other proton blazar models.

We obtain the maximum diffuse neutrino intensity predicted by hadronic photoproduction models of the type which have been applied to the jets of active galactic nuclei or gamma ray bursts. For this, we compare the proton and gamma ray fluxes associated with hadronic photoproduction in extragalactic neutrino sources with the present experimental upper limit on cosmic ray protons and the extragalactic gamma ray background, employing a transport calculation of energetic protons traversing cosmic photon backgrounds. We take into account the effects of the photon spectral shape in the sources on the photoproduction process, cosmological source evolution, the optical depth for cosmic ray ejection, and discuss the possible effects of magnetic fields in the vicinity of the sources. For photohadronic neutrino sources which are optically thin to the emission of neutrons we find that the cosmic ray flux imposes a stronger bound than the extragalactic gamma ray background in the energy range between 10^5 GeV and 10^11 GeV, as previously noted by Waxman & Bahcall (1999). We also determine the maximum contribution from the jets of active galactic nuclei, using constraints set to their neutron opacity by gamma-ray observations. This present upper limit is consistent with the jets of active galactic nuclei producing the extragalactic gamma ray background hadronically, but we point out future observations in the GeV-to-TeV regime could lower this limit. We also briefly discuss the contribution of gamma ray bursts to ultra-high energy cosmic rays as it can be inferred from possible observations or limits on their correlated neutrino fluxes.

The April 1997 giant flare of Mkn 501 is modelled by a stationary Synchrotron-Proton-Blazar model. Our derived model parameters are consistent with X-ray-to-TeV-data in the flare state and diffusive shock acceleration of electron and proton in a Kolmogorov/Kraichnan turbulence sprectrum. While the emerging pair-synchrotron cascade spectra initiated by photons from neutral pi-decay and electrons from pi+- --> mu+- --> e+- decay turn out to be relatively featureless, muon and proton synchrotron radiation and their cascade radiation produce a double-humped spectral energy distribution. For the present model we find proton synchrotron radiation to dominate the TeV emission, while the contribution from the synchrotron radiation of the pairs, produced by the high energy hump, is only minor.

Some gamma ray bursts may be produced by supernovae exploding in close massive binary systems (type Ib/c supernovae) as suggested by the recent observation of SN 1998bw/GRB 980425. We propose that high energy radiation observed in such gamma ray bursts may be produced by synchrotron radiation of electrons accelerated by 1st order Fermi acceleration at a quasi-stationary shock in the high velocity SN ejecta colliding with the companion star or some other nearby massive object. Nuclei would also be accelerated, and could give rise to an observable fluence of high energy neutrinos at Earth. This paper has been withdrawn by the authors following evidence that SN 1998bw and GRB 980425 are two distinct sources.

The recent detection of gamma radiation from Mkn 501 at energies as high as 25 TeV suggests stringent upper bounds on the diffuse, far infrared, extragalactic radiation density. The production of electron-positron pairs through photon-photon collisions would prevent gamma photons of substantially higher energies from reaching us across distances of order 100 Mpc. However, coherently arriving TeV or sub-TeV gammas - Bose-Einstein condensations of photons at these energies - could mimic the Cherenkov shower signatures of extremely energetic gammas. To better understand such events, we describe their observational traits and discuss how they might be generated.

In a recent paper, Waxman and Bahcall (PRD 59, 023002) claimed that the present data on ultra-high energy cosmic rays imply a model-independent upper bound on extragalactic neutrino fluxes of 2*10^-8 GeV cm^-2 s^-1 ster^-1, independent of neutrino energy. Mannheim, Protheroe and Rachen (astro-ph/9812398) have repeated this calculation and confirmed the WB-bound, within a factor of 2, only at a very limited energy range of E_\nu \sim 10^(16-18) eV, while at other energies the neutrino flux is mainly limited by the extragalactic gamma ray background to a level about two orders of magnitude higher than the WB bound. In this paper we present a simple, (almost) no-math approach to the problem, and discuss under which astrophysical assumptions the WB-bound and the MPR-bound, respectively, apply. Then we discuss to which respect these assumptions apply to presently discussed models of extragalactic neutrino production. We note that, averaged over the observed luminosity function, blazars are sufficiently opaque to ultra-high energy neutrons that there is no conflict of the predicted neutrino fluxes with the cosmic ray data, and that these models are rather constrained by their contributions to the extragalactic gamma ray background. At present, no modifications are implied to the predicted neutrino events from these models in active or planned neutrino detectors.

Nature produces cosmic ray particles, probably protons, with energies well above $10^{20}$ eV -- how are they produced? Where do they come from? Gamma rays with energies above $10^{13}$ eV are produced in jets of active galaxies -- are these produced by energetic electrons or protons? What is the correct model of Gamma Ray Bursts? These are just some of the fundamental questions in high energy astrophysics to be answered by observations made with large area neutrino telescopes.

We determine the physical parameters (magnetic field and Doppler factor) of the homogeneous synchrotron self-Compton model allowed by the observed X-ray to gamma-ray spectra and variability of Markarian~501 during the 15-16 April 1997 flaring activity. We find that magnetic fields between 0.07 G and 0.6 G and Doppler factors between 12 and 36 could fit (depending on observed variability time scale) these observations. We take account of photon-photon pair production interactions of gamma-ray photons occurring both inside the emission region and during propagation to Earth and find these to be extremely important in correctly determining the allowed model parameters. Previous estimates of the allowed parameter space have neglected this effect. Future multi-wavelength campaigns during strong flaring activity, including observations from optical to TeV gamma-rays, should enable the physical parameters to be further constrained.

Photomeson production is the main energy loss for relativistic nucleons in dense radiation fields like the cosmic microwave background and the radiation fields in Gamma Ray Bursts (GRB) and jets of Active Galactic Nuclei (AGN). In this paper we study photomeson production in typical GRB and AGN jet radiation fields by using the recently developed Monte Carlo event generator SOPHIA (see these proceedings). We discuss processes that are relevant for the physics of cosmic ray acceleration and the production of neutrinos and gamma rays. We compare our results with widely used approximations, and find significant deviations, particularly for GRBs. The photoproduction of antibaryons as a so far not considered effect in astrophysics is briefly discussed.

A new Monte Carlo program for photohadronic interactions of relativistic nucleons with an ambient photon radiation field is presented. The event generator is designed to fulfil typical astrophysical requirements, but can also be used for radiation and background studies at high energy colliders such as LEP2 and HERA, as well as for simulations of photon induced air showers. We consider the full photopion production cross section from the pion production threshold up to high energies. It includes resonance excitation and decay, direct single pion production and diffractive and non-diffractive multiparticle production. The cross section of each individual process is calculated by fitting experimental data, while the kinematics is determined by the underlying particle production process. We demonstrate that our model is capable of reproducing known accelerator data over a wide energy range.

Radiative deceleration of relativistic jets in active galactic nuclei as the result of inverse Compton scattering of soft photons from accretion discs is discussed. The Klein-Nishina (KN) cross section is used in the calculation of the radiation force due to inverse Compton scattering. Our result shows that deceleration due to scattering in the KN regime is important only for jets starting with a bulk Lorentz factor larger than 1000. When the bulk Lorentz factor satisfies this condition, particles scattering in the Thomson regime contribute positively to the radiation force (acceleration), but those particles scattering in the KN regime are dominant and the overall effect is deceleration. In the KN limit, the drag due to Compton scattering, though less severe than in the Thomson limit, strongly constrains the bulk Lorentz factor. Most of the power from the deceleration goes into radiation and hence the ability of the jet to transport significant power (in particle kinetic energy) out of the subparsec region is severely limited. The deceleration efficiency decreases significantly if the jet contains protons and the proton to electron number density ratio satisfies the condition $n_p/n_{e0}>2\gamma_{\rm min}/\mu_p$ where $\gamma_{\rm min}$ is the minimum Lorentz factor of relativistic electrons (or positrons) in the jet frame and $\mu_p$ is the proton to electron mass ratio.

In this chapter I give an overview of shock acceleration, including a discussion of the maximum energies possible and the shape of the spectrum near cut-off, interactions of high energy cosmic rays with, and propagation through, the background radiation, and the resulting electron-photon cascade. Possible sources of the highest energy cosmic rays are discussed including active galaxies, gamma ray bursts and topological defects. I argue that while the origin of the highest energy cosmic rays is still uncertain, it is not necessary to invoke exotic models such as emission by topological defects to explain the existing data. It seems likely that shock acceleration at Fanaroff-Riley Class II radio galaxies can account for the existing data. However, new cosmic ray data, as well as better estimates of the extragalactic radiation fields and magnetic fields will be necessary before we will be certain of the origin of the highest energy particles occurring in nature.

I give a brief discussion of possible sources of high energy neutrinos of astrophysical origin over the energy range from $\sim 10^{12}$ eV to $\sim 10^{25}$ eV. In particular I shall review predictions of the diffuse neutrino intensity. Neutrinos from interactions of galactic cosmic rays with interstellar matter are guaranteed, and the intensity can be reliably predicted to within a factor of 2. Somewhat less certain are intensities in the same energy range from cosmic rays escaping from normal galaxies or active galactic nuclei (AGN) and interacting with intracluster gas. At higher energies, neutrinos will definitely be produced by interactions of extragalactic cosmic rays with the microwave background. With the discovery that gamma ray bursts (GRB) are extragalactic, and therefore probably the most energetic phenomena in the Universe, it seems likely that they will be copious sources of high energy neutrinos. Other sources, such as AGN and topological defects, are more speculative. However, searches for neutrinos from all of these potential sources should be made because their detection would have important implications for high energy astrophysics and cosmology.

We discuss the first applications of our newly developed Monte Carlo event generator SOPHIA to multiparticle photoproduction of relativistic protons with thermal and power law radiation fields. The measured total cross section is reproduced in terms of excitation and decay of baryon resonances, direct pion production, diffractive scattering, and non-diffractive multiparticle production. Non--diffractive multiparticle production is described using a string fragmentation model. We demonstrate that the widely used `$\Delta$--approximation' for the photoproduction cross section is reasonable only for a restricted set of astrophysical applications. The relevance of this result for cosmic ray propagation through the microwave background and hadronic models of active galactic nuclei and gamma-ray bursts is briefly discussed.

We have examined cutoffs and pile-ups due to various processes in the spectra of particles produced by shock acceleration, and found that, even in the absence of energy losses, the shape of the spectrum of accelerated particles at energies well below the nominal maximum energy depends strongly on the energy dependence of the diffusion coefficient. This has implications in many areas, for example, in fitting the observed cosmic ray spectrum with models based on power-law source spectra and rigidity dependent diffusive escape from the galaxy. With continuous energy losses, prominent pile-ups may arise, and these should be included when modelling synchrotron X-ray and inverse Compton gamma-ray spectra from a shock-accelerated electron population. We have developed a Monte Carlo/numerical technique to model the shape of the spectrum for the case of non-continuous energy losses such as inverse Compton scattering in the Klein-Nishina regime. We find that the shapes of the resulting cut-offs differ substantially from those arising from continuous processes, and we suggest that such differences could be observable through their effect on the spectrum of radiation emitted by a population of recently accelerated electrons as, for example, may exist in young supernova remnants.

We discuss the consequences of almost rectilinear acceleration of protons to extremely high energies in a reconnection region on the surface of an accretion disk which surrounds a central black hole in an active galaxy. The protons produce $\gamma$-rays and neutrinos in interactions with the disk radiation as considered in several previous papers. However, in this model the secondary $\gamma$-rays can initiate cascades in the magnetic and radiation fields above the disk. We compute the spectra of $\gamma$-rays and neutrinos emerging from regions close to the disk surface. Depending on the parameters of the reconnection regions, this model predicts the appearance of $\gamma$-ray and neutrino flares if protons takes most of the energy from the reconnection region. In contrast, if leptons take most of the energy, they produce pure $\gamma$-ray flares. The $\gamma$-ray spectrum expected in the case of hadronic cascading is compared with the spectrum observed during the flare in June 1991 from 3C 279. The neutrino flares which should accompany these gamma-ray flares may be detected by future large scale neutrino telescopes sensitive at $\sim 10^{5}$ TeV.

It is shown that for moderately hot polar caps (with effective temperature of $\sim 10^6$ K), the efficiency of polar gap acceleration is lower compared to the case in which the polar caps are relatively cool and inverse Compton scattering plays no role in controlling the gap. For young pulsars with superstrong magnetic fields ($\geq10^9 T$) and hot polar caps (with temperature of $\ge 5\times10^6 K$), because of the energy loss of electrons or positrons due to resonant inverse Compton scattering in the vicinity of polar caps, pair cascades occur at distances further away from the polar cap, and in this case, we have a relatively high acceleration efficiency with ions carrying most of the particle luminosity.

During remarkable flaring activity in 1997 Markarian 501 was the brightest source in the sky at TeV energies, outshining the Crab Nebula by a factor of up to 10. Periods of flaring activity each lasting a few days were observed simultaneously by several gamma ray telescopes. Contemporaneous multiwavelength observations in April 1997 show that a substantial fraction of the total AGN power is at TeV energies indicating that high energy processes dominate the energetics of this object. Rapid variability, on time scales of less than a day, and a flat spectrum extending up to at least 10 TeV characterize this object. Results of 1997 observations by 6 telescopes are summarized and some of the implications of these results are discussed.

We consider the result of acceleration of heavy ions in the slot gap potential of a very young pulsar with a hot polar cap. Photodisintegration of the heavy ions in the radiation field of the polar cap and pulsar surface gives rise to a flux of energetic neutrons. Some fraction of these neutrons interact with target nuclei in the supernova shell to produce neutrino and gamma-ray signals which should be observable from very young supernova remnants in our galaxy for a range of pulsar parameters.

The bulk of the cosmic rays up to about 100 TeV are thought to be accelerated by the 1st order Fermi mechanism at supernova shocks, producing a power-law spectrum. Both electrons and protons should be accelerated, but their ratio on acceleration is not well known. Recently, the EGRET instrument on the Compton Gamma Ray Observatory has observed supernova remnants IC 443 and $\gamma$ Cygni at GeV energies. On the assumption that the observed gamma-rays are produced by accelerated particles in the remnants (rather than, for example, from a central compact object) we model the contributions due to pion production, bremsstrahlung, and inverse Compton scattering on the cosmic microwave, diffuse galactic radiation, and locally produced radiation fields. In the case of the same spectral index for both electrons and nuclei, and a cut-off at 80 TeV, we find that a spectral index of accelerated particles close to 2.4, and a ratio of electrons to protons in the range 0.2 to 0.3, gives a good fit to the observed spectra. For lower cut-off energies flatter spectra are possible. We also investigate the case where the electron spectrum is steeper than that of nuclei. We discuss the implications of our results for observations at air shower energies, and for the propagation of cosmic rays.

We investigate the consequences of the acceleration of heavy nuclei (e.g. iron nuclei) by the Crab pulsar. Accelerated nuclei can photodisintegrate in collisions with soft photons produced in the pulsar's outer gap, injecting energetic neutrons which decay either inside or outside the Crab Nebula. The protons from neutron decay inside the nebula are trapped by the Crab Nebula magnetic field, and accumulate inside the nebula producing gamma-rays and neutrinos in collisions with the matter in the nebula. Neutrons decaying outside the Crab Nebula contribute to the Galactic cosmic rays. We compute the expected fluxes of gamma-rays and neutrinos, and find that our model could account for the observed emission at high energies and may be tested by searching for high energy neutrinos with future neutrino telescopes currently in the design stage.

Based on the detected variability time scales of X-ray and TeV gamma-ray emission, and the observed multiwavelength photon spectrum, of Mrk 421 we place constraints on the allowed parameter space (magnetic field and Doppler factor of the emission region) for the homogeneous synchrotron self-Compton model. The spectra calculated for the allowed parameters are marginally consistent with the available spectral information above $\sim 1$ TeV reported by the Whipple Observatory in the case of a 1 day flare time scale. However, for the recently reported very short duration flares varying on a time scale of 15 min, the calculated spectra are significantly steeper, suggesting that the homogeneous synchrotron self Compton model has problems in describing the relatively flat observed spectra extending above a few TeV. We determine the maximum ratio of TeV gamma-ray luminosity to X-ray luminosity during flaring which is allowed by the homogeneous synchrotron self-Compton model for the case of no significant photon-photon absorption in the source.

The bulk of the diffuse galactic gamma-ray emission above a few tens of GeV has been conventionally ascribed to the decay of neutral pions produced in cosmic-ray interactions with interstellar matter. Cosmic-ray electrons may, however, make a significant contribution to the gamma-ray spectrum at high energies, and even dominate at TeV-PeV energies depending on their injection spectral index and acceleration cut-off energy. If the injection spectrum is flat, the highest energy electrons will also contribute a diffuse hard X-ray/soft gamma-ray flux via synchrotron emission, and this may offer an explanation for the OSSE observation of a steep spectrum below a few MeV from the inner Galaxy. We perform a propagation calculation for cosmic-ray electrons, and use the resulting interstellar electron spectrum to obtain the gamma-ray spectrum due to inverse Compton, synchrotron and bremsstrahlung interactions consistently from MeV to PeV energies. We compare our results with available observations from satellite-borne telescopes, optical Cherenkov telescopes and air shower arrays and place constraints on the injection spectrum of cosmic-ray electrons. With future observations at TeV-PeV energies it should be possible to determine the average interstellar spectrum of cosmic-ray electrons, and hence estimate their spectrum on acceleration.

Using the propagation of ultra high energy nucleons, photons and electrons in the universal radiation backgrounds, we obtain limits on the luminosity of topological defects scenarios for the origin of the highest energy cosmic rays. The limits are set as function of the mass of the X-particles emitted by the cosmic strings or other defects, the cosmological evolution of the topological defects, and the strength of the extragalactic magnetic fields. The existing data on the cosmic ray spectrum and on the isotropic 100 MeV gamma-ray background limit significantly the parameter space in which topological defects can generate the flux of the highest energy cosmic rays.

In this lecture I give an overview of shock acceleration, interactions of high energy cosmic rays with, and propagation through, the background radiation, and the resulting electron-photon cascade. I argue that while the origin of the highest energy cosmic rays is still uncertain, it is not necessary to invoke exotic models such as emission by topological defects to explain the existing data. It seems likely that shock acceleration at Fanaroff-Riley Class II radio galaxies can account for the existing data. However, new cosmic ray data, as well as better estimates of the extragalactic radiation fields and magnetic fields will be necessary before we will be certain of the origin of the highest energy particles occurring in nature.

I give a brief critical review of the predicted intensity of diffuse high energy neutrinos of astrophysical origin over the energy range from 10^12 eV to 10^24 eV. Neutrinos from interactions of galactic cosmic rays with interstellar matter are guaranteed, and the intensity can be reliably predicted to within a factor of about 2 up to 10^17 eV. Somewhat less certain are intensities in the same energy range from cosmic rays escaping from normal galaxies or active galactic nuclei (AGN) and interacting with intracluster gas. At higher energies, neutrinos will be produced by interactions of extragalactic cosmic rays with the microwave background. Other sources, such as AGN, in particular blazars, and topological defects, are more speculative. However, searches for neutrinos from all of these potential sources should be made because their observation would have important implications for high energy astrophysics and cosmology.

We have developed a model for the variability of gamma ray emission in jets of active galactic nuclei in which the variability arises as a result of photon-photon pair production interactions with X-rays emitted by a hot spot in the inner part of the accretion disk. As the hot spot orbits around the central engine, the amount of absorption varies periodically. Our model may account for the observed variability of TeV emission from Markarian 421 and other blazars detected by the EGRET instrument on the Compton Gamma Ray Observatory, as well as correlated variations of X-ray and TeV gamma flux from Markarian 421. Quasi-periodic variations observed in Markarian 421 enable us to place an upper limit on the black hole mass of (2-30) times 10^8 solar masses.

We have developed a model for gamma ray emission in jets of active galactic nuclei in which particle acceleration takes place at a shock in the relativistic jet plasma due to a massive star in the central region of the host galaxy moving through the jet. The gamma rays are produced in a pair-Compton cascade in the radiation field of the star initiated by accelerated electrons. Our model may account for the observed GeV to TeV gamma ray spectrum and variability of Markarian 421 and other blazars detected by the EGRET instrument on the Compton Gamma Ray Observatory.

Gamma rays have been observed from two blazars at TeV energies. One of these, Markarian 421, has been observed also at GeV energies and has roughly equal luminosity per decade at GeV and TeV energies. Photon-photon pair production on the infrared background radiation is expected to prevent observation above about 1 TeV. However, the infrared background is not well known and it may be possible to observe the nearest blazars up to energies somewhat below about 100 TeV where absorption on the cosmic microwave background will give a sharp cut-off. Blazars are commonly believed to correspond to low power radio galaxies, seen down along a relativistic jet; as such they are all expected to have the nuclear activity encircled by a dusty molecular torus, which subtends an angle of 90 degrees or more in width as seen from the central source. Photon-photon pair production can also take place on the infrared radiation produced at the AGN by this molecular torus and surrounding outer disk. We calculate the optical depth for escaping gamma-rays produced near the central black hole and at various points along the jet axis for the case of blazars where the radiation is observed in a direction closely aligned with the jet. We find that the TeV emission site must be well above the top of the torus. For example, if the torus has an inner radius of 0.1 pc and an outer radius of 0.2 pc, then the emission site in Mrk 421 would have be at least 0.25 pc above the upper surface of the torus, and if Mrk 421 is observed above 50 TeV in the future, the emission site would have to be at least 0.5 pc above the upper surface. This has important implications for models of gamma-ray emission in active galactic nuclei.

I describe a model of energetic gamma ray and neutrino emission in blazars which is consistent with most of the gamma-ray observations, and use this model to estimate the diffuse intensity of high energy neutrinos from blazars.

We make a new estimate of the extragalactic radio background down to kHz frequencies based on the observed luminosity functions and radio spectra of normal galaxies and radio galaxies. We have constructed models for the spectra of these two classes of objects down to low frequencies based on observations of our Galaxy, other normal galaxies and radio galaxies. We check that the models and evolution of the luminosity functions give source counts consistent with data and calculate the radio background expected from kHz to GHz frequencies. The motivation for this calculation is that the propagation of ultra-high energy gamma-rays in the universe is limited by photon-photon pair production on the radio background. Electromagnetic cascades involving photon-photon pair production and subsequent synchrotron radiation in the intergalactic magnetic field may develop. Such gamma-rays may be produced in acceleration sites of ultra-high energy cosmic rays, as a result of interactions with the microwave background, or emitted as a result of decay or annihilation of topological defects. We find that photon-photon pair production on the radio background remains the dominant attenuation process for gamma-rays from $3 \times 10^{10}$ GeV up to GUT scale energies.

We use of a hybrid matrix--Monte Carlo method to simulate the cascade through the cosmic background radiation initiated by UHE particles and radiation emitted by topological defects. We follow the cascade over cosmological distances and calculate the intensities of hadrons, gamma-rays and neutrinos produced. We compare our results with the observed cosmic ray intensity at 300 EeV and lower energies, and conclude that topological defects are most unlikely to be the origin of the most energetic cosmic ray events.

We describe the results of a hybrid matrix--Monte Carlo calculation of cascading of UHE cosmic rays and $\gamma$-rays through the cosmic background radiation fields over cosmological distances. We calculate the $\gamma$-ray and neutrino emission that results from the cascade, as well as the effect of cascading on the primary spectrum. We discuss the results for various cosmic ray injection spectra and primary species. Certain models for the production of the highest energy cosmic rays are ruled out.

We describe a calculation of electromagnetic cascading in radiation and matter in the early universe initiated by the decay of massive particles or by some other process. We have used a combination of Monte Carlo and numerical techniques which enables us to use exact cross sections, where known, for all the relevant processes. In cascades initiated after the epoch of big bang nucleosynthesis $\gamma$-rays in the cascades will photodisintegrate $^4$He, producing $^3$He and deuterium. Using the observed $^3$He and deuterium abundances we are able to place constraints on the cascade energy deposition as a function of cosmic time. In the case of the decay of massive primordial particles, we place limits on the density of massive primordial particles as a function of their mean decay time, and on the expected intensity of decay neutrinos.

We consider the production of high energy neutrinos and cosmic rays in radio-quiet active galactic nuclei (AGN) or in the central regions of radio-loud AGN. We use a model in which acceleration of protons takes place at a shock in an accretion flow onto a supermassive black hole, and follow the cascade that results from interactions of the accelerated protons in the AGN environment. We use our results to estimate the diffuse high energy neutrino intensity and cosmic ray intensity due to AGN. We discuss our results in the context of high energy neutrino telescopes under construction, and measurements of the cosmic ray composition in the region of the ``knee'' in the energy spectrum at $\sim 10^7$ GeV.

We calculate the spectrum of photons resulting from electromagnetic cascades through thermal radiation, and examine the consequences of including triplet production in these cascades. We assume that the cascade is one-dimensional, and we find that this approximation is justified in the present work for thermal radiation with temperature less than $10^{-3} mc^2$. Results are obtained for both monoenergetic and power-law primary spectra, and for a variety of path lengths. We find that triplet production is particularly important in electron--photon cascades through thermal radiation when the primary energy exceeds $10^5 m^2c^4/kT$ for propagation over small path lengths. The importance of triplet production decreases as the path length increases, and it has no effect on saturated cascades.

We consider the emission of high energy to very high energy $\gamma$-rays in radio-quiet active galactic nuclei (AGN) or the central regions of radio-loud AGN. We use our results to estimate the $\gamma$-ray flux from the central regions of nearby AGN, and then to calculate the contribution to the diffuse $\gamma$-ray flux from unresolved AGN.