Activation of MHD reconnection on ideal timescales

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.Comment: Accepted for publication in Plasma Physics and Controlled Fusio

Onset and Evolution of the Oblique, Resonant Electron Firehose Instability in the Expanding Solar Wind Plasma

A double adiabatically expanding solar wind would quickly develop large parallel to perpendicular temperature anisotropies in electrons and ions that are not observed. One reason is that firehose instabilities would be triggered, leading to an ongoing driving/saturation evolution mechanism. We verify this assumption here for the first time for the electron distribution function and the electron firehose instability (EFI), using fully kinetic simulations with the Expanding Box Model. This allows the self-consistent study of onset and evolution of the oblique, resonant EFI in an expanding solar wind. We characterize how the competition between EFI and adiabatic expansion plays out in high- A nd low-beta cases, in high- A nd low-speed solar wind streams. We observe that, even when competing against expansion, the EFI results in perpendicular heating and parallel cooling. These two concurrent processes effectively limit the expansion-induced increase in temperature anisotropy and parallel electron beta. We show that the EFI goes through cycles of stabilization and destabilization: When higher wave number EFI modes saturate, lower wave number modes are destabilized by the effects of the expansion. We show how resonant wave/particle interaction modifies the electron velocity distribution function after the onset of the EFI. The simulations are performed with the fully kinetic, semi-implicit expanding box code EB-iPic3D

Particle-in-cell simulations of Alfv\'en wave parametric decay in a low-beta plasma

We study the parametric decay instability of parallel propagating Alfv\'en wave in a low-beta plasma using one-dimensional fully kinetic simulations. We focus for the first time on the conversion of the energy stored in the initial Alfv\'en wave into particle internal energy, and on its partition between particle species. We show that compressible fluctuations generated by the decay of the pump wave into a secondary ion-acoustic mode and a reflected Alfv\'en wave contribute to the gain of internal energy via two distinct mechanisms. First, the ion-acoustic mode leads nonlinearly to proton trapping and proton phase space mixing, in agreement with previous work based on hybrid simulations. Second, during the nonlinear stage, a compressible front of the fast type develops at the steepened edge of the backward Alfv\'en wave leading to a field-aligned proton beam propagating backwards at the Alfv\'en speed. We find that parametric decay heats preferentially protons, which gain about 50% of the pump wave energy in the form of internal energy. However, we find that electrons are also energized and that they contribute to the total energy balance by gaining 10% of the pump wave energy. By investigating energy partition and particle heating during parametric decay, our results contribute to determine the role of compressible and kinetic effects in wave-driven models of the solar wind

The role of parametric instabilities in turbulence generation and proton heating: Hybrid simulations of parallel propagating Alfv\'en waves

Large amplitude Alfv\'en waves tend to be unstable to parametric instabilities which result in a decay process of the initial wave into different daughter waves depending upon the amplitude of the fluctuations and the plasma beta. The propagation angle with respect to the mean magnetic field of the daughter waves plays an important role in determining the type of decay. In this paper, we revisit this problem by means of multi-dimensional hybrid simulations. In particular, we study the decay and the subsequent nonlinear evolution of large-amplitude Alfv\'en waves by investigating the saturation mechanism of the instability and its final nonlinear state reached for different wave amplitudes and plasma beta conditions. As opposed to one-dimensional simulations where the Decay instability is suppressed for increasing plasma beta values, we find that the decay process in multi-dimensions persists at large values of the plasma beta via the filamentation/magnetosonic decay instabilities. In general, the decay process acts as a trigger both to develop a perpendicular turbulent cascade and to enhance mean field-aligned wave-particle interactions. We find indeed that the saturated state is characterized by a turbulent plasma displaying a field-aligned beam at the Alfv\'en speed and increased temperatures that we ascribe to the Landau resonance and pitch angle scattering in phase space

Local proton heating at magnetic discontinuities in Alfvenic and non-Alfvenic solar wind

