Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemistry 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 469 records · Page 26

Interplay of coulomb and exciton–phonon coupling controls singlet fission dynamics in two pentacene polymorphs

Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.

14 SOLAR ENERGY↗

The promise of optical NMR spectroscopy for experimental aqueous geochemistry

Here, we discuss new quantum technologies being adapted to detect extraordinarily small magnetic fields. The signals are derived from the stimulated emission of light from single-atom defects in diamonds. The intensity of light reports the spin state of the electrons in the color center. The physics is relevant to geochemists because it points to a means of conducting NMR experiments on solutions in diamond-anvil cells. Because the method is so sensitive, it can detect signals in subnanoliter volumes, making work at elevated temperatures and pressures possible, and also on solutions that are dangerous in larger volumes.

58 GEOSCIENCES↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Neutrons in Structural Biology: Challenges and Opportunities (Workshop Report)

Gaining a thorough understanding of biological systems requires building our knowledge about biological processes from the level of atoms and electrons, and up to whole organisms. Such comprehensive knowledge will allow for a predictive understanding of complex biological systems behavior. It will guide us in the design and development of novel therapeutics and vaccines to tackle existing health threats and to prepare for future pandemics, and it will provide information necessary to create new biomaterials and bio-inspired technologies through manipulation of biological macromolecules, their assemblies, single cells and even microorganisms. Reaching these goals will require a synergistic combination of multiple experimental techniques with molecular calculations and predictive simulations, and the design and development of new techniques and capabilities that bridge current knowledge and technology gaps. Neutron scattering provides unique information about the biomacromolecular structure and function and can play a major role in achieving these goals. A workshop was held to engage the scientific community in identifying pressing challenges in biochemistry, structural biology, enzymology and structure-guided drug design not solved with the current neutron scattering technologies or utilizing other structural biology techniques such as X-ray crystallography, NMR, and cryo-EM. The workshop brought together structural biology, biochemistry and computational experts, as well as early career researchers and students, creating a forum for discussing scientific advancement and collaboration. The workshop included a one-day satellite training workshop where graduate students and postdoctoral researchers were educated in the application of neutron crystallography and small-angle scattering in structural biology. Furthermore, the Instrument Scientific Advisory Board (ISAB) for the development of a macromolecular neutron diffractometer at ORNL’s Second Target Station was introduced at the workshop. The major outcome was that neutrons can provide atomic-level understanding of biomacromolecular structure, function and dynamics which is of paramount importance for addressing the identified challenges. Neutron crystallography, in particular, can resolve long-standing biochemical issues regarding enzyme function by delineating the underlying chemistry and can have a major impact on the design of small-molecule therapeutics, especially in combination with molecular computation (quantum chemistry and molecular dynamics simulations) and the emerging artificial intelligence (AI)-assisted drug design technologies. The unique properties of neutrons, including their high sensitivity to hydrogen and their non-destructive nature, make them ideal probes of biological matter. There is a palpable need in the scientific community to expand and enhance the impact of neutron sciences on biology. Neutron crystallography is the only structural biology method capable of determining positions of all hydrogen atoms in proteins, nucleic acids and their complexes at near-physiological temperatures and of unstable species at cryogenic temperatures. Moreover, neutron analysis is non-ionizing, non-destructive and does not perturb the structure or redox chemistry of active site metal centers and clusters in proteins, which can be invaluable for studying radiation-sensitive metalloprotein complexes. Further, neutron energies used in scattering applications are similar to atomic motions, permitting neutron spectroscopies to characterize the dynamics of biomacromolecules on the picosecond to microsecond timescales. The different sensitivities of neutrons to protium (H) and deuterium (D) isotopes of hydrogen allow enhanced visibility of specific parts of biological complexes through isotopic labeling. The impact of neutrons will be most powerful when neutron scattering is combined with complementary experimental techniques that use photons and electrons, and with high-performance computing. The interconnection and mutuality of the experimental and theoretical capabilities will drive discoveries in biological and health sciences to generate more complete picture of complex biological systems. The major limitation in the field of biological neutron crystallography has been signal-to-noise, demanding large samples that are difficult to produce for the majority of biomacromolecules and limiting the applicability of this technique in biological sciences. A neutron crystallography instrument at the Second Target Station will revolutionize biological science with neutrons by engaging a large scientific community of structural biologists, enabling successful neutron diffraction experiments from radically smaller biomacromolecular crystals, resolving unanswered biochemical questions, and meaningfully contributing to rational drug design. The meeting highlighted 10 grand challenges that will be addressed with this advanced capability over the next decade and beyond, and the recommendations required to help address them are given below.

