Evolution of electron energy spectrum during solar flares

Particle acceleration by direct current electric field in the current sheet has been extensively studied, in which an electric and a magnetic field are generally prescribed, and a power law distribution of the electron energy is obtained. Based on MHD numerical simulations of flares, this paper aims at investigating the time evolution of the electron energy spectrum during solar flares. It turns out that the model reproduces the soft–hard–hard spectral feature which was observed in some flares.     

Modeling of “gradual” solar energeticparticle events using a stochastic differential equation method

We have modeled “gradual” solar energetic particle events through numerical simulations using a StochasticDifferential Equation (SDE) method. We consider that energetic particle events are roughly divided into two groups: (1) where the shock was driven by coronal mass ejections (CMEs) associated with large solar flares, and (2) where they have no related solar events apart from the CMEs. (The detailed classification of energetic particle events was discussed in our previous paper.) What we call “gradual” solar energetic particle events belong to the former group. Particles with energies greater than 10 MeV are observed within several hours after the occurrence of flares and CMEs in many gradual events. By applying the SDE method coupled with particle splitting to diffusive acceleration, we found that an injection of high energy particles is necessary for early enhancement of such a high-energy proton flux and that it should not be presumed that the solar wind particles act as the seed population.     

利用质子强度变化对太阳质子事件进行短期预报方法初步研究

采用GOSE-10卫星4～9 MeV(P2),9～15 MeV(P3),15～40 MeV(P4),40～80 MeV(P5)能段上的质子通量数据,结合质子能谱,对太阳质子事件发生前各能谱参数的变化特征进行分析,详细介绍利用能谱参数的变化特征及能量E>10 MeV的质子通量数据对太阳质子事件进行预报的新方法,并运用这种方法对2002-2006年期间太阳质子事件进行了预报.预报结果显示,预报提前量最多达到100 h以上,对质子事件的报准率达97.5%,预报方法具备一定的有效性和实用性.      

太阳纯电子事件和质子-电子事件的电子能谱的机制

本文计算、分析了太阳耀斑加速电子在日冕中传输时激发的等离子体尾场的效应,认为耀斑电子的高能成份激发的尾场,能够加速低能耀斑电子,低能耀斑电子的能量增值可达几十keV至上百keV,这种尾场加速将软化约100keV以下的能量范围内(探测阈之上)的耀斑电子能谱。结合考虑尾场效应,本文提出了太阳耀斑加速电子从加速区到形成电子事件之间的能谱演化模式,说明了太阳纯电子事件的双幂律电子能谱和太阳质子-电子事件的单幂律电子能谱的形成,认为两类事件的电子能谱差异为耀斑电子日冕传输中不同程度的尾场效应所致,前者尾场效应弱,电子能谱呈双幂律,后者尾场效应较强,电子能谱为单幂律谱。      

First-order Fermi shock acceleration in solar flares

First order Fermi shock acceleration of electrons, protons and alpha particles is compared to observations of energetic particle events. For each event, a unique shock compression ratio produces spectra in good agreement with observation. The simple model predicts that the acceleration time to a given energy will be approximately equal for electrons and protons and, for reasonable solar parameters, can be less than 1 second to ～ 100 MeV.     

一个太阳硬X射线爆的观测特性

本文分析了1981年4月27日的一个特大高能爆发。它在硬X射线(HXR)、γ射线和微波(MW)记录上显示一一对应的多脉冲结构。文中讨论了这些脉冲的寿命、时延与HXR的能量、MW的波长之间的关系,发现时延量与湍动加速所预期的值相符合。提出了各脉冲期间的谱呈软→硬→软的演化,可能是由高能电子受到加速的看法。还对产生HXR的高能电子的谱指数、电子总数和MW源区的磁场也作了估计。      

Variability of spectra in large solar energeticparticle events

Using data from the Solar Isotope Spectrometer on the Advanced Composition Explorer obtained during 36large solar energetic particle events (SEPs) that occurred during 1997–2002 we have examined the spectral characteristics of oxygen and iron. Based on the shape of the oxygen spectrum during the decay phase following the peak in particle intensity, each SEP event was categorized as either exponential (7 events) or power law (29 events). We find that the exponential events were typically the larger events (in terms of peak oxygen intensity) and had Fe/0 ratios that strongly decreased with increasing energy.Event-averaged Fe/0 ratios (integrated over 12 to 60 MeV/nucleon) were at or below coronal abundances for nearly all these events, while the ratios obtained in the power law events were typically enhanced over coronal values. The majority of the power law events had the same spectral index for both 0 and Fe resulting in an Fe/0 ratio independent of energy. However 6 of the 29 power law events had Fe/0 ratios that increased with increasing energy due to an Fe spectral index less negative than that of 0. We consider simple diffusion theory in an effort to understand the nature of these events and obtain preliminary but promising results.     

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.     

太阳质子事件中短期预报模型研究

