SEARCH · Search NASA

Results for “Local density approximations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗

DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗

DOE OSTI2025

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗

DOE OSTI2026

Unraveling Grain Boundary Instability in Dense Proton-Conducting Oxides

The long-term stability of protonic ceramic electrolysis cell (PCEC) materials under high-steam operating conditions remains a critical barrier to device commercialization. Here, we investigate the fundamental degradation mechanisms of dense BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3-δ (BCZYYb) electrolytes operated at 550 °C, 50% H 2 O in air. Over 1,000 h, the total electrolyte conductivity decreases by 11.1%, driven primarily by a >130% increase in grain-boundary resistivity. Post-mortem analyses reveal that damage is localized to near-surface grain boundaries extending ∼50 μm into the dense electrolyte pellet. This surface localization indicates that degradation is likely to be severe in thin, device-level electrolytes. Degradation is primarily attributed to chemo-mechanical grain-boundary weakening arising from hydration-induced chemical expansion, culminating in the formation of intergranular cracks oriented parallel to the pellet surface. These internal cracks subsequently react with steam and/or CO 2 , leading to the formation of nanoscale insulating phases, including Ba(OH) 2 , nanocrystalline BaCO 3 , and amorphous Ce/Zr/Y/Yb-containing oxides or hydroxycarbonates. After an initial degradation period of approximately 200 h, the overall conductivity stabilizes. Incorporating NiO sintering aids reduces grain-boundary density by an order of magnitude under identical sintering conditions. Although addition of NiO increases the initial resistivity by >160% at 550 °C, it substantially suppresses grain-boundary instability and mitigates chemical degradation. These findings underscore the urgent need for chemical and/or physical stabilization of BCZYYb electrolytes and offer design guidelines to enable durable, high-performance PCECs.

08 HYDROGEN↗

DOE OSTI2024

ELM-free H-mode phase and decoupling of peeling–ballooning stability boundary in the MAST Upgrade tokamak

A linear magnetohydrodynamic (MHD) peeling–ballooning stability analysis of the edge-localized mode (ELM)-free phase of a MAST Upgrade (MAST-U) H-mode plasma is presented. In contrast to other similar discharges, #47018 is found to have a significantly higher and wider pedestal during its ELM-free H-mode phase that lasts for approximately 80 ms; this is made possible by the reduced core MHD mode activity on the q = 2 surface. During this period, there is sustained decoupling of peeling and ballooning branches of the stability boundary on J–α space, opening an access channel to the second stability regime with higher peaks in pedestal current density $J_{\mathrm{N,ped}}$ and pressure gradient (α). Such decoupling of the stability boundary has not previously been observed in MAST-U H-modes, and if such a condition can be readily reproduced, it opens a wide range of opportunities for MAST-U to explore low-collisionality peeling-limited pedestal regimes as well as advanced scenarios such as quiescent H-modes that are relevant to future reactors such as STEP and ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Observation of Orbital-Selective Dual Modulations in an Anisotropic Antiferromagnetic Kagome Metal TbTi 3 ⁢Bi 4

Orbital selectivity is pivotal in dictating the phase diagrams of multiorbital systems, with prominent examples including the orbital-selective Mott phase and superconductivity. The intercalation of anisotropic layers represents an effective method for enhancing orbital selectivity and thereby shaping the low-energy physics of multiorbital systems. Despite its potential, related experimental studies, especially those elucidating the correlation between orbital selectivity and magnetism, remain limited. In this work, we systematically examine the interplay between orbital selectivity and magnetism in the newly discovered anisotropic kagome TbTi 3 ⁢Bi 4 single crystal, and report the coexistence of orbital-selective dual-band modulations (𝑞 1 ∼ 1/3⁢𝑎*, 𝑞 2 ∼ 0.28⁢𝑏*) within the antiferromagnetic (AFM) state. By combining soft x-ray and vacuum ultraviolet angle-resolved photoemission spectroscopy measurements, neutron powder diffraction, scanning tunneling microscopy, and density-functional-theory calculations, we identify these dual-band reconstructions as manifestations of the AFM order driven by a (approximately 1/3, 0.28, 0) nesting instability of the intercalated Tb 5⁢𝑑 𝑥⁢𝑧 orbitals. These orbital-selective modulations induce unusual momentum-dependent band folding and lead to the emergence of Dirac cones only at the $\bar{M}$ 1 point, signaling a topological phase transition in the AFM state. Importantly, the discovery of orbital-selective (approximately 1/3, 0.28, 0) AFM order offers crucial insights into the mechanism underlying the fractional magnetization plateau in this kagome AFM metal. Our findings not only underscore the essential role of both conducting and localized electrons in determining the magnetic orders of Ln⁢Ti 3 ⁢Bi 4 (Ln = lanthanide) kagome metals but also offer a pathway for manipulating magnetism through selective control of anisotropic electronic structures.