59 BASIC BIOLOGICAL SCIENCES↗

Characterizing the Excited States and Electronic Absorption Spectra of Small Alkylperoxy (RO 2 •) and Hydroperoxy (•QOOH) Radicals

Organic peroxy (ROO•) and hydroperoxy (•QOOH) radicals are key reactive intermediates that are formed via the oxidation of volatile organic compounds during combustion or in the Earth’s atmosphere. Their primary fate is continued unimolecular decay or bimolecular chemistry, the relative branching for which is heavily structure- and temperature-dependent. This article outlines a combined single- and multi-reference quantum chemical study to characterize the near-UV accessible electronically excited states of the prototypical ROO• and •QOOH intermediates, tert-butyl peroxy and hydroperoxy-tert-butyl radicals—the ground-state chemistries of which have been well studied both experimentally and computationally. Additionally, we simulate the electronic absorption profiles of these ROO• and •QOOH intermediates with a variety of multi- and single-reference methods. The results show an interesting conformer dependence on the electronically excited-state character and electronic absorption maxima of •QOOH. The results show promise for electronic absorption spectroscopy to be used as a selected probe for determining •QOOH conformers. Additionally, electronic absorption may contribute to the daytime removal of long-lived •QOOH intermediates formed in the troposphere. We expect that our studies will motivate experiments on the electronic absorption spectra of experimentally achievable ROO• and •QOOH.

excited states↗

Vibrational Spectra of HNIW and its Isotopologues: A Combined Experimental and Computational Study

Incorporating isotopically labelled materials in degradation experiments could help unravel the mechanism(s) of decomposition through use of the kinetic isotope effect. Characterizing synthesized isotopologues however requires an understanding of what observable signals are affected by the isotopic substitution. As vibrational spectroscopy can distinguish between isotopologues, it is an ideal characterization technique to evaluate isotopic variants. To this end, the vibrational spectra of HNIW and its deuterated ( 2 H), 13 C, 15 N (all), 15 N (nitro), and 18 O isotopologues have been computationally predicted in the gas phase using density functional theory. These results are compared to experimentally measured FTIR/ATR and Raman spectra of both unsubstituted HNIW and 15 N-labeled HNIW in which the six nitro groups were synthetically tagged with 15 N atoms ( 15 N nitro -HNIW). The experimental isotopic frequency shift for the -NO 2 asymmetric stretching frequencies agrees with that theoretically calculated (~35.7 cm-1 vs. 36.5cm-1, respectively). Furthermore, analysis of the theoretically predicted frequency shifts for all isotopologues suggest the -NO 2 bending modes are lower in frequency than previously reported. This assignment is supported by the experimentally measured isotopic shift of ~10.1 cm -1 for these features (consistent with the predicted shift of ~13.1 cm -1 ). This work expands our current understanding of the vibrational modes in HNIW as well as provides a method for future work on similar systems.

36 MATERIALS SCIENCE↗

QForte: An Efficient State-Vector Emulator and Quantum Algorithms Library for Molecular Electronic Structure

Here, we introduce a novel open-source software package QForte, a comprehensive development tool for new quantum simulation algorithms. QForte incorporates functionality for handling molecular Hamiltonians, fermionic encoding, ansatz construction, time evolution, and state-vector emulation, requiring only a classical electronic structure package as a dependency. QForte also contains black-box implementations of a wide variety of quantum algorithms, including variational and projective quantum eigensolvers, adaptive eigensolvers, quantum imaginary time evolution, and quantum Krylov methods. We highlight two features of QForte: (i) how the Python class structure of QForte enables the facile implementation of new algorithms, and (ii) how existing algorithms can be executed in just a few lines of code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗

Impact of Humidity and Temperature on the Stability of the Optical Properties and Structure of MAPbI3, MA0.7FA0.3PbI3 and (FAPbI3)0.95(MAPbBr3)0.05 Perovskite Thin Films

