1
|
Ding Y, Cheng M, Gao H. Imaging the Collision Energy-Dependent Charge-Transfer Dynamics between the Spin-Orbit Ground Ar +( 2P 3/2) Ion and CO. J Phys Chem Lett 2024; 15:6038-6044. [PMID: 38819134 DOI: 10.1021/acs.jpclett.4c01233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
The collisional charge-transfer reaction between Ar+(2P3/2,1/2) and CO represents one of the most studied ion-molecule systems; many controversies persist among different studies, and the detailed quantum state-to-state charge-transfer dynamics remains unknown. Here, differential cross sections of the charge-transfer process between the spin-orbit ground Ar+(2P3/2) ion and CO are reported at three center-of-mass collision energies of 1.02, 0.72, and 0.40 eV using a home-built three-dimensional velocity-map imaging-based ion-molecule crossed beam setup. At all three collision energies, the direct energy resonant charge-transfer mechanism dominates the reaction, featuring predominantly forward scattering with the CO+ product population peaking at the v' = 6 and v' = 7 vibrational levels. Only at the lowest collision energy of 0.40 eV is the significant backward peaked scattering product observed, with CO+ populated from v' = 4 to v' = 8. There is no obvious evidence for the formation of CO+ in excited electronic state A2Π+, in qualitative accord with previous theoretical predictions.
Collapse
Affiliation(s)
- Yufan Ding
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Min Cheng
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hong Gao
- Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
2
|
Vodeb J, Diego M, Vaskivskyi Y, Logaric L, Gerasimenko Y, Kabanov V, Lipovsek B, Topic M, Mihailovic D. Non-equilibrium quantum domain reconfiguration dynamics in a two-dimensional electronic crystal and a quantum annealer. Nat Commun 2024; 15:4836. [PMID: 38844460 PMCID: PMC11156939 DOI: 10.1038/s41467-024-49179-z] [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: 06/13/2023] [Accepted: 05/23/2024] [Indexed: 06/09/2024] Open
Abstract
Relaxation dynamics of complex many-body quantum systems trapped into metastable states is a very active field of research from both the theoretical and experimental point of view with implications in a wide array of topics from macroscopic quantum tunnelling and nucleosynthesis to non-equilibrium superconductivity and energy-efficient memory devices. In this work, we investigate quantum domain reconfiguration dynamics in the electronic superlattice of a quantum material using time-resolved scanning tunneling microscopy and unveil a crossover from temperature to noisy quantum fluctuation dominated dynamics. The process is modeled using a programmable superconducting quantum annealer in which qubit interconnections correspond directly to the microscopic interactions between electrons in the quantum material. Crucially, the dynamics of both the experiment and quantum simulation is driven by spectrally similar pink noise. We find that the simulations reproduce the emergent time evolution and temperature dependence of the experimentally observed electronic domain dynamics.
Collapse
Affiliation(s)
- Jaka Vodeb
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia.
- Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.
- Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52425, Jülich, Germany.
| | - Michele Diego
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia
| | - Yevhenii Vaskivskyi
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia
- Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia
| | - Leonard Logaric
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia
| | | | - Viktor Kabanov
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia
| | - Benjamin Lipovsek
- Faculty for Electrical Engineering, University of Ljubljana, Tržaška 25, 1000, Ljubljana, Slovenia
| | - Marko Topic
- Faculty for Electrical Engineering, University of Ljubljana, Tržaška 25, 1000, Ljubljana, Slovenia
| | - Dragan Mihailovic
- Jozef Stefan Institute, Jamova 39, 1000, Ljubljana, Slovenia.
- Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.
- CENN Nanocenter, Jamova 39, 1000, Ljubljana, Slovenia.
| |
Collapse
|
3
|
Pant R, Verma PK, Rangi C, Mondal E, Bhati M, Srinivasan V, Wüster S. Universal Measure for the Impact of Adiabaticity on Quantum Transitions. PHYSICAL REVIEW LETTERS 2024; 132:126903. [PMID: 38579224 DOI: 10.1103/physrevlett.132.126903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 01/25/2024] [Accepted: 02/26/2024] [Indexed: 04/07/2024]
Abstract
Adiabaticity is crucial for our understanding of complex quantum dynamics and thus for advancing fundamental physics and technology, but its impact cannot yet be quantified in complex but common cases where dynamics is only partially adiabatic, several eigenstates are simultaneously populated and transitions between noneigenstates are of key interest. We construct a universally applicable measure that can quantify the adiabaticity of quantum transitions in an arbitrary basis. Our measure distinguishes transitions that occur due to the adiabatic change of populated system eigenstates from transitions that occur due to beating between several eigenstates and can handle nonadiabatic events. While all quantum dynamics fall within the scope of the measure, we demonstrate its usage and utility through two important material science problems-energy and charge transfer-where adiabaticity could be effected by nuclear motion and its quantification will aid not only in unraveling mechanisms but also in system design, for example, of light harvesting systems.
Collapse
Affiliation(s)
- R Pant
- Department of Physics, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - P K Verma
- Department of Chemistry, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - C Rangi
- Department of Physics, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - E Mondal
- Department of Chemistry, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - M Bhati
- Department of Chemistry, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - V Srinivasan
- Department of Chemistry, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| | - S Wüster
- Department of Physics, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462 066, India
| |
Collapse
|
4
|
López X, Sánchez-Mansilla A, Sousa C, Straatsma TP, Broer R, de Graaf C. Comparison of Computational Strategies for the Calculation of the Electronic Coupling in Intermolecular Energy and Electron Transport Processes. J Phys Chem A 2023; 127:10717-10731. [PMID: 38084088 PMCID: PMC10749449 DOI: 10.1021/acs.jpca.3c05998] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 11/29/2023] [Accepted: 11/29/2023] [Indexed: 12/22/2023]
Abstract
Electronic couplings in intermolecular electron and energy transfer processes calculated by six different existing computational techniques are compared to nonorthogonal configuration interaction for fragments (NOCI-F) results. The paper addresses the calculation of the electronic coupling in diketopyrrolopyrol, tetracene, 5,5'-difluoroindigo, and benzene-Cl for hole and electron transport, as well as the local exciton and singlet fission coupling. NOCI-F provides a rigorous computational scheme to calculate these couplings, but its computational cost is rather elevated. The here-considered ab initio Frenkel-Davydov (AIFD), Dimer projection (DIPRO), transition dipole moment coupling, Michl-Smith, effective Hamiltonian, and Mulliken-Hush approaches are computationally less demanding, and the comparison with the NOCI-F results shows that the NOCI-F results in the couplings for hole and electron transport are rather accurately predicted by the more approximate schemes but that the NOCI-F exciton transfer and singlet fission couplings are more difficult to reproduce.
Collapse
Affiliation(s)
- Xavier López
- Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili, C. Marcel·lí Domingo 1, 43007 Tarragona, Spain
| | - Aitor Sánchez-Mansilla
- Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili, C. Marcel·lí Domingo 1, 43007 Tarragona, Spain
| | - Carmen Sousa
- Departament
de Ciència de Materials i Química Física and
Institut de Química Teòrica i Computacional, Universitat de Barcelona, C. Martí i Franquès, 08028 Barcelona, Spain
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | - Ria Broer
- Zernike
Institute of Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, Netherlands
| | - Coen de Graaf
- Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili, C. Marcel·lí Domingo 1, 43007 Tarragona, Spain
- ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
5
|
Skinner KC, Kammeraad JA, Wymore T, Narayan ARH, Zimmerman PM. Simulating Electron Transfer Reactions in Solution: Radical-Polar Crossover. J Phys Chem B 2023; 127:10097-10107. [PMID: 37976536 PMCID: PMC11135460 DOI: 10.1021/acs.jpcb.3c06120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Single-electron transfer (SET) promotes a wide variety of interesting chemical transformations, but modeling of SET requires a careful treatment of electronic and solvent effects to give meaningful insight. Therefore, a combined constrained density functional theory and molecular mechanics (CDFT/MM) tool is introduced specifically for SET-initiated reactions. Mechanisms for two radical-polar crossover reactions involving the organic electron donors tetrakis(dimethylamino)ethylene (TDAE) and tetrathiafulvalene (TTF) were studied with the new tool. An unexpected tertiary radical intermediate within the TDAE system was identified, relationships between kinetics and substitution in the TTF system were explained, and the impact of the solvent environments on the TDAE and TTF reactions were examined. The results highlight the need for including solvent dynamics when quantifying SET kinetics and thermodynamics, as a free energy difference of >20 kcal/mol was observed. Overall, the new method informs mechanistic analysis of SET-initiated reactions and therefore is envisioned to be useful for studying reactions in the condensed phase.
Collapse
Affiliation(s)
- Kevin C Skinner
- Life Sciences Institute, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Josh A Kammeraad
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Troy Wymore
- Laufer Center, Stony Brook University, Stony Brook, New York 11794, United States
| | - Alison R H Narayan
- Life Sciences Institute, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
6
|
Zeitter N, Hippchen N, Weidlich A, Jäger P, Ludwig P, Rominger F, Dreuw A, Freudenberg J, Bunz UHF. Hexakis-TIPS-Alkynylated Nonacenes: Persistent and Processible. Chemistry 2023; 29:e202302323. [PMID: 37490332 DOI: 10.1002/chem.202302323] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 07/25/2023] [Accepted: 07/25/2023] [Indexed: 07/26/2023]
Abstract
Four substituted nonacenes were prepared and characterized by UV-vis and EPR spectroscopy and X-ray crystallography. The compounds are the most stable and soluble nonacenes to date - due to six strategically placed triisopropylsilyl(TIPS)-ethynyl groups. They are stable for several weeks in the solid state. In dilute solution their half-life is 5-9 h. Crystal structure analyses of two nonacenes prove their structures. A nonacene derivative was tested in a solution-processed transistor and exhibits ambipolar charge transport (μe =0.007 cm2 /Vs; μh =0.023 cm2 /Vs).
Collapse
Affiliation(s)
- Nico Zeitter
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Nikolai Hippchen
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Anna Weidlich
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 205, 69120, Heidelberg, Germany
| | - Patrick Jäger
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Philipp Ludwig
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Frank Rominger
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 205, 69120, Heidelberg, Germany
| | - Jan Freudenberg
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Uwe H F Bunz
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| |
Collapse
|
7
|
Bazant MZ. Unified quantum theory of electrochemical kinetics by coupled ion-electron transfer. Faraday Discuss 2023; 246:60-124. [PMID: 37676178 DOI: 10.1039/d3fd00108c] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
A general theory of coupled ion-electron transfer (CIET) is presented, which unifies Marcus kinetics of electron transfer (ET) with Butler-Volmer kinetics of ion transfer (IT). In the limit of large reorganization energy, the theory predicts normal Marcus kinetics of "electron-coupled ion transfer" (ECIT). In the limit of large ion transfer energies, the theory predicts Butler-Volmer kinetics of "ion-coupled electron transfer" (ICET), where the charge transfer coefficient and exchange current are connected to microscopic properties of the electrode/electrolyte interface. In the ICET regime, the reductive and oxidative branches of Tafel's law are predicted to hold over a wide range of overpotentials, bounded by the ion-transfer energies for oxidation and reduction, respectively. The probability distribution of transferring electron energies in CIET smoothly interpolates between a shifted Gaussian distribution for ECIT (as in the Gerischer-Marcus theory of ET) to an asymmetric, fat-tailed Meixner distribution centered at the Fermi level for ICET. The latter may help interpret asymmetric line shapes in x-ray photo-electron spectroscopy (XPS) and Auger electron spectroscopy (AES) for metal surfaces in terms of shake-up relaxation of the ionized atom and its image polaron by ICET. In the limit of large overpotentials, the theory predicts a transition to inverted Marcus ECIT, leading to a universal reaction-limited current for metal electrodes, dominated by barrierless quantum transitions. Uniformly valid, closed-form asymptotic approximations are derived that smoothly transition between the limiting rate expressions for ICET and ECIT for metal electrodes, using simple but accurate mathematical functions. The theory is applied to lithium intercalation in lithium iron phosphate (LFP) and found to provide a consistent description of the observed current dependence on overpotential, temperature and concentration. CIET theory thus provides a critical bridge between quantum electrochemistry and electrochemical engineering, which may find many other applications and extensions.
Collapse
Affiliation(s)
- Martin Z Bazant
- Department of Chemical Engineering and Department of Mathematics, Massachusetts Institute of Technology, Cambridge 02139, MA, USA.
| |
Collapse
|
8
|
Mandal A, Taylor MA, Weight BM, Koessler ER, Li X, Huo P. Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics. Chem Rev 2023; 123:9786-9879. [PMID: 37552606 PMCID: PMC10450711 DOI: 10.1021/acs.chemrev.2c00855] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 08/10/2023]
Abstract
When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | - Michael A.D. Taylor
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| | - Braden M. Weight
- Department
of Physics and Astronomy, University of
Rochester, Rochester, New York 14627, United
States
| | - Eric R. Koessler
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Xinyang Li
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
9
|
Zhang X, Liu L, Li W, Wang C, Wang J, Fang WH, Chen X. Extended Single-Electron Transfer Model and Dynamically Associated Energy Transfer Event in a Dual-Functional Catalyst System. JACS AU 2023; 3:1452-1463. [PMID: 37234115 PMCID: PMC10206599 DOI: 10.1021/jacsau.3c00098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 04/20/2023] [Accepted: 04/21/2023] [Indexed: 05/27/2023]
Abstract
Organic photocatalysis has been developed flourishingly to rely on bimolecular energy transfer (EnT) or oxidative/reductive electron transfer (ET), promoting a variety of synthetic transformations. However, there are rare examples to merge EnT and ET processes rationally within one chemical system, of which the mechanistic investigation still remains in its infancy. Herein, the first mechanistic illustration and kinetic assessments of the dynamically associated EnT and ET paths were conducted for realizing the C-H functionalization in a cascade photochemical transformation of isomerization and cyclization by using the dual-functional organic photocatalyst of riboflavin. An extended single-electron transfer model of transition-state-coupled dual-nonadiabatic crossings was explored to analyze the dynamic behaviors in the proton transfer-coupled cyclization. This can also be used to clarify the dynamic correlation with the EnT-driven E → Z photoisomerization that has been kinetically evaluated by using Fermi's golden rule with the Dexter model. The present computational results of electron structures and kinetic data contribute to a fundamental basis for understanding the photocatalytic mechanism of the combined operation of EnT and ET strategies, which will guide the design and manipulation for the implementation of multiple activation modes based on a single photosensitizer.
Collapse
Affiliation(s)
- Xiaorui Zhang
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| | - Lin Liu
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| | - Weijia Li
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| | - Chu Wang
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| | - Juanjuan Wang
- College
of Nuclear Science and Technology, Beijing
Normal University, Xin-wai-da-jie
No. 19, Beijing 100875, China
- Laboratory
of Beam Technology and Energy Materials, Advanced Institute of Natural
Science, Beijing Normal University, Zhuhai 519087, China
| | - Wei-Hai Fang
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| | - Xuebo Chen
- Department
of Chemistry, Beijing Normal University, Xin-wai-da-jie No. 19, Beijing 100875, China
| |
Collapse
|
10
|
Zhang J, Tang Z, Zhang X, Zhu H, Zhao R, Lu Y, Gao J. Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States. J Chem Theory Comput 2023; 19:1777-1789. [PMID: 36917687 DOI: 10.1021/acs.jctc.2c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Zhen Tang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Xiaoyong Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Hong Zhu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China
| | - Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,Institute of Theoretical Chemistry, Jilin University, Changchun, 130023 Jilin, P. R. China
| | - Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
11
|
Kim TI, Lee IS, Kim H, Min SK. Calculation of exciton couplings based on density functional tight-binding coupled to state-interaction state-averaged ensemble-referenced Kohn-Sham approach. J Chem Phys 2023; 158:044106. [PMID: 36725518 DOI: 10.1063/5.0132361] [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/07/2023] Open
Abstract
We introduce the combination of the density functional tight binding (DFTB) approach, including onsite correction (OC) and long-range corrected (LC) functional and the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method with extended active space involving four electrons and four orbitals [LC-OC-DFTB/SSR(4,4)], to investigate exciton couplings in multichromophoric systems, such as organic crystals and molecular aggregates. We employ the LC-OC-DFTB/SSR(4,4) method to calculate the excitonic coupling in anthracene and tetracene. As a result, the LC-OC-DFTB/SSR(4,4) method provides a reliable description of the locally excited (LE) state in a single chromophore and the excitonic couplings between chromophores with reasonable accuracy compared to the experiment and the conventional SSR(4,4) method. In addition, the thermal fluctuation of excitonic couplings from dynamic nuclear motion in an anthracene crystal with LC-OC-DFTB/SSR(4,4) shows a similar fluctuation of excitonic coupling and spectral density with those of first-principle calculations. We conclude that LC-OC-DFTB/SSR(4,4) is capable of providing reasonable features related to LE states, such as Frenkel exciton with efficient computational cost.
Collapse
Affiliation(s)
- Tae In Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Hwon Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| |
Collapse
|
12
|
Niemeyer N, Eschenbach P, Bensberg M, Tölle J, Hellmann L, Lampe L, Massolle A, Rikus A, Schnieders D, Unsleber JP, Neugebauer J. The subsystem quantum chemistry program
Serenity. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1647] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Niklas Niemeyer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Moritz Bensberg
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Johannes Tölle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lars Hellmann
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lukas Lampe
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anja Massolle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anton Rikus
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - David Schnieders
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Jan P. Unsleber
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| |
Collapse
|
13
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
14
|
de Sousa LE, de Silva P. Diabatic Decomposition Perspective on the Role of Charge Transfer and Local Excitations in Thermally Activated Delayed Fluorescence. J Chem Theory Comput 2022; 18:5459-5470. [PMID: 36007256 DOI: 10.1021/acs.jctc.2c00310] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Thermally activated delayed fluorescence (TADF) is a phenomenon that relies on the upconversion of triplet excitons to singlet excitons by means of reverse intersystem crossing (rISC). It has been shown both experimentally and theoretically that the TADF mechanism depends on the interplay between charge transfer and local excitations. However, the difference between the diabatic and adiabatic character of the involved excited states is rarely discussed in the literature. Here we develop a diabatization procedure to implement a four-state model Hamiltonian to a set of TADF molecules. We provide physical interpretations of the Hamiltonian elements and show their dependence on the electronic state of the equilibrium geometry. We also demonstrate how vibrations affect the TADF efficiency by modifying the diabatic decomposition of the molecule. Finally, we provide a simple model that connects the diabatic Hamiltonian to the electronic properties relevant to TADF and show how this relationship translates into different optimization strategies for rISC, fluorescence, and overall TADF performance.
Collapse
Affiliation(s)
- Leonardo Evaristo de Sousa
- Department of Energy Conversion and Storage, Technical University of Denmark, Anker Engelunds Vej 301, 2800 Kongens Lyngby, Denmark
| | - Piotr de Silva
- Department of Energy Conversion and Storage, Technical University of Denmark, Anker Engelunds Vej 301, 2800 Kongens Lyngby, Denmark
| |
Collapse
|
15
|
Hong Y, Schlosser F, Kim W, Würthner F, Kim D. Ultrafast Symmetry-Breaking Charge Separation in a Perylene Bisimide Dimer Enabled by Vibronic Coupling and Breakdown of Adiabaticity. J Am Chem Soc 2022; 144:15539-15548. [PMID: 35951363 DOI: 10.1021/jacs.2c03916] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Perylene bisimides (PBIs) have received great attention in their applicability to optoelectronics. Especially, symmetry-breaking charge separation (SB-CS) in PBIs has been investigated to mimic the efficient light capturing and charge generation in natural light-harvesting systems. However, unlike ultrafast CS dynamics in donor-acceptor heterojunction materials, ultrafast SB-CS in a stacked homodimer has still been challenging due to excimer formation in the absence of rigidifying surroundings such as a special pair in the natural systems. Herein, we present the detailed mechanism of ultrafast photoinduced SB-CS occurring in a 1,7-bis(N-pyrrolidinyl) PBI dimer within a cyclophane. Through narrow-band and broad-band transient absorption spectroscopy, we demonstrate that ultrafast SB-CS in the dimer is enabled by the combination of (1) vibrationally coherent charge-transfer resonance-enhanced excimer formation and (2) breakdown of adiabaticity (formation of SB-CS diabats) in the excimer state via structural and solvent fluctuation. Quantum chemical calculations also underpin that the participation of strong electron-donating substituents in overall vibrational modes plays a crucial role in triggering the ultrafast SB-CS. Therefore, our work provides an alternative route to facilitate ultrafast SB-CS in PBIs and thereby establishes a novel strategy for the design of optoelectronic materials.
Collapse
Affiliation(s)
- Yongseok Hong
- Department of Chemistry, Yonsei University, Seoul 03722, Korea
| | - Felix Schlosser
- Institut für Organische Chemie & Center for Nanosystems Chemistry, Universitat Würzburg, Am Hubland, Würzburg 97074, Germany
| | - Woojae Kim
- Department of Chemistry, Yonsei University, Seoul 03722, Korea
| | - Frank Würthner
- Institut für Organische Chemie & Center for Nanosystems Chemistry, Universitat Würzburg, Am Hubland, Würzburg 97074, Germany
| | - Dongho Kim
- Department of Chemistry, Yonsei University, Seoul 03722, Korea.,Division of Energy Materials, Pohang University of Science and Technology (POSTECH), Pohang 37673, Republic of Korea
| |
Collapse
|
16
|
Ten Brink M, Gräber S, Hopjan M, Jansen D, Stolpp J, Heidrich-Meisner F, Blöchl PE. Real-time non-adiabatic dynamics in the one-dimensional Holstein model: Trajectory-based vs exact methods. J Chem Phys 2022; 156:234109. [PMID: 35732530 DOI: 10.1063/5.0092063] [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/14/2022] Open
Abstract
We benchmark a set of quantum-chemistry methods, including multitrajectory Ehrenfest, fewest-switches surface-hopping, and multiconfigurational-Ehrenfest dynamics, against exact quantum-many-body techniques by studying real-time dynamics in the Holstein model. This is a paradigmatic model in condensed matter theory incorporating a local coupling of electrons to Einstein phonons. For the two-site and three-site Holstein model, we discuss the exact and quantum-chemistry methods in terms of the Born-Huang formalism, covering different initial states, which either start on a single Born-Oppenheimer surface, or with the electron localized to a single site. For extended systems with up to 51 sites, we address both the physics of single Holstein polarons and the dynamics of charge-density waves at finite electron densities. For these extended systems, we compare the quantum-chemistry methods to exact dynamics obtained from time-dependent density matrix renormalization group calculations with local basis optimization (DMRG-LBO). We observe that the multitrajectory Ehrenfest method, in general, only captures the ultrashort time dynamics accurately. In contrast, the surface-hopping method with suitable corrections provides a much better description of the long-time behavior but struggles with the short-time description of coherences between different Born-Oppenheimer states. We show that the multiconfigurational Ehrenfest method yields a significant improvement over the multitrajectory Ehrenfest method and can be converged to the exact results in small systems with moderate computational efforts. We further observe that for extended systems, this convergence is slower with respect to the number of configurations. Our benchmark study demonstrates that DMRG-LBO is a useful tool for assessing the quality of the quantum-chemistry methods.
Collapse
Affiliation(s)
- M Ten Brink
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - S Gräber
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M Hopjan
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - D Jansen
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - J Stolpp
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - F Heidrich-Meisner
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - P E Blöchl
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
17
|
Axelrod S, Shakhnovich E, Gómez-Bombarelli R. Excited state non-adiabatic dynamics of large photoswitchable molecules using a chemically transferable machine learning potential. Nat Commun 2022; 13:3440. [PMID: 35705543 PMCID: PMC9200747 DOI: 10.1038/s41467-022-30999-w] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Accepted: 05/23/2022] [Indexed: 12/31/2022] Open
Abstract
Light-induced chemical processes are ubiquitous in nature and have widespread technological applications. For example, photoisomerization can allow a drug with a photo-switchable scaffold such as azobenzene to be activated with light. In principle, photoswitches with desired photophysical properties like high isomerization quantum yields can be identified through virtual screening with reactive simulations. In practice, these simulations are rarely used for screening, since they require hundreds of trajectories and expensive quantum chemical methods to account for non-adiabatic excited state effects. Here we introduce a diabatic artificial neural network (DANN), based on diabatic states, to accelerate such simulations for azobenzene derivatives. The network is six orders of magnitude faster than the quantum chemistry method used for training. DANN is transferable to azobenzene molecules outside the training set, predicting quantum yields for unseen species that are correlated with experiment. We use the model to virtually screen 3100 hypothetical molecules, and identify novel species with high predicted quantum yields. The model predictions are confirmed using high-accuracy non-adiabatic dynamics. Our results pave the way for fast and accurate virtual screening of photoactive compounds.
Collapse
Affiliation(s)
- Simon Axelrod
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, 02138, USA.,Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Eugene Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, 02138, USA
| | - Rafael Gómez-Bombarelli
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA.
| |
Collapse
|
18
|
Cederbaum LS. Cooperative molecular structure in polaritonic and dark states. J Chem Phys 2022; 156:184102. [DOI: 10.1063/5.0090047] [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/14/2022] Open
Abstract
The interaction of quantum light with matter is known to give rise to mixed light-matter states. An ensemble of identical molecules is discussed. The resulting hybrid light-matter states exhibit complex structure even if only a single vibrational coordinate per molecule is considered. Starting from the uniform situation where all molecules possess the same value of this coordinate, polaritons and dark states follow like in atoms, but are functions of this coordinate. It is proven that any point on a resulting polariton energy curve is a (local) minimum or maximum for distorting molecules perpendicular to this curve. It is shown how to explicitly compute the impact of distortion solely based on the data of a free molecule. The structure of the dark states and their behavior upon distortion is analyzed as well. Useful techniques are introduced and general results on, for example, minimum energy path, symmetry breaking and restoration, are obtained. The developed strategy is transferred to include several or even many nuclear degrees of freedom per molecule and it is demonstrated that the interplay of several vibrational degrees of freedom in a single molecule of the ensemble is expected to lead to qualitatively different physics. General consequences are discussed.
Collapse
|
19
|
Dakua KK, Rajak K, Mishra S. Spin–vibronic coupling in the quantum dynamics of a Fe(III) trigonal-bipyramidal complex. J Chem Phys 2022; 156:134103. [DOI: 10.1063/5.0080611] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The presence of a high density of excited electronic states in the immediate vicinity of the optically bright state of a molecule paves the way for numerous photo-relaxation channels. In transition-metal complexes, the presence of heavy atoms results in a stronger spin–orbit coupling, which enables spin forbidden spin-crossover processes to compete with the spin-allowed internal conversion processes. However, no matter how effectively the states cross around the Franck–Condon region, the degree of vibronic coupling, of both relativistic and non-relativistic nature, drives the population distribution among these states. One such case is demonstrated in this work for the intermediate-spin Fe(III) trigonal-bipyramidal complex. A quantum dynamical investigation of the photo-deactivation mechanism in the Fe(III) system is presented using the multi-configurational time-dependent Hartree approach based on the vibronic Hamiltonian whose coupling terms are derived from the state-averaged complete active space self-consistent field/complete active space with second-order perturbation theory (CASPT2) calculations and spin–orbit coupling of the scalar-relativistic CASPT2 states. The results of this study show that the presence of a strong (non-relativistic) vibronic coupling between the optically bright intermediate-spin state and other low-lying states of the same spin-multiplicity overpowers the spin–orbit coupling between the intermediate-spin and high-spin states, thereby lowering the chances of spin-crossover while exhibiting ultrafast relaxation among the intermediate-spin states. In a special case, where the population transfer pathway via the non-relativistic vibronic coupling is blocked, the probability of the spin-crossover is found to increase. This suggests that a careful modification of the complex by incorporation of heavier atoms with stronger relativistic effects can enhance the spin-crossover potential of Fe(III) intermediate-spin complexes.
Collapse
Affiliation(s)
- Kishan Kumar Dakua
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, India
| | - Karunamoy Rajak
- Centre for Theoretical Studies, Indian Institute of Technology Kharagpur, Kharagpur, India
| | - Sabyashachi Mishra
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, India
- Centre for Computational and Data Sciences, Indian Institute of Technology Kharagpur, Kharagpur, India
| |
Collapse
|
20
|
Nishio S, Kurashige Y. Importance of dynamical electron correlation in diabatic couplings of electron-exchange processes. J Chem Phys 2022; 156:114107. [PMID: 35317578 DOI: 10.1063/5.0075978] [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/14/2022] Open
Abstract
We demonstrate the importance of the dynamical electron correlation effect in diabatic couplings of electron-exchange processes in molecular aggregates. To perform a multireference perturbation theory with large active space of molecular aggregates, an efficient low-rank approximation is applied to the complete active space self-consistent field reference functions. It is known that kinetic rates of electron-exchange processes, such as singlet fission, triplet-triplet annihilation, and triplet exciton transfer, are not sufficiently explained by the direct term of the diabatic couplings but efficiently mediated by the low-lying charge transfer states if the two molecules are in close proximity. It is presented in this paper, however, that regardless of the distance of the molecules, the direct term is considerably underestimated by up to three orders of magnitude without the dynamical electron correlation, i.e., the diabatic states expressed in the active space are not adequate to quantitatively reproduce the electron-exchange processes.
Collapse
Affiliation(s)
- Soichiro Nishio
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
21
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
22
|
Tölle J, Neugebauer J. The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory. J Phys Chem Lett 2022; 13:1003-1018. [PMID: 35061387 DOI: 10.1021/acs.jpclett.1c04023] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
23
|
Cao K, Yin S, Wang Y. Insightful understanding of charge transfer processes in metalated phthalocyanines. Phys Chem Chem Phys 2022; 24:7635-7641. [DOI: 10.1039/d2cp00680d] [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
Marcus electron transfer theory coupling with quantum-mechanics (QM) calculations was applied to study the hole mobilities of a series of metalated phthalocyanine molecular crystals. The effect of metal on the...
Collapse
|
24
|
Carmona-Espíndola J, Gázquez JL. Charge transfer excitations and constrained density functional theory. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02860-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
25
|
Changala PB. Vibronic mean-field and perturbation theory for Jahn-Teller and pseudo-Jahn-Teller molecules. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1962556] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
26
|
Ko HY, Santra B, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles. J Chem Theory Comput 2021; 17:7789-7813. [PMID: 34775753 DOI: 10.1021/acs.jctc.0c01194] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H2O)128) and ≈6.8 s ((H2O)256) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
27
|
Ren M, Zhang L, Jiao Y, Chen Z, Wu W. Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer. J Chem Theory Comput 2021; 17:6861-6875. [PMID: 34605634 DOI: 10.1021/acs.jctc.1c00603] [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/29/2022]
Abstract
A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH2+ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.
Collapse
Affiliation(s)
- Mingxing Ren
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yang Jiao
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Zhenhua Chen
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
28
|
Bai X, Guo X, Wang L. Machine Learning Approach to Calculate Electronic Couplings between Quasi-diabatic Molecular Orbitals: The Case of DNA. J Phys Chem Lett 2021; 12:10457-10464. [PMID: 34672582 DOI: 10.1021/acs.jpclett.1c03053] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Diabatization of one-electron states in flexible molecular aggregates is a great challenge due to the presence of surface crossings between molecular orbital (MO) levels and the complex interaction between MOs of neighboring molecules. In this work, we present an efficient machine learning approach to calculate electronic couplings between quasi-diabatic MOs without the need of nonadiabatic coupling calculations. Using MOs of rigid molecules as references, the MOs that can be directly regarded to be quasi-diabatic in molecular dynamics are selected out, state tracked, and phase corrected. On the basis of this information, artificial neural networks are trained to characterize the structure-dependent onsite energies of quasi-diabatic MOs and the intermolecular electronic couplings. A representative sequence of DNA is systematically studied as an illustration. Smooth time evolution of electronic couplings in all base pairs is obtained with quasi-diabatic MOs. In particular, our method can calculate electronic couplings between different quasi-diabatic MOs independently, and thus, this possesses unique advantages in many applications.
Collapse
Affiliation(s)
- Xin Bai
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
29
|
Chen Z, Song J, Chen X, Zhou C, Wu W. N-Body Reduced Density Matrix-Based Valence Bond Theory and Its Applications in Diabatic Electronic-Structure Computations. Acc Chem Res 2021; 54:3895-3905. [PMID: 34615356 DOI: 10.1021/acs.accounts.1c00421] [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/30/2022]
Abstract
Valence bond (VB) theory, as a helpful complement to the more popular molecular orbital theory, is a fundamental electronic-structure theory that aims at interpreting molecular structure and chemical reactions in a lucid way. Both theoretical and experimental chemists have shown great interest in VB theory because of its capability of providing intuitive insight into the nature of chemical bonding and the mechanism of chemical reaction in a clear and comprehensible language rooted in Lewis structure. Therefore, there is a great call for the renaissance of VB theory. Nevertheless, this is possible only after a series of methods and algorithms were developed and efficiently implemented in user-friendly programs so as to serve computational chemists for general applications. In the past three decades, we have devoted a great amount of scientific enthusiasm toward this goal. In this Account, we will concisely summarize and briefly but insightfully discuss recent developments in ab initio VB theory, especially the N-body reduced density matrices (RDM)-based approach and its applications in diabatic electronic-structure computations, which is very useful for the vivid interpretation of many fundamental chemical processes such as electron and energy transfers. Furthermore, because of the fundamentally important role that the diabatic state plays in electron and energy transfers, which are two frontier research topics in both molecular and biochemical sciences, there are a broad range of applications that VB theory can handle.We start by briefly reviewing the general feature of ab initio VB wave functions. In particular, we focus on the multistructural ab initio VB theory that uses strictly localized orbitals, including the fundamental VB self-consistent field (VBSCF) and two post-SCF methods, VBCI and VBPT2, that use the VBSCF wave function as reference. We then allot a section to describing the recent developments of the RDM-based VB approach in the second quantization language. In this section, the enhanced Wick theorem is first outlined, followed by a brief discussion of its applications in evaluating VBSCF energy gradients and a Hessian with respect to the orbital expansion coefficients, together with a short review of the implementation of an automatic formula and code generator (AFCG) designed for many-body methods with nonorthogonal orbitals. Then, we introduce the application of the RDM-based approach in implementing the post-SCF method that addresses dynamic electronic correlation via perturbation theory, viz., the icVBPT2 method that adopts an internal contraction technique naturally. We finish this section by incorporating VB theory with the concept of seniority number, in which the tensor analysis technique is carefully exploited with the RDM-based approach, resulting in significant improvements in both the number of the active electrons/orbitals and in the speedup of the computational efficiency, thus pushing VB theory to its new limit. With these achievements available, we present the applications of VB theory in diabatic electronic-structure computations by using the intuitive insight rendered by VB theory. Therefore, we believe that there is a bright future in VB theory with true opportunities and new challenges coexisting both for theoretical developments and computational applications.
Collapse
Affiliation(s)
- Zhenhua Chen
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Jinshuai Song
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Xun Chen
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Chen Zhou
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
30
|
Dunnett AJ, Gowland D, Isborn CM, Chin AW, Zuehlsdorff TJ. Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue. J Chem Phys 2021; 155:144112. [PMID: 34654312 DOI: 10.1063/5.0062950] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.
Collapse
Affiliation(s)
- Angus J Dunnett
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
31
|
Hou S, Ji X, Gaskell K, Wang PF, Wang L, Xu J, Sun R, Borodin O, Wang C. Solvation sheath reorganization enables divalent metal batteries with fast interfacial charge transfer kinetics. Science 2021; 374:172-178. [PMID: 34618574 DOI: 10.1126/science.abg3954] [Citation(s) in RCA: 102] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
[Figure: see text].
Collapse
Affiliation(s)
- Singyuk Hou
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| | - Xiao Ji
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| | - Karen Gaskell
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA
| | - Peng-Fei Wang
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| | - Luning Wang
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA
| | - Jijian Xu
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| | - Ruimin Sun
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| | - Oleg Borodin
- Battery Science Branch, Sensors and Electron Devices Directorate, US Army Combat Capabilities Development Command Army Research Laboratory, Adelphi, MD 20783, USA
| | - Chunsheng Wang
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20740, USA
| |
Collapse
|
32
|
Rumble CA, Vauthey E. Molecular Dynamics Simulations of Bimolecular Electron Transfer: the Distance-Dependent Electronic Coupling. J Phys Chem B 2021; 125:10527-10537. [PMID: 34519508 DOI: 10.1021/acs.jpcb.1c05013] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Understanding the distance dependence of the parameters underpinning Marcus theory is imperative when interpreting the results of experiments on electron transfer (ET). Unfortunately, most of these parameters are difficult or impossible to access directly with experiments, necessitating the use of computer simulations to model them. In this work, we use molecular dynamics simulations in conjunction with constrained density functional theory calculations to study the distance dependence of the electronic coupling matrix element, |HRP|, for bimolecular ET. Contrary to what is typically assumed for such intermolecular reactions, we find that the magnitude of |HRP| does not decay exponentially with the center-of-mass separation of the reactants, rCOM. The addition of other simple measures of donor/acceptor (D/A) orientation did not improve the correlation of |HRP| with rCOM. Using the minimum distance separation, rmin, of the reactants as the structural descriptor allowed the system to be partitioned into high-coupling/close-contact and low-coupling/non-contact regimes, but large fluctuations of |HRP| were still found for the close-contact reactant pairs. Despite the persistent large fluctuations of |HRP|, its mean value was found to decay piecewise exponentially with increasing rmin, which was attributed to significant changes in the average D/A pair structure. The results herein advise one to use caution when interpreting the experimental results derived from spherical reactant models of bimolecular ET.
Collapse
Affiliation(s)
- Christopher A Rumble
- Départment de Chimie Physique, Université de Genève, 30 Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Eric Vauthey
- Départment de Chimie Physique, Université de Genève, 30 Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
33
|
Weight BM, Mandal A, Huo P. Ab initio symmetric quasi-classical approach to investigate molecular Tully models. J Chem Phys 2021; 155:084106. [PMID: 34470343 DOI: 10.1063/5.0061934] [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/14/2022] Open
Abstract
We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
34
|
Carmona-Espíndola J, Gázquez JL. Perturbation approach to constrained electron transfer in density functional theory. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02798-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
35
|
Lee Y, Kolesov G, Yao X, Kaxiras E, Cho K. Nonadiabatic dynamics of cobalt tricarbonyl nitrosyl for ligand dissociation induced by electronic excitation. Sci Rep 2021; 11:8997. [PMID: 33903609 PMCID: PMC8076320 DOI: 10.1038/s41598-021-88243-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 03/23/2021] [Indexed: 11/24/2022] Open
Abstract
We utilize real-time time-dependent density functional theory and Ehrenfest dynamics scheme to investigate excited-state nonadiabatic dynamics of ligand dissociation of cobalt tricarbonyl nitrosyl, Co(CO)3NO, which is a precursor used for cobalt growth in advanced technologies, where the precursor's reaction is enhanced by electronic excitation. Based on the first-principles calculations, we demonstrate two dissociation pathways of the NO ligand on the precursor. Detailed electronic structures are further analyzed to provide an insight into dynamics following the electronic excitations. This study sheds light on computational demonstration and underlying mechanism of the electronic-excitation-induced dissociation, especially in molecules with complex chemical bonds such as the Co(CO)3NO.
Collapse
Affiliation(s)
- Yeonghun Lee
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, TX, 75080, USA
| | - Grigory Kolesov
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Xiaolong Yao
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, TX, 75080, USA
| | - Efthimios Kaxiras
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Department of Physics, Harvard University, Cambridge, MA, 02138, USA
| | - Kyeongjae Cho
- Department of Materials Science and Engineering, University of Texas at Dallas, Richardson, TX, 75080, USA.
| |
Collapse
|
36
|
Pedron FN, Issoglio F, Estrin DA, Scherlis DA. Electron transfer pathways from quantum dynamics simulations. J Chem Phys 2021; 153:225102. [PMID: 33317287 DOI: 10.1063/5.0023577] [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/14/2022] Open
Abstract
This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, -1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant.
Collapse
Affiliation(s)
- F N Pedron
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Química Inorgánica, Analítica y Química Física, Buenos Aires, Argentina
| | - F Issoglio
- CONICET-Universidad de Buenos Aires, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), Buenos Aires, Argentina
| | - D A Estrin
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Química Inorgánica, Analítica y Química Física, Buenos Aires, Argentina
| | - D A Scherlis
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Química Inorgánica, Analítica y Química Física, Buenos Aires, Argentina
| |
Collapse
|
37
|
Camilo SRG, Curtolo F, Galassi VV, Arantes GM. Tunneling and Nonadiabatic Effects on a Proton-Coupled Electron Transfer Model for the Q o Site in Cytochrome bc1. J Chem Inf Model 2021; 61:1840-1849. [PMID: 33793213 DOI: 10.1021/acs.jcim.1c00008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Cytochrome bc1 is a fundamental enzyme for cellular respiration and photosynthesis. This dimeric protein complex catalyzes a proton-coupled electron transfer (PCET) from the reduced coenzyme-Q substrate (Q) to a bimetallic iron-sulfur cluster in the Qo active site. Herein, we combine molecular dynamics simulations of the complete cytochrome bc1 protein with electronic-structure calculations of truncated models and a semiclassical tunneling theory to investigate the electron-proton adiabaticity of the initial reaction catalyzed in the Qo site. After sampling possible orientations between the Q substrate and a histidine side chain that functions as hydrogen acceptor, we find that a truncated model composed by ubiquinol-methyl and imidazole-iron(III)-sulfide captures the expected changes in oxidation and spin states of the electron donor and acceptor. Diabatic electronic surfaces obtained for this model with multiconfigurational wave function calculations demonstrate that this reaction is electronic nonadiabatic, and proton tunneling is faster than mixing of electronic configurations. These results indicate the formalism that should be used to calculate vibronic couplings and kinetic parameters for the initial reaction in the Qo site of cytochrome bc1. This framework for molecular simulation may also be applied to investigate other PCET reactions in the Q-cycle or in various metalloproteins that catalyze proton translocation coupled to redox processes.
Collapse
Affiliation(s)
- Sofia R G Camilo
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| | - Felipe Curtolo
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| | - Vanesa V Galassi
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| | - Guilherme M Arantes
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| |
Collapse
|
38
|
Yin Z, Braams BJ, Fu B, Zhang DH. Neural Network Representation of Three-State Quasidiabatic Hamiltonians Based on the Transformation Properties from a Valence Bond Model: Three Singlet States of H3+. J Chem Theory Comput 2021; 17:1678-1690. [DOI: 10.1021/acs.jctc.0c01336] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zhengxi Yin
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Bastiaan J. Braams
- Centrum Wiskunde & Informatica (CWI), the Dutch National Center for Mathematics and Computer Science, 1098 XG Amsterdam, Netherlands
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| |
Collapse
|
39
|
Abstract
Recently, there has been much interest in the chirality-induced spin selectivity effect, whereby electron spin polarization, which is dependent on molecular chirality, is produced in electrode-molecule electron transfer processes. Naturally, one might consider if a similar effect can be observed in simple molecular charge transfer reactions, for example, in light-induced electron transfer from an electron donor to an electron acceptor. In this work, I explore the effect of electron transfer on spins in chiral single radicals and chiral radical pairs using Nakajima-Zwanzig theory. In these cases, chirality, in conjuction with spin-orbit coupling, does not lead to spin polarization, but instead, the electron transfer generates quantum coherence between spins states. In principle, this chirality-induced spin coherence could manifest in a range of experiments, and in particular, I demonstrate that the out of phase electron spin echo envelope modulation pulse electron paramagnetic resonance experiment would be able to detect this effect in oriented radical pairs.
Collapse
Affiliation(s)
- Thomas P Fay
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, U.K
| |
Collapse
|
40
|
Paz ASP, Glover WJ. Diabatic Many-Body Expansion: Development and Application to Charge-Transfer Reactions. J Chem Theory Comput 2021; 17:1497-1511. [DOI: 10.1021/acs.jctc.0c01231] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amiel S. P. Paz
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - William J. Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
41
|
Wang YC, Feng S, Liang W, Zhao Y. Electronic Couplings for Photoinduced Charge Transfer and Excitation Energy Transfer Based on Fragment Particle-Hole Densities. J Phys Chem Lett 2021; 12:1032-1039. [PMID: 33470827 DOI: 10.1021/acs.jpclett.0c03514] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
A new scheme is proposed to calculate the electronic couplings for photoinduced charge transfer and excitation energy transfer for both singlet and triplet states. In this scheme, the locally excited and charge-transfer states are constructed from the adiabatic ones by maximally localizing the particle (i.e., electron) and hole densities in terms of predefined molecular fragments. The construction process, after which the electronic couplings are directly obtained, is highly efficient and can be combined with various kinds of preliminary electronic structure calculations as long as the adiabatic excitation energies and transition densities are available. The method also applies to the systems with multiple charge or excitation centers. Its validity is demonstrated by the applications to the 6,13-dichloropentacene dimer and tetramer and the C60-Zn porphyrin dyad. The results reveal that the environment has a strong impact on the electronic couplings and can even enlarge those for long-range charge transfer by several orders of magnitude.
Collapse
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shishi Feng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
42
|
Abstract
Singlet fission (SF) is a photophysical downconversion pathway, in which a singlet excitation transforms into two triplet excited states. As such, it constitutes an exciton multiplication generation process, which is currently at the focal point for future integration into solar energy conversion devices. Beyond this, various other exciting applications were proposed, including quantum cryptography or organic light emitting diodes. Also, the mechanistic understanding evolved rapidly during the last year. Unfortunately, the number of suitable SF-chromophores is still limited. This is per se problematic, considering the wide range of envisaged applicability. With that in mind, we emphasize uncommon SF-scaffolds and outline requirements as well as strategies to expand the chromophore pool of SF-materials.
Collapse
Affiliation(s)
- Tobias Ullrich
- Friedrich-Alexander-Universität (FAU) Erlangen-Nürnberg, Department für Chemie und Pharmazie, Egerlandstr. 1-3, 91058 Erlangen, Germany.
| | | | | |
Collapse
|
43
|
Ramos P, Pavanello M. Nonadiabatic couplings from a variational excited state method based on constrained DFT. J Chem Phys 2021; 154:014110. [DOI: 10.1063/5.0028872] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
44
|
Barragan AM, Soudackov AV, Luthey-Schulten Z, Hammes-Schiffer S, Schulten K, Solov'yov IA. Theoretical Description of the Primary Proton-Coupled Electron Transfer Reaction in the Cytochrome bc1 Complex. J Am Chem Soc 2021; 143:715-723. [PMID: 33397104 DOI: 10.1021/jacs.0c07799] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The cytochrome bc1 complex is a transmembrane enzymatic protein complex that plays a central role in cellular energy production and is present in both photosynthetic and respiratory chain organelles. Its reaction mechanism is initiated by the binding of a quinol molecule to an active site, followed by a series of charge transfer reactions between the quinol and protein subunits. Previous work hypothesized that the primary reaction was a concerted proton-coupled electron transfer (PCET) reaction because of the apparent absence of intermediate states associated with single proton or electron transfer reactions. In the present study, the kinetics of the primary bc1 complex PCET reaction is investigated with a vibronically nonadiabatic PCET theory in conjunction with all-atom molecular dynamics simulations and electronic structure calculations. The computed rate constants and relatively high kinetic isotope effects are consistent with experimental measurements on related biomimetic systems. The analysis implicates a concerted PCET mechanism with significant hydrogen tunneling and nonadiabatic effects in the bc1 complex. Moreover, the employed theoretical framework is shown to serve as a general strategy for describing PCET reactions in bioenergetic systems.
Collapse
Affiliation(s)
- Angela M Barragan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Zaida Luthey-Schulten
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States.,Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Klaus Schulten
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Ilia A Solov'yov
- Department of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Strasse 9-11, Oldenburg 26111, Germany
| |
Collapse
|
45
|
Koch T, Bachmann J, Lettmann T, Doltsinis NL. Multiscale modelling of charge transport in P3HT:DIPBI bulk heterojunction organic solar cells. Phys Chem Chem Phys 2021; 23:12233-12250. [PMID: 34009221 DOI: 10.1039/d1cp00674f] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Charge transport properties of a P3HT:DIPBI bulk heterojunction solar cell are modelled by kinetic Monte Carlo simulations based on a morphology obtained from coarse-grained molecular dynamics. Different methods for calculating the hopping integrals entering the charge transfer rates are compared and calibrated for hole transport in amorphous P3HT. The influence of intermolecular and intramolecular charge transfer on the total charge carrier mobility and hence the power conversion efficiency is investigated in detail. An analysis of the most probable pathways with low resistance for hole transport is performed, establishing a connection between charge mobility and morphology.
Collapse
Affiliation(s)
- Tobias Koch
- Institut für Festkörpertheorie, Westfälische Wilhelms-Universität Münster and Center for Multiscale Theory & Computation, Wilhelm-Klemm-Str. 10, 48149 Münster, Germany.
| | - Jim Bachmann
- Institut für Festkörpertheorie, Westfälische Wilhelms-Universität Münster and Center for Multiscale Theory & Computation, Wilhelm-Klemm-Str. 10, 48149 Münster, Germany.
| | - Tobias Lettmann
- Institut für Festkörpertheorie, Westfälische Wilhelms-Universität Münster and Center for Multiscale Theory & Computation, Wilhelm-Klemm-Str. 10, 48149 Münster, Germany.
| | - Nikos L Doltsinis
- Institut für Festkörpertheorie, Westfälische Wilhelms-Universität Münster and Center for Multiscale Theory & Computation, Wilhelm-Klemm-Str. 10, 48149 Münster, Germany.
| |
Collapse
|
46
|
Bao P, Hettich CP, Shi Q, Gao J. Block-Localized Excitation for Excimer Complex and Diabatic Coupling. J Chem Theory Comput 2020; 17:240-254. [PMID: 33370101 DOI: 10.1021/acs.jctc.0c01015] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We describe a block-localized excitation (BLE) method to carry out constrained optimization of block-localized orbitals for constructing valence bond-like, diabatic excited configurations using multistate density functional theory (MSDFT). The method is an extension of the previous block-localized wave function method through a fragment-based ΔSCF approach to optimize excited determinants within a molecular complex. In BLE, both the number of electrons and the electronic spin of different fragments in a whole system can be constrained, whereas electrostatic, exchange, and polarization interactions among different blocks can be fully taken into account of. To avoid optimization collapse to unwanted states, a ΔSCF projection scheme and a maximum overlap of wave function approach have been presented. The method is illustrated by the excimer complex of two naphthalene molecules. With a minimum of eight spin-adapted configurational state functions, it was found that the inversion of La- and Lb- states near the optimal structure of the excimer complex is correctly produced, which is in quantitative agreement with DMRG-CASPT2 calculations and experiments. Trends in the computed transfer integrals associated with excited-state energy transfer both in the singlet and triplet states are discussed. The results suggest that MSDFT may be used as an efficient approach to treat intermolecular interactions in excited states with a minimal active space (MAS) for interpretation of the results and for dynamic simulations, although the selection of a small active space is often system dependent.
Collapse
Affiliation(s)
- Peng Bao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Christian P Hettich
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of the Chinese Academy of Sciences, Beijing 100049, China
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.,Beijing University Shenzhen Graduate School, Shenzhen 518055, China
| |
Collapse
|
47
|
Mao Y, Montoya-Castillo A, Markland TE. Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer. J Chem Phys 2020; 153:244111. [PMID: 33380087 DOI: 10.1063/5.0035593] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
48
|
Tölle J, Cupellini L, Mennucci B, Neugebauer J. Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization. J Chem Phys 2020; 153:184113. [PMID: 33187428 DOI: 10.1063/5.0022677] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| | - Lorenzo Cupellini
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
49
|
Shizu K, Adachi C, Kaji H. Visual Understanding of Vibronic Coupling and Quantitative Rate Expression for Singlet Fission in Molecular Aggregates. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2020. [DOI: 10.1246/bcsj.20200159] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Katsuyuki Shizu
- Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Chihaya Adachi
- Center for Organic Photonics and Electronics Research (OPERA), Kyushu University, 744 Motooka, Nishi, Fukuoka 819-0395, Japan
- International Institute for Carbon Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Motooka, Nishi, Fukuoka 819-0395, Japan
- Japan Science and Technology Agency (JST), ERATO, Adachi Molecular Exciton Engineering Project, 744 Motooka, Nishi, Fukuoka 819-0395, Japan
| | - Hironori Kaji
- Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
50
|
Ghan S, Kunkel C, Reuter K, Oberhofer H. Improved Projection-Operator Diabatization Schemes for the Calculation of Electronic Coupling Values. J Chem Theory Comput 2020; 16:7431-7443. [DOI: 10.1021/acs.jctc.0c00887] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Simiam Ghan
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Christian Kunkel
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
- Fritz-Haber Institute of the Max-Planck Society, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| |
Collapse
|