Prospects for observing dark-matter remnants with GLAST

The satellite-based experiment, GLAST (Gamma-ray Large Area Space Telescope), is under construction and is planned to measure the cosmic γ-ray flux in the energy range 20 MeV to >300 GeV, with supporting measurements for γ-ray bursts from 10 keV to 25 MeV. With its launch in 2007, GLAST will open a new and important window on a wide variety of high-energy phenomena, including exotic relics from the Big Bang. Among these may be the decay/annihilation products of the hypothesized super symmetric image of the known particles. Single-photon energy thresholds for channels leading to such final states have been excluded in a model-dependent manner by accelerator searches to energies greater than 50 GeV. The ability of GLAST to set limits on this important component of cosmological evolution is presented along with an update on the present status of this mission.     

Relationship of solar flare accelerated particles to solar energetic particles (SEPs) observed in the interplanetary medium

Observations of hard X-ray (HXR)/γ-ray continuum and γ-ray lines produced by energetic electrons and ions, respectively, colliding with the solar atmosphere, have shown that large solar flares can accelerate ions up to many GeV and electrons up to hundreds of MeV. Solar energetic particles (SEPs) are observed by spacecraft near 1 AU and by ground-based instrumentation to extend up to similar energies as in large SEP events, but it appears that a different acceleration process, one associated with fast coronal mass ejections is responsible. Much weaker SEP events are observed that are generally rich in electrons, ³He, and heavy elements. The energetic particles in these events appear to be similar to those accelerated in flares. The Ramaty high energy solar spectroscopic imager (RHESSI) mission provides high-resolution spectroscopy and imaging of flare HXRs and γ-rays. Such observations can provide information on the location, energy spectra, and composition of the flare accelerated energetic particles at the Sun. Here, preliminary comparisons of the RHESSI observations with observations of both energetic electron and ion near 1 AU are reviewed, and the implications for the particle acceleration and escape processes are discussed.     

Gamma-ray observations of explosive nucleosynthesis products

In this review I discuss the various γ-ray emission lines that can be expected and, in some cases have been observed, from radioactive explosive nucleosynthesis products. The most important γ-ray lines result from the decay chains of ⁵⁶Ni, ⁵⁷Ni, and ⁴⁴Ti. ⁵⁶Ni is the prime explosive nucleosynthesis product of Type Ia supernovae, and its decay determines to a large extent the Type Ia light curves. ⁵⁶Ni is also a product of core-collapse supernovae, and in fact, γ-ray line emission from its daughter product, ⁵⁶Co, has been detected from SN1987A by several instruments. The early occurrence of this emission was surprising and indicates that some fraction of ⁵⁶Ni, which is synthesized in the innermost supernova layers, must have mixed with the outermost supernova ejecta.Special attention is given to the γ-ray line emission of the decay chain of ⁴⁴Ti (⁴⁴Ti → ⁴⁴Sc → ⁴⁴Ca), which is accompanied by line emission at 68, 78, and 1157 keV. As the decay time of ⁴⁴Ti is ∼86 yr, one expects this line emission from young supernova remnants. Although the ⁴⁴Ti yield (typically 10⁻⁵–10⁻⁴M_⊙) is not very high, its production is very sensitive to the energetics and asymmetries of the supernova explosion, and to the mass cut, which defines the mass of the stellar remnant. This makes ⁴⁴Ti an ideal tool to study the inner layers of the supernova explosion. This is of particular interest in light of observational evidence for asymmetric supernova explosions.The γ-ray line emission from ⁴⁴Ti has so far only been detected from the supernova remnant Cas A. I discuss these detections, which were made by COMPTEL (the 1157 keV line) and BeppoSAX (the 68 and 78 keV lines), which, combined, give a flux of (2.6 ± 0.4 ± 0.5) × 10⁻⁵ ph cm⁻² s⁻¹ per line, suggesting a ⁴⁴Ti yield of (1.5 ± 1.0) × 10⁻⁴M_⊙. Moreover, I present some preliminary results of Cas A observations by INTEGRAL, which so far has yielded a 3σ detection of the 68 keV line with the ISGRI instrument with a flux that is consistent with the BeppoSAX detections. Future observations by INTEGRAL-ISGRI will be able to constrain the continuum flux above 90 keV, as the uncertainty about the continuum shape, is the main source of systematic error for the 68 and 78 keV line flux measurements. Moreover, with the INTEGRAL-SPI instrument it will be possible to measure or constrain the line broadening of the 1157 keV line. A preliminary analysis of the available data indicates that narrow line emission (i.e., Δv < 1000 km s⁻¹) can be almost excluded at the 2σ level, for an assumed line flux of 1.9 × 10⁻⁵ ph cm⁻² s⁻¹.     

