Effect of ion dynamics on the evolution of electron phase-space holes

The evolution of two-dimensional (2D) electron phase-space holes (electron holes) has been previously investigated with electrostatic Particle-in-Cell (PIC) simulations, which neglect ion dynamics. The electron holes are found to be unstable to the transverse instability, and their evolution is determined by the combined action between the transverse instability and the stabilization by the background magnetic field. In this paper, the effect of ion dynamics on the evolution of an electron hole is studied. In weakly magnetized plasma (Ω_e < ω_pe, where Ω_e and ω_pe are electron gyrofrequency and plasma frequency, respectively), the electron hole is still unstable to the transverse instability. However, it evolves a little faster and is destroyed in a shorter time when ion dynamics is considered. In strongly magnetized plasma (Ω_e > ω_pe), the electron hole is broken due to the lower hybrid waves, and its evolution is much faster.      

From solar nanoflares to stellar giant flares: Scaling laws and non-implications for coronal heating

In this study we explore physical scaling laws applied to solar nanoflares, microflares, and large flares, as well as to stellar giant flares. Solar flare phenomena exhibit a fractal volume scaling, V(L)  L^1.9, with L being the flare loop length scale, which explains the observed correlation between the total emission measure EM_p and flare peak temperature T_p in both solar and stellar flares. However, the detected stellar flares have higher emission measures EM_p than solar flares at the same flare peak temperature T_p, which can be explained by a higher electron density that is caused by shorter heating scale height ratios s_H/L ≈ 0.04–0.1. Using these scaling laws we calculate the total radiated flare energies E_X and thermal flare energies E_T and find that the total counts C are a good proxy for both parameters. Comparing the energies of solar and stellar flares we find that even the smallest observed stellar flares exceed the largest solar flares, and thus their observed frequency distributions are hypothetically affected by an upper cutoff caused by the maximum active region size limit. The powerlaw slopes fitted near the upper cutoff can then not reliably be extrapolated to the microflare regime to evaluate their contribution to coronal heating.     

Storm-time formation of a relativistic electron belt and some relevant phenomena in other magnetospheric plasma domains

An empirical formula relating the strength of a storm given by its |Dst|_max with the L-coordinate of the peak of storm-injected relativistic electrons is one of a few well-confirmed quantitative relations found in the magneto-spheric physics. We successively extended a dataset of the formula’s basic storms with several events of high Dst-amplitude up to the highest observed |Dst|_max = 600 nT. Possible applying of the formula to the predicting of the ring-current plasma-pressure distribution and the lowest westward electrojet position for a storm are discussed. We have also analyzed the 2000–2001 years’ data on relativistic electrons from our instruments installed on EXPRESS-A (geosynchronous orbit; E_e = 0.8–6 MeV), Molniya-3 (h = 500 × 40 000 km, i = 63°; E_e = 0.8–5.5 MeV) and GLONASS (h = 20 000 km, i = 64°; E_e  l MeV) along with other correlated measurements: GOES series (E_e > 2 MeV), geomagnetic indices (Dst, AE, AL) and interplanetary parameters (solar wind, IMF). The goal is to investigate which outer conditions are most responsible for the high/low output of the storm-injected relativistic electrons. For the geosynchronous orbit, two factors are found as the necessary condition of the highest electron output: high and long-lasting substorm activity on a storm recovery phase and high velocity of solar wind. On the contrary, extremely low substorm activity surely observed during whole the storm recovery phase constitutes a sufficient condition of the non-increased after-storm electron intensity. For the first time found cases of the increased after-storm electron intensity observed at the inner L-shells with no simultaneously seen increase in the geosynchronous distances are presented.     

Electrostatic Structure of the Electron Phase-space Holes Generated by the Electron Two-stream Instability with a Finite Width

