Simulations of Galactic Cosmic Ray impacts onto the Herschel/PACS photoconductor arrays with Geant4 code

We present results of simulations performed with the Geant4 software code of the effects of Galactic Cosmic Ray impacts on the photoconductor arrays of the PACS instrument. This instrument is part of the ESA-Herschel payload, which will be launched in 2008 and will operate at the Lagrangian L2 point of the Sun-Earth system. Both the Satellite plus the cryostat (the shield) and the detector act as source of secondary events, affecting the detector performance. Secondary event rates originated within the detector and from the shield are of comparable intensity. The impacts deposit energy on each photoconductor pixel but do not affect the behaviour of nearby pixels. These latter are hit with a probability always lower than 7%. The energy deposited produces a spike which can be hundreds times larger than the noise. We then compare our simulations with proton irradiation tests carried out for one of the detector modules and follow the detector behaviour under ‘real’ conditions.     

On the secondary production of galactic cosmic ray electrons

The differential energy spectrum and charge ratio of primary cosmic ray electrons produced by collisions of primary cosmic ray particles with the interstellar medium is calculated by means of the two temperature statistical model of high-energy interactions. Two realistic models for the primary cosmic ray flux are considered. Contributions to the primary cosmic ray electron intensity from both the pion and kaon decay modes have been included. The distribution of matter in the galaxy and energy loss of produced secondaries and electrons are considered. The results are compared to recent experimental data.     

Effects of active solar regions on the galactic cosmic ray intensity

During the year 1969 two long-lived centres were active on the Sun at Carrington longitudes 50° < L < 100° and 250° < L < 300°. About 80% of the flares of importance 1B, produced during this period, originated in these zones.The solar modulation of galactic cosmic ray intensity during 1969 was dominated by effects resulting from the activity in the two zones. In fact all the decreases can be related to the passage at the central meridian of the active centres. Persistence of the effects connected to solar regions is found also during rotations in which they do not produce flares in front of the Earth.Seventeen among the twenty-six intensity decreases, observed during this period, can also be correlated to individual flares belonging to the region at central meridian (longitudes ± 40° with respect to the CM).The data suggest that two phenomena are operative to produce decreases of the cosmic ray flux: the passage of the interplanetary corotating stream associated with the active region near the central meridian and the blast wave produced by the flares in front of the Earth.     

Modulation of galactic cosmic ray anisotropy in heliomagnetosphere: Influence of cosmic ray scattering on sidereal daily variation

In previous papers, the present authors have shown that the galactic anisotropy is modulated due to cosmic ray orbital deflection in the heliomagnetosphere, and that the sidereal time daily variations of galactic origin can be expressed using the basic vectors, which have been obtained by calculating trajectories of cosmic rays in a model magnetosphere having Parker's Archimedian spiral structure with a flat or a wavy neutral sheet. In the present paper, the magnetic irregularities superposed on the Parker's spiral field have been taken into account, which cause the scattering of cosmic rays and disturb their orbits. We examined the fluctuations of asymptotic directions calculating their orbits by the Monte-Carlo simulation, based on the theory of the multiple scattering process. It is shown that the dispersion of the projected deviation angle is determined mainly by the scattering mean free path and by the structure of the order magnetic field, e.g. the polarity state of the heliomagnetosphere and the extent of the neutral sheet. We investigated also the influence of the fluctuations of asymptotic directions on the sidereal daily variation. It is found that, under some conditions, the scattering causes only the attenuation of the amplitude of the basic vector, and does not change its phase. The attenuation is negligibly small at high rigidities larger than ~ 1000 GV, but becomes more serious with decreasing rigidity. The rigidity dependence curve of the attenuation rate was calculated for various cases. A simple and approximate method is also presented for the derivation of those curves for any value of the magnitude of the mean free path and for various model magnetospheres. It is noted, however, that the lower limiting rigidity below which the present method is not applicable is relatively high in the Positive polarity state.     

The insensitivity of the cosmic ray galactic anisotropy to heliomagnetic polarity reversals

