Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure method”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Prediction of impact sensitivity, heat of formation and heat of explosion using atomic connectivity

In these proceedings we revisit a large collection of explosives and explosive descriptors with the goal of predicting impact sensitivity using only local atomic environments that can be deciphered from molecular SMILES strings as descriptors without utilizing empirically measured values or computationally expensive electronic structure calculations. From the original database of nearly 500 descriptors, removing empirically measured and electronic structure values decreased the number of descriptors to 135, which we reduced to 18 the most important descriptors using Random Forests. The condensed model predicted impact sensitivity with essentially the same accuracy as the existing, more complex model (R 2 = 0.788 and RSME = 0.312), while remaining applicable to all types of explosives (Peroxides, azides, C-Nitros, Nitroamines, Nitrate Esters, etc.). In addition to impact sensitivity, we proposed similar models to accurately predict values heat of formation (ΔH f ) and heat of explosion (Q), with R 2 = 0.966 and 0.916, respectively. In conclusion, the work in these proceedings allows for prediction of explosive performance and sensitivity with only chemical structure information and an estimate of density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

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

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

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

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

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

Excitons: Energetics and spatiotemporal dynamics

The concept of an exciton as a quasiparticle that represents collective excited states was originally adapted from solid-state physics and has been successfully applied to molecular aggregates by relying on the well-established limits of the Wannier exciton and the Frenkel exciton. However, the study of excitons in more complex chemical systems and solid materials over the past two decades has made it clear that simple concepts based on Wannier or Frenkel excitons are not sufficient to describe detailed excitonic behavior, especially in nano-structured solid materials, multichromophoric macromolecules, and complex molecular aggregates. In addition, important effects such as vibronic coupling, the influence of charge-transfer (CT) components, spin-state interconversion, and electronic correlation, which had long been studied but not fully understood, have turned out to play a central role in many systems. This has motivated new experimental approaches and theoretical studies of increasing sophistication. This article provides an overview of works addressing these issues that were published for A Special Topic of the Journal of Chemical Physics on ``Excitons: Energetics and spatio-temporal dynamics" and discusses their implications.

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

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

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

Mass-selected ion–molecule cluster beam apparatus for ultrafast photofragmentation studies

We describe an apparatus for investigating the excited-state dissociation dynamics of mass-selected ion–molecule clusters by mass-resolving and detecting photofragment-ions and neutrals, in coincidence, using an ultrafast laser operating at high repetition rates. The apparatus comprises a source that generates ion–molecule clusters, a time-of-flight spectrometer, and a mass filter that selects the desired anions, and a linear-plus-quadratic reflectron mass spectrometer that discriminates the fragment anions after the femtosecond laser excites the clusters. The fragment neutrals and anions are then captured by two channeltron detectors. The apparatus performance is tested by measuring the photofragments: I – , CF 3 I – , and neutrals from photoexcitation of the ion–molecule cluster CF 3 I·I – using femtosecond UV laser pulses with a wavelength of 266 nm. Here, the experimental results are compared with our ground state and excited state electronic structure calculations as well as the existing results and calculations, with particular attention to the generation mechanism of the anion fragments and dissociation channels of the ion–molecule cluster CF 3 I·I – in the charge-transfer excited state.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Hybrid programming-model strategies for GPU offloading of electronic structure calculation kernels

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. In this study, we demonstrate the portability of this approach and its performance in benchmark problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring direct photodetachment and photodissociation–photodetachment dynamics of platinum iodide anions (PtI n - , n = 2–5) using cryogenic photoelectron spectroscopy

The direct photodetachment and two-photon photodissociation–photodetachment processes of a series of PtI n - (n = 2–5) anions were systematically studied using cryogenic anion photoelectron spectroscopy and first-principles electronic structure calculations. The adiabatic/vertical detachment energies (ADEs/VDEs) of these anions were determined from their 193 nm photoelectron (PE) spectra, i.e., 3.54/3.63, 4.04/4.09, 4.33/4.36, and 4.37/4.41 eV for n = 2–5, respectively, and well reproduced by B3LYP-D3(BJ)/aug-cc-pVTZ-pp calculations. As the coordination number increases, the electron affinity (EA) of PtI n • (n = 2–5) neutrals (equivalent to the corresponding anion’s ADE) gradually increases, exceeding the EA of Cl at n = 3 and exhibiting superhalogen characteristics for n ≥ 3. Meanwhile, the ground state transition contributed from detaching electrons in the highest occupied molecular orbital gradually evolves from the central metal Pt to the iodine ligands. For the PtI 3 - anion, besides one-photon direct detachment, four distinct two-photon photodissociation–photodetachment channels were identified, and the competition between them was discussed.

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

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

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

Symmetry breaking as predicted by a phase space Hamiltonian with a spin Coriolis potential

