DCMB that combines divide‐and‐conquer and mixed‐basis set methods for accurate geometry optimizations, total energies, and vibrational frequencies of large molecules

We present a method, named DCMB, for the calculations of large molecules. It is a combination of a parallel divide‐and‐conquer (DC) method and a mixed‐basis (MB) set scheme. In this approach, atomic forces, total energy and vibrational frequencies are obtained from a series of MB calculations, which are derived from the target system utilizing the DC concept. Unlike the fragmentation based methods, all DCMB calculations are performed over the whole target system and no artificial caps are introduced so that it is particularly useful for charged and/or delocalized systems. By comparing the DCMB results with those from the conventional method, we demonstrate that DCMB is capable of providing accurate prediction of molecular geometries, total energies, and vibrational frequencies of molecules of general interest. We also demonstrate that the high efficiency of the parallel DCMB code holds the promise for a routine geometry optimization of large complex systems. © 2012 Wiley Periodicals, Inc.

[1]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[2]  Peter Pulay,et al.  Second-order Møller–Plesset calculations with dual basis sets , 2003 .

[3]  Yihan Shao,et al.  Dual-basis second-order Moller-Plesset perturbation theory: A reduced-cost reference for correlation calculations. , 2006, The Journal of chemical physics.

[4]  Michael A Collins,et al.  Approximate ab initio energies by systematic molecular fragmentation. , 2005, The Journal of chemical physics.

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

[6]  Masato Kobayashi,et al.  Implementation of divide‐and‐conquer method including Hartree‐Fock exchange interaction , 2007, J. Comput. Chem..

[7]  L. Curtiss,et al.  COMPUTATIONAL MATERIALS CHEMISTRY: METHODS AND APPLICATIONS , 2011 .

[8]  Anan Wu,et al.  XO: An extended ONIOM method for accurate and efficient geometry optimization of large molecules , 2010 .

[9]  Masato Kobayashi,et al.  Is the divide-and-conquer Hartree–Fock method valid for calculations of delocalized systems? , 2007 .

[10]  Hans Horn,et al.  Prescreening of two‐electron integral derivatives in SCF gradient and Hessian calculations , 1991 .

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

[12]  Masato Kobayashi,et al.  Second-order Møller-Plesset perturbation energy obtained from divide-and-conquer Hartree-Fock density matrix. , 2006, The Journal of chemical physics.

[13]  Yuto Komeiji,et al.  Fragment molecular orbital method: analytical energy gradients , 2001 .

[14]  H. Schlegel,et al.  Combining Synchronous Transit and Quasi-Newton Methods to Find Transition States , 1993 .

[15]  John W. Mintmire,et al.  Electronic and structural properties of carbon nanotubes , 1995 .

[16]  Martin Head-Gordon,et al.  Sparsity of the Density Matrix in Kohn-Sham Density Functional Theory and an Assessment of Linear System-Size Scaling Methods , 1997 .

[17]  Ryan P A Bettens,et al.  A new algorithm for molecular fragmentation in quantum chemical calculations. , 2006, The journal of physical chemistry. A.

[18]  Dieter Cremer,et al.  The Concept of the Chemical Bond , 1990 .

[19]  D. B. Chesnut,et al.  The use of locally dense basis sets in correlated NMR chemical shielding calculations , 1996 .

[20]  Wei Li,et al.  Geometry optimizations and vibrational spectra of large molecules from a generalized energy-based fragmentation approach. , 2008, The journal of physical chemistry. A.

[21]  Jiali Gao,et al.  Hybrid Quantum and Molecular Mechanical Simulations: An Alternative Avenue to Solvent Effects in Organic Chemistry , 1996 .

[22]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[23]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[24]  Ron Shepard,et al.  Elimination of the diagonalization bottleneck in parallel Direct-SCF methods , 1993 .

[25]  Wei Li,et al.  Divide-and-conquer local correlation approach to the correlation energy of large molecules. , 2004, The Journal of chemical physics.

[26]  Qingshi Zhu,et al.  Linear scaling calculation of band edge states and doped semiconductors. , 2007, The Journal of chemical physics.

[27]  Martin Head-Gordon,et al.  Approaching the Basis Set Limit in Density Functional Theory Calculations Using Dual Basis Sets without Diagonalization , 2004 .

[28]  M. Terrones.  Science and Technology of the Twenty-First Century: Synthesis, Properties, and Applications of Carbon Nanotubes , 2003 .

[29]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

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

[31]  Jan Řezáč,et al.  Multilevel Fragment-Based Approach (MFBA): A Novel Hybrid Computational Method for the Study of Large Molecules. , 2010, Journal of chemical theory and computation.

[32]  Paul G. Mezey,et al.  Ab Initio-Quality Electrostatic Potentials for Proteins: An Application of the ADMA Approach , 2002 .

