Symmetry and broken symmetry in molecular orbital description of unstable molecules IV

Comparison between single- and multi-reference computational results for antiaromtic molecules

Toru Saito, Satomichi Nishihara, Shusuke Yamanaka, Yasutaka Kitagawa, Takashi Kawakami, Satoru Yamada, Hiroshi Isobe, Mitsutaka Okumura, Kizashi Yamaguchi

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

First principle calculations of the effective exchange integrals (J) in the Heisenberg model for diradical species are presented for both symmetry-adapted multi-reference (MR) and single-reference broken-symmetry (BS) methods. The Mukherjee-type state-specific MR coupled cluster singles and doubles (MkCCSD) method with several different reference orbitals including BS natural orbitals is used to calculate the singlet-triplet energy gaps (S-T energy gap or 2J) and diradical characters for the antiaromatic molecules [a cyclopropenyl anion (CPA), b cyclobutadiene (CBD), and c cyclopentadienyl cation (CPC)], the cyclobutadiene derivatives with polar substitutents [d aminocyclobutadiene (ACBD), e formylcyclobutadiene (FCBD), and f 1-amino-2-formyl-cyclobutadiene (AFCBD)] and finally the cyclobutadine derivatives with radical substitutents [g 1,2-bis(methylene)cyclobutadiene (1,2-BMCBD) and h 1,3-bis(methylene)cyclobutadiene (1,3-BMCBD)]. For the BS methods, the spin-unrestricted Hartree-Fock based CCSD (UHF-CCSD), the CCD with the spin-unrestricted Brueckner determinant (UBD), and BS density functional theory (UDFT) computations are performed. Comparison between MkCCSD and the UHF-CCSD results indicates that spin-contamination of UHF-CCSD solutions still remains. In comparison with UHF-CCSD, the UBD results show that spin-contamination involved in BS solutions is greatly suppressed. To eliminate the spin contamination, an approximate spin-projection (AP) scheme is applied to the BS solutions. The AP procedure with the use of the expectation value of the total-spin operator corresponding to UHF-CCSD and UBD results yields good agreement with the MkCCSD results. As for the AP correction of the UDFT methods, three different computational schemes for predicting the expectation value of the total-spin operator are examined. Systematic comparisons between these methods are presented for the S-T energy gaps (2J). Implications of the present computational results have been discussed in relation to the design of magnetic oligomers and polymers.

Original languageEnglish
Pages (from-to)749-763
Number of pages15
JournalTheoretical Chemistry Accounts
Volume130
Issue number4-6
DOIs
Publication statusPublished - Dec 2011
Externally publishedYes

Fingerprint

Molecular orbitals
broken symmetry
molecular orbitals
Energy gap
Contamination
Molecules
Density functional theory
symmetry
Derivatives
molecules
Oligomers
Charge coupled devices
Anions
Cations
determinants
Polymers
contamination
projection
methylene
density functional theory

Keywords

  • Approximate spin projection procedure
  • Diradical character
  • Multireference coupled cluster
  • Spin contamination

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Symmetry and broken symmetry in molecular orbital description of unstable molecules IV : Comparison between single- and multi-reference computational results for antiaromtic molecules. / Saito, Toru; Nishihara, Satomichi; Yamanaka, Shusuke; Kitagawa, Yasutaka; Kawakami, Takashi; Yamada, Satoru; Isobe, Hiroshi; Okumura, Mitsutaka; Yamaguchi, Kizashi.

In: Theoretical Chemistry Accounts, Vol. 130, No. 4-6, 12.2011, p. 749-763.

Research output: Contribution to journalArticle

