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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

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

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD)

Here, we present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-MD classical and nuclear quantum dynamical simulations in which the energies and forces of the full system are approximated by the 2-, 3-, and 4-body terms of the MBE expansion and compare the average potential and the vibrational density of states (VDOS) with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the 4-body term converges with near identical behavior to the one from the full simulation. The 3-body makes a substantial contribution (~20%) to the energy, whereas the 4- body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (~2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm -1 ) at the 4-body level, and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a 4-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N 4 ,N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

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

Validation of Density Functionals for Adsorption Energies on Transition Metal Surfaces

The quantitative prediction of adsorption energies of radicals and molecules on surfaces is essential for the design and understanding of heterogeneous catalytic processes. Here, a recent paper by Wellendorff et al. collected an experimental database of 39 reaction energies involving adsorption energies on transition metal surfaces that can be used as benchmarks for testing quantum mechanical electronic structure methods, and we compared the experimental data to Kohn–Sham density functional calculations with six exchange–correlation functionals. In this paper, we rearranged the data into two categories: open-shell radical adsorption reactions and closed-shell molecular adsorption reactions. We recalculated the adsorption energies with PBE, and we also calculated them with three functionals, M06-L, GAM, and MN 15 -L, that were not studied in the Wellendorff et al. paper; then we compared our results to the benchmark data. Of the nine functionals that have been compared to the databases, we find that BEEF-vdW, GAM, and RPBE perform best for the open-shell radical adsorption reactions, and MN15-L performs best for the closed-shell molecular adsorption, followed by BEEF-vdW and M06-L.

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

Examining the Effect of Exchange-Correlation Approximations in First-Principles Dynamics Simulation of Interfacial Charge Transfer

We examine the extent to which the exchange-correlation (XC) approximation influences modeling interfacial charge transfer using fewest-switches surface hopping (FSSH) simulations within the single-particle description. A heterogeneous interface between a lithium ion and an extended boron-nitride sheet was considered here, being an extreme case in which wave function localization and energy level alignments are highly sensitive to the XC approximation. The PBE0 hybrid XC approximation yields nonadiabatic couplings (NACs) that are significantly smaller than the values obtained from the PBE-GGA approximation by an order of magnitude for localized electronic states. This difference between the two XC functionals for the calculated NACs was found to derive mainly from the wave function characteristics rather than from the lattice movement although first-principles molecular dynamics trajectories, along which NACs are obtained, differ noticeably between the two XC functionals. Using the NACs and single-particle energy level alignments at different levels of theory, FSSH simulations were performed to model the electron transfer dynamics at the interface. Here, the electron transfer time scale was found to vary as much as, but not more than, 1 order of magnitude. The time scale was found to be quite sensitive to both NACs and energy level alignments. While the order of magnitude consistency for the charge transfer rate is encouraging even for this rather extreme model of heterojunction interface, continued advancement in electronic structure methods is required for quantitatively accurate determination of the transfer rate.

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

Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems

Nonadiabatic Molecular Dynamics (NAMD) of excited states has been widely used in the simulation of photoinduced phenomena. However, the inability to treat bond breaking and forming processes with single-reference electronic structure methods limits their application in photochemistry for extended molecular systems. In this work, the extension of excited-state NAMD for open-shell systems is developed and implemented in the NEXMD software. We present the spin-unrestricted CIS and TDSCF formalism for the ground and excited states, analytical derivatives, and nonadiabatic derivative couplings for the respective potential energy surfaces. Next, this methodology is employed to study the photochemical reaction of three model molecules. The results demonstrate the advantage of the open-shell approach in modeling photochemical reactions, especially involving bond breaking processes. We find that the open-shell method lowers the reaction barrier at the bond-breaking limits resulting in larger calculated photochemical quantum yields compared to the respective closed-shell results. We also address problems related to spin contamination in the open-shell method, especially when molecular geometries are far from equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accurate Atomic Correlation and Total Energies for Correlation Consistent Effective Core Potentials

Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches, we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly exact total energy calculations for these ccECPs. The calculations cover H–Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n = 6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of ≈1–10 mHa for K–Zn atoms and ≈0.1–0.3 mHa for all other elements—within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high-accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.

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

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

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

