2026-11-06 | | Total: 9
Several closed-form non-axisymmetric magnetic fields with nested toroidal flux surfaces have recently been discovered as exact benchmarks for verifying three-dimensional magnetohydrodynamic equilibrium codes. However, all of these solutions have at most two field periods ($N \le 2$), i.e., identical segments that repeat toroidally. Here, we derive closed-form non-axisymmetric magnetohydrodynamic equilibria for an arbitrary integer number of field periods $N \geq 1$. The equilibria have non-constant pressure and closed magnetic field lines with either zero rotational transform ($ι= 0$) for $N \geq 1$ or integer rotational transform ($ι= \pm N$) for $N \ge 3$. For the reflection-symmetric $ι= 0$ equilibria, we also construct a matched exterior vacuum field with continuous magnetic field across the plasma boundary, providing a benchmark for free-boundary equilibrium codes.
We present an analytical stability analysis of the laser hosing instability in a parabolic channel in the weakly relativistic regime. We employ a variational principle approach to derive a set of coupled equations describing laser dynamics in a transversely inhomogeneous plasma. The quadratic Lagrangian formalism is applied to perform a stability analysis and to derive the growth rates for hosing in both long- and short-wavelength regimes. We find that a parabolic channel suppresses the long-wavelength hosing mode completely, for any channel depth, and that the unstable band of the short-wavelength mode is bounded below by a cutoff wavenumber set by the ratio of the channel frequency to the hosing coupling. A steeper channel narrows the unstable band and lowers the short-wavelength growth rate, but cannot stabilize the mode, since the dispersion relation forbids damping identically. The analysis is accompanied by two-dimensional particle-in-cell simulations, which reproduce the prediction in the nonlinear regime: the long-wavelength mode is absent in a channel, the surviving modulation shortens as the channel steepens, and its amplitude falls with channel depth. Further, our simulation results demonstrate the sensitivity of off-axis laser propagation to the plasma channel and distinguish the hosing instability from laser oscillations due to off-axis laser propagation in a plasma channel.
It has been known for some decades that the plasma in the Solar System is invariably in a state away from thermal equilibrium. In several plasma environments, be it in the Solar Wind or within planetary magnetospheres, the plasma has been found in a quasi-stationary nonthermal state, with a substantial fraction of its high-energy particles fitting along power-law distributions, instead of the exponential decay of thermal systems. Translating the observations in energy to velocity distribution functions (VDFs), the Kappa plasma model has been employed by the space plasma community, with the (standard) Kappa VDF (SKD) being the preferred choice. However, recent works have called attention to the limitations of the SKD model, proposing as a viable alternative the regularized Kappa VDF (RKD) model, which corrects the deficiencies of the SKD. The RKD model introduces a Gaussian-like dependency on the very high-energy part of the VDF that effectively cut-offs, via a regularization parameter, the Lorentzian profile of the intermediate-energy range. In this work, we employ a generalized RKD that supports anisotropies both in the thermal speed parameters and in the regularization parameters, including as well the possibility of drifts parallel to the magnetic field. This generalized RKD (called the bi-regularized Kappa with drift distribution, bRKd), is employed to derive for the first time the components of the susceptibility tensor for a suprathermal plasma species/population modelled by the bRKd. The tensor components are valid for waves with any propagation angle, frequency range and polarization. During the derivation, new special functions are defined and analyzed, which generalize previous expressions restricted to the parallel propagation limit.
The shear Alfvén continuum has been studied extensively in tokamaks, but its gap structure in modern optimized stellarators has received less attention, despite its importance for energetic-particle confinement in fusion pilot plants. We compute and compare the continua of five reactor-scale stellarator equilibria: two quasi-axisymmetric (QA), two quasi-helically symmetric (QH), and one quasi-isodynamic (QI). We present selective Fourier mode schemes that provide a sufficient resolution of eigenfrequencies within the frequency range of interest. In all five configurations, the widest gap is the HAE(2,$N_{\text{fp}}$) gap, which is driven by rotating ellipticity. This gap lies in the frequency range of strong resonance with prompt fusion-born alphas in the QA configurations but not in the QH configurations, suggesting that QH configurations may be easier to optimize for alpha confinement. Near-axis theory predicts the HAE(2,$N_{\text{fp}}$) and EAE gap widths well but overestimates the MAE gap widths and the gap widths of the QI configuration. These discrepancies likely arise from interactions between gaps, which depend on the rotational transform profile and are not captured by near-axis theory.
Quench propagation in cable-in-conduit conductors is thermo-hydraulically driven, yet existing 1D models treat a conductor as one straight adiabatic cable with prescribed pressure and mass-flow boundaries, so secondary quenches and flow redistribution between windings cannot be resolved. We extend the open-source code OPENSC$^2$ with transverse heat conduction between turns, a 3D field resolving the current-sharing temperature, induced losses, and a lumped hydraulic network coupled to the conductors by a mass-flow formulation. Validated against the 2006 W7-X coil tests, it predicts the quench site, cuts the outlet-temperature error from 55 K to 12 K, and reproduces the neighboring windings' heating as flow reversal through the coil inlet.
We construct a divergence-free velocity field $\mathbf{u}(\mathbf{x},t)$ on $\mathbb{T}^{n}$ that is a genuine function of position and carries both defining properties of turbulent transport: Richardson pair dispersion and a dissipative anomaly. No synthetic field has carried both, because the frame that gives the first destroys the field the second requires. The key is a dichotomy: deforming a spectrally localised pattern smears it across scales, whereas pushing its generating centres forward and re-convolving leaves the spectral envelope exactly invariant --- so transporting the eddy centres preserves both cascade and frame. In the field the finite-size Lyapunov exponent obeys $ζ(R)=h(R)-1$ threshold by threshold within the scaling range, with a sweeping deficit $Δ_{\rm sw}=+0.038\pm0.019$, consistent with zero, against $+0.091\pm0.019$ for the frozen field on the same data, a gap that doubles at Navier--Stokes intermittency. Companion papers prove that this family shows Richardson dispersion and a dissipative anomaly (Fefferman and Mayboroda 2026), and that a scalar advected by this flow obeys Yaglom's relation (De Lellis 2026), at small Kubo number. Trajectories forget their starting point as the smallest eddy is removed: the median separation of pairs released from one point with independent Brownian kicks varies as $κ_{p}^{0.0333\pm0.0038}$ over 22 decades of particle diffusivity, against $κ_{p}^{1/2}$ for a smooth control. A scalar in the same field closes the exact scale-by-scale variance budget, four independently measured terms balancing the dissipation to 15\% over 1.2 decades with no adjustable constant, and reaches its inertial--convective asymptote $S_{L}=-χr$ to $0.9995\pm0.0018$ once the measured diffusive transfer is restored, while under an imposed mean gradient $χ$ changes by 14\% over a factor 5.8 in $κ$.
The numerical simulation of kinetic plasma models in six-dimensional phase space poses significant computational and memory challenges, particularly when considering modern high-performance computing architectures. In this work, we propose an adaptive dynamical low-rank approach for the efficient solution of the 6D Vlasov--Poisson system. Our strategy is an extension of the projector-splitting dynamical low-rank integrator. The algorithm changes the rank based on physical observables (e.g., the electric field and energy) as opposed to controlling the error in the distribution function. Our approach is memory efficient in the sense that only rank $r$ objects are required to implement the scheme. This makes the algorithm suitable for performing high-dimensional computations on graphic processing units. Comparisons with rank adaptive Basis-Update & Galerkin integrators exhibit both significant memory savings and reduction in run time. We show, by presenting a variety of 6D benchmark experiments implemented in the C++/CUDA low-rank framework Ensign, that our technique allows us to run simulations with a significantly smaller rank, while still preserving the physical fidelity. Thus, the approach substantially reduces memory footprint and computational cost and enables simulations that are prohibitively expensive for full-rank or less memory-efficient low-rank schemes.
Sample-based quantum subspace methods construct a small subspace from measured basis states and solve the target problem classically within it. We study these methods for heat conduction in magnetized plasmas, where the anisotropy $χ_\parallel/χ_\perp$ reaches $10^8$ to $10^{10}$. We prove two results for the energy-norm Galerkin solution. First, if the measurement basis is aligned with the flux surfaces and the subspace contains the poloidal-mode-zero block, the error is bounded by $\|r_0\|_{L^+}(χ_\parallel-χ_\perp)^{-1/2}$, where the numerator does not depend on $χ_\parallel$. The bound is sharp, so the accuracy improves as the anisotropy grows. Second, in the computational basis, a subspace that misses a single node of a flux surface suppresses the flux-function component there as $χ_\parallel^{-1/2}$, which gives an error floor when the limiting solution is nonzero on that surface. We verify both results on a shaped cross-section and on a KSTAR-like Miller equilibrium with ten qubits. At $χ_\parallel/χ_\perp=10^9$, the computational basis gives an error of order unity for every shot budget, whereas a quantum Fourier transform on the poloidal register gives $1.3\times10^{-5}$. The same barrier limits sample-based quantum diagonalization. The aligned measurement needs no time evolution, so the circuit is only state preparation and a five-qubit quantum Fourier transform. On a ten-qubit superconducting processor, the aligned basis gives a median error of $7.0\times10^{-6}$ with 768 shots, whereas the computational basis gives errors above 0.94 up to $1.2\times10^4$ shots. A trapped-ion processor shows the same behavior.
The integration of a theory regarding the dynamics of averaged magnetic moment in a turbulent setting with the concept of a self-optimized cascade loop of helical fluxes, either spontaneously or intentionally generated near the separatrix, enables the derivation of spectral laws for the floating potential and ion saturation current, which align with findings from various experiments conducted on tokamaks, stellarators, and RFX-mod reversed field pinches. The notion of distributed chaos enables a quantitative evaluation of the randomness levels of chaotic/turbulent states observed inside and outside the separatrix, linking them to self-optimized helical fluxes.