Showing 1–44 of 44 results for author: Hellinger, P

Anisotropic electron heating in turbulence-driven magnetic reconnection in the near-Sun solar wind

Authors: Luca Franci, Emanuele Papini, Alfredo Micera, Giovanni Lapenta, Petr Hellinger, Daniele Del Sarto, David Burgess, Simone Landi

Abstract: We perform a high-resolution two-dimensional fully-kinetic numerical simulation of a turbulent plasma system with observation-driven conditions, in order to investigate the interplay between turbulence, magnetic reconnection, and particle heating from ion to sub-electron scales in the near-Sun solar wind. We find that the power spectra of the turbulent plasma and electromagnetic fluctuations show… ▽ More We perform a high-resolution two-dimensional fully-kinetic numerical simulation of a turbulent plasma system with observation-driven conditions, in order to investigate the interplay between turbulence, magnetic reconnection, and particle heating from ion to sub-electron scales in the near-Sun solar wind. We find that the power spectra of the turbulent plasma and electromagnetic fluctuations show multiple power-law intervals down to scales smaller than the electron gyroradius. Magnetic reconnection is observed to occur in correspondence of current sheets with a thickness of the order of the electron inertial length, which form and shrink due to interacting ion-scale vortexes. In some cases, both ion and electron outflows are observed (the classic reconnection scenario), while in others -- typically for the shortest current sheets -- only electron jets are presents ("electron-only reconnection"). At the onset of reconnection, the electron temperature starts to increase and a strong parallel temperature anisotropy develops. This suggests that in strong turbulence electron-scale coherent structures may play a significant role for electron heating, as impulsive and localized phenomena such as magnetic reconnection may transfer energy from the electromagnetic fields to particles more efficiently than damping mechanisms related to interactions with wave-like fluctuations. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 18 pages, 6 figures

Ion-scale transition of plasma turbulence: Pressure-strain effect

Authors: Petr Hellinger, Victor Montagud-Camps, Luca Franci, Lorenzo Matteini, Emanuele Papini, Andrea Verdini, Simone Landi

