Search NASA⌕ Search

SEARCH · Search NASA

Results for “multireference”

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 OSTI2024

The Localized Active Space Method with Unitary Selective Coupled Cluster

Here, we introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2 , (H 2 ) 4 , and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

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

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2025

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

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

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium–Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D e and bond length r e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15–0.19. For the single-reference methods, Kohn–Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D e and r e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core–valence terms; this results in highly accurate dissociation energies and good bond lengths. In conclusion, averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

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

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
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↗
DOE OSTI2026

Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster (ΔCC) ansatz is introduced that computes the energies and wave functions of states from any degenerate or nondegenerate Slater-determinant references with any numbers of α- and β-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Møller–Plesset perturbation (ΔMP) theory. For ionized and electron-attached references, it is a coupled-cluster Green’s function, although the present theory is convergent toward the full-configuration-interaction limits, while the Feynman–Dyson many-body Green’s function (MBGF) theory generally is not. Its single-excitation instance is a projection Hartree–Fock theory as per the Thouless theorem, which may be useful for core ionizations, high-spin states, and possibly electron affinities. Additionally, a new multireference coupled-cluster theory for a general model space is developed. This quasidegenerate coupled-cluster (QCC) theory is exactly converging, but not black-box, and intended for strong correlation. Determinant-based, general-order algorithms of ΔCC and QCC theories are implemented and compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with ΔMP and MBGF theories up to the nineteenth order. An algebraic, optimal-scaling algorithm of the ΔCC theory is computer-synthesized at the levels of single excitations (ΔCCS) and of single and double excitations (ΔCCSD). As a result, the order of performance is QCC ≈ ΔCC > EOM-CC > CI at the same order or QCC ≈ ΔCC > ΔMP > MBGF at the same cost scaling.

Hirata, So [University of Illinois at Urbana-Champ↗
DOE OSTI2024

Fragment-based initialization for quantum subspace methods

Here, we present a novel quantum-classical algorithm called LAS-QKSD for multireference systems, by combining a classical localized active space (LAS) fragment-based multireference algorithm with the quantum Krylov subspace diagonalization (QKSD) method for quantum computers. The algorithm uses wave function information from a LAS self-consistent field (LASSCF) calculation to prepare an initial state with better overlap with the target ground state than the Hartree-Fock state. This is coupled with the use of QKSD to ultimately converge to the exact energy, providing faster convergence than starting from the Hartree-Fock state. Fragmentation has the two-fold benefit of fewer configurations on the classical side of the algorithm as well as fewer state preparation gates on the quantum side. First, we compare the LAS-QKSD method to the classical LASSCF method and to QKSD with a Hartree-Fock initial state. We then examine ways to load the LASSCF wave function using direct initialization and a QKSD-motivated spectral filtering approach. Finally, using a bimetallic complex, we show that the LAS-QKSD method is an efficient alternative to highly expensive complete active space SCF (CASSCF) calculations on strongly correlated systems.

D'Cunha, Ruhee↗
DOE OSTI2023

Advanced Electronic Structure Theories for Strongly Correlated Ground and Excited States

The aim of this project was to create widely applicable, numerically robust, and systematically improvable multireference theories based on the Driven Similarity Renormalization Group (DSRG). The DSRG is a many-body formalism recently developed in our lab that maps a complex problem involving strongly and weakly interacting electrons to a simpler one in which a few electrons interact strongly via "renormalized" interactions. During this project, we focused our efforts on developing a multireference version of the DSRG (MR-DSRG) applicable to a wide range of problems in theoretical chemistry, including the computation of global ground and excited potential energy surfaces, spin-splittings in transition metal complexes and predicting the interaction of near-degenerate electronically excited states.

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

Formation of phosphorus monoxide through the $$\mathbf {P}(^{4}S)+\mathbf {O}_{\mathbf {2}}(^{3}\Sigma ^{-})\rightarrow \mathbf {O}(^{3}P)+\mathbf {PO}(^{2}\Pi )$$ reaction

Phosphorus is a key and vital element for a diverse set of important biological molecules, being indispensable for life as we know. A deeper comprehension of its role in astrochemistry and atmospheric chemistry may aid in finding answers to how this element became available on Earth. The PO molecule is one of the main reservoirs of phosphorus in the interstellar medium (ISM), and a better understanding of the mechanisms and rate coefficients for its formation in the ISM is important for modelling its abundances. Here in this work, we perform multireference configuration interaction calculations on the formation of PO via the P( 4 S) + O 2 (3 Σ - ) reaction, analyzing its potential energy surface and rate coefficients for the global reaction on both doublet and quartet states. We also perform DFT (M06-2X) and CCSD(T) calculations, in order to compare the results. We found that the OPO system possesses a high multiconfigurational character, making DFT and CCSD methodologies not suitable for its potential energy landscape calculation. The rate coefficients have been calculated using the master equation system solver (MESS) package, and the results compared to recent experimental data. It is shown that the quartet state contributes for temperatures higher than 700K. The computed rate coefficient can be described by a modified Arrhenius equation [α(T/300) β exp (-γ/T)] with α = 1.44 x 10 -12 cm 3 s -1 , β= -1.66 and γ= 704 K.

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

Molecular Geometry and Electronic Structure of Copper Corroles

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree–Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S o ) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T o ), which is approximately 5 kcal/mol higher in energy compared to the S o state for both the bare and substituted copper corroles. Notably, the planarity of the T o geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S o and T o geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D 3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

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

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast Jahn‐Teller Photoswitching in Cobalt Single‐Ion Magnets

