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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

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

Binding and excitations in Si x H y molecular systems using quantum Monte Carlo

In this work, we present high-accuracy correlated calculations of small SixHy molecular systems in both the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within ≈0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between 1% and 3.5% with absolute values at most ≈0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations.

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

Bringing chemical structures to life with augmented reality, machine learning, and quantum chemistry

Visualizing 3D molecular structures is crucial to understanding and predicting their chemical behavior. However, static 2D hand-drawn skeletal structures remain the preferred method of chemical communication. Here, we combine cutting-edge technologies in augmented reality (AR), machine learning, and computational chemistry to develop MolAR, an open-source mobile application for visualizing molecules in AR directly from their hand-drawn chemical structures. Users can also visualize any molecule or protein directly from its name or protein data bank ID and compute chemical properties in real time via quantum chemistry cloud computing. MolAR provides an easily accessible platform for the scientific community to visualize and interact with 3D molecular structures in an immersive and engaging way.

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

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

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

Nuclear Fuel and Pu Redox Studies from The Glenn T. Seaborg Institute at Idaho National Laboratory

The Glenn T. Seaborg Institute at Idaho National Laboratory (INL-GTSI) focuses on advancing fundamental research in the actinide sciences by providing unique opportunities to early career scientists and engineers to gain experience studying the actinide elements and their associated systems. The INL-GTSI is built from three focus areas that are based on the expertise and supporting infrastructure at INL and include solid state chemistry and physics, solution phase chemistry and physics, and forensic and isotope science. INL is the lead Laboratory for nuclear energy research and development in the U. S. and the research on nuclear fuels performed under the INL-GTSI gives good examples of solid state studies. Uranium-Molybdenum (U-Mo) alloys are leading fuel candidates for conversion of high performance research and test reactors to low-enriched fuels. During irradiation, generated fission gas accumulates into bubbles and self-organizes into a gas bubble superlattice (GBS) that effectively stores fission gases and inhibits fuel swelling. A study on the early self- organizing behavior of the GBS shows that not only grain boundaries but the interfaces between the U-Mo matrix and uranium carbide (UC) impurities are important to GBS formation.[1] In solution, understanding the complex redox behavior of plutonium in aqueous environments is critical for establishing optimized nuclear waste reprocessing solvent systems and storage tank environments. INL-GTSI researchers have produced an experimentally validated multi-scale model of the gamma radiation induced behavior of plutonium ions in concentrated aqueous HNO3 solutions.[2] Here, gamma radiation effected only minimal steady state changes in the redox distribution of the plutonium oxidation states. The redox cycling between Pu(IV) and Pu(III) is demonstrated to be mediated by the •OH/NO3• radical oxidation of Pu(III) and the H2O2/HNO3 driven reduction of Pu(IV). The INL-GTSI offers young researchers the unique chance to work directly with actinide bearing materials in a U. S. National Laboratory environment. Further topical areas of interest to the INL-GTSI include, but are not limited to, fundamental actinide properties, structure/property (electronic, magnetic, thermal) relations, actinide quantum criticality, f- electron interactions, electron correlations, computational studies, new phases, defect effects, interface interactions, isotope production and separation, forensic analytical chemistry, structure and dynamic properties of actinides in non-aqueous media, separations chemistry and kinetics for advanced nuclear fuel cycles, radiation effects, and innovative and advanced ligand design for complexation of the actinides.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Ab Initio Molecular Dynamics Spectra for Characterization of Hydrated Al 2 O 3 Supported MoO x

The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.

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

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Coulomb Interaction-Driven Entanglement of Electrons on Helium

