Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Oppenheimer”

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

Attosecond response of molecules to impulsive ionization

When matter interacts with energetic radiation it can undergo sudden, or impulsive, ionization. This process can drive chemical change and occurs widely in space and planetary atmospheres, yet its comprehensive description challenges our current theoretical and computational capabilities as it requires advanced treatment of electron correlation and nonadiabatic dynamics beyond the Born–Oppenheimer approximation. Here, in this study, we measure the response of the para-aminophenol molecule to sudden ionization. Using attosecond X-ray absorption spectroscopy, we resolve the ultrafast dynamics of the ionized molecule with atomic precision. A subfemtosecond decay corresponds to states undergoing non-radiative decay, whereas few-femtosecond oscillatory signatures are associated with electronic wavepacket motion in stable cation states that later couple to nuclear motion. We compare our measurement with state-of-the-art computational modelling, qualitatively reproducing the observed response across multiple timescales. These results provide a benchmark for computational models of sudden ionization and ultrafast charge motion in matter.

Driver, Taran [SLAC National Accelerator Laborator↗
DOE OSTI2020

Real-time degradation dynamics of hydrated per- and polyfluoroalkyl substances (PFASs) in the presence of excess electrons

Per- and polyfluoroalkyl substances (PFASs) are synthetic chemicals that are harmful to both the environment and human health. In this work, using self-interaction-corrected Born–Oppenheimer molecular dynamics simulations, we provide the first real-time assessment of PFAS degradation in the presence of excess electrons. In particular, we show that the initial phase of the degradation involves the transformation of an alkane-type C–C bond into an alkene-type C[double bond, length as m-dash]C bond in the PFAS molecule, which is initiated by the trans elimination of fluorine atoms bonded to these adjacent carbon atoms.

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

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

We report nonadiabatic dynamics, which goes beyond the Born–Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

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

Polaron absorption in aligned conjugated polymer films: breakdown of adiabatic treatments and going beyond the conventional mid-gap state model

This study provides the first experimental polarized intermolecular and intramolecular optical absorption components of field-induced polarons in regioregular poly(3-hexylthiophene-2,5-diyl), rr-P3HT, a polymer semiconductor. Highly aligned rr-P3HT thin films were prepared by a high temperature shear-alignment process that orients polymer backbones along the shearing direction. rr-P3HT in-plane molecular orientation was measured by electron diffraction, and out-of-plane orientation was measured through series of synchrotron X-ray scattering techniques. Then, with molecular orientation quantified, polarized charge modulation spectroscopy was used to probe mid-IR polaron absorption in the ℏω = 0.075 — 0.75 eV range and unambiguously assign intermolecular and intramolecular optical absorption components of hole polarons in rr-P3HT. This data represents the first experimental quantification of these polarized components and allowed long-standing theoretical predictions to be compared to experimental results. The experimental data is discrepant with predictions of polaron absorption based on an adiabatic framework that works under the Born–Oppenheimer approximation, but the data is entirely consistent with a more recent nonadiabatic treatment of absorption based on a modified Holstein Hamiltonian. This nonadiabatic treatment was used to show that both intermolecular and intramolecular polaron coherence break down at length scales significantly smaller than estimated structural coherence in either direction. Importantly, this strongly suggests that polaron delocalization is fundamentally limited by energetic disorder in rr-P3HT.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

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

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

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

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. In this study, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

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

Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane

Ab Initio Multiple Spawning (AIMS) simulates the excited-state dynamics of molecular systems by representing nuclear wavepackets in a basis of coupled traveling Gaussian functions, called trajectory basis functions (TBFs). New TBFs are spawned when nuclear wavepackets enter regions of strong nonadiabaticity, permitting the description of non-Born–Oppenheimer processes. The spawning algorithm is simultaneously the blessing and the curse of the AIMS method: it allows for an accurate description of the transfer of nuclear amplitude between different electronic states, but it also dramatically increases the computational cost of the AIMS dynamics as all TBFs are coupled. Recently, a strategy coined stochastic-selection AIMS (SSAIMS) was devised to limit the ever-growing number of TBFs and tested on simple molecules. In this work, we use the photodynamics of three different molecules—cyclopropanone, fulvene, and 1,2-dithiane—to investigate (i) the potential of SSAIMS to reproduce reference AIMS results for challenging nonadiabatic dynamics, (ii) the compromise achieved by SSAIMS in obtaining accurate results while using the smallest average number of TBFs as possible, and (iii) the performance of SSAIMS in comparison to the mixed quantum/classical method trajectory surface hopping (TSH)—both in terms of its accuracy and computational cost. We show that SSAIMS can accurately reproduce the AIMS results for the three molecules considered at a much cheaper computational cost, often close to that of TSH. We deduce from these tests that an overlap-based criterion for the stochastic-selection process leads to the best agreement with the reference AIMS dynamics for the smallest average number of TBFs.

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

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

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

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Classical-quantum scattering

We analyze the framework recently proposed by Oppenheim et al (2023 Nat. Commun. 14; 2023 Phys. Rev. X 13 041040; arXiv:2302.07283 [gr-qc]; 2023 J. High Energy Phys. JHEP08(2023)163) to model relativistic quantum fields coupled to relativistic, classical, stochastic fields (in particular, as a model of quantum matter coupled to ‘classical gravity’). Perhaps surprisingly, we find that we can define and calculate scattering probabilities which are Lorentz-covariant and conserve total probability, at least at tree level. As a concrete example, we analyze 2→2 scattering of quantum matter mediated by a classical Yukawa field. Mapping this to a gravitational coupling in the non-relativistic limit, and assuming that we can treat large objects as point masses, we find that the simplest possible ‘classical-quantum’ gravity theory constructed this way gives predictions for 2→2 gravitational scattering which are inconsistent with simple observations of, e.g. spacecraft undergoing slingshot maneuvers. We comment on lessons learned for attempts to couple quantum matter to ‘non-quantum’ gravity, or more generally, for attempts to couple relativistic quantum and classical systems.

Carney, Daniel↗
DOE OSTI2026

Direct inference of nuclear equation-of-state parameters from gravitational-wave observations

The observation of neutron star mergers with gravitational waves (GWs) has provided a new method to constrain the dense-matter equation of state (EOS) and to better understand its nuclear physics. However, inferring nuclear microphysics from GW observations necessitates the sampling of EOS model parameters that serve as input for each EOS used during the GW data analysis. The sampling of the EOS parameters requires solving the Tolman–Oppenheimer–Volkoff (TOV) equations a large number of times—a process that slows down each likelihood evaluation in the analysis on the order of a few seconds. Here, we employ emulators for the TOV equations built using multilayer perceptron neural networks to enable direct inference of nuclear EOS parameters from GW strain data. Our emulators allow us to rapidly solve the TOV equations, taking in EOS parameters and outputting the associated tidal deformability of a neutron star in only a few tens of milliseconds. We implement these emulators in PyCBC to directly infer the EOS parameters using the event GW170817, providing posteriors on these parameters informed solely by GWs. We benchmark these runs against analyses performed using the full TOV solver and find that the emulators achieve speed ups of nearly two orders of magnitude, with negligible differences in the recovered posteriors. Additionally, we constrain the slope and curvature of the symmetry energy at the 90% upper credible interval to be $L$ sym ≲ 106 MeV and $K$ sym ≲ 26 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2021

Superconductivity and electron-phonon interaction in Sr x Bi 2 Se 3 under pressure

Pressure-induced superconductivity has been widely explored and observed in Bi 2 Se 3 -based topological materials to hunt for topological superconductors. Although their triggered superconductivity has a close connection to their pressure-induced structural phase transitions, the quest for the electron pairing mechanism of these superconducting semiconductors in both their initial rhombohedral and high-pressure phases remains unknown. In this work, we present a systematic study of the pressure effect on superconducting properties and lattice dynamics using a combination of electrical transport, Raman-scattering, and synchrotron x-ray diffraction measurements using diamond anvil cells. One key finding is our observation of a cooperative connection between the strength of the electron–phonon interaction (EPI) generated by optical branches and the pressure-tunable superconductivity in rhombohedral Sr x Bi 2 Se 3 crystal. The underlying suppression mechanism of the T c by pressure is ascribed to the weakening of the electrons' interaction with the optical phonon modes in the rhombohedral phase. In the intermediate monoclinic phase, the T c value underwent a sharp increase with carrier density accumulation accompanying the concurrent enhanced EPI. This is intuitively unusual since it is expected that the EPI shall be weakened by inducing more conducting carriers in a normal metal. In the tetragonal phase, the superconductivity is interpreted within BCS theory, since it is fully metallized and obeys the adiabatic Born–Oppenheimer approximation well. Our findings are important to fully understand unconventional superconductivity and the unusual pairing mechanism in the layered rhombohedral Bi 2 Se 3 -based superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dynamical instability of collapsed dark matter halos

A self-interacting dark matter halo can experience gravothermal collapse, resulting in a central core with an ultrahigh density. It can further contract and collapse into a black hole, a mechanism proposed to explain the origin of supermassive black holes. We study dynamical instability of the core in general relativity. Here we use a truncated Maxwell-Boltzmann distribution to model the dark matter distribution and solve the Tolman-Oppenheimer-Volkoff equation. For given model parameters, we obtain a series of equilibrium configurations and examine their dynamical instability based on considerations of total energy, binding energy, fractional binding energy, and adiabatic index. Our numerical results indicate that the core can collapse into a black hole when the fractional binding energy reaches 0.035 with a central gravitational redshift of 0.5. We further show for the instability to occur in the classical regime, the boundary temperature of the core should be at least 10% of the mass of dark matter particles; for a 10 9 M ⊙ seed black hole, the particle mass needs to be larger than a few keV. These results can be used to constrain different collapse models, in particular, those with dissipative dark matter interactions.

79 ASTRONOMY AND ASTROPHYSICS↗