Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELECTRON ORBIT”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Forecasting Megaelectron-Volt Electrons Inside Earth's Outer Radiation Belt: PreMevE 2.0 Based on Supervised Machine Learning Algorithms

We present the recent progress in upgrading a predictive model for megaelectron-volt (MeV) electrons inside the Earth's outer Van Allen belt. This updated model, called PreMevE 2.0, provides improved forecasts, particularly at outer L-shells, by adding upstream solar wind speeds to the model's input parameter list that originally includes precipitating electrons observed at low Earth orbits and MeV electron fluxes in situ measured by a geosynchronous satellite. Furthermore, based on several kinds of linear and artificial neural networks algorithms, a list of models was constructed, trained, validated, and tested with 42-month MeV electron observations from Van Allen Probes. Out-of-sample test results from these models show that, with optimized model hyperparameters and input parameter combinations, the top performer from each category of models has the similar capability of making reliable 1-day (2-day) forecasts of 1-MeV electron flux distributions with performance efficiency values ~0.87 (~0.82) averaged over the L-shell range of 2.8–6.6, significantly outperforming the previous version of PreMevE particularly at L-shells > ~4.5. Interestingly, the linear regression model is often the most successful when compared to other models, which suggests the relationship between dynamics of trapped 1-MeV electrons and precipitating electrons is dominated by linear components. Results also show that PreMevE 2.0 can reasonably well predict the onsets of MeV electron events in 2-day forecasts. PreMevE 2.0 is designed to be driven by observations from longstanding space infrastructure to make high-fidelity forecasts for MeV electrons, thus can be an invaluable space weather forecasting tool for the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗
DOE OSTI2022

Comparison of Energetic Electron Fluxes Measured by GPS and THEMIS Spacecraft in the Inner Magnetosphere

Key elements of space weather models are energetic electron fluxes in the inner magnetosphere and the outer radiation belt. Flux depletion is driven by various loss processes: scattering into atmosphere, magnetopause shadowing. Flux enhancement is driven by various acceleration processes: local wave-particle interactions, radial transport, plasma sheet injections. Many of these processes operate on ~ hour timescales. Such mesoscale flux variations are not well traced by equatorial spacecraft with much longer orbits. Energetic electron detectors onboard the Global Positioning System (GPS) constellation provide a unique opportunity for probing such ~ hour-scale flux variations. Measurements from up to 23 identically instrumented GPS satellites cover a wide energy and L-shell range with a subhour time resolution. However, their orbits are inclined and thus all measurements at L-shell >4.3 are off-equatorial. Here, in this report, we present a comparison of equatorial THEMIS and nonequatorial GPS measurements of omnidirectional ≤600 keV electron fluxes. Such a comparison allows us to derive coefficients for using off-equatorial GPS fluxes to infer the equatorial values. These coefficients depend on particle energy and L-shell. We demonstrate a new data set derived from GPS measurements and discuss how it can be used to investigate mesoscale dynamics of energetic electron fluxes in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

ORBIT CORRECTION STUDIES ON THE ELECTRON TRANSPORT LINE FROM RCS TO ESR

A dedicated electron transfer line from Rapid Cycling Synchrotron (RCS) to Electron Storage Ring (ESR), referred to as the RTE line has been designed for the Electron-Ion Collider (EIC). The beamline follows a straight-line geometry, with a length of 133 m, and is consists with two matching sections and a FODO section for beam diagnostics. Imperfections with magnet alignments introduce orbit distortions, making orbit correction scheme a critical component in the design. To facilitate orbit correction, each quadrupole magnet is equipped with a pair of beam position monitors (BPMs) and kickers. The Singular Value Decomposition (SVD) algorithm is used for orbit correction and tolerance studies. This paper presents the ongoing progress in the optics design and error correction scheme of the RTE line.

