Recent Publications

Assessing the Projector Augmented-Wave Method for Stopping Power Calculations

Bryn Lloyd, Dirk O. Gericke, Gilles Rodway-Gant, Gianluca Gregori

The stopping power of charged particles is investigated using time-dependent density functional theory (TDDFT). Such simulations are made possible by recent advances in computational resources and numerical implementations of this first-principles method. In practice, DFT simulations widely employ the projector augmented-wave (PAW) method to approximate all-electron behaviour, but the implications of the PAW approximation for non-adiabatic stopping simulations remain insufficiently explored. Here, the suitability of the PAW method for stopping power simulations is evaluated. A workflow for generating and selecting PAW datasets tailored to these simulations is developed, enabling systematic optimisation of augmentation radii and projector constructions. The approach is applied to proton stopping in FCC aluminium, demonstrating how dataset design influences stopping predictions, and enabling an investigation of crystal channelling effects on charged-particle transport.

Proton Radiography with Isotrajectory Optics

Dmitry Varentsov, Maksym Miski-Oglu, Zs. Major, Paul Neumayer, Martin Schanz, John Schmidt

We consider a compact isotrajectory-optics scheme for proton radiography. The calculation uses a 30-60 MeV proton interval as a reference example, typical of high-energy short-pulse laser acceleration in the TNSA regime. The idea is to use synchronized time-dependent fields such that, for one selected energy-arrival-time branch, the proton energy changes mainly the arrival time and not the transverse image position. The reference lattice consists of four pulsed electric quadrupoles with horizontal sign pattern + - + -, magnification Mx = My = -3.0 at a detector plane of 2.57 m, a common Fourier plane at zF = 0.962 m, and equal outer and equal inner lens lengths. The ideal closure residual is calculated from second central moments. Fourth central moments are kept separately as a diagnostic of non-Gaussian tails. The quoted residual widths are numerical closure tests of the ideal map, not experimental resolution predictions. Real electrode fields, waveform errors, scattering, detector blur, and space charge define the next level of the problem.

Microcoulomb-level electron beam and multi-Joule hard X-rays driven by a high-efficiency laser-plasma accelerator

B. Mahieu, L. Ribotte, W. Cayzac, G. Boutoux, R. Parreault, J. Gastineau, E. Lamoine, F. Audo, R. Babjak, D. Batani, et al.

We report on the production of ultrahigh-charge relativistic electron beams and the development of a laser-wakefield acceleration platform at the LMJ facility. Making use of the kilojoule-class, sub-picosecond PETAL laser pulse focused onto a supersonic helium gas jet, electron beams carrying a total charge beyond 1 $μ$C were generated, with energies up to $\sim$500 MeV. Given the ps-scale laser pulse duration, an on-target intensity approaching $10^{19}~\mathrm{W/cm^2}$, and a plasma density reaching 2% of the critical density, electron energisation arises from a combination of self-modulated laser wakefield acceleration (SMLWFA) and direct laser acceleration (DLA). The resulting electron spectrum exhibits a Maxwellian-like distribution, characteristic of this mixed SMLWFA/DLA regime. The total energy carried by the electron beam is estimated to be up to 17 J, within a sub-ps duration. A broadband Joule-level photon beam was also produced by Bremsstrahlung, demonstrating the potential for future applications. Experimental results are supported by start-to-end numerical simulations, including 3-D particle-in-cell and Monte-Carlo particle transport calculations. These findings pave the way for applications requiring high-charge electron beams, including the generation of high-power secondary radiation or particle sources. The use of these beams to probe matter in high-energy density states driven by the nanosecond-duration LMJ beams represents another promising avenue.

Data-Driven Reconstruction of Spatially Resolved Electron and Ion Energy Distributions from Macroscopic Plasma Quantities with Deep Neural Networks

Libin Varghese, Kaushik Prajapati, Bhaskar Chaudhury

