Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure theory”

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 73 records · Page 4

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

Li 5 VF 4 (SO 4 ) 2 : A Prototype High-Voltage Li-Ion Cathode.

A Li-rich polyanionic compound based on V 3+ with a previously unknown structure, Li 5 VF 4 (SO 4 ) 2 , has been developed as a high-voltage cathode material for Li-ion batteries. The solvothermal preparation of this material, crystal structure solution, and initial electrochemical characterization are presented. An analysis based on density functional theory electronic structure calculations suggests that a high voltage close to 5 V is required to extract two Li ions and to reach the oxidation state of V 5+ . However, the use of conventional carbonate-based electrolytes, which exhibit increasing degradation above a potential of 4.3 V, does not permit the full capacity of this compound to be achieved at this time.

25 ENERGY STORAGE↗

Roadmap on data-centric materials science

Science is and always has been based on data, but the terms ‘data-centric’ and the ‘4th paradigm’ of materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of artificial intelligence and its subset machine learning, has become pivotal in addressing all these challenges. This Roadmap on Data-Centric Materials Science explores fundamental concepts and methodologies, illustrating diverse applications in electronic-structure theory, soft matter theory, microstructure research, and experimental techniques like photoemission, atom probe tomography, and electron microscopy. While the roadmap delves into specific areas within the broad interdisciplinary field of materials science, the provided examples elucidate key concepts applicable to a wider range of topics. The discussed instances offer insights into addressing the multifaceted challenges encountered in contemporary materials research.

36 MATERIALS SCIENCE↗

False metals, real insulators, and degenerate gapped metals

This paper deals with a significant family of compounds predicted by simplistic electronic structure theory to be metals but are, in fact, insulators. This false metallic state has been traditionally attributed in the literature to reflect the absence of proper treatment of electron-electron correlation (“Mott insulators”) whereas, in fact, even mean-field like density functional theory describes the insulating phase correctly if the restrictions posed on the simplistic theory are avoided. Such unwarranted restrictions included different forms of disallowing symmetry breaking described in this article. As the science and technology of conductors have transitioned from studying simple elemental metals such as Al or Cu to compound conductors such as binary or ternary oxides and pnictides, a special class of degenerate but gapped metals has been noticed. Their presumed electronic configurations show the Fermi level inside the conduction band or valence band, yet there is an “internal band gap” between the principal band edges. The significance of this electronic configuration is that it might be unstable toward the formation of states inside the internal band gap when the formation of such states costs less energy than the energy gained by transferring carriers from the conduction band to these lower energy acceptor states, changing the original (false) metal to an insulator. The analogous process also exists for degenerate but gapped metals with the Fermi level inside the valence band, where the energy gain is defined by transfer of electrons from the donor level to the unoccupied part of the valence band. We focus here on the fact that numerous electronic structure methodologies have overlooked some physical factors that could stabilize the insulating alternative, predicting instead false metals that do not really exist (note that this is in general not a physical phase transition, but a correction of a previous error in theory that led to a false prediction of a metal). Such errors include: (i) ignoring spin symmetry breaking, such as disallowing magnetic spin ordering in CuBi2O4 or disallowing the formation of polymorphous spin networks in paramagnetic LaTiO 3 and YTiO 3 ; (ii) ignoring structural symmetry breaking, e.g., not enabling energy-lowering bond disproportionation (Li-doped TiO 2 , SrBiO 3 , and rare-earth nickelates), or not exploring pseudo-Jahn–Teller-like distortions in LaMnO 3 , or disallowing spontaneous formation of ordered vacancy compounds in Ba 4 As 3 and Ag 3 Al 22 O 34 ; and (iii) ignoring spin–orbit coupling forcing false metallic states in CaIrO 3 and Sr 2 IrO 4 . The distinction between false metals vs real insulators is important because (a) predicting theoretically that a given compound is metal even though it is found to be an insulator often creates the temptation to invoke high order novel physical effects (such as correlation in d-electron Mott insulators) to explain what was in effect caused by a more mundane artifact in a lower-level mean-field band theory, (b) recent prediction of exotic physical effects such as topological semimetals were unfortunately based on the above compounds that were misconstrued by theory to be metal, but are now recognized to be stable insulators not hosting exotic effects, and (c) practical technological applications based on stable degenerate but gapped metals such as transparent conductors or electrides for catalysis must rely on the systematically correct and reliable theoretical classification of metals vs insulators.

36 MATERIALS SCIENCE↗

Impact of Noncovalent Interactions on the Structural Chemistry of Thorium(IV)-Aquo-Chloro Complexes

