Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive quantum 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.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

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

Excitons in two-dimensional metal-halide perovskites

Georgia Tech School of Chemistry and Biochemistry Special Summer Seminar Series, Prof Carlos Silva; While polarons --- charges bound to a lattice deformation induced by electron-phonon coupling --- are primary photoexcitations at room temperature in bulk metal-halide hybrid organic-inorganic perovskites (HOIP), excitons --- Coulomb-bound electron-hole pairs --- are the stable quasi-particles in their two-dimensional (2D) analogues. Here we address the fundamental question: are polaronic effects consequential for excitons in 2D-HIOPs? Based on our recent work, we argue that polaronic effects are manifested intrinsically in the exciton spectral structure, which is comprised of multiple non-degenerate resonances with constant inter-peak energy spacing. We highlight measurements of population and dephasing dynamics that point to the apparently deterministic role of polaronic effects in excitonic properties. We contend that an interplay of long-range and short-range exciton-lattice couplings give rise to exciton polarons, a character that fundamentally establishes their effective mass and radius, and consequently, their quantum dynamics. Given this complexity, a fundamentally far-reaching issue is how Coulomb-mediated many-body interactions---elastic scattering such as excitation-induced dephasing, inelastic exciton bimolecular scattering, and multi-exciton binding---depend upon the specific exciton-lattice coupling within the structured excitation lineshape. We measure the intrinsic and density-dependent exciton dephasing rates of the multiple excitons and their dependence on temperature by means of two-dimensional coherent excitation spectroscopy. We find that diverse excitons display distinct intrinsic dephasing rates mediated by phonon scattering involving different effective phonons, and contrasting rates of exciton-exciton elastic scattering. These findings establish specifically the consequence of distinct lattice dressing on exciton many-body quantum dynamics, which critically define fundamental optical properties that underpin photonics and quantum optoelectronics.

solar cells, metal-halide perovskites, excitons, p↗
DOE OSTI2025

Alkali Metal Control of Triplet-Mediated C–H Activation in Iron-Mediated Coupling of Dinitrogen and Benzene

Direct coupling of N 2 with abundant feedstocks like benzene to form N-containing organic compounds is a promising strategy for N 2 fixation pathways. The challenge of coupling N 2 activation and C–H bond oxidative addition was recently solved by introducing a reversible benzene C–H bond activation process mediated by a low-valent Fe(0) complex, which gave an organometallic product that could couple with partially reduced N 2 . Interestingly, the energetics of the C–H oxidative addition/reductive elimination step depends on the choice of alkali metal. However, the reason why the alkali metal influences the C–H bond activation remained elusive. In this work, we present a comprehensive study on this Fe(0)-mediated reversible C–H activation. Through density functional theory combined with high-level coupled cluster calculations, we discovered that the intermediate-spin triplet (S = 1) controls the energy of the transition state for C–H cleavage, while the high-spin quintet (S = 2) controls the position of the equilibrium. Na + drives the equilibrium toward oxidative addition due to an electrostatic effect, while K + and Rb + are dominated by a steric effect that favors the iron(0) species. The key role played by nonbonding interactions in the Fe-mediated C–H activation provides a conceptual model for alkali control over organometallic transformations.

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

Shining light on the microscopic resonant mechanism responsible for cavity-mediated chemical reactivity

Strong light–matter interaction in cavity environments is emerging as a promising approach to control chemical reactions in a non-intrusive and efficient manner. The underlying mechanism that distinguishes between steering, accelerating, or decelerating a chemical reaction has, however, remained unclear, hampering progress in this frontier area of research. We leverage quantum-electrodynamical density-functional theory to unveil the microscopic mechanism behind the experimentally observed reduced reaction rate under cavity induced resonant vibrational strong light-matter coupling. We observe multiple resonances and obtain the thus far theoretically elusive but experimentally critical resonant feature for a single strongly coupled molecule undergoing the reaction. While we describe only a single mode and do not explicitly account for collective coupling or intermolecular interactions, the qualitative agreement with experimental measurements suggests that our conclusions can be largely abstracted towards the experimental realization. Specifically, we find that the cavity mode acts as mediator between different vibrational modes. In effect, vibrational energy localized in single bonds that are critical for the reaction is redistributed differently which ultimately inhibits the reaction.

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

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

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

Crossover from hydrogen to chemical bonding

Hydrogen bonds (H-bonds) can be interpreted as a classical electrostatic interaction or as a covalent chemical bond if the interaction is strong enough. As a result, short strong H-bonds exist at an intersection between qualitatively different bonding descriptions, with few experimental methods to understand this dichotomy. The [F-H-F] − ion represents a bare short H-bond, whose distinctive vibrational potential in water is revealed with femtosecond two-dimensional infrared spectroscopy. It shows the superharmonic behavior of the proton motion, which is strongly coupled to the donor-acceptor stretching and disappears on H-bond bending. In combination with high-level quantum-chemical calculations, we demonstrate a distinct crossover in spectroscopic properties from conventional to short strong H-bonds, which identify where hydrogen bonding ends and chemical bonding begins.

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

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

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

Imaging the spin chirality of ferrimagnetic Néel skyrmions stabilized on topological antiferromagnetic Mn 3 Sn

