Search NASA⌕ Search

SEARCH · Search NASA

Results for “self-consistent field”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Accurate energies for ππ* excited states via exchange scaling: the XS-CASSCF method

The state-averaged complete-active space self-consistent field method (SA-CASSCF) is a widely employed electronic structure method used for studying photochemistry and dynamics owing to its ability to provide a reliable description even of complicated cases while still retaining computational efficiency. However, SA-CASSCF suffers from one Achilles heel, related to the description of ionic ππ* excited states, whose energy is often overestimated by 1–2 eV. In light of this challenge, we present the XS-CASSCF method, a new approach based on the idea of exchange scaling (XS) that screens the involved energy terms to improve the excitation energies of singlet ionic ππ* states. First, we illustrate the power of the XS-CASSCF method using hexatriene and para-quinodimethane as examples, showing that it corrects the targeted ionic states while leaving the other states largely unaffected, giving root-mean-square errors (RMSE) below 0.2 eV for the four lowest states in both cases. Subsequently, XS-CASSCF vertical excitation energies are tested against theoretical best estimates for a set of 11 molecules and 56 excited states. XS-CASSCF performs exceptionally well for the ππ* states of hydrocarbons, reducing the RMSE over 21 excitation energies from 0.96 to 0.27 eV. In the challenging subset of molecules with heteroatoms and a larger number of ππ* and nπ* states, we find that improvements can also be obtained, albeit not as pronounced. We conclude with an outlook into more realistic molecular materials focusing on their singlet–triplet (S 1 /T 1 ) gaps, finding that significant improvements can be obtained along the whole range of S 1 /T 1 gaps studied, going from 0.1 eV to more than 1.5 eV. Owing to notable improvements across significant classes of molecules combined with its conceptual simplicity, we believe that XS-CASSCF is a promising addition to the electronic structure toolbox, serving both as a standalone electronic structure method and as a starting point for further correlated treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2026

Radiation-induced segregation in dilute Fe–Cr: A rate-theory framework for the Cr enrichment–depletion transition at the grain boundary

Radiation-induced segregation (RIS) poses a significant challenge for ferritic Fe–Cr alloys under irradiation, as it can compromise mechanical integrity and increase susceptibility to intergranular corrosion. Yet, the mechanisms governing Cr segregation remain incompletely understood. Here, in this study, we present a physics-based rate-theory model parameterized using self-consistent mean field theory-based Onsager transport coefficients to investigate RIS at the grain boundary in dilute Fe-(0.1 at. %) Cr. Under equal production rates of vacancies and self-interstitial atoms (SIA), and their equal absorption rates by bulk dislocations, the model simulates the experimentally observed transition from Cr enrichment at low temperatures to depletion at higher temperatures. Under these unbiased conditions, systematic investigation reveals that while temperature-dependent transport properties dictate the segregation direction, dose rate, grain size, and dislocation density only influence the magnitude and spatial extent of Cr segregation. However, under more realistic conditions of preferential vacancy production within damage cascade and/or preferential SIA absorption by bulk dislocations, the enrichment-to-depletion transition shifts to lower temperatures. Our findings demonstrate that RIS predictions based solely on transport coefficients are valid only under symmetric point defect flux conditions, and that biases in defect production and absorption must be considered for accurate predictions. This work provides a mechanistic framework for understanding RIS in ferritic alloys and informs alloy design for advanced nuclear systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Structure of domain walls in chiral spin liquids

The chiral spin liquid is a canonical state of quantum spins combining topological and symmetry-breaking order, and possible experimental realizations have attracted growing interest. We examine the physics at interfaces between chiral spin liquid domains of opposite chirality. We show that a self-consistent mean-field description of spinons remains possible in the vicinity of a domain wall and use this to formulate a Ginzburg–Landau theory of the domain wall. The bulk of a chiral spin liquid contains gapped spinon excitations and gauge fluctuations, set by a finite spinon mass and a nonzero spinon Chern number. A third class of excitations consists of amplitude fluctuations of the spinon hoppings, which admit a geometric interpretation in terms of effective vielbein fields. These fluctuations are usually neglected because they are irrelevant for a homogeneous chiral spin liquid and are suppressed in standard large- N treatments. Going beyond the purely topological Chern–Simons limit, we incorporate these fluctuations into an effective field theoretic framework and show that they generate momentum-dependent corrections, including Chern–Simons-like linear terms and higher-order contributions, beyond the universal topological limit. We then analyze nontopological properties, including domain wall tension and edge velocity, and explain how they modify observables relative to the uniform case. These results connect measurable, nonuniversal quantities such as domain wall width, domain wall tension, and edge velocity to microscopic parameters and provide concrete targets for experiments.

domain walls↗
DOE OSTI2026

Chiral superconductivity from a parent Chern band and its non-Abelian generalization

