Status of the GAMMA-400 project

The preliminary design of the new space gamma-ray telescope GAMMA-400 for the energy range 100 MeV–3 TeV is presented. The angular resolution of the instrument, 1–2° at E_γ ∼ 100 MeV and ∼0.01° at E_γ > 100 GeV, its energy resolution ∼1% at E_γ > 100 GeV, and the proton rejection factor ∼10⁶ are optimized to address a broad range of science topics, such as search for signatures of dark matter, studies of Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission, gamma-ray bursts, as well as high-precision measurements of spectra of cosmic-ray electrons, positrons, and nuclei.     

Using gamma-ray burst prompt emission to probe relativistic shock acceleration

It is widely accepted that the prompt transient signal in the 10 keV–10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broad-band spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).     

The ground-based large-area wide-angle γ-ray and cosmic-ray experiment HiSCORE

The question of the origin of cosmic rays and other questions of astroparticle and particle physics can be addressed with indirect air-shower observations above 10 TeV primary energy. We propose to explore the cosmic ray and γ-ray sky (accelerator sky) in the energy range from 10 TeV to 1 EeV with the new ground-based large-area wide angle (ΔΩ ∼ 0.85 sterad) air-shower detector HiSCORE (Hundred^∗i Square-km Cosmic ORigin Explorer). The HiSCORE detector is based on non-imaging air-shower Cherenkov light-front sampling using an array of light-collecting stations. A full detector simulation and basic reconstruction algorithms have been used to assess the performance of HiSCORE. First prototype studies for different hardware components of the detector array have been carried out. The resulting sensitivity of HiSCORE to γ-rays will be comparable to CTA at 50 TeV and will extend the sensitive energy range for γ-rays up to the PeV regime. HiSCORE will also be sensitive to charged cosmic rays between 100 TeV and 1 EeV.     

Study of the variability of blazars gamma-ray emission

The γ-ray emission of blazar jets shows a pronounced variability and this feature provides limits to the size and to the speed of the emitting region. We study the γ-ray variability of bright blazars using data from the first 18 months of activity of the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. From the daily light-curves of the blazars characterized by a remarkable activity, we firstly determine the minimum variability time-scale, giving an upper limit for the size of the emitting region of the sources, assumed to be spheroidal blobs in relativistic motion. These regions must be smaller than ∼10⁻³ parsec. Another interesting time-scale is the duration of the outbursts. We conclude that they cannot correspond to radiation produced by a single blob moving relativistically along the jet, but they are either the signature of emission from a standing shock extracting energy from a modulated jet, or the superposition of a number of flares occurring on a shorter time-scale. We also derive lower limits on the bulk Lorentz factor needed to make the emitting region transparent for gamma-rays interacting through photon–photon collisions.     

Acceleration of charged particles by fluctuating and steady electric fields in a X-type magnetic field

Release of stored magnetic energy via particle acceleration is a characteristic feature of astrophysical plasmas. Magnetic reconnection is one of the mechanisms for releasing energy from magnetized plasmas. Collisionless magnetic reconnection could provide both the energy release mechanism and the particle accelerator in space plasmas. Here we studied particle acceleration when fluctuating (in-time) electric fields are superposed on an static X-type magnetic field in collisionless hot solar plasma. This system is chosen to mimic the reconnective dissipation of a linear MHD disturbance. Our results are compared to particle acceleration from constant electric field superposed on an X-type magnetic field. The constant electric field configuration represents the effects of steady state magnetic reconnection. Time evolution of ion and electron distributions are obtained by numerically integrating particle trajectories. The frequencies of the electric field represent a turbulent range of waves. Depending on the frequency and amplitude of the electric field, electrons and ions are accelerated to different degrees and have energy distributions of bimodal form consisting of a lower energy part and a high energy tail. For frequencies (ω in dimensioless units) in the range 0.5 ? ω ? 1.0 a substantial fraction (20%–30%) of the proton distribution is accelerated to gamma-ray producing energies. For frequencies in the range 1 ? ω ? 100.0 the bulk of the electron distribution is accelerated to hard X-ray producing energies. The acceleration mechanism is important for solar flares and solar noise storms but it could be applicable to all collisionless astrophysical plasmas.     

Spectral analysis of time-integrated Konus–Wind GRBs: Implication on radiative mechanisms