Space satellite observations in an electron phase-space hole (electron hole) have shown that bipolar structures are discovered at the parallel cut of parallel electric field, while unipolar structures spring from the parallel cut of perpendicular electric field. Particle-in-cell (PIC) simulations have demonstrated that the electron bi-stream instability induces several electron holes during its nonlinear evolution. However, how the unipolar structure of the parallel cut of the perpendicular electric field formed in these electron holes is still an unsolved problem, especially in a strongly magnetized plasma (Ω_e > ω_pe, where Ω_e is defined as electron gyrofrequency and ω_pe is defined as plasma frequency, respectively). In this paper, with two-dimensional (2D) electrostatic PIC simulations, the evolution of the electron two-stream instability with a finite width in strongly magnetized plasma is investigated. Initially, those conditions lead to monochromatic electrostatic waves, and these waves coalesce with each other during their nonlinear evolution. At last, a solitary electrostatic structure is formed. In such an electron hole, a bipolar structure is formed in the parallel cut of parallel electric field, while a unipolar structure presents in the parallel cut of perpendicular electric field.      

High accuracy satellite drag model (HASDM)

The dominant error source in force models used to predict low-perigee satellite trajectories is atmospheric drag. Errors in operational thermospheric density models cause significant errors in predicted satellite positions, since these models do not account for dynamic changes in atmospheric drag for orbit predictions. The Air Force Space Battlelab’s High Accuracy Satellite Drag Model (HASDM) estimates and predicts (out three days) a dynamically varying global density field. HASDM includes the Dynamic Calibration Atmosphere (DCA) algorithm that solves for the phases and amplitudes of the diurnal and semidiurnal variations of thermospheric density near real-time from the observed drag effects on a set of Low Earth Orbit (LEO) calibration satellites. The density correction is expressed as a function of latitude, local solar time and altitude. In HASDM, a time series prediction filter relates the extreme ultraviolet (EUV) energy index E_10.7 and the geomagnetic storm index a_p, to the DCA density correction parameters. The E_10.7 index is generated by the SOLAR2000 model, the first full spectrum model of solar irradiance. The estimated and predicted density fields will be used operationally to significantly improve the accuracy of predicted trajectories for all low-perigee satellites.     

The collective gyration of a heavy ion cloud in a magnetized plasma

In both the ionospheric barium injection experiments CRIT I and CRIT II, a long-duration oscillation was seen with a frequency close to the gyro frequency of barium and a time duration of about one second. A model for the phenomenon which was proposed for the CRIT 1 experiment is here compared to the results from CRIT II which made a much more complete set of measurements. The model follows the motion of a low-β ion cloud through a larger ambient plasma. The internal field of the model is close to antiparallel to the injection direction v_i but slightly tilted towards the self-polarization direction E_P = −v_ixB. As the ions move across the magnetic field, the space charge is continuously neutralized by magnetic-field aligned electron currents from the ambient ionosphere, drawn by the divergence in the perpendicular electric field. These currents give a perturbation of the magnetic field related to the electric field perturbation by ΔE/ΔB ≈ V_a. The model predictions agree quite well with the observed vector directions, field strengths, and decay times of the electric and magnetic fields in CRIT II. The possibility to extend the model to the active region, where the ions are produced in this type of self-ionizing injection experiments, is discussed.     

电离层Es离子分布数值建模

从离子连续性方程和动量方程出发,比较全面地考虑了太阳光电离、星光电离、地冕电离、星际背景电离和流星离子流等电离源,综合分析Es层的主要动力学过程和光化学过程,以风剪切理论为基础,研究中性成分碰撞、电场作用、金属离子作用和分子离子作用等机制对离子分布的影响,建立了一维时变电离层Es物理模型.对离子产生率、垂直方向的离子速度在高度和时间上的变化进行仿真和计算,得到了在电场及风剪切作用下的离子密度剖面24h内的变化规律.根据建立的模型,以昆明相干散射雷达观测的径向风结果作为模型输入参考,仿真并得到了电离层Es电子密度的时空分布,据此反演出电离层Es临界频率f₀E_s.与昆明站电离层测高仪同时间观测的情况进行对比,结果比较一致,验证了建立的Es层物理模型正确性.      

Diagnosing the storm-time ring current and other magnetospheric plasma domains using relativistic electron data