Néel skyrmions are generally realized in asymmetric multilayers made of heavy metals (HMs) and ultrathin ferromagnets possessing strong interfacial Dzyaloshinskii-Moriya interactions (iDMIs). Depending on the relative strengths of iDMIs at the interfaces, various types of Néel skyrmions have been suggested, which are typified with characteristically different topological properties and current-driven dynamics. This suggests the importance of a precise quantification of their spin chiralities. In this paper, we explore the possibility of realizing Néel skyrmions in magnetic multilayers without the direct usage of standard HMs. Specifically, through depositing a thin layer of ferrimagnetic (FIM) CoTb layer on top of an antiferromagnetic (AFM) quantum material of composition Mn 3 Sn, the AFM exchange interaction at the asymmetric interface provides an equivalent iDMI for stabilizing FIM Néel skyrmions. Secondly, through using advanced four-dimensional Lorentz scanning transmission electron microscopy (4D LSTEM), in combination with x-ray magnetic circular dichroism photoemission electron microscopy (XMCD-PEEM), we can directly determine the spin chirality of FIM Néel skyrmions. Finally, the present findings not only broaden the phase space for chiral interfacial magnetism but also provide a possibility for future applications of heavy-metal-free skyrmionic devices.

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

Colossal Terahertz Emission with Ultrafast Tunability Based on Van der Waals Ferroelectric NbOI 2

Abstract Terahertz (THz) technology is critical for quantum material physics, biomedical imaging, ultrafast electronics, and next‐generation wireless communications. However, standing in the way of widespread applications is the scarcity of efficient ultrafast THz sources with on‐demand fast modulation and easy on‐chip integration capability. Here the discovery of colossal THz emission is reported from a van der Waals (vdW) ferroelectric semiconductor NbOI 2 . Using THz emission spectroscopy, a THz generation efficiency an order of magnitude higher than that of ZnTe, a standard nonlinear crystal for ultrafast THz generation is observed. The underlying generation mechanisms associated are further uncovered with its large ferroelectric polarization by studying the THz emission dependence on excitation wavelength, incident polarization, and fluence. Moreover, the ultrafast coherent amplification and annihilation of the THz emission and associated coherent phonon oscillations by employing a double‐pump scheme are demonstrated. These findings combined with first‐principles calculations, inform a new understanding of the THz light–matter interaction in emergent vdW ferroelectrics and pave the way to develop high‐performance THz devices on them for quantum materials sensing and ultrafast electronics.

Subedi, Sujan [Department of Materials Science and↗
DOE OSTI2022

Mesoscale interplay between phonons and crystal electric field excitations in quantum spin liquid candidate CsYbSe 2

In this report CsYbSe 2 , a recently identified quantum spin liquid (QSL) candidate, exhibits strong crystal electric field (CEF) excitations. Here, we identify phonon and CEF modes with Raman spectroscopy and observe strong CEF-phonon coupling resulting in a vibronic bound state. Complex, mesoscale interplay between phonon modes and CEF modes is observed in real space. Additionally, an unexpected resonant Raman excitation condition is satisfied, yielding up to third-order combination modes, with a total of 17 modes identified in the spectra. This study paves the way to coherent control of possible QSL ground states with optically accessible CEF-phonon manifolds and mesoscale engineering of CEF-phonon interactions.

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

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. In conclusion, our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Prediction of correlation energies using variational subspace valence bond

In the variational subspace valence bond (VSVB) [G. D. Fletcher, J. Chem. Phys. 142, 134112 (2015)] method, the electronic orbitals comprising the wave function correspond to chemically meaningful objects, such as bonds, lone pairs, atomic cores, and so on. Selected regions of a molecule (for example, a single chemical bond, as opposed to all the bonds) can be modeled with different levels of basis set and possible methods for modeling correlation from the other regions. The interactions between the components of a molecule (say, a bond and a neighboring orbital) can then be studied in detail for their impact on a chemical phenomenon while avoiding the expense of necessarily applying the higher levels and methods to the entire molecule. Here, this work presents the theoretical basis for modeling correlation effects between specific electron pairs by incorporating terms in the inter-electronic coordinates (“r 12 ”) into VSVB. The approach is validated with calculations on small systems using single-reference wave functions.

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

Dependency of solvation effects on metal identity in surface reactions

Solvent interactions with adsorbed moieties involved in surface reactions are often believed to be similar for different metal surfaces. However, solvents alter the electronic structures of surface atoms, which in turn affects their interaction with adsorbed moieties. To reveal the importance of metal identity on aqueous solvent effects in heterogeneous catalysis, we studied solvent effects on the activation free energies of the O–H and C–H bond cleavages of ethylene glycol over the (111) facet of six transition metals (Ni, Pd, Pt, Cu, Ag, Au) using an explicit solvation approach based on a hybrid quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface. A significant metal dependence on aqueous solvation effects was observed that suggests solvation effects must be studied in detail for every reaction system. The main reason for this dependence could be traced back to a different amount of charge-transfer between the adsorbed moieties and metals in the reactant and transition states for the different metal surfaces.

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

Melting of charge order in the low-temperature state of an electronic ferroelectric-like system

Strong electronic interactions can drive a system into a state with a symmetry breaking. Lattice frustration or competing interactions tend to prevent symmetry breaking, leading to quantum disordered phases. In spin systems frustration can produce a spin liquid state. Frustration of a charge degree of freedom also can result in various exotic states, however, experimental data on these effects is scarce. In this work we demonstrate how in a Mott insulator on a weakly anisotropic triangular lattice a charge ordered state melts on cooling down to low temperatures. Raman scattering spectroscopy finds that κ -(BEDT-TTF) 2 Hg(SCN) 2 Cl enters an insulating “dipole solid” state at T =30K, but below T =15K the order melts, while preserving the insulating energy gap. Based on these observations, we suggest a phase diagram relevant to other quantum paraelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

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

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na 2 + ) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that for the Na 2 + molecule, solute--solvent interactions lead to two stable, chemically-distinct coordination states (Na(THF) 4 -Na(THF) 5 + and Na(THF) 5 -Na(THF) 5 + ) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Lastly, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

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

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