Particle acceleration and transport at an oblique CME-driven shock

In gradual solar energetic particle (SEP) events, protons and heavy ions are often accelerated to >100 MeV/nucleon at a CME-driven shock. In this work, we study particle acceleration at an oblique shock by extending our earlier particle acceleration and transport in heliosphere (PATH) code to include shocks with arbitrary θ_BN, where θ_BN is the angle between the upstream magnetic field and the shock normal. Instantaneous particle spectra at the shock front are obtained by solving the transport equation using the total diffusion coefficient κ, which is a function of the parallel diffusion coefficient κ_∥ and the perpendicular diffusion coefficient κ_⊥. In computing κ_∥ and κ_⊥, we use analytic expressions derived previously. The particle maximum energy at the shock front as a function of time, the time intensity profiles and particle spectra at 1 AU for five θ_BN’s are calculated for an example shock.     

Modeling the effects of pitch-angle scattering processes on the transport of solar energetic particles along the interplanetary magnetic field

We present a Monte-Carlo technique to study the time-dependent transport of energetic particles in the interplanetary medium. We use the guiding center approximation between discrete finite pitch-angle scatterings to quantify the competing effects of focusing and pitch-angle scattering on energetic particles propagating along a Parker spiral magnetic field. We consider that the pitch-angle scattering process is produced by small-scale magnetic field irregularities frozen in the expanding solar wind. We also include the effects of both solar wind convection and adiabatic deceleration. We use a joint probability distribution P(h, μ′) = p(h; μ′)q(μ′; μ) to describe both the distance traveled by the particle between two scattering processes and the change in the particle pitch-angle after a scattering process. Here, p(h; μ′) is the conditional probability that the particle travels a distance h along the field line before the next scattering if it had a pitch-angle cosine μ′ after the previous scattering, and q(μ′; μ) is the conditional probability for the pitch-angle cosine μ′ if the pitch-angle cosine was μ before the scattering. We consider several functional forms to describe the processes of pitch-angle scattering, such as an isotropic scattering without any memory of the initial particle’s pitch-angle or processes in which the scattering result depends upon the initial particle’s pitch-angle. The results of our simulations are pitch-angle distributions and time-intensity profiles that can be directly compared to spacecraft observations. Comparison of our simulations with near-relativistic (45–290 keV) electron events observed by the Electron, Proton and Alpha Monitor on board the Advanced Composition Explorer allows us to estimate both the time dependence of the injection of near-relativistic electrons into the interplanetary medium and the conditions for electron propagation along the interplanetary magnetic field.     

Hard X-ray footpoint motions in solar flares: Comparing magnetic reconnection models with observations

Hard X-ray observations from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of the October 29, 2003 GOES X10 two-ribbon flare are used together with magnetic field observations from the Michelson Doppler Imager (MDI) onboard SoHO to compare footpoint motions with predictions from magnetic reconnection models. The temporal variations of the velocity v of the hard X-ray footpoint motions and the photospheric magnetic field strength B in footpoints are investigated. The underlying photospheric magnetic field strength is generally higher (B ∼ 700–1200 G) in the slower moving (v ∼ 20–50 km s⁻¹) western footpoint than in the faster (v ∼ 20–100 km s⁻¹) moving eastern source (∼100–600 G). Furthermore, a rough temporal correlation between the HXR flux and the product vB² is observed.     

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

Scaling and universality of the thermal conductivity of liquid 4He near the superfluid transition and in restricted geometries

We present measurements of the thermal conductivity λ(t, P, L) = l/ρ(t, P, L) near the superfluid transition of ⁴He at saturated vapor pressure and confined in cylindrical geometries with radii L = 0.5 and 1.0 μm (t ≡ T/T_λ(P) − 1). For L = 1.0 μm measurements at six pressures P are presented. At and above T_λ the data are consistent with a universal scaling function F(X) = (L/ξ_o)^x/ν(ρ/ρ₀), X = (L/ξ_o)^1/νt valid for all P (ρ₀ and x are the pressure-dependent amplitude and effective exponent of the bulk resistivity ρ(t, P, ∞) = ρ₀t^x and ξ = ξ₀t^−ν is the correlation length). Indications of breakdown of scaling and universality are observed below T_λ.     

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.     

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.     