太阳质子事件通量的预测对航天器抗辐射加固设计和航天员出舱活动具有非常重要的意义.针对一年以下的航天任务,利用经验统计方法,确认太阳活跃年和太阳平静年期间,1——365天不同时间段内 > 10MeV,> 30MeV和 > 60MeV的太阳质子事件积分通量符合对数正态分布,且通量对数的标准偏差σ和期望值μ随任务期时间的变化满足对数函数形式.以此为基础,构建太阳质子通量的中短期预报模型.该模型能够针对太阳活跃年和太阳平静年,给出一定置信度下1——365天不同时间内 > 10MeV,> 30MeV和 > 60MeV的质子事件通量分布,从而为执行中短期航天任务提供太阳质子事件通量的预测,以规避不必要的风险.      

Detection of solar neutrons on a long duration balloon flight

The observation of large solar flares on high altitude balloons requires long duration balloon flights because large flares are infrequent and cannot be predicted with enough reliability and lead time to allow a conventional balloon to be launched and reach altitude before the flare occurs. With the many weeks at float altitude expected for a long duration flight, the probability of “catching” a large flare during solar maximum becomes reasonably high and the study of phenomena which heretofore have required a satellite become accessible to a balloon platform. One example of this type of experiment is the observation of neutrons produced by the interaction of flare accelerated nucleons with the solar atmosphere. Because the neutrons are produced immediately by the flare accelerated particles and are unaffected by their transmission through the upper solar atmosphere and the intervening magnetic fields, their observation at 1 A.U. will provide direct information on the flare acceleration process. Specifically, a measurement of the neutron energy and time spectra will yield the energy spectrum of the charged nucleons in the interval 50 to 500 MeV/amu, the charged particle anisotropy, the height of the acceleration region for limb flares, and information on the two-stage acceleration process. Because the γ-ray spectrum is also sensitive to these factors, a combined neutron and γ-ray measurement will provide a much more stringent test of flare models than either done separately. CWRU and the University of Melbourne have designed the EOSCOR (Extended Observation of Solar and Cosmic Radiation) detector to have the necessary sensitivity to detect neutrons from a flare 0.1 the size of the 4 Aug. 1972 event and to be compatible with the constraints of the long duration balloon system. The detector has been test flown on short duration balloon flights and calibrated at E_n = 38, 58, and 118 MeV. It is planned to launch it on a long duration balloon flight from Australia in December 1982 when simultaneous γ-ray observations will be possible with the SMM and/or HINTORI satellites.     

FEM code simulation of the magnetospheric proton fluxes

We report a study of the numeric solution to the diffusive transport equation for energetic protons magnetically trapped in the Earth's equatorial magnetosphere. The analysis takes into account the pertinent physical processes in this region, including deceleration of protons by Coulomb collisional interactions with free and bound electrons, the charge exchange process, cosmic ray albedo neutron decay source, and electric and magnetic radial diffusion. These results were obtained using the Finite Element Method with magnetic moment and geomagnetic L-shell as free variables. Steady state boundary conditions were imposed at L=1 as zero distribution function and at L=7 with proton distribution function extracted from ATS 6 satellite observations. The FEM-code yields unidirectional proton flux in the energy range of 0.1–1000 MeV at the equatorial top of the geomagnetic lines, and the results are found to be in satisfactorily agreement with the empirical NASA AP-8 model proton flux within the energy range of 0.5–100 MeV. Below 500 keV, the empirical AP-8 model proton fluxes are several orders of magnitude greater than those computed with the FEM-code at L<3. This discrepancy is difficult to explain by uncertainties of boundary spectrum parameters or transport coefficients.     

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.     

Change in the rigidity dependence of the galactic cosmic ray modulation in 2008–2009

Observations of charged particle fluxes in the stratosphere of the polar regions represent the cosmic rays variations with energy above 100 MeV. At the end of 2009 these fluxes reached the highest level for the time of observations from mid 1957 and were by 17% higher than the previous extremum value of May 1965. In the mean time the ground-based neutron monitors showed the remarkably less count rate enhancement. These results argue for the significant change in the energy spectrum of incoming particles in 2008–2009 in the energy range of ∼100–1500 MeV/n.     

Solar flare particle injection spectra and spectra of flux maxima at the point of observation

Numerical models of impulsive solar flare particle events usually assume the radial diffusion coefficient to be independent of energy per nucleon, T, although the observations indicate a T^0.5 dependence (constant mean free path). The assumption of a constant diffusion coefficient results in a preservation of a power law injection spectrum at all radial distances throughout the event. We investigate the effect of an energy dependent diffusion coefficient on the spectrum of flux maxima at a fixed point in interplanetary space. This spectrum is harder than that of initial differential number densities close to the sun. Furthermore, the spectrum hardens with increasing radial distance which seems to be at variance with observations.     

（20～400）MeV准单能中子参考辐射场的建立方法CSCD