On the HCN – HNC Energy Difference

The value for the HCN → HNC 0 K isomerization energy has been investigated by combining state-of-the-art electronic structure methods with the Active Thermochemical Tables (ATcT) approach. The directly computed energy difference between HCN and HNC at the HEAT-456QP level of theory is 5236 ± 50 cm–1. This is substantially lower (by ~470 cm–1 or ~1.3 kcal/mol) than the recently proposed high-level multireference configuration interaction value of 5705 ± 20 cm–1 of Barber et al. (Mon. Not. R. Astron. Soc.2014, 437, 1828–1835). The discrepancy was analyzed by the ATcT approach, using several distinct steps, which (a) independently corroborated the current single-reference HEAT-456QP result, (b) independently found that the recent multireference-based value is highly unlikely to be correct within its originally stated uncertainty, and (c) produced a recommended value of 5212 ± 30 cm–1 for the HCN → HNC isomerization energy at 0 K, based on all currently available knowledge. The ATcT standard enthalpies of formation at 0 and 298 K for HCN, HNC, and their cations and anions are also presented.

Active Thermochemical Tables↗
DOE OSTI2017

Accurate anharmonic zero-point energies for some combustion-related species from diffusion Monte Carlo

Full dimensional analytic potential energy surfaces based on CCSD(T)/cc-pVTZ calculations have been determined for 48 small combustion related molecules. The analytic surfaces have been used in Diffusion Monte Carlo calculations of the anharmonic, zero point energies. Here, the resulting anharmonicity corrections are compared to vibrational perturbation theory results based both on the same level of electronic structure theory and on lower level electronic structure methods (B3LYP and MP2).

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

Calculating Free Energy Changes in Continuum Solvation Models

We recently showed for a large dataset of pK a s and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents (Phys. Chem. Chem. Phys. 2015, 17, 2859). In this paper, we significantly expand the scope of our study to examine the suitability of this approach for the calculation of general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMDM062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems and reaction types. This includes cluster-continuum schemes for pK a calculations, as well as various neutral, radical and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and bimolecular SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar – the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 kJ mol -1 and 25 kJ mol -1 respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment. The results indicate that when appropriate combinations of electronic structure methods are employed, the direct approach provides a reliable alternative to the thermodynamic cycle calculations of solution phase kinetics and thermodynamics across a broad range of organic reactions.

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

Relevance of the Nuclear Quantum Effects on the Proton/Deuteron Transmission through Hexagonal Boron Nitride and Graphene Monolayers

According to recent experiments, atomically thin hexagonal boron nitride and graphene are permeable to protons and deuterons (and not to other atomic species), and the experimental estimates of the activation energy are lower than the theoretical values by about 0.5 eV for the isolated proton-membrane transfer model. Our analysis of the electronic potential energy surfaces along the normal to the transmission direction, obtained using correlated electronic structure methods, suggests that the aqueous environment is essential to stabilize the proton { as opposed to the hydrogenatom { transmission. Therefore, the process is examined within a molecular model of H 2 O { H(D) + { material { H 2 O. Exact quantum-mechanical scattering calculations are performed to assess the relevance of the nuclear quantum eects, such as tunneling factors and the kinetic isotope eect (KIE). Deuteration is found to aect the thermal reaction rate constants (KIE of 3-4 for hexagonal boron nitride and 20-30 for the graphene) and to eectively lower the barriers to the proton transfer by 0.2 and 0.4 eV for the two membranes, respectively. This lowering eect is reduced for the deuteron by approximately a factor of three. A more comprehensive description of the proton transmission is likely to require an extended explicit aqueous environment.

36 MATERIALS SCIENCE↗
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↗