Plasmoid Formation Research Articles

Extreme-ultraviolet (EUV) observables of simulated plasmoid-mediated reconnection in the solar corona

Context. Understanding the role of magnetic reconnection in the heating and dynamics of the solar atmosphere requires detailed observational data of any observable aspect of the reconnection process, including small-scale features such as plasmoids. Aims. Here, we examine the capability of active and upcoming instruments to detect plasmoids generated by reconnection in the corona including low-density regimes. Methods. We used the Bifrost code to perform simulations of plasmoid-mediated reconnection in the corona with a 2D idealized setup: a fan-spine topology with uniform density including thermal conduction. Through a forward-modeling of extreme-ultraviolet (EUV) observables, we checked whether our simulated plasmoids could be detected with the instruments of Solar Dynamics Observatory (SDO) and Solar Orbiter (SO), as well as the upcoming Multi-Slit Solar Explorer (MUSE) and Solar-C missions. Results. Short-lived (∼10 − 20 s) small-scale (∼0.2 − 0.5 Mm) coronal plasmoids are not resolvable with the Atmospheric Imaging Assembly (AIA) on board SDO. In contrast, they could be captured with the EUV High-Resolution Imager at the Extreme Ultraviolet Imager (EUI-HRIEUV) of SO. The spatial and temporal high-resolution planned for the MUSE spectrograph (SG) is adequate to obtain full spectral information of these plasmoids. To achieve a sufficient signal-to-noise ratio (S/N) for ∼0.8 MK plasmoids in the MUSE/SG 171 Å channel, full-raster images are attainable for regions with electron densities above 109 cm−3, while sit-and-stare observations are recommended for lower-density regions. The future Solar-C mission could also capture these coronal plasmoids using the EUV High-Throughput Spectroscopic Telescope (EUVST), considering rapid changes in Doppler shift and line widths in different EUV lines caused by plasmoid motions along the current sheet. Conclusions. With the combined spectra of MUSE/SG and Solar-C/EUVST in multiple emission lines, along with high-resolution images from SO/EUI-HRIEUV and MUSE/CI, it should be possible to gain new insights about plasmoid formation in the corona.

Plasmoid identification and statistics in two-dimensional Harris sheet and GRMHD simulations

ABSTRACT Magnetic reconnection is a ubiquitous phenomenon for magnetized plasma and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These so-called plasmoids could therefore observationally manifest themselves as hotspots, which are associated with flaring behaviour in supermassive black hole systems, such as Sagittarius A*. We have developed a novel algorithm for identifying plasmoid structures, which incorporates watershed and custom closed contouring steps. From the identified structures, we determine the plasmoids’ plasma characteristics and energetics in magnetohydrodynamical simulations. The algorithm’s performance is showcased for a high-resolution suite of axisymmetric ideal and resistive magnetohydrodynamical simulations of turbulent accretion discs surrounding a supermassive black hole. For validation purposes, we also evaluate several Harris current sheets that are well-investigated in the literature. We recover the characteristic power-law distribution of plasmoid sizes for both the black hole and Harris sheet simulations. This indicates that while the dynamics are vastly different, with different dominant plasma instabilities, the plasmoid creation characteristics are similar. Plasmoid formation rates for resistive general relativistic magnetohydrodynamical simulations are significantly higher than for their ideal counterpart. Moreover, the largest identified plasmoids are consistent with sizes typically assumed for semi-analytical interpretation of observations. We recover a positive correlation between the plasmoid formation rate and decreases in black-hole-horizon-penetrating magnetic flux, during which the accretion flow is temporarily halted. These results demonstrate the efficacy of the newly developed algorithm which has enabled an extensive quantitative analysis of plasmoid formation in black hole accretion simulations.

Plasmoid Formation and Strong Radiative Cooling in a Driven Magnetic Reconnection Experiment.

