Search NASA⌕ Search

SEARCH · Search NASA

Results for “multireference 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The study of molecular spectroscopy by ab initio methods

This review illustrates the potential of theory for solving spectroscopic problems. The accuracy of approximate techniques for including electron correlation have been calibrated by comparison with full configuration-interaction calculations. Examples of the application of ab initio calculations to vibrational, rotational, and electronic spectroscopy are given. It is shown that the state-averaged, complete active space self-consistent field, multireference configuration-interaction procedure provides a good approach for treating several electronic states accurately in a common molecular orbital basis.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2021

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

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

Coupled-cluster theory for strong entanglement in nuclei

Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shake-Down Spectroscopy as State- and Site-Specific Probe of Ultrafast Chemical Dynamics

Tracking the multifarious ultrafast electronic and structural changes occurring in a molecule during a photochemical transformation is a challenging endeavor that benefits from recent experimental and computational progress in time-resolved techniques. Measurements of valence electronic states, which provide a global picture of the bonding structure of the molecule, and core electronic states, which provide insight into the local environment, traditionally require different approaches and are often studied separately. Here, we demonstrate that X-ray pulses from a seeded free-electron laser (FEL) enable the measurement of high-resolution, time-resolved X-ray photoelectron spectra (XPS) that capture weak satellite states resulting from shake-down processes in a valence-excited molecule. This approach effectively combines the advantages of both valence- and core-state investigations. We applied this method to investigate photoexcited CS2 molecules, where the role of internal conversion (IC) and intersystem crossing (ISC) in determining the predissociation dynamics is controversial. We present XPS spectra from photoexcited CS2, obtained at the FERMI FEL. High-resolution measurements, compared to the corresponding spectra obtained from accurate multireference quantum chemical calculations, reveal that shake-down satellite channels are highly sensitive to both valence electronic and geometric changes. Previous studies of the predissociation dynamics have led to uncertain assignments of the branching between singlet and triplet excited states. We derive a propensity rule that demonstrates the spin-selectivity of the shake-downs. This selectivity allows us to unequivocally assign contributions from the bright and dark singlet excited states, with populations tracked along the predissociation dynamic pathway.

