CONQUEST

CONQUEST is a linear scaling, or O(N), density functional theory (DFT) electronic structure open-source code.[1][2] The code is designed to perform DFT calculations on very large systems containing many thousands of atoms. It can be run at different levels of precision ranging from ab initio tight binding up to full DFT with plane wave accuracy. It has been applied to the study of three-dimensional reconstructions formed by Ge on Si(001), containing over 20,000 atoms.[3] Tests on the UK's national supercomputer HECToR in 2009 demonstrated the capability of the code to perform ground-state calculations on systems of over 1,000,000 atoms.

Methodology

Instead of solving for the Kohn-Sham eigenstates as normal DFT codes do, CONQUEST solves for the one particle density matrix, ${\displaystyle \rho (\mathbf {r} ,\mathbf {r} ^{\prime })}$. To make the problem computationally tractable, the density matrix is written in separable form:
${\displaystyle \rho (\mathbf {r} ,\mathbf {r} ^{\prime })=\sum _{i\alpha j\beta }\phi _{i\alpha }(\mathbf {r} )K_{i\alpha j\beta }\phi _{j\beta }(\mathbf {r} ^{\prime })}$, where ${\displaystyle \phi _{i\alpha }(\mathbf {r} )}$ is a support function centred on atom i (with support functions on the same atom notated by ${\displaystyle \alpha }$) and ${\displaystyle K_{i\alpha j\beta }}$ is the density matrix in the basis of the support functions.[4] The ground state is found as a series of nested loops:

• Minimise the energy with respect to the density matrix for fixed charge density and support functions
• Find self-consistency between charge density and potential
• Minimise the energy with respect to the support functions

The support functions are confined within spheres of given cutoff radius and the density matrix is forced to zero beyond a given range: ${\displaystyle K_{ij}=0,|\mathbf {R} _{i}-\mathbf {R} _{j}|>R_{c}}$. These approximations give linear scaling behaviour, and as the radii are increased tend to the exact result.

Developers

CONQUEST is jointly developed at the Department of Physics and Astronomy and London Centre for Nanotechnology, University College London in the UK and at the Computational Materials Science Centre, National Institute for Materials Science, Tsukuba, Japan. In the UK, the development team includes Dr. David Bowler, Dr. Veronika Brazdova, Prof. Mike Gillan, Dr. Andrew Horsfield, Mr. Alex Sena, Mr. Lianheng Tong, Mr. Jack Baker and Mr. Shereif Mujahed who are all members of the Thomas Young Centre; in Japan, the development team includes Dr. Tsuyoshi Miyazaki, Dr. Takahisa Ohno, Dr. Takao Ohtsuka, Dr. Milica Todorovic and Dr. Antonio Torralba.

Previous developers include Ian Bush, Rathin Choudhury, Chris Goringe and Eduardo Hernandez.

See also

• Density functional theory
• Quantum chemistry computer programs

External links

• CONQUEST official website
• National Institute for Materials Science website
• London Centre for Nanotechnology website

References

1. D. R. Bowler; T. Miyazaki; M. J. Gillan (2002). "Recent progress in linear-scaling ab initio electronic structure techniques". J. Phys. Condens. Matter. 14 (11): 2781–2798. Bibcode:2002JPCM...14.2781B. doi:10.1088/0953-8984/14/11/303.
2. D. R. Bowler; R. Choudhury; M. J. Gillan; T. Miyazaki (2006). "Recent progress with large-scale ab initio calculations: The CONQUEST code". Physica Status Solidi. 243 (5): 989–1000. arXiv:cond-mat/0603063. Bibcode:2006PSSBR.243..989B. doi:10.1002/pssb.200541386.
3. T. Miyazaki; D. R. Bowler; M. J. Gillan; T. Ohno (2008). "The energetics of hut-cluster self-assembly in Ge/Si(001) from linear-scaling DFT calculations". J. Phys. Soc. Jpn. 77 (12): 123706. Bibcode:2008JPSJ...77l3706M. doi:10.1143/JPSJ.77.123706.
4. E. Hernandez; M. J. Gillan (1995). "Self-consistent first-principles technique with linear scaling". Phys. Rev. B. 51 (15): 10157–10160. arXiv:mtrl-th/9501005. Bibcode:1995PhRvB..5110157H. doi:10.1103/PhysRevB.51.10157. PMID 9977692.