Spatially resolved EEDFs/IEDFs provide essential kinetic information about low-temperature plasmas (LTPs) and play a central role in determining transport, chemical reaction rates, and plasma surface interactions. While kinetic simulations directly resolve these distributions, experimental measurements remain challenging and are often invasive, spatially limited, or require assumptions regarding the distribution shape such as a Maxwellian. However, several macroscopic plasma observables can be measured non-invasively using advanced diagnostic techniques, providing spatially resolved information about the plasma state. An important inverse problem is therefore whether readily measurable macroscopic plasma quantities contain sufficient information to reconstruct the underlying kinetic state. In this work, we investigate this problem by learning a nonlinear mapping from spatially resolved macroscopic plasma observables to the corresponding spatially resolved EEDFs/IEDFs using a deep learning framework. Paired datasets comprising 2D macroscopic observables and spatially resolved EDFs are generated using 2D-3V PIC-MCC simulations. Three representative learning paradigms, a U-Net, a FNO, and a MeshGraphNet, are employed in this study to learn this inverse mapping. The predicted EDFs reproduce both bulk plasma and sheath characteristics with good agreement to the PIC-MCC reference data, with the FNO providing the best overall performance. Beyond conventional metrics, physics-based validation demonstrates that the reconstructed EDFs accurately recover the corresponding density and temperature, and rate coefficients. These results demonstrate that macroscopic plasma observables encode sufficient information to infer important kinetic properties in LTPs, providing a potential foundation for surrogate kinetic modeling and next-generation plasma diagnostics.

Modelling of Flowing Plasma in the Magnetic Field of the Small Volume Plasma System Experiment

Subhasish Bag, Ashish Ganguli, Vikrant Saxena, Ramesh Narayanan, Debaprasad Sahu

A flow model for a magnetized plasma has been developed to investigate the flow dynamics in the small volume plasma system (SVPS) experiment. The SVPS experimental conditions require the model to describe a stationary, collisional, quasineutral, axisymmetric plasma. Also, the ions are cold while the electrons are isothermal and in thermal equilibrium, obeying the Boltzmann relation. In a plasma flowing along a magnetic field, the velocity of ions along the magnetic field lines is much greater than the velocity perpendicular to the field. The latter feature permits a unique ordering of the relevant variables, when the flow equations are transformed to the magnetic coordinate system (MCS), where the coordinate axes are parallel and perpendicular to the field lines. The ordering of the flow variables in the MCS allows a further simplification of the flow equations, by permitting their splitting into set of reduced, simplified equations. The SVPS experimental data are used to provide the requisite boundary conditions for initializing and solving the reduced flow equations on a magnetic coordinate grid along the different lines of the MCS. An important aspect of the present work is the validation of the splitting scheme used to derive the simplified and reduced flow equations. This is achieved by an in-depth comparison of the predictions from the model equations with the experimental data. The obtained numerical results compare favourably with the SVPS observations and have been discussed rigorously. The model developed here provides a framework for exploring magnetized plasma dynamics in the given cylindrically symmetric magnetic field configuration and can be further extended to more complex configurations.

Minimum-q induced alternation between infernal modes and EP-driven modes in advanced tokamak configurations

Shiwei Xue, Ping Zhu, Haolong Li

For an advanced tokamak configuration in the presence of energetic particles (EPs), the dominant instability is found to alternate between infernal modes and Alfvén eigenmodes with the variation of the minimum safety factor $q_{\min}$. For relatively high $q_{\min}$, the mode is identified as a reversed-shear Alfvén eigenmode (RSAE), characterized by its finite Alfvénic frequency and radial localization near the minimum of safety-factor profile. As $q_{\min}$ is further reduced, the dominant branch sequentially transitions through a low-frequency infernal-mode interval, then an energetic-particle-mode (EPM) regime, and finally another low-frequency infernal-mode interval. Increasing the EP beta fraction $β_h$ tends to destabilize the RSAE and EPM branches but to stabilize the infernal modes. Phase-space diagnostics further indicate that the destabilizing effects of EPs on the RSAE and EPM branches are mainly associated with trapped-particle drive, whereas the stabilizing effects of EPs on the infernal modes is dominated by passing particles.

