2026-07-27 | | Total: 7
We present the first direct experimental observation of fishbone-driven zonal flows in the core of the EAST tokamak. In contrast to the global pattern predicted by previous models and simulations based on the energetic-particle-expulsion mechanism, the observed flows exhibit a fine-scale, radially reversed structure inside the q = 1 rational surface. The flow rises faster and saturates earlier than the fishbone within a single burst, indicating that a beat-driven nonlinear process dominates the early stage rather than the energetic-particle-expulsion mechanism. Global nonlinear gyrokinetic simulations quantitatively reproduce the observed radial profile and reveal that this structure arises from the cancellation of comparable but opposite contributions from thermal ions and electrons. This cancellation mechanism is not captured in previous theoretical frameworks. These findings establish that the fishbone can generate sheared flows with a distinct radial topology, offering a promising pathway for regulating turbulence and improving core confinement.
Results of a full 3D electromagnetic plasma fusion simulator are presented. It will enable investigations into the influence of structured electromagnetic potentials, anomalous magnetic moments, or charge cluster structures on the fusion yield. The simulator calculates fusion yields similar to those found in real world fusors. In addition, the ion energy and density are accurately simulated. Given that the ion properties arise from a full Maxwell solver there is ample evidence that this part of the simulation is accurate. Concluded with a recommendation for continued research.
Discrepancies between experimental data and radiation-hydrodynamic models of burning plasmas at the National Ignition Facility have been attributed, in part, to possible deviations of reactant ion distributions from Maxwellian equilibrium. In particular, it has been hypothesized that the collisional relaxation of energetic fusion products with the bulk plasma may generate suprathermal ion populations not captured by reduced models. To assess this hypothesis in a fully kinetic setting, we present a multi-species reactive--elastic kinetic framework, in which fusion sources and sinks are evaluated directly from the reactant velocity distribution function with a reactive Boltzmann collision operator. For computational efficiency, elastic interactions are modeled in the small-angle collision limit using the Landau and Lenard--Bernstein collision operators, although the framework admits large-angle elastic scattering generalizations. The integro-differential collision operators are discretized with a fast Fourier spectral method, enabling efficient, high-order accuracy evaluation of reactive processes in three-dimensional velocity space and providing a deterministic alternative to Monte Carlo collision methods. Numerical experiments are performed for spatially homogeneous two- and three-species D--D fusion systems in weakly coupled, early-time regimes. We find that fusion product heating does not generate a significant suprathermal population in the reactant ions: deviation from the corresponding Maxwellian distribution remains small in both the Frobenius norm and the fusion reactivity.
The blowout regime of plasma wakefield acceleration is a strongly nonlinear plasma response in which an intense relativistic beam or laser pulse expels plasma electrons forming a nearly uniform ion column. Existing self-consistent theories describe this structure primarily in the axisymmetric limit, reducing the bubble dynamics to a single radial degree of freedom. This restriction is inadequate for flat beams with highly asymmetric transverse emittances, as required in future linear colliders to control luminosity and beamstrahlung. Using a reduced Lagrangian formulation, we generalize the existing theory beyond axisymmetry and derive a coupled system of equations governing the wake potential driven by an asymmetric beam. This theory self-consistently predicts the focusing fields and their variation along the co-moving coordinate, as well as the accelerating field and its transverse nonuniformity. The resulting theoretical predictions are benchmarked against particle-in-cell simulations, demonstrating the theory's ability to describe and interpret wakefields driven by nonaxisymmetric sources.
Rotational symmetry protects the topological charge of on-axis electron vortices but not of off-axis vortices. We identify an additional SU(1,1) dynamical invariant that guarantees conservation of their intrinsic orbital angular momentum within the near-axis approximation. First-principles simulations of an off-axis electron vortex traversing a Glaser lens confirm this prediction, establishing a robust transport mechanism in axisymmetric nonuniform magnetic fields.
We derive relativistic resistive magnetohydrodynamics framework for a two-component ultrarelativistic plasma of massless, oppositely charged particles directly from the Boltzmann-Vlasov equation using the 14-moment approximation. The resulting second-order equations couple the net-charge diffusion current to the shear-stress tensor through the electric field, with all transport coefficients given in closed microscopic form. In the homogeneous limit, the charge-current dynamics is well described by relaxation-type Ohm's law for moderate field strengths, while large viscosity drives the system into an underdamped oscillatory regime absent from the standard Israel-Stewart formulation. Most strikingly, a purely electric field generates sizable momentum anisotropy even without any underlying flow gradient. Under Bjorken expansion this field-induced anisotropy persists but becomes subleading to the hydrodynamic expansion source.
An elastica knot is defined in terms of the Frenet-Serret curvature $κ(s,t)$ as a function of the arclength $s$ along the spatial curve ${\bf r}(s,t)$ at a fixed time $t$, which is a solution of the curvature differential equation $\partial^{2}_{s}κ(s,t) = -\;κ^{3}/2 + k_{0}^{4}τ_{0}^{2}\;κ^{-3} + λ\,k_{0}^{2}κ/2$ that is obtained from a variational principle that minimizes the bending energy of the spatial curve under the constraint of a constant curve length. Here, the Frenet-Serret torsion $τ(s,t)$ satisfies the conservation law $κ^{2}(s,t)\,τ(s,t) \equiv k_{0}^{2}\,τ_{0}$, while $λ$ is a constant of integration. After briefly reviewing the Hasimoto transformation from a space curve ${\bf r}(s,t)$ to the nonlinear Schrödinger equation (NLSE) $-\,iD^{-1}\partial_{t}ψ= \partial^{2}_{s}ψ+ \frac{1}{2}\,|ψ|^{2}ψ$, where the constant $D$ has units of fluid circulation (m$^{2}$/sec), we show how the traveling-wave solution $ψ(s,t) = Ψ(s_{t} \equiv s - c\,t) \equiv κ(s_{t})\;\exp[iθ(s_{t})]$ is mapped onto the curvature equation for an elastica knot, with $θ^{\prime}(s_{t}) \equiv c/(2D) + k_{0}^{2}τ_{0}/κ^{2}(s_{t})$ and the elastica-knot constant $k_{0}^{2}λ= -\frac{1}{2}\,(c/D)^{2}$ expressed in terms of the traveling-wave NLSE parameters $(c,D)$. The constraint of a closed 3D elastica knot imposes spatial periodicity conditions that introduce a unique set of knot parameters for which the NLSE traveling wave can exist. The present work shows that the traveling-wave solution on a closed elastica knot requires an extension of the classical elastica-knot parameter space.