The accuracy of molecular bond lengths computed by multireference electronic structure methods

Ron Shepard, Gary S. Kedziora, Hans Lischka, Isaiah Shavitt, Thomas Müller, P. Szalay, M. Kállay, Michael Seth

Research output: Contribution to journalArticle

21 Citations (Scopus)

Abstract

We compare experimental Re values with computed Re values for 20 molecules using three multireference electronic structure methods, MCSCF, MR-SDCI, and MR-AQCC. Three correlation-consistent orbital basis sets are used, along with complete basis set extrapolations, for all of the molecules. These data complement those computed previously with single-reference methods. Several trends are observed. The SCF Re values tend to be shorter than the experimental values, and the MCSCF values tend to be longer than the experimental values. We attribute these trends to the ionic contamination of the SCF wave function and to the corresponding systematic distortion of the potential energy curve. For the individual bonds, the MR-SDCI Re values tend to be shorter than the MR-AQCC values, which in turn tend to be shorter than the MCSCF values. Compared to the previous single-reference results, the MCSCF values are roughly comparable to the MP4 and CCSD methods, which are more accurate than might be expected due to the fact that these MCSCF wave functions include no extra-valence electron correlation effects. This suggests that static valence correlation effects, such as near-degeneracies and the ability to dissociate correctly to neutral fragments, play an important role in determining the shape of the potential energy surface, even near equilibrium structures. The MR-SDCI and MR-AQCC methods predict Re values with an accuracy comparable to, or better than, the best single-reference methods (MP4, CCSD, and CCSD(T)), despite the fact that triple and higher excitations into the extra-valence orbital space are included in the single-reference methods but are absent in the multireference wave functions. The computed Re values using the multireference methods tend to be smooth and monotonic with basis set improvement. The molecular structures are optimized using analytic energy gradients, and the timings for these calculations show the practical advantage of using variational wave functions for which the Hellmann-Feynman theorem can be exploited.

Original languageEnglish
Pages (from-to)37-57
Number of pages21
JournalChemical Physics
Volume349
Issue number1-3
DOIs
Publication statusPublished - Jun 16 2008

Fingerprint

Bond length
chemical bonds
Wave functions
Electronic structure
wave functions
electronic structure
valence
self consistent fields
potential energy
Hellmann-Feynman theorem
trends
Electron correlations
orbitals
Potential energy surfaces
Molecules
Potential energy
Extrapolation
complement
Molecular structure
extrapolation

Keywords

  • Bond length
  • CBS
  • Complete basis set limit
  • Configuration interaction
  • GVB
  • MCSCF
  • MR-AQCC
  • MR-SDCI
  • Multiconfiguration self-consistent field
  • Multiference
  • R

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Spectroscopy
  • Atomic and Molecular Physics, and Optics

Cite this

The accuracy of molecular bond lengths computed by multireference electronic structure methods. / Shepard, Ron; Kedziora, Gary S.; Lischka, Hans; Shavitt, Isaiah; Müller, Thomas; Szalay, P.; Kállay, M.; Seth, Michael.

In: Chemical Physics, Vol. 349, No. 1-3, 16.06.2008, p. 37-57.

Research output: Contribution to journalArticle

Shepard, Ron ; Kedziora, Gary S. ; Lischka, Hans ; Shavitt, Isaiah ; Müller, Thomas ; Szalay, P. ; Kállay, M. ; Seth, Michael. / The accuracy of molecular bond lengths computed by multireference electronic structure methods. In: Chemical Physics. 2008 ; Vol. 349, No. 1-3. pp. 37-57.
@article{77719a8ccf0b4c64b5a75c17f1dd4198,
title = "The accuracy of molecular bond lengths computed by multireference electronic structure methods",
abstract = "We compare experimental Re values with computed Re values for 20 molecules using three multireference electronic structure methods, MCSCF, MR-SDCI, and MR-AQCC. Three correlation-consistent orbital basis sets are used, along with complete basis set extrapolations, for all of the molecules. These data complement those computed previously with single-reference methods. Several trends are observed. The SCF Re values tend to be shorter than the experimental values, and the MCSCF values tend to be longer than the experimental values. We attribute these trends to the ionic contamination of the SCF wave function and to the corresponding systematic distortion of the potential energy curve. For the individual bonds, the MR-SDCI Re values tend to be shorter than the MR-AQCC values, which in turn tend to be shorter than the MCSCF values. Compared to the previous single-reference results, the MCSCF values are roughly comparable to the MP4 and CCSD methods, which are more accurate than might be expected due to the fact that these MCSCF wave functions include no extra-valence electron correlation effects. This suggests that static valence correlation effects, such as near-degeneracies and the ability to dissociate correctly to neutral fragments, play an important role in determining the shape of the potential energy surface, even near equilibrium structures. The MR-SDCI and MR-AQCC methods predict Re values with an accuracy comparable to, or better than, the best single-reference methods (MP4, CCSD, and CCSD(T)), despite the fact that triple and higher excitations into the extra-valence orbital space are included in the single-reference methods but are absent in the multireference wave functions. The computed Re values using the multireference methods tend to be smooth and monotonic with basis set improvement. The molecular structures are optimized using analytic energy gradients, and the timings for these calculations show the practical advantage of using variational wave functions for which the Hellmann-Feynman theorem can be exploited.",
keywords = "Bond length, CBS, Complete basis set limit, Configuration interaction, GVB, MCSCF, MR-AQCC, MR-SDCI, Multiconfiguration self-consistent field, Multiference, R",
author = "Ron Shepard and Kedziora, {Gary S.} and Hans Lischka and Isaiah Shavitt and Thomas M{\"u}ller and P. Szalay and M. K{\'a}llay and Michael Seth",
year = "2008",
month = "6",
day = "16",
doi = "10.1016/j.chemphys.2008.03.009",
language = "English",
volume = "349",
pages = "37--57",
journal = "Chemical Physics",
issn = "0301-0104",
publisher = "Elsevier",
number = "1-3",

}

