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.     

Possible consequences of a disk around B0656+14 on e–e and near-1 Hz gravitational wave production

Nearby pulsars B0656+14 and Geminga were proposed in the literature as the main sources of cosmic-ray positrons observed near Earth above 10 GeV. B0656+14 has comparable distance from Earth, similar magnetic field and period of Geminga. However, observations in the R and I bands indicate the presence of a disk of approximately 10⁻⁴ M_⊙ around B0656+14. Radio and pulsed γ-ray flux observations from this pulsar are also consistent with supernova fallback material and disk entering the light cylinder and partially quenching the development of electromagnetic showers in the magnetosphere. If this is the case, B0656+14 has unlikely given any contribution to e⁺ and e⁻ observed near Earth. Absolute flux measurements and the level of anisotropy in the high energy electron and positron arrival directions above 50 GeV will help in revealing if none, one of both nearby pulsars are sources of these particles observed near Earth.     

Time dependent charging of layer clouds in the global electric circuit

There is much observational data consistent with the hypothesis that the ionosphere-earth current density (J_z) in the global electric circuit, which is modulated by both solar activity and thunderstorm activity, affects atmospheric dynamics and cloud cover. One candidate mechanism involves J_z causing the accumulation of space charge on droplets and aerosol particles, that affects the rate of scavenging of the latter, notably those of Cloud Condensation Nuclei (CCN) and Ice Forming Nuclei (IFN) (  and ). Space charge is the difference, per unit volume, between total positive and total negative electrical charge that is on droplets, aerosol particles (including the CCN and IFN) and air ions. The cumulative effects of the scavenging in stratiform clouds and aerosol layers in an air mass over the lifetime of the aerosol particles of 1–10 days affects the concentration and size distribution of the CCN, so that in subsequent episodes of cloud formation (including deep convective clouds) there can be effects on droplet size distribution, coagulation, precipitation processes, and even storm dynamics.     

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.     

XMM-Newton observations across the Cygnus Loop from northeastern rim to southwestern rim

We have observed the Cygnus Loop from the northeast (NE) to the southwest (SW) with XMM-Newton. We extracted spectra from 3′-spaced annular regions across the Loop and fitted them either with a one-kT_e-component non-equilibrium ionization (NEI) model or with two-kT_e-component NEI model. We found that the two-kT_e-component model yields significantly better fits in almost all the spectra than the one-kT_e-component model. Judging from the abundances, the high-kT_e-component in the two-temperature model must be fossil ejecta while the low-kT_e-component comes from the swept-up interstellar medium (ISM). The distributions of Ne, Mg, Si, and S for fossil ejecta appear to retain the onion-skin structure at the time of a supernova explosion, suggesting that significant overturning of the ejecta has not occurred yet. Comparing the relative abundances of fossil ejecta estimated for the entire Cygnus Loop with those from theoretical calculations for Type-II SN, the mass of the progenitor star is likely to be ∼13 M_⊙. The trends of the radial profiles of kT_e and emission integral for the swept-up ISM are adequately described by the Sedov model, suggesting that the swept-up ISM is concentrated in a shell-like structure. Comparing our data with the Sedov model, we found the ambient medium density to be ∼0.7 cm⁻³. Then, we estimated the total mass of the swept-up ISM and the age of the remnant to be ∼130 M_⊙ and 13,000 years, respectively. Note that if the explosion occurred within a stellar wind cavity, then the density in the cavity, the total swept-up mass in the cavity, and the age of the remnant are estimated to be ∼0.14 cm⁻³, ∼25 M_⊙, and ∼10,000 years, respectively.     

Plasma electron temperature variability in lunar surface potential and in electric field under average solar wind conditions

