Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Rotational constants of linear and/or bent Cn+1H and CnN(n = 1–6): A DFT study

The geometries, dipole moments, and rotational constants for the linear and/or bent cations, C_n+1H⁺ and C_nN⁺(n = 1–6), were studied by the B3LYP method with the modest basis sets. For C_nH⁺(n = odd; 3, 5, 7) and C_nN⁺(n = even; 2, 4, 6), the theoretical rotational constants (B_es) of closed-shell singlet C₃H⁺, C₅H⁺, C₇H⁺, CCN⁺, C₄N⁺, and C₆N⁺ were calculated to be about 11,244, 2420, 885.2, 11,970, 2439, and 880.8 MHz, respectively. By contrast, the triplets are stable than the corresponding singlets for C_nH⁺(n = odd; 2, 4, 6) and C_nN⁺(n = even; 3, 5) except CN⁺.     

Model calculations of the dayside ionosphere of Venus

Model calculations of the dayside ionosphere of Venus are presented. The coupled continuity and momentum equations were solved for O₂⁺, O⁺, CO₂⁺, C⁺, N⁺, He⁺, and H⁺ density distributions, which are compared with measurements from the Pioneer Venus ion mass spectrometer. The agreement between the model results and the measurements is good for some species, such as O⁺, and rather poor for others, such as N⁺, indicating that our understanding of the dayside ion composition of Venus is incomplete. The coupled heat conduction equations for ions and electrons were solved and the calculated temperatures compared with Pioneer Venus measurements. It is shown that fluctuations in the magnetic field have a significant effect on the energy balance of the ionosphere.     

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.     

Relativistic electron fluxes and dose rate variations during April–May 2010 geomagnetic disturbances in the R3DR data on ISS

Space radiation has been monitored successfully using the Radiation Risks Radiometer-Dosimeter (R3D) installed at the ESA EXPOSE-R (R3DR) facility outside of the Russian Zvezda module of the International Space Station (ISS) between March 2009 and January 2011. R3DR is a Liulin type spectrometer–dosimeter with a single Si PIN detector 2 cm² of area and 0.3 mm thick. The R3DR instrument accumulated about 2 million measurements of the absorbed dose rate and flux of 10 s resolution. The total external and internal shielding before the detector of R3DR device is 0.41 g cm⁻². The calculated stopping energy of normally incident particles to the detector is 0.78 MeV for electrons and 15.8 MeV for protons. After the Coronal Mass Ejection (CME) at 09:54 UTC on 3 April 2010, a shock was observed at the ACE spacecraft at 0756 UTC on 5 April, which led to a sudden impulse on Earth at 08:26 UTC. Nevertheless, while the magnetic substorms on 5 and 6 of April were moderate; the second largest in history of GOES fluence of electrons with energy >2 MeV was measured. The R3DR data show a relatively small amount of relativistic electrons on 5 April. The maximum dose rate of 2323 μGy day⁻¹ was reached on 7 April; by 9 April, a dose of 6600 μGy was accumulated. By the end of the period on 7 May 2010 a total dose of 11,587 μGy was absorbed. Our data were compared with AE-8 MIN, CRESS and ESA-SEE1 models using SPENVIS and with similar observations on American, Japanese and Russian satellites.     

Influence of solar-geomagnetic disturbances on SABER measurements of 4.3 μm emission and the retrieval of kinetic temperature and carbon dioxide

Thermospheric infrared radiance at 4.3 μm is susceptible to the influence of solar-geomagnetic disturbances. Ionization processes followed by ion-neutral chemical reactions lead to vibrationally excited NO⁺ (i.e., NO⁺(v)) and subsequent 4.3 μm emission in the ionospheric E-region. Large enhancements of nighttime 4.3 μm emission were observed by the TIMED/SABER instrument during the April 2002 and October–November 2003 solar storms. Global measurements of infrared 4.3 μm emission provide an excellent proxy to observe the nighttime E-region response to auroral dosing and to conduct a detailed study of E-region ion-neutral chemistry and energy transfer mechanisms. Furthermore, we find that photoionization processes followed by ion-neutral reactions during quiescent, daytime conditions increase the NO⁺ concentration enough to introduce biases in the TIMED/SABER operational processing of kinetic temperature and CO₂ data, with the largest effect at summer solstice. In this paper, we discuss solar storm enhancements of 4.3 μm emission observed from SABER and assess the impact of NO⁺(v) 4.3 μm emission on quiescent, daytime retrievals of Tk/CO₂ from the SABER instrument.     

An analytical model for the prediction of a micro-dosimeter response function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (δ ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCRs) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 s/ion. The ionizing radiation environment at LEO is represented by O’Neill’s GCR model (2004), covering charged particles in the 1 ? Z ? 28 range. O’Neill’s free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.     

A detailed FLUKA-2005 Monte-Carlo simulation for the ATIC detector

