Showing 1–50 of 73 results for author: Loureiro, N F

Electron-only reconnection and inverse magnetic-energy transfer at sub-ion scales

Authors: Zhuo Liu, Caio Silva, Lucio M. Milanese, Muni Zhou, Noah R. Mandell, Nuno F. Loureiro

Abstract: We derive, and validate numerically, an analytical model for electron-only magnetic reconnection applicable to strongly magnetized (low-beta) plasmas. Our model predicts sub-ion-scale reconnection rates significantly higher than those pertaining to large-scale reconnection, aligning with recent observations and simulations. We apply this reconnection model to the problem of inverse magnetic-energy… ▽ More We derive, and validate numerically, an analytical model for electron-only magnetic reconnection applicable to strongly magnetized (low-beta) plasmas. Our model predicts sub-ion-scale reconnection rates significantly higher than those pertaining to large-scale reconnection, aligning with recent observations and simulations. We apply this reconnection model to the problem of inverse magnetic-energy transfer at sub-ion scales. We derive time-dependent scaling laws for the magnetic energy decay and the typical magnetic structure dimensions that differ from those previously found in the MHD regime. These scaling laws are validated via two- and three-dimensional simulations, demonstrating that sub-ion scale magnetic fields can reach large, system-size scales via successive coalescence. △ Less

Submitted 8 July, 2024; originally announced July 2024.

A transport-like approach to turbulence-mediated magnetic reconnection

Authors: Nuno F. Loureiro

Abstract: A transport-like framework for the study of magnetic reconnection mediated by self-driven turbulence is proposed, based on timescale separation between the reconnection time and the characteristic timescale of the turbulent fluctuations which arise in the reconnection layer. We argue that the mean fields remain on MHD scales even in collisionless cases. These observations provide theoretical justi… ▽ More A transport-like framework for the study of magnetic reconnection mediated by self-driven turbulence is proposed, based on timescale separation between the reconnection time and the characteristic timescale of the turbulent fluctuations which arise in the reconnection layer. We argue that the mean fields remain on MHD scales even in collisionless cases. These observations provide theoretical justification for an efficient computational approach to this problem, which we discuss. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 7 pages, 1 figure

A Two-dimensional Numerical Study of Ion-Acoustic Turbulence

Authors: Zhuo Liu, Ryan White, Lucio M. Milanese, Nuno F. Loureiro

Abstract: We investigate the linear and nonlinear evolution of the ion-acoustic instability in a collisionless plasma via two-dimensional (2D2V) Vlasov-Poisson numerical simulations. We initialize the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding super-critical initial conditions that are physically unrealizable. The nonlinea… ▽ More We investigate the linear and nonlinear evolution of the ion-acoustic instability in a collisionless plasma via two-dimensional (2D2V) Vlasov-Poisson numerical simulations. We initialize the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding super-critical initial conditions that are physically unrealizable. The nonlinear evolution of ion-acoustic turbulence (IAT) is characterized in detail, including the particles' distribution functions, particle heating, (two-dimensional) wave spectrum, and the resulting anomalous resistivity. An important result is that no steady saturated nonlinear state is ever reached in our simulations: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived. Electron-acoustic waves (EAWs) are triggered during the late nonlinear evolution of the system, caused by strong modifications to the particle distribution induced by IAT. △ Less

Submitted 2 January, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Magnetogenesis in a collisionless plasma: from Weibel instability to turbulent dynamo

Authors: Muni Zhou, Vladimir Zhdankin, Matthew W. Kunz, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We report on a first-principles numerical and theoretical study of plasma dynamo in a fully kinetic framework. By applying an external mechanical force to an initially unmagnetized plasma, we develop a self-consistent treatment of the generation of ``seed'' magnetic fields, the formation of turbulence, and the inductive amplification of fields by the fluctuation dynamo. Driven large-scale motions… ▽ More We report on a first-principles numerical and theoretical study of plasma dynamo in a fully kinetic framework. By applying an external mechanical force to an initially unmagnetized plasma, we develop a self-consistent treatment of the generation of ``seed'' magnetic fields, the formation of turbulence, and the inductive amplification of fields by the fluctuation dynamo. Driven large-scale motions in an unmagnetized, weakly collisional plasma are subject to strong phase mixing, which leads to the development of thermal pressure anisotropy. This anisotropy triggers the Weibel instability, which produces filamentary ``seed'' magnetic fields on plasma-kinetic scales. The plasma is thereby magnetized, enabling efficient stretching and folding of the fields by the plasma motions and the development of Larmor-scale kinetic instabilities such as the firehose and mirror. The scattering of particles off the associated microscale magnetic fluctuations provides an effective viscosity, regulating the field morphology and turbulence. During this process, the seed field is further amplified by the fluctuation dynamo until they reach energy equipartition with the turbulent flow. By demonstrating that equipartition magnetic fields can be generated from an initially unmagnetized plasma through large-scale turbulent flows, this work has important implications for the origin and amplification of magnetic fields in the intracluster and intergalactic mediums. △ Less

Submitted 28 July, 2023; originally announced August 2023.

Comments: 16 pages, 10 figures

A quantum algorithm for the linear Vlasov equation with collisions

Authors: Abtin Ameri, Paola Cappellaro, Hari Krovi, Nuno F. Loureiro, Erika Ye

Abstract: The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability. In this work, we present a quantum algorithm that simulates the linearized Vlasov equation with and without collisions, in the one-dimensional, electrostatic l… ▽ More The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability. In this work, we present a quantum algorithm that simulates the linearized Vlasov equation with and without collisions, in the one-dimensional, electrostatic limit. Rather than solving this equation in its native spatial and velocity phase-space, we adopt an efficient representation in the dual space yielded by a Fourier-Hermite expansion. The Fourier-Hermite representation is exponentially more compact, thus yielding a classical algorithm that can match the performance of a previously proposed quantum algorithm for this problem. This representation results in a system of linear ordinary differential equations which can be solved with well-developed quantum algorithms: Hamiltonian simulation in the collisionless case, and quantum ODE solvers in the collisional case. In particular, we demonstrate that a quadratic speedup in system size is attainable. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Observation of electromagnetic filamentary structures produced by the Weibel instability in laser-driven plasmas

Authors: G. D. Sutcliffe, P. J. Adrian, J. A. Pearcy, T. M. Johnson, J. Kunimune, B. Pollock, J. D. Moody, N. F. Loureiro, C. K. Li

Abstract: We present experimental observations of electron-scale structures in an expanding high-energy-density (HED) plasma generated with a modest intensity $\sim 2 \times 10^{14}$ W/cm$^2$, $\sim1$ ns laser. The observed structures have wavelengths ($\sim 150-220 μみゅーm$) and growth rates ($\sim 0.4-1.0$ ns$^{-1}$) consistent with an electron-driven Weibel instability where the anisotropy in the electron dis… ▽ More We present experimental observations of electron-scale structures in an expanding high-energy-density (HED) plasma generated with a modest intensity $\sim 2 \times 10^{14}$ W/cm$^2$, $\sim1$ ns laser. The observed structures have wavelengths ($\sim 150-220 μみゅーm$) and growth rates ($\sim 0.4-1.0$ ns$^{-1}$) consistent with an electron-driven Weibel instability where the anisotropy in the electron distribution is small, $A\sim 0.002$. This instability is found to be a better match to the observed phenomena than other typical field-generation mechanisms found in HED plasmas, including counter-streaming ion Weibel and magnetothermal instabilities. These observations experimentally demonstrate for the first time that the electron Weibel instability must be considered alongside other magnetic field generation and amplification mechanisms in expanding ablation plasmas, which are ubiquitous in HED research. They also provide physics insight into the generation of magnetic fields in large-scale astrophysical plasmas. Additionally, inspection of the magnetic power spectrum shows a possible scaling match to analytic gyrokinetic predictions, $|B_k|^2\propto k^{-16/3}$, at scales below the electron Larmor radius. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 6 pages, 5 figures

