Ab initio correlation functionals from second-order perturbation theory
J. Chem. Phys. 125, 104108 (2006); doi:10.1063/1.2212936
Published 12 September 2006
You are not logged in to this journal. Log in
Orbital-dependent exchange-correlation functionals are not limited by the explicit dependence on the density and present an attractive alternative to conventional functionals. With the successful implementation of the exact orbital-dependent exchange functional, the challenge lies in developing orbital-dependent approximations for the correlation functional. Ab initio many-body methods can provide such approximations. In particular, perturbation theory with the Kohn-Sham model as the reference [Görling and Levy, Phys. Rev. A 50, 196 (1994)] defines the exact correlation functional via an infinite perturbation series. The second-order term of these series gives the lowest-order approximation to the correlation functional. However, it has been suggested [Bartlett et al., J. Chem. Phys. 122, 034104 (2005)] that the Kohn-Sham Hamiltonian is not the optimal choice for the perturbation expansion and a different reference Hamiltonian may lead to an improved perturbation series and more accurate second-order approximation. Here, we demonstrate explicitly that the modified series can be used to define superior functional and potential. We present results of atomic and molecular calculations with both second-order functionals. Our results demonstrate that the modified functional offers a significantly improved description of the correlation effects as it does not suffer from convergence problems and results in energies and densities that are more accurate than those obtained with second-order Møller-Plesset perturbation theory or generalized-gradient approximation functionals.
©2006 American Institute of Physics
| History: | Received 8 March 2006; accepted 17 May 2006; published 12 September 2006 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/125/104108/1 |
Supplemental Material
- 305625jcp.pdf (213 kB) 23-Aug-2006 9:9
- README.TXT (1 kB) 15-Sep-2006 16:12
KEYWORDS and PACS
RELATED DATABASES
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (56)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- P. Hohenberg and W. Kohn,
Phys. Rev. 136, B864 (1964) . - W. Kohn and L. J. Sham,
Phys. Rev. 140, 1133 (1965) . - S. H. Vosko, L. Wilk, and M. Nusair,
Can. J. Phys. 58, 1200 (1980) . - J. P. Perdew, Phys. Rev. B 33, 8822 (1986).
- A. D. Becke, Phys. Rev. A 38, 3098 (1988).
- C. T. Lee, W. T. Yang, and R. G. Parr, Phys. Rev. B 37, 785 (1988).
- J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
- A. Becke, J. Chem. Phys. 98, 5648 (1993).
- J. M. Tao, J. P. Perdew, V. N. Staroverov, and G. E. Scuseria, Phys. Rev. Lett. 91, 146401 (2003).
- J. P. Perdew and K. Schmidt, in Density Functional Theory and Its Applications to Materials, edited by V. V. Doren, C. V. Alsenoy, and P. Geerlings (American Institute of Physics, New York, 2001).
- E. Engel and R. M. Dreizler,
J. Comput. Chem. 20, 31 (1999) . - A. Görling and M. Levy, Phys. Rev. A 50, 196 (1994).
- R. T. Sharp and G. K. Horton,
Phys. Rev. 90, 317 (1953) . - J. D. Talman and W. F. Shadwick, Phys. Rev. A 14, 36 (1976).
- V. Sahni, J. Gruenebaum, and J. P. Perdew, Phys. Rev. B 26, 4371 (1982).
- M. R. Norman and D. D. Koelling, Phys. Rev. B 30, 5530 (1984).
- J. B. Krieger, Y. Li, and G. J. Iafrate,
Phys. Lett. A 148, 470 (1990) . - M. Gruning, O. V. Gritsenko, and E. J. Baerends, J. Chem. Phys. 116, 6435 (2002).
- F. Della Sala and A. Görling, J. Chem. Phys. 115, 5718 (2001).
- S. Ivanov, S. Hirata, and R. J. Bartlett, Phys. Rev. Lett. 83, 5455 (1999).
- A. Görling, Phys. Rev. Lett. 83, 5459 (1999).
- S. Kummel and J. P. Perdew, Phys. Rev. Lett. 90, 043004 (2003).
- It is a common viewpoint that the GGA exchange and correlation components compensate the deficiencies of each other, providing more accurate approximation for the total exchange-correlation energy than for the exchange and correlation separately. An illustration of this is the PBE correlation potential of He atom which is almost exactly the negative of the exact correlation potential.
- A. Gorling and M. Levy, Phys. Rev. B 47, 13105 (1993).
- E. Engel, A. F. Bonetti, S. Keller, I. Andrejkovics, and R. M. Dreizler, Phys. Rev. A 58, 964 (1998).
- I. Grabowski, S. Hirata, S. Ivanov, and R. J. Bartlett, J. Chem. Phys. 116, 4415 (2002).
- A. F. Bonetti, E. Engel, R. N. Schmid, and R. M. Dreizler, Phys. Rev. Lett. 86, 2241 (2001).
- M. Warken,
Chem. Phys. Lett. 237, 256 (1995) . - M. Ernzerhof,
Chem. Phys. Lett. 263, 499 (1996) . - M. Seidl, J. P. Perdew, and S. Kurth, Phys. Rev. Lett. 84, 5070 (2000).
- P. Bour,
J. Comput. Chem. 21, 8 (2000) . - R. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005).
- R. J. Bartlett, V. F. Lotrich, and I. V. Schweigert, J. Chem. Phys. 123, 074106 (2005).
- P. Mori-Sanchez, Q. Wu, and W. T. Yang, J. Chem. Phys. 123, 062204 (2005).
- Note that when both q and p are occupied, the terms under the summation change sign with respect to the interchange of q and p, and therefore vanish.
- A. Görling and M. Levy, Int. J. Quantum Chem. 29, 93 (1995).
- Note that the definition of the KS potential via the functional derivative is just a consequence of the genuine condition that the model density must reproduce the exact one. Combining the density condition with the fact that the ground-state density minimizes the functional, one obtains the functional derivative definition.
- See EPAPS Document No. E-JCPSA6-125-305625 for derivation of the final expressions for the PT2SC functional and potential. This document can be reached via a direct link in the online article's HTML reference section or via the EPAPS homepage (http://www.aip.org/pubservs/epaps.html). [EPAPS]
- If the products of occupied and unoccupied orbitals form a linearly independent set, then satisfying Eq. (23) at Nocc×Nunocc different points in space requires that

