Interactions of solar wind plasma with dust grains: Effects of strong plasma anisotropy

In this paper we re-examined the fundamental physics of charging of a dust particle in the moon environment by tenuous anisotropic solar wind plasma. The majority of work on dusty (complex) plasmas is largely concerns with laboratory plasmas, in which charging process of dust grains is very fast, thus making practical the working concept of dynamically equilibrium floating potential and grain charge. However, solar wind plasma parameters are considerably different at the moon orbit, and we found the characteristic charging time of lunar dust grains to be considerably longer, ranging from 3 to 4.6 min for micron size particles, and up to 7.6 h for 10-nm grains, depending on the value of plasma streaming velocity. These findings make it clear that the transient stage of charging process is important in the moon environment, and equilibrium floating potential and grain charge could be considered as long time asymptotic values. For this reason we re-formulated the moon dust charging process as an inherently time-dependent problem and derived the time-dependent charging equation for the grain potential for general case of anisotropic solar wind plasma. Using the results of our kinetics analysis we found that the distribution of charge density over grain surface submerged into solar wind plasma is highly anisotropic, thus making the OML model, which is based on the assumption of isotropic distribution of surface charge density, not applicable to the grain charging problem by the solar wind plasma.     

Observations of magnetohydrodynamic turbulence in the 3D heliosphere

The solar wind is a high Reynolds’ number plasma flow of solar origin that permeates the whole heliosphere. It is also the only accessible medium in which to study collisionless magnetohydrodynamic turbulence performing direct measurements. This represents a topic of fundamental importance to both plasma physics and astrophysics. During the past decades, in situ observations on the ecliptic and at high heliographic latitudes have been very valuable to shed some light on the intricate nature of space plasma turbulence. In this brief review, we will mainly describe the evolution experienced by the turbulence as the solar wind expands into the interplanetary space. We will also address implications due to different processes of local generation of turbulence which might be at work on the ecliptic and at high latitude. Moreover, the fact that solar wind fluctuations are not isotropic and poorly single scale-invariant, two of the fundamental hypotheses at the basis of Kolmogorov’s theory (K41), will give us the possibility to discuss also the relevance of intermittency in the study of space plasma turbulence.     

Charging of dust grains by anisotropic solar wind multi-component plasma

In this paper we study the charging process of small grain particles by anisotropic multi-component solar wind plasmas (electrons, protons and heavy ions), versus two-component (electron/proton) plasmas. We are focusing attention on the important characteristics of the charging process, namely the charging time, floating potential and current content as functions of plasma parameters such as He⁺⁺/H⁺ (α/p) number density and T_α/T_p temperature ratios of alpha particles to protons, as well as plasma streaming velocity v₀. Measured statistical properties of solar wind plasma parameters at 1 AU show considerable variations in α/p-temperature ratios from 1 to 10, in α/p-number density ratio from 0.01 to 0.35, as well as in values of streaming velocity v₀ from 200 km/s to 1000 km/s and more. Periods of these variations could last for several days each, leading to significant variability in the charging process, according to newly derived general analytical expressions. Numerical calculations performed for protons/alphas plasmas showed large disparity in the charging characteristics. For example, in anisotropic plasma, grain charging time varies up to 90% depending on α/p-particles temperature and number density ratios, whereas changes in floating potential are up to 40%. In contrast, in isotropic plasma, charging characteristic for grains do not change very much for the same plasma parameters variations, with charging time varying about 12% and floating potential only varying about 4%. It is also shown that in highly anisotropic plasma, with all ballistic electrons and ions, dust grains could not hold their charges, and characteristic discharged time is calculated. We note that the analysis is equally applicable to any sized body immersed in solar wind plasma.     

Ionospheric losses of Venus in the solar wind

The nature of ionospheric losses from Venus is of essential importance for understanding the ionosphere dynamics of this unmagnetized planet. A plausible mechanism that can explain the escape of charged particles involves the solar wind interaction with the upper atmospheric layers of Venus. The hydrodynamic approach proposed for plasma expansion in the present study comprises two populations of positive ions and the neutralizing electrons, which interact with the solar wind electrons and protons. The fluid equations describing the plasma are solved numerically using a self-similar approach. The behavior of plasma density, velocity, and electric potential, as well as their reliance upon solar wind parameters have been examined. It is found that for noon midnight sites, the oxygen ion-to-electron relative density may be the main factor to enhance the ionic loss. However, the other parameters, like hydrogen density and solar wind density and velocity seem to do not stimulate the runaway ions. For lower dawn-dusk region, the plasma are composed of hydrogen and oxygen ions as well as electrons, but for higher altitudes only hydrogen ions and electrons are encountered. All ionic densities play an important role either to reduce or boost the ionic loss. The streaming solar wind velocity has no effect on the plasma escaping for lower altitudes, but it reduces the expansion at higher altitudes.     

