ABSTRACT It has been established that the Sweet-Parker current layer in high Lundquist number reconnection is unstable to the super-Alfvénic plasmoid instability. Past two-dimensional ...magnetohydrodynamic simulations have demonstrated that the plasmoid instability leads to a new regime where the Sweet-Parker current layer changes into a chain of plasmoids connected by secondary current sheets, and the averaged reconnection rate becomes nearly independent of the Lundquist number. In this work, a three-dimensional simulation with a guide field shows that the additional degree of freedom allows plasmoid instabilities to grow at oblique angles, which interact and lead to self-generated turbulent reconnection. The averaged reconnection rate in the self-generated turbulent state is of the order of a hundredth of the characteristic Alfvén speed, which is similar to the two-dimensional result but is an order of magnitude lower than the fastest reconnection rate reported in recent studies of externally driven three-dimensional turbulent reconnection. Kinematic and magnetic energy fluctuations both form elongated eddies along the direction of the local magnetic field, which is a signature of anisotropic magnetohydrodynamic turbulence. Both energy fluctuations satisfy power-law spectra in the inertial range, where the magnetic energy spectral index is in the range from −2.3 to −2.1, while the kinetic energy spectral index is slightly steeper, in the range from −2.5 to −2.3. The anisotropy of turbulence eddies is found to be nearly scale-independent, in contrast with the prediction of the Goldreich-Sridhar theory for anisotropic turbulence in a homogeneous plasma permeated by a uniform magnetic field.
The scaling of the plasmoid instability maximum linear growth rate with respect to the Lundquist number S in a Sweet-Parker current sheet, , indicates that at high S, the current sheet will break ...apart before it approaches the Sweet-Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We carry out a series of two-dimensional simulations and develop diagnostics to separate fluctuations from an evolving background. It is found that the fluctuation amplitude starts to grow only when the linear growth rate is sufficiently high to overcome advection loss and the stretching effect due to the outflow. The linear growth rate continues to rise until the sizes of plasmoids become comparable to the inner layer width of the tearing mode. At this point, the current sheet is disrupted and the instability enters the early nonlinear regime. The growth rate suddenly decreases, but the reconnection rate starts to rise rapidly, indicating that current sheet disruption triggers the onset of fast reconnection. We identify important timescales of the instability development, as well as scalings for the linear growth rate, current sheet width, and dominant wavenumber at disruption. These scalings depend not only on the Lundquist number, but also on the noise amplitude. A phenomenological model that reproduces scalings from simulation results is proposed. The model incorporates the effect of reconnection outflow, which is crucial for yielding a critical Lundquist number Sc below which disruption does not occur. The critical Lundquist number Sc is not a constant value, but has a weak dependence on the noise amplitude.
Magnetic reconnection and particle acceleration in relativistic Harris sheets in low-density electron-positron plasmas with no guide field have been studied by means of two-dimensional ...particle-in-cell simulations. Reconnection rates are of the order of one when the background density in a Harris sheet is of the order of 1% of the density in the current sheet, which is consistent with previous results in the non-relativistic regime. It has been demonstrated that the increase of the Lorentz factors of accelerated particles significantly enhances the collisionless resistivity needed to sustain a large reconnection electric field. It is shown analytically and numerically that the energy spectrum of accelerated particles near the X-line is the product of a power law and an exponential function of energy, gamma super(-1/4) exp(-a gamma super(1/2), where gamma is the Lorentz factor and a is a constant. However, in the low-density regime, while the most energetic particles are produced near X-lines, many more particles are energized within magnetic islands. Particles are energized in contracting islands by multiple reflection, but the mechanism is different from Fermi acceleration in magnetic islands for magnetized particles in the presence of a guide field. In magnetic islands, strong core fields are generated and plasma beta values are reduced. As a consequence, the fire-hose instability condition is not satisfied in most of the island region, and island contraction and particle acceleration can continue. In island coalescence, reconnection between two islands can accelerate some particles, however, many particles are decelerated and cooled, which is contrary to what has been discussed in the literature on particle acceleration due to reconnection in non-relativistic hydrogen plasmas.
The distribution function f(ψ) of magnetic flux ψ in plasmoids formed in high-Lundquist-number current sheets is studied by means of an analytic phenomenological model and direct numerical ...simulations. The distribution function is shown to follow a power law f(ψ)∼ψ(-1), which differs from other recent theoretical predictions. Physical explanations are given for the discrepant predictions of other theoretical models.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We propose a new mechanism for a turbulent mean-field dynamo in which the magnetic fluctuations resulting from a small-scale dynamo drive the generation of large-scale magnetic fields. This is in ...stark contrast to the common idea that small-scale magnetic fields should be harmful to large-scale dynamo action. These dynamos occur in the presence of a large-scale velocity shear and do not require net helicity, resulting from off-diagonal components of the turbulent resistivity tensor as the magnetic analogue of the "shear-current" effect. Given the inevitable existence of nonhelical small-scale magnetic fields in turbulent plasmas, as well as the generic nature of velocity shear, the suggested mechanism may help explain the generation of large-scale magnetic fields across a wide range of astrophysical objects.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Plasmoid Instability in Forming Current Sheets Comisso, L.; Lingam, M.; Huang, Y.-M. ...
Astrophysical journal/The Astrophysical journal,
12/2017, Volume:
850, Issue:
2
Journal Article
Peer reviewed
Open access
The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is ...crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most previous studies have focused on Sweet-Parker current sheets, which are unattainable in typical astrophysical systems. Here we derive a general set of scaling laws for the plasmoid instability in resistive and visco-resistive current sheets that evolve over time. Our method relies on a principle of least time that enables us to determine the properties of the reconnecting current sheet (aspect ratio and elapsed time) and the plasmoid instability (growth rate, wavenumber, inner layer width) at the end of the linear phase. After this phase the reconnecting current sheet is disrupted and fast reconnection can occur. The scaling laws of the plasmoid instability are not simple power laws, and they depend on the Lundquist number (S), the magnetic Prandtl number (Pm), the noise of the system ( ), the characteristic rate of current sheet evolution ( ), and the thinning process. We also demonstrate that previous scalings are inapplicable to the vast majority of astrophysical systems. We explore the implications of the new scaling relations in astrophysical systems such as the solar corona and the interstellar medium. In both of these systems, we show that our scaling laws yield values for the growth rate, wavenumber, and aspect ratio that are much smaller than the Sweet-Parker-based scalings.