The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

The properties of cosmic rays with energies above 10**6 GeV have to be deduced from the spacetime structure and particle content of the air showers which they initiate. In this review we summarize the phenomenology of these giant air showers. We describe the hadronic interaction models used to extrapolate results from collider data to ultra high energies, and discuss the prospects for insights into forward physics at the LHC. We also describe the main electromagnetic processes that govern the longitudinal shower evolution, as well as the lateral spread of particles. Armed with these two principal shower ingredients and motivation from the underlying physics, we provide an overview of some of the different methods proposed to distinguish primary species. The properties of neutrino interactions and the potential of forthcoming experiments to isolate deeply penetrating showers from baryonic cascades are also discussed. We finally venture into a terra incognita endowed with TeV-scale gravity and explore anomalous neutrino-induced showers.

Using data from the SUGAR and the AGASA experiments taken during a 10 yr period with nearly uniform exposure to the entire sky, we search for large-scale anisotropy patterns in the arrival directions of cosmic rays with energies > 10^{19.6} eV. We determine the angular power spectrum from an expansion in spherical harmonics for modes out to \ell=5. Based on available statistics, we find no significant deviation from isotropy. We compare the rather modest results which can be extracted from existing data samples with the results that should be forthcoming as new full-sky observatories begin operation.

With the aim of testing recent claims for a particularly strong correlation between ultra-high energy cosmic rays (UHECRs), observed with the AGASA and the Yakutsk experiments, and a sample of BL Lacertae (BL Lacs), we here conduct a blind statistical assessment. We search for associations between the same set of BL Lac objects and the arrival directions of 33 relevant UHECRs observed with the Haverah Park and the Volcano Ranch experiments. Within the accuracy of angle determination, there are no positional coincidences. The probability that this null result arises as a statistical fluctuation from the strongly correlated case is less than 5%. This implies that the possible correlation between the arrival directions of UHECRs and BL Lacs is not statistically sustained. We discuss the impact of our findings on the propose additional connection among UHECRs, BL Lacs, and EGRET gamma-ray blazars. Recently, such an association was used as classification technique for EGRET sources. Here we show that its main underlying hypothesis, i.e., the EGRET angular uncertainty is twice that quoted in the Third EGRET Catalog, underestimates the goodness of existing gamma-ray data.

A numerical likelihood approach to the determination of cosmic ray anisotropies is presented which offers many advantages over other approaches. It allows a wide range of statistically meaningful hypotheses to be compared even when full sky coverage is unavailable, can be readily extended in order to include measurement errors, and makes maximum unbiased use of all available information.

We report progress on applying technologies developed for LHC-era experiments to cosmic ray detection, using scintillating tiles with embedded wavelength-shifting fibers and avalanche photodiode readouts as parts of a robust, inexpensive cosmic air shower detector. We are planning to deploy such detectors in high schools as part of an outreach effort able to search for long-distance correlations between airshowers.

In this review we discuss the important progress made in recent years towards understanding the experimental data on cosmic rays with energies $\agt 10^{19}$ eV. We begin with a brief survey of the available data, including a description of the energy spectrum, mass composition, and arrival directions. At this point we also give a short overview of experimental techniques. After that, we introduce the fundamentals of acceleration and propagation in order to discuss the conjectured nearby cosmic ray sources. We then turn to theoretical notions of physics beyond the Standard Model where we consider both exotic primaries and exotic physical laws. Particular attention is given to the role that TeV-scale gravity could play in addressing the origin of the highest energy cosmic rays. In the final part of the review we discuss the potential of future cosmic ray experiments for the discovery of tiny black holes that should be produced in the Earth's atmosphere if TeV-scale gravity is realized in Nature.

We examine the correlation between compact radio quasars (redshifts in the range $z = 0.3 - 2.2$) and the arrival direction of ultrahigh energy cosmic rays forming clusters. Our Monte Carlo simulation reveals a statistically significant correlation on the AGASA sample: the chance probability of this effect being less than 1%. The implications of this result on the origin of ultrahigh energy cosmic rays are discussed.

In this article we examine the hypothesis that the highest energy cosmic rays are complex nuclei from extragalactic sources. Under reasonable physical assumptions, we show that the nearby metally rich starburst galaxies (M82 and NGC 253) can produce all the events observed above the ankle. This requires diffusion of particles below $10^{20}$ eV in extragalactic magnetic fields $B \approx 15$ nG. Above $10^{19}$ eV, the model predicts the presence of significant fluxes of medium mass and heavy nuclei with small rate of change of composition. Notwithstanding, the most salient feature of the starburst-hypothesis is a slight anisotropy induced by iron debris just before the spectrum-cutoff.

