Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONFIGURATION INTERACTION”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A new set of potential energy surfaces for HCO: Influence of Renner-Teller coupling on the bound and resonance vibrational states

It is commonly understood that the Renner-Teller effect can strongly influence the spectroscopy of molecules through coupling of electronic states. Here we investigate the vibrational bound states and low-lying resonances of the formyl radical treating the Renner-Teller coupled $ \tilde{X} $ 2 A' and à 2 A" states using the MultiConfiguration Time Dependent Hartree (MCTDH) method. The calculations were performed using the improved relaxation method for the bound states and a recently published extension to compute resonances. A new set of accurate global potential energy surfaces were computed at the explicitly correlated multireference configuration interaction (MRCI-F12) level and yielded remarkably close agreement with experiment in this application and thus enable future studies including photodissociation and collisional dynamics. Here, the results show the necessity of including the large contribution from a Davidson correction in the electronic structure calculations in order to appreciate the relatively small effect of the Renner-Teller coupling on the states considered here.

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

Calculating potential energy curves with fixed-node diffusion Monte Carlo: CO and N 2

This study reports on the prospect for the routine use of Quantum Monte Carlo (QMC) for the electronic structure problem, applying fixed-node Diffusion Monte Carlo (DMC) to generate highly accurate Born-Oppenheimer potential energy curves (PECs) for small molecular systems. The singlet ground electronic states of CO and N 2 were used as test cases. The PECs obtained by DMC employing multiconfigurational trial wavefunctions were compared with those obtained by conventional high-accuracy electronic structure methods such as multireference configuration interaction and/or the best available empirical spectroscopic curves. The goal was to test whether a straightforward procedure using available QMC codes could be applied robustly and reliably. Results obtained with DMC codes were found to be in close agreement with the benchmark PECs, and the n 3 scaling with the number of electrons (compared with n 7 or worse for conventional high-accuracy quantum chemistry) could be advantageous depending on the system size. Due to a large pre-factor in the scaling, for the small systems tested here, it is currently still much more computationally intensive to compute PECs with QMC. Nevertheless, QMC algorithms are particularly well-suited to large-scale parallelization and are therefore likely to become more relevant for future massively parallel hardware architectures.

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

Merging symmetry projection methods with coupled cluster theory: Lessons from the Lipkin model Hamiltonian

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail in different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. In conclusion, despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate in the weakly and strongly correlated limits, as well as the recoupling regime.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Driven similarity renormalization group: Third-order multireference perturbation theory

Here, a third-order multireference perturbation theory based on the driven similarity renormalization group (DSRG-MRPT3) approach is presented. The DSRG-MRPT3 method has several appealing features: (a) it is intruder free, (b) it is size consistent, (c) it leads to a non-iterative algorithm with O(N 6 ) scaling, and (d) it includes reference relaxation effects. The DSRG-MRPT3 scheme is benchmarked on the potential energy curves of F 2 , H 2 O 2 , C 2 H 6 , and N 2 along the F–F, O–O, C–C, and N–N bond dissociation coordinates, respectively. The nonparallelism errors of DSRG-MRPT3 are consistent with those of complete active space third-order perturbation theory and multireference configuration interaction with singles and doubles and show significant improvements over those obtained from DSRG second-order multireference perturbation theory. Our efficient implementation of the DSRG-MRPT3 based on factorized electron repulsion integrals enables studies of medium-sized open-shell organic compounds. This point is demonstrated with computations of the singlet-triplet splitting (Δ ST = E T –E S ) of 9,10-anthracyne. At the DSRG-MRPT3 level of theory, our best estimate of the adiabatic Δ ST is 3.9 kcal mol –1 , a value that is within 0.1 kcal mol –1 from multireference coupled cluster results.

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

Ortho-para interconversion in cation-water complexes: The case of V + (H 2 O) and Nb + (H 2 O) clusters

Vanadium and niobium cation-water complexes, V + (H 2 O) and Nb + (H 2 O), are produced by laser vaporization in a pulsed supersonic expansion, mass selected in a time-of-flight spectrometer, and studied with infrared photodissociation spectroscopy using rare gas atom (Ar, Ne) complex predissociation. The vibrational bands measured in the O–H stretching region contain K-type rotational sub-band structure, which provides insight into the structures of these complexes. However, rotational sub-bands do not exhibit the simple patterns seen previously for other metal ion-water complexes. We report the $A$ rotational constants are smaller than expected and the normal 1:3 intensity ratios for K = even:odd levels for independent ortho:para nuclear spin states are missing for some complexes. We relied on highly correlated internally contracted Multi-Reference Configuration Interaction (icMRCI) and Coupled Cluster [CCSD(T)] electronic structure calculations of those complexes with and without the rare gas atoms to investigate these anomalies. Rare gas atoms were found to bind via asymmetric motifs to the hydrated complexes undergoing large amplitude motions that vibrationally average to quasi-C 2v symmetry with significant probability off the C 2 axis, thus explaining the reduced $A$ values. Both vanadium and iobium cations exhibit unusually strong nuclear spin coupling to the hydrogen atoms of water, the values of which vary with their electronic state. This catalyzes ortho-para interconversion in some complexes and explains the rotational patterns. The rate of ortho-para relaxation in the equilibrated complexes must therefore be greater than the collisional cooling rate in the supersonic expansion (about 106 sec-1).

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

