Linear-scaling density-functional theory with tens of thousands of atoms: Expanding the scope and scale of calculations with ONETEP

ONETEP is an ab initio electronic structure package for total energy calculations within density-functional theory. It combines 'linear scaling', in that the total computational effort scales only linearly with system size, with 'plane-wave' accuracy, in that the convergence of the total energy is systematically improvable in the manner typical of conventional plane-wave pseudopotential methods. We present recent progress on improving the performance, and thus in effect the feasible scope and scale, of calculations with ONETEP on parallel computers comprising large clusters of commodity servers. Our recent improvements make calculations of tens of thousands of atoms feasible, even on fewer than 100 cores. Efficient scaling with number of atoms and number of cores is demonstrated Lip to 32,768 atoms on 64 cores. (C) 2008 Elsevier B.V. All rights reserved.