Fragment molecular orbital
The fragment molecular orbital method (FMO) is a computational method that can be used to calculate very large molecular systems with thousands of atoms using ab initio quantum-chemical wave functions.
History of FMO and related methods
The fragment molecular orbital method (FMO) was developed by Kazuo Kitaura and coworkers in 1999.[1] FMO is deeply interconnected with the energy decomposition analysis (EDA) by Kazuo Kitaura and Keiji Morokuma, developed in 1976. The main use of FMO is to compute very large molecular systems by dividing them into fragments and performing ab initio or density functional quantum-mechanical calculations of fragments and their dimers, whereby the Coulomb field from the whole system is included. The latter feature allows fragment calculations without using caps.
The mutually consistent field (MCF) method [2] had introduced the idea of self-consistent fragment calculations in their embedding potential, which was later used with some modifications in various methods including FMO. There had been other methods related to FMO including the incremental correlation method by H. Stoll (1992).[3]
Later, other methods closely related to FMO were proposed including the kernel energy method of L. Huang[4] and the electrostatically embedded many-body expansion by E. Dahlke,[5] S. Hirata[6] and later M. Kamiya[7] suggested approaches also very closely related to FMO. Effective fragment molecular orbital (EFMO) method combines some features of the effective fragment potentials (EFP) and FMO. A detailed perspective on the fragment-based method development can be found in a review.[8]
Introduction to FMO
In addition to the calculation of the total properties, such as the energy, energy gradient, dipole moment etc., an interaction energy is obtained for each pair of fragments. This pair interaction energy can be further decomposed into electrostatic, exchange, charge transfer and dispersion contributions. This analysis is known as the pair interaction energy decomposition analysis (PIEDA) and it can be thought of as FMO-based EDA. Alternatively, configuration analysis for fragment interaction (CAFI) and fragment interaction analysis based on local MP2 (FILM) were suggested within the FMO framework.
In FMO, various wave functions can be used for ab initio calculations of fragments and their dimers, such as Hartree–Fock, Density functional theory (DFT), Multi-configurational self-consistent field (MCSCF), time-dependent DFT (TDDFT), configuration interaction (CI), second order Møller–Plesset perturbation theory (MP2), coupled cluster (CC), and density functional tight binding (DFTB). Excited states can be computed with Configuration interaction (CI), Time-dependent density functional theory (TDDFT), time-dependent DFTB (TD-DFTB), equation-of-motion coupled-cluster (EOM-CC), and GW approximation. The solvent effects can be treated with the Polarizable continuum model (PCM). The FMO code is very efficiently parallelized utilising the generalized distributed data interface (GDDI) and hundreds of CPUs can be used with nearly perfect scaling.
In the FMO book published in 2009,[9] one can find 10 illustrated chapters written by the experts in the FMO development and applications, as well as a CDROM with annotated samples of input and output files, Facio modelling software and video tutorials (AppliGuide movies, showing mouse clicks) for treating difficult PDB files with Facio. In addition to this book, there are several chapters published in other books.[10][11][12][13]
In 2013–2014, a Japanese journal, CICSJ Bulletin, published a series of FMO papers in Japanese (about 100 pages in total), which give a representative summary of the recent FMO development and applications done in Japan, including papers on the GAMESS/FMO interface in Facio and developing an OpenMP version of GAMESS/FMO on the K computer.[14] Later, FMO was ported on the Fugaku [15] and Summit supercomputers.[16]
There are several reviews of FMO published in 2007,[17] 2012,[18] 2014,[19] 2017,[20] and 2022.[21] Two new FMO books were published in 2021[22] and 2023.[23]
The largest systems computed with FMO using DFTB are (a) a slab of fullerite surface containing 1,030,440 atoms, whose geometry was fully optimized and (b) a 10.7 μm white graphene nano material containing 1,180,800 atoms, for which Molecular dynamics simulations were performed.
To facilitate applications of FMO to drug discovery, FMO consortium was established.[24]
Applications of FMO
There are two main application fields of FMO: biochemistry and molecular dynamics of chemical reactions in solution. In addition, there is an emerging field of inorganic applications.
In 2005, an application of FMO to the calculation of the ground electronic state of photosynthetic protein with more than 20,000 atoms was distinguished with the best technical paper award at Supercomputing 2005. A number of applications of FMO to biochemical problems has been published, for instance, to Drug design, quantitative structure-activity relationship (QSAR) as well as the studies of excited states and chemical reactions of biological systems. The adaptive frozen orbital (AFO) treatment of the detached bonds was developed for FMO, making it possible to study solids, surfaces and nano systems, such as silicon nanowrires. FMO-TDDFT was applied to the excited states of molecular crystals (quinacridone).
Among inorganic systems, silica-related materials (zeolites, mesoporous nanoparticles and silica surfaces) were studied with FMO, as well as ionic liquids and boron nitride ribbons. There are other applications of FMO.[25]
Software for FMO
The FMO method is implemented in GAMESS (US), ABINIT-MP,[26] PAICS,[27] and OpenFMO [28] software packages, distributed free of charge.
Fu,[29] is a general open-source GUI that can generate input files for FMO. Another graphical user interface Facio[30] developed by M. Suenaga has a very convenient specialised support of FMO (in addition to other features), with which an automatic fragmentation of molecular clusters, proteins, nucleotides, saccharides and any combination thereof (e.g., DNA and protein complexes in explicit solvent) can be done in a few minutes, and a manual fragmentation of solids and surfaces can be accomplished by clicking the bonds to be detached. Facio can also visualise results of FMO calculations, such as the pair interactions.
FMO implementation in GAMESS
(E - energy, G - gradient, H - Hessian; bold - can be used with PCM)
|
See also
References
- K. Kitaura; E. Ikeo; T. Asada; T. Nakano; M. Uebayasi (1999). "Fragment molecular orbital method: an approximate computational method for large molecules". Chem. Phys. Lett. 313 (3–4): 701–706. Bibcode:1999CPL...313..701K. doi:10.1016/S0009-2614(99)00874-X.
- P. Otto; J. Ladik (1975). "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". Chem. Phys. 8 (1–2): 192–200. Bibcode:1975CP......8..192O. doi:10.1016/0301-0104(75)80107-8.
- Stoll, Hermann (1992). "Correlation energy of diamond". Physical Review B. 46 (11): 6700–6704. Bibcode:1992PhRvB..46.6700S. doi:10.1103/PhysRevB.46.6700. PMID 10002370.
- Huang, Lulu; Massa, Lou; Karle, Jerome (2005). "Kernel energy method illustrated with peptides". International Journal of Quantum Chemistry. 103 (6). Wiley: 808–817. Bibcode:2005IJQC..103..808H. doi:10.1002/qua.20542. ISSN 0020-7608.
- Dahlke, Erin E.; Truhlar, Donald G. (2006-11-04). "Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters". Journal of Chemical Theory and Computation. 3 (1). American Chemical Society (ACS): 46–53. CiteSeerX 10.1.1.186.8799. doi:10.1021/ct600253j. ISSN 1549-9618. PMID 26627150.
- Hirata, So; Valiev, Marat; Dupuis, Michel; Xantheas, Sotiris S.; Sugiki, Shinichiro; Sekino, Hideo (2005-08-10). "Fast electron correlation methods for molecular clusters in the ground and excited states". Molecular Physics. 103 (15–16). Informa UK Limited: 2255–2265. Bibcode:2005MolPh.103.2255H. doi:10.1080/00268970500083788. ISSN 0026-8976. S2CID 95428135.
- Kamiya, Muneaki; Hirata, So; Valiev, Marat (2008-02-21). "Fast electron correlation methods for molecular clusters without basis set superposition errors". The Journal of Chemical Physics. 128 (7). AIP Publishing: 074103. Bibcode:2008JChPh.128g4103K. doi:10.1063/1.2828517. ISSN 0021-9606. PMID 18298136.
- Gordon, Mark S.; Fedorov, Dmitri G.; Pruitt, Spencer R.; Slipchenko, Lyudmila V. (2011-08-26). "Fragmentation Methods: A Route to Accurate Calculations on Large Systems". Chemical Reviews. 112 (1). American Chemical Society (ACS): 632–672. doi:10.1021/cr200093j. ISSN 0009-2665. PMID 21866983.
- The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems, edited by D. G. Fedorov, K. Kitaura, CRC Press, Boca Raton, Florida, 2009 ISBN 978-1-4200-7848-0
- "(a) D. G. Fedorov, K. Kitaura, Theoretical development of the fragment molecular orbital (FMO) method and (b) T. Nakano, Y. Mochizuki, K. Fukuzawa, S. Amari, S. Tanaka, Developments and applications of ABINIT-MP software based on the fragment molecular orbital method in Modern methods for theoretical physical chemistry of biopolymers, edited by E. Starikov, J. Lewis, S. Tanaka, Elsevier, Amsterdam, 2006, ISBN 978-0-444-52220-7
- T. Nagata, D. G. Fedorov, K. Kitaura (2011). "Mathematical Formulation of the fragment molecular orbital method" in Linear-Scaling Techniques in Computational Chemistry and Physics. R. Zalesny, M. G. Papadopoulos, P. G. Mezey, J. Leszczyński (Eds.), Springer, New York, pp. 17-64.
- Y. Komeiji, Y. Mochizuki, T. Nakano, H. Mori (2012). "Recent advances in fragment molecular orbital-based molecular dynamics (FMO-MD) simulations", in Molecular Dynamics - Theoretical Developments and Applications in Nanotechnology and Energy, L. Wang (Ed.), Intech, pp. 3-24.
- Quantum mechanics in drug discovery., edited by A. Heifetz, Springer, New York, 2020 ISBN 978-1-0716-0281-2
- "日本化学会情報化学部会誌".
- "FMO@Fugaku". YouTube.
- "FMO@Summit". 18 November 2022. pp. 1–14.
- D. G. Fedorov; et al. (2007). "Extending the Power of Quantum Chemistry to Large Systems with the Fragment Molecular Orbital Method". J. Phys. Chem. A. 111 (30): 6904–6914. Bibcode:2007JPCA..111.6904F. doi:10.1021/jp0716740. PMID 17511437.
- Fedorov, Dmitri G.; Nagata, Takeshi; Kitaura, Kazuo (2012). "Exploring chemistry with the fragment molecular orbital method". Physical Chemistry Chemical Physics. 14 (21). Royal Society of Chemistry (RSC): 7562–7577. Bibcode:2012PCCP...14.7562F. doi:10.1039/c2cp23784a. ISSN 1463-9076. PMID 22410762.
- Tanaka, Shigenori; Mochizuki, Yuji; Komeiji, Yuto; Okiyama, Yoshio; Fukuzawa, Kaori (2014). "Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems". Phys. Chem. Chem. Phys. 16 (22). Royal Society of Chemistry (RSC): 10310–10344. Bibcode:2014PCCP...1610310T. doi:10.1039/c4cp00316k. ISSN 1463-9076. PMID 24740821.
- Fedorov, Dmitri G. (2017). "The fragment molecular orbital method: theoretical development, implementation in GAMESS, and applications". WIREs Comput. Mol. Sci. 7 (6): e1322. doi:10.1002/wcms.1322. S2CID 103065280.
- Fukuzawa, Kaori; Tanaka, Shigenori (2022). "Fragment molecular orbital calculations for biomolecules". Current Opinion Structural Biology. 72: 127–134. doi:10.1016/j.sbi.2021.08.010. PMID 34656048. S2CID 239010379.
- Recent advances of the fragment molecular orbital method, edited by Y. Mochizuki, S. Tanaka, K. Fukuzawa, Springer, Singapore, 2021 ISBN 978-981-15-9234-8
- Complete guide to the fragment molecular orbital method in GAMESS., D. G. Fedorov, World Scientific, Singapore, 2023 ISBN 978-981-126-362-0
- "FMO consortium".
- http://staff.aist.go.jp/d.g.fedorov/fmo/fmoref.txt
- "ABINIT-MP (in Japanese)".
- "PAICS".
- "OpenFMO". GitHub.
- "fu". GitHub.
- "Facio the 3D-Molecular Modeling Software".