A study on emission geometry and the separation of peaks in the γ-ray pulse profiles of pulsar

In the inner annular gap (IAG) model the γ-ray radiation sources are suggested to be located close to the null charge surface (NCS). A method to explore two geometric parameters (κ, λ) of the IAG model is given in this paper. Pulsar population statistics are proposed to test the radiation locations of γ-ray pulsars. Within the IAG model, predictions for the GLAST observation of γ-ray pulsars are also made.     

What gamma-ray deexcitation lines reveal about solar-flares

Gamma-ray emission from solar flares reveals information about the nature of the accelerated particles and about the physical conditions of the medium through which the accelerated particles are transported. In this paper, we present the gamma-ray line-production and loop transport models used in our calculations of high-energy emission. We discuss the calculated interaction time history, the depth distribution, the interacting-particle angular distribution, and fluence ratios of the narrow gamma-ray lines. We show the relationship between the γ-ray observables and the parameters of the transport and line-production models. For illustration, we use calculations of 4.44 MeV ¹²C nuclear deexcitation line-production. Applications of the calculations to flare observations by both SMM and RHESSI are also presented.     

The neutral gas in the environs of the Geminga gamma-ray pulsar

We present a high-resolution (24″) study of the HI interstellar gas distribution around the radio-quiet neutron star Geminga. Based on very large array and MPIfR Effelsberg telescope data, we analyzed a 40′ × 40′ field around Geminga. These observations have revealed the presence of a neutral gas shell, 0.4 pc in radius, with an associated HI mass of 0.8M_⊙, which surrounds Geminga at a radial velocity compatible with the kinematical distance of the neutron star. In addition, morphological agreement is observed between the internal face of the HI shell and the brightest structure of Geminga’s tail observed in X-rays. We explore the possibility that this morphological agreement is the result of a physical association.     

Highly significant detection of solar neutrons on 2005 September 7

We have successfully detected solar neutrons at ground level in association with the X17.0 solar flare that occurred on 2005 September 7. Observations were made with the solar neutron telescopes and neutron monitors located in Bolivia and Mexico. In this flare, large fluxes of hard X-rays and γ-rays were observed by the GEOTAIL and the INTEGRAL satellites. The INTEGRAL observations include the 4.4 MeV line γ-rays of ¹²C. The data suggest that solar neutrons were produced at the same time as these hard electromagnetic radiations. We have however found an apparent discrepancy between the observed and the expected time profiles. This fact suggests a possible extended neutron emission.     

Analysis of electromagnetic and electrostatic effects of particle impacts on spacecraft

Particle impacts on spacecraft can cause considerable damage, even leading to complete failure. A theory for the resulting vehicle potential changes and the electromagnetic radiation from impact-induced plasma has been published by Close et al. (2010). Here we compare this theory to impacts registered by the Radio and Plasma Wave Science instrumentation on the Cassini spacecraft. We study both low-velocity (16 km/s) large particles (2.6 μm radius) detected in Saturn’s rings and high-velocity (450 km/s) small particles (1 nm radius) in the solar wind. The agreement with the theory is quite good. We also apply these results to earth orbit and conclude that both Electrostatic Discharge and Electromagnetic Pulse radiation are likely and could lead to spacecraft failure.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

