Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “excited state methods”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

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

Nonequilibrium information entropy approach to ternary fission of actinides

Ternary fission of actinides probes the state of the nucleus at scission. Light clusters are produced in space and time very close to the scission point. Within the nonequilibrium statistical operator method, a generalized Gibbs distribution is constructed from the information given by the observed yields of isotopes. Using this relevant statistical operator, yields are calculated taking excited states and continuum correlations into account, in accordance with the virial expansion of the equation of state. Furthermore, clusters with mass number A≤10 are well described using the nonequilibrium generalizations of temperature and chemical potentials. Improving the virial expansion, in-medium effects may become of importance in determining the contribution of weakly bound states and continuum correlations to the intrinsic partition function. Yields of larger clusters, which fail to reach this quasiequilibrium form of the relevant distribution, are described by nucleation kinetics, and a saddle-to-scission relaxation time of about 7000 fm/c is inferred. Light-charged particle emission, described by reaction kinetics and virial expansions, may therefore be regarded as a very important tool to probe the nonequilibrium time evolution of actinide nuclei during fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

A massively parallel and scalable multi-CPU material point method

Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid (G2P2G) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array (AoSoA) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 10243 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Wang, Xinlei↗
DOE OSTI2020

Ab Initio Methods for L-edge X-ray Absorption Spectroscopy

The theoretical prediction of X-ray absorption spectra (XAS) has become common- place in electronic structure theory. The ability to better model and understand L-edge spectra is of great interest in the study of transition metal complexes and a wide variety of solid state materials. However, until recently few rst-principles works have mod- eled L-edge XAS due to the presence of strong spin-orbit coupling in the 2p orbitals which splits the observed peaks into multiple groups of features. Therefore, a proper description of spin-orbit coupling is vital for the successful prediction of L-edge spectra. A number of new approaches that incorporate spin-orbit coupling have recently made advances in the computation of L-edge spectra. In this review, we describe recent work in computational L-edge XAS and how these methods may continue to improve in the future. Comparison of the advantages and disadvantages of the various approaches are considered, with special attention to not only the computational cost of the level of theory, but also the various approaches that can be used to compute the absorption spectra with a large number of high energy excited states.

Kasper, Joseph M.↗
DOE OSTI2021

Hot-carrier dynamics in InAs/AlAsSb multiple-quantum wells

A type-II InAs/AlAs 0.16 Sb 0.84 multiple-quantum well sample is investigated for the photoexcited carrier dynamics as a function of excitation photon energy and lattice temperature. Time-resolved measurements are performed using a near-infrared pump pulse, with photon energies near to and above the band gap, probed with a terahertz probe pulse. The transient terahertz absorption is characterized by a multi-rise, multi-decay function that captures long-lived decay times and a metastable state for an excess-photon energy of >100 meV. For sufficient excess-photon energy, excitation of the metastable state is followed by a transition to the long-lived states. Excitation dependence of the long-lived states map onto a nearly-direct band gap (E g ) density of states with an Urbach tail below E g . As temperature increases, the long-lived decay times increase < E g , due to the increased phonon interaction of the unintentional defect states, and by phonon stabilization of the hot carriers > E g . Additionally, Auger (and/or trap-assisted Auger) scattering above the onset of the plateau may also contribute to longer hot-carrier lifetimes. Meanwhile, the initial decay component shows strong dependence on excitation energy and temperature, reflecting the complicated initial transfer of energy between valence-band and defect states, indicating methods to further prolong hot carriers for technological applications.

14 SOLAR ENERGY↗
DOE OSTI2023

Coherent detection of hidden spin–lattice coupling in a van der Waals antiferromagnet

