4.7 Article

Can orbital-free density functional theory simulate molecules?

Journal

JOURNAL OF CHEMICAL PHYSICS
Volume 136, Issue 8, Pages -

Publisher

AIP Publishing
DOI: 10.1063/1.3685604

Keywords

ab initio calculations; bond lengths; density functional theory; dissociation energies; molecular electronic states; vibrational states

Funding

  1. Office of Naval Research
  2. National Science Foundation (NSF)
  3. Direct For Mathematical & Physical Scien
  4. Division Of Chemistry [0910563] Funding Source: National Science Foundation

Ask authors/readers for more resources

Orbital-free density functional theory (OFDFT), with its attractive linearly scaling computation cost and low prefactor, is one of the most powerful first principles methods for simulating large systems (similar to 10(4)-10(6) atoms). However, approximating the electron kinetic energy with density functionals limits the accuracy and generality of OFDFT compared to Kohn-Sham density functional theory (KSDFT). In this work, we test whether the Huang-Carter (HC) kinetic energy density functional (KEDF), which contains the physics to properly describe covalently bonded semiconductor materials, can also be used to describe covalent bonds in molecules. In particular, we calculate a variety of homonuclear diatomic molecules with the HC functional within OFDFT. The OFDFT bond dissociation energy, equilibrium bond length, and vibrational frequency of these dimers are in remarkably good agreement with benchmark KSDFT results, given the lack of orbitals in the calculation. We vary the two parameters lambda (controlling the reduced density gradient contribution to the nonlocal kernel) and beta (the exponent of the density in the nonlocal term) present in the HC KEDF and find that the optimal lambda correlates with the magnitude of the highest occupied molecular orbital - lowest unoccupied molecular orbital energy gap. Although the HC KEDF represents a significant improvement over previous KEDFs in describing covalent systems, deficiencies still exist. Despite the similar overall shape of the KSDFT and OFDFT ground state electron densities, the electron density within the bonding region is still quite different. Furthermore, OFDFT is not yet able to give reasonable description of magnetic states. The energy orderings of the triplet and singlet states of Si-2 and Al family dimers are not consistent with KSDFT or experimental results and the spin polarization distributions also differ widely between the two theories. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3685604]

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available