We have analysed a sample of 328 time-integrated GRB prompt emission spectra taken via the Konus instrument on board the US GGS-Wind spacecraft between 2002 and 2004 using a couple of two-components models, Cut-off Power Law (CPL) + Power Law (PL) and blackbody (BB) + PL. The spectra show clear deviation from the Band function. The PL term is interpreted as the low energy tail of a nonthermal emission mechanism. The distributions of corresponding index β give values β < −2/3 consistent with synchrotron and synchrotron self-Compton mechanisms. The distribution of low energy index α associated with the CPL term shows clear discordance with synchrotron models for 31.4% of the analysed GRBs with values exceeding that for the line of death, α = −2/3. Then, a set of nonthermal radiation mechanisms producing harder slopes, i.e., α > −2/3, are presented and discussed. For the remaining majority (68.6%) of GRBs with CPL index α < −2/3, we show that optically thin synchrotron produced by a power law electron distribution of type, N(γ) ∼ γ^−p, γ₁ < γ < γ₂, for finite energy range (γ₂ ≠ ∞) is a likely emission mechanism with α ∼−(p + 1)/2 in the frequency range ν₁ ? ν ? ν₂ (where ν₂ = η²ν₁ with η = γ₂/γ₁), such that for p > 1/3, one gets α < −2/3. We also show that corresponding spectra in terms of F_ν and νF_ν functions are peaked around frequency ν₂ instead of ν₁, respectively for p < 1 and p < 3. Besides, thermal emission is examined taking a single Planck function for fitting the low energy range. It can be interpreted as an early emission from the GRB fireball photosphere with observed mean temperature, kT′ ∼ 16.8 keV. Furthermore, we have performed a statistical comparison between the CPL + PL and BB + PL models finding comparable χ²-values for an important fraction of GRBs, which makes it difficult to distinguish which model and specific radiation mechanism (possible thermal or nonthermal γ-ray emissions) are best suitable for describing the reported data. Therefore, additional information for those bursts, such as γ-ray polarization, would be highly desirable in future determinations of GRBs observational data.     

Electron temperature anisotropy effects on tearing mode in ion-scale current sheets

Recent two-dimensional (2-D) particle-in-cell (PIC) simulations have shown that there is a critical thickness of a current sheet, above which no significant saturation amplitude of the 2-D tearing (TI) mode can be expected. Here, we have introduced the initial electron temperature anisotropy (α_e0 = T_e⊥/T_e|| > 1), which is known to raise significantly the linear growth rates, and inspected if α_e0 > 1 can change the saturation level of the TI in a super-critical current sheet. Varying α_e0 and D (D: the current sheet half-thickness) systematically, we have found that while α_e0 boosts up the linear growth rate in both sub- and super-critical current sheets, macroscopic effects are obtained only in sub-critical current sheets, that is, energy transfer from the fastest growing short wavelength modes to longer wavelength modes are available only in the sub-critical regime. Since the critical thickness is a fraction of the ion inertial length, the tearing mode assisted by the electron temperature anisotropy alone, despite its significant boost in the linear growth rate, cannot be the agent for reconnection triggering in a current sheet of ion-scale thickness.     

Solar cycle and latitude dependence of high-beta suprathermal plasma conditions in interplanetary space between 1.3 and 5.4 AU

The analysis of energetic particles and magnetic field measurements from the Ulysses spacecraft has shown that in a series of events, the energy density contained in the suprathermal tail particle distribution is comparable to or larger than that of the magnetic field, creating conditions of high-beta plasma. In this work we analyze periods of high-beta suprathermal plasma occurrences (β_ep > 1) in interplanetary space, using the ratio (β_ep) of the energetic particle (20 keV to ∼5 MeV) and magnetic field energy densities from measurements covering the entire Ulysses mission lifetime (1990–2009) in order to reveal new or to reconfirm some recently defined interesting characteristics. The main key-results of the work are summarized as follows: (i) we verify that high-beta events are detected within well identified regions corresponding mainly to the vicinity of shock surfaces and magnetic structures, and associated with energetic particle intensity enhancements due to (a) reacceleration at shock-fronts and (b) unusually large magnetic field depressions. (ii) We define three considerable features for the high-beta events, concentrated on the next points: (a) there is an appreciable solar-activity influence on the high-beta events, during the maximum and middle solar-cycle phase, (b) the annual peak magnitude and the number of occurrences of high events are well correlated with the sunspot number, (c) the high-beta suprathermal plasma events present a spatial distribution in heliographic latitudes (HL) up to ∼±80°, and a specific important concentration on the low (−25° ? HL < −6°, 6° < HL ? 25°) and median (−45° ? HL < −25°, 25° < HL ? 45°) latitudes. We also reconfirm by a statistical analysis the results of Marhavilas and Sarris (2011), that the high-beta suprathermal plasma (β_ep > 1) events are characterized by a very large parameter β_ep (up to 1732.5), a great total duration (406 days) and a large percentage of the Ulysses-mission lifetime (which is equal to 6.34% of the total duration with usable measurements, and 11.3% of the duration in presence of suprathermal particles events).     