On the 27-day variations of the galactic cosmic ray anisotropy and intensity for different periods of solar magnetic cycle

The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy calculated based on the neutron monitors experimental data is larger in the qA > 0 period than in the qA < 0 period of solar magnetic cycle. The amplitudes of the 27-day variation of the galactic cosmic rays anisotropy do not depend on the tilt angles of the heliospheric neutral sheet for different the qA > 0 and the qA < 0 periods of solar magnetic cycle. A good correlation has been revealed between the changes of the amplitudes of the 27-day variations of the galactic cosmic ray anisotropy and intensity versus the qA > 0 and the qA < 0 periods of solar magnetic cycle.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

Detection of meteors and sub-relativistic dust grains by the fluorescence detectors of ultra high energy cosmic rays

Fluorescence detectors of ultra high energy cosmic rays (UHECR) allow to record not only the extensive air showers, initiated by the UHECR particles, but also to detect light, produced by meteors and by the fast dust grains. It is shown that the fluorescence detector operated at the mountain site can register signals from meteors with kinetic energy threshold of about 25 J (meteor mass ∼ 5 × 10⁻⁶ g, velocity ∼ 3 × 10⁶ cm/s). The same detector might be used for recording of the dust grains of smaller mass (of about 10⁻¹⁰ g) but with velocity 10⁹ cm/s, close to the light velocity (sub-relativistic dust grains). The light signal from a sub-relativistic dust grain is expected in much shorter time scale (∼0.001 s), in comparison with the meteor signal (∼0.1–1 s), and much longer than duration of the UHECR signals (tens of μs). The fluorescence detector capable to register various phenomena: from meteors to UHECR – should have a variable pixel and selecting system integration time. A study of the new phenomenon of sub-relativistic grains will help to understand the mechanism of particle and dust grain acceleration in astrophysical objects (in SN explosions, for example).     

Modeling of effect of polarization on UV sky radiance during twilight

The spatial distribution of the vector of the Stokes parameters characterizing the radiance intensity and the radiance polarization describes the radiation field in the atmosphere. A simplified treatment of light as the scalar has only restricted application. A few studies compared previously results of the vector and scalar radiative transfer models and showed that scalar models are in error by up to 10% for many cases. Though several observational conditions were exploited, an effect of polarization on modeling of UV radiance has not been investigated yet for twilight. The paper presents a preliminary study of modeled UV radiance during twilight taking into account polarization. The intensity and the degree of linear polarization of the scattered UV radiance for two cases of the ground-based observations are discussed. In the first case, radiation incoming from the zenith for the solar zenith angles (SZA) from 90° to 98° is under investigation. Radiation in the solar principal plane for the beginning of twilight (SZA = 90.1°) was calculated in the second case. The study showed that the UV radiation field in the twilight atmosphere can be handled correctly only using the vector theory. The errors of scalar radiative transfer strongly depend on wavelength, line of an observation and solar position. The revealed distortion of the zenith radiance caused by using of the scalar approximation reaches maximum of 15% at 340 nm for the solar zenith angle (SZA) equal to 98°. The shorter wavelengths have the smaller errors, about 5% at 305 nm for SZA = 98°, due to the larger part of the single scattered radiance. The error of the scalar modeling may be as large as −17% for radiance incoming from the horizon for SZA = 90.1°. Scalar radiative transfer models underestimate the integral intensity in the principal plane up to 3–4% ± 0.5% at SZA = 90.1° for wavelengths from 320 to 340 nm. This should be taken into account in problems of radiative budget estimation and remote sensing of the atmosphere exploiting the twilight period.     

Effect of the troposphere on surface neutron counter measurements