Zhang, Renjie [Shanghai Jiao Tong University (Chin↗

DOE OSTI2025

A charged particle transport approximation for thick and thin plasmas

In this work, we will construct a simple method to be utilized in the interpretation of the properties of a thin plasma from a study of the transport of charged particles through it. We do so by first demonstrating how charged particle fluxes encode local plasma information in their spectra for sufficiently thick plasmas and then propose how to extend that analysis to the thin plasma limit where the extensive geometry cannot be neglected through arguments from scale separation. We provide a numerical treatment of the transport problem and demonstrate that it is in good agreement with the approximation we present. Finally, we utilize both the numerical solution and our novel approximation to study the impact of the extensive scale of the plasma on the yield and flux normalized high-energy neutron spectra resulting from the upscattering of charged fuel ions. Using this analysis, we show that (after controlling for fusion yield) for the same uniform densities and temperatures, larger plasmas have higher magnitude but softer reaction-in-flight (RIF) neutron spectra relative to smaller plasmas. This is because larger plasmas retain more knocked-on suprathermal ions within their bulk and can downscatter them to lower average energies, while smaller plasmas allow a larger fraction of high energy particles to “range out” of the system prior to substantial downscattering or inducing an RIF reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Libration of hydroxyl groups in layered aluminum (oxy)hydroxides and other material analogs: insights from inelastic neutron scattering and theory

We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.

Wang, Hsiu-Wen [Oak Ridge National Laboratory (ORN↗

DOE OSTI2025

Cosmological quasiparticles and the cosmological collider

The interplay between cosmology and strongly coupled dynamics can yield transient spectral features that vanish at late times, but which may leave behind phenomenological signatures in the spectrum of primordial fluctuations. Of particular interest are strongly coupled extensions of the standard model featuring approximate conformal invariance. In flat space, the spectral density for a scalar operator in a conformal field theory is characterized by a continuum with scaling law governed by the dimension of the operator, and is otherwise featureless. Anti–de Sitter/conformal field theory arguments suggest that for large 𝑁, in an inflationary background with Hubble rate 𝐻, this continuum is gapped. We demonstrate that there can be additional peak structures that become sharp and particlelike at phenomenologically interesting regions in parameter space, and we estimate their contribution to cosmological observables. We find phenomena that are potentially observable in future experiments that are unique to these models, including displaced oscillatory features in the squeezed limit of the bispectrum. These particles can be either fundamental, and localized to a UV brane, or composite at the Hubble scale, 𝐻, and bound to a horizon in the bulk of the 5D geometry. We comment on how stabilization of conformal symmetry breaking vacua can be correlated with these spectral features and their phenomenology.

conformal field theory↗

DOE OSTI2025

How Many Bursts Does It Take to Form a Core at the Center of a Galaxy?

We present a novel method for systematically assessing the impact of the central potential fluctuations associated with bursty outflows on the structures of dark matter halos for classical and ultrafaint dwarf (UFD) galaxies. Specifically, we use dark-matter-only simulations augmented with a manually added massive particle that modifies the central potential and approximately accounts for a centrally concentrated baryonic component. This approach enables precise control over the magnitude, frequency, and timing of rapid outflow events. We demonstrate that this method can reproduce the established result of core formation for systems that undergo multiple episodes of bursty outflows. In contrast, we also find that equivalent models involving only single (or a small number of) burst episodes do not form cores with the same efficacy. This is important because many UFDs in the local Universe are observed to have tightly constrained star formation histories that are best described by a single early burst of star formation. Using a suite of cosmological zoom-in simulations, we identify the regimes in which single bursts can and cannot form a cored density profile. Our results suggest that it may be difficult to form cores in UFD-mass systems with a single early burst, regardless of its magnitude.

cold dark matter↗

DOE OSTI2025

Simulations of Quantum Approximate Optimization Algorithm on HPC-QC Integrated Systems

The Quantum Approximate Optimization Algorithm (QAOA) has emerged as a promising tool for accelerating optimization processes in the Noisy Intermediate-Scale Quantum (NISQ) era. Compared to classical methods, QAOA efficiently solves optimization problems, often formulated as Quadratic Unconstrained Binary Optimization (QUBO) problems. Classical quantum simulators are crucial for evaluating quantum algorithms due to limited quantum resources. However, QAOA's performance can vary with different simulation methods. This study analyzes QAOA's performance using various quantum simulators (e.g., density _matrix, statevector, and matrix_product_state) and demonstrates the benefits of HPC-QC integrated systems in solving QUBO problems on an active learning workflow. By simulating QAOA on dense, large-matrix QUBO problems, we evaluate accuracy and problem-solving time. We also assess QAOA's performance on local computers and HPC-QC inte-grated systems, using Oak Ridge Leadership Computing Facility (OLCF)'s Frontier supercomputer with local Qiskit Aer and remote IBM Quantum simulators.

Kim, Seongmin [ORNL] (ORCID:0000000159063004)↗

DOE OSTI2025

Ultrasonic characterization of material heterogeneities in stainless steel components produced by laser powder bed fusion

We introduce pulse-echo ultrasound as a method for characterizing the impact of powder bed fusion parameters on the properties of additively manufactured stainless-steel components, their material anisotropy, and location-dependent heterogeneity. Our results indicate that accurate characterization requires careful selection of ultrasonic propagation paths, which must consider the direction of additive layering, variations in processing parameters, and the component's geometry. We employed two distinct methods to estimate material properties from ultrasonic data: One assumes isotropy, while the other accounts for anisotropic interactions during the propagation of elastic waves. When applied to samples fabricated with laser energy densities ranging from 24 to 42 J/mm³ , these methods revealed transverse isotropy and weak anisotropy (quantified by small Thomsen parameters, ε = 0.0651 and γ = 0.0092) and less than a ∼6 % change in acoustic impedance. The assumption of isotropy, in this case, leads to small errors (less than 4 % or 1 % for Young's modulus in the build or transverse directions) when estimating orthotropic material properties using ultrasonic data measured along just two orthogonal directions, one of which must align with the build direction. By comparing ultrasonic measurements — which aggregate the spatial variability in material properties along the length of elastic wave propagation into a single value — with localized measurements obtained from surface nanoindentation, we uncovered and spatially profiled significant differences between the surface and interior properties. Specifically, the surface Young's modulus decreased from approximately 210 GPa to 180 GPa within a depth of about 3 mm. We attribute this surface-localized heterogeneity in PBF-fabricated components to distinct thermal histories experienced by the surface and interior regions. Collectively, the results of this study establish a framework for the ultrasonic characterization of material heterogeneity and anisotropy in material properties and demonstrate its application in additively manufactured metal components.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Electron transport barrier and high confinement in configurations with internal islands close to the plasma edge of W7-X

The low magnetic shear in the Wendelstein 7-X (W7-X) stellarator makes it feasible to shape the separatrix by the large islands constituting an island-divertor, and this can be exploited to access various magnetic configurations, including samples of different internal island sizes and locations. To investigate the configuration effects on the plasma confinement, a configuration scan was performed by changing the coil currents to vary the rotational transform between values 5/4 and 5/6 at the plasma boundary with different power levels (2, 4, 6 MW) of electron cyclotron resonance heating (ECRH) at a maximum plasma density of 8 × 10 19 m -2 . neutral beam injection (NBI) heating was also applied during some configurations of the scan to create a density ramp and access high densities beyond the X2 ECRH cutoff. For the magnetic configurations, where the 5/5 and 5/6 island chains were moved closer to separatrix but remaining inside the last closed flux surface, the electron cyclotron emission shows that an electron temperature, T e , pedestal develops already during ECRH heated plasma buildup phase indicating a transport barrier, and the barrier sustains irrespective of changed plasma heating conditions such as NBI in the later part of discharge. The transport barrier is broken by subsequent fast crashes, observed with multiple plasma diagnostics with characteristics such as tokamak edge localized modes, and the corresponding crash amplitude and frequency vary with plasma pressure. The impact of the transport barrier on plasma confinement can be seen through the increased core T e profile, which could be responsible for the overall increase in the stored diamagnetic energy by approximately 10% for these configurations. After the plasma heating is terminated, a backwards transition to a degraded confinement state is also observed. These observations indicate a configuration triggered high confinement mode in low shear W7-X. This work focuses on the occurrence of this transport barrier for different magnetic configurations and its relation to internal magnetic islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Mapping and Characterizing the Visual Impacts of the Existing US Wind Turbine Fleet

Visual impacts of wind turbines have been a persistent concern for wind energy development in the United States (US) for decades and remain a major source of project delays and cancellations. Assessments of visual impacts are frequently performed at a local scale for individual projects, but a comprehensive understanding of broader geographic patterns in visual impacts across the US is lacking. This paper presents a visual impact assessment of the existing land-based wind turbine fleet of the contiguous United States (CONUS). The assessment combines geographic information systems and 3D simulation methods to account for key factors driving the visual magnitude of impacts from the installed turbines. The results indicate that, despite the deployment of approximately 70,000 turbines and over 144 gigawatts of land-based wind in the CONUS, the visual impacts are very small when measured as a proportion of land area, population, and sensitive visual resources. Nonetheless, visual impacts are not evenly distributed: people experience a concentrated share in a small number of natural settings, primarily including plains, prairies, and deserts. Finally, we find that although increased density of wind development consistently leads to visual impacts across a greater proportion of land, it does not always lead to impacts to a greater share of the population. These findings suggest that visual impacts from wind energy are generally well-mitigated across the CONUS to date but also highlight the need for a deeper understanding of landscape sensitivity and individual perceptions of wind turbines in the most heavily impacted natural settings.

3D simulation↗

DOE OSTI2024

Hole doping and electronic correlations in Cr substituted BaFe$_{2}$As$_{2}$

For a significant composition range, the suppression of the spin density wave transition temperature (T SDW ) in Cr- and Mn-substituted BaFe 2 As 2 (CrBFA and MnBFA, respectively) coincides as a function of Cr/Mn content, despite the distinct electronic effects of these substitutions. Additionally, for any Cr/Mn content superconductivity (SC) is absent and this topic is particularly less explored in the case of CrBFA. Here, in this work, we employ angle-resolved photoemission spectroscopy (ARPES) and combined density functional theory plus dynamical mean field theory (DFT+DMFT) to address the evolution of the Fermi surface (FS) and electronic correlations in CrBFA. Our findings reveal that incorporating Cr leads to an effective hole doping of the states near the FS, which is well described within the virtual crystal approximation (VCA). Moreover, analysis of the ARPES spectra of the bands with main d yz -orbital character reveals a fractional scaling of the imaginary part of self-energy as a function of the binding energy, a signature property of Hund's correlations. Our DFT+DMFT calculations support these experimental findings. We conclude that CrBFA is a correlated electron system for which the changes in the FS as a function of Cr are unrelated to the suppression of T SDW . In addition, we suggest that the absence of SC is primarily due to the competition between Cr local moments and the Fe-derived itinerant spin fluctuations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