Modulation of cosmic rays: Perpendicular diffusion and drifts in a heliosphere with a solar wind termination shock

In this study the roles of polar perpendicular diffusion and drifts are illustrated in a model containing a heliosheath and diffusive shock acceleration as applied to the solar wind termination shock. Of particular interest is the relation of polar perpendicular diffusion to particle drifts and how the effectiveness of the termination shock acceleration of galactic and anomalous protons is influenced by this relation. We found that drifts have a more prominent effect than the polar enhancement of perpendicular diffusion so that its omission from termination shock models would produce unrealistically large shock acceleration and consequently also larger modulation effects throughout the heliosphere. The computed spectra at a heliolatitude of 35° are almost similar for the two polarity epochs indicating that the two Voyager spacecraft might not observe differences between the two cycles in future.     

Microscopic approach to analyze solar-sail space-environment effects

Near-sun space-environment effects on metallic thin films solar sails as well as hollow-body sails with inflation fill gas are considered. Analysis of the interaction of the solar radiation with the solar-sail materials is presented. This analysis evaluates worst-case solar radiation effects during solar-radiation-pressure acceleration. The dependence of the thickness of solar sail on temperature and on wavelength of the electromagnetic spectrum of solar radiation is investigated. Physical processes of the interactions of photons, electrons, protons and α-particles with sail material atoms and nuclei, and inflation fill gas molecules are analyzed. Calculations utilized conservative assumptions with the highest values for the available cross sections for interactions of solar photons, electrons and protons with atoms, nuclei and hydrogen molecules. It is shown that for high-energy photons, electrons and protons the beryllium sail is mostly transparent. Sail material will be partially ionized by solar UV and low-energy solar electrons. For a hollow-body photon sail effects including hydrogen diffusion through the solar-sail walls, and electrostatic pressure is considered. Electrostatic pressure caused by the electrically charged sail’s electric field may require mitigation since sail material tensile strength decreases with elevated temperature. It also can substitute inflation-gas pressure loss due to gas diffusion and perforation by micrometeoroids impact to keep the sail inflated.     

Energetic positive ions in the cometary “foreshock” region

Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading); hence, they can energize the incoming ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind. It is suggested that Comet Halley may be the most efficient “cosmic ray shock” in the solar system.     

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.     

2008年4月24日行星际激波事件相关联的超热电子90°投掷角的增强

利用测试粒子数值模拟的方法研究了与STEREO-A卫星观测到的2008年4月24日行星际激波事件相关联的超热电子90°投掷角的增强.根据激波到达前给定时刻超热电子的观测分布,拟合得到不同投掷角的初始分布函数;在给定的激波参数下,采用时间向后的方法计算特定能道上激波下游超热电子的投掷角分布.由于超热电子具有较高的共振频率,模拟采用的磁场湍流谱包含了低能电子发生共振的耗散区.对以215.76,151.67,106.63,eV为中心的三个能道进行了模拟.结果表明,不同能道上超热电子在激波下游的投掷角分布均在90°投掷角附近出现峰值,呈现出明显的90°投掷角增强,这与观测结果符合得很好.可以认为在激波对电子的加速过程中,电子与湍流耗散区的共振对90°投掷角的增强具有重要作用.      

Energetic particles in the comet Halley environment

