Articles

The Fly's Eye (1981-1993)

A "New" Experiment (1981)

After the successful test at Volcano Ranch, The University of Utah Cosmic Ray group began the construction of a full-scale observatory based on the same basic design. The new group was led by George Cassiday. Other Utah faculty members working on the new project included Haven Bergeson, Gene Loh, and Pierre Sokolsky. The new Utah experiment began observations in 1981 and was operated until 1993. A second detector site was completed in 1986.


"New" Home: Dugway Proving Ground

The new experiment was located in the West Desert of Utah, within the United States Army Dugway Proving Ground (DPG). The detectors sit atop Little Granite Mountain - also known as Five Mile Hill (FMH). Dugway is located 160 km southwest of Salt Lake City. The Proving Ground was chosen because it featured a unique combination of available infra-structure (electricity, communications) and relative isolation. The presence of the military installation also offers excellent security for personnel and equipment.

The West Desert provides an ideal location for fluorescence observations. The floor of the desert lies at an altitude of ~4,500 feet or a vertical depth of ~860 g/cm2, which is well suited for observing cosmic ray air showers in the UHE regime. The nearest population centers are more than 30 miles away, and light pollution is mostly blocked by the surrounding mountains. The average annual rainfall and snowfall are 7.5 and 15.5 inches, respectively. There is a complete absence of surface water within the boundary of DPG, and the average humidity is 48%.

For 347 days per year, the visibility at DPG is better than 10 miles. With an average wind velocity of 6 miles per hour, and its isolation from major urban centers, there is a minimum of aerosol contamination. These conditions allow the greatest number of clear, moonless nights possible for fluorescence observations, and also give the largest possible detection range for air showers to be found anywhere in the world.


The Detector: a Multi-Faceted "Eye"

Between 1981-1986, the new Utah detector array at Dugway comprised of 67 modules constructed on the same design as the three prototypes used at Volcano Ranch in 1976. The new modules are housed on corrugated steel barrels on motor-driven rotary mounts. During down-time, the barrels are rotated such that the mirrors are pointed downward to protect the module from sunlight and the elements. During run-time on clear, moonless nights they are oriented as shown in the photograph of Little Granite Mountain.

As in the Volcano Ranch prototype, each detector module contains a spherical mirror of 1.6 m diameter, and 1.5 m.

During observation, the mirrors and the PMTs are oriented to divide the sky into 880 pixels. Each pixel covers 0.0066 steradians (~5° by 5°). The arrangement of the PMT-pixels is shown in the figure below. Because of the light cone, the actual shape of the field-of-view (FOV) for each pixel is hexagonal, leaving no gaps between PMTs. As an air shower crosses the sky, its trajectory is imaged onto the PMTs firing a succession of pixels, shown by the shaded hexagons in the same figure.

The coverage of the sky with hexagonal pixels is reminiscent of the compound eye of a fly. For this reason, the Utah experiment was named "The Fly's Eye".Extending the analogy, the read-out electronic (amplifier-discriminator) boards were referred to as "ommatidial boards" (OMB), as the neural tissues behind the retina are known as ommatidial nerves. The major parameters of the Fly's Eye Experiment are listed in the table below.

The Cornell detector is triggered by requiring a coincidence between any two adjoining pixels. The signals are piped to a bank of 3" cathode ray tube displays, and recorded on 70 mm film. This detector operated for several years but was not sensitive enough to detect UHE cosmic rays reliably. In particular, the 0.1 m2 lenses are too small to collect sufficient light, and the atmosphere in Up-state New York is too contaminated with water vapor and aerosols.

Fly's Eye Detector Parameters 
Number of mirrors 67
Diameter of mirrors 1.575 m
Focal length 1.500 m
Number of PMT (and Winston cones) 880
Mirror Obscuration by PMT cluster 13%
Mirror-cone efficiency product ~0.7
PMT type EMI 9861 B
Peak PMT quantum efficiency at 360 nm 0.21
Angular Aperture per PMT 91.5 mr
Solid Angle per PMT 6.57 msr
Number of electronic channels 2640
Charge dynamic range 105 linear
Time Resolution

