Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

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

A most misunderstood conditionally-solvable quantum-mechanical model

Highlights: • The Schrödinger equation for some quantum-mechanical models is separable in cylindrical coordinates. • The radial equation exhibits harmonic, linear and Coulomb-like interactions. • The Frobenius method leads to three-term recurrence relations. • Some particular energies are obtained from truncation of the recurrence relation. • Many authors misunderstood these results. In this paper we show that several authors have derived wrong physical conclusions from a gross misunderstanding of the exact eigenvalues and eigenfunctions of a conditionally-solvable quantum-mechanical model. It consists of an eigenvalue equation with seemingly Coulomb, linear and harmonic terms. Here we compare the results derived by those authors with the actual eigenvalues of the models calculated by means of the Ritz variational method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Exploring the elastic and electronic properties of chromium molybdenum diboride alloys

We perform first-principles calculations to study the structural, mechanical, thermal, electronic, and magnetic properties of Cr 1-x Mo x B 2 for x = 0.25, 0.33, 0.50, 0.67 and 0.75. Based on structural search methods, we determine the ground-state structure for each concentration. The ternaries are either monoclinic (x = 0.25, 0.75) or trigonal (x = 0.33, 0.50, 0.67). The calculated mechanical properties reveal that the strength of Cr 1-x Mo x B 2 is maximized for x = 0.50. Cr 0.5 Mo 0.5 B 2 exhibits excellent mechanical properties (B = 298 GPa, Y = 558 GPa, G = 235 Gpa, ν = 0.19, H v =27 GPa), surpassing those of β-MoB 2 at a lower cost. All of these ternaries are hard alloys with Vickers hardness greater than 24 GPa. Chemical bonding analysis demonstrates that the strength of the new compounds is related to the alternating planar and buckled B-B layers, as well as the strong TM-B bonds. Finally, the enhanced strength of Cr 0.5 Mo 0.5 B 2 is a consequence of the high density of strong interlayer Cr-Mo metallic bonds around the Fermi level.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two-Dimensional Moiré Polaronic Electron Crystals

Two-dimensional moiré materials have emerged as the most versatile platform for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe 2 /WS 2 moiré superlattices. We find that ultrafast electronic excitation leads to partial melting of the Mott states on timescales 5 times longer than predictions from the charge hopping integrals and that the melting rates are thermally activated, with activation energies of 1⁢8±3 and 1⁢3±2 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. A density functional theory calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Interaction effects in graphene in a weak magnetic field

A weak perpendicular magnetic field, B, breaks the chiral symmetry of each valley in the electron spectrum of graphene, preserving the overall chiral symmetry in the Brillouin zone. Here we explore the consequences of this symmetry breaking for the interaction effects in graphene. In particular, we demonstrate that the electron-electron interaction lifetime acquires an anomalous B dependence. Also, the ballistic zero-bias anomaly, δv(ω), where ω is the energy measured from the Fermi level, emerges at a weak B and has the form δv(B) ~ B 2 /ω 2 . Temperature dependence of the magnetic-field corrections to the thermodynamic characteristics of graphene is also anomalous. We discuss experimental manifestations of the effects predicted. The microscopic origin of the B-field sensitivity is an extra phase acquired by the electron wave function resulting from the chirality-induced pseudospin precession.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

All about the neutron from lattice QCD

I describe how simulations of lattice QCD using the path integral formulation provide the two basic quantum mechanical properties of QCD, its ground state in which correlation functions are calculated, and Fock state wavefunctions between which matrix elements of operators are calculated. Both constructs are stochastic, so unfortunately one gets no intuitive picture or even qualitative understanding of what they look like, nevertheless they contain and display all the subtleties of the quantum field theory. Today, these simulations provide many quantities that are impacting phenomenology and experiments. I illustrate the methods and the steps in the analysis using, as examples, three observables: the isovector charges of the nucleon, the contribution of the quark’s intrinsic spin to the nucleon spin, and the pion-nucleon sigma term.

Gupta, Rajan (ORCID:0000000317843058)↗
DOE OSTI2026

Pb-Apatite Framework as a Generator of Novel Flat-Band CuO-Based Physics

