We investigate the self-consistent particle acceleration physics associated with the development of the kink instability (KI) in nonrelativistic, electron-ion plasma jets. Using 3D fully kinetic particle-in-cell simulations, we show that the KI efficiently converts the initial toroidal magnetic field energy into energetic ions. The accelerated ions form a nonthermal power-law tail in the energy spectrum, containing of the initial magnetic field energy, with the maximum ion energy extending to the confinement energy of the jet. We find that the ions are efficiently accelerated by the concerted action of the motional electric field and the highly tangled magnetic field that develop in the nonlinear phase of the KI: fast curvature drift motions of ions across magnetic field lines enable their acceleration along the electric field. We further investigate the role of Coulomb collisions in the ion acceleration efficiency and identify the collisional threshold above which nonthermal ion acceleration is suppressed. Our results reveal how energetic ions may result from unstable nonrelativistic plasma jets in space and astrophysics and provide constraints on the plasma conditions required to reproduce this acceleration mechanism in laboratory experiments.
I. INTRODUCTION
Jets of plasma threaded by magnetic fields are ubiquitous in space and astrophysical environments. Throughout their evolution, these jets can be subject to current-driven magnetohydrodynamic (MHD) instabilities, with the most common being the sausage (m = 0) and kink (m = 1) modes.1 The development of such instabilities is believed to be a candidate mechanism to explain the conversion of the jet's magnetic energy into energetic particles and radiation in a wide range of space physics and astrophysics scenarios. These include the flaring activity and energetic particle generation in jets emanating from the solar atmosphere,2–4 the synchrotron emission from bright knots, and cosmic ray acceleration in jets from active galaxies,5–7 among others. These instabilities are also associated with violent disruptions and particle acceleration in laboratory fusion devices, such as Tokamaks8 and Z-pinches.9–11 Yet, the detailed physics underlying how the development of such instabilities in jets accelerates particles remains poorly understood.
MHD simulations play an important role in understanding where and when sausage and kink modes develop over the course of the formation and propagation of astrophysical6,12 and laboratory13 plasma jets and how their nonlinear development impacts the structural integrity of jets.14 However, the kinetic physics that underpins the acceleration of particles is not captured by MHD models. Simulations of test particle dynamics in the fields computed by MHD simulations have shown nonthermal acceleration15 but do not provide a self-consistent picture of the injection mechanisms or feedback of energetic particles on the plasma fields. Fully kinetic simulations are necessary to address this central problem in a self-consistent manner. Using 3D particle-in-cell (PIC) simulations, we have recently demonstrated the self-consistent acceleration of nonthermal particles resulting from the development of the kink instability (KI) in relativistic magnetized pair plasma and pair-proton jets.7 These simulations revealed that the KI provides a viable means of producing radiating nonthermal particle distributions and of accelerating ultrahigh-energy cosmic rays in the jets of active galactic nuclei. These findings encourage further exploration of the KI's ability to efficiently accelerate particles in nonrelativistic electron-ion plasma jets of relevance to a range of space and laboratory conditions.
In this paper, we report on 3D fully kinetic simulations of the KI in nonrelativistic, electron-ion plasma jets and discuss the associated particle acceleration physics. We demonstrate that the nonlinear development of the KI in collisionless plasma results in the efficient dissipation of the jet's magnetic field energy, which is preferentially transferred to high-energy ions. Approximately 10% of the initial magnetic field energy is converted into nonthermal ions. The ion energy spectrum develops a power-law tail , with p approximately fixed at over the parameter range explored in this work. We discuss the ion acceleration mechanism and how it depends on the size and magnetization of the jet and how the mass ratio between electrons and ions dictates the energy partition between species. In addition, we explore how the plasma collisionality impacts particle acceleration, which is relevant for a variety of laboratory experiments. Previous experiments using plasma guns,16,17 modified plasma Z-pinch platforms,11,18–20 and high-power lasers21 have successfully studied MHD instabilities in jets. In some cases, the observation of nonthermal ions has been reported,11 yet without a clear understanding of the underlying acceleration mechanism. Our work helps establish the conditions for which the present acceleration mechanism can be studied in laboratory experiments.
This paper is organized as follows. In Sec. II, we describe the initial physical configuration of the plasma jet and our simulation setup. We then discuss the evolution of the electric and magnetic fields produced by the KI and the associated particle acceleration dynamics in the collisionless regime in Sec. III. The effects of Coulomb collisions on the particle acceleration dynamics are discussed in Sec. IV, and conclusions are drawn in Sec. V.
II. SIMULATION SETUP
We simulate the particle acceleration dynamics associated with the development of the KI in a nonrelativistic electron-ion plasma jet using the 3D PIC code OSIRIS 4.0.22,23 We consider a purely toroidal magnetic field profile of the form , where B0 is the peak amplitude of the magnetic field and R is the characteristic radius of the jet. The current density is supported by symmetrically streaming electrons and ions along the jet axis. The thermal pressure profile of the plasma (P) is chosen to achieve hydromagnetic equilibrium, ; we consider plasma with uniform number density and a nonuniform temperature profile with such that . Note that we have tested different initial conditions with jet current being carried solely by the electrons and no net flow of ions and have found that the particle acceleration dynamics remains unaffected. This is to be expected since the particle thermal velocities are significantly higher than the flow velocities required to support the electric current.
The jet's initial conditions can thus be characterized by the following parameters. A dimensionless measure of the jet's magnetization is , where is the ion cyclotron frequency and is the ion plasma frequency; the dimensionless jet radius is , with being the ion inertial length; and the ion to electron mass ratio is mi/me. Note that due to pressure balance , we have that the typical Larmor radius of thermal ions at the core of the jet is given by .
We explore jet magnetizations in the range of [corresponding to nonrelativistic magnetic energy densities in the range of 0.006–0.06], and we simulate jet radii in the range of R/di = 5–10. Given the high computational cost of simulating realistic ions to electron mass ratios, we perform simulations with reduced mass ratios mi/me between 4 and 36. By progressively increasing the mass ratio, we uncover how the particle acceleration physics scales with mi/me, allowing us to infer the behavior for realistic mass-ratio conditions.
The dimensions of our simulation domain are , with the jet placed at the center of the domain and oriented along . The grid resolution is chosen to resolve the gyroradius of thermal electrons at the core of the jet with >2 points, and we use 12 particles per cell per species with quadratic splines for the particle shapes. Finer spatial resolutions and more particles per cell were tested to ensure numerical convergence. Periodic boundary conditions are imposed in all directions; the transverse ( and ) dimensions of the domain are sufficiently large to avoid artificial recirculation of particles during the time scale of interest. We follow the evolution of the system until the particle acceleration dynamics terminates, which occurs at is the characteristic Alfvén speed).
III. PARTICLE ACCELERATION IN COLLISIONLESS ELECTRON-ION JETS
A. Structure and evolution of the KI-induced E and B fields
The general dynamics of the KI in collisionless electron-ion plasma is illustrated in Fig. 1 with the case of , and ; all other simulations explored in this work exhibited similar behavior. The onset of the KI is triggered by small transverse displacements of the jet about its equilibrium. Those displacements produce imbalances in the magnetic pressure across the jet, inducing motions that reinforce the initial distortion. This ultimately manifests as a growing helical modulation of the jet structure (left column of Fig. 1), with a wavelength of and developing at a rate of , in agreement with linear MHD theory.1
The radial motions of the magnetized plasma jet give rise to an inductive electric field , which at early times is also shown to have a helical and harmonic structure in Fig. 1(c1). Similar to the behavior found in the relativistic pair plasma regime,7,24 the nonlinear distortions of the jet (distortions comparable to the jet radius R) result in the formation of a coherent structure in the axial component of the electric field Ez [Fig. 1(c2)]; the mean amplitude of the electric field along the jet axis is found to be . The nonlinear distortions of the jet current density also result in a highly tangled magnetic field structure [Fig. 1(b2)]. It is this configuration of electric and magnetic fields that mediates the efficient conversion of magnetic energy [≃65% of initial magnetic energy ] into plasma kinetic energy [Fig. 2(a)]. This efficient accelerating field structure persists during the transit time of the kink perturbations across the diameter of the jet, . After this period, the fields decay and particle acceleration ceases.
B. Efficient acceleration of nonthermal ions
As shown in Fig. 2(a), the dissipated magnetic energy is preferentially transferred to the jet ions, gaining as much energy as the electron population. More interestingly, this process results in the acceleration of nonthermal ions, forming a high-energy power-law tail in the ion energy spectrum with p = 2.6 [Fig. 2(b1)]. We find that of the total initial magnetic field energy is transferred to nonthermal ions with . We further verify that ions are accelerated up to the confinement energy of the jet (valid for subrelativistic ion energies), i.e., the energy beyond which the ion Larmor radius ρi exceeds the system size R. In the astrophysical context, this limiting energy is also known as the Hillas energy or the Hillas constraint.25 For the simulated parameters, we find that the cutoff energy of the ion spectrum occurs at . This is consistent with the maximum energy gain of an ion accelerating freely in the coherent axial electric field structure during the time scale τKI, .
The jet electrons, however, are not accelerated to nonthermal energies by the KI [Fig. 2(b2)]. In fact, the electrons within the jet (within r < R) cool down as shown by the receding spectral tail in Fig. 2(b2). It is the electrons at the periphery of the jet () that absorb a fraction of the dissipated magnetic energy but are not accelerated efficiently to many times their initial thermal energy. The difference between the acceleration dynamics of ions and electrons will be further discussed in Sec. III C.
C. Ion acceleration mechanism
We find that the mechanism by which nonthermal ions are accelerated is similar to that identified in the regime of relativistic pair plasma.7 By following the trajectories of a representative sample of nonthermal ions [gray curves in Fig. 3(a)], we have verified that their acceleration is due to the work of the motional electric field . This is shown for both the ion that attained the highest energy in the sample (red curve) and the mean energy gain of the sample (blue curve) in Fig. 3(a). This indicates that nonideal electric fields (), commonly associated with reconnecting current layers, play a negligible role in the acceleration dynamics of ions.
The efficient acceleration of ions by the motional electric field implies that ions must move efficiently transverse to the local magnetic field. Indeed, we find that the highly inhomogeneous and highly tangled structure of the magnetic field that develops in the nonlinear phase of the KI [Fig. 1(b2)] facilitates the displacement of ions across magnetic field lines via guiding center drift motions. In particular, we find that curvature drift motions of the ions play a dominant role. This is verified in Fig. 3(b), which shows that the energy gain of a random subset of 3 particle tracks (from the same sample) is well described by the work performed by the motional electric field along their guiding center curvature drift trajectories, i.e., , where is the curvature drift velocity and is the magnetic curvature vector field. The guiding center description of the ion trajectories breaks down when they achieve a large fraction of the confinement energy of the jet, . Beyond these energies, ions become effectively unmagnetized, moving nearly freely along the motion electric field in meandering type orbits until they escape the acceleration region or until the accelerating fields decay.
Figure 3(c) shows the evolution of the distribution of magnetic curvature κB during the nonlinear development of the KI. Enhancement in the distribution of magnetic curvature is observed in between the scales of the jet radius and Larmor radius of thermal electrons . However, this enhancement is significantly larger at the Larmor scale of thermal ions when compared to that of thermal electrons. This result indicates that thermal ions are more likely to encounter the magnetic field curvature that approaches the scale of their Larmor radii than electrons, allowing them to be more easily injected into a rapid accelerating phase via curvature drifts. Moreover, ions can become locally unmagnetized (when they experience ) with higher probability than electrons. It is likely that this distribution of magnetic curvature is key to explaining the difference in acceleration efficiency between the two species. A more in detailed analysis of the electron dynamics and acceleration efficiency for different initial conditions will be the subject of future work.
D. Particle acceleration dependence on magnetization (), system size , and mass ratio mi/me
We have performed a set of simulations of varying magnetization (), system sizes (), and ion to electron mass ratios () to probe how these parameters influence the particle acceleration dynamics of the KI. For the range of parameters explored in this work, we have verified that all simulations exhibit identical growth and saturation of the KI upon normalizing space and time to R and , respectively. We systematically observe the conversion of of the initial magnetic energy in the system into plasma kinetic energy after in all our simulations. The magnitude of the accelerating electric field structure that develops in the nonlinear phase of the KI is also similar across our simulations when normalized to the characteristic value of ().
We observe the persistent acceleration of nonthermal ions in all our simulations, while electrons are always found to remain thermal. The spectra of accelerated ions for varying mi/me and at a fixed system size are shown in Fig. 4(a). Interestingly, we find that the power-law index of the nonthermal tail remains approximately constant at −2.6 ± 0.1 for varying mi/me and (within the explored parameter range), and thus, these results may be extrapolated to realistic mass ratios. We further verify in all simulations that the nonthermal ion tail extends up to the confinement energy of the jet [Fig. 4(a)]. For a fixed normalized jet radius , we find that , which is independent of the ion to electron mass ratio and jet magnetization, as seen in Fig. 4(a). In all the simulated cases, we find that the cutoff energy of the spectrum is .
While the fraction of initial magnetic energy that is dissipated is approximately constant in all our simulations, the ratio of the electron to ion energy gain () is shown to depend on both the mass ratio and the magnetization, as illustrated in Fig. 4(b). We observe that an increasing fraction of the dissipated magnetic energy is transferred to the ions with both the increasing mass ratio and increasing magnetization. We note that a similar dependence of on mi/me has been reported in 3D PIC simulations of relativistic magnetized plasma turbulence.26
IV. ROLE OF COULOMB COLLISIONS IN ION ACCELERATION EFFICIENCY
There is a significant experimental effort to study the MHD stability of plasma jets in the laboratory in conditions relevant to astrophysical environments.11,16–21 However, laboratory-produced plasma jets are generally far more collisional than those naturally occurring in space or astrophysical settings. Thus, in order to assess if the ion acceleration mechanism discussed in Sec. III can be studied in laboratory plasma experiments, it is important to consider the effects of Coulomb collisions. One naturally anticipates that if the jet plasma is too collisional, nonthermal particle distributions will not be allowed to develop as they will rapidly relax to thermal equilibrium.
We have investigated the impact of a finite Coulomb collision frequency on the particle acceleration dynamics of the KI, with the aim of determining the collisional threshold above which nonthermal particle acceleration is suppressed. We include the physics of binary Coulomb collisions (the Spitzer-Harm model for weakly collisional plasmas) in our PIC simulations using a Monte Carlo approach.27–30 This method randomly pairs particles locally within a cell and scatters their momenta such that energy and momentum are conserved on each collision (equal particle weights are used). We use this method to model electron-ion and ion-ion collisions in our simulations of the KI, and we vary the collision frequency by artificially varying the Coulomb logarithm. In this way, all other physical parameters of the system can be fixed (plasma density, temperature, and magnetic field strength) and the effect of Coulomb collisions can be varied independently.
We have examined the effect of Coulomb collisions on the illustrative case explored in Sec. III, i.e., a jet with , and (note that is close to the jet radii produced in recent experiments11,19,21). We progressively increased the ion collision frequency νi relative to the dynamical time of the KI (which also corresponds to the acceleration time scale of the KI) from to 0.4; note that can also be conveniently interpreted as the ratio of the jet radius to the ion mean-free-path, , since (as imposed by the initial equilibrium configuration). The resulting spectra of accelerated ions are shown in Fig. 5(a) for varying νi values. We find that the spectrum of accelerated ions remains unaffected for , relative to the collisionless case: the slope of the nonthermal tail, the number of nonthermal particles, and the maximum energy reached remain unchanged by Coulomb collisions below . At , the ion collision frequency is no longer negligible compared to the acceleration time scale of the KI. As ions gain energy through the KI-induced electric field, they also lose energy through collisions and heat up the background particles. This effect leads to a strong reduction in the acceleration efficiency of nonthermal ions as shown in Fig. 5(a): the number of nonthermal particles is strongly reduced and the slope of the tail of the distribution is hardened. At , nonthermal ion acceleration is nearly fully suppressed.
Note that as the collision frequency was varied, no significant changes were observed in the morphology or temporal evolution of the KI-induced electric and magnetic fields relative to the collisionless regime. This is illustrated Figs. 5(b1)–5(b3) by the snapshots of the cross section of Ez for varying collisionality. These snapshots were all taken at , revealing that the temporal development and spatial structure of the KI-induced electric field are nearly unaffected by collisions; only in the highly collisional case of does the amplitude of Ez become noticeably lower, mainly as a consequence of the increased magnetic diffusivity that lowers the peak amplitude of the magnetic field and hence of the motional electric field.
The laboratory plasma jets produced by the radial wire array Z-pinch19 are characterized by low-temperature and high-Z plasmas, achieving . The acceleration of nonthermal ions by the mechanism reported here is thus expected to be significantly suppressed within the dense body of the jet. However, high-energy ions have indeed been reported, resulting from the development of MHD instabilities in the jets produced by this platform.11 It is likely these high-energy ions are accelerated outside of the dense jet in the low-density ambient plasma where the collisionality is reduced. Our simulations have not considered the role of a low-density plasma background surrounding the jet, which could allow ambient particles to interact with the fields produced by the KI of the jet while remaining in a weakly collisional environment. Indeed, we do see in our simulations that particles outside of the jet, within , do interact with the KI-induced fields and absorb a significant fraction of the dissipated magnetic energy. It is therefore possible that the mechanism reported here can participate in the acceleration of nonthermal particles in the low-density plasma surrounding the dense jet. A more detailed exploration of our work in the conditions of these experiments will be pursued in the near future.
Laser-driven high-energy-density plasmas are more likely to produce weakly collisional jets in the laboratory. Recent laser-driven plasma experiments produced near kilo-electron-volt-temperature plasma jets, with , aimed at investigating the development of the KI in astrophysically relevant conditions.21 Future experiments using a more energetic laser drive may be able to produce jets with , making it possible to directly probe the efficient particle acceleration physics of the KI within the body of the jet, as described in this work.
V. CONCLUSIONS
We have shown via 3D PIC simulations that the development of the KI in nonrelativistic, electron-ion plasma results in the efficient acceleration of nonthermal ions. Approximately 10% of initial magnetic energy is transferred into nonthermal ions over the course of a few 10s of dynamical times of the KI (), forming a power-law tail in the energy spectrum with index . We showed that the power-law index of the nonthermal tail remains nearly constant over the range of jet magnetizations explored in our work and that the maximum ion energy systematically reaches the confinement energy of the jet .
The ion acceleration mechanism is similar to the mechanism found to operate in relativistic pair plasma jets.7 The nonlinear development of the KI produces a coherent motional electric field along the axis of the jet that is embedded in a highly tangled magnetic field. Ions experience fast curvature drift motions in the highly tangled magnetic field, which permit their efficient displacement along the motional electric field and hence their efficient acceleration. Our results indicate that this is a viable mechanism to explain efficient acceleration of nonthermal protons in nonrelativistic jets in space and astrophysical environments.
The similarity between the particle acceleration physics in the relativistic and nonrelativistic regimes of the KI should also motivate the development of laboratory experimental platforms capable of investigating the particle acceleration mechanisms relevant to astrophysical jets. Our work indicates that in order to study this acceleration mechanism in the laboratory, the plasma collisionality needs to be significantly reduced, such that the ion mean free path becomes greater than the jet radius. These conditions can likely be produced in laser-driven high-energy-density plasma experiments.31
ACKNOWLEDGMENTS
This work was supported by the U.S. Department of Energy SLAC Contract No. DE-AC02-76SF00515, by the U.S. DOE Office of Science, Fusion Energy Sciences under FWP 100237, and by the U.S. DOE Early Career Research Program under FWP 100331. The authors acknowledge the OSIRIS Consortium, consisting of UCLA and IST (Portugal) for the use of the OSIRIS 4.0 framework and the visXD framework. Simulations were run on Mira (ALCF) through an ALCC award.