Effect of Charge accumulation on Magnetic Rayleigh-Taylor Instability

The intuitive physical description of magnetic Rayleigh-Taylor instability in some textbooks is generally considered to be: a small perturbation causes current discontinuity, which produce charge accumulation, the electric field produced by the accumulated charge amplify the initial perturbation. However, in calculating the linear growth rate of magnetic Rayleigh-Taylor instability (MRTI), the displacement current term in the Maxwell’s equations is ignored, which means the contribution of charge accumulation to the growth of MRTI is totally ignored. In this article, we calculated the linear growth rate of MRTI with the displacement current term in Maxwell’s equations retained. We show that the contribution of charge accumulation to the growth of MRTI is negligible only when the nominal Alfvén speed is much smaller than the light speed. For space plasma whose nominal Alfvén speed is generally much smaller than the light speed, the linear growth rate previous calculated is right but the intuitive physical description of MRTI is wrong. For laboratory plasma whose nominal Alfvén speed maybe comparable to light speed, the intuitive physical description of MRTI is also inaccurate and the linear growth rate of MRTI is undervalued.


Effect of Charge accumulation on Magnetic Rayleigh-Taylor Instability
Kangkang Liu 1,2 The intuitive physical description of magnetic Rayleigh-Taylor instability in some textbooks is generally considered to be: a small perturbation causes current discontinuity, which produce charge accumulation, the electric field produced by the accumulated charge amplify the initial perturbation. However, in calculating the linear growth rate of magnetic Rayleigh-Taylor instability (MRTI), the displacement current term in the Maxwell's equations is ignored, which means the contribution of charge accumulation to the growth of MRTI is totally ignored. In this article, we calculated the linear growth rate of MRTI with the displacement current term in Maxwell's equations retained. We show that the contribution of charge accumulation to the growth of MRTI is negligible only when the nominal Alfvén speed is much smaller than the light speed. For space plasma whose nominal Alfvén speed is generally much smaller than the light speed, the linear growth rate previous calculated is right but the intuitive physical description of MRTI is wrong. For laboratory plasma whose nominal Alfvén speed maybe comparable to light speed, the intuitive physical description of MRTI is also inaccurate and the linear growth rate of MRTI is undervalued.
When a heavier fluid is supported by a lighter fluid against gravity, the equilibrium is unstable to any perturbations of the interface. For if a parcel of heavier fluid is displaced downward with an equal volume of lighter fluid displaced upwards, the potential energy of the configuration is lower than that in the initial state, and the process goes on. This instability is called the Rayleigh-Taylor Instability (RTI) 1 . In addition to neutral fluids, RTI also plays an important role in space and laboratory plasmas [2][3][4][5][6] .
For the neutral fluid, the momentum equation is , where ρ, V, g, p is the plasma density, the plasma bulk velocity, the gravity acceleration and the thermal pressure, respectively. In equilibrium state ρg = ∇p. The heavier fluids are supported by lighter fluids through pressure gradient in equilibrium state. For plasma in magnetic field, the momentum equation is , where c, J, B is the light speed, electric current and magnetic field strength, respectively. In equilibrium The heavier fluids are supported by lighter fluids through pressure gradient and the J × B force in equilibrium state. The intuitive physical description of magnetic Rayleigh-Taylor instability (MRTI) in some textbooks 7 is generally considered to be ( Fig. 1): In equilibrium state a net current flows in in the horizontal direction and the current is proportional to the plasma density. There is thus a divergence, and charge will pile up on the edges of the small initial perturbation (charge accumulation). As a result, perturbation electric fields build up in the directions shown. These fields in turn cause an upward (downward) drift in the region where the density is low (high). Lower (higher) density plasma is therefore advected upward (downward), creating a larger perturbation, and the system is unstable. In the above physical description charge accumulation is the cause for the growth of RTI. However, in calculating the linear growth rate of magnetic Rayleigh-Taylor instability (MRTI), the displacement current term in the Maxwell's equations is ignored 8,9 . When the displacement current term in the Maxwell's equations is ignored, take the divergence of the Ampère's circuital law without Maxwell's addition shows that the current is always continuous. Which means the contribution of charge accumulation to the growth of MRTI is totally ignored. In this paper, we derived the expression of the linear growth rate of MRTI when charge accumulation is considered and discussed under what circumstances the displacement current can be ignored.
Assume a bundle of stratified plasma ρ = ρ(y) immersed in a magnetic field B (0, 0, B) and a gravity field g (0, −g, 0) (Fig. 1). To analyze the instability of the system, the relevant magnetohydrodynamic (MHD) equations and Maxwell's equations written in Gaussian unit are: where V, ρ, c, J, B, p, E, ρ c are the velocity of fluid element, plasma mass density, light speed, electric current density, magnetic field, thermal pressure, electric field, charge density, respectively. The following derivation process of the linear growth rate of stratified plasma in a horizontal magnetic field in this manuscript was similar to that of Sharma et al. 10 and Wu et al. 11 , except that charge conservation equation (Eq. (6)) instead of the current continuity equation was used in the derivation process. Form the momentum equation (Eq. (1)), one can see that the three term on the right side of the equation was equivalent. If charge accumulation effects considered, the J × B term will affect the linear growth rate of magnetic Rayleigh-Taylor instability.
To examine the stability of the system, we assume the following perturbation in physical quantities In the above equations, quantities with superscripts 0 and 1 denote equilibrium and perturbation quantities, respectively.
Assuming perturbations are in the form i kx t ( ) where ω is the frequency of the perturbation, k is the wave number. Linearizing the Eq. (1), we get  www.nature.com/scientificreports www.nature.com/scientificreports/ where z is the unit vector in z direction. From Eqs (5) and (6), we get is the charge accumulation in (outside) the fluid element. Substituting Eq. (10) in Eq. (9) we get Supposing that the plasma is incompressible we get the following equation x y From the linearized continuity equation 1 0 We get x y 0 Substituting the Eqs (13) and (16)  where s is a constant, Eq. (17) can be written as   www.nature.com/scientificreports www.nature.com/scientificreports/ The full expression of the linear growth rate of MRTI is To investigate the effect of charge accumulation on the linear growth rate of RTI, we nondimensionalize Eq. (21) with the following expressions where ω pe is the plasma frequency. We get  15 , and the growth rate is the same as that of the classical RTI. When V A > 0, the growth rate is larger than that of the classical RTI, and increases with the increase of V /c A 2 2 . In planetary plasma, in general  V c A , the growth rate is basically the same as that of the classical RTI. The contribution of charge accumulation to the growth of MRTI can be neglected. Ignore the displacement current term in Maxwell's equations when calculate the linear growth rate of MRTI is acceptable. However, the physical description of the MRTI in planetary plasma which attribute the growth of MRTI to charge accumulation is wrong. The physical description of the growth of MRTI is similar to that of neutral fluid 16,17 . In laboratory plasma, the nominal Alfvén speed is not subject to the constraint  V c A , the growth rate can be much larger than that of the classical RTI. Under such condition, the charge accumulation effect and the displacement current term cannot be ignored. However, the physical description in Fig. 1 which attribute the growth of MRTI to charge accumulation is not accurate, the mechanism for the growth of RTI in neutral fluid works too 16,17 . Also from Eq. (22) one can see that when the wave number k is very large, the linear growth rate reduce to that of classical RTI, and the effect of charge accumulation can be neglected.
Time scale analysis showed that everything involving charge separation happens on time scales of the inverse plasma frequency, and in such short time scales, the displacement current term cannot be ignored 18,19 . The MRTI process involves charge accumulation, so the displacement current term should not be ignored during the calculation of the growth rate of MRTI. The displacement current term is neglected because it's generally small compared to J and curl B term. However, as we can see from Eq. (3) that the ∂ ∂t E has the same order of magnitude as J, and the changing electric field have significant effects on the dynamics of the plasma [20][21][22] . Ignore the displacement www.nature.com/scientificreports www.nature.com/scientificreports/ current term exclude the effects of electric field on the dynamics of the plasma. It is generally considered that the displacement current, which is closely related to charge separation, must be become negligible under the assumption of quasi-neutrality. However, studies show that when the nominal Alfvén speed is comparable to or larger than the light speed, the displacement current cannot be neglected, regardless of the quasi-neutrality considerations 16,23,24 , and this is consistent with our result discussed above.