At Auger energies only model predictions enable us to extract primary cosmic ray features. The simulation of the shower evolution depends sensitively on the first few interactions, necessarily related to the quality of our understanding of high energy hadronic collisions. Distortions of the standard ``soft semi-hard'' scenario include novel large compact dimensions and a string or quantum gravity scale not far above the electroweak scale. Na\"{\i}vely, the additional degrees of freedom yield unification of all forces in the TeV range. In this article we study the influence of such precocious unification during atmospheric cascade developments by analyzing the most relevant observables in proton induced showers.

One of the possible consequences of the existence of extra degrees of freedom beyond the electroweak scale is the increase of neutrino-nucleon cross sections ($\sigma_{\nu N}$) beyond Standard Model predictions. At ultra-high energies this may allow the existence of neutrino-initiated extensive air showers. In this paper, we examine the most relevant observables of such showers. Our analysis indicates that the future Pierre Auger Observatory could be potentially powerful in probing models with large compact dimensions.

In this article we examine the question of whether the highest energy cosmic ray primaries could be ultra relativistic magnetic monopoles. The analysis is performed within the framework of large compact dimensions and TeV scale quantum gravity. Our study indicates that while this hypothesis must be regarded as highly speculative it cannot be ruled out with present data.

A re-examination of the atmospheric cascade profile of the highest energy cosmic ray is presented. The study includes air-shower simulations considering different cross sections, particle multiplicity and variation of the hadronic-event-generator to model interactions above 200 GeV. The analysis provides evidence that a medium mass nucleus primary reproduces the shower profile quite well. This result does not support the idea, increasingly popular at present, that the highest energy particles are protons, derived from the decay of supermassive relic particles. On the other hand, we show that debris of relativistic super-heavy nuclei, which can survive a 100 Mpc journey through the primeval radiation are likely to generate such a kind of cascade.

We present an evaluation of a simulated cosmic ray shower, based on {\sc geant4} and {\sc top-c}, which tracks all the particles in the shower. {\sc top-c} (Task Oriented Parallel C) provides a framework for parallel algorithm development which makes tractable the problem of following each particle. This method is compared with a simulation program which employs the Hillas thinning algorithm.

We study the self-similar structure of electromagnetic showers and introduce the notion of the fractal dimension of a shower. Studies underway of showers in various materials and at various energies are presented, and the range over which the fractal scaling behaviour is observed is discussed. Applications to fast shower simulations and identification, particularly in the context of extensive air showers, are also discussed.

{\sc top-c} (Task Oriented Parallel C) is a freely available package for parallel computing. It is designed to be easy to learn and to have good tolerance for the high latencies that are common in commodity networks of computers. It has been successfully used in a wide range of examples, providing linear speedup with the number of computers. A brief overview of {\sc top-c} is provided, along with recent experience with cosmic ray physics simulations.

We present a detailed simulation of the performance of water Cerenkov detectors suitable for use in the Pierre Auger Observatory. Using {\sc geant4}, a flexible object-oriented simulation program, including all known physics processes, has been developed. The program also allows interactive visualization, and can easily be modified for any experimental setup.

We critically review the common belief that ultrahigh energy cosmic rays are protons or atomic nuclei with masses not exceeding that of iron. We find that heavier nuclei are indeed possible, and discuss possible sources and acceleration mechanisms for such primaries. We also show detailed simulations of extensive air showers produced by ``superheavy'' nuclei, and discuss prospects for their detection in future experiments.

A re-examination of the energy cosmic ray spectrum above $10^{20}$ eV is presented. The overall data-base provides evidence, albeit still statistically limited, that non-nucleon primaries could be present at the end of the spectrum. In particular, the possible appearance of superheavy nuclei (seldom discussed in the literature) is analysed in detail.

We describe a scalable parallelization of Geant4 using commodity hardware in a collaborative effort between the College of Computer Science and the Department of Physics at Northeastern University. The system consists of a Beowulf cluster of 32 Pentium II processors with 128 MBytes of memory each, connected via ATM and fast Ethernet. The bulk of the parallelization is done using TOP-C (Task Oriented Parallel C), software widely used in the computational algebra community. TOP-C provides a flexible and powerful framework for parallel algorithm development, is easy to learn, and is available at no cost. Its task oriented nature allows one to parallelize legacy code while hiding the details of interprocess communications. Applications include fast interactive simulation of computationally intensive processes such as electromagnetic showers. General results motivate wider applications of TOP-C to other simulation problems as well as to pattern recognition in high energy physics.