带电粒子在近地球区运动的计算方法

基于磁层粒子动力学理论,根据单粒子方法和偶极磁场模型,在近地球区(L＜10),详细讨论了带电粒子的运动特征以及适合不同能量的带电粒子的计算方法,定量分析了磁层中的电场和磁场对各种能量的带电粒子漂移运动的影响.结果表明,对于能量低于105 eV的电子和102 eV的质子,宜采用引导中心近似方法;对于能量介于105～108 eV之间的电子和102～108 eV之间的质子,可在部分区域内采用引导中心近似方法,若运动区在10 Re内,则只能采用变步长的轨道法;而对于引导中心方法,粒子能量高于105 eV时可以忽略电场漂移的影响,粒子能量低于103 eV时不必考虑磁场漂移的影响,从而简化了引导中心方程组,提高了数值计算效率.      

磁扰期间的三维漂移壳分离

基于磁层粒子动力学理论,首先对比了计算漂移壳分离的引导中心法和磁力线追踪法,计算表明两种方法的计算结果一致.然后分别采用T89c和T96磁层磁场模式,用磁力线追踪法数值计算了不同初始位置(≤9Re)、不同初始投掷角、不同Kp指数和不同太阳风压力下,带电粒子的漂移壳分离.计算结果揭示了漂移壳分离随初始位置、投掷角、Kp指数和太阳风压力的变化.其具体特征如下. (1)随着径向距离的增大,漂移壳分离效应愈加显著,由正午出发的粒子将被稳定捕获,而午夜出发的径向距离≥7Re的部分大投掷角粒子将沿磁层顶逃逸. (2)正午出发的粒子,漂移到午夜时其漂移壳随投掷角减小向外排列;午夜出发的粒子,漂移到正午时其漂移壳随投掷角增大排列; 90°投掷角粒子在磁赤道面的漂移壳沿着磁场等值线排列. (3)漂移壳分离随Kp指数和太阳风压力增大变得显著,且随这两种扰动参数的变化特征和趋势是基本相似的.      

In situ observations of magnetic field fluctuations

It is well known that the irregularities of the magnetic field are intimately related to the motion of charged particles. Although transport theories need the spatial and time variations of the magnetic field as input, in situ observations are very limited. Ulysses observations have provided a major step forward by entering the unexplored high latitude regions of the heliosphere, the knowledge of which is vital to interpret particle flux measurements, even at the ecliptic. We analyze the magnetic field data of Ulysses during the mission to study the waves and discontinuities in the heliosphere at different locations, covering a total sunspot cycle. Various tools are employed, including power spectral and structure function analysis. A remarkable difference was found between the fluctuations in the fast and slow solar wind. We argue that the latitudinal extent of the high speed solar wind contributes significantly to the latitudinal variation of the transport parameters, which should also affect the 11 (and 22) year modulation cycle.     

The devil is in the details: Compact structures in pulsar wind nebulae

The large-scale structure of pulsar wind nebulae (PWNe) tells us a considerable amount about their average magnetic fields, the total particle input from the pulsar winds, and the confining pressure at their outer boundaries. However, the details of the pulsar outflow, the sites of shocks and particle acceleration, the effects of instabilities in the magnetic field, and the interaction between the relativistic wind and the surrounding ejecta are contained in small-scale structures, where we observe jets and toroidal structures, time-varying emission from compact clumps, and filaments in both the inner and outer regions of the nebulae. Here, I review recent observational studies of compact structures in PWNe and present current scenarios (and questions) regarding their origin.     

星间洛仑兹力编队飞行的平衡点及零速度曲面

提出了一种新的利用星间洛仑兹力控制卫星相对运动的方法:使主星产生自旋磁场,副星带电,通过控制副星所受的星间洛仑兹力进行编队。假设主星产生的人造磁场表现为偶极子并且运动在一条开普勒圆形轨道上,副星恒定带电并在主星附近运动,同时假设星间洛仑兹力只影响副星的运动而不影响主星的运动。推导了副星在主星HCW(当地垂直当地地平坐标系)坐标系下的相对运动动力学方程。针对偶极子与HCW坐标轴X轴重合的情况下,推导了动力系统下的平衡点(并采用稳定性分析方法分析其线性化意义下的稳定性)、积分常数和零速度曲面,证明了有界相对运动的存在。最后用数值仿真验证了上述结论。     