An overview is presented of electrons, protons and heavier ions (E > 20 keV) recorded by the energetic particle detector EPONA in the Comet Halley environment, 12–15 March, 1986. Pick-up ions were detected at distances of up to at least 7.5 × 10⁶ km from the nucleus. Estimates of the energies that typical cometary ions may be expected to acquire from the solar wind pertaining at Encounter show that the pick-up process is insufficient to account for the energies of the particles detected. An additional mechanism must thus be postulated to account for the observed particle signatures. Preliminary correlations with magnetic and plasma wave data from other instruments suggest that the presence of MHD turbulence at several million kilometers upstream of the bowshock may have contributed to the acceleration of the first pick-up ions observed. The bowshock boundary (inbound) does not appear to have constituted a location where particle acceleration to high energies took place. Downstream of the shock boundary, hardening of the energy spectrum and the development of less anisotropic particle streaming was observed to occur when the spacecraft was in a turbulent environment 1 × 10⁶ km from the nucleus. The waxing influence of mass loading as a mechanism for reducing energetic particle fluxes as well as the depletion of energetic ions due to their escape along open field lines and to charge exchange collision with neutrals in a progressively more stagnant solar wind, may be inferred in a regime (seen on the magnetometer data to be largely non-turbulent) traversed by the spacecraft from 5 × 10⁵ km from the nucleus to within the magnetic pile-up region. A major burst of ions and electrons (not yet established to be of cometary origin) occurred when the spacecraft was close to the Contact Surface. A population of high energy electrons (from 180 keV to at least 300 keV) was detected for about one hour before Closest Approach and for several hours thereafter. Also an energetic beam of electrons was identified exiting from a location at about 1 × 10⁶ km from the nucleus (outbound). Finally, differences between inbound and outbound particle signatures are described.     

Interplanetary acceleration of low energy protons in association with shock waves: The 11 February 1979 event

The problem of interplanetary acceleration of low energy protons in association with shock waves is examined in the context of the specific event observed on 11 February 1979 on board the ISEE-3 spacecraft. This event has been selected for special study as it apparently was not associated with a solar flare event. The low energy proton telescope system on ISEE-3 measures the proton distribution function with good spectral, directional and temporal resolution from E_p = 35 keV. The evolution of the anisotropies and of the energy spectrum during the event are consistent with particle acceleration taking place in the vicinity of the shock wave.     

Energy content of accelerated electrons and ions in intense solar flares

The energy content of nonthermal particles in solar flares is shared between accelerated electrons and ions. It isimportant for understanding the particle acceleration mechanism in solar flares. Yohkoh observed a few intense flares which produced both strong gamma-ray lines and electron bremsstrahlung continuum. We analyze energy spectra of X-class solar flares on October 27, 1991(X6.1), November 6, 1997 (X9.4), July 14, 2000 (X5.7) and November 24, 2000 (X2.3). The accelerated electron and proton spectra are derived from a spectral analysis of their high-energy photon emission and the energy contents in >1 MeV electrons and >10 MeV protons are estimated to be 6×l0²⁸ – 4×10³⁰ and 2×10²⁸ – 5×10²⁹ erg, respectively. We study the flare to flare variation in the energy content of >1 MeV electrons and >10 MeV protons for the four Yohkoh gamma-ray flares. Ratios of >1 MeV electron energy content to >10 MeV proton energy content are roughly within an order of magnitude.     

Simulation of ASTROD I test mass charging due to solar energetic particles and interplanetary electrons

As ASTROD I travels through space, its test mass will accrue charge due to exposure of the spacecraft to high-energy particles. This test mass charge will result in Coulomb forces between the test mass and the surrounding electrodes. In earlier work, we have used the GEANT 4 toolkit to simulate charging of the ASTROD test mass due to cosmic-ray protons of energies between 0.1 and 1000 GeV at solar maximum and at solar minimum. Here we use GEANT 4 to simulate the charging process due to solar energetic particle events and interplanetary electrons. We then estimate the test mass acceleration noise due to these fluxes. The predicted charging rates range from 2247 e⁺/s to 47,055 e⁺/s, at peak intensity, for the four largest SEP events in September and October 1989. Although the noise due to charging exceeds the ASTROD I budget for the two larger events, it can be suppressed through continuous discharging. The acceleration noise during the two small events is well below the design target. The charging rate of the ASTROD I test mass due to interplanetary electrons in this simulation is about −11% of the cosmic-ray protons at solar minimum, and over −37% at solar maximum. In addition to the Monte Carlo uncertainty, an error of ±30% in the net charging rates should be added to account for uncertainties in the spectra, physics models and geometry implementations.     

The role of radial perpendicular diffusion and latitude dependent acceleration along the solar wind termination shock

A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E^−1.38 to E^−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E^−1.38 to E^−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.     

Observations of solar radio bursts from meter to kilometer wavelengths