The Moon is immersed in plasma environment. The most interesting challenge of the lunar plasma– field environment is that it is alternatively dominated by the extended but variable outer atmosphere of the Earth – the magnetosphere – and by the extended but highly variable solar atmosphere – the solar wind. Understanding the plasma environment and its interaction with the lunar surface will be beneficial to both manned and robotic surface exploration activities and to scientific investigations. Presented is a preliminary map of variations of lunar surface electric potential over the day side and night side using probe equations and a discussion on dust dynamics in this E-field structure using the data from Electron Reflectometer in Lunar Prospector spacecraft during 1998–1999. On the day side, potential is around 5 V and on the night side it reaches up to −82 V. On the night side region, only highly energetic electrons can overcome this large negative potential. The variation at electron temperature (T_e) strongly reflects in the surface potential. The potential reaches to a value of −82 V for T_e = 58 eV. Surface charging causes the electrostatic transport of charged dust grains. Dust grain size of 0.1 μm shows a levitation height of 4.92 m on lunar day side, 748 m on terminator region and 3.7 km on the night side. The radius of maximum sized grain to be lofted, R_max, peaks at the terminator region (R_max = 0.83 μm). At the transition region dust levitation is almost absent. This region is most suited for exploration activities as the region is free from hazards caused by lunar dust.     

Particle acceleration by stochastic fluctuations and dawn-dusk electric field in the Earth’s magnetotail

This work is devoted to investigate the interaction between protons and stochastic time-dependent electromagnetic fields generated by oscillating clouds of finite size, randomly positioned in the x–y  -plane. The geometry of the system is two-dimensional and, beside the time-dependent electromagnetic fluctuations, a steady-state, dawn-dusk electric field, _Ey$E_y$, has been added along the y-direction. The simultaneous presence of the stochastic time-dependent fluctuations and of the constant electric field component in the same system gives rise to two types of acceleration mechanisms operating on test particles: a second order Fermi-like process and a direct acceleration. By performing a parametric study, we extensively study the contribution of the two processes to proton acceleration. The energy values reached by test particles in this simple model are in good agreement with those observed in the Earth’s magnetotail region. Possible applications to the problem of particle acceleration in the terrestrial magnetosphere are widely discussed, and guidelines for future works are drawn.     

The optical emission lines of type 1 X-ray bright Active Galactic Nuclei

A strong X-ray emission is one of the defining signatures of nuclear activity in galaxies. According to the Unified Model for Active Galactic Nuclei (AGN), both the X-ray radiation and the prominent broad emission lines, characterizing the optical and UV spectra of Type 1 AGNs, are originated in the innermost regions of the sources, close to the Super Massive Black Holes (SMBH), which power the central engine. Since the emission is concentrated in a very compact region (with typical size r?0.1$r?0.1$ pc) and it is not possible to obtain resolved images of the source, spectroscopic studies of this radiation represent the only valuable key to constrain the physical properties of matter and its structure in the center of active galaxies. Based on previous studies on the physics of the Broad Line Region (BLR) and on the X-ray spectra of broad (FWHM_Hβ ? 2000 km s⁻¹) and narrow line (1000 km s⁻¹ ?FWHM_Hβ ? 2000 km s⁻¹) emitting objects, it has been observed that the kinematic and ionization properties of matter close to the SMBHs are related together, and, in particular, that ionization is higher in narrow line sources. Here we report on the study of the optical and X-ray spectra of a sample of Type 1 AGNs, selected from the Sloan Digital Sky Survey (SDSS) database, within an upper redshift limit of z=0.35$z=0.35$, and detected at X-ray energies. We present analysis of the broad emission line fluxes and profiles, as well as the properties of the X-ray continuum and Fe Kα emission and we use these parameters to assess the consistency of our current AGN understanding.     

Theoretical properties of offset bipolar electric field solitary structures in space plasmas

The bipolar electric field solitary (EFS) structures have been frequently observed in the near Earth plasma regions, such as auroral zone, magnetopause, cusp regions, and magneto-tail. Sometimes these structures are observed as offset bipolar structures. In this paper, the properties of the offset bipolar EFS structures parallel to the magnetic field are studied with an ion fluid model in a cylindrical symmetry by considering electrostatic condition. The model results show that the offset bipolar EFS structures can develop from both ion-acoustic waves and ion cyclotron waves, and propagate along the magnetic field line in the space plasmas if plasma satisfies some conditions. The offset bipolar EFS structures can have both polarities. It will be first negative pulse and then positive pulse if the initial electric field E₀ < 0 or reverse in order if E₀ > 0. The amplitude of the offset bipolar EFS structures first decreases and then increases with the wave propagation velocity. Some results from our model are consistent with the observations.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

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.     

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.     

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.     