Hybrid kinetic-MHD simulation on the formation of runaway electron current plateau during a current quench process

Weimin Fu, Ping Zhu

Runaway electron (RE) generation is of great concern for high-current tokamak operations. A full-f particle-in-cell (PIC) model for RE dynamics has been developed and coupled with the 3D nonlinear extended magnetohydrodynamic (MHD) model implemented in the NIMROD code. Our model accounts for RE generation using analytical source terms and advances the RE motion along guiding-center (GC) orbits. The model was employed to simulate the formation of RE current plateau during a disruption, in which the plasma current becomes dominated by the RE current. The model agrees well with several RE codes based on fluid models when the RE GC drifts are ignored. For highly relativistic REs, the grad-$B$ and curvature drifts can play a significant role in the RE generation and motion due to the increase in the safety factor profile during the current quench process.

Electron Injection and Beam Dynamics in a Laser Wakefield Acceleration Driven by Laser Pulses Carrying Orbital Angular Momentum

Q. Qian, M. Li, Y. Ma, P. T. Campbell, C-H. Chang, E. Denis, N. Ernst, R. Fedosejevs, B. Hou, A. James, et al.

Laser pulses carrying orbital angular momentum (OAM) provide a new degree of freedom for controlling plasma-based accelerators. Here, we demonstrate experimentally laser wakefield acceler- ation driven by OAM laser pulses, producing electron beams with distinct broad-band energy spectra and angularly structured phase space. The measured spectra exhibit a two-beamlet structure with correlated angular dispersion, indicating injection occurring at multiple azimuthal phases of the plasma wake. Particle-in-cell simulations reproduced the observed spectral features and revealed helical electron trajectories driven by the OAM-induced wakefields. These results show that the phase structure of the laser driver can shape electron injection and acceleration dynamics, opening a route toward optimal control of beam structure in compact laser-driven accelerators.

A self-consistent solar coronal heating model by Alfvenic waves

Hidetaka Kuniyoshi, Richard J. Morton

Alfvenic waves are prevalent throughout the solar atmosphere and are believed to play an essential role in coronal heating, classified as alternating current (AC) heating in contrast to direct current (DC) heating associated with quasi-static magnetic field line braiding. The relative importance of AC versus DC heating depends on the details of the photospheric driver and on the configuration of the magnetic field. Moreover, even if AC heating prevails, several wave dissipation mechanisms have been proposed, and which of them dominates remains unclear, as its efficiency depends on plasma compressibility and density inhomogeneity. We address these issues by performing three-dimensional radiative magnetohydrodynamic (MHD) simulations of a coronal loop spanning from the upper convection zone to the corona, which self-consistently capture many relevant physical processes. We find that the corona is predominantly heated by AC heating, with Alfven wave turbulence providing the primary contribution, accounting for at least 80% of the entire coronal heating in the present simulation. Our results strongly support the use of Alfven wave turbulence-based models employed in space weather and stellar activity research, such as the Alfven Wave Solar Model (AWSoM) and the Magnetohydrodynamic Algorithm outside a Sphere (MAS).

AlgoPlasma: Open Algorithms for Plasma Modeling

Yinjian Zhao, Zhongping Zhao, Zhe Liu, Baisheng Wang, Zilong Peng, Xin Luo, Lihuan Xie, Xi Chen, Zhijun Zhou, Kunpeng Zhong, et al.