On the electrodynamics of counter stream in pulsar’s inner annular cap region

The electric properties of pulsar’s inner annular gap are explored in this paper. Under two main assumptions, (1) the pulsar is alive, (2) the total charge of pulsar should not vary with time, the condition for the acceleration of negative particle in the annular region is derived. The acceleration condition is j₋ ? 0.5j₊, i.e., the current carried by negative particles is greater than or equal to 0.5 times of the current carried by positive particles. This condition holds even when the backward flow of positive particles exists in the annular region. It is noted that the outflow of negative particles offers good opportunities to understand the current closure problem of pulsar as well as wide radiation beam of pulsar observed at high energy band.     

Results from the study of an M7.6 flare and its associated CME

An M7.6 flare was well observed on October 24, 2003 in active region 10486 by a few instruments and satellites, including GOES, TRACE, SOHO, RHESSI and NoRH. Multi-wavelength study shows that this flare underwent two episodes. During the first episode, only a looptop source of <40 keV was observed in reconstructed RHESSI images, which showed shrinkage with a velocity of 12–14 km s⁻¹ in a period of about 12 min. During the second process, in addition to the looptop source, two footpoint sources were observed in energy channel of as high as ∼200 keV. One of them showed fast propagation along one of the two TRACE 1600 Å flare ribbons and the 195 Å loop footpoints, which could be explained by successive magnetic reconnection. The associated CME showed a mass pickup process with decreasing center-of-mass velocity. The decrease of the CME kinetic energy and the increase of its potential energy lead to an almost constant total energy during the CME propagation. Our results reveal that the flare and its associated CME have comparable energy content, and the flare is of non-thermal property.     

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

Statistical study of the detection of solar protons of highest energies at 20 January 2005

On January 20, 2005, 7:02–7:04 UT the Aragats Multichannel Muon Monitor (AMMM) registered enhancement of the high energy secondary muon flux (energy threshold ∼5 GeV). The enhancement, lasting 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES satellite and the ground level enhancement detected by the world-wide network of neutron monitors and by muon detectors. The most probable proton energy corresponding to the measured 5 GeV muon flux is within 23–30 GeV. Due to upmost importance of the detection of solar particles of highest energies in presented paper we perform detailed statistical analysis of the detected peak. The statistical technique introduced in the paper is also appropriate for the searches of sources of ultra-high energy cosmic rays.     

Jovian periodicities (∼10 h, ∼40 min) on Ulysses’ Distant Jupiter Encounter observations around the Halloween CIR events

We analyzed data from four different instruments (HI-SCALE, URAP, SWOOPS, VHM/FGM) onboard Ulysses spacecraft (s/c) and we searched for possible evidence of Jovian emissions when the s/c approached Jupiter during the times of Halloween events (closest time approach/position to Jupiter: February 5, 2004/R = 1683 R_J,θ = ∼49°). In particular, we analyzed extensively the low energy ion measurements obtained by the HI-SCALE experiment in order to examine whether low energy ion/electron emissions show a symmetry, and whether they are observed at north high latitudes upstream from the jovian bow shock, as is known to occur in the region upstream from the south bow shock as well ( Marhavilas et al., 2001). We studied the period from October 2003 to March 2004, as Ulysses moved at distances 0.8–1.2 AU from the planet at north Jovicentric latitudes <75°, and we present here an example of characteristic Jovian periodicities in the measurements around a CIR observed by Ulysses on days ∼348–349/2003 (R = 1894 R_J,θ = 72°). We show that Ulysses observed low energy ion (∼0.055–4.7 MeV) and electron (>∼40 keV) flux and/or spectral modulation with the Jupiter rotation period (∼10 h) as well as variations with the same period in solar wind parameters, radio and magnetic field directional data. In addition, characteristic strong ∼40 min periodic variations were found superimposed on the ∼10 h ion spectral modulation. Both the ∼10 h and ∼40 min ion periodicities in HI-SCALE measurements were present in several cases during the whole period examined (October 2003 to March 2004) and were found to be more evident during some special conditions, for instance during enhanced fluxes around the start (forward shock) and the end (reverse shock) of CIRs. We infer that the Jovian magnetosphere was triggered by the impact of the CIRs, after the Halloween events, and it was (a) a principal source of forward and reverse shock-associated ion flux structures and (b) the cause of generation of ∼10 h quasi-periodic magnetic field and plasma modulation observed by Ulysses at those times.     

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.     

Quantifying the role of solar radiative forcing over the 20th century

The 20th century temperature anomaly record is reproduced using an energy balance model, with a diffusive deep ocean. The model takes into account all the standard radiative forcings, and in addition the possibility of a non-thermal solar component. The model is parameterized and then optimized to produce the most likely values for the climate parameters and radiative forcings which reproduce the 20th century global warming. We find that the best fit is obtained with a negligible net feedback. We also show that a non-thermal solar component is necessarily present, indicating that the total solar contribution to the 20th century global warming, of ΔT_solar = 0.27 ± 0.07 °C, is much larger than can be expected from variation in the total solar irradiance alone. However, we also find that the largest contribution to the 20th century warming comes from anthropogenic sources, with ΔT_man = 0.42 ± 0.11 °C.     

Energy conversion of the flare due to direct electric fields from the sheared reconnection

In this paper we present a new mechanism of the main energy conversion of the solar flare. Since a flare inducing prominence (flux tube) rises V_z ? 300 km s⁻¹, the plasmas below it cannot continuously eject with Alfvén speeds of V_A = 3000 km s⁻¹ but probably with V_z ≈ ±100 km s⁻¹. Plasma up and downflows with V_A will within a short duration be blocked between the chromosphere where reconnected flux tubes are piling up, and the slowly rising flux rope. Hence the Petschek slow shock mechanism is difficult to be realized as a major energy converting mechanism.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

The hydroxyl radical as an indicator of SEP fluxes in the high-latitude terrestrial atmosphere

The low background values at nighttime of the mesospheric hydroxyl (OH) radical make it easier to single out the atmospheric response to the external solar forcing in Polar Regions. Because of the short lifetime of HO_x, it is possible to follow the trails of Solar Energetic Particle (SEP) events in the terrestrial atmosphere, as shown by Storini and Damiani (2008). The sensitivity of this indicator makes discernible not only extreme particle events with a flux peak of several thousand pfu [1 pfu = 1 particle/(cm² s sr)] at energies >10 MeV, but also those with lower flux up to about 300 pfu. Using data from the Microwave Limb Sounder (MLS) on board the EOS AURA satellite, we examined the correlation of OH abundance vs. solar proton flux for almost all the identified SEP events spanning from November 2004 to December 2006 (later on no more SEP events occurred during Solar Cycle no. 23). The channels at energies greater than 5 MeV and 10 MeV showed the best correlation values (r ∼ 0.90–0.95) at altitudes around 65–75 km whereas, as expected, the most energetic channels were most highly correlated at lower altitudes. Therefore, it is reasonably possible to estimate the solar proton flux from values of mesospheric OH (and viceversa) and it could be useful in studying periods with gaps in the records of solar particles.     

X-ray study of accretion flow in narrow-line Seyfert 1 galaxies

We present the results of a systematic study of narrow-line Seyfert 1 galaxies (NLS1s) observed with XMM-Newton. The 2–12 keV X-ray spectra of NLS1s are well represented by a single power law with a photon index Γ ∼ 2. When this hard power law continuum is extrapolated into the low energy band, we found that all objects in our sample show prominent soft excess emission. This excess emission is well parameterized by the thermal emission expected from an optically thick accretion disk, and we found the following three peculiar features: (1) The derived disk temperatures are significantly higher than the expectation from a standard Shakura-Sunyaev accretion disk, if we assume a central mass of a black hole to be 10^6–8M_⊙. (2) The temperatures are distributed within narrow range (ΔkT ∼ 0.08 keV) with an average temperature of 0.18 keV in spite of the range of four orders of magnitude in luminosity (10^41–45 erg s⁻¹). (3) We found a peculiar temperature–luminosity relation, where the luminosity seems to be almost saturated in spite of the significant change in temperature, during the observations of the most luminous NLS1 PKS 0558-504. These results strongly suggest that the standard accretion disk picture is no longer appropriate in the nuclei of NLS1s. We discuss a possible origin for the soft excess component, and suggest that a slim disk may be able to explain the observational results, if the photon trapping effect is properly taken into account.