Embedded divide-and-conquer algorithm on hierarchical real-space grids: parallel molecular dynamics simulation based on linear-scaling density functional theory

A linear-scaling algorithm has been developed to perform large-scale molecular-dynamics (MD) simulations, in which interatomic forces are computed quantum mechanically in the framework of the density functional theory. A divide-and-conquer algorithm is used to compute the electronic structure, where non-additive contribution to the kinetic energy is included with an embedded cluster scheme. Electronic wave functions are represented on a real-space grid, which is augmented with coarse multigrids to accelerate the convergence of iterative solutions and adaptive fine grids around atoms to accurately calculate ionic pseudopotentials. Spatial decomposition is employed to implement the hierarchical-grid algorithm on massively parallel computers. A converged solution to the electronic-structure problem is obtained for a 32,768-atom amorphous CdSe system on 512 IBM POWER4 processors. The total energy is well conserved during MD simulations of liquid Rb, showing the applicability of this algorithm to first principles MD simulations. The parallel efficiency is 0.985 on 128 Intel Xeon processors for a 65,536-atom CdSe system.

[1]  Rajiv K. Kalia,et al.  Scalable and portable implementation of the fast multipole method on parallel computers , 2003 .

[2]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[3]  William H. Press,et al.  Numerical recipes , 1990 .

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

[5]  Frederick H. Streitz,et al.  Electrostatic potentials for metal-oxide surfaces and interfaces. , 1994 .

[6]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .

[7]  Peter S. Lomdahl,et al.  Recent advances in large-scale atomistic materials simulations , 1999, Comput. Sci. Eng..

[8]  Arieh Warshel,et al.  Ab Initio Free Energy Perturbation Calculations of Solvation Free Energy Using the Frozen Density Functional Approach , 1994 .

[9]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[10]  A. Nakano,et al.  Multiresolution molecular dynamics algorithm for realistic materials modeling on parallel computers , 1994 .

[11]  White,et al.  Detonations at nanometer resolution using molecular dynamics. , 1993, Physical review letters.

[12]  Rajiv K. Kalia,et al.  Linear-scaling density-functional-theory calculations of electronic structure based on real-space grids: design, analysis, and scalability test of parallel algorithms , 2001 .

[13]  T. Beck.  Real-space mesh techniques in density-functional theory , 2000, cond-mat/0006239.

[14]  Timothy Campbell,et al.  A scalable molecular-dynamics algorithm suite for materials simulations: design-space diagram on 1024 Cray T3E processors , 2000, Future Gener. Comput. Syst..

[15]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[16]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[17]  N. Govind,et al.  Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems , 1999 .

[18]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[20]  Jean-Luc Fattebert,et al.  Linear scaling first-principles molecular dynamics with controlled accuracy , 2004, Comput. Phys. Commun..

[21]  Rajiv K. Kalia,et al.  Short- and intermediate-range structural correlations in amorphous silicon carbide: a molecular dynamics study , 2004 .

[22]  James R. Chelikowsky,et al.  Electronic Structure Methods for Predicting the Properties of Materials: Grids in Space , 2000 .

[23]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[24]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

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

[26]  Rajiv K. Kalia,et al.  Massively parallel algorithms for computational nanoelectronics based on quantum molecular dynamics , 1994 .

[27]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[28]  Robert Walkup,et al.  Simulating materials failure by using up to one billion atoms and the world's fastest computer: Brittle fracture , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[29]  J. Tersoff,et al.  New empirical approach for the structure and energy of covalent systems. , 1988, Physical review. B, Condensed matter.

[30]  J. Fattebert,et al.  Towards grid-based OÑNÖ density-functional theory methods: Optimized nonorthogonal orbitals and multigrid acceleration , 2000 .

[31]  Donald W. Brenner,et al.  The Art and Science of an Analytic Potential , 2000 .

[32]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. II. Operators for fast iterative diagonalization. , 1991, Physical review. B, Condensed matter.

[33]  Lin,et al.  Real-space implementation of nonlocal pseudopotentials for first-principles total-energy calculations. , 1991, Physical review. B, Condensed matter.

[34]  Rajiv K. Kalia,et al.  DYNAMICS OF OXIDATION OF ALUMINUM NANOCLUSTERS USING VARIABLE CHARGE MOLECULAR-DYNAMICS SIMULATIONS ON PARALLEL COMPUTERS , 1999 .

[35]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[36]  S. L. Dixon,et al.  Fast, accurate semiempirical molecular orbital calculations for macromolecules , 1997 .

[37]  David E. Keyes,et al.  Domain-Based Parallelism and Problem Decomposition Methods in Computational Science and Engineering , 1995 .

[38]  A. V. van Duin,et al.  Shock waves in high-energy materials: the initial chemical events in nitramine RDX. , 2003, Physical review letters.

[39]  Rajiv K. Kalia,et al.  Molecular dynamics simu-lations of Coulombic systems on ditributed-memory MIMD machines , 1993 .

[40]  Tomoya Ono,et al.  Timesaving Double-Grid Method for Real-Space Electronic-Structure Calculations , 1999 .

[41]  Weitao Yang,et al.  A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules , 1995 .

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

[43]  Rajiv K. Kalia,et al.  Hybrid quantum mechanical/molecular dynamics simulation on parallel computers: density functional theory on real-space multigrids , 2002 .

[44]  Sullivan,et al.  Real-space multigrid-based approach to large-scale electronic structure calculations. , 1996, Physical review. B, Condensed matter.

[45]  Subhash Saini,et al.  Scalable atomistic simulation algorithms for materials research , 2002 .

[46]  Rajiv K. Kalia,et al.  Environmental effects of H2O on fracture initiation in silicon: A hybrid electronic-density-functional/molecular-dynamics study , 2004 .

[47]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[48]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[49]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[50]  Rajiv K. Kalia,et al.  Hybrid finite-element/molecular-dynamics/electronic-density-functional approach to materials simulations on parallel computers , 2001 .

[51]  Aiichiro Nakano,et al.  Parallel multilevel preconditioned conjugate-gradient approach to variable-charge molecular dynamics , 1997 .

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

[53]  Wu,et al.  Higher-order finite-difference pseudopotential method: An application to diatomic molecules. , 1994, Physical review. B, Condensed matter.