Dust in magnetised plasmas: Basic theory and some applications

In this paper we review the theory of charged test particle motion in magnetic fields. This theory is then extended to charged $\text{dust}$ particles, for which gravity and charge fluctuations play an important role. It is shown that systematic drifts perpendicular to the magnetic field and stochastic transport effects may then have to be considered - none of which occur in the case of atomic particles (with the exception of charge exchange reactions). Some applications of charged dust particle transport theory to planetary rings are then briefly discussed.     

Field-aligned currents and related phenomena in the cusp

This review deals with observations of field-aligned currents (FAC) based on measurements of magnetic and electrical fields, particle fluxes (0.2–15 keV) and plasma parameters, provided by the INTERCOSMOS BULGARIA-1300 satellite. Special attention is paid to the cusp sector (10–14 MLT). The current regions inferred from the satellite magnetometer are compared with the charecteristics of the dayside oval and the cusp derived from a variety of parameters. Our topics are as follows: penetration of region 1 in the cusp; FAC reversals and plasma convection; thermal plasma structure and collocation of FAC; Associations with the IMF components etc.     

Magnetic connectivity and solar energetic proton event intensity profiles at deka-MeV energy

We present an analysis of the time-intensity profiles of 25 solar energetic proton events at 18.2 MeV, modelled by fitting an analytical function form (a modified Weibull function) to the observed intensities. Additionally relying on previous work that characterized the magnetic connectivity between the event-related solar flare and the observer in these events with three angular parameters, we investigate the fit function parameters, the connectivity parameters, and the iron-to-carbon ratio of the events for dependencies and correlations. We find that the fit parameter controlling the basic shape of the profile (parameter a) is not clearly dependent on the connectivity parameters or the Fe/C ratio, suggesting that the profile shapes of neither well and weakly connected nor generally “impulsive” and “gradual” events differ systematically during the early stages of the event at 1 AU. In contrast, the time scaling of the fit function (parameter b) is at least moderately correlated with both the magnetic connectivity parameters and the Fe/C ratio, in that well-connected and iron-rich events are typically shorter in relative duration than weakly connected and nominal-abundance events; intensity rise times display a similar correlation with the connectivity parameters. We interpret the former result as following from the combined effect of various transport processes acting on the particles in interplanetary space, while the latter is essentially consistent with established knowledge regarding the observed dependence of the time-intensity profile shapes of solar energetic particle events on their magnetic connectivity and heavy ion abundances. The desirability of modelling the particle transport effects in detail and extending the analysis to cover higher energies is indicated.     

World grid of cosmic ray vertical cut-off rigidity for the last decade

Cosmic ray cut-off rigidity tables and maps over the world concerning the epochs 2010, 2015 and the current one 2020 have been constructed. These maps display the effective cut-off rigidity in every five degrees in latitude and in longitude at the altitude of 20 km above the surface of the international reference ellipsoid. The values of the geomagnetic cut-off rigidity were calculated in every 5° in latitude and in every 15° in longitude applying the well-known method of particle trajectory calculations resulted from the theory of the particle motion in the Earth's magnetic field. The applied software employed the 12th Generation of the International Geomagnetic Reference Field (IGRF 12) and trajectories were calculated at 0.01 GV intervals in order to determine the vertical cut-off rigidity for each location. Beyond the use of the calculated cut-off rigidity values as a basic reference of charged particle access to different geographical locations during quiet and/or more intense geomagnetic periods, these results can be used for a long- term forecasting of the geomagnetic conditions variations.     

Energization of charged particle in a time-dependent chaotic magnetic field with an implication of the production of seed particles in solar energetic particle events

We investigate the acceleration of charged particles in a time-dependent chaotic magnetic field in this work. In earlier works, it has been demonstrated that in an asymmetric wire-loop current systems (WLCSs), the magnetic field is of chaotic in nature. Furthermore, observations also showed that there exist time-varying current loops and current filaments in solar corona. It is therefore natural to conceive that the magnetic field on the solar surface is chaotic and time-dependent. Here, we develop a numerical model to study the acceleration process of charged particles in a time-varying chaotic magnetic field that is generated by an ensemble of 8 WLCSs. We found that the motion of energetic particles in the system is of diffusive in nature and a power law spectrum can quickly develop. The mechanism examined here may serve as an efficient pre-acceleration mechanism that generates the so-called seed particles for diffusive shock acceleration at a coronal mass ejection (CME) driven shock in large solar energetic particle (SEP) events.     

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.     

