太阳高能粒子事件强度与相关双CME事件关系的两个事例

太阳高能粒子事件常伴随太阳耀斑和日冕物质抛射事件(Coronal Mass Ejections,CME)出现,由于太阳高能粒子事件的关键因素是双CME的相互作用,利用SOHO卫星观测的高能粒子强度、耀斑强度以及CME的相对高度与时间,通过高度与时间拟合得到的速度,分析了2001年4月15日和2005年1月20日的太阳高能粒子事件强度与相关双CME事件的关系,发现这两个太阳高能粒子事件中E ≥ 10MeV质子的强度与双CME事件无关.因此在这两次太阳高能粒子事件早期,E ≥ 10MeV质子的强度只与相关太阳耀斑和CME有关.      

空间高能粒子与器件布线层核反应后次级粒子LET分布研究

空间高能质子和重离子是导致元器件发生单粒子效应的根本原因,为准确评估元器件在轨遭遇的单粒子效应风险,必须清楚高能质子、重离子与器件材料发生核反应的物理过程及生成的次级重离子LET（Line EnergyTransfer）分布规律。针对典型CMOS工艺器件模拟计算了不同能量质子和氦核粒子在器件灵敏单元内产生的反冲核、平均能量及线性能量转移值,并分析了半导体器件金属布线层中重金属对次级重离子LET分布的影响规律。计算结果表明：高能粒子与器件相互作用后产生大量次级重离子,且高能质子作用后产生的次级粒子的LET值主要分布为0～25MeV·cm²/mg;高能氦核粒子作用后产生的次级粒子的LET值主要分布为0～35 MeV·cm²/mg;有重金属钨（W）存在时能提高次级粒子的LET值,增加了半导体器件发生单粒子效应的概率,该研究结果可为元器件单粒子效应风险分析、航天器抗单粒子效应指标确定提供重要依据。     

太阳高能粒子事件上升时间统计研究

选取1997-2006年共66个较大的缓变型太阳高能粒子(SEP)事件, 分析了不同条件下太阳高能粒子通量廓线上升时间与源区日面经向分布之间的相关关系, 研究了日冕物质抛射(CME)和耀斑在SEP上升阶段的作用特点.统计结果表明,大SEP事件的源区主要分布在太阳西半球, 特别是磁足点东西两侧45°范围内; 在高速太阳风条件下, 低能通道的通量上升时间与日面相对经度有较好的相关性,即离磁足点越远, 上升时间越长,而高能通道相关性则不明显; 全晕状CME产生的SEP事件对应的上升时间与源区位置没有明显的相关性, 而部分晕状CME伴随的SEP事件则与二次拟合曲线符合很好.分析表明,在缓变型SEP事件的通量上升阶段, 耀斑加速过程起着重要作用,这在部分晕状CME伴随的SEP事件中尤为显著.      

空间粒子环境对单粒子效应影响的比较

空间的高能重离子和高能质子都能引起单粒子效应。几种典型环境中的重离子和质子对不同器件单粒子效应的影响的比较表明，在地磁捕获质子强度较大的区域，捕获质子对单粒子效应的作用必须考虑；而在该区域以外，对单粒子效应的发生有显著影响的是宇宙射线的重离子而不是占绝大多数的质子。     

2 AU以内的“渐进型”太阳高能粒子事件模拟

太阳高能粒子（Solar Energetic Particle,SEP）事件是影响地球空间以及深空辐射环境的主要因素之一。“渐进型”太阳高能粒子事件中的高能粒子主要来自于日冕物质抛射（Coronal Mass Ejection,CME）所驱动的激波扩散加速（Diffusive Shock Acceleration,DSA）过程。CME驱动的激波在行星际的传播过程中,其结构不断演化,进而影响到高能粒子的加速过程。本文利用二维太阳高能粒子加速和传播模型,对发生于2014年4月18日的太阳高能粒子事件实例进行了数值模拟。模型考察了黄道面上2 AU的距离以内包含地球所在位置的4个不同点,分别计算了每个点上高能粒子的通量。数值模拟的结果表明：黄道面内不同位置的观察点,与激波波前的磁力线连接不同,从而导致观察点处高能粒子的通量有着显著的差异。该模型的计算结果可以为深空探测计划开展辐射环境研究提供必要的输入。     

高速太阳风中α粒子分布的形成

以多级波粒相互作用模型和1维电磁混合模拟方法研究了高速太阳风中α粒子分布的形成机制．认为α粒子从其产生区到1AU,经历4个特征区域：在冕洞下方的小尺度磁重联区,α粒子产生并得到加速：在约十几个太阳半径日冕区,α粒子形成球壳分布函数：在几十个太阳半径处,在Alfven波和温度各向异性不稳定性的作用下,α粒子分布函数演化为漂移Maxwell分布：当α粒子进入0．3－1AU时,Alfven波的作用使α粒子速度超过太阳风速度约1个局域Alfven速度．     