Neththikumara, I. [Thomas Jefferson National Accel↗
DOE OSTI2024

First-principles investigation of elastic, vibrational, and thermodynamic properties of kagome metals CsM 3 Te 5 (M = Ti, Zr, Hf)

Kagome metals are a unique class of quantum materials characterized by their distinct atomic lattice arrangement, featuring interlocking triangles and expansive hexagonal voids. These lattice structures impart exotic properties, including superconductivity, interaction-driven topological many-body phenomena, and magnetism, among others. The kagome metal CsM 3 ⁢Te 5 (where M = Ti, Zr, or Hf) exhibits both superconductivity and nontrivial topological electronic properties, offering a promising platform for exploring topological superconductivity. This study employs first-principles density functional theory calculations to systematically analyze the elastic, mechanical, vibrational, thermodynamic, and electronic properties of CsM 3 ⁢Te 5 (M = Ti, Zr, Hf). Our calculations reveal that the studied compounds—CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 —are ductile metals with elastic properties akin to the hexagonal Bi and Sb, with average elastic constants, including a bulk modulus of 27 GPa, a shear modulus of 11 GPa, and Young's modulus of 29 GPa. We observe peculiar dispersionless, flat, phonon branches in the vibrational spectra of these metals. Additionally, we thoroughly analyze the symmetries of the zone-center phonon eigenvectors and predict vibrational fingerprints of the Raman- and infrared-active phonon modes. The analysis of thermodynamic properties reveals the Einstein temperature for CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 to be 66, 54, and 53 K, respectively. Our orbital-decomposed electronic structure calculations reveal significant in-plane steric interactions and multiple Dirac band crossings near the Fermi level. We further investigate the role of spin-orbit coupling effect on the studied properties. Furthermore, this theoretical investigation sheds light on the intriguing quantum behavior of kagome metals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

Effects of the surface termination and the oxygen vacancy position on LaNiO 3 ultra-thin films: First-principles study

In this study, while ultra-thin layers of the LaNiO 3 film exhibit a remarkable metal-insulator transition as the film thickness becomes smaller than a few unit cell (u.c.), the formation of possible oxygen vacancies and their effects on the correlated electronic structure have been rarely studied using first principles. Here, we investigate the effects of the surface termination and the oxygen vacancy position on the electronic properties and vacancy energetics of LaNiO 3 ultra-thin films under the compressive strain using density functional theory plus U (DFT + U). We find that oxygen vacancies can be easily formed in the Ni layers with the NiO 2 terminated surface (0.5 u.c. and 1.5 u.c. thickness) compared to the structures with the LaO terminated surface and the in-plane vacancy is energetically favored than the out-of-plane vacancy. When two vacancy sites are allowed, the Ni square plane geometry is energetically more stable in most cases as two oxygen vacancies tend to stay near a Ni ion. Strong anisotropy between the in-plane and out-of-plane vacancy formation as well as the layer and orbital dependent electronic structure occur due to strain, surface termination, charge reconstruction, and quantum confinement effects. The in-plane vacancy of the NiO 2 terminated structure is favored since the released charge due to the oxygen vacancy can be easily accommodated in the d x 2 -y 2 orbital, which is less occupied than the d z 2 orbital. Remarkably, the oxygen vacancy structure containing the Ni square-plane geometry becomes an insulating state in DFT + U with a sizable band gap of 1.2 eV because the large crystal field splitting between d z 2 and d x 2 -y 2 orbitals in the square-plane favors an insulating state and the Mott insulating state is induced in other Ni sites due to strong electronic correlations. In the thin-film structure without oxygen vacancies, the strong correlation effect in DFT + U drives a pseudogap ground state at the Fermi energy, similarly as the experimental photo-emission spectra; however, the variation of the DFT + U electronic structure depending on the surface termination becomes weaker compared to those obtained in DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Role of orbital selectivity on crystal structures and electronic states in BiMnO 3 and LaMnO 3 perovskites

Correlated oxides, such as BiMnO 3 and LaMnO 3 , show a complex interplay of electronic correlations and crystal structure exhibiting multiple first-order phase transitions, some without a clear order parameter. The quantitative theoretical description of this temperature-dependent electronic-structural interplay in the vicinity of a Mott transition is still a challenge. Here we address this issue by simultaneously considering both structural and electronic degrees of freedom, within a self-consistent density functional theory with embedded dynamical mean field theory. Our results show the existence of novel electronic states characterized by the coexistence of insulating, semimetallic, and metallic orbitals. This state is in agreement with resonant x-ray scattering. Here we also show that electronic entropy plays a decisive role in both electronic and structural phase transitions. By self-consistent determination of both the electronic state and the corresponding crystal structure, we show that the temperature evolution of these phases can be quantitatively explained from first principles, thus demonstrating the predictive power of the theoretical method for both the structural and the electronic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Electron-Ion Collider - A machine that will unlock the secrets of the strongest force in nature!

The computers and smartphones we use every day depend on what we learned about the atom in the last century. All information technology – and much of our economy today – relies on understanding the electromagnetic force between the atomic nucleus and the electrons that orbit it. The science of that force is well understood, but we still know little about the microcosm within the protons and neutrons that make up the atomic nucleus. That’s where Brookhaven National Laboratory (BNL) comes in. Brookhaven National Laboratory (located in Suffolk County, NY, about 60 miles east of midtown Manhattan) was recently chosen as the building site for an Electron-Ion Collider (EIC), a one-of-a-kind nuclear physics research facility. The EIC will be a discovery machine for unlocking the secrets of the “glue” that binds the building blocks of visible matter in the universe. The machine design will take advantage of the existing and highly optimized Relativistic Heavy Ion Collider (RHIC) that’s been operating at Brookhaven Lab since 2000. Beyond sparking scientific discoveries in a new frontier of fundamental physics, the Electron-Ion Collider will trigger technological breakthroughs that have broad-ranging impact on human health and national challenges.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

EIC ESR beam position monitors thermal and thermomechanical study preliminary design

The Electron Ion Collider (EIC) Electron Storage Ring (ESR) Beam Position Monitors (BPM) are designed to monitor the electron beam orbit position on a turn-by-turn and bunch-by-bunch basis in the vacuum chamber. Large variations in ESR fill pattern and bunch intensity imposes a large dynamic range on the BPM system design. And given the large ESR beam current - up to 2.5 A with 1160 x 20 nC bunches – the heating by synchrotron radiation (SR) and beam impedance is a significant concern. The high energy of the electrons – up to 18 GeV – leads to scat tering of the SR that produces high radiation doses and as sociated heating of structures away from the SR impact area. Because of this intense heating, and stringent position stability requirements, an in-depth engineering analysis was carried out. This paper will review the analysis made of the preliminary design and its main outcomes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Hund’s coupling mediated multi-channel quantum phase transition of a single magnetic impurity in Fe(Se, Te)

Understanding the interplay between individual magnetic impurities and superconductivity is crucial for bottom-up construction of novel phases of matter. Sub-gap bound states that are used in this endeavor are typically considered as independent entities that each result from the exchange scattering between the respective impurity orbitals and electrons of the superconducting condensate. Here we present experimental evidence of individual multi-spin impurities where the sub-gap states are not independent. Specifically, we find that by tuning the energy of the state closest to zero through zero, all other sub-gap states change particle-hole asymmetry as well. We show that this can be understood by including Hund’s coupling, which favors high-spin configurations, into a multi-orbital Anderson model. Unlike for the case of independent spins, the transition we observe signals the simultaneous departure of more than one quasiparticle from the impurity, while the parity of the ground state may remain unchanged. Our results show that Hund’s coupling is not only crucial in generating high-spin impurities, but also to understand the transition between two distinct ground states, and should therefore be taken into account for e.g. impurity-based band-structure engineering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

CoSn-type NiIn1–xSbx (0 ≤ x ≤ 0.17): Site-Selective Substitution, Electronic Structure, Chemical Bonding, and Structural Transformation

CoSn-type intermetallic compounds have emerged as a model platform for Kagome-derived flat-band physics, where subtle chemical perturbations can strongly influence electronic structure and phase stability. Here, we present a combined experimental and theoretical study of Sb-substitution in CoSn-type NiIn1–xSbx (0 ≤ x ≤ 0.17) to elucidate the interplay between site selectivity, solubility limit, chemical bonding, and electronic structure. Rietveld refinements on Neutron powder diffraction data confirmed the selective Sb-substitution at the electron-rich In2 (2d) site forming the honeycomb substructure, while the In1 (1a) site within the Kagome layer remains exclusively occupied by In. Density functional theory (DFT) calculations revealed that pristine CoSn-type NiIn hosts Ni 3d-dominated flat bands near the Fermi level (EF), originating from the Kagome-like Ni substructure. Partial replacement of In by Sb within the honeycomb layer alters these flat-band features below EF, reducing the density of states and suppressing the flat-band topology near the Fermi level. Orbital-resolved electronic structure and chemical-bonding analyses show that Sb-substitution enhances Ni-p-block (In/Sb) covalency and optimizes charge compensation, stabilizing the CoSn-type structure up to the solubility limit x ≈ 0.17. Beyond the limit, the higher-Sb compositions show satellite reflections consistent with an incommensurately modulated phase. These results establish a link between site-selective chemical substitution, bonding optimization, and flat-band electronic structure evolution, providing fundamental insights into how chemical substitution influences the electronic properties of Kagome-based intermetallic compounds.

Roy, Nilanjan [National Institute of Technology Si↗
DOE OSTI2022

Disorder driven variations in magnetoresistance and planar Hall effect in Bi 2 Te 3 thin films

The degenerately doped sesquichalcogenide Bi 2 Te 3 , a well-known thermoelectric material, has attracted much attention in recent years as a topological insulator in which conducting surface states lead to exotic electronic transport. Here, in this study, we report the observation of anisotropic magnetoresistance (AMR) and planar Hall effect (PHE) in polycrystalline thin films of Bi 2 Te 3 over a broad range of temperature and magnetic field strength. The values of AMR, PHE and the Hall mobility (μ H ) are compared with results of similar measurements on single crystals and epitaxial films of Bi 2 Te 3 , where these properties are understood in the framework of electronic transport in the topologically protected conducting surface states. Our textured polycrystalline films show a robust metallic character with a carrier density (n s ) and mobility μ H of 9×10 20 cm -3 and 16.7 cm 2 V -1 s -1 respectively. This carrier density is higher and the mobility lower by an order of magnitude compared to the reported values of n s and μ H in epitaxial films and single crystals. However, we see a comparable AMR and PHE in our films for which the concept of conducting surface states is ill-defined due to their polycrystallinity and disordered structure. We suggest that an anisotropic backscattering of charge carriers in the strongly spin – orbit coupled electronic states of this heavy-metal-based non-magnetic compound may lead to the measured resistivity anisotropy. The observation of AMR and PHE in polycrystalline films may be important for scalable production of Bi 2 Te 3 based AMR and PHE sensors of magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

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

Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation

A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