Here, we perform electronic structure calculations for a set of molecules with degenerate spin-dependent ground states ( 3 CH 2 , 2 CH$^{•}_{3}$, 3 O 2 ) going beyond the Born–Oppenheimer approximation and accounting for nuclear motion. According to a phase space approach that parameterizes electronic states (|Φ⟩) and electronic energies (E) by nuclear position and momentum [i.e., |Φ(R, P)⟩ and E(R, P)], we find that the presence of degenerate spin degrees of freedom leads to broken symmetry ground states. More precisely, rather than a single degenerate minimum at (R, P) = (R min , 0), the ground state energy has two minima at (R,P)=(R' min ,±P min ) (where R' min is close to R min ), dramatically contradicting the notion that the total energy of the system can be written in separable form as E = $\frac{P^2}{2M}$ + V el . Although we find that the broken symmetry solutions have small barriers between them for the small molecules, we hypothesize that the barriers should be macroscopically large for metallic solids, thus offering up a new phase-space potential energy surface for simulating the Einstein–de Haas effect.

Berry connection↗
DOE OSTI2021

Variational quantum eigensolver for approximate diagonalization of downfolded Hamiltonians using generalized unitary coupled cluster ansatz

In this paper we discuss the utilization of Variational Quantum Solver (VQE) and recently introduced Generalized Unitary Coupled Cluster (GUCC) formalism for the diagonalization of downfolded/effective Hamiltonians in active spaces. In addition to effective Hamiltonians defined by the downfolding of a subset of virtual orbitals we also consider their form defined by freezing core orbitals, which enables us to deal with larger systems. Furthermore, we also consider various solvers to identify solutions of the GUCC equations. We use N 2 , H 2 O, and C 2 H 4 , and benchmark systems to illustrate the performance of the combined framework.

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

Correlated purification for restoring 𝑁-representability in quantum simulation

Experimentally measured reduced density matrices (RDMs) often violate constraints that ensure they represent N-electron states—known as N-representability conditions—because of statistical and hardware noise. In this work, we present a correlated purification framework based on semidefinite programming to restore the accuracy of a noisy, unphysical two-electron RDM (2-RDM). The method performs a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the correction to the measured 2-RDM. The nuclear norm, often employed in matrix completion, promotes low-rank corrections, while the energy term acts as a regularization term that can improve the purity of the ground state. While the method is particularly effective for ground states, it can also be applied to excited and nonstationary states by decreasing the weight of the energy relative to the error norm. In an application to fermionic shadow tomography of large hydrogen chains, correlated purification yields substantial reductions in both energy and 2-RDM error, achieving chemical accuracy across dissociation curves. This framework provides a robust strategy for tomography in many-body quantum simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gauge invariance of light-matter interactions in first-principle tight-binding models

In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.

2-dimensional systems↗
DOE OSTI2021

Twisted symmetric trilayer graphene. II. Projected Hartree-Fock study

The Hamiltonian of the magic-angle twisted symmetric trilayer graphene (TSTG) can be decomposed into a twisted-bilayer-graphene- (TBG-) like flat band Hamiltonian and a high-velocity Dirac fermion Hamiltonian. We use Hartree-Fock mean field approach to study the projected Coulomb interacting Hamiltonian of TSTG developed in Calugaru et al. [Phys. Rev. B 103, 195411 (2021)] at integer fillings ν = -3, -2, -1, and 0 measured from charge neutrality. We study the phase diagram with w 0 /w 1 , the ratio of AA and AB interlayer hoppings, and the displacement field, which introduces an interlayer potential U and hybridizes the TBG-like bands with the Dirac bands. At small U, we find the ground states at all fillings ν are in the same phases as the tensor products of a Dirac semimetal with the filling ν TBG insulator ground states, which are spin-valley polarized at ν = -3, and fully (partially) intervalley coherent at ν = -2, 0(ν = -1) in the flat bands. An exception is ν = -3 with w 0 /w 1 ≳ 0.7 , which possibly becomes a metal with competing orders at small U due to charge transfers between the Dirac and flat bands. At strong U where the bandwidths exceed interactions, all the fillings ν enter a metal phase with small or zero valley polarization and intervalley coherence. Lastly, at intermediate U, semimetal or insulator phases with zero intervalley coherence may arise for ν = -2, -1, 0. Finally, our results provide a simple picture for the electron interactions in TSTG systems, and reveal the connection between the TSTG and TBG ground states.

2-dimensional systems↗
DOE OSTI2022

Interaction-expansion inchworm Monte Carlo solver for lattice and impurity models

Multiorbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications, including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. We generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. So, the problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive time stepping for the two-time integro-differential Kadanoff-Baym equations

The nonequilibrium Green's function gives access to one-body observables for quantum systems. Of particular interest are quantities such as density, currents, and absorption spectra which are important for interpreting experimental results in quantum transport and spectroscopy. We present an integration scheme for the Green's function's equations of motion, the Kadanoff-Baym equations (KBE), which is both adaptive in the time integrator step size and method order as well as the history integration order. We analyze the importance of solving the KBE self-consistently and show that adapting the order of history integral evaluation is important for obtaining accurate results. To examine the efficiency of our method, we compare runtimes to a state-of-the-art fixed time step integrator for several test systems and show an order of magnitude speedup at similar levels of accuracy. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