磁场重联的二维粒子模拟研究

利用二维全粒子模拟方法研究了无碰撞等离子体中的磁场重联过程，得到了不同区域的离子和电子速度分布．计算结果表明，电子和离子在扩散区中的不同动力学特性产生的Hall电流使磁场的y分量By呈现四极形分布．离子和电子的速度分布偏离了初态时的Maxwell分布，呈现非局域的多重分布．同时由于磁场重联而产生的电场使电子在X点附近得到加速和加热，因而在电子的能谱分布中形成-高能尾。     

关于磁尾高能粒子脉冲的一些特性

利用行星际监察卫星IMP-J取得的高能粒子探测数据(质子能档P4:230keV>E>160keV)与极光电激流指数AE作相关分析,在地心太阳磁层坐标下,按照Fairfield关于中性片对地心太阳磁层"赤道面"的偏离模式,把磁尾分成三个区域:中性片区域、低纬区域和高纬区域。结果表明:(1)高能粒子脉冲的平均强度在中性片区域最强,低纬次之,高纬最弱,表明高能粒子脉冲源区在中性片区域;(2)中性片附近,粒子脉冲和AE指数相关最好,达0.59,低纬次之,高纬几乎无相关,表明粒子脉冲与亚暴事件有关,它是磁尾中性片附近磁力线重联产生的感应电场加速的结果;(3)粒子加速区局限于中性片附近的薄层内,与国外结果完全一致。      

质子单粒子翻转截面计算方法

在分析质子与硅反应的基础上,建立了质子单粒子翻转截面理论计算模型,提出了模拟计算方法。计算得到了不同能量的高能质子与硅反应产生的次级粒子种类、截面、能谱和双微分截面。采用Monte Carlo方法模拟质子与硅的反应;应用TRIM程序计算次级粒子的射程;计算得到次级粒子在存储单元的灵敏区内沉积的能量,产生的电荷。通过与临界电荷的比较,判断是否导致单粒子翻转,从而得到单粒子翻转截面。计算得到的单粒子翻转截面与实验数据符合较好。     

太阳高能粒子观测特征经向分布统计

基于多卫星联合观测数据,筛选了2006年12月至2017年10月期间122个太阳高能粒子(SEP)事件及其伴随的日冕物质抛射(CME),分析了SEP事件属性随相对经度的变化、与CME属性之间相关性的经向分布以及与Fe/O比值的关联.研究结果显示:低Fe/O类事件的峰值通量Ip通常更高,伴随CME更大,但通量上升速度较慢...     

Direct observation of charge state abundances of energetic He,C, O,and Fe emitted in solar flares

The ionic charge states of helium, carbon, oxygen, and iron have been determined for three solar particle enhancements by an electrostatic deflection analyzer, which is combined with a thin window proportional counter dE/dX vs. E system. The observations are obtained during the periods September 23 to 29, 1978, June 6 to 8, 1979, and September 15 to 26, 1979, with an instrument onboard the ISEE-3 spacecraft. The mean ionic charge states for He, C, and oxygen exhibit a high degree of ionization with values of Q = 2, 6, and 7.2, respectively. The charge state of iron is near 13 charge units. Variations from flare to flare and within the September 23, 1978 flare are small. The most surprising feature of the charge state measurement is the observation of a small (～10%) but finite contribution of singly ionized helium.     

Propagation of energetic particles in the solar wind

The characteristics of solar energetic particle events as observed in interplanetary space depend on many physical processes acting at the source and during the transport from the source to the observer. These processes, such as acceleration at the Sun and the propagation near the Sun and in interplanetary space depend, in general, on both the particle velocity and rigidity. Thus, the evaluation of both the nuclear charge and/or atomic mass and the ionic charge of heavy ions turns out to be essential for the interpretation of the physical parameters observed, such as the energy spectra and the compositional variations during individual solar energetic particle events. In this paper recent results on the direct determination of the charge states of He, C, O, and Fe will be summarized. Using these results the compositional variations during individual solar particle events will be discussed. It will be shown that ratio changes by a factor of ～ 10 during the onset phase of solar particle events, as frequently observed, could be explained not only by rigidity dependent interplanetary propagation, but also by rigidity dependent diffusive propagation in the corona. However, there is now increasing experimental evidence that also other processes, such as compositional variations at the source and discontinuities of the interplanetary magnetic field, separating two different particle populations, may be important. Thus the picture emerges that these variations do not have a unique explanation but rather that each event has to be investigated individually.     

Composition of energetic particles from solar flares.

