Partitioning scheme for density functional calculations of extended systems.
We show that, at least for the ground electronic state of systems treated using semilocal density functionals (like in local density or generalized gradient approximations), a calculation of the entire extended nonperiodic system consisting of several well distinguished parts (e.g., a collection of molecules) can be replaced with a finite set of calculations on specifically chosen smaller subsystems that overlap with each other. Every subsystem is terminated with link (or pseudo) atoms (or groups of atoms) to reduce the effect of the termination. However, because of the particular choice of the subsystems, the effect of the link atoms is largely compensated in the final total energy if the subsystems are chosen sufficiently large. In fact, we prove that the proposed method should result in nearly the same total energy, electronic density and atomic forces as a single (considered as a reference) density functional calculation on the entire system. Our method, however, should be much more efficient due to unfavorable scaling of the modern electronic structure methods with the system size. The method is illustrated on examples of serine water, lysine-water and lysine dimer systems. We also discuss possible approximate applications of our method for quantum-classical calculations of extended systems, when, as compared to widely used quantum-mechanical/molecular-mechanical methods, the problem of the quantum cluster boundary can be eliminated to a large degree.
Department of Physics, King's College London, London WC2R 2LS, United Kingdom.
Jiangping
J
He
C
C
Di Paola
L
L
Kantorovich