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 55 records · Page 3

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↗

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

Herein, we present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling Multireference Calculations on Multimetallic Systems with Graphic Processing Units

Modeling multimetallic systems efficiently enables faster prediction of desirable chemical properties and the design of new materials. This work describes an initial implementation for performing multireference wave function method localized active-space self-consistent field (LASSCF) calculations through the use of multiple graphics processing units (GPUs) to accelerate time-to-solution. Density fitting is leveraged to reduce memory requirements, and we demonstrate the ability to fully utilize multi-GPU compute nodes. Performance improvements of 5–10x in total application runtime were observed in LASSCF calculations for multimetallic catalyst systems up to 1200 AOs and an active space of (22e,40o) using up to four NVIDIA A100 GPUs. Furthermore, written with performance portability in mind, a comparable performance is also observed in early runs on the Aurora exascale system using Intel Max Series GPUs.

Algorithms↗

Reaction Profiles and Kinetics for Radical–Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C2H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol(-1) and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio Potential-Energy Surfaces and Electron-Spin-Exchange Cross Sections for H-O2 Interactions

Accurate quartet- and doublet-state potential-energy surfaces for the interaction of a hydrogen atom and an oxygen molecule in their ground states have been determined from an ab initio calculation using large-basis sets and the internally contracted multireference configuration interaction method. These potential surfaces have been used to calculate the H-O2 electron-spin-exchange cross section; the square root of the cross section (in a(sub 0)), not taking into account inelastic effects, can be obtained approximately from the expressions 2.390E(sup -1/6) and 5.266-0.708 log10(E) at low and high collision energies E (in E(sub h)), respectively. These functional forms, as well as the oscillatory structure of the cross section found at high energies, are expected from the nature of the interaction energy. The mean cross section (the cross section averaged over a Maxwellian velocity distribution) agrees reasonably well with the results of measurements.

Stallcop, James R.↗

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗

The structures, binding energies and vibrational frequencies of Ca3 and Ca4: An application of the CCSD(T) method

The Ca3 and Ca4 metallic clusters have been investigated using state-of-the-art ab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca3 and Ca4 have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both the e' mode of Ca3 and the t2 mode of Ca4 are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.

Lee, Timothy J.↗

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗

Photodissociation of CH2. I - Potential energy surfaces of the dissociation into CH and H

The possible photodissociation pathways of the CH2 radical are studied using ab initio multireference configuration-interaction methods, and accurate photodissociation cross sections and branching ratios for the production of CH + H and C + H2 are obtained. Potential energy surfaces were calculated using the Wuppertal-Bonn self-consistent field plus a multireference single and double-excitation configuration interaction package of programs. Two-dimensional potential energy surfaces of the ten lowest triplet states correlating with the seven lowest states of CH were calculated as functions of bond angle and one C-H bond distance, keeping the other C-H bond distance fixed at the equilibrium CH2 value.

Bearda, Robert A.↗

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Praseodymium in the formal +5 oxidation state

Praseodymium in the +5 oxidation state is a long-sought connection between lanthanide, early-transition and actinide metal redox chemistries. Unique among the lanthanide series, evidence for molecular pentavalent praseodymium species has been observed in the gas phase and noble gas matrix isolation conditions. Here we report the low-temperature synthesis and characterization of a molecular praseodymium complex in the formal +5 oxidation state, [Pr 5+ (NP t Bu 3 ) 4 ][X − ] (where t Bu = tert -butyl and X − = tetrakis(pentafluorophenyl)borate or hexafluorophosphate). Single-crystal X-ray diffraction, solution-state spectroscopic, solution magnetometric, density functional theory and multireference wavefunction-based methods indicate a highly multiconfigurational singlet ground state. Finally, an inverted ligand field drives this unique electronic structure, which establishes a critical link in understanding the bonding of high-valent metal complexes across the periodic table.

chemical bonding↗

Theoretical study of the NO beta system

A theoretical determination of the transition moment functions (TMFs) for the beta system of NO is presented. High levels of correlation treatment are required to show the changing degree of Rydberg character in the B2II with decreasing r values. The state-averaged complete-active-space self-consistent-field multireference configuration-interaction method is used for the determination. Previous lifetime measurements made with laser-induced fluorescence, varying between 2 and 0.85 microns, are discussed in terms of the calculated lifetimes for v-prime values 0-6, which vary from 2.12-1.17 microns. When larger r values are used for the transition moment function, the calculated lifetimes correlate with experimental lifetimes. The Einstein coefficients agree with experimental results, although limitations in the calibration of the spectrometer can account for systematic differences. The correlation with earlier experimental results suggests that radiative lifetimes are in the range of 1-2 microns.