Using the Clark Lake Radioheliograph data we present direct evidence that type III electron streams propagate in dense coronal streamers. We also present imaging observations of meter-decameter microbursts, which appear to be similar to those observed in hard X-rays. At meter-decameter wavelengths, these microbursts appear to be due to plasma radiation. From observations made with ISSE-3, we discuss the characteristics of hectometer and kilometer wavelength radio bursts. In particular, we show that from studies of type III storms that the exciter electrons propagate along spiral structures, where the density is enhanced and that there is an acceleration of the solar wind. We discuss type II bursts at kilometer wavelengths, compare them with meter type II bursts and discuss their association with interplanetary shocks. We show that the interaction between type III electron streams and shocks at kilometer wavelengths can provide information on the interplanetary shock geometry. Finally, we discuss the possibility that some shock associated (SA) events may be emissions caused by electrons accelerated lower in the atmosphere rather than high in the corona in type II shocks.Recent advances in solar research have resulted from new work on plasma radiation theory, new observations of active regions and flares across the electromagnetic spectrum and the availability of spacecraft in situ measurements of solar ejecta. In this paper, we review some results obtained with the Clark Lake multifrequency radioheliograph at meter-decameter wavelengths and from satellite multifrequency directive observations at hectometer and kilometer wavelengths. We present evidence that type III electrons propagate in dense coronal streamers, and that frequently observed microbursts (presumably of type III) at meter-decameter wavelengths are due to plasma radiation. We discuss observations of hectometer and kilometer type III radio storms which reveal information about active region structures, interplanetary magnetic field configuration, and solar wind acceleration. We also discuss kilometer type II bursts, interactions between type III electrons and interplanetary shocks, and present some new results on shock associated (SA) events.     

Enhancement in low-energy region ofa proton flux associated with interplanetary shock waves

We study time evolution of an energy spectrum of a proton flux in the range of 47 – 4750 keV for the energeticparticle event occurred on 255 DOY in 1999, which we consider as one of typical diffusive acceleration events associated with interplanetary shocks and irrespective of large X-ray solar flares. Fast enhancement during evolution is found in the range of less than about 0.5 MeV. Our previous numerical simulations using Stochastic Differential Equation method could not show this behavior, although we obtained results showing a power law energy spectrum, which suggesting that energetic particles are accelerated diffusively by shock waves, the first-order Fermi acceleration. We consider that less than 0.5 MeV protons need to exist to explain behavior of the observational energy spectrum and perform numerical simulations in order to investigate proper injection models for this event.     

Electron–ion temperature equilibration at collisionless shocks in supernova remnants

The topic of this review is the current state of our knowledge about the degree of initial equilibration between electrons, protons and ions at supernova remnant (SNR) shocks. Specifically, the question has been raised as to whether there is an inverse relationship between the shock velocity and the equilibration. This review aims to compile every method that has been used to measure the equilibration and every SNR on which they have been tested. I review each method, its problems and uncertainties and how those would effect the degree of equilibration (or velocity) inferred. The final compilation of observed electron to proton temperature ratios as a function of shock velocity gives an accurate, conservative picture of the state of our knowledge and the avenues we need to pursue to make progress in our understanding of the relation between the velocity of a shock and the degree of equilibration.     

Effects of the position of the solar wind termination shock and the heliopause on the heliospheric modulation of cosmic rays

The effects of changing the position of the solar wind termination shock and the position of the heliopause, and therefore the extent of the heliosheath, on the modulation of cosmic ray protons are illustrated. An improved numerical model with diffusive termination shock acceleration, a heliosheath and drifts is used. The modulation is computed in the equatorial plane and at 35 heliolatitude using recently derived diffusion coefficients applicable to a number of cosmic ray species during both magnetic polarity cycles of the Sun. It was found that qualitatively the modulation results for the different heliopause positions are similar although they differ quantitatively, e.g., clearly different radial gradients are predicted for the regions beyond the termination shock compared to inside the shock. The difference between the modulation for the two solar polarity cycles are less significant at a heliolatitude of 35° than in the equatorial plane. We found that moving the termination shock from 90 to 100 AU, with the heliopause fixed at 120 AU, caused only quantitative differences so that the exact position of the TS in the outer heliosphere seems not crucially important to global modulation. Moving the heliopause outwards, to represent the modulation in the tail region of the heliosphere, causes overall decreases in the cosmic ray intensities but not linearly as a function of energy, e.g., at 1 GeV the effect is insignificant. We conclude from this modelling that the modulation of protons in the heliospheric nose and tail regions are qualitatively similar although, clear quantitative and interesting differences occur.     

