Recent Publications

Particle-resolved pathways to energetic-ion formation in a fluctuating low-current hollow-cathode plume

Baisheng Wang, Zilong Peng, Yinjian Zhao, Zhongxi Ning

Energetic-ion formation in a low-current hollow-cathode plume is investigated using experiments, self-consistent electrostatic particle-in-cell (PIC) simulation, and particle-resolved analysis. Retarding potential analyzer measurements show a substantial energetic-ion population over discharge currents of 0.8-3.5 A, while probe measurements reveal broadband plume fluctuations. Two-point phase-derived frequency-wavenumber measurements do not resolve a continuous ion-acoustic dispersion branch within the principal apparent-wavenumber interval. Because the inferred wavenumber is obtained from a cross-spectral phase defined modulo 2pi, the fluctuation diagnostics do not provide an unambiguous modal attribution for the energetic-ion population. A representative PIC plume, used as a qualitative kinetic reference, likewise develops broadband time-dependent electrostatic fluctuations together with a nonthermal energetic-ion population. Particle-resolved analysis shows that the energetic outflow is dominated by ions generated through ionization inside the plume, while source localization biases access to distinct trajectory and escape families. Matched field controls further show that time-averaged and frozen fields strongly suppress access to high-energy trajectories relative to the full time-dependent field over the analyzed interval. At the single-particle level, ion kinetic-energy gain is determined by electrostatic-field work accumulated along the actual trajectory, with different escape families exhibiting distinct radial and axial work contributions. These results establish a source-trajectory-field-work pathway for energetic-ion formation that can be identified without first assigning the fluctuating plume to a unique resolved plasma mode.

Modeling of plasma transport during edge-localized mode in tokamak using a kinetic Vlasov-Poisson code

Ce Wang, Sven Van Loo, Geert Verdoolaege

A kinetic parallel transport code KOBRA based on a finite-volume method is developed to study edge localized mode (ELM) plasma transport from the mid-plane to divertor targets. The large scale separation between the Debye length (~cm) and the connection length (~10-20 m) leads to prohibitive computational cost in full 6D simulations. To alleviate this, an adaptive-mesh refinement (AMR) strategy is employed. Comparisons with uniform-grid simulations show that AMR accurately reproduces the characteristic ELM dynamics, including the rapid rise and slow decay of divertor fluxes, as well as the early-time peak induced by fast electrons. Analysis of the electron distribution and self-consistent electric field reveals that AMR efficiency is closely linked to phase-space evolution. Overall, AMR achieves comparable physical accuracy while reducing memory usage by 30-40% and accelerating computations by up to a factor of two, demonstrating its effectiveness for high-dimensional kinetic ELM simulations.

kobra: a new Vlasov code intended for plasma-wall modeling

Sebastian Konewko, Nathan Maestracci, Sven Van Loo

In a fusion device plasma-wall interactions \edit{on the sheath scale} can be modeled as a collisionless problem. When modeling these regions particle-in-cell codes suffer from statistical error originating from undersampling the velocity space. On the other hand, Vlasov codes do not have this issue as they evolve the full distribution function. Here, we present a new finite-volume Vlasov code, kobra, equipped with adaptive-mesh refinement to reduce computational effort. Currently, the code solves the Vlasov-Poisson equations. We validate our code in 1d1v and 1d2v using established benchmarks, i.e. the two-stream instability, Landau damping, the Dory-Guest-Harris instability, and also a classical electrostatic plasma sheath. We find that the code reproduces the theoretical properties of these problems well. More importantly, the adaptive grid provides a computational gain that is likely to scale to higher dimensional, plasma-wall simulations.

Physics-Informed Neural Networks to Infer the Perpendicular Energy Conductivity in the Scrape-Off Layer of Stellarator Devices

J. Gallego, P. Protopapas, A. Bustos, A. Alonso, S. Barquero, A. Baciero, I. Rivera, J. A. Moríñigo, R. Mayo-García