Five novel tetravalent thorium (Th) compounds that consist of Th(H 2 O) x Cl y structural units were isolated from acidic aqueous solutions using a series of nitrogen-containing heterocyclic hydrogen (H) bond donors. Taken together with three previously reported phases, the compounds provide a series of monomeric Th IV complexes wherein the effects of noncovalent interactions (and H-bond donor identity) on Th structural chemistry can be examined. Seven distinct structural units of the general formulas [Th(H 2 O) x Cl 8–x ] x-4 (x = 2, 4) and [Th(H 2 O) x Cl 9–x] x-5 (x = 5–7) are described. The complexes range from chloride-deficient [Th(H 2 O) 7 Cl 2 ] 2+ to chloride-rich [Th(H 2 O) 2 Cl 6 ] 2– species, and theory was used to understand the relative energies that separate complexes within this series via the stepwise chloride addition to an aquated Th cation. Electronic structure theory predicted the reaction energies of chloride addition and release of water through a series of transformations, generally highlighting an energetic driving force for chloride complexation. To probe the role of the counterion in the stabilization of these complexes, electrostatic potential (ESP) surfaces were calculated. Finally, the ESP surfaces indicated a dependence of the chloride distribution about the Th metal center on the pK a of the countercation, highlighting the directing effects of noncovalent interactions (e.g., Hbonding) on Th speciation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simultaneous Optimization of Nuclear–Electronic Orbitals

Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear–electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. In conclusion, benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory-Based Mechanism for Fluoromethane Combustion I: Thermochemistry and Abstraction Reactions

A new detailed chemical kinetic mechanism is presented for small fluorinated hydrocarbons. Ab initio electronic structure theory is used to provide heats of formation with subchemical accuracy. The ANL0 method is extended to include fluorine. The resulting heats of formation at 0 K are in excellent agreement with 36 benchmark species in the Active Thermochemical Tables, with a mean error of μ = –0.02 kJ/mol and a standard deviation of σ = 0.91 kJ/mol. The thermophysical properties for 92 small-molecule H/C/O/F species are computed. The rate coefficients for 40+ H-abstraction reactions involving H, O, F, OH, OF, HO 2 , and various methyl radicals with CH 4 , CH 3 F, CH 2 F 2 , CHF 3 , CH 2 O, and CHFO are discussed. Here, the computed rate constants are in excellent agreement with the available literature. Additionally, 30+ rate constants are provided for F abstraction, which are several orders of magnitude smaller than H abstraction. The thermophysical properties and rate constants are provided in a mechanism. This mechanism is the first in a series of theory-based investigations into the thermal destruction of per- and polyfluorinated species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Studies of surface adsorbate electronic structure and femtochemistry at the fundamental length and time scales. Final report

The electronic structure and ultrafast (10-15 s-femtosecond timescale) electron dynamics were investigated for clean and atom/molecule covered metal surfaces. The studies were performed by scanning tunneling microscopy (STM) to measure the structure of adsorbed atoms and molecules on metal surfaces, and to investigate their electronic properties. The electronic structure of the observed molecular networks was calculated by electronic structure theory in collaboration with Prof. Jin Zhao, who is a long-time collaborator, a Professor at the University of Science and Technology of China, and holds an Adjunct Professorship at the University of Pittsburgh. We also investigated the electronic properties of C60 molecules when they are templated by corrugated black phosphorous surfaces. We found unexpected charge delocalization that is enabled by the templating. This research was done in collaboration with Professor Min Feng at the Wuhan University, and who also holds an Adjunct Professorship at the University of Pittsburgh. Moreover, the electronic structure and electron dynamics in metal surfaces were investigated by time-resolved photoemission electron spectroscopy. The focus of ultrafast spectroscopy has been on the plasmonic response of silver surfaces. One direction has been to develop multidimensional (energy, momentum, and time) photoelectron spectroscopy of the coherent response of solid surfaces. This method was applied to study the collective electron excitations known generally as plasmons, which screen optical fields from penetration into metals. Although this collective response has been known for more than 60 years and is used extensively to deposit optical energy into metals, how this happens is poorly known. We investigated the plasmonic response of silver at the point where the dielectric response passes through zero and bulk plasmon is excited by light. We discovered that the plasmon excitation decays by exciting electrons from the Fermi level of a metal, which is contrary to what is believed in the plasmonic science community. This research has been performed in collaboration with Dr. Marcel Reutzel, who was a postdoctoral fellow working on this research at the University of Pittsburgh, and now has a faculty position at the University of Göttingen in Germany. Prof. Branko Gumhalter from the Institute of Physics in Zagreb contributed on the theory of plasmonic decay processes. Furthermore, we investigated the Floquet engineering of electronic bands in metals leading to multiphoton photoemission and above threshold photoemission. Finally, we demonstrated that it is possible to change the electronic structure of metals by application of optical fields. Our studies indicated that this happens on subfemtosecond time scale and could potentially be used in ultrafast information processing and quantum computation. Related document information

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pseudogap in elemental plutonium

