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.
Protonation changes the molecular architecture of the crystal structures, and thus modifies the intermolecular interactions in the studied structures.
Not Available
Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.
The quantum confinement effect has driven the development and exploration of colloidal semiconductor nanocrystals over approximately the past three decades. These semiconductor nanostructures are typically synthesized as stable colloidal dispersions in solution-phase chemical batch reactions at relatively low temperatures. Researchers have developed a wide range of reaction conditions to produce a variety of shapes, compositions, and structures with well-controlled sizes and low polydispersities. The shape can be controlled to produce confinement in one, two, or all three spatial dimensions, the latter defining a “quantum dot.” Combining two or more semiconductors within the various nanostructures produces additional optical and electronic degrees of freedom. These mixed compositional nanostructures have properties not accessible in bulk semiconductor systems. Further, researchers have fabricated novel nanoheterostructures, such as “dot-in-rods,” quantum wells, core/shell structures, giant-core/shell tetrapods where the core is one material and the arm is another, Janus structures that consist of spherical structures with well-defined internal heterostructures, nano-dumbbells, hollow-core structures, and many other novel and interesting structures. The structural variations are seemingly only limited by the creativity of the various researchers engaged in these endeavors.
Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.
We incorporate explicit, non-perturbative treatment of spin–orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin–orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br 2 and I 2 , and solid Bi.
The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.
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.
NASA has undertaken a concerted effort to generate a complete model for the chemical reactions and other collisional processes that take place in the bow shock layer and wake of spacecraft entering the atmospheres of Earth and Mars. This new model is based on first-principles physics calculations: interaction potentials between colliding species obtained from accurate quantum mechanics calculations and reaction rate coefficients and cross sections based on classical mechanics simulations of heavy particle collisions.This presentation will give an introduction to the underlying theory and computational methods used in this work and describe the range of chemical and physical processes that take place in the close vicinity of spacecraft during atmospheric entry. I will also present some of our recent results for reactions involving nitrogen, carbon monoxide and CO2.
State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space–unitary coupled cluster (LAS–UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. Here, in this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C–C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.
Quantum entangled photons provide a sensitive probe of many-body interactions and offer a unique experimental portal for quantifying many-body correlations in a material system. In this paper, we present a theoretical demonstration of how photon–photon entanglement can be generated via interactions between coupled qubits. Here, we develop a model for the scattering of an entangled pair of photons from a molecular dimer. We develop a diagrammatic theory for the scattering matrix and show that one can correlate the von Neumann entropy of the outgoing bi-photon wave function with exciton exchange and repulsion interactions. Finally, we conclude by discussing possible experimental scenarios for realizing these ideas.
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.
Ab initio quantum chemistry calculations for model molecules can be used to parameterize force fields for molecular dynamics simulations of polymers. Emphasis in our research group is on using quantum chemistry-based force fields for molecular dynamics simulations of organic polymers in the melt and glassy states, but the methodology is applicable to simulations of small molecules, multicomponent systems and solutions. Special attention is paid to deriving reliable descriptions of the non-bonded and electrostatic interactions. Several procedures have been developed for deriving and calibrating these parameters. Our force fields for aromatic polyimide simulations will be described. In this application, the intermolecular interactions are the critical factor in determining many properties of the polymer (including its color).
SparcleQC is a Python package that, given a protein:ligand complex in the Protein Data Bank (PDB) file format, can create quantum mechanics/molecular mechanics (QM/MM)-like input files for the electronic structure theory packages PSI4, QChem, and NWChem. The resulting input files include quantum mechanical representations of the ligand and a small section of the protein, surrounded by point charges that represent the rest of the protein. Creation of these QM/MM input files includes cutting and capping the QM subregion, obtaining point charges for the protein, and adjusting charges at the QM/MM boundary; and each of these tasks are automated by the software. In this article, we describe the details of SparcleQC’s procedure, show examples of the Python API, and explain additional features that are helpful in protein:ligand interaction studies. Finally, we show that SparcleQC enables automated preparation of input files for QM/MM calculations, which can return can return accurate interaction energies in minutes, while a fully quantum mechanical computation on the protein:ligand complex could take days, if it is even possible.
We report solids irradiated by energetic photons can be eroded in two modes, depending on the radiation intensity. High average, low-peak power sources, e.g., synchrotron radiation and high-order harmonics, induce desorption of the material at a low etch rate. In contrast, high-peak-power radiation from extreme ultraviolet and x-ray lasers usually causes a massive removal of the material even by a single shot. In this contribution, an effective material erosion is reported in PMMA exposed to multiple accumulated pulses generated by the free-electron x-ray-laser Linac Coherent Light Source (LCLS, tuned at a photon energy of 830 eV in this study, operated in Menlo Park at Stanford, CA, USA) at a fluence below the single-pulse ablation threshold. The effect is caused by polymer-chain scissions initiated by single photons carrying enough energy to break the C—C bounds. High efficiency of the erosion is supposed to occur due to a correlation of the single-photon effects. The subthreshold damage exhibits a nonlinear dose dependence resulting from a competition between chain scissions and cross-linking processes. The cross-linking is proven by Raman spectroscopy of the irradiated polymer. Two theoretical models of the x-ray free-electron-laser-induced erosion are suggested, which provide an excellent agreement with the experimental results.
When light interacts with matter by means of scattering and absorption, we observe the resulting color. Light also probes the symmetry of matter and the result is encoded in its polarization. In the special case of circularly-polarized light, which is especially relevant in nonlinear optics, quantum photonics, and physical chemistry, a critical dimension of symmetry is along the longitudinal direction. We examine recent advances in controlling circularly-polarized light and reveal that the commonality in these advances is in judicious control of longitudinal symmetry. In particular, in the use of high quality-factor modes in dielectric metasurfaces, the finite thickness can be used to tune the modal profile. These symmetry considerations can be applied in multiplexed optical communication schemes, deterministic control of quantum emitters, and sensitive detection of the asymmetry of small molecules.
We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.
Developing multimetallic complexes with tunable metal–metal interactions has long been a target of synthetic inorganic chemistry efforts due to the unique properties that such compounds can exhibit. However, understanding relationships between metal–metal bonding and chemical properties is challenging due to system-dependent factors that influence metal–metal and metal–ligand interactions, including ligand identity, coordination geometry, and metal–metal distance. In this work, we apply X-ray absorption and emission spectroscopy and quantum chemical calculations to describe electronic structure and bonding in a series of dicobalt complexes. The compounds with silane ligands and pseudo-octahedral coordination geometry exhibit Co–Co σ and multicentered bonding character, which we characterize from both the occupied and vacant perspectives via their contributions to the Co X-ray emission and absorption spectra, respectively. In contrast, the dicobalt complexes with a pseudotetrahedral coordination environment do not exhibit Co–Co bonding due to symmetry constraints on orbital overlap. We extend these insights to the theoretical evaluation of related dicobalt complexes to explain how ligand coordination and symmetry dictate the presence or absence of a Co–Co bond. In conclusion, this work highlights how fundamental insights into electronic structure and bonding through X-ray spectroscopy uncover important factors governing metal–metal interactions and guide the rational design of multimetallic complexes with tunable metal–metal bonds.