1
|
Bhat V, Ganapathysubramanian B, Risko C. Rapid Estimation of the Intermolecular Electronic Couplings and Charge-Carrier Mobilities of Crystalline Molecular Organic Semiconductors through a Machine Learning Pipeline. J Phys Chem Lett 2024; 15:7206-7213. [PMID: 38973725 DOI: 10.1021/acs.jpclett.4c01309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Organic semiconductors (OSC) offer tremendous potential across a wide range of (opto)electronic applications. OSC development, however, is often limited by trial-and-error design, with computational modeling approaches deployed to evaluate and screen candidates through a suite of molecular and materials descriptors that generally require hours to days of computational time to accumulate. Such bottlenecks slow the pace and limit the exploration of the vast chemical space comprising OSC. When considering charge-carrier transport in OSC, a key parameter of interest is the intermolecular electronic coupling. Here, we introduce a machine learning (ML) model to predict intermolecular electronic couplings in organic crystalline materials from their three-dimensional (3D) molecular geometries. The ML predictions take only a few seconds of computing time compared to hours by density functional theory (DFT) methods. To demonstrate the utility of the ML predictions, we deploy the ML model in conjunction with mathematical formulations to rapidly screen the charge-carrier mobility anisotropy for more than 60,000 molecular crystal structures and compare the ML predictions to DFT benchmarks.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| | - Baskar Ganapathysubramanian
- Department of Mechanical Engineering & Translational AI Center, Iowa State University, Ames, Iowa 50010, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
2
|
Sokolov M, Hoffmann DS, Dohmen PM, Krämer M, Höfener S, Kleinekathöfer U, Elstner M. Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex. Phys Chem Chem Phys 2024; 26:19469-19496. [PMID: 38979564 DOI: 10.1039/d4cp02116a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.
Collapse
Affiliation(s)
- Monja Sokolov
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - David S Hoffmann
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - Philipp M Dohmen
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Mila Krämer
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - Sebastian Höfener
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | | | - Marcus Elstner
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| |
Collapse
|
3
|
Kermarrec M, Dumont E, Gillet N. What tunes guanine ionization potential in a nucleosome? An all-in-one systematic QM/MM assessment. Biophys J 2024:S0006-3495(24)00448-X. [PMID: 38988071 DOI: 10.1016/j.bpj.2024.07.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2024] [Revised: 05/31/2024] [Accepted: 07/05/2024] [Indexed: 07/12/2024] Open
Abstract
Guanine radical cations are precursors to oxidatively induced DNA lesions, and the determination of oxidative DNA hot spots beyond oligonucleotides remains a current challenge. In order to rationalize the finetuned ionization properties of the ∼60 guanines in a nucleosome core particle, we report a robust molecular dynamics-then-FO-DFTB/MM (fragment-orbital tight-binding density functional theory/molecular mechanics) simulation protocol spanning 20 μs. Our work allows us to identify several factors governing guanine ionization potential and map oxidative hotspots. Our results highlight the predominant role of the proximity of positively charged histone residues in the modulation of the guanine ionization potential up to 0.6 eV. Consequently, fast, long-range hole transfer in nucleosomal DNA could be tuned by the proximity of histone tails, which differs, from a biological point of view, on the chromatin state.
Collapse
Affiliation(s)
- Maxime Kermarrec
- ENS de Lyon, CNRS, Université Claude Bernard Lyon 1, Laboratoire de Chimie UMR5182, Lyon, France
| | - Elise Dumont
- Université Côte d'Azur, CNRS, Institut de Chimie de Nice, UMR 7272, Nice, France; Institut Universitaire de France, Paris, France
| | - Natacha Gillet
- ENS de Lyon, CNRS, Université Claude Bernard Lyon 1, Laboratoire de Chimie UMR5182, Lyon, France.
| |
Collapse
|
4
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
5
|
Lin HH, Wang CI, Yang CH, Secario MK, Hsu CP. Two-Step Machine Learning Approach for Charge-Transfer Coupling with Structurally Diverse Data. J Phys Chem A 2024; 128:271-280. [PMID: 38157315 DOI: 10.1021/acs.jpca.3c04524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
Electronic coupling is important in determining charge-transfer rates and dynamics. Coupling strength is sensitive to both intermolecular, e.g., orientation or distance, and intramolecular degrees of freedom. Hence, it is challenging to build an accurate machine learning model to predict electronic coupling of molecular pairs, especially for those derived from the amorphous phase, for which intermolecular configurations are much more diverse than those derived from crystals. In this work, we devise a new prediction algorithm that employs two consecutive KRR models. The first model predicts molecular orbitals (MOs) from structural variation for each fragment, and coupling is further predicted by using the overlap integral included in a second model. With our two-step procedure, we achieved mean absolute errors of 0.27 meV for an ethylene dimer and 1.99 meV for a naphthalene pair, much improved accuracy amounting to 14-fold and 3-fold error reductions, respectively. In addition, MOs from the first model can also be the starting point to obtain other quantum chemical properties from atomistic structures. This approach is also compatible with a MO predictor with sufficient accuracy.
Collapse
Affiliation(s)
- Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Molecular Science and Digital Innovation Center, Genetics Generation Advancement Corp, No. 28, Ln. 36, Xinhu First Rd., Neihu, Taipei 114, Taiwan
| | - Chun-I Wang
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Chou-Hsun Yang
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
| | - Muhammad Khari Secario
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Taiwan International Graduate Program on Sustainable Chemical Science & Technology, Academia Sinica Institute of Chemistry, 128 Academia Road Sec.2, Nankang, Taipei 115, Taiwan
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 University Road, Hsinchu 300, Taiwan
| | - Chao-Ping Hsu
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Division of Physics, National Center for Theoretical Sciences, 1, Section 4, Roosevelt Road, Taipei 106, Taiwan
| |
Collapse
|
6
|
Salerno KM, Domenico J, Le NQ, Balakrishnan K, McQuillen RJ, Stiles CD, Solov'yov IA, Martino CF. Long-Time Oxygen and Superoxide Localization in Arabidopsis thaliana Cryptochrome. J Chem Inf Model 2023; 63:6756-6767. [PMID: 37874902 DOI: 10.1021/acs.jcim.3c00325] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
Cryptochromes are proteins that are highly conserved across species and in many instances bind the flavin adenine dinucleotide (FAD) cofactor within their photolyase-homology region (PHR) domain. The FAD cofactor has multiple redox states that help catalyze reactions, and absorbs photons at about 450 nm, a feature linked to the light-related functions of cryptochrome proteins. Reactive oxygen species (ROS) are produced from redox reactions involving molecular oxygen and are involved in a myriad of biological processes. Superoxide O2•- is an exemplary ROS that may be formed through electron transfer from FAD to O2, generating an electron radical pair. Although the formation of a superoxide-FAD radical pair has been speculated, it is still unclear if the required process steps could be realized in cryptochrome. Here, we present results from molecular dynamics (MD) simulations of oxygen interacting with the PHR domain of Arabidopsis thaliana cryptochrome 1 (AtCRY1). Using MD simulation trajectories, oxygen binding locations are characterized through both the O2-FAD intermolecular distance and the local protein environment. Oxygen unbinding times are characterized through replica simulations of the bound oxygen. Simulations reveal that oxygen molecules can localize at certain sites within the cryptochrome protein for tens of nanoseconds, and superoxide molecules can localize for significantly longer. This relatively long-duration molecule binding suggests the possibility of an electron-transfer reaction leading to superoxide formation. Estimates of electron-transfer rates using the Marcus theory are performed for the identified potential binding sites. Molecular oxygen binding results are compared with recent results demonstrating long-time oxygen binding within the electron-transfer flavoprotein (ETF), another FAD binding protein.
Collapse
Affiliation(s)
- K Michael Salerno
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Janna Domenico
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Nam Q Le
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Krithika Balakrishnan
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Ryan J McQuillen
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Christopher D Stiles
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| | - Ilia A Solov'yov
- Institute of Physics, Carl von Ossietzky University Oldenburg, Carl von Ossietzky Straße 9-11, 26129 Oldenburg, Germany
- Centre for Nanoscale Dynamics (CENAD), Carl von Ossietzky Universität Oldenburg, Carl von Ossietzky Straße 9-11, 26129 Oldenburg, Germany
- Centre for Neurosensory Science, Carl von Ossietzky Universität Oldenburg, Carl von Ossietzky Straße 9-11, 26129 Oldenburg, Germany
| | - Carlos F Martino
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, Maryland 20723, United States
| |
Collapse
|
7
|
Zahabi N, Baryshnikov G, Linares M, Zozoulenko I. Charge carrier dynamics in conducting polymer PEDOT using ab initio molecular dynamics simulations. J Chem Phys 2023; 159:154801. [PMID: 37843059 DOI: 10.1063/5.0169363] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 10/01/2023] [Indexed: 10/17/2023] Open
Abstract
As conducting polymers become increasingly important in electronic devices, understanding their charge transport is essential for material and device development. Various semi-empirical approaches have been used to describe temporal charge carrier dynamics in these materials, but there have yet to be any theoretical approaches utilizing ab initio molecular dynamics. In this work, we develop a computational technique based on ab initio Car-Parrinello molecular dynamics to trace charge carrier temporal motion in archetypical conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT). Particularly, we analyze charge dynamics in a single PEDOT chain and in two coupled chains with different degrees of coupling and study the effect of temperature. In our model we first initiate a positively charged polaron (compensated by a negative counterion) at one end of the chain, and subsequently displace the counterion to the other end of the chain and trace polaron dynamics in the system by monitoring bond length alternation in the PEDOT backbone and charge density distribution. We find that at low temperature (T = 1 K) the polaron distortion gradually disappears from its initial location and reappears near the new position of the counterion. At the room temperature (T = 300 K), we find that the distortions induced by polaron, and atomic vibrations are of the same magnitude, which makes tracking the polaron distortion challenging because it is hidden behind the temperature-induced vibrations. The novel approach developed in this work can be used to study polaron mobility along and between the chains, investigate charge transport in highly doped polymers, and explore other flexible polymers, including n-doped ones.
Collapse
Affiliation(s)
- Najmeh Zahabi
- Laboratory of Organic Electronics (LOE), Department of Science and Technology (ITN), Campus Norrköping, Linköping University, SE-60174 Norrköping, Sweden
| | - Glib Baryshnikov
- Laboratory of Organic Electronics (LOE), Department of Science and Technology (ITN), Campus Norrköping, Linköping University, SE-60174 Norrköping, Sweden
| | - Mathieu Linares
- Group of Scientific Visualization, Department of Science and Technology (ITN), Campus Norrköping, Linköping University, SE-60174 Norrköping, Sweden
- Swedish e-Science Center (SeRC), Linköping University, SE-581 83 Linköping, Sweden
| | - Igor Zozoulenko
- Laboratory of Organic Electronics (LOE), Department of Science and Technology (ITN), Campus Norrköping, Linköping University, SE-60174 Norrköping, Sweden
| |
Collapse
|
8
|
Kubař T, Elstner M, Cui Q. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annu Rev Biophys 2023; 52:525-551. [PMID: 36791746 PMCID: PMC10810093 DOI: 10.1146/annurev-biophys-111622-091140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.
Collapse
Affiliation(s)
- T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - Q Cui
- Department of Chemistry, Boston University, Boston, Massachusetts, USA;
- Department of Physics, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| |
Collapse
|
9
|
Tully JC. Ehrenfest dynamics with quantum mechanical nuclei. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
|
10
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, France
| |
Collapse
|
11
|
Pradhan CS, Jain A. Detailed Balance and Independent Electron Surface-Hopping Method: The Importance of Decoherence and Correct Calculation of Diabatic Populations. J Chem Theory Comput 2022; 18:4615-4626. [PMID: 35880817 DOI: 10.1021/acs.jctc.2c00320] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We benchmark and improve the independent electron surface-hopping (IESH) method developed by J. C. Tully's group for nonadiabatic simulations near metal surfaces. We have incorporated decoherence within the IESH method as well as implemented a scheme for the accurate calculation of diabatic populations. We benchmark the original IESH method with the above inclusions for a model system to calculate rate constants and long-time populations. The original IESH method fails to capture the detailed balance for some of the parameters, which is corrected with the inclusion of decoherence and accurate calculation of diabatic populations. Total rate constants are well captured both within the original IESH method as well as within our modified IESH.
Collapse
Affiliation(s)
- Chinmay S Pradhan
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
12
|
Giannini S, Blumberger J. Charge Transport in Organic Semiconductors: The Perspective from Nonadiabatic Molecular Dynamics. Acc Chem Res 2022; 55:819-830. [PMID: 35196456 PMCID: PMC8928466 DOI: 10.1021/acs.accounts.1c00675] [Citation(s) in RCA: 32] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
Organic semiconductors (OSs) are an exciting
class of materials
that have enabled disruptive technologies in this century including
large-area electronics, flexible displays, and inexpensive solar cells.
All of these technologies rely on the motion of electrical charges
within the material and the diffusivity of these charges critically
determines their performance. In this respect, it is remarkable that
the nature of the charge transport in these materials has puzzled
the community for so many years, even for apparently simple systems
such as molecular single crystals: some experiments would better fit
an interpretation in terms of a localized particle picture, akin to
molecular or biological electron transfer, while others are in better
agreement with a wave-like interpretation, more akin to band transport
in metals. Exciting recent progress in the theory and simulation
of charge
carrier transport in OSs has now led to a unified understanding of
these disparate findings, and this Account will review one of these
tools developed in our laboratory in some detail: direct charge carrier
propagation by quantum-classical nonadiabatic molecular dynamics.
One finds that even in defect-free crystals the charge carrier can
either localize on a single molecule or substantially delocalize over
a large number of molecules depending on the relative strength of
electronic couplings between the molecules, reorganization, or charge
trapping energy of the molecule and thermal fluctuations of electronic
couplings and site energies, also known as electron–phonon
couplings. Our simulations predict that in molecular OSs exhibiting
some of
the highest measured charge mobilities to date, the charge carrier
forms “flickering” polarons, objects that are delocalized
over 10–20 molecules on average and that constantly change
their shape and extension under the influence of thermal disorder.
The flickering polarons propagate through the OS by short (≈10
fs long) bursts of the wave function that lead to an expansion of
the polaron to about twice its size, resulting in spatial displacement,
carrier diffusion, charge mobility, and electrical conductivity. Arguably
best termed “transient delocalization”, this mechanistic
scenario is very similar to the one assumed in transient localization
theory and supports its assertions. We also review recent applications
of our methodology to charge transport in disordered and nanocrystalline
samples, which allows us to understand the influence of defects and
grain boundaries on the charge propagation. Unfortunately, the
energetically favorable packing structures of
typical OSs, whether molecular or polymeric, places fundamental constraints
on charge mobilities/electronic conductivity compared to inorganic
semiconductors, which limits their range of applications. In this
Account, we review the design rules that could pave the way for new
very high-mobility OS materials and we argue that 2D covalent organic
frameworks are one of the most promising candidates to satisfy them. We conclude that our nonadiabatic dynamics method is a powerful
approach for predicting charge carrier transport in crystalline and
disordered materials. We close with a brief outlook on extensions
of the method to exciton transport, dissociation, and recombination.
This will bring us a step closer to an understanding of the birth,
survival, and annihiliation of charges at interfaces of optoelectronic
devices.
Collapse
Affiliation(s)
- Samuele Giannini
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
13
|
Ziogos OG, Blumberger J. Ultrafast estimation of electronic couplings for electron transfer between pi-conjugated organic molecules. II. J Chem Phys 2021; 155:244110. [PMID: 34972358 DOI: 10.1063/5.0076555] [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
The development of highly efficient methods for the calculation of electronic coupling matrix elements between the electron donor and acceptor is an important goal in theoretical organic semiconductor research. In Paper I [F. Gajdos, S. Valner, F. Hoffmann, J. Spencer, M. Breuer, A. Kubas, M. Dupuis, and J. Blumberger, J. Chem. Theory Comput. 10, 4653 (2014)], we introduced the analytic overlap method (AOM) for this purpose, which is an ultrafast electronic coupling estimator parameterized to and orders of magnitude faster than density functional theory (DFT) calculations at a reasonably small loss in accuracy. In this work, we reparameterize and extend the AOM to molecules containing nitrogen, oxygen, fluorine, and sulfur heteroatoms using 921 dimer configurations from the recently introduced HAB79 dataset. We find again a very good linear correlation between the frontier orbital overlap, calculated ultrafast in an optimized minimum Slater basis, and DFT reference electronic couplings. The new parameterization scheme is shown to be transferable to sulfur-containing polyaromatic hydrocarbons in experimentally resolved dimeric configurations. Our extension of the AOM enables high-throughput screening of very large databases of chemically diverse organic crystal structures and the application of computationally intense non-adiabatic molecular dynamics methods to charge transport in state-of-the-art organic semiconductors, e.g., non-fullerene acceptors.
Collapse
Affiliation(s)
- Orestis George Ziogos
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
14
|
Aggarwal A, Sahoo AK, Bag S, Kaliginedi V, Jain M, Maiti PK. Fine-tuning the DNA conductance by intercalation of drug molecules. Phys Rev E 2021; 103:032411. [PMID: 33862831 DOI: 10.1103/physreve.103.032411] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 03/08/2021] [Indexed: 11/07/2022]
Abstract
In this work we study the structure-transport property relationships of small ligand intercalated DNA molecules using a multiscale modeling approach where extensive ab initio calculations are performed on numerous MD-simulated configurations of dsDNA and dsDNA intercalated with two different intercalators, ethidium and daunomycin. DNA conductance is found to increase by one order of magnitude upon drug intercalation due to the local unwinding of the DNA base pairs adjacent to the intercalated sites, which leads to modifications of the density of states in the near-Fermi-energy region of the ligand-DNA complex. Our study suggests that the intercalators can be used to enhance or tune the DNA conductance, which opens new possibilities for their potential applications in nanoelectronics.
Collapse
Affiliation(s)
- Abhishek Aggarwal
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Anil Kumar Sahoo
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Saientan Bag
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Veerabhadrarao Kaliginedi
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| | - Manish Jain
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Prabal K Maiti
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
15
|
Yang ZH. Speed-dependent adaptive partitioning in QM/MM MD simulations of displacement damage in solid-state systems. Phys Chem Chem Phys 2021; 23:3417-3428. [PMID: 33506833 DOI: 10.1039/d0cp05149g] [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/21/2022]
Abstract
Solids undergo displacement damage (DD) when interacting with energetic particles, which may happen during the fabrication of semiconductor devices, in harsh environments and in certain analysis techniques. Simulations of DD generation are usually carried out using classical molecular dynamics (MD), but classical MD does not account for all the effects in DD, as demonstrated by ab initio calculations of model systems in the literature. A complete ab initio simulation of DD generation is impractical due to the large number of atoms involved. In my previous paper [Yang, Phys. Chem. Chem. Phys., 2020, 22, 19307], I developed an adaptive-center (AC) method for the adaptive-partitioning (AP) of quantum mechanics/molecular mechanics (QM/MM) simulations, allowing the active region centers and the QM/MM partition to be determined on-the-fly for energy-conserving AP-QM/MM methods. I demonstrated that the AC-AP-QM/MM is applicable to the simulation of DD generation, so that the active regions can be treated using an ab initio method. The AC method could not be used to identify the fast-moving recoil ions in DD generation as active region centers, however, and the accuracy is negatively affected by the rapid change in the QM/MM partition of the system. In this paper, I extend the AC method and develop a speed-dependent adaptive-center (SDAC) method for accurate AP-QM/MM simulations of DD. The SDAC method is applicable to general problems with speed-dependent active regions, and is compatible with all existing energy-conserving partitioning-by-distance AP-QM/MM methods. The artifact due to the speed-dependent potential energy surface can be made small by choosing suitable criteria. I demonstrate the SDAC method by simulations of DD generation in bulk silicon.
Collapse
Affiliation(s)
- Zeng-Hui Yang
- Microsystem and Terahertz Research Center, China Academy of Engineering Physics, Chengdu, 610200, China.
| |
Collapse
|
16
|
Aggarwal A, Vinayak V, Bag S, Bhattacharyya C, Waghmare UV, Maiti PK. Predicting the DNA Conductance Using a Deep Feedforward Neural Network Model. J Chem Inf Model 2020; 61:106-114. [PMID: 33320660 DOI: 10.1021/acs.jcim.0c01072] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Double-stranded DNA (dsDNA) has been established as an efficient medium for charge migration, bringing it to the forefront of the field of molecular electronics and biological research. The charge migration rate is controlled by the electronic couplings between the two nucleobases of DNA/RNA. These electronic couplings strongly depend on the intermolecular geometry and orientation. Estimating these electronic couplings for all the possible relative geometries of molecules using the computationally demanding first-principles calculations requires a lot of time and computational resources. In this article, we present a machine learning (ML)-based model to calculate the electronic coupling between any two bases of dsDNA/dsRNA and bypass the computationally expensive first-principles calculations. Using the Coulomb matrix representation which encodes the atomic identities and coordinates of the DNA base pairs to prepare the input dataset, we train a feedforward neural network model. Our neural network (NN) model can predict the electronic couplings between dsDNA base pairs with any structural orientation with a mean absolute error (MAE) of less than 0.014 eV. We further use the NN-predicted electronic coupling values to compute the dsDNA/dsRNA conductance.
Collapse
Affiliation(s)
- Abhishek Aggarwal
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Vinayak Vinayak
- Undergraduate Program, Indian Institute of Science, Bangalore 560012, India
| | - Saientan Bag
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Chiranjib Bhattacharyya
- Department of Computer Science and Automation, Indian Institute of Science, Bangalore 560012, India
| | - Umesh V Waghmare
- Theoretical Sciences Unit, Jawaharlal Nehru Center for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| | - Prabal K Maiti
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
17
|
Kröncke S, Herrmann C. Toward a First-Principles Evaluation of Transport Mechanisms in Molecular Wires. J Chem Theory Comput 2020; 16:6267-6279. [PMID: 32886502 DOI: 10.1021/acs.jctc.0c00667] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding charge transport through molecular wires is important for nanoscale electronics and biochemistry. Our goal is to establish a simple first-principles protocol for predicting the charge transport mechanism in such wires, in particular the crossover from coherent tunneling for short wires to incoherent hopping for longer wires. This protocol is based on a combination of density functional theory with a polarizable continuum model introduced by Kaupp et al. for mixed-valence molecules, which we had previously found to work well for length-dependent charge delocalization in such systems. We combine this protocol with a new charge delocalization measure tailored for molecular wires, and we show that it can predict the tunneling-to-hopping transition length with a maximum error of one subunit in five sets of molecular wires studied experimentally in molecular junctions at room temperature. This suggests that the protocol is also well suited for estimating the extent of hopping sites as relevant, for example, for the intermediate tunneling-hopping regime in DNA.
Collapse
Affiliation(s)
- Susanne Kröncke
- Department of Chemistry, University of Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| | - Carmen Herrmann
- Department of Chemistry, University of Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
| |
Collapse
|
18
|
Shuai Z, Li W, Ren J, Jiang Y, Geng H. Applying Marcus theory to describe the carrier transports in organic semiconductors: Limitations and beyond. J Chem Phys 2020; 153:080902. [PMID: 32872875 DOI: 10.1063/5.0018312] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.
Collapse
Affiliation(s)
- Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Yuqian Jiang
- Laboratory for Nanosystem and Hierarchy Fabrication, National Center for Nanoscience and Technology, Chinese Academy of Sciences, 100084 Beijing, People's Republic of China
| | - Hua Geng
- Department of Chemistry, Capital Normal University, 100048 Beijing, People's Republic of China
| |
Collapse
|
19
|
Nematiaram T, Troisi A. Modeling charge transport in high-mobility molecular semiconductors: Balancing electronic structure and quantum dynamics methods with the help of experiments. J Chem Phys 2020; 152:190902. [DOI: 10.1063/5.0008357] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Tahereh Nematiaram
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Alessandro Troisi
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| |
Collapse
|
20
|
Plasser F, Mai S, Fumanal M, Gindensperger E, Daniel C, González L. Strong Influence of Decoherence Corrections and Momentum Rescaling in Surface Hopping Dynamics of Transition Metal Complexes. J Chem Theory Comput 2019; 15:5031-5045. [DOI: 10.1021/acs.jctc.9b00525] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Felix Plasser
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K
| | - Sebastian Mai
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| | - Maria Fumanal
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Etienne Gindensperger
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Chantal Daniel
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Leticia González
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|
21
|
Holub D, Lamparter T, Elstner M, Gillet N. Biological relevance of charge transfer branching pathways in photolyases. Phys Chem Chem Phys 2019; 21:17072-17081. [PMID: 31313765 DOI: 10.1039/c9cp01609k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The repair of sun-induced DNA lesions by photolyases is driven by a photoinduced electron transfer from a fully reduced FAD to the damaged DNA. A chain of several aromatic residues connecting FAD to solvent ensures the prior photoreduction of the FAD cofactor. In PhrA, a class III CPD photolyase, two branching tryptophan charge transfer pathways have been characterized. According to previous experiments, both pathways play a role in the FAD photoreduction. To provide a molecular insight to the charge transfer abilities of both pathways, we perform multiscales simulations where the protein motion and the positive charge are simultaneously propagated. Our computational approach reveals that one pathway drives a very fast charge transfer whereas the other pathway provides a very good thermodynamic stabilization of the positive charge. During the simulations, the positive charge firstly moves on the fast triad, while a reorganization of the close FAD˙- environment occurs. Then, backward transfers can lead to the propagation of the positive charge on the second pathway. After one nanosecond, we observe a nearly equal probability to find the charge at ending tryptophan of either pathway; eventually the charge distribution will likely evolve towards a charge stabilization on the last tryptophan of the slowest pathway. Our results highlight the role the protein environment, which manages the association of a kinetic and a thermodynamic pathways to trigger a fast and efficient FAD photoreduction.
Collapse
Affiliation(s)
- Daniel Holub
- Department for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute for Technology, Kaiserstr. 12, 76131, Karlsruhe, Germany.
| | - Tilman Lamparter
- Botanical Institute, Karlsruhe Institute of Technology, Fritz Haber Weg 4, 76131, Karlsruhe, Germany
| | - Marcus Elstner
- Department for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute for Technology, Kaiserstr. 12, 76131, Karlsruhe, Germany. and Institute of Biological Interfaces (IBG2), Karlsruhe Institute for Technology, Kaiserstr. 12, 76131, Karlsruhe, Germany
| | - Natacha Gillet
- Department for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute for Technology, Kaiserstr. 12, 76131, Karlsruhe, Germany.
| |
Collapse
|
22
|
Sjulstok E, Lüdemann G, Kubař T, Elstner M, Solov'yov IA. Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome. Biophys J 2019; 114:2563-2572. [PMID: 29874607 DOI: 10.1016/j.bpj.2018.04.014] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2017] [Revised: 04/04/2018] [Accepted: 04/05/2018] [Indexed: 01/23/2023] Open
Abstract
Cryptochrome proteins are activated by the absorption of blue light, leading to the formation of radical pairs through electron transfer in the active site. Recent experimental studies have shown that once some of the amino acid residues in the active site of Xenopus laevis cryptochrome DASH are mutated, radical-pair formation is still observed. In this study, we computationally investigate electron-transfer pathways in the X. laevis cryptochrome DASH by extensively equilibrating a previously established homology model using molecular dynamics simulations and then mutating key amino acids involved in the electron transfer. The electron-transfer pathways are then probed by using tight-binding density-functional theory. We report the alternative electron-transfer pathways resolved at the molecular level and, through comparison of amino acid sequences for cryptochromes from different species, we demonstrate that one of these alternative electron-transfer pathways could be general for all cryptochrome DASH proteins.
Collapse
Affiliation(s)
- Emil Sjulstok
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Gesa Lüdemann
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Ilia A Solov'yov
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark, Odense, Denmark.
| |
Collapse
|
23
|
Peng J, Xie Y, Hu D, Lan Z. Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics. J Chem Phys 2019; 150:164126. [PMID: 31042919 DOI: 10.1063/1.5079426] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.
Collapse
Affiliation(s)
- Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
24
|
Carof A, Giannini S, Blumberger J. Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping. J Chem Phys 2018; 147:214113. [PMID: 29221382 DOI: 10.1063/1.5003820] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.
Collapse
Affiliation(s)
- Antoine Carof
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Samuele Giannini
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
25
|
Giannini S, Carof A, Blumberger J. Crossover from Hopping to Band-Like Charge Transport in an Organic Semiconductor Model: Atomistic Nonadiabatic Molecular Dynamics Simulation. J Phys Chem Lett 2018; 9:3116-3123. [PMID: 29787275 PMCID: PMC6077769 DOI: 10.1021/acs.jpclett.8b01112] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Accepted: 05/22/2018] [Indexed: 05/25/2023]
Abstract
The mechanism of charge transport (CT) in a 1D atomistic model of an organic semiconductor is investigated using surface hopping nonadiabatic molecular dynamics. The simulations benefit from a newly implemented state tracking algorithm that accounts for trivial surface crossings and from a projection algorithm that removes decoherence correction-induced artificial long-range charge transfer. The CT mechanism changes from slow hopping of a fully localized charge to fast diffusion of a polaron delocalized over several molecules as electronic coupling between the molecules exceeds the critical threshold V ≥ λ/2 (λ is the reorganization energy). With increasing temperature, the polaron becomes more localized and the mobility exhibits a "band-like" power law decay due to increased site energy and electronic coupling fluctuations (local and nonlocal electron-phonon coupling). Thus, reducing both types of electron-phonon coupling while retaining high mean electronic couplings should be part of the strategy toward discovery of new organics with high room-temperature mobilities.
Collapse
Affiliation(s)
- Samuele Giannini
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Antoine Carof
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
- Institute
for Advanced Study, Technische Universität
München, Lichtenbergstrasse
2 a, D-85748 Garching, Germany
| |
Collapse
|
26
|
Meng Z, Kubar T, Mu Y, Shao F. A Molecular Dynamics-Quantum Mechanics Theoretical Study of DNA-Mediated Charge Transport in Hydrated Ionic Liquids. J Chem Theory Comput 2018; 14:2733-2742. [PMID: 29570288 DOI: 10.1021/acs.jctc.7b01201] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Charge transport (CT) through biomolecules is of high significance in the research fields of biology, nanotechnology, and molecular devices. Inspired by our previous work that showed the binding of ionic liquid (IL) facilitated charge transport in duplex DNA, in silico simulation is a useful means to understand the microscopic mechanism of the facilitation phenomenon. Here molecular dynamics simulations (MD) of duplex DNA in water and hydrated ionic liquids were employed to explore the helical parameters. Principal component analysis was further applied to capture the subtle conformational changes of helical DNA upon different environmental impacts. Sequentially, CT rates were calculated by a QM/MM simulation of the flickering resonance model based upon MD trajectories. Herein, MD simulation illustrated that the binding of ionic liquids can restrain dynamic conformation and lower the on-site energy of the DNA base. Confined movement among the adjacent base pairs was highly related to the increase of electronic coupling among base pairs, which may lead DNA to a CT facilitated state. Sequentially combining MD and QM/MM analysis, the rational correlations among the binding modes, the conformational changes, and CT rates illustrated the facilitation effects from hydrated IL on DNA CT and supported a conformational-gating mechanism.
Collapse
Affiliation(s)
- Zhenyu Meng
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences , Nanyang Technological University , 21 Nanyang Link , Singapore 637371
| | - Tomas Kubar
- Institute of Physical Chemistry & Center for Functional Nanostructures , Karlsruhe Institute of Technology , 76131 Karlsruhe , Germany
| | - Yuguang Mu
- School of Biological Sciences , Nanyang Technological University , 60 Nanyang Drive , Singapore 637551
| | - Fangwei Shao
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences , Nanyang Technological University , 21 Nanyang Link , Singapore 637371
| |
Collapse
|
27
|
Holub D, Ma H, Krauß N, Lamparter T, Elstner M, Gillet N. Functional role of an unusual tyrosine residue in the electron transfer chain of a prokaryotic (6-4) photolyase. Chem Sci 2017; 9:1259-1272. [PMID: 29675172 PMCID: PMC5887102 DOI: 10.1039/c7sc03386a] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2017] [Accepted: 12/09/2017] [Indexed: 11/21/2022] Open
Abstract
Cryptochromes and photolyases form a flavoprotein family in which the FAD chromophore undergoes light induced changes of its redox state. During this process, termed photoreduction, electrons flow from the surface via conserved amino acid residues to FAD. The bacterial (6-4) photolyase PhrB belongs to a phylogenetically ancient group. Photoreduction of PhrB differs from the typical pattern because the amino acid of the electron cascade next to FAD is a tyrosine (Tyr391), whereas photolyases and cryptochromes of other groups have a tryptophan as direct electron donor of FAD. Mutagenesis studies have identified Trp342 and Trp390 as essential for charge transfer. Trp342 is located at the periphery of PhrB while Trp390 connects Trp342 and Tyr391. The role of Tyr391, which lies between Trp390 and FAD, is however unclear as its replacement by phenylalanine did not block photoreduction. Experiments reported here, which replace Tyr391 by Ala, show that photoreduction is blocked, underlining the relevance of Tyr/Phe at position 391 and indicating that charge transfer occurs via the triad 391-390-342. This raises the question, why PhrB positions a tyrosine at this location, having a less favourable ionisation potential than tryptophan, which occurs at this position in many proteins of the photolyase/cryptochrome family. Tunnelling matrix calculations show that tyrosine or phenylalanine can be involved in a productive bridged electron transfer between FAD and Trp390, in line with experimental findings. Since replacement of Tyr391 by Trp resulted in loss of FAD and DMRL chromophores, electron transfer cannot be studied experimentally in this mutant, but calculations on a mutant model suggest that Trp might participate in the electron transfer cascade. Charge transfer simulations reveal an unusual stabilization of the positive charge on site 391 compared to other photolyases or cryptochromes. Water molecules near Tyr391 offer a polar environment which stabilizes the positive charge on this site, thereby lowering the energetic barrier intrinsic to tyrosine. This opens a second charge transfer channel in addition to tunnelling through the tyrosine barrier, based on hopping and therefore transient oxidation of Tyr391, which enables a fast charge transfer similar to proteins utilizing a tryptophan-triad. Our results suggest that evolution of the first site of the redox chain has just been possible by tuning the protein structure and environment to manage a downhill hole transfer process from FAD to solvent.
Collapse
Affiliation(s)
- Daniel Holub
- Department for Theoretical Chemical Biology , Institute for Physical Chemistry , Karlsruhe Institute for Technology , Kaiserstr. 12 , 76131 , Karlsruhe , Germany .
| | - Hongju Ma
- Botanical Institute , Karlsruhe Institute for Technology , Fritz Haber Weg 4 , 76131 , Karlsruhe , Germany
| | - Norbert Krauß
- Botanical Institute , Karlsruhe Institute for Technology , Fritz Haber Weg 4 , 76131 , Karlsruhe , Germany
| | - Tilman Lamparter
- Botanical Institute , Karlsruhe Institute for Technology , Fritz Haber Weg 4 , 76131 , Karlsruhe , Germany
| | - Marcus Elstner
- Department for Theoretical Chemical Biology , Institute for Physical Chemistry , Karlsruhe Institute for Technology , Kaiserstr. 12 , 76131 , Karlsruhe , Germany . .,Institute of Biological Interfaces (IGB2) , Karlsruhe Institute for Technology , Kaiserstr. 12 , 76131 , Karlsruhe , Germany
| | - Natacha Gillet
- Department for Theoretical Chemical Biology , Institute for Physical Chemistry , Karlsruhe Institute for Technology , Kaiserstr. 12 , 76131 , Karlsruhe , Germany .
| |
Collapse
|
28
|
Voityuk AA. Electronic Couplings for Photoinduced Electron Transfer and Excitation Energy Transfer Computed Using Excited States of Noninteracting Molecules. J Phys Chem A 2017; 121:5414-5419. [DOI: 10.1021/acs.jpca.7b03924] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexander A. Voityuk
- Institució Catalana de Recerca i Estudis Avançats, 08010 Barcelona, Spain
- Institut de Química Computacional i Catàlisi (IQCC), Universitat de Girona, 17071 Girona, Spain
| |
Collapse
|
29
|
Oberhofer H, Reuter K, Blumberger J. Charge Transport in Molecular Materials: An Assessment of Computational Methods. Chem Rev 2017. [PMID: 28644623 DOI: 10.1021/acs.chemrev.7b00086] [Citation(s) in RCA: 191] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.
Collapse
Affiliation(s)
- Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom.,Institute for Advanced Study, Technische Universität München , Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
30
|
Nakano H, Sato H. Introducing the mean field approximation to CDFT/MMpol method: Statistically converged equilibrium and nonequilibrium free energy calculation for electron transfer reactions in condensed phases. J Chem Phys 2017; 146:154101. [DOI: 10.1063/1.4979895] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Hiroshi Nakano
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| |
Collapse
|
31
|
Wolter M, Elstner M, Kleinekathöfer U, Kubař T. Microsecond Simulation of Electron Transfer in DNA: Bottom-Up Parametrization of an Efficient Electron Transfer Model Based on Atomistic Details. J Phys Chem B 2017; 121:529-549. [PMID: 28045546 DOI: 10.1021/acs.jpcb.6b11384] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The transfer of electrons over long distances in complex molecular systems is a phenomenon of significance in both biochemistry and technology. In recent years, we have been developing efficient models to study ET in complex systems, including DNA as a prominent example. Ab initio and model approaches have been combined in an "on-the-fly" calculation of ET parameters, which can be used to propagate nuclear and electronic degrees of freedom simultaneously. These previous efforts have aimed at deriving an efficient nonadiabatic quantum mechanical-molecular mechanical (QM/MM) simulation scheme for ET, making nanosecond simulations of ET in realistic systems possible. This, however, is still insufficient for the treatment of large donor-bridge-acceptor systems, like the ET in DNA, overcoming long adenine bridges. Therefore, we have constructed a theoretical model in a bottom-up manner. All quantum-chemical as well as force-field calculations are substituted by theoretical models of the involved phenomena on a molecular level, including polarization and relaxation of the molecular environment, which are often omitted in other recently developed theoretical models of ET. A nonadiabatic simulation scheme is employed, and no assumptions regarding the ET mechanism are needed. Thus, the predictive power of the simulations is preserved, while pushing the limits of the accessible time scales beyond microseconds. This model-based simulation scheme is applied to ET in various DNA species. Good agreement with the "full" atomistic nonadiabatic QM/MM scheme is observed for the archetypal DNA ET systems, the polyA sequence, as well as the sequences GTnGGG, containing adenines as bridge sites. Furthermore, ET in larger, more complex DNA sequences is simulated, and the results are discussed.
Collapse
Affiliation(s)
| | | | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen , 28759 Bremen, Germany
| | | |
Collapse
|
32
|
Marquetand P, Nogueira JJ, Mai S, Plasser F, González L. Challenges in Simulating Light-Induced Processes in DNA. Molecules 2016. [PMCID: PMC6155660 DOI: 10.3390/molecules22010049] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.
Collapse
|
33
|
Kubař T, Elstner M, Popescu B, Kleinekathöfer U. Polaron Effects on Charge Transport through Molecular Wires: A Multiscale Approach. J Chem Theory Comput 2016; 13:286-296. [DOI: 10.1021/acs.jctc.6b00879] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Tomáš Kubař
- Institute of Physical Chemistry & Center for Functional Nanostructures, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Bogdan Popescu
- Department
of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department
of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
34
|
Spencer J, Gajdos F, Blumberger J. FOB-SH: Fragment orbital-based surface hopping for charge carrier transport in organic and biological molecules and materials. J Chem Phys 2016. [DOI: 10.1063/1.4960144] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
|
35
|
Kranz JJ, Elstner M. Simulation of Singlet Exciton Diffusion in Bulk Organic Materials. J Chem Theory Comput 2016; 12:4209-21. [DOI: 10.1021/acs.jctc.6b00235] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Julian J. Kranz
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Marcus Elstner
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| |
Collapse
|
36
|
Heck A, Kranz JJ, Elstner M. Simulation of Temperature-Dependent Charge Transport in Organic Semiconductors with Various Degrees of Disorder. J Chem Theory Comput 2016; 12:3087-96. [DOI: 10.1021/acs.jctc.6b00215] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alexander Heck
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
- HEiKA - Heidelberg Karlsruhe Research Partnership, Heidelberg University, Karlsruhe Institute of Technology (KIT), Germany
| | - Julian J. Kranz
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
- HEiKA - Heidelberg Karlsruhe Research Partnership, Heidelberg University, Karlsruhe Institute of Technology (KIT), Germany
| |
Collapse
|
37
|
Cailliez F, Müller P, Firmino T, Pernot P, de la Lande A. Energetics of Photoinduced Charge Migration within the Tryptophan Tetrad of an Animal (6–4) Photolyase. J Am Chem Soc 2016; 138:1904-15. [DOI: 10.1021/jacs.5b10938] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Fabien Cailliez
- Laboratoire
de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Pavel Müller
- Institute
for Integrative Biology of the Cell (I2BC), CEA, CNRS, University
Paris-Sud, University Paris-Saclay, 91198 Gif-sur-Yvette
cedex, France
| | - Thiago Firmino
- Laboratoire
de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Pascal Pernot
- Laboratoire
de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Aurélien de la Lande
- Laboratoire
de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
38
|
Spencer J, Scalfi L, Carof A, Blumberger J. Confronting surface hopping molecular dynamics with Marcus theory for a molecular donor–acceptor system. Faraday Discuss 2016; 195:215-236. [DOI: 10.1039/c6fd00107f] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We investigate the performance of fewest switches surface hopping (SH) in describing electron transfer (ET) for a molecular donor–acceptor system. Computer simulations are carried out for a wide range of reorganisation energy (λ), electronic coupling strength (Hab) and driving force using our recently developed fragment orbital-based SH approach augmented with a simple decoherence correction. This methodology allows us to compute SH ET rates over more than four orders of magnitude, from the sub-picosecond to the nanosecond time regime. We find good agreement with semi-classical ET theory in the non-adiabatic ET regime. The correct scaling of the SH ET rate with electronic coupling strength is obtained and the Marcus inverted regime is reproduced, in line with previously reported results for a spin-boson model. Yet, we find that the SH ET rate falls below the semi-classical ET rate in the adiabatic regime, where the free energy barrier is in the order of kBT in our simulations. We explain this by first signatures of non-exponential population decay of the initial charge state. For even larger electronic couplings (Hab = λ/2), the free energy barrier vanishes and ET rates are no longer defined. At this point we observe a crossover from ET on the vibronic time scale to charge relaxation on the femtosecond time scale that is well described by thermally averaged Rabi oscillations. The extension of the analysis from the non-adiabatic limit to large electronic couplings and small or even vanishing activation barriers is relevant for our understanding of charge transport in organic semiconductors.
Collapse
Affiliation(s)
- Jacob Spencer
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Laura Scalfi
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
- Department of Chemistry
| | - Antoine Carof
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Jochen Blumberger
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| |
Collapse
|
39
|
Blumberger J. Recent Advances in the Theory and Molecular Simulation of Biological Electron Transfer Reactions. Chem Rev 2015; 115:11191-238. [DOI: 10.1021/acs.chemrev.5b00298] [Citation(s) in RCA: 238] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jochen Blumberger
- Department of Physics and
Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
40
|
Heck A, Kranz JJ, Kubař T, Elstner M. Multi-Scale Approach to Non-Adiabatic Charge Transport in High-Mobility Organic Semiconductors. J Chem Theory Comput 2015; 11:5068-82. [DOI: 10.1021/acs.jctc.5b00719] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexander Heck
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
- Heidelberg
Karlsruhe Research Partnership (HEiKA), Heidelberg University, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Julian J. Kranz
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Department
of Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
- Heidelberg
Karlsruhe Research Partnership (HEiKA), Heidelberg University, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
41
|
Ramos P, Papadakis M, Pavanello M. Performance of Frozen Density Embedding for Modeling Hole Transfer Reactions. J Phys Chem B 2015; 119:7541-57. [DOI: 10.1021/jp511275e] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Markos Papadakis
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
42
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 301] [Impact Index Per Article: 33.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
43
|
Computational modeling of photoexcitation in DNA single and double strands. Top Curr Chem (Cham) 2015; 356:89-122. [PMID: 24647841 DOI: 10.1007/128_2014_533] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The photoexcitation of DNA strands triggers extremely complex photoinduced processes, which cannot be understood solely on the basis of the behavior of the nucleobase building blocks. Decisive factors in DNA oligomers and polymers include collective electronic effects, excitonic coupling, hydrogen-bonding interactions, local steric hindrance, charge transfer, and environmental and solvent effects. This chapter surveys recent theoretical and computational efforts to model real-world excited-state DNA strands using a variety of established and emerging theoretical methods. One central issue is the role of localized vs delocalized excitations and the extent to which they determine the nature and the temporal evolution of the initial photoexcitation in DNA strands.
Collapse
|
44
|
Beratan DN, Liu C, Migliore A, Polizzi NF, Skourtis SS, Zhang P, Zhang Y. Charge transfer in dynamical biosystems, or the treachery of (static) images. Acc Chem Res 2015; 48:474-81. [PMID: 25307316 PMCID: PMC4333612 DOI: 10.1021/ar500271d] [Citation(s) in RCA: 124] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
![]()
The image is not the thing. Just as a pipe rendered
in an oil painting cannot be smoked, quantum mechanical coupling pathways
rendered on LCDs do not convey electrons. The aim of this Account
is to examine some of our recent discoveries regarding biological
electron transfer (ET) and transport mechanisms that emerge when one
moves beyond treacherous static views to dynamical frameworks. Studies over the last two decades introduced both atomistic detail
and macromolecule dynamics to the description of biological ET. The
first model to move beyond the structureless square-barrier tunneling
description is the Pathway model, which predicts how protein secondary
motifs and folding-induced through-bond and through-space tunneling
gaps influence kinetics. Explicit electronic structure theory is applied
routinely now to elucidate ET mechanisms, to capture pathway interferences,
and to treat redox cofactor electronic structure effects. Importantly,
structural sampling of proteins provides an understanding of how dynamics
may change the mechanisms of biological ET, as ET rates are exponentially
sensitive to structure. Does protein motion average out tunneling
pathways? Do conformational fluctuations gate biological ET? Are transient
multistate resonances produced by energy gap fluctuations? These questions
are becoming accessible as the static view of biological ET recedes
and dynamical viewpoints take center stage. This Account introduces
ET reactions at the core of bioenergetics, summarizes our team’s
progress toward arriving at an atomistic-level description, examines
how thermal fluctuations influence ET, presents metrics that characterize
dynamical effects on ET, and discusses applications in very long (micrometer
scale) bacterial nanowires. The persistence of structural effects
on the ET rates in the face of thermal fluctuations is considered.
Finally, the flickering resonance (FR) view of charge transfer is
presented to examine how fluctuations control low-barrier transport
among multiple groups in van der Waals contact. FR produces exponential
distance dependence in the absence of tunneling; the exponential character
emerges from the probability of matching multiple vibronically broadened
electronic energies within a tolerance defined by the rms coupling
among interacting groups. FR thus produces band like coherent transport
on the nanometer length scale, enabled by conformational fluctuations.
Taken as a whole, the emerging context for ET in dynamical biomolecules
provides a robust framework to design and interpret the inner workings
of bioenergetics from the molecular to the cellular scale and beyond,
with applications in biomedicine, biocatalysis, and energy science.
Collapse
|
45
|
Lüdemann G, Solov’yov IA, Kubař T, Elstner M. Solvent Driving Force Ensures Fast Formation of a Persistent and Well-Separated Radical Pair in Plant Cryptochrome. J Am Chem Soc 2015; 137:1147-56. [DOI: 10.1021/ja510550g] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Gesa Lüdemann
- Department
for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Ilia A. Solov’yov
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark
| | - Tomáš Kubař
- Department
for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Department
for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
46
|
Bashir A, Heck A, Narita A, Feng X, Nefedov A, Rohwerder M, Müllen K, Elstner M, Wöll C. Charge carrier mobilities in organic semiconductors: crystal engineering and the importance of molecular contacts. Phys Chem Chem Phys 2015; 17:21988-96. [DOI: 10.1039/c5cp03171k] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have conducted a combined experimental and theoretical study on the packing optimization of hexa-peri-hexabenzocoronene (HBC) as organic semiconductor.
Collapse
Affiliation(s)
- Asif Bashir
- Max-Planck-Institut für Eisenforschung GmbH
- 40237 Düsseldorf
- Germany
| | - Alexander Heck
- Institute of Physical Chemistry
- Karlsruhe Institute of Technology
- 76131 Karlsruhe
- Germany
- HEiKA – Heidelberg Karlsruhe Research Partnership
| | | | - Xinliang Feng
- Max-Planck Institute for Polymer Research
- 55128 Mainz
- Germany
| | - Alexei Nefedov
- Institute of Functional Interfaces
- Karlsruhe Institute of Technology
- 76344 Eggenstein-Leopoldshafen
- Germany
| | | | - Klaus Müllen
- Max-Planck Institute for Polymer Research
- 55128 Mainz
- Germany
| | - Marcus Elstner
- Institute of Physical Chemistry
- Karlsruhe Institute of Technology
- 76131 Karlsruhe
- Germany
- HEiKA – Heidelberg Karlsruhe Research Partnership
| | - Christof Wöll
- Institute of Functional Interfaces
- Karlsruhe Institute of Technology
- 76344 Eggenstein-Leopoldshafen
- Germany
| |
Collapse
|
47
|
Kubas A, Hoffmann F, Heck A, Oberhofer H, Elstner M, Blumberger J. Electronic couplings for molecular charge transfer: benchmarking CDFT, FODFT, and FODFTB against high-level ab initio calculations. J Chem Phys 2014; 140:104105. [PMID: 24628150 DOI: 10.1063/1.4867077] [Citation(s) in RCA: 127] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
We introduce a database (HAB11) of electronic coupling matrix elements (H(ab)) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H(ab) values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.
Collapse
Affiliation(s)
- Adam Kubas
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Felix Hoffmann
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Alexander Heck
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 6, 76131 Karlsruhe, Germany
| | - Harald Oberhofer
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 6, 76131 Karlsruhe, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
48
|
Watanabe HC, Kubař T, Elstner M. Size-Consistent Multipartitioning QM/MM: A Stable and Efficient Adaptive QM/MM Method. J Chem Theory Comput 2014; 10:4242-52. [DOI: 10.1021/ct5005593] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hiroshi C. Watanabe
- Center
for Biological Resources and Informatics, Tokyo Institute of Technology, B-62 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501 Japan
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstraße
12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstraße
12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstraße
12, 76131 Karlsruhe, Germany
| |
Collapse
|
49
|
Gajdos F, Valner S, Hoffmann F, Spencer J, Breuer M, Kubas A, Dupuis M, Blumberger J. Ultrafast Estimation of Electronic Couplings for Electron Transfer between π-Conjugated Organic Molecules. J Chem Theory Comput 2014; 10:4653-60. [DOI: 10.1021/ct500527v] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Fruzsina Gajdos
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Siim Valner
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Felix Hoffmann
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität Bochum, Universitätsstrasse
150, 44801 Bochum, Germany
| | - Jacob Spencer
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Marian Breuer
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Adam Kubas
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Michel Dupuis
- Pacific Northwest
National Laboratory, Richland, Washington 99354, United States
| | - Jochen Blumberger
- Department
of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| |
Collapse
|
50
|
Cui Q, Elstner M. Density functional tight binding: values of semi-empirical methods in an ab initio era. Phys Chem Chem Phys 2014; 16:14368-77. [PMID: 24850383 PMCID: PMC4836871 DOI: 10.1039/c4cp00908h] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.
Collapse
|