Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemistry methods”

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

Quantum error mitigation by hidden inverses protocol in superconducting quantum devices

We present a method to improve the convergence of variational algorithms based on hidden inverses (HIs) to mitigate coherent errors. In the context of error mitigation, this means replacing the hardware implementation of certain Hermitian gates with their inverses. Doing so results in noise cancellation and a more resilient quantum circuit. This approach improves performance in a variety of two-qubit error models where the noise operator also inverts with the gate inversion. We apply the mitigation scheme on superconducting quantum processors running the variational quantum eigensolver (VQE) algorithm to find the H 2 ground-state energy. When implemented on superconducting hardware we find that the mitigation scheme effectively reduces the energy fluctuations in the parameter learning path in VQE, reducing the number of iterations for a converged value. We also provide a detailed numerical simulation of VQE performance under different noise models and explore how HIs and randomized compiling affect the underlying loss landscape of the learning problem. These simulations help explain our experimental hardware outcomes, helping to connect lower-level gate performance to application-specific behavior in contrast to metrics like fidelity which often do not provide an intuitive insight into observed high level performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hilbert space multireference coupled cluster tailored by matrix product states

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants—the state universal one and the Brillouin–Wigner and Mukherjee’s state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO–LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Chemistry↗
DOE OSTI2024

Accessing metal-specific orbital interactions in C–H activation with resonant inelastic X-ray scattering

Photochemically prepared transition-metal complexes are known to be effective at cleaving the strong C–H bonds of organic molecules in room temperature solutions. There is also ample theoretical evidence that the two-way, metal to ligand (MLCT) and ligand to metal (LMCT), charge-transfer between an incoming alkane C–H group and the transition metal is the decisive interaction in the C–H activation reaction. What is missing, however, are experimental methods to directly probe these interactions in order to reveal what determines reactivity of intermediates and the rate of the reaction. Here, using quantum chemical simulations we predict and propose future time-resolved valence-to-core resonant inelastic X-ray scattering (VtC-RIXS) experiments at the transition metal L-edge as a method to provide a full account of the evolution of metal–alkane interactions during transition-metal mediated C–H activation reactions. For the model system cyclopentadienyl rhodium dicarbonyl (CpRh(CO) 2 ), we demonstrate, by simulating the VtC-RIXS signatures of key intermediates in the C–H activation pathway, how the Rh-centered valence-excited states accessible through VtC-RIXS directly reflect changes in donation and back-donation between the alkane C–H group and the transition metal as the reaction proceeds via those intermediates. We benchmark and validate our quantum chemical simulations against experimental steady-state measurements of CpRh(CO) 2 and Rh(acac)(CO)2 (where acac is acetylacetonate). Our study constitutes the first step towards establishing VtC-RIXS as a new experimental observable for probing reactivity of C–H activation reactions. More generally, the study further motivates the use of time-resolved VtC-RIXS to follow the valence electronic structure evolution along photochemical, photoinitiated and photocatalytic reactions with transition metal complexes.

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

Development of Lipid-Coated Semiconductor Nanosensors for Recording of Membrane Potential in Neurons

In the past decade, optical imaging methods have significantly improved our understanding of the information processing principles in the brain. Although many promising tools have been designed, sensors of membrane potential are lagging behind the rest. Semiconductor nanoparticles are an attractive alternative to classical voltage indicators, such as voltage-sensitive dyes and proteins. Such nanoparticles exhibit high sensitivity to external electric fields via the quantum-confine Stark effect. Here we report the development of semiconductor voltage-sensitive nanorods (vsNRs) that self-insert into the neuronal membrane. To facilitate interaction of the nanorods with the membrane, we functionalized their surface with the lipid mixture derived from brain extract. We describe a workflow to detect and process the photoluminescent signal of vsNRs after widefield time-lapse recordings. We also present data indicating that vsNRs are feasible for sensing membrane potential in neurons at a single-particle level. This shows the potential of vsNRs for the detection of neuronal activity with unprecedentedly high spatial and temporal resolution.

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

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗
DOE OSTI2023

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

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

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

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

State predictive information bottleneck

We report the ability to make sense of the massive amounts of high-dimensional data generated from molecular dynamics simulations is heavily dependent on the knowledge of a low-dimensional manifold (parameterized by a reaction coordinate or RC) that typically distinguishes between relevant metastable states, and which captures the relevant slow dynamics of interest. Methods based on machine learning and artificial intelligence have been proposed over the years to deal with learning such low-dimensional manifolds, but they are often criticized for a disconnect from more traditional and physically interpretable approaches. To deal with such concerns, in this work we propose a deep learning based state predictive information bottleneck approach to learn the RC from high-dimensional molecular simulation trajectories. We demonstrate analytically and numerically how the RC learnt in this approach is connected to the committor in chemical physics and can be used to accurately identify transition states. A crucial hyperparameter in this approach is the time delay or how far into the future the algorithm should make predictions about. Through careful comparisons for benchmark systems, we demonstrate that this hyperparameter choice gives useful control over how coarse-grained we want the metastable state classification of the system to be. We thus believe that this work represents a step forward in systematic application of deep learning based ideas to molecular simulations.

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