Centre-to-limb variation of photospheric facular radiance and image resolution

We study the effect of the angular resolution on the determination of the angular properties of the facular radiance. We analyze photospheric intensity in the continuum, around the Ni 676.8 nm line, and longitudinal magnetic field along the line of sight, measured by the MDI instrument aboard SOHO with two spatial resolutions, 4″ and 1.2″ (2″ and 0.6″ pixels, respectively). The effect of the limited photometric sensitivity of the instrument and the limited information on the angular structure of the magnetic field tubes are considered. Our study of the high-resolution data shows that intensity contrast of magnetic features between 80 and 600 Gauss increases from centre to limb up to a maximum that occurs at higher heliocentric angles (θ) when obtained with higher resolution data than for lower resolution data. There is a suggestion that at heliocentric angles below about 75° there is only a monotonic increase in the contrast as one goes from cos (θ) = 1 to cos (θ) = 0.2.     

Low-frequency electric field and density fluctuation measurements on Solar Orbiter

Solar Orbiter will orbit the Sun down to a distance of 0.22 AU allowing detailed in situ studies of important but unexplored regions of the solar wind in combination with coordinated remote sensing of the Sun. In-situ measurements require high quality measurements of particle distributions and electric and magnetic fields. We show that such important scientific topics as the identification of coronal heating remnants, solar wind turbulence, magnetic reconnection and shock formation within coronal mass ejections all require electric field and plasma density measurements in the frequency range from DC up to about 100 Hz. We discuss how such measurements can be achieved using the double-probe technique. We sketch a few possible antenna design solutions.     

The spatial distribution of galactic and anomalous cosmic rays in the heliosphere at solar minimum

The spatial distributions of galactic and anomalous cosmic rays in the heliosphere at the solar minima of Cycles 20/22 (qA > 0) and of Cycle 21 (qA < 0) are studied, using data from IMP 8, Voyagers 1/2 and Pioneer 10. It is found that the radial dependences of intensities J can be approximated by a power of radial distance r as J ∝ r^α with a different value of a constant in the inner and outer heliosphere with a transition at a radial distance of 10–15 AU. To study the physical meaning of these radial intensity profiles we examined the rigidity dependences of the intensity gradients by determining the particle mean free paths, using a simple one-dimensional modulation model. The particle mean free path λ was assumed to be a separable function of distance of the form r^γ and rigidity R of R^δ over the range of 0.5–3.0 GV in the inner and outer heliosphere. It was shown that λ of rigidity dependence of R^1.6 determined for Cycle 20/22 (qA > 0) with anomalous He is about 4–5 times larger than that of Cycle 21 (qA < 0) with R^0.9 at around 1 GV in the outer heliosphere, and that the radial dependences are r^1.4 and r^1.1, respectively, for Cycles 20/22 and for Cycle 21.     

A Gemini observation of the anomalous X-ray pulsar 1RXS J170849-400910

The anomalous X-ray pulsars (AXPs) represent a growing class of neutron stars discovered at X-ray energies. While the nature of their multi-wavelength emission mechanism is still under debate, evidence has been recently accumulating in favor of their magnetar nature. Their study in the optical and infrared (IR) wavelengths has recently opened a new window to constrain the proposed models. We here present a brief overview of AXPs and our Gemini-South observation of 1RXS J170849-400910, which is a relatively bright AXP discovered with ROSAT and later found to be an 11 s X-ray pulsar by ASCA. The observation was taken with the near-IR imager Flamingos in J (1.25 μm), H (1.65 μm), and K_s (2.15 μm). We confirm the recent detection by (ApJ, 589, L93–L96) of a source coincident with the CHANDRA source (candidate ‘A’). Our derived magnitudes of J = 20.6 (0.2), H = 18.6 (0.2), and K_s = 17.1 (0.2) are consistent with those derived by (ApJ, 589, L93–L96), and indicate that if this source is indeed the IR counterpart to 1RXS J170849-400910, then there is no evidence of variability from this AXP. However, given the lack of IR variability and the relatively high IR to X-ray flux of this source when compared to the other AXPs, we conclude that this source is unlikely the counterpart of the AXP, and that the other source (candidate ‘B’) within the CHANDRA error circle should not be ruled out as the counterpart. Further monitoring of these sources and a deep observation of this complex field are needed to confirm the nature of these sources and their association with the AXP.     

