Search NASA⌕ Search

SEARCH · Search NASA

Results for “ORBIT PERTURBATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Relativistic Exact Two-Component Theory in the Generalized Pseudospectral Representation

We present a formulation and implementation of exact two-component (X2C) relativistic theory in the generalized pseudospectral representation. When combined with the Hartree-Fock-Slater framework, this approach enables efficient and accurate treatments of scalar-relativistic and spin-orbit effects in atomic electronic structures without the computational overhead of four-component methods. Benchmark calculations across light and heavy elements demonstrate that the our X2C scheme yields substantially more accurate relativistic corrections to core-electron binding energies than perturbational Breit-Pauli treatments while converging more rapidly with respect to basis size. The method provides an improved computational framework for modeling ultrafast x-ray-induced processes in heavy-element systems.

Wang, Xubo↗
DOE OSTI2026

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2026

A chemical bonding based descriptor for predicting the role of anharmonicity induced by quantum nuclear effects in hydride superconductors

Quantum nuclear effects (QNEs) can significantly alter a material’s crystal structure and phonon spectra, impacting properties such as thermal conductivity and superconductivity. However, predicting a priori whether these effects will enhance or suppress superconductivity, or destabilize a structure, remains a grand challenge. Herein, we address this unresolved problem by introducing two possible descriptors, based upon the integrated crystal orbital bonding index (iCOBI) or the bond valence function, to predict the influence of QNEs on a crystal lattice’s dynamic stability, phonon spectra and superconducting properties. We find that structures with atoms in symmetric chemical bonding environments exhibit greater resilience to structural perturbations induced by QNEs, while those with atoms in asymmetric bonding environments are more susceptible to structural alterations, resulting in enhanced superconducting critical temperatures.

Physical chemistry↗
DOE OSTI2024

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗
DOE OSTI2026

Collapse of Jahn-Teller phonons in La 1−x Sr x MnO 3 with weak magnetoresistance

Perovskite manganites are quantum materials exhibiting competing interactions inducing colossal magnetoresistance (CMR). The prevailing theory of CMR highlights the essential role of electron-phonon coupling (EPC), but mounting evidence suggests the underlying mechanism is more complicated. Here, we investigate phonons and spin-phonon coupling in ferromagnetic CMR manganites La 1−x Sr x MnO 3 (x=0.2,0.3) with relatively small CMR associated with melting of the magnetic order above room temperature. High-resolution neutron scattering experiments combined with density functional theory (DFT) show that the low-temperature ferromagnetic phase is conventional: neutron scattering from phonons agrees with DFT predictions and magnons follow sinusoidal dispersions. Fluctuating magnetic moments and low-energy phonons remain conventional in the high-temperature paramagnetic phase, indicating the Mn and La/Sr sublattices are not strongly perturbed by melting of ferromagnetism. In contrast, the Jahn–Teller-active optical oxygen vibrations collapse entirely above the Curie temperature, despite low CMR in these compositions, with some of the lost spectral weight reappearing as quasielastic scattering. We attribute this highly anomalous behavior to giant EPC in the charge and/or orbital channel. It drives cooperative diffusive motion of quasistatic carrier-trapping oxygen sublattice distortions once ferromagnetism disappears. We hypothesize the magnitude of magnetoresistance correlates with the rate of diffusion rather than with the strength of Jahn–Teller EPC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗
DOE OSTI2025

Effects of chemical disorder and spin-orbit coupling on the electronic-structure and Fermi-surface topology of YbSb-based monopnictides

Here, in this work, we study the influence of disorder on the electronic structure of YbSb—a rare-earth monopnictide featuring a simple rocksalt (B1) crystal structure and a well-defined Fermi-surface topology—by employing first-principles density-functional theory. We focus on chemical disorder introduced through Te and Al doping, selected based on their thermodynamic stability in alloyed configurations, to understand how such perturbations modify the electronic states of YbSb. Our results indicate that Te doping predominantly introduces electronlike states at the 𝑋 and 𝐿 points, while Al doping leads to a suppression of holelike states at Γ, effectively driving the system from a semimetallic state to one characterized by very narrow-gap behavior at Γ. This modulation of the Fermi surface, particularly the reduction of central hole pockets at Γ, plays a central role in altering interpocket scattering—a mechanism critical for tuning quantum transport properties, including superconductivity. This disorder-driven modulation of the Fermi surface, particularly the suppression of central hole pockets at Γ, controls interpocket scattering, which is essential for optimizing quantum transport properties, including superconductivity. Our results show that disorder can be effectively used for engineering band topology, thereby tuning quantum related response through a tailored electronic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

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

Near-half-metallic state in the half-Heusler PtMnSb film on a III-V substrate