Abstract Single‐ion magnets (SIMs) constitute the ultimate size limit in the quest for miniaturizing magnetic materials. Several bottlenecks currently hindering breakthroughs in quantum information and communication technologies could be alleviated by new generations of SIMs displaying multifunctionality. Here, ultrafast optical absorption spectroscopy and X‐ray emission spectroscopy are employed to track the photoinduced spin‐state switching of the prototypical complex [Co(terpy) 2 ] 2+ (terpy = 2,2′:6′,2″‐terpyridine) in solution phase. The combined measurements and their analysis supported by density functional theory (DFT), time‐dependent‐DFT (TD‐DFT) and multireference quantum chemistry calculations reveal that the complex undergoes a spin‐state transition from a tetragonally elongated doublet state to a tetragonally compressed quartet state on the femtosecond timescale, i.e., it sustains ultrafast Jahn‐Teller (JT) photoswitching between two different spin multiplicities. Adding new Co‐based complexes as possible contenders in the search for JT photoswitching SIMs will greatly widen the possibilities for implementing magnetic multifunctionality and eventually controlling ultrafast magnetization with optical photons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probing Multiconfigurational States by Spectroscopy: The Cerium XAS L 3 ‐edge Puzzle

Abstract The Ce L 3 edge XAS spectra of CeO 2 and cerocene [Ce(C 8 H 8 ) 2 ] were calculated with relativistic ab‐initio multireference wavefunction approaches capable of reproducing the observed spectra accurately. The study aims to resolve the decades‐long puzzle regarding the relationship between the number and relative intensities of the XAS peaks and the 4f electron occupation in the ground state (GS) versus the core‐excited states (ESs). CeO 2 and cerocene exemplify the different roles of covalent bonding and wavefunction configurational composition in the observed intensity patterns. Good agreement is found between the calculated GS 4f‐shell occupations and the value derived from XAS measurements using peak areas ( n f ). The identity of the two‐peaked Ce L 3 edge is fully rationalized from the perspective of the relaxed wavefunctions for the GS and core ESs. The states underlying the different peaks differ from each other in a surprisingly simple way that can be associated with 4f 1 vs. 4f 0 sub‐configurations. Furthermore, part of one of the cerocene spectral peaks is associated with 4f 2 sub‐configurations. The pattern therefore reveals excited states that can be interpreted in terms of Ce IV and III oxidation numbers, as long assumed, with Ce II states additionally appearing in the cerocene spectrum. While this work demonstrates the rough accuracy of the conventional approach to determining n f from Ce L 3 ‐edge XAS, limitations are highlighted in terms of the ultimate accuracy of this approach and the potential of observing new types of excited states. The need to determine the sources of n f by calculations, is stressed.

Sergentu, Dumitru‐Claudiu↗
DOE OSTI2023

Electronic structure of Li 1,2,3 +,0,– and nature of the bonding in Li 2,3 +,0,–

Abstract The current study of the small lithium molecules Li 2 +,0,− and Li 3 +,0,− focuses on the nature of the bonding in these molecules as well as their structures and energetics (bond energies, ionization energies, and electron affinities). Valence CASSCF (2s,2p) calculations incorporate nondynamical electron correlation in the calculations, while the corresponding multireference configuration interaction and coupled cluster calculations incorporate dynamical electron correlation. Treatment of nondynamical correlation is critical for properly describing the Li 2,3 +,0,− molecules as well as the Li − anion with dynamical correlation, in general, only fine‐tuning the predictions. All lithium molecules and ions are bound, with the Li 3 + and Li 2 + ions being the most strongly bound, followed by Li 3 − , Li 2 , Li 2 − and Li 3 . The minimum energy structures of Li 3 +,0,− are, respectively, an equilateral triangle, an isosceles triangle, and a linear structure. The results of SCGVB calculations are analyzed to obtain insights into the nature of the bonding in these molecules. An important finding of this work is that interstitial orbitals, a concept first put forward by McAdon and Goddard in 1985, play an essential role in the bonding of all lithium molecules considered here except for Li 2 . The interstitial orbitals found in the Li 3 +,0 molecules likely give rise to the non‐nuclear attractors/maxima observed in these molecules.

Chemistry↗