We present ONETEP (order-N electronic total energy package), a density functional program for parallel computers whose computational cost scales linearly with the number of atoms and the number of processors. ONETEP is based on our reformulation of the plane wave pseudopotential method which exploits the electronic localization that is inherent in systems with a nonvanishing band gap. We summarize the theoretical developments that enable the direct optimization of strictly localized quantities expressed in terms of a delocalized plane wave basis. These same localized quantities lead us to a physical way of dividing the computational effort among many processors to allow calculations to be performed efficiently on parallel supercomputers. We show with examples that ONETEP achieves excellent speedups with increasing numbers of processors and confirm that the time taken by ONETEP as a function of increasing number of atoms for a given number of processors is indeed linear. What distinguishes our approach is that the localization is achieved in a controlled and mathematically consistent manner so that ONETEP obtains the same accuracy as conventional cubic-scaling plane wave approaches and offers fast and stable convergence. We expect that calculations with ONETEP have the potential to provide quantitative theoretical predictions for problems involving thousands of atoms such as those often encountered in nanoscience and biophysics.

We present a reformulation of the plane-wave pseudopotential method for insulators. This new approach allows us to perform density-functional calculations by solving directly for ''nonorthogonal generalized Wannier functions'' rather than extended Bloch states. We outline the theory on which our method is based and present test calculations on a variety of systems. Comparison of our results with a standard plane-wave code shows that they are equivalent. Apart from the usual advantages of the plane-wave approach such as the applicability to any lattice symmetry and the high accuracy, our method also benefits from the localization properties of our functions in real space. The localization of all our functions greatly facilitates the future extension of our method to linear-scaling schemes or calculations of the electric polarization of crystalline insulators.

No abstract

Linear-scaling electronic structure methods are essential for calculations on large systems. Some of these approaches use a systematic basis set, the completeness of which may be tuned with an adjustable parameter similar to the energy cut-off of plane-wave techniques. The search for the electronic ground state in such methods suffers from an ill-conditioning which is related to the kinetic contribution to the total energy and which results in unacceptably slow convergence. We present a general preconditioning scheme to overcome this ill-conditioning and implement it within our own first-principles linear-scaling density functional theory method. The scheme may be applied in either real space or reciprocal space with equal success. The rate of convergence is improved by an order of magnitude and is found to be almost independent of the size of the basis.

In this work we study the solvatochromic shift of a selected low-energy excited state of alizarin in water by using a linear-scaling implementation of large-scale time-dependent density functional theory (TDDFT). While alizarin, a small organic dye, is chosen as a simple example of solute-solvent interactions, the findings presented here have wider ramifications for the realistic modeling of dyes, paints, and pigment-protein complexes. We find that about 380 molecules of explicit water need to be considered in order to yield an accurate representation of the solute-solvent interaction and a reliable solvatochromic shift. By using a novel method of constraining the TDDFT excitation vector, we confirm that the origin of the slow convergence of the solvatochromic shift with system size is due to two different effects. The first factor is a strong redshift of the excitation due to an explicit delocalization of a small fraction of the electron and the hole from the alizarin onto the water, which is mainly confined to within a distance of 7 Å from the alizarin molecule. The second factor can be identified as long-range electrostatic influences of water molecules beyond the 7 Å region on the ground-state properties of alizarin. We also show that these electrostatic influences are not well reproduced by a QM/MM model, suggesting that full QM studies of relatively large systems may be necessary in order to obtain reliable results.