为了满足对高能中子辐射场测量仪器进行校准的需要,必需建立20MeV以上的准单能中子参考辐射场.对国际上（20～400） MeV准单能中子参考辐射场的建立情况及建立方法进行了调研,可通过7Li（p,n）7Be反应获得该能区准单能中子场,中子束流轮廓可通过转换体和影像板进行测量,辐射场中子能谱可利用液体闪烁探测器或者238U裂变电离室通过飞行时间法进行测量,可建立高效率的反冲质子望远镜作为该能区的中子注量测量初级标准.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

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.     

Monte Carlo calibration of the response of the University of Chicago’s Cosmic Ray Nuclei Experiment (CRNE) on IMP-8 to electrons above 0.5 MeV

Modern instrument-simulation techniques offer the possibility of increasing the scientific yield from archival space datasets. In this paper, we report on a simulation of the electron response of the University of Chicago’s Cosmic Ray Nuclei Experiment (CRNE) instrument on the IMP-8 satellite. IMP-8/CRNE returned data from 1973 to 2006. The CRNE particle telescope was designed to measure the isotopic composition of Galactic cosmic-ray (GCR) nuclei and has also been used in many studies of protons and ions above 10 MeV/nucleon from solar energetic particle (SEP) events. But CRNE also functions as a highly-capable detector for solar electrons above 0.5 MeV, an energy range that has not been extensively studied. Utilization of the CRNE electron data has heretofore been limited by the fact that CRNE was never calibrated for electrons. We have therefore used the GEANT4 Monte Carlo simulation package to model the CRNE response to electrons and (separately) protons for multiple energies and incident angles. The results were used to compute the energy- and angle-dependence of the effective area and the energy-dependence of the geometric factor. The response to protons, which was already well understood, was used to verify the mass model, the simulation settings, and the post-processing software. Our simulation of the IMP-8/CRNE electron response now allows analysis of hundreds of relativistic solar electron events observed by CRNE over the years, including studies of evolution of electron energy spectra with high time resolution. We show examples of these results and briefly discuss potential applications to future scientific investigations.     

Solar particle precipitation into the polar atmosphere and their dependence on hemisphere and local time

The precipitation of solar energetic particles, protons as well as electrons, at high latitudes is commonly assumed to be homogeneous across both polar caps. Using Low-Earth Orbit POES (Polar Orbiting Environmental Satellites) we determine particle penetration ratios into the polar atmosphere for protons ranging from about 0.1 MeV to 500 MeV and for electrons spanning about one order of magnitude in energy with a maximum of 0.3 MeV. Based on power law fits for the POES spectrum we show, that for energies interesting for middle and lower atmosphere chemistry, particle flux over the poles is comparable in magnitude to flux at the geostationary orbit or at L1 in interplanetary space. The time period under study are the solar energetic particle (SEP) event series of October/November 2003 and January 2005.     

Particle acceleration and gamma rays in solar flares: Recent observations and new modeling

Experiments on SMM, GAMMA, Yohkoh, GRANAT, Compton GRO, INTEGRAL, RHESSI and CORONAS-F satellites over the past three decades have provided copious data for fundamental research relating to particle acceleration, transport and energetics of flares and to the ambient abundance of the solar corona, chromosphere and photosphere. We summarize main results of solar gamma-astronomy (including some results of several joint Russian–Chinese projects) and try to appraise critically a real contribution of those results into modern understanding of solar flares, particle acceleration at the Sun and some properties of the solar atmosphere. Recent findings based on the RHESSI, INTEGRAL and CORONAS-F measurements (source locations, spectrum peculiarities, ³He abundance etc.) are especially discussed. Some unusual features of extreme solar events (e.g., 28 October 2003 and 20 January 2005) have been found in gamma-ray production and generation of relativistic particles (solar cosmic rays, or SCR). A number of different plausible assumptions are considered concerning the details of underlying physical processes during large flares: (1) existence of a steeper distribution of surrounding medium density as compared to a standard astrophysical model (HSRA) for the solar atmosphere; (2) enhanced content of the ³He isotope; (3) formation of magnetic trap with specific properties; (4) prevailing non-uniform (e.g., fan-like) velocity (angular) distributions of secondary neutrons, etc. It is emphasized that real progress in this field may be achieved only by combination of gamma-ray data in different energy ranges with multi-wave and energetic particle observations during the same event. We especially note several promising lines for the further studies: (1) resonant acceleration of the ³He ions in the corona; (2) timing of the flare evolution by gamma-ray fluxes in energy range above 90 MeV; (3) separation of gamma-ray fluxes from different sources at/near the Sun (e.g., different acceleration sources/episodes during the same flare, contribution of energetic particles accelerated by the CME-driven shocks etc.); (4) asymmetric magnetic geometry and new magnetic topology models of the near-limb flares; (5) modeling of self-consistent time scenario of the event.