Ring current ions and relativistic electrons simultaneously measured on board MOLNIYA-1 are analyzed in comparison with the ground-based magnetometer data for the period of a strong magnetic storm (|D_st|max≈230 nT). Injection of >500 keV electrons into the slot region (L≈3) near equatorial plane is occurred on time scale ≈1 hour, when, during the magnetic storm maximum, the extreme low-latitude position of auroral electrojets is reached and ring current becomes more symmetrical. Positions of both the ring current maximum and electron intensity maximum (L_max) are consistent to our previous result: |D_st|_max = 2.75 • 10⁴/L⁴_max. An extreme storm-time low-latitude position of the west electrojet center (for amplitudes of |D_st|_max up to 600 nT) is shown to be in a good consistence with this empirical dependence. It is supposed the trapped radiation boundary collapses down to L≈L_max in the course of the storm main phase.     

Energy deposition in the upper atmosphere in the EXCEDE III experiment

The EXCEDE III sounding rocket flight of April 27, 1990 used a 18 Ampere 2.5 keV electron beam to produce an artificial aurora in the region 90 to 115 km. A “daughter” sensor payload remotely monitored the low-energy X-ray spectrum while scanning photometers measured the spatial profile of prompt emissions of N₂⁺ (1N) and N₂ (2P) transitions (3914Å and 3805Å, respectively). Two Ebert-Fastie spectrometers measured the spectral region from 1800 to 8000Å. On the “mother” accelerator payload, the return current electron differential energy spectra were monitored by an electrostatic analyzer (up to 10 keV) and by a retarding potential analyzer (0 eV to 100 eV). We present an overview of the results from this experiment.     

稳态等离子体推进器羽流场数值模拟

采用二维轴对称模型,使用粒子网格法(PIC)和直接模拟蒙特卡洛法(DSMC)相结合的方法,对稳态等离子体推进器(SPT)羽流场进行了数值模拟.采用DSMC方法中的随机取样频率法(RSF)求解粒子碰撞过程,并对比了不同的分配电荷方式、电子运动模型及SPT出口条件时的羽流场.将不同条件计算得到的羽流场中距SPT出口0.1?m,0.5?m及1.0?m处的离子电流密度和电荷密度与实验结果进行了对比,得出在采用面积权重法分配电荷、等熵模型描述电子运动和用实验值设定发动机出口参数时对SPT羽流场数值模拟的电流密度和轴向附近的电荷密度结果与实验结果符合程度较好的结论.     

Substorms in the inner plasma sheet

Thin Current Sheets (TCS) are regularly formed prior to substorm breakup, even in the near-Earth plasma sheet, as close as the geostationary orbit. A self-consistent kinetic theory describing the response of the plasma sheet to an electromagnetic perturbation is given. This perturbation corresponds to an external forcing, for instance caused by the solar wind (not an internal instability). The equilibrium of the configuration of this TCS in the presence of a time varying perturbation is shown to produce a strong parallel thermal anisotropy (T T) of energetic electrons and ions (E>50keV) as well as an enhanced diamagnetic current carried by low energy ions (E<50keV). Both currents tend to enhance the confinement of this current sheet near the magnetic equator. These results are compared with data gathered by GEOS-2 at the geostationary orbit, where the magnetic signatures of TCS, and parallel anisotropics are regularly observed prior to breakup. By ensuring quasi-neutrality everywhere we find, when low frequency electromagnetic perturbations are applied, that although the magnetic field line remains an equipotential to the lowest order in T_e/T_i, a field-aligned potential drop exists to the next order in (T_e/T_i). Thus the development of a TCS implies the formation of a field-aligned potential drop ( few hundred volts) to ensure the quasi-neutrality everywhere. For an earthward directed pressure gradient, a field-aligned electric field, directed towards the ionosphere, is obtained, on the western edge of the perturbation (i.e. western edge of the current sheet). Thus field aligned beams of electrons are expected to flow towards the equatorial region on the western edge of the current sheet. We study the stability of these electron beams and show that they are unstable to “High Frequency” (HF) waves. These “HF” waves are regularly observed at frequencies of the order of the proton gyrofrequency (f_H+) just before, or at breakup. The amplitude of these HF waves is so large that they can produce a strong pitch-angle diffusion of energetic ions and a spatial diffusion that leads to a reduction of the diamagnetic current. The signature of a fast ion diffusion is indeed regularly observed during the early breakup; it coincides with the sudden development of large amplitude transient fluctuations, ballooning modes, observed at much lower frequencies (fH+). These results suggest that the HF waves, generated by field-aligned electron beams, provide the dissipation which is necessary to destabilize low frequency (ballooning) modes.     