We present the first experimental study of plasmoid formation in a magnetic reconnection layer undergoing rapid radiative cooling, a regime relevant to extreme astrophysical plasmas. Two exploding aluminum wire arrays, driven by the Z machine, generate a reconnection layer (S_{L}≈120) in which the cooling rate far exceeds the hydrodynamic transit rate (τ_{hydro}/τ_{cool}>100). The reconnection layer generates a transient burst of >1 keV x-ray emission, consistent with the formation and subsequent rapid cooling of the layer. Time-gated x-ray images show fast-moving (up to 50 km s^{-1}) hotspots in the layer, consistent with the presence of plasmoids in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (around 1.6keV) prior to the onset of cooling, and exhibit temperatures (170eV) much greater than that of the plasma inflows and the rest of the reconnection layer, thus providing insight into the generation of high-energy radiation in radiatively cooled reconnection events.

2.5D Magnetohydrodynamic Simulation of the Formation and Evolution of Plasmoids in Coronal Current Sheets

In the present paper, using MPI-AMRVAC, we perform a 2.5D numerical magnetohydrodynamic simulation of the dynamics and associated thermodynamical evolution of an initially force-free Harris current sheet subjected to an external velocity perturbation under the condition of uniform resistivity. The amplitude of the magnetic field is taken to be 10 G, typical of the solar corona. We impose a Gaussian velocity pulse across this current sheet that mimics the interaction of fast magnetoacoustic waves with a current sheet in the corona. This leads to a variety of dynamics and plasma processes in the current sheet, which is initially quasi-static. The initial pulse interacts with the current sheet and splits into a pair of counterpropagating wavefronts, which form a rarefied region that leads to an inflow and a thinning of the current sheet. The thinning results in Petschek-type magnetic reconnection followed by a tearing instability and plasmoid formation. The reconnection outflows containing outward-moving plasmoids have accelerated motions with velocities ranging from 105 to 303 km s−1. The average temperature and density of the plasmoids are found to be 8 MK and twice the background density of the solar corona, respectively. These estimates of the velocity, temperature, and density of the plasmoids are similar to values reported from various solar coronal observations. Therefore, we infer that the external triggering of a quasi-static current sheet by a single-velocity pulse is capable of initiating magnetic reconnection and plasmoid formation in the absence of a localized enhancement of resistivity in the solar corona.

Characteristics of Foreshock Subsolar Compressive Structures

AbstractThe turbulent foreshock region upstream of the quasi‐parallel bow shock is dominated by waves and reflected particles that interact with each other and create a large number of different foreshock transients. The structures with the enhanced magnetic field, Short Large Amplitude Magnetic Structures, and density spikes named plasmoids are frequently observed. They are one of the suggested sources of transient flux enhancements or jets in the magnetosheath. Using measurements of the Magnetospheric Multiscale Spacecraft (MMS) and OMNI solar wind database between the 2015 and 2018 years, we have found that there is a category of events exhibiting both magnetic field and density enhancements simultaneously and we introduce the term “mixed structures” for them. Consequently, we divided our set of observations into three groups of events and present their comparative statistical analysis in the subsolar foreshock. Based on our results and previous research, we discuss the properties, possible origin and occurrence rates of these events under different upstream conditions and their possible relation to the jet and plasmoid formation in the magnetosheath.

Species dependence of the impurity injection induced poloidal flow and magnetic island rotation in a tokamak

Recent experiments have demonstrated the species dependence of the impurity poloidal drift direction along with the magnetic island rotation in the poloidal plane. Our resistive MHD simulations have reproduced such a dependence of the impurity poloidal flow, which is found mainly determined by a local plasmoid formation due to the impurity injection. The synchronized magnetic island rotation is dominantly driven by the electromagnetic torque produced by the impurity radiation primarily through the modification to the axisymmetric components of current density.

Two-fluid reconnection jets in a gravitationally stratified atmosphere