AlgoPlasma is an open-source library in which core numerical algorithms for plasma modeling are implemented as modular, well-documented, and independently testable components. Rather than offering a complete simulation code, it allows researchers to select, adapt, and assemble the required components into application-specific workflows. The current release is centered on particle-based simulation, while AlgoPlasma is designed to encompass a broader range of approaches to plasma modeling. It provides components for particle initialization and advancement, particle--grid coupling, field solution, collision modeling, parallel data exchange, input/output, and selected fluid updates. Documentation links mathematical formulations to source implementations, interfaces, and usage, while verification and validation cases evaluate numerical accuracy and physical behavior. AlgoPlasma thus establishes a shared algorithmic foundation for plasma modeling, transforming repeatedly reimplemented numerical methods into open, reusable, tested, and explainable components for research, verification, education, and collaborative development.

Relative and absolute dosimetric commissioning of the ELIMAIA--ELIMED laser-driven proton beamline at 23.45 MeV

G. A. P. Cirrone, G. Petringa, A. Kurmanova, R. Catalano, A. Amato, A. Pappalardo, M. G. Guarrera, F. Schillaci, P. Blaha, F. Grepl, et al.

The development of laser-driven proton beamlines for biomedical and radiobiological applications requires traceable dosimetry and reliable online monitoring at the irradiation point. In this work, we report the relative and absolute dosimetric commissioning of the ELIMAIA--ELIMED laser-driven proton beamline at ELI Beamlines using an energy-selected proton beam with an average energy of about 24~MeV. Radiochromic-film measurements were used to characterize the transverse dose distribution, the depth--dose profile, and the proton energy spectrum at the irradiation point. The reconstructed spectrum was centred at $23.45~\mathrm{MeV}$ with a FWHM of $2.60~\mathrm{MeV}$, while the transverse dose distribution showed an approximately $5.5~\mathrm{mm}$ field size with a millimetre-scale homogeneous region. The Faraday Cup was used as the absolute reference detector for dose to water determination and for cross-calibrating the Dual-Gap Ionization Chamber, which was operated as the primary online dose monitor. The Integrating Current Transformer and Secondary Electron Monitor were evaluated as upstream relative beam monitors. An independent RCF--FC dose comparison performed over fifty consecutive shots yielded $30.45 \pm 3.5~\mathrm{cGy}$ from the FC and $36.46 \pm 1.8~\mathrm{cGy}$ from the EBT3 film. Dedicated G4ELIMED Monte Carlo simulations showed that the RCF positioned upstream of the FC perturbs the FC measurement through proton losses outside the FC acceptance and through an increase in the effective beam area at the FC entrance. Applying the resulting correction factor, $C_{\mathrm{MC}}=1.109$, reduced the residual difference between the FC- and RCF-derived doses to about $7\%$. These results establish the dosimetric chain of the ELIMAIA--ELIMED beamline under the present low-fluence commissioning conditions and identify the operational limits of the online monitoring system.

Enhancement of axial magnetic field generation during relativistic self-channeling of laser radiation propagating along thin films

Vitalia A. Kuleshova, Artem V. Korzhimanov

Within the framework of a stationary model of relativistic self-focusing, it is shown that the addition of a thin layer of denser plasma on the propagation axis of a circularly polarized beam increases the axial magnetic field generated by it via the inverse Faraday effect by more than a factor of 10 compared to the case of a homogeneous plasma. The magnitude of the axial field in this case can exceed the magnetic field of the laser wave, which makes it possible, at radiation intensities $\sim 10^{26}$~W/cm$^2$ potentially achievable in the near future, to approach the level of teragauss quasi-stationary fields.

Physics of Circular Polarized Ion-Scale Waves in Hybrid Simulations of Alfvénic Fluctuations

Hai Yang Harry Qian, Trevor A. Bowen, Carlos A. Gonzalez, Nikos Sioulas, Alfred Mallet, Kristopher G. Klein, Daniel Verscharen, Stuart D. Bale