Saito, Toru ; Nishihara, Satomichi ; Yamanaka, Shusuke ; Kitagawa, Yasutaka ; Kawakami, Takashi ; Yamada, Satoru ; Isobe, Hiroshi ; Okumura, Mitsutaka ; Yamaguchi, Kizashi. / Symmetry and broken symmetry in molecular orbital description of unstable molecules IV : Comparison between single- and multi-reference computational results for antiaromtic molecules. In: Theoretical Chemistry Accounts. 2011 ; Vol. 130, No. 4-6. pp. 749-763.
@article{145c1823bbcc478eb1c04232f39512f7,
title = "Symmetry and broken symmetry in molecular orbital description of unstable molecules IV: Comparison between single- and multi-reference computational results for antiaromtic molecules",
abstract = "First principle calculations of the effective exchange integrals (J) in the Heisenberg model for diradical species are presented for both symmetry-adapted multi-reference (MR) and single-reference broken-symmetry (BS) methods. The Mukherjee-type state-specific MR coupled cluster singles and doubles (MkCCSD) method with several different reference orbitals including BS natural orbitals is used to calculate the singlet-triplet energy gaps (S-T energy gap or 2J) and diradical characters for the antiaromatic molecules [a cyclopropenyl anion (CPA), b cyclobutadiene (CBD), and c cyclopentadienyl cation (CPC)], the cyclobutadiene derivatives with polar substitutents [d aminocyclobutadiene (ACBD), e formylcyclobutadiene (FCBD), and f 1-amino-2-formyl-cyclobutadiene (AFCBD)] and finally the cyclobutadine derivatives with radical substitutents [g 1,2-bis(methylene)cyclobutadiene (1,2-BMCBD) and h 1,3-bis(methylene)cyclobutadiene (1,3-BMCBD)]. For the BS methods, the spin-unrestricted Hartree-Fock based CCSD (UHF-CCSD), the CCD with the spin-unrestricted Brueckner determinant (UBD), and BS density functional theory (UDFT) computations are performed. Comparison between MkCCSD and the UHF-CCSD results indicates that spin-contamination of UHF-CCSD solutions still remains. In comparison with UHF-CCSD, the UBD results show that spin-contamination involved in BS solutions is greatly suppressed. To eliminate the spin contamination, an approximate spin-projection (AP) scheme is applied to the BS solutions. The AP procedure with the use of the expectation value of the total-spin operator corresponding to UHF-CCSD and UBD results yields good agreement with the MkCCSD results. As for the AP correction of the UDFT methods, three different computational schemes for predicting the expectation value of the total-spin operator are examined. Systematic comparisons between these methods are presented for the S-T energy gaps (2J). Implications of the present computational results have been discussed in relation to the design of magnetic oligomers and polymers.",
keywords = "Approximate spin projection procedure, Diradical character, Multireference coupled cluster, Spin contamination",
author = "Toru Saito and Satomichi Nishihara and Shusuke Yamanaka and Yasutaka Kitagawa and Takashi Kawakami and Satoru Yamada and Hiroshi Isobe and Mitsutaka Okumura and Kizashi Yamaguchi",
year = "2011",
month = "12",
doi = "10.1007/s00214-011-0941-9",
language = "English",
volume = "130",
pages = "749--763",
journal = "Theoretical Chemistry Accounts",
issn = "1432-881X",
publisher = "Springer New York",
number = "4-6",

}

TY - JOUR

T1 - Symmetry and broken symmetry in molecular orbital description of unstable molecules IV

T2 - Comparison between single- and multi-reference computational results for antiaromtic molecules

AU - Saito, Toru

AU - Nishihara, Satomichi

AU - Yamanaka, Shusuke

AU - Kitagawa, Yasutaka

AU - Kawakami, Takashi

AU - Yamada, Satoru

AU - Isobe, Hiroshi

AU - Okumura, Mitsutaka

AU - Yamaguchi, Kizashi

PY - 2011/12

Y1 - 2011/12