Correlation consistent basis sets for actinides. II. The atoms Ac and Np–Lr

New correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are reported for the actinide elements Ac and Np through Lr. These complete the series of sets reported previously for Th–U [K. A. Peterson, J. Chem. Phys. 142, 074105 (2015); M. Vasiliu et al., J. Phys. Chem. A 119, 11422 (2015)]. The new sets range in size from double- to quadruple-zeta and encompass both those optimized for valence (6s6p5f7s6d) and outer-core electron correlations (valence + 5s5p5d). The final sets have been contracted for both the DKH3 and eXact 2-component (X2C) Hamiltonians, yielding cc-pVnZ-DK3/cc-pVnZ-X2C sets for valence correlation and cc-pwCVnZ-DK3/cc-pwCVnZ-X2C sets for outer-core correlation (n = D, T, Q in each case). In order to test the effectiveness of the new basis sets, both atomic and molecular benchmark calculations have been carried out. In the first case, the first three atomic ionization potentials (IPs) of all the actinide elements Ac–Lr have been calculated using the Feller-Peterson-Dixon (FPD) composite approach, primarily with the multireference configuration interaction (MRCI) method. Excellent convergence towards the respective complete basis set (CBS) limits is achieved with the new sets, leading to good agreement with experiment, where these exist, after accurately accounting for spin-orbit effects using the 4-component Dirac-Hartree-Fock method. For a molecular test, the IP and atomization energy (AE) of PuO 2 have been calculated also using the FPD method but using a coupled cluster approach with spin-orbit coupling accounted for using the 4-component MRCI. The present calculations yield an IP 0 for PuO 2 of 159.8 kcal/mol, which is in excellent agreement with the experimental electron transfer bracketing value of 162 ± 3 kcal/mol. Likewise, the calculated 0 K AE of 305.6 kcal/mol is in very good agreement with the currently accepted experimental value of 303.1 ± 5 kcal/mol. The ground state of PuO 2 is predicted to be the 5 Σ$+\atop{0g}$ state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

One-particle many-body Green’s function theory: Algebraic recursive definitions, linked-diagram theorem, irreducible-diagram theorem, and general-order algorithms

A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green’s function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green’s function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green’s function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green’s function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green’s function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams

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

Charge-transfer excited states: Seeking a balanced and efficient wave function ansatz in variational Monte Carlo

We present a simple and efficient wave function ansatz for the treatment of excited charge-transfer states in real-space quantum Monte Carlo methods. Using the recently-introduced variation-after-response method, this ansatz allows a crucial orbital optimization step to be performed beyond a configuration interaction singles expansion, while only requiring calculation of two Slater determinant objects. As a result, we demonstrate this ansatz for the illustrative example of the stretched LiF molecule, for a range of excited states of formaldehyde, and finally for the more challenging ethylene-tetrafluoroethylene molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Theoretical investigation of intersystem crossing in the cyanonitrene molecule, 1 NCN → 3 NCN

The NCN diradical is an important intermediate of prompt nitric oxide formation in flames. The mechanism of intersystem crossing (ISC) in the NCN molecule formed via pyrolysis or photolysis of NCN 3 is of relevance to the interpretation of experiments that utilize NCN 3 as a precursor for laboratory studies of NCN kinetics. This mechanism has been investigated by means of multi-reference configuration interaction calculations. From the potential energy surfaces for NCN 3 dissociation, it was inferred that both thermal and photo-chemical decomposition initially lead to NCN in its lowest singlet state, $\tilde{a}^{-1}$$Δ_g$, with a possible contribution from the $\tilde{b}^{-1}\Sigma_g^+$ state at low photolysis wavelengths. Direct formation of the triplet ground state $\tilde{X}^{-3}\Sigma_g^-$ is also feasible for the photolytic pathway. Ananalysis of surface crossings between $\tilde{a}$ or $\tilde{b}$ and the triplet ground state $\tilde{X}^{-3}\Sigma_g^-$ in the absence and presence of a helium atom revealed an ISC channel 1 NCN ($\tilde{a}$) → 3 NCN($\tilde{X}$) via a strongly bent structure. However, its barrier of 38 kcal mol -1 relative to the singlet minimum turned out to be much too high to explain the fast ISC observed in experiments. A rigid-bender model including Renner-Teller interactions was used to examine the occurrence of mixed-multiplicity rovibrational states-so-called gateway states-that could enhance collision-induced ISC. The results of this study indicate that a gateway mechanism is probably not operative in the case of the $\tilde{a}$/$\tilde{X}$pair of states in NCN.

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

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods

In this work, we report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states’ character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. Finally, the numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

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

Guided ion beam and theoretical studies of the bond energy of SmS +

Previous work has shown that atomic samarium cations react with carbonyl sulfide to form SmS + + CO in an exothermic and barrierless process. To characterize this reaction further, the bond energy of SmS + is determined in the present study using guided ion beam tandem mass spectrometry. Reactions of SmS + with Xe, CO, and O 2 are examined. Results for collision-induced dissociation processes with all three molecules along with the endothermicity of the SmS + + CO → Sm + + COS exchange reaction are combined to yield D 0 (Sm + –S) = 3.37 ± 0.20 eV. The CO and O 2 reactions also yield a SmSO + product, with measured endothermicities that indicate D 0 (SSm + –O) = 3.73 ± 0.16 eV and D 0 (OSm + –S) = 1.38 ± 0.27 eV. The SmS + bond energy is compared with theoretical values characterized at several levels of theory, including CCSD(T) complete basis set extrapolations using all-electron basis sets. Multireference configuration interaction calculations with explicit spin-orbit calculations along with composite thermochemistry using the Feller-Peterson-Dixon method and all-electron basis sets were also explored for SmS+, and for comparison, SmO, SmO + , and EuO

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

Driven similarity renormalization group for excited states: A state-averaged perturbation theory

The multireference driven similarity renormalization group (MRDSRG) approach [C. Li and F. A. Evangelista, J. Chem. Theory Comput. 11, 2097 (2015)] is generalized to treat quasi-degenerate electronic excited states. The new scheme, termed state-averaged (SA) MRDSRG, is a state-universal approach that considers an ensemble of quasi-degenerate states on an equal footing. Using the SA-MRDSRG framework, we implement second- (SA-DSRG-PT2) and third-order (SA-DSRG-PT3) perturbation theories. These perturbation theories can treat a manifold of near-degenerate states at the cost of a single state-specific computation. At the same time, they have several desirable properties: (1) they are intruder-free and size-extensive, (2) their energy expressions can be evaluated non-iteratively and require at most the three-body density cumulant of the reference states, and (3) the reference states are allowed to relax in the presence of dynamical correlation effects. Numerical benchmarks on the potential energy surfaces of lithium fluoride, ammonia, and the penta-2,4-dieniminium cation reveal that the SA-DSRG-PT2 method yields results with accuracy similar to that of other second-order quasi-degenerate perturbation theories. The SA-DSRG-PT3 results are instead consistent with those from multireference configuration interaction with singles and doubles (MRCISD). Finally, we compute the vertical excitation energies of (E,E)-1,3,5,7-octatetraene. As a result, the ordering of the lowest three states is predicted to be 2 1 A – g <1 1 B + u <1 1 B – u by both SA-DSRG-PT2 and SA-DSRG-PT3, in accordance with MRCISD plus Davidson correction.

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

Normal and resonant Auger spectroscopy of isocyanic acid, HNCO

In this study, we investigate HNCO by resonant and nonresonant Auger electron spectroscopy at the K-edges of carbon, nitrogen, and oxygen, employing soft X-ray synchrotron radiation. In comparison with the isosteric but linear CO 2 molecule, spectra of the bent HNCO molecule are similar but more complex due to its reduced symmetry, wherein the degeneracy of the π-orbitals is lifted. Resonant Auger electron spectra are presented at different photon energies over the first core-excited 1s → 10a' resonance. All Auger electron spectra are assigned based on ab initio configuration interaction computations combined with the one-center approximation for Auger intensities and moment theory to consider vibrational motion. The calculated spectra were scaled by a newly introduced energy scaling factor, and generally, good agreement is found between experiment and theory for normal as well as resonant Auger electron spectra. A comparison of resonant Auger spectra with nonresonant Auger structures shows a slight broadening as well as a shift of the former spectra between –8 and –9 eV due to the spectating electron. Since HNCO is a small molecule and contains the four most abundant atoms of organic molecules, the reported Auger electron decay spectra will provide a benchmark for further theoretical approaches in the computation of core electron spectra.

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

Perspective: Multireference coupled cluster theories of dynamical electron correlation

Predicting the electronic structure and properties of molecular systems that display strong electron correlation effects continues to remain a fundamental theoretical challenge. This perspective investigates the recent progress and current challenges in multireference wave function methods for dynamical electron correlation, honing in on systematically improvable methods that go beyond the limitations of configuration interaction and perturbation theory

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

