TY - JOUR
T1 - Grid-based energy density analysis
T2 - Implementation and assessment
AU - Imamura, Yutaka
AU - Takahashi, Asuka
AU - Nakai, Hiromi
PY - 2007/1/29
Y1 - 2007/1/29
N2 - Grid-based energy density analysis (grid-EDA) that decomposes the total energy into atomic energies by a space-partitioning function is proposed. The kinetic energy, nuclear attraction, and exchange-correlation functional are evaluated on grid points and are split into atomic contributions. To reduce numerical errors in the conventional scheme of numerical integration, the electronic Coulomb and HF exchange interactions are evaluated by the pseudospectral method, which was first applied to an ab initio method by Friesner [Chem. Phys. Lett. 116, 39 (1985)], and are decomposed into atomic contributions. Grid-EDA using the pseudospectral method succeeds in ensuring less than 1 kcalmol error in total energies for small molecules and providing reliable atomic energy contributions for the problematic lithium cluster, which exhibits a strong basis-set dependence for Mulliken-type EDA. Also, site-dependent atomization energies are estimated by grid-EDA for cluster models such as Li48, C41 H60, and Mg32 O32. Grid-EDA reveals that these models imitate crystal environments reasonably because atomization energies estimated from the inner atoms of the models are close to the experimental cohesive energies.
AB - Grid-based energy density analysis (grid-EDA) that decomposes the total energy into atomic energies by a space-partitioning function is proposed. The kinetic energy, nuclear attraction, and exchange-correlation functional are evaluated on grid points and are split into atomic contributions. To reduce numerical errors in the conventional scheme of numerical integration, the electronic Coulomb and HF exchange interactions are evaluated by the pseudospectral method, which was first applied to an ab initio method by Friesner [Chem. Phys. Lett. 116, 39 (1985)], and are decomposed into atomic contributions. Grid-EDA using the pseudospectral method succeeds in ensuring less than 1 kcalmol error in total energies for small molecules and providing reliable atomic energy contributions for the problematic lithium cluster, which exhibits a strong basis-set dependence for Mulliken-type EDA. Also, site-dependent atomization energies are estimated by grid-EDA for cluster models such as Li48, C41 H60, and Mg32 O32. Grid-EDA reveals that these models imitate crystal environments reasonably because atomization energies estimated from the inner atoms of the models are close to the experimental cohesive energies.
UR - http://www.scopus.com/inward/record.url?scp=33846415416&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=33846415416&partnerID=8YFLogxK
U2 - 10.1063/1.2428290
DO - 10.1063/1.2428290
M3 - Article
AN - SCOPUS:33846415416
VL - 126
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 3
M1 - 034103
ER -