An updated theory of the polar wind

The ‘classical’ polar wind is an ambipolar outflow of thermal plasma from the terrestrial ionosphere at high latitudes. As the plasma escapes along diverging geomagnetic flux tubes, it undergoes four major transitions, including a transition from chemical to diffusion dominance, a transition from subsonic to supersonic flow, a transition from collision-dominated to collisionless regimes, and a transition from a heavy to a light ion. A further complication arises because of horizontal convection of the flux tubes owing to magnetospheric electric fields. Recent modelling predictions indicate that the polar wind has the following characteristics: (1) The ion and electron distributions are anisotropic and asymmetric in the collisionless regime; (2) Elevated electron temperatures ( ∼ 10,000 K) act to produce significant escape fluxes of suprathermal O⁺ ions; (3) The interaction of the hot magnetospheric and cold ionospheric electron populations leads to a localized (double layer) electric field which accelerates the polar wind ions; (4) A time-dependent expansion produces suprathermal ions; and (5) Large perturbations lead to the formation of forward and reverse shocks. These and other results are reviewed.     

Numerical model of time-dependent hydrodynamic flows with mass loading

We study the influence of mass loading on the formation of shocks employing a time-dependent hydrodynamic model. Numerical examples illustrate in which way the formation and propagation of shocks in the plasma is affected by variations of the ionization rate (due to changes of solar wind parameters or of neutral gas production). It is found that mass loading may be an efficient mechanism to prevent the formation of discontinuities in hydrodynamic flows.     

Study of the solar Slow Sonic,Alfvén and Fast Magnetosonic transition surfaces

In this paper we study the shape, extend and time variations of the solar wind transition surfaces using the Lima and Priest (1993) hydrodynamic model adequately adapted for the case of the solar wind flow. The transition surfaces, namely the Slow (Sonic), the Alfvén, and the Fast Magnetosonic surface, are important boundaries around the Sun and play a crucial role in the development of the solar wind and the structure of the inner heliosphere. We determine the shape and dimension of these surfaces as a function of heliographic latitude using measurements from Ulysses spacecraft, and we also study their temporal variation using data from spacecrafts at 1 AU (OMNI database). Furthermore, we establish their dependence with the solar activity, demonstrating their shape and location for the last two solar cycles. From this we noticed that the temporal variation of all transition surfaces follows the 11-year solar cycle. Finally, from the OMNI database, we have studied the temporal variation over the past 40 years of the plasma β parameter, the kinetic to magnetic and the kinetic to thermal energy ratios, at a distance of 1 AU from the Sun.     

太阳风中航天器带电与尾迹效应的模拟

航天器充电和尾迹效应会对周围等离子体造成扰动,影响测量装置结果的准确性.利用SPIS （Spacecraft Plasma Interaction Software）分别模拟了航天器与太阳风的相互作用,考察了光电效应以及航天器尺度对表面充电情况和尾迹效应的影响.结果表明:太阳风环境下,等离子体密度稀薄,电子电流比光电子电流小得多,航天器表面为正电势,航天器后部有清晰的尾迹结构,尾迹带负电;光电效应可改变尾迹结构,与无光电效应相比,光电效应使得航天器尾迹尺度变大;由于太阳风定向运动动能大于航天器表面势能,航天器的尾迹结构与其几何尺寸有关,航天器尺寸越大,尾迹尺度越大.      

Dependence of ESP intensity on collisionless shock wave characteristics

Flux variations of 1 – 5 MeV protons are studied in energetic storm particle events with respect to the preshock solar wind plasma parameters and to the thickness of the collisionless interplanetary shock wave. It is found that the peak intensity in ESP events depends on pre-shock plasma density and on the thickness of the transition region. These relations predict, in agreement with recent observations, the increase of ESP events at larger heliocentric distances.     

The Compton–Getting effect of energetic particles with an anisotropic pitch-angle distribution

