1
|
De PK, Jain A. Exciton energy transfer inside cavity-A benchmark study of polaritonic dynamics using the surface hopping method. J Chem Phys 2024; 161:054117. [PMID: 39105549 DOI: 10.1063/5.0216787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 07/22/2024] [Indexed: 08/07/2024] Open
Abstract
Strong coupling between the molecular system and photon inside the cavity generates polaritons, which can alter reaction rates by orders of magnitude. In this work, we benchmark the surface hopping method to simulate non-adiabatic dynamics in a cavity. The comparison is made against a numerically exact method (the hierarchical equations of motion) for a model system investigating excitonic energy transfer for a broad range of parameters. Surface hopping captures the effects of the radiation mode well, both at resonance and off-resonance. We have further investigated parameters that can increase or decrease the rate of population transfer, and we find that surface hopping in general can capture both effects well. Finally, we show that the dipole self-energy term within our parameter regime does not significantly affect the system's dynamics.
Collapse
Affiliation(s)
- Priyam Kumar De
- Department of Chemistry, Indian Institute of Technology, Mumbai 400076, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology, Mumbai 400076, India
| |
Collapse
|
2
|
Dines A, Ellis M, Blumberger J. Stabilized coupled trajectory mixed quantum-classical algorithm with improved energy conservation: CTMQC-EDI. J Chem Phys 2023; 159:234118. [PMID: 38117021 DOI: 10.1063/5.0183589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 11/26/2023] [Indexed: 12/21/2023] Open
Abstract
Coupled trajectory mixed quantum-classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularization procedures in certain circumstances. Such procedures in the latter can cause instabilities, leading to undesirable effects, such as energy drift and spurious population transfer, which is expected to become increasingly prevalent when the system gets larger as such events would happen more frequently. We propose a further modification to CTMQC-E, which includes a redefinition of the quantum momentum, CTMQC-EDI (double intercept), such that it has no formal divergences. We then show for Tully models I-III and the double arch model that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI, therefore, shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.
Collapse
Affiliation(s)
- Aaron Dines
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Matthew Ellis
- 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
|
3
|
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
|
4
|
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] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 03/01/2023] [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
|
5
|
Ten Brink M, Gräber S, Hopjan M, Jansen D, Stolpp J, Heidrich-Meisner F, Blöchl PE. Real-time non-adiabatic dynamics in the one-dimensional Holstein model: Trajectory-based vs exact methods. J Chem Phys 2022; 156:234109. [PMID: 35732530 DOI: 10.1063/5.0092063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We benchmark a set of quantum-chemistry methods, including multitrajectory Ehrenfest, fewest-switches surface-hopping, and multiconfigurational-Ehrenfest dynamics, against exact quantum-many-body techniques by studying real-time dynamics in the Holstein model. This is a paradigmatic model in condensed matter theory incorporating a local coupling of electrons to Einstein phonons. For the two-site and three-site Holstein model, we discuss the exact and quantum-chemistry methods in terms of the Born-Huang formalism, covering different initial states, which either start on a single Born-Oppenheimer surface, or with the electron localized to a single site. For extended systems with up to 51 sites, we address both the physics of single Holstein polarons and the dynamics of charge-density waves at finite electron densities. For these extended systems, we compare the quantum-chemistry methods to exact dynamics obtained from time-dependent density matrix renormalization group calculations with local basis optimization (DMRG-LBO). We observe that the multitrajectory Ehrenfest method, in general, only captures the ultrashort time dynamics accurately. In contrast, the surface-hopping method with suitable corrections provides a much better description of the long-time behavior but struggles with the short-time description of coherences between different Born-Oppenheimer states. We show that the multiconfigurational Ehrenfest method yields a significant improvement over the multitrajectory Ehrenfest method and can be converged to the exact results in small systems with moderate computational efforts. We further observe that for extended systems, this convergence is slower with respect to the number of configurations. Our benchmark study demonstrates that DMRG-LBO is a useful tool for assessing the quality of the quantum-chemistry methods.
Collapse
Affiliation(s)
- M Ten Brink
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - S Gräber
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M Hopjan
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - D Jansen
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - J Stolpp
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - F Heidrich-Meisner
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - P E Blöchl
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
6
|
Giannini S, Peng WT, Cupellini L, Padula D, Carof A, Blumberger J. Exciton transport in molecular organic semiconductors boosted by transient quantum delocalization. Nat Commun 2022; 13:2755. [PMID: 35589694 PMCID: PMC9120088 DOI: 10.1038/s41467-022-30308-5] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Accepted: 04/26/2022] [Indexed: 11/09/2022] Open
Abstract
Designing molecular materials with very large exciton diffusion lengths would remove some of the intrinsic limitations of present-day organic optoelectronic devices. Yet, the nature of excitons in these materials is still not sufficiently well understood. Here we present Frenkel exciton surface hopping, an efficient method to propagate excitons through truly nano-scale materials by solving the time-dependent Schrödinger equation coupled to nuclear motion. We find a clear correlation between diffusion constant and quantum delocalization of the exciton. In materials featuring some of the highest diffusion lengths to date, e.g. the non-fullerene acceptor Y6, the exciton propagates via a transient delocalization mechanism, reminiscent to what was recently proposed for charge transport. Yet, the extent of delocalization is rather modest, even in Y6, and found to be limited by the relatively large exciton reorganization energy. On this basis we chart out a path for rationally improving exciton transport in organic optoelectronic materials.
Collapse
Affiliation(s)
- Samuele Giannini
- Department of Physics and Astronomy and Thomas Young Centre, University College London, WC1E 6BT, London, UK.
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc 20, 7000, Mons, Belgium.
| | - Wei-Tao Peng
- Department of Physics and Astronomy and Thomas Young Centre, University College London, WC1E 6BT, London, UK
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via G. Moruzzi 13, 56124, Pisa, Italy
| | - Daniele Padula
- Dipartimento di Biotecnologie, Chimica e Farmacia, Universitá di Siena, Via A. Moro 2, 53100, Siena, Italy
| | - Antoine Carof
- Laboratoire de Physique et Chimie Théoriques, CNRS, UMR No. 7019, Université de Lorraine, BP 239, 54506, Vandoeuvre-lés-Nancy Cedex, France
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, WC1E 6BT, London, UK.
| |
Collapse
|
7
|
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
|
8
|
Ahart CS, Rosso KM, Blumberger J. Electron and Hole Mobilities in Bulk Hematite from Spin-Constrained Density Functional Theory. J Am Chem Soc 2022; 144:4623-4632. [PMID: 35239359 PMCID: PMC9097473 DOI: 10.1021/jacs.1c13507] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Transition metal oxide materials have attracted much attention for photoelectrochemical water splitting, but problems remain, e.g. the sluggish transport of excess charge carriers in these materials, which is not well understood. Here we use periodic, spin-constrained and gap-optimized hybrid density functional theory to uncover the nature and transport mechanism of holes and excess electrons in a widely used water splitting material, bulk-hematite (α-Fe2O3). We find that upon ionization the hole relaxes from a delocalized band state to a polaron localized on a single iron atom with localization induced by tetragonal distortion of the six surrounding iron-oxygen bonds. This distortion is responsible for sluggish hopping transport in the Fe-bilayer, characterized by an activation energy of 70 meV and a hole mobility of 0.031 cm2/(V s). By contrast, the excess electron induces a smaller distortion of the iron-oxygen bonds resulting in delocalization over two neighboring Fe units. We find that 2-site delocalization is advantageous for charge transport due to the larger spatial displacements per transfer step. As a result, the electron mobility is predicted to be a factor of 3 higher than the hole mobility, 0.098 cm2/(V s), in qualitative agreement with experimental observations. This work provides new fundamental insight into charge carrier transport in hematite with implications for its photocatalytic activity.
Collapse
Affiliation(s)
- Christian S Ahart
- Department of Physics and Astronomy, University College London, London WC1E 6BT, U.K
| | - Kevin M Rosso
- 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
|
9
|
Ellis M, Yang H, Giannini S, Ziogos OG, Blumberger J. Impact of Nanoscale Morphology on Charge Carrier Delocalization and Mobility in an Organic Semiconductor. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2021; 33:e2104852. [PMID: 34558126 DOI: 10.1002/adma.202104852] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Revised: 08/11/2021] [Indexed: 06/13/2023]
Abstract
A central challenge of organic semiconductor research is to make cheap, disordered materials that exhibit high electrical conductivity. Unfortunately, this endeavor is hampered by the poor fundamental understanding of the relationship between molecular packing structure and charge carrier mobility. Here a novel computational methodology is presented that fills this gap. Using a melt-quench procedure it is shown that amorphous pentacene spontaneously self-assembles to nanocrystalline structures that, at long quench times, form the characteristic herringbone layer of the single crystal. Quantum dynamical simulations of electron hole transport show a clear correlation between the crystallinity of the sample, the quantum delocalization, and the mobility of the charge carrier. Surprisingly, the long-held belief that charge carriers form relatively localized polarons in disordered OS is only valid for fully amorphous structures-for nanocrystalline and crystalline samples, significant charge carrier delocalization over several nanometers occurs that underpins their improved conductivities. The good agreement with experimentally available data makes the presented methodology a robust computational tool for the predictive engineering of disordered organic materials.
Collapse
Affiliation(s)
- Matthew Ellis
- Department of Physics and Astronomy and Thomas Young Centre London, University College London, Gower Street, London, WC1E 6BT, UK
| | - Hui Yang
- Department of Physics and Astronomy and Thomas Young Centre London, University College London, Gower Street, London, WC1E 6BT, UK
| | - Samuele Giannini
- Department of Physics and Astronomy and Thomas Young Centre London, University College London, Gower Street, London, WC1E 6BT, UK
| | - Orestis G Ziogos
- Department of Physics and Astronomy and Thomas Young Centre London, University College London, Gower Street, London, WC1E 6BT, UK
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre London, University College London, Gower Street, London, WC1E 6BT, UK
| |
Collapse
|
10
|
Ghan S, Kunkel C, Reuter K, Oberhofer H. Improved Projection-Operator Diabatization Schemes for the Calculation of Electronic Coupling Values. J Chem Theory Comput 2020; 16:7431-7443. [DOI: 10.1021/acs.jctc.0c00887] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Simiam Ghan
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Christian Kunkel
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
- Fritz-Haber Institute of the Max-Planck Society, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technical University of Munich, Lichtenbergstraße 4, D-85747 Garching, Germany
| |
Collapse
|
11
|
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
|
12
|
Ahart CS, Blumberger J, Rosso KM. Polaronic structure of excess electrons and holes for a series of bulk iron oxides. Phys Chem Chem Phys 2020; 22:10699-10709. [DOI: 10.1039/c9cp06482f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
With the use of a gap-optimized hybrid functional and large supercells, it is found that while the electron hole polaron generally localises onto a single iron site, the electron polaron localises across two iron sites of the same spin layer.
Collapse
Affiliation(s)
- Christian S. Ahart
- 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
|
13
|
Ghosh S, Giannini S, Lively K, Blumberger J. Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping. Faraday Discuss 2020; 221:501-525. [DOI: 10.1039/c9fd00046a] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Exploring effects of quantizing nuclei in non-adiabatic dynamics for simulating charge transfer in a dimer of “ethylene-like-molecules” at different temperatures.
Collapse
Affiliation(s)
- Soumya Ghosh
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Samuele Giannini
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Kevin Lively
- 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
- Institute for Advanced Study
| |
Collapse
|
14
|
Kim T, Kim W, Vakuliuk O, Gryko DT, Kim D. Two-Step Charge Separation Passing Through the Partial Charge-Transfer State in a Molecular Dyad. J Am Chem Soc 2019; 142:1564-1573. [DOI: 10.1021/jacs.9b12016] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Taeyeon Kim
- Department of Chemistry and Spectroscopy Laboratory for Functional π-Electronic Systems, Yonsei University, Seoul 03722, Korea
| | - Woojae Kim
- Department of Chemistry and Spectroscopy Laboratory for Functional π-Electronic Systems, Yonsei University, Seoul 03722, Korea
| | - Olena Vakuliuk
- Institute of Organic Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Daniel T. Gryko
- Institute of Organic Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Dongho Kim
- Department of Chemistry and Spectroscopy Laboratory for Functional π-Electronic Systems, Yonsei University, Seoul 03722, Korea
| |
Collapse
|
15
|
Giannini S, Carof A, Ellis M, Yang H, Ziogos OG, Ghosh S, Blumberger J. Quantum localization and delocalization of charge carriers in organic semiconducting crystals. Nat Commun 2019; 10:3843. [PMID: 31451687 PMCID: PMC6710274 DOI: 10.1038/s41467-019-11775-9] [Citation(s) in RCA: 87] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Accepted: 08/02/2019] [Indexed: 02/06/2023] Open
Abstract
Charge carrier transport in organic semiconductors is at the heart of many revolutionary technologies ranging from organic transistors, light-emitting diodes, flexible displays and photovoltaic cells. Yet, the nature of charge carriers and their transport mechanism in these materials is still unclear. Here we show that by solving the time-dependent electronic Schrödinger equation coupled to nuclear motion for eight organic molecular crystals, the excess charge carrier forms a polaron delocalized over up to 10-20 molecules in the most conductive crystals. The polaron propagates through the crystal by diffusive jumps over several lattice spacings at a time during which it expands more than twice its size. Computed values for polaron size and charge mobility are in excellent agreement with experimental estimates and correlate very well with the recently proposed transient localization theory.
Collapse
Affiliation(s)
- Samuele Giannini
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Antoine Carof
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Matthew Ellis
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Hui Yang
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Orestis George Ziogos
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Soumya Ghosh
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, WC1E 6BT, UK.
- Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748, Garching, Germany.
| |
Collapse
|
16
|
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
|
17
|
Sifain AE, Wang L, Tretiak S, Prezhdo OV. Numerical tests of coherence-corrected surface hopping methods using a donor-bridge-acceptor model system. J Chem Phys 2019; 150:194104. [DOI: 10.1063/1.5092999] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew E. Sifain
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089-0485, USA
- Theoretical Division, Center for Nonlinear Studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Angeles, New Mexico 87545, USA
| | - Linjun Wang
- Center for Chemistry of Novel and High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Angeles, New Mexico 87545, USA
- Skolkovo Institute of Science and Technology, Moscow 143026, Russia
| | - Oleg V. Prezhdo
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089-0485, USA
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| |
Collapse
|
18
|
Jensen KT, Benson RL, Cardamone S, Thom AJW. Modeling Electron Transfers Using Quasidiabatic Hartree–Fock States. J Chem Theory Comput 2018; 14:4629-4639. [DOI: 10.1021/acs.jctc.8b00379] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kristopher T. Jensen
- Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1TN, UK
| | - Raz L. Benson
- Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1TN, UK
| | - Salvatore Cardamone
- Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1TN, UK
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1TN, UK
| |
Collapse
|
19
|
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
|
20
|
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
|
21
|
Mattiat J, Richardson JO. Effects of tunnelling and asymmetry for system-bath models of electron transfer. J Chem Phys 2018; 148:102311. [PMID: 29544261 DOI: 10.1063/1.5001116] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We apply the newly derived nonadiabatic golden-rule instanton theory to asymmetric models describing electron-transfer in solution. The models go beyond the usual spin-boson description and have anharmonic free-energy surfaces with different values for the reactant and product reorganization energies. The instanton method gives an excellent description of the behaviour of the rate constant with respect to asymmetry for the whole range studied. We derive a general formula for an asymmetric version of the Marcus theory based on the classical limit of the instanton and find that this gives significant corrections to the standard Marcus theory. A scheme is given to compute this rate based only on equilibrium simulations. We also compare the rate constants obtained by the instanton method with its classical limit to study the effect of tunnelling and other quantum nuclear effects. These quantum effects can increase the rate constant by orders of magnitude.
Collapse
Affiliation(s)
- Johann Mattiat
- Laboratory of Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | | |
Collapse
|
22
|
Sato K, Pradhan E, Asahi R, Akimov AV. Charge transfer dynamics at the boron subphthalocyanine chloride/C60 interface: non-adiabatic dynamics study with Libra-X. Phys Chem Chem Phys 2018; 20:25275-25294. [DOI: 10.1039/c8cp03841d] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The Libra-X software for non-adiabatic molecular dynamics is reported. It is used to comprehensively study the charge transfer dynamics at the boron subphtalocyanine chloride (SubPc)/fullerene (C60) interface.
Collapse
Affiliation(s)
- Kosuke Sato
- Toyota Central Research and Development Laboratories, Inc
- Nagakute
- Japan
| | - Ekadashi Pradhan
- Department of Chemistry
- University at Buffalo
- The State University of New York
- New York 14260-3000
- USA
| | - Ryoji Asahi
- Toyota Central Research and Development Laboratories, Inc
- Nagakute
- Japan
| | - Alexey V. Akimov
- Department of Chemistry
- University at Buffalo
- The State University of New York
- New York 14260-3000
- USA
| |
Collapse
|
23
|
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
|
24
|
Steiner F, Poelking C, Niedzialek D, Andrienko D, Nelson J. Influence of orientation mismatch on charge transport across grain boundaries in tri-isopropylsilylethynyl (TIPS) pentacene thin films. Phys Chem Chem Phys 2017; 19:10854-10862. [DOI: 10.1039/c6cp06436a] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present a multi-scale model for charge transport across grain boundaries in molecular electronic materials that incorporates packing disorder, electrostatic and polarisation effects.
Collapse
Affiliation(s)
- Florian Steiner
- Department of Physics and Centre for Plastic Electronics
- Imperial College London
- London SW7 2AZ
- UK
| | - Carl Poelking
- Max Planck Institute for Polymer Research
- 55128 Mainz
- Germany
| | - Dorota Niedzialek
- Department of Physics and Centre for Plastic Electronics
- Imperial College London
- London SW7 2AZ
- UK
| | | | - Jenny Nelson
- Department of Physics and Centre for Plastic Electronics
- Imperial College London
- London SW7 2AZ
- UK
| |
Collapse
|