Density-Functional Tight-Binding Combined with the Fragment Molecular Orbital Method.

We developed the energy and its gradient for the self-consistent-charge density-functional tight-binding (DFTB) method, combined with the fragment molecular orbital (FMO) approach, FMO-DFTB, including an optional a posteriori treatment for dispersion interaction, and evaluated its accuracy as well as computational efficiency for a set of representative systems: polypeptides, a DNA segment, and a small protein. The error in the total energy of FMO-DFTB versus full SCC-DFTB was below 1 kcal/mol for the polyalanine system consisting of about 2000 atoms partitioned into fragments containing 2 residues, and the optimized structures had root-mean-square deviations below 0.1 Å. The scaling of FMO-DFTB with the system size N is only marginally larger than linear [O(N(1.2)) in the worst case]. A parallelization efficiency of 94% was achieved using 128 CPU cores, and we demonstrate the applicability of FMO-DFTB for systems containing more than one million atoms by performing a geometry optimization of a fullerite cluster.

[1]  Serguei Patchkovskii,et al.  An Efficient a Posteriori Treatment for Dispersion Interaction in Density-Functional-Based Tight Binding. , 2005, Journal of chemical theory and computation.

[2]  K. Kitaura,et al.  Derivatives of the approximated electrostatic potentials in unrestricted Hartree–Fock based on the fragment molecular orbital method and an application to polymer radicals , 2014, Theoretical Chemistry Accounts.

[3]  Masami Uebayasi,et al.  The fragment molecular orbital method for geometry optimizations of polypeptides and proteins. , 2007, The journal of physical chemistry. A.

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

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

[6]  K. Morokuma,et al.  Molecular and electronic structures of endohedral fullerenes, Sc2C2@C3v–C82 and Sc2@C3v–C82: Benchmark for SCC‐DFTB and proposal of new inner cluster structures , 2012 .

[7]  Kazuo Kitaura,et al.  Importance of the hybrid orbital operator derivative term for the energy gradient in the fragment molecular orbital method , 2010 .

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

[9]  Tai-Sung Lee,et al.  A variational linear-scaling framework to build practical, efficient next-generation orbital-based quantum force fields. , 2013, Journal of chemical theory and computation.

[10]  Miguel Fuentes-Cabrera,et al.  The application of approximate density functionals to complex systems , 1998 .

[11]  James J. P. Stewart,et al.  Application of the PM6 method to modeling proteins , 2009, Journal of molecular modeling.

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

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

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

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

[16]  G. Seifert,et al.  Calculations of molecules, clusters, and solids with a simplified LCAO-DFT-LDA scheme , 1996 .

[17]  Michael A Collins,et al.  Systematic fragmentation of large molecules by annihilation. , 2012, Physical chemistry chemical physics : PCCP.

[18]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

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

[20]  Yuichi Inadomi,et al.  Fragment molecular orbital study of the electronic excitations in the photosynthetic reaction center of Blastochloris viridis , 2009, J. Comput. Chem..

[21]  Lou Massa,et al.  Kernel energy method applied to vesicular stomatitis virus nucleoprotein , 2009, Proceedings of the National Academy of Sciences.

[22]  Gotthard Seifert,et al.  Density‐functional tight binding—an approximate density‐functional theory method , 2012 .

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

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

[25]  Anders M N Niklasson,et al.  Energy conserving, linear scaling Born-Oppenheimer molecular dynamics. , 2012, The Journal of chemical physics.

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

[27]  T. Nakano,et al.  Examination of numerical accuracy on fragment-DFT calculations with integral values of total electron density functions , 2006 .

[28]  K. Kitaura,et al.  Time-dependent density functional theory with the multilayer fragment molecular orbital method , 2007 .

[29]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[30]  Pekka Koskinen,et al.  Density-functional tight-binding for beginners , 2009, 0910.5861.

[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]  Hélio A. Duarte,et al.  Density-functional based tight-binding: an approximate DFT method , 2009 .

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

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

[35]  D. Bowler,et al.  Linear Scaling Constrained Density Functional Theory in CONQUEST. , 2011, Journal of chemical theory and computation.

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

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

[38]  Steve Kaminski,et al.  Improved electronic properties from third-order SCC-DFTB with cost efficient post-SCF extensions. , 2012, The journal of physical chemistry. A.

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

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

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

[42]  Lou Massa,et al.  Quantum kernel applications in medicinal chemistry. , 2012, Future medicinal chemistry.

[43]  Kazuo Kitaura,et al.  Binding of influenza A virus hemagglutinin to the sialoside receptor is not controlled by the homotropic allosteric effect. , 2010, The journal of physical chemistry. B.

[44]  Thomas Frauenheim,et al.  Atomistic simulations of complex materials: ground-state and excited-state properties , 2002 .

[45]  T. Frauenheim,et al.  A density-functional based tight-binding approach to GaAs surface reconstructions , 1997 .

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

