Density-functional method for very large systems with LCAO basis sets

✍ Scribed by Daniel Sánchez-Portal; Pablo Ordejón; Emilio Artacho; José M. Soler

• Publisher: John Wiley and Sons
• Year: 1997
• Tongue: English
• Weight: 181 KB
• Volume: 65
• Category: Article
• ISSN: 0020-7608
• DOI: 10.1002/(sici)1097-461x(1997)65:5<453::aid-qua9>3.0.co;2-v

✦ Synopsis

We have implemented a linear scaling, fully self-consistent densityfunctional method for performing first-principles calculations on systems with a large number of atoms, using standard norm-conserving pseudopotentials and flexible linear Ž . combinations of atomic orbitals LCAO basis sets. Exchange and correlation are treated within the local-spin-density or gradient-corrected approximations. The basis functions and the electron density are projected on a real-space grid in order to calculate the Hartree and exchange᎐correlation potentials and matrix elements. We substitute the customary diagonalization procedure by the minimization of a modified energy functional, which gives orthogonal wave functions and the same energy and density as the Kohn᎐Sham energy functional, without the need of an explicit orthogonalization. The additional restriction to a finite range for the electron wave functions allows the Ž . computational effort time and memory to increase only linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, allowing structural relaxation and molecular dynamics simulations. We present test calculations beginning with small molecules and ending with a piece of DNA. Using double-z, polarized bases, geometries within 1% of experiments are obtained.