25 ns


Geometrical Reconstruction with a Single "Eye"

For each detected air shower, the Fly's Eye records the pattern of fired PMTs. Each hit PMT also stores the time of arrival of the light pulse from the shower, and the integrated amount of light received. From the pointing directions and the arrival times of the PMT hits, it is possible to reconstruct the trajectory of the air shower. The procedure is traditionally divided into two steps:

  1. Detector-shower plane fit
  2. Timing fit

For the "plane fit", we treat the air shower trajectory as a straight line. Strictly speaking, the shower has a lateral size of the order of ~100 m, depending on the stage of shower development. At sufficiently large distances (more than a few km), the light from the shower would appear to come from a line source. We also treat the detector itself as a single point, even though it is spread over a ~100 m wide area. Again, this size is negligible at distances of a few km or more. The point of the detector and line of the shower together defines a plane which is referred to as the "detector-shower plane", as illustrated in the figure.

The parameters of the shower-detector plane can be determined by fitting the direction cosines of the hit pixels to a plane. It is customary in this procedure to choose the detector site as the origin and fit for the normal unit vector of the shower-detector plane.

Next, the exact trajectory of the shower within the shower-detector plane can be specified by two additional parameters: (a) the impact parameter Rp, and (b) the incline angle c0, or its complementary angle y. It turns out that there exists a simple trigonometric relationship between the arrival times ti and the angles ci of the hit pixels. The angle ci for a given tube is measured from the position where the shower strikes the ground (called the "core" location) to the point of observation P. This relationship is given in the diagram below. This analytic formula involves Rp and c0 as parameters, plus an arbitrary time offset t0, and relies on the approximation that both the air shower particles and the light emitted travel at the speed of light in vacuum, c.

The parameters Rp, c0, and t0 are determined from the fit to the above formula. This step is usually referred to as the "timing fit".


A Typical Fly's Eye Air Shower Event

The first figure below shows a display of a typical air shower event recorded by Fly's Eye. The x- and y-direction cosines of the hit pixels are plotted and numbered (1-36) in time-order. This event in fact passes almost directly over the detector. The solid curve below shows the fitted shower-detector plane.

The timing-fit for the same event is shown in the plot above. The pixels are numbered exactly as before. Those pixels which were excluded from the fit are shown with a box around the number. These are typically random noise hits unassociated with the air shower, or from pixels receiving very little light and hence lacking a reliable time measurement.


The Shower Profile

Once the geometry of the event has been determined, the amount of light seen by the detector can be corrected for 1/r2 spread and attenuation losses to yield a shower size (in number of electrons, since most of the shower is dominated by the electromagnetic component.) as a function of the atmospheric depth penetrated. This function is referred to as a "shower profile".

Two important quantities can be extracted from the shower profile. First, the energy of the primary cosmic ray is proportional to the shower size. Second, the depth of the shower maximum can be read off the plot and gives a statistical measure of the composition of the primary cosmic ray. For the above event, there are ~6x108 particles at shower maximum. This translates to an energy of ~8x1017 eV for the primary particle. The shower maximum is measured at just over 800 g/cm2, which is more consistent with light composition for this energy. One should not attach too much significance to the measured value of Xmax, because the expected fluctuations are very large and conclusions should only be drawn from the average values over many events.

The actual shower profile fitting procedure is more complicated than the simplified description above would lead you to believe. The light observed by the detector includes several components in addition to the scintillation light. This is illustrated for a simulated 1018 eV event in the figure below. Here the shower development is parameterized in terms of altitude (top) and the "viewing angle" qi. The angle qi, introduced in the diagram above illustrating the timing fit, is the angle between the reverse direction of the shower and the pointing direction of the pixel in question. For example, those pixels which view the shower head-on would have qi values near zero.

