Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method.

Analytic second derivatives of the energy with respect to nuclear coordinates have been developed for spin restricted density functional theory (DFT) based on the fragment molecular orbital method (FMO). The derivations were carried out for the three-body expansion (FMO3), and the two-body expressions can be obtained by neglecting the three-body corrections. Also, the restricted Hartree-Fock (RHF) Hessian for FMO3 can be obtained by neglecting the density-functional related terms. In both the FMO-RHF and FMO-DFT Hessians, certain terms with small magnitudes are neglected for computational efficiency. The accuracy of the FMO-DFT Hessian in terms of the Gibbs free energy is evaluated for a set of polypeptides and water clusters and found to be within 1 kcal/mol of the corresponding full (non-fragmented) ab initio calculation. The FMO-DFT method is also applied to transition states in SN2 reactions and for the computation of the IR and Raman spectra of a small Trp-cage protein (PDB: 1L2Y). Some computational timing analysis is also presented.

[1]  Axel D. Becke,et al.  Density‐functional thermochemistry. IV. A new dynamical correlation functional and implications for exact‐exchange mixing , 1996 .

[2]  K. Kitaura,et al.  Analytic energy gradient for second-order Møller-Plesset perturbation theory based on the fragment molecular orbital method. , 2011, The Journal of chemical physics.

[3]  Yuto Komeiji,et al.  Three-body expansion and generalized dynamic fragmentation improve the fragment molecular orbital-based molecular dynamics (FMO-MD)☆ , 2010 .

[4]  Stephan Irle,et al.  Density-Functional Tight-Binding Combined with the Fragment Molecular Orbital Method. , 2014, Journal of chemical theory and computation.

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

[6]  John D. Head,et al.  COMPUTATION OF VIBRATIONAL FREQUENCIES FOR ADSORBATES ON SURFACES , 1997 .

[7]  Martin Karplus,et al.  Molecular properties from combined QM/MM methods. I. Analytical second derivative and vibrational calculations , 2000 .

[8]  Ye Mei,et al.  QUANTUM CALCULATION OF PROTEIN SOLVATION AND PROTEIN–LIGAND BINDING FREE ENERGY FOR HIV-1 PROTEASE/WATER COMPLEX , 2009 .

[9]  Noriyuki Kurita,et al.  Fragment molecular orbital method with density functional theory and DIIS convergence acceleration , 2003 .

[10]  Chris-Kriton Skylaris,et al.  Hybrid MPI-OpenMP Parallelism in the ONETEP Linear-Scaling Electronic Structure Code: Application to the Delamination of Cellulose Nanofibrils. , 2014, Journal of chemical theory and computation.

[11]  Jan H. Jensen,et al.  FragIt: A Tool to Prepare Input Files for Fragment Based Quantum Chemical Calculations , 2012, PloS one.

[12]  Ian M. Mills,et al.  Force Constants and Dipole-Moment Derivatives of Molecules from Perturbed Hartree-Fock Calculations. I , 1968 .

[13]  Shinichiro Nakamura,et al.  Efficient vibrational analysis for unrestricted Hartree–Fock based on the fragment molecular orbital method , 2014 .

[14]  T Verstraelen,et al.  Vibrational modes in partially optimized molecular systems. , 2007, The Journal of chemical physics.

[15]  M. Gordon,et al.  Accurate first principles model potentials for intermolecular interactions. , 2013, Annual review of physical chemistry.

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

[17]  Peng Bai,et al.  Water 26-mers Drawn from Bulk Simulations: Benchmark Binding Energies for Unprecedentedly Large Water Clusters and Assessment of the Electrostatically Embedded Three-Body and Pairwise Additive Approximations. , 2014, The journal of physical chemistry letters.

[18]  Hui Li,et al.  Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy , 2002 .

[19]  Robert Zaleśny,et al.  Linear-Scaling Techniques in Computational Chemistry and Physics , 2011 .

[20]  K. Kitaura,et al.  Unrestricted Hartree-Fock based on the fragment molecular orbital method: energy and its analytic gradient. , 2012, The Journal of chemical physics.

[21]  Filipp Furche,et al.  An efficient implementation of second analytical derivatives for density functional methods , 2002 .

[22]  Michael J. Frisch,et al.  Analytic second derivatives of the gradient-corrected density functional energy. Effect of quadrature weight derivatives , 1993 .

[23]  Kazuo Kitaura,et al.  Derivatives of the approximated electrostatic potentials in the fragment molecular orbital method , 2009 .

