skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces

Authors:
; ; ORCiD logo
Publication Date:
Research Org.:
SLAC National Accelerator Lab., Menlo Park, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1369418
DOE Contract Number:
AC02-76SF00515; ACI-1429830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Physical Chemistry Letters; Journal Volume: 8; Journal Issue: 11
Country of Publication:
United States
Language:
English

Citation Formats

Snyder, James W., Parrish, Robert M., and Martínez, Todd J. α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces. United States: N. p., 2017. Web. doi:10.1021/acs.jpclett.7b00940.
Snyder, James W., Parrish, Robert M., & Martínez, Todd J. α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces. United States. doi:10.1021/acs.jpclett.7b00940.
Snyder, James W., Parrish, Robert M., and Martínez, Todd J. 2017. "α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces". United States. doi:10.1021/acs.jpclett.7b00940.
@article{osti_1369418,
title = {α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces},
author = {Snyder, James W. and Parrish, Robert M. and Martínez, Todd J.},
abstractNote = {},
doi = {10.1021/acs.jpclett.7b00940},
journal = {Journal of Physical Chemistry Letters},
number = 11,
volume = 8,
place = {United States},
year = 2017,
month = 5
}
  • (12/11)CASSCF and (12/11)CASPT2 ab initio electronic structure calculations with both the cc-pVDZ and cc-pVTZ basis sets find that there is a barrier to the very exothermic hydrogen shift that converts singlet methylnitrene, CH{sub 3}N, to methyleneimine, H{sub 2}C{double{underscore}bond}NH. These two energy minima are connected by a transition structure of C{sub s} symmetry, which is computed to lie 3.8 kcal/mol above the reactant at the (12/11)CASPT2/cc-pVTZ//(12/11)CASSCF/cc-pVTZ level of theory. The (12/11)CASSCF/cc-pVTZ value for the lowest frequency vibration in the transition structure is 854 cm{sup {minus}1}, and CASPT2 calculations concur that this a{double{underscore}prime} vibration does indeed have a positive force constant. Thus,more » there is no evidence that this geometry is actually a mountain top, rather than a transition structure, on the global potential energy surface or that a C{sub 1} pathway of lower energy connects the reactant to the product. Therefore, computational results indicate that the bands seen for singlet methylnitrene in the negative ion photoelectron spectrum of CH{sub 3}N{sup {minus}} are due to singlet methylnitrene being an energy minimum, rather than a transition state. These results also lead to the prediction that, at least in principle, singlet methylnitrene should be an observable intermediate in the formation of methyleneimine.« less
  • Low-energy (0-3 eV) ballistic electrons originated during the irradiation of biological material can interact with DNA/RNA nucleobases yielding transient-anion species which undergo decompositions. Since the discovery that these reactions can eventually lead to strand breaking of the DNA chains, great efforts have been dedicated to their study. The main fragmentation at the 0-3 eV energy range is the ejection of a hydrogen atom from the specific nitrogen positions. In the present study, the methodological approach introduced in a previous work on uracil [I. González-Ramírez et al., J. Chem. Theory Comput. 8, 2769-2776 (2012)] is employed to study the DNA canonicalmore » nucleobases fragmentations of N–H bonds induced by low-energy electrons. The approach is based on minimum energy path and linear interpolation of internal coordinates computations along the N–H dissociation channels carried out at the complete-active-space self-consistent field//complete-active-space second-order perturbation theory level. On the basis of the calculated theoretical quantities, new assignations for the adenine and cytosine anion yield curves are provided. In addition, the π{sub 1}{sup −} and π{sub 2}{sup −} states of the pyrimidine nucleobases are expected to produce the temporary anions at electron energies close to 1 and 2 eV, respectively. Finally, the present theoretical results do not allow to discard neither the dipole-bound nor the valence-bound mechanisms in the range of energies explored, suggesting that both possibilities may coexist in the experiments carried out with the isolated nucleobases.« less
  • We report reaction paths starting from N({sup 2}D) + H{sub 2}O for doublet spin states, D{sub 0} and D{sub 1}. The potential energy surfaces are explored in an automated fashion using the global reaction route mapping strategy. The critical points and reaction paths have been fully optimized at the complete active space second order perturbation theory level taking all valence electrons in the active space. In addition to direct dissociation pathways that would be dominant, three roaming processes, two roaming dissociation, and one roaming isomerization: (1) H{sub 2}ON → H–O(H)N → H–HON → NO({sup 2}Π) + H{sub 2}, (2) cis-HNOHmore » → HNO–H → H–HNO → NO + H{sub 2}, (3) H{sub 2}NO → H–HNO → HNO–H → trans-HNOH, are confirmed on the D{sub 0} surface.« less
  • Complete active space MCSCF followed by MRSDCI (multireference singles and doubles configuration interaction) calculations are carried out on the electronic states of PtH/sub 2/. Spin--orbit interaction is introduced using a relativistic configuration interaction scheme on PtH/sup +/ whose d orbital Mulliken population is close to that of the d population of PtH/sub 2/ and thus enables calculation of spin--orbit splittings for the electronic states of PtH/sub 2/. The bending potential energy surfaces of the /sup 1/A/sub 1/ and /sup 3/A/sub 1/ states are obtained. The /sup 1/A/sub 1/ surface has a bent minimum and dissociates almost without a barrier intomore » Pt(/sup 1/S/sub 0/) and H/sub 2/, while the /sup 3/A/sub 1/ state has a large (--55 kcal/mol) barrier to dissociation into Pt(/sup 3/D/sub 3/)+H/sub 2/. The ground state of PtH/sub 2/ is a bent /sup 1/A/sub 1/ state (theta = 85/sup 0/).« less
  • Complete active space MCSCF (CASSCF) and second-order configuration interaction (SOCI) calculations of low-lying electronic states (/sup 2/B/sub 2/,/sup 2/A/sub 1/) of Au/sub 3/ as well as the /sup 1/..sigma../sup +//sub g/ state of Au/sub 2/ are carried out. The bending potential energy surfaces of /sup 2/A/sub 1/ and /sup 2/B/sub 2/ states are also presented. A barrier is found in the potential energy surface of the /sup 2/A/sub 1/ state in moving from the linear to bent structure. Two nearly-degenerate structures are found for the ground state. The /sup 2/..sigma../sup +//sub u/ state arising from the linear structure with anmore » Au--Au bond length of 2.66 A is only 3.2 kcal/mol below the /sup 2/A/sub 1/ bent state. The equilibrium geometry of the /sup 2/A/sub 1/ state is an isosceles triangle with an apex angle of 54/sup 0/. The Au/sub 3/ cluster is found to be more stable than the gold dimer. The effect of d correlation is studied on Au/sub 2/ by carrying out MRSDCI (multireference singles and doubles CI) calculations on the /sup 1/..sigma../sup +//sub g/ state of Au/sub 2/ which include excitations from the d orbitals.« less