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

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

[1]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[2]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[3]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[4]  Steven G. Louie,et al.  Nonlinear ionic pseudopotentials in spin-density-functional calculations , 1982 .

[5]  Leonard Kleinman,et al.  Efficacious Form for Model Pseudopotentials , 1982 .

[6]  Wolfram Saenger,et al.  Principles of Nucleic Acid Structure , 1983 .

[7]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[8]  H. Schaefer,et al.  Extensive theoretical studies of the hydrogen‐bonded complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2 , 1986 .

[9]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[10]  O. Sankey,et al.  Ab initio multicenter tight-binding model for molecular-dynamics simulations and other applications in covalent systems. , 1989, Physical review. B, Condensed matter.

[11]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[12]  Galli,et al.  Large scale electronic structure calculations. , 1992, Physical review letters.

[13]  Dennis R. Salahub,et al.  Gaussian density functional calculations on hydrogen-bonded systems , 1992 .

[14]  Car,et al.  Orbital formulation for electronic-structure calculations with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[15]  Barnett,et al.  Born-Oppenheimer molecular-dynamics simulations of finite systems: Structure and dynamics of (H2O)2. , 1993, Physical review. B, Condensed matter.

[16]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[17]  Soler,et al.  Augmented-plane-wave calculations on small molecules. , 1993, Physical review. B, Condensed matter.

[18]  Li,et al.  Density-matrix electronic-structure method with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[19]  Martin,et al.  Unconstrained minimization approach for electronic computations that scales linearly with system size. , 1993, Physical review. B, Condensed matter.

[20]  Stechel,et al.  Order-N methods in self-consistent density-functional calculations. , 1994, Physical review. B, Condensed matter.

[21]  Peter A. Kollman,et al.  Theoretical Investigation of the Hydrogen Bond Strengths in Guanine-Cytosine and Adenine-Thymine Base Pairs , 1994 .

[22]  Hernández,et al.  Self-consistent first-principles technique with linear scaling. , 1995, Physical review. B, Condensed matter.

[23]  Kim,et al.  Total-energy global optimizations using nonorthogonal localized orbitals. , 1995, Physical review. B, Condensed matter.

[24]  Martin,et al.  Linear system-size scaling methods for electronic-structure calculations. , 1995, Physical review. B, Condensed matter.

[25]  Mixed Approach to Incorporate Self-Consistency into Order-N LCAO Methods , 1995 .

[26]  J. Šponer,et al.  Structures and Energies of Hydrogen-Bonded DNA Base Pairs. A Nonempirical Study with Inclusion of Electron Correlation , 1996 .

[27]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[28]  Soler,et al.  Self-consistent order-N density-functional calculations for very large systems. , 1996, Physical review. B, Condensed matter.

[29]  D. Sánchez-Portal,et al.  Analysis of atomic orbital basis sets from the projection of plane-wave results , 1995, cond-mat/9509053.

[30]  Otto F. Sankey,et al.  ELECTRONIC-STRUCTURE-BASED MOLECULAR-DYNAMICS METHOD FOR LARGE BIOLOGICAL SYSTEMS : APPLICATION TO THE 10 BASEPAIR POLY(DG).POLY(DC) DNA DOUBLE HELIX , 1997 .