Using the cosmic ray sidereal and anti-sidereal diurnal variations observed underground in London and Hobart during the period 1958–1983, it is demonstrated that: (1) the phase changes of the apparent sidereal diurnal variation observed only in the Northern Hemisphere cannot be attributed to the change of the heliomagnetospheric modulation of galactic cosmic ray anisotropy caused by the polarity reversal of the solar magnetic field, but that they are due to the fluctuation of the spurious sidereal variation produced from the anisotropy responsible for the solar semi-diurnal variation; (2) the spurious sidereal variation can be eliminated from the apparent variation by using the observed anti-sidereal diurnal variation; and (3) after the elimination, the sidereal diurnal variations in the Northern and Southern Hemispheres almost coincide with each other and are stationary throughout the period, regardless of the polarity reversal of the heliomagnetosphere. The origin of the corrected sidereal variation is discussed.     

Spatial distribution of high energy cosmic ray electrons perpendicular to the galactic plane

With the help of empirical data concerning the latitudinal distribution of galactic gamma rays the contribution of inverse Compton scattered gamma rays is calculated using various models concerning the distribution of high energy cosmic ray electrons perpendicular to the galactic plane. It is shown that gamma ray astronomy from regions with vanishing stellar and interstellar matter densities at energies greater than 100 MeV provides instructive information on the cosmic ray electron density. We find evidence for the existence of a broad galactic electron disk with a total thickness of at least 6.4 kpc. The uncertainties of the cosmic ray electron spectrum measurements above 100 GeV imply an additional uncertainty in the inverse Compton source function of at least a factor 6.     

A model of possible variations of the galactic cosmic ray intensity over the recent billion years

Based on the analysis of published data on exposure ages of iron meteorites determined with the ⁴⁰K/K method (T _K) and ages calculated using short-lived cosmogenic radionuclides (with the half-life T _1/2 < 1 Myr) in combination with stable cosmogenic isotopes of noble gases (T_RS), the following results have been obtained. (1) The distribution of T _RS ages (106 values) has an exponential shape, similar to that for ordinary chondrites, but different from the distribution of T _K ages (80 values). The difference is most likely due to small amounts of data for meteorites with low T _K ages (less than ~200–300 Myr). The latter can be ascribed to the difficulty of measurement of small concentrations of cosmogenic potassium isotopes. This circumstance makes the selection of meteorites with ⁴⁰K/K ages nonrepresentative and casts doubt on the correctness of conclusions about the variations of the intensity of galactic cosmic rays (GCR) based on the analysis of distribution of these ages. (2) The magnitude of the known effect (systematic overestimation of T _K ages in comparison with T _RS ages) has been refined. The value k = T _K/T _RS = 1.51 ± 0.03 is acquired for the whole population of data. We have shown the inefficiency of the explanation of this effect on account of an exponential change in the GCR intensity (I _T ) with time (T) according to the relation I _T = I ₀exp(–γT) over the whole range of ages of iron meteorites. (3) In order to explain the overestimation of T _K ages in comparison with T _RS ages, a model has been proposed, according to which the GCR intensity has exponentially increased in the interval of 0–1500 Myr governed by the relation: I _T = I _{T = 1500} (1 + αexp(–βT)). For one of the variants of this model, the GCR intensity has exponentially increased by a factor of two only over the recent ~300 Myr, remaining approximately constant for the rest of the time. The data acquired with the use of this model indicate that the measured T _K ages are close to the actual time that the meteorites existed in space; the data are in agreement with the observed exponential distribution of T _RS ages.     

Effect of the heliosheath and standing termination shock on galactic cosmic ray propagation in a stationary heliosphere model

We consider a stationary model of the propagation of galactic cosmic rays (GCR) in the heliosphere and adjacent interstellar space. The heliosphere is assumed to be a two-layer medium consisting of two adjacent regions that are spherically symmetric relative to the sun. The solar wind velocity is supersonic in the inner heliosphere bounded by the standing termination shock, and this velocity is subsonic in the outer heliosphere bounded by the heliosheath. The GCR scattering in these regions is due to different factors characterized by relevant diffusion coefficients. The solar wind velocity is assumed to be zero in the interstellar medium, where the scattering becomes weaker. No particle sources are presumed to exist at the boundaries between the layers. An exact analytical solution of the corresponding mathematical problem can be obtained without essential difficulties, although it is extremely cumbersome. Analytical expressions for the GCR spectra of particles with very high energies (>2500 MeV) and very low energies (<1400 MeV) are obtained for each region of particle propagation. The low-energy particle distribution corresponds to the data obtained by the Voyager spacecraft. It is shown that the low-energy particle density continuously increases from the sun toward the heliospheric boundary, regardless of the scattering mode in the inner and outer parts of the heliosphere.     