In this work, we develop an inverse Physics-Informed Neural Network (PINN) framework to infer the dependence of the scrape-off layer (SOL) perpendicular heat conductivity on plasma density and temperature, $κ_\perp(n,T)$. The method combines radial profile measurements of electron density and temperature with the residual of a reduced one-dimensional SOL transport equation, so that the inferred conductivity is constrained by both the measurements and the underlying transport model. Three neural networks are trained simultaneously: two reconstruct the temperature and density profiles as functions of the radial coordinate and transported power, while a third represents the effective conductivity as a function of the local density and temperature. The framework is first validated using synthetic data generated from a prescribed conductivity function, allowing the inferred $κ_\perp(n,T)$ to be compared directly with the ground truth. The model recovers the imposed functional dependence with errors below $10~\%$ in the data-constrained region. Bootstrap resampling is shown to provide a practical indicator of prediction reliability and consistency. A scan in the number of plasma profiles used for training and the number of radial measurement positions per profile identifies a practical trade-off between reconstruction accuracy and data availability. Finally, the method is applied to an experimental dataset from the TJ-II stellarator obtained with the helium-beam diagnostic. This exploratory application provides an initial estimate of the effective SOL conductivity and illustrates the potential of inverse PINNs for extracting transport information from plasma edge measurements.

Physics-Informed Drift Diagnosis for Laser-Plasma Accelerator Operations

Ou Labun, Calin Hojbota, Mara Klebonas, Mike Downer, Rafal Zgadzaj, Phil Franke, Lance Labun

Laser-plasma accelerators (LPAs) sustain accelerating gradients of order $100\,\mathrm{GV/m}$, but routine operation remains difficult: electron beam metrics drift over an operating shift, and the root physical cause is often invisible to the available diagnostics. We formulate LPA operation as a latent state-space model in which three effective interaction-point variables, the normalized laser amplitude $a_0$, the normalized plasma electron density $\tilde n_e$ and the residual pulse chirp $\mathcal{C}$, are inferred from routine electron beam observations by an extended Kalman filter. The emission model, which maps the latent state to the diagnostics, is kept structurally separate from the {transition} model, which describes how the latent state evolves between shots. The separation supports diagnosis in two stages, one asking which latent variable moved and one asking what moved it. The implemented emission model is a toy model, yielding an expected performance in line with current facilities and using 3D blow-out regime dependencies where relevant. We conduct synthetic sessions to test the effectiveness of the detection and attribution protocols, finding that attribution is limited by excitation rather than by shot count or diagnostic resolution. Because the construction needs only a set of physical latent variables, an emission model and a family of hardware-derived transition models, it transfers to other drift-prone subsystems. We argue that the accuracy of the whole procedure is limited by the emission model rather than by the inference method.

Effect of pedestal current on the density window for ELM suppression using n = 4 RMP in EAST

Xuemin Wu, Youwen Sun, Qun Ma, Shuai Gu, Manni Jia, Yueqiang Liu, Yifeng Wang, Cheng Ye, Pengcheng Xie, Alberto Loarte, et al.

Existence of operational window in both edge safety factor and line averaged plasma density for suppression of ELMs using n=4 Resonant Magnetic Perturbations in low input torque plasmas has been observed in EAST experiment, in which q95 and plasma normalized beta (βN) close to that required in ITER high-Q operation. Here, n is toroidal mode number of the magnetic perturbation. In contrast to previous reports from other tokamaks, there is not only an upper density limit but also a lower one for accessing ELM suppression. Modelling results using the MARS-F code show that the RMP with linear plasma response has a peak at an intermediate density and decays as the density increases or decreases, which results in a minimal RMP field penetration threshold at the intermediate density. In this experiment, the observed lower density limit operationally manifests a sensitivity of the q-profile: different densities alter the edge current profile, which change the alignment of the eigenmode structure with the RMP coil configuration, causing a reduction of the resonant field in both low- and high-density cases, and hence making field penetration more difficult. The modelled window of the strongest resonant plasma response in terms of [⟨ne⟩, q95] agrees well with the observed ELM suppression in EAST. Peeling-ballooning modes stability analysis using the ELITE code shows that plasmas gradually approach peeling instability boundary caused by increase of edge bootstrap current as the plasma density decreases, which is consistent with the observation that ELMs come back again in lower density plasmas for fixed q95. These results indicate that linear modelling with full toroidal geometry can well predict the optimized RMP configuration for ELM suppression and reveal the important roles of pedestal plasma current, which need to be carefully considered in the application of high n RMPs for ELM suppression in future ITER.

