Gamma-ray astronomy and cosmic ray origin problem

We analyse the results of the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) in order to describe their relevant properties: the remnant dynamics and the characteristics of nonthermal emission produced by CRs. It is shown that the theory fits the existing data in a satisfactory way and that the magnetic field in SNRs is significantly amplified due to efficient acceleration of the nuclear CR component. From the fact that magnetic field amplification occurs in all the young SNRs for which relevant data exist, and given the strong theoretical connection between magnetic field amplification and efficiently accelerating the nuclear CRs, we tentatively conclude that the Galactic SNRs are the source population of the Galactic CRs. Due to high interior magnetic fields in young SNRs the π⁰-decay γ-rays generated by the nuclear CR component as a rule dominate over γ-rays generated by the electron CR component, and the calculated γ-ray flux fits existing data.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

Radiation from relativistic shocks in turbulent magnetic fields

Using our new 3-D relativistic particle-in-cell (PIC) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron–positron jet propagating in an unmagnetized ambient electron–positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic shock compression. In the jet (reverse) shock behind the bow (forward) shock the strongest electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to calculate radiation from first principles that goes beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations. Initially we calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We then used the technique to calculate emission from electrons in a small simulation system. From these simulations we obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger simulation system may generate a jitter/synchrotron spectrum.     

New insights into SNR evolution revealed by the discovery of recombining plasmas

We report the discovery of recombining plasmas in three supernova remnants (SNRs) with the Suzaku X-ray astronomy satellite. During SNR’s evolution, the expanding supernova ejecta and the ambient matter are compressed and heated by the reverse and forward shocks to form an X-ray emitting hot plasma. Since ionization proceeds slowly compared to shock heating, most young or middle-aged SNRs have ionizing (underionized) plasmas. Owing to high sensitivity of Suzaku, however, we have detected radiative recombination continua (RRCs) from the SNRs IC 443, W49B, and G359.1–0.5. The presence of the strong RRC is the definitive evidence that the plasma is recombining (overionized). As a possible origin of the overionization, an interaction between the ejecta and dense circumstellar matter is proposed; the highly ionized gas was made at the initial phase of the SNR evolution in dense regions, and subsequent rapid adiabatic expansion caused sudden cooling of the electrons. The analysis on the full X-ray band spectrum of IC 443, which is newly presented in this paper, provides a consistent picture with this scenario. We also comment on the implications from the fact that all the SNRs having recombining plasmas are correlated with the mixed-morphology class.     

Magnetic reconnection in a high-temperature plasma of solar flares

In the frame of a simple self-consistent model for high-temperature turbulent current sheet (HTCS) /1/, three effects are considered. (i) Gradient instabilities create anamalous plasma diffusion across magnetic field and increase the power of energy release in HTCS. (ii) Penetration of a small transverse component of magnetic field into HTCS also can significantly increase an energy output of HTCS. (iii) There appears electric current circulating around a current sheet at a compression of longitudinal magnetic field. This current induces a Joule heat; however, a total flux of the longitudinal field remains constant.     

Magnetogasdynamic shock wave generated by a moving piston in a rotational axisymmetric isothermal flow of perfect gas with variable density

The propagation of a strong cylindrical shock wave in an ideal gas with azimuthal magnetic field, and with or without axisymmetric rotational effects, is investigated. The shock wave is driven out by a piston moving with time according to power law. The ambient medium is assumed to have radial, axial and azimuthal component of fluid velocities. The fluid velocities, the initial density and the initial magnetic field of the ambient medium are assumed to be varying and obey power laws. Solutions are obtained, when the flow between the shock and the piston is isothermal. The gas is assumed to have infinite electrical conductivity and the angular velocity of the ambient medium is assumed to be decreasing as the distance from the axis increases. It is expected that such an angular velocity may occur in the atmospheres of rotating planets and stars. The shock wave moves with variable velocity and the total energy of the wave is non-constant. The effects of variation of the initial density and the Alfven-Mach number on the flow-field are obtained. A comparison is also made between rotating and non-rotating cases.     

Solar wind interaction with lunar crustal magnetic anomalies