Specific features of electron distributions at altitudes of 400 km

New experimental data obtained on the orbital station ‘MIR’ in 1991 during solar maximum are discussed. Electron fluxes with E_e>75 keV were registered for three different directions as well as for electrons with E_e>300 and 600 keV. Spatial and time distributions of electron fluxes in the trapping region are presented. In the inner radiation belt an additional maximum is observed at L=1.25–1.35, and the fluxes in the 22-05^h MLT interval are 2–3 orders of magnitude smaller, than during other local times. In this region a flattening of the electron spectrum is observed. The results obtained were compared with the AE-8 model.     

The effect of the longitudinal propagating waves on the electron acceleration in the reconnecting current sheet

Based upon the most efficient electron acceleration near the midplane of 3D non-neutral driven reconnecting current sheet (RCS) and the electrostatic wave excitation by the drift Maxwellian distribution of electrons in Vlasov simulation, we assume that the electrostatic waves mainly propagate opposite to the reconnecting electric field and investigated how these waves affect the electron acceleration. The main results are: (1) when the electron’s velocity equals to the phase speed of the waves, they will be trapped and have the different accelerating characteristics from the untrapped electrons through solving the momentum equations of electrons analytically; (2) the test particle simulations further prove that the number of the energetic electrons decreases with the increasing intensity of unstable waves, and the distribution of the energetic electrons takes on the double power-law.     

Numerical simulation of three-dimensional reconnection due to the instability of collisionless current sheets

Based on analytical calculations we have currently argued that spontaneous reconnection through thin collisionless current sheets is an essentially three-dimensional (3 D) process (Büchner, 1996 a, b). Since 3 D kinetic PIC codes have become available, the three dimensional nature of the collisionless current sheet decay are now illustrated by numerical simulations (Büchner and Kuska, 1996; Pritchett and Coroniti, 1996; Zhu and Winglee, 1996). While the latter two claim a coupling to a longer wavelength kink mode as a main factor, destabilizing thin current sheets in 3 D, our simulations have revealed that even shorter scale perturbations in the current direction suffice to destabilize thin sheets very quickly. Since past simulation runs, however, were limited to mass ratios near unity, the influence of the electrons was not treated adequately. We have now investigated the stability of thin collisionless current sheets including 64 times lighter negatively charged particles. We can now show that while the two-dimensional tearing instability slows down for M = M_p/m_e = 64, the three-dimensional current sheet decay is a much faster process — practically as fast as the mass ratio M = 1 3 D sheet decay, even without kinking the sheet. We further conclude that, unlike the two-dimensional tearing instability, the three-dimensional decay of thin current sheets is not controlled by the electrons. For a sheet width comparable with the ion inertial length, we also recovered signatures of the Hall effect as predicted by Vasyliunas (1975) in the mass ratio M = 64 case. The ion inertial length seems to be the critical scale at which the sheet starts to decay.     

释放SF6气体扰动电离层时的人工气辉研究

在电离层高度释放SF₆气体能够显著扰动电离层.根据SF₆分子在电离层中的扩散方程,同时考虑其在电离层中主要的离子化学反应,研究了SF₆气体释放后电离层各粒子浓度的时空变化,计算了产生人工气辉的体发射系数和发射强度.结果表明,SF₆气体在电离层高度释放后,电子和O⁺的密度均有大幅度下降,主要的负离子成分由电子转变成SF₅^-;在释放过程中,主要产生777.4 nm和135.6 nm两种气辉,且前者的气辉强度远小于后者;电离层温度对气辉的强度有很大的影响.本文的数值计算与美国IMS/SF₆实验观测数据进行比较,结果近似,且通过数据比较还能准确推断出实验时当地的电离层温度.     