Five distinct sets of curves are shown in the figure to the right. The curve marked Ne shows the shower size (relative to the scale on the right) as a function of shower development. The other four curves show the fractional contribution (left-hand side scale) to the total amount of light seen by the detector at each stage of shower development: (a) Cv: direct Cerenkov light for small viewing angles, (b) Sc: scintillation (fluorescence) light, (c) R: Cerenkov light from molecular (Rayleigh) scattering, and (d) M: Cerenkov light from particulate (Mie) scattering. All four components must be included when the observed light signal is fitted to a shower profile.


Stereoscopic Observations (1986-1993)

In 1986, a second Fly's Eye detector site was completed with 36 mirrors. The new detector was referred to as Fly's Eye 2 (FE2), and was located 3.4 km from the original site, which was henceforth denoted as FE1. The new mirrors are identical in design and construction to the original FE1 mirrors.

The addition of the FE2 detector allowed the trajectory of the air shower to be reconstructed by a process of triangulation. This "stereoscopic" reconstruction technique was much more accurate than the timing fit described previously, and hence gave a much more reliable energy measurement than was afforded by monocular measurement.

Both FE1 and FE2 sites were operated until 1993. The total monocular exposure of the original FE1 detector reached ~930 km2-sr-yrs for cosmic rays with energy 5x1019 eV (the aperture of the detector, measured in km2-sr, is energy-dependent). The final stereo exposure accumulated was ~150 km2-sr-yrs.


{multithumb popup_type=lightbox thumb_width=300}

Results From The Fly's Eye

At the time of it shutdown in 1993, the Fly's Eye had compiled the world's largest UHE cosmic ray data set. This yielded a high precision monocular energy spectrum, shown in the figure below. In this plot, the differential flux is multiplied by E3. This multiplication is customary when showing the details of UHE cosmic ray spectra.

The Fly's Eye data set was also used in an extensive search for both large- and small-scale arrival direction anisotropy. No significant excesses or deficits were observed. It turns out that the superior resolution of the stereo data set yielded some very interesting results. In particular, the stereo energy sepctrum, shown in the figure below (filled squares with error bars) revealed an apparent "dip" in the spectrum near ~3x1018 eV. This structure was not seen in the monocular spectrum with its inferior energy resolution, but have been observed by other experiments (e.g. Yakutsk and AGASA), and is usually referred to as the "ankle". It should be noted that there is siginificant disagreement between the different experiments in the actual energy of the "ankle".

In addition toe the "dip" in the energy spectrum, the stereo data also allowed a very accurate measurement of the evolution of the average Xmax with energy, as shown in the figure below (filled circles with error bars).

The slope of this semi-log plot, dXmax/dlogE, is referred to as the "elongation rate". The prediction for the absolute Xmax values is very sensitive to the hadronic models used in shower simulations. However, the elongation rates are essentially model-independent. For a single component (light element or heavy element), the value of dXmax/dlogE is predicted to be about 55 g/cm2/decade. In the figure above, the predictions for iron and protons are shown by the open circles and open aquares, respectively. The measured elongation rate above 3x1017 eV is ~70 g/cm2/decade. This slope is suggestive of a shift from a heavier to a lighter compositoin across the ankle region.

To explore this possibility further, both the energy spectrum and the Xmax data are fitted to a two-component model. One component is assumed to consist of iron, and the other protons. The iron component is assumed to have a steeper spetral index than the protons. It turns out that this model is able to simultansouly fit both sets of observations. In the energy spectrum plot above, the spectra of the two components are shown by the two dashed lines. The combined spectrum, shown by the open diamonds, reproduces the observed energy spectrum quite well. Using the same proportions of the two components extratced from this fit, one can make a prediction for the average Xmax values from this varying mix. The results are plotted in the Xmax figure above shown with open diamonds. It is clear that the predictions are in excellent agreement with the measured points.

The two component model used above is probably much too simple-minded to represent the actual situation. However, its success in describing the observations should be seen as a strong indication of a shift towards a lighter composition for cosmic rays near the "ankle" region. A possible explanation for such an observation would be the shift from a magnetically confined, heavier, intra-galatic component to a lighter, extra-galatic component.


The Highest Energy Particle Ever Recorded