Using Lunar Prospector data, we review the magnetic field and electron signatures of solar wind interaction with lunar crustal magnetic sources. Magnetic field amplifications, too large to represent direct measurements of crustal fields, appear in the solar wind over strong crustal sources, with the chance of observing these amplifications depending on upstream solar wind parameters. We often observe increases in low-energy (?100 eV) electron energy fluxes simultaneously with large magnetic field amplifications, consistent with an increase in plasma density across a shock surface. We also often observe low frequency wave activity in the magnetic field data (both broadband turbulence and monochromatic waves), often associated with electron energization, sometimes up to keV energies. Electron energization appears to be correlated more closely with wave activity than with magnetic amplifications. Detailed studies of the interaction region will be necessary in order to understand the physics of the Moon–solar wind interaction. At present, the Moon represents the only natural laboratory available to us to study solar wind interaction with small-scale crustal magnetic fields, though simulation results and theoretical work can also help us understand the physical processes at work.     

强磁场结构在行星际空间膨胀的数值模拟

本文研究了强磁场结构在18—240R_s的行星际空间的膨胀效应。模拟结果表明, 强磁场结构的膨胀能够形成典型的磁云剖面结构;膨胀速度明显依赖于磁场强度的大小;强磁场结构运动的速度直接影响到它的动力学特性。      

Geomagnetic activity associated with magnetic clouds,magnetic cloud-like structures and interplanetary shocks for the period 1995–2003

Using nine years (1995–2003) of solar wind plasma and magnetic field data, solar sunspot number, and geomagnetic activity data, we investigated the geomagnetic activity associated with magnetic clouds (MCs), magnetic cloud-like structures (MCLs), and interplanetary shock waves. Eighty-two MCs and one hundred and twenty-two MCLs were identified by using solar wind and magnetic field data from the WIND mission, and two hundred and sixty-one interplanetary shocks were identified over the period of 1995–2003 in the vicinity of Earth. It is found that MCs are typically more geoeffective than MCLs or interplanetary shocks. The occurrence frequency of MCs is not well correlated with sunspot number. By contrast, both occurrence frequency of MCLs and sudden storm commencements (SSCs) are well correlated with sunspot number.     

Interplanetary coronal mass ejections during the descending cycle 23: Sheath and ejecta properties comparison

We have used Omniweb data in order to identify the sheath and the ejecta boundaries of 67 shock-driving interplanetary coronal mass ejections during the time period 2003–2006. We examine and compare their statistical properties (speed, magnetic field strength, proton density and temperature, proton plasma beta), with those of the typical solar wind. We also calculate their passage time and radial width. We study the correlation between the ejecta and sheath characteristics.     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

Evolution of solar wind structures from 0.72 to 1 AU

Understanding the evolution of solar wind structures in the inner heliosphere as they approach the Earth is important to space weather prediction. From the in situ solar wind plasma and magnetic field measurements of Pioneer Venus Orbiter (PVO) at 0.72 AU (1979–1988), and of Wind/Advanced Composition Explorer (ACE) missions at 1 AU (1995–2004), we identify and characterize two major solar wind structures, stream interaction regions (SIRs) and interplanetary coronal mass ejections (ICMEs). The average percentage of SIRs occurring with shocks increases significantly from 3% to 24% as they evolve from 0.72 to 1 AU. The average occurrence rate, radial extent, and bulk velocity variation of SIRs do not change from 0.72 to 1 AU, while peak pressure and magnetic field strength both decrease with the radial evolution of SIRs. Within the 0.28 AU distance from the orbit of Venus to that of Earth, the average fraction of ICMEs with shocks increases from 49% to 66%, and the typical radial extent of ICMEs expands by about a fraction of 1.4, with peak pressure and magnetic field strength decreasing significantly. The mean occurrence rate and expansion velocity of ICMEs do not change from 0.72 to 1 AU.     

Dynamics of the convection in the inner magnetosphere by observations of the diffuse aurora,SAR arc and ionospheric drift

The first results of the comparison of subauroral luminosity dynamics in 557,7 and 630,0 nm emission with simultaneous measurements of the ionospheric drift in the F2 region with a digisonde DPS-4 at the Yakutsk meridian (CGMC: 55–60N, 200°E) at Kp = 2–6 are presented. It is shown from the analysis of individual events that during the magnetospheric convection intensification after the turn of the IMF B_z – component to the south the equatorward extension of diffuse aurora takes place. At the same time the westward ionospheric drift velocity increases both in the diffuse aurora region and much equatorward of it due to the occurrence of the northward polarization electric field. We suppose that the generation of polarization field can be associated with the development of the region 2 FAC during the intensification of magnetospheric convection. The comparison of ground-based observations with measurements of the plasma drift aboard the DMSP-F15 satellite has been carried out.     

The evolution of mirror type magnetic fluctuations in the magnetosheath based on multipoint observations