Auroral signatures of transient processes in the outer magnetosphere

We present an analysis of sporadic and recurrent injections of magnetospheric ions in the midnight auroral oval during substorms and of the associated ionospheric ion outflows. The source of plasma sheet precipitating ions is determined using a simple method, based on the measured relation between the ion inverse velocity and time (l = v × t). This method is applied here to two typical passes of the Interball-Auroral (IA) satellite at distances of 3 R_E above the auroral regions. Substorm related ion injections are shown to be mainly due to time of flight effects. In contrast with particle trajectory computations (Sauvaud et al., 1999), the inverse velocity method does not require magnetic and electric field models and can thus be used systematically for the detection of time of flight dispersed ion structures (TDIS). This allowed us to build a large database of TDIS events and to perform a statistical analysis of their spatial distribution. For the cases presented here the source region of the injected ions is found at radial distances from 18 to 30 R_E near the equatorial magnetosphere. At Interball altitudes ( 3 R_E), ion injections detected at the poleward boundary of the nighside auroral oval are associated with shear Alfvén waves superimposed over large-scale quasi-static current structures. We show that the most poleward TDIS are collocated with a large outflow of ionospheric H⁺ and O⁺ displaying pitch-angle distributions peaked in the pitch-angle range 90°–120°. These ions are thus accelerated perpendicularly to the magnetic field not only in the main auroral acceleration region but also up to at least 3 R_E. The expanding auroral bulge thus constitutes a significant source of H⁺ and O⁺ ions for the mid-tail magnetosphere.     

Solitary and periodic waves in magneto–rotating plasmas: Effect of the Kappa–Cairns distributed electrons

The nonlinear propagation of ion–acoustic (IA) waves in a magneto–rotating plasma is studied by considering the Kappa-Cairns electron distribution. Employing the perturbation scheme, Korteweg–de Vries equation is derived. It is seen that both positive and negative potential solitons can be supported in the present plasma model. The numerical results reveal that the Kappa-Cairns distributed electrons modify features of the electrostatic waves significantly. The effects of non–thermal parameters, plasma rotation frequency, ion temperature, and the wave propagation angle on electrostatic solitary wave structures are also discussed here. It is found that the plasma parameters considerably influence the propagation of IA waves in rotating plasmas. Furthermore, using the bifurcation theory of planar dynamical systems to the K-dV equation, we have presented the existence of solitary and periodic traveling waves. Our study may be helpful to understand the behavior of ion–acoustic wave in the rotating plasma.     

基于拉曼努金和的非均匀多载波调制系统

针对高速移动通信,多普勒效应和多径效应导致信道非均匀,提出一种基于拉曼努金和的频谱可调的非均匀调制多载波系统.首先,根据拉曼努金和的正交性和周期性,推导了拉曼努金傅里叶正反变换的完全重建条件,进而建立基于拉曼努金和的多载波调制系统（RFMT,Ramanujan Fourier Multi-Tone system）.由于拉曼努金和的非均匀频谱分布性质及频率共振性质,RFMT在不同载波通道的误比特率不同,可实现对数据的不均等保护.在非均匀信道下RFMT具有优于OFDM的抗多径能力,仿真验证了RFMT的抗多径有效性.在使用迫零均衡算法时,10-5误比特率下,根据信道情况设计的非均匀多载波系统RFMT对E_b/N₀的要求可以比OFDM低4 dB.      

混合球磨对YAl_2颗粒形貌及表面改性的影响

采用扫描电镜、能谱分析、透射电镜等分析方法研究了球磨对YAl₂及YAl₂/Mg颗粒的组织形貌、YAl₂颗粒在YAl₂/Mg颗粒中的分散性及YAl₂/Mg界面结合状况的影响.结果表明:YAl₂颗粒经球磨20 h后,其平均粒径约为1.5 μm,颗粒边缘尖角数量明显减少,呈圆滑的近球形.加入Mg颗粒进行混合球磨20 h后,YAl₂颗粒在混合颗粒中呈弥散分布,且表面包裹一层单质Mg,YAl₂/Mg的界面为直接结合型,未发现有界面反应.