The tearing mode instability is one important mechanism that may explain the triggering of fast magnetic reconnection in astrophysical plasmas such as the solar corona. In this paper, the linear stability analysis of the tearing mode is carried out for a current sheet in the presence of a guide field, and including the Hall effect. We show that the presence of a strong guide field does not modify the most unstable mode in the two-dimensional wave vector space orthogonal to the current gradient direction, which remains the fastest growing parallel mode. With the Hall effect, the inclusion of a guide field turns the non-dispersive propagation along the guide field direction to a dispersive one. The oblique modes have a wave-like structure along the normal direction of the current sheet and a strong guide field suppresses this structure while making the eigen-functions asymmetric.

In many astrophysical environments the plasma is only partially ionized, and therefore the interaction of charged and neutral particles may alter both the triggering of reconnection and its subsequent dynamical evolution. We derive the tearing mode maximum growth rate for partially ionized plasmas in the cases of weak and strong coupling between the plasma and the neutrals. In addition, critical scalings for current sheet aspect ratios are presented in terms of Lundquist number and ion-neutral collision frequencies. In the decoupled regime the standard tearing mode is recovered with a small correction depending on the ion-neutral collision frequency; in the intermediate regime collisions with neutrals are shown to stabilize current sheets, resulting in larger critical aspect ratios for ideal tearing to occur. Nonetheless, the additional electron-neutral collisions, hidden in the definition of the Lundquist number, can shrink the critical aspect ratios below the fully ionized case. In the coupled regime, the growth rate depends on the density ratio between ions and neutrals through the collision frequency between these two species. These provide critical aspect ratios for which the tearing mode instability transitions from slow to ideal, that depend on the neutral-ion density ratio.

In contrast with the fast solar wind, that originates in coronal holes, the source of the slow solar wind is still debated. Often intermittent and enriched with low FIP elements -- akin to what is observed in closed coronal loops -- the slow wind could form in bursty events nearby helmet streamers. Slow winds also exhibit density perturbations which have been shown to be periodic and could be associated with flux ropes ejected from the tip of helmet streamers, as shown recently by the WISPR white light imager onboard Parker Solar Probe (PSP). In this work, we propose that the main mechanism controlling the release of flux ropes is a flow-modified tearing mode at the heliospheric current sheet (HCS). We use MHD simulations of the solar wind and corona to reproduce realistic configurations and outflows surrounding the HCS. We find that this process is able to explain long ($\sim 10-20$h) and short ($\sim 1-2$h) timescales of density structures observed in the slow solar wind. This study also sheds new light on the structure, topology and composition of the slow solar wind, and could be, in the near future, compared with white light and in situ PSP observations.

During Parker Solar Probe's first orbit, the solar wind plasma has been observed in situ closer than ever before, the perihelion on November 6th 2018 revealing a flow that is constantly permeated by large amplitude Alfv\'enic fluctuations. These include radial magnetic field reversals, or switchbacks, that seem to be a persistent feature of the young solar wind. The measurements also reveal a very strong, unexpected, azimuthal velocity component. In this work, we numerically model the solar corona during this first encounter, solving the MHD equations and accounting for Alfv\'en wave transport and dissipation. We find that the large scale plasma parameters are well reproduced, allowing the computation of the solar wind sources at Probe with confidence. We try to understand the dynamical nature of the solar wind to explain both the amplitude of the observed radial magnetic field and of the azimuthal velocities.

We carry out two-dimensional magnetohydrodynamic (MHD) simulations of an ensemble of Alfv\'enic fluctuations propagating in a structured, expanding solar wind including the presence of fast and slow solar wind streams. Using an appropriate expanding box model, the simulations incorporate the effects of fast-slow stream shear and compression and rarefaction self-consistently. We investigate the radial and longitudinal evolution of the cross-helicity, the total and residual energies and the power spectra of outward and inward Alfv\'enic fluctuations. The stream interaction is found to strongly affect the radial evolution of Alfv\'enic turbulence. The total energy in the Alfv\'en waves is depleted within the velocity shear regions, accompanied by the decrease of the normalized cross-helicity. The presence of stream-compression facilitates this process. Residual energy fluctuates around zero due to the correlation and de-correlation between the inward/outward waves but no net growth or decrease of the residual energy is observed. The radial power spectra of the inward/outward Alfv\'en waves show significant longitudinal variations. Kolmogorov-like spectra are developed only inside the fast and slow streams and when both the compression and shear are present. On the other hand, the spectra along the longitudinal direction show clear Kolmogorov-like inertial ranges in all cases.

Parker Solar Probe (PSP) observations during its first encounter at 35.7 $R_\odot$ have shown the presence of magnetic field lines which are strongly perturbed to the point that they produce local inversions of the radial magnetic field, known as switchbacks. Their counterparts in the solar wind velocity field are local enhancements in the radial speed, or jets, displaying (in all components) the velocity-magnetic field correlation typical of large amplitude Alfv\'en waves propagating away from the Sun. Switchbacks and radial jets have previously been observed over a wide range of heliocentric distances by Helios, WIND and Ulysses, although they were prevalent in significantly faster streams than seen at PSP. Here we study via numerical MHD simulations the evolution of such large amplitude Alfv\'enic fluctuations by including, in agreement with observations, both a radial magnetic field inversion and an initially constant total magnetic pressure. Despite the extremely large excursion of magnetic and velocity fields, switchbacks are seen to persist for up to hundreds of Alfv\'en crossing times before eventually decaying due to the parametric decay instability. Our results suggest that such switchback/jet configurations might indeed originate in the lower corona and survive out to PSP distances, provided the background solar wind is sufficiently calm, in the sense of not being pervaded by strong density fluctuations or other gradients, such as stream or magnetic field shears, that might destabilize or destroy them over shorter timescales.

We use 2.5D Magnetohydrodynamic simulations to investigate the spectral signatures of the nonlinear disruption of a tearing unstable current sheet via the generation of multiple secondary current sheets and magnetic islands. During the nonlinear phase of tearing mode evolution, there develops a regime in which the magnetic energy density shows a spectrum with a power-law close to $B(k)^2\sim k^{-0.8}$. Such an energy spectrum is found in correspondence of the neutral line, within the diffusion region of the primary current sheet, where energy is conveyed towards smaller scales via a ``recursive'' process of fast tearing-type instabilities. Far from the neutral line we find that magnetic energy spectra evolve towards slopes compatible with the ``standard'' Kolmogorov spectrum. Starting from a self-similar description of the nonlinear stage at the neutral line, we provide a model that predicts a reconnecting magnetic field energy spectrum scaling as $k^{-4/5}$, in good agreement with numerical results. An extension of the predicted power-law to generic current sheet profiles is also given and possible implications for turbulence phenomenology are discussed. These results provide a step forward to understand the ``recursive'' generation of magnetic islands (plasmoids), which has been proposed as a possible explanation for the energy release during flares, but which, more in general, can have an impact on the subsequent turbulent evolution of unstable sheets that naturally form in the high-Lundquist number and collisionless plasmas found in most of the astrophysical environments.

Magnetohydrodynamic (MHD) theory and simulations have shown that reconnection is triggered via a fast "ideal" tearing instability in current sheets whose inverse aspect ratio decreases to $a/L\sim S^{-1/3}$, with $S$ is the Lundquist number defined by the half-length $L$ of the current sheet (of thickness $2a$). Ideal tearing, in 2D sheets, triggers a hierarchical collapse via stretching of X-points and recursive instability. At each step, the local Lundquist number decreases, until the subsequent sheet thickness either approaches kinetic scales or the Lundquist number becomes sufficiently small. Here we carry out a series of Hall-MHD simulations to show how the Hall effect modifies recursive reconnection once the ion inertial scale is approached. We show that as the ion inertial length becomes of the order of the inner, singular layer thickness at some step of the recursive collapse, reconnection transits from the plasmoid-dominant regime to an intermediate plasmoid+Hall regime and then to the Hall-dominant regime. The structure around the X-point, the reconnection rate, the dissipation property and the power spectra are also modified significantly by the Hall effect.

Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal' - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.

Magnetohydrodynamic simulations have shown that a non-unique critical Lundquist number $S_c$ exists, hovering around $S_c \sim 10^4$, above which threshold Sweet-Parker type stationary reconnecting configurations become unstable to a fast tearing mode dominated by plasmoid generation. It is known that the flow along the sheet plays a stabilizing role, though a satisfactory explanation of the non-universality and variable critical Lundquist numbers observed is still lacking. Here we discuss this question using 2D linear MHD simulations and linear stability analyses of Sweet-Parker type current sheets in the presence of background stationary inflows and outflows at low Lundquist numbers ($S\le 10^4$). Simulations show that the inhomogeneous outflow stabilizes the current sheet by stretching the growing magnetic islands and at the same time evacuating the magnetic islands out of the current sheet. This limits the time during which fluctuations which begin at any given wave-length can remain unstable, rendering the instability non-exponential. We find that the linear theory based on the expanding-wavelength assumption works well for $S$ larger than $\sim 1000$. However we also find that the inflow and location of the initial perturbation also affect the stability threshold.

In this paper we study the scaling relations for the triggering of the fast, or "ideal", tearing instability starting from equilibrium configurations relevant to astrophysical as well as laboratory plasmas that differ from the simple Harris current sheet configuration. We present the linear tearing instability analysis for equilibrium magnetic fields which a) go to zero at the boundary of the domain and b) contain a double current sheet system (the latter previously studied as a cartesian proxy for the m=1 kink mode in cylindrical plasmas). More generally, we discuss the critical aspect ratio scalings at which the growth rates become independent of the Lundquist number $S$, in terms of the dependence of the $\Delta'$ parameter on the wavenumber $k$ of unstable modes. The scaling $\Delta'(k)$ with $k$ at small $k$ is found to categorize different equilibria broadly: the critical aspect ratios may be even smaller than $L/a \sim S^{\alpha}$ with $\alpha=1/3$ originally found for the Harris current sheet, but there exists a general lower bound $\alpha\ge1/4$.

Alfv\'enic fluctuations in the solar wind display many properties reflecting an ongoing nonlinear cascade, e.g. a well-defined spectrum in frequency, together with some characteristics more commonly associated with the linear propagation of waves from the Sun, such as the variation of fluctuation amplitude with distance, dominated by solar wind expansion effects. Therefore both nonlinearities and expansion must be included simultaneously in any successful model of solar wind turbulence evolution. Because of the disparate spatial scales involved, direct numerical simulations of turbulence in the solar wind represent an arduous task, especially if one wants to go beyond the incompressible approximation. Indeed, most simulations neglect solar wind expansion effects entirely. Here we develop a numerical model to simulate turbulent fluctuations from the outer corona to 1 AU and beyond, including the sub-Alfv\'enic corona. The accelerating expanding box (AEB) extends the validity of previous expanding box models by taking into account both the acceleration of the solar wind and the inhomogeneity of background density and magnetic field. Our method incorporates a background accelerating wind within a magnetic field that naturally follows the Parker spiral evolution using a two-scale analysis in which the macroscopic spatial effect coupling fluctuations with background gradients becomes a time-dependent coupling term in a homogeneous box. In this paper we describe the AEB model in detail and discuss its main properties, illustrating its validity by studying Alfv\'en wave propagation across the Alfv\'en critical point.

One of the main questions in magnetic reconnection is the origin of triggering behavior with on/off properties that accounts, once it is activated, for the fast magnetic energy conversion to kinetic and thermal energies at the heart of explosive events in astrophysical and laboratory plasmas. Over the past decade progress has been made on the initiation of fast reconnection via the plasmoid instability and what has been called "ideal" tearing, which sets in once current sheets thin to a critical inverse aspect ratio $(a/L)_c$: as shown by Pucci and Velli (2014), at $(a/L)_c \sim S^{-1/3}$ the time scale for the instability to develop becomes of the order of the Alfv\'en time and independent of the Lundquist number (here defined in terms of current sheet length $L$). However, given the large values of $S$ in natural plasmas, this transition might occur for thicknesses of the inner resistive singular layer which are comparable to the ion inertial length $d_i$. When this occurs, Hall currents produce a three-dimensional quadrupole structure of magnetic field, and the dispersive waves introduced by the Hall effect accelerate the instability. Here we present a linear study showing how the "ideal" tearing mode critical aspect ratio is modified when Hall effects are taken into account, including more general scaling laws of the growth rates in terms of sheet inverse aspect ratio: the critical inverse aspect ratio is amended to $a/L \simeq (di/L)^ {0.29} (1/S)^{0.19}$, at which point the instability growth rate becomes Alfv\'enic and does not depend on either of the (small) parameters $d_i/L, 1/S$. We discuss the implications of this generalized triggering aspect ratio for recently developed phase diagrams of magnetic reconnection.

Magnetic reconnection in laboratory, space and astrophysical plasmas is often invoked to explain explosive energy release and particle acceleration. However, the timescales involved in classical models within the macroscopic MHD regime are far too slow to match the observations. Here we revisit the tearing instability by performing visco-resistive two-dimensional numerical simulations of the evolution of thin current sheets, for a variety of initial configurations and of values of the Lunquist number $S$, up to $10^7$. Results confirm that when the critical aspect ratio of $S^{1/3}$ is reached in the reconnecting current sheets, the instability proceeds on ideal (Alfv\'enic) macroscopic timescales, as required to explain observations. Moreover, the same scaling is seen to apply also to the local, secondary reconnection events triggered during the nonlinear phase of the tearing instability, thus accelerating the cascading process to increasingly smaller spatial and temporal scales. The process appears to be robust, as the predicted scaling is measured both in inviscid simulations and when using a Prandtl number $P=1$ in the viscous regime.

Magnetic reconnection is thought to be the dynamical mechanism underlying many explosive phenomena observed both in space and in the laboratory, though the question of how fast magnetic reconnection is triggered in such high Lundquist ($S$) number plasmas has remained elusive. It has been well established that reconnection can develop over timescales faster than those predicted traditionally once kinetic scales are reached. It has also been shown that, within the framework of resistive Magnetohydrodynamics (MHD), fast reconnection is achieved for thin enough sheets via the onset of the so-called plasmoid instability. The latter was discovered in studies specifically devoted to the Sweet-Parker current sheet, either as an initial condition or an apparent transient state developing in nonlinear studies. On the other hand, a fast tearing instability can grow on an ideal, i.e., $S$-independent, timescale (dubbed "ideal" tearing) within current sheets whose aspect ratio scales with the macroscopic Lundquist number as $L/a\sim S^{1/3}$ -- much smaller than the Sweet-Parker one -- suggesting a new way to approach to the initiation of fast reconnection in collapsing current configurations. Here we present an overview of what we have called "ideal" tearing in resistive MHD, and discuss how the same reasoning can be extended to other plasma models commonly used that include electron inertia and kinetic effects. We then discuss a scenario for the onset of "ideal" fast reconnection via collapsing current sheets and describe a quantitative model for the interpretation of the nonlinear evolution of "ideally" unstable sheets in two dimensions.