Based upon density functional theory (DFT) calculations, we present the basic electronic structure of CuPb 9 (PO 4 ) 6 O (Cu-doped lead apatite, aka LK-99), in two scenarios: (1) where the structure is constrained to the P3 symmetry and (2) where no symmetry is imposed. At the DFT level, the former is predicted to be metallic while the latter is found to be a charge-transfer insulator. In both cases the filling of these states is nominally d 9 , consistent with the standard Cu 2+ valence state, and Cu with a local magnetic moment of order 0.7 μ B . In the metallic case we find these states to be unusually flat (∼0.2 eV dispersion), giving a very high density of electronic states (DOS) at the Fermi level that we argue can be a host for novel electronic physics. The flatness of the bands is the likely origin of symmetry-lowering gapping possibilities that would remove the spectral weight from E F . Motivated by some initial experimental observations of metallic or semiconducting behavior, we propose that disorder (likely structural) is responsible for closing the gap. Here, we consider a variety of possibilities that could possibly close the charge-transfer gap but limit consideration to kinds of disorder that preserve electron count. Of the possible kinds we considered (spin disorder, O populating vacancy sites, and Cu on less energetically favorable Pb sites), the local Cu moment, and consequently the charge-transfer gap, remains robust. We conclude that disorder responsible for metallic behavior entails some kind of doping where the electron count changes. Further, we claim that the emergence of the flat bands should be due to weak wave function overlap between the orbitals on Cu and O sites, owing to the directional character of the constituent orbitals. Therefore, finding an appropriate host structure for minimizing hybridization between Cu and O while allowing them to still weakly interact should be a promising route for generating flat bands at E F which can lead to interesting electronic phenomena, regardless of whether LK-99 is a superconductor.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

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

Characterizing the Mechanisms of Ca and Mg Carbonate Ion-Pair Formation with Multi-Level Molecular Dynamics/Quantum Mechanics Simulations

The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Furthermore, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.

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

Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H 2 → SH + H Reaction

We study the spin-dependent stereodynamics of the S + H 2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3 A" states and singlet 1 A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet–triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Furthermore, our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.

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

Anion Photoelectron Spectroscopy and Ab Initio Studies of the UF – Anion

A synergistic anion photoelectron spectroscopic and ab initio computational study of photodetachment of UF – is reported. The measurement determined a vertical detachment energy of 0.63(03) eV, which is consistent with a spinor-based relativistic coupled-cluster CCSD(T) value of 0.61 eV. The complex spectral features due to excited electronic states and vibrational progressions of UF are analyzed and assigned with the help of spin–orbit-coupled multireference perturbation theory and spinor-based relativistic coupled-cluster calculations. UF and UF – are confirmed to be dominated by ionic bonding. Furthermore, the usefulness of the spinor CCSD(T) approach is demonstrated.

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

Overlap-Driven Splitting of Triplet Pairs in Singlet Fission

Here, we analyze correlated-triplet-pair (TT) singlet-fission intermediates toward two-triplet separation (T...T) using spin-state-averaged density matrix renormalization group electronic-structure calculations. Specifically, we compare the triplet–triplet exchange (J) for tetracene dimers, bipentacene, a subunit of the benzodithiophene–thiophene dioxide polymer, and a carotenoid (neurosporene). Exchange-split energy gaps of J and 3J separate a singlet from a triplet and a singlet from a quintet, respectively. We draw two new insights: (a) the canonical tetracene singlet-fission unit cell supports precisely three low-lying TT intermediates with order-of-magnitude differences in J, and (b) the separable TT intermediate in carotenoids emanates from a pair of excitations to the second triplet state. Therefore, unlike with tetracenes, carotenoid fission requires above-gap excitations. In all cases, the distinguishability of the molecular triplets—that is, the extent of orbital overlap—determines the splitting within the spin manifold of TT states. Consequently, J represents a spectroscopic observable that distnguishes the resemblance between TT intermediates and the T...T product.

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

Geometrical picture of the electron–electron correlation at the large- D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron–electron correlations and might give a bound on the correlation energy. Furthermore, the results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron–electron correlations.

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

Charge transfer as a mechanism for chlorophyll fluorescence concentration quenching

Highly concentrated solutions of chlorophyll display rapid fluorescence quenching. The same devastating energy loss is not seen in photosynthetic light-harvesting antenna complexes, despite the need for chromophores to be in close proximity to facilitate energy transfer. A promising, though unconfirmed mechanism for the observed quenching is energy transfer from an excited chlorophyll monomer to a closely associated chlorophyll pair that subsequently undergoes rapid nonradiative decay to the ground state via a short-lived intermediate charge-transfer state. In this work, we make use of newly emerging fast methods in quantum chemistry to assess the feasibility of this proposed mechanism. We calculate rate constants for the initial charge separation, based on Marcus free-energy surfaces extracted from molecular dynamics simulations of solvated chlorophyll pairs, demonstrating that this pathway will compete with fluorescence (i.e., drive quenching) at experimentally measured quenching concentrations. We show that the rate of charge separation is highly sensitive to interchlorophyll distance and the relative orientations of chromophores within a quenching pair. We discuss possible solvent effects on the rate of charge separation (and consequently the degree of quenching), using the light-harvesting complex II (LH2) protein from rps. acidophila as a specific example of how this process might be controlled in a protein environment. Crucially, we reveal that the LH2 antenna protein prevents quenching, even at the high chlorophyll concentrations required for efficient energy transfer, by restricting the range of orientations that neighboring chlorophyll pairs can adopt.

59 BASIC BIOLOGICAL SCIENCES↗