Strong interactions between different degrees of freedom lead to exotic phases of matter with complex order parameters and emergent collective excitations. Conventional techniques, such as scattering and transport, probe the amplitudes of these excitations, but they are typically insensitive to phase. Therefore, novel methods with phase sensitivity are required to understand ground states with phase modulations and interactions that couple to the phase of collective modes. Here, by performing phase-resolved coherent phonon spectroscopy (CPS), we reveal a hidden spin–lattice coupling in a vdW antiferromagnet FePS 3 that eluded other phase-insensitive conventional probes, such as Raman and X-ray scattering. With comparative analysis and analytical calculations, we directly show that the magnetic order in FePS 3 selectively couples to the trigonal distortions through partially filled t 2 g orbitals. This magnetoelastic coupling is linear in magnetic order and lattice parameters, rendering these distortions inaccessible to inelastic scattering techniques. Our results not only capture the elusive spin–lattice coupling in FePS 3 but also establish phase-resolved CPS as a tool to investigate hidden interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Dodecane Radiolysis Yields by Time-Resolved and Steady-State Methods

Liquid organic molecules are present as solvents, complexing ligands, and additives in both used nuclear fuel reprocessing solvent systems and in their subsequent nuclear waste streams. Under these extreme environments, these organic molecules are constantly exposed to ionizing radiation which promotes their radiolysis, forming a variety of short-lived, highly energetic, excited state and radical species.1-4 Here, we demonstrate new experimental results for the steady-state and time-resolved irradiations of dodecane (C12H26), a long chain, liquid, aliphatic hydrocarbon that is the prototypical solvent used for benchtop studies of aqueous-organic solvent extraction systems. When ionizing radiation interacts with neat dodecane, the energy transfer can result in molecular ionization, to give the dodecane radical cation (C12H26+•) and the solvated electron (eS–), and molecular electronic excitation (C12H26*), which rapidly produces transient carbon-centered radical fragments (CxHy•) and hydrogen atoms (H•).1-4 Studies on the initial yields of the ionization and excitation products were performed using time-resolved picosecond electron pulse radiolysis with the use of molecular probes. Using steady-state cobalt-60 gamma irradiations, the suite of products formed by dodecane radiolysis in aerated and deaerated solutions was determined. Then, using iodine as an alkyl radical scavenger, the loss of molecular iodine with dose was quantified, and by correlating with the molecular hydrogen yields of the system, the initial yields of the various carbon-centered radicals were also determined. Finally, the rates of reactions of the C12H26+• and eS– with ligands proposed for use in spent nuclear fuel reprocessing were studied as a function of temperature from 10 – 40 °C.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2023

Analytic gradients for relativistic exact-two-component equation-of-motion coupled-cluster singles and doubles method

A first implementation of analytic gradients for spinor-based relativistic equation-of-motion coupled-cluster singles and doubles method using an exact two-component Hamiltonian augmented with atomic mean-field spin–orbit integrals is reported. To demonstrate its applicability, we present calculations of equilibrium structures and harmonic vibrational frequencies for the electronic ground and excited states of the radium mono-amide molecule (RaNH 2 ) and the radium mono-methoxide molecule (RaOCH 3 ). Spin–orbit coupling is shown to quench Jahn–Teller effects in the first excited state of RaOCH 3 , resulting in a C 3v equilibrium structure. Furthermore, the calculations also show that the radium atoms in these molecules serve as efficient optical cycling centers.

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

Decay spectroscopy of 160 Eu: Quasiparticle configurations of excited states and structure of K π = 4 + bandheads in 160 Gd