[24]  Jan H. Jensen,et al.  Analytic gradient for the adaptive frozen orbital bond detachment in the fragment molecular orbital method , 2009 .

[25]  Yuichi Inadomi,et al.  DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations. , 2007, The journal of physical chemistry. B.

[26]  Shinichiro Nakamura,et al.  Analytic second derivatives of the energy in the fragment molecular orbital method. , 2013, The Journal of chemical physics.

[27]  Kazuo Kitaura,et al.  The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems , 2009 .

[28]  Leo Radom,et al.  Harmonic Vibrational Frequencies: An Evaluation of Hartree−Fock, Møller−Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors , 1996 .

[29]  Kaori Fukuzawa,et al.  Fragment molecular orbital method: use of approximate electrostatic potential , 2002 .

[30]  Mark S Gordon,et al.  Structure and dynamics of the 1-hydroxyethyl-4-amino-1,2,4-triazolium nitrate high-energy ionic liquid system. , 2012, The journal of physical chemistry. B.

[31]  Masato Kobayashi,et al.  An effective energy gradient expression for divide-and-conquer second-order Møller-Plesset perturbation theory. , 2013, The Journal of chemical physics.

[32]  Yuto Komeiji,et al.  Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems. , 2014, Physical chemistry chemical physics : PCCP.

[33]  Kazuo Kitaura,et al.  Fully analytic energy gradient in the fragment molecular orbital method. , 2011, The Journal of chemical physics.

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

[35]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[36]  Mark S. Gordon,et al.  Parallel coupled perturbed CASSCF equations and analytic CASSCF second derivatives , 2006, J. Comput. Chem..

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

[38]  Michael W. Schmidt,et al.  Efficient molecular dynamics simulations of multiple radical center systems based on the fragment molecular orbital method. , 2014, The journal of physical chemistry. A.

[39]  Trygve Helgaker,et al.  Coupled cluster energy derivatives. Analytic Hessian for the closed‐shell coupled cluster singles and doubles wave function: Theory and applications , 1990 .

[40]  J. Ladik,et al.  Investigation of the interaction between molecules at medium distances: I. SCF LCAO MO supermolecule, perturbational and mutually consistent calculations for two interacting HF and CH2O molecules , 1975 .

[41]  J. Pople,et al.  Derivative studies in configuration–interaction theory , 1980 .

[42]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

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

[44]  Kenneth M Merz,et al.  Divide-and-Conquer Hartree-Fock Calculations on Proteins. , 2010, Journal of chemical theory and computation.

[45]  Mark S. Gordon,et al.  Parallel algorithm for SCF analytic Hessians. I. Small scale algorithm , 1993 .

[46]  Jiali Gao,et al.  Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations , 1997 .

[47]  Michael A Collins.  Molecular forces, geometries, and frequencies by systematic molecular fragmentation including embedded charges. , 2014, The Journal of chemical physics.

[48]  Jacob Kongsted,et al.  Ligand Affinities Estimated by Quantum Chemical Calculations. , 2010, Journal of chemical theory and computation.

[49]  Trygve Helgaker,et al.  Molecular Hessians for large‐scale MCSCF wave functions , 1986 .

[50]  Stephan Irle,et al.  Analytical second-order geometrical derivatives of energy for the self-consistent-charge density-functional tight-binding method. , 2004, The Journal of chemical physics.

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

[52]  Mark S. Gordon,et al.  Advances in electronic structure theory , 2005 .

[53]  Kazuo Kitaura,et al.  Exploring chemistry with the fragment molecular orbital method. , 2012, Physical chemistry chemical physics : PCCP.

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

[55]  K. Hirao,et al.  A long-range correction scheme for generalized-gradient-approximation exchange functionals , 2001 .

[56]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[57]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .

[58]  A. Becke.  Perspective: Fifty years of density-functional theory in chemical physics. , 2014, The Journal of chemical physics.

[59]  Y. Yamaguchi,et al.  A New Dimension to Quantum Chemistry: Analytic Derivative Methods in AB Initio Molecular Electronic Structure Theory , 1994 .

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

[61]  Ayako Nakata,et al.  Efficient Calculations with Multisite Local Orbitals in a Large-Scale DFT Code CONQUEST. , 2014, Journal of chemical theory and computation.

[62]  S. Asher,et al.  UV-resonance raman thermal unfolding study of Trp-cage shows that it is not a simple two-state miniprotein. , 2005, Journal of the American Chemical Society.