Context. Density decreases exponentially with height in the gravitationally stratified solar atmosphere, and therefore collisional coupling between the ionized plasma and the neutrals also decreases. Reconnection is a process observed at all heights in the solar atmosphere. Aims. Here, we investigate the role of collisions between ions and neutrals in the reconnection process occurring at various heights in the atmosphere. Methods. We performed simulations of magnetic reconnection induced by a localized resistivity in a gravitationally stratified atmosphere, in which we varied the height of the initial reconnection X-point. We compared a magnetohydrodynamic (MHD) model and two two-fluid configurations: one in which the collisional coupling was calculated from local plasma parameters, and another in which the coupling was decreased so that collisional effects would be enhanced. The latter setup has a more representative solar collisionality regime. Results. Simulations in a stratified atmosphere show similar structures in MHD and two-fluid simulations, with strong coupling. However, when collisional effects are increased to attain representative parameter regimes, we find a nonlinear runaway instability, which separates the plasma-neutral densities across the current sheet (CS). With increased collisional effects, the initial decoupling in velocity heats the neutrals and this sets up a nonlinear feedback loop, according to which neutrals migrate outside the CS, replacing charged particles that accumulate toward the center of the CS. Conclusions. The reconnection rate has a maximum value of around 0.1 for both reconnection heights, and is consistent with the locally enhanced resistivity used in all three models. The early-stage plasmoid formation observed near the end of our simulations is influenced by the outflow from the primary reconnection point, rather than by collisions. We synthesized optically thin emission for both MHD and two-fluid models, which can show a very different evolution when the charged-particle density is used instead of the total density. Our simulations have relevance for observed plasmoid features associated with chromospheric to low-coronal flare events.

Viscous Effects on Nonlinear Double Tearing Mode and Plasmoid Formation in Adjacent Harris Sheets

In this paper, we study the effects of viscosity on the evolution of the double tearing mode (DTM) in a pair of adjacent Harris sheets based on the resistive MHD model in the NIMROD code. Similar to the tearing mode in the conventional single Harris sheet, a transition is observed in the generation of both normal and monster plasmoids at Prandtl number Pr=1. In the Pr<1 regime of the DTM, normal plasmoids (small plasmoids) are generated along with monster plasmoid, whereas in the single tearing mode (STM) cases, such a generation is not observed. When Pr is above the critical value, the generation of monster plasmoid is halted. Correspondingly, in the Pr<1 regime, a quadrupolar flow advects along the poloidal direction, but in the Pr>1 regime this flow advection is inhibited.

Two-temperature GRMHD simulations of black hole accretion flows with multiple magnetic loops

ABSTRACT We have performed a series of two-dimensional two-temperature general relativistic magnetohydrodynamic simulations of magnetized accretion flows initiated from tori with different sizes and poloidal magnetic loop polarities. In these two-temperature simulations, we trace the process of heating electrons through turbulence and reconnection, most of the time these electrons are trapped in plasmoids. We found that the accretion process strongly depends on the size of the magnetic loops. The accretion flows never reach the magnetically arrested (MAD) regime in small loop cases. Interaction between magnetic field with different polarities dissipates and decreases the efficiency of magnetorotational instability. The dependence on the wavelength of the loops places a lower limit on the loop size. In the large loop cases, after reaching a quasi-steady phase, a transition from standard and normal evolution flow to MAD flow is observed. The transition of the accretion state and the transition time depends on the initial loop wavelength. The formation of plasmoids strongly depends on the size of the magnetic loops. The frequent magnetic reconnection between the magnetic loops is responsible for the formation of most of the plasmoids. For some plasmoids, Kelvin–Helmholtz and tearing instabilities are coexisting, showing another channel of plasmoid formation. The simulations present that electrons in the plasmoids are well-heated up by turbulent and magnetic reconnection. Different properties of plasmoid formation in different magnetic field configurations provide new insights for the understanding of flaring activity and electron thermodynamics in Sagittarius A*.

Closed field line vortices in planetary magnetospheres

