1
|
Nishimura R, Yoshikawa T, Sakata K, Nakai H. Excitation configuration analysis for divide-and-conquer excited-state calculation method using dynamical polarizability. J Chem Phys 2024; 160:244103. [PMID: 38913842 DOI: 10.1063/5.0207935] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 05/23/2024] [Indexed: 06/26/2024] Open
Abstract
The authors previously developed a divide-and-conquer (DC)-based non-local excited-state calculation method for large systems using dynamical polarizability [Nakai and Yoshikawa, J. Chem. Phys. 146, 124123 (2017)]. This method evaluates the excitation energies and oscillator strengths using information on the dynamical polarizability poles. This article proposes a novel analysis of the previously developed method to obtain further configuration information on excited states, including excitation and de-excitation coefficients of each excitation configuration. Numerical applications to simple molecules, such as ethylene, hydrogen molecule, ammonia, and pyridazine, confirmed that the proposed analysis could accurately reproduce the excitation and de-excitation coefficients. The combination with the DC scheme enables both the local and non-local excited states of large systems with an excited nature to be treated.
Collapse
Affiliation(s)
- Ryusei Nishimura
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Ken Sakata
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
2
|
Fujita T, Hoshi T. Ab Initio Study of Charge Separation Dynamics and Pump-Probe Spectroscopy in the P3HT/PCBM Blend. J Phys Chem B 2023; 127:7615-7623. [PMID: 37639551 DOI: 10.1021/acs.jpcb.3c02458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
We develop a bottom-up computational method for excited-state dynamics and time-resolved spectroscopy signals in molecular aggregates, on the basis of ab initio excited-state calculations. As an application, we consider the charge separation dynamics and pump-probe spectroscopy in the amorphous P3HT/PCBM blend. To simulate quantum dynamics and time-resolved spectroscopy, the model Hamiltonian for single-excitation and double-excitation manifolds was derived on the basis of fragment-based excited-state calculations within the GW approximation and the Bethe-Salpeter equation. After elucidating the energetics of the electron-hole separation and examining linear absorption spectrum, we investigated the quantum dynamics of exciton and charge carriers in comparison with the pump-probe transient absorption spectra. In particular, we introduced the pump-probe excited-state absorption (ESA) anisotropy as a spectroscopic signature of charge carrier dynamics after exciton dissociation. We found that the charge separation dynamics can be probed by the pump-probe ESA anisotropy dynamics after charge-transfer excitations. The present study provides the fundamental information for understanding the experimental spectroscopy signals, by elucidating the relationship between the excited states, the exciton and charge carrier dynamics, and time-resolved spectroscopy.
Collapse
Affiliation(s)
- Takatoshi Fujita
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, Chiba 263-8555, Japan
| | - Takeo Hoshi
- Department of Mechanical and Physical Engineering, Faculty of Engineering, Tottori University, Tottori-shi 680-8552, Tottori, 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
|
Einsele R, Hoche J, Mitrić R. Long-range corrected fragment molecular orbital density functional tight-binding method for excited states in large molecular systems. J Chem Phys 2023; 158:044121. [PMID: 36725509 DOI: 10.1063/5.0136844] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Herein, we present a new method to efficiently calculate electronically excited states in large molecular assemblies, consisting of hundreds of molecules. For this purpose, we combine the long-range corrected tight-binding density functional fragment molecular orbital method (FMO-LC-DFTB) with an excitonic Hamiltonian, which is constructed in the basis of locally excited and charge-transfer configuration state functions calculated for embedded monomers and dimers and accounts explicitly for the electronic coupling between all types of excitons. We first evaluate both the accuracy and efficiency of our fragmentation approach for molecular dimers and aggregates by comparing it with the full LC-TD-DFTB method. The comparison of the calculated spectra of an anthracene cluster shows a very good agreement between our method and the LC-TD-DFTB reference. The effective computational scaling of our method has been explored for anthracene clusters and for perylene bisimide aggregates. We demonstrate the applicability of our method by the calculation of the excited state properties of pentacene crystal models consisting of up to 319 molecules. Furthermore, the participation ratio of the monomer fragments to the excited states is analyzed by the calculation of natural transition orbital participation numbers, which are verified by the hole and particle density for a chosen pentacene cluster. The use of our FMO-LC-TDDFTB method will allow for future studies of excitonic dynamics and charge transport to be performed on complex molecular systems consisting of thousands of atoms.
Collapse
Affiliation(s)
- Richard Einsele
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Strasse 42, 97074 Würzburg, Germany
| | - Joscha Hoche
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Strasse 42, 97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Strasse 42, 97074 Würzburg, Germany
| |
Collapse
|
5
|
Nakai H, Kobayashi M, Yoshikawa T, Seino J, Ikabata Y, Nishimura Y. Divide-and-Conquer Linear-Scaling Quantum Chemical Computations. J Phys Chem A 2023; 127:589-618. [PMID: 36630608 DOI: 10.1021/acs.jpca.2c06965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.
Collapse
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo, Hokkaido060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo, Hokkaido001-0021, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba274-8510, Japan
| | - Junji Seino
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Yasuhiro Ikabata
- Information and Media Center, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan.,Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| |
Collapse
|
6
|
Chakraborty A, Tribedi S, Maitra R. A double exponential coupled cluster theory in the fragment molecular orbital framework. J Chem Phys 2022; 156:244117. [DOI: 10.1063/5.0090115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Fragmentation-based methods enable electronic structure calculations for large chemical systems through partitioning them into smaller fragments. Here, we have developed and benchmarked a dual exponential operator-based coupled cluster theory to account for high-rank electronic correlation of large chemical systems within the fragment molecular orbital (FMO) framework. Upon partitioning the molecular system into several fragments, the zeroth order reference determinants for each fragment and fragment pair are constructed in a self-consistent manner with two-body FMO expansion. The dynamical correlation is induced through a dual exponential ansatz with a set of fragment-specific rank-one and rank-two operators that act on the individual reference determinants. While the single and double excitations for each fragment are included through the conventional rank-one and rank-two cluster operators, the triple excitation space is spanned via the contraction between the cluster operators and a set of rank-two scattering operators over a few optimized fragment-specific occupied and virtual orbitals. Thus, the high-rank dynamical correlation effects within the FMO framework are computed with rank-one and rank-two parametrization of the wave operator, leading to significant reduction in the number of variables and associated computational scaling over the conventional methods. Through a series of pilot numerical applications on various covalent and non-covalently bonded systems, we have shown the quantitative accuracy of the proposed methodology compared to canonical, as well as FMO-based coupled-cluster single double triple. The accuracy of the proposed method is shown to be systematically improvable upon increasing the number of contractible occupied and virtual molecular orbitals employed to simulate triple excitations.
Collapse
Affiliation(s)
- Anish Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Soumi Tribedi
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
7
|
Orimoto Y, Hisama K, Aoki Y. Local electronic structure analysis by ab initio elongation method: A benchmark using DNA block polymers. J Chem Phys 2022; 156:204114. [DOI: 10.1063/5.0087726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ab initio elongation (ELG) method based on a polymerization concept is a feasible way to perform linear-scaling electronic structure calculations for huge aperiodic molecules while maintaining computational accuracy. In the method, the electronic structures are sequentially elongated by repeating (1) the conversion of canonical molecular orbitals (CMOs) to region-localized MOs (RLMOs), that is, active RLMOs localized onto a region close to an attacking monomer or frozen RLMOs localized onto the remaining region, and the subsequent (2) partial self-consistent-field calculations for an interaction space composed of the active RLMOs and the attacking monomer. For each ELG process, one can obtain local CMOs for the interaction space and the corresponding local orbital energies. Local site information, such as the local highest-occupied/lowest-unoccupied MOs, can be acquired with linear-scaling efficiency by correctly including electronic effects from the frozen region. In this study, we performed a local electronic structure analysis using the ELG method for various DNA block polymers with different sequential patterns. This benchmark aimed to confirm the effectiveness of the method toward the efficient detection of a singular local electronic structure in unknown systems as a future practical application. We discussed the high-throughput efficiency of our method and proposed a strategy to detect singular electronic structures by combining with a machine learning technique.
Collapse
Affiliation(s)
- Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Keisuke Hisama
- Department of Interdisciplinary Engineering Sciences, Chemistry and Materials Science, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| |
Collapse
|
8
|
Friedl C, Fedorov DG, Renger T. Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches. Phys Chem Chem Phys 2022; 24:5014-5038. [PMID: 35142765 PMCID: PMC8865841 DOI: 10.1039/d1cp03566e] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0–0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang–Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang–Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0–0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling. The importance of reaction field and screening effects on the excitonic couplings is demonstrated, and from quantum-chemical calculations a single scaling factor is derived that can be used to improve simple models based on the Poisson equation.![]()
Collapse
Affiliation(s)
- Christian Friedl
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria.
| | - 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.
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria.
| |
Collapse
|
9
|
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
|
10
|
Uratani H, Yoshikawa T, Nakai H. Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding. J Chem Theory Comput 2021; 17:1290-1300. [PMID: 33577323 DOI: 10.1021/acs.jctc.0c01155] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
11
|
Paz ASP, Glover WJ. Diabatic Many-Body Expansion: Development and Application to Charge-Transfer Reactions. J Chem Theory Comput 2021; 17:1497-1511. [DOI: 10.1021/acs.jctc.0c01231] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amiel S. P. Paz
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - William J. Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
12
|
MOCHIZUKI Y, NAKANO T, SATO S, SAKAKURA K, WATANABE H, OKUWAKI K, OHSHIMA S, KATAGIRI T. Development Status of ABINIT-MP in 2021. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2021. [DOI: 10.2477/jccj.2022-0001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Yuji MOCHIZUKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
| | - Tatsuya NAKANO
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| | - Shinya SATO
- NEC Solution Innovators, Ltd. (Shiromi 1-4-24, NEC Kansai Building, Chuo-ku, Osaka, 540-8551, Japan)
| | - Kota SAKAKURA
- Foundation for Computational Science (Minatojima Minamicho 7-1-28 Computational Science Center Building 1F, Chuo-ku, Kobe, 650-0047, Japan)
| | - Hiromasa WATANABE
- HPC Systems Inc. (Kaigan 3-9-15 LOOP-X 8F, Minato-ku, Tokyo, 108-0022, Japan)
| | - Koji OKUWAKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
| | - Satoshi OHSHIMA
- Information Technology Center, Nagoya University (Furo-cho, Chikusa-ku, Nagoya, 464-8601, Japan)
| | - Takahiro KATAGIRI
- Information Technology Center, Nagoya University (Furo-cho, Chikusa-ku, Nagoya, 464-8601, Japan)
| |
Collapse
|
13
|
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
|
14
|
Yoshikawa T, Yoshihara J, Nakai H. Large-scale excited-state calculation using dynamical polarizability evaluated by divide-and-conquer based coupled cluster linear response method. J Chem Phys 2020; 152:024102. [PMID: 31941302 DOI: 10.1063/1.5124909] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This study attempted to propose an efficient scheme at the coupled cluster linear response (CCLR) level to perform large-scale excited-state calculations of not only local excitations but also nonlocal ones such as charge transfers and transitions between delocalized orbitals. Although standard applications of fragmentation techniques to the excited-state calculations brought about the limitations that could only deal with local excitations, this study solved the problem by evaluating the excited states as the poles of dynamical polarizability. Because such an approach previously succeeded at the time-dependent density functional theory level [H. Nakai and T. Yoshikawa, J. Chem. Phys. 146, 124123 (2017)], this study was considered as an extension to the CCLR level. To evaluate the dynamical polarizability at the CCLR level, we revisited three equivalent formulas, namely, coupled-perturbed self-consistent field (CPSCF), random phase approximation (RPA), and Green's function (GF). We further extended these formulas to the linear-scaling methods based on the divide-and-conquer (DC) technique. We implemented the CCLR with singles and doubles (CCSDLR) program for the six schemes, i.e., the standard and DC-type CPSCF, RPA, and GF. Illustrative applications of the present methods demonstrated the accuracy and efficiency. Although the standard three treatments could exactly reproduced the conventional frequency-domain CCSDLR results, their computational costs were commonly higher than that of the conventional ones due to large amount of computations for individual frequencies of the external electric field. The DC-type treatments, which approximately reproduced the conventional results, could achieve quasilinear scaling computational costs. Among them, DC-GF was found to exhibit the best performance.
Collapse
Affiliation(s)
- Takeshi Yoshikawa
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Jyunya Yoshihara
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| |
Collapse
|
15
|
MOCHIZUKI Y, SAKAKURA K, WATANABE H, OKUWAKI K, KATO K, WATANABE N, OKIYAMA Y, FUKUZAWA K, NAKANO T. Development Status of ABINIT-MP in 2020. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2020. [DOI: 10.2477/jccj.2021-0015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Yuji MOCHIZUKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
| | - Kota SAKAKURA
- Foundation for Computational Science (Minatojima Minamicho 7-1-28 Computational Science Center Building 1F, Chuo-ku, Kobe, 650-0047, Japan)
| | - Hiromasa WATANABE
- HPC Systems Inc. (Kaigan 3-9-15 LOOP-X 8F, Minato-ku, Tokyo, 108-0022, Japan)
| | - Koji OKUWAKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
| | - Koichiro KATO
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University (744 Moto-oka, Nishi-ku, Fukuoka, 819-0395, Japan)
| | - Naoki WATANABE
- Mizuho Information & Research Institute, Inc. (Kanda-Nishikicho 2-3, Chiyoda-ku, Tokyo, 101-8443, Japan)
| | - Yoshio OKIYAMA
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| | - Kaori FUKUZAWA
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University (Ebara 2-4-41, Shinagawa-ku, Tokyo, 142-8501, Japan)
| | - Tatsuya NAKANO
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| |
Collapse
|
16
|
Yoshikawa T, Komoto N, Nishimura Y, Nakai H. GPU-Accelerated Large-Scale Excited-State Simulation Based on Divide-and-Conquer Time-Dependent Density-Functional Tight-Binding. J Comput Chem 2019; 40:2778-2786. [PMID: 31441083 DOI: 10.1002/jcc.26053] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Revised: 08/04/2019] [Accepted: 08/07/2019] [Indexed: 01/09/2023]
Abstract
The present study implemented the divide-and-conquer time-dependent density-functional tight-binding (DC-TDDFTB) code on a graphical processing unit (GPU). The DC method, which is a linear-scaling scheme, divides a total system into several fragments. By separately solving local equations in individual fragments, the DC method could reduce slow central processing unit (CPU)-GPU memory access, as well as computational cost, and avoid shortfalls of GPU memory. Numerical applications confirmed that the present code on GPU significantly accelerated the TDDFTB calculations, while maintaining accuracy. Furthermore, the DC-TDDFTB simulation of 2-acetylindan-1,3-dione displays excited-state intramolecular proton transfer and provides reasonable absorption and fluorescence energies with the corresponding experimental values. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Takeshi Yoshikawa
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan
| | - Nana Komoto
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan.,Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto, 615-8520, Japan
| |
Collapse
|
17
|
Fujita T, Noguchi Y, Hoshi T. Charge-transfer excited states in the donor/acceptor interface from large-scale GW calculations. J Chem Phys 2019; 151:114109. [PMID: 31542033 DOI: 10.1063/1.5113944] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Predicting the charge-transfer (CT) excited states across the donor/acceptor (D/A) interface is essential for understanding the charge photogeneration process in an organic solar cell. Here, we present a fragment-based GW implementation that can be applied to a D/A interface structure and thus enables accurate determination of the CT states. The implementation is based on the fragmentation approximation of the polarization function and the combined GW and Coulomb-hole plus screened exchange approximations for self-energies. The fragment-based GW is demonstrated by application to the pentacene/C60 interface structure containing more than 2000 atoms. The CT excitation energies were estimated from the quasiparticle energies and electron-hole screened Coulomb interactions; the computed energies are in reasonable agreement with experimental estimates from the external quantum efficiency measurements. We highlight the impact of the induced polarization effects on the electron-hole energetics. The proposed fragment-based GW method offers a first-principles tool to compute the quasiparticle energies and electronic excitation energies of organic materials.
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
|
18
|
Mester D, Kállay M. Reduced-Scaling Approach for Configuration Interaction Singles and Time-Dependent Density Functional Theory Calculations Using Hybrid Functionals. J Chem Theory Comput 2019; 15:1690-1704. [DOI: 10.1021/acs.jctc.8b01199] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box
91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box
91, H-1521 Budapest, Hungary
| |
Collapse
|
19
|
MOCHIZUKI Y, AKINAGA Y, SAKAKURA K, WATANABE H, KATO K, WATANABE N, OKUWAKI K, NAKANO T, FUKUZAWA K. Current Status of ABINIT-MP Open Series. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2019. [DOI: 10.2477/jccj.2019-0016] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Yuji MOCHIZUKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
| | - Yoshinobu AKINAGA
- VINAS Co.,Ltd.(Dojima 2-1-31, Keihan-Dojima Build. 8F, Kita-ku, Osaka, 530-003, Japan)
| | - Kota SAKAKURA
- NEC Corporation (Shiba 5-7-1, Minato-ku, Tokyo, 108-8001, Japan)
| | - Hiromasa WATANABE
- HPC Systems Inc. (Kaigan 3-9-15 LOOP-X 8F, Minato-ku, Tokyo, 108-0022, Japan)
| | - Kouichiro KATO
- Mizuho Information & Research Institute, Inc. (Kanda-Nishikicho 2-3, Chiyoda-ku, Tokyo, 101-8443, Japan)
| | - Naoki WATANABE
- Mizuho Information & Research Institute, Inc. (Kanda-Nishikicho 2-3, Chiyoda-ku, Tokyo, 101-8443, Japan)
| | - Koji OKUWAKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
| | - Tatsuya NAKANO
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| | - Kaori FUKUZAWA
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University (Ebara 2-4-41, Shinagawa-ku, Tokyo, 142-8501, Japan)
| |
Collapse
|
20
|
Fujita T, Alam MK, Hoshi T. Thousand-atom ab initio calculations of excited states at organic/organic interfaces: toward first-principles investigations of charge photogeneration. Phys Chem Chem Phys 2018; 20:26443-26452. [PMID: 30306163 DOI: 10.1039/c8cp05574b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Predicting electronically excited states across electron-donor/electron acceptor interfaces is essential for understanding the charge photogeneration process in organic solar cells. However, organic solar cells are large and disordered systems, and their excited states cannot be easily accessed by conventional quantum chemistry approaches. Moreover, a large number of excited states must be obtained to fully understand the charge separation mechanism. Recently, we have developed a novel fragment-based excited state method which can efficiently calculate a large number of states in molecular aggregates. In this article, we demonstrate the large-scale excited-state calculations by investigating interfacial charge transfer (ICT) states across the electron-donor/electron acceptor interfaces. As the model systems, we considered the face-on and edge-on configurations of pentacene/C60 bilayer heterojunction structures. These model structures contain approximately 1.8 × 105 atoms, and their local interface regions containing 2000 atoms were treated quantum mechanically, embedded in the electrostatic potentials from the remaining parts. Therefore, the charge delocalization effect, structural disorder, and the resulting heterogeneous electrostatic and polarizable environments were taken into account in the excited-state calculations. The computed energies of the low-lying ICT states are in reasonable agreement with experimental estimates. By comparing the edge-on and face-on configurations of the pentacene/C60 interfaces, we discuss the influence of interfacial morphologies on the energetics and charge delocalization of ICT states. In addition, we present the detailed characterization of excited states and highlight the importance of hybridization effects between pentacene excited states and ICT states. The large-scale ab initio calculations for the interface systems enabled the exploration of the ICT states, leading to first-principles investigation of the charge separation mechanism in organic solar cells.
Collapse
Affiliation(s)
- Takatoshi Fujita
- Institute for Molecular Science, Okazaki, Aichi 444-0865, Japan.
| | - Md Khorshed Alam
- Department of Physics, University of Barisal, Barisal-8200, Bangladesh
| | - Takeo Hoshi
- Department of Applied Mathmatics and Physics, Tottori University, Tottori 680-8550, Japan
| |
Collapse
|
21
|
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
|
22
|
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
|
23
|
Ge Q, Mao Y, White AF, Epifanovsky E, Closser KD, Head-Gordon M. Simulating the absorption spectra of helium clusters (N = 70, 150, 231, 300) using a charge transfer correction to superposition of fragment single excitations. J Chem Phys 2017; 146:044111. [DOI: 10.1063/1.4973611] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qinghui Ge
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alec F. White
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Evgeny Epifanovsky
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Kristina D. Closser
- Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
24
|
Fujita T, Atahan-Evrenk S, Sawaya NPD, Aspuru-Guzik A. Coherent Dynamics of Mixed Frenkel and Charge-Transfer Excitons in Dinaphtho[2,3-b:2'3'-f]thieno[3,2-b]-thiophene Thin Films: The Importance of Hole Delocalization. J Phys Chem Lett 2016; 7:1374-1380. [PMID: 27011327 DOI: 10.1021/acs.jpclett.6b00364] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Charge-transfer states in organic semiconductors play crucial roles in processes such as singlet fission and exciton dissociation at donor/acceptor interfaces. Recently, a time-resolved spectroscopy study of dinaphtho[2,3-b:2'3'-f]thieno[3,2-b]-thiophene (DNTT) thin films provided evidence for the formation of mixed Frenkel and charge-transfer excitons after the photoexcitation. Here, we investigate optical properties and excitation dynamics of the DNTT thin films by combining ab initio calculations and a stochastic Schrödinger equation. Our theory predicts that the low-energy Frenkel exciton band consists of 8-47% CT character. The quantum dynamics simulations show coherent dynamics of Frenkel and CT states in 50 fs after the optical excitation. We demonstrate the role of charge delocalization and localization in the mixing of CT states with Frenkel excitons as well as the role of their decoherence.
Collapse
Affiliation(s)
- Takatoshi Fujita
- Department of Chemistry, Graduate School of Science, Kyoto University , Kyoto 606-8502, Japan
| | - Sule Atahan-Evrenk
- TOBB University of Economics and Technology , Sogutozu, Ankara 06560, Turkey
| | - Nicolas P D Sawaya
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States
| |
Collapse
|
25
|
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
|
26
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
27
|
SAKAGUCHI M, MOCHIZUKI Y, WATANABE C, FUKUZAWA K. Effects of Water Molecules and Configurations of Neighboring Amino Acid Residues Surrounding DsRed Chromophore on Its Excitation Energy. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2015. [DOI: 10.2477/jccj.2015-0033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Masataka SAKAGUCHI
- Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yuji MOCHIZUKI
- Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Chiduru WATANABE
- Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Kaori FUKUZAWA
- Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
- School of Dentistry at Matsudo, Nihon University, 2-870-1 Sakaemachi-Nishi, Matsudo 271-8571, Japan
| |
Collapse
|
28
|
A Mini-review on Chemoinformatics Approaches for Drug Discovery. JOURNAL OF COMPUTER AIDED CHEMISTRY 2015. [DOI: 10.2751/jcac.16.15] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
29
|
Pruitt SR, Bertoni C, Brorsen KR, Gordon MS. Efficient and accurate fragmentation methods. Acc Chem Res 2014; 47:2786-94. [PMID: 24810424 DOI: 10.1021/ar500097m] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Conspectus Three novel fragmentation methods that are available in the electronic structure program GAMESS (general atomic and molecular electronic structure system) are discussed in this Account. The fragment molecular orbital (FMO) method can be combined with any electronic structure method to perform accurate calculations on large molecular species with no reliance on capping atoms or empirical parameters. The FMO method is highly scalable and can take advantage of massively parallel computer systems. For example, the method has been shown to scale nearly linearly on up to 131 000 processor cores for calculations on large water clusters. There have been many applications of the FMO method to large molecular clusters, to biomolecules (e.g., proteins), and to materials that are used as heterogeneous catalysts. The effective fragment potential (EFP) method is a model potential approach that is fully derived from first principles and has no empirically fitted parameters. Consequently, an EFP can be generated for any molecule by a simple preparatory GAMESS calculation. The EFP method provides accurate descriptions of all types of intermolecular interactions, including Coulombic interactions, polarization/induction, exchange repulsion, dispersion, and charge transfer. The EFP method has been applied successfully to the study of liquid water, π-stacking in substituted benzenes and in DNA base pairs, solvent effects on positive and negative ions, electronic spectra and dynamics, non-adiabatic phenomena in electronic excited states, and nonlinear excited state properties. The effective fragment molecular orbital (EFMO) method is a merger of the FMO and EFP methods, in which interfragment interactions are described by the EFP potential, rather than the less accurate electrostatic potential. The use of EFP in this manner facilitates the use of a smaller value for the distance cut-off (Rcut). Rcut determines the distance at which EFP interactions replace fully quantum mechanical calculations on fragment-fragment (dimer) interactions. The EFMO method is both more accurate and more computationally efficient than the most commonly used FMO implementation (FMO2), in which all dimers are explicitly included in the calculation. While the FMO2 method itself does not incorporate three-body interactions, such interactions are included in the EFMO method via the EFP self-consistent induction term. Several applications (ranging from clusters to proteins) of the three methods are discussed to demonstrate their efficacy. The EFMO method will be especially exciting once the analytic gradients have been completed, because this will allow geometry optimizations, the prediction of vibrational spectra, reaction path following, and molecular dynamics simulations using the method.
Collapse
Affiliation(s)
- Spencer R. Pruitt
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
- Argonne Leadership Computing Facility, Argonne, Illinois 60439, United States
| | - Colleen Bertoni
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Kurt R. Brorsen
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Mark S. Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
30
|
Fedorov DG, Asada N, Nakanishi I, Kitaura K. The use of many-body expansions and geometry optimizations in fragment-based methods. Acc Chem Res 2014; 47:2846-56. [PMID: 25144610 DOI: 10.1021/ar500224r] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Conspectus Chemists routinely work with complex molecular systems: solutions, biochemical molecules, and amorphous and composite materials provide some typical examples. The questions one often asks are what are the driving forces for a chemical phenomenon? How reasonable are our views of chemical systems in terms of subunits, such as functional groups and individual molecules? How can one quantify the difference in physicochemical properties of functional units found in a different chemical environment? Are various effects on functional units in molecular systems additive? Can they be represented by pairwise potentials? Are there effects that cannot be represented in a simple picture of pairwise interactions? How can we obtain quantitative values for these effects? Many of these questions can be formulated in the language of many-body effects. They quantify the properties of subunits (fragments), referred to as one-body properties, pairwise interactions (two-body properties), couplings of two-body interactions described by three-body properties, and so on. By introducing the notion of fragments in the framework of quantum chemistry, one obtains two immense benefits: (a) chemists can finally relate to quantum chemistry, which now speaks their language, by discussing chemically interesting subunits and their interactions and (b) calculations become much faster due to a reduced computational scaling. For instance, the somewhat academic sounding question of the importance of three-body effects in water clusters is actually another way of asking how two hydrogen bonds affect each other, when they involve three water molecules. One aspect of this is the many-body charge transfer (CT), because the charge transfers in the two hydrogen bonds are coupled to each other (not independent). In this work, we provide a generalized view on the use of many-body expansions in fragment-based methods, focusing on the general aspects of the property expansion and a contraction of a many-body expansion in a formally two-body series, as exemplified in the development of the fragment molecular orbital (FMO) method. Fragment-based methods have been very successful in delivering the properties of fragments, as well as the fragment interactions, providing insights into complex chemical processes in large molecular systems. We briefly review geometry optimizations performed with fragment-based methods and present an efficient geometry optimization method based on the combination of FMO with molecular mechanics (MM), applied to the complex of a subunit of protein kinase 2 (CK2) with a ligand. FMO results are discussed in comparison with experimental and MM-optimized structures.
Collapse
Affiliation(s)
- Dmitri G. Fedorov
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Naoya Asada
- Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Isao Nakanishi
- Department of Pharmaceutical Sciences, Kinki University, 3-4-1,
Kowakae, Higashi-Osaka, Osaka 577-8502, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
31
|
Structures and electronic properties of metal organic frameworks: DFT and ab initio FMO calculations for model systems. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.08.012] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
32
|
Fujita T, Huh J, Saikin SK, Brookes JC, Aspuru-Guzik A. Theoretical characterization of excitation energy transfer in chlorosome light-harvesting antennae from green sulfur bacteria. PHOTOSYNTHESIS RESEARCH 2014; 120:273-289. [PMID: 24504540 DOI: 10.1007/s11120-014-9978-7] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/10/2013] [Accepted: 01/23/2014] [Indexed: 06/03/2023]
Abstract
We present a theoretical study of excitation dynamics in the chlorosome antenna complex of green photosynthetic bacteria based on a recently proposed model for the molecular assembly. Our model for the excitation energy transfer (EET) throughout the antenna combines a stochastic time propagation of the excitonic wave function with molecular dynamics simulations of the supramolecular structure and electronic structure calculations of the excited states. We characterized the optical properties of the chlorosome with absorption, circular dichroism and fluorescence polarization anisotropy decay spectra. The simulation results for the excitation dynamics reveal a detailed picture of the EET in the chlorosome. Coherent energy transfer is significant only for the first 50 fs after the initial excitation, and the wavelike motion of the exciton is completely damped at 100 fs. Characteristic time constants of incoherent energy transfer, subsequently, vary from 1 ps to several tens of ps. We assign the time scales of the EET to specific physical processes by comparing our results with the data obtained from time-resolved spectroscopy experiments.
Collapse
Affiliation(s)
- Takatoshi Fujita
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, 02138, USA,
| | | | | | | | | |
Collapse
|
33
|
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
|
34
|
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
|
35
|
Yoshikawa T, Kobayashi M, Fujii A, Nakai H. Novel Approach to Excited-State Calculations of Large Molecules Based on Divide-and-Conquer Method: Application to Photoactive Yellow Protein. J Phys Chem B 2013; 117:5565-73. [DOI: 10.1021/jp401819d] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Masato Kobayashi
- Waseda Institute for Advanced Study, Waseda University, Tokyo 169-8050, Japan
| | | | - Hiromi Nakai
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama
332-0012, Japan
- ESICB, Kyoto University, Kyotodaigaku-Katsura, Nishigyoku,
Kyoto 615-8520, Japan
| |
Collapse
|
36
|
Nakata H, Nagata T, Fedorov DG, Yokojima S, Kitaura K, Nakamura S. Analytic second derivatives of the energy in the fragment molecular orbital method. J Chem Phys 2013; 138:164103. [DOI: 10.1063/1.4800990] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Nakata H, Fedorov DG, Nagata T, Yokojima S, Ogata K, Kitaura K, Nakamura S. Unrestricted Hartree-Fock based on the fragment molecular orbital method: Energy and its analytic gradient. J Chem Phys 2012; 137:044110. [DOI: 10.1063/1.4737860] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Nagata T, Fedorov DG, Li H, Kitaura K. Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method. J Chem Phys 2012; 136:204112. [DOI: 10.1063/1.4714601] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Nagata T, Fedorov DG, Kitaura K. Analytic gradient and molecular dynamics simulations using the fragment molecular orbital method combined with effective potentials. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1136-8] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
40
|
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
|
41
|
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]
|
42
|
Pruitt SR, Addicoat MA, Collins MA, Gordon MS. The fragment molecular orbital and systematic molecular fragmentation methods applied to water clusters. Phys Chem Chem Phys 2012; 14:7752-64. [DOI: 10.1039/c2cp00027j] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
43
|
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
|
44
|
Brorsen K, Fedorov DG. Fully analytic energy gradient in the fragment molecular orbital method. J Chem Phys 2011; 134:124115. [DOI: 10.1063/1.3568010] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
|
46
|
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]
|
47
|
Ishikawa T, Yamamoto N, Kuwata K. Partial energy gradient based on the fragment molecular orbital method: Application to geometry optimization. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.09.071] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
48
|
Chiba M, Koido T. Electronic excitation energy calculation by the fragment molecular orbital method with three-body effects. J Chem Phys 2010; 133:044113. [DOI: 10.1063/1.3462247] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Nagata T, Fedorov DG, Kitaura K. Importance of the hybrid orbital operator derivative term for the energy gradient in the fragment molecular orbital method. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.04.043] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
50
|
Fedorov DG, Kitaura K. The role of the exchange in the embedding electrostatic potential for the fragment molecular orbital method. J Chem Phys 2010; 131:171106. [PMID: 19894991 DOI: 10.1063/1.3250349] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have examined the role of the exchange in describing the electrostatic potential in the fragment molecular orbital method and showed that it should be included in the total Fock matrix to obtain an accurate one-electron spectrum; however, adding it to the Fock matrices of individual fragments and pairs leads to very large errors. For the error analysis we have used the virial theorem; numerical tests have been performed for solvated phenol at the Hartree-Fock level with the 6-31G( *) and 6-311G( * *) basis sets.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan.
| | | |
Collapse
|