BLUEMIRA: a modular, open-source framework for designing tokamak fusion reactors

Matti Coleman, James Cook, Fabrizio Franza, Ivan Maione, Simon Mcintosh, Hudson Baker, Alexander Ian Blair, Shail Desai, Oliver Funk, Georgina Graham, et al.

As efforts across the world to deliver fusion power enter their respective conceptual design stages, many begin to encounter the difficulties of generating complete designs, evaluating them, performing trade-off studies, and comparing different alternatives. Many organisations successfully use fusion reactor systems codes with low-fidelity formulations of the critical physics and engineering aspects to inform the initial design stage and steer them towards a region of the design space they consider attractive, but struggle to explore the design space in more depth. This work is an overview of the open-source tokamak reactor design framework, BLUEMIRA, built from the BLUEPRINT [1, 2] and MIRA [3, 4] codes, which goes beyond the “classical” systems code paradigm and enables the user to parameterise their reactor designs at higher fidelity. Reactor design workflows ranging from 0-D “radial build” design, equilibrium design, to 3-D CAD and automated neutronics analyses can be run within minutes. We present the functionality implemented in the BLUEMIRA code and demonstrate its application to the conceptual design of a conventional aspect ratio fusion reactor, performing an indicative design study over aspect ratio and number of toroidal field coils.

A polar-harmonic unified gas-kinetic scheme for magnetized ion dynamics from cyclotron kinetics to the Hall-Pedersen constitutive limit

Yixiao Wang, Zhigang Pu, Xing Ji, Kun Xu

Magnetized ion transport in weakly ionized plasmas ranges from gyroangle-dependent kinetics to Hall-Pedersen drift-diffusion as collisionality and magnetization vary. We develop a polar-harmonic unified gas-kinetic scheme (PH-UGKS) for the ion Vlasov-BGK equation in a uniform magnetic field. The scheme evolves the full ion distribution by coupling a conservative density update to exponential evolution of its nonequilibrium component. Exact collision-rotation integration in gyroangle Fourier space is combined with a time-averaged kinetic flux that incorporates spatial transport and electric acceleration, together with a compact Hall-Pedersen correction to the density flux. The scheme conserves ion number, and analysis establishes second-order temporal consistency and asymptotic preservation of the Hall-Pedersen density limit at fixed magnetization. Numerical tests reproduce ion Bernstein dispersion and Dory-Guest-Harris growth rates and resolve changes in the gyroharmonic spectrum as the collision-to-gyrofrequency ratio varies. The driven ion-flux response agrees with an independent characteristic-Volterra reference, including finite-frequency departures from the instantaneous Hall-Pedersen relation. In collisional tests, accurate responses are obtained with time steps far larger than both the collision time and the gyroperiod. Fixed-resolution density tests confirm convergence to the corresponding Hall-Pedersen discretization. The same kinetic formulation thus connects kinetic response and macroscopic transport without switching to a fluid solver or subcycling microscopic time scales.

Connecting Dynamo Theory with DNS Data: A Computational Analysis of α and \b{eta} Effects

Kiwan Park

We investigate the influence of current helicity on the turbulent magnetic diffusivity $β$ using three complementary derivations of the $α$ and $β$ coefficients, based on the large-scale magnetic field $\overline{\mathbf{B}}$, the turbulent velocity $\mathbf{u}$, and the turbulent magnetic field $\mathbf{b}$. Applying these coefficients to raw DNS data, we reconstruct $\overline{\mathbf{B}}$ and compare the results with the original simulations. In the kinematic regime all models agree well with the DNS data. In the nonlinear regime, however, $β_{\mathrm{vv-vw}}$ alone produces unbounded growth of $\overline{\mathbf{B}}$. Including the contribution from turbulent magnetic fields ($β_{\mathrm{bb+jb}}$) suppresses this unphysical growth and restores agreement with the DNS results. We find that kinetic helicity drives $β$ more negative, while current helicity shifts it back toward zero. Weighted combinations of the coefficients further show that the $β$ effect dominates the evolution of $\overline{\mathbf{B}}$ throughout, whereas the $α$ effect becomes important mainly for sustaining the field in the nonlinear regime. The corresponding IDL analysis scripts are provided to facilitate practical implementation of the theoretical models.

