Search NASA⌕ Search

SEARCH · Search NASA

Results for “electronic structure 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 OSTI2021

Quantum Machine-Learning for Eigenstate Filtration in Two-Dimensional Materials

Quantum machine-learning algorithms have emerged to be a promising alternative to their classical counterparts as they leverage the power of quantum computers. Such algorithms have been developed to solve problems like electronic structure calculations of molecular systems and spin models in magnetic systems. However, the discussion in all these recipes focuses specifically on targeting the ground state. Herein we demonstrate a quantum algorithm that can filter any energy eigenstate of the system based on either symmetry properties or a predefined choice of the user. The workhorse of our technique is a shallow neural network encoding the desired state of the system with the amplitude computed by sampling the Gibbs–Boltzmann distribution using a quantum circuit and the phase information obtained classically from the nonlinear activation of a separate set of neurons. We show that the resource requirements of our algorithm are strictly quadratic. To demonstrate its efficacy, we use state filtration in monolayer transition metal dichalcogenides which are hitherto unexplored in any flavor of quantum simulations. We implement our algorithm not only on quantum simulators but also on actual IBM-Q quantum devices and show good agreement with the results procured from conventional electronic structure calculations. We thus expect our protocol to provide a new alternative in exploring the band structures of exquisite materials to usual electronic structure methods or machine-learning techniques that are implementable solely on a classical computer.

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

Low Thermal Conductivity in Heteroanionic Materials with Layers of Homoleptic Polyhedra

Although BiAgOSe, an analogue of a well-studied thermoelectric material BiCuOSe, is thermodynamically stable, its synthesis is complicated by the low driving force of formation from the stable binary and ternary intermediates. Here we have developed a “subtraction strategy” to suppress byproducts and produce pure phase BiAgOSe using hydrothermal methods. Electronic structure calculations and optical characterization show that BiAgOSe is an indirect bandgap semiconductor with a bandgap of 0.95 eV. The prepared sample exhibits lower lattice thermal conductivities (0.61 W·m –1 ·K –1 at room temperature and 0.35 W·m –1 ·K –1 at 650 K) than BiCuOSe. Further, lattice dynamical simulations and variable temperature diffraction measurements demonstrate that the low lattice thermal conductivity arises from both the low sound velocity and high phonon–phonon scattering rates in BiAgOSe. These in turn result primarily from the soft Ag–Se bonds in the edge-sharing AgSe 4 tetrahedra and large sublattice mismatch between the quasi-two-dimensional [Bi 2 O 2 ] 2+ and [Ag 2 Se 2 ] 2– layers. These results highlight the advantages of manipulating the chemistry of homoleptic polyhedra in heteroanionic compounds for electronic structure and phonon transport control.

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

Modeling Methane Adsorption in Interpenetrating Porous Polymer Networks

Porous polymer networks (PPNs) are a class of porous materials of particular interest in a variety of energy-related applications because of their stability, high surface areas, and gas uptake capacities. Computationally derived structures for five recently synthesized PPN frameworks, PPN-2, -3, -4, -5, and -6, were generated for various topologies, optimized using semiempirical electronic structure methods, and evaluated using classical grand canonical Monte Carlo simulations. We show that a key factor in modeling the methane uptake performance of these materials is whether, and how, these material frameworks interpenetrate and demonstrate a computational approach for predicting the presence, degree, and nature of interpenetration in PPNs that enables the reproduction of experimental adsorption data.

Martin, RL↗
DOE OSTI2013

Dielectric Properties of Water Ice, the Ice Ih/XI Phase Transition, and an Assessment of Density Functional Theory

The dielectric properties of the hydrogen disordered hexagonal phase (Ih) of water ice have been computed using density functional theory (DFT) based Monte Carlo simulations in the isobaric–isothermal ensemble. Temperature dependent data yield a fit for the Curie–Weiss law of the system and hence a prediction of the temperature of the phase transition from the Ih phase to the hydrogen ordered ice XI phase. Direct simulations around the phase transition temperature confirm and refine the predicted phase transition temperatures and provide data for further properties, such as the linear thermal expansion coefficient. Results have been obtained with both hybrid and semilocal density functionals, which yields insight in the performance of the electronic structure method. In particular, the hybrid functional yields significantly more realistic dielectric constants than the semilocal variant, namely ε ≈ 116 as opposed to ε ≈ 151 at 273 K (ε experiment = 95). This can be attributed to the tendency of semilocal functionals to be biased to configurations with a large dipole moment, and their overestimation of the dipole moments of these configurations. This is also reflected in the estimates of the Ih/XI transition temperature, which is 70–80 and 90–100 K for the hybrid and semilocal functional respectively. Finally, DFT based sampling of the millions of configurations necessary for this work has been enabled by a Tree Monte Carlo algorithm, designed for massively parallel computers.

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

Electronic structure with direct diagonalization on a D-wave quantum annealer

Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Computational simulation and interpretation of the low-lying excited electronic states and electronic spectrum of thioanisole