In October of 1991, The FE1 detector observed an air shower with an energy of 3.2x1020 eV. This corresponds to ~50 joules or ~12 calories, or roughly the kinetic energy of a well-pitched baseball. As of the year 2012, this remains the highest energy particle ever recorded from any source. A display of the event is shown below, where the x- and z-direction cosines of the hit pixels are circled.

The 3.2x1020 eV event was not seen by the FE2 detector. The absence of a corresponding trigger in FE2 is consist with the FE1 monocular reconstruction, which places the shower outside the field of view of FE2. The plot below shows the reconstrucetd shower profile. The observed Xmax of ~800 g/cm2 is perfectly consistent with a hadron-initiated air shower of the measured 3.2x1020 eV energy.

The most important aspect of this 3.2x1020 eV event is that it is ~5 times above the theoretical Greisen-Zaptsepin-K'uzmin (GZK) cut-off at ~6x1019 eV. Super-GZK events (those above 6x1019 eV) had previously been reported by a number of ground array experiments, including Volcano Ranch, Yakutsk, and Haverah Park. The Fly's Eye event is the first such event seen with the fluorescence technique and is siginificantly higher in energy than the other candidates. This discovery was seen to give strong evidence for a non-zero flux of cosmic rays above the GZK cut-off. Since the Fly's Eye event, the AGASA group has reported a number of super-GZK events. A discussion of the GZK cut-off is given in the Physics of HiRes section of this web-site.

History of the Air Fluorescence Technique

Air Fluorescence

The term "fluorescence" refers to the process by which atoms absorb photons of one wavelength and emits photons at a longer wavelength. A common application of this effect is in nearly every household around the world: fluorescent lights. Fluorescent lights were first introduced in 1939 at the World's Fair in New York City by Interelectric corporation and offered to the public in 1940. An electric current passes through an elongated bulb, colliding with mercury atoms. The collision process excites the mercury atoms, which then emits ultra-violet (UV) light. This emission is actually referred to as "luminescence" or "scintillation". These UV photons are then absorbed by the phosphor coating of the bulbs, which re-emits in the visible. It is of course the re-emission process which is properly called "fluorescence".

The passage of charged particles in an extensive air shower through the atmosphere results in the ionization and excitation of the gas molecules (mostly nitrogen). Some of this excitation energy is emitted in the form of visible and UV radiation.

Rigorously speaking, this is a "luminescence" process analogous to the emission by mercury in a fluorescent light. Much to the horror of the optical physicists, the name "Air Fluorescence" has been adopted by the astrophysics community to describe the scintillation light from extensive air showers. This misuse of the term is in part due to the apparent similarity to the workings of a fluorescent light. On the positive side, this usage makes it easy to distinguish between a fluorescence detector from a scintillation detector (the latter is the name commonly used for desk-top particle detection devices made from inorganic salts or organic plastics).


The Origin of the Fluorescence Technique

Air fluorescence was studied in the early 1960's by the Los Alamos National Laboratory (LANL), then called the Los Alamos Scientific Laboratory. The technique was investigated as a method for detecting the yield of nuclear explosions in atmospheric tests.

Many charged particles are expelled from a nuclear explosion, and these particles will also produce scintillation light as they pass through air. The amount of light collected can then be use to estimate the total amount of energy released from the device.


The Cornell Experiment

{multithumb popup_type=slimbox thumb_width=200}

The first attempts to observe extensive air showers by the fluorescent (more correctly: luminescent) emissions were made by a group led by Kenneth Greisen at Cornell University in the middle 1960's. This group included Dr. Alan Bunner, working then as a graduate student under Greisen. Greisen himself was the first graduate student of Bruno Rossi, one of the foremost cosmic ray physicists of the 20th Century. Rossi and Greisen both worked on the Manhattan Project in Los Alamos during World War II. Greisen was in fact an eyewitness at the Trinity test and filed an official report of his observations.

The Cornell detector images the nigh-sky using 500 photo-multiplier tubes (PMT). Each PMT corresponds to a pixel covering a solid angle of 0.01 steradian (~6° by 6°). The 500 PMT's are divided into 10 modules. Each module is equipped with a 0.1 m2 Fresnel lens seen in the foreground of the photograph above (Note: the "s" in "Fresnel" is silent, and the stress is on the second syllable).