We present a model for composition of heavy ions in the solar energetic particles (SEP). The SEP composition in a typical large solar particle event reflects the composition of the Sun, with adjustments due to fractionation effects which depend on the first ionization potential (FIP) of the ion and on the ratio of ionic charge to mass (Q/M). Flare-to-flare variations in composition are represented by parameters describing these fractionation effects and the distributions of these parameters are presented.     

Towards multi-dimensional space weather modeling for energetic oxygen ions in the earth''s inner magnetosphere: Equilibrium configuration

We report the first 3+1 dimensional model development for energetic atomic oxygen ions in the Earth's radiation belts. Energetic Oxygen ions cans be supplied to the Earth's Inner magnetosphere from the sun (as a component of solar wind and solar energetic particles), from anomalous cosmic rays, and from acceleration processes acting on ionospheric atomic oxygen ions. We have built a multi-dimensional oxygen ion model in the following free parameters: geomagnetic L-shell, the magnetic moment, the second adiabatic invariant, and the discrete charge state number. Quiet time, steady state oxygen ion distributions have been obtained numerically from an assumed outer radiation zone boundary condition at L=7, average values of the radial diffusion coefficients, and standard values for the exospheric neutral densities due to the MSIS-86 upper atmosphere and exosphere neutral thermal particle density model. Average distributions of free electrons in the plasmasphere were also assumed with a mean plasmapause location just beyond L=4. We included the six lowest ionic charge states of atomic oxygen (¹⁶O) based on an existing charge exchange cross section compilation by Spjeldvik and Fritz (1978). Computed oxygen ion distributions include the resulting equilibrium structure of energy oxygen ions between 10 KeV and 100 MeV.     

耀斑脉冲相离子电离的弛豫过程

观测表明:在耀斑前活动区上空的日冕中普遍存在低温的磁环结构, 其中一些这种磁环结构的温度在105K以下;同时, 耀斑脉冲相存在非常迅速的加热过程.本文在这些观测事实的基础上, 提出了一个太阳耀斑爆发时离子电离的模型, 并由这一模型研究了耀斑时离子电离的弛豫过程.在太阳宇宙线的源物质来自高色球层的条件下, 得到了与观测结果符合得较好的太阳耀斑能量离子的平均电荷及其电离态分布.      

Long-term modulation of Galactic Cosmic Radiation and its model for space exploration.

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop in exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertain ties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa (r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa (-r,t)=K(t)kappa 0 beta P kappa 1(r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that kappa(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and 14C measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.     

Particle charge interchange during acceleration in flare regions

We have examined the conditions for the establishment of charge equilibrium of solar particles during their acceleration. We derive criteria for charge interchange with the atomic and ionized hydrogen at the particles'sources, for two different acceleration mechanisms. It is found that charge interchange is established whenever a particle event is produced. The implications related to mass and charge spectra of particles are discussed. The measured charge state of solar particles cannot in general be directly used for diagnosis of the source temperature, so we suggest another alternative based on the emitted radiation from electron capture.     

太阳宇宙线的有效电荷及计算

本文讨论了太阳宇宙线成份的高电离态的产生, 提出了计算平均有效电荷的公式, 结果与观测资料对比能很好地符合.      

On compositional variations of heavy ions during solar particle events

Intensity-time profiles of protons, alpha particles, and heavy ions (C, O, Fe) in the MeV/nucleon energy range have been analyzed for one solar particle event following the solar flare on September 23, 1978. The data have been obtained with the wide angle double dE/dx-E sensor of the Max-Planck-Institut/University of Maryland experiment onboard ISEE-3. We found time variations in the iron to helium ratio of up to 2 orders of magnitude and a significant variation of the O/He ratio during this event, whereas the C/O-ratio at the same energy/nucleon appears to be time independent. We investigated the influence of a rigidity dependent mean free path in interplanetary space and of rigidity dependent coronal propagation on heavy ion ratios during solar particle events. We found that both the magnitude and time scale of the ratio changes during the September 23 event cannot be explained by rigidity dependent interplanetary or coronal propagation alone. These ratio changes are probably caused by multiple injection at the sun.     

Energetic particles during the first two years of SOHO science mission

Energetic particle intensities observed by ERNE instrument onboard SOHO spacecraft during the first two years of SOHO science mission have been analysed and compared to observations of IMP 8 satellite around two earlier sunspot minima. During an eight month period around the latest sunspot minimum, which occurred in October 1996, energetic particle intensities stayed at a lower level than during any equivalent period around solar minima of July 1976 and June 1986. During the period from March to October, 1996, there was not a single day, when the daily averaged intensity of 1.6–3 MeV protons exceeded the level of 1 proton/(cm² sr s MeV). Also monthly counts of grouped solar flares and mean monthly sunspot number were lower than during minima of 1976 and 1986.