i|nlx
x|
a
are exactly zero. However, in virtually any finite basis, these products are linearly dependent and Eq. (23) can be viewed as a weighted fit such that 
i|nlx
x|
a
are minimized. - J. F. Stanton, J. Gauss, J. D. Watts et al., ACES II, a program product of the Quantum Theory Project, University of Florida.
- S. Ivanov, S. Hirata, and R. J. Barlett, J. Chem. Phys. 116, 1269 (2002).
- P. O. Widmark, P. A. Malmqvist, and B. O. Roos,
Theor. Chim. Acta 77, 291 (1990) . - P. O. Widmark, B. J. Persson, and B. O. Roos,
Theor. Chim. Acta 79, 419 (1991) . - T. H. Dunning, J. Chem. Phys. 90, 1007 (1989).
- S. Hirata, S. Ivanov, I. Grabowski, R. J. Bartlett, K. Burke, and J. D. Talman, J. Chem. Phys. 115, 1635 (2001).
- S. Hirata, S. Ivanov, I. Grabowski, and R. J. Bartlett, J. Chem. Phys. 116, 6468 (2002).
- J. C. Slater,
Phys. Rev. 81, 385 (1951) . - E. Engel, H. Jiang, and A. F. Bonetti, Phys. Rev. A 72, 052503 (2005).
- Y. M. Niquet, M. Fuchs, and X. Gonze, J. Chem. Phys. 118, 9504 (2003).
- I. Grabowski and V. Lotrich,
Mol. Phys. 103, 2085 (2005) . - Note that we did not include the PBE correlation energies in Table II because the absolute values of GGA energies can be quite different from wave function values. However, it is the ab initio nature of the PT2 and PT2SC functionals that allows us to compare the absolute values of the correlation energy. See also Ref. 56.
- L. Curtiss, K. Raghavachari, G. W. Trucks, and J. A. Pople, J. Chem. Phys. 94, 7221 (1991).
- NIST Computational Chemistry Comparison and Benchmark Database, NIST Standard Reference Database No. 101, Release 12, edited by R. D. Johnson III, 2005 (http://srdata.nist.gov/cccbdb).
- R. J. Bartlett, in Chemistry for the 21st Century, edited by E. Keinan and I. Schechter (Wiley-VCH, Weinheim, 2001), pp. 271286.
- H. Jiang and E. Engel, J. Chem. Phys. 123, 224101 (2005).
- E. J. Baerends and O. V. Gritsenko, J. Chem. Phys. 123, 062202 (2005).
P. Taylor, VPROPS;
and T. Helgaker, H. J. A. Jensen, P. Jörgensen, J. Olsen, and P. R. Taylor, ABACUS.








