Exact nonadditive kinetic potentials for embedded density functional theory
Abstract
We describe an embedded density functional theory (DFT) protocol in which the nonadditive kinetic energy component of the embedding potential is treated exactly. At each iteration of the Kohn–Sham equations for constrained electron density, the Zhao–Morrison–Parr constrained search method for constructing Kohn–Sham orbitals is combined with the King-Handy expression for the exact kinetic potential. We use this formally exact embedding protocol to calculate ionization energies for a series of three- and four-electron atomic systems, and the results are compared to embedded DFT calculations that utilize the Thomas–Fermi (TF) and the Thomas–Fermi–von Weisacker approximations to the kinetic energy functional. These calculations illustrate the expected breakdown due to the TF approximation for the nonadditive kinetic potential, with errors of 30%–80% in the calculated ionization energies; by contrast, the exact protocol is found to be accurate and stable. To significantly improve the convergence of the new protocol, we introduce a density-based switching function to map between the exact nonadditive kinetic potential and the TF approximation in the region of the nuclear cusp, and we demonstrate that this approximation has little effect on the accuracy of the calculated ionization energies. Finally, we describe possible extensions of the exact protocol to perform accurate embedded DFT calculations in large systems with strongly overlapping subsystem densities.
Additional Information
© 2010 American Institute of Physics. Received 26 April 2010; accepted 12 July 2010; published online 23 August 2010. This work was partially supported by a Camille and Henry Dreyfus Foundation New Faculty Award and an Alfred P. Sloan Foundation Research Fellowship (T.F.M.).Attached Files
Published - Goodpaster2010p11497J_Chem_Phys.pdf
Files
Name | Size | Download all |
---|---|---|
md5:eb3ca6e0f323662f491be6c71f8c3cb9
|
428.4 kB | Preview Download |
Additional details
- Eprint ID
- 20382
- Resolver ID
- CaltechAUTHORS:20101011-131128803
- Camille and Henry Dreyfus Foundation
- Alfred P. Sloan Foundation
- Created
-
2010-11-29Created from EPrint's datestamp field
- Updated
-
2021-11-08Created from EPrint's last_modified field