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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

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

Density functional study of relaxation of adsorbate vibration modes: Dominance of anharmonic interaction

Formulation and density functional workflow for calculating the lifetime of vibrational modes of molecular adsorbates on solid surfaces due to vibration–phonon coupling are presented. The anharmonic coupling is invoked to give the correct description of the origin of temperature dependence. Using pyrrolidine (C 4 H 9 N) absorbed on the Cu(001) surface as a concrete example, here we show that the anharmonic coupling can be one to two orders more significant than the harmonic interaction for the broadening of vibrational spectra, especially as temperature increases. These results challenge the common assumption that the anharmonic interaction is weak and call for attention of considering its effect in quantum relaxation and related problems.

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

Aggregation of solutes in bosonic versus fermionic quantum fluids

Quantum fluid droplets made of helium-3 ( 3 He) or helium-4 ( 4 He) isotopes have long been considered as ideal cryogenic nanolabs, enabling unique ultracold chemistry and spectroscopy applications. The droplets were believed to provide a homogeneous environment in which dopant atoms and molecules could move and react almost as in free space but at temperatures close to absolute zero. Here, we report ultrafast x-ray diffraction experiments on xenon-doped 3 He and 4 He nanodroplets, demonstrating that the unavoidable rotational excitation of isolated droplets leads to highly anisotropic and inhomogeneous interactions between the host matrix and enclosed dopants. Superfluid 4 He droplets are laced with quantum vortices that trap the embedded particles, leading to the formation of filament-shaped clusters. In comparison, dopants in 3 He droplets gather in diffuse, ring-shaped structures along the equator. The shapes of droplets carrying filaments or rings are direct evidence that rotational excitation is the root cause for the inhomogeneous dopant distributions.

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

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Modeling Intermolecular Interactions in Nanotubes, Fullerenes and Graphite using a New Long-Range Potential

The cohesive energy and compressibility of strands of a single-wall nanotube rope has been computed using a new long-range potential energy function derived from accurate ab initio quantum chemistry calculations of the benzene dimer and calibrated for energetic and mechanical properties of graphite (at pressures up to 12 GPa). We also use this potential to calculate a variety of properties of carbon nanotubes (both single- and multi-wall) and fullerenes. Extensive comparisons are made with previously published potentials.

Jaffe, Richard↗
DOE OSTI2024

Microscopic Derivation of Transition-state Theory for Complex Quantum Systems

The decay of quantum complex systems through a potential barrier is often described with transition-state theory, also known as RRKM theory in chemistry. Here we derive the basic formula for transition-state theory based on a generic Hamiltonian as might be constructed in a configuration-interaction basis. Two reservoirs of random Hamiltonians from Gaussian orthogonal ensembles are coupled to intermediate states representing the transition states at a barrier. Under the condition that the decay of the reservoirs to open channels is large, an analytic formula for reaction rates is derived. Here, the transition states act as independent Breit–Wigner resonances which contribute additively to the total transition probability, as is well known for electronic conductance through resonant tunneling states. It is also found that the transition probability is independent of the decay properties of the states in the second reservoir over a wide range of decay widths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Binding of noble gas atoms by superhalogens

Because of their closed shells, noble gas (Ng) atoms (Ng = Ne, Ar, Kr, and Xe) seldom take part in chemical reactions, yet finding such mechanisms not only is of scientific interest but also has practical significance. Following a recent work by Mayer et al. [Proc. Natl. Acad. Sci. U. S. A. 116, 8167–8172 (2019)] on the room temperature binding of Ar to a superelectrophilic boron site embedded in a negative ion complex, B 12 (CN) 11 - , we have systematically studied the effect of cluster size and terminal ligands on the interaction of Ng by focusing on B 12 X 11 (Ng) (X = H, CN, and BO) and B 12 X 10 (Ng) 2 (X = CN and BO) whose stabilities are governed by the Wade–Mingos rule and on C 5 BX 5 (Ng) (X = H, F, and CN) and C 4 B 2 (CN) 4 (Ng) 2 whose stabilities are governed by the Huckel’s aromaticity rule. Our conclusion, based on density functional theory, is that both the cluster size and the terminal ligands matter—the interaction between the cluster and the Ng atoms becomes stronger with increasing cluster size and the electron affinity of the terminal ligands. Our studies also led to a counter-intuitive finding—removing multiple terminal ligands can enable electrophilic centers to bind multiple Ng atoms simultaneously without compromising their binding strength.

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