The drawing on the right, taken from Dr. Bunner's thesis, shows a cut-away section of a detector module. The Fresnel lens is shown on the left, and the PMT's are arranged at the focal surface (roughly spherical). An optical filter is placed before the lens at the entrace aperture to reduce night-sky background and eliminate contamination from filament lamps visible near the horizon.

The Cornell detector is triggered by requiring a coincidence between any two adjoining pixels. The signals are piped to a bank of 3" cathode ray tube displays, and recorded on 70 mm film. This detector operated for several years but was not sensitive enough to detect UHE cosmic rays reliably. In particular, the 0.1 m2 lenses are too small to collect sufficient light, and the atmosphere in Up-state New York is too contaminated with water vapor and aerosols.


University of Utah: 1976

In 1976, Physicists from University of Utah were the first to detect fluorescence light from cosmic ray air showers. Three prototype modules were used in a test at Volcano Ranch near Albuquerque, New Mexico.

Each prototype module contained a 1.8m diameter mirror for light collection, with 14 PMT's at the focal plane. Each PMT covers a solid angle about 0.008 steradians (~5° by 5°) in the sky. The large mirrors provided a 20-fold increase in the light collection area over that of the lenses used in the Cornell detector. The clear desert air also provided much improved visibility over the Cornell experiment.

Volcano Ranch was chosen for the test because the site also hosted a large ground array which had been operated since 1958 by a group from MIT under the leadership of John Linsley and Livio Scarsi. On Thanksgiving Night in 1976, the prototype fluorescence detectors were able to observe air showers in coincidence with the ground array.

It should be noted that the Volcano Ranch Ground Array itself was one of the great pioneering efforts in UHE cosmic ray physics. It was known as the "Desert Queen" within MIT. In 1961, the Volcano Ranch Array reported the detection of a cosmic ray air shower with a measured energy of 1020 eV. The significance of such a cosmic ray shower was not realized until after the discovery of the Cosmic Microwave Background Radiation (CMBR) in 1965 by Penzias and Wilson.

Selected References

  1. Bunner, A. N., Cosmic Ray Detection by Atmospheric Fluorescence, Ph. D. Thesis, Cornell University (1967).

Ultra High Energy Cosmic Rays

Energy Spectrum of Cosmic Rays

Cosmic rays have been observed with energies from 109 eV to over 1020 eV. Over this range, the "flux" of cosmic rays (the number of arriving particles per unit area, per solid angle, per unit time) appears to follow a single power law ~E-3. The variation of the flux with energy is referred to as the "Energy Spectrum". This spectrum is shown in the figure to the right. Remarkably, it appears to be a smooth curve over 10 decades of energy with only a few noticeable structures. The most prominent of these are small, abrupt changes in the local spectral index (the power exponent ~3) just above 1015 eV and again just above 1018 eV. Within the Cosmic Ray field, these features are usually referred to as the "knee" and the "ankle", respectively.


The Mystery of Ultra-High Energy Cosmic Rays

Cosmic Rays with energies above ~1018 eV are referred to as "Ultra-High Energy Cosmic Rays" (UHECR). These are microscopic particles with a macroscopic amount of energy-about a joule (for comparison: an electron-volt is about 1.6x10-19 joules) or more. The existence of such energetic particles remains a mystery. The three main questions concerning UHECR's are:

  • How are UHE cosmic rays accelerated to such extreme energies?
  • Where do UHE cosmic rays come from?
  • What is the composition of the UHE cosmic rays?

To study the acceleration mechanism, one must make careful measurements of the energy spectrum of UHE cosmic rays to compare to the predictions from different acceleration models. To understand where the UHE cosmic rays come from, one needs to make a careful survey of the arrival directions, and search for both small- and large-scale anisotropies in their distribution.

Composition one of the most difficult measurements because UHE cosmic rays cannot be detected directly using conventional particle detectors (see Detection of Ultra-High Energy Cosmic Ray). Consequently, the composition must be inferred from auxiliary measurements. These are discussed in the other sections.


