A new path for the generation of a sub-ion scale cascade in collisionless space and astrophysical plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of powerlaw energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ∼ −2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of −5/3 in the MHD range and of ∼ −3 at sub-ion scales.

We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of twodimensional numerical simulations, within the framework of compressible, resistive MHD. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number S as S 1 3-. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth that are independent of S provided this is high enough, as it should be natural having the ideal case as a limit for S  ¥. Our simulations confirm that the tearing instability growth rate can be as fast aswhere A t is the ideal Alfvénic time set by the macroscopic scales, for our least diffusive case with S 10 7 = . The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of S explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with S, here calculated on the local, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin, and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvénic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.

Context. Local helioseismology has so far relied on semi-analytical methods to compute the spatial sensitivity of wave travel times to perturbations in the solar interior. These methods are cumbersome and lack flexibility. Aims. Here we propose a convenient framework for numerically solving the forward problem of time-distance helioseismology in the frequency domain. The fundamental quantity to be computed is the cross-covariance of the seismic wavefield. Methods. We choose sources of wave excitation that enable us to relate the cross-covariance of the oscillations to the Green's function in a straightforward manner. We illustrate the method by considering the 3D acoustic wave equation in an axisymmetric reference solar model, ignoring the effects of gravity on the waves. The symmetry of the background model around the rotation axis implies that the Green's function can be written as a sum of longitudinal Fourier modes, leading to a set of independent 2D problems. We use a high-order finite-element method to solve the 2D wave equation in frequency space. The computation is 'embarrassingly parallel', with each frequency and each azimuthal order solved independently on a computer cluster. Results. We compute travel-time sensitivity kernels in spherical geometry for flows, sound speed, and density perturbations under the first Born approximation. Convergence tests show that travel times can be computed with a numerical precision better than one millisecond, as required by the most precise travel-time measurements. Conclusions. The method presented here is computationally efficient and will be used to interpret travel-time measurements in order to infer, e.g., the large-scale meridional flow in the solar convection zone. It allows the implementation of (full-waveform) iterative inversions, whereby the axisymmetric background model is updated at each iteration.

We investigate the properties of plasma turbulence by means of two-dimensional Hallmagnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations. We find that HMHD simulations exhibit spectral properties that are in most cases in agreement with the results of the HPIC simulations and with solar wind observations. The energy spectra of magnetic fluctuations exhibit a double power-law with spectral index −5/3 at MHD scales and −3 at kinetic scales, while for velocity fluctuations the spectral index is −3/2 at MHD scales. The break between the MHD and the kinetic scales occurs at the same scale in both simulations. In the MHD range the slopes of the total energy and residual energy spectra satisfy a fast Alfvén-dynamo balance. The development of a turbulent cascade is concurrently characterized by magnetic reconnection events taking place in thin current sheets that form between large eddies. A statistical analysis reveals that reconnection is qualitatively the same and fast in both the HMHD and HPIC models, characterized by inverse reconnection rates much smaller than the characteristic large-eddy nonlinear time. The agreement extends to other statistical properties, such us the kurtosis of the magnetic field. Moreover, the observation of a direct energy transfer from the large vortices to the small sub-ion scales triggered by magnetic reconnection, further supports the existence of a reconnection-mediated turbulent regime at kinetic scales. We conclude that the Hall-MHD fluid description captures to a large extent the transition of the turbulent cascade between the large MHD scales and the sub-ion scales.

Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequent nonlinear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics, as appropriate in situations where the Alfvén velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S −1/2 , with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S −1/3 are considered, the so-called ideal tearing regime is retrieved, with modes growing independently on S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfvén velocity. In the fully nonlinear stage we observe an inverse cascade towards the fundamental mode, with Petschektype supersonic jets propagating at the external Alfvén speed from the X-point, and a fast reconnection rate at the predicted value R ∼ (ln S) −1 .

We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particlein-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.

We investigate the transition of the solar wind turbulent cascade from MHD to sub‐ion range by means of a detailed comparison between in situ observations and hybrid numerical simulations. In particular, we focus on the properties of the magnetic field and its component anisotropy in Cluster measurements and hybrid 2D simulations. First, we address the angular distribution of wave vector in the kinetic range between ion and electron scales by studying the variance anisotropy of the magnetic field components. When taking into account a single-direction sampling, like that performed by spacecraft in the solar wind, the main properties of the fluctuations observed in situ are also recovered in our numerical description. This result confirms that solar wind turbulence in the sub‐ion range is characterized by a quasi-2D gyrotropic distribution of k-vectors around the mean field. We then consider the magnetic compressibility associated with the turbulent cascade and its evolution from large-MHD to sub‐ion scales. The ratio of field aligned to perpendicular fluctuations, typically low in the MHD inertial range, increases significantly when crossing ion scales and its value in the sub‐ion range is a function of the total plasma beta only, as expected from theoretical predictions, with higher magnetic compressibility for higher beta. Moreover, we observe that this increase has a gradual trend from low to high beta values in the in situ data; this behavior is well captured by the numerical simulations. The level of magnetic field compressibility that is observed in situ and in the simulations is in fairly good agreement with theoretical predictions, especially at high beta, suggesting that, in the kinetic range explored, the turbulence is supported by low-frequency and highly oblique fluctuations in pressure balance, like kinetic Alfvén waves or other slowly evolving coherent structures. The resulting scaling properties as a function of the plasma beta and the main differences between numerical and theoretical expectations and in situ observations are also discussed.

No abstract