We present and compare the predictions of various cosmic-ray Monte Carlo models for the energy (dE/deta) and particle (dN/deta) flows in p-p, p-Pb and Pb-Pb collisions at sqrt(s) = 14, 8.8, and 5.5 TeV respectively, in the range covered by forward LHC detectors like CASTOR or TOTEM (5.2<|eta|<6.6) and ZDC or LHCf (|eta|>8.1 for neutrals).

In the course of its operation, the EGRET experiment detected high-energy gamma ray sources at energies above 100 MeV over the whole sky. In this communication, we search for large-scale anisotropy patterns among the catalogued EGRET sources using an expansion in spherical harmonics, accounting for EGRET's highly non-uniform exposure. We find significant excess in the quadrupole and octopole moments. This is consistent with the hypothesis that, in addition to the galactic plane, a second mid-latitude (5^{\circ} < |b| < 30^{\circ}) population, perhaps associated with the Gould belt, contributes to the gamma ray flux above 100 MeV.

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.

Linsley used the Volcano Ranch array to collect data on the lateral distribution of showers produced by cosmic rays at energies above 10^17 eV. Very precise measurements of the steepness of the lateral distribution function were made on 366 events. The current availability of sophisticated hadronic interaction models has prompted an interpretation of the measurements. In this analysis we use the AIRES Monte Carlo code to generate showers, together with GEANT4 to simulate the detector response to ground particles. The results show that, with the assumption of a bi-modal proton and iron mix, iron is the dominant component of cosmic rays between 5x10^17 and 10^19 eV, assuming that hadronic interactions are well-described by QGSJET at this energy range.

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.

Accurate measurements of the lateral distribution of extensive air showers produced by cosmic rays of energy greater than $10^{17}$ eV were made in the 1970s by Linsley. At the time, the state of knowledge about the best hadronic interaction model to use to describe such data prevented conclusions about the cosmic ray mass from being drawn. We have used a modern model, {\sc qgsjet98}, to infer the primary mass from these data, using the very careful records left by Linsley. We find that at a median energy of $10^{18}$ eV, the data are well described by an iron-dominated composition (88 $\pm$ 6(stat) $\pm$ 18 (syst))%. We discuss the systematic errors in this estimate that arise from model uncertainties and from the range of energies used in the work of Linsley. These data are used with the permission of the late John Linsley to whom this work is dedicated.

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.

The magnetospheres of accreting neutron stars develop electrostatic gaps with huge potential drops. Protons and ions, accelerated in these gaps along the dipolar magnetic field lines to energies greater than 100 TeV, can impact onto a surrounding accretion disc. A proton-induced cascade so develops, and $\nu$-emission is produced from charged pion decays. Using GEANT4, a computer code that tracks particles produced in high energy collisions, we have calculated the resulting $\nu$-spectrum with extensive disc shower simulations. We show that the $\nu$-spectrum produced out of the proton beam is a power law. We use this result to propose accretion-powered X-ray binaries (with highly magnetized neutron stars) as a new population of point-like $\nu$-sources for km-scale detectors, such as ICECUBE. As a particular example we discuss the case of A0535+26. We show that ICECUBE should find A0535+26 to be a periodic $\nu$-source: one for which the formation and loss of its accretion disc can be fully detected. Finally, we briefly comment on the possibility that smaller telescopes, like AMANDA, could also detect A0535+26 by folding observations with the orbital period.

In the course of its operation, the Volcano Ranch array collected data on the lateral distribution of showers produced by cosmic rays at energies above $10^{17}$ {\rm eV}. From these data very precise measurements of the steepness of the lateral distribution function, characterized by the $\eta$ parameter, were made. The current availability of sophisticated hadronic interaction models has prompted a reinterpretation of the measurements. We use the interaction models {\sc qgsjet} and {\sc sibyll} in the {\sc aires} Monte Carlo code to generate showers together with {\sc geant4} to simulate the response of the detectors to ground particles. As part of an effort to estimate the primary mass composition of cosmic rays at this energy range, we present the results of our preliminary analysis of the distribution of $\eta$.

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.