Car-Parrinello density matrix search with a first principles fictitious electron mass method for electronic wave function optimization.

In spite of its success in molecular dynamics and the advantage of being a first order propagation technique, the Car-Parrinello method and its variations have not been successful in self-consistent-field (SCF) wave function optimization due to slow convergence. In this article, we introduce a first principles fictitious mass scheme to weigh each individual density element differently and instantaneously. As an alternative to diagonalization in SCF, the Car-Parrinello scheme is implemented as a density matrix search method. Not only does the fictitious mass scheme developed herein allow a very fast SCF convergence, but also the Car-Parrinello density matrix search (CP-DMS) exhibits linear scaling with respect to the system size for alanine helical chain test molecules. The excellent performance of CP-DMS holds even for very challenging compact three-dimensional quantum particles. While the conventional diagonalization based SCF method has difficulties optimizing electronic wave functions for CdSe quantum dots, CP-DMS shows both smooth and faster convergence.

[1]  Todd J. Martinez,et al.  Graphical Processing Units for Quantum Chemistry , 2008, Computing in Science & Engineering.

[2]  D. Gamelin,et al.  Investigation of pure and Co2+-doped ZnO quantum dot electronic structures using the density functional theory: choosing the right functional , 2008 .

[3]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation. , 2008, Journal of chemical theory and computation.

[4]  Michael J Frisch,et al.  Efficient evaluation of short-range Hartree-Fock exchange in large molecules and periodic systems. , 2006, The Journal of chemical physics.

[5]  S. Iyengar.  Dynamical effects on vibrational and electronic spectra of hydroperoxyl radical water clusters. , 2005, The Journal of chemical physics.

[6]  M. Frisch,et al.  Effect of time-dependent basis functions and their superposition error on atom-centered density matrix propagation (ADMP): connections to wavelet theory of multiresolution analysis. , 2004, The Journal of chemical physics.

[7]  Gustavo E Scuseria,et al.  Efficient hybrid density functional calculations in solids: assessment of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional. , 2004, The Journal of chemical physics.

[8]  Gustavo E Scuseria,et al.  Assessment and validation of a screened Coulomb hybrid density functional. , 2004, The Journal of chemical physics.

[9]  Michael J. Frisch,et al.  Density matrix search using direct inversion in the iterative subspace as a linear scaling alternative to diagonalization in electronic structure calculations , 2003 .

[10]  G. Scuseria,et al.  Hybrid functionals based on a screened Coulomb potential , 2003 .

[11]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with gaussian orbitals. IV. Formal analysis of the deviations from born‐oppenheimer dynamics , 2002 .

[12]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. III. Comparison with Born–Oppenheimer dynamics , 2002 .

[13]  G. Scuseria,et al.  A black-box self-consistent field convergence algorithm: One step closer , 2002 .

[14]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy conservation and choice of initial conditions , 2001 .

[15]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals , 2001 .

[16]  Gustavo E. Scuseria,et al.  Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations , 1997 .

[17]  Michael J. Frisch,et al.  A linear scaling method for Hartree–Fock exchange calculations of large molecules , 1996 .

[18]  Eric Schwegler,et al.  Linear scaling computation of the Hartree–Fock exchange matrix , 1996 .

[19]  Benny G. Johnson,et al.  Linear scaling density functional calculations via the continuous fast multipole method , 1996 .

[20]  Martin Head-Gordon,et al.  Derivation and efficient implementation of the fast multipole method , 1994 .

[21]  E. Kryachko.  On the Car–Parrinello computational scheme: Rigorous treatment , 1993 .

[22]  E. Kryachko.  On the formalism of many-electron dynamics , 1993 .

[23]  Kryachko.  Formalism for molecular-dynamics calculations of many-electron systems. , 1993, Physical review. B, Condensed matter.

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

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

[26]  R. Mcweeny.  Some Recent Advances in Density Matrix Theory , 1960 .

[27]  G. Golub.  Matrix computations , 1983 .

[28]  Per-Olov Löwdin,et al.  On the Nonorthogonality Problem , 1970 .