F3 layer feature under low and medium solar activity observed at a Chinese low latitude station Fuke

Ionogram observations from the ionosonde at Fuke (9.5°N geomagnetic latitude), a Chinese low latitude station, in 2010–2012 are analyzed to present the features of F3 layer under low and moderate solar activity conditions. Structure of the ionogram, displaying the F3 layer, was more distinct and clear during MSA than LSA periods especially during spring to summer. Start time of occurrence of the F3 layer is about at 0830–0900 LT and is approximately the same for LSA and MSA conditions. The average duration time of the F3 layer occurrence was 181 min per day under F_10.7 = 75 condition, 263 min in F_10.7 = 99 and 358 min in F_10.7 = 125, respectively. The differences of h′F2 and h′F3 exhibited obvious semiannual variation observed at Fuke from March 2010 to June 2012 and increased with increasing solar activity. The difference of foF2 and foF3 in the months February, March, September, October and November is less evident in the middle solar activity period 2011–2012 than the low solar activity 2010 and in the other period it shows a slight increase (0.5 MHz) or keeps constant. The results show that the solar activity dependence of the F3 layer occurrence at low latitude away from the magnetic equator is different from that at near the magnetic equator.     

Interpretation of the SOHO/UVCS observations of two CME-driven shocks

We report on the analysis of two fast CME-driven shocks observed with the UltraViolet Coronagraph Spectrometer (UVCS) on board the Solar and Heliospheric Observatory (SOHO). The first event, detected on 2002 March 22 at 4.1 R_⊙ with the UVCS slit placed in correspondence with the flank of the expanding CME surface, represents the highest UV detection of a shock obtained so far with the UVCS instrument in the corona. The second one, detected on 2002 July 23 at 1.6 R_⊙ with the UVCS slit placed in correspondence with the front of the expanding CME surface, shows an anomalous deficiency of ion heating with respect to what observed in previous CME/shocks observed by UVCS, possibly reflecting the effect of different coronal plasma conditions over the solar cycle. From the two different sets of observations we derived an estimate for the shock compression ratio X, which turns out to be X = 2.4 ± 0.2 and X = 2.2 ± 0.1, respectively, for the first and second event. Comparison between the two events provides complementary perspectives on the dynamical evolution of CME-driven shocks.     

X-ray redshifts with the International X-ray Observatory (IXO)

We explore the capabilities of the future space science mission IXO (International X-ray Observatory) for obtaining cosmological redshifts of distant Active Galactic Nuclei (AGNs) using the X-ray data only. We first find in which regions of the X-ray luminosity (L_X) versus redshift (z) plane the weak but ubiquitous Fe Kα narrow emission line can deliver an accurate redshift (δz < 5%) as a function of exposure time, using a CCD-based Wide Field Imager (IXO/WFI) as the one baselined for IXO. Down to a 2–10 keV X-ray flux of 10⁻¹⁴ erg cm⁻² s⁻¹ IXO/WFI exposures of 100 ks, 300 ks and 1 Ms will deliver 20%, 40% and 60% of the redshifts. This means that in a typical 18′ × 18′ IXO/WFI field of view, 4, 10 and 25 redshifts will be obtained for free from the X-ray data alone, spanning a wide range up to z ∼ 2–3 and fairly sampling the real distribution. Measuring redshifts of fainter sources will indeed need spectroscopy at other wavebands.     

High energy charged particle experiment – A payload proposal to KuaFu-B mission

High Energy Charged Particle Experiment (HECPE) is to measure the fluxes of MeV electrons and tens of MeV protons. The two satellites of KuaFu-B are in the same polar orbit with apogee 7.0R_E, perigee 1.8R_E. They can sweep large L values and pass through the inner and outer radiation belts. The high energy electrons and protons in the radiation belts are principal sources for failures of satellites and spacecrafts in the Earth orbits. The enhancements of the high energy electrons and protons, so-called energetic particle events, are important phenomena of the Space Weather. The energy ranges monitored by HECPE are 0.3–0.5 MeV, 0.5–1.0 MeV, 1.0–2.0 MeV, and E > 2.0 MeV for electrons, 5–10 MeV, 10–20 MeV, 20–40 MeV, and 40–80 MeV for protons.     