Error mitigation, optimization, and extrapolation on a trapped-ion testbed

Current noisy intermediate-scale quantum (NISQ) trapped-ion devices are subject to errors which can significantly impact the accuracy of calculations if left unchecked. A form of error mitigation called zero noise extrapolation (ZNE) can decrease an algorithm’s sensitivity to these errors without increasing the number of required qubits. Here we explore different methods for integrating this error mitigation technique into the Variational Quantum Eigensolver (VQE) algorithm for calculating the ground state of the HeH + molecule at 0.8 Å in the presence of experimental noise. Using the Quantum Scientific Computing Open User Testbed (QSCOUT) trapped-ion device, we test three methods of scaling noise for extrapolation: time stretching the two-qubit gates, scaling the sideband detuning parameter, and inserting two-qubit gate identity operations into the ansatz circuit. We find that time stretching and sideband detuning scaling fail to scale the noise on our particular hardware in a way that can be extrapolated to zero noise. Scaling our noise with global gate identity insertions and extrapolating after variational optimization, we achieve error suppression of 96.8%, resulting in an energy estimate within –0.004 ± 0.04 hartree of the ground state energy. This is an improvement, but still outside the chemical accuracy threshold of 0.0016 hartree. Furthermore, our results show that the efficacy of this error mitigation technique depends on choosing the correct implementation for a given device architecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Exploring the elastic and electronic properties of chromium molybdenum diboride alloys

We perform first-principles calculations to study the structural, mechanical, thermal, electronic, and magnetic properties of Cr 1-x Mo x B 2 for x = 0.25, 0.33, 0.50, 0.67 and 0.75. Based on structural search methods, we determine the ground-state structure for each concentration. The ternaries are either monoclinic (x = 0.25, 0.75) or trigonal (x = 0.33, 0.50, 0.67). The calculated mechanical properties reveal that the strength of Cr 1-x Mo x B 2 is maximized for x = 0.50. Cr 0.5 Mo 0.5 B 2 exhibits excellent mechanical properties (B = 298 GPa, Y = 558 GPa, G = 235 Gpa, ν = 0.19, H v =27 GPa), surpassing those of β-MoB 2 at a lower cost. All of these ternaries are hard alloys with Vickers hardness greater than 24 GPa. Chemical bonding analysis demonstrates that the strength of the new compounds is related to the alternating planar and buckled B-B layers, as well as the strong TM-B bonds. Finally, the enhanced strength of Cr 0.5 Mo 0.5 B 2 is a consequence of the high density of strong interlayer Cr-Mo metallic bonds around the Fermi level.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods like the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite like GAMESS is maintained and developed is considered. Future developments are briefly summarized.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

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

Using resonant energy X-ray diffraction to extract chemical order parameters in ternary semiconductors

II–IV–V 2 materials, ternary analogs to III–V materials, are emerging for their potential applications in devices such as LEDs and solar cells. Controlling cation ordering in II–IV–V 2 materials offers the potential to tune properties at nearly fixed compositions and lattice parameters. While tuning properties at a fixed lattice constant through ordering has the potential to be a powerful tool used in device fabrication, cation ordering also creates challenges with characterization and quantification of ordering. In this work, we investigate two different methods to quantify cation ordering in ZnGeP 2 thin films: a stretching parameter calculated from lattice constants , and an order parameter determined from the cation site occupancies ( S ). We use high resolution X-ray diffraction (HRXRD) to determine and resonant energy X-ray diffraction (REXD) to extract S . REXD is critical to distinguish between elements with similar Z -number ( e.g. Zn and Ge). We found that samples with a corresponding to the ordered chalcopyrite structure had only partially ordered S values. The optical absorption onset for these films occurred at lower energy than expected for fully ordered ZnGeP 2 , indicating that S is a more accurate descriptor of cation order than the stretching parameter. Since disorder is complex and can occur on many length scales, metrics for quantifying disorder should be chosen that most accurately reflect the physical properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonradiative relaxation dynamics in the [Au 25- n Ag n (SH) 18 ] -1 ( n = 1, 12, 25) thiolate-protected nanoclusters

Evaluation of the electron-nuclear dynamics and relaxation mechanisms of gold and silver nanoclusters and their alloys is important for future photocatalytic, light harvesting, and photoluminescence applications of these systems. In this work, the effect of silver doping on the nonradiative excited state relaxation dynamics of the atomically precise thiolate-protected gold nanocluster [Au 25- n Ag n (SH) 18 ] –1 ( n = 1, 12, 25) is studied theoretically. Time-dependent density functional theory is used to study excited states lying in the energy range 0.0–2.5 eV. The fewest switches surface hopping method with decoherence correction was used to investigate the dynamics of these states. The HOMO–LUMO gap increases significantly upon doping of 12 silver atoms but decreases for the pure silver nanocluster. Doped clusters show a different response for ground state population increase lifetimes and excited state population decay times in comparison to the undoped system. The ground state recovery times of the S 1 –S 6 states in the first excited peak were found to be longer for [Au 13 Ag 12 (SH) 18 ] –1 than the corresponding recovery times of other studied nanoclusters, suggesting that this partially doped nanocluster is best for preserving electrons in an excited state. The decay time constants were in the range of 2.0–20 ps for the six lowest energy excited states. Among the higher excited states, S 7 has the slowest decay time constant although it occurs more quickly than S 1 decay. Overall, these clusters follow common decay time constant trends and relaxation mechanisms due to the similarities in their electronic structures.

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