Ion cyclotron waves (ICW) and fast magnetosonic/whistler waves (FMW) are fundamental electromagnetic modes at ion kinetic scales, yet their generation mechanisms and roles in plasma evolution remain poorly understood. We analyze a 2.5D hybrid simulation of broadband Alfvénic fluctuations, where the proton velocity distribution is modeled as a sum of two bi-Maxwellian components: a thermal core and a drifting beam. Using wavelet-based wave identification, bi-Maxwellian VDF fitting, and the PLUME linear dispersion solver, we find that ICW behave as linear modes. Growth is intermittent, occurring when core temperature anisotropy builds up, and is driven mainly by the core (the beam contributes negligibly). Poynting flux analysis shows that ICW are predominantly forward-propagating, with a net energy flux ratio of $+1$ across all frequencies, consistent with the initial condition. FMW present a stark contrast: PLUME solutions often yield very small (near-zero) linear growth/damping rates. The species decomposition breaks down when $|γ/ω_r| \gtrsim 0.368$, indicating that linear theory predicts these waves to be strongly damped and not describable by linear eigenmodes. Nevertheless, FMW are clearly observed in the wavelet helicity spectrogram, indicating that they are generated by nonlinear processes (e.g., parametric decay or phase steepening) and persist despite linear damping. The net energy flux ratio for FMW is close to $+1$ at low frequencies but decreases at higher frequencies, yet never reaches zero (net energy flow remains forward). These results demonstrate that ICW are linear, core-driven waves that transfer energy to the plasma, while FMW are heavily damped, nonlinearly generated waves.

Scalable production of nuclear battery alpha emitters using fusion neutrons

J. F. Parisi

Nuclear batteries powered by alpha decay have been deployed successfully for over 60 years, on a worldwide $^{238}$Pu supply of kilograms per year. We show that the 14 MeV neutrons of a single deuterium-tritium fusion plant can produce alpha emitter battery fuels up to tons per year, in three classes: fuels with completely new production pathways ($^{236}$Pu, $^{227}$Ac, $^{210}$Pb), fuels previously proposed whose scarce feedstock the same pathways now breed at scale ($^{232}$U, $^{228}$Th), and the established $^{238}$Pu. OpenMC simulations of actinide channels in a tokamak blanket give, per GW yr of fusion: 11 to 57 kg of $^{236}$Pu, whose chain releases 18 GJ per gram over a century, ending at stable $^{208}$Pb, plus up to 5.2 t of co-product $^{238}$Pu; up to 1.4 t of $^{231}$Pa from thorium, and, from channel with $^{231}$Pa feedstock, up to $\sim$15 t of $^{232}$U or $\sim$122 kg of $^{210}$Pb, with $^{227}$Ac produced at 21 g/yr per tonne of $^{231}$Pa. Neutron capture also upgrades $^{241}$Am to a $^{242}$Cm/$^{242m}$Am/$^{241}$Am/$^{238}$Pu blend with up to 10 times higher power density. The same $^{236}$Pu and $^{232}$U also serve as proliferation safeguards: the $^{237}$Np, $^{232}$Th, and $^{231}$Pa channel products are self-protecting, the plutonium by $^{236}$Pu and $^{238}$Pu decay heat and the 2.6 MeV gammas from $^{208}$Tl content, and similarly the uranium from its $^{232}$U. Many of these fuels ($^{236}$Pu, $^{232}$U, $^{228}$Th, $^{227}$Ac) have an order of magnitude higher power and energy density than current alpha emitters, and at human spaceflight-relevant doses the $^{227}$Ac and $^{210}$Pb chains need less shield mass than $^{238}$Pu or $^{241}$Am above a few hundred watts. Fusion neutrons could therefore enable nuclear batteries at the kilowatt to megawatt scale and unlock new possibilities for power sources requiring exceptionally high energy density.

Time-Domain Benchmark Solutions for Bernstein Waves

M. Pelkner, K. Hallatschek