ABSTRACTIn a rotation-dominated magnetosphere, there is a region where closed field lines rotate around the planet, and also a region where the open field lines stretch away from the planet, forming the lobes of the magnetotail. This paper shows that there could be a third, significantly different region, where the closed field lines form twisted vortex structures anchored in the magnetotail. Such patterns form when there are significant plasma sources inside the magnetosphere and the time-scale of the plasmoid formation process is substantially larger than the planetary rotation period. In the presence of vortices, the Dungey and Vasyliunas cycles act differently. The Dungey flow does not penetrate the central region of the polar cap. Tail reconnection events are rare, thus leaving the plasma time enough to participate in the essentially 3D vortex-forming plasma motion. The above conditions are fulfilled for Saturn. We discovered vortex-like patterns in the plasma and magnetic field data measured by the Cassini spacecraft in the nightside magnetosphere of Saturn. The plasma whirling around in these vortices never reaches the dayside, instead, it performs a retrograde motion in the high latitude regions of the magnetotail. Low-energy plasma data suggest that the observed patterns correspond to the closed field line vortices.

Kilotesla plasmoid formation by a trapped relativistic laser beam.

A strong quasistationary magnetic field is generated in hollow targets with curved internal surface under the action of a relativistically intense picosecond laser pulse. Experimental data evidence the formation of quasistationary strongly magnetized plasma structures decaying on a hundred picoseconds timescale, with the magnetic field strength of the kilotesla scale. Numerical simulations unravel the importance of transient processes during the magnetic field generation and suggest the existence of fast and slow regimes of plasmoid evolution depending on the interaction parameters. The proposed setup is suited for perspective highly magnetized plasma application and fundamental studies.

MHD turbulence: a biased review

This review of scaling theories of magnetohydrodynamic (MHD) turbulence aims to put the developments of the last few years in the context of the canonical time line (from Kolmogorov to Iroshnikov–Kraichnan to Goldreich–Sridhar to Boldyrev). It is argued that Beresnyak's (valid) objection that Boldyrev's alignment theory, at least in its original form, violates the Reduced-MHD rescaling symmetry can be reconciled with alignment if the latter is understood as an intermittency effect. Boldyrev's scalings, a version of which is recovered in this interpretation, and the concept of dynamic alignment (equivalently, local 3D anisotropy) are thus an example of a physical theory of intermittency in a turbulent system. The emergence of aligned structures naturally brings into play reconnection physics and thus the theory of MHD turbulence becomes intertwined with the physics of tearing, current-sheet disruption and plasmoid formation. Recent work on these subjects by Loureiro, Malletet al.is reviewed and it is argued that we may, as a result, finally have a reasonably complete picture of the MHD turbulent cascade (forced, balanced, and in the presence of a strong mean field) all the way to the dissipation scale. This picture appears to reconcile Beresnyak's advocacy of the Kolmogorov scaling of the dissipation cutoff (as$\mathrm {Re}^{3/4}$) with Boldyrev's aligned cascade. It turns out also that these ideas open the door to some progress in understanding MHD turbulence without a mean field – MHD dynamo – whose saturated state is argued to be controlled by reconnection and to contain, at small scales, a tearing-mediated cascade similar to its strong-mean-field counterpart (this is a new result). On the margins of this core narrative, standard weak-MHD-turbulence theory is argued to require some adjustment – and a new scheme for such an adjustment is proposed – to take account of the determining part that a spontaneously emergent 2D condensate plays in mediating the Alfvén-wave cascade from a weakly interacting state to a strongly turbulent (critically balanced) one. This completes the picture of the MHD cascade at large scales. A number of outstanding issues are surveyed: imbalanced turbulence (for which a new, tentative theory is proposed), residual energy, MHD turbulence at subviscous scales, and decaying MHD turbulence (where there has been dramatic progress recently, and reconnection again turned out to feature prominently). Finally, it is argued that the natural direction of research is now away from the fluid MHD theory and into kinetic territory – and then, possibly, back again. The review lays no claim to objectivity or completeness, focusing on topics and views that the author finds most appealing at the present moment.