First results with the boloSource() algorithm: photometry of faint standard stars observed by Herschel/PACS

The boloSource() algorithm is a tool to separate the signal of compact sources from that of the diffuse background in the timeline of far-infrared measurements performed by the PACS camera of the Herschel Space Observatory. An important characteristic and quality indicator of this method is that how well it can reproduce the flux of faint standard stars which have reliable flux estimates. For this propose we selected a few calibrator targets and constructed light curves by extracting point source flux for each repetition of the measurements independently using standard aperture photometry methods. These were compared with the light curves obtained using the boloSource() method on the same dataset. The results indicate that boloSource() provides a similar level of photometric accuracy and reproducibility as the usual flux extraction and photometry methods. This new technique will be developed further and also tested against other methods in more complex fields with the goal to make it usable for large-scale studies in the future.     

Preparation for the solar system observations with Herschel: Simulation of Jupiter observations with PACS

Observations of the water inventory as well as other chemically important species on Jupiter will be performed in the frame of the guaranteed time key project of the Herschel Space Observatory entitled “Water and related chemistry in the Solar system”. Among other onboard instruments, PACS (Photodetector Array Camera and Spectrometer) will provide new data of the spectral atlas in a wide region covering the far-infrared and submillimetre domains, with an improved spectral resolution and a higher sensitivity compared to previous observations carried out by Cassini/CIRS (Composite InfraRed Spectrometer) and by ISO (Infrared Space Observatory).In order to optimise the observational plan and to prepare for the data analysis, we have simulated the expected spectra of PACS Jupiter observations. Our simulation shows that PACS will promisingly detect several H₂O emission lines. As PACS is capable of spatially resolving the Jovian disk, we will be able to discern the external oxygen sources in the giant planets by exploring the horizontal distribution of water. In addition to H₂O lines, some absorption lines due to tropospheric CH₄, HD, PH₃ and NH₃ lines will be observed with PACS. Furthermore, owing to the high sensitivity of the instrument, the current upper limit on the abundance of hydrogen halides such as HCl will be also improved.     

From cosmic ray physics to cosmic ray astronomy: Bruno Rossi and the opening of new windows on the universe

Bruno Rossi is considered one of the fathers of modern physics, being also a pioneer in virtually every aspect of what is today called high-energy astrophysics. At the beginning of 1930s he was the pioneer of cosmic ray research in Italy, and, as one of the leading actors in the study of the nature and behavior of the cosmic radiation, he witnessed the birth of particle physics and was one of the main investigators in this fields for many years. While cosmic ray physics moved more and more towards astrophysics, Rossi continued to be one of the inspirers of this line of research. When outer space became a reality, he did not hesitate to leap into this new scientific dimension. Rossi’s intuition on the importance of exploiting new technological windows to look at the universe with new eyes, is a fundamental key to understand the profound unity which guided his scientific research path up to its culminating moments at the beginning of 1960s, when his group at MIT performed the first in situ measurements of the density, speed and direction of the solar wind at the boundary of Earth’s magnetosphere, and when he promoted the search for extra-solar sources of X rays. A visionary idea which eventually led to the breakthrough experiment which discovered Scorpius X-1 in 1962, and inaugurated X-ray astronomy.     

Measurement of the cosmic ray composition at the knee with the SPASE-2/AMANDA-B10 detectors

The mass composition of high-energy cosmic rays at energies above 10¹⁵ eV can provide crucial information for the understanding of their origin. Air showers were measured simultaneously with the SPASE-2 air shower array and the AMANDA-B10 Cherenkov telescope at the South Pole. This combination has the advantage to sample almost all high-energy shower muons and is thus a new approach to the determination of the cosmic ray composition. The change in the cosmic ray mass composition was measured versus existing data from direct measurements at low energies. Our data show an increase of the mean log atomic mass lnA by about 0.8 between 500 TeV and 5 PeV. This trend of an increasing mass through the “knee” region is robust against a variety of systematic effects.     

Operations and performance of the PACS instrument 3He sorption cooler on board of the Herschel space observatory