In situ real-time spectroscopic ellipsometry (RTSE) measurements have been conducted on MAPbI3, MA0.7FA0.3PbI3, and (FAPbI3)0.95(MAPbBr3)0.05 perovskite thin films when exposed to different levels of relative humidity at given temperatures over time. Analysis of RTSE measurements track changes in the complex dielectric function spectra and structure, which indicate variations in stability influenced by the underlying material, preparation method, and perovskite composition. MAPbI3 and MA0.7FA0.3PbI3 films deposited on commercial fluorine-doped tin oxide coated glass are more stable than corresponding films deposited on soda lime glass directly. (FAPbI3)0.95(MAPbBr3)0.05 films on soda lime glass showed improved stability over the other compositions regardless of the substrate, and this is attributed to the preparation method as well as the final composition.

36 MATERIALS SCIENCE↗

Atomic Resolution Cryogenic 4D-STEM Imaging via Robust Distortion Correction

Cryogenic four-dimensional scanning transmission electron microscopy (4D-STEM) imaging is a useful technique for studying quantum materials and their interfaces by simultaneously probing charge, lattice, spin, and chemistry on the atomic scale with the sample held at temperatures ranging from room to cryogenic. However, its applications are currently limited by the instabilities of cryo-stages and electronics. To overcome this challenge, we develop an algorithm to effectively correct the complex distortions present in atomic resolution cryogenic 4D-STEM data sets. This method uses nonrigid registration to identify localized distortions in a 4D-STEM and relate them to an undistorted experimental STEM image, followed by a series of affine transformations for distortion corrections. This method allows a minimum loss of information in both reciprocal and real spaces, enabling the reconstruction of sample information from 4D-STEM data sets. This method is computationally cheap, fast, and applicable for on-the-fly data analysis in future in situ cryogenic 4D-STEM experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A quantum theoretical study of polyimides

One of the most important contributions of theoretical chemistry is the correct prediction of properties of materials before any costly experimental work begins. This is especially true in the field of electrically conducting polymers. Development of the Valence Effective Hamiltonian (VEH) technique for the calculation of the band structure of polymers was initiated. The necessary VEH potentials were developed for the sulfur and oxygen atoms within the particular molecular environments and the explanation explored for the success of this approximate method in predicting the optical properties of conducting polymers.

Burke, Luke A.↗

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↗

A Physics-Based Approach to Hypersonic Nonequilibrium Chemistry

NASA has undertaken a concerted effort to generate a complete model for the chemical reactions and other collisional processes that take place in the bow shock layer and wake of spacecraft entering the atmospheres of Earth and Mars. This new model is based on first-principles physics calculations: interaction potentials between colliding species obtained from accurate quantum mechanics calculations and reaction rate coefficients and cross sections based on classical mechanics simulations of heavy particle collisions.This presentation will give an introduction to the underlying theory and computational methods used in this work and describe the range of chemical and physical processes that take place in the close vicinity of spacecraft during atmospheric entry. I will also present some of our recent results for reactions involving nitrogen, carbon monoxide and CO2.

Jaffe, Richard↗

High Resolution Infrared Spectroscopy of Highly Reactive Chemical Intermediates: Berkeley Inspiration and a C.B. Moore Retrospective

One of the long-standing paradoxes in chemistry is that the molecular species of greatest interest are often highly reactive chemical intermediates and thus present only at vanishingly small concentrations under typical steady state conditions. Here, this has proven both a key frustration and powerful motivation for physical chemists in probing such species, particularly for high resolution laser spectroscopic methods in the infrared. Over the past 4-5 decades, this fundamental sensitivity challenge has stimulated the development of many novel techniques for efficient generation, cooling, and infrared probing of such transient chemical intermediates, achieving increasingly sophisticated levels of spectroscopic detail, chemical insight, and a testbed for comparison with first principles ab initio quantum calculations. Two of the pioneering contributors to this scientific arena have been George Pimentel and C. Bradley Moore, both responsible for nourishing DJN’s own enduring fascination with spectroscopic IR study of elusive chemical intermediates. This chapter is not intended to provide a comprehensive review of this large and enormously successful field, but it rather more simply attempts to capture a few selective “Moore-centric” snapshots of the Nesbitt group scientific evolution, specifically in i) advances by Pimentel for the original development of rapid scan flash kinetic spectroscopy based on spectrometers with rapidly rotating IR gratings, ii) the modernization/extension of these methods in the Moore group toward the first high resolution IR laser spectroscopy of singlet/triplet methylene 1,3 CH 2 diradical, and finally to iii) the development of slit supersonic discharge expansion methods at JILA for sub-Doppler infrared laser study of multiple radicals, jet cooled molecular ions, and highly reactive chemical intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

Non-Adiabatic Excited State Molecular Dynamics Methodologies: comparison and convergence

Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