This paper presents a simulation of anisotropy measurements by the low-energy charged particle (LECP) experiment on Voyager 1 for cases when the particle pitch-angle distribution function in the solar wind plasma reference frame is not isotropic. The model includes both the Compton–Getting anisotropy and perpendicular diffusion anisotropy that possibly exists in the upstream region of the termination shock. The results show that the Voyager 1 data cannot rule out either the model with zero solar wind speed or the one with a finite speed on qualitative basis. The determination of solar wind speed using the Compton–Getting effect is affected by the assumption of the magnetic field direction and perpendicular diffusion anisotropy. Because the pitch-angle distribution anisotropy is so large, a small uncertainty in the magnetic field direction can produce very different solar wind speeds ranging from ∼0 to >400 km/s. In fact, if the magnetic field is chosen to be in the Parker spiral direction, which is consistent with the magnetometer measurement on Voyager 1, the derived solar wind speed is still close to the supersonic value. Only the two lowest-energy channels of the LECP instrument may give a definitive answer to the solar wind speed. However, because these channels contain a very high level of cosmic ray background, an uncertainty of just a few percent in the background can entirely hamper the estimate of solar wind speed.     

膨胀磁场和太阳风粒子相互作用三维数值模拟

采用三维模型,使用混合网格质点法HPIC(Hybrid Particle-in-Cell)对膨胀的磁场和太阳风相互作用过程进行数值模拟.研究了线圈产生的偶极子磁场在注入等离子体后和太阳风粒子的相互作用过程,并对以不同速度入射的等离子体引起的太阳风粒子的变化和磁场变化进行了比较.研究结果表明,偶极子磁场和太阳风作用时会产生弓形激波,此时磁压等于太阳风粒子的动压,当向线圈产生的偶极子磁场中注入高能等离子体时引起磁场膨胀,膨胀的磁场将会排斥太阳风粒子向外运动,从而引起弓形激波的变化,增大与太阳风相互作用的面积,并且粒子入射速度越大,磁场膨胀越明显,与太阳风相互作用愈强.      

Bow shock: Power aspects

It is clear that the primary energy source for magnetospheric processes is the solar wind, but the process of energy transfer from the solar wind into the magnetosphere, or rather, to convecting magnetospheric plasma, appears to be rather complicated. Bow shock is a powerful transformer of the solar wind kinetic energy into the gas dynamic and electromagnetic energy. A jump of the magnetic field tangential component at front crossing means that the front carries an electric current. The solar wind kinetic energy partly transforms to gas kinetic and electromagnetic energy during its passage through the bow shock front. The transition layer (magnetosheath) can use part of this energy for accelerating of plasma, but can conversely spend part its kinetic energy on the electric power generation, which afterwards may be used by the magnetosphere. Thereby, transition layer can be both consumer (sink) and generator (source) of electric power depending upon special conditions. The direction of the current behind the bow shock front depends on the sign of the IMF Bz-component. It is this electric current which sets convection of plasma in motion.     

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.     

月表磁异常影响的太阳风电子反射运动

模拟太阳风电子向月表运动的轨迹, 研究由于月表磁异常的存在造成的电子反射运动. 首先设定行星际磁场Bsw 指向月球并与月表垂直, 将月表的磁异常区看成是一个磁偶极子, 偶极矩大小为Mcb; 然后分别考察该偶极矩与行星际磁场方向平行, 反平行以及±90° 的情形, 通过计算发现, 被反射的电子数目会随着磁偶极矩和行星际磁场的方向改变而改变. 在偶极矩与行星际磁场平行的情况下, 反射率最大; 随着夹角的增大, 反射率减小. 这些结果为利用电子反射法高精度遥测月表磁场提供了很重要的信息.      

The plasma sheet and boundary layers under northward IMF: A multi-point and multi-instrument perspective

During conditions of northward interplanetary magnetic field (IMF), the near-tail plasma sheet is known to become denser and cooler, and is described as the cold-dense plasma sheet (CDPS). While its source is likely the solar wind, the prominent penetration mechanisms are less clear. The two main candidates are solar wind direct capture via double high-latitude reconnection on the dayside and Kelvin–Helmholtz/diffusive processes at the flank magnetopause. This paper presents a case study on the formation of the CDPS utilizing a wide variety of space- and ground-based observations, but primarily from the Double Star and Polar spacecraft on December 5th, 2004. The pertinent observations can be summarized as follows: TC-1 observes quasi-periodic (∼2 min period) cold-dense boundary layer (compared to a hot-tenuous plasma sheet) signatures interspersed with magnetosheath plasma at the dusk flank magnetopause near the dawn-dusk terminator. Analysis of this region suggests the boundary to be Kelvin–Helmholtz unstable and that plasma transport is ongoing across the boundary. At the same time, IMAGE spacecraft and ground based SuperDARN measurements provide evidence of high-latitude reconnection in both hemispheres. The Polar spacecraft, located in the southern hemisphere afternoon sector, sunward of TC-1, observes a persistent boundary layer with no obvious signature of boundary waves. The plasma is of a similar appearance to that observed by TC-1 inside the boundary layer further down the dusk flank, and by TC-2 in the near-Earth magnetotail. We present comparisons of electron phase space distributions between the spacecraft. Although the dayside boundary layer at Polar is most likely formed via double high-altitude reconnection, and is somewhat comparable to the flank boundary layer at Double Star, some differences argue in favour of additional transport that augment solar wind plasma entry into the tail regions.     

