2026-05-06 | | Total: 15
The anisotropic diffusion of electrons with various rigidity and the anomalous conductivity of a collisionless plasma in the presence of Weibel-generated quasi-static turbulent and uniform external magnetic fields are examined. Using an original code based on the Boris algorithm, the electron diffusion coefficients and the longitudinal, transverse, and Hall mobility factors are determined for a representative set of plasma parameters. It is shown that these values and their anisotropy depend strongly on the electron temperature, external magnetic field, average level of magnetic turbulence, and its spectrum. The physical origin and expected limits of such dependencies are indicated. Applications of the results are discussed in the case of coronal plasma, where the anomalous resistivity prevails over the collisional one and can be responsible for the redistribution of large-scale currents in magnetic loops and small-scale currents in the regions of reconnection of magnetic-field lines.
Systematic discrepancies between electron temperature (Te) measurements from radially viewing electron cyclotron emission (ECE) and Thomson scattering (TS) diagnostics have been observed in multiple tokamaks and are widely attributed to non-Maxwellian features in the electron velocity distribution function (EVDF). As the International Thermonuclear Experimental Reactor (ITER) is expected to operate at much higher core temperatures than present devices, accurate Te measurements from ECE become increasingly critical, particularly in the presence of non-Maxwellian EVDFs. This work presents ECE spectra simulations performed to assess the diagnostic capability of the ITER oblique view under ITER conditions where non-Maxwellian distributions are present. The results show that above a critical oblique angle, Doppler broadening becomes the dominant effect, masking fine-scale spectral signatures of non-thermal distortions across a wide range of conditions. Furthermore, at higher EC harmonics in either polarization, the spectra remain largely unaffected by non-thermal emissions, enabling reliable reconstruction of Te-profiles. These findings demonstrate that the ITER oblique ECE view retains sufficient sensitivity to non-thermal electrons while providing robust and accurate Te-profile measurements under reactor-relevant conditions. The results presented here also have direct implications for ITER ECE system operation and channel selection, particularly in the presence of non-Maxwellian electron populations.
Tokamaks and stellarators are the leading magnetic-confinement concepts for fusion, but they rely on complementary design principles. Tokamaks use simple axisymmetric coils and plasma current, whereas stellarators use externally generated three-dimensional fields for steady-state operation. Here, we propose a programmable stellarator--tokamak hybrid that uses a fixed set of simple planar coils to access a broad magnetic-configuration space. The device adds 288 dipole-field coils to a tokamak-like coil set, with only six independent coil geometries required by symmetry. By programming coil currents, the same hardware generates more than 1.66 million optimized stellarator configurations spanning quasi-axisymmetry, quasi-helical symmetry, and quasi-isodynamicity, as well as tokamak-relevant three-dimensional perturbations. Representative configurations exhibit nested magnetic surfaces, low neoclassical transport, and favorable energetic-particle confinement. This approach enables rapid magnetic-configuration discovery without hardware redesign.
Low-pressure radio-frequency capacitively coupled plasmas operated in Ar/O$_2$ gas mixtures are widely adopted in critical semiconductor manufacturing processes. O($^3$P) and O($^1$D) are key highly reactive species for oxidation or as oxygen sources for deposited thin films. Optimizing external parameters to realize efficient generation of these species under limited energy deposition is essential for improving process yield.Based on a one-dimensional (1D) fluid/electron Monte Carlo (EMC) hybrid model, this study investigates the energy cost of O($^1$D) and O($^3$P) generation driven by sawtooth up-type voltage waveforms at a fixed peak-to-peak voltage, focusing on the effects of the harmonic number ($N$) and the O$_2$ ratio. The results show that O($^3$P) generation is consistently more efficient than that of O($^1$D). The generation energy cost decreases with increasing O$_2$ ratio, yet increases as $N$ increases. However, in the specific scenario of 10% O$_2$, an inflection point can be observed at $N = 2$. As $N$ increases from 1 to 2, the discharge mode shifts from the DA mode to the $α$-DA hybrid mode, expanding the effective spatio-temporal range of the ionization rate and enhancing its peak, which increases electron density. Consequently, the generation rates are significantly enhanced, leading to a reduction in the generation energy cost.Moreover, as discussed above, monotonically increasing the harmonic number $N$ does not reduce the generation energy cost of O($^1$D) and O($^3$P) associated with medium-energy (8-20 eV) electrons. Only by selecting the appropriate $N$ to sustain the discharge in the hybrid $α$-DA mode, thereby increasing the electron density and promoting the generation of these species, can the generation energy cost be reduced.
The idea of an energy cascade in the inertial range is often invoked in turbulent space plasmas to estimate the energy dissipation rate. Laws governing the behavior of third-order structure functions in the inertial range, so-called third-order laws, are among the few rigorous theoretical results quantifying cross-scale energy transfer. The widely used third-order-law derived rate assumes isotropy, which fundamentally conflicts with the anisotropic nature of space plasmas. Elementary questions persist regarding how such anisotropic energy cascades can be quantified using multi-spacecraft constellations. As the heliospheric community increasingly progresses towards multi-spacecraft, multi-scale constellations, such as Plasma Observatory and HelioSwarm, we revisit these crucial issues pertinent to accurately measuring the inertial-range energy transfer. Here we make a systematic comparison between two methods: direction-averaging (DA) and lag polyhedral derivative ensemble (LPDE) to determine the full three-dimensional (3D) dependence of cross-scale energy transfer. We find that DA exhibits both polar and azimuthal dependence, but is insensitive to spacecraft configuration. By contrast, LPDE is strongly affected by spacecraft separation and tetrahedral shape, while being comparatively insensitive to the sampling trajectory. Our findings have direct implications for current and future multi-spacecraft missions. Both DA and LPDE will provide crucial information on the nature of turbulence in space and astrophysics.
We present iGENE, a fully-differentiable TensorFlow implementation of the electromagnetic local nonlinear gyrokinetic model, which allows us to compute gradients of any simulation output with respect to any input via automatic differentiation. We show that even if the stochastic nature of turbulence prevents the exact evaluation of gradients of nonlinear quantities of interest, they can still be successfully used to perform outer-loop tasks, such as profile predictions. This work enables the integration of gyrokinetics into automated parameter optimization, uncertainty quantification, sensitivity analysis, and AI workflows.
Building on recent progress in the understanding of particle transport in magnetized plasmas, we derive a scaling law for the formation of nonthermal spectral tails in mildly and strongly magnetized turbulent environments. We validate this scaling using driven-turbulence particle-in-cell simulations that incorporate particle escape, allowing the system to reach a steady state. The simulation results show good agreement with our theoretical predictions. We then discuss the astrophysical implications of these findings, focusing on proton acceleration in the coronae of supermassive black holes and the resulting high-energy neutrino emission.
Pre-plasma conditions strongly influence laser-plasma interactions in the multi-petawatt (MPW) regime, increasing the need for reliable early-time plasma evolution diagnostics. Among available pre-plasma diagnostics, interferometry remains the most direct method for measuring the spatially resolved electron density of pre-formed plasmas. However, its implementation becomes increasingly challenging at MPW scale dueto steep density gradients, phase-recovery difficulties, strong electromagnetic pulses (EMP), debris accumulation, and high-repetition-rate operation. Compounding these technical challenges, many large-scale facilities lack permanent probe-line architecture and trained diagnostic support, reducing experimental reproducibility and consuming limited beamtime. Future MPW facilities should standardize probe-line architecture, adopt off frequency probing strategies, improve phase-recovery methods for non-symmetric plasmas, integrate emerging real-time analysis capabilities, and engineer diagnostic systems resilient to EMP and high-repetition-rate environments. These advances will enable the user community to reliably characterize pre-plasma formation and laser-plasma dynamics at next-generation MPW facilities.
Hydrogen and some of its derivatives (such as e-methanol, e-methane, and e-ammonia) are promising energy carriers that have the potential to replace conventional fuels, thereby eliminating their harmful environmental impacts. An innovative use of hydrogen as a zero-emission fuel is forming weakly ionized plasma by seeding the combustion products of hydrogen with a small amount of an alkali metal vapor (cesium or potassium). This formed plasma can be used as a working fluid in supersonic open-cycle magnetohydrodynamic (OCMHD) power generators. In these OCMHD generators, direct-current (DC) electricity is generated straightforwardly without rotary turbogenerators. In the current study, we quantitatively and qualitatively explore the levels of electric conductivity and the resultant volumetric electric output power density in a typical OCMHD supersonic channel, where thermal equilibrium plasma is accelerated at a Mach number of two (Mach 2) while being subject to a strong applied magnetic field (applied magnetic-field flux density) of five teslas (5 T), and a temperature of 2300 K (2026.85 °C). We varied the total pressure of the pre-ionization seeded gas mixture between 1/16 atm and 16 atm. We also varied the seed level between 0.0625% and 16% (pre-ionization mole fraction). We also varied the seed type between cesium and potassium. We also varied the oxidizer type between air (oxygen-nitrogen mixture, 21-79% by mole) and pure oxygen. Our results suggest that the ideal power density can reach exceptional levels beyond 1000 MW/m3 (or 1 kW/cm3) provided that the total absolute pressure can be reduced to about 0.1 atm only and cesium is used for seeding rather than potassium.
The centroid oscillation of an offset laser pulse propagating in a preformed plasma channel is investigated through theoretical analysis and three-dimensional particle-in-cell simulations. For non-relativistic laser pulses, the mode leakage of a finite channel and the temporal walk-off between the fundamental and high order modes of a finite-duration laser induce a decay in the laser centroid oscillation. An analytical model characterizing these decay mechanisms is derived and validated by simulations. For relativistic laser pulses, the slice-based centroid oscillation frequency develops an axial chirp due to relativistic channel modification and photon deceleration. This chirp leads to phase mixing across different axial slices of the pulse, resulting in a rapid damping of the overall centroid oscillation. Understanding this oscillation damping is crucial for mitigating electron beam pointing jitter and maintaining beam quality in high-energy, channel-guided laser wakefield accelerators.
Due to computational efficiency and numerical stability limitations, coronal simulations constrained by static magnetograms are typically performed first and then used to drive inner-heliosphere (IH) models. In this paper, we calculate the Sun-to-Earth coronal and wind evolutions using a single time-evolving MHD model, showing that implicit MHD models have the potential to meaningfully simplify and improve the overall Sun-to-Earth modelling pipeline. We extend the implicit time-evolving coronal MHD model COCONUT out to 1 AU, and utilise it to investigate solar coronal and wind evolutions around a solar maximum Carrington rotation (CR). We compare quasi-steady-state and time-evolving Sun-to-Earth simulations to evaluate the impact of the inner-boundary magnetic field evolution, which is neglected in steady-state simulations, on background plasma parameters. Comparisons with commonly used coupled Sun-to-Earth simulations are also conducted to further validate and assess the Sun-to-Earth model COCONUT. The results show that the time-evolving implicit MHD modelling approach yields noticeable differences compared to oversimplified steady-state simulations, and is efficient enough for practical applications. Modelling the solar corona and wind using a single MHD model simplifies the modelling pipeline and avoids uncertainties associated with coupling different coronal and IH models. The noticeable differences in the temporal evolution of plasma parameters at the L1 and L5 points highlight the need to use continuously evolving, synchronised magnetic field observations to improve global coronal and solar wind simulations, whereas the overall consistent evolutionary trend reveals the reliability of using L5 observations to forecast solar wind conditions near Earth about four days in advance.
Spectroscopy of laser-produced plasmas offers an avenue for real-time, standoff and non-preparatory sensing of rare-earth elements (REEs) within a mineralogical context with applications spanning exploration geology to ore body mapping to ore sorting. Demonstrations of laser-induced breakdown spectroscopy (LIBS) in rock samples have employed both atomic and molecular detection for REE sensors. In this work we evaluate a complementary technique of absorption spectroscopy, realized with dual-frequency combs. This approach provides multi-THz (nm) spectral coverage with simultaneous sub-GHz (pm) resolution. It can improve accuracy and line identification confidence in congested multi-species spectra, which makes it ideal for multi-species evaluations present within mineralogical samples. We analyze REE signatures in calibrated reference materials (CRMs) and a synthesized, REE-containing alloy for atomic, ionic and molecular (oxide) absorptions across three spectral windows. We identify lines from rare-earth and matrix elements, compare absorption line strengths and investigate their temporal evolution. For La I, Sm I and Ce I, preliminary limits of detection from 54-583 ppm are estimated for CRMs, using univariate analysis of selected transitions. Comparing the CRM signatures to those of REEs synthesized in a copper alloy, we observe that most all REE lines appear earlier and disappear faster in the CRM samples. We attribute these dynamics to matrix effects: Among other elements, the increased oxygen content in the CRM could favor molecular formation. For rock samples, observations will once again differ due to grain sizes and bonding mechanisms. Compared to LIBS, we can resolve individual REE and matrix lines with minimal spectral overlap. These proof-of-principle results form a foundation for further development of this laser-based method as a mining sensor.
Optical skyrmions have been widely explored in recent years. Among them, Stokes skyrmions require sophisticated wave engineering or photonic devices for their generation. We show that Stokes skyrmions can emerge spontaneously in second-order harmonic generation in underdense plasmas driven by optical vortices. The nonlinear response produces a structured frequency-doubled field whose polarization texture maps the Poincaré sphere. When plasma inhomogeneities are taken into account, the electron density gradient deforms the skyrmionic texture, enabling topological diagnosis of plasma density.
Reconnection-driven nanoflares are widely considered a leading mechanism for coronal-loop heating, but their direct fingerprints in the tenuous coronal plasma remain elusive. The recently discovered coronal nanojets offer a potential probe of reconnection dynamics, but their extreme collimation, directionality and multi-wavelength visibility are not fully understood. Here we present a 3D rMHD simulation that unprecedentedly reproduces the key properties of nanojets, offering a viable model to explain their nature. These results provide a unified picture in which nanojet morphology, dynamics and detectability are contingent on the thermodynamic environment of reconnection. Together, our results point to a cool origin of coronal nanojets, where cool and dense material permits narrow, multi-band jet signatures to emerge from reconnection.
Long-period cyclic reversals of the large-scale magnetic field are a prominent feature of the dynamo driven by the magnetorotational instability (MRI) in accretion disks, but their physical origin remains unclear. We develop a quasilinear theory (QLT) of the MRI dynamo where the electromotive force (emf) is computed from the linear eigenfunctions under the WKB approximation. The emf depends on the mean field $\mathbf{B}$ more generally than standard mean-field closures allow. In the unstratified case, the leading order contribution to the large-scale dynamo is the shear-current effect: the emf depends on the current $\mathbf{J}$ as $\pmb{\varepsilon} = \pmbβ\cdot\mathbf{J}$, with a tensor $\pmbβ(\mathbf{B},t)$ that oscillates with time $t$ and whose off-diagonal components generate the mean field. The oscillations arise from beats between the two branches of MRI eigenfrequencies. Since the beat frequency varies only weakly with wavenumber, the beats remain coherent and drive the long-period butterfly cycle seen in local shearing box simulations. We predict a dominant cycle period $\sim 30{\left(1+a^2\right)}^{1/2}\,t_{\rm orb}$, with $a$ the vertical-to-radial aspect ratio and $t_{\rm orb}$ the orbital period, and an amplitude scaling $\sim a^2$ before saturation at $a\gtrsim 5$. Both trends agree with zero-net-flux unstratified shearing box simulations with Athena++. A carrier-envelope analysis of the simulation spectra shows that the same interference mechanism extends beyond strict QLT, through higher-order linear combinations of the eigenfrequencies, with observed cycles arising from pairwise beats within this spectral network. These results identify coherent interference between nearly degenerate eigenfrequencies as a key mechanism behind large-scale cyclic dynamos, with implications for magnetic variability in protoplanetary disks, X-ray binaries, and AGNs.