A ³He sorption cooler produced the operational temperature of 285 mK for the bolometer arrays of the Photodetector Array Camera and Spectrometer (PACS) instrument of the Herschel Space Observatory. This cooler provided a stable hold time between 60 and 73 h, depending on the operational conditions of the instrument. The respective hold time could be determined by a simple functional relation established early on in the mission and reliably applied by the scientific mission planning for the entire mission. After exhaustion of the liquid ³He due to the heat input by the detector arrays, the cooler was recycled for the next operational period following a well established automatic procedure. We give an overview of the cooler operations and performance over the entire mission and distinguishing in-betweenthe start conditions for the cooler recycling and the two main modes of PACS photometer operations. As a spin-off, the cooler recycling temperature effects on the Herschel cryostat ⁴He bath were utilized as an alternative method to dedicated Direct Liquid Helium Content Measurements in determining the lifetime of the liquid Helium coolant.     

The effect of solar wind high-speed streams on the galactic cosmic rays intensity

The effect of high-speed recurrent solar wind streams from coronal holes on the galactic cosmic rays intensity is investigated. The distribution of galactic cosmic rays for different solar cycles is considered based on the data of the world network of neutron monitors. Within the inhomogeneous model, which includes a homogeneous background and regions of high-speed streams (HSS’s), the transport equation has been solved and the effect of HSS’s on the spatial distribution of galactic cosmic rays is estimated. It is shown that theoretical calculations are agreed with the experimental results obtained for 2000–2014 under different assumptions about the mean free path of cosmic rays in the corresponding period of HSS’s.     

The synergy of the cosmic ray and high energy atmospheric physics: Particle bursts observed by arrays of particle detectors

Particle bursts detected on the earth's surface during thunderstorms by various particle detectors originated from the relativistic runaway electron avalanches (RREAs) initiated by free electrons accelerated in the strong atmospheric electric fields. Two oppositely directed dipoles in the thundercloud accelerate electrons in the direction of the earth's surface, and to the open space. The particle bursts observed by orbiting gamma ray observatories are called terrestrial gamma ray flashes (TGFs, with energies of several MeV, only sometimes reaching tens of MeV); ones registered by particle detectors located on the ground – are called thunderstorm ground enhancements (TGEs, with energies, usually reaching 40-50 MeV). Balloons and aircraft in the troposphere register gamma ray glows (with energies of several MeV). Recently, high-energy atmospheric physics includes also, so-called, downward TGFs (DTGFs), intense particle bursts with a duration of a few milliseconds.Well-known extensive air showers (EASs) originate from the interactions of galactic protons and fully-stripped nuclei with the atmosphere atoms. EAS particles have very dense cores around the shower axes. However, high-energy particles in the EAS cores comprise a very thin disc of (a few tens of ns), and a particle detector traversed by an EAS core will not register a particle burst, but only one very large pulse. Only neutron monitor, by collecting delayed thermal neutrons from EAS core particle interactions with soil, can register particle bursts. We discuss the relation between short particle bursts available from the largest particle arrays with EAS phenomena. We demonstrate that the neutron monitors can extend the EAS “lifetime” up to a few milliseconds, a time comparable with DTGFs duration. The possibility to use the network of neutron monitors for high-energy cosmic ray research is also deliberated.Plain Language Summary: Short and extended particle bursts are registered in space, the troposphere, and the earth's surface. Coordinated monitoring of the particle fluxes, near-surface electric fields, and lightning flashes makes it possible to formulate a hypothesis on the origin of intense bursts and their relation to extensive air showers and atmospheric discharges. Analysis of the observational data and possible origination scenarios of particle bursts allows us to conclude that the bursts can be explained by the electron acceleration in the thunderous atmosphere and by gigantic showers developed in the terrestrial atmosphere by high-energy protons and fully-stripped nuclei accelerated in Galaxy.     

Solar cosmic ray composition above 10 MeV/nucleon and its energy dependence in the 4 August 1972 event