Acceleration Mechanisms

One of the earliest theories on the acceleration of cosmic rays was proposed by Enrico Fermi in 1949 [1]. It became known as the "Second Order Fermi Mechanism". In this model, particles collide stochastically with magnetic clouds in the interstellar medium. Those particles involved in head-on collisions will gain energy (similar to a sling-shot process used to accelerate spacecrafts around planets), and those involved in tail-end collisions will lose energy. On average, however, head-on collisions are more probable. In this way, particles gain energy over many collisions.

This mechanism naturally predict a power law energy spectrum, but the power index depends on the local details of the model and would not give rise to a universal power law for cosmic rays arriving from all directions. This mechanism is also too slow and too inefficient to account for the observed UHE cosmic rays.

A more efficient version of Fermi Acceleration was proposed independently by a number of workers in the late 1970's [2-5]. In this model, particles are accelerated by a strong shock propagating through interstellar space. The following gives a schematic of the process as described in Prof. Longhair's book [6]:

Consider the case of a strong shock propagating at a supersonic, but non-relativistic speed U through a stationary interstellar gas. Figure (a) at left shows the situation in the rest frame of the gas: the density, pressure, and temperature of the gas upstream and downstream of the shock front are r2, p2, T2 and r1, p1, T1, respectively.

When viewed in the rest frame of the shock front as in figure (b) below, particles are arriving from downstream with speed v1=U and exiting upstream at speed v2. Conservation of the number of particles implies the relation: r1v1=r2v2. In the case of strong shock we expect r2/r1=(g+1)/(g-1), where g is the usual ratio of heat capacities. For a fully ionized plasma, one expects g=5/3, leading to a velocity ratio of v1/v2=4.

First order Fermi acceleration naturally predicts a power law spectrum of DNA(E)/dE ~ E-2. While the power index of 2 does not agree with the measured index of ~3, this model predicts, for the first time, a power law spectrum with a unique spectral index that is independent of the details of the local environment. The mechanism requires only the presence of strong shocks, which are quite plausibly present in the suspected sources of cosmic rays.


Possible Sources of UHE Cosmic Rays

The leading candidates for the source of UHE cosmic rays are large, energetic structures where strong shocks are expected to be found. The most well known of these are supernova remnants, which have long been suspected to generate cosmic rays. In 1995, Japan's ASCA X-ray Satellite, reported positive observation non-thermal X-ray emissions from the Supernova Remnant SN1006. The observed emission spectrum is consistent with synchrotron emission by accelerated charged particles. This report is widely seen as confirmation of supernova remnants as a known source of cosmic rays.

The observed emission from SN1006, with some fine tuning of the emission models, can explain the existence of cosmic rays up to ~1015 eV. However, it is difficult to explain the existence of cosmic rays above 1018 eV, because supernovae are simply not large enough to maintain acceleration to the UHE regime. Furthermore, no positive correlation has been observed between the arrival directions of UHE cosmic rays and supernova remnants.

There are many larger objects in the sky where strong shocks are expected. For example, strong shocks are possible around colliding galaxies such as NGC 4038/9. However, there is no evidence to indicate these objects are sources of UHE cosmic rays.

Another class of objects which are candidate sources of UHE cosmic rays are active galactic nuclei (AGN). AGN is the generic name given to a class of galaxies which are suspected to have at their center a super massive black-holes. AGN's are typically accompanied by jets which can extend 50-100 thousand light-years. Roughly one of every ten known galaxy is an AGN. It is therefore always possible to find an AGN within error of the arrival direction of a UHE cosmic ray. Even so, there is still no evidence to indicate that AGN's actually produce UHE cosmic rays.


Exotic Mechanisms

Other ideas for explaining the existence of UHE cosmic rays include:

  • Top-Down" Models:
    Decay or annihilation of some super-heavy particles or cosmological relics {e.g. topological defects, relic magnetic monopoles.)
  • Acceleration in Catastrophic events
  • Gamma ray bursts
  • New Physics?