The generation and evolution of entanglement in many-body systems is an active area of research that spans multiple fields, from quantum information science to the simulation of quantum many-body systems encountered in condensed matter, subatomic physics, and quantum chemistry. Motivated by recent experiments exploring quantum information processing systems with electrons trapped above the surface of cryogenic noble gas substrates, we theoretically investigate the generation of entanglement between two electrons via their unscreened Coulomb interaction. The model system consists of two electrons confined in separate electrostatic traps that establish microwave-frequency quantized states of their motion. We compute the motional energy spectra of the electrons, as well as their entanglement, by diagonalizing the model Hamiltonian with respect to a single-particle Hartree product basis. We also compare our results with the predictions of an effective Hamiltonian. The computational procedure outlined here can be employed for device design and guidance of experimental implementations. In particular, the theoretical tools developed here can be used for fine-tuning and optimization of control parameters in future experiments with electrons trapped above the surface of superfluid helium or solid neon. Published by the American Physical Society 2024

Physics↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

chimes_calculator

The Chebyshev Interaction Model for Efficient Simulation (ChIMES) is a machine-learned interatomic potential targeting chemistry in condensed phase systems. ChIMES models are able to approach quantum-accuracy through a systematically improvable explicitly many-bodied basis comprised of linear combinations of Chebyshev polynomials. Though originally developed as to enable description of organic molecular materials, ChIMES has successfully been applied to systems spanning ambient water to molten carbon, and leveraged as correction for density functional based tight binding simulations. The ChIMES calculator comprises a flexible toolset for evaluating ChIMES interactions (e.g. in simulations, single point calculations, etc). Users have the option of directly imbedding the ChIMES calculator within their codes, or evaluating interactions through the beginner-friendly serial interface, each of which have Python, C++, C, and Fortran API's.

Goldman, Nir↗
DOE OSTI2024

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

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

Dynamic Symmetry Breaking in Chiral Magnetic Systems

The Dzyaloshinskii–Moriya interaction (DMI) in magnetic systems stabilizes spin textures with preferred chirality, applicable to next-generation memory and computing architectures. In perpendicularly magnetized heavy-metal/ferromagnet films, the interfacial DMI originating from structural inversion asymmetry and strong spin-orbit coupling favors chiral Néel-type domain walls (DWs) whose energetics and mobility remain at issue. Thus, a new effect is characterized in which domains expand unidirectionally in response to a combination of out-of-plane and in-plane magnetic fields, with the growth direction controlled by the in-plane field strength. These growth directionalities and symmetries with applied fields cannot be understood from static treatments alone. The authors theoretically demonstrate that perpendicular field torques stabilize steady-state magnetization profiles highly asymmetric in elastic energy, resulting in a dynamic symmetry breaking consistent with the experimental findings. This phenomenon sheds light on the mechanisms governing the dynamics of Néel-type DWs and expands the utility of field-driven DW motion to probe and control chiral DWs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2024

Small tensor product distributed active space (STP-DAS) framework for relativistic and non-relativistic multiconfiguration calculations: Scaling from 10 9 on a laptop to 10 12 determinants on a supercomputer

Despite the power and flexibility of configuration interaction (CI) based methods in computational chemistry, their broader application is limited by an exponential increase in both computational and storage requirements, particularly due to the substantial memory needed for excitation lists that are crucial for scalable parallel computing. Here, the objective of this work is to develop a new CI framework, namely, the small tensor product distributed active space (STP-DAS) framework, aimed at drastically reducing memory demands for extensive CI calculations on individual workstations or laptops, while simultaneously enhancing scalability for extensive parallel computing. Moreover, the STP-DAS framework can support various CI-based techniques, such as complete active space (CAS), restricted active space, generalized active space, multireference CI, and multireference perturbation theory, applicable to both relativistic (two- and four-component) and non-relativistic theories, thus extending the utility of CI methods in computational research. We conducted benchmark studies on a supercomputer to evaluate the storage needs, parallel scalability, and communication downtime using a realistic exact-two-component CASCI (X2C-CASCI) approach, covering a range of determinants from 10 9 to 10 12 . Additionally, we performed large X2C-CASCI calculations on a single laptop and examined how the STP-DAS partitioning affects performance.

Complete-active space self-consistent field↗