1
|
Stojanovic L, Giannini S, Blumberger J. Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics. J Chem Theory Comput 2024; 20:6241-6252. [PMID: 38967252 PMCID: PMC11270823 DOI: 10.1021/acs.jctc.4c00605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Revised: 06/11/2024] [Accepted: 06/12/2024] [Indexed: 07/06/2024]
Abstract
Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.
Collapse
Affiliation(s)
- Ljiljana Stojanovic
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, U.K.
| | - Samuele Giannini
- Institute
of Chemistry of OrganoMetallic Compounds, National Research Council (ICCOM-CNR), Pisa I-56124, Italy
| | - Jochen Blumberger
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, U.K.
| |
Collapse
|
2
|
Walker DW, Sing CE. Effect of Hydrodynamic Interactions and Flow on Charge Transport in Redox-Active Polymer Solutions. J Phys Chem B 2024; 128:1796-1811. [PMID: 38330099 DOI: 10.1021/acs.jpcb.3c07657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Redox-active polymers (RAPs) are a subclass of polyelectrolytes that can store charge and undergo redox self-exchange reactions. RAPs are of great interest in the field of redox flow batteries (RFBs) due to their ability to quickly charge and discharge, their chemical modularity, and their molecular size. However, designing RAPs for efficient charge transport at the molecular level requires a fundamental understanding of the charge transport mechanisms that occur in RFBs. Previous work from our group has explored these mechanisms, and in this paper, we seek to improve upon the previous model by incorporating both hydrodynamic interactions (HIs) and out-of-equilibrium dynamics, which are both highly pertinent to flow battery systems. We use a hybrid Brownian dynamics and Monte Carlo simulation to model redox-active polymer chains in both dilute and semidilute solutions. This model is used to show that HI is an important feature when charge hopping is not the major mechanism for charge displacement and leads to more rapid segmental and translational motion of polymer chains that expedites charge transport at low polymer concentrations. We demonstrate that strong extensional flows may result in either enhanced or decreased transport depending on the fraction of charges present on the RAP chain. We show that flow not only can promote charge transport by extending polymer conformations but can also suppress nonadjacent charge hopping processes that are important for transport at high charge fractions. Shear flows can similarly enhance charge transport through chain extension, but tumbling dynamics lead to oscillatory displacements that become dominant features with high charge fractions and strong flows.
Collapse
Affiliation(s)
- Dejuante W Walker
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Charles E Sing
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
3
|
Li W, Giannini S, Quarti C, Hou Z, Prezhdo OV, Beljonne D. Interlayer Charge Transport in 2D Lead Halide Perovskites from First Principles. J Chem Theory Comput 2023; 19:9403-9415. [PMID: 38048307 DOI: 10.1021/acs.jctc.3c00904] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/06/2023]
Abstract
We report on the implementation of a versatile projection-operator diabatization approach to calculate electronic coupling integrals in layered periodic systems. The approach is applied to model charge transport across the saturated organic spacers in two-dimensional (2D) lead halide perovskites. The calculations yield out-of-plane charge transfer rates that decay exponentially with the increasing length of the alkyl chain, range from a few nanoseconds to milliseconds, and are supportive of a hopping mechanism. Most importantly, we show that the charge carriers strongly couple to distortions of the Pb-I framework and that accounting for the associated nonlocal dynamic disorder increases the thermally averaged interlayer rates by a few orders of magnitude compared to the frozen-ion 0 K-optimized structure. Our formalism provides the first comprehensive insight into the role of the organic spacer cation on vertical transport in 2D lead halide perovskites and can be readily extended to functional π-conjugated spacers, where we expect the improved energy alignment with the inorganic layout to speed up the charge transfer between the semiconducting layers.
Collapse
Affiliation(s)
- Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, China
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, B-7000 Mons, Belgium
| | - Samuele Giannini
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, B-7000 Mons, Belgium
| | - Claudio Quarti
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, B-7000 Mons, Belgium
| | - Zhufeng Hou
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, China
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, B-7000 Mons, Belgium
| |
Collapse
|
4
|
Giannini S, Di Virgilio L, Bardini M, Hausch J, Geuchies JJ, Zheng W, Volpi M, Elsner J, Broch K, Geerts YH, Schreiber F, Schweicher G, Wang HI, Blumberger J, Bonn M, Beljonne D. Transiently delocalized states enhance hole mobility in organic molecular semiconductors. NATURE MATERIALS 2023; 22:1361-1369. [PMID: 37709929 DOI: 10.1038/s41563-023-01664-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 08/14/2023] [Indexed: 09/16/2023]
Abstract
Evidence shows that charge carriers in organic semiconductors self-localize because of dynamic disorder. Nevertheless, some organic semiconductors feature reduced mobility at increasing temperature, a hallmark for delocalized band transport. Here we present the temperature-dependent mobility in two record-mobility organic semiconductors: dinaphtho[2,3-b:2',3'-f]thieno[3,2-b]-thiophene (DNTT) and its alkylated derivative, C8-DNTT-C8. By combining terahertz photoconductivity measurements with atomistic non-adiabatic molecular dynamics simulations, we show that while both crystals display a power-law decrease of the mobility (μ) with temperature (T) following μ ∝ T -n, the exponent n differs substantially. Modelling reveals that the differences between the two chemically similar semiconductors can be traced to the delocalization of the different states that are thermally accessible by charge carriers, which in turn depends on their specific electronic band structure. The emerging picture is that of holes surfing on a dynamic manifold of vibrationally dressed extended states with a temperature-dependent mobility that provides a sensitive fingerprint for the underlying density of states.
Collapse
Affiliation(s)
- Samuele Giannini
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium.
| | | | - Marco Bardini
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium
| | - Julian Hausch
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | | | - Wenhao Zheng
- Max Planck Institute for Polymer Research, Mainz, Germany
| | - Martina Volpi
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Jan Elsner
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, UK
| | - Katharina Broch
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | - Yves H Geerts
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
- International Solvay Institutes for Physics and Chemistry, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Frank Schreiber
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | - Guillaume Schweicher
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Hai I Wang
- Max Planck Institute for Polymer Research, Mainz, Germany.
- Nanophotonics, Debye Institute for Nanomaterials Science, Utrecht University, Utrecht, The Netherlands.
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, UK
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Mainz, Germany.
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium.
| |
Collapse
|
5
|
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
|
6
|
Craven GT, Nitzan A. Electron hopping heat transport in molecules. J Chem Phys 2023; 158:2887563. [PMID: 37125714 DOI: 10.1063/5.0144248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Accepted: 04/10/2023] [Indexed: 05/02/2023] Open
Abstract
The realization of single-molecule thermal conductance measurements has driven the need for theoretical tools to describe conduction processes that occur over atomistic length scales. In macroscale systems, the principle that is typically used to understand thermal conductivity is Fourier's law. At molecular length scales, however, deviations from Fourier's law are common in part because microscale thermal transport properties typically depend on the complex interplay between multiple heat conduction mechanisms. Here, the thermal transport properties that arise from electron transfer across a thermal gradient in a molecular conduction junction are examined theoretically. We illustrate how transport in a model junction is affected by varying the electronic structure and length of the molecular bridge in the junction as well as the strength of the coupling between the bridge and its surrounding environment. Three findings are of note: First, the transport properties can vary significantly depending on the characteristics of the molecular bridge and its environment; second, the system's thermal conductance commonly deviates from Fourier's law; and third, in properly engineered systems, the magnitude of electron hopping thermal conductance is similar to what has been measured in single-molecule devices.
Collapse
Affiliation(s)
- Galen T Craven
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
7
|
Schütze Y, Gayen D, Palczynski K, de Oliveira Silva R, Lu Y, Tovar M, Partovi-Azar P, Bande A, Dzubiella J. How Regiochemistry Influences Aggregation Behavior and Charge Transport in Conjugated Organosulfur Polymer Cathodes for Lithium-Sulfur Batteries. ACS NANO 2023; 17:7889-7900. [PMID: 37014093 PMCID: PMC10141565 DOI: 10.1021/acsnano.3c01523] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 03/29/2023] [Indexed: 06/19/2023]
Abstract
For lithium-sulfur (Li-S) batteries to become competitive, they require high stability and energy density. Organosulfur polymer-based cathodes have recently shown promising performance due to their ability to overcome common limitations of Li-S batteries, such as the insulating nature of sulfur. In this study, we use a multiscale modeling approach to explore the influence of the regiochemistry of a conjugated poly(4-(thiophene-3-yl)benzenethiol) (PTBT) polymer on its aggregation behavior and charge transport. Classical molecular dynamics simulations of the self-assembly of polymer chains with different regioregularity show that a head-to-tail/head-to-tail regularity can form a well-ordered crystalline phase of planar chains allowing for fast charge transport. Our X-ray diffraction measurements, in conjunction with our predicted crystal structure, confirm the presence of crystalline phases in the electropolymerized PTBT polymer. We quantitatively describe the charge transport in the crystalline phase in a band-like regime. Our results give detailed insights into the interplay between microstructural and electrical properties of conjugated polymer cathode materials, highlighting the effect of polymer chain regioregularity on its charge transport properties.
Collapse
Affiliation(s)
- Yannik Schütze
- Research
Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Theoretical
Chemistry, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Diptesh Gayen
- Applied Theoretical
Physics - Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder-Straße 3, 79104 Freiburg, Germany
| | - Karol Palczynski
- Research
Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Ranielle de Oliveira Silva
- Department
Electrochemical Energy Storage, Helmholtz-Zentrum
Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Yan Lu
- Department
Electrochemical Energy Storage, Helmholtz-Zentrum
Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Institute
of Chemistry, University of Potsdam, Am Neuen Palais 10, 14469 Potsdam, Germany
| | - Michael Tovar
- Department
Structure and Dynamics of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Pouya Partovi-Azar
- Institute
for Chemistry, Martin Luther Universität
Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Annika Bande
- Theory of
Electron Dynamics and Spectroscopy, Helmholtz-Zentrum
Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Joachim Dzubiella
- Research
Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
- Applied Theoretical
Physics - Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder-Straße 3, 79104 Freiburg, Germany
| |
Collapse
|
8
|
Willson JT, Liu W, Balzer D, Kassal I. Jumping Kinetic Monte Carlo: Fast and Accurate Simulations of Partially Delocalized Charge Transport in Organic Semiconductors. J Phys Chem Lett 2023; 14:3757-3764. [PMID: 37044057 DOI: 10.1021/acs.jpclett.3c00388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Developing devices using disordered organic semiconductors requires accurate and practical models of charge transport. In these materials, charge transport occurs through partially delocalized states in an intermediate regime between localized hopping and delocalized band conduction. Partial delocalization can increase mobilities by orders of magnitude compared to those with conventional hopping, making it important for the design of materials and devices. Although delocalization, disorder, and polaron formation can be described using delocalized kinetic Monte Carlo (dKMC), it is a computationally expensive method. Here, we develop jumping kinetic Monte Carlo (jKMC), a model that approaches the accuracy of dKMC for modest amounts of delocalization (such as those found in disordered organic semiconductors), with a computational cost comparable to that of conventional hopping. jKMC achieves its computational performance by modeling conduction using identical spherical polarons, yielding a simple delocalization correction to the Marcus hopping rate that allows polarons to jump over their nearest neighbors. jKMC can be used in regimes of partial delocalization inaccessible to dKMC to show that modest delocalization can increase mobilities by as much as 2 orders of magnitude.
Collapse
Affiliation(s)
- Jacob T Willson
- School of Chemistry, University of Sydney, Sydney, NSW 2006, Australia
| | - William Liu
- School of Chemistry, University of Sydney, Sydney, NSW 2006, Australia
| | - Daniel Balzer
- School of Chemistry, University of Sydney, Sydney, NSW 2006, Australia
| | - Ivan Kassal
- School of Chemistry, University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
9
|
Tsai H, Ghosh D, Kinigstein E, Dryzhakov B, Driscoll H, Owczarek M, Hu B, Zhang X, Tretiak S, Nie W. Light-Induced Structural Dynamics and Charge Transport in Layered Halide Perovskite Thin Films. NANO LETTERS 2023; 23:429-436. [PMID: 36603204 DOI: 10.1021/acs.nanolett.2c03403] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The dynamic nature of the metal halide perovskite lattice upon photoexcitation plays a vital role in their properties. Here we report an observation of light-induced structure dynamics in quasi-2D Ruddlesden-Popper phase perovskite thin films and its impact on the carrier transport properties. By a time-resolved X-ray scattering technique, we observe a rapid lattice expansion upon photoexcitation, followed by a slow relaxation over the course of 100 ns in the dark. Theoretical modeling suggests that the expansion originates from the lattice's thermal fluctuations caused by photon energy deposition. Power dependent optical spectroscopy and photoconductivity indicate that high laser powers triggered a strong local structural disorder, which increased the charge dissociation activation energy that results in localized transport. Our study investigates the impact of laser energy deposition on the lattices and the subsequent carrier transport properties, that are relevant to device operations.
Collapse
Affiliation(s)
- Hsinhan Tsai
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
- Department of Chemistry, University of California, Berkeley, Berkeley, California94720, United States
| | - Dibyajyoti Ghosh
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
- Department of Materials Science and Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi110016, India
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi110016, India
| | - Eli Kinigstein
- X-ray Science Division, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Bogdan Dryzhakov
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee37996, United States
| | - Honora Driscoll
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Magdalena Owczarek
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Bin Hu
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee37996, United States
| | - Xiaoyi Zhang
- X-ray Science Division, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Sergei Tretiak
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Wanyi Nie
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| |
Collapse
|
10
|
Balzer D, Kassal I. Even a little delocalization produces large kinetic enhancements of charge-separation efficiency in organic photovoltaics. SCIENCE ADVANCES 2022; 8:eabl9692. [PMID: 35960797 PMCID: PMC9374333 DOI: 10.1126/sciadv.abl9692] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Accepted: 06/28/2022] [Indexed: 05/25/2023]
Abstract
In organic photovoltaics, charges can separate efficiently even if their Coulomb attraction is an order of magnitude greater than the available thermal energy. Delocalization has been suggested to explain this fact, because it could increase the initial separation of charges in the charge-transfer (CT) state, reducing their attraction. However, understanding the mechanism requires a kinetic model of delocalized charge separation, which has proven difficult because it involves tracking the correlated quantum-mechanical motion of the electron and the hole in large simulation boxes required for disordered materials. Here, we report the first three-dimensional simulations of charge-separation dynamics in the presence of disorder, delocalization, and polaron formation, finding that even slight delocalization, across less than two molecules, can substantially enhance the charge-separation efficiency, even starting with thermalized CT states. Delocalization does not enhance efficiency by reducing the Coulomb attraction; instead, the enhancement is a kinetic effect produced by the increased overlap of electronic states.
Collapse
|
11
|
Peng WT, Brey D, Giannini S, Dell’Angelo D, Burghardt I, Blumberger J. Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree. J Phys Chem Lett 2022; 13:7105-7112. [PMID: 35900333 PMCID: PMC9376959 DOI: 10.1021/acs.jpclett.2c01928] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 07/19/2022] [Indexed: 05/20/2023]
Abstract
Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Dominik Brey
- Institute
of Physical and Theoretical Chemistry, Goethe
University Frankfurt, Max-von-Laue-Strasse 7, 60438 Frankfurt am Main, Germany
| | - Samuele Giannini
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - David Dell’Angelo
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Irene Burghardt
- Institute
of Physical and Theoretical Chemistry, Goethe
University Frankfurt, Max-von-Laue-Strasse 7, 60438 Frankfurt am Main, Germany
| | - Jochen Blumberger
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
12
|
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
|
13
|
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
|
14
|
Qiu J, Lu Y, Wang L. Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States. J Chem Theory Comput 2022; 18:2803-2815. [PMID: 35380833 DOI: 10.1021/acs.jctc.2c00130] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We present a multilayer subsystem surface hopping (MSSH) method to deal with nonadiabatic dynamics in large-scale systems. A small subsystem instead of the full system is adopted for surface hopping and is updated on-the-fly to achieve a reliable description of important adiabatic states and the wave function evolution. Additional subsystems for molecular dynamics and statistical description are introduced to further improve the simulation reliability. The global flux hopping probabilities with optimal state assignments are utilized to treat the complex surface crossings. As demonstrated in a series of one- and two-dimensional Holstein models with up to hundreds of thousands of states, MSSH shows weak parameter dependence in all investigated systems. Especially, the computational costs are reduced by 2-6 orders of magnitude compared to traditional surface hopping simulations in full systems, and size-independent results are achieved with a large time-step size of 2-5 fs. The new method is compatible with different decoherence correction strategies and achieves a much better balance between efficiency and reliability, thus promising for applications in general charge and exciton dynamics simulations.
Collapse
Affiliation(s)
- Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Yao Lu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
Dilmurat R, Prodhan S, Wang L, Beljonne D. Thermally activated intra-chain charge transport in high charge-carrier mobility copolymers. J Chem Phys 2022; 156:084115. [DOI: 10.1063/5.0082569] [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/15/2022] Open
Abstract
Disordered or even seemingly amorphous, donor–acceptor type, conjugated copolymers with high charge-carrier mobility have emerged as a new class of functional materials, where transport along the conjugated backbone is key. Here, we report on non-adiabatic molecular dynamics simulations of charge-carrier transport along chains of poly (indacenodithiophene-co-benzothiadiazole), within a model Hamiltonian parameterized against first-principles calculations. We predict thermally activated charge transport associated with a slightly twisted ground-state conformation, on par with experimental results. Our results also demonstrate that the energy mismatch between the hole on the donor vs the acceptor units of the copolymer drives localization of the charge carriers and limits the intra-chain charge-carrier mobility. We predict that room-temperature mobility values in excess of 10 cm2 V−1 s−1 can be achieved through proper chemical tuning of the component monomer units.
Collapse
Affiliation(s)
- Rishat Dilmurat
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| | - Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| |
Collapse
|
18
|
Wang Z, Dong J, Qiu J, Wang L. All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning. ACS APPLIED MATERIALS & INTERFACES 2022; 14:22929-22940. [PMID: 35100503 DOI: 10.1021/acsami.1c22181] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.
Collapse
Affiliation(s)
- Zedong Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiawei Dong
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
19
|
Wu R, Matta M, Paulsen BD, Rivnay J. Operando Characterization of Organic Mixed Ionic/Electronic Conducting Materials. Chem Rev 2022; 122:4493-4551. [PMID: 35026108 DOI: 10.1021/acs.chemrev.1c00597] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Operando characterization plays an important role in revealing the structure-property relationships of organic mixed ionic/electronic conductors (OMIECs), enabling the direct observation of dynamic changes during device operation and thus guiding the development of new materials. This review focuses on the application of different operando characterization techniques in the study of OMIECs, highlighting the time-dependent and bias-dependent structure, composition, and morphology information extracted from these techniques. We first illustrate the needs, requirements, and challenges of operando characterization then provide an overview of relevant experimental techniques, including spectroscopy, scattering, microbalance, microprobe, and electron microscopy. We also compare different in silico methods and discuss the interplay of these computational methods with experimental techniques. Finally, we provide an outlook on the future development of operando for OMIEC-based devices and look toward multimodal operando techniques for more comprehensive and accurate description of OMIECs.
Collapse
Affiliation(s)
- Ruiheng Wu
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Micaela Matta
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Bryan D Paulsen
- Department of Biomedical Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Jonathan Rivnay
- Department of Biomedical Engineering, Northwestern University, Evanston, Illinois 60208, United States.,Simpson Querrey Institute, Northwestern University, Chicago, Illinois 60611, United States
| |
Collapse
|
20
|
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
|
21
|
Han Y, Nickle C, Maglione MS, Karuppannan SK, Casado‐Montenegro J, Qi D, Chen X, Tadich A, Cowie B, Mas‐Torrent M, Rovira C, Cornil J, Veciana J, del Barco E, Nijhuis CA. Bias-Polarity-Dependent Direct and Inverted Marcus Charge Transport Affecting Rectification in a Redox-Active Molecular Junction. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2021; 8:e2100055. [PMID: 34145786 PMCID: PMC8292891 DOI: 10.1002/advs.202100055] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 04/25/2021] [Indexed: 05/11/2023]
Abstract
This paper describes the transition from the normal to inverted Marcus region in solid-state tunnel junctions consisting of self-assembled monolayers of benzotetrathiafulvalene (BTTF), and how this transition determines the performance of a molecular diode. Temperature-dependent normalized differential conductance analyses indicate the participation of the HOMO (highest occupied molecular orbital) at large negative bias, which follows typical thermally activated hopping behavior associated with the normal Marcus regime. In contrast, hopping involving the HOMO dominates the mechanism of charge transport at positive bias, yet it is nearly activationless indicating the junction operates in the inverted Marcus region. Thus, within the same junction it is possible to switch between Marcus and inverted Marcus regimes by changing the bias polarity. Consequently, the current only decreases with decreasing temperature at negative bias when hopping is "frozen out," but not at positive bias resulting in a 30-fold increase in the molecular rectification efficiency. These results indicate that the charge transport in the inverted Marcus region is readily accessible in junctions with redox molecules in the weak coupling regime and control over different hopping regimes can be used to improve junction performance.
Collapse
Affiliation(s)
- Yingmei Han
- Department of ChemistryNational University of Singapore3 Science Drive 3Singapore117543Singapore
| | - Cameron Nickle
- Department of PhysicsUniversity of Central FloridaOrlandoFL32816USA
| | - Maria Serena Maglione
- Institut de Ciència de Materials de Barcelona (ICMAB‐CSIC)/CIBER‐BBNCampus de la UABBellaterra08193Spain
| | | | - Javier Casado‐Montenegro
- Institut de Ciència de Materials de Barcelona (ICMAB‐CSIC)/CIBER‐BBNCampus de la UABBellaterra08193Spain
| | - Dong‐Chen Qi
- Centre for Materials ScienceSchool of Chemistry and PhysicsQueensland University of TechnologyBrisbaneQueensland4001Australia
| | - Xiaoping Chen
- Department of ChemistryNational University of Singapore3 Science Drive 3Singapore117543Singapore
| | - Anton Tadich
- Australian Synchrotron ClaytonVictoria3168Australia
| | - Bruce Cowie
- Australian Synchrotron ClaytonVictoria3168Australia
| | - Marta Mas‐Torrent
- Institut de Ciència de Materials de Barcelona (ICMAB‐CSIC)/CIBER‐BBNCampus de la UABBellaterra08193Spain
| | - Concepció Rovira
- Institut de Ciència de Materials de Barcelona (ICMAB‐CSIC)/CIBER‐BBNCampus de la UABBellaterra08193Spain
| | - Jérôme Cornil
- Laboratory for Chemistry of Novel MaterialsUniversity of MonsPlace du Parc 20MonsB‐7000Belgium
| | - Jaume Veciana
- Institut de Ciència de Materials de Barcelona (ICMAB‐CSIC)/CIBER‐BBNCampus de la UABBellaterra08193Spain
| | | | - Christian A. Nijhuis
- Department of ChemistryNational University of Singapore3 Science Drive 3Singapore117543Singapore
- Centre for Advanced 2D Materials and Graphene Research CenterNational University of Singapore6 Science Drive 2Singapore117546Singapore
- Hybrid Materials for Opto‐Electronics GroupDepartment of Molecules and MaterialsMESA+ Institute for Nanotechnology and Center for Brain‐Inspired Nano SystemsFaculty of Science and TechnologyUniversity of TwenteP.O. Box 217EnschedeAE 7500The Netherlands
| |
Collapse
|
22
|
Wei YC, Shen SW, Wu CH, Ho SY, Zhang Z, Wu CI, Chou PT. Through-Space Exciton Delocalization in Segregated HJ-Crystalline Molecular Aggregates. J Phys Chem A 2021; 125:943-953. [DOI: 10.1021/acs.jpca.0c09075] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yu-Chen Wei
- Department of Chemistry, National Taiwan University, Taipei, 10617 Taiwan, ROC
| | - Shin-Wei Shen
- Graduate Institute of Photonics and Optoelectronics, National Taiwan University, Taipei, 10617 Taiwan, ROC
| | - Cheng-Ham Wu
- Department of Chemistry, National Taiwan University, Taipei, 10617 Taiwan, ROC
| | - Ssu-Yu Ho
- Department of Chemistry, National Taiwan University, Taipei, 10617 Taiwan, ROC
| | - Zhiyun Zhang
- Key Laboratory for Advanced Materials and Institute of Fine Chemicals, East China University of Science & Technology, Shanghai 200237, P. R. China
| | - Chih-I Wu
- Graduate Institute of Photonics and Optoelectronics, National Taiwan University, Taipei, 10617 Taiwan, ROC
| | - Pi-Tai Chou
- Department of Chemistry, National Taiwan University, Taipei, 10617 Taiwan, ROC
| |
Collapse
|
23
|
Termine R, Golemme A. Charge Mobility in Discotic Liquid Crystals. Int J Mol Sci 2021; 22:E877. [PMID: 33467214 PMCID: PMC7830985 DOI: 10.3390/ijms22020877] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2020] [Revised: 01/11/2021] [Accepted: 01/12/2021] [Indexed: 12/12/2022] Open
Abstract
Discotic (disk-shaped) molecules or molecular aggregates may form, within a certain temperature range, partially ordered phases, known as discotic liquid crystals, which have been extensively studied in the recent past. On the one hand, this interest was prompted by the fact that they represent models for testing energy and charge transport theories in organic materials. However, their long-range self-assembling properties, potential low cost, ease of processability with a variety of solvents and the relative ease of tailoring their properties via chemical synthesis, drove the attention of researchers also towards the exploitation of their semiconducting properties in organic electronic devices. This review covers recent research on the charge transport properties of discotic mesophases, starting with an introduction to their phase structure, followed by an overview of the models used to describe charge mobility in organic substances in general and in these systems in particular, and by the description of the techniques most commonly used to measure their charge mobility. The reader already familiar or not interested in such details can easily skip these sections and refer to the core section of this work, focusing on the most recent and significant results regarding charge mobility in discotic liquid crystals.
Collapse
Affiliation(s)
- Roberto Termine
- LASCAMM CR-INSTM, CNR-NANOTEC SS di Rende, Dipartimento di Fisica, Università Della Calabria, 87036 Rende, Italy;
| | | |
Collapse
|
24
|
Balzer D, Smolders TJAM, Blyth D, Hood SN, Kassal I. Delocalised kinetic Monte Carlo for simulating delocalisation-enhanced charge and exciton transport in disordered materials. Chem Sci 2020; 12:2276-2285. [PMID: 34163994 PMCID: PMC8179315 DOI: 10.1039/d0sc04116e] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Charge transport is well understood in both highly ordered materials (band conduction) or highly disordered ones (hopping conduction). In moderately disordered materials—including many organic semiconductors—the approximations valid in either extreme break down, making it difficult to accurately model the conduction. In particular, describing wavefunction delocalisation requires a quantum treatment, which is difficult in disordered materials that lack periodicity. Here, we present the first three-dimensional model of partially delocalised charge and exciton transport in materials in the intermediate disorder regime. Our approach is based on polaron-transformed Redfield theory, but overcomes several computational roadblocks by mapping the quantum-mechanical techniques onto kinetic Monte Carlo. Our theory, delocalised kinetic Monte Carlo (dKMC), shows that the fundamental physics of transport in moderately disordered materials is that of charges hopping between partially delocalised electronic states. Our results reveal why standard kinetic Monte Carlo can dramatically underestimate mobilities even in disordered organic semiconductors, where even a little delocalisation can substantially enhance mobilities, as well as showing that three-dimensional calculations capture important delocalisation effects neglected in lower-dimensional approximations. The first three-dimensional model of transport in moderately disordered materials shows that a little delocalisation can dramatically enhance mobilities.![]()
Collapse
Affiliation(s)
- Daniel Balzer
- School of Chemistry and University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| | - Thijs J A M Smolders
- School of Chemistry and University of Sydney Nano Institute, University of Sydney NSW 2006 Australia .,Institute for Molecules and Materials, Radboud University 6525 AJ Nijmegen The Netherlands
| | - David Blyth
- School of Mathematics and Physics, University of Queensland St. Lucia QLD 4072 Australia
| | - Samantha N Hood
- School of Mathematics and Physics, University of Queensland St. Lucia QLD 4072 Australia
| | - Ivan Kassal
- School of Chemistry and University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| |
Collapse
|
25
|
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
|
26
|
Abstract
When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
- Jiabo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
27
|
Gütlein P, Blumberger J, Oberhofer H. An Iterative Fragment Scheme for the ACKS2 Electronic Polarization Model: Application to Molecular Dimers and Chains. J Chem Theory Comput 2020; 16:5723-5735. [PMID: 32701273 DOI: 10.1021/acs.jctc.0c00151] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The treatment of electrostatic interactions is a key ingredient in the force field-based simulation of condensed phase systems. Most approaches used fixed, site-specific point charges. Yet, it is now clear that many applications of force fields (FFs) demand more sophisticated treatments, prompting the implementation of charge equilibration methods in polarizable FFs to allow the redistribution of charge within the system. One approach allowing both, charge redistribution and site-specific polarization, while at the same time solving methodological shortcomings of earlier methods, is the first-principles-derived atom-condensed Kohn-Sham density functional theory method approximated to the second order (ACKS2). In this work, we present two fragment approaches to ACKS2, termed f-ACKS2 and a self-consistent version, scf-ACKS2, that treat condensed phase systems as a collection of electronically polarizable molecular fragments. The fragmentation approach to ACKS2 not only leads to a more transferable and less system-specific collection of electronic response parameters but also opens up the method to large condensed phase systems. We validate the accuracies of f-ACKS2 and scf-ACKS2 by comparing polarization energies and induced dipole moments for a number of charged hydrocarbon dimers against DFT reference calculations. Finally, we also apply both fragmented ACKS2 variants to calculate the polarization energy for electron-hole pair separation along a chain of anthracene molecules and find excellent agreement with reference DFT calculations.
Collapse
Affiliation(s)
- Patrick Gütlein
- 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, London WC1E 6BT, U.K.,Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| |
Collapse
|
28
|
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
|
29
|
Prodhan S, Qiu J, Ricci M, Roscioni OM, Wang L, Beljonne D. Design Rules to Maximize Charge-Carrier Mobility along Conjugated Polymer Chains. J Phys Chem Lett 2020; 11:6519-6525. [PMID: 32692920 DOI: 10.1021/acs.jpclett.0c01793] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The emergence of polymeric materials displaying high charge-carrier mobility values despite poor interchain structural order has spawned a renewal of interest in the identification of structure-property relationships pertaining to the transport of charges along conjugated polymer chains and the subsequent design of optimized architectures. Here, we present the results of intrachain charge transport simulations obtained by applying a robust surface hopping algorithm to a phenomenological Hamiltonian parametrized against first-principles simulations. Conformational effects are shown to provide a clear signature in the temperature-dependent charge-carrier mobility that complies with recent experimental observations. We further contrast against molecular crystals the evolution with electronic bandwidth and electron-phonon interactions of the room-temperature mobility in polymers, showing that intrachain charge-carrier mobility values in excess of 100 cm2/(V s) can be achieved through a proper chemical engineering of the backbones.
Collapse
Affiliation(s)
- Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| | - Jing Qiu
- Center for Chemistry of Novel & High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | | | | | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| |
Collapse
|
30
|
Giannini S, Ziogos OG, Carof A, Ellis M, Blumberger J. Flickering Polarons Extending over Ten Nanometres Mediate Charge Transport in High‐Mobility Organic Crystals. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000093] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Samuele Giannini
- 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
| | - Antoine Carof
- Laboratoire de Physique et Chimie Théoriques, CNRS, UMR No. 7019 Université de Lorraine BP 239 Vandœuvre‐lès‐Nancy Cedex 54506 France
| | - Matthew Ellis
- 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
| |
Collapse
|
31
|
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
|
32
|
Fratini S, Nikolka M, Salleo A, Schweicher G, Sirringhaus H. Charge transport in high-mobility conjugated polymers and molecular semiconductors. NATURE MATERIALS 2020; 19:491-502. [PMID: 32296138 DOI: 10.1038/s41563-020-0647-2] [Citation(s) in RCA: 261] [Impact Index Per Article: 65.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Accepted: 02/20/2020] [Indexed: 06/11/2023]
Abstract
Conjugated polymers and molecular semiconductors are emerging as a viable semiconductor technology in industries such as displays, electronics, renewable energy, sensing and healthcare. A key enabling factor has been significant scientific progress in improving their charge transport properties and carrier mobilities, which has been made possible by a better understanding of the molecular structure-property relationships and the underpinning charge transport physics. Here we aim to present a coherent review of how we understand charge transport in these high-mobility van der Waals bonded semiconductors. Specific questions of interest include estimates for intrinsic limits to the carrier mobilities that might ultimately be achievable; a discussion of the coupling between charge and structural dynamics; the importance of molecular conformations and mesoscale structural features; how the transport physics of conjugated polymers and small molecule semiconductors are related; and how the incorporation of counterions in doped films-as used, for example, in bioelectronics and thermoelectric devices-affects the electronic structure and charge transport properties.
Collapse
Affiliation(s)
| | - Mark Nikolka
- Cavendish Laboratory, University of Cambridge, Cambridge, UK
| | - Alberto Salleo
- Department of Materials Science and Engineering, Stanford University, Stanford, CA, USA
| | | | | |
Collapse
|
33
|
Cheng CY, Campbell JE, Day GM. Evolutionary chemical space exploration for functional materials: computational organic semiconductor discovery. Chem Sci 2020; 11:4922-4933. [PMID: 34122948 PMCID: PMC8159259 DOI: 10.1039/d0sc00554a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2020] [Accepted: 04/21/2020] [Indexed: 11/26/2022] Open
Abstract
Computational methods, including crystal structure and property prediction, have the potential to accelerate the materials discovery process by enabling structure prediction and screening of possible molecular building blocks prior to their synthesis. However, the discovery of new functional molecular materials is still limited by the need to identify promising molecules from a vast chemical space. We describe an evolutionary method which explores a user specified region of chemical space to identify promising molecules, which are subsequently evaluated using crystal structure prediction. We demonstrate the methods for the exploration of aza-substituted pentacenes with the aim of finding small molecule organic semiconductors with high charge carrier mobilities, where the space of possible substitution patterns is too large to exhaustively search using a high throughput approach. The method efficiently explores this large space, typically requiring calculations on only ∼1% of molecules during a search. The results reveal two promising structural motifs: aza-substituted naphtho[1,2-a]anthracenes with reorganisation energies as low as pentacene and a series of pyridazine-based molecules having both low reorganisation energies and high electron affinities.
Collapse
Affiliation(s)
- Chi Y Cheng
- Computational Systems Chemistry, School of Chemistry, University of Southampton Highfield Southampton SO17 1NX UK
| | - Josh E Campbell
- Computational Systems Chemistry, School of Chemistry, University of Southampton Highfield Southampton SO17 1NX UK
| | - Graeme M Day
- Computational Systems Chemistry, School of Chemistry, University of Southampton Highfield Southampton SO17 1NX UK
| |
Collapse
|
34
|
Xie W, Holub D, Kubař T, Elstner M. Performance of Mixed Quantum-Classical Approaches on Modeling the Crossover from Hopping to Bandlike Charge Transport in Organic Semiconductors. J Chem Theory Comput 2020; 16:2071-2084. [PMID: 32176844 DOI: 10.1021/acs.jctc.9b01271] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
In the present study, several mixed quantum-classical (MQC) methods are applied to on-the-fly nonadiabatic molecular dynamics simulations of hole transport in molecular organic semiconductors (OSCs). The tested MQC methods contain the mean-field Ehrenfest (MFE), trajectory surface hopping (TSH) approaches based on Tully's fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH), the latter in the diabatic/adiabatic representation, and a Landau-Zener type trajectory surface hopping (LZSH). We also tested several correction schemes which were proposed to identify trivial crossings and to remove unphysical long-range charge transfers due to decoherence corrections. In addition, several cost-effective approaches for the nuclear velocity adjustment after an energy-allowed/energy-forbidden hop are investigated with respect to detailed balance and internal consistency conditions. To model a broad spectrum of OSCs with different charge transport characteristics, we derived from the anthracene structural model the construction of two additional models by uniformly scaling down the electronic couplings by the factors of 0.1 and 0.5. Anthracene shows a bandlike charge transport mechanism, characterized by slightly delocalized charge carriers 'diffusing' through the crystal. For smaller couplings, the mechanism changes to a hopping type, characteristically differing in the charge delocalization and temperature dependence. The MFE and corrected adiabatic TSH approaches are able to quantitatively reproduce the expected behavior, while the diabatic LZSH method fails for large couplings, as do approaches which are based on the hopping of localized charge between neighboring sites. Moreover, we find that while the hole mobility of the anthracene crystal simulated using the celebrated Marcus theory is in good agreement with the experimental value, its agreement has to be regarded as an accident due to the overestimation of the prefactor in the Marcus rate equation.
Collapse
Affiliation(s)
- Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Daniel Holub
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.,Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| |
Collapse
|
35
|
Krysko ID, Freidzon AY, Bagaturyants AA. Hole hopping in dimers of N,N' di(1-naphthyl)-N,N'-diphenyl-4,4'-diamine (α-NPD): a theoretical study. Phys Chem Chem Phys 2020; 22:3539-3544. [PMID: 31994567 DOI: 10.1039/c9cp06455a] [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/21/2022]
Abstract
Hole-hopping parameters for Marcus-like charge transport, Marcus hole hopping rates, and hole mobilities are calculated for a series of model dimers of a typical hole-transporting material α-NPD using multireference quantum chemistry. The parameters are extracted from the two-state energy profiles built for charge hopping between two states with a hole localized on each of the monomers. The dependence of the hopping integral on the intermolecular arrangement in the dimer is studied. It is shown that at short intermolecular distances strong orbital interactions between molecules cause a drastic increase in the hopping integral and, therefore, in the hopping rate.
Collapse
Affiliation(s)
- Ilya Dmitrievich Krysko
- National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoye Shosse 31, Moscow 115409, Russia. and Federal Research Center "Crystallography and Photonics" Photochemistry Center, Russian Academy of Sciences, ul. Novatorov 7a, Moscow 119421, Russia
| | - Alexandra Yakovlevna Freidzon
- National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoye Shosse 31, Moscow 115409, Russia. and Federal Research Center "Crystallography and Photonics" Photochemistry Center, Russian Academy of Sciences, ul. Novatorov 7a, Moscow 119421, Russia
| | - Alexander Alexandrovich Bagaturyants
- National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoye Shosse 31, Moscow 115409, Russia. and Federal Research Center "Crystallography and Photonics" Photochemistry Center, Russian Academy of Sciences, ul. Novatorov 7a, Moscow 119421, Russia
| |
Collapse
|
36
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
37
|
Hsu CP. Reorganization energies and spectral densities for electron transfer problems in charge transport materials. Phys Chem Chem Phys 2020; 22:21630-21641. [DOI: 10.1039/d0cp02994g] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Various contributions to the outer reorganization energy of an electron transfer system and their theoretical and computational aspects have been discussed.
Collapse
Affiliation(s)
- Chao-Ping Hsu
- 128 Academia Road Section 2
- Institute of Chemistry
- Academia Sinica
- Taipei
- Taiwan
| |
Collapse
|
38
|
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
|
39
|
Zhong YJ, Lan CF, Lin BC, Hu CD, Cheng YC, Hsu CP. The anisotropy and temperature dependence in the mobility of rubrene. ADVANCES IN QUANTUM CHEMISTRY 2020. [DOI: 10.1016/bs.aiq.2020.04.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
40
|
Carof A, Giannini S, Blumberger J. How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm. Phys Chem Chem Phys 2019; 21:26368-26386. [PMID: 31793569 DOI: 10.1039/c9cp04770k] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.
Collapse
Affiliation(s)
- Antoine Carof
- 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.
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK. and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
41
|
Wang CI, Braza MKE, Claudio GC, Nellas RB, Hsu CP. Machine Learning for Predicting Electron Transfer Coupling. J Phys Chem A 2019; 123:7792-7802. [PMID: 31429287 DOI: 10.1021/acs.jpca.9b04256] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electron transfer coupling is a critical factor in determining electron transfer rates. This coupling strength can be sensitive to details in molecular geometries, especially intermolecular configurations. Thus, studying charge transporting behavior with a full first-principle approach demands a large amount of computation resources in quantum chemistry (QC) calculation. To address this issue, we developed a machine learning (ML) approach to evaluate electronic coupling. A prototypical ML model for an ethylene system was built by kernel ridge regression with Coulomb matrix representation. Since the performance of the ML models highly dependent on their building strategies, we systematically investigated the generality of the ML models, the choice of features and target labels. The best ML model trained with 40 000 samples achieved a mean absolute error of 3.5 meV and greater than 98% accuracy in predicting phases. The distance and orientation dependence of electronic coupling was successfully captured. Bypassing QC calculation, the ML model saved 10-104 times the computation cost. With the help of ML, reliable charge transport models and mechanisms can be further developed.
Collapse
Affiliation(s)
- Chun-I Wang
- Institute of Chemistry , Academia Sinica , Taipei 115 , Taiwan
| | - Mac Kevin E Braza
- Institute of Chemistry, College of Science , University of the Philippines Diliman , Quezon City 1101 , Philippines
| | - Gil C Claudio
- Institute of Chemistry, College of Science , University of the Philippines Diliman , Quezon City 1101 , Philippines
| | - Ricky B Nellas
- Institute of Chemistry, College of Science , University of the Philippines Diliman , Quezon City 1101 , Philippines
| | - Chao-Ping Hsu
- Institute of Chemistry , Academia Sinica , Taipei 115 , Taiwan
| |
Collapse
|
42
|
Jiang X, Futera Z, Blumberger J. Ergodicity-Breaking in Thermal Biological Electron Transfer? Cytochrome C Revisited. J Phys Chem B 2019; 123:7588-7598. [PMID: 31405279 DOI: 10.1021/acs.jpcb.9b05253] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
It was recently suggested that certain redox proteins operate in an ergodicity-breaking regime to facilitate biological electron transfer (ET). A signature for this is a large variance reorganization free energy (several electronvolts) but a significantly smaller Stokes reorganization free energy due to incomplete protein relaxation on the time scale of the ET event. Here we investigate whether this picture holds for oxidation of cytochrome c in aqueous solution, at various levels of theory including classical molecular dynamics with two additive and one electronically polarizable force field, and QM/MM calculations with the QM region treated by full electrostatic DFT embedding and by the perturbed matrix method. Sampling the protein and energy gap dynamics over more than 250 ns, we find no evidence for ergodicity-breaking effects. In particular, the inclusion of electronic polarizability of the heme group at QM/MM levels did not induce nonergodic effects, contrary to previous reports by Matyushov et al. The well-known problem of overestimation of reorganization free energies with additive force fields is cured when the protein and solvent are treated as electronically polarizable. Ergodicity-breaking effects may occur in other redox proteins, and our results suggest that long simulations, ideally on the ET time scale, with electronically polarizable force fields are required to obtain strong numerical evidence for them.
Collapse
Affiliation(s)
- Xiuyun Jiang
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Zdenek Futera
- 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
|
43
|
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: 81] [Impact Index Per Article: 16.2] [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
|
44
|
Gütlein P, Lang L, Reuter K, Blumberger J, Oberhofer H. Toward First-Principles-Level Polarization Energies in Force Fields: A Gaussian Basis for the Atom-Condensed Kohn-Sham Method. J Chem Theory Comput 2019; 15:4516-4525. [PMID: 31276382 DOI: 10.1021/acs.jctc.9b00415] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The last 20 years of force field development have shown that even well parametrized classical models need to at least approximate the dielectric response of molecular systems-based, e.g., on atomic polarizabilities-in order to correctly render their structural and dynamic properties. Yet, despite great advances most approaches tend to be based on ad hoc assumptions and often insufficiently capture the dielectric response of the system to external perturbations, such as, e.g., charge carriers in semiconducting materials. A possible remedy was recently introduced with the atom-condensed Kohn-Sham density-functional theory approximated to second order (ACKS2), which is fully derived from first principles. Unfortunately, specifically its reliance on first-principles derived parameters so far precluded the widespread adoption of ACKS2. Opening up ACKS2 for general use, we here present a reformulation of the method in terms of Gaussian basis functions, which allows us to determine many of the ACKS2 parameters analytically. Two sets of parameters depending on exchange-correlation interactions are still calculated numerically, but we show that they could be straightforwardly parametrized owing to the smoothness of the new basis. Our approach exhibits three crucial benefits for future applications in force fields: i) efficiency, ii) accuracy, and iii) transferability. We numerically validate our Gaussian augmented ACKS2 model for a set of small hydrocarbons which shows a very good agreement with density-functional theory reference calculations. To further demonstrate the method's accuracy and transferability for realistic systems, we calculate polarization responses and energies of anthracene and tetracene, two major building blocks in organic semiconductors.
Collapse
Affiliation(s)
- Patrick Gütlein
- Chair for Theoretical Chemistry and Catalysis Research Center , Technische Universität München , Lichtenbergstrasse 4 , D-85747 Garching , Germany
| | - Lucas Lang
- 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 , London WC1E 6BT , U.K.,Institute for Advanced Study , Technische Universität München , Lichtenbergstrasse 2 a , D-85748 Garching , Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center , Technische Universität München , Lichtenbergstrasse 4 , D-85747 Garching , Germany
| |
Collapse
|
45
|
Wang L, Qiu J, Bai X, Xu J. Surface hopping methods for nonadiabatic dynamics in extended systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1435] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Linjun Wang
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jing Qiu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Xin Bai
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jiabo Xu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| |
Collapse
|
46
|
Lian M, Wang YC, Ke Y, Zhao Y. Non-Markovian stochastic Schrödinger equation in k-space toward the calculation of carrier dynamics in organic semiconductors. J Chem Phys 2019; 151:044115. [DOI: 10.1063/1.5096219] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Man Lian
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yaling Ke
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
47
|
Abstract
We present a subspace surface hopping strategy to deal with complex surface crossings in nonadiabatic dynamics. By focusing on only important adiabatic states, we make subspace crossing correction (SCC) in the framework of the standard fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH). The resulting SCC-FSSH and SCC-GFSH approaches show much better performance than the counterparts using all adiabatic states for surface hopping. As demonstrated in a series of Holstein models with up to over 1000 molecular sites, both SCC-FSSH and SCC-GFSH show excellent size independence with a large time step size of 1 fs. Especially, SCC-GFSH does not refer to nonadiabatic couplings at all and gives a more proper description of superexchange, and thus, it is promising for realistic applications with complex potential energy surfaces.
Collapse
Affiliation(s)
- Jing Qiu
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Xin Bai
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials , and Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| |
Collapse
|
48
|
Gryn'ova G, Lin KH, Corminboeuf C. Read between the Molecules: Computational Insights into Organic Semiconductors. J Am Chem Soc 2018; 140:16370-16386. [PMID: 30395466 PMCID: PMC6287891 DOI: 10.1021/jacs.8b07985] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
The
performance and key electronic properties of molecular organic
semiconductors are dictated by the interplay between the chemistry
of the molecular core and the intermolecular factors of which manipulation
has inspired both experimentalists and theorists. This Perspective
presents major computational challenges and modern methodological
strategies to advance the field. The discussion ranges from insights
and design principles at the quantum chemical level, in-depth atomistic
modeling based on multiscale protocols, morphological prediction and
characterization as well as energy-property maps involving data-driven
analysis. A personal overview of the past achievements and future
direction is also provided.
Collapse
Affiliation(s)
- Ganna Gryn'ova
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Kun-Han Lin
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland.,Laboratory for Computational Molecular Design and National Center for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland.,Laboratory for Computational Molecular Design and National Center for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| |
Collapse
|
49
|
Rice B, Guilbert AAY, Frost JM, Nelson J. Polaron States in Fullerene Adducts Modeled by Coarse-Grained Molecular Dynamics and Tight Binding. J Phys Chem Lett 2018; 9:6616-6623. [PMID: 30380880 DOI: 10.1021/acs.jpclett.8b02320] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Strong electron-phonon coupling leads to polaron localization in molecular semiconductor materials and influences charge transport, but it is expensive to calculate atomistically. Here, we propose a simple and efficient model to determine the energy and spatial extent of polaron states within a coarse-grained representation of a disordered molecular film. We calculate the electronic structure of the molecular assembly using a tight-binding Hamiltonian and determine the polaron state self-consistently by perturbing the site energies by the dielectric response of the surrounding medium to the charge. When applied to fullerene derivatives, the method shows that polarons extend over multiple molecules in C60 but localize on single molecules in higher adducts of phenyl-C61-butyric-acid-methyl-ester (PCBM) because of packing disorder and the polar side chains. In PCBM, polarons localize on single molecules only when energetic disorder is included or when the fullerene is dispersed in a blend. The method helps to establish the conditions under which a hopping transport model is justified.
Collapse
Affiliation(s)
- Beth Rice
- Department of Physics , Imperial College London , London SW7 2BZ , U.K
| | - Anne A Y Guilbert
- Department of Physics , Imperial College London , London SW7 2BZ , U.K
| | - Jarvist M Frost
- Department of Physics , Imperial College London , London SW7 2BZ , U.K
- Department of Physics , King's College London , London WC2R 2LS , U.K
| | - Jenny Nelson
- Department of Physics , Imperial College London , London SW7 2BZ , U.K
| |
Collapse
|
50
|
Qiu J, Bai X, Wang L. Crossing Classified and Corrected Fewest Switches Surface Hopping. J Phys Chem Lett 2018; 9:4319-4325. [PMID: 30011207 DOI: 10.1021/acs.jpclett.8b01902] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
In the traditional fewest switches surface hopping (FSSH), trivial crossings between uncoupled or weakly coupled states have highly peaked nonadiabatic couplings and thus are difficult to deal with in the preferred, adiabatic representation. Here, we classify surface crossings into four general types and propose a parameter-free crossing corrected FSSH (CC-FSSH) algorithm, which could treat multiple trivial crossings within a time interval. As examples, Holstein Hamiltonians with different parameters are adopted to mimic electron dynamics in tens to hundreds of molecules, which suffer from severe trivial crossing problems. Using existed surface hopping approaches as references, we show that CC-FSSH exhibits significantly fast time interval convergence and weak system size dependence. In all cases, a reliable description is achieved with a large time interval of 1 fs. With a simple formalism and the ability to describe complex surface crossings, CC-FSSH could potentially simulate general nonadiabatic dynamics in nanoscale materials with a high efficiency.
Collapse
Affiliation(s)
- Jing Qiu
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Xin Bai
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| | - Linjun Wang
- Department of Chemistry , Zhejiang University , Hangzhou 310027 , China
| |
Collapse
|