A search for radio emission from the young 16-ms X-ray pulsar PSR J0537−6910

PSR J0537−6910 is a young, energetic, rotation-powered X-ray pulsar with a spin period of 16 ms located in the Large Magellanic Cloud. We have searched for previously undetected radio pulsations (both giant and standard) from this pulsar in a 12-h observation taken at 1400 MHz with the Parkes 64-m radio telescope. The very large value of the magnetic field at the light cylinder radius suggests that this pulsar might be emitting giant radio pulses like those seen in other pulsars with similar field strengths. No radio emission of either kind was detected from the pulsar, and we have established an upper limit of ∼25 mJy kpc² for the average 1400-MHz radio luminosity of PSR J0537−6910. The 5σ single-pulse detection threshold was ∼750 mJy for a single 80-μs sample. These limits are likely to be the best obtainable until searches with greatly improved sensitivity can be made with next-generation radio instruments.     

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⁻¹.     

XMM-Newton observations of the Mouse,SLX 1744–299 and SLX 1744–300

We observed the radio and X-ray source G359.23–0.82, also known as “the Mouse”, with XMM-Newton. The X-ray image of this object shows a point-like source at the Mouse’s “head”, accompanied by a “tail” that extends for about 40″ westward. The morphology is consistent with that observed recently with Chandra [Gaensler, B.M., van der Swaluw, E., Camilo, F., et al. The Mouse that soared: high resolution X-ray imaging of the pulsar-powered bow shock G359.23–0.82, ApJ 616, 383–402, 2004]. The spectrum of the head can be described by a power-law model with a photon index Γ ≃ 1.9. These results confirm that the Mouse is a bow-shock pulsar wind nebula (PWN) powered by PSR J1747–2958. We found that the hydrogen column density toward the Mouse, N_H = (2.60 ± 0.09) × 10²² cm⁻², is 20%–40% lower than those toward two serendipitously detected X-ray bursters, SLX 1744–299 and SLX 1744–300. At a plausible distance of 5 kpc, the X-ray luminosity of the Mouse, L(0.5–10 keV) = 3.7 × 10³⁴ erg s⁻¹, is 1.5% of the pulsar’s spin-down luminosity. We detected a Type I X-ray burst from SLX 1744–300 and found a possible decrease of N_H and persistent luminosity for this source, in comparison with those observed with ROSAT in 1992.     

Supernova remnant 1987A: The latest report from the Chandra X-ray Observatory

We continue monitoring supernova remnant (SNR) 1987A with the Chandra X-ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the southwest are now evident in addition to the bright eastern ring. The overall X-ray spectrum, since 2002 December, can be described by a planar shock with an electron temperature of ∼2.1 keV. The soft X-ray flux is now 8 × 10⁻¹³ ergs cm⁻² s⁻¹, which is about five times higher than four years ago. This flux increase rate is consistent with our prediction based on an exponential density distribution along the radius of the SNR between the HII region and the inner ring. We still have no direct evidence of a central point source, and place an upper limit of L_X = 1.3 × 10³⁴ ergs s⁻¹ on the 3–10 keV band X-ray luminosity.     

The radio signature of a kink unstable loop

Based upon multi-wavelength observations outlined by Huang et al. [Solar Phys. 213 (2003) 341], especially the dynamic spectrum at 4.5–7.5 GHz, we study the physical nature of the radio fine structures (FS) during the flare on August 25, 1999 in AR 8674 (S28E21). The main results are: (1) the helical loop of the event related to the FS is unstable for m = 1 kink mode; (2) the time interval between the beginning of reconnection and relaxation of the unstable loop, inferred from the observation, is quantitatively consistent with the results of the numerical simulations on kink unstable loop; (3) the magnetic field strength estimated from the fast kink standing wave is basically of the same order as that estimated from the photospheric magnetic field, which provides strong support to our analysis.     