[47]  Jiali Gao,et al.  Communication: variational many-body expansion: accounting for exchange repulsion, charge delocalization, and dispersion in the fragment-based explicit polarization method. , 2012, The Journal of chemical physics.

[48]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[49]  Johnson,et al.  Modified Broyden's method for accelerating convergence in self-consistent calculations. , 1988, Physical review. B, Condensed matter.

[50]  N. Mel’nik,et al.  Metastable crystalline and amorphous carbon phases obtained from fullerite C 60 by high-pressure–high-temperature treatment , 1997 .

[51]  Hui Li,et al.  Energy gradients in combined fragment molecular orbital and polarizable continuum model (FMO/PCM) calculation , 2009, J. Comput. Chem..

[52]  D. York,et al.  Extension of the self-consistent-charge density-functional tight-binding method: third-order expansion of the density functional theory total energy and introduction of a modified effective coulomb interaction. , 2007, The journal of physical chemistry. A.

[53]  Klaus Schulten,et al.  Mature HIV-1 capsid structure by cryo-electron microscopy and all-atom molecular dynamics , 2013, Nature.

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

[55]  Foulkes,et al.  Tight-binding models and density-functional theory. , 1989, Physical review. B, Condensed matter.

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

[57]  Paul G. Mezey,et al.  A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions , 1989 .

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

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

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

[61]  Shridhar R. Gadre,et al.  Appraisal of molecular tailoring approach for large clusters. , 2013, The Journal of chemical physics.

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

[63]  Stephan Irle,et al.  The C60 formation puzzle "solved": QM/MD simulations reveal the shrinking hot giant road of the dynamic fullerene self-assembly mechanism. , 2006, The journal of physical chemistry. B.

[64]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[65]  Thomas Frauenheim,et al.  Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment , 2001 .

[66]  Michael Gaus,et al.  DFTB3: Extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB). , 2011, Journal of chemical theory and computation.

[67]  Vladimir B. Sulimov,et al.  Semiempirical calculations of binding enthalpy for protein-ligand complexes , 2004 .

[68]  K. Kitaura,et al.  Use of an auxiliary basis set to describe the polarization in the fragment molecular orbital method , 2014 .

[69]  Donald G Truhlar,et al.  X-Pol Potential: An Electronic Structure-Based Force Field for Molecular Dynamics Simulation of a Solvated Protein in Water. , 2009, Journal of chemical theory and computation.

[70]  Hao Hu,et al.  Simulating water with the self-consistent-charge density functional tight binding method: from molecular clusters to the liquid state. , 2007, The journal of physical chemistry. A.

[71]  Wei Li,et al.  Generalized Energy-Based Fragmentation CCSD(T)-F12a Method and Application to the Relative Energies of Water Clusters (H2O)20. , 2014, Journal of chemical theory and computation.

[72]  T. Hirano,et al.  A third-generation density-functional-theory-based method for calculating canonical molecular orbitals of large molecules. , 2014, Physical chemistry chemical physics : PCCP.

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

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

[75]  Kaori Fukuzawa,et al.  Development of the four-body corrected fragment molecular orbital (FMO4) method , 2012 .

[76]  Stephan Irle,et al.  Performance of the DFTB method in comparison to DFT and semiempirical methods for geometries and energies of C20–C86 fullerene isomers , 2005 .

[77]  N. Hush,et al.  Density-functional geometry optimization of the 150,000-atom photosystem-I trimer. , 2006, The Journal of chemical physics.

[78]  R. Parr,et al.  Absolute hardness: companion parameter to absolute electronegativity , 1983 .

[79]  Yuri Alexeev,et al.  Geometry Optimization of the Active Site of a Large System with the Fragment Molecular Orbital Method , 2011 .

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

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

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

[83]  Michael Gaus,et al.  Density functional tight binding: application to organic and biological molecules , 2014 .

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

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

[86]  Seifert,et al.  Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon. , 1995, Physical review. B, Condensed matter.

[87]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

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

[89]  K. Kitaura,et al.  Systematic study of the embedding potential description in the fragment molecular orbital method. , 2010, The journal of physical chemistry. A.

[90]  Michael Methfessel,et al.  Crystal structures of zirconia from first principles and self-consistent tight binding , 1998 .

[91]  Kaori Fukuzawa,et al.  Large-scale FMO-MP3 calculations on the surface proteins of influenza virus, hemagglutinin (HA) and neuraminidase (NA) , 2010 .

[92]  Shigenori Tanaka,et al.  Intra‐ and intermolecular interactions between cyclic‐AMP receptor protein and DNA: Ab initio fragment molecular orbital study , 2006, J. Comput. Chem..

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

[94]  Takeshi Ishikawa,et al.  How does an S(N)2 reaction take place in solution? Full ab initio MD simulations for the hydrolysis of the methyl diazonium ion. , 2008, Journal of the American Chemical Society.

[95]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .