A parallel implementation of linear-scaling first-principles total-ene
rgy calculations is presented. The theoretical basis is density functi
onal theory (DFT) and the pseudopotential approximation. The linear-sc
aling method is the one due to the present authors, but the parallelis
ation techniques are also relevant to other linear-scaling DFT methods
. The theoretical and computational framework of the linear-scaling me
thod is summarised, in order to identify the main classes of computer
operation required. For each class of operation, the issues involved i
n distributing tasks and data between processors are discussed and a s
olution is proposed. Practical tests of the proposed implementation on
a Gray T3D are presented, and it is shown that the scaling with respe
ct to both the number of atoms and the number of processors is excelle
nt for systems containing up to over 6000 atoms.