Background: Detailed spectroscopy of neutron-rich, heavy, deformed nuclei is of broad interest for nuclear astrophysics and nuclear structure. Nuclei in the r-process path and following freeze-out region impact the resulting r-process abundance distribution, and the structure of nuclei midshell in both proton and neutron number helps to understand the evolution of subshell gaps and large deformation in these nuclei. Purpose: We aim to improve the understanding of the nuclear structure of 160 Gd, specifically the K π = 4 + bands, as well as study the β decay of 160 Eu into 160 Gd. Methods: High-statistics decay spectroscopy of 160 Gd resulting from the β-decay of 160 Eu was collected using the GRIFFIN spectrometer at the TRIUMF-ISAC facility. Results: Two new excited states and ten new transitions were observed in 160 Gd. The β-decaying half-lives of the low- and high-spin isomers in 160 Eu were determined, and the low-spin state's half-life was measured to be t 1/2 = 26.0 (8) s, ≈ 16% shorter than previous measurements. Lifetimes of the two K π = 4 + bandheads in 160 Gd were measured for the first time, as well as γ – γ angular correlations and mixing ratios of intense transitions out of those bandheads. Conclusions: Lifetimes and mixing ratios suggest that the hexadecapole phonon model of the K π = 4 + bandheads in 160 Gd is preferred over a simple two-state strong mixing scenario, although further theoretical calculations are needed to fully understand these states. Additionally, the 1999.0-keV state in 160 Gd heavily populated in β decay is shown to have positive parity, which raises questions regarding the structure of the high-spin β-decaying state in 160 Eu.

150 ≤ A ≤ 189↗
DOE OSTI2021

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2024

Thermonuclear 28 P(p, γ ) 29 S reaction rate and astrophysical implication in ONe nova explosion

An accurate 28 P(p, γ) 29 S reaction rate is crucial to defining the nucleosynthesis products of explosive hydrogen burning in ONe novae. Using the recently released nuclear mass of 29 S, together with a shell model and a direct capture calculation, we reanalyzed the 28 P(p, γ) 29 S thermonuclear reaction rate and its astrophysical implication. We focus on improving the astrophysical rate for 28 P(p, γ) 29 S based on the newest nuclear mass data. Our goal is to explore the impact of the new rate and associated uncertainties on the nova nucleosynthesis. We evaluated this reaction rate via the sum of the isolated resonance contribution instead of the previously used Hauser-Feshbach statistical model. The corresponding rate uncertainty at different energies was derived using a Monte Carlo method. Nova nucleosynthesis is computed with the 1D hydrodynamic code SHIVA. The contribution from the capture on the first excited state at 105.64 keV in 28 P is taken into account for the first time. We find that the capture rate on the first excited state in 28 P is up to more than 12 times larger than the ground-state capture rate in the temperature region of 2.5 × 10 7 K to 4 × 10 8 K, resulting in the total 28 P(p, γ) 29 S reaction rate being enhanced by a factor of up to 1.4 at ~1 × 10 9 K. In addition, the rate uncertainty has been quantified for the first time. It is found that the new rate is smaller than the previous statistical model rates, but it still agrees with them within uncertainties for nova temperatures. The statistical model appears to be roughly valid for the rate estimation of this reaction in the nova nucleosynthesis scenario. Using the 1D hydrodynamic code SHIVA, we performed the nucleosynthesis calculations in a nova explosion to investigate the impact of the new rates of 28 P(p, γ) 29 S. Our calculations show that the nova abundance pattern is only marginally affected if we use our new rates with respect to the same simulations but statistical model rates. Finally, the isotopes whose abundance is most influenced by the present 28 P(p, γ) 29 S uncertainty are 28 Si, 33,34 S, 35,37 Cl, and 36 Ar, with relative abundance changes at the level of only 3% to 4%.

Astronomy & Astrophysics↗
DOE OSTI2025

Quantum computation of mass gap in an asymptotically free theory

In relativistic field theories, the mass spectrum is given by the difference between the energy of the vacuum and the excited states. Near the continuum limit, the cancellation between these two values leads to loss of precision. We propose a method to extract the mass gap directly using quantum computers and apply it to a particular version of the nonlinear $σ$-model with the correct continuum limit and perform calculations in quantum hardware (at strong coupling) and simulation in classical computers (at weak coupling).

Bedaque, Paulo F. [Maryland U.] (ORCID:00000001521↗
DOE OSTI2025

Strongly Correlated States of Transition Metal Spin Defects: The Case of an Iron Impurity in Aluminum Nitride

We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. Here, we consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.

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

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