2026-09-04 | | Total: 9
We develop a quasi-cylindrical direction-splitting (QDS) dispersionless Maxwell solver for the open-source EPOCH particle-in-cell framework. The method preserves the exact axial phase velocity for every retained azimuthal mode by advecting transverse transport variables with the $cΔt=Δx$ shift. Building on the same transport viewpoint, we formulate an axion solver that advances the Klein--Gordon equation and the axion-regenerated electromagnetic fields within the QDS solver, so that axion--photon coupling is handled self-consistently. Benchmarks demonstrate that the cylindrical QDS solver reproduces the standard Yee wakefield structure while removing the axial group-velocity error and reducing computational cost by orders of magnitude. In the resonant axion generation test via two-color laser mixing, the QDS solver eliminates the spurious vacuum source term and maintains phase matching, whereas the Yee solver suffers from rapid dephasing due to numerical dispersion. The proposed algorithm therefore enables quasi-3D simulations of axion--photon physics in laser-plasma interaction with dispersion-free laser transport.
The nonlinear evolution and mutual interaction of two-dimensional electron acoustic (EA) nonlinear structures in superthermal plasma environment are studied. The plasma model consists of inertial cold electrons, superthermal hot electrons described by kappa ($κ$) distribution, and stationary ions providing overall charge neutrality. Using the extended Poincaré-Lighthill-Kuo (PLK) reductive perturbation technique, a pair of two-sided Kadomtsev-Petviashvili (KP) equations governing right- and left-propagating EA solitary waves (EASWs) is derived. Exact analytical solutions of the KP equations, including single soliton, multisoliton, breather, and lump structures, are obtained via the Hirota bilinear method. The effects of key plasma parameters such as hot electron concentration, temperature ratio, and superthermality index on the characteristics of these nonlinear excitations are examined. Particular attention is devoted to the head-on collision dynamics between solitons, breather-soliton, and breather-breather interactions. The results reveal quasi-elastic collisions accompanied by phase shifts, transient amplitude modulation, and localized energy concentration, with clear distinctions between oscillatory and non-oscillatory mode interactions. The present study provides new insights into multidimensional electron acoustic wave (EAW) dynamics and energy redistribution mechanisms in superthermal space plasmas, with direct relevance to planetary magnetospheric environments such as Saturn's ring region.
Laser wakefield acceleration (LWFA) is an advanced method of high gradient acceleration of charged particles with wakefield excited in plasma by laser pulse. A distinctive feature of this method is the ability to create and accelerate so-called self-injected bunches with unique parameters - charge, energy, emittance, and geometric dimensions - without the need for an external injector of the required bunches for their subsequent acceleration to higher energies. Self-injected bunches emerge due to the plasma electrons trapping by the excited wakefield (self-injection phenomenon). For many applications, self-injected bunches can be used directly in relevant experiments. In this paper the dependence of the self-injected bunch parameters on the laser pulse amplitude and longitudinal gradient of plasma density in tapered plasma channel is investigated using numerical simulation with the WarpX code. At the laser amplitude $a_0=E_0(m_{e0}cω/e)^{-1}=3.6$ in a conical channel with the radius decreasing from $4.0\,c/ω_{pe}$ to $2.16\,c/ω_{pe}$ and with the plasma density increasing linearly from $n_e=2.61\cdot10^{19}$ cm$^{-3}$ at the channel entrance to $3n_e$ at the exit, a self-injected bunch is obtained with the charge $32.1\,μ$C/m, the mean longitudinal momentum $111.9\,m_ec$, the length $5.10\,μ$m, the area $5.3\,μ$m$^2$, the transverse emittance $1.6\cdot10^{-2}$ mm$\cdot$mrad.
Global gyrokinetic particle simulations remain computationally expensive, as they demand both adequate marker statistics and three-dimensional field solvers. In this work, we present a hybrid spectral method within the particle-in-Fourier (PIF) framework and implement it in the electrostatic model of GTC. Charge scatter and field gather are performed between particles and fields on a two-dimensional poloidal mesh, while the corresponding Poisson solver is discretized using radial finite differences and poloidal $m$-harmonics. Truncated spectral transforms are employed to connect multiple representations for fields, avoiding costly particle-grid operations for each individual $m$-harmonic within the particle loop. Benchmarks against conventional particle-in-cell (PIC) simulations successfully reproduce single-$n$ ion temperature gradient (ITG) mode structures and dispersion relations, as well as multi-$n$ nonlinear ITG transport and its regulation by zonal flows. Compared to conventional PIC, the proposed method reduces the effective problem size by more than a factor of 48 and achieves a speedup of over two orders of magnitude for single-$n$ cases. A 2000-step single-$n$ simulation with approximately 2 million markers completes in 78.2 seconds on a laptop GPU, while multi-$n$ turbulence simulation also completes within minutes. Furthermore, the elimination of toroidal particle-shift communication yields promising preliminary scaling performance on multiple NVIDIA A100 GPUs. The numerical scheme is broadly applicable for accelerating particle simulations on platforms ranging from laptops to supercomputers.
The origin of magnetic fields in turbulent astrophysical systems has long been a central problem in plasma astrophysics. Fluctuation dynamos are a class of field amplification mechanisms that occur in turbulent magnetohydrodynamics whereby stochastically forced motions of plasma at sufficiently high magnetic Reynolds numbers exponentially amplify magnetic energy. For steady forcing, such dynamos saturate with magnetic energies at a sizable fraction of the turbulent kinetic energy. Although fluctuation dynamo is widely studied within the framework of resistive magnetohydrodynamics (MHD) and driven-turbulence numerical simulations, this work explores the difference when the Hall term is included in the magneto-fluid's generalized Ohm's law. The inclusion is motivated in part by recent high energy-density plasma experiments studying fluctuation dynamo that are governed by an extended magnetohydrodynamics (xMHD) ansatz, which includes the Hall term. We first discuss the details of the Hall-MHD implementation in the FLASH code, the tool we use to model xMHD fluctuation dynamo. We then investigate the influence of the Hall term on the fluctuation dynamo in a three-dimensional periodic box, driven with stochastic forcing at the box scale. We compare cases with a Hall term of varying magnitude to no-Hall cases with respect to the magnetic field growth rate, saturation level, and magnetic field structure. The Hall-MHD fluctuation dynamo is found to saturate at lower magnetic energies and with fewer small-scale magnetic structures than the no-Hall cases. Both findings are consistent with the interpretation that the Hall term acts as an additional, non-linear transport term, akin to an enhanced turbulent diffusivity.
Expensive fusion simulations are commonly preserved in code-specific formats that limit discovery, comparison, and reuse. We present the Multiscale GyroKinetic DataBase (MGKDB), an open-source software framework and curated archive that converts heterogeneous simulation campaigns into traceable scientific records. Each record links code-native inputs and outputs to provenance and quality metadata, an IMAS-aligned physics representation, and derived diagnostics, preserving model-specific evidence while enabling common-field queries. Production pathways support linear and nonlinear GENE and CGYRO calculations and reduced quasilinear TGLF evaluations. At the September 1, 2026 snapshot, MGKDB contained 1,068,089 records, nearly all of which included a populated gyrokinetics IMAS branch. The software is openly available, while access to the NERSC-hosted production records is managed. Three demonstrations show how these linked representations support scientific reuse. Standardized quantities stored in the Diagnostics branch enable population-scale analysis of archived linear modes; common input coordinates reveal coverage, redundancy, and campaign-driven sampling structure across a multicode collection; and record-level retrieval of native CGYRO inputs drives matched TGLF calculations and produces a traceable dataset for exploratory surrogate modeling. Together, these examples demonstrate how MGKDB supports archive characterization, candidate cross-code and cross-fidelity comparisons, campaign planning, and reproducible data-driven modeling without treating different models as automatically equivalent.
Physics-informed neural networks for high-order partial differential equations require mixed input derivatives and their gradients with respect to the network parameters. Standard implementations obtain an order-$K$ derivative by repeated automatic differentiation. We instead organize the forward recursion of the multivariate Faà di Bruno formula and its explicit backpropagation over a prescribed downward-closed set of multi-indices, with the Bell-polynomial convolutions tabulated once. The forward pass carries only the derivatives the differential operator requires. The backward pass propagates gradients from losses formed from any subset of them, including nonlinear products and coupled fields. Both recursions are exact up to roundoff and avoid nested computational graphs. An independent Taylor-jet implementation, symbolic checks of the test problems, and finite differences verify derivatives and loss gradients through order seven. On one CPU core, the method evaluates 330 mixed derivatives with respect to four inputs through order seven and the corresponding loss gradient without the memory failures observed for several nested implementations. Numerical tests include third-, fifth-, and seventh-order dispersive equations, incompressible flow, and a manufactured five-field electrohydrodynamic system. The seventh-order Zakharov-Kuznetsov test in $3+1$ dimensions has a relative solution error of $6\times10^{-4}$.
For a simple model fluid, the Yukawa fluid, the condition for dynamic crossover, known as the Frenkel line, is defined. The condition is related to the fact that the roton minimum in the dispersion relation of longitudinal acoustic-like excitations exists only when the collective vibrational dynamics of particles dominates in the liquid. Based on the self-consistent relaxation theory for the Yukawa fluid, thermodynamic states are determined in which the roton minimum disappears. The obtained values of the state parameters for the Frenkel line are consistent with the results of studies in which the position of this line on the phase diagram of the Yukawa fluid was determined using molecular dynamics simulations. It is shown that the Frenkel line in this system can be determined directly from the structural characteristic -- the static structure factor. A physical interpretation of the roton minimum frequency for the simple liquids near the Frenkel line is proposed.
Learning solution operators across broad parameter ranges can require substantial coverage of both input functions and physical parameters, particularly for purely data-driven parametric models. In addition, the resulting models may fail silently outside the training distribution. We introduce equation recast, which reformulates parametric operator learning as the learning of a single canonical operator. Parameter-induced operator variations are derived analytically from the governing equation and absorbed into effective sources, enabling zero-shot prediction across new parameter regimes. Across multi-parameter, nonlinear, and singular PDE settings, equation recast supports extrapolation, integrates sparse heterogeneous datasets in a shared canonical representation, and uses loss of convergence as an internal warning signal for failure of the recast iteration. In high-fidelity tokamak simulations for nuclear fusion, the framework unifies electron-temperature data across four device geometries through canonical-domain mapping within one jointly trained operator. Equation recast provides a route toward reusable neural PDE solvers combining equation-guided transfer, data efficiency, and monitorable inference.