2026-09-25 | | Total: 16
Sustained ICRF operation in a fusion power plant may require low edge densities to mitigate plasma-wall interactions, a regime which was recently achieved in WEST with very little impurity sputtering. Cold plasma theory, however, predicts singular radiofrequency electric fields in this regime, both at the lower hybrid resonance and along the resonance cones, raising the question of whether standard collisional cold plasma models suffice to describe low-density edge ICRF at all. Collisions in principle remove these singularities, replacing them with finite but sharply peaked fields. We derive these peak length scales analytically and confirm them with a 2D finite-element simulation using exponential mesh refinement, achieving micrometer resolution where needed. We conclude that edge collisions in cold plasma do not remove the need to resolve length scales ordinarily associated with hot-plasma and Bernstein-wave physics.
Slow magnetoacoustic pulsations (SMAPs) are found in MHD simulations of a tokamak plasma whose safety factor $q$ near the center is flat and slightly above unity ($q \gtrsim 1$). SMAPs are located on the central plateau of the slow magnetoacoustic continuum $ω_{\rm S} = k_\parallel c_{\rm S}$, where $c_{\rm S}$ is the speed of sound and $k_\parallel$ the wavenumber parallel to the magnetic field. In our model, SMAPs exist when the ion viscosity and thermal diffusivity are sufficiently low. They can be driven unstable by a pressure gradient in the $q \sim 1$ region when the electric resistivity is sufficiently high. Exponentially growing SMAPs consist of standing slow waves on magnetic surfaces that are radially synchronized into a quasi-interchange structure with poloidal/toroidal mode numbers $m/n=1/1$. After free energy depletion, saturated quasi-linear pulsations (alternating $n=1$ and $0$) in the central $q\sim 1$ region couple to distant $q\geq 2$ rational surfaces that undergo "reversible" magnetic reconnection: as the magnetic islands wax and wane with period $2π/ω_{\rm S}$, their X- and O-points alternate. These results show how the MHD model facilitates non-local coupling of slow waves, pressure-driven resistive interchange and tearing. This motivates further study in kinetic models, where collisionless mechanisms for fast reversible reconnection exist and proper treatment of parallel dynamics will allow to assess the role of Landau damping as well as the question whether the ${\mathbf B}$ field's weak ergodicity in the $q\sim 1$ region allows the waves to outpace the ion's parallel streaming to maintain the thermal misbalance underlying SMAPs. Also of interest are pulsations closer to the Alfvénic branch, satisfying $ω\approx k_\parallel v_{\rm A}$ with Alfvén speed $v_{\rm A}$, which require no resistivity and are less dependent on thermal misbalance.
Reliable correction of nonaxisymmetric error fields (EFs) is essential to the safety and performance of tokamak operation. Although resonant error field correction (EFC) is well established, supported by an improved understanding of 3D plasma response, the residual fields left after resonant EFC remain an open question: they are expected to be predominantly nonresonant, yet can still degrade both confinement and stability. Here we introduce a systematic EFC scheme that minimizes resonant and nonresonant EF effects simultaneously, based on neoclassical torque response computed from self-consistent perturbed equilibria. The scheme extends the method developed for designing quasi-symmetric magnetic perturbations to the case where actual error fields are present. Application to standard ITER target plasmas with a range of intrinsic EF scenarios demonstrates the advantages of this quasi-symmetric (QS) EFC scheme for controlling residual EFs. QS EFC consistently yields low-torque solutions while strongly suppressing resonant response, outperforming single-mode resonant overlap EFC in most cases and approaching the performance of multimodal resonant EFC. We also show that the QS EFC solution varies only tolerably between half- and full-$I_p$ ITER scenarios, despite the greater sensitivity expected from its higher-order nature.
Direct numerical simulation of magnetic reconnection is limited by the scale separation of resistive magnetohydrodynamics. At large Lundquist numbers $S$ the current layer thins as $S^{-1/2}$, forcing fine grids and short time steps that make parameter scans prohibitively expensive. Deep learning neural operators offer an alternative by learning the map between function spaces rather than individual solutions, so that a single trained model returns the state for any parameter and time at inference cost. We have developed a Fourier Neural Operator (FNO) based Physics-Informed Neural Operator (PINO) surrogate for two-dimensional compressible, viscous, resistive reconnection in a wall bounded domain. We condition on the initial state, the Lundquist number, and a continuous query time. Predicting the magnetic flux function makes $\grad\!\cdot\!\bB=0$ exact, direct time queries remove autoregressive error accumulation, and parity-aware spectral differentiation matches the wall treatment of the reference solver. Trained across $10^{3}\le S\le2\times10^{5}$ with two values withheld, the surrogate recovers density, pressure, and guide field to better than $0.2\%$, velocity to $3.2$ to $3.9\%$, the flux function to $1.0\%$, and the spectral current density to $7.5\%$. It resolves the thin, narrow current sheet, Alfvénic jets, and spectral resolution beyond mode cutoff. Our surrogate reproduces the layer-averaged reconnection rate to within $3\%$ through $S=1.1\times10^{4}$, recovering a reconnection time exponent of $0.497$ consistent with Sweet-Parker scaling. When queried zero-shot on the $2049^{2}$ solver mesh, a full trajectory run is on the order of seconds on one GPU against several hours for the DNS, over two orders of magnitude faster, establishing learned operators as a practical route to surveying reconnection across parameter space.
Aluminum is a reference standard in high-energy-density research and serves as a liner material in megampere Z-pinch facilities. Using mixed deterministic-stochastic finite-temperature density functional theory, we compute the equation of state and transport coefficients of liquid aluminum up to temperatures of 1000 eV and compare the results with model-based approaches. We find that for T<200 eV the density-functional results differ from the models by more than 10% at high density, with a maximum deviation of over 30% (at T = 10 eV), leading to significant discrepancies in the Hugoniot curve at high compression; the calculated electrical and thermal conductivities deviate from model predictions by 26% to 63% at T~100 eV. We further employ the Drude model and the Epperlein-Haines framework to examine the magnetic-field dependence of the conductivities. These results can provide valuable input data for radiation-hydrodynamics codes.
Integrated codes simulating interactions between Alfvén waves and fast ions in tokamak plasmas use perturbative models for the relatively slow processes of instability growth, saturation, chirping and bursting, and transport. Faster processes by which an Alfvén mode's spatiotemporal structure forms are assumed to have been completed within the mode's oscillation period, $τ_0 \equiv 2π/ω_0$. This separation of time scales underlies the computational efficiency of perturbative models, where the Alfvén mode's time-dependence is reduced to that of a scalar signal $s(t) = A(t)\sin(-ω_0 t - φ(t))$ with variable amplitude $A(t)$ and phase $φ(t)$. For this, accurate input data in the form of a mode's spatial structure $δΦ({\mathbf x})$, damping rate $γ_{\rm d}$, and initial frequency $ω_0$ are required. For modes residing in dense or continuous spectra, $δΦ$ and $γ_{\rm d}$ could be estimated from the form of the continua and fast ion orbits, but it is difficult to guess the seed frequency $ω_0$. Here, we report results of numerical experiments showing that it is possible to find $ω_0$ using a prompt frequency shift that occurs during the first few $100$ time steps of a simulation. Restarts with the shifted frequency iteratively converge to a value of $ω_0$ that seems to maximize the resonant drive, suggesting an auto-optimization process. The need for iteration is attributed to the fact that the terms required for rapid frequency adjustments were truncated when deriving the perturbative model. Meanwhile, the fact that partial auto-optimization is possible at all is attributed to the fact that series truncation alone (without filter) does not strictly enforce slowness. Remnants of and cross-talk with faster dynamics still occur in numerical implementations. This entails potential for both uncertainty and utility.
Quasilinear models make gyrokinetic turbulent-transport predictions sufficiently fast for integrated modelling, but their predictive capability is limited by two factors: the physical and geometrical applicability of the linear solver, and the validity of the saturation rule used to close the model. We present the Predictive Quasilinear Solver (PQLS), a quasi- linear gyrokinetic transport solver formulated in general magnetic geometry. Its implementation as an eigenvalue solver retains electromagnetic and collisional effects, provides access to dominant and subdominant modes and is differen- tiable with respect to all plasma parameters. Linear benchmarks against GENE reproduce the growth rates, frequencies, and eigenfunctions. We additionally formulate the saturation-rule closure as a Bayesian inference problem that distin- guishes uncertainty in its fitted coefficients from the residual model-form uncertainty. The approach is demonstrated by calibrating the SAT3 rule on PQLS quasilinear weights against published nonlinear CGYRO cases. In addition to improving the robustness of the calibration, the new method also quantifies the uncertainty in each of the fit coefficients. Such uncertainty is propagated through transport calculations to produce error-aware profiles that are compared to the ones obtained from the full gyrokinetic simulation, showing excellent agreement.
Discharge ignition events in dielectric barrier discharges (DBDs) self-organise into spatio-temporal patterns with varying degrees of order. The complex dynamics of a DBD and intricate structure of occurring patterns complicate the formulation of predictive, mechanistic descriptions. We present the formulation of a reduced-order model that describes the re-ignition dynamics between consecutive discharges appearing at the same position inside a DBD arrangement. The model is derived from an equivalent electric circuit and validated against fluid-Poisson simulations and experiments performed with a multi-filament arrangement in air-like gas mixtures at atmospheric pressure driven by sinusoidal high-voltage waveforms. The experimental scenarios include a highly ordered regime where discharges ignite at regular time and space intervals generating a pattern stable over several periods, and an unstable regime with discharges appearing at seemingly random positions and times. The model accuracy is assessed in both regimes and it is found that the associated prediction uncertainty provides a quantitative measure of the spatial order of the discharge pattern.
We develop and implement the quadratic density response function of non-interacting fermions at arbitrary temperature, frequencies, and wave vectors. Starting from a Green's function formulation, we derive the quadratic response and demonstrate its equivalence to the result obtained from the Wigner equation. We further derive the classical limit through a perturbative expansion of the Vlasov equation and demonstrate that the quantum and classical formulations agree in the high-temperature limit. We analyse the limiting behaviour with respect to wavenumber and derive the zeroth harmonic response. Two independent implementations are provided and extensively benchmarked against density-functional theory, canonical path integral Monte Carlo (PIMC), and grand canonical PIMC simulations. As the density response of the interacting electron gas is commonly modelled through the ideal response functions and approximate models for the local field correction, the presented formulation will also allow for more complete explorations of interacting systems. Especially, our efficient implementation, which evaluates the ideal static and dynamic quadratic response functions in less than 0.5 ms on a 1.3 GHz processor, will enable evaluation of quadratic corrections to integrated quantities such as interaction potentials and stopping powers in warm dense matter.
We present molecular dynamics simulation results for an ultracold $^{40}$Ca plasma in a constant, homogeneous magnetic field. The magnetic field induces a significant spatial separation of charges. Furthermore, a significant fraction of the electrons escapes from the region bounded by the ions, in contrast to the case of plasma expansion into a vacuum. This leads to the formation of a quasineutral plasma core containing the remaining electrons and surrounded by an outer, thick ion shell, where the ion density drops sharply. Anisotropy in the ion kinetic energy is also observed, depending on the direction of ion motion relative to the magnetic field lines. The simulation results are in agreement with experimental data on ion dynamics.
Previous work built AI-based surrogates for a nonlinear gyrokinetic simulation code with the goal of using them for fast, direct calculations of turbulent ion heat flux in stellarator design optimizations and scenario planning for experiments. These AI surrogates were trained on data from >200k nonlinear, adiabatic electron gyrokinetic simulations with the gyrokinetic flux-tube code GX, using a wide range of stellarator magnetic configurations ($\sim$23k), positions in the plasma, and gradient scale lengths. In this paper, we demonstrate the use of the AI-based turbulence surrogate in the optimization of stellarator magnetic equilibrium and to speed up stellarator transport solvers. Due to its speed ($\sim$ms), the AI-based surrogate enables previously unattainable optimization objectives, such as full radial profiles of ion turbulent heat flux, or directly optimizing to maximize the turbulent critical gradient at multiple locations across the plasma. These direct calculations provide a potentially more accurate optimization target and reduce reliance on ad-hoc heuristics that may not accurately capture the variation of turbulent transport with magnetic configuration. By including the AI-based surrogate for turbulent heat flux in a transport solver, we can quickly postprocess and confirm the improved ion temperature resulting from the optimized equilibrium. Finally, we demonstrate the use of AI agents with strong reasoning AI models to automate the outer loop, exploring many objective and hyperparameter configurations with this AI-based turbulence surrogate to discover improved turbulence optimized magnetic equilibria.
Resistive switching in oxide-based devices is widely governed by stochastic defect processes, yet a predictive link between fabrication conditions and functional behavior remains elusive. Here, we establish a multiscale framework connecting plasma-defined deposition conditions to macroscopic device functionality in sputtered SiO$_x$/Cu/SiO$_x$-based systems. By combining large-scale statistical analysis of more than 50,000 experimentally characterized devices with physics-based plasma and atomistic simulations, we show that device behavior does not emerge from deterministic process-to-performance mappings, but from a probabilistic cascade spanning defect formation, defect-state evolution, and functional-regime emergence. Data-driven clustering reveals a continuous functional state space composed of operational switching types, while inverse modeling identifies the reconstructed oxygen-vacancy density as an effective latent descriptor capturing the combined influence of structural disorder and defect topology. This latent descriptor is strongly coupled to both Cu redistribution and electrical response, linking otherwise hidden material properties to observable device characteristics. Furthermore, macroscopic switching behavior is argued to arise from ensemble integration across spatially heterogeneous subdomains, providing a physical explanation for the pronounced variability of large-area devices. These findings shift the perspective from deterministic defect engineering toward probabilistic defect-state design and establish a physically grounded framework for understanding and controlling functional variability in such oxide-based systems, such as memristive or resistive-switching devices.
The production of multiply charged titanium ions is of significant interest for accelerator-based experiments, surface modification technologies, and applications in nuclear physics. Due to the low saturated vapor pressure of titanium at temperatures below 1000 °C, the generation of titanium ion beams from an electron cyclotron resonance ion source (ECRIS) remains a challenging task. The choice of titanium vapor injection methods into the ECR plasma plays a key role in achieving high ionization efficiency, plasma stability, and reliable long-term source operation. In this work, alternative approaches for the injection of neutral titanium atoms into the ECR plasma are investigated, focusing on the use of fluorine-containing compounds. Titanium fluorides (TiF3, TiF4) possess a relatively high saturated vapor pressure at moderate temperatures, allowing to use standard resistively heated ovens operating up to 1000 °C. Another method involves the in-situ formation of titanium fluorides inside the plasma chamber via chemical reactions between metal titanium and the dissociation products of sulfur hexafluoride (SF6). These approaches enable controlled and efficient titanium injection into the ECRIS plasma with satisfactory extracted ion beam stability.
Monte Carlo sampling of strong-field quantum electrodynamics processes underpins simulations of high-intensity laser experiments and of astrophysical compact-object magnetospheres. Sampling an event requires the total rate of the process together with the cumulative probability that determines how energy is partitioned between the produced particles. Simulations typically tabulate both in advance and invert the tabulated probability numerically. Here we replace this procedure with elementary-function approximations for synchrotron radiation and the nonlinear Breit--Wheeler process. For each process, we approximate the auxiliary function that sets the total rate, as well as the cumulative probability, with Padé approximants chosen so that the inversion reduces to a quartic equation. This yields the sampled quantum parameter---electron $χ_e$ or photon $χ_γ$---in closed form. The approximations and the particle spectra sampled from them agree with the exact results to within $1\%$. The procedure requires no lookup tables, no interpolation, and no numerical root finding, and can be inserted directly into radiative particle-in-cell codes.
Artificially generated muon beams, with tunable parameters, are highly valuable in imaging objects that are opaque to any other scanning technique. Lower-energy muon beams are currently produced at accelerator facilities, while higher energies have only recently been observed as a result of the interaction of laser-wakefield-accelerated electrons and solid targets. This work focuses on using muography imaging detectors to measure the profile of a muon beam generated via the Bethe-Heitler process and to perform muon imaging in an experimental setup extending up to 42 meters from the muon source. By comparing the data with Monte Carlo simulations, we confirm that the measured beam profile and object imaging are consistent with artificially produced muons of energy of a few GeV. To our knowledge, this represents the first demonstration of imaging dominated by an artificial, laser-driven muon beam.
The axion coupling to electromagnetism may leave signatures of uncertain nature in astrophysical spectra; one is thought to be a single line centered at the axion frequency $ω= \sqrt{\left(m_a c^2/\hbar\right)^2 + c^2 k_a^2}$, with $m_a c^2/\hbar$ the axion Compton frequency and $\hbar k_a$ its momentum. If a parametric plasma instability is excited, this signature is transformed into a set of lines, which may be proportionally spaced and linearly polarized. A Crab pulsar signal at a frequency near recent cosmological predictions for the axion Compton frequency has some qualitative similarity to this signature, though this similarity is not evidence that the signal is caused by axions.