The interplay between half-metallic ferromagnetism and spin-orbit coupling within the inversion-symmetry-broken structure of half Heuslers provides an ideal platform for various spintronics functionalities. Taking advantage of good lattice matching, it is highly desired to epitaxially integrate promising Heuslers into III-V semiconductor-based devices. PtMnSb is one of the first half Heuslers predicted to be an above-room-temperature half-metal with large spin orbit coupling, however, its half-metallicity and potential as a spintronics material has remained elusive due to lack of high quality samples. Here, in this study, we demonstrate epitaxial growth of single crystal PtMnSb(001) film on GaSb(001) substrates using molecular beam epitaxy. Direct observation of the band structure via angle-resolved photoemission spectroscopy and many-body perturbation theory within the quasi-particle self-consistent GW approximation (QPGW) reveal that PtMnSb hosts rather a near-half-metallic state with both spin bands crossing the Fermi level and with high spin polarization over 90%. Temperature dependence of magnetization also shows an anomalous enhancement below 60 K, which can be associated with the development of such a near-half-metallic state at low temperatures. Epitaxial growth of high crystalline PtMnSb on a III-V paves the way for systematic clarification of its spin transport properties with fine-tuning of strain in heterostructure devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Search for 3D topological superconductors using laser-based spectroscopy (Final Technical Report)

The three-dimensional topological superconductor (3D TSC) is a novel quantum phase of matter that is predicted to exhibit exotic thermal and electrical properties, which may serve applications ranging from precision thermal management and sensing technologies to fault-tolerant quantum information processing. Over the past decade, there has been an intensive effort to realize 1D and 2D TSCs by artificially engineering them from combinations of already existing materials. For example, coupling conventional superconductors to spin-orbit coupled 1D nanowires or to 2D topological insulator surfaces by the proximity effect provides a route to realizing 1D and 2D TSCs respectively. However no artificial engineering strategy is available for a 3D TSC. Therefore, current searches are limited to intrinsically superconducting compounds whose topological properties are difficult to predict. The overarching goal of this project was to realize and to identify 3D TSCs in bulk single crystals using novel laser-based spectroscopic techniques. One major thrust was to develop ultralow temperature and ultra-high-resolution angle-resolved photoemission spectroscopy (ARPES) to directly measure the dispersion of Majorana excitations that are predicted to exist on the surfaces of 3D TSCs. Another major thrust was to develop a rational route to realize 3D TSCs by identifying their precursor phases. It is theoretically proposed that the critical fluctuations of certain ordered electronic phases can mediate Cooper pairing in odd-parity channels and lead to topological superconductivity. Therefore, 3D TSCs may potentially emerge upon suppressing these inversion symmetry broken precursor phases to a critical point with external perturbations such as pressure, strain or even light. To test this hypothesis, we studied candidate materials using a suite of symmetry-sensitive optical probes including rotational anisotropy optical second harmonic generation (RA-SHG) and time-resolved coherent phonon spectroscopy, all operable under high-pressure and low temperature environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗
DOE OSTI2025

Bridging hadronic and vacuum structure by heavy quarkonia

We discuss the central and, mostly, spin-dependent potentials in heavy quarkonia $\overline{b}$$b$,$\overline{c}$$c$, with two goals in mind. The first is phenomenological: using the splitting between the 1S and 2S pairs, as well as the 1P and 2P quartet masses, we obtain pertinent matrix elements of all these potentials. The second is theoretical: using standard wave functions, we compute these matrix elements, quantifying the perturbative, nonperturbative and “string” contributions. The model for nonperturbative effects is a “dense instanton liquid model,” in which the QCD vacuum is made mostly of (strongly color correlated) instanton-anti-instanton pairs or “molecules.” Their contributions are calculated via standard Wilson lines, with or without extra powers of gauge fields. We find a reasonable description of all central, spin-spin, and spin-orbit forces at distances 𝑟 = 0–0.7 fm, relevant for $\overline{b}$$b$ and $\overline{c}$$c$ quarkonia.

Bound states↗
DOE OSTI2024

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗
DOE OSTI2024

Transport and losses of energetic particles in tokamaks in the presence of Alfvén activity using the new full orbit TAPaS code coupled to FAR3d

Recent developments and tools integrated into the TAPaS code are presented, enabling realistic scenario simulations of particle dynamics within experimental tokamak magnetic equilibria. In particular, the enhanced capabilities of TAPaS enable seamless coupling with external simulations, provided the metric and equilibrium magnetic field of the external code are known. Coupling TAPaS with the gyro-fluid code FAR3d, the transport and losses of energetic particles in the presence Alfvén eigenmodes (AEs) in DIII-D plasma discharge #159243 were investigated. Detailed analyses of prompt losses with and without collisions were performed. Then, further analysis was performed in the presence of electromagnetic perturbations resulting from AEs activity. The results indicate that, for the energies and the initial conditions considered here, the presence of AEs enhances the particle losses.

Alfven eigenmodes↗