High-charge, highly polarized positron beams generated from a laser-driven nanowire-array target

De-Sheng Zhang, Cui-Wen Zhang, Kun Xue, Feng Wan, Xue-Ren Hong, Jian-Xing Li, Bai-Song Xie

The generation of high-charge, highly polarized positron beams in the interaction of a linearly polarized laser pulse with a nanowire-array target is investigated. Here, laser-driven electrons emit high-energy photons through nonlinear Compton scattering (NCS), which subsequently produce electron--positron pairs through the nonlinear Breit--Wheeler (NBW) process. We model this interaction using two-dimensional spin-resolved quantum electrodynamics particle-in-cell (QED-PIC)} simulations. At positron birth, the sign of $S_z$ is statistically correlated with that of the local $B_z$. The spatiotemporal field structure arising from the laser--nanowire interaction strengthens the correlation between the birth spin sign and the direction of the subsequent transverse Lorentz impulse, thereby limiting cancellation between opposite-spin contributions at a given angle. The results show that the average polarization degree reaches $|\bar S_z|\approx0.46$, and the positron charge satisfying $|\bar S_z|>0.3$ is approximately $308\,\mathrm{nC}$. Parameter scans reveal that the high-polarization positron charge is maximized at intermediate target densities and nanowire periods. Such a source could enable polarization-sensitive studies of strong-field QED and spin-dependent phenomena in high-energy and materials physics.

JET neutron emissivity reconstruction using the Minimum Fisher Information method for 1 ms temporal resolution

Katarzyna Mikszuta-Michalik, Daniele Marocco, Basilio Esposito, Marco Riva, Gianluca Pucella

A tomography code based on the Minimum Fisher Regularisation method with a geometry description based on the field of view has been applied to reconstruct the neutron emissivity profiles from measurements performed in 2021 during the second JET deuterium-tritium campaign (DTE2). The Neutron Profile Monitor Upgrade installed at JET offers a unique opportunity for investigating the neutron emissivity in a full poloidal cross-section with unprecedented time resolution (~1 ms), thanks to the availability of digital data reprocessing and the high neutron emission that enables milisecond line-integrated raw data with low statistical error. Measurements of 14.1 MeV neutrons provided by the Bicron BC418 plastic scintillators installed in the JET neutron camera have been used.The analysis focuses on the effect of sawtooth crashes on the neutron emission profiles. Tomography results were compared with 1D reconstructions using two different methods and 10 ms time resolution. The use of the neutron camera data sampled at 1 ms enables the observation of hollowing of the neutron emissivity profiles during the sawtooth crash and their subsequent flattening.

Machine learning tomography for sparse-view soft X-ray diagnostics in MAST-U

Brian Steward, Marco Cecconello

Understanding magnetohydrodynamic (MHD) instabilities is crucial for advancing magnetic confinement fusion. Soft X-Ray (SXR) tomography enables reconstruction of the spatial structure of these instabilities from one dimensional line integrated measurements, however traditional tomography methods struggle with sparse diagnostic arrays. The MAST-U tokamak has a sparse camera array of 28 non-intersecting lines of sight (LoS), making traditional tomography techniques challenging. We present a machine learning (ML) tomographic method using a 50 layer residual neural network (ResNet) trained on synthetic data simulated by the transport code TRANSP. ML tomography is compared with traditional minimum Fisher information tomography, and with second derivative smoothing tomography. The ML technique shows a 78 to 94 times median improvement in the mean square error (MSE) over traditional 28 LoS reconstructions of synthetic data, at a computational speed which is 144 to 294 times greater (25 ms on GPU hardware). Notably, ML tomography with 28 LoS outperforms traditional methods even when provided with three times as much data in the form of 84 intersecting synthetic LoS. Validation on experimental data from MAST-U demonstrates accurate reconstructions of peaked and broad plasma SXR profiles, and the evolution of sawtooth instabilities are visualized with ML tomography. The sub-30 ms reconstruction time enables real-time plasma diagnostics, opening new possibilities for MHD monitoring in fusion devices.