Thermally enhanced tearing in solar current sheets: Explosive reconnection with plasmoid-trapped condensations

Context. Thermal instability plays a major role in condensation phenomena in the solar corona (e.g., for coronal rain and prominence formation). In flare-relevant current sheets, tearing instability may trigger explosive reconnection and plasmoid formation. However, how the two instabilities influence the disruption of current concentrations in the solar corona has received less attention to date. Aims. We explore how the thermal and tearing modes reinforce each other in the fragmentation of a current sheet in the solar corona through an explosive reconnection process, characterized by the formation of plasmoids which interact and trap condensing plasma. Methods. We used a resistive magnetohydrodynamic simulation of a 2D current layer, incorporating the non-adiabatic effects of optically thin radiative energy loss and background heating using the open-source code MPI-AMRVAC. Multiple levels of adaptive mesh refined grids are used to achieve a high resolution to resolve the fine structures during the evolution of the system. Results. Our parametric survey explores different resistivities and plasma-β to quantify the instability growth rate in the linear and nonlinear regimes. We note that for dimensionless resistivity values within 10−4 − 5 × 10−3, we obtain explosive behavior where thermal instability and tearing behavior reinforce each other. This is clearly below the usual critical Lundquist number range of pure resistive explosive plasmoid formation. We calculate the mean growth rate for the linear phase and different nonlinear phases of the evolution. The nonlinear growth rates follow a weak power-law dependency with resistivity. The fragmentation of the current sheet and the formation of the plasmoids in the nonlinear phase of the evolution due to the thermal and tearing instabilities are obtained. The formation of plasmoids is found for the Lundquist number (SL) range 4.6 × 103 − 2.34 × 105. We quantify the temporal variation in the plasmoid numbers and the density filling factor of the plasmoids for different physical conditions. We also find that the maximum plasmoid numbers scale as SL0.223.Within the nonlinearly coalescing plasmoid chains, localized cool condensations gather, realizing density and temperature contrasts similar to coronal rain or prominences.

Magnetic reconnection and plasmoid formation in three-dimensional accretion flows around black holes

ABSTRACT Magnetic reconnection is thought to be one of the main energy-dissipation mechanisms fueling energy to the plasma in the vicinity of a black hole. Indeed, plasmoids formed through magnetic reconnection may play a key role in γ-ray, X-ray and near-infrared flares from the black hole at the centre of our galaxy, SgrA*. We report the results of 3D general-relativistic ideal and resistive magnetohydrodynamics simulations modelling magnetic reconnection in accretion flows around astrophysical black holes. As an important difference with similar works, our accretion discs have an initial dipolar magnetic-field configuration with loops of alternating polarity. We show that current sheets are formed and destroyed rapidly in the turbulent environment of black hole accretion. Plasmoids are formed from current sheets close to the event horizon, in a region of ∼2–15 gravitational radii. We further quantify the magnetic dissipation and the process of energy transfer to the plasmoids, reporting the reconnection rate, the relative current density with respect to the local magnetic field, and the size of the plasmoids. We find that plasmoids gain energy through reconnection and heat up to relativistic temperatures, with the largest ones being sufficiently energetic to leave the black hole near the polar regions. During their evolution, plasmoids are stretched and elongated, becoming disrupted when the shear is sufficiently large, although some plasmoids survive as well-distinguished structures at distances of ∼30–40 gravitational radii from the black hole. Finally, we find that in some cases the plasmoids acquire a super-Keplerian azimuthal velocity, as suggested by recent observations of flares from Sgr A*.

Intra-pulse variability induced by plasmoid formation in pulsar magnetospheres