In previous work, we introduced a semi-analytical method for computing time-domain solutions of linearized Vlasov problems. Rather than representing the plasma response as a sum of residues associated with Landau poles, the method constructs a regularized frequency-domain response spectrum that is subsequently inverted numerically. Explicit applications have so far been limited to unmagnetized plasmas. In this work, we extend the construction to plasmas with a uniform background magnetic field and a Maxwellian equilibrium distribution. We apply the resulting formulation to the electrostatic ion-Bernstein density response of a plasma with kinetic ions and adiabatic electrons, and provide error estimates for truncating both the spectral integration domain and the cyclotron-harmonic expansion. The solution serves as a time-domain reference for damped finite-$k_\parallel$ regimes, in which dispersion-relation roots alone are insufficient for pointwise-in-time verification of simulation codes. Finally, we indicate how the framework can be extended to fully electromagnetic problems.

Flow-thermochemistry coupling governs pressure-dependent CO$_2$ conversion in vortex-stabilized microwave plasma reactors: Insights from three-dimensional CFD modeling

Qinghao Shen, Cas van Deursen, Pieter Willem Groen, Lex Kuijpers, Mauritius C. M. van de Sanden

In this work, a three-dimensional computational fluid dynamics model is developed for a vortex-stabilized microwave CO$_2$ plasma reactor operating over the pressure range of 100-400 mbar. The model combines experimentally constrained, emission-based plasma sizes and volumetric heat-source distributions with thermally dominated finite-rate heavy-particle chemistry for a multi-component mixture. Turbulent flow and transport are described using the SST k-omega model. The model reproduces the measured radial gas-temperature profiles in the plasma core and the non-monotonic pressure dependence of CO$_2$ conversion, including a maximum at 150 mbar and a pronounced decrease at 400 mbar. A vortex-driven recirculation region redistributes gas upstream. Turbulent mixing and cooling are strongest near the upper reactor boundary, but their contribution decreases as pressure increases. The pressure dependence of conversion is determined by the competition between CO$_2$ dissociation and CO recombination. CO$_2$ direct dissociation reaction dominates in the high-temperature plasma core, whereas O-assisted conversion reaction contributes near the plasma edges and in the surrounding hot region. At 150 mbar, enhanced CO$_2$ dissociation is accompanied by limited CO loss, resulting in the highest conversion. With pressure increasing to 400 mbar, slower cooling and more frequent three-body collisions promote CO recombination in the afterglow, causing more than 60% of the CO formed near the plasma to be lost downstream. Moreover, additional CO loss occurs in the upper region of the reactor at higher pressures because of the reduced cooling rate.

The Production of Electron-Capture Elements in Thermonuclear Supernovae: Theory vs. Observations

S. Shiber, P. Hoeflich, T. Mera, E. Fereidouni, Z. Levy, D. Maci, C. Ashall, K. Medler, J. M. DerKacy, E. Baron, et al.

Type Ia supernovae (SNe Ia) explosively destroy carbon-oxygen white dwarfs (WDs) in multiple stellar systems. They produce approximately 50% of the iron-group elements in the Universe, synthesize electron-capture (EC) elements, drive nuclear physics experiments, and underpin high-precision cosmology. To first order, the outcome is governed by nuclear physics, a property often described as stellar amnesia. Recently, this stellar amnesia has begun to be broken by the nearly universal detection of EC elements with JWST. These elements trace high-density burning, largely ruling out the currently popular helium-triggered, sub-Mch detonation models as the dominant channel. Instead, the ubiquitous presence of EC is shifting back the focus to dynamical and secular mergers, and near-Mch explosions similar to the deflagration model W7, but in which the nuclear flame undergoes a deflagration-to-detonation transition. The early deflagration phase is especially important because spherical simulations identify the central WD density, and thus the WD mass, as a key parameter governing the explosion. Here, we present detailed magneto-hydrodynamical simulations. We find that small-scale, pre-existing turbulence expected from the pre-explosion smoldering phase is essential for overcoming the fundamental challenges imposed by the intrinsic 3D physics. This turbulence systematically reduces the production of EC elements by about a factor of two, implying the need for WD central densities closer to those associated with accretion-induced collapse to a neutron star. We also demonstrate the effect of magnetic fields near the saturation field strength and highlight the need for higher-precision EC rates at low Ye.