Interplanetary coronal mass ejections during the descending cycle 23: Sheath and ejecta properties comparison

We have used Omniweb data in order to identify the sheath and the ejecta boundaries of 67 shock-driving interplanetary coronal mass ejections during the time period 2003–2006. We examine and compare their statistical properties (speed, magnetic field strength, proton density and temperature, proton plasma beta), with those of the typical solar wind. We also calculate their passage time and radial width. We study the correlation between the ejecta and sheath characteristics.     

太阳风电子在月表电磁场中的运动

月表磁异常区的分布是月球探测工程的重要内容. 但是由于月表电磁环境错综复 杂, 通常认为月球表面在特殊的空间天气条件下会带有数千伏电压. 以往的空间研究已经证实, 表面的带电与放电容易造成卫星仪器的异常或失联. 月表电场对电子 反射法有重要影响, 研究分析不同电磁条件下太阳风电子的运动轨迹,对月表环境 (电磁环境, 太阳风条件, 等离子体参数等)的研究可以更加深入细致. 通过模拟向月表运动的太阳风电子的运动轨迹, 分析了月表电磁环境的改变对太阳风电子反射 的影响, 并着重研究了月表电场对电子反射法遥感探测月表磁异常的影响, 为探测 月表电磁环境提供了重要的信息.      

Formation and characteristics of low latitude boundary layer

Characteristics of low latitude boundary layer (LLBL) of the Earth’s magnetosphere are investigated using data of Interball/Tail probe observations. The role of different processes of LLBL formation is discussed. The high level of magnetosheath turbulence is taken into account. It is shown that the turbulent nature of magnetic field and plasma fluctuations in the magnetosheath is one of the main factors determining the structure of LLBL. The results of Interball/Tail probe observations of the event 9 March 1996 are analyzed. The thickness of LLBL is determined for the number of cases. The change of LLBL thickness under the influence of the changes of solar wind parameters is investigated. It is shown that variability of solar wind conditions can be the important factor controlling LLBL thickness. Results of observations are compared with the theory which explains the value of LLBL thickness as the result of plasma transport inside the magnetosphere. It is shown that the theory gives the qualitative explanation of the observed dependence of LLBL thickness on solar wind parameters.     

Gas-dynamic approach to the theory of non-linear ion-acoustic waves in plasma with Kaniadakis’ distributed species

The non-linear theory of ion-acoustic waves in two-spesies isotropic collisionless plasma is developed. Both light electron and heavy ion species in the plasma are distributed with Kaniadakis’ statistics. Kaniadakis’ gas law is derived. The exact formula for the Sagdeev pseudopotential in parametric form is derived by the method based on the integration of the inverse function. The pseudopotential is analyzed. It is shown that periodic ion-acoustic waves and solitons are possible in the studied plasma. The differential equation describing the profile of the ion concentration in the wave is derived. The profiles of this concentration in the periodic ion-acoustic wave and soliton are calculated.     

Remote upgrading of a space-borne instrument

The main purposes of experiment “Obstanovka” (“Environment” in Russian) consisting of several instruments are to measure a set of electromagnetic and plasma phenomena characterizing the space weather conditions, and to evaluate how such a big and highly energy consuming body as the International Space Station disturbs the surrounding plasma, and how the station itself is charged due to the operation of so many instruments, solar batteries, life supporting devices, etc. Two identical Langmuir electrostatic probes are included in the experiment “Obstanovka”. In this paper the Langmuir probes for “Obstanovka” experiment are described, including the choice of geometry (spherical or cylindrical), a more reliable method for the sweep voltage generation, an adaptive algorithm for the probe’s operation. Special attention is paid to the possibility for remote upgrading of the instrument from the ground using the standard communication channels.