Here, we propose a minimal model starting from a parent Chern band with quartic dispersion that can describe the spin-valley polarized electrons in rhombohedral tetralayer graphene. The interplay between repulsive and attractive interactions on top of that parent Chern band is studied. We conduct standard self-consistent mean-field calculations, and find a rich phase diagram that consists of metal, quantum anomalous Hall crystal, chiral topological superconductor, as well as trivial gapped Bose-Einstein condensate. In particular, there exists a topological phase transition from the chiral superconductor to the Bose-Einstein condensate at zero temperature. Motivated by the recent experimental and theoretical studies of composite Fermi liquid in rhombohedral stacked multilayer graphene, we further generalize the physical electron model to its composite fermion counterpart based on a field theory analysis. The chiral superconductor phase of the composite fermion becomes the non-abelian Moore-Read quantum Hall phase. We argue that a chiral (pseudo-)spin liquid phase can emerge in the vicinity of this Moore-Read quantum Hall phase. Our work suggests rhombohedral multilayer graphene as a potential platform for rich correlated topological phases.

Wang, Yan-Qi [University of Maryland, College Park↗
DOE OSTI2026

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2026

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Instant folded strings, dark energy and a cyclic bouncing universe

We present a wholly self-consistent, complete cyclic bouncing cosmology based on components drawn from string theory and constructed in a way that is under perturbative control throughout (e.g., with temperature much less than the string scale and string coupling g s ≪ 1 at all times). The cyclic evolution is governed by standard dilaton-gravity in (3+1)-dimensions with a perturbatively generated potential and a coupling between the dilaton and a second field that becomes massless at ϕ = ϕ ESP , resulting in an enhanced symmetry point (ESP) that prevents the dilaton from running all the way to zero coupling. A central role is played by instant folded strings (IFSs) — fundamental strings with the unusual property of being much lighter than the string mass while extending far beyond the string length, and violating the Null Energy Condition (NEC). IFSs are produced classically when the string coupling grows with time, which occurs at two critical points in each cycle. In turn, they fulfill a dual function: Enabling cosmological bounces and initiating transient epochs of dark-energy domination that naturally transition into slow contraction. The resulting cosmology eliminates the cosmic singularity and multiverse problems of big bang inflationary models and robustly predicts time-varying IFS-induced dark energy and the absence of primordial B-mode polarization in the cosmic microwave background.

alternatives to inflation↗
DOE OSTI2026

Cyclotron breaking: a mechanism for parallel ion cyclotron waves to heat the fast solar wind

The Parker Solar Probe mission has observed near-continuous power in parallel ion cyclotron waves (PICWs) in the young, fast solar wind. These waves are unlikely to be directly produced by the turbulent cascade and are likely born of a local instability; yet, they are observed to both cool – and heat – the plasma. We propose that these observations can be self-consistently explained as the natural consequence of PICWs propagating in the inhomogeneous solar wind after they have been driven unstable. In this work, we argue that strong proton heating by a turbulent cascade of oblique ICWs will result in PICWs being driven unstable in a process known as quasi-linear focusing. Because the power in the turbulent cascade is concentrated at scales above the turbulent transition region, PICWs will be driven unstable within a range of wavenumbers parallel to the background magnetic field, 𝑘 ∥ , that is bounded from above by 𝑘$^{∗}_{∥P}$, corresponding to the start of the transition region. As unstable PICWs propagate away from the Sun to regions of lower proton density, their 𝑘 ∥ , multiplied by the proton inertial length 𝑑 p , increases. Eventually, 𝑘$^{∗}_{∥P}$ of the PICWs becomes larger than 𝑘$^{∗}_{∥P}$⁢𝑑 p and the waves damp, heating the solar wind. We call this effect ‘cyclotron breaking’, in analogy with ocean waves breaking on the shore. We then discuss the testable predictions of the theory, including a distinct heating signature in which PICWs cool fast protons and heat slow protons at any given heliocentric distance 𝑟. Finally, we conjecture that cyclotron breaking can lead to net heating by PICWs if the power emitted as PICWs decreases sufficiently rapidly with 𝑟 that local emission of PICWs is overwhelmed by the local damping of PICWs generated closer to the Sun.

plasma heating↗
DOE OSTI2026

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗
DOE OSTI2026

Muon-induced fission of actinide nuclei

A negative muon captured by an actinide cascades down through the muonic atomic levels; deeply bound transitions can proceed via inverse internal conversion, depositing the muonic transition energy directly into the nucleus and, when the deposited energy exceeds the fission barrier, inducing prompt fission. Because the muon mean lifetime exceeds the saddle-to-scission timescale by orders of magnitude, the muon can survive the entire fission event as a 1⁢𝑠 spectator and ultimately attach to one or both of the emerging fragments. Its postscission attachment probability to the light fragment, 𝑃 𝐿 , can be used as a direct electromagnetic probe of fission dynamics on a timescale of 10 −21 s. In previous work, we introduced a three-dimensional lattice solution of the time-dependent Dirac equation coupled to the electromagnetic field generated by a fissioning nucleus and reported 𝑃 𝐿 for several actinides at a single dissipation strength. In this work, we extend that framework to a systematic survey of 232 Th , 238 U , and 240 Pu and implement a more realistic fission model which incorporates dynamic pairing correlations. We find that 𝑃 𝐿 falls steeply with the fragment charge asymmetry, a robust structural fingerprint of the fissioning system, while its dependence on nuclear dissipation is secondary and sensitive to the phenomenological friction prescription. These results establish 𝑃 𝐿 as a clean electromagnetic probe of fragment charge asymmetry and motivate a self-consistent, coordinate- and time-dependent treatment of nuclear dissipation as the natural next step.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