Intermittency and electron heating in kinetic-Alfvén-wave turbulence

Authors: Muni Zhou, Zhuo Liu, Nuno F. Loureiro

Abstract: We report analytical and numerical investigations of sub-ion-scale turbulence in low-beta plasmas, focusing on the spectral properties of the fluctuations and electron heating. In the isothermal limit, the numerical results strongly support a description of the turbulence as a critically-balanced Kolmogorov-like cascade of kinetic Alfvén wave fluctuations, as amended by Boldyrev & Perez (Astrophys… ▽ More We report analytical and numerical investigations of sub-ion-scale turbulence in low-beta plasmas, focusing on the spectral properties of the fluctuations and electron heating. In the isothermal limit, the numerical results strongly support a description of the turbulence as a critically-balanced Kolmogorov-like cascade of kinetic Alfvén wave fluctuations, as amended by Boldyrev & Perez (Astrophys. J. Lett. 758, L44 (2012)) to include intermittent effects. When the constraint of isothermality is removed (i.e., with the inclusion of electron kinetic physics), the energy spectrum is found to steepen due to electron Landau damping, which is enabled by the local weakening of advective nonlinearities around current sheets, and yields significant energy dissipation via a velocity-space cascade. The use of a Hermite-polynomial representation to express the velocity-space dependence of the electron distribution function allows us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 19 pages, 10 figures

Experimental Evidence of Plasmoids in High-$βべーた$ Magnetic Reconnection

Authors: J. A. Pearcy, M. J. Rosenberg, T. M. Johnson, G. D. Sutcliffe, B. L. Reichelt, J. D. Hare, N. F. Loureiro, R. D. Petrasso, C. K. Li

Abstract: Magnetic reconnection is a ubiquitous and fundamental process in plasmas by which magnetic fields change their topology and release magnetic energy. Despite decades of research, the physics governing the reconnection process in many parameter regimes remains controversial. Contemporary reconnection theories predict that long, narrow current sheets are susceptible to the tearing instability and spl… ▽ More Magnetic reconnection is a ubiquitous and fundamental process in plasmas by which magnetic fields change their topology and release magnetic energy. Despite decades of research, the physics governing the reconnection process in many parameter regimes remains controversial. Contemporary reconnection theories predict that long, narrow current sheets are susceptible to the tearing instability and split into isolated magnetic islands (or plasmoids), resulting in an enhanced reconnection rate. While several experimental observations of plasmoids in the regime of low- to intermediate-$βべーた$ (where $βべーた$ is the ratio of plasma thermal pressure to magnetic pressure) have been made, there is a relative lack of experimental evidence for plasmoids in the high-$βべーた$ reconnection environments which are typical in many space and astrophysical contexts. Here, we report the observation of strong experimental evidence for plasmoid formation and dynamics in laser-driven high-$βべーた$ reconnection experiments. △ Less

Submitted 8 July, 2022; originally announced July 2022.

A quantum-inspired method for solving the Vlasov-Poisson equations

Authors: Erika Ye, Nuno F. G. Loureiro

Abstract: Kinetic simulations of collisionless (or weakly collisional) plasmas using the Vlasov equation are often infeasible due to high resolution requirements and the exponential scaling of computational cost with respect to dimension. Recently, it has been proposed that matrix product state (MPS) methods, a quantum-inspired but classical algorithm, can be used to solve partial differential equations wit… ▽ More Kinetic simulations of collisionless (or weakly collisional) plasmas using the Vlasov equation are often infeasible due to high resolution requirements and the exponential scaling of computational cost with respect to dimension. Recently, it has been proposed that matrix product state (MPS) methods, a quantum-inspired but classical algorithm, can be used to solve partial differential equations with exponential speed-up, provided that the solution can be compressed and efficiently represented as an MPS within some tolerable error threshold. In this work, we explore the practicality of MPS methods for solving the Vlasov-Poisson equations in 1D1V, and find that important features of linear and nonlinear dynamics, such as damping or growth rates and saturation amplitudes, can be captured while compressing the solution significantly. Furthermore, by comparing the performance of different mappings of the distribution functions onto the MPS, we develop an intuition of the MPS representation and its behavior in the context of solving the Vlasov-Poisson equations, which will be useful for extending these methods to higher dimensional problems. △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: main text 14 pages, 7 figures; SI appended

Hybrid Quantum-Classical Algorithm for Hydrodynamics

Authors: Julien Zylberman, Giuseppe Di Molfetta, Marc Brachet, Nuno F. Loureiro, Fabrice Debbasch

Abstract: A new model of nonlinear charged quantum relativistic fluids is presented. This model can be discretized into Discrete Time Quantum Walks (DTQWs), and a new hybrid (quantum-classical) algorithm for implementing these walks on NISQ devices is proposed. High resolution (up to $N=2^{17}$ grid points) hybrid numerical simulations of relativistic and non-relativistic hydrodynamical shocks on current IB… ▽ More A new model of nonlinear charged quantum relativistic fluids is presented. This model can be discretized into Discrete Time Quantum Walks (DTQWs), and a new hybrid (quantum-classical) algorithm for implementing these walks on NISQ devices is proposed. High resolution (up to $N=2^{17}$ grid points) hybrid numerical simulations of relativistic and non-relativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm and shown to reproduce equivalent simulations on classical computers. This work demonstrates that nonlinear fluid dynamics can be simulated on NISQs, and opens the door to simulating other, quantum and non-quantum fluids, including plasmas, with more general quantum walks and quantum automata. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 11 pages, 5 figures

Spontaneous magnetization of collisionless plasma through the action of a shear flow

Authors: Muni Zhou, Vladimir Zhdankin, Matthew W. Kunz, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the se… ▽ More We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated seed fields on key parameters (e.g., ratio of system scale to electron skin depth, the forcing amplitude) are confirmed by 3D and 2D particle-in-cell simulations using an electron-positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow, and therefore has important implications for magnetogenesis in dilute astrophysical systems. △ Less

Submitted 3 October, 2021; originally announced October 2021.

Comments: 20 pages, 16 figures

Statistical description of coalescing magnetic islands via magnetic reconnection

Authors: Muni Zhou, David H. Wu, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath, and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral… ▽ More The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath, and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in 2D. We solve our island kinetic equation numerically for three different types of initial distribution: delta-distribution, Gaussian and power-law distribution. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde t^{-1}$, the number of islands decreases as $\tilde t^{-1}$, and the averaged area of islands grows as $\tilde t$, where $\tilde t$ is the time normalized to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. △ Less

Submitted 30 October, 2021; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 38 pages; 10 figures; Accepted by Journal of Plasma Physics

Dynamic Phase Alignment in Navier-Stokes Turbulence

Authors: Lucio M. Milanese, Nuno F. Loureiro, Stanislav Boldyrev

Abstract: In Navier-Stokes turbulence, energy and helicity injected at large scales are subject to a joint direct cascade, with both quantities exhibiting a spectral scaling $\propto k^{-5/3}$. We demonstrate via direct numerical simulations that the two cascades are compatible due to the existence of a strong scale-dependent phase alignment between velocity and vorticity fluctuations, with the phase alignm… ▽ More In Navier-Stokes turbulence, energy and helicity injected at large scales are subject to a joint direct cascade, with both quantities exhibiting a spectral scaling $\propto k^{-5/3}$. We demonstrate via direct numerical simulations that the two cascades are compatible due to the existence of a strong scale-dependent phase alignment between velocity and vorticity fluctuations, with the phase alignment angle scaling as $\cosαあるふぁ_k\propto k^{-1}$. △ Less

Submitted 5 November, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Laminar and Turbulent Plasmoid Ejection in a Laboratory Parker Spiral Current Sheet

Authors: Ethan E. Peterson, Douglass A. Endrizzi, Michael Clark, Jan Egedal, Kenneth Flanagan, Nuno F. Loureiro, Jason Milhone, Joseph Olson, Carl R. Sovinec, John Wallace, Cary B. Forest

Abstract: Quasi-periodic plasmoid formation at the tip of magnetic streamer structures is observed to occur in experiments on the Big Red Ball as well as in simulations of these experiments performed with the extended-MHD code, NIMROD. This plasmoid formation is found to occur on a characteristic timescale dependent on pressure gradients and magnetic curvature in both experiment and simulation. Single mode,… ▽ More Quasi-periodic plasmoid formation at the tip of magnetic streamer structures is observed to occur in experiments on the Big Red Ball as well as in simulations of these experiments performed with the extended-MHD code, NIMROD. This plasmoid formation is found to occur on a characteristic timescale dependent on pressure gradients and magnetic curvature in both experiment and simulation. Single mode, or laminar, plasmoids exist when the pressure gradient is modest, but give way to turbulent plasmoid ejection when the system drive is higher, producing plasmoids of many sizes. However, a critical pressure gradient is also observed, below which plasmoids are never formed. A simple heuristic model of this plasmoid formation process is presented and suggested to be a consequence of a dynamic loss of equilibrium in the high-$βべーた$ region of the helmet streamer. This model is capable of explaining the periodicity of plasmoids observed in the experiment and simulations and produces plasmoid periods of 90 minutes when applied to 2D models of solar streamers with a height of $3R_\odot$. This is consistent with the location and frequency at which periodic plasma blobs have been observed to form by LASCO and SECCHI instruments. △ Less

Submitted 22 June, 2021; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: 17 pages, 10 figures, 1 movie

Efficient quantum algorithm for dissipative nonlinear differential equations

Authors: Jin-Peng Liu, Herman Øie Kolden, Hari K. Krovi, Nuno F. Loureiro, Konstantina Trivisa, Andrew M. Childs

Abstract: Nonlinear differential equations model diverse phenomena but are notoriously difficult to solve. While there has been extensive previous work on efficient quantum algorithms for linear differential equations, the linearity of quantum mechanics has limited analogous progress for the nonlinear case. Despite this obstacle, we develop a quantum algorithm for dissipative quadratic $n$-dimensional ordin… ▽ More Nonlinear differential equations model diverse phenomena but are notoriously difficult to solve. While there has been extensive previous work on efficient quantum algorithms for linear differential equations, the linearity of quantum mechanics has limited analogous progress for the nonlinear case. Despite this obstacle, we develop a quantum algorithm for dissipative quadratic $n$-dimensional ordinary differential equations. Assuming $R < 1$, where $R$ is a parameter characterizing the ratio of the nonlinearity and forcing to the linear dissipation, this algorithm has complexity $T^2 q~\mathrm{poly}(\log T, \log n, \log 1/εいぷしろん)/εいぷしろん$, where $T$ is the evolution time, $εいぷしろん$ is the allowed error, and $q$ measures decay of the solution. This is an exponential improvement over the best previous quantum algorithms, whose complexity is exponential in $T$. While exponential decay precludes efficiency, driven equations can avoid this issue despite the presence of dissipation. Our algorithm uses the method of Carleman linearization, for which we give a novel convergence theorem. This method maps a system of nonlinear differential equations to an infinite-dimensional system of linear differential equations, which we discretize, truncate, and solve using the forward Euler method and the quantum linear system algorithm. We also provide a lower bound on the worst-case complexity of quantum algorithms for general quadratic differential equations, showing that the problem is intractable for $R \ge \sqrt{2}$. Finally, we discuss potential applications, showing that the $R < 1$ condition can be satisfied in realistic epidemiological models and giving numerical evidence that the method may describe a model of fluid dynamics even for larger values of $R$. △ Less

Submitted 15 October, 2021; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 36 pages, 1 figure. Published in PNAS

Journal ref: Proceedings of the National Academy of Sciences 118, e2026805118 (2021)

Dynamic Phase Alignment in Inertial Alfven Turbulence

Authors: Lucio M. Milanese, Nuno F. Loureiro, Maximilian Daschner, Stanislav Boldyrev

Abstract: In weakly-collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin-depth, $k_\perp d_e\gtrsim 1$. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence… ▽ More In weakly-collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin-depth, $k_\perp d_e\gtrsim 1$. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as $\cosαあるふぁ_k\propto k_{\perp}^{-1}$. The kinetic and magnetic energy spectra scale as $\propto k_{\perp}^{-5/3}$ and $\propto k_{\perp}^{-11/3}$, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as $\propto k^{-5/3}_{\perp}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the $k_{\perp} d_e \gg 1$ range. Turbulent eddies exhibit a phase-space anisotropy $k_{\parallel} \propto k_{\perp}^{5/3}$, consistent with critically-balanced inertial Alfvén fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Journal ref: Phys. Rev. Lett. 125, 265101 (2020)

Inverse energy transfer in decaying, three dimensional, nonhelical magnetic turbulence due to magnetic reconnection

Authors: Pallavi Bhat, Muni Zhou, Nuno F. Loureiro

Abstract: It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, nonhelical, magnetically dominated, magnetohydrodynamic turbulence in 3-dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to s… ▽ More It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, nonhelical, magnetically dominated, magnetohydrodynamic turbulence in 3-dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and the 3D cases, i.e., the magnetic energy evolves as $t^{-1}$, and the magnetic power spectrum follows a slope of $k^{-2}$. We show that on normalizing time by the magnetic reconnection timescale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse onto one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as $ t^{-1.4}$ and, interestingly, a dynamo effect is observed. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: 13 pages, 18 figures

Multi-scale dynamics of magnetic flux tubes and inverse magnetic energy transfer

Authors: Muni Zhou, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the mag… ▽ More We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in critically-balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field. △ Less

Submitted 13 June, 2020; v1 submitted 20 January, 2020; originally announced January 2020.

Comments: 33 pages, 19 figures, accepted for publication in Journal of Plasma Physics

Tearing instability in Alfvén and kinetic-Alfvén turbulence

Authors: Stanislav Boldyrev, Nuno F. Loureiro

Abstract: Recently, it has been realized that magnetic plasma turbulence and magnetic field reconnection are inherently related phenomena. Turbulent fluctuations generate regions of a sheared magnetic field that become unstable to the tearing instability and reconnection, thus modifying turbulence at the corresponding scales. In this contribution, we give a brief discussion of some recent results on tearing… ▽ More Recently, it has been realized that magnetic plasma turbulence and magnetic field reconnection are inherently related phenomena. Turbulent fluctuations generate regions of a sheared magnetic field that become unstable to the tearing instability and reconnection, thus modifying turbulence at the corresponding scales. In this contribution, we give a brief discussion of some recent results on tearing-mediated magnetic turbulence. We illustrate the main ideas of this rapidly developing field of study by concentrating on two important examples -- magnetohydrodynamic Alfvén turbulence and small-scale kinetic-Alfvén turbulence. Due to various potential applications of these phenomena in space physics and astrophysics, we specifically try not to overload the text by heavy analytical derivations, but rather present a qualitative discussion accessible to a non-expert in the theories of turbulence and reconnection. △ Less

Submitted 27 October, 2019; originally announced October 2019.

Comments: A book chapter in AGU Book "Solar and Heliospheric Plasma Structures" (to appear)

Nonlinear Reconnection in Magnetized Turbulence

Authors: Nuno F. Loureiro, Stanislav Boldyrev

Abstract: Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those co… ▽ More Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those conditions are met, a new turbulence regime is accessed where reconnection-driven energy dissipation becomes common, rather than the rare feature that it must be when they are not. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: Submitted to Astrophysical Journal

Interactions of magnetized plasma flows in pulsed-power driven experiments

Authors: L G Suttle, G C Burdiak, C L Cheung, T Clayson, J W D Halliday, J D Hare, S Rusli, D Russell, E Tubman, A Ciardi, N F Loureiro, J Li, A Frank, S V Lebedev

Abstract: A supersonic flow of magnetized plasma is produced by the application of a 1 MA-peak, 500 ns current pulse to a cylindrical arrangement of parallel wires, known as an inverse wire array. The plasma flow is produced by the JxB acceleration of the ablated wire material, and a magnetic field of several Tesla is embedded at source by the driving current. This setup has been used for a variety of exper… ▽ More A supersonic flow of magnetized plasma is produced by the application of a 1 MA-peak, 500 ns current pulse to a cylindrical arrangement of parallel wires, known as an inverse wire array. The plasma flow is produced by the JxB acceleration of the ablated wire material, and a magnetic field of several Tesla is embedded at source by the driving current. This setup has been used for a variety of experiments investigating the interactions of magnetized plasma flows. In experiments designed to investigate magnetic reconnection, the collision of counter-streaming flows, carrying oppositely directed magnetic fields, leads to the formation of a reconnection layer in which we observe ions reaching temperatures much greater than predicted by classical heating mechanisms. The breakup of this layer under the plasmoid instability is dependent on the properties of the inflowing plasma, which can be controlled by the choice of the wire array material. In other experiments, magnetized shocks were formed by placing obstacles in the path of the magnetized plasma flow. The pile-up of magnetic flux in front of a conducting obstacle produces a magnetic precursor acting on upstream electrons at the distance of the ion inertial length. This precursor subsequently develops into a steep density transition via ion-electron fluid decoupling. Obstacles which possess a strong private magnetic field affect the upstream flow over a much greater distance, providing an extended bow shock structure. In the region surrounding the obstacle the magnetic pressure holds off the flow, forming a void of plasma material, analogous to the magnetopause around planetary bodies with self-generated magnetic fields. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 46th EPS conference on Plasma Physics, Milan, Italy 2019

[Plasma 2020 Decadal] Disentangling the Spatiotemporal Structure of Turbulence Using Multi-Spacecraft Data

Authors: J. M. TenBarge, O. Alexandrova, S. Boldyrev, F. Califano, S. S. Cerri, C. H. K. Chen, G. G. Howes, T. Horbury, P. A. Isenberg, H. Ji, K. G. Klein, C. Krafft, M. Kunz, N. F. Loureiro, A. Mallet, B. A. Maruca, W. H. Matthaeus, R. Meyrand, E. Quataert, J. C. Perez, O. W. Roberts, F. Sahraoui, C. S. Salem, A. A. Schekochihin, H. Spence , et al. (4 additional authors not shown)

Abstract: This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and… ▽ More This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

Role of reconnection in inertial kinetic-Alfven turbulence

Authors: Stanislav Boldyrev, Nuno F. Loureiro

Abstract: In a weakly collisional, low-electron-beta plasma, large-scale Alfvén turbulence transforms into inertial kinetic-Alfvén turbulence at scales smaller than the ion microscale (gyroscale or inertial scale). We propose that at such kinetic scales, the nonlinear dynamics tend to organize turbulent eddies into thin current sheets, consistent with the existence of two conserved integrals of the ideal eq… ▽ More In a weakly collisional, low-electron-beta plasma, large-scale Alfvén turbulence transforms into inertial kinetic-Alfvén turbulence at scales smaller than the ion microscale (gyroscale or inertial scale). We propose that at such kinetic scales, the nonlinear dynamics tend to organize turbulent eddies into thin current sheets, consistent with the existence of two conserved integrals of the ideal equations, energy and helicity. The formation of strongly anisotropic structures is arrested by the tearing instability that sets a critical aspect ratio of the eddies at each scale $a$ in the plane perpendicular to the guide field. This aspect ratio is defined by the balance of the eddy turnover rate and the tearing rate, and varies from $(d_e/a)^{1/2}$ to $d_e/a$ depending on the assumed profile of the current sheets. The energy spectrum of the resulting turbulence varies from $k^{-8/3}$ to $k^{-3}$, and the corresponding spectral anisotropy with respect to the strong background magnetic field from $k_z\lesssim k_\perp^{2/3}$ to $k_z\lesssim k_\perp$. △ Less

Submitted 27 October, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: published version

Journal ref: Phys. Rev. Research 1, 012006 (2019)

Magnetic island merger as a mechanism for inverse magnetic energy transfer

Authors: Muni Zhou, Pallavi Bhat, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magneti… ▽ More Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magnetic energy decays as $\tilde t^{-1}$, where $\tilde t$ is time normalized to the (appropriately defined) reconnection timescale; and the correlation length of the field grows as $\tilde t^{1/2}$. The magnetic energy spectrum is self-similar, and evolves as $\propto \tilde t ^{-3/2}k^{-2}$, where the $k$-dependence is imparted by the formation of thin current sheets. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: 6 pages, 5 figures, submitted for publication

Journal ref: Phys. Rev. Research 1, 012004 (2019)

Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality

Authors: R. Jorge, P. Ricci, S. Brunner, S. Gamba, V. Konovets, N. F. Loureiro, L. M. Perrone, N. Teixeira

Abstract: The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength… ▽ More The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard-Bernstein and Dougherty collision operators, and the electron-ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode. △ Less

Submitted 5 March, 2019; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 23 pages, 10 figures, accepted for publication on Journal of Plasma Physics

Journal ref: J. Plasma Phys. 85 (2019) 905850211

Inertial Kinetic-Alfvén Turbulence

Authors: Vadim Roytershteyn, Stanislav Boldyrev, Gian Luca Delzanno, Christopher H. K. Chen, Daniel Grošelj, Nuno F. Loureiro

Abstract: Recent observational and analytical studies suggested that a new regime of kinetic turbulence may exist in plasma environments with low electron beta (Chen and Boldyrev, 2017). Such a regime, termed inertial kinetic-Alfvén turbulence, is relevant for the solar corona, Earth's magnetosheath, and other astrophysical systems where the electron and ion plasma beta parameters satisfy the condition… ▽ More Recent observational and analytical studies suggested that a new regime of kinetic turbulence may exist in plasma environments with low electron beta (Chen and Boldyrev, 2017). Such a regime, termed inertial kinetic-Alfvén turbulence, is relevant for the solar corona, Earth's magnetosheath, and other astrophysical systems where the electron and ion plasma beta parameters satisfy the condition $βべーた_e\ll βべーた_i\lesssim 1$. In this paper we present kinetic numerical simulations that confirm existence of the iKAW regime. Specifically, the simulations demonstrate a transition at scales below electron inertial length $d_e$ when $βべーた_e\ll βべーた_i\lesssim 1$. Spectral slopes and other statistical properties of turbulence at sub-$d_e$ scales are consistent with the phenomenological theory of inertial kinetic-Alfvén turbulence proposed by Chen and Boldyrev (2017) and with the recent observations in the Earth's magnetosheath. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Dependence of alpha-particle-driven Alfvén eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Authors: E. A. Tolman, N. F. Loureiro, P. Rodrigues, J. W. Hughes, E. S. Marmar

Abstract: Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conventional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfvén eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabiliz… ▽ More Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conventional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfvén eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently-developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvénic or nearly sub-Alfvénic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities. △ Less

Submitted 16 February, 2019; v1 submitted 19 September, 2018; originally announced September 2018.

Comments: Accepted to Nuclear Fusion

Calculations in the theory of tearing instability

Authors: Stanislav Boldyrev, Nuno F. Loureiro

Abstract: Recent studies have suggested that the tearing instability may play a significant role in magnetic turbulence. In this work, we review the theory of the magnetohydrodynamic tearing instability in the general case of an arbitrary tearing parameter, which is relevant for applications in turbulence. We discuss a detailed derivation of the results for the standard Harris profile and accompany it by th… ▽ More Recent studies have suggested that the tearing instability may play a significant role in magnetic turbulence. In this work, we review the theory of the magnetohydrodynamic tearing instability in the general case of an arbitrary tearing parameter, which is relevant for applications in turbulence. We discuss a detailed derivation of the results for the standard Harris profile and accompany it by the derivation of the results for a lesser known sine-shaped profile. We devote special attention to the exact solution of the inner equation, which is the central result in the theory of tearing instability. We also briefly discuss the influence of shear flows on tearing instability in magnetic structures. Our presentation is self-contained; we expect it to be accessible to researchers in plasma turbulence who are not experts in magnetic reconnection. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: Proceedings of the 17th Annual International Astrophysics Conference (2018)

Theory of the Drift-Wave Instability at Arbitrary Collisionality

Authors: R. Jorge, P. Ricci, N. F. Loureiro

Abstract: A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis, to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show… ▽ More A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis, to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show that our framework allows retrieving established collisional and collisionless limits. At the intermediate collisionalities relevant for present and future magnetic nuclear fusion devices, deviations with respect to collision operators used in state-of-the-art turbulence simulation codes show the need for retaining the full Coulomb operator in order to obtain both the correct instability growth rate and eigenmode spectrum, which, for example, may significantly impact quantitative predictions of transport. The exponential convergence of the spectral representation that we propose makes the representation of the velocity space dependence, including the full collision operator, more efficient than standard finite difference methods. △ Less

Submitted 18 September, 2018; v1 submitted 18 June, 2018; originally announced June 2018.

Comments: 7 pages, 3 figures, accepted for publication on Physical Review Letters

Journal ref: Phys. Rev. Lett. 121, 165001 (2018)

Turbulence in Magnetized Pair Plasmas

Authors: N. F. Loureiro, S. Boldyrev

Abstract: Alfvénic-type turbulence in strongly magnetized, low-beta pair plasmas is investigated. A coupled set of equations for the evolution of the magnetic and flow potentials are derived, covering both fluid and kinetic scales. In the fluid (MHD) range those equations are the same as for electron-ion plasmas, so turbulence at those scales is expected to be of the Alfvénic nature, exhibiting critical bal… ▽ More Alfvénic-type turbulence in strongly magnetized, low-beta pair plasmas is investigated. A coupled set of equations for the evolution of the magnetic and flow potentials are derived, covering both fluid and kinetic scales. In the fluid (MHD) range those equations are the same as for electron-ion plasmas, so turbulence at those scales is expected to be of the Alfvénic nature, exhibiting critical balance, dynamic alignment, and transition to a tearing mediated regime at small scales. The critical scale at which a transition to a tearing-mediated range occurs is derived, and the spectral slope in that range is predicted to be $k_\perp^{-8/3}$ (or $k_\perp^{-3}$, depending on details of the reconnecting configuration assumed). At scales below the electron (and positron) skin depth, it is argued that turbulence is dictated by a cascade of the inertial Alfvén wave, which we show to result in the magnetic energy spectrum $\propto k_\perp^{-11/3}$. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: submitted for publication

Ion heating and magnetic flux pile-up in a magnetic reconnection experiment with super-Alfvenic plasma inflows

Authors: L. G. Suttle, J. D. Hare, S. V. Lebedev, A. Ciardi, N. F. Loureiro, G. C. Burdiak, J. P. Chittenden, T. Clayson, J. W. D. Halliday, N. Niasse, D. Russell, F. Suzuki Vidal, E. Tubman, T. Lane, J. Ma, T. Robinson, R. A. Smith, N. Stuart

Abstract: This work presents a magnetic reconnection experiment in which the kinetic, magnetic and thermal properties of the plasma each play an important role in the overall energy balance and structure of the generated reconnection layer. Magnetic reconnection occurs during the interaction of continuous and steady flows of super-Alfvenic, magnetized, aluminum plasma, which collide in a geometry with two-d… ▽ More This work presents a magnetic reconnection experiment in which the kinetic, magnetic and thermal properties of the plasma each play an important role in the overall energy balance and structure of the generated reconnection layer. Magnetic reconnection occurs during the interaction of continuous and steady flows of super-Alfvenic, magnetized, aluminum plasma, which collide in a geometry with two-dimensional symmetry, producing a stable and long-lasting reconnection layer. Optical Thomson scattering measurements show that when the layer forms, ions inside the layer are more strongly heated than electrons, reaching temperatures of Ti~ZTe>300 eV - much greater than can be expected from strong shock and viscous heating alone. Later in time, as the plasma density in the layer increases, the electron and ion temperatures are found to equilibrate, and a constant plasma temperature is achieved through a balance of the heating mechanisms and radiative losses of the plasma. Measurements from Faraday rotation polarimetry also indicate the presence of significant magnetic field pile-up occurring at the boundary of the reconnection region, which is consistent with the super-Alfvenic velocity of the inflows. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: 15 pages, 8 figures, 1 table

Fully kinetic large scale simulations of the collisionless Magnetorotational instability

Authors: Giannandrea Inchingolo, Thomas Grismayer, Nuno F. Loureiro, Ricardo A. Fonseca, Luis O. Silva

Abstract: We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $βべーた$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the… ▽ More We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $βべーた$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI, together with ensuing magnetic reconnection, activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady state, the magnetic energy has an MHD-like spectrum with a slope of $k^{-5/3}$ for $kρろー<1$ and $k^{-3}$ for sub-Larmor scale ($kρろー>1$). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. In particular, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of $εいぷしろん^{-2}$ in the high-energy region. △ Less

Submitted 29 April, 2018; v1 submitted 25 January, 2018; originally announced January 2018.

Comments: 13 pages, 9 figures, accepted to Astrophysical Journal

Development of Tearing Instability in a Current Sheet Forming by Sheared Incompressible Flow

Authors: Elizabeth A. Tolman, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; R… ▽ More Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time-dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet. △ Less

Submitted 28 February, 2018; v1 submitted 24 November, 2017; originally announced November 2017.

Journal ref: Journal of Plasma Physics 84.1 (2018)

An Experimental Platform for Pulsed-Power Driven Magnetic Reconnection

Authors: J. D. Hare, L. G. Suttle, S. V. Lebedev, N. F. Loureiro, A. Ciardi, J. P. Chittenden, T. Clayson, S. J. Eardley, C. Garcia, J. W. D. Halliday, T. Robinson, R. A. Smith, N. Stuart, F. Suzuki-Vidal, E. R. Tubman

Abstract: We describe a versatile pulsed-power driven platform for magnetic reconnection experiments, based on exploding wire arrays driven in parallel [Suttle, L. G. et al. PRL, 116, 225001]. This platform produces inherently magnetised plasma flows for the duration of the generator current pulse (250 ns), resulting in a long-lasting reconnection layer. The layer exists for long enough to allow evolution o… ▽ More We describe a versatile pulsed-power driven platform for magnetic reconnection experiments, based on exploding wire arrays driven in parallel [Suttle, L. G. et al. PRL, 116, 225001]. This platform produces inherently magnetised plasma flows for the duration of the generator current pulse (250 ns), resulting in a long-lasting reconnection layer. The layer exists for long enough to allow evolution of complex processes such as plasmoid formation and movement to be diagnosed by a suite of high spatial and temporal resolution laser-based diagnostics. We can access a wide range of magnetic reconnection regimes by changing the wire material or moving the electrodes inside the wire arrays. We present results with aluminium and carbon wires, in which the parameters of the inflows and the layer which forms are significantly different. By moving the electrodes inside the wire arrays, we change how strongly the inflows are driven. This enables us to study both symmetric reconnection in a range of different regimes, and asymmetric reconnection. △ Less

Submitted 12 January, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: 14 pages, 9 figures. Version revised to include referee's comments. Submitted to Physics of Plasmas

Fully Kinetic Simulation of 3D Kinetic Alfven Turbulence

Authors: D. Groselj, A. Mallet, N. F. Loureiro, F. Jenko

Abstract: We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation… ▽ More We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation $k_{\parallel}(k_{\perp})$, where $k_{\parallel}$ is a characteristic wave number along the local mean magnetic field at perpendicular scale $l_{\perp}\sim 1/k_{\perp}$. The subion range anisotropy is scale dependent with $k_{\parallel}<k_{\perp}$ and the ratio of linear to nonlinear time scales is of order unity, suggesting that the kinetic cascade is close to a state of critical balance. Our results compare favorably against a number of \emph{in situ} solar wind observations and demonstrate---from first principles---the feasibility of plasma turbulence models based on a critically balanced cascade of kinetic Alfv\' en waves. △ Less

Submitted 7 March, 2018; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: revised manuscript; Fig. 4 was replaced to account for a small correction in the normalization of $k_{\parallel}$

Journal ref: Phys. Rev. Lett. 120, 105101 (2018)

A Drift-Kinetic Analytical Model for SOL Plasma Dynamics at Arbitrary Collisionality

Authors: R. Jorge, P. Ricci, N. F. Loureiro

Abstract: A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition… ▽ More A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously-derived fluid models. △ Less

Submitted 27 November, 2017; v1 submitted 4 September, 2017; originally announced September 2017.

Journal ref: Journal of Plasma Physics, 83(6) (2017)

The fully kinetic Biermann battery and associated generation of pressure anisotropy

Authors: K. M. Schoeffler, N. F. Loureiro, L. O. Silva

Abstract: The dynamical evolution of a fully kinetic, collisionless system with imposed background density and temperature gradients is investigated analytically. The temperature gradient leads to the generation of temperature anisotropy, with the temperature along the gradient becoming larger than that in the direction perpendicular to it. This causes the system to become unstable to pressure anisotropy dr… ▽ More The dynamical evolution of a fully kinetic, collisionless system with imposed background density and temperature gradients is investigated analytically. The temperature gradient leads to the generation of temperature anisotropy, with the temperature along the gradient becoming larger than that in the direction perpendicular to it. This causes the system to become unstable to pressure anisotropy driven instabilities, dominantly to electron Weibel. When both density and temperature gradients are present and non-parallel to each other, we obtain a Biermann-like linear in time magnetic field growth. Accompanying particle in cell numerical simulations are shown to confirm our analytical results. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 5 pages, 2 figures, + Supplementary materials (4 pages, 2 figures )

Journal ref: Phys. Rev. E 97, 033204 (2018)

Collisionless Reconnection in Magnetohydrodynamic and Kinetic Turbulence

Authors: N. F. Loureiro, S. Boldyrev

Abstract: It has recently been proposed (Loureiro & Boldyrev 2017; Mallet et al. 2017) that the inertial interval in magnetohydrodynamic (MHD) turbulence is terminated at small scales not by a Kolmogorov-like dissipation region, but rather by a new sub-inertial interval mediated by tearing instability. However, many astrophysical plasmas are nearly collisionless so that the MHD approximation is not applicab… ▽ More It has recently been proposed (Loureiro & Boldyrev 2017; Mallet et al. 2017) that the inertial interval in magnetohydrodynamic (MHD) turbulence is terminated at small scales not by a Kolmogorov-like dissipation region, but rather by a new sub-inertial interval mediated by tearing instability. However, many astrophysical plasmas are nearly collisionless so that the MHD approximation is not applicable to turbulence at small scales. In this Letter, we propose the extension of the theory of reconnection-mediated turbulence to plasmas which are so weakly collisional that the reconnection occurring in the turbulent eddies is caused by electron inertia rather than by resistivity. We find that the transition scale to reconnection-mediated turbulence depends on the plasma beta and on the assumptions of the plasma turbulence model. However, in all cases analyzed, the energy spectra in the reconnection-mediated interval range from $E(k)dk\propto k^{-8/3}dk$ to $E(k)dk\propto k^{-3}dk$. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: submitted for publication

Magnetohydrodynamic turbulence mediated by reconnection

Authors: Stanislav Boldyrev, Nuno F. Loureiro

Abstract: Magnetic field fluctuations in MHD turbulence can be viewed as current sheets that are progressively more anisotropic at smaller scales. As suggested by Loureiro & Boldyrev (2017) and Mallet et al (2017), below a certain critical thickness $λらむだ_c$ such current sheets become tearing-unstable. We propose that the tearing instability changes the effective alignment of the magnetic field lines in such a… ▽ More Magnetic field fluctuations in MHD turbulence can be viewed as current sheets that are progressively more anisotropic at smaller scales. As suggested by Loureiro & Boldyrev (2017) and Mallet et al (2017), below a certain critical thickness $λらむだ_c$ such current sheets become tearing-unstable. We propose that the tearing instability changes the effective alignment of the magnetic field lines in such a way as to balance the eddy turnover rate at all scales smaller than $λらむだ_c$. As a result, turbulent fluctuations become progressively less anisotropic at smaller scales, with the alignment angle increasing as $θしーた\sim (λらむだ/λらむだ_*)^{-4/5+βべーた}$, where $λらむだ_*\sim L_0 S_0^{-3/4}$ is the resistive dissipation scale. Here $L_0$ is the outer scale of the turbulence, $S_0$ is the corresponding Lundquist number, and {$0\leq βべーた<4/5$} is a parameter. The resulting Fourier energy spectrum is $E(k_\perp)\propto k_\perp^{-11/5+2βべーた/3}$, where $k_\perp$ is the wavenumber normal to the local mean magnetic field, and the critical scale is $λらむだ_c\sim S_L^{-(4-5βべーた)/(7-{20βべーた/3})}$. The simplest model corresponds to $βべーた=0$, in which case the predicted scaling formally agrees with one of the solutions obtained in (Mallet et al 2017) from a discrete hierarchical model of abruptly collapsing current sheets, an approach different and complementary to ours. We also show that the reconnection-mediated interval is non-universal with respect to the dissipation mechanism. Hyper-resistivity of the form ${\tilde ηいーた}k^{2+2s}$ leads (in the simplest case of $βべーた=0$) to the different transition scale $λらむだ_c\sim L_0{\tilde S}_0^{-4/(7+9s)}$ and the energy spectrum $E(k_\perp)\propto k_\perp^{-(11+9s)/(5+3s)}$, where ${\tilde S}_0$ is the corresponding hyper-resistive Lundquist number. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: submitted for publication

Fully kinetic versus reduced-kinetic modelling of collisionless plasma turbulence

Authors: D. Groselj, S. S. Cerri, A. Banon Navarro, C. Willmott, D. Told, N. F. Loureiro, F. Califano, F. Jenko

Abstract: We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully-kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta ($βべーた_i$) re… ▽ More We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully-kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta ($βべーた_i$) regimes are considered: 0.1 and 0.5. For $βべーた_i=0.5$, good agreement between the GK and FK models is found at scales ranging from the ion to the electron gyroradius, thus providing firm evidence for a kinetic Alfvén cascade scenario. In the same range, the HK model produces shallower spectral slopes, presumably due to the lack of electron Landau damping. For $βべーた_i=0.1$, a detailed analysis of spectral ratios reveals a slight disagreement between the GK and FK descriptions at kinetic scales, even though kinetic Alfvén fluctuations likely still play a significant role. The discrepancy can be traced back to scales above the ion gyroradius, where the FK and HK results seem to suggest the presence of fast magnetosonic and ion Bernstein modes in both plasma beta regimes, but with a more notable deviation from GK in the low-beta case. The identified practical limits and strengths of reduced-kinetic approximations, compared here against the fully-kinetic model on a case-by-case basis, may provide valuable insight into the main kinetic effects at play in turbulent collisionless plasmas, such as the solar wind. △ Less

Submitted 20 September, 2017; v1 submitted 8 June, 2017; originally announced June 2017.

Journal ref: Astrophys. J. 847, 28 (2017)

Formation and Structure of a Current Sheet in Pulsed-Power Driven Magnetic Reconnection Experiments

Authors: J. D. Hare, S. V. Lebedev, L. G. Suttle, N. F. Loureiro, A. Ciardi, G. C. Burdiak, J. P. Chittenden, T. Clayson, S. J. Eardley, C. Garcia, J. W. D. Halliday, N. Niasse, T. Robinson, R. A. Smith, N. Stuart, F. Suzuki-Vidal, G. F. Swadling, J. Ma, J. Wu

Abstract: We describe magnetic reconnection experiments using a new, pulsed-power driven experimental platform in which the inflows are super-sonic but sub-Alfvénic.The intrinsically magnetised plasma flows are long lasting, producing a well-defined reconnection layer that persists over many hydrodynamic time scales.The layer is diagnosed using a suite of high resolution laser based diagnostics which provid… ▽ More We describe magnetic reconnection experiments using a new, pulsed-power driven experimental platform in which the inflows are super-sonic but sub-Alfvénic.The intrinsically magnetised plasma flows are long lasting, producing a well-defined reconnection layer that persists over many hydrodynamic time scales.The layer is diagnosed using a suite of high resolution laser based diagnostics which provide measurements of the electron density, reconnecting magnetic field, inflow and outflow velocities and the electron and ion temperatures.Using these measurements we observe a balance between the power flow into and out of the layer, and we find that the heating rates for the electrons and ions are significantly in excess of the classical predictions. The formation of plasmoids is observed in laser interferometry and optical self-emission, and the magnetic O-point structure of these plasmoids is confirmed using magnetic probes. △ Less

Submitted 13 September, 2017; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: 14 pages, 12 figures. Accepted for publication in Physics of Plasmas

Role of magnetic reconnection in MHD turbulence

Authors: N. F. Loureiro, S. Boldyrev

Abstract: The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale… ▽ More The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale $λらむだ\sim L S_L^{-4/7}$, where $S_L$ is the outer-scale ($L$) Lundquist number and $λらむだ$ is the smallest of the field-perpendicular eddy dimensions. This scale is larger than the scale set by the resistive diffusion of eddies, therefore implying a fundamentally different route to energy dissipation than that predicted by the Kolmogorov-like phenomenology. In particular, our analysis predicts the existence of the sub-inertial, reconnection interval of MHD turbulence, with the Fourier energy spectrum $E(k_\perp)\propto k_\perp^{-5/2}$, where $k_\perp$ is the wave number perpendicular to the local mean magnetic field. The same calculation is also performed for high (perpendicular) magnetic Prandtl number plasmas ($Pm$), where the reconnection scale is found to be $λらむだ/L\sim S_L^{-4/7}Pm^{-2/7}$. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: submitted for publication

Journal ref: Phys. Rev. Lett. 118, 245101 (2017)

Anomalous Heating and Plasmoid Formation in a Driven Magnetic Reconnection Experiment

Authors: J. D. Hare, L. Suttle, S. V. Lebedev, N. F. Loureiro, A. Ciardi, G. C. Burdiak, J. P. Chittenden, T. Clayson, C. Garcia, N. Niasse, T. Robinson, R. A. Smith, N. Stuart, F. Suzuki-Vidal, G. F. Swadling, J. Ma, J. Wu, Q. Yang

Abstract: We present a detailed study of magnetic reconnection in a quasi-two-dimensional pulsed-power driven laboratory experiment. Oppositely directed magnetic fields $(B=3$ T), advected by supersonic, sub-Alfvénic carbon plasma flows $(V_{in}=50$ km/s), are brought together and mutually annihilate inside a thin current layer ($δでるた=0.6$ mm). Temporally and spatially resolved optical diagnostics, including i… ▽ More We present a detailed study of magnetic reconnection in a quasi-two-dimensional pulsed-power driven laboratory experiment. Oppositely directed magnetic fields $(B=3$ T), advected by supersonic, sub-Alfvénic carbon plasma flows $(V_{in}=50$ km/s), are brought together and mutually annihilate inside a thin current layer ($δでるた=0.6$ mm). Temporally and spatially resolved optical diagnostics, including interferometry, Faraday rotation imaging and Thomson scattering, allow us to determine the structure and dynamics of this layer, the nature of the inflows and outflows and the detailed energy partition during the reconnection process. We measure high electron and ion temperatures $(T_e=100$ eV, $T_i=600$ eV), far in excess of what can be attributed to classical (Spitzer) resistive and viscous dissipation. We observe the repeated formation and ejection of plasmoids, which we interpret as evidence of two-fluid effects in our experiment. △ Less

Submitted 25 January, 2017; v1 submitted 29 September, 2016; originally announced September 2016.

Comments: Accepted for publication in PRL, expected 17th Febuary 2017. 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 118, 085001 (2017)

Plasma Turbulence in the Scrape-off Layer of the ISTTOK Tokamak

Authors: Rogerio Jorge, Paolo Ricci, Federico D. Halpern, Nuno F. Loureiro, Carlos Silva

Abstract: The properties of plasma turbulence in a poloidally limited scrape-off layer (SOL) are addressed, with focus on ISTTOK, a large aspect ratio tokamak with a circular cross section. Theoretical investigations based on the drift-reduced Braginskii equations are carried out through linear calculations and non-linear simulations, in two- and three-dimensional geometries. The linear instabilities drivin… ▽ More The properties of plasma turbulence in a poloidally limited scrape-off layer (SOL) are addressed, with focus on ISTTOK, a large aspect ratio tokamak with a circular cross section. Theoretical investigations based on the drift-reduced Braginskii equations are carried out through linear calculations and non-linear simulations, in two- and three-dimensional geometries. The linear instabilities driving turbulence and the mechanisms that set the amplitude of turbulence as well as the SOL width are identified. A clear asymmetry is shown to exist between the low-field and the high-field sides of the machine. A comparison between experimental measurements and simulation results is presented. △ Less

Submitted 2 October, 2016; v1 submitted 30 June, 2016; originally announced June 2016.

Comments: 24 pages, 15 figures, Accepted for publication - Physics of Plasmas

Journal ref: Physics of Plasmas 23, 102511 (2016)

Comprehensive evaluation of the linear stability of Alfvén eigenmodes driven by alpha particles in an ITER baseline scenario

Authors: A. C. A. Figueiredo, P. Rodrigues, D. Borba, R. Coelho, L. Fazendeiro, J. Ferreira, N. F. Loureiro, F. Nabais, S. D. Pinches, A. R. Polevoi, S. E. Sharapov

Abstract: The linear stability of Alfvén eigenmodes in the presence of fusion-born alpha particles is thoroughly assessed for two variants of an ITER baseline scenario, which differ significantly in their core and pedestal temperatures. A systematic approach is used that considers all possible eigenmodes for a given magnetic equilibrium and determines their growth rates due to alpha-particle drive and Landa… ▽ More The linear stability of Alfvén eigenmodes in the presence of fusion-born alpha particles is thoroughly assessed for two variants of an ITER baseline scenario, which differ significantly in their core and pedestal temperatures. A systematic approach is used that considers all possible eigenmodes for a given magnetic equilibrium and determines their growth rates due to alpha-particle drive and Landau damping on fuel ions, helium ashes and electrons. This extensive stability study is efficiently conducted through the use of a specialized workflow that profits from the performance of the hybrid MHD drift-kinetic code $\mbox{CASTOR-K}$ (Borba D. and Kerner W. 1999 J. Comput. Phys. ${\bf 153}$ 101; Nabais F. ${\it et\,al}$ 2015 Plasma Sci. Technol. ${\bf 17}$ 89), which can rapidly evaluate the linear growth rate of an eigenmode. It is found that the fastest growing instabilities in the aforementioned ITER scenario are core-localized, low-shear toroidal Alfvén eigenmodes. The largest growth-rates occur in the scenario variant with higher core temperatures, which has the highest alpha-particle density and density gradient, for eigenmodes with toroidal mode numbers $n\approx30$. Although these eigenmodes suffer significant radiative damping, which is also evaluated, their growth rates remain larger than those of the most unstable eigenmodes found in the variant of the ITER baseline scenario with lower core temperatures, which have $n\approx15$ and are not affected by radiative damping. △ Less

Submitted 25 January, 2016; originally announced January 2016.

Sensitivity of alpha-particle-driven Alfven eigenmodes to q-profile variation in ITER scenarios

Authors: P. Rodrigues, A. C. A. Figueiredo, D. Borba, R. Coelho, L. Fazendeiro, J. Ferreira, N. F. Loureiro, F. Nabais, S. D. Pinches, A. R. Polevoi, S. E. Sharapov

Abstract: A perturbative hybrid ideal-MHD/drift-kinetic approach to assess the stability of alpha-particle-driven Alfven eigenmodes in burning plasmas is used to show that certain foreseen ITER scenarios, namely the Ip = 15 MA baseline scenario with very low and broad core magnetic shear, are sensitive to small changes in the background magnetic equilibrium. Slight variations (of the order of 1%) of the saf… ▽ More A perturbative hybrid ideal-MHD/drift-kinetic approach to assess the stability of alpha-particle-driven Alfven eigenmodes in burning plasmas is used to show that certain foreseen ITER scenarios, namely the Ip = 15 MA baseline scenario with very low and broad core magnetic shear, are sensitive to small changes in the background magnetic equilibrium. Slight variations (of the order of 1%) of the safety-factor value on axis are seen to cause large changes in the growth rate, toroidal mode number, and radial location of the most unstable eigenmodes found. The observed sensitivity is shown to proceed from the very low magnetic shear values attained throughout the plasma core, raising issues about reliable predictions of alpha-particle transport in burning plasmas. △ Less

Submitted 21 April, 2016; v1 submitted 7 January, 2016; originally announced January 2016.

Comments: 10 pages, 4 figures; Accepted manuscript

Journal ref: Nucl. Fusion 56 (2016) 112006

The generation of magnetic fields by the Biermann battery and the interplay with the Weibel instability

Authors: K. M. Schoeffler, N. F. Loureiro, R. A. Fonseca, L. O. Silva

Abstract: An investigation of magnetic fields generated in an expanding bubble of plasma with misaligned temperature and density gradients (driving the Biermann battery mechanism) is performed. With gradient scales $L$, large-scale magnetic fields are generated by the Biermann battery mechanism with plasma $βべーた\sim 1$, as long as $L$ is comparable to the ion inertial length $d_i$. For larger system sizes,… ▽ More An investigation of magnetic fields generated in an expanding bubble of plasma with misaligned temperature and density gradients (driving the Biermann battery mechanism) is performed. With gradient scales $L$, large-scale magnetic fields are generated by the Biermann battery mechanism with plasma $βべーた\sim 1$, as long as $L$ is comparable to the ion inertial length $d_i$. For larger system sizes, $L/d_e > 100$ (where $d_e$ is the electron inertial length), the Weibel instability generates magnetic fields of similar magnitude but with wavenumber $k d_e \sim 0.2$. In both cases, the growth and saturation of these fields have a weak dependence on mass ratio $m_i/m_e$, indicating electron mediated physics. A scan in system size is performed at $m_i/m_e = 2000$, showing agreement with previous results with $m_i/m_e = 25$. In addition, the instability found at large system sizes is quantitatively demonstrated to be the Weibel instability. Furthermore, magnetic and electric energy spectra at scales below the electron Larmor radius are found to exhibit power law behavior with spectral indices $-16/3$ and $-7/3$, respectively. △ Less

Submitted 5 May, 2016; v1 submitted 16 December, 2015; originally announced December 2015.

Journal ref: Phys. Plasmas 23, 056304 (2016)

Magnetic reconnection: from the Sweet-Parker model to stochastic plasmoid chains

Authors: N. F. Loureiro, D. A. Uzdensky

Abstract: (abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergo… ▽ More (abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergone a major paradigm shift. The steady-state model of reconnection described by the famous Sweet-Parker (SP) theory, which dominated the field for ~50 years, has been replaced with an essentially time-dependent, bursty picture of the reconnection layer, dominated by the continuous formation and ejection of multiple secondary islands (plasmoids). Whereas in the SP model reconnection was predicted to be slow, a major implication of this new paradigm is that reconnection in fluid systems is fast (i.e., independent of the Lundquist number), provided that the system is large enough. This conceptual shift hinges on the realization that SP-like current layers are violently unstable to the plasmoid instability - implying, therefore, that such current sheets are super-critically unstable and thus can never form in the first place. This suggests that the formation of a current sheet and the subsequent reconnection process cannot be decoupled, as is commonly assumed. This paper provides an introductory-level overview of the recent developments in reconnection theory and simulations that led to this essentially new framework. We briefly discuss the role played by the plasmoid instability in selected applications, and describe some of the outstanding challenges that remain at the frontier of this subject. Amongst these are the analytical and numerical extension of the plasmoid instability to (i) 3D and (ii) non-MHD regimes. New results are reported in both cases. △ Less

Submitted 28 July, 2015; originally announced July 2015.

Comments: Submitted to Plasma Phys. Control. Fusion, EPS 2015 Special Issue

Viriato: a Fourier-Hermite spectral code for strongly magnetised fluid-kinetic plasma dynamics

Authors: N. F. Loureiro, W. Dorland, L. Fazendeiro, A. Kanekar, A. Mallet, M. S. Vilelas, A. Zocco

Abstract: We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihi… ▽ More We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes. △ Less

Submitted 11 May, 2015; originally announced May 2015.

Comments: 42 pages, 15 figures, submitted to J. Comp. Phys

Kinetic microtearing modes and reconnecting modes in strongly magnetised slab plasmas

Authors: A Zocco, N F Loureiro, D Dickinson, R Numata, C M Roach

Abstract: The problem of the linear microtearing mode in a slab magnetised plasma, and its connection to kinetic reconnecting modes, is addressed. Electrons are described using a novel hybrid fluid-kinetic model that captures electron heating, ions are gyrokinetic. Magnetic reconnection can occur as a result of either electron conductivity and inertia, depending on which one predominates. We eschew the use… ▽ More The problem of the linear microtearing mode in a slab magnetised plasma, and its connection to kinetic reconnecting modes, is addressed. Electrons are described using a novel hybrid fluid-kinetic model that captures electron heating, ions are gyrokinetic. Magnetic reconnection can occur as a result of either electron conductivity and inertia, depending on which one predominates. We eschew the use of an energy dependent collision frequency in the collisional operator model, unlike previous works. A model of the electron conductivity that matches the weakly collisional regime to the exact Landau result at zero collisionality and gives the correct electron isothermal response far from the reconnection region is presented. We identify in the breaking of the constant-$A_{\parallel}$ approximation the necessary condition for microtearing instability in the collisional regime. Connections with the theory of collisional non-isothermal (or semicollisional) and collisionless tearing-parity electron temperature gradient driven (ETG) modes are elucidated. △ Less

Submitted 20 November, 2014; originally announced November 2014.

Comments: Submitted for publication