Relativistic proton production at the Sun in the 20 January 2005 solar event

Based on the concept of multiple acceleration of solar energetic particles (SEP) we analyzed the super-event of 20 January 2005 by the data of ground level, balloon and spacecraft observations. The main characteristics of relativistic solar protons (energy spectra, anisotropy directions and pitch-angle distributions) are derived and their dynamics during the event is studied. It is shown that the flux of relativistic solar protons may consist of two distinct components, the so-called prompt and delayed ones. Within a two-source model of particle generation, one of which is associated with an expanding magnetic loop, we solved the transport equation in energy phase space, including adiabatic losses simultaneously with the stochastic acceleration process, and calculate the expected spectra of the delayed component at the source. The confrontation of experimental spectra with theoretical ones shows that the delayed component may be correctly described by stochastic acceleration, but not the prompt component. The required acceleration efficiencies turned out to be rather high, so that, for this particular event, adiabatic cooling is practically negligible. Our results provide a new support to the existence of two populations of relativistic solar protons in some SEP events.     

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

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

A model for proton interaction in the atmosphere including fluxes of generated electrons

In this paper we analyse from a new point of view the energy deposition due to precipitated protons when they interact with the atmospheric components. The method described presents a different way of calculating the proton interaction and allows us to obtain the production rate and fluxes of the generated electrons as function of height and energy. Also the model gives the possibility of estimating independently the effects of protons and resulting secondary and tertiary electrons in protons events.     

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.     

HXMT卫星空间环境监测器及初步观测结果

硬X射线调制望远镜（HXMT）卫星是中国首个专门进行天文探测的空间科学实验卫星,运行于高度约550km、倾角约43°的低地球轨道.星载空间环境监测器为星上科学任务开展提供背景辐射实测资料.该监测器采用固体探测器望远镜系统和扇形阵列全新组合设计,可获取轨道空间高能质子和高能电子能谱、方向综合动态结果,给出更为全面的粒子辐射分布图像.初步探测结果显示,卫星运行轨道遭遇的带电粒子辐射集中分布在经度80°W-20°E,纬度0°-40°S的南大西洋异常区,粒子辐射在该区域表现出不同程度的方向差异分布,高能电子方向差异分布显著强于高能质子.2017年9月空间环境扰动期间,爆发的太阳质子事件并未对该轨道粒子辐射产生影响,而地磁活动导致该轨道穿越经度120°W-60°E,纬度40°-43°N的北美上空和经度60°-120°E,纬度43°-40°S的澳大利亚西南区域时遭遇增强粒子辐射影响,增强的粒子辐射表现出极强的方向分布.      

Energetic ion acceleration and transport in the upstream region of Jupiter: Voyager 1 and 2

Long-lived upstream energetic ion events at Jupiter appear to be very similar in nearly all respects to upstream ion events at earth. A notable difference between the two planetary systems is the enhanced heavy ion compositional signature reported for the Jovian events. This compositional feature has suggested that ions escaping from the Jovian magnetosphere play an important role in forming upstream ion populations at Jupiter. In contrast, models of energetic upstream ions at earth emphasize $\text{in situ}$ acceleration of reflected solar wind ions within the upstream region itself. Using Voyager 1 and 2 energetic (? 30 keV) ion measurements near the magnetopause, in the magnetosheath, and immediately upstream of the bow shock, we examine the compositional patterns together with typical energy spectra in each of these regions. We find characteristic spectral changes late in ion events observed upstream of the bow shock at the same time that heavy ion fluxes are enhanced and energetic electrons are present. A model involving upstream Fermi acceleration early in events and emphasizing energetic particle escape in the prenoon part of the Jovian magnetosphere late in events is presented to explain many of the features in the upstream region of Jupiter.