Molecular tailoring approach for geometry optimization of large molecules: energy evaluation and parallelization strategies.

A linear-scaling scheme for estimating the electronic energy, gradients, and Hessian of a large molecule at ab initio level of theory based on fragment set cardinality is presented. With this proposition, a general, cardinality-guided molecular tailoring approach (CG-MTA) for ab initio geometry optimization of large molecules is implemented. The method employs energy gradients extracted from fragment wave functions, enabling computations otherwise impractical on PC hardware. Further, the method is readily amenable to large scale coarse-grain parallelization with minimal communication among nodes, resulting in a near-linear speedup. CG-MTA is applied for density-functional-theory-based geometry optimization of a variety of molecules including alpha-tocopherol, taxol, gamma-cyclodextrin, and two conformations of polyglycine. In the tests performed, energy and gradient estimates obtained from CG-MTA during optimization runs show an excellent agreement with those obtained from actual computation. Accuracy of the Hessian obtained employing CG-MTA provides good hope for the application of Hessian-based geometry optimization to large molecules.

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

[2]  Shridhar R. Gadre,et al.  Tailoring approach for exploring electron densities and electrostatic potentials of molecular crystals , 2004 .

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

[4]  Paul G. Mezey,et al.  The Field-Adapted ADMA Approach: Introducing Point Charges , 2004 .

[5]  Donald E. Knuth,et al.  The Art of Computer Programming, Volume 4, Fascicle 2: Generating All Tuples and Permutations (Art of Computer Programming) , 2005 .

[6]  Kotoko Nakata,et al.  Ab initio quantum mechanical study of the binding energies of human estrogen receptor α with its ligands: An application of fragment molecular orbital method , 2005, J. Comput. Chem..

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

[8]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[9]  Wei Li,et al.  A localized molecular-orbital assembler approach for Hartree-Fock calculations of large molecules. , 2005, The Journal of chemical physics.

[10]  Mark S. Gordon,et al.  A new hierarchical parallelization scheme: Generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO) , 2004, J. Comput. Chem..

[11]  Steven L. Dixon,et al.  Parallel implementation of a divide and conquer semiempirical algorithm , 1998 .

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

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

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

[15]  Thom Vreven,et al.  Geometry optimization with QM/MM, ONIOM, and other combined methods. I. Microiterations and constraints , 2003, J. Comput. Chem..

[16]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[17]  Gustavo E. Scuseria,et al.  A quantitative study of the scaling properties of the Hartree–Fock method , 1995 .

[18]  Ching-Hsing Yu,et al.  Ab Initio Geometry Determinations of Proteins. 1. Crambin , 1998 .

[19]  Francesco Luigi Gervasio,et al.  Electronic structure of wet DNA. , 2002, Physical review letters.

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

[21]  Shridhar R Gadre,et al.  Many-body interaction analysis: algorithm development and application to large molecular clusters. , 2004, The Journal of chemical physics.

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

[23]  C. L. Liu.  Elements of Discrete Mathematics , 1985 .

[24]  Robert Wieczorek,et al.  Comparison of fully optimized alpha- and 3(10)-helices with extended beta-strands. An ONIOM density functional theory study. , 2004, Journal of the American Chemical Society.

[25]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

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

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

[28]  Henry F. Schaefer,et al.  Quantum chemistry in the 21st century (Special topic article) , 2000 .

[29]  Ranbir Singh,et al.  J. Mol. Struct. (Theochem) , 1996 .

[30]  Kazuo Kitaura,et al.  Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method. , 2004, The Journal of chemical physics.

[31]  Yuichi Inadomi,et al.  Fragment molecular orbital method: application to molecular dynamics simulation, ‘ab initio FMO-MD’ , 2003 .

[32]  Shridhar R. Gadre,et al.  Ab initio quality one‐electron properties of large molecules: Development and testing of molecular tailoring approach , 2003, J. Comput. Chem..