We survey the underlying theory behind the large-scale and linear scaling DFT code, Conquest, which shows excellent parallel scaling and can be applied to thousands of atoms with diagonalisation, and millions of atoms with linear scaling. We give details of the representation of the density matrix and the approach to finding the electronic ground state, and discuss the implementation of molecular dynamics with linear scaling. We give an overview of the performance of the code, focussing in particular on the parallel scaling, and provide examples of recent developments and applications.

Given the widespread use of density functional theory (DFT), there is an increasing need for the ability to model large systems (beyond 1,000 atoms). We present a brief overview of the large-scale DFT code Conquest, which is capable of modelling such large systems, and discuss approaches to the generation of consistent, well-converged pseudo-atomic basis sets which will allow such large scale calculations. We present tests of these basis sets for a variety of materials, comparing to fully converged plane wave results using the same pseudopotentials and grids.

We have used density functional theory to study the energetics and electronic structure of aluminium dopants in crystalline silicon. We present data regarding the atomic and electronic structure and properties of pairs of substitutional aluminium dopants. We find that pairs of dopants, when occupying nearest neighbouring subsitutional sites in a high spin state, can bond to form aluminium pairs. This suggests that such a configuration of dopants will be electrically active when made to occupy a high spin state, whereas in the low spin state the neighbouring dopant pairs are found to be self compensating.