Quantum Monte Carlo and Density Functional Theory Study of Strain and Magnetism in 2D 1T-VSe 2 with Charge Density Wave States

Two-dimensional (2D) 1T-VSe 2 has prompted significant interest due to the discrepancies regarding alleged ferromagnetism (FM) at room temperature, charge density wave (CDW) states, and the interplay between the two. We employed a combined Diffusion Monte Carlo (DMC) and density functional theory (DFT) approach to accurately investigate the magnetic properties, CDW states, and their responses to strain in monolayer 1T-VSe 2 . Our calculations show the delicate competition between various phases, revealing critical insights into the relationship between their energetic and structural properties. Here, we performed classical Monte Carlo simulations informed by our DMC and DFT results and found the magnetic transition temperature (T c ) of the undistorted (non-CDW) FM phase to be 228 K and the distorted (CDW) phase to be 68 K. Additionally, we studied the response of biaxial strain on the energetic stability and magnetic properties of various phases of 2D 1T-VSe 2 and found that small amounts of strain can increase the T c , suggesting a promising route for engineering and enhancing magnetic behavior. Finally, we synthesized 1T-VSe 2 and performed Raman spectroscopy measurements, which were in close agreement with our calculated results, validating our computational approach. Our work emphasizes the role of highly accurate DMC methods in advancing the understanding of monolayer 1T-VSe 2 and provides a robust framework for future studies of 2D magnetic materials.

2D magnets↗
DOE OSTI2022

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO 2 reduction (CO 2 R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO 2 R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. Furthermore, this work illustrates one possible origin of the facet dependence of CO 2 R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO 2 R electrocatalysts can be controlled by the surface morphology.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Monitoring Redox Processes in Lithium-Ion Batteries by Laboratory-Scale Operando X-ray Emission Spectroscopy

Tracking changes in the chemical state of transition metals in alkali-ion batteries is crucial to understanding the redox chemistry during operation. X-ray absorption spectroscopy (XAS) is often used to follow the chemistry through observed changes in the chemical state and local atomic structure as a function of the state-of-charge (SoC) in batteries. In this study, we utilize an operando X-ray emission spectroscopy (XES) method to observe changes in the chemical state of active elements in batteries during operation. Operando XES and XAS were compared by using a laboratory-scale setup for four different battery systems: LiCoO 2 (LCO), Li[Ni 1/3 Co 1/3 Mn 1/3 ]O 2 (NMC111), Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 (NMC811), and LiFePO 4 (LFP) under a constant current charging the battery in 10 h (C/10 charge rate). We show that XES, despite narrower chemical shifts in comparison to XAS, allows us to fingerprint the battery SOC in real time. We further demonstrate that XES can be used to track the change in net spin of the probed atoms by analyzing changes in the emission peak shape. As a test case, the connection between net spin and the local chemical and structural environment was investigated by using XES and XAS in the case of electrochemically delithiated LCO in the range of 2–10% lithium removal.

25 ENERGY STORAGE↗
DOE OSTI2026

Nanosecond Structure of Radical Pair Intermediates from High-Frequency Quantum Oscillations: Insight into the Q A •– to Q B Electron Transfer Step in Purple Bacterial Photosynthesis

We demonstrate the validity of our approach to deduce, from the anisotropy of quantum oscillations, the geometry of short-lived radical pair intermediates in photosynthesis. A global fit of a two-dimensional W-band (94 GHz) electron paramagnetic resonance (EPR) experiment provides the same global minimum values for the geometry of the A-side radical pair P 700 •+ A 1A •− in photosystem I (PSI) as observed in a previous Q-band (34 GHz) EPR study, yet with a significantly increased convergence rate of 62%. This demonstrates that the global fit yields the correct radical pair geometry even at Q-band frequencies. With this information, we revisit our previous Q-band study of the cofactor arrangement of P 865 •+ Q A •− , the stabilized charge-separated state in purple bacterial reaction centers (RCs). Analysis of calculated two-dimensional data sets of P 865 •+ Q A •− reveals that the quantum oscillation technique is unaffected by a mirror ambiguity in disordered solids and thus can provide unambiguous solutions for all five Euler angles of the radical pair geometry. This enables us to elucidate the Q A •− to Q B electron transfer step in purple bacterial photosynthesis, the subject of controversial discussions for more than 25 years. Our results show that this electron transfer step involves a gating mechanism requiring a 60° rotation of the headgroup of Q A •− in its binding pocket.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