We investigate the local proton energization at magnetic discontinuities/intermittent structures and the corresponding kinetic signatures in velocity phase space in Alfv\'enic and non-Alfv\'enic wind streams observed by Parker Solar Probe. By means of the Partial Variance of Increments method, we find that the hottest proton populations are localized around compressible, kinetic-scale magnetic structures in both types of wind. Furthermore, the Alfv\'enic wind shows preferential enhancements of $T_\parallel$ as smaller scale structures are considered, whereas the non-Alfvenic wind shows preferential $T_\bot$ enhancements. Although proton beams are present in both types of wind, the proton velocity distribution function displays distinct features. Hot beams, i.e., beams with beam-to-core perpendicular temperature up to three times larger than the total distribution anisotropy, are found in the non-Alfv\'enic wind, whereas colder beams in the Alfv\'enic wind. Our data analysis is complemented by 2.5D hybrid simulations in different geometrical setups, which support the idea that proton beams in Alfv\'enic and non-Alfv\'enic wind have different kinetic properties and origins. The development of a perpendicular nonlinear cascade, favored in balanced turbulence, allows a preferential relative enhancement of the perpendicular plasma temperature and the formation of hot beams. Cold field-aligned beams are instead favored by Alfv\'en wave steepening. Non-Maxwellian distribution functions are found near discontinuities and intermittent structures, pointing to the fact that the nonlinear formation of small-scale structures is intrinsically related to the development of highly non-thermal features in collisionless plasmas

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.Comment: Improved version, accepted to JPP Lecture Notes. Some parts were shortened and some parts were expanded. The text now contains Conclusion

First Solar Orbiter observation of the Alfvénic slow wind and identification of its solar source

Context: Turbulence dominated by large amplitude nonlinear Alfvén-like fluctuations mainly propagating away from the Sun is ubiquitous in high speed solar wind streams. Recent studies have shown that also slow wind streams may show strong Alfvénic signatures, especially in the inner heliosphere. Aims: The present study focuses on the characterisation of an Alfvénic slow solar wind interval observed by Solar Orbiter on July 14-18, 2020 at a heliocentric distance of 0.64 AU. Methods: Our analysis is based on plasma moments and magnetic field measurements from the Solar Wind Analyser (SWA) and Magnetometer (MAG) instruments, respectively. We compare the behaviour of different parameters to characterise the stream in terms of the Alfvénic content and magnetic properties. We perform also a spectral analysis to highlight spectral features and waves signature using power spectral density and magnetic helicity spectrograms, respectively. Moreover, we reconstruct the Solar Orbiter magnetic connectivity to the solar sources via both a ballistic and a Potential Field Source Surface (PFSS) model. Results: The Alfvénic slow wind stream described in this paper resembles in many respects a fast wind stream. Indeed, at large scales, the time series of the speed profile shows a compression region, a main portion of the stream and a rarefaction region, characterised by different features. Moreover, before the rarefaction region, we pinpoint several structures at different scales recalling the spaghetti-like flux-tube texture of the interplanetary magnetic field. Finally, we identify the connections between Solar Orbiter in situ measurements, tracing them down to coronal streamer and pseudostreamer configurations. Conclusions. The characterisation of the Alfvénic slow wind stream observed by Solar Orbiter and the identification of its solar source are extremely important aspects to understand possible future observations of the same solar wind regime, especially as solar activity is increasing toward a maximum, where a higher incidence of this solar wind regime is expected

Magnetic field spectral evolution in the inner heliosphere

Parker Solar Probe and Solar Orbiter data are used to investigate the radial evolution of magnetic turbulence between $0.06 ~ \lesssim R ~\lesssim 1$ au. The spectrum is studied as a function of scale, normalized to the ion inertial scale $d_{i}$. In the vicinity of the Sun, the inertial range is limited to a narrow range of scales and exhibits a power-law exponent of, $\alpha_{B} = -3/2$, independent of plasma parameters. The inertial range grows with distance, progressively extending to larger spatial scales, while steepening towards a $\alpha_{B} =-5/3$ scaling. It is observed that spectra for intervals with large magnetic energy excesses and low Alfv\'enic content steepen significantly with distance, in contrast to highly Alfv\'enic intervals that retain their near-Sun scaling. The occurrence of steeper spectra in slower wind streams may be attributed to the observed positive correlation between solar wind speed and Alfv\'enicity.Comment: Accepted to APJ letters with minor revision