[33]  Williams,et al.  N-scaling algorithm for density-functional calculations of metals and insulators. , 1994, Physical review. B, Condensed matter.

[34]  P. Jørgensen,et al.  Walking on potential energy surfaces , 1983 .

[35]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[36]  Tai-Sung Lee,et al.  A pseudobond approach to combining quantum mechanical and molecular mechanical methods , 1999 .

[37]  Shridhar R. Gadre,et al.  Molecular tailoring approach : towards PC-based ab initio treatment of large molecules , 2006 .

[38]  John Z H Zhang,et al.  New Advance in Computational Chemistry:  Full Quantum Mechanical ab Initio Computation of Streptavidin-Biotin Interaction Energy. , 2003, The journal of physical chemistry. B.

[39]  V Ganesh,et al.  Molecular tailoring approach for geometry optimization of large molecules: energy evaluation and parallelization strategies. , 2006, The Journal of chemical physics.

[40]  Kazuo Kitaura,et al.  The importance of three-body terms in the fragment molecular orbital method. , 2004, The Journal of chemical physics.

[41]  K. Hirao,et al.  A dual-level approach to density-functional theory. , 2006, Journal of Chemical Physics.

[42]  M. Frisch,et al.  Using redundant internal coordinates to optimize equilibrium geometries and transition states , 1996, J. Comput. Chem..

[43]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[44]  Mark S Gordon,et al.  Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies. , 2009, The journal of physical chemistry. A.

[45]  Gustavo E. Scuseria,et al.  Linear Scaling Density Functional Calculations with Gaussian Orbitals , 1999 .

[46]  Z. Maksić,et al.  Theoretical Models of Chemical Bonding , 1991 .

[47]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[48]  J. Almlöf,et al.  Dual basis sets in calculations of electron correlation , 1991 .

[49]  M. Head‐Gordon,et al.  Dual-basis self-consistent field methods: 6-31G* calculations with a minimal 6-4G primary basis , 2007 .

[50]  J. B. Collins,et al.  Self‐consistent molecular orbital methods. XVII. Geometries and binding energies of second‐row molecules. A comparison of three basis sets , 1976 .

[51]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[52]  G. A. Petersson,et al.  A complete basis set model chemistry. I. The total energies of closed‐shell atoms and hydrides of the first‐row elements , 1988 .

[53]  Manuel F. Ruiz-López,et al.  MODELING OF PEPTIDE HYDROLYSIS BY THERMOLYSIN. A SEMIEMPIRICAL AND QM/MM STUDY , 1998 .

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

[55]  L. Girifalco,et al.  Energy of Cohesion, Compressibility, and the Potential Energy Functions of the Graphite System , 1956 .

[56]  Wei Li,et al.  Generalized energy-based fragmentation approach for computing the ground-state energies and properties of large molecules. , 2007, The journal of physical chemistry. A.

[57]  Weitao Yang,et al.  Analytical energy gradients and geometry optimization in the divide‐and‐conquer method for large molecules , 1995 .

[58]  D. York,et al.  Linear‐scaling semiempirical quantum calculations for macromolecules , 1996 .

[59]  Philippe Sautet,et al.  Computational methods in catalysis and materials science , 2009 .

[60]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

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

[62]  Feng Long Gu,et al.  Efficiency and accuracy of the elongation method as applied to the electronic structures of large systems , 2006, J. Comput. Chem..

[63]  Wei Pan,et al.  Parallel implementation of divide‐and‐conquer semiempirical quantum chemistry calculations , 1998 .

[64]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

[65]  M. Head‐Gordon,et al.  Dual-basis analytic gradients. 1. Self-consistent field theory. , 2006, The journal of physical chemistry. A.

[66]  Thom Vreven,et al.  Chapter 3 Hybrid Methods: ONIOM(QM:MM) and QM/MM , 2006 .

[67]  Wei Li,et al.  An efficient fragment-based approach for predicting the ground-state energies and structures of large molecules. , 2005, Journal of the American Chemical Society.

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

[69]  J. Pople,et al.  Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements , 2002 .

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

[71]  Masato Kobayashi,et al.  Implementation of Surjan's density matrix formulae for calculating second-order Møller-Plesset energy , 2006 .

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

[73]  Yi Luo,et al.  An elongation method for first principle simulations of electronic structures and electron transport properties of finite nanostructures. , 2006, The Journal of chemical physics.

[74]  Anuja P. Rahalkar,et al.  Enabling ab initio Hessian and frequency calculations of large molecules. , 2008, The Journal of chemical physics.

[75]  Tai-Sung Lee,et al.  Linear-scaling quantum mechanical calculations of biological molecules: The divide-and-conquer approach , 1998 .