Three singlet states, namely a closed-shell ground state and two excited states with 1 ππ* and 1 nσ* character, have been suggested to be responsible for the radiationless decay or photochemical reaction of photoexcited thioanisole. The correct interpretation of the electronic spectrum is critical for understanding the character of these low-lying excited states, but the experimental spectrum is yet to be fully interpreted. In the work reported here, we investigated the nature of those three states and a fourth singlet state of thioanisole using electronic structure calculations by multireference perturbation theory, by completely-renormalized equation-of-motion coupled cluster theory with single and double excitations and noniterative inclusion of connected triples (CR-EOM-CCSD(T)), and by linear-response time-dependent density functional theory (TDDFT). We clarified the assignment of the electronic spectrum by simulating it using a normal-mode sampling approach combined with TDDFT in the Tamm–Dancoff approximation (TDA). The understanding of the electronic states and of the accuracy of the electronic structure methods lays the foundation of our future work of constructing potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2026

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗
DOE OSTI2026

Strong field ionization and dissociation dynamics of vinyl bromide (C 2 H 3 Br) initiated by few-cycle pulses

Photoelectron-photoion (PEPICO) and photoion-photoion (PIPICO) coincidence measurements, coupled with a few-cycle intense laser field, were employed to investigate strong field ionization/dissociation dynamics of vinyl bromide (C₂H₃Br). The angular streaking technique was used to map the recoil-frame angle-dependent ionization rates to identify ionizing orbitals. Dissociative single ionization was mainly attributed to ionizing the lower lying molecular orbitals, suggesting that direct ionization dominates with few-cycle pulses. Three dissociative double ionization channels, both prompt and delayed fragmentation, were identified. Theoretical analysis using time-dependent configuration interaction with complex absorbing potential (TDCI-CAP) and high-level electronic structure methods was performed to elucidate the underlying ionization and dissociation mechanisms.

Olowolafe, Temitayo A. [Wayne State University, De↗
DOE OSTI2017

Investigation of the ozone formation reaction pathway: Comparisons of full configuration interaction quantum Monte Carlo and fixed-node diffusion Monte Carlo with contracted and uncontracted MRCI

The association/dissociation reaction path for ozone (O 2 + O ↔ O 3 ) is notoriously difficult to describe accurately using ab initio electronic structure theory, due to the importance of both strong and dynamic electron correlations. Experimentally, spectroscopic studies of the highest lying recorded vibrational states combined with the observed negative temperature dependence of the kinetics of oxygen isotope exchange reactions confirm that the reaction is barrierless, consistent with the latest potential energy surfaces. Previously reported potentials based on Davidson-corrected internally contracted multireference configuration interaction (MRCI) suffer from a spurious reef feature in the entrance channel even when extrapolated towards the complete basis set limit. In this work, we report an analysis of comparisons between a variety of electronic structure methods including internally contracted and uncontracted MRCI (with and without Davidson corrections), as well as full configuration interaction quantum Monte Carlo, fixed-node diffusion Monte Carlo, and density matrix renormalization group.

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

A theoretical and experimental benchmark study of core-excited states in nitrogen

The high resolution near edge X-ray absorption fine structure spectrum of nitrogen displays the vibrational structure of the core-excited states. This makes nitrogen well suited for assessing the accuracy of different electronic structure methods for core excitations. We report high resolution experimental measurements performed at the SOLEIL synchrotron facility. These are compared with theoretical spectra calculated using coupled cluster theory and algebraic diagrammatic construction theory. The coupled cluster singles and doubles with perturbative triples model known as CC3 is shown to accurately reproduce the experimental excitation energies as well as the spacing of the vibrational transitions. In conclusion, the computational results are also shown to be systematically improved within the coupled cluster hierarchy, with the coupled cluster singles, doubles, triples, and quadruples method faithfully reproducing the experimental vibrational structure.

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

Non-additive non-interacting kinetic energy of rare gas dimers

We show approximations of the non-additive non-interacting kinetic energy (NAKE) as an explicit functional of the density are the basis of several electronic structure methods that provide improved computational efficiency over standard Kohn-Sham calculations. However, within most fragment-based formalisms, there is no unique exact NAKE, making it difficult to develop general, robust approximations for it. When adjustments are made to the embedding formalisms to guarantee uniqueness, approximate functionals may be more meaningfully compared to the exact unique NAKE. We use numerically accurate inversions to study the exact NAKE of several rare-gas dimers within partition density functional theory, a method that provides the uniqueness for the exact NAKE. We find that the NAKE decreases nearly exponentially with atomic separation for the rare-gas dimers. We compute the logarithmic derivative of the NAKE with respect to the bond length for our numerically accurate inversions as well as for several approximate NAKE functionals. We show that standard approximate NAKE functionals do not reproduce the correct behavior for this logarithmic derivative and propose two new NAKE functionals that do. The first of these is based on a re-parametrization of a conjoint Perdew-Burke-Ernzerhof (PBE) functional. The second is a simple, physically motivated non-decomposable NAKE functional that matches the asymptotic decay constant without fitting.

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

On simulation of local fluxes in molecular junctions

We present a pedagogical review of the current density simulation in molecular junction models indicating its advantages and deficiencies in analysis of local junction transport characteristics. In particular, we argue that current density is a universal tool which provides more information than traditionally simulated bond currents, especially when discussing inelastic processes. However, current density simulations are sensitive to the choice of basis and electronic structure method. We note that while discussing the local current conservation in junctions, one has to account for the source term caused by the open character of the system and intra-molecular interactions. Our considerations are illustrated with numerical simulations of a benzenedithiol molecular junction.

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

Gaussian process based optimization of molecular geometries using statistically sampled energy surfaces from quantum Monte Carlo

Optimization of atomic coordinates and lattice parameters remains a significant challenge to the wide use of stochastic electronic structure methods such as quantum Monte Carlo (QMC). Measurements of the forces and stress tensor by these methods contain statistical errors, challenging conventional gradient-based numerical optimization methods that assume deterministic results. Additionally, forces are not yet available for some methods, wavefunctions, and basis sets and when available may be expensive to compute to sufficiently high statistical accuracy near energy minima, where the energy surfaces are flat. Here, we explore the use of Gaussian process based techniques to sample the energy surfaces and reduce sensitivity to the statistical nature of the problem. We utilize Latin hypercube sampling, with the number of sampled energy points scaling quadratically with the number of optimized parameters. Furthermore, we show these techniques may be successfully applied to systems consisting of tens of parameters, demonstrating QMC optimization of a benzene molecule starting from a randomly perturbed, broken symmetry geometry.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

High-accuracy extrapolated ab initio thermochemistry. IV. A modified recipe for computational efficiency

A number of economical modifications to the high-accuracy extrapolated ab initio thermochemistry (HEAT) model chemistry are evaluated. The two resulting schemes, designated as mHEAT and mHEAT+, are designed for efficient and pragmatic evaluation of molecular energies in systems somewhat larger than can be practically studied by the unapproximated HEAT scheme. It is found that mHEAT+ produces heats of formation with nearly subchemical ± 1 kJ/mol) accuracy at a substantially reduced cost relative to the full scheme. Total atomization energies calculated using the new thermochemical recipes are compared to the results of the HEAT-345(Q) model chemistry, and enthalpies of formation for the three protocols are also compared to Active Thermochemical Tables. Finally, a small selection of transition states is studied using mHEAT and mHEAT+, which illuminates some interesting features of reaction barriers and serves as an initial benchmark of the performance of these model chemistries for chemical kinetics applications.

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