Kinetic modeling of lithium impurity transport in the HIDRA stellarator during controlled lithium evaporation

Steven Gula, Nina Mihajlov, Giovanni Diaz, Arnav Goyal, R Maingi, D Andruczyk, Davide Curreli

Post-operational imaging of the Hybrid Illinois Device for Research and Applications (HIDRA) stellarator following controlled lithium evaporation reveals distinct streak patterns of lithium deposition along the vacuum vessel wall. These localized structures exhibit strong alignment with magnetic field topology, suggesting that impurity transport in the scrape-off layer (SOL) is governed primarily by the three-dimensional magnetic geometry and associated electrostatic fields. In this work, large-scale kinetic simulations of trace lithium impurity transport are used to investigate the formation of these deposition structures. The magnetic field topology is reconstructed using a Biot-Savart law approach, from which a normalized flux-surface parameter is computed and used to prescribe the background helium plasma profiles. Electron density and plasma potential are modeled as normalized flux-surface-dependent fields scaled by representative scalar measurements and analytical sheath-based estimates, with the electric field obtained from ∇V p . These reconstructed fields are used to compute lithium ion trajectories and generate spatial maps of deposition on the vessel wall. The simulations reproduce streak-like deposition features aligned with magnetic flux tubes intersecting the wall, and the simulated deposition locations and orientations show qualitative agreement with experimental observations. These results demonstrate that magnetic field structure plays a dominant role in shaping impurity transport and deposition in HIDRA and establish a predictive framework for modeling impurity transport and surface deposition in plasma-surface interaction experiments involving liquid-metal plasma-facing components.

Positive sheath generation in front of fishscaled divertor plates

Choong-Seock Seock Chang, Seung-Hoe Ku, Xin Zhang, Trenton Brewer, Nicolas Lopez, Chris Marsden

Fish-scale divertor plates are becoming popular in the design of next-generation magnetic fusion reactors to make the edge of one divertor tile sits safely behind the shadowed profile of the preceding tile, hence to protect the leading edge from extreme heat deposition. Recently, ST40 tokamak has observed an extremely narrow and peaked divertor heat-load footprint near the separatrix strike point on the fish-scale plates, on top of the usual ion-drift-width scale footprint. This observation raised concern over the severely localized burn even when the well-known ion-drift-width scale burn issue can be resolved [X. Zhang et al., Nucl. Mater. Energy 41, 101772 (2024)]. In this report we demonstrate that the lost ion gyro-orbits to the tilted tile edges can lead to a significant positive-potential sheath at a practical magnetic field incidence angle that is much greater than the usual incidence angle for the ion gyro-sheath formation on flat strike-surface. This could make the kinetic electron heat-flow spilled over from the confined region down the divertor legs to gain a significant kinetic energy amplification and allow a highly localized heat-load peak under a low edge turbulence condition. Possible intrinsic mitigation mechanism is discussed. This issue should be considered carefully in the design of next generation magnetic fusion reactors.

RF-Specific Tungsten Erosion and Global Transport in ITER under Neon Seeding

Atul Kumar, Dhyanjyoti Nath, Wouter Tierens, Jeremy D. Lore, Andrei Pshenov, Tom Wauters, Andrea Galvan, Davide Curreli, Syun'ichi Shiraiwa, Nicola Bertelli, et al.

Ion cyclotron radio-frequency heating (ICRH) is a key auxiliary heating system in ITER, but high-power RF operation can enhance plasma-material interactions through rectified RF sheath potentials on antenna structures and nearby plasma-facing components. We present the first predictive application of the STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to assess RF sheath-driven tungsten (W) erosion and global impurity transport from the ITER ICRH antenna under ITER-relevant neon-seeded conditions. STRIPE couples SOLPS-ITER plasma backgrounds, full-wave RF sheath calculations, geometry-specific ion energy-angle distributions, sputtering physics, and three-dimensional impurity transport. Simulations predict RF sheath potentials of 1 to 3 kV on antenna limiter sidewalls, increasing gross W erosion by about a factor of 64 relative to thermal sheath conditions and producing a gross source of 3.34e18 W atoms per second. Erosion is governed by RF-modified ion energy-angle distributions together with local plasma flux rather than sheath voltage alone. About 10 percent of sputtered W is locally redeposited, giving a net source of 3.01e18 W atoms per second. The RF-induced antenna source remains about three orders of magnitude smaller than the thermal divertor source and more than two orders of magnitude smaller than the integrated thermal main-chamber source. After 100 ms, about 22 percent of the mobile W inventory resides within the SOLPS-covered confined-plasma region, corresponding to an annular W concentration of 1.70e-6. These results indicate that the ITER ICRH antenna is unlikely to dominate the total W source budget under the conditions considered and demonstrate the need for coupled modeling of RF waves, sheaths, sputtering, redeposition, and global impurity transport.

