1
|
Ohno S, Uratani H, Nakai H. Implementation of Nonadiabatic Molecular Dynamics for Intersystem Crossing Based on a Time-Dependent Density-Functional Tight-Binding Method. J Phys Chem A 2024. [PMID: 38990848 DOI: 10.1021/acs.jpca.4c02422] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
Intersystem crossing (ISC) and internal conversion (IC) are types of nonadiabatic transitions that play important roles in a wide range of fields, including photochemistry, photophysics, and photobiology. The nonadiabatic molecular dynamics (NA-MD) method is a powerful tool for computational simulations of dynamic phenomena involving nonadiabatic transitions. In this study, we implemented the NA-MD method, which treats ISC and IC on an equal footing, where the electronic structure is treated at the level of the time-dependent (TD) density-functional tight-binding (DFTB) method, a low-cost semiempirical analog of TD density functional theory (DFT). In particular, the spin-orbit coupling calculation algorithm was implemented in the TD-DFTB framework, and the results showed trends similar to those obtained using TD-DFT. In addition, the NA-MD method successfully reproduced ultrafast ISC of 2-nitronaphthalene.
Collapse
Affiliation(s)
- Shota Ohno
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiroki Uratani
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto Daigaku-Katsura, Nishikyo, Kyoto 615-8510, Japan
- PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, 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
|
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
|
3
|
Sakhraoui T, Karlický F. Prediction of induced magnetism in 2D Ti 2C based MXenes by manipulating the mixed surface functionalization and metal substitution computed by xTB model Hamiltonian of the DFTB method. Phys Chem Chem Phys 2024; 26:12862-12868. [PMID: 38623885 DOI: 10.1039/d3cp05665a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
We employed the recently developed density functional tight binding (DFTB) method's Hamiltonian, GFN1-xTB, for modeling the mixed termination in Ti2C MXenes, namely three types of termination by combining -O and -OH, -O and -F, and -F and -OH. We demonstrated that the approach yields reliable predictions for the electronic and magnetic properties of such MXenes. The first highlighted result is that the mixed surface functionalization in Ti2CAxBy MXenes induces spin polarization with diverse magnetic alignments, including ferromagnetism and two types of antiferromagnetism. We further identified the magnetic alignment for the investigated MXene in terms of the compositions of the terminal groups. Moreover, the effect of the transition metal (Ti) substituted by the Sc atom on the electronic and magnetic properties was also investigated. We found that the studied systems maintain the magnetism and the metallic characteristics. A magnetic transition from antiferromagnetic (AFM) to ferrimagnetic (FiM) ordering was found for ScTi15C8F8(OH)8 and ScTi15C8F12(OH)4 compounds. Finally, we proved that incorporating the Sc atom into the lattice of Ti2CO2 and the mixed surface termination in Ti2CAxBy is an effective strategy to induce magnetism. Our study may provide a new potential application for designing MXene-based spintronics.
Collapse
Affiliation(s)
- Taoufik Sakhraoui
- Department of Physics, Faculty of Science, University of Ostrava, 701 03 Ostrava, Czech Republic.
| | - František Karlický
- Department of Physics, Faculty of Science, University of Ostrava, 701 03 Ostrava, Czech Republic.
| |
Collapse
|
4
|
Wickramasinghe S, Hoehn A, Wetthasinghe ST, Lin H, Wang Q, Jakowski J, Rassolov V, Tang C, Garashchuk S. Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes. J Phys Chem B 2023; 127:10129-10141. [PMID: 37972315 DOI: 10.1021/acs.jpcb.3c04118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. The effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.
Collapse
Affiliation(s)
- Sachith Wickramasinghe
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Alexandria Hoehn
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Shehani T Wetthasinghe
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Huina Lin
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Qi Wang
- Department of Mathematics, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Jacek Jakowski
- Center for Nanophase Materials Science, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Vitaly Rassolov
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Chuanbing Tang
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Sophya Garashchuk
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| |
Collapse
|
5
|
Putro PA, Sakti AW, Ahmad F, Nakai H, Alatas H. Quantum mechanical assessment on the optical properties of capsanthin conformers. J Comput Chem 2023; 44:2319-2331. [PMID: 37548072 DOI: 10.1002/jcc.27199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Revised: 06/25/2023] [Accepted: 07/15/2023] [Indexed: 08/08/2023]
Abstract
As optical properties, the ultraviolet-visible (UV-Vis) absorption spectra of capsanthin-based red natural dye are a decisive parameter for their usage in various applications. Thus, accurately predicting the maximum UV-Vis wavelength (λ max ) values is critical in designing dye-conjugated material. Extensive metadynamics simulations were carried out to generate capsanthin conformers at various levels of the extended tight-binding method. Benchmarking the time-dependent density-functional theory (TD-DFT) methods help understand the results of a particular functional and allows a comparison between results obtained with different functional. The long-range correction (LC) scheme in LC-TD-DFT-D4/ωB97X/def2-SVP has been found to reproduce the experimentalλ max , and exhibited the effect of conformational changes to the calculated wavelengths. On the other hand, an inexpensive yet efficient LC-TD-DFTB method reproduced the experimentalλ max insensitive to conformational changes.
Collapse
Affiliation(s)
- Permono Adi Putro
- Theoretical Physics Division, Department of Physics, Faculty of Mathematics and Natural Sciences, IPB University, Bogor, Indonesia
- Department of Physics, Faculty of Science, Universitas Mandiri, Subang, Indonesia
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor, Indonesia
| | - Aditya Wibawa Sakti
- Theoretical Physics Division, Department of Physics, Faculty of Mathematics and Natural Sciences, IPB University, Bogor, Indonesia
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor, Indonesia
- Global Center for Science and Engineering, Waseda University, Tokyo, Japan
- Department of Chemistry, Faculty of Science and Computer, Universitas Pertamina, Jakarta, Indonesia
| | - Faozan Ahmad
- Theoretical Physics Division, Department of Physics, Faculty of Mathematics and Natural Sciences, IPB University, Bogor, Indonesia
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor, Indonesia
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, Waseda University, Tokyo, Japan
| | - Husin Alatas
- Theoretical Physics Division, Department of Physics, Faculty of Mathematics and Natural Sciences, IPB University, Bogor, Indonesia
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor, Indonesia
- Center for Transdisciplinary & Sustainability Sciences (CTSS), IPB University, Bogor, Indonesia
| |
Collapse
|
6
|
Budiutama G, Li R, Manzhos S, Ihara M. Hybrid Density Functional Tight Binding (DFTB)─Molecular Mechanics Approach for a Low-Cost Expansion of DFTB Applicability. J Chem Theory Comput 2023. [PMID: 37450317 DOI: 10.1021/acs.jctc.3c00310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
The density functional-based tight binding (DFTB) method has seen a rise in adoption for materials modeling, as it offers significant improvement in scalability with accuracy comparable to the density functional theory (DFT) when good parameterizations exist. The cost reduction in DFTB compared to DFT is achieved by the pre-parameterization of the elements of the Hamiltonian matrix as well as the repulsion potential between all pairs of atoms. Parameterization for new systems with accuracies competitive with DFT in specific applications requires specialized manpower and computational resources. This prevents the application of the DFTB method to systems for which it was not parameterized. In this work, we explore an approach to address the problem of missing parameters of DFTB by modeling the interactions with missing Slater-Koster parameters with an interatomic interaction potential. When the distance between two atoms modeled at the force-field level is sufficiently large, the approach results in accurate structural and electronic properties. The resulting calculation is therefore a hybrid between DFTB and molecular mechanics, a pure DFTB for atoms with a complete set of interatomic parameterizations, and a mix between DFTB and molecular mechanics for atoms with a missing interatomic parameterization. The approach is expected to be particularly useful for hybrid materials and interfaces. The method is tested on the examples of 2D materials, mixed oxides, and a large-scale calculation of an oxide-oxide interface.
Collapse
Affiliation(s)
- Gekko Budiutama
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Ruicheng Li
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| |
Collapse
|
7
|
Kubař T, Elstner M, Cui Q. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annu Rev Biophys 2023; 52:525-551. [PMID: 36791746 PMCID: PMC10810093 DOI: 10.1146/annurev-biophys-111622-091140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.
Collapse
Affiliation(s)
- T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - Q Cui
- Department of Chemistry, Boston University, Boston, Massachusetts, USA;
- Department of Physics, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| |
Collapse
|
8
|
Uratani H, Nakai H. Nanoscale and Real-Time Nuclear-Electronic Dynamics Simulation Study of Charge Transfer at the Donor-Acceptor Interface in Organic Photovoltaics. J Phys Chem Lett 2023; 14:2292-2300. [PMID: 36827224 DOI: 10.1021/acs.jpclett.2c03808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Charge-transfer (CT) processes in donor-acceptor interfaces of organic photovoltaics have been challenging targets for computational chemistry owing to their nanoscale and ultrafast nature. Herein, we report real-time nuclear-electronic dynamics simulations of CT processes in a nanometer-scale donor-acceptor interface model composed of a donor poly(3-hexylthiophene-2,5-diyl) crystal and an acceptor [6,6]-phenyl-C61-butyric acid methyl ester aggregate. The simulations were realized using our original reduced-scaling computational technique, namely, patchwork-approximation-based Ehrenfest dynamics. The results illustrated the CT pathway with atomic resolution, thereby rationalizing the observed excitation-energy dependence of the quantity of CT. Further, nuclear motion, which is affected by the electronic dynamics, was observed to play a significant role in the CT process by modulating molecular orbital energies. The present study suggests that microscopic CT processes strongly depend on local structures of disordered donor-acceptor interfaces as well as coupling between nuclear and electronic dynamics.
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
| | - 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
| |
Collapse
|
9
|
Nishimura Y, Nakai H. Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations. J Chem Phys 2023; 158:054106. [PMID: 36754823 DOI: 10.1063/5.0132573] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.
Collapse
Affiliation(s)
- 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
| |
Collapse
|
10
|
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
|
11
|
Ono J, Okada C, Nakai H. Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin. Chemphyschem 2022; 23:e202200109. [PMID: 35818319 DOI: 10.1002/cphc.202200109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 06/21/2022] [Indexed: 11/10/2022]
Abstract
In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (~50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.
Collapse
Affiliation(s)
- Junichi Ono
- Kyoto University: Kyoto Daigaku, Elements Strategy Initiative for Catalysts & Batteries (ESICB), 1-30 Goryo-Ohara, 615-8245, Nishi-ku, JAPAN
| | - Chika Okada
- Waseda University: Waseda Daigaku, Department of Chemistry and Biochemistry, 3-4-1 Okubo, 169-8555, Shinjuku, JAPAN
| | - Hiromi Nakai
- Waseda University Faculty of Science and Engineering: Waseda Daigaku Riko Gakujutsuin, Department of Chemistry and Biochemistry, 3-4-1 Okubo, 169-8555, Shinjuku, JAPAN
| |
Collapse
|
12
|
Peng L, Peng D, Gu FL, Yang W. Regularized Localized Molecular Orbitals in a Divide-and-Conquer Approach for Linear Scaling Calculations. J Chem Theory Comput 2022; 18:2975-2982. [PMID: 35416665 PMCID: PMC9972215 DOI: 10.1021/acs.jctc.2c00142] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Non-orthogonal localized molecular orbitals (NOLMOs) have been employed as building blocks for the divide-and-conquer (DC) linear scaling method. The NOLMOs are calculated from subsystems and used for constructing the density matrix (DM) of the entire system, instead of the subsystem DM in the original DC approach. Also, unlike the original DC method, the inverse electronic temperature parameter β is not needed anymore. Furthermore, a new regularized localization approach for NOLMOs has been developed, in which the localization cost function is a sum of the spatial spread function, as in the Boys method, and the kinetic energy, as a regularization measure to limit the oscillation of the NOLMOs. The optimal weight of the kinetic energy can be determined by optimization with analytical gradients. The resulting regularized NOLMOs have enhanced smoothness and better transferability because of reduced kinetic energies. Compared with the original DC, while NOLMO-DC has a similar computational linear scaling cost, the accuracy of NOLMO-DC is better by several orders of magnitude for large conjugated systems and by about 1 order of magnitude for other systems. The NOLMO-DC method is thus a promising development of the DC approach for linear scaling calculations.
Collapse
Affiliation(s)
- Liang Peng
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Daoling Peng
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Environment, South China Normal University, Guangzhou 510006, People’s Republic of China
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708-0346, United States
| |
Collapse
|
13
|
Multiple protonation states in ligand-free SARS-CoV-2 main protease revealed by large-scale quantum molecular dynamics simulations. Chem Phys Lett 2022; 794:139489. [PMID: 35221345 PMCID: PMC8863314 DOI: 10.1016/j.cplett.2022.139489] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2022] [Revised: 02/16/2022] [Accepted: 02/20/2022] [Indexed: 12/16/2022]
Abstract
The main protease (Mpro) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) catalyzes the cleavage of polyproteins for viral replication. Here, large-scale quantum molecular dynamics and metadynamics simulations for ligand-free Mpro were performed, where all the atoms were treated quantum-mechanically, focusing on elucidation of the controversial active-site protonation state. The simulations clarified that the interconverting multiple protonation states exist in unliganded Mpro, and the catalytically relevant ion-pair state is more stable than the neutral state, which is consistent with neutron crystallography. The results highlight the importance of the ion-pair state for repurposing or discovering antiviral drugs that target Mpro.
Collapse
|
14
|
Chan B, Dawson W, Nakajima T. Searching for a Reliable Density Functional for Molecule-Environment Interactions, Found B97M-V/def2-mTZVP. J Phys Chem A 2022; 126:2397-2406. [PMID: 35390254 DOI: 10.1021/acs.jpca.2c02032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present study, we have examined density functional theory methods for the calculation of the interaction energy between a small molecule and its environment. For simple systems such as a neutral solute in a neutral solvent, good accuracy can be attained using low-cost "3c" methods, in particular r2SCAN-3c. When part(s) of the system is charged, the accurate computation of the interactions is more challenging. In these cases, we find the B97M-V/def2-mTZVP method to agree well with reference values; it also shows good accuracy for the more straightforward neutral systems. Thus, B97M-V/def2-mTZVP provides a means for accurate and low-cost computation of interaction energies, notably the binding between a substrate or a drug molecule and an enzyme, which may facilitate rational drug design.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - William Dawson
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| |
Collapse
|
15
|
Chan B, Dawson W, Nakajima T. Modeling the Conformational Preference of the Lignocellulose Interface and Its Interaction with Weak Acids. J Phys Chem A 2022; 126:2119-2126. [PMID: 35349294 DOI: 10.1021/acs.jpca.2c00308] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
We have examined the conformational space of model systems for the hydrogen-bonded and covalent linkages between the sugar and lignin components of lignocellulose. Specifically, glucose and paracoumaryl alcohol moieties are used in our models. Multistage screening protocols are used to identify and validate a set of lowest-energy isomers. We found a cost-effective screening process involving an initial screening with DFTB3/3ob using a 20 kJ mol-1 threshold, a refinement with the SCANh/6-31+G(2d,p) method with a 10 kJ mol-1 cutoff, a third step at the DSD-PBEP86/ma-def2-TZVP//MS1-D3/6-31+G(2d,p) level with the same 10 kJ mol-1 threshold, and a last step at the CCSD(T)/CBS//B3LYP/cc-pVTZ level with a tighter 5 kJ mol-1 threshold. The use of machine learning (with the London and Axilrod-Teller-Muto potential) can further accelerate the screening process. In general, all low-energy conformers are characterized by hydrogen bonding between the sugar and lignin moieties. We examined the interactions of covalently bonded sugar-lignin models with weak acids (HSO3-, H2PO3-, HSeO3-, H2citrate-, etc.) and found that they interacted strongly with the oxygen of the sugar-O-lignin linkage. Our results suggest that acids such as dihydrogen citrate may be attractive alternatives to the commonly used HSO3- for lignocellulose processing.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki, Nagasaki 852-8521, Japan
| | - William Dawson
- RIKEN Center for Computational Science, Minatojima-minami 7-1-26, Kobe, Hyogo 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, Minatojima-minami 7-1-26, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
16
|
Sakti AW, Wahyudi ST, Ahmad F, Darmawan N, Hardhienata H, Alatas H. Effects of Salt Concentration on the Water and Ion Self-Diffusion Coefficients of a Model Aqueous Sodium-Ion Battery Electrolyte. J Phys Chem B 2022; 126:2256-2264. [PMID: 35271293 DOI: 10.1021/acs.jpcb.1c09619] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The aqueous sodium-ion battery is a promising alternative to the well-known lithium-ion battery owing to the large abundance of sodium ion resources. Although it is safer than the lithium-ion battery, the voltage window of the sodium-ion battery is narrower than that of the lithium-ion battery, thus limiting its practical implementation. Therefore, a highly concentrated electrolyte is required to address this issue. In the present work, the effect of the salt concentration on the transport properties of water molecules is investigated via theoretical analyses at the quantum mechanical level. A molecular dynamics simulation at the quantum mechanical level revealed that as the salt concentration increases, the ion-water interactions became stronger, leading to a lower diffusivity and a lower electronic band gap. These imply that the superconcentrated aqueous-based electrolytes have high potentials for the sodium-ion battery applications.
Collapse
Affiliation(s)
- Aditya Wibawa Sakti
- Department of Chemistry, Faculty of Science and Computer, Universitas Pertamina, Jakarta 12220, Indonesia.,Theoretical Physics Division, Department of Physics, IPB University, Bogor 16680, Indonesia.,Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan.,Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia
| | - Setyanto Tri Wahyudi
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia.,Biophysics Division, Department of Physics, IPB University, Bogor 16680, Indonesia
| | - Faozan Ahmad
- Theoretical Physics Division, Department of Physics, IPB University, Bogor 16680, Indonesia.,Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia
| | - Noviyan Darmawan
- Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia.,Inorganic Chemistry Division, Department of Chemistry, IPB University, Bogor 16680, Indonesia
| | - Hendradi Hardhienata
- Theoretical Physics Division, Department of Physics, IPB University, Bogor 16680, Indonesia.,Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia
| | - Husin Alatas
- Theoretical Physics Division, Department of Physics, IPB University, Bogor 16680, Indonesia.,Indonesia Computational-Research Consortium on Renewable Energy (ICRC-RE), IPB University, Bogor 16680, Indonesia
| |
Collapse
|
17
|
Kitheka M, Redington M, Zhang J, Yao Y, Goyal P. BENCHMARKS OF THE DENSITY FUNCTIONAL TIGHT-BINDING METHOD FOR REDOX, PROTONATION AND ELECTRONIC PROPERTIES OF QUINONES. Phys Chem Chem Phys 2022; 24:6742-6756. [DOI: 10.1039/d1cp05333g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Organic materials with controllable molecular design and sustainable resources are promising electrode materials. Crystalline quinones have been investigated in a variety of rechargeable battery chemistries due to their ubiquitous nature,...
Collapse
|
18
|
Khanal R, Irle S. Quantum chemical investigation of the effect of alkali metal ions on the dynamic structure of water in aqueous solutions. RSC Adv 2022; 12:25500-25510. [PMID: 36275866 PMCID: PMC9480497 DOI: 10.1039/d2ra04563j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 08/26/2022] [Indexed: 11/21/2022] Open
Abstract
Cation size primarily governs dynamic correlations in aqueous solutions. However, drawing a clear line between chaotropic and cosmotrope ions is complicated by the fact that all cations delay dynamic correlation decay.
Collapse
Affiliation(s)
- Rabi Khanal
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 37831, USA
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 37831, USA
| |
Collapse
|
19
|
Uratani H, Nakai H. Scalable Ehrenfest Molecular Dynamics Exploiting the Locality of Density-Functional Tight-Binding Hamiltonian. J Chem Theory Comput 2021; 17:7384-7396. [PMID: 34860019 DOI: 10.1021/acs.jctc.1c00950] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
To explore the science behind excited-state dynamics in high-complexity chemical systems, a scalable nonadiabatic molecular dynamics (MD) technique is indispensable. In this study, by treating the electronic degrees of freedom at the density-functional tight-binding level, we developed and implemented a reduced scaling and multinode-parallelizable Ehrenfest MD method. To achieve this goal, we introduced a concept called patchwork approximation (PA), where the effective Hamiltonian for real-time propagation of the electronic density matrix is partitioned into a set of local parts. Numerical results for giant icosahedral fullerenes, which comprise up to 6000 atoms, suggest that the scaling of the present PA-based method is less than quadratic, which yields a significant advantage over the conventional cubic scaling method in terms of computational time. The acceleration by the parallelization on multiple nodes was also assessed. Furthermore, the electronic and structural dynamics resulting from the perturbation by the external electric field were accurately reproduced with the PA, even when the electronic excitation was spatially delocalized.
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
| | - 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
|
20
|
Matsumoto RA, Thompson MW, Vuong VQ, Zhang W, Shinohara Y, van Duin ACT, Kent PRC, Irle S, Egami T, Cummings PT. Investigating the Accuracy of Water Models through the Van Hove Correlation Function. J Chem Theory Comput 2021; 17:5992-6005. [PMID: 34516134 DOI: 10.1021/acs.jctc.1c00637] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present molecular-simulation-based calculations of the Van Hove correlation function (VHF) of water using multiple modeling approaches: classical molecular dynamics with simple three-site nonpolarizable models, with a polarizable model, and with a reactive force field; density functional tight-binding molecular dynamics; and ab initio molecular dynamics. Due to the many orders of magnitude difference in the computational cost of these approaches, we investigate how small and short the simulations can be while still yielding sufficiently accurate and interpretable results for the VHF. We investigate the accuracy of the different models by comparing them to recently published inelastic X-ray scattering measurements of the VHF. We find that all of the models exhibit qualitative agreement with the experiments, and in some models and for some properties, the agreement is quantitative. This work lays the foundation for future simulation approaches to calculating the VHF for aqueous solutions in bulk and under nanoconfinement.
Collapse
Affiliation(s)
- Ray A Matsumoto
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| | - Matthew W Thompson
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| | - Van Quan Vuong
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Weiwei Zhang
- Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Yuya Shinohara
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Adri C T van Duin
- Department of Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Paul R C Kent
- Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Takeshi Egami
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States.,Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37916, United States.,Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Peter T Cummings
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37212, United States.,Multiscale Modeling and Simulation Center, Vanderbilt University, Nashville, Tennessee 37212, United States
| |
Collapse
|
21
|
Nakai H, Takemura T, Ono J, Nishimura Y. Quantum-Mechanical Molecular Dynamics Simulations on Secondary Proton Transfer in Bacteriorhodopsin Using Realistic Models. J Phys Chem B 2021; 125:10947-10963. [PMID: 34582194 DOI: 10.1021/acs.jpcb.1c06231] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Bacteriorhodopsin (BR) transports a proton from intracellular to extracellular (EC) sites through five proton transfers. The second proton transfer is the release of an excess proton stored in BR into the EC medium, and an atomistic understanding of this whole process has remained unexplored due to its ubiquitous environment. Here, fully quantum mechanical (QM) molecular dynamics (MD) and metadynamics (MTD) simulations for this process were performed at the divide-and-conquer density-functional tight-binding level using realistic models (∼50000 and ∼20000 atoms) based on the time-resolved photointermediate structures from an X-ray free electron laser. Regarding the proton storage process, the QM-MD/MTD simulations confirmed the Glu-shared mechanism, in which an excess proton is stored between Glu194 and Glu204, and clarified that the activation occurs by localizing the proton at Glu204 in the photocycle. Furthermore, the QM-MD/MTD simulations elucidated a release pathway from Glu204 through Ser193 to the EC water molecules and clarified that the proton release starts at ∼250 μs. In the ubiquitous proton diffusion in the EC medium, the transient proton receptors predicted experimentally were assigned to carboxylates in Glu9 and Glu74. Large-scale QM-MD/MTD simulations beyond the conventional sizes, which provided the above findings and confirmations, were possible by adopting our Dcdftbmd program.
Collapse
Affiliation(s)
- 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 (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Toshiaki Takemura
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
22
|
Takenaka N, Bouibes A, Yamada Y, Nagaoka M, Yamada A. Frontiers in Theoretical Analysis of Solid Electrolyte Interphase Formation Mechanism. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2021; 33:e2100574. [PMID: 34338349 DOI: 10.1002/adma.202100574] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Revised: 05/13/2021] [Indexed: 06/13/2023]
Abstract
Solid electrolyte interphase (SEI) is an ion conductive yet electron-insulating layer on battery electrodes, which is formed by the reductive decomposition of electrolytes during the initial charge. The nature of the SEI significantly impacts the safety, power, and lifetime of the batteries. Hence, elucidating the formation mechanism of the SEI layer has become a top priority. Conventional theoretical calculations reveal initial elementary steps of electrolyte reductive decomposition, whereas experimental approaches mainly focus on the characterization of the formed SEI in the final form. Moreover, both theoretical and experimental methodologies could not approach intermediate or transient steps of SEI growth. A major breakthrough has recently been achieved through a novel multiscale simulation method, which has enriched the understanding of how the reduction products are aggregated near the electrode and influence the SEI morphologies. This review highlights recent theoretical achievements to reveal the growth mechanism and provides a clear guideline for designing a stable SEI layer for advanced batteries.
Collapse
Affiliation(s)
- Norio Takenaka
- Graduate School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
- ESICB, Kyoto University, Kyodai Katsura, Nishikyo-ku, Kyoto, 615-8520, Japan
| | - Amine Bouibes
- Graduate School of Informatics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8601, Japan
| | - Yuki Yamada
- Graduate School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
- ESICB, Kyoto University, Kyodai Katsura, Nishikyo-ku, Kyoto, 615-8520, Japan
| | - Masataka Nagaoka
- ESICB, Kyoto University, Kyodai Katsura, Nishikyo-ku, Kyoto, 615-8520, Japan
- Graduate School of Informatics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8601, Japan
| | - Atsuo Yamada
- Graduate School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-8656, Japan
- ESICB, Kyoto University, Kyodai Katsura, Nishikyo-ku, Kyoto, 615-8520, Japan
| |
Collapse
|
23
|
Sakti AW, Chou CP, Nishimura Y, Nakai H. Is Oxygen Diffusion Faster in Bulk CeO2 or on a (111)-CeO2 Surface? A Theoretical Study. CHEM LETT 2021. [DOI: 10.1246/cl.200895] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Aditya Wibawa Sakti
- Department of Chemistry, Faculty of Science and Computer, Universitas Pertamina, Jakarta 12220, Indonesia
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
| | - Chien-Pin Chou
- JSR Coorporation Yokkaichi Research Center, Yokkaichi, Mie 510-8552, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
- Element Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyotodaigaku-Katsura, Kyoto 615-8520, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| |
Collapse
|
24
|
|
25
|
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
|
26
|
Sokolov M, Bold BM, Kranz JJ, Höfener S, Niehaus TA, Elstner M. Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies. J Chem Theory Comput 2021; 17:2266-2282. [PMID: 33689344 DOI: 10.1021/acs.jctc.1c00095] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The absorption and emission of light is a ubiquitous process in chemical and biological processes, making a theoretical description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of analytical nuclear gradients for time-dependent long-range corrected density functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 organic molecules taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivatives based on the analytical gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.
Collapse
Affiliation(s)
- Monja Sokolov
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Julian J Kranz
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Thomas A Niehaus
- Université Claude Bernard Lyon 1, Université Lyon, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany.,Institute of Biological Interfaces (IBG2), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| |
Collapse
|
27
|
NAKAI H. Commentary toward the 20th Anniversary of the Society ofComputer Chemistry, Japan. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2021. [DOI: 10.2477/jccj.2021-0020] [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)
- Hiromi NAKAI
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University,3-4-1 Okubo, Shinjuku, Tokyo 169-8555, JAPAN
| |
Collapse
|
28
|
Uratani H, Morioka T, Yoshikawa T, Nakai H. Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings. J Chem Theory Comput 2020; 16:7299-7313. [PMID: 33197192 DOI: 10.1021/acs.jctc.0c00936] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.
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
| | - Toshiki Morioka
- 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
|
29
|
Ono J, Imai M, Nishimura Y, Nakai H. Hydroxide Ion Carrier for Proton Pumps in Bacteriorhodopsin: Primary Proton Transfer. J Phys Chem B 2020; 124:8524-8539. [DOI: 10.1021/acs.jpcb.0c05507] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Minori Imai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
30
|
Ono J, Nakai H. Weighted histogram analysis method for multiple short-time metadynamics simulations. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137384] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
31
|
Yoshikawa T, Doi T, Nakai H. Finite-temperature-based time-dependent density-functional theory method for static electron correlation systems. J Chem Phys 2020; 152:244111. [PMID: 32610978 DOI: 10.1063/1.5144527] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this study, we developed a time-dependent density-functional theory (TDDFT) with a finite-temperature (FT) scheme, denoted as FT-TDDFT. We introduced the concept of fractional occupation numbers for random phase approximation equation and evaluated the excited-state electronic entropy terms with excited-state occupation number. The orbital occupation numbers for the excited state were evaluated from the change in the ground-state electron configuration with excitation and deexcitation coefficients. Furthermore, we extended the FT formulation to the time-dependent density-functional tight-binding (TDDFTB) method for larger systems, denoted as FT-TDDFTB. Numerical assessment for the FT-(TD)DFT method showed smooth potential curves for double-bond rotation of ethylene in both ground and excited states. Excited-state calculations based on the FT-TDDFTB method were applied to the uniform π-stacking columns composed of trioxotriangulene, possessing neutral radicals in strong correlation systems.
Collapse
Affiliation(s)
- Takeshi Yoshikawa
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Toshiki Doi
- 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
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| |
Collapse
|
32
|
Uratani H, Nakai H. Non-adiabatic molecular dynamics with divide-and-conquer type large-scale excited-state calculations. J Chem Phys 2020; 152:224109. [DOI: 10.1063/5.0006831] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
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
| | - 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
|
33
|
Uratani H, Nakai H. Simulating the Coupled Structural-Electronic Dynamics of Photoexcited Lead Iodide Perovskites. J Phys Chem Lett 2020; 11:4448-4455. [PMID: 32418430 DOI: 10.1021/acs.jpclett.0c01028] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Motivated by the optoelectronic applications of lead halide perovskites (LHPs), researchers have paid considerable attention to their photoexcited-state dynamics, where the coupling between the electronic and nuclear dynamics is pronounced. Here, we present simulations of the photoexcited-state dynamics of representative lead iodide perovskites, CsPbI3 and MAPbI3 (MA = CH3NH3), by adopting nonadiabatic molecular dynamics combined with the linear-response time-dependent density-functional tight-binding (LR-TD-DFTB) method, an efficient excited-state calculation framework. In the calculations, the electronic wave function and the nuclear coordinates were propagated in a mutually dependent manner. The results suggest that the excited LHPs undergo exciton dissociation, hot carrier cooling, and polaron formation on similar time scales. In particular, the decay of the carrier energy is attributed to not only the relaxation toward the band edge but also the change in orbital energy originating from the structural deformation, highlighting the importance of coupling between the electronic and nuclear degrees of freedom.
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
| | - 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
|
34
|
Nishimura Y, Nakai H. Hierarchical parallelization of divide‐and‐conquer density functional tight‐binding molecular dynamics and metadynamics simulations. J Comput Chem 2020; 41:1759-1772. [DOI: 10.1002/jcc.26217] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Revised: 04/15/2020] [Accepted: 04/20/2020] [Indexed: 11/08/2022]
Affiliation(s)
- Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry School of Advanced Science and Engineering, Waseda University Tokyo Japan
- Elements Strategy Initiative for Catalysts and Batteries Kyoto University Kyoto Japan
| |
Collapse
|
35
|
Sakti A, Chou CP, Nakai H. Density-Functional Tight-Binding Study of Carbonaceous Species Diffusion on the (100)-γ-Al 2O 3 Surface. ACS OMEGA 2020; 5:6862-6871. [PMID: 32258922 PMCID: PMC7114690 DOI: 10.1021/acsomega.0c00203] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Accepted: 03/05/2020] [Indexed: 05/17/2023]
Abstract
Carbonaceous or oxy-carbon species are intermediates formed during C x H y combustion on a Pt n /Al2O3 catalyst, which contain carbon, hydrogen, and oxygen atoms. The accumulation of the carbonaceous species, arguably, leads to catalytic deactivation; therefore, their removal is of importance. As the diffusion process is occasionally the rate-determining step in the growth of carbonaceous species, the present study aims to reveal the diffusion mechanisms. The free energy barriers of acetate, formate, and methoxy diffusion on the (100)-γ-Al2O3 surface were evaluated through extensive metadynamics simulations at the density-functional tight-binding level. The present work deduces that each adopted carbonaceous species exhibits different diffusion mechanisms and supports experimental evidence that the acetate species exhibits the slowest diffusivity among the adopted carbonaceous species.
Collapse
Affiliation(s)
- Aditya
W. Sakti
- Element
Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyotodaigaku-Katsura, Kyoto 615-8520, Japan
- Waseda
Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
| | - Chien-Pin Chou
- Waseda
Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Element
Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Kyotodaigaku-Katsura, Kyoto 615-8520, Japan
- Waseda
Research Institute for Science and Engineering (WISE), Waseda University, Tokyo 169-8555, Japan
- Department
of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- E-mail: . Phone: +81 3-5286-3452. Fax: +81 3-3205-2504
| |
Collapse
|
36
|
Inamori M, Yoshikawa T, Ikabata Y, Nishimura Y, Nakai H. Spin‐flip approach within time‐dependent density functional tight‐binding method: Theory and applications. J Comput Chem 2020; 41:1538-1548. [DOI: 10.1002/jcc.26197] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2020] [Revised: 03/07/2020] [Accepted: 03/10/2020] [Indexed: 12/25/2022]
Affiliation(s)
- Mayu Inamori
- Department of Chemistry and Biochemistry, School of Advanced Science and EngineeringWaseda University Tokyo Japan
| | - Takeshi Yoshikawa
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Yasuhiro Ikabata
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and EngineeringWaseda University Tokyo Japan
- Waseda Research Institute for Science and EngineeringWaseda University Tokyo Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB)Kyoto University Kyoto Japan
| |
Collapse
|
37
|
Komoto N, Yoshikawa T, Nishimura Y, Nakai H. Large-Scale Molecular Dynamics Simulation for Ground and Excited States Based on Divide-and-Conquer Long-Range Corrected Density-Functional Tight-Binding Method. J Chem Theory Comput 2020; 16:2369-2378. [DOI: 10.1021/acs.jctc.9b01268] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- 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
| | - Takeshi Yoshikawa
- Waseda Research Institute for 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
- 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, 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
|
38
|
Spiegelman F, Tarrat N, Cuny J, Dontot L, Posenitskiy E, Martí C, Simon A, Rapacioli M. Density-functional tight-binding: basic concepts and applications to molecules and clusters. ADVANCES IN PHYSICS: X 2020; 5:1710252. [PMID: 33154977 PMCID: PMC7116320 DOI: 10.1080/23746149.2019.1710252] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Accepted: 12/19/2019] [Indexed: 06/10/2023] Open
Abstract
The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.
Collapse
Affiliation(s)
- Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Nathalie Tarrat
- CEMES, Université de Toulouse (UPS), CNRS, UPR8011, Toulouse, Toulouse, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Leo Dontot
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Evgeny Posenitskiy
- Laboratoire Collisions Agrégats et Réactivité LCAR/IRSAMC, UMR5589, Université de Toulouse (UPS) and CNRS, Toulouse, France
| | - Carles Martí
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
- Laboratoire de Chimie, UMR5182, Ecole Normale Supérieure de Lyon, Université de Lyon and CNRS, Lyon, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS)and CNRS, Toulouse, France
| |
Collapse
|
39
|
Uratani H, Chou CP, Nakai H. Quantum mechanical molecular dynamics simulations of polaron formation in methylammonium lead iodide perovskite. Phys Chem Chem Phys 2020; 22:97-106. [DOI: 10.1039/c9cp04739e] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Polaron formation in a halide perovskite is analyzed via nanometre-scale quantum mechanical molecular dynamics simulations.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry
- School of Advanced Science and Engineering
- Waseda University
- Shinjuku-ku
- Japan
| | - Chien-Pin Chou
- Waseda Research Institute for Science and Engineering (WISE)
- Tokyo 169-8555
- Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry
- School of Advanced Science and Engineering
- Waseda University
- Shinjuku-ku
- Japan
| |
Collapse
|
40
|
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
|
41
|
Finite-temperature-based linear-scaling divide-and-conquer self-consistent field method for static electron correlation systems. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.04.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|