Efficient and Accurate Methods for the Geometry Optimization of Water Clusters: Application of Analytic Gradients for the Two-Body:Many-Body QM:QM Fragmentation Method to (H2O)n, n = 3-10.

The structures of more than 70 low-lying water clusters ranging in size from (H2O)3 to (H2O)10 have been fully optimized with several different quantum mechanical electronic structure methods, including second-order Møller-Plesset perturbation theory (MP2) in conjunction with correlation consistent triple-ζ basis sets (aug-cc-pVTZ for O and cc-pVTZ for H, abbreviated haTZ). Optimized structures obtained with less demanding computational procedures were compared to the MP2/haTZ ones using both MP2/haTZ single point energies and the root-mean-square (RMS) deviations of unweighted Cartesian coordinates. Based on these criteria, B3LYP/6-31+G(d,2p) substantially outperforms both HF/haTZ and MP2/6-31G*. B3LYP/6-31+G(d,2p) structures never deviate from the MP2/haTZ geometries by more than 0.44 kcal mol(-1) on the MP2/haTZ potential energy surface, whereas the errors associated with the HF/haTZ and MP2/6-31G* structures grow as large as 12.20 and 2.98 kcal mol(-1), respectively. The most accurate results, however, were obtained with the two-body:many-body QM:QM fragmentation method for weakly bound clusters, in which all one- and two-body interactions are calculated at the high-level, while a low-level calculation is performed on the entire cluster to capture the cooperative effects (nonadditivity). With the haTZ basis set, the MP2:HF two-body:many-body fragmentation method generates structures that deviate from the MP2/haTZ ones by 0.01 kcal mol(-1) on average and not by more than 0.03 kcal mol(-1).

[1]  T. Steiner.  The hydrogen bond in the solid state. , 2002, Angewandte Chemie.

[2]  Brian W. Hopkins,et al.  Ab initio studies of π...π interactions: The effects of quadruple excitations , 2004 .

[3]  Gregory S. Tschumper,et al.  Anchoring the potential energy surface of the cyclic water trimer. , 2004, The Journal of chemical physics.

[4]  Yuriko Aoki,et al.  A theoretical synthesis of polymers by using uniform localization of molecular orbitals: Proposal of an elongation method , 1991 .

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

[6]  Yutaka Akiyama,et al.  Fragment molecular orbital method: application to polypeptides , 2000 .

[7]  Robert J. Harrison,et al.  Development of transferable interaction models for water. II. Accurate energetics of the first few water clusters from first principles , 2002 .

[8]  Gregory S. Tschumper.  Multicentered integrated QM:QM methods for weakly bound clusters: An efficient and accurate 2-body:many-body treatment of hydrogen bonding and van der Waals interactions , 2006 .

[9]  Gregory S. Tschumper,et al.  Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations. , 2008, The journal of physical chemistry. B.

[10]  Roger E. Miller,et al.  Molecular clusters: Structure and dynamics of weakly bound systems , 1996 .

[11]  Brian W. Hopkins,et al.  Anchoring the potential energy surface of the diacetylene dimer , 2009 .

[12]  Mark S. Gordon,et al.  Accurate methods for large molecular systems. , 2009, The journal of physical chemistry. B.

[13]  Gregory S. Tschumper,et al.  Analytic gradients for the multicentred integrated QM:QM method for weakly bound clusters: efficient and accurate 2-body:many-body geometry optimizations , 2007 .

[14]  Peter Pulay,et al.  Local configuration interaction: An efficient approach for larger molecules , 1985 .

[15]  Jing Ma,et al.  Linear scaling local correlation approach for solving the coupled cluster equations of large systems , 2002, J. Comput. Chem..

[16]  Shridhar R. Gadre,et al.  Molecular Tailoring Approach for Simulation of Electrostatic Properties , 1994 .

[17]  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.

[18]  Brian W. Hopkins,et al.  A multicentered approach to integrated QM/QM calculations. Applications to multiply hydrogen bonded systems , 2003, J. Comput. Chem..

[19]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[20]  Donald G Truhlar,et al.  Benchmark Databases for Nonbonded Interactions and Their Use To Test Density Functional Theory. , 2005, Journal of chemical theory and computation.