Thompson, Henry J. [Univ. of Southampton (United K↗
DOE OSTI2024

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

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

Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry

The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t sub 1 vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appropriate. This diagnostic, T sub 1, is defined for use with self consistent field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T sub 1 is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of nondynamical electron correlation and is far superior to C sub 0 from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T sub 1 (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

Lee, Timothy J.↗
NASA NTRS1989

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t(sub 1) vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appopriate. This diagnostic, T(sub 1) is defined for use with self-consistent-field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T(sub 1) is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of non-dynamical electron correlation and is far superior to C(sub 0) from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T(sub 1) (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

Lee, Timothy J.↗
DOE OSTI2026

Interlayer coupling in two-dimensional MoS 2 and phosphorene bilayers: Benchmark quantum Monte Carlo study of interaction energies and quasiparticle band gaps

Using high-accuracy many-body quantum Monte Carlo (QMC) methods, we study the effect of interlayer coupling on the properties of two-dimensional freestanding bilayers (BLs) of MoS 2 and phosphorene. The properties of the two BL-materials are very different and largely determined by the interlayer interaction, which is purely van der Waals in MoS 2 and partially electronic/chemical in phosphorene, resulting in a modest layer-dependent property modulation in MoS 2 and strong modulation in phosphorene. Multireference and symmetry considerations are used to construct state-of-the-art accuracy QMC trial wave functions. We determine the quasiparticle band gaps for both materials, $Δ^{\textrm{qp}}_{Γ→\textrm{K}}$ = 2.45 ± 0.05 eV in BL-MoS 2 and $Δ^{\textrm{qp}}_{Γ→Γ}$=1.59 ±0.1eV in BL-phosphorene. In conclusion, these benchmark band-gap values make it possible to consolidate the interpretation of the widely scattered experimental and theory data.

Huang, Yongda [Slovak Academy of Sciences (SAS), B↗
DOE OSTI2021

DMRG on Top of Plane-Wave Kohn–Sham Orbitals: A Case Study of Defected Boron Nitride

In this paper, we analyze the numerical aspects of the inherent multi reference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

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

Quantum algorithm for polaritonic chemistry based on an exact ansatz

Abstract Cavity-modified chemistry uses strong light-matter interactions to modify the electronic properties of molecules in order to enable new physical phenomena such as novel reaction pathways. As cavity chemistry often involves critical regions where configurations become nearly degenerate, the ability to treat multireference problems is crucial to understanding polaritonic systems. In this Letter, we show through the use of a unitary ansatz derived from the anti-Hermitian contracted Schrödinger equation that cavity-modified systems with strong correlation, such as the deformation of rectangular H 4 coupled to a cavity mode, can be solved efficiently and accurately on a quantum device. In contrast, while our quantum algorithm can be made formally exact, classical-computing methods as well as other quantum-computing algorithms often yield answers that are both quantitatively and qualitatively incorrect. Additionally, we demonstrate the current feasibility of the algorithm on near intermediate-scale quantum hardware by computing the dissociation curve of H 2 strongly coupled to a bosonic bath.

Warren, Samuel (ORCID:0000000157134454)↗
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

Linearized Pair-Density Functional Theory with Spin–Orbit Coupling

Here, we include spin–orbit coupling (SOC) effects in linearized pair-density functional theory (L-PDFT), which is a multistate extension of multiconfiguration pair-density functional theory (MC-PDFT). Both 1-electron and 2-electron SOC integrals are computed using Breit-Pauli and Douglas–Kroll–Hess Hamiltonians in the atomic mean-field approximation. SO-L-PDFT removes the unphysical J-symmetry breaking observed in MC-PDFT. The accuracy of SO-L-PDFT is validated by calculations of zero-field splittings, fine-structure excitation energies, and low-energy excited-state spectra for a diverse group of atoms and molecules spanning the whole range of the periodic table, including atoms of groups 3, 11, and 13–17, the Ce 3+ and U 5+ ions, group 16 monohydrides, group 17 monoxides, lanthanide hexachlorides ([CeCl 6 ] 3− , [PrCl 6 ] 3− , and [NdCl 6 ] 3− ), actinyl ions ([UO 2 ] + , [NpO 2 ] 2+ ), and tricarbonatoactinyl complexes ([UO 2 (CO 3 ) 3 ] 5− , [NpO 2 (CO 3 ) 3 ] 4− ). We also compare the results to new spin–orbit-inclusive calculations by single-state and multistate multireference perturbation theory.

Hamiltonians↗
DOE OSTI2025

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗
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 OSTI2025

Bayesian model mixing with multireference energy density functional

Reliably predicting nuclear properties across the entire chart of isotopes is important for applications ranging from nuclear astrophysics to superheavy science to nuclear technology. To this day, however, all the theoretical models that can scale at the level of the chart of isotopes remain semiphenomenological. Because they are fitted locally, their predictive power can vary significantly; different versions of the same theory provide different predictions. Bayesian model mixing takes advantage of such imperfect models to build a local mixture of a set of models to make improved predictions. Earlier attempts to use Bayesian model mixing for mass table calculations relied on models treated at single-reference energy density functional level, which fail to capture some of the correlations caused by configuration mixing or the restoration of broken symmetries. In this study we have applied Bayesian model mixing techniques within a multireference energy density functional (MR-EDF) framework. We considered predictions of two-particle separation energies from particle number projection or angular momentum projection with four different energy density functionals—a total of eight different MR-EDF models. We used a hierarchical Bayesian stacking framework with a Dirichlet prior distribution over weights together with an inverse log-ratio transform to enable positive correlations between different models. We found that Bayesian model mixing provides significantly improved predictions compared to the participating models. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Bond Dissociation Energies in Heavy Element Chalcogen and Halogen Small Molecules

Thermodynamic properties including bond dissociation energies, heats of formation, and gas phase acidities for the hydrides and dimers of chalcogens and halogens: H 2 Y, HX, Y 2 and X 2 for Y = Se, Te, At and X = Br, I, and At have been predicted using the Feller-Peterson-Dixon composite correlated molecular orbital theory approach. A full 4-component CCSD(T) approach was used to calculate the spin-orbit effects on thermodynamic properties, except for Se 2 where the AoC-DHF value was used due to strong multireference effects in Se 2 for the SO calculations. The calculated results show that the At 2 bond dissociation energy is quite small, 19.5 kcal/mol, with much of the low bond energy due to spin orbit effects. H 2 Po is not predicted to be stable to dehydrogenation to Po + H 2 in terms of the free energy at 298 K. In the gas phase, HAt is predicted to be a stronger acid than H 2 SO 4 . Here, the current results provide insights into potential difficulties in the actual experimental observation of such species for heavy elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