Generalized Similarity Theory for Plasmas

Yangyang Fu

A generalized theory of plasma similarity is established based on the scaling of the Boltzmann equation coupled with the full set of Maxwell's equations. The predicted similarity scalings are demonstrated through first-principles particle-in-cell and fluid simulations across diverse operational regimes, including collisionless plasmas with two-stream instabilities, pure electron diodes spanning the classical to relativistic regimes, electromagnetic-wave-driven plasmas, and discharge plasmas ranging from low to high ionization regimes. The generalized theory reveals the intrinsic scale-invariant nature of plasmas under specified conditions, rooted in the symmetry and scaling transformations of the governing equations.

On the modeling of oblique firehose instabilities in regularized Kappa plasmas using ALPS

D. L. Schröder, M. Lazar, H. Fichtner, K. G. Klein, D. Verscharen

In-situ measurements in space plasmas indicate that the velocity distributions of charged particles are not in thermal equilibrium, deviating from a standard Maxwellian mainly due to anisotropies and suprathermal populations which enhance high-energy tails. Although the Standard \k{appa}-Distribution (SKD) is well established in modeling these non-equilibrium distributions, its application is often viewed controversially due to certain unphysical implications, in particular divergent velocity moments. To address these issues, the Regularized \k{appa}-Distribution (RKD) was introduced. Such advanced, in general anisotropic RKDs are invoked here for the first time to investigate oblique firehose instabilities, including those induced by the temperature anisotropy of electrons and protons, the dominant species in space plasmas. In weakly collisional plasmas, both of these instabilities are expected to play significant roles in the self-regulation of the macroscopic properties of space plasmas (e.g. the expanding solar wind) reported by observations. Is not yet possible to resort to a general dispersion tensor related to RKD plasmas (whose derivation is still a challenge), instead the Arbitrary Linear Plasma Solver (ALPS) is exploited here. The unstable solutions obtained for the already established Maxwellian & SKDs successfully validate the capability of ALPS. For RKDs, the instabilities confirm the stimulating effect of suprathermal populations, with lower \k{appa} values generally enhancing the firehose growth rates. For the oblique electron firehose instability, RKDs substantially modify the competition between periodic & aperiodic branches only at intermediate angles & can sustain significantly increased growth rates of especially aperiodic modes at highly oblique propagation beyond both the Maxwellian unstable regime & the one that is accessible with SKDs.

An Interactive Three-Dimensional Volume Visualization Pipeline for Fusion Simulation Data in Unreal Engine and Omniverse

Darian Figuera-Michal, Eui-Hyun Jung, Mintae Na, Sun-Jeong Kim, Jong-Hyun Kim, Jung Lee, Yu-Seop Kim, Jungpyo Lee, Jaemin Seo, Il-Hwan Kim, et al.

Study on X-ray Energy Deposition in Slab and Droplet Cases Using FLASH

Abinash Manikandan, Sophie L. Sharpe, Kirk A. Flippo, Lars Pedersen, Eric Cervi, Gustav Blankenberg, S. Woodruff

System-Level Tritium Transport Modeling for the Outer Fuel Cycle of the SCYLLA Breeding Blanket and Relevance for UNITY-2

Luigi Candido, Luke Taylor-King, Colin Baus, Christian Day, Italo Godoy-Morison, John McGrady, Minoru Jimma, Kyosuke Namba, Jonas Schwenzer, Tim Teichmann, et al.