[21]  Nicholas J Mayhall,et al.  Molecules-in-Molecules: An Extrapolated Fragment-Based Approach for Accurate Calculations on Large Molecules and Materials. , 2011, Journal of chemical theory and computation.

[22]  Kazuya Ishimura,et al.  Accuracy of the three‐body fragment molecular orbital method applied to Møller–Plesset perturbation theory , 2007, J. Comput. Chem..

[23]  Wei Li,et al.  Local correlation calculations using standard and renormalized coupled-cluster approaches. , 2009, The Journal of chemical physics.

[24]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Correlation Energy, with Applications to the Calculation of Accurate Second-Order Møller-Plesset Perturbation Theory Energies for Large Water Clusters. , 2007, Journal of chemical theory and computation.

[25]  Pavel Hobza,et al.  2 Theoretical studies of hydrogen bonding , 2004 .

[26]  Shogo Sakai,et al.  Ab initio integrated multi-center molecular orbitals method for large cluster systems: total energy and normal vibration. , 2005, The journal of physical chemistry. A.

[27]  Brian W. Hopkins,et al.  Integrated quantum mechanical approaches for extended π systems: Multicentered QM/QM studies of the cyanogen and diacetylene trimers , 2005 .

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

[29]  Edward F. Valeev,et al.  Estimates of the Ab Initio Limit for π−π Interactions: The Benzene Dimer , 2002 .

[30]  A. W. Castleman,et al.  Clusters: Structure, Energetics, and Dynamics of Intermediate States of Matter , 1996 .

[31]  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.

[32]  Donald G Truhlar,et al.  Evaluation of the Electrostatically Embedded Many-Body Expansion and the Electrostatically Embedded Many-Body Expansion of the Correlation Energy by Application to Low-Lying Water Hexamers. , 2008, Journal of chemical theory and computation.

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

[34]  Ping Qian,et al.  Ab initio investigation of water clusters (H2O)n (n = 2–34) , 2010 .

[35]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters in the ground and excited states , 2005 .

[36]  M. Quack,et al.  Explicitly correlated coupled cluster calculations of the dissociation energies and barriers to concerted hydrogen exchange of (HF)n oligomers (n=2,3,4,5) , 1998 .

[37]  Brian W. Hopkins,et al.  Multicentred QM/QM methods for overlapping model systems , 2005 .

[38]  G. Shields,et al.  Accurate predictions of water cluster formation, (H₂O)(n=2-10). , 2010, The journal of physical chemistry. A.

[39]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters. , 2007, Journal of chemical theory and computation.

[40]  Nan Jiang,et al.  Electrostatic field-adapted molecular fractionation with conjugated caps for energy calculations of charged biomolecules. , 2006, The Journal of chemical physics.

[41]  P. Rossky,et al.  Nonadiabatic mixed quantum-classical dynamic simulation of pi-stacked oligophenylenevinylenes. , 2009, The journal of physical chemistry. A.

[42]  Ralf Ludwig,et al.  Water: From Clusters to the Bulk. , 2001, Angewandte Chemie.

[43]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[44]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[45]  Hannah R. Leverentz,et al.  Assessment of the accuracy of density functionals for prediction of relative energies and geometries of low-lying isomers of water hexamers. , 2008, The journal of physical chemistry. A.

[46]  Matthew L. Leininger,et al.  Anchoring the water dimer potential energy surface with explicitly correlated computations and focal point analyses , 2002 .

[47]  Guntram Rauhut,et al.  Integral transformation with low‐order scaling for large local second‐order Møller–Plesset calculations , 1998 .

[48]  M. Schütz,et al.  Ab Initio Calculations of the Binding Energies of Small (H2O)n Clusters (n = 1…4) , 1995 .

[49]  Wei Yang,et al.  A divide-and-conquer strategy to improve diffusion sampling in generalized ensemble simulations. , 2008, The Journal of chemical physics.

[50]  W. Klopper,et al.  Basis set limit coupled cluster study of h-bonded systems and assessment of more approximate methods. , 2007, The journal of physical chemistry. A.

[51]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[52]  Joshua R. Smith,et al.  Development of a 3-body:many-body integrated fragmentation method for weakly bound clusters and application to water clusters (H2O)(n = 3-10, 16, 17). , 2011, The Journal of chemical physics.