Observations of the proton, helium, (C, N, O) and Fe-group nuclei fluxes made during the large 4 August 1972 solar particle event are presented. The results show a small, but significant variation of the composition of multiply-charged nuclei as a function of energy in the energy region above 10 MeV nucleon^–1. In particular, the He/(C, N, O) abundance ratio varies by a factor 2 between 10 and 50 MeV nucleon^–1 and the Fe-group/(C, N, O) ratio suggests a similar variation. Abundance ratios from the 4 August 1972 event are compared as a function of energy with ratios measured in other solar events to show that several of the earlier results are consistent with an energy variation like that observed in August 1972, while certain other events must have had a substantially different dependence of composition on energy. At energies 50 MeV nucleon^–1, the He/(C, N, O) abundance ratio for August 1972 is consistent with all earlier measurements made above that energy which suggests that variations may vanish at high energies.NASA/NAS Senior Resident Research Associate, on leave from TATA Institute of Fundamental Research, Bombay.     

On the protection of extrasolar Earth-like planets around K/M stars against galactic cosmic rays

Previous studies have shown that extrasolar Earth-like planets in close-in habitable zones around M-stars are weakly protected against galactic cosmic rays (GCRs), leading to a strongly increased particle flux to the top of the planetary atmosphere. Two main effects were held responsible for the weak shielding of such an exoplanet: (a) For a close-in planet, the planetary magnetic moment is strongly reduced by tidal locking. Therefore, such a close-in extrasolar planet is not protected by an extended magnetosphere. (b) The small orbital distance of the planet exposes it to a much denser stellar wind than that prevailing at larger orbital distances. This dense stellar wind leads to additional compression of the magnetosphere, which can further reduce the shielding efficiency against GCRs. In this work, we analyse and compare the effect of (a) and (b), showing that the stellar wind variation with orbital distance has little influence on the cosmic ray shielding. Instead, the weak shielding of M star planets can be attributed to their small magnetic moment. We further analyse how the planetary mass and composition influence the planetary magnetic moment, and thus modify the cosmic ray shielding efficiency. We show that more massive planets are not necessarily better protected against galactic cosmic rays, but that the planetary bulk composition can play an important role.     

The intensity variations of solar and galactic cosmic rays with azimuthal angle in the polar region

A balloon-born multidirectional detector is used to measure the intensity variation of galactic and solar cosmic rays with the azimuthal angle, the zenith angle being maintained at 60°. In polar regions, the intensity towards the north is found to be 20% larger than that towards the south. It is shown that this anisotropy does not originate in interplanetary space and is not produced by a magnetospheric source. It is suggested that the effect is due to propagation effects within the magnetosphere.     

Simulation of the interaction of galactic cosmic ray protons with meteoroids: On the production of 3H and light noble gas isotopes in isotropically irradiated thick gabbro and iron targets

Abstract— Thick spherical targets, one made of gabbro (R = 25 cm) and one made of iron (R = 10 cm), were irradiated isotropically with 1.6 GeV protons at Laboratoire National Saturne (LNS)/Saclay to simulate the interactions of galactic cosmic ray protons with meteoroids in space. At various depths, both artificial meteoroids contained a large number of high‐purity, single‐element target foils and chemical compounds of up to 28 target elements. In these individual target foils, the elemental production rates of radionuclides and noble gas isotopes were measured. Here, we report the results for the light noble gas isotopes ^3,4He, ^20,21,22Ne, and ^36,38,39Ar for the most cosmochemically relevant target elements as well as for some meteoritic material from Jilin, Farmington, and Cape York. From ³He analyses done several years apart, ³H diffusive losses during sample storage have been obtained, and direct as well as cumulative ³He production rates for O, Mg, Al, Si, Fe, Ni, and the meteoritic material are given. Losses by diffusion of tritium from metallic Mg and Fe are found to occur on time scales of months, while metallic Al, Si, and stone meteorites are much more retentive. The production rate ratios P(³H)/P(³He)_d obtained in the simulation experiments are 0.73, 1.28, and 1.16 for O, Al, and Si, respectively. These rates are based on our best knowledge about the ³H and ³He production rates and should, therefore, replace data published earlier (Leya et al. 2000a). The earlier calculations for ⁴He, ^20,21,22Ne, and ^36,38,39Ar remain valid. The new modeled correlation ³He_cum/²¹Ne versus ²²Ne/²¹Ne for chondrites exposed to cosmic rays with an energy spectrum characterized by a modulation parameter of φ = 650 MeV is in fair agreement with the empirical relationship (“Berne plot”). However, for small meteorites and little shielding in larger ones, there are systematic differences that most likely are due to an underestimation of the spallogenic ²²Ne/²¹Ne ratio by ?2%.