Coupled Cluster Downfolding Methods: the effect of double commutator terms on the accuracy of ground-state energies

Downfolding coupled cluster (CC) techniques have recently been introduced into quantum chemistry as a tool for the dimensionality reduction of the many-body quantum problem. As opposed to earlier formulations in physics and chemistry based on the concept of effective Hamiltonians, the appearance of the downfolded Hamiltonians is a natural consequence of the single-reference exponential parametrization of the wave function. In this paper, we discuss the impact of higher-order terms originating in double commutators. In analogy to previous studies, we consider the case when only one- and two-body interactions are included in the downfolded Hamiltonians. We demonstrate the efficiency of the many-body expansions involving single and double commutators for the unitary extension of the downfolded Hamiltonians on the example of the beryllium atom, and bond-breaking processes in the Li2 and H2O molecules. For the H2O system, we also analyze energies obtained with downfolding procedures as functions of the active space size.

quantum computing, quantum algorithms, electronic ↗
DOE OSTI2022

Perturbative Quantum Monte Carlo Method for Nuclear Physics

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. Here we demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to 16 O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second-order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Entangled Electrons Drive a Non-superexchange Mechanism in a Cobalt Quinoid Dimer Complex

A central theme in chemistry is the understanding of the mechanisms that drive chemical transformations. A well-known, highly cited mechanism in organometallic chemistry is the superexchange mechanism in which unpaired electrons on two or more metal centers interact through an electron pair of the bridging ligand. We use a combination of novel synthesis and computation to show that such interactions may in fact occur by a more direct mechanism than superexchange that is based on direct quantum entanglement of the two metal centers. Specifically, we synthesize and experimentally characterize a novel cobalt dimer complex with benzoquinoid bridging ligands and investigate its electronic structure with the variational two-electron reduced density matrix method using large active spaces. The result draws novel connections between inorganic mechanisms and quantum entanglement, thereby opening new possibilities for the design of strongly correlated organometallic compounds whose magnetic and spin properties have applications in superconductors, energy storage, thermoelectrics, and spintronics.

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

Trigonometric continuous-variable gates and hybrid quantum simulations of the sine-Gordon model

Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present an ancilla-based framework for implementing trigonometric gates with arguments given by arbitrary Hermitian functions of qumode quadratures. The protocol yields unitary gates deterministically, and non-unitary gates through probabilistic post-selection. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
NASA NTRS1989

Accurate quantum chemical calculations

An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2020

Understanding lead halide perovskites

It has been ten years since research on lead halide perovskites (LHPs) took off spectacularly and expanded into one of the most active fields today in chemistry, physics, materials science, and engineering. There is no doubt that the rapid growth of this field has been fueled by the exceptional performance of this class of materials in solar cells, light emitting devices, and other electronic, optoelectronic, and photonic devices. Such application oriented advance has been accompanied by major development in mechanistic understanding. Growing mechanistic studies on LHPs have revealed structure–property relationships including, among others, structural instability, unique defect chemistry, dynamic disorder and anharmonicity of phonons, local symmetry breaking, efficient screening of charge carriers in large polaron formation, and the Rashba effect due to spin–orbital coupling and the presence of local electric fields. This special issue collects the cutting edge of our understanding of why LHPs work and what lessons LHPs have taught us about structure and dynamics, in general.

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

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

We report recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry– Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