Selected References

  1. Fermi, E. Phys. Rev., 75, 1169 (1949).
  2. Axelford, W. I., Leer, E., and Skadron, G. Proc. 15th International Cosmic Ray Conference, 11, 132 (1977).
  3. Krymsky, G. F. Dok. Acad. Nauk. USSR, 234, 1306 (1977).
  4. Bell, A. R. Mon. Not. R. Astr. Soc., 182, 147 (1978).
  5. Blandford, R.D. and Ostriker, J. P. Astrophys. J., 221, L29 (1978)
  6. Longhair, M. S., Chapter 21 of High Energy Astrophysics Vol. 2, 2nd Ed. Cambridge University Press (1994).

Detection of UHECRs

The Energy Spectrum and Detection of Cosmic Rays

The techniques by which cosmic rays in a given energy range are detected depend critically on the rate of arrival. In the energy spectrum to the left (which you already saw in the section on Ultra-High Energy Cosmic Rays), you see that in the energy range of 1011-1012 eV, the flux of cosmic rays is roughly 1 particle per square meter (one square meter is roughly 10 square feet). This rate is high enough to allow direct detection.

 The atmosphere absorbs most of the cosmic rays, as was demonstrated by Hess's original experiments. Radiation detected at ground level are actually secondary particles produced from interactions between primary cosmic rays and the air. To measure the primary cosmic rays directly, the detection equipment must be placed above the atmosphere. This is accomplished usually by carrying the instrument aboard high-altitude balloons flying at above 100,000 feet, on Earth-orbit satellites, or in the future aboard the International Space Station (ISS). An example of a detector scheduled for deployment on the ISS is the Alpha Magnetic Spectrometer (AMS), which is designed to search for nuclear antimatter in cosmic rays.

 At above 1015-1016 eV, the flux of cosmic rays drops to below one particle per square meter per year. This rate makes direct measurements impractical, as it would require flying very large detectors in order to collect sufficient number of particles. A different method is required.


<h3>Extensive Air Showers</h3>

Over the last 70 years, physicists have studied cosmic rays with energies in excess of ~1014 eV by using the Earth's atmosphere itself as part of the detection equipment .This technique takes advantage of the fact that interaction between high-energy cosmic rays and the air produces a correlated cascade of secondary particles. The process begins with the collision of the primary cosmic ray with a nucleus near the top of the atmosphere. This first collision typically produces more than 50 secondary particles, a majority of which are pi-mesons* (usually referred to as pions).

Pions come in three different flavors: positively charged, negatively charged, and neutral. All pions are unstable, but the charged pions are relatively long-lived and will most probably collide with another nucleus before decaying. The subsequent collisions are similar in nature to the primary collision. This process then leads to a cascade of particles which is referred to as a "hadronic shower".

One third of the pions produced are neutral. The neutral pions are very short-lived and will almost all decay into a pair of photons (gamma rays) before interacting with nuclei in the atmosphere. The photons interact with the nuclei in the air to produce electron-positron pairs, which in turn will produce photons via the "bremsstrahlung" process. This cascading process leads to the formation of an "electromagnetic shower". The hadronic shower itself is continuously producing neutral pions and thus initiating secondary electromagnetic showers along its path.

While high-energy cosmic rays are believed to consist mostly of charged nuclei. Gamma rays have been observed with energies as high as ~1012 eV. In the case of a gamma-ray primary particle, the particle shower produced would be purely electromagnetic. Generically, both types of cascades are called "extensive air showers" (EAS).

Extensive air showers were discovered in the 1930's by French physicist Pierre Victor Auger. In addition to his contributions to the field of cosmic rays, Pierre Auger was most well known for his discovery in the 1920's of a spontaneous process by which an atom with a vacancy in the K-shell achieves a more stable state by the emission of an electron instead of an X-ray photon, commonly known as the Auger Effect. This process forms the basis for the technique of Auger Electron Spectroscopy developed in the 1960's for characterizing surface properties of materials. Pierre Auger held the position of the Head of the Natural Science Sector at UNESCO during the years 1948-1959. Between 1962-1967, he served as Director General for the European Space Research Organization, the forerunner of the European Space Agency (ESA).