On the quality of mutations in mammalian cells induced by high LET radiations

The deleterious effects of accelerated heavy ions as component of the space radiation environment on living cells are of increasing importance for long duration human space flight activities. The most important aspect of such densely ionizing particle radiation is attributed to the type and quality of biological damage induced by them. This issue is addressed by investigating cell inactivation and mutation induction at the Hprt locus (coding for hypoxanthine-guanine-phosphoribosyl-transferase) of cultured V79 Chinese hamster cells exposed to densely ionizing radiation (accelerated heavy ions with different LETs from oxygen to gold, specific energies ranging from 1.9 to 69.7 MeV/u, corresponding LET values range from 62 to 13,223 keV/μm) and to sparsely ionizing radiation (200 kV X-rays). 30 spontaneous, 40 X-ray induced and 196 heavy ion induced 6-thioguanine resistant Hprt mutant colonies were characterized by Southern technique using the restriction enzymes EcoRI, PstI and BglII and a full length Hprt cDNA probe isolated from the plasmid pHPT12. Restriction patterns of the spontaneous Hprt mutants were indistinguishable from the wild type pattern, as these mutants probably contain only small deletions or even point mutations in the Hprt locus. In contrast, the overall spectrum of heavy ion induced mutations revealed a majority of partial or total deletions of the Hprt gene. With constant particle fluence (3 × 10⁶ particles/cm²) the quality of heavy ion induced mutations in the Hprt locus depends on physical parameters of the beam (atomic number, specific energy, LET). This finding suggests a relationship between the type of DNA damage and track structure. The fraction of mutants with severe deletions in the Hprt locus after exposure to oxygen ions increases from 65% at 60 keV/μm up to a maximum (100%) at 300 keV/μm and declines with higher LET values to 75% at 750 keV/μm. With heavier ions (Ca- and Au-ions) and even higher LET-values this mutant fraction decreases to 58% at 13,200 keV/μm. Heavy ion induced DNA break points in the Hprt locus are not randomly distributed.     

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.     

Active galactic nuclei imaging programs of the RadioAstron mission

Imaging relativistic jets in active galactic nuclei (AGN) at angular resolution significantly surpassing that of the ground-based VLBI at centimetre wavelengths is one of the key science objectives of the RadioAstron space-VLBI mission. There are three RadioAstron imaging key science programs that target both nearby radio galaxies and blazars, with one of the programs specifically focusing on polarimetry of the jets. The first images from these programs reach angular resolution of a few tens of microarcseconds and reveal unprecedented details about the jet collimation profile, magnetic field configuration, and Kelvin-Helmholtz instabilities along the flow in some of the most studied AGN (3C 84, BL Lac, 3C 273, S5 0836 + 710). Here we give an overview of the goals and strategy of these three ongoing programs, highlight their early results, and discuss the challenges of space-VLBI imaging.     

Spectral analysis of a two-ribbon microflare

Using high-resolution Hα, CaII 8542 Å and FeI 6302.5 Å Stokes spectral data obtained simultaneously with THEMIS in 2002 September, we have analyzed the spectra and the characteristics of a two-ribbon microflare (MF). The hard X-ray emission provides evidence of non-thermal particle acceleration in the microflare. The two-ribbons are located on either sides of the magnetic polarity inversion line. The non-thermal characteristics mainly appeared at the outer edges of the flare ribbons. It indicates that the instantaneous magnetic reconnection and the particle acceleration mainly took place at the outer edges of the flare ribbons. Using the Hα and CaII 8542 Å line profiles and the non-LTE calculation, we obtain the semi-empirical atmospheric model for the bright kernel of the MF. The result indicates that the temperature enhancement in the chromosphere is about 2000–2500 K.