1
|
Fedorov DG. Use of caps in the auxiliary basis set formulation of the fragment molecular orbital method. J Comput Chem 2024; 45:1540-1551. [PMID: 38490813 DOI: 10.1002/jcc.27345] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 02/19/2024] [Accepted: 02/29/2024] [Indexed: 03/17/2024]
Abstract
An auxiliary polarization formulation of the fragment molecular orbital (FMO) method is developed, combining a basis set correction computed for capped isolated fragments with a polarization obtained from uncapped fragments. For a set of organic and inorganic test systems, it is shown that the total energy and atomic charges are accurately reproduced with respect to full unfragmented calculations. It is demonstrated that the method is accurate for computing electronic excited states. The developed approach is applied to rank the isomers of chignolin from experimental NMR data (PDB: 1UAO) according to their relative energy. Contributions of polarization and basis set effects to pair interactions between fragments are elucidated.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan
| |
Collapse
|
2
|
Fedorov DG. Site-Specific Ionization Potentials and Electron Affinities in Large Molecular Systems at Coupled Cluster Level. J Phys Chem A 2023; 127:9357-9364. [PMID: 37782030 DOI: 10.1021/acs.jpca.3c04847] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
A many-body expansion of ionization potentials and electron affinities is developed based on a combination of the fragment molecular orbital method and equation-of-motion coupled-cluster (EOM-CC). In addition to site-specific values, obtained as one-body properties, pair and triple corrections are added to account for nonlocal EOM-CC contributions of the molecular environment of a chromophore. The developed method is applied to carboxylic acids, alkyl cations, a protein ubiquitin (Protein Data Bank ID 1UBQ), and a nano ribbon of white graphene elucidating the effect of environment on ionization potential and electron affinity.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
3
|
Nakata H, Fedorov DG. Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2023; 19:1276-1285. [PMID: 36753486 DOI: 10.1021/acs.jctc.2c01177] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).
Collapse
Affiliation(s)
- Hiroya Nakata
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
4
|
Fujita T, Noguchi Y. Fragment-Based Excited-State Calculations Using the GW Approximation and the Bethe-Salpeter Equation. J Phys Chem A 2021; 125:10580-10592. [PMID: 34871000 DOI: 10.1021/acs.jpca.1c07337] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Herein, we present a fragment-based approach for calculating the charged and neutral excited states in molecular systems, based on the many-body Green's function method within the GW approximation and the Bethe-Salpeter equation (BSE). The implementation relies on the many-body expansion of the total irreducible polarizability on the basis of fragment molecular orbitals. The GW quasi-particle energies in complex molecular environments are obtained by the GW calculation for the target fragment plus induced polarization contributions of the surrounding fragments at the static Coulomb-hole plus screened exchange level. In addition, we develop a large-scale GW/BSE method for calculating the delocalized excited states of molecular aggregates, based on the fragment molecular orbital method and the exciton model. The accuracy of fragment-based GW and GW/BSE methods was evaluated on molecular clusters and molecular crystals. We found that the accuracy of the total irreducible polarizability can be improved systematically by including two-body correction terms, and the fragment-based calculations can reasonably reproduce the results of the corresponding unfragmented calculations with a relative error of less than 100 meV. The proposed approach enables efficient excited-state calculations for large molecular systems with reasonable accuracy.
Collapse
Affiliation(s)
- Takatoshi Fujita
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, Tokai, Ibaraki 319-1106, Japan
| | - Yoshifumi Noguchi
- Department of Applied Chemistry and Biochemical Engineering, Graduate School of Engineering, Shizuoka University, Hamamatsu, Shizuoka 432-8561, Japan
| |
Collapse
|
5
|
Fujita T, Noguchi Y, Hoshi T. Revisiting the Charge-Transfer States at Pentacene/C 60 Interfaces with the GW/Bethe-Salpeter Equation Approach. MATERIALS (BASEL, SWITZERLAND) 2020; 13:E2728. [PMID: 32560127 PMCID: PMC7345661 DOI: 10.3390/ma13122728] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Revised: 06/11/2020] [Accepted: 06/13/2020] [Indexed: 11/16/2022]
Abstract
Molecular orientations and interfacial morphologies have critical effects on the electronic states of donor/acceptor interfaces and thus on the performance of organic photovoltaic devices. In this study, we explore the energy levels and charge-transfer states at the organic donor/acceptor interfaces on the basis of the fragment-based GW and Bethe-Salpeter equation approach. The face-on and edge-on orientations of pentacene/C60 bilayer heterojunctions have employed as model systems. GW+Bethe-Salpeter equation calculations were performed for the local interface structures in the face-on and edge-on bilayer heterojunctions, which contain approximately 2000 atoms. Calculated energy levels and charge-transfer state absorption spectra are in reasonable agreements with those obtained from experimental measurements. We found that the dependence of the energy levels on interfacial morphology is predominantly determined by the electrostatic contribution of polarization energy, while the effects of induction contribution in the edge-on interface are similar to those in the face-on. Moreover, the delocalized charge-transfer states contribute to the main absorption peak in the edge-on interface, while the face-on interface features relatively localized charge-transfer states in the main absorption peak. The impact of the interfacial morphologies on the polarization and charge delocalization effects is analyzed in detail.
Collapse
Affiliation(s)
| | - Yoshifumi Noguchi
- Department of Applied Chemistry and Biochemical Engineering, Graduate School of Engineering, Shizuoka University, Hamamatsu, Shizuoka 432-8561, Japan;
| | - Takeo Hoshi
- Department of Applied Mathematics and Physics, Tottori University, Tottori 680-8550, Japan;
| |
Collapse
|
6
|
Fujita T, Mochizuki Y. Development of the Fragment Molecular Orbital Method for Calculating Nonlocal Excitations in Large Molecular Systems. J Phys Chem A 2018; 122:3886-3898. [DOI: 10.1021/acs.jpca.8b00446] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute for Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
7
|
Fedorov DG. The fragment molecular orbital method: theoretical development, implementation in
GAMESS
, and applications. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1322] [Citation(s) in RCA: 71] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD‐FMat)National Institute of Advanced Industrial Science and Technology (AIST)TsukubaJapan
| |
Collapse
|
8
|
Tanaka S, Mochizuki Y, Komeiji Y, Okiyama Y, Fukuzawa K. Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems. Phys Chem Chem Phys 2015; 16:10310-44. [PMID: 24740821 DOI: 10.1039/c4cp00316k] [Citation(s) in RCA: 189] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.
Collapse
Affiliation(s)
- Shigenori Tanaka
- Graduate School of System Informatics, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan.
| | | | | | | | | |
Collapse
|
9
|
Nakata H, Fedorov DG, Yokojima S, Kitaura K, Sakurai M, Nakamura S. Unrestricted density functional theory based on the fragment molecular orbital method for the ground and excited state calculations of large systems. J Chem Phys 2014; 140:144101. [DOI: 10.1063/1.4870261] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|
10
|
Tsuchiya T, Shrestha K, Jakubikova E. Orbital Analysis and Excited-State Calculations in an Energy-Based Fragmentation Method. J Chem Theory Comput 2013; 9:3350-63. [DOI: 10.1021/ct400025a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Takashi Tsuchiya
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United
States
| | - Kushal Shrestha
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United
States
| | - Elena Jakubikova
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United
States
| |
Collapse
|
11
|
Ma Y, Ma H. Calculating Excited States of Molecular Aggregates by the Renormalized Excitonic Method. J Phys Chem A 2013; 117:3655-65. [DOI: 10.1021/jp401168s] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yingjin Ma
- Key Laboratory of Mesoscopic Chemistry of MOE, School
of Chemistry and Chemical Engineering, Institute of Theoretical and
Computational Chemistry, Nanjing University, Nanjing 210093, China
| | - Haibo Ma
- Key Laboratory of Mesoscopic Chemistry of MOE, School
of Chemistry and Chemical Engineering, Institute of Theoretical and
Computational Chemistry, Nanjing University, Nanjing 210093, China
| |
Collapse
|
12
|
Ma Y, Liu Y, Ma H. A new fragment-based approach for calculating electronic excitation energies of large systems. J Chem Phys 2012; 136:024113. [DOI: 10.1063/1.3675915] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Fedorov DG, Nagata T, Kitaura K. Exploring chemistry with the fragment molecular orbital method. Phys Chem Chem Phys 2012; 14:7562-77. [DOI: 10.1039/c2cp23784a] [Citation(s) in RCA: 290] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
14
|
Wu F, Liu W, Zhang Y, Li Z. Linear-Scaling Time-Dependent Density Functional Theory Based on the Idea of “From Fragments to Molecule”. J Chem Theory Comput 2011; 7:3643-60. [DOI: 10.1021/ct200225v] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Fangqin Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Wenjian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Yong Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Zhendong Li
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People’s Republic of China
| |
Collapse
|
15
|
Gordon MS, Fedorov DG, Pruitt SR, Slipchenko LV. Fragmentation Methods: A Route to Accurate Calculations on Large Systems. Chem Rev 2011; 112:632-72. [DOI: 10.1021/cr200093j] [Citation(s) in RCA: 836] [Impact Index Per Article: 64.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Spencer R. Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
16
|
Taguchi N, Mochizuki Y, Nakano T. Fragment molecular orbital calculations for excitation energies of blue- and yellow-fluorescent proteins. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.01.054] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|