Reinterpretation of the electronic absorption spectrum of the methylene amidogen radical (H 2 CN)

The odd electronic absorption spectrum of H 2 CN has been of great interest to experiment. Herein, this system is reported on extensively by applying theoretical methods to the ground and low-lying excited electronic states. Employing a large breadth of high-level ab initio computations, including coupled cluster [CCSD(T) and CCSDT(Q)] and multireference configuration interaction [MRCISD+Q] methods, we comprehensively demonstrate that the most recent experimental and theoretical interpretations of the electronic spectrum of H 2 CN are in error. The previous assignments of the two broad features in the spectrum as the origin 0$^{0}_{0}$ (~35 050 cm –1 ) and 4$^{2}_{0}$ (~35 600 cm –1 )$\widetilde{B}$ 2 A 1 ←$\widetilde{X}$ 2 B 2 transitions are both found to be incorrect. The presently reported transition energies suggest that the higher energy band near 35 600 cm –1 is the true origin band. Moreover, from the computed anharmonic vibrational frequencies of the $\widetilde{X}$ 2 B 2 and $\widetilde{B}$ 2 A 1 states, we demonstrate that this ~550 cm –1 band spacing cannot be attributed to a simple vibronic transition, as claimed by the 4$^{2}_{0}$ assignment. Possible alternative explanations for the appearance of the lower intensity band near 35 050 cm –1 are discussed.

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

Communication: A mean field platform for excited state quantum chemistry

We present a mean field theory for excited states that is broadly analogous to ground state Hartree-Fock theory. Like Hartree-Fock, our approach is deterministic, state-specific, applies a variational principle to a minimally correlated ansatz, produces energy stationary points, relaxes the orbital basis, has a Fock-build cost-scaling, and can serve as the foundation for correlation methods such as perturbation theory and coupled cluster theory. To emphasize this last point, we pair our mean field approach with an excited state analog of second order Møller-Plesset theory and demonstrate that in water, formaldehyde, neon, and stretched lithium fluoride, the resulting accuracy far exceeds that of configuration interaction singles and rivals that of equation of motion coupled cluster.

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

Cluster perturbation theory. III. Perturbation series for coupled cluster singles and doubles excitation energies

The cluster perturbation series, CPS(D), for coupled cluster singles and doubles excitation energies is considered. It is demonstrated that the second-order model CPS(D-2) is identical to the configuration interaction singles with perturbative doubles, CIS(D) model. The third-order model, CPS(D-3), provides excitation energies of coupled cluster singles and doubles (CCSD) quality in the sense that the difference between CPS(D-3) and CCSD excitation energies is of the same size or smaller than the effect of adding triples corrections to CCSD excitation energies. We further show that the third-order corrections can be efficiently implemented, in particular, when the resolution of the identity approximation is used for integrals. Here, we also show that the CPS(D-3) excitation energies can be determined for system sizes that are far beyond what can be considered in conventional CCSD excitation energy calculations.

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

Bond dissociation energy and electronic spectroscopy of Cr + (NH 3 ) and its isotopomers

The electronic spectra of Cr + (NH 3 ), Cr + (ND 3 ), and Cr + ( 15 NH 3 ) have been measured from 14 200 to 17 400 cm -1 using photodissociation spectroscopy. Transitions are predominantly observed from the 6 A 1 ground state, in which the Cr + has a 3$d$ 5 electronic configuration, to the $|tilde{B}$ 6 E ($Π$) state (3$d$ 4 4s). There is extensive vibronic structure in the spectrum due to a long progression in the Cr–N stretch and transitions to all six spin-orbit levels in the upper state. The spin-orbit splitting in the excited state is observed to be A so $'$ = 39 cm -1 . For the lowest spin-orbit level, the Cr–N stretching frequency in the excited state is 343 cm -1 , with an anharmonicity of 4.2 cm -1 . The 6 E ($Π$) origin is predicted to lie at T 0 = 14 697 cm -1 . The first peak observed is due to v$'$ = 1, so the observed photodissociation onset is thermodynamic rather than spectroscopic, giving D 0 (Cr + -NH 3 ) = 14 830 ± 100 cm -1 (177.4 ± 1.2 kJ/mol) and D 0 (Cr + –ND 3 ) = 15 040 ± 30 cm -1 (179.9 ± 0.4 kJ/mol). The 6 E ($Π$) state of Cr + (NH 3 ) is ~2740 cm -1 less strongly bound than the ground state, and the Cr–N bond length increases by 0.23 ± 0.03 Å upon electronic excitation. Calculations at the time-dependent density functional theory (M06) and equations of motion coupled cluster, with single and double excitations (EOM-CCSD) level fairly accurately predict the energy and vibrational frequency of the excited state. Multi-reference configuration interaction calculations show how the spin-orbit states of Cr + (NH 3 ) evolve into those of Cr + + NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