Solar activity dependence of total electron content derived from GPS observations over Mbarara

Vertical total electron content (VTEC) observed at Mbarara (geographic co-ordinates: 0.60°S, 30.74°E; geomagnetic coordinates: 10.22°S, 102.36°E), Uganda, for the period 2001–2009 have been used to study the diurnal, seasonal and solar activity variations. The daily values of the 10.7 cm radio flux (F_10.7) and sunspot number (R) were used to represent Solar Extreme Ultraviolet Variability (EUV). VTEC is generally higher during high solar activity period for all the seasons and increases from 0600 h LT and reaches its maximum value within 1400 h–1500 h LT. All analysed linear and quadratic fits demonstrate positive VTEC-F_10.7 and positive VTEC-R correlation, with all fits at 0000 h and 1400 h LT being significant with a confidence level of 95% when both linear and quadratic models are used. All the fits at 0600 h LT are insignificant with a confidence level of 95%. Generally, over Mbarara, quadratic fit shows that VTEC saturates during all seasons for F_10.7 more than 200 units and R more than 150 units. The result of this study can be used to improve the International Reference Ionosphere (IRI) prediction of TEC around the equatorial region of the African sector.     

Propagation of normal and faster CMEs in the interplanetary medium

We have analyzed 101 Coronal Mass Ejection (CME) events and their associated interplanetary CMEs (ICMEs) and interplanetary (IP) shocks observed during the period 1997–2005 from the list given by Mujiber Rahman et al. (2012). The aim of the present work is to correlate the interplanetary parameters such as, the speeds of IP shocks and ICMEs, CME transit time and their relation with CME parameters near the Sun. Mainly, a group of 10 faster CME events (V_INT > 2200 km/s) are compared with a list of 91 normal events of Manoharan et al. (2004). From the distribution diagrams of CME, ICME and IP shock speeds, we note that a large number of events tends to narrow towards the ambient (i.e., background) solar wind speed (∼500 km/s) in agreement with the literature. Also, we found that the IP shock speed and the average ICME speed measured at 1 AU are well correlated. In addition, the IP shock speed is found to be slightly higher than the ICME speed. While the normal events show CME travel time in the range of ∼40–80 h with a mean value of 65 h, the faster events have lower transit time with a mean value of 40 h. The effect of solar wind drag is studied using the correlation of CME acceleration with interplanetary (IP) acceleration and with other parameters of ICMEs. While the mean acceleration values of normal and faster CMEs in the LASCO FOV are 1 m/s², 18 m/s², they are −1.5 m/s² and −14 m/s² in the interplanetary medium, respectively. The relation between CME speed and IP acceleration for normal and faster events are found to agree with that of  and  except slight deviations for the faster events. It is also seen that the faster events with less travel time face higher negative acceleration (>−10 m/s²) in the interplanetary medium up to 1 AU.     

Modeling PSBL high speed ion beams observed by Cluster and Double Star

On October 8, 2004, the Cluster and Double Star spacecraft crossed the near-Earth (12–19 R_E) magnetotail neutral sheet during the recovery phase of a small, isolated substorm. Although they were separated in distance by ∼7 R_E and in time by ∼30 min, both Cluster and Double Star observed steady, but highly structured Earthward moving >1000 km/s high speed H⁺ beams in the PSBL. This paper utilizes a global magnetohydrodynamic (MHD) simulation driven by Wind spacecraft solar wind input to model the large-scale structure of the PSBL and large-scale kinetic (LSK) particle tracing calculations to investigate the similarities and differences in the properties of the observed beams. This study finds that the large-scale shape of the PSBL is determined by the MHD configuration. On smaller scales, the LSK calculations, in good qualitative agreement with both Cluster and Double Star observations, demonstrated that the PSBL is highly structured in both time and space, on time intervals of less than 2 min, and spatial distances of the order of 0.2–0.5 R_E. This picture of the PSBL is different from the ordered and structured region previously reported in observations.