SEARCH · Engineering Papers
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.
NWChem is a computational chemistry software package which includes quantum chemical and molecular dynamics functionality. It aims to be scalable both in its ability to treat large problems efficiently, and in its usage of available parallel computing resources.
NWChem is a computational chemistry software package which includes quantum chemical and molecular dynamics functionality. It aims to be scalable both in its ability to treat large problems efficiently, and in its usage of available parallel computing resources.
We introduce a general setup for the analog quantum simulation of the dynamics of open quantum systems based on semiconductor quantum dots electrically connected to a chain of quantum RLC electronic circuits. The dots are chosen to be in the regime of spin-charge hybridization to enhance their sensitivity to the RLC circuits while mitigating the detrimental effects of unwanted noise. In this context, we establish an experimentally realizable map between the hybrid system and a qubit coupled to thermal harmonic environments of arbitrary complexity that enables the analog quantum simulation of open quantum systems. We assess the utility of the simulator by numerically exact emulations that indicate that the experimental setup can faithfully mimic the intended target even in the presence of its natural inherent noise. We further provide a detailed analysis of the physical requirements on the quantum dots and the RLC circuits needed to experimentally realize this proposal that indicates that the simulator can be created with existing technology. The approach can exactly capture the effects of highly structured non-Markovian quantum environments typical of photosynthesis and chemical dynamics, and offer clear potential advantages over conventional and even quantum computation. The proposal opens a general path for effective quantum dynamics simulations based on semiconductor quantum dots.
Explore the source record for details and available documents.
Chemical reactions, in which bonds break and form, are highly dynamic quantum processes. A fundamental question is whether coherence can be preserved in chemical reactions and then harnessed to generate entangled products. Here we investigated this question by studying the 2KRb → K 2 + Rb 2 reaction at 500 nanokelvins, focusing on the nuclear spin degrees of freedom. We prepared the initial nuclear spins in KRb (potassium-rubidium) in an entangled state by lowering the magnetic field to where the spin-spin interaction dominates and characterized the preserved coherence in nuclear spin wave function after the reaction. We observed an interference pattern that is consistent with full coherence at the end of the reaction, suggesting that entanglement prepared within the reactants could be redistributed through the atom-exchange process.
Solute-solvent interactions are often evaluated to understand and explain gas solubilities. This information is regularly used to interpret gas solubility and separation within ionic liquids (ILs). Here, we use several different quantum chemical and previous molecular dynamics simulation results to show that strong solute-solvent interactions neither lead to high solubility nor to high selectivity, with respect to other solutes. Instead, the solute-solvent interaction strength with different isolated IL ions is affected by the electrostatic extremes. Also, the solute-solvent interaction with different IL molecules is strongly correlated to the molecular volumes or polarizabilities (within a series of same-anion based ILs).
Isolated cobalt ions on nitrogen-doped carbon (Co-N-C) catalyze CO oxidation at temperatures as low as 196 K, but the active site and mechanism for this reaction remain elusive. In this work, steady-state CO oxidation around 273 K over Co-N-C revealed nearly first order behavior in both CO and O 2 as well as a negative apparent activation energy. Isotopic transient analysis of the reaction confirmed a rapid turnover frequency and low surface coverage of adsorbed intermediates leading to CO 2 (< 10% of the Co). Results from kinetics experiments combined with quantum chemical calculations and molecular dynamics simulations are consistent with a reaction path involving weak adsorption of CO onto Co ions followed by a low barrier for the CO-assisted activation of weakly adsorbed O 2 . Furthermore, this proposed mechanism for dioxygen activation does not involve a redox cycle with the transition metal ion and may be important in other low temperature catalytic reactions involving O 2 .
Abstract In this paper, we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models, we demonstrate how reactivity in the collective VSC regime does not exhibit altered rate behavior in equilibrium but may exhibit resonant cavity modification of reactivity when the system is explicitly out of equilibrium. Our results suggest experimental protocols that may be used to modify reactivity in the collective regime and point to features not included in the models studied, which demand further scrutiny.
In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.
Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.
Cation size primarily governs dynamic correlations in aqueous solutions. However, drawing a clear line between chaotropic and cosmotrope ions is complicated by the fact that all cations delay dynamic correlation decay.
We report nonadiabatic dynamics, which goes beyond the Born–Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.
Data from quantum mechanical calculations were generated to train a NN model for predicting phosphorescence energies. Model uses localization layer, which learns to resemble the triplet spin density without explicit training.
The calcium monofluoride (CaF) molecule has emerged as a promising candidate for precision measurements, quantum simulation, and ultracold chemistry experiments. Inelastic and reactive collisions of laser cooled CaF molecules in optical tweezers have recently been reported and collisions of cold Li atoms with CaF are of current experimental interest. In this paper, we report ab initio electronic structure and full-dimensional quantum dynamical calculations of the Li + CaF → LiF + Ca chemical reaction. The electronic structure calculations are performed using the internally contracted multi-reference configuration-interaction method with Davidson correction (MRCI + Q). An analytic fit of the interaction energies is obtained using a many-body expansion method. Furthermore, a coupled-channel quantum reactive scattering approach implemented in hyperspherical coordinates is adopted for the scattering calculations under cold conditions. Results show that the Li + CaF reaction populates several low-lying vibrational levels and many rotational levels of the product LiF molecule and that the reaction is inefficient in the 1–100 mK regime allowing sympathetic cooling of CaF by collisions with cold Li atoms.
In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.
Data from quantum mechanical calculations were generated to train a NN model for predicting phosphorescence energies. Model uses localization layer, which learns to resemble the triplet spin density without explicit training. As part of this allocation, the HIPPYNN software package was tested and openly released though the Feynman Center for Innovation.
This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.
Non-classical features like interference are already being harnessed to control the output of chemical reactions. However, quantum entanglement which is an equally enigmatic many-body quantum correlation can also be used as a powerful resource yet has eluded explicit attention. In this report, an experimental scheme under the crossed beam molecular dynamical setup, with the F + HD reaction, is proposed aiming to study the possible influence of entanglement within reactant pairs on the angular features of the product distribution. The aforesaid reaction has garnered interest recently, as an unusual horseshoe shape pattern in the product (HF) distribution was observed, which has been attributed to the coupling of spin and orbital degrees of freedom. An experimental scheme is proposed aiming to study the possible influence of entanglement on the necessity for the inclusion of such spin–orbit characteristics, under circumstances wherein the existence of entanglement and spin–orbit interaction is simultaneously detectable. The attainable results are further numerically simulated highlighting specific patterns corresponding to various possibilities. Finally, such studies if extended can provide unforeseen mechanistic insight into analogous reactions, too, from the lens of quantum information.