Context. Pulsars show irregularities in their pulsed radio emission that originate from propagation effects and the intrinsic activity of the source. Aims. In this work, we investigate the role played by magnetic reconnection and the formation of plasmoids in the pulsar wind current sheet as a possible source of intrinsic pulse-to-pulse variability in the incoherent, high-energy emission pattern. Methods. We used a two-dimensional particle-in-cell simulation of an orthogonal pulsar magnetosphere restricted to the plane perpendicular to the star spin axis. We evolved the solution for several tens of pulsar periods to gather a statistically significant sample of synthetic pulse profiles. Results. The formation of plasmoids leads to strong pulse-to-pulse variability in the form of multiple short, bright subpulses, which appear only on the leading edge of each main pulse. These secondary peaks of emission are dominated by the dozen plasmoids that can grow up to macroscopic scales. They emerge from the high end of the hierarchical merging process occurring along the wind current layer. The flux of the subpulses is correlated with their width in phase. Although the full-scale separation is not realistic, we argue that the simulation correctly captures the demographics and the properties of the largest plasmoids, and therefore of the brightest subpulses. Conclusions. The prediction of subpulses at specific pulse phases provides a new observational test of the magnetic reconnection scenario as the origin of the pulsed incoherent emission. High-time-resolution observations of the Crab pulsar in the optical range may be the most promising source to target for this purpose.

Coexistence of Plasmoid and Kelvin–Helmholtz Instabilities in Collisionless Plasma Turbulence

The plasmoid formation in collisionless plasmas, where magnetic reconnection within turbulence may take place driven by the electron inertia, is analyzed. We find a complex situation in which, due to the presence of strong velocity shears, the typical plasmoid formation, observed to influence the energy cascade in the magnetohydrodynamic context, has to coexist with the Kelvin–Helmholtz (KH) instability. We find that the current density layers may undergo the plasmoid or the KH instability depending on the local values of the magnetic and velocity fields. The competition among these instabilities affects not only the evolution of the current sheets, that may generate plasmoid chains or KH-driven vortices, but also the energy cascade, that is different for the magnetic and kinetic spectra.

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

Effect of Magnetic and Electric Fields on the Dynamics of Plasmoid Formation in the Gatchina Discharge

Effect of Magnetic and Electric Fields on the Dynamics of Plasmoid Formation in the Gatchina Discharge

Influence of aspect ratio, plasma viscosity, and radial position of the resonant surfaces on the plasmoid formation in the low resistivity plasma in Tokamak

In the present paper, we systematically investigate the nonlinear evolution of the resistive kink mode in the low resistivity plasma in Tokamak geometry. We find that the aspect ratio of the initial equilibrium can significantly influence the critical resistivity for plasmoid formation. With the aspect ratio of 3/1, the critical resistivity can be one magnitude larger than that in cylindrical geometry due to the strong mode–mode coupling. We also find that the critical resistivity for plasmoid formation η crit decreases with increasing plasma viscosity in the moderately low resistivity regime. Due to the geometry of Tokamaks, the critical resistivity for plasmoid formation increases with the increasing radial location of the resonant surface.

Connecting theory of plasmoid-modulated reconnection to observations of solar flares

Abstract The short timescale of the solar flare reconnection process has long proved to be a puzzle. Recent studies suggest the importance of the formation of plasmoids in the reconnecting current sheet, with quantifying the aspect ratio of the width to length of the current sheet in terms of a negative power $ \alpha $ of the Lundquist number, that is, $ {S}^{-\alpha } $ , being key to understanding the onset of plasmoids formation. In this paper, we make the first application of theoretical scalings for this aspect ratio to observed flares to evaluate how plasmoid formation may connect with observations. For three different flares that show plasmoids we find a range of $ \alpha $ values of $ \alpha =0.26 $ to $ 0.31 $ . The values in this small range implies that plasmoids may be forming before the theoretically predicted critical aspect ratio ( $ \alpha =1/3 $ ) has been reached, potentially presenting a challenge for the theoretical models.