Electronic correlations associated with incipient magnetism have long been recognized as an important factor in stabilizing the largest atomic volume δ phase of plutonium, although their strength compared to those in the rare earths and neighboring actinides in the periodic table has largely remained a mystery. We show here using calorimetry measurements, together with prior detailed measurements of the phonon dispersion, that the 5f electrons of the δ phase reside in a pseudogapped state, accompanied by reductions in various physical properties below a characteristic temperature T* ≈ 100 K. The small characteristic energy scale of the pseudogapped state implies that the 5f electrons in plutonium are much closer to the threshold for localization and magnetic order than has been suggested by state-of-the-art electronic structure theory, revealing plutonium to be arguably the most strongly correlated of the elements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantifying Spin-Mixed States in Ferromagnets

We quantify the presence of spin-mixed states in ferromagnetic 3D transition metals by precise measurement of the orbital moment. While central to phenomena such as Elliot-Yafet scattering, quantification of the spin-mixing parameter has hitherto been confined to theoretical calculations. We demonstrate that this information is also available by experimental means. Comparison of ferromagnetic resonance spectroscopy with x-ray magnetic circular dichroism results show that Kittel’s original derivation of the spectroscopic g factor requires modification, to include spin mixing of valence band states. Our results are supported by ab initio relativistic electronic structure theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wavelength-Dependent Excitonic Properties of Covalent Organic Frameworks Explored by Theory and Experiments

Many aspects of the correlation between physical structure, light harvesting, and excitonic properties of covalent organic frameworks (COFs) remain unclear despite being key properties determining their photocatalytic function. One area of COF research that could bring clarity is through using both electronic structure theory and time-resolved spectroscopic analysis over a series of isomeric COFs. Here, we show structure-property relationships between four imine COFs built from a combination of ditopic and tritopic monomers using transient absorption spectroscopy together with time-dependent density functional theory. We find that monomer selection only moderately affects the charge transfer (CT) behavior of the COFs. Instead, we infer that imine chemistry profoundly impacts CT by acting as a CT mediator. Moreover, we discover two distinct valence bands arising from varying degrees of locally excited/CT mixing, which is responsible for energy dependent exciton dynamics. Finally, we use theory to hypothesize that interlayer interactions can modify excitonic properties that we correlate with tail states commonly observed, but rarely investigated in COFs. Furthermore, these results reveal that imine chemistry should be recognized as a very important factor to consider in the development of COF photocatalysts and the correlation of their structural environment with light harvesting and CT properties that should ultimately determine their photocatalytic function.

14 SOLAR ENERGY↗

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO 2 C 2 H 7 + and PO 2 CH 4 + , oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO 2 C 2 H 7 + and PO 2 CH 4 + . The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

59 BASIC BIOLOGICAL SCIENCES↗

Vibrational and Rovibrational Spectroscopy Applied to Astrochemistry

The detection of molecules in astrophysical environments almost always requires remote sensing. While radioastronomical observation and associated rotational spectroscopy are powerful astronomical tools, infrared spectral analysis provides a unique means of examining the observable universe, especially for molecules where permanent dipole moments are small or even non-existent. The molecular vibrations of small molecules are now able to be modeled via quantum chemistry and electronic structure theory conjoined to vibrational analysis to within spectroscopic accuracy in many cases. This chapter will showcase this success and build upon it to show how such advances are now being leveraged to describe molecular vibrations for molecules containing dozens of atoms, electronically excited states, "hot bands," exoplanetary atmospheric opacity data, and even emission of polycyclic aromatic hydrocarbons. All of these are required to prune the interstellar spectral garden of its "weeds" in search of "flowers" that will provide the necessary fingerprints for astronomers to be able to probe the heavens for its past, present, and future secrets.

Ryan C. Fortenberry↗

Structure, electronic and magnetic characterization, and calculated electronic structures of two oxyhalide hexagonal perovskites

Herein we report the crystal structures, initial magnetic and charge transport characterization, and calculated electronic structures of the hexagonal oxyhalide perovskites, Ba 7 Ru 4 O 15 Cl 2 and Ba 7 Ru 4 O 15 Br 2 . The experimental information is obtained through the study of single crystals. Face-sharing Ru O 6 octahedra form Ru 2 O 9 dimers in a layered triangular geometry in these materials; minor amounts of off-magnetic-site structural disorder are present in our crystals. Both are magnetically isotropic, with 2.5 μB/mol-Ru, Curie-Weiss theta –185 K, and 2.9 μB/mol-Ru, Curie-Weiss theta –168 K, respectively. Broad features in the magnetic susceptibility and heat capacity associated with magnetic ordering (at 35 and 37 K, respectively) are observed. The charge transport band gaps are 0.03 eV for the oxychloride and 0.006 eV for the oxybromide. There is no gap at the Fermi level for either of these compounds in density-functional theory electronic structure calculations.

36 MATERIALS SCIENCE↗