Surface neutron counter data are often used as a proxy for atmospheric ionisation from cosmic rays in studies of extraterrestrial effects on climate. Neutron counter instrumentation was developed in the 1950s and relationships between neutron counts, ionisation and meteorological conditions were investigated thoroughly using the techniques available at the time; the analysis can now be extended using modern data. Whilst surface neutron counts are shown to be a good proxy for ionisation rate, the usual meteorological correction applied to surface neutron measurements, using surface atmospheric pressure, does not completely compensate for tropospheric effects on neutron data. Residual correlations remain between neutron counts, atmospheric pressure and geopotential height, obtained from meteorological reanalysis data. These correlations may be caused by variations in the height and temperature of the atmospheric layer at ∼100 hPa. This is where the primary cosmic rays interact with atmospheric air, producing a cascade of secondary ionising particles.     

On some large solar energetic particle events: Direct impulsive component and broad associated coronal mass ejections

Using the proton intensity and X-ray flux data from the GOES, combined with the observations of the associated solar eruptions by the Large Angle and Spectrometric Coronagraph Experiment (LASCO) on board the Solar and Heliospheric Observatory (SOHO), 14 large SEP events occurring in the period 2000 January–2002 April have been studied. It is found that: (1) events with the SEPs increasing shortly after the maximum of their parent flares (<1 h; hereafter prompt events) have rapid and great (up to four orders of magnitude) SEP increments in high-energy channels (> ∼100 MeV); however, for events whose onset of the SEP injection lags the flare maximum for a long time (>3 h; hereafter delayed events), the high-energy SEPs show no obvious enhancements (within one order of magnitude); (2) peak intensity of the prompt events is distinctly larger than that of the delayed events; (3) CMEs associated with the poorly magnetically connected events (source region <W30°) in our survey are all halo CMEs. From these observational differences, we propose a special scenario of the production of the largest SEP events: both CMEs and flares are induced in the same coronal process; high-energy particles accelerated in the reconnection region can escape easily from the open field lines and/or be transported by fast CMEs into interplanetary space, indicating a direct impulsive component in large gradual SEP events. Meanwhile, the broad width of the associated CMEs implies that the CME width is more important in SEP events production than previously considered.     

Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

Launching of a 500,000 cubic meter balloon with the semi-dynamic launching method

Launching a large balloon in a limited launching field is a long standing problem in Japan. The largest balloon ever launched successfully was 200,000 m³ in volume. It was launched in 1973. A larger balloon with a volume of 500,000 m³ was tried later, but it burst during the ascending phase. For launching balloons with a large lift exceeding 500 kg, the conventional static launching method had the most serious problem with possible damage to the polyethylene film of the balloon caused by the holding mechanism. After that, we had developed a new static launching method to launch balloons with a total lift of 1.0 ton. For launching a large balloon with a total lift above 1.5 ton, the new static launching method had a weak point in that if there was an air bubble in the folded part of the balloon, it may puncture the balloon as it is pushed by a spool. To avoid this problem, we developed a semi-dynamic launching method in 1999 using a launcher fixed to the ground leaving a freedom of rotation around the vertical axis. We have launched some balloons using the method and have gradually enriched our experience in using this system.In 2003, we successfully launched a balloon with a volume of 500,000 m³ by using the method. This balloon was made of polyethylene films with a thickness of 20 μm and it is the largest balloon ever launched in Japan.     

A radio study of the mouse,G359.23 − 0.82

The recent detection of a young pulsar powering “the Mouse”, G359.23 − 0.82, as well as detailed imaging of surrounding nebular X-ray emission, have motivated us to investigate the structural details and polarization characteristics of the radio emission from this axisymmetric source with a supersonic bow shock. Using polarization data at 3.6 and 6 cm, we find that the magnetic field wraps around the bow-shock structure near the apex of the system, but downnstream runs parallel to the inferred direction of the pulsar’s motion. The rotation measure (RM) distribution of the Mouse also suggests that the low degree of polarization combined with a high RM ahead of the pulsar result from internal plasma within the bow-shock region. In addition, using sub-arcsecond radio image of the Mouse, we identify modulations in the brightness distribution of the Mouse that may be associated with the unshocked pulsar wind behind the pulsar. Lastly, we discuss the relationship between the Mouse and its neighboring shell-type supernova remnant G359.1 − 0.5 and argue that these two sources could potentially have the same origin.