Changes in cosmic ray propagation induced by corotating interaction regions

By simulating the trajectories for scatter free and diffusive propagation of relativistic cosmic rays in a model of the heliospheric magnetic fields containing a representation of corotating interaction regions (CIR's) we find that, due to the large gradients associated with these compression regions, the motion is strongly affected and differs substantially from the predictions of current modulation theory. For positive (outward) northern hemisphere polarity particles do not stream purely from high latitudes but can come from almost any latitude in the outer heliosphere; for negative polarity many particles come along the current sheet (as predicted) but a second, equally important, poipulation exists comprising particles that do not start on the current sheet but are brought to low latitudes by their interaction with CIR's. Thus, we conclude that CIR's (and other large scale structures) cannot be ignored in analyses of cosmic ray modulation.     

Modeling the radiation belts: what are the important physical processes to be taken into account in models?

During active periods, physical processes acting on particle dynamics are well known to have several origins and can play a role at different times, different locations, and with different time scales. As their effects can be opposite, it is necessary to identify them and their relative importance:

• the particle sources (plasmasheet or direct solar wind entry)

• the particle losses (precipitation or drift loss)

• the particle transport and acceleration (convection access, magnetic or electric pulse or variation and recirculation).

The various phenomena are explained and results are presented. In particular, we demonstrate that classical diffusion models like the Salammbo code can account for all these phenomena.     

带电粒子在中性线磁场中运动的解析轨道

本文我们计算了带电粒子在中性线磁场中运动的解析轨道。其结果是:(1)带电粒子在中性片磁场中的运动是粒子在中性线磁场或在具有北向分量的中性片磁场中的第一级近似形式。(2)带电粒子在中性片磁场中的解析轨道的第三级近似形式与电子计算机计算的数值轨道基本相同。它们仅仅在小扰动区与非小扰动区的交界线上出现一些偏差。(3)带电粒子在整个中性片磁场的运动可以分成三种形式。粒子一方面在垂直于磁场的平面上作闭合的周期性轨道运动, 同时闭合轨道的中心还沿着垂直于磁场平行于中性线方向漂移。另一方面粒子还沿磁力线方向做等速运动。(4)在小扰动区中粒子的闭合轨道是一个圆轨道, 但在非小扰动区中却是一个“8”字形轨道, 其漂移速度与小扰动区漂移方向相反, 其大小也比小扰动区漂移大很多。以上结果本文都给出一个完整的解析形式。      

Emission heights of coronal bright points on Fe XII radiance map

The study of coronal bright points (BPs) is important for understanding coronal heating and the origin of the solar wind. Previous studies indicated that coronal BPs have a highly significant tendency to coincide with magnetic neutral lines in the photosphere. Here we further studied the emission heights of the BPs above the photosphere in the bipolar magnetic loops that are apparently associated with them. As BPs are seen in projection against the disk their true emission heights are unknown. The correlation of the BP locations on the Fe XII radiance map from EIT with the magnetic field features (in particular neutral lines) was investigated in detail. The coronal magnetic field was determined by an extrapolation of the photospheric field (derived from 2-D magnetograms obtained from the Kitt Peak observatory) to different altitudes above the disk. It was found that most BPs sit on or near a photospheric neutral line, but that the emission occurs at a height of about 5 Mm. Some BPs, while being seen in projection, still seem to coincide with neutral lines, although their emission takes place at heights of more than 10 Mm. Such coincidences almost disappear for emissions above 20 Mm. We also projected the upper segments of the 3-D magnetic field lines above different heights, respectively, on to the tangent x–y plane, where x is in the east–west and y in the south–north direction. The shape of each BP was compared with the respective field-line segment nearby. This comparison suggests that most coronal BPs are actually located on the top of their associated magnetic loops. Finally, we calculated for each selected BP region the correlation coefficient between the Fe XII intensity enhancement and the horizontal component of the extrapolated magnetic field vector at the same x–y position in planes of different heights, respectively. We found that for almost all the BP regions we studied the correlation coefficient, with increasing height, increases to a maximal value and then decreases again. The height corresponding to this maximum was defined as the correlation height, which for most bright points was found to range below 20 Mm.     

耀斑后环珥系速度场的定量研究

本文研究了在太阳重力、偶极磁压强梯度和气体压力梯度联合作用下环珥系内物质的下落运动,并在等温、定常假设下用数值方法计算了该环系的二维视向速度场。通过将速度场的理论计算跟文[1]中给出的观测速度场拟合,导出了该环系的有关参数及分布。计算表明环系内物质密度和磁场强度对物质下落运动影响较显著。     

Flux emergence,flux imbalance,magnetic free energy and solar flares

Emergence of complex magnetic flux in the solar active regions lead to several observational effects such as a change in sunspot area and flux embalance in photospheric magnetograms. The flux emergence also results in twisted magnetic field lines that add to free energy content. The magnetic field configuration of these active regions relax to near potential-field configuration after energy release through solar flares and coronal mass ejections. In this paper, we study the relation of flare productivity of active regions with their evolution of magnetic flux emergence, flux imbalance and free energy content. We use the sunspot area and number for flux emergence study as they contain most of the concentrated magnetic flux in the active region. The magnetic flux imbalance and the free energy are estimated using the HMI/SDO magnetograms and Virial theorem method. We find that the active regions that undergo large changes in sunspot area are most flare productive. The active regions become flary when the free energy content exceeds 50% of the total energy. Although, the flary active regions show magnetic flux imbalance, it is hard to predict flare activity based on this parameter alone.     

Analytic and numerical treatment of motion of dust grain particle around triangular equilibrium points with post-AGB binary star and disc

This paper investigates the motion around the triangular equilibrium points, of a passively gravitating dust particle in the gravitational field of a low-mass post-AGB binary system, surrounded by circumbinary disc. The two bodies of the binary are modeled as a triaxial star and a radiating-oblate star. Due to small deviation of disc stars on circular orbits, we have assumed that the Coriolis and centrifugal forces of the stars are slightly perturbed. The triangular equilibrium points of the particle are found. These points are defined by, triaxiality of the primary star, oblateness and radiation of the secondary one and the gravitational potential from the disc mass. Further, when the disc mass increases, the particle moves nearer to the stars and farther away from the disc. In general, these equilibrium points are linearly stable when μ < μ_C; where μ is the mass ratio and μ_C is the critical mass function, defined by the parameters of the system. The effects of each of these parameters on the size of the stability region are stated, and the periodic motion around the stable points is examined. It is seen that the orbits are ellipses, and the orientation, eccentricities, lengths of the semi-major and semi-minor axes are influenced by the parameters of the problem. In particular, for our numerical linear stability analysis, we have taken an extremely depleted pulsating star, IRAS 11472-0800 as the post-AGB triaxial star, with a weakly-radiating young white dwarf star; G29-38 as the secondary. For this system, the stability result of the triangular points comes out different. Here, μ_C < μ throughout the entire range of the mass ratio and the critical mass function. Hence, the triangular equilibrium points are unstable. The stability of the orbits is tested using the Poincaré surfaces of section (Pss). The region of stability is controlled by the introduced parameters and the Jacobi constant.     

Particle transport in a vortex medium

We study numerically particle transport in a two-dimensional coherent vortex field. Reasonable agreement exists between previously derived radial transport coefficients for energetic particles (Verkhoglyadova, O.P., le Roux, J.A. Particle diffusion on vortices in nearly incompressible magnetohydrodynamics. Astrophys. J. 602, 1002–1005, 2004a; Verkhoglyadova, O.P., le Roux, J.A. Cosmic ray transport in a vortex flow. IGPP-UCR Conf. Physics of the Outer Heliosphere (Riverside), AIP Conf. Proc., pp. 243–248, 2004b) and results of numerical simulations. Different physical factors controlling particle momentum change and drifts are analysed. It is shown that the vortex electric field is the main cause of trapped particle motion. Drift due to magnetic field inhomogeneity predominantly disturbs free particle gyroorbit along background magnetic field in the vicinity of a vortex. Our simulations show the development of a subdiffusion regime.