Journal
COMPUTER PHYSICS COMMUNICATIONS
Volume 183, Issue 8, Pages 1664-1673Publisher
ELSEVIER
DOI: 10.1016/j.cpc.2012.03.004
Keywords
Linear scaling algorithm; Real-space implementation; Density functional theory; Divide-and-conquer; Electronic structure; Embedding effect
Funding
- JST Core Research for Evolutional Science and Technology (CREST)
- JSPS [23310074]
- MEXT of Japan
- MEXT Computational Materials Science Initiative (CMSI) of Japan
- MEXT Open Advanced Research Facilities Initiative
- Fujitsu FX1 at Information Technology Center of Nagoya University
- Grants-in-Aid for Scientific Research [23310074] Funding Source: KAKEN
Ask authors/readers for more resources
The real-space grid based implementation of the Kohn-Sham density functional theory of electrons using the finite difference method for derivatives of variables, has attractive features of parallelizability and applicability to various boundary conditions in addition to universality in target materials. Following the divide-and-conquer strategy, we propose a linear scaling algorithm of it by advancing the algorithm in [F. Shimojo et al., Comput. Phys. Comm. 167 (2005) 151]. In the Kohn-Sham-type equation for a domain, we introduce (i) the density-template potential for density continuity with simple stepwise weight functions and (ii) the embedding potential to take into account all the quantum correlation effects with other overlapping domains in addition to the classical effects of ionic and electronic Coulomb potentials. We thereby realize reasonably high accuracies in atomic forces with relatively small numbers of buffer ions irrespective of the electronic characters of materials. The timing tests on parallel machines demonstrate the linear scaling of the code with little communication time between the domains. (C) 2012 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available