1
|
Bhat V, Ganapathysubramanian B, Risko C. Rapid Estimation of the Intermolecular Electronic Couplings and Charge-Carrier Mobilities of Crystalline Molecular Organic Semiconductors through a Machine Learning Pipeline. J Phys Chem Lett 2024:7206-7213. [PMID: 38973725 DOI: 10.1021/acs.jpclett.4c01309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Organic semiconductors (OSC) offer tremendous potential across a wide range of (opto)electronic applications. OSC development, however, is often limited by trial-and-error design, with computational modeling approaches deployed to evaluate and screen candidates through a suite of molecular and materials descriptors that generally require hours to days of computational time to accumulate. Such bottlenecks slow the pace and limit the exploration of the vast chemical space comprising OSC. When considering charge-carrier transport in OSC, a key parameter of interest is the intermolecular electronic coupling. Here, we introduce a machine learning (ML) model to predict intermolecular electronic couplings in organic crystalline materials from their three-dimensional (3D) molecular geometries. The ML predictions take only a few seconds of computing time compared to hours by density functional theory (DFT) methods. To demonstrate the utility of the ML predictions, we deploy the ML model in conjunction with mathematical formulations to rapidly screen the charge-carrier mobility anisotropy for more than 60,000 molecular crystal structures and compare the ML predictions to DFT benchmarks.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| | - Baskar Ganapathysubramanian
- Department of Mechanical Engineering & Translational AI Center, Iowa State University, Ames, Iowa 50010, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
2
|
Nakata H, Kitoh-Nishioka H, Sakai W, Choi CH. Toward Accurate Prediction of Ion Mobility in Organic Semiconductors by Atomistic Simulation. J Chem Theory Comput 2023; 19:1517-1528. [PMID: 36757219 DOI: 10.1021/acs.jctc.2c01221] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
A multiscale scheme (MLMS: Multi-Level Multi-Scale) to predict the ion mobility (μ) of amorphous organic semiconductors is proposed, which was successfully applied to the hole mobility predictions of 14 organic systems. An inverse relationship between μ and reorganization energy is observed due to local polaronic distortions. Another moderate inverse correlation between μ and distribution of site energy change exists, representing the effects of geometric flexibility. The former and the latter represent the intramolecular and intermolecular geometric effects, respectively. In addition, a linear correlation between transfer coupling and μ is observed, showing the importance of orbital overlaps between monomers. Especially, the highest hole mobility of C6-2TTN is due to its large transfer coupling. On the other hand, another high hole mobility of CBP turned out to come from the high first neighbor density (ρFND) of its first self-solvation, emphasizing the proper description of amorphous structural configurations with a sufficiently large number of monomers. In general, systems with either unusually high transfer coupling or high first neighbor density can potentially have high μ regardless of geometric effects. Especially, the newly suggested design parameter, ρFND, is pointing to a new direction as opposed to the traditional π-conjugation strategy.
Collapse
Affiliation(s)
- Hiroya Nakata
- Research Institute for Advanced Materials and Devices, Kyocera Corporation, 3-5-3 Hikaridai Seika-cho, Soraku-gun, Kyoto 619-0237, Japan
| | - Hirotaka Kitoh-Nishioka
- Department of Energy and Materials, Faculty of Science and Engineering, Kindai University, 3 Chome-4-1 Kowakae, Higashiosaka, Osaka 577-8502, Japan
| | - Wakana Sakai
- Research Institute for Advanced Materials and Devices, Kyocera Corporation, 3-5-3 Hikaridai Seika-cho, Soraku-gun, Kyoto 619-0237, Japan
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
3
|
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
|
4
|
Kitoh-Nishioka H, Shigeta Y, Ando K. Tunneling matrix element and tunneling pathways of protein electron transfer calculated with a fragment molecular orbital method. J Chem Phys 2020; 153:104104. [PMID: 32933280 DOI: 10.1063/5.0018423] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Practical ways to calculate the tunneling matrix elements and analyze the tunneling pathways for protein electron-transfer (ET) reactions with a fragment molecular orbital (FMO) method are presented. The straightforward use of minimal basis sets only for the atoms involved in the covalent bond detachment in FMO can properly describe the ETs through the protein main-chains with the cost-effective two-body corrections (FMO2) without losing the quality of double-zeta basis sets. The current FMO codes have been interfaced with density functional theory, polarizable continuum model, and model core potentials, with which the FMO-based protein ET calculations can consider the effects of electron correlation, solvation, and transition-metal redox centers. The reasonable performance of the FMO-based ET calculations is demonstrated for three different sets of protein-ET model molecules: (1) hole transfer between two tryptophans covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, (2) ET between two plastoquinones covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, and (3) hole transfer between ruthenium (Ru) and copper (Cu) complexes covalently bridged by a stretch of a polyglycine linker as a model for Ru-modified derivatives of azurin.
Collapse
Affiliation(s)
- Hirotaka Kitoh-Nishioka
- Japan Science and Technology Agency (JST), PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Koji Ando
- Department of Information and Sciences, Tokyo Woman's Christian University, 2-6-1 Zenpukuji, Suginami-ku, Tokyo 167-8585, Japan
| |
Collapse
|
5
|
Sato R, Kitoh-Nishioka H, Ando K, Yamato T. Electron Transfer Pathways of Cyclobutane Pyrimidine Dimer Photolyase Revisited. J Phys Chem B 2018; 122:6912-6921. [PMID: 29890068 DOI: 10.1021/acs.jpcb.8b04333] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The photoinduced electron transfer (ET) reaction of cyclobutane pyrimidine dimer (CPD) photolyase plays an essential role in its DNA repair reaction, and the molecular mechanism of the ET reaction has attracted a large number of experimental and theoretical studies. We investigated the quantum mechanical nature of their ET reactions, characterized by multiple ET pathways of the CPD photolyase derived from Anacystis nidulans. Using the generalized Mulliken-Hush (GMH) method and the bridge green function (GF) methods, we estimated the electronic coupling matrix element, TDA, to be 36 ± 30 cm-1 from the donor (FADH-) to the acceptor (CPD). The estimated ET time was 386 ps, in good agreement with the experimental value (250 ps) in the literature. Furthermore, we performed the molecular dynamics (MD) simulations and ab initio molecular orbital (MO) calculations, and explored the electron tunneling pathway. We examined 20 different structures during the MD trajectory and quantitatively evaluated the electron tunneling currents for each of them. As a result, we demonstrated that the ET route via Asn349 was the dominant pathway among the five major routes via (Adenine/Asn349), (Adenine/Glu283), (Adenine/Glu283/Asn349/Met353), (Met353/Asn349), and (Asn349), indicating that Asn349 is an essential amino acid residue in the ET reaction.
Collapse
Affiliation(s)
- Ryuma Sato
- Department of Physics, Graduate School of Science , Nagoya University , Furo-cho, Chikusa-ku , Nagoya 464-8602 , Japan
| | - Hirotaka Kitoh-Nishioka
- Center for Computational Sciences , University of Tsukuba 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Koji Ando
- Department of Information and Sciences , Tokyo Woman's Christian University , 2-6-1 Zempukuji, Suginami-ku , Tokyo 167-8585 , Japan
| | - Takahisa Yamato
- Department of Physics, Graduate School of Science , Nagoya University , Furo-cho, Chikusa-ku , Nagoya 464-8602 , Japan.,Institute of Genetics and Molecular and Cellular Biology , University of Strasbourg , 1 rue Laurent Fries Parc d'Innovation 67404 Illkirch, Cedex, France
| |
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, Kitaura K. Many-body expansion of the Fock matrix in the fragment molecular orbital method. J Chem Phys 2017; 147:104106. [DOI: 10.1063/1.5001018] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
8
|
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
|
9
|
Kitoh-Nishioka H, Ando K. FMO3-LCMO study of electron transfer coupling matrix element and pathway: Application to hole transfer between two tryptophans through cis- and trans-polyproline-linker systems. J Chem Phys 2016. [DOI: 10.1063/1.4962626] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hirotaka Kitoh-Nishioka
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan
| | - Koji Ando
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
10
|
Fedorov DG, Kitaura K. Subsystem Analysis for the Fragment Molecular Orbital Method and Its Application to Protein-Ligand Binding in Solution. J Phys Chem A 2016; 120:2218-31. [PMID: 26949816 DOI: 10.1021/acs.jpca.6b00163] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A subsystem analysis is derived incorporating interfragment interactions into the fragment properties, such as energies or charges. The relative stabilities of three alanine isomers, the α-helix, the β-turn, and the extended form are studied and the differences in fragment properties are elucidated. The analysis is further elaborated for studies of binding energies. The binding of the Trp-cage protein (PDB: 1L2Y ) to two ligands is studied in detail. Binding energies defined for each fragment can be used as a convenient descriptor for analyzing contributions to binding in solution.
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
| | - Kazuo Kitaura
- Graduate School of System Informatics, Kobe University , 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
11
|
Otsuka T, Okimoto N, Taiji M. Assessment and acceleration of binding energy calculations for protein-ligand complexes by the fragment molecular orbital method. J Comput Chem 2015; 36:2209-18. [DOI: 10.1002/jcc.24055] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2015] [Revised: 07/17/2015] [Accepted: 07/30/2015] [Indexed: 12/16/2022]
Affiliation(s)
- Takao Otsuka
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| | - Noriaki Okimoto
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| | - Makoto Taiji
- Laboratory for Computational Molecular Design, Computational Biology Research Core, RIKEN Quantitative Biology Center (QBiC); 1-6-5 Minatojima Minamimachi, Chuo-Ku, Kobe Hyogo 650-0047 Japan
| |
Collapse
|
12
|
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
|
13
|
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
|
14
|
Kitoh-Nishioka H, Ando K. Charge-transfer matrix elements by FMO-LCMO approach: Hole transfer in DNA with parameter tuned range-separated DFT. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2014.12.057] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
15
|
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
|
16
|
Shimojo F, Hattori S, Kalia RK, Kunaseth M, Mou W, Nakano A, Nomura KI, Ohmura S, Rajak P, Shimamura K, Vashishta P. A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations. J Chem Phys 2014; 140:18A529. [DOI: 10.1063/1.4869342] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
17
|
Fukuzawa K, Watanabe C, Kurisaki I, Taguchi N, Mochizuki Y, Nakano T, Tanaka S, Komeiji Y. Accuracy of the fragment molecular orbital (FMO) calculations for DNA: Total energy, molecular orbital, and inter-fragment interaction energy. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.02.002] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|