Random forest machine learning models for interpretable X-ray absorption near-edge structure spectrum-property relationships
Not Available
36 MATERIALS SCIENCE↗
SEARCH · Search NASA
Results for “Theoretical Chemistry”
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.
Vibrationally resolved optical excitations of the nitrogen-vacancy center in diamond
A comprehensive description of the optical cycle of spin defects in solids requires the understanding of the electronic and atomistic structure of states with different spin multiplicity, including singlet states which are particularly challenging from a theoretical standpoint. We present a general framework, based on spin-flip time-dependent density function theory, to determine the excited state potential energy surfaces of the many-body singlet states of spin defects; we then predict the vibrationally resolved absorption spectrum between singlet shelving states of a prototypical defect, the nitrogen-vacancy center in diamond. Our results, which are in very good agreement with experiments, provide an interpretation of the measured spectra and reveal the key role of specific phonons in determining absorption processes, and the notable influence of non-adiabatic interactions. The insights gained from our calculations may be useful in defining strategies to improve infrared-absorption-based magnetometry and optical pumping schemes. The theoretical framework developed here is general and applicable to a variety of other spin defects and materials.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Circumventing the phonon bottleneck by multiphonon-mediated hot exciton cooling at the nanoscale
Abstract Nonradiative processes govern efficiencies of semiconductor nanocrystal (NC)-based devices. A central process is hot exciton cooling, or the nonradiative relaxation of a highly excited electron/hole pair to form a band-edge exciton. Due to quantum confinement effects, the timescale and mechanism of cooling are not well understood. A mismatch between electronic energy gaps and phonon frequencies has led to the hypothesis of a phonon bottleneck and extremely slow cooling, while enhanced electron-hole interactions have suggested ultrafast cooling. Experimental measurements of the cooling timescale range six orders of magnitude. Here, we develop an atomistic approach to describe phonon-mediated exciton dynamics and simulate cooling in NCs of experimentally relevant sizes. We find that cooling occurs on ~30 fs timescales in CdSe NCs, in agreement with the most recent measurements, and that the phonon bottleneck is circumvented through a cascade of multiphonon-mediated relaxation events. Furthermore, we identify NC handles for tuning the cooling timescale.
36 MATERIALS SCIENCE↗
Candidate ferroelectrics via ab initio high-throughput screening of polar materials
Ferroelectrics are a class of polar and switchable functional materials with diverse applications, from microelectronics to energy conversion. Computational searches for new ferroelectric materials have been constrained by accurate prediction of the polarization and switchability with electric field, properties that, in principle, require a comparison with a nonpolar phase whose atomic-scale unit cell is continuously deformable from the polar ground state. For most polar materials, such a higher-symmetry nonpolar phase does not exist or is unknown. Here, we introduce a general high-throughput workflow that screens polar materials as potential ferroelectrics. We demonstrate our workflow on 1978 polar structures in the Materials Project database, for which we automatically generate a nonpolar reference structure using pseudosymmetries, and then compute the polarization difference and energy barrier between polar and nonpolar phases, comparing the predicted values to known ferroelectrics. Focusing on a subset of 182 potential ferroelectrics, we implement a systematic ranking strategy that prioritizes candidates with large polarization and small polar-nonpolar energy differences. To assess stability and synthesizability, we combine information including the computed formation energy above the convex hull, the Inorganic Crystal Structure Database id number, a previously reported machine learning-based synthesizability score, and ab initio phonon band structures. To distinguish between previously reported ferroelectrics, materials known for alternative applications, and lesser-known materials, we combine this ranking with a survey of the existing literature on these candidates through Google Scholar and Scopus databases, revealing ~130 promising materials uninvestigated as ferroelectric. Our workflow and large-scale high-throughput screening lays the groundwork for the discovery of novel ferroelectrics, revealing numerous candidates materials for future experimental and theoretical endeavors.
36 MATERIALS SCIENCE↗
P SI4 1.4: Open-source software for high-throughput quantum chemistry
PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4 's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. We announce a rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Pulsed power accelerator surface Joule heating models
Understanding the effects of contaminant plasmas generated within the Z machine at Sandia is critical to understanding current loss mechanisms. The plasmas are generated at the accelerator electrode surfaces and include desorbed species found in the surface and substrate of the walls. These desorbed species can become ionized. The timing and location of contaminant species desorbed from the wall surface depend non-linearly on the local surface temperature. For accurate modeling, it is necessary to utilize wall heating models to estimate the amount and timing of material desorption. One of these heating mechanisms is Joule heating. In this paper, we propose several extended semi-analytic magnetic diffusion heating models for computing surface Joule heating and demonstrate their effects for several representative current histories. We quantitatively assess under what circumstances these extensions to classical formulas may provide a validatable improvement to the understanding of contaminant desorption timing.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ultrafast Metal Electrodeposition Revealed by In Situ Optical Imaging and Theoretical Modeling towards Fast–Charging Zn Battery Chemistry
Metallic Zn is a preferred anode material for rechargeable aqueous batteries towards a smart grid and renewable energy storage. Importantly, understanding how the metal nucleates and grows at the aqueous Zn anode is a critical and challenging step to achieve full reversibility of Zn battery chemistry, especially under fast-charging conditions. Here, by combining in situ optical imaging and theoretical modeling, we uncover the critical parameters governing the electrodeposition stability of the metallic Zn electrode, that is, the competition among crystallographic thermodynamics, kinetics, and Zn 2+ -ion diffusion. Moreover, steady-state Zn metal plating/ stripping with Coulombic efficiency above 99 % is achieved at 10-100 mA cm -2 in a reasonably high concentration (3 M) ZnSO 4 electrolyte. Significantly, a long-term cycling-stable Zn metal electrode is realized with a depth of discharge of 66.7% under 50 mA cm -2 in both Zn || Zn symmetrical cells and MnO 2 || Zn full cells.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Light-induced electron dynamics of molecules in cavities: comparison of model Hamiltonians
The rapidly expanding field of polaritonic chemistry requires accurate theoretical simulations to understand new phenomena at the atomic scale. Computing the optoelectronic properties of molecules using established electronic structure methods is a careful balance of accuracy and computational expense, and expanding these methods to quantum electrodynamics to describe coupled cavity-molecule systems is an active topic of development. Key to these methods are the Hamiltonian operators representing the photon cavity modes. The recently introduced quantum electrodynamics time-dependent configuration interaction (QED-TDCI) method allows for the combination of electron dynamics simulations with quantum electrodynamics, enabling the simulation of time-dependent optoelectronic properties of cavity-molecule systems. Using this method, a comparison of two many-state QED Hamiltonians—the Pauli-Fierz and quantum Rabi model Hamiltonians—is presented, with a particular focus on time-dependent properties in applied electric fields.
Peyton, Benjamin G. (ORCID:0000000163190751)↗
The ezSpectra suite: An easy‐to‐use toolkit for spectroscopy modeling
Abstract A molecule's spectrum encodes information about its structure and electronic properties. It is a unique fingerprint that can serve as a molecular ID. Quantum chemistry calculations provide key ingredients for interpreting spectra, but modeling the spectra rarely ends there; it requires additional steps that entail combined treatments of electronic and nuclear degrees of freedom and account for specifics of the experimental setup (light energy, polarization, averaging over molecular orientations, temperature, etc.). This Software Focus article describes the ezSpectra suite, which currently comprises two stand‐alone open‐source codes: ezFCF and ezDyson . ezFCF calculates Franck–Condon factors, which yield vibrational progressions for polyatomic molecules, within the double‐harmonic approximation. ezDyson calculates absolute cross‐sections for photodetachment/photoionization processes and photoelectron angular distributions using Dyson orbitals computed by a quantum chemistry program. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods
Gozem, Samer↗
On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond
Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.
Chemistry↗
Spectroscopic and theoretical studies of UN and UN +
The low-energy electronic states of UN and UN+ have been examined using high-level electronic structure calculations and two-color photoionization techniques. The experimental measurements provided an accurate ionization energy for UN (IE = 50 802 +/- 5 cm -1 ). Spectra for UN+ yielded ro-vibrational constants and established that the ground state has the electronic angular momentum projection omega = 4. Ab initio calculations were carried out using the spin-orbit state interacting approach with the complete active space second-order perturbation theory method. A series of correlation consistent basis sets were used in conjunction with small-core relativistic pseudopotentials on U to extrapolate to the complete basis set limits. The results for UN correctly obtained an omega = 3.5 ground state and demonstrated a high density of configurationally related excited states with closely similar ro-vibrational constants. Similar results were obtained for UN + , with reduced complexity owing to the smaller number of outer-shell electrons. The calculated IE for UN was in excellent agreement with the measured value. Improved values for the dissociation energies of UN and UN + , as well as their heats of formation, were obtained using the Feller-Peterson-Dixon composite thermochemistry method, including corrections up through coupled cluster singles, doubles, triples and quadruples. An analysis of the ab initio results from the perspective of the ligand field theory shows that the patterns of electronic states for both UN and UN + can be understood in terms of the underlying energy level structure of the atomic metal ion.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape
Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
First-principles modeling of chemistry in mixed solvents: Where to go from here?
Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Theoretical insights into the surface physics and chemistry of redox-active oxides
Redox active oxides are ubiquitous in materials applications, including catalysis, photovoltaics, self-cleaning glasses, chemical sensors, and electronic components. Their utility derives from their unique ability to access multiple metal charge states within a finite energy window. However, this property also confounds our ability to study reducible oxides, because it leads to structural, compositional, and electronic complexities that elude simplistic models of materials structure and function. Oxygen vacancies play a critical role in shaping the functional properties of such oxides; most notably, they lead to mobile charge imbalances that impact surface processes at significant distances from the originating defect. Atomistic simulations are inherently equipped to illuminate these phenomena at a fundamental level; however, reducible oxides pose great challenges due to the high level of electron correlation needed to correctly describe them. Knowing how defects form, couple, propagate, agglomerate, or repel each other and influence the surface properties of reducible oxides is only now coming into the grasp of modern theory and simulation capabilities. This knowledge is also key to discovering and controlling emergent materials processes at nanometer scale and beyond. The manuscript was written with support of all three authors from U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. Work by RR and VG was performed at Pacific Northwest National Laboratory (PNNL), which is a multi-program laboratory operated by Battelle for DOE under Contract DE AC05 76RL01830. AS was supported under Award DE-SC0007347.
Rousseau, Roger J.↗
MXene-Derived Potassium-Preintercalated Bilayered Vanadium Oxide Nanostructures for Cathodes in Nonaqueous K-Ion Batteries
Bilayered vanadium oxides (BVOs) are promising cathode materials for beyond-Li-ion batteries due to their tunable chemistries and high theoretical capacities. However, the large size of beyond-Li + ions limits electrochemical cycling and rate capability of BVO electrodes. Recent reports of MXene-derived BVOs with nanoscale flower-like morphology have shown improved electrochemical stability at high rates up to 5C in nonaqueous lithium-ion batteries. Here, we report how morphological stabilization can lead to improved rate capability in potassium-ion batteries (PIBs) through the synthesis and electrochemical characterization of MXene-derived K-preintercalated BVOs (MD-KVOs), which were derived from two V 2 CT x precursor materials prepared using two different etching protocols. We show that the etching conditions affect the surface chemistry of the MXene, which plays a role in the MXene-to-oxide transformation process. MXene derived from a milder etchant transformed into a nanoflower MD-KVO with two-dimensional (2D) nanosheet petals (KVO-DMAE) while a more aggressive etchant produced a MXene that transformed into a MD-KVO with one-dimensional (1D) nanorod morphology (KVO-CMAE). Electrochemical cycling of the produced MD-KVOs after drying at 200 °C under vacuum (KVO-DMAE-200 and KVO-CMAE-200) in PIBs showed that electrochemical stability of MD-KVO at high rates improved through the morphological stabilization of 2D particles combined with the control of interlayer water and K + ion content. Structure refinement of KVO-DMAE-200 further corroborates the behavior observed during K + ion cycling, connecting structural and compositional characteristics to the improved rate capability. This work demonstrates how proper synthetic methodology can cause downstream effects in the control of structure, chemical composition, and morphology of nanostructured layered oxide materials, which is necessary for development of future materials for beyond-Li-ion battery technologies.
25 ENERGY STORAGE↗
Understanding flavin electronic structure and spectra
Abstract Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Theoretical and Physical Chemistry > Spectroscopy
59 BASIC BIOLOGICAL SCIENCES↗
Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules
Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy
Chemistry↗
Robust Inversion of Time-Resolved Data via Forward-Optimization in a Trajectory Basis
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