We propose a simple linear scaling expression in reciprocal space for evaluating the ion--electron potential of crystalline solids. The expression replaces the long-range ion--electron potential with an equivalent localized charge distribution and corresponding boundary conditions on the unit cell. Given that no quadratic scaling structure factor is required---as used in traditional methods---the expression shows inherent linear behavior, and is well suited to simulating large-scale systems within orbital-free density functional theory. The scheme is implemented in the ATLAS software package and benchmarked by using a solid Mg bcc lattice containing tens of thousands of atoms in the unit cell. The test results show that the method can efficiently model large crystals with high computational accuracy.
تحميل البحث