Premium
Density‐functional method for very large systems with LCAO basis sets
Author(s) -
SánchezPortal Daniel,
Ordejón Pablo,
Artacho Emilio,
Soler José M.
Publication year - 1997
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/(sici)1097-461x(1997)65:5<453::aid-qua9>3.0.co;2-v
Subject(s) - orthogonalization , linear scale , linear combination of atomic orbitals , density functional theory , basis (linear algebra) , basis function , atomic orbital , orbital free density functional theory , wave function , density matrix , chemistry , statistical physics , quantum mechanics , scaling , hybrid functional , time dependent density functional theory , physics , computational chemistry , electron , basis set , mathematics , algorithm , geometry , geodesy , quantum , geography
We have implemented a linear scaling, fully self‐consistent density‐functional 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. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65 : 453–461, 1997