Premium
Implementation of linear‐scaling plane wave density functional theory on parallel computers
Author(s) -
Skylaris ChrisKriton,
Haynes Peter D.,
Mostofi Arash A.,
Payne Mike C.
Publication year - 2006
Publication title -
physica status solidi (b)
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.51
H-Index - 109
eISSN - 1521-3951
pISSN - 0370-1972
DOI - 10.1002/pssb.200541328
Subject(s) - computer science , scaling , computation , subroutine , parallel computing , basis (linear algebra) , range (aeronautics) , message passing interface , linear scale , computational science , set (abstract data type) , plane (geometry) , message passing , algorithm , theoretical computer science , mathematics , materials science , geometry , geodesy , composite material , programming language , geography , operating system
We describe the algorithms we have developed for linear‐scaling plane wave density functional calculations on parallel computers as implemented in the onetep program. We outline how onetep achieves plane wave accuracy with a computational cost which increases only linearly with the number of atoms by optimising directly the single‐particle density matrix expressed in a psinc basis set. We describe in detail the novel algorithms we have developed for computing with the psinc basis set the quantities needed in the evaluation and optimisation of the total energy within our approach. For our parallel computations we use the general Message Passing Interface (MPI) library of subroutines to exchange data between processors. Accordingly, we have developed efficient schemes for distributing data and computational load to processors in a balanced manner. We describe these schemes in detail and in relation to our algorithms for computations with a psinc basis. Results of tests on different materials show that onetep is an efficient parallel code that should be able to take advantage of a wide range of parallel computer architectures. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)