^{1}

^{1}

^{1}

^{*}

The time evolution of the magnetic field aligned super-small-scale irregularities in the concentration of charged particles, existing in the near-Earth rarefied plasma, is studied with the help of the model simulation. A new version of the two-dimensional mathematical model, developed earlier in the Polar Geophysical Institute, is utilized to investigate the temporal history of the irregularity with circular cross section, created initially in the near-Earth plasma. The utilized model is based on a numerical solution of the Vlasov-Poisson system of equations, with the Vlasov equations describing the distribution functions of charged particles and the Poisson equation governing the self-consistent electric field. The results of simulation indicate that the mobility of the positive ions ought to influence essentially on the time evolution of the super-small-scale irregularities in the concentration of charged particles, existing in the near-Earth rarefied plasma.

Irregularities in the concentration of charged particles are often observed in the near-Earth plasma. These irregularities are a natural phenomenon in the Earth’s ionosphere and magnetosphere. The irregularities have a wide range of spatial scales, ranging from thousands of kilometers to a few Debye lengths. The charged particles density increases and depletions inside irregularities can lie in the range from some tens of percentages to a few portions. In the ionosphere, the well-known equatorial anomaly at F-layer altitudes is the example of large-scale irregularities [

Not large-scale irregularities are predominately magnetic field aligned. Usually, there are three generic types of their structures: rods, wings, and sheets. Rods are isotropic in the plane perpendicular to the geomagnetic field. Wings and sheets are elongated not only along geomagnetic field but also in the perpendicular plane along a certain direction [

It is known that the irregularities either are naturally present or may be artificially produced as a result of active experiments in the near-Earth plasma. For example, both a large-scale irregularity and super-small-scale irregularities in the concentration of charged particles may be formed by high-power high-frequency radio waves, pumped into the ionosphere by ground based ionospheric heaters. To investigate the behavior of the artificially created irregularities in the near-Earth plasma not only the experimental and theoretical but also computational studies may be applied. The large-scale F-layer modification by powerful high frequency waves was investigated with the help of mathematical models in some studies (for example, see [

The purpose of the present paper is to investigate numerically the time evolution of the magnetic field aligned super-small-scale irregularities in the concentration of charged particles, existing in the near-Earth rarefied plasma, whose cross section is much less than the mean free path of particles between successive collisions and commensurable with a Debye length.

The near-Earth plasma at altitudes of the ionospheric F layer and magnetosphere is a rarefied compound consisting of electrons and positive ions in the presence of a strong, external, magnetic field. The parameters of the magnetic field aligned irregularities, studied in the present paper, have been described in details by Wong et al. [

The studied irregularities are geomagnetic field-aligned. In the vicinity of the irregularity, gradients of the plasma parameters in the longitudinal direction are supposed to be much less than those in a plane perpendicular to a magnetic field. Therefore, plasma parameters inside and beyond the irregularity may be considered as independent on the longitudinal coordinate. Hence, it is sufficient to consider a two dimensional flow of plasma in a plane perpendicular to a magnetic field line.

In the present paper, we study the magnetic field aligned irregularities, having circular cross sections, that is, the rodlike irregularities, with their diametrical sizes being much less than the mean free path of particles between successive collisions. To investigate the time evolution of the studied irregularities a two-dimensional mathematical model, developed earlier in the Polar Geophysical Institute, is utilized. In this model, kinetic processes in the plasma are simulated by using the Vlasov-Poisson system of equations, with the Vlasov equations describing the distribution functions of charged particles and the Poisson equation governing the self-con- sistent electric field. The Vlasov equations are numerically solved applying a macroparticle method. The Poisson equation is solved using a finite-difference method. The utilized system of equations has been earlier considered, for example, in the studies by Hockney and Eastwood [

where

In the latter studies, the mathematical model has been utilized for numerical simulation of the behavior of super-small-scale rodlike and sheetlike irregularities existing in the near-Earth rarefied plasma. The time evolutions of plasma irregularities, having initial cross-section dimension commensurable with a Debye length, have been simulated during the period sufficient for the irregularities to decay completely. In these studies, a non-sta- tionary process started from a situation when, at the initial moment, the spatial distribution of the electron concentration contained an irregularity at the center of the simulation region, with the spatial distribution of the positive ion concentration having been homogeneous in all simulation region. Thus, an electrical neutrality of the plasma was disturbed at the initial moment.

Unlike, in the present study, a non-stationary process starts from a situation when, at the initial moment, the spatial distribution of the electron concentration coincides with the spatial distribution of the positive ion concentration in all simulation region. Thus, the process starts from the completely electrically neutral state in all simulation region. However, the spatial distributions of the electron and positive ion concentrations may be inhomogeneous and can contain an irregularity at the center of the simulation region at the initial moment.

Moreover, in the present study, a new version of the mathematical model is applied in which positive ions can move. In the mathematical model, utilized in the studies of Mingalev et al. [

Also, in the new version of the mathematical model, some parameters of the numerical method were improved. In particular, the quantity of the grid cells was enlarged from 768 × 768 to 1024 × 1024, the average number of macro-particles in the Debye cell for the model plasma was enlarged from 2^{13} to 2^{15}. A two-dimen- sional simulation region, lain in the plane perpendicular to the magnetic field line, is a square and its side length became equal to 128 Debye lengths of the plasma (instead of 96 Debye lengths of the plasma in the previous version of the mathematical model).

The utilized mathematical model can describe the time evolution of the magnetic field aligned super-small-scale irregularities in the concentration of charged particles, having various parameters at the initial moment. The results of calculations to be presented in this paper were obtained for the following initial state of the irregularity. The spatial distribution of the electron concentration (equal to the positive ion concentration at the initial moment), contains a circular irregularity at the center of the simulation region. Beyond this irregularity, the electron concentration is homogeneous and coincides with the non-disturbed electron concentration, n_{0}, which must be given in the simulation region at the initial moment. The initially created irregularity has the cross-section diameter of

The applied mathematical model can describe the behavior of the near-Earth plasma under various conditions. The results of calculations to be presented in this paper were obtained using the input parameters of the model typical for the nocturnal ionospheric plasma at the altitude of 300 km . At this altitude, the plasma contains only positive ions, with the oxygen ion, O^{+}, being the bulk of ion content (99%). At the altitude of 300 km , the value of the non-disturbed electron concentration (equal to the positive ion concentration), n_{0}, is equal to 10^{11} m^{−3}. The electron and ion temperatures are supposed to be equal to 1213 K and 930 K, respectively. The bulk flow velocities of electrons and positive ions are assumed to be zero. The value of the magnetic field, B_{0}, is 4.4 ´ 10^{−}^{5} T.

The above pointed out values yield the following quantities of some physically significant parameters. The electron thermal velocity, ^{7} s^{−1}. The Debye length of the plasma, ^{−}^{3} m, the electron gyro radius, ^{−}^{3} m. The equilibrium period of Langmuir oscillations of electrons, ^{−7} s. The period of cyclotron oscillations of electrons, ^{−7} s, that is, approximately a factor of 2.3 larger than the equilibrium period of Langmuir oscillations of electrons (

We have calculated the time evolution of the distribution functions of charged particles as well as self-con- sisting electric field, taking the input parameters of the mathematical model typical for the nocturnal ionosphere at the level of 300 km , for two distinct cases. The previous version of the mathematical model was utilized for the first case, whereas, the new version of the mathematical model was used for the second case.

Simulation results, obtained for the first case when positive ions were supposed to be immovable, indicate that, after initial moment, the spatial distribution of the electron concentration changes essentially while the positive ion concentration is retained practically invariable. As a consequence, after initial moment, the electrical neutrality of the plasma was broken and an electric charge arises at the center of the simulation region. The configuration of the disturbed area in the spatial distribution of the electric charge density was similar to the configuration of the disturbed area in the spatial distribution of the concentration of charged particles at the initial moment, which contained a circular irregularity at the center of the simulation region. Further calculations indicated that the changes in the spatial distribution of the electric charge density were continued. After a short period, the disturbed area in the spatial distribution of the electric charge density vanished almost completely, with the plasma having become almost electrically neutral in all simulation region. It turned out that, after a short period, the disturbed area in the spatial distribution of the electric charge density was recovered. Later, the cycle of vanishing and recovering of the disturbed area in the spatial distribution of the electric charge density was repeated again and again. In the course of time, the disturbed area in the spatial distribution of the electric charge density was enlarged, and additional rings were arisen around the position of the circular initial irregularity situated in the center of the simulation region. The spatial distribution of the normalized electric charge density,

Let us compare the results, obtained in the present paper for the first case, when the process starts from the completely electrically neutral state in all simulation region, with the results, obtained in the studies of Mingalev et al. [

of the fluctuations are very similar. During the processes of evolutions, around the initial irregularities, additional almost symmetrical alternate rings with an excess of charge of different sign began to appear. In the course of time, these additional rings filled up all simulation region.

The time evolutions of the magnetic field aligned super-small-scale irregularities display the following dis- tinctions caused by different initial conditions. It may be recalled that the time interval of about 35 periods of Langmuir oscillations of electrons was sufficient for the irregularity to decay completely, when the electric neutrality of the plasma was broken inside the irregularity at the initial moment [

exists during all period of calculations. As a consequence, an electrical neutrality of the plasma can not be reached because of a mobility of the electrons whose concentration deviate from the concentration of positive ions after the initial moment.

Let us consider the results, obtained in the present paper for the second case, when the process starts from the completely electrically neutral state in all simulation region and when a motion of the positive ions is taken into account. The spatial distribution of the normalized electric charge density,

The simulation results, obtained in the present study for two different cases, have essential distinctions. In the vicinity of the center of the simulation region, for the first case, when positive ions were supposed to be immovable, the normalized electric charge density stays negative during all period of calculations, whereas, for the second case, when positive ions can move, the sign of the normalized electric charge density can change in the course of time (

The time evolution of the magnetic field aligned super-small-scale irregularities in the concentration of charged particles, existing in the near-Earth rarefied plasma, was studied with the help of the model simulation. The two- dimensional mathematical model, developed earlier in the Polar Geophysical Institute, has been improved by

taking into account a motion of the positive ions and a new version of the model has been utilized to study a behavior of initially created irregularities. The utilized model is based on a numerical solution of the Vlasov- Poisson system of equations, with the Vlasov equations describing the distribution functions of charged particles and the Poisson equation governing the self-consistent electric field. The Vlasov equations are numerically solved applying a macroparticle method. The Poisson equation is solved using a finite-difference method. In the new version of the mathematical model, some parameters of the numerical method have been improved.

The dynamics of the magnetic field aligned irregularity, having circular cross sections, was studied, with the initial cross-section diameter having been equal to twelve Debye lengths of the plasma. Inside the irregularity, the electron concentration is equal to the positive ion concentration at the initial moment, with the concentrations of charged particles being disturbed. Nevertheless, the process started from the completely electrically neutral state in all simulation regions. Calculations were made using the input parameters of the model typical for the nocturnal ionospheric plasma at the altitude of 300 km .

Firstly, calculations were performed for the case when positive ions were supposed to be immovable. Simulation results indicated that the initially created irregularity vanishes and recovers periodically, with its parameters fluctuating. The fluctuations possess of two main periods, namely, the period of Langmuir oscillation of electrons and period of cyclotron oscillations of electrons. It turned out that the irregularity exists during all period of calculations (100 periods of Langmuir oscillations of electrons). The obtained results were compared with

simulation results, obtained in earlier studies devoted to investigation of the dynamics of the irregularity which had not electrical neutrality at the initial moment. For such irregularity, the time interval of about 35 periods of Langmuir oscillations of electrons was sufficient for the irregularity to decay completely. The distinction of the periods of existing may be explained by difference of the distributions of the positive ion concentration inside the irregularities at the initial moment.

Secondly, calculations were performed for the case when the process started from the completely electrically neutral state in all simulation regions and when a motion of the positive ions was taken into account. Obtained simulation results are appeared to be similar to the results obtained for the case when positive ions were supposed to be immovable. Nevertheless, essential distinctions between these simulation results exist. In particular, the behavior of some fluctuating parameters of the plasma, such as X-component of the electric field in the vicinity of the center of the simulation region and the normalized potential energy of the plasma filling up all simulation region, may be different during the initial time interval of about 50 equilibrium periods of Langmuir oscillations of electrons. Thus, the mobility of the positive ions ought to influence essentially on the time evolution of the magnetic field aligned super-small-scale irregularities in the concentration of charged particles, existing in the near-Earth rarefied plasma.

This work was partly supported by Grant No. 13-01-00063 from the Russian Foundation for Basic Research and

by the Division of the Physical Sciences of the RAS through the program “Plasma processes in space and laboratory”.