TY - JOUR

T1 - The accuracy of molecular bond lengths computed by multireference electronic structure methods

AU - Shepard, Ron

AU - Kedziora, Gary S.

AU - Lischka, Hans

AU - Shavitt, Isaiah

AU - Müller, Thomas

AU - Szalay, P.

AU - Kállay, M.

AU - Seth, Michael

PY - 2008/6/16

Y1 - 2008/6/16

N2 - We compare experimental Re values with computed Re values for 20 molecules using three multireference electronic structure methods, MCSCF, MR-SDCI, and MR-AQCC. Three correlation-consistent orbital basis sets are used, along with complete basis set extrapolations, for all of the molecules. These data complement those computed previously with single-reference methods. Several trends are observed. The SCF Re values tend to be shorter than the experimental values, and the MCSCF values tend to be longer than the experimental values. We attribute these trends to the ionic contamination of the SCF wave function and to the corresponding systematic distortion of the potential energy curve. For the individual bonds, the MR-SDCI Re values tend to be shorter than the MR-AQCC values, which in turn tend to be shorter than the MCSCF values. Compared to the previous single-reference results, the MCSCF values are roughly comparable to the MP4 and CCSD methods, which are more accurate than might be expected due to the fact that these MCSCF wave functions include no extra-valence electron correlation effects. This suggests that static valence correlation effects, such as near-degeneracies and the ability to dissociate correctly to neutral fragments, play an important role in determining the shape of the potential energy surface, even near equilibrium structures. The MR-SDCI and MR-AQCC methods predict Re values with an accuracy comparable to, or better than, the best single-reference methods (MP4, CCSD, and CCSD(T)), despite the fact that triple and higher excitations into the extra-valence orbital space are included in the single-reference methods but are absent in the multireference wave functions. The computed Re values using the multireference methods tend to be smooth and monotonic with basis set improvement. The molecular structures are optimized using analytic energy gradients, and the timings for these calculations show the practical advantage of using variational wave functions for which the Hellmann-Feynman theorem can be exploited.

AB - We compare experimental Re values with computed Re values for 20 molecules using three multireference electronic structure methods, MCSCF, MR-SDCI, and MR-AQCC. Three correlation-consistent orbital basis sets are used, along with complete basis set extrapolations, for all of the molecules. These data complement those computed previously with single-reference methods. Several trends are observed. The SCF Re values tend to be shorter than the experimental values, and the MCSCF values tend to be longer than the experimental values. We attribute these trends to the ionic contamination of the SCF wave function and to the corresponding systematic distortion of the potential energy curve. For the individual bonds, the MR-SDCI Re values tend to be shorter than the MR-AQCC values, which in turn tend to be shorter than the MCSCF values. Compared to the previous single-reference results, the MCSCF values are roughly comparable to the MP4 and CCSD methods, which are more accurate than might be expected due to the fact that these MCSCF wave functions include no extra-valence electron correlation effects. This suggests that static valence correlation effects, such as near-degeneracies and the ability to dissociate correctly to neutral fragments, play an important role in determining the shape of the potential energy surface, even near equilibrium structures. The MR-SDCI and MR-AQCC methods predict Re values with an accuracy comparable to, or better than, the best single-reference methods (MP4, CCSD, and CCSD(T)), despite the fact that triple and higher excitations into the extra-valence orbital space are included in the single-reference methods but are absent in the multireference wave functions. The computed Re values using the multireference methods tend to be smooth and monotonic with basis set improvement. The molecular structures are optimized using analytic energy gradients, and the timings for these calculations show the practical advantage of using variational wave functions for which the Hellmann-Feynman theorem can be exploited.

KW - Bond length

KW - CBS

KW - Complete basis set limit

KW - Configuration interaction

KW - GVB

KW - MCSCF

KW - MR-AQCC

KW - MR-SDCI

KW - Multiconfiguration self-consistent field

KW - Multiference

KW - R

UR - http://www.scopus.com/inward/record.url?scp=44749090595&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=44749090595&partnerID=8YFLogxK

U2 - 10.1016/j.chemphys.2008.03.009

DO - 10.1016/j.chemphys.2008.03.009

M3 - Article

VL - 349

SP - 37

EP - 57

JO - Chemical Physics

JF - Chemical Physics

SN - 0301-0104

IS - 1-3

ER -