Molecular symmetry. II. Gradient of electronic energy with respect to nuclear coordinates

Michel Dupuis, Harry F. King

Research output: Contribution to journalArticle

275 Citations (Scopus)

Abstract

Symmetry methods employed in the HONDO ab initio polyatomic SCF program are extended to the analytic computation of energy gradients. Validity of the Hellmann-Feynman theorem is not assumed, i.e., all two-electron contributions to the gradient are included explicitly. The method is geared to the efficient computation of entire blocks of two-electron integrals. Just one of a set of symmetrically related blocks must be computed. The gradient contribution from each unique block is multiplied by q4, the number of equivalent blocks, and added into a "skeleton gradient vector," all other blocks are simply omitted. After processing molecular integrals, the true gradient vector is generated by projecting the symmetric component out of the skeleton vector. The analysis is based on Eqs. (26) and (33) which are valid for many variational wavefunctions including restricted closed shell and unrestricted open shell self-consistent field functions. We also extend the use of translational symmetry proposed previously by Morokuma et al. To illustrate the method, the gradient of the restricted SCF energy is computed for eclipsed ethane using a Pople-type 631G** basis and D3h symmetry. The same calculation is repeated using various subgroups of D3h. Computation times for SCF and for the gradient are each roughly inversely proportional to the order of the group, and for a given symmetry, the gradient computation takes about two and a half times as long as SCF.

Original languageEnglish
Pages (from-to)3998-4004
Number of pages7
JournalJournal of Chemical Physics
Volume68
Issue number9
Publication statusPublished - 1978

Fingerprint

gradients
self consistent fields
symmetry
electronics
energy
musculoskeletal system
Ethane
Electrons
Wave functions
Hellmann-Feynman theorem
subgroups
ethane
Processing
electrons

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Molecular symmetry. II. Gradient of electronic energy with respect to nuclear coordinates. / Dupuis, Michel; King, Harry F.

In: Journal of Chemical Physics, Vol. 68, No. 9, 1978, p. 3998-4004.

Research output: Contribution to journalArticle

@article{3bd7b8ff10824b9a97432860291ca019,
title = "Molecular symmetry. II. Gradient of electronic energy with respect to nuclear coordinates",
abstract = "Symmetry methods employed in the HONDO ab initio polyatomic SCF program are extended to the analytic computation of energy gradients. Validity of the Hellmann-Feynman theorem is not assumed, i.e., all two-electron contributions to the gradient are included explicitly. The method is geared to the efficient computation of entire blocks of two-electron integrals. Just one of a set of symmetrically related blocks must be computed. The gradient contribution from each unique block is multiplied by q4, the number of equivalent blocks, and added into a {"}skeleton gradient vector,{"} all other blocks are simply omitted. After processing molecular integrals, the true gradient vector is generated by projecting the symmetric component out of the skeleton vector. The analysis is based on Eqs. (26) and (33) which are valid for many variational wavefunctions including restricted closed shell and unrestricted open shell self-consistent field functions. We also extend the use of translational symmetry proposed previously by Morokuma et al. To illustrate the method, the gradient of the restricted SCF energy is computed for eclipsed ethane using a Pople-type 631G** basis and D3h symmetry. The same calculation is repeated using various subgroups of D3h. Computation times for SCF and for the gradient are each roughly inversely proportional to the order of the group, and for a given symmetry, the gradient computation takes about two and a half times as long as SCF.",
author = "Michel Dupuis and King, {Harry F.}",
year = "1978",
language = "English",
volume = "68",
pages = "3998--4004",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "9",

}

TY - JOUR

T1 - Molecular symmetry. II. Gradient of electronic energy with respect to nuclear coordinates

AU - Dupuis, Michel

AU - King, Harry F.

PY - 1978

Y1 - 1978

N2 - Symmetry methods employed in the HONDO ab initio polyatomic SCF program are extended to the analytic computation of energy gradients. Validity of the Hellmann-Feynman theorem is not assumed, i.e., all two-electron contributions to the gradient are included explicitly. The method is geared to the efficient computation of entire blocks of two-electron integrals. Just one of a set of symmetrically related blocks must be computed. The gradient contribution from each unique block is multiplied by q4, the number of equivalent blocks, and added into a "skeleton gradient vector," all other blocks are simply omitted. After processing molecular integrals, the true gradient vector is generated by projecting the symmetric component out of the skeleton vector. The analysis is based on Eqs. (26) and (33) which are valid for many variational wavefunctions including restricted closed shell and unrestricted open shell self-consistent field functions. We also extend the use of translational symmetry proposed previously by Morokuma et al. To illustrate the method, the gradient of the restricted SCF energy is computed for eclipsed ethane using a Pople-type 631G** basis and D3h symmetry. The same calculation is repeated using various subgroups of D3h. Computation times for SCF and for the gradient are each roughly inversely proportional to the order of the group, and for a given symmetry, the gradient computation takes about two and a half times as long as SCF.

AB - Symmetry methods employed in the HONDO ab initio polyatomic SCF program are extended to the analytic computation of energy gradients. Validity of the Hellmann-Feynman theorem is not assumed, i.e., all two-electron contributions to the gradient are included explicitly. The method is geared to the efficient computation of entire blocks of two-electron integrals. Just one of a set of symmetrically related blocks must be computed. The gradient contribution from each unique block is multiplied by q4, the number of equivalent blocks, and added into a "skeleton gradient vector," all other blocks are simply omitted. After processing molecular integrals, the true gradient vector is generated by projecting the symmetric component out of the skeleton vector. The analysis is based on Eqs. (26) and (33) which are valid for many variational wavefunctions including restricted closed shell and unrestricted open shell self-consistent field functions. We also extend the use of translational symmetry proposed previously by Morokuma et al. To illustrate the method, the gradient of the restricted SCF energy is computed for eclipsed ethane using a Pople-type 631G** basis and D3h symmetry. The same calculation is repeated using various subgroups of D3h. Computation times for SCF and for the gradient are each roughly inversely proportional to the order of the group, and for a given symmetry, the gradient computation takes about two and a half times as long as SCF.

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

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

M3 - Article

AN - SCOPUS:0011398687

VL - 68

SP - 3998

EP - 4004

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 9

ER -