Langhoff, Stephen R.↗

Theoretical study of the spectroscopy of the alkali oxides LiO, NaO, and KO

A state-averaged complete-active-space self-consistent-field multireference configuration-interaction method is presented to characterize the bound-bound emission from the CPi-2 state into the two lowest ionic states of LiO, NaO, and KO. Ab initio calculations use the experimental results obtained by Woodward et al. (1989) of the emitting state as CPi-2, but indicate that the tentative experimental band assignments are incorrect.

Langhoff, Stephen R.↗

The Low-Lying States of AlCu and AlAg

The singlet and triplet states of AlCu and AlAg below about 32 000/cm are studied using the internally contracted multireference configuration-interaction method. A more elaborate study of the X(sup 1)Sum(sup +) ground state of AlCu is undertaken using extended Gaussian basis sets, including the effect of inner-shell correlation and including a perturbational estimate of relativistic effects. Our best estimate of the spectroscopic constants (r(sub 0), DeltaG(sub 1/2), and D(sub 0)) for the X(sup 1)Sum(sup+) state with the experimental values in parentheses are: 4.416(4.420) a(sub 0), 295 (294) /cm, and 2.318 (2.315) eV. The calculations definitively assign the upper state in the observed transition at 14 892/cm to the lowest (sup 1)Prod state. The calculated spectroscopic constants and radiative lifetime for the (sup 1)Prod state are in good agreement with experiment. The calculations support the tentative assignments of Behm et al. for three band systems observed in the visible region between 25 000 and 28 000 / cm. However, the computed spectroscopic constants are in very poor agreement with those deduced from an analysis of the spectra. Analogous theoretical results for AlAg suggest that the (2)(sup 3)Prod, (3)(sup 3)Prod, and (3)(sup 1)Sum(sup +) states account for the bands observed, but not assigned, by Duncan and co-workers.

Bauschlicher, Charles W., Jr.↗

An Ab Initio Study of the Low-Lying Doublet States of AgO and AgS

Spectroscopic constants (D(sub o), r(sub e), mu(sub e), T(sub e)) are determined for the doublet states of AgO and AgS below approx. = 30000/cm. Large valence basis sets are employed in conjunction with relativistic effective core potentials (RECPs). Electron correlation is included using the modified coupled-pair functional (MCPF) and multireference configuration interaction (MRCI) methods. The A(sup 2)Sigma(sup +) - X(sup 2)Pi band system is found to occur in the near infrared (approx. = 9000/cm) and to be relatively weak with a radiative lifetime of 900 microns for A(sup 2)Sigma(sup +) (upsilon = 0). The weakly bound C(sup 2)Pi state (our notation), the upper state of the blue system, is found to require high levels of theoretical treatment to determine a quantitatively accurate potential. The red system is assigned as a transition from the C(sup 2)Pi state to the previously unobserved A(sup 2)Sigma(sup +) state. Several additional transitions are identified that should be detectable experimentally. A more limited study is performed for the vertical excitation spectrum of AgS. In addition, a detailed all-electron study of the X(sup 2)Pi and A(sup 2)Sigma(sup +) states of AgO is carried out using large atomic natural orbital (ANO) basis sets. Our best calculated D(sub o) value for AgO is significantly less than the experimental value, which suggests that there may be some systematic error in the experimental determination.

Bauschlicher, Charles W., Jr.↗

The Theoretical Transition Probabilities Between the B(sup 3)Pi(sub g) and the A(sup 3)Sigma(Sup +, sub u), W(sup 3)Delta(sub u), B'(sup 3)Sigma(sup -, sub u) States of N2

The electronic transition moment functions between the B(sup 3)Pi(sub g) and the A(sup 3)Sigma(sup +, sub u), W(sup 3)Delta(sub u), B'(sup 3)Sigma(sup -, sub u) states of N2 are studied using the internally contracted multireference configuration interaction (ICMRCI) method based upon complete active space SCF (CASSCF) reference wave-functions. The dependence of the moments on both the one and n-particle basis sets has been investigated in detail. The calculated radiative lifetimes for the vibrational levels of B(sup 3)Pi(sub g) are in excellent agreement with the most recent measurement of Euler and Pipkin (1983)

Thuemmel, Helmar T.↗