Abstract: We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (ex… ▽ More We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (extended by considering the plasma pressure as a tensor quantity) to the numerical results. The KHM analysis indicates that the transition from MHD to ion scales (the so called ion break in the power spectrum) results from a combination of an onset of Hall physics and of an effective dissipation owing to the pressure-strain energy-exchange channel and resistivity. We discuss the simulation results in the context of the solar wind. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 7 pages, 3 figures

Spectrum of kinetic plasma turbulence at 0.3-0.9 astronomical units from the Sun

Authors: Olga Alexandrova, Vamsee Krishna Jagarlamudi, Petr Hellinger, Milan Maksimovic, Yuri Shprits, Andre Mangeney

Abstract: We investigate spectral properties of turbulence in the solar wind that is a weakly collisional astrophysical plasma, accessible to in-situ observations. Using the Helios search coil magnetometer measurements in the fast solar wind, in the inner heliosphere, we focus on properties of the turbulent magnetic fluctuations at scales smaller than the ion characteristic scales, the so-called kinetic pla… ▽ More We investigate spectral properties of turbulence in the solar wind that is a weakly collisional astrophysical plasma, accessible to in-situ observations. Using the Helios search coil magnetometer measurements in the fast solar wind, in the inner heliosphere, we focus on properties of the turbulent magnetic fluctuations at scales smaller than the ion characteristic scales, the so-called kinetic plasma turbulence. At such small scales, we show that the magnetic power spectra between 0.3 and 0.9 AU from the Sun have a generic shape $\sim f^{-8/3}\exp{(-f/f_d)}$ where the dissipation frequency $f_d$ is correlated with the Doppler shifted frequency $f_{ρe}$ of the electron Larmor radius. This behavior is statistically significant: all the observed kinetic spectra are well described by this model, with $f_d = f_{ρe}/1.8$. Our results indicate that the electron gyroradius plays the role of the dissipation scale and marks the end of the electromagnetic cascade in the solar wind. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Comments: Published in PhysRevE in June 2021 https://ui.adsabs.harvard.edu/abs/2021PhRvE.103f3202A/abstract. arXiv admin note: substantial text overlap with arXiv:2004.01102

Journal ref: Physical Review E, Volume 103, Issue 6, article id.063202, June 2021

Spacetime Hall-MHD turbulence at sub-ion scales: structures or waves?

Authors: Emanuele Papini, Antonio Cicone, Luca Franci, Mirko Piersanti, Simone Landi, Petr Hellinger, Andrea Verdini

Abstract: Spatiotemporal properties of two-dimensional (2D) Hall-magnetohydrodynamic turbulence at intermediate plasma $β=2$ are studied by means of Fast Iterative Filtering, a new technique for the decomposition of nonstationary nonlinear signals. Results show that the magnetic energy at sub-ion scales is concentrated in perturbations with frequencies smaller than the ion-cyclotron (IC) frequency and with… ▽ More Spatiotemporal properties of two-dimensional (2D) Hall-magnetohydrodynamic turbulence at intermediate plasma $β=2$ are studied by means of Fast Iterative Filtering, a new technique for the decomposition of nonstationary nonlinear signals. Results show that the magnetic energy at sub-ion scales is concentrated in perturbations with frequencies smaller than the ion-cyclotron (IC) frequency and with polarization properties that are incompatible with both kinetic Alfvén waves (KAWs) and IC waves. At higher frequencies, we clearly identify signatures of both whistler waves and KAWs, however their energetic contribution to the magnetic power spectrum is negligible. We conclude that the dynamics of 2D Hall-MHD turbulence at sub-ion scales is mainly driven by localized intermittent structures, with no significant contribution of wavelike fluctuations. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 10 pages, 5 figures. Accepted for publication on The Astrophysical Journal Letters

Spectral transfer and Kármán-Howarth-Monin equations for compressible Hall magnetohydrodynamics

Authors: Petr Hellinger, Emanuele Papini, Andrea Verdini, Simone Landi, Luca Franci, Lorenzo Matteini, Victor Montagud-Camps

Abstract: We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to… ▽ More We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to the periodic boundary conditions and have complementary cumulative behavior. They are used here to analyze the onset of turbulence and its properties when it is fully developed. These approaches give equivalent results characterizing: the decay of the kinetic + magnetic energy at large scales, the MHD and Hall cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. The Hall cascade appears when the MHD one brings the energy close to the ion inertial range and is related to the formation of reconnecting current sheets. At later times, the pressure-dilation energy-exchange rate oscillates around zero with no net effect on the cross-scale energy transfer when averaged over a period of its oscillations. A reduced one-dimensional analysis suggests that all three methods may be useful to estimate the energy cascade rate from in situ observations. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 11 pages, 7 figues

Proton energization by phase-steepening of parallel propagating Alfvénic fluctuations

Authors: C. A. González, A. Tenerani, L. Matteini, P. Hellinger, M. Velli

Abstract: Proton energization at magnetic discontinuities generated by phase-steepened fronts of parallel propagating, large-amplitude Alfvénic fluctuation is studied using hybrid simulations. We find that dispersive effects yield to the collapse of the wave via phase steepening and the subsequent generation of compressible fluctuations that mediate an efficient local energy transfer from the wave to the pr… ▽ More Proton energization at magnetic discontinuities generated by phase-steepened fronts of parallel propagating, large-amplitude Alfvénic fluctuation is studied using hybrid simulations. We find that dispersive effects yield to the collapse of the wave via phase steepening and the subsequent generation of compressible fluctuations that mediate an efficient local energy transfer from the wave to the protons. Proton scattering at the steepened edges causes non-adiabatic proton perpendicular heating. Furthermore, the parallel electric field at the propagating fronts mediates the acceleration of protons along the mean field. A steady-state is achieved where proton distribution function displays a field-aligned beam at the Alfvén speed, and compressible fluctuations are largely damped. We discuss the implications of our results in the context of Alfvénic solar wind. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Scale dependence and cross-scale transfer of kinetic energy in compressible hydrodynamic turbulence at moderate Reynolds numbers

Authors: P. Hellinger, A. Verdini, S. Landi, E. Papini, L. Franci, L. Matteini

Abstract: We investigate properties of the scale dependence and cross-scale transfer of kinetic energy in compressible three-dimensional hydrodynamic turbulence, by means of two direct numerical simulations of decaying turbulence with initial Mach numbers M = 1/3 and M = 1, and with moderate Reynolds numbers, R_lambda ~ 100. The turbulent dynamics is analyzed using compressible and incompressible versions o… ▽ More We investigate properties of the scale dependence and cross-scale transfer of kinetic energy in compressible three-dimensional hydrodynamic turbulence, by means of two direct numerical simulations of decaying turbulence with initial Mach numbers M = 1/3 and M = 1, and with moderate Reynolds numbers, R_lambda ~ 100. The turbulent dynamics is analyzed using compressible and incompressible versions of the dynamic spectral transfer (ST) and the Karman-Howarth-Monin (KHM) equations. We find that the nonlinear coupling leads to a flux of the kinetic energy to small scales where it is dissipated; at the same time, the reversible pressure-dilatation mechanism causes oscillatory exchanges between the kinetic and internal energies with an average zero net energy transfer. While the incompressible KHM and ST equations are not generally valid in the simulations, their compressible counterparts are well satisfied and describe, in a quantitatively similar way, the decay of the kinetic energy on large scales, the cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. There exists a simple relationship between the KHM and ST results through the inverse proportionality between the wave vector k and the spatial separation length l as k l ~ 3^1/2. For a given time the dissipation and pressure-dilatation terms are strong on large scales in the KHM approach whereas the ST terms become dominant on small scales; this is owing to the complementary cumulative behavior of the two methods. The effect of pressure dilatation is weak when averaged over a period of its oscillations and may lead to a transfer of the kinetic energy from large to small scales without a net exchange between the kinetic and internal energies. Our results suggest that for large-enough systems there exists an inertial range for the kinetic energy cascade ... △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 14 pages, 10 figures

Evidence of a "current-mediated" turbulent regime in space and astrophysical plasmas

Authors: Luca Franci, Daniele Del Sarto, Emanuele Papini, Alice Giroul, Julia E. Stawarz, David Burgess, Petr Hellinger, Simone Landi, Stuart D. Bale

Abstract: How the turbulent energy cascade develops below the magnetohydrodynamic scales in space and astrophysical plasmas is a major open question. Here, we measure the power spectrum of magnetic fluctuations in Parker Solar Probe's observations close to the Sun and in state-of-the-art numerical simulations of plasma turbulence. Both reveal a power-law behavior with a slope compatible with $-11/3$ at scal… ▽ More How the turbulent energy cascade develops below the magnetohydrodynamic scales in space and astrophysical plasmas is a major open question. Here, we measure the power spectrum of magnetic fluctuations in Parker Solar Probe's observations close to the Sun and in state-of-the-art numerical simulations of plasma turbulence. Both reveal a power-law behavior with a slope compatible with $-11/3$ at scales smaller than the ion characteristic scales, steeper than what is typically observed in the solar wind and in the Earth's magnetosheath. We explain such behavior by developing a simple two-fluid model which does not require any kinetic processes nor electron-inertia effects. This is characterized by a significant contribution of the ion kinetic energy to the total turbulent energy cascade at sub-ion scales, although the dynamics is driven by the magnetic field through the current density. We expect that this regime may be relevant for a broad class of low-beta plasmas, e.g. the solar corona, non-relativistic magnetized jets and disks, and laboratory plasmas. △ Less

Submitted 10 October, 2020; originally announced October 2020.

Magnetic field turbulence in the solar wind at sub-ion scales: in situ observations and numerical simulations

Authors: L. Matteini, L. Franci, O. Alexandrova, C. Lacombe, S. Landi, P. Hellinger, E. Papini, A. Verdini

Abstract: We investigate the transition of the solar wind turbulent cascade from MHD to sub-ion range by means of a detail 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-vectors in the kine… ▽ More We investigate the transition of the solar wind turbulent cascade from MHD to sub-ion range by means of a detail 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-vectors in the kinetic range between ion and electron scales by studying the variance anisotropy of the magnetic field components. When taking into account the single-direction sampling 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 behaviour 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. △ Less

Submitted 30 August, 2020; originally announced August 2020.

Comments: Manuscript submitted to Frontiers Astronomy and Space Sciences, Research Topic: Improving the Understanding of Kinetic Processes in Solar Wind and Magnetosphere: From CLUSTER to MMS

The role of parametric instabilities in turbulence generation and proton heating: Hybrid simulations of parallel propagating Alfvén waves

Authors: C. A. González, A. Tenerani, M. Velli, P. Hellinger

Abstract: Large amplitude Alfvén 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 thi… ▽ More Large amplitude Alfvén 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én 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én speed and increased temperatures that we ascribe to the Landau resonance and pitch angle scattering in phase space. △ Less

Submitted 28 September, 2020; v1 submitted 21 June, 2020; originally announced June 2020.

Multidimensional Iterative Filtering: a new approach for investigating plasma turbulence in numerical simulations

Authors: Emanuele Papini, Antonio Cicone, Mirko Piersanti, Luca Franci, Petr Hellinger, Simone Landi, Andrea Verdini

Abstract: Turbulent space and astrophysical plasmas exhibit a complex dynamics, which involves nonlinear coupling across different temporal and spatial scales. There is growing evidence that impulsive events, such as magnetic reconnection instabilities, lead to a spatially localized enhancement of energy dissipation, thus speeding up the energy transfer at small scales. Capturing such a diverse dynamics is… ▽ More Turbulent space and astrophysical plasmas exhibit a complex dynamics, which involves nonlinear coupling across different temporal and spatial scales. There is growing evidence that impulsive events, such as magnetic reconnection instabilities, lead to a spatially localized enhancement of energy dissipation, thus speeding up the energy transfer at small scales. Capturing such a diverse dynamics is challenging. Here, we employ the Multidimensional Iterative Filtering (MIF) method, a novel technique for the analysis of nonstationary multidimensional signals. Unlike other traditional methods (e.g., based on Fourier or wavelet decomposition), MIF does not require any previous assumption on the functional form of the signal to be identified. Using MIF, we carry out a multiscale analysis of Hall-magnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations of decaying plasma turbulence. The results assess the ability of MIF to spatially identify and separate the different scales (the MHD inertial range, the sub-ion kinetic, and the dissipation scales) of the plasma dynamics. Furthermore, MIF decomposition allows to detect localized current structures and to characterize their contribution to the statistical and spectral properties of turbulence. Overall, MIF arises as a very promising technique for the study of turbulent plasma environments. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 19 pages, 9 figures, submitted to Journal of Plasma Physics

On cascade of kinetic energy in compressible hydrodynamic turbulence

Authors: Petr Hellinger, Andrea Verdini, Simone Landi, Luca Franci, Emanuele Papini, Lorenzo Matteini

Abstract: Properties of the turbulent cascade of kinetic energy are studied using direct numerical simulations of three-dimensional hydrodynamic decaying turbulence with a moderate Reynolds number and the initial Mach number $M=1$. Compressible and incompressible versions of the Karman-Howarth-Monin (KHM) and low-pass filtering/coarse-graining approaches are compared. In the simulation the total energy is w… ▽ More Properties of the turbulent cascade of kinetic energy are studied using direct numerical simulations of three-dimensional hydrodynamic decaying turbulence with a moderate Reynolds number and the initial Mach number $M=1$. Compressible and incompressible versions of the Karman-Howarth-Monin (KHM) and low-pass filtering/coarse-graining approaches are compared. In the simulation the total energy is well conserved; the scale dependent KHM and coarse-grained energy equations are also well conserved; the two approaches show similar results, the system does not have an inertial range for the cascade of kinetic energy, the region where this cascade dominates also have a non-negligible contribution of the kinetic-energy decay, dissipation, and pressure-dilatation effects. While the two approaches give semi-quantitatively similar results for the kinetic energy cascade, dissipation and pressure-dilatation rates, they differ in the increment separation and filtering scales; these scales are not simply related. The two approaches may be used to find the inertial range and to determine the cascade/dissipation rate of the kinetic energy. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 13 pages, 5 figures

Kinetic turbulence in space plasmas observed in the near-Earth and near-Sun solar wind

Authors: Olga Alexandrova, Vamsee Krishna Jagarlamudi, Claudia Rossi, Milan Maksimovic, Petr Hellinger, Yuri Shprits, André Mangeney

Abstract: Turbulence develops in any stressed flow when the scales of the forcing are much larger than those of the dissipation. In neutral fluids, it consists of chaotic motions in physical space but with a universal energy spectrum in Fourier space. Intermittency (non-Gaussian statistics of fluctuations) is another general property and it is related to the presence of coherent structures. Space plasmas ar… ▽ More Turbulence develops in any stressed flow when the scales of the forcing are much larger than those of the dissipation. In neutral fluids, it consists of chaotic motions in physical space but with a universal energy spectrum in Fourier space. Intermittency (non-Gaussian statistics of fluctuations) is another general property and it is related to the presence of coherent structures. Space plasmas are turbulent as well. Here, we focus on the kinetic plasma scales, which are not yet well understood. We address the following fundamental questions: (1) Do the turbulent fluctuations at kinetic scales form a universal spectrum? and (2) What is the nature of the fluctuations? Using measurements in the solar wind we show that the magnetic spectra of kinetic turbulence at 0.3, 0.6 and 0.9 AU from the Sun have the same shape as the ones close to the Earth orbit at 1 AU, indicating universality of the phenomenon. The fluctuations, which form this spectrum, are typically non-linearly interacting eddies that tend to generate magnetic filaments. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 12 pages, 7 figures, submitted to Nature Communications on March 2, 2020

Modeling Kelvin-Helmholtz instability-driven turbulence with hybrid simulations of Alfvénic turbulence

Authors: Luca Franci, Julia E. Stawarz, Emanuele Papini, Petr Hellinger, Takuma Nakamura, David Burgess, Simone Landi, Andrea Verdini, Lorenzo Matteini, Robert Ergun, Olivier Le Contel, Per-Arne Lindqvist

Abstract: Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and… ▽ More Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfvén wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath. △ Less

Submitted 17 November, 2019; originally announced November 2019.

Comments: 11 pages, 6 figures, submitted to The Astrophysical Journal

Fast acceleration of transrelativistic electrons in astrophysical turbulence

Authors: Domenico Trotta, Luca Franci, David Burgess, Petr Hellinger

Abstract: Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the ene… ▽ More Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfvénic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible. △ Less

Submitted 17 April, 2020; v1 submitted 21 October, 2019; originally announced October 2019.

Comments: Accepted for publication in The Astrophysical Journal

Turbulence vs. fire hose instabilities: 3-D hybrid expanding box simulations

Authors: P. Hellinger, L. Matteini, S. Landi, L. Franci, A. Verdini, E. Papini

Abstract: The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-hel… ▽ More The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 12 pages, 10 figures

In-situ observation of Hall Magnetohydrodynamic Cascade in Space Plasma

Authors: Riddhi Bandyopadhyay, Luca Sorriso-Valvo, Alexandros Chasapis, Petr Hellinger, William H. Matthaeus, Andrea Verdini, Simone Landi, Luca Franci, Lorenzo Matteini, Barbara L. Giles, Daniel J. Gershman Craig J. Pollock, Christopher T. Russell, Robert J. Strangeway, Roy B. Torbert, Thomas E. Moore, James L. Burch

Abstract: We present estimates of the turbulent energy cascade rate, derived from a Hall-MHD third-order law. We compute the contribution from the Hall term and the MHD term to the energy flux. We use MMS data accumulated in the magnetosheath and the solar wind, and compare the results with previously established simulation results. We find that in observation, the MHD contribution is dominant at inertial s… ▽ More We present estimates of the turbulent energy cascade rate, derived from a Hall-MHD third-order law. We compute the contribution from the Hall term and the MHD term to the energy flux. We use MMS data accumulated in the magnetosheath and the solar wind, and compare the results with previously established simulation results. We find that in observation, the MHD contribution is dominant at inertial scales, as in the simulations, but the Hall term becomes significant in observations at larger scales than in the simulations. Possible reasons are offered for this unanticipated result. △ Less

Submitted 2 May, 2020; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: Accepted for Publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 124, 225101 (2020)

Spectral anisotropies and intermittency of plasma turbulence at ion kinetic scales

Authors: Simone Landi, Luca Franci, Emanuele Papini, Andrea Verdini, Lorenzo Matteini, Petr Hellinger

Abstract: By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories a… ▽ More By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales. △ Less

Submitted 17 April, 2019; v1 submitted 8 April, 2019; originally announced April 2019.

Comments: 9 pages, 4 figures

Can Hall Magnetohydrodynamics explain plasma turbulence at sub-ion scales?

Authors: Emanuele Papini, Luca Franci, Simone Landi, Andrea Verdini, Lorenzo Matteini, Petr Hellinger

Abstract: We investigate the properties of plasma turbulence by means of two-dimensional Hall-magnetohydrodynamic (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 p… ▽ More We investigate the properties of plasma turbulence by means of two-dimensional Hall-magnetohydrodynamic (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. △ Less

Submitted 23 April, 2019; v1 submitted 4 October, 2018; originally announced October 2018.

Comments: 12 pages, 8 figures

Journal ref: The Astrophysical Journal, 870:52, 2019

Electron mirror instability: Particle-in-cell simulations

Authors: Petr Hellinger, Stepan Stverak

Abstract: Properties of the electron mirror instability and its competition with the usually dominant whistler (electron cyclotron) instability driven by the electron perpendicular temperature anisotropy are investigated on the linear level using a Vlasov linear solver and on the nonlinear level using a two-dimensional full particle code. The simulation results show that the linearly subdominant electron mi… ▽ More Properties of the electron mirror instability and its competition with the usually dominant whistler (electron cyclotron) instability driven by the electron perpendicular temperature anisotropy are investigated on the linear level using a Vlasov linear solver and on the nonlinear level using a two-dimensional full particle code. The simulation results show that the linearly subdominant electron mirror instability may compete on the nonlinear level with the whistler instability and may even become eventually the dominant mode that generates robust non-propagating sub-ion-scale coherent structures in the form of magnetic peaks. △ Less

Submitted 13 June, 2018; originally announced June 2018.

Comments: 12 pages, 7 figures, J. Plasma Phys

Von Karman-Howarth equation for Hall magnetohydrodynamics: Hybrid simulations

Authors: Petr Hellinger, Andrea Verdini, Simone Landi, Luca Franci, Lorenzo Matteini

Abstract: A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is… ▽ More A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: 3 Figures, 5 pages, ApJL, submitted

Density Fluctuations in the Solar Wind Driven by Alfvén Wave Parametric Decay

Authors: Trevor A. Bowen, Samuel Badman, Petr Hellinger, Stuart D. Bale

Abstract: Measurements and simulations of inertial compressive turbulence in the solar wind are characterized by anti-correlated magnetic fluctuations parallel to the mean field and density structures. This signature has been interpreted as observational evidence for non-propagating pressure balanced structures (PBS), kinetic ion acoustic waves, as well as the MHD slow-mode. Given the high damping rates of… ▽ More Measurements and simulations of inertial compressive turbulence in the solar wind are characterized by anti-correlated magnetic fluctuations parallel to the mean field and density structures. This signature has been interpreted as observational evidence for non-propagating pressure balanced structures (PBS), kinetic ion acoustic waves, as well as the MHD slow-mode. Given the high damping rates of parallel propagating compressive fluctuations, their ubiquity in satellite observations is surprising, and suggestive of a local driving process. One possible candidate for the generation of compressive fluctuations in the solar wind is Alfvén wave parametric instability. Here we test the parametric decay process as a source of compressive waves in the solar wind by comparing the collisionless damping rates of compressive fluctuations with the growth rates of the parametric decay instability daughter waves. Our results suggest that generation of compressive waves through parametric decay is overdamped at 1 AU, but that the presence of slow-mode like density fluctuations is correlated with the parametric decay of Alfvén waves. △ Less

Submitted 5 February, 2018; v1 submitted 22 December, 2017; originally announced December 2017.

Three-dimensional simulations of solar wind turbulence with the hybrid code CAMELIA

Authors: L. Franci, P. Hellinger, M. Guarrasi, C. H. K. Chen, E. Papini, A. Verdini, L. Matteini, S. Landi

Abstract: 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 particle-in-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 ag… ▽ More 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 particle-in-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. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 13 pages, 10 figures, submitted to the Proceedings of the 12th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2017)

The parametric instability of Alfvén waves: effects of temperature anisotropy

Authors: Anna Tenerani, Marco Velli, Petr Hellinger

Abstract: We study the stability of large amplitude, circularly polarized Alfvén waves in an anisotropic plasma described by the double-adiabatic/CGL closure, and in particular the effect of a background thermal pressure anisotropy on the well-known properties of Alfvén wave parametric decay in Magnetohydrodynamics (MHD). Anisotropy allows instability over a much wider range of values of parallel plasma bet… ▽ More We study the stability of large amplitude, circularly polarized Alfvén waves in an anisotropic plasma described by the double-adiabatic/CGL closure, and in particular the effect of a background thermal pressure anisotropy on the well-known properties of Alfvén wave parametric decay in Magnetohydrodynamics (MHD). Anisotropy allows instability over a much wider range of values of parallel plasma beta ($β_\parallel$) when $ξ=p_{0\bot}/p_{0\parallel} > 1$. When the pressure anisotropy exceeds a critical value, $ξ\geqξ^*$ with $ξ^*\simeq2.7$, there is a new regime in which the parametric instability is no longer quenched at high $β_\parallel$ and in the limit $β_\parallel\gg1$ the growth rate becomes independent of $β_\parallel$. In the opposite case of $ξ< ξ^*$, the instability is strongly suppressed with increasing parallel plasma beta, similarly to the MHD case. We analyze marginal stability conditions for parametric decay in the $(ξ,\,β_\parallel)$ parameter space, and discuss possible implications for Alfvénic turbulence in the solar wind. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 9 pages, 9 figures. Under review

Solar wind turbulent cascade from MHD to sub-ion scales: large-size 3D hybrid particle-in-cell simulations

Authors: Luca Franci, Simone Landi, Andrea Verdini, Lorenzo Matteini, Petr Hellinger

Abstract: Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic fiel… ▽ More Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, both at magnetohydrodynamic (MHD) and at kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a two-dimensional (2D) simulation with similar physical parameters. This provides a partial, a-posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 11 pages, 8 figures, submitted to The Astrophysical Journal

Magnetic reconnection as a driver for a sub-ion scale cascade in plasma turbulence

Authors: L. Franci, S. S. Cerri, F. Califano, S. Landi, E. Papini, A. Verdini, L. Matteini, F. Jenko, P. Hellinger

Abstract: A new path for the generation of a sub-ion scale cascade in collisionless 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 power-l… ▽ More A new path for the generation of a sub-ion scale cascade in collisionless 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 power-law 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. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 6 pages, 5 figures

Mirror instability in the turbulent solar wind

Authors: P. Hellinger, S. Landi, L. Matteini, A. Verdini, L. Franci

Abstract: The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition betwe… ▽ More The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy. △ Less

Submitted 17 March, 2017; originally announced March 2017.

Comments: ApJ, 7 pages, 9 figures. arXiv admin note: text overlap with arXiv:1508.03159

Plasma turbulence at ion scales: a comparison between PIC and Eulerian hybrid-kinetic approaches

Authors: S. S. Cerri, L. Franci, F. Califano, S. Landi, P. Hellinger

Abstract: Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid.… ▽ More Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning about two decades in wavenumber - from MHD inertial range to scales well below the ion gyroradius - with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell (HPIC) simulations of freely-decaying Alfvénic turbulence, whereas the other consists of Eulerian hybrid Vlasov-Maxwell (HVM) simulations of turbulence continuously driven with partially-compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_\perpρ_i\gtrsim1$ are observed. The system indeed self-consistently "reprocesses" the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter. Small-scale turbulence has been found to be mainly populated by kinetic Alfvén wave (KAW) fluctuations for $β\geq1$, whereas KAW fluctuations are only sub-dominant for low-$β$. △ Less

Submitted 7 March, 2017; originally announced March 2017.

Comments: 18 pages, 4 figures, accepted for publication in J. Plasma Phys. (Collection: "The Vlasov equation: from space to laboratory plasma physics")

Proton fire hose instabilities in the expanding solar wind

Authors: Petr Hellinger

Abstract: Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For quasi radial ambient magnetic field the expansion drives $T_{\mathrm{p}\|}>T_{\mathrm{p}\perp}$ and the system becomes eventually unstable with respect to the dominant parallel fi… ▽ More Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For quasi radial ambient magnetic field the expansion drives $T_{\mathrm{p}\|}>T_{\mathrm{p}\perp}$ and the system becomes eventually unstable with respect to the dominant parallel fire hose instability. This instability is generally unable to counteract the induced anisotropization and the system typically becomes unstable with respect to the oblique fire hose instability later on. The oblique instability efficiently reduces the anisotropy and the system rapidly stabilizes while a significant part of the generated electromagnetic fluctuations is damped to protons. As long as the magnetic field is in the quasi radial direction, this evolution repeats itself and the electromagnetic fluctuations accumulate. For sufficiently oblique magnetic field the expansion drives $T_{\mathrm{p}\perp}>T_{\mathrm{p}\|}$ and brings the system to the stable region with respect to the fire hose instabilities. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: Journal of Plasma Physics, 14 pages, 9 figures

Journal ref: JPP, 83, 705830105, 2017

Plasma beta dependence of the ion-scale spectral break of solar wind turbulence: high-resolution 2D hybrid simulations

Authors: Luca Franci, Simone Landi, Lorenzo Matteini, Andrea Verdini, Petr Hellinger

Abstract: We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, di… ▽ More We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to -5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with $β$ (from around -3.6 for small values to around -2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, $d_i$, for $β\ll 1$ and to the ion gyroradius, $ρ_i = d_i\sqrtβ$, for $β\gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $β\sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $β$ that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: 8 pages, 4 figures, accepted for publication on The Astrophysical Journal

Proton heating by pick-up ion driven cyclotron waves in the outer heliosphere: Hybrid expanding box simulations

Authors: Petr Hellinger, Pavel M. Travnicek

Abstract: Using one-dimensional hybrid expanding box model we investigate properties of the solar wind in the outer heliosphere. We assume a proton-electron plasma with a strictly transverse ambient magnetic field and, beside the expansion, we take into account influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of inter… ▽ More Using one-dimensional hybrid expanding box model we investigate properties of the solar wind in the outer heliosphere. We assume a proton-electron plasma with a strictly transverse ambient magnetic field and, beside the expansion, we take into account influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of interstellar origin. The injected cold pick-up protons form a ring distribution function that rapidly becomes unstable and generate Alfvén cyclotron waves. The Alfvén cyclotron waves scatter pick-up protons to a spherical shell distribution function that thickens over that time owing to the expansion-driven cooling. The Alfven cyclotron waves heat solar wind protons in the perpendicular direction (with respect to the ambient magnetic field) through the cyclotron resonance. At later times, the Alfvén cyclotron waves become parametrically unstable and the generated ion acoustic waves heat protons in the parallel direction through the Landau resonance. The resulting heating of the solar wind protons is efficient on the expansion time scale. △ Less

Submitted 7 September, 2016; originally announced September 2016.

Comments: ApJ, 7 pages, 7 figures

Journal ref: ApJ, 832, 32, 2016

Ion collisional transport coefficients in the solar wind at 1 AU

Authors: Petr Hellinger

Abstract: Proton and alpha particle collisional transport coefficients (isotropization, relative deceleration frequencies and heating rates) at 1 AU are quantified using the WIND/SWE data. In agreement with previous studies the ion-ion Coulomb collisions are generally important for slow solar wind streams and tend to reduce the temperature anisotropies, the differential streaming and the differences between… ▽ More Proton and alpha particle collisional transport coefficients (isotropization, relative deceleration frequencies and heating rates) at 1 AU are quantified using the WIND/SWE data. In agreement with previous studies the ion-ion Coulomb collisions are generally important for slow solar wind streams and tend to reduce the temperature anisotropies, the differential streaming and the differences between proton and alpha particle temperatures. In slow solar wind streams the Coulomb collisions between protons and alpha particles are important for the overall proton energetics as well as for the relative deceleration between the two species. It is also shown that ion temperature anisotropies and differential streaming need to be generally taken into account for evaluation of the collisional transport coefficients. △ Less

Submitted 3 May, 2016; originally announced May 2016.

Comments: ApJ, 5 pages, 4 figures

Two-dimensional Hybrid Simulations of Kinetic Plasma Turbulence: Current and Vorticity vs Proton Temperature

Authors: Luca Franci, Petr Hellinger, Lorenzo Matteini, Andrea Verdini, Simone Landi

Abstract: Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate… ▽ More Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Comments: 4 pages, 2 figures, Proceedings of the Fourteenth International Solar Wind Conference

Journal ref: AIP Conf. Proc. 1720, 040003 (2016)

Microstructure in two- and three-dimensional hybrid simulations of perpendicular collisionless shocks

Authors: David Burgess, Petr Hellinger, Peter W. Gingell, Pavel M. Trávníček

Abstract: Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been… ▽ More Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfvén speed of the ramp, and motion in the sense of gyration of the reflected ions. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Submitted to Journal of Plasma Physics

Fire Hose instability driven by alpha particle temperature anisotropy

Authors: Lorenzo Matteini, Petr Hellinger, Steven Schwartz, Simone Landi

Abstract: We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated… ▽ More We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma. △ Less

Submitted 23 August, 2015; originally announced August 2015.

Comments: 10 pages, 9 figures, Accepted for publication in The Astrophysical Journal

Plasma turbulence and kinetic instabilities at ion scales in the expanding solar wind

Authors: Petr Hellinger, Lorenzo Matteini, Simone Landi, Andrea Verdini, Luca Franci, Pavel M. Travnicek

Abstract: The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition… ▽ More The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime. △ Less

Submitted 29 September, 2015; v1 submitted 13 August, 2015; originally announced August 2015.

Comments: 6 pages, 5 figures, in press in ApJL

Journal ref: Astrophys. J. Lett. (2015), 811, L32

High-resolution hybrid simulations of kinetic plasma turbulence at proton scales

Authors: Luca Franci, Simone Landi, Lorenzo Matteini, Andrea Verdini, Petr Hellinger

Abstract: We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbu… ▽ More We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters. △ Less

Submitted 19 June, 2015; originally announced June 2015.

Comments: 15 pages, 13 figures, submitted to ApJ

Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations

Authors: Luca Franci, Andrea Verdini, Lorenzo Matteini, Simone Landi, Petr Hellinger

Abstract: We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several propert… ▽ More We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law. △ Less

Submitted 18 April, 2015; v1 submitted 18 March, 2015; originally announced March 2015.

Comments: 5 pages, 5 figures; introduction and conclusions changed, references updated, accepted for publication in ApJL

Anisotropy of third-order structure functions in MHD turbulence

Authors: Andrea Verdini, Roland Grappin, Petr Hellinger, Simone Landi, Wolf-Christian Müller

Abstract: The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (… ▽ More The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfvén mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion. △ Less

Submitted 2 April, 2015; v1 submitted 16 February, 2015; originally announced February 2015.

Comments: Some text editing and typos corrected, 13 pages, 6 figures, to be published in ApJ

Proton temperature-anisotropy-driven instabilities in weakly collisional plasmas: Hybrid simulations

Authors: Petr Hellinger, Pavel M. Travnicek

Abstract: Kinetic instabilities in weakly collisional, high beta plasmas are investigated using two-dimensional hybrid expanding box simulations with Coulomb collisions modeled through the Langevin equation (corresponding to the Fokker-Planck one). The expansion drives a parallel or perpendicular temperature anisotropy (depending on the orientation of the ambient magnetic field). For the chosen parameters t… ▽ More Kinetic instabilities in weakly collisional, high beta plasmas are investigated using two-dimensional hybrid expanding box simulations with Coulomb collisions modeled through the Langevin equation (corresponding to the Fokker-Planck one). The expansion drives a parallel or perpendicular temperature anisotropy (depending on the orientation of the ambient magnetic field). For the chosen parameters the Coulomb collisions are important with respect to the driver but are not strong enough to keep the system stable with respect to instabilities driven by the proton temperature anisotropy. In the case of the parallel temperature anisotropy the dominant oblique fire hose instability efficiently reduces the anisotropy in a quasilinear manner. In the case of the perpendicular temperature anisotropy the dominant mirror instability generates coherent compressive structures which scatter protons and reduce the temperature anisotropy. For both the cases the instabilities generate temporarily enough wave energy so that the corresponding (anomalous) transport coefficients dominate over the collisional ones and their properties are similar to those in collisionless plasmas. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: in press in J. Plasma Phys

Journal ref: J. Plasma Phys. 81 (2015) 305810103

Solar wind protons at 1 AU: trends and bounds, constraints and correlations

Authors: Petr Hellinger, Pavel M. Travnicek

Abstract: The proton temperature anisotropy in the solar wind exhibits apparent bounds which are compatible with the theoretical constraints imposed by temperature-anisotropy driven kinetic instabilities. Recent statistical analyses based on conditional averaging indicate that near these theoretical constraints the solar wind protons have typically enhanced temperatures and a weaker collisionality. Here we… ▽ More The proton temperature anisotropy in the solar wind exhibits apparent bounds which are compatible with the theoretical constraints imposed by temperature-anisotropy driven kinetic instabilities. Recent statistical analyses based on conditional averaging indicate that near these theoretical constraints the solar wind protons have typically enhanced temperatures and a weaker collisionality. Here we carefully analyze the solar wind data and show that these results are a consequence of superposition of multiple correlations in the solar wind, namely, they mostly result from the correlation between the proton temperature and the solar wind velocity and from the superimposed anti-correlation between the proton temperature anisotropy and the proton parallel beta in the fast solar wind. Colder and more collisional data are distributed around temperature isotropy whereas hotter and less collisional data have a wider range of the temperature anisotropy anti-correlated with the proton parallel beta with signatures of constraints owing to the temperature-anisotropy driven instabilities. However, most of the hot and weakly collisional data, including the hottest and least collisional ones, lies far from the marginal stability regions. Consequently, we conclude that there is no clear relation between the enhanced temperatures and instability constraints and that the conditional averaging used for these analyses must be used carefully and need to be well tested. △ Less

Submitted 19 February, 2014; originally announced February 2014.

Comments: Accepted for publication in ApJL

Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: from fluid to kinetic modeling

Authors: P. Henri, S. S. Cerri, F. Califano, F. Pegoraro, C. Rossi, M. Faganello, O. Šebek, P. M. Trávníček, P. Hellinger, J. T. Frederiksen, Å. Nordlund, S. Markidis, R. Keppens, G. Lapenta

Abstract: The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understand… ▽ More The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime. △ Less

Submitted 29 October, 2013; originally announced October 2013.

Journal ref: Phys. Plasmas 20, 102118 (2013)

Cross-Scale: Multi-Scale Coupling in Space Plasma, Assessment Study Report

Authors: Steve Schwartz, Stuart D. Bale, Masaki Fujimoto, Petr Hellinger, Mona Kessel, Guan Le, William Liu, Philippe Louarn, Ian Mann, Rumi Nakamura, Chris Owen, Jean-Louis Pinçon, Luca Sorriso-Valvo, Andris Vaivads, Robert F. Wimmer-Schweingruber

Abstract: Driven by the support and interest of the international space plasma community to examine simultaneous physical plasma scales and their interactions, the Cross-Scale Mission concept was submitted and accepted as an ESA Cosmic Vision M-class candidate mission. This report presents an overview of the assessment study phase of the 7 ESA spacecraft Cross-Scale mission. Where appropriate, discussion… ▽ More Driven by the support and interest of the international space plasma community to examine simultaneous physical plasma scales and their interactions, the Cross-Scale Mission concept was submitted and accepted as an ESA Cosmic Vision M-class candidate mission. This report presents an overview of the assessment study phase of the 7 ESA spacecraft Cross-Scale mission. Where appropriate, discussion of the benefit of international collaboration with the SCOPE mission, as well as other interested parties, is included. △ Less

Submitted 4 December, 2009; originally announced December 2009.

Comments: Assessment Study Report as part of the ESA Cosmic Vision M-Class mission selection process

Report number: ESA/SRE[2009]1

Quantum-double for a nonstandard deformation of a Borel subalgebra $sl(2,C)$

Authors: C. Burdik, P. Hellinger

Abstract: We give a construction of Drienfeld's quantum double for a nonstandard deformation of Borel subalgebra of $sl(2)$. We construct explicitly some simple representations of this quantum algebra and from the universal R-matrix we obtain the explicit solutions of the Yang-Baxter equation in those cases. We give a construction of Drienfeld's quantum double for a nonstandard deformation of Borel subalgebra of $sl(2)$. We construct explicitly some simple representations of this quantum algebra and from the universal R-matrix we obtain the explicit solutions of the Yang-Baxter equation in those cases. △ Less

Submitted 5 March, 1993; originally announced March 1993.

Comments: 9 pages

Report number: PRA-HEP-93/2