As an EAS develops into the atmosphere, more and more particles are produced. A small fraction of the kinetic energy of the primary particle is converted into mass energy. The remaining kinetic energy is then distributed over the shower. The process of multiplication continues until the average energy of the shower particles is insufficient to produce more particles in subsequent collisions. This point of the EAS development is called the "shower maximum". Beyond the maximum, the shower particles are gradually absorbed with an attenuation length of ~200 g/cm2 (rigorously this is a measure of the depth of material penetrated by the shower, which will be explained below).

 *Properties of elementary particles can also be found at the web-site of the Particle Data Group.


Properties of the Shower Maximum

Two properties of the shower maximum are important to note. First, at maximum, an EAS typically contains ~1-1.6 particles for every GeV (109 eV) of energy carried by the primary cosmic ray. Second, the average "slant depth"at which the shower maximum occurs, varies logarithmically with the energy of the primary cosmic ray.{multithumb thumb_width=250}

The "slant depth" refers to the amount of materials penetrated by the shower at a given point in its development, and is customarily denoted by the symbol "X". The value of X is calculated by integrating the density of air from the point of entry of the air shower at the top of the atmosphere, along the trajectory of the shower, to the point in question. Hence X has units given as density (g/cm3) multiplied by distance (cm). An air shower traveling along an exactly vertical, downward trajectory traverses ~1,000 g/cm2 in reaching sea-level. This value of 1,000 g/cm2 can also be interpreted as an atmospheric pressure. Obviously, an inclined shower will traverse more than 1,000 g/cm2 to reach sea-level.

Following the above convention, the depth of shower maximum is denoted "Xmax". With a value of about 500 g/cm2 at 1015 eV, the average Xmax for cosmic ray showers increases by 60-70 g/cm2 for every decade of energy.

Various hadronic shower models tend to predict significantly different absolute values for average Xmax. This makes direct comparison of measured Xmax to theoretical predictions somewhat problematic as a means of studying composition. However, nearly all the models predict (a) the same slope dXmax/dLog10E ~ 55 g/cm2 for any single element, and (b) roughly the same separation in Xmax between heavier and lighter elements. A deviation in the measured slope dXmax/dLog10E, referred to as the "Elongation Rate", from the canonical single-component value would be a clear indication of an evolution in the composition mix with energy.


Ground Arrays

For air showers with energies in excess of 1015 eV, the shower maximum penetrates to half the vertical atmospheric depth or more. There is also sufficient number of particles in the cascade such that the remnant of the shower can be detected as a correlated event by an array of individual particle detectors on the ground. The threshold (the lowest energy detectable by an instrument) of such a "Ground Array" depends on the altitude of the array. Typically it is difficult to measure cosmic rays with energies below 1014 eV with ground arrays.

The figure to the left shows the schematic of a ground array. Each station of the array samples the density of particles in its neighborhood of the shower. The footprint of air showers typically extends hundreds of meters. The particles in the air shower arrive in the form of a thin pancake traveling at essentially the speed of light. By measuring the time of arrival of the shower front at the individual stations, the direction of the primary cosmic rays can be calculated to about one degree accuracy. Conventionally, the energy of the shower is deduced from r(600), the density measured at 600 meters from the core of the shower at ground level. This density was chosen because it was the quantity which showed least amount of variations between different shower models.

An early example of a ground array is the Haverah Park array operated by University of Leeds between 1967-1987. A more recent experiment is the Akeno Giant Air Shower Array (AGASA) operated by University of Tokyo. When coupled with an underground muon array, it is also possible to measure the composition of the primary particle with a ground array. An example of such a combination used to search for very-high energy gamma rays was the CASA-MIA array.

CASA-MIA, depicted in this photograph, was operated between 1989-1997 at the U.S. Army Dugway Proving Ground. It is commonly agreed in the astrophysics community that CASA the definitive search for very-high energy gamma ray sources with energies above ~1014 eV. Their results effectively put to rest the controversy in the 1980's over reports of gamma ray detections from Cygnus X-3 and Hercules X-1.