Macroscopic Stability of a Rapidly Rotating Theta Pinch

Richard Fitzpatrick

The macroscopic ideal-MHD stability of an axisymmetric mirror device with sonic levels of plasma rotation is analyzed by approximating the plasma equilibrium as a rotating theta pinch possessing an artificial gravity. An eigenmode equation is derived that governs the stability of the equilibrium to small perturbations in the case of an arbitrary plasma angular velocity profile. The stability of the m=1 and m=2 modes is investigated. The plasma is found to be stable to these two modes provided that it is sufficiently short in the axial direction. The critical axial length of the device below which the modes are stabilized first decreases with increasing plasma rotation, attains a minimum value when the rotation is roughly sonic, and then increases with increasing plasma rotation. The value of the plasma rotation off the magnetic axis is found to have a significantly stronger effect on the stability of the modes than the value on the magnetic axis.

Large-Scale Dynamos Driven by Shear-Flow-Induced Jets

B. Tripathi, A. E. Fraser, P. W. Terry, E. G. Zweibel, M. J. Pueschel, R. Fan

University of Wisconsin-Madison, Columbia University, University of Colorado, Ruhr-Universität Bochum, Dutch Institute for Fundamental Energy Research

At every scale they occupy, magnetic fields affect various phenomena, including star formation, cosmic ray transport, charged particle acceleration, space weather, transport in planetary atmospheres, and laboratory plasmas. These fields are often generated and sustained by turbulent flows in a process called the dynamo. In 1955, E. N. Parker parameterized the effects of small-scale turbulence to propose a mean-field dynamo theory. The widely used theory reproduces observed large-scale fields but suffers from difficulty in tuning parameters as they are not justified from first principles: Studies of turbulent flows show tangled magnetic fields, which are folded and fragmented into small-scale structures due to shear-flow straining. Here, considering a shear flow that is unstable and driven, we develop analytic theory and perform three-dimensional (3D), advanced computer simulations of turbulence with up to 4096 x 4096 x 8192 grid points, showing ab initio generation of quasi-periodic, large-scale magnetic fields. The generation occurs via the mean-vorticity effect---an additional mean-field dynamo process postulated in 1990. Crucial to this dynamo is the prior generation of large-scale 3D jets, robustly produced as topologically protected and exact nonlinear solutions of the magnetohydrodynamic equations. The jet-driven dynamo applies to shear-driven laboratory and astrophysical systems. These include binary neutron star mergers, where the reported dynamo likely operates on microsecond timescales to produce in milliseconds some of the strongest magnetic fields in the Universe, providing signals for multimessenger astronomy.

Monoenergetic acceleration of charge-neutralized ion bunches to GeV-scale energies by the combination of a high-current electron beam and an ionization front

Jiyuan Chen, Jihoon Kim, Roopendra Singh Rajawat, Gennady Shvets

Compact heavy ion accelerators have numerous applications, ranging from heavy ion fusion to carbon ion radiotherapy, and testing radiation-hardened electronics. The demand could be met by developing high-gradient traveling wave plasma accelerators of high-charge ($\simμ\mathrm{C}$) relativistic ion beams. We will discuss a novel ion acceleration regime -- Counter-propagating ionization Front Acceleration (CFA) -- utilizing counter-propagating Ionization Front (IF) and high-current Relativistic Electron Beam (REB). Theoretical modeling and 3D PIC simulations demonstrate the possibility of using typical REBs produced by induction voltage adders propagating through a gas-filled tube undergoing laser ionization to achieve acceleration gradients in excess of $\sim 250 {\rm MeV/m}$ while accelerating micro-Coulombs of ions over meters distance. A unique energy conversion mechanism -- from the REB to electromagnetic fields to the ions is discussed, as well as the limits on the accelerated ions charge and the degree of its neutralization, acceleration gradient, and ion energy spread.