Two orbits were selected in January–February 2006 when the separation between the Cluster spacecraft was large and mirror type magnetic field fluctuations were observed by all spacecraft in different regions of the terrestrial magnetosheath. Minimum variance analysis was applied to find the mirror type fluctuations, and the amplitude of the fluctuations was determined individually. Mirror mode structures are moving along the streamlines frozen in the plasma. A model was developed for the calculation of plasma flowtime from the bow shock to the observation point. The growth rate of the field strength perturbations was estimated by comparing the amplitudes of fluctuations observed simultaneously at distant locations (∼10,000 km) based on the assumption that δB ∼ exp(γt). The obtained growth rate values were about an order of magnitude smaller than those provided by linear models and they decreased in the inner regions of the magnetosheath, indicating some saturation in the growth of the waves when proceeding towards the magnetopause. The results of these two case studies suggest that mirror type fluctuations originate from the compression region downstream of the quasi-perpendicular bow shock, and the growth of the fluctuations cannot be described by linear approximations.     

Observational appearances of isolated stellar-mass black hole accretion – Theory and observations

General properties of accretion onto isolated stellar-mass black holes in the Galaxy are discussed. An analysis of plasma internal energy growth during the infall is performed. Adiabatic heating of collisionless accretion flow due to magnetic adiabatic invariant conservation is 25% more efficient than in the standard non-magnetized gas case. It is shown that magnetic field line reconnections in discrete current sheets lead to significant nonthermal electron component formation, which leads to a formation of a hard (UV, X-ray, up to gamma), highly variable spectral component in addition to the standard synchrotron optical component first derived by Shvartsman generated by thermal electrons in the magnetic field of the accretion flow. Properties of accretion flow emission variability are discussed. Observation results of two single black hole candidates – gravitational lens MACHO-1999-BLG-22 and radio-loud X-ray source with featureless optical spectrum J1942+10 – in optical band with high temporal resolution are presented and interpreted in the framework of the proposed model.     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

Cosmic ray acceleration by supernova shocks

We analyse the results of recent measurements of nonthermal emission from individual supernova remnants (SNRs) and their correspondence to the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is shown that the theory fits these data in a satisfactory way and provides the strong evidences for the efficient CR production in SNRs accompanied by significant magnetic field amplification. Magnetic field amplification leads to considerable increase of CR maximum energy so that the spectrum of CRs accelerated in SNRs is consistent with the requirements for the formation of Galactic CR spectrum up to the energy ∼10¹⁷ eV.     

无碰撞激波上游低频波的激发及其对激波耗散机制的影响

由于很大一部分来自激波上游的粒子被激波面所反射,因此在准平行无碰撞激波的上游存在着等离子体束流．通过一维混合模拟方法,计算了束流密度较小(nb／n0=0．02)和较大(nb／n0=0．2)两种不同情形下的等离子体束流不稳定性．结果表明,在束流密度较小时,束流激发的主要是平行于背景磁场方向传播的右旋波,此波动只能对束流粒子产生影响,包括减速和加热．在束流密度较大时,束流可同时激发平行和反平行于背景磁场方向传播的右旋波,除能对束流粒子产生影响外,还可通过非共振作用加速和加热背景粒子．文中对准平行无碰撞激波耗散机制的影响也进行了讨论。     

Kink-like mode of a double gradient instability in a compressible plasma current sheet

A linear MHD instability of the electric current sheet, characterized by a small normal magnetic field component, varying along the sheet, is investigated. The tangential magnetic field component is modeled by a hyperbolic function, describing Harris-like variations of the field across the sheet. For this problem, which is formulated in a 3D domain, the conventional compressible ideal MHD equations are applied. By assuming Fourier harmonics along the electric current, the linearized 3D equations are reduced to 2D ones. A finite difference numerical scheme is applied to examine the time evolution of small initial perturbations of the plasma parameters. This work is an extended numerical study of the so called “double gradient instability”, – a possible candidate for the explanation of flapping oscillations in the magnetotail current sheet, which has been analyzed previously in the framework of a simplified analytical approach for an incompressible plasma. The dispersion curve is obtained for the kink-like mode of the instability. It is shown that this curve demonstrates a quantitative agreement with the previous analytical result. The development of the instability is investigated also for various enhanced values of the normal magnetic field component. It is found that the characteristic values of the growth rate of the instability shows a linear dependence on the square root of the parameter, which scales uniformly the normal component of the magnetic field in the current sheet.