[63]  Yuri Alexeev,et al.  GAMESS as a free quantum-mechanical platform for drug research. , 2012, Current topics in medicinal chemistry.

[64]  T. Exner,et al.  Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 2. Level of Theory, Basis Set, and Solvents Model Dependence. , 2012, Journal of chemical theory and computation.

[65]  Kazuo Kitaura,et al.  Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method. , 2012, The Journal of chemical physics.

[66]  Kazuo Kitaura,et al.  Role of the key mutation in the selective binding of avian and human influenza hemagglutinin to sialosides revealed by quantum-mechanical calculations. , 2010, Journal of the American Chemical Society.

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

[68]  M. Gordon,et al.  Analytic Gradient for Density Functional Theory Based on the Fragment Molecular Orbital Method. , 2014, Journal of chemical theory and computation.

[69]  B. Brooks,et al.  Probing the local dynamics of nucleotide-binding pocket coupled to the global dynamics: myosin versus kinesin. , 2005, Biophysical journal.

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

[71]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[72]  M. Malagoli,et al.  The effect of grid quality and weight derivatives in density functional calculations of harmonic vibrational frequencies , 2003 .

[73]  Mark S. Gordon,et al.  Diffusion energy profiles in silica mesoporous molecular sieves modelled with the fragment molecular orbital method , 2013 .

[74]  Gregory S Tschumper,et al.  N-body:Many-body QM:QM vibrational frequencies: application to small hydrogen-bonded clusters. , 2013, The Journal of chemical physics.

[75]  K. Kitaura,et al.  Analytic gradient and molecular dynamics simulations using the fragment molecular orbital method combined with effective potentials , 2012, Theoretical Chemistry Accounts.

[76]  David E. Bernholdt,et al.  High performance computational chemistry: An overview of NWChem a distributed parallel application , 2000 .

[77]  Andrew Komornicki,et al.  Molecular gradients and hessians implemented in density functional theory , 1993 .

[78]  Kaori Fukuzawa,et al.  Modeling of peptide–silica interaction based on four-body corrected fragment molecular orbital (FMO4) calculations , 2013 .

[79]  J. Gauss,et al.  Analytic UHF-CCSD(T) second derivatives: implementation and application to the calculation of the vibration-rotation interaction constants of NCO and NCS , 1998 .

[80]  Mark Whittaker,et al.  Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method , 2011, J. Cheminformatics.

[81]  Kazuo Kitaura,et al.  On the accuracy of the 3-body fragment molecular orbital method (FMO) applied to density functional theory , 2004 .

[82]  Kazuo Kitaura,et al.  The three-body fragment molecular orbital method for accurate calculations of large systems , 2006 .

[83]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[84]  J. S. Binkley,et al.  Derivative studies in hartree-fock and møller-plesset theories , 2009 .

[85]  J. Gauss,et al.  Analytic CCSD(T) second derivatives , 1997 .

[86]  Kazuo Kitaura,et al.  Analytic gradient for the embedding potential with approximations in the fragment molecular orbital method , 2012 .

[87]  K. Kitaura,et al.  Simulations of Raman Spectra Using the Fragment Molecular Orbital Method. , 2014, Journal of chemical theory and computation.

[88]  Jon Baker,et al.  The effect of grid quality and weight derivatives in density functional calculations , 1994 .

[89]  Benny G. Johnson,et al.  A standard grid for density functional calculations , 1993 .

[90]  Lucas Visscher,et al.  Quantum-Chemical Electron Densities of Proteins and of Selected Protein Sites from Subsystem Density Functional Theory. , 2013, Journal of chemical theory and computation.

[91]  Stefan Goedecker,et al.  Linear scaling electronic structure methods in chemistry and physics , 2003, Comput. Sci. Eng..

[92]  Yoshihiro Matsumoto,et al.  Intrinsic Edge Asymmetry in Narrow Zigzag Hexagonal Heteroatomic Nanoribbons Causes their Subtle Uniform Curvature , 2012 .

[93]  Kwang S. Kim,et al.  Theory and applications of computational chemistry : the first forty years , 2005 .

[94]  C. W. Murray,et al.  Quadrature schemes for integrals of density functional theory , 1993 .

[95]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[96]  Mark S. Gordon,et al.  A parallel distributed data CPHF algorithm for analytic Hessians , 2007, J. Comput. Chem..

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