N2 - First principle calculations of the effective exchange integrals (J) in the Heisenberg model for diradical species are presented for both symmetry-adapted multi-reference (MR) and single-reference broken-symmetry (BS) methods. The Mukherjee-type state-specific MR coupled cluster singles and doubles (MkCCSD) method with several different reference orbitals including BS natural orbitals is used to calculate the singlet-triplet energy gaps (S-T energy gap or 2J) and diradical characters for the antiaromatic molecules [a cyclopropenyl anion (CPA), b cyclobutadiene (CBD), and c cyclopentadienyl cation (CPC)], the cyclobutadiene derivatives with polar substitutents [d aminocyclobutadiene (ACBD), e formylcyclobutadiene (FCBD), and f 1-amino-2-formyl-cyclobutadiene (AFCBD)] and finally the cyclobutadine derivatives with radical substitutents [g 1,2-bis(methylene)cyclobutadiene (1,2-BMCBD) and h 1,3-bis(methylene)cyclobutadiene (1,3-BMCBD)]. For the BS methods, the spin-unrestricted Hartree-Fock based CCSD (UHF-CCSD), the CCD with the spin-unrestricted Brueckner determinant (UBD), and BS density functional theory (UDFT) computations are performed. Comparison between MkCCSD and the UHF-CCSD results indicates that spin-contamination of UHF-CCSD solutions still remains. In comparison with UHF-CCSD, the UBD results show that spin-contamination involved in BS solutions is greatly suppressed. To eliminate the spin contamination, an approximate spin-projection (AP) scheme is applied to the BS solutions. The AP procedure with the use of the expectation value of the total-spin operator corresponding to UHF-CCSD and UBD results yields good agreement with the MkCCSD results. As for the AP correction of the UDFT methods, three different computational schemes for predicting the expectation value of the total-spin operator are examined. Systematic comparisons between these methods are presented for the S-T energy gaps (2J). Implications of the present computational results have been discussed in relation to the design of magnetic oligomers and polymers.

AB - First principle calculations of the effective exchange integrals (J) in the Heisenberg model for diradical species are presented for both symmetry-adapted multi-reference (MR) and single-reference broken-symmetry (BS) methods. The Mukherjee-type state-specific MR coupled cluster singles and doubles (MkCCSD) method with several different reference orbitals including BS natural orbitals is used to calculate the singlet-triplet energy gaps (S-T energy gap or 2J) and diradical characters for the antiaromatic molecules [a cyclopropenyl anion (CPA), b cyclobutadiene (CBD), and c cyclopentadienyl cation (CPC)], the cyclobutadiene derivatives with polar substitutents [d aminocyclobutadiene (ACBD), e formylcyclobutadiene (FCBD), and f 1-amino-2-formyl-cyclobutadiene (AFCBD)] and finally the cyclobutadine derivatives with radical substitutents [g 1,2-bis(methylene)cyclobutadiene (1,2-BMCBD) and h 1,3-bis(methylene)cyclobutadiene (1,3-BMCBD)]. For the BS methods, the spin-unrestricted Hartree-Fock based CCSD (UHF-CCSD), the CCD with the spin-unrestricted Brueckner determinant (UBD), and BS density functional theory (UDFT) computations are performed. Comparison between MkCCSD and the UHF-CCSD results indicates that spin-contamination of UHF-CCSD solutions still remains. In comparison with UHF-CCSD, the UBD results show that spin-contamination involved in BS solutions is greatly suppressed. To eliminate the spin contamination, an approximate spin-projection (AP) scheme is applied to the BS solutions. The AP procedure with the use of the expectation value of the total-spin operator corresponding to UHF-CCSD and UBD results yields good agreement with the MkCCSD results. As for the AP correction of the UDFT methods, three different computational schemes for predicting the expectation value of the total-spin operator are examined. Systematic comparisons between these methods are presented for the S-T energy gaps (2J). Implications of the present computational results have been discussed in relation to the design of magnetic oligomers and polymers.

KW - Approximate spin projection procedure

KW - Diradical character

KW - Multireference coupled cluster

KW - Spin contamination

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

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

U2 - 10.1007/s00214-011-0941-9

DO - 10.1007/s00214-011-0941-9

M3 - Article

VL - 130

SP - 749

EP - 763

JO - Theoretical Chemistry Accounts

JF - Theoretical Chemistry Accounts

SN - 1432-881X

IS - 4-6

ER -