A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigenfunctions?
J. Chem. Phys. 118, 9084 (2003); doi:10.1063/1.1568735
Issue Date: 22 May 2003
You are not logged in to this journal. Log in
Spin-complete versions of the spin-flip configuration-interaction-singles (SF-CIS) approach have been investigated to determine the impact of making the wave function an eigenfunction of
2. The method has been implemented within an extended restricted active space configuration interaction formalism. Spin-complete results are presented for excitation energies, equilibrium geometries, and potential energy curves for dissociation of a single bond in several small molecules. The effect of different orbital choices has also been investigated. The spin-complete results are compared both to results using the original spin-flip method and to more computationally expensive benchmarks. Using spin eigenfunctions dramatically improves upon the accuracy of the SF-CIS approach. ©2003 American Institute of Physics.
| History: | Received 2 January 2003; accepted 27 February 2003 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/118/9084/1 |
KEYWORDS and PACS
RELATED DATABASES
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (26)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- J. R. Thomas, B. J. DeLeeuw, G. Vacek, T. D. Crawford, Y. Yamaguchi, and H. F. Schaefer, J. Chem. Phys. 99, 403 (1993).
- T. Helgaker, J. Gauss, P. Jørgensen, and J. Olsen, J. Chem. Phys. 106, 6430 (1997).
- D. Feller and K. A. Peterson, J. Chem. Phys. 108, 154 (1998).
- R. J. Bartlett and J. F. Stanton, in Computational Chemistry, edited by K. B. Lipkowitz and D. B. Boyd (VCH, New York, 1994), Vol. 5, pp. 65169.
- B. O. Roos, P. R. Taylor, and P. E. M. Siegbahn,
Chem. Phys. 48, 157 (1980) . - K. Ruedenberg, M. Schmidt, M. Gilbert, and S. Elbert,
Chem. Phys. 71, 41 (1982) . - M. Schmidt and M. Gordon,
Annu. Rev. Phys. Chem. 49, 233 (1998) . - A. I. Krylov, C. D. Sherrill, E. F. C. Byrd, and M. Head-Gordon, J. Chem. Phys. 109, 10669 (1998).
- A. I. Krylov,
Chem. Phys. Lett. 338, 375 (2001) . - A. I. Krylov,
Chem. Phys. Lett. 350, 522 (2001) . - A. I. Krylov and C. D. Sherrill, J. Chem. Phys. 116, 3194 (2002).
- Y. Shao, M. Head-Gordon, and A. Krylov, J. Chem. Phys. 118, 4807 (2003).
- A. I. Krylov, J. Chem. Phys. 113, 6052 (2000).
- C. D. Sherrill and H. F. Schaefer, The Configuration Interaction Method: Advances in Highly Correlated Approaches, in Advances in Quantum Chemistry, Vol. 34, edited by P.-O. Löwdin (Academic, New York, 1999), pp. 143269.
- J. Olsen, B. O. Roos, P. Jørgensen, and H. J. Aa. Jensen, J. Chem. Phys. 89, 2185 (1988).
- C. D. Sherrill and H. F. Schaefer,
J. Phys. Chem. 100, 6069 (1996) . - J. Stanton, J. Chem. Phys. 101, 8928 (1994).
- In case of two nonequivalent fragments, such as Be and Ne, or two H2 molecules with one of the bonds being stretched, the lowest energy HartreeFock triplet state is naturally localized on Be or on the stretched H2 molecule, respectively. However, in case of two identical fragments in a symmetric configuration, the lowest triplet state cannot be described by a single Slater determinant using symmetry-adapted orbitals. Therefore, it is necessary to use a symmetry-broken HartreeFock wave function to describe localization in the case of symmetrically arranged, noninteracting equivalent fragments.
- H. Koch, H. Jørgen Aa. Jensen, and P. Jørgensen, J. Chem. Phys. 93, 3345 (1990).
- Since the target states are found by diagonalization of the Hamiltonian matrix, the block-diagonal structure of the latter ensures that the eigenvalues of the subsystems A and B are also eigenvalues of the combined system AB. The total SF energies, which are the sum of the high-spin reference energy and the corresponding transition energy, are thus also size-consistent due to size-consistency of the HartreeFock model. Note that the matrix of the Hamiltonian (2) has eigenvalues corresponding to the simultaneous excitations of both fragments, or to the charge transfer states, i.e., those from the HAB,AB block. The interpretation of these roots requires caution as described in footnote 22 of Ref. 17. For example, due to the truncated nature of the excitation operators, the EOM and CI models fail to describe situations when both fragments are excited with the same accuracy as the excitations on individual fragments. Due to similar considerations, Eq. (1) does not mean that for two noninteracting fragments, the total SF energy of the composite system equals the sum of total SF energy of individual fragments; rather, the total SF energy is equal to the sum of SF energy of fragment A and HartreeFock energy for the fragment B.
- J. B. Foresman, M. Head-Gordon, J. A. Pople, and M. J. Frisch,
J. Phys. Chem. 96, 135 (1992) . - S. Huzinaga, J. Chem. Phys. 42, 1293 (1965).
- T. H. Dunning, J. Chem. Phys. 53, 2823 (1970).
- W. D. Laidig, P. Saxe, and R. J. Bartlett, J. Chem. Phys. 86, 887 (1987).
- L. Salem and C. Rowland,
Angew. Chem., Int. Ed. Engl. 11, 92 (1972) . - L. Slipchenko and A. Krylov, J. Chem. Phys. 117, 4694 (2002).