Stochastic embedding DFT: Theory and application to p -nitroaniline in water

Over this past decade, we combined the idea of stochastic resolution of identity with a variety of electronic structure methods. In our stochastic Kohn-Sham density functional theory (DFT) method, the density is an average over multiple stochastic samples, with stochastic errors that decrease as the inverse square root of the number of sampling orbitals. Here in this paper, we develop a stochastic embedding density functional theory method (se-DFT) that selectively reduces the stochastic error (specifically on the forces) for a selected subsystem(s). The motivation, similar to that of other quantum embedding methods, is that for many systems of practical interest, the properties are often determined by only a small subsystem. In stochastic embedding DFT, two sets of orbitals are used: a deterministic one associated with the embedded subspace and the rest, which is described by a stochastic set. The method agrees exactly with deterministic calculations in the limit of a large number of stochastic samples. We apply se-DFT to study a p-nitroaniline molecule in water, where the statistical errors in the forces on the system (the p-nitroaniline molecule) are reduced by an order of magnitude compared with nonembedding stochastic DFT.

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

Interaction between water and carbon nanostructures: How good are current density functional approximations?

Due to their current and future technological applications, including realization of water filters and desalination membranes, water adsorption on graphitic sp2-bonded carbon is of overwhelming interest. However, these systems are notoriously challenging to model, even for electronic structure methods such as density functional theory (DFT), because of the crucial role played by London dispersion forces and noncovalent interactions, in general. Recent efforts have established reference quality interactions of several carbon nanostructures interacting with water. Here, we compile a new benchmark set (dubbed WaC18), which includes a single water molecule interacting with a broad range of carbon structures and various bulk (3D) and two-dimensional (2D) ice polymorphs. The performance of 28 approaches, including semilocal exchange-correlation functionals, nonlocal (Fock) exchange contributions, and long-range van der Waals (vdW) treatments, is tested by computing the deviations from the reference interaction energies. The calculated mean absolute deviations on the WaC18 set depend crucially on the DFT approach, ranging from 135 meV for local density approximation (LDA) to 12 meV for PBE0-D4. We find that modern vdW corrections to DFT significantly improve over their precursors. Within the 28 tested approaches, we identify the best performing within the functional classes of generalized gradient approximated (GGA), meta-GGA, vdW-DF, and hybrid DF, which are BLYP-D4, TPSS-D4, rev-vdW-DF2, and PBE0-D4, respectively.

Brandenburg, Jan Gerit (ORCID:0000000292192948)↗