We have performed a detailed Monte-Carlo (MC) simulation for the Advanced Thin Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2005 which is capable of simulating particles up to 10 PeV. The ATIC detector has completed two successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 days. ATIC is designed as a multiple, long duration balloon flight, investigation of the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully active Bismuth Germanate (BGO) calorimeter. It is equipped with a large mosaic of silicon detector pixels capable of charge identification, and, for particle tracking, three projective layers of x–y scintillator hodoscopes, located above, in the middle and below a 0.75 nuclear interaction length graphite target. Our simulations are part of an analysis package of both nuclear (A) and energy dependences for different nuclei interacting in the ATIC detector. The MC simulates the response of different components of the detector such as the Si-matrix, the scintillator hodoscopes and the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2005 MC calculations with GEANT calculations and with the ATIC CERN data.     

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).     

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

AVS-F observations of γ-ray emission during January 20, 2005 solar flare up to 140 MeV

The solar flare of January 20, 2005 (X7.1, 06:36–07:26 UT, maximum at 07:01 UT by the GOES soft X-ray data) was the most powerful one in January 2005 series. The AVS-F apparatus onboard CORONAS-F registered γ-emission during soft X-ray rising phase of this flare in two energy ranges of 0.1–20 MeV and 2–140 MeV. The highest γ-ray energy registered during this flare was ∼140 MeV. Six spectral features were registered in energy spectrum of this solar flare: annihilation + αα (0.4–0.6 MeV), ²⁴Mg + ²⁰Ne + ²⁸Si + neutron capture (1.7–2.3 MeV), ²¹Ne + ²²Ne + ¹⁶O + ¹²С (3.2–5.0 MeV), ¹⁶O (5.3–6.9 MeV), one from neutral pions decay (25–110 MeV) and one in energy band 15–21 MeV. Four of them contain typical for solar flares lines – annihilation, nuclear de-excitation and neutron capture at ¹H. Spectral feature caused by neutral pions decay was registered during several flares too. Some spectral peculiarities in the region of 15–21 MeV were first observed in solar flare energy spectrum.     

Stochastic particle acceleration by the forced interaction of relativistic current sheets

We address the problem of interacting relativistic current sheets in self-consistent kinetic plasma simulations within the framework of the Particle-In-Cell model. The interaction is enforced in head-on collisions of up to 10 current sheets at relativistic bulk speeds. The simulations are motivated by the general problem of Poynting flux dissipation in ‘striped wind’ configurations presumably governing the relativistic outflows pervasive in pulsar winds and gamma-ray bursts. We identify the generation of non-thermal particles and formation of a stable power-law shape in the particle energy distributions f(γ) dγ ∝ γ^−s dγ. In 1D, a spectral index s ∼ 2 is observed and attributed to a stochastic Fermi-type acceleration mechanism. In 2D, the generic index of s ∼ 3–4 is retained as in previous simulations of individual current sheets. Whereas in 2D the high energy cut-off is constrained by the limited dissipation of magnetic energy, in 1D the process converts the bulk motion of current sheets towards directed particle momentum of an exclusive class of non-thermal particles.     

F-Region ion composition modeling

Measurements of the principal ion species of the F₁- and F₂- regions have been used to develop an empirical model of the ion composition for altitudes between 150 and 500 km. The species measured by the S3-1 satellite include N⁺, O⁺, N₂⁺, NO⁺ and O₂⁺. The data were obtained near the minimum of the solar cycle, thus limited information on the ionospheric variation with solar flux is available. However, the range of latitude, altitude, local time and geomagnetic activity does provide a useful basis for modeling the F-region. The ion composition measurements have been used to provide a model for relative ion composition which is compatible with the total ion density from the International Reference Ionosphere model.     

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.     

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.     

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.     

Feasibility of developing an ionospheric E-region electron density storm model using TIMED/SABER measurements

We present a new technique for improving ionospheric models of nighttime E-region electron densities under geomagnetic storm conditions using TIMED/SABER measurements of broadband 4.3 μm limb radiance. The response of E-region electron densities to geomagnetic activity is characterized by SABER-derived NO⁺(v) 4.3 μm Volume Emission Rates (VER). A storm-time E-region electron density correction factor is defined as the ratio of storm-enhanced NO⁺(v) VER to a quiet-time climatological average NO⁺(v) VER, which will be fit to a geomagnetic activity index in a future work. The purpose of this paper is to demonstrate the feasibility of our technique in two ways. One, we compare storm-to-quiet ratios of SABER-derived NO⁺(v) VER with storm-to-quiet ratios of electron densities measured by Incoherent Scatter Radar. Two, we demonstrate that NO⁺(v) VER can be parameterized by widely available geomagnetic activity indices. The storm-time correction derived from NO⁺(v) VER is applicable at high-latitudes.