1
|
Liu D, Wang B, Vasenko AS, Prezhdo OV. Decoherence ensures convergence of non-adiabatic molecular dynamics with number of states. J Chem Phys 2024; 161:064104. [PMID: 39120030 DOI: 10.1063/5.0222557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Accepted: 07/24/2024] [Indexed: 08/10/2024] Open
Abstract
Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems.
Collapse
Affiliation(s)
| | - Bipeng Wang
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Andrey S Vasenko
- HSE University, 101000 Moscow, Russia
- Donostia International Physics Center (DIPC), 20018 San Sebastián-Donostia, Euskadi, Spain
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
2
|
Ghalami F, Dohmen PM, Krämer M, Elstner M, Xie W. Nonadiabatic Simulation of Exciton Dynamics in Organic Semiconductors Using Neural Network-Based Frenkel Hamiltonian and Gradients. J Chem Theory Comput 2024; 20:6160-6174. [PMID: 38976696 DOI: 10.1021/acs.jctc.4c00220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
In this study, we present a multiscale method to simulate the propagation of Frenkel singlet excitons in organic semiconductors (OSCs). The approach uses neural network models to train a Frenkel-type Hamiltonian and its gradient, obtained by the long-range correction version of density functional tight-binding with self-consistent charges. Our models accurately predict site energies, excitonic couplings, and corresponding gradients, essential for the nonadiabatic molecular dynamics simulations. Combined with the fewest switches surface hopping algorithm, the method was applied to four representative OSCs: anthracene, pentacene, perylenediimide, and diindenoperylene. The simulated exciton diffusion constants align well with experimental and reported theoretical values and offer valuable insights into exciton dynamics in OSCs.
Collapse
Affiliation(s)
- Farhad Ghalami
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
- Institute of Nano Technology (INT), Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Philipp M Dohmen
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Mila Krämer
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
- Institute of Nano Technology (INT), Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Weiwei Xie
- Frontiers Science Center for New Organic Matter, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), State Key Laboratory of Advanced Chemical Power Sources, College of Chemistry, Nankai University, Tianjin 300071, China
| |
Collapse
|
3
|
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
|
4
|
Bhat V, Ganapathysubramanian B, Risko C. Rapid Estimation of the Intermolecular Electronic Couplings and Charge-Carrier Mobilities of Crystalline Molecular Organic Semiconductors through a Machine Learning Pipeline. J Phys Chem Lett 2024; 15:7206-7213. [PMID: 38973725 DOI: 10.1021/acs.jpclett.4c01309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Organic semiconductors (OSC) offer tremendous potential across a wide range of (opto)electronic applications. OSC development, however, is often limited by trial-and-error design, with computational modeling approaches deployed to evaluate and screen candidates through a suite of molecular and materials descriptors that generally require hours to days of computational time to accumulate. Such bottlenecks slow the pace and limit the exploration of the vast chemical space comprising OSC. When considering charge-carrier transport in OSC, a key parameter of interest is the intermolecular electronic coupling. Here, we introduce a machine learning (ML) model to predict intermolecular electronic couplings in organic crystalline materials from their three-dimensional (3D) molecular geometries. The ML predictions take only a few seconds of computing time compared to hours by density functional theory (DFT) methods. To demonstrate the utility of the ML predictions, we deploy the ML model in conjunction with mathematical formulations to rapidly screen the charge-carrier mobility anisotropy for more than 60,000 molecular crystal structures and compare the ML predictions to DFT benchmarks.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| | - Baskar Ganapathysubramanian
- Department of Mechanical Engineering & Translational AI Center, Iowa State University, Ames, Iowa 50010, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
5
|
Sokolov M, Hoffmann DS, Dohmen PM, Krämer M, Höfener S, Kleinekathöfer U, Elstner M. Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex. Phys Chem Chem Phys 2024; 26:19469-19496. [PMID: 38979564 DOI: 10.1039/d4cp02116a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.
Collapse
Affiliation(s)
- Monja Sokolov
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - David S Hoffmann
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - Philipp M Dohmen
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Mila Krämer
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | - Sebastian Höfener
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | | | - Marcus Elstner
- Institute of Physical Chemistry (IPC), Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| |
Collapse
|
6
|
Zhang Q, Shao X, Li W, Mi W, Pavanello M, Akimov AV. Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:385901. [PMID: 38866023 DOI: 10.1088/1361-648x/ad577d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 06/12/2024] [Indexed: 06/14/2024]
Abstract
In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.
Collapse
Affiliation(s)
- Qingxin Zhang
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| | - Xuecheng Shao
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China
| | - Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, People's Republic of China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People's Republic of China
| | - Michele Pavanello
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| |
Collapse
|
7
|
Pavalamuthu M, Navamani K. Entropy-ruled nonequilibrium charge transport in thiazolothiazole-based molecular crystals: a quantum chemical study. Phys Chem Chem Phys 2024; 26:16488-16504. [PMID: 38751327 DOI: 10.1039/d3cp05739a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
The charge and energy fluctuations in molecular solids are crucial factors for a better understanding of charge transport (CT) in organic semiconductors. The energetic disorder-coupled molecular charge transport is still not well-established. Moreover, the conventional Einstein's diffusion (D)-mobility (μ) relation fails to explain the quantum features of organic semiconductors, including nonequilibrium and degenerate transport systems, where kB is the Boltzmann constant, T is the temperature and q is the electric charge. To overcome this issue, a unified version of the entropy-ruled D/μ relation was proposed by Navamani (J. Phys. Chem. Lett., 2024, 15, 2519-2528) for hopping and band transport systems as where d, η and heff are the dimension (d = 1, 2, 3), chemical potential and effective entropy, respectively. Within this context, we investigate the CT properties of 2,5-bis(4-methoxyphenyl)thiazolo[5,4-d]thiazole (MOP-TZTZ) and 2,5-bis(2,4,5 trifluorophenyl)-thiazolo[5,4-d]thiazole (TFP-TZTZ) molecular solids using electronic structure calculations and the entropy-ruled method. The CT key parameters such as charge transfer integral and site energy are computed by matrix elements of the Kohn-Sham Hamiltonian. Using Marcus theory, the charge transfer rate is numerically calculated for MOP-TZTZ and TFP-TZTZ molecular crystals under different site energy disorder (ΔEij(E⃑)) situations. Using our entropy-ruled method, the exact diffusion-mobility (D/μ) and other transport quantities such as thermodynamic density of states, conductivity, and current density are calculated for these derivatives at different applied electric field values via the site energy disorder. The theoretical results show that the molecule TFP-TZTZ has good hole mobility (∼0.012 cm2 V-1 s-1) at a site energy disorder value of 90 meV. The obtained ideality factor from the Navamani-Shockley diode current density equation categorizes the typical transport as either the Langevin-type or Shockley-Read-Hall mechanism in the studied molecular solids. Our analysis clearly shows that both the electron and hole transport in these MOP-TZTZ and TFP-TZTZ molecules follow the trap-free Langevin mechanism, which is indeed ideal for designing charge-transporting molecular devices.
Collapse
Affiliation(s)
- M Pavalamuthu
- Department of Physics, Centre for Research and Development (CFRD), KPR Institute of Engineering and Technology, Coimbatore-641407, India.
| | - K Navamani
- Department of Physics, Centre for Research and Development (CFRD), KPR Institute of Engineering and Technology, Coimbatore-641407, India.
| |
Collapse
|
8
|
Liu Z, Song Z, Sun X. All-Atom Photoinduced Charge Transfer Dynamics in Condensed Phase via Multistate Nonlinear-Response Instantaneous Marcus Theory. J Chem Theory Comput 2024; 20:3993-4006. [PMID: 38657208 PMCID: PMC11099976 DOI: 10.1021/acs.jctc.4c00010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/30/2024] [Accepted: 04/08/2024] [Indexed: 04/26/2024]
Abstract
Photoinduced charge transfer (CT) in the condensed phase is an essential component in solar energy conversion, but it is challenging to simulate such a process on the all-atom level. The traditional Marcus theory has been utilized for obtaining CT rate constants between pairs of electronic states but cannot account for the nonequilibrium effects due to the initial nuclear preparation. The recently proposed instantaneous Marcus theory (IMT) and its nonlinear-response formulation allow for incorporating the nonequilibrium nuclear relaxation to electronic transition between two states after the photoexcitation from the equilibrium ground state and provide the time-dependent rate coefficient. In this work, we extend the nonlinear-response IMT method for treating photoinduced CT among general multiple electronic states and demonstrate it in the organic photovoltaic carotenoid-porphyrin-fullerene triad dissolved in explicit tetrahydrofuran solvent. All-atom molecular dynamics simulations were employed to obtain the time correlation functions of energy gaps, which were used to generate the IMT-required time-dependent averages and variances of the relevant energy gaps. Our calculations show that the multistate IMT could capture the significant nonequilibrium effects due to the initial nuclear state preparation, and this is corroborated by the substantial differences between the population dynamics predicted by the multistate IMT and the Marcus theory, where the Marcus theory underestimates the population transfer. The population dynamics by multistate IMT is also shown to have a better agreement with the all-atom nonadiabatic mapping dynamics than the Marcus theory does. Because the multistate nonlinear-response IMT is straightforward and cost-effective in implementation and accounts for the nonequilibrium nuclear effects, we believe this method offers a practical strategy for studying charge transfer dynamics in complex condensed-phase systems.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Zailing Song
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
9
|
Lin HH, Wang CI, Yang CH, Secario MK, Hsu CP. Two-Step Machine Learning Approach for Charge-Transfer Coupling with Structurally Diverse Data. J Phys Chem A 2024; 128:271-280. [PMID: 38157315 DOI: 10.1021/acs.jpca.3c04524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
Electronic coupling is important in determining charge-transfer rates and dynamics. Coupling strength is sensitive to both intermolecular, e.g., orientation or distance, and intramolecular degrees of freedom. Hence, it is challenging to build an accurate machine learning model to predict electronic coupling of molecular pairs, especially for those derived from the amorphous phase, for which intermolecular configurations are much more diverse than those derived from crystals. In this work, we devise a new prediction algorithm that employs two consecutive KRR models. The first model predicts molecular orbitals (MOs) from structural variation for each fragment, and coupling is further predicted by using the overlap integral included in a second model. With our two-step procedure, we achieved mean absolute errors of 0.27 meV for an ethylene dimer and 1.99 meV for a naphthalene pair, much improved accuracy amounting to 14-fold and 3-fold error reductions, respectively. In addition, MOs from the first model can also be the starting point to obtain other quantum chemical properties from atomistic structures. This approach is also compatible with a MO predictor with sufficient accuracy.
Collapse
Affiliation(s)
- Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Molecular Science and Digital Innovation Center, Genetics Generation Advancement Corp, No. 28, Ln. 36, Xinhu First Rd., Neihu, Taipei 114, Taiwan
| | - Chun-I Wang
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Chou-Hsun Yang
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
| | - Muhammad Khari Secario
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Taiwan International Graduate Program on Sustainable Chemical Science & Technology, Academia Sinica Institute of Chemistry, 128 Academia Road Sec.2, Nankang, Taipei 115, Taiwan
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 University Road, Hsinchu 300, Taiwan
| | - Chao-Ping Hsu
- Institute of Chemistry, Academia Sinica, 128 Section 2 Academia Road, Nankang, Taipei 115, Taiwan
- Division of Physics, National Center for Theoretical Sciences, 1, Section 4, Roosevelt Road, Taipei 106, Taiwan
| |
Collapse
|
10
|
Carey RL, Giannini S, Schott S, Lemaur V, Xiao M, Prodhan S, Wang L, Bovoloni M, Quarti C, Beljonne D, Sirringhaus H. Spin relaxation of electron and hole polarons in ambipolar conjugated polymers. Nat Commun 2024; 15:288. [PMID: 38177094 PMCID: PMC10767019 DOI: 10.1038/s41467-023-43505-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Accepted: 11/09/2023] [Indexed: 01/06/2024] Open
Abstract
The charge-transport properties of conjugated polymers have been studied extensively for opto-electronic device applications. Some polymer semiconductors not only support the ambipolar transport of electrons and holes, but do so with comparable carrier mobilities. This opens the possibility of gaining deeper insight into the charge-transport physics of these complex materials via comparison between electron and hole dynamics while keeping other factors, such as polymer microstructure, equal. Here, we use field-induced electron spin resonance spectroscopy to compare the spin relaxation behavior of electron and hole polarons in three ambipolar conjugated polymers. Our experiments show unique relaxation regimes as a function of temperature for electrons and holes, whereby at lower temperatures electrons relax slower than holes, but at higher temperatures, in the so-called spin-shuttling regime, the trend is reversed. On the basis of theoretical simulations, we attribute this to differences in the delocalization of electron and hole wavefunctions and show that spin relaxation in the spin shuttling regimes provides a sensitive probe of the intimate coupling between charge and structural dynamics.
Collapse
Affiliation(s)
- Remington L Carey
- Cavendish Laboratory, University of Cambridge, Cambridge, CB3 0HE, UK
| | - Samuele Giannini
- Laboratory for Chemistry of Novel Materials, University of Mons, 7000, Mons, Belgium
- Institute of Chemistry of OrganoMetallic Compounds, National Research Council (ICCOM-CNR), I-56124, Pisa, Italy
| | - Sam Schott
- Cavendish Laboratory, University of Cambridge, Cambridge, CB3 0HE, UK
| | - Vincent Lemaur
- Laboratory for Chemistry of Novel Materials, University of Mons, 7000, Mons, Belgium
| | - Mingfei Xiao
- Cavendish Laboratory, University of Cambridge, Cambridge, CB3 0HE, UK
| | - Suryoday Prodhan
- Department of Chemistry, University of Liverpool, Liverpool, L69 3BX, UK
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou, 310058, China
| | - Michelangelo Bovoloni
- Laboratory for Chemistry of Novel Materials, University of Mons, 7000, Mons, Belgium
| | - Claudio Quarti
- Laboratory for Chemistry of Novel Materials, University of Mons, 7000, Mons, Belgium
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, 7000, Mons, Belgium
| | | |
Collapse
|
11
|
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
|
12
|
Dines A, Ellis M, Blumberger J. Stabilized coupled trajectory mixed quantum-classical algorithm with improved energy conservation: CTMQC-EDI. J Chem Phys 2023; 159:234118. [PMID: 38117021 DOI: 10.1063/5.0183589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 11/26/2023] [Indexed: 12/21/2023] Open
Abstract
Coupled trajectory mixed quantum-classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularization procedures in certain circumstances. Such procedures in the latter can cause instabilities, leading to undesirable effects, such as energy drift and spurious population transfer, which is expected to become increasingly prevalent when the system gets larger as such events would happen more frequently. We propose a further modification to CTMQC-E, which includes a redefinition of the quantum momentum, CTMQC-EDI (double intercept), such that it has no formal divergences. We then show for Tully models I-III and the double arch model that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI, therefore, shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.
Collapse
Affiliation(s)
- Aaron Dines
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Matthew Ellis
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
13
|
Mostaghimi M, Pacheco Hernandez H, Jiang Y, Wenzel W, Heinke L, Kozlowska M. On-off conduction photoswitching in modelled spiropyran-based metal-organic frameworks. Commun Chem 2023; 6:275. [PMID: 38110545 PMCID: PMC10728195 DOI: 10.1038/s42004-023-01072-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 11/22/2023] [Indexed: 12/20/2023] Open
Abstract
Materials with photoswitchable electronic properties and conductance values that can be reversibly changed over many orders of magnitude are highly desirable. Metal-organic framework (MOF) films functionalized with photoresponsive spiropyran molecules demonstrated the general possibility to switch the conduction by light with potentially large on-off-ratios. However, the fabrication of MOF materials in a trial-and-error approach is cumbersome and would benefit significantly from in silico molecular design. Based on the previous proof-of-principle investigation, here, we design photoswitchable MOFs which incorporate spiropyran photoswitches at controlled positions with defined intermolecular distances and orientations. Using multiscale modelling and automated workflow protocols, four MOF candidates are characterized and their potential for photoswitching the conductivity is explored. Using ab initio calculations of the electronic coupling between the molecules in the MOF, we show that lattice distances and vibrational flexibility tremendously modulate the possible conduction photoswitching between spiropyran- and merocyanine-based MOFs upon light absorption, resulting in average on-off ratios higher than 530 and 4200 for p- and n-conduction switching, respectively. Further functionalization of the photoswitches with electron-donating/-withdrawing groups is demonstrated to shift the energy levels of the frontier orbitals, permitting a guided design of new spiropyran-based photoswitches towards controlled modification between electron and hole conduction in a MOF.
Collapse
Affiliation(s)
- Mersad Mostaghimi
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Helmy Pacheco Hernandez
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Yunzhe Jiang
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Lars Heinke
- Institute of Functional Interfaces (IFG), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany.
| | - Mariana Kozlowska
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131, Karlsruhe, Germany.
| |
Collapse
|
14
|
Lei Y, Xie X, Ma H, Ma J. Vitality of Intralayer Vibration in hBN for Effective Long-Range Interlayer Hole Transfer across High Barriers in MoSe 2/hBN/WSe 2 Heterostructures. J Phys Chem Lett 2023:11190-11199. [PMID: 38055859 DOI: 10.1021/acs.jpclett.3c03040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
Introducing the two-dimensional (2D) hexagonal boron nitride (hBN) between 2D transition metal dichalcogenide (TMD) layers promises convenient manipulation of the interlayer exciton (IX) and interlayer charge transfer in TMD/hBN/TMD heterostructures, while the role of inserted hBN layers during IX formation is controversial. Employing ab initio nonadiabatic molecular dynamics (NAMD) simulations and the electron-phonon coupling model, we systematically investigate interlayer hole transfer in MoSe2/WSe2 bilayers intercalated by hBN layers with various thicknesses. The conventional direct hole transfer from MoSe2 to WSe2 is decelerated by 2-3 orders of magnitude after the hBN insertion. Meanwhile, a novel channel intermediated by a deeper hole of WSe2 becomes dominant, where the intralayer shear mode of hBN plays a crucial role by reducing the energy barriers for this new channel. The unique role of hBN layers is revealed for the first time, enriching the knowledge of the underlying microscopic mechanisms and providing instructive guidance to practical van der Waals optoelectronic devices.
Collapse
Affiliation(s)
- Yuli Lei
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xiaoyu Xie
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
15
|
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
|
16
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
17
|
Kohn JT, Gildemeister N, Grimme S, Fazzi D, Hansen A. Efficient calculation of electronic coupling integrals with the dimer projection method via a density matrix tight-binding potential. J Chem Phys 2023; 159:144106. [PMID: 37818996 DOI: 10.1063/5.0167484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Accepted: 09/04/2023] [Indexed: 10/13/2023] Open
Abstract
Designing organic semiconductors for practical applications in organic solar cells, organic field-effect transistors, and organic light-emitting diodes requires understanding charge transfer mechanisms across different length and time scales. The underlying electron transfer mechanisms can be efficiently explored using semiempirical quantum mechanical (SQM) methods. The dimer projection (DIPRO) method combined with the recently introduced non-self-consistent density matrix tight-binding potential (PTB) [Grimme et al., J. Chem. Phys. 158, 124111 (2023)] is used in this study to evaluate charge transfer integrals important for understanding charge transport mechanisms. PTB, parameterized for the entire Periodic Table up to Z = 86, incorporates approximate non-local exchange, allowing for efficient and accurate calculations for large hetero-organic compounds. Benchmarking against established databases, such as Blumberger's HAB sets, or our newly introduced JAB69 set and comparing with high-level reference data from ωB97X-D4 calculations confirm that DIPRO@PTB consistently performs well among the tested SQM approaches for calculating coupling integrals. DIPRO@PTB yields reasonably accurate results at low computational cost, making it suitable for screening purposes and applications to large systems, such as metal-organic frameworks and cyanine-based molecular aggregates further discussed in this work.
Collapse
Affiliation(s)
- J T Kohn
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - N Gildemeister
- Department of Chemistry, Greinstrasse 4-6, 50939 Köln, Germany
| | - S Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - D Fazzi
- Dipartimento di Chimica "Giacomo Ciamician," Via Selmi 2, 40126 Bologna, Italy
| | - A Hansen
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| |
Collapse
|
18
|
Zhang J, Peng K, Xu ZK, Xiong Y, Liu J, Cai C, Huang X. A comprehensive review on the behavior and evolution of oil droplets during oil/water separation by membranes. Adv Colloid Interface Sci 2023; 319:102971. [PMID: 37562248 DOI: 10.1016/j.cis.2023.102971] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 07/01/2023] [Accepted: 07/31/2023] [Indexed: 08/12/2023]
Abstract
Membrane separation technology has significant advantages for treating oil-in-water emulsions. Understanding the evolution of oil droplets could reveal the interfacial and colloidal interactions, facilitate the design of advanced membranes, and improve the separation performances. This review on the characteristic behavior and evolution of oil droplets focuses on the advanced analytical techniques, and the subsequent fouling as well as demulsification effects during membrane separation. A detailed introduction is provided on microscopic observations and numerical simulations of the dynamic evolution of oil droplets, featuring real-time in-situ visualization and accurate reconstruction, respectively. Characteristic behaviors of these oil droplets include attachment, pinning, wetting, spreading, blockage, intrusion, coalescence, and detachment, which have been quantified by specific proposed parameters and criteria. The fouling process can be evaluated using Hermia and resistance models. The related adhesion force and intrusion pressure as well as droplet-droplet/membrane interfacial interactions can be accurately quantified using various force analysis methods and advanced force measurement techniques. It is encouraging to note that oil coalescence has been achieved through various effects such as electrostatic interactions, mechanical actions, Laplace pressure/surface free energy gradients, and synergistic effects on functional membranes. When oil droplets become destabilized and coalesce into larger ones, the functional membranes can overcome the limitations of size-sieving effect to attain higher separation efficiency. This not only bypasses the trade-off between permeability and rejection, but also significantly reduces membrane fouling. Finally, the challenges and potential research directions in membrane separation are proposed. We hope this review will support the engineering of advanced materials for oil/water separation and research on interface science in general.
Collapse
Affiliation(s)
- Jialu Zhang
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China
| | - Kaiming Peng
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China; Institute of Carbon Neutrality, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China.
| | - Zhi-Kang Xu
- MOE Key Laboratory of Macromolecular Synthesis and Functionalization, and Key Lab of Adsorption and Separation Materials & Technologies of Zhejiang Province, Department of Polymer Science and Engineering, Zhejiang University, No.38 Zheda Road, Hangzhou 310027, PR China
| | - Yongjiao Xiong
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China
| | - Jia Liu
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China; Institute of Carbon Neutrality, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China
| | - Chen Cai
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China; Institute of Carbon Neutrality, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China
| | - Xiangfeng Huang
- State Key Laboratory of Pollution Control and Resource Reuse, Ministry of Education Key Laboratory of Yangtze River Water Environment, Shanghai Institute of Pollution Control and Ecological Security, College of Environmental Science and Engineering, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China; Institute of Carbon Neutrality, Tongji University, No.1239 Siping Road, Shanghai 200092, PR China.
| |
Collapse
|
19
|
Shao C, Shi Z, Xu J, Wang L. Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics. J Phys Chem Lett 2023; 14:7680-7689. [PMID: 37606199 DOI: 10.1021/acs.jpclett.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
20
|
Gumber S, Prezhdo OV. Zeno and Anti-Zeno Effects in Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:7274-7282. [PMID: 37556319 PMCID: PMC10440816 DOI: 10.1021/acs.jpclett.3c01831] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Accepted: 08/01/2023] [Indexed: 08/11/2023]
Abstract
Decoherence plays an important role in nonadiabatic (NA) molecular dynamics (MD) simulations because it provides a physical mechanism for trajectory hopping and can alter transition rates by orders of magnitude. Generally, decoherence effects slow quantum transitions, as exemplified by the quantum Zeno effect: in the limit of infinitely fast decoherence, the transitions stop. If the measurements are not sufficiently frequent, an opposite quantum anti-Zeno effect occurs, in which the transitions are accelerated with faster decoherence. Using two common NA-MD approaches, fewest switches surface hopping and decoherence-induced surface hopping, combined with analytic examination, we demonstrate that including decoherence into NA-MD slows down NA transitions; however, many realistic systems operate in the anti-Zeno regime. Therefore, it is important that NA-MD methods describe both Zeno and anti-Zeno effects. Numerical simulations of charge trapping and relaxation in graphitic carbon nitride suggest that time-dependent NA Hamiltonians encountered in realistic systems produce robust results with respect to errors in the decoherence time, a favorable feature for NA-MD simulations.
Collapse
Affiliation(s)
- Shriya Gumber
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department
of Physics and Astronomy, University of
Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
21
|
Bhat V, Callaway CP, Risko C. Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials. Chem Rev 2023. [PMID: 37141497 DOI: 10.1021/acs.chemrev.2c00704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.
Collapse
Affiliation(s)
- Vinayak Bhat
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Connor P Callaway
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| | - Chad Risko
- Department of Chemistry & Center for Applied Energy Research, University of Kentucky, Lexington, Kentucky 40506-0055, United States
| |
Collapse
|
22
|
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
|
23
|
Ambrosio F, Wiktor J, Landi A, Peluso A. Charge Localization in Acene Crystals from Ab Initio Electronic Structure. J Phys Chem Lett 2023; 14:3343-3351. [PMID: 36994951 PMCID: PMC10084468 DOI: 10.1021/acs.jpclett.3c00191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Accepted: 03/10/2023] [Indexed: 06/19/2023]
Abstract
The performance of Koopmans-compliant hybrid functionals in reproducing the electronic structure of organic crystals is tested for a series of acene crystals. The calculated band gaps are found to be consistent with those achieved with the GW method at a fraction of the computational cost and in excellent accord with the experimental results at room temperature, when including the thermal renormalization. The energetics of excess holes and electrons reveals a struggle between polaronic localization and band-like delocalization. The consequences of these results on the transport properties of acene crystals are discussed.
Collapse
Affiliation(s)
- Francesco Ambrosio
- Dipartimento
di Chimica e Biologia Adolfo Zambelli, Università
di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
- Dipartimento
di Scienze, Università degli Studi
della Basilicata, Viale
dell’Ateneo Lucano, 10-85100 Potenza, Italy
| | - Julia Wiktor
- Department
of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Alessandro Landi
- Dipartimento
di Chimica e Biologia Adolfo Zambelli, Università
di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
| | - Andrea Peluso
- Dipartimento
di Chimica e Biologia Adolfo Zambelli, Università
di Salerno, Via Giovanni Paolo II, I-84084 Fisciano (SA), Italy
| |
Collapse
|
24
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, France
| |
Collapse
|
25
|
Abstract
With the advent of a new era of smart-technology, the demand for more economic optoelectronic materials that do not compromise with efficiency is gradually on the rise. Organic semiconductors provide greener alternatives to the conventional inorganic ones, but encounter the challenge of balancing charge carrier mobility with processability in devices. Discotic liquid crystals (DLCs), a class of self-assembling soft organic materials, possess the perfect degree of order and dynamics to address this challenge. Providing unidimensional charge carrier pathways through their nanoscale columnar architecture, DLCs can behave as efficient charge transport systems across a wide range of optoelectronic devices. Moreover, DLCs are solution-processable, thus reducing the fabrication cost. In this article, we have discussed the approaches towards developing DLCs as semiconductors, focusing on their molecular design concepts, supramolecular structures and electronic properties in the context of their charge carrier mobilities.
Collapse
Affiliation(s)
- Ritobrata De
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector-81, SAS Nagar, Knowledge City, Manauli-140306, India.
| | - Santanu Kumar Pal
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Sector-81, SAS Nagar, Knowledge City, Manauli-140306, India.
| |
Collapse
|
26
|
Knepp ZJ, Masso GB, Fredin LA. Efficiently predicting directional carrier mobilities in organic materials with the Boltzmann transport equation. J Chem Phys 2023; 158:064704. [PMID: 36792516 DOI: 10.1063/5.0128125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Describing charge carrier anisotropy in crystalline organic semiconductors with ab initio methods is challenging because of the weak intermolecular interactions that lead to both localized and delocalized charge transport mechanisms. Small polaron hopping models (localized) are generally used to describe materials with small charge carrier mobilities, while periodic band models (delocalized) are used to describe materials with high charge carrier mobilities. Here, we prove the advantage of applying the constant relaxation time approximation of the Boltzmann transport equation (BTE) to efficiently predict the anisotropic hole mobilities of several unsubstituted (anthracene, tetracene, pentacene, and hexacene) and substituted (2,6-diphenylanthracene, rubrene, and TIPS-pentacene) high-mobility n-acene single crystals. Several density functionals are used to optimize the crystals, and the composite density functional PBEsol0-3c/sol-def2-mSVP predicts the most experimentally similar geometries, adequate indirect bandgaps, and the theoretically consistent n-acene charge transport mobility trend. Similarities between BTE and Marcus mobilities are presented for each crystal. BTE and Marcus charge carrier mobilities computed at the same geometry result in similar mobility trends, differing mostly in materials with more substitutions or structurally complex substituents. By using a reduced number of calculations, BTE is able to predict anisotropic carrier mobilities efficiently and effectively for a range of high-mobility organic semiconductors.
Collapse
Affiliation(s)
- Zachary J Knepp
- Department of Chemistry, Lehigh University, 6 E. Packer Ave., Bethlehem, Pennsylvania 18015, USA
| | - Gabriel B Masso
- Department of Chemistry, Lehigh University, 6 E. Packer Ave., Bethlehem, Pennsylvania 18015, USA
| | - Lisa A Fredin
- Department of Chemistry, Lehigh University, 6 E. Packer Ave., Bethlehem, Pennsylvania 18015, USA
| |
Collapse
|
27
|
Zhu Y, Jiang Q, Zhang J, Ma Y. Recent Progress of Organic Semiconductor Materials in Spintronics. Chem Asian J 2023; 18:e202201125. [PMID: 36510771 DOI: 10.1002/asia.202201125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 12/12/2022] [Accepted: 12/12/2022] [Indexed: 12/15/2022]
Abstract
Spintronics, a new discipline focusing on the spin-dependent transport process of electrons, has been developing rapidly. Spin valves are the most significant carriers of spintronics utilizing the spin freedom of electrons. It is expected to pierce "Moore's Law" and become the core component in processors of the next generation. Organic semiconductors advance in their adjustable band gap, weak spin-orbit coupling and hyperfine interaction, excellent film-forming property, having enormous promise for spin valves. Here, the principle of spin valves is introduced, and the history and progress in organic spin injection and transport materials are summarized. Then we analyze the influence of spinterface on device performance and introduce reliable methods of constructing organic spin valves. Finally, the challenges for spin valves are discussed, and the future is proposed. We aim to draw the attention of researchers to organic spin valves and promote further research in spintronics through this paper.
Collapse
Affiliation(s)
- Yanuo Zhu
- Institute of Polymer Optoelectronic Materials and Devices, State Key Laboratory of Luminescent Materials and Devices, South China University of Technology, 381 Wushan Road, Guangzhou, Guangdong, 510640, P. R. China
| | - Qinglin Jiang
- Institute of Polymer Optoelectronic Materials and Devices, State Key Laboratory of Luminescent Materials and Devices, South China University of Technology, 381 Wushan Road, Guangzhou, Guangdong, 510640, P. R. China
| | - Jiang Zhang
- Department of Physics, South China University of Technology 381 Wushan Road, Guangzhou, Guangdong, 510640, P. R. China
| | - Yuguang Ma
- Institute of Polymer Optoelectronic Materials and Devices, State Key Laboratory of Luminescent Materials and Devices, South China University of Technology, 381 Wushan Road, Guangzhou, Guangdong, 510640, P. R. China
| |
Collapse
|
28
|
Futera Z, Wu X, Blumberger J. Tunneling-to-Hopping Transition in Multiheme Cytochrome Bioelectronic Junctions. J Phys Chem Lett 2023; 14:445-452. [PMID: 36622944 DOI: 10.1021/acs.jpclett.2c03361] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Multiheme cytochromes (MHCs) have attracted much interest for use in nanobioelectronic junctions due to their high electronic conductances. Recent measurements on dry MHC junctions suggested that a coherent tunneling mechanism is operative over surprisingly long long distances (>3 nm), which challenges our understanding of coherent transport phenomena. Here we show that this is due to (i) a low exponential distance decay constant for coherent conduction in MHCs (β = 0.2 Å-1) and (ii) a large density of protein electronic states which prolongs the coherent tunneling regime to distances that exceed those in molecular wires made of small molecules. Incoherent hopping conduction is uncompetitive due to the large energy level offset at the protein-electrode interface. Removing this offset, e.g., by gating, we predict that the transport mechanism crosses over from coherent tunneling to incoherent hopping at a protein size of ∼7 nm, thus enabling transport on the micrometer scale with a shallow polynomial (∼1/r) distance decay.
Collapse
Affiliation(s)
- Zdenek Futera
- Faculty of Science, University of South Bohemia, Branisovska 1760, 370 05 Ceske Budejovice, Czech Republic
| | - Xiaojing Wu
- University College London, Department of Physics and Astronomy, Gower Street, London WC1E 6BT, U.K
| | - Jochen Blumberger
- University College London, Department of Physics and Astronomy, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
29
|
Morab S, Sundaram MM, Pivrikas A. Time-Dependent Charge Carrier Transport with Hall Effect in Organic Semiconductors for Langevin and Non-Langevin Systems. NANOMATERIALS (BASEL, SWITZERLAND) 2022; 12:4414. [PMID: 36558267 PMCID: PMC9782042 DOI: 10.3390/nano12244414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Revised: 12/07/2022] [Accepted: 12/08/2022] [Indexed: 06/17/2023]
Abstract
The time-dependent charge carrier transport and recombination processes in low-mobility organic semiconductor diodes are obtained through numerical simulations using the finite element method (FEM). The application of a Lorentz force across the diode alters the charge transport process leading to the Hall effect. In this contribution, the Hall effect parameters, such as the Hall voltage and charge carrier concentration with varying magnetic fields, are computed for both Langevin and non-Langevin type recombination processes. The results indicate the charge carrier concentration within the diode for the Langevin system is about seven and fourteen times less while the maximum amount of extracted charge is nearly five and ten times less than that in the non-Langevin system of 0.01 and 0.001, respectively. The Hall voltage values obtained for the steady-state case are similar to the non-Langevin system of ββL=0.01. However, the values obtained for the Langevin and non-Langevin systems of ββL=1 and 0.001 exhibit anomalies. The implications of these findings advance the understanding of the charge transport and Hall effect measurements in organic semiconductors that underpins the device's performance.
Collapse
|
30
|
Shi R, Fang Q, Vasenko AS, Long R, Fang WH, Prezhdo OV. Structural Disorder in Higher-Temperature Phases Increases Charge Carrier Lifetimes in Metal Halide Perovskites. J Am Chem Soc 2022; 144:19137-19149. [PMID: 36206144 DOI: 10.1021/jacs.2c08627] [Citation(s) in RCA: 29] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Solar cells and optoelectronic devices are exposed to heat that degrades performance. Therefore, elucidating temperature-dependent charge carrier dynamics is essential for device optimization. Charge carrier lifetimes decrease with temperature in conventional semiconductors. The opposite, anomalous trend is observed in some experiments performed with MAPbI3 (MA = CH3NH3+) and other metal halide perovskites. Using ab initio quantum dynamics simulation, we establish the atomic mechanisms responsible for nonradiative electron-hole recombination in orthorhombic-, tetragonal-, and cubic MAPbI3. We demonstrate that structural disorder arising from the phase transitions is as important as the disorder due to heating in the same phase. The carrier lifetimes grow both with increasing temperature in the same phase and upon transition to the higher-temperature phases. The increased lifetime is rationalized by structural disorder that induces partial charge localization, decreases nonadiabatic coupling, and shortens quantum coherence. Inelastic and elastic electron-vibrational interactions exhibit opposite dependence on temperature and phase. The partial disorder and localization arise from thermal motions of both the inorganic lattice and the organic cations and depend significantly on the phase. The structural deformations induced by thermal fluctuations and phase transitions are on the same order as deformations induced by defects, and hence, thermal disorder plays a very important role. Since charge localization increases carrier lifetimes but inhibits transport, an optimal regime maximizing carrier diffusion can be designed, depending on phase, temperature, material morphology, and device architecture. The atomistic mechanisms responsible for the enhanced carrier lifetimes at elevated temperatures provide guidelines for the design of improved solar energy and optoelectronic materials.
Collapse
Affiliation(s)
- Ran Shi
- Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing100875, People's Republic of China
| | - Qiu Fang
- Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing100875, People's Republic of China
| | | | - Run Long
- Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing100875, People's Republic of China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing100875, People's Republic of China
| | - Oleg V Prezhdo
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, California90089, United States
| |
Collapse
|
31
|
Hu Z, Sun X. All-Atom Nonadiabatic Semiclassical Mapping Dynamics for Photoinduced Charge Transfer of Organic Photovoltaic Molecules in Explicit Solvents. J Chem Theory Comput 2022; 18:5819-5836. [PMID: 36073792 DOI: 10.1021/acs.jctc.2c00631] [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/30/2022]
Abstract
Direct all-atom simulation of nonadiabatic dynamics in disordered condensed phases like liquid solutions and amorphous solids has been challenging. The first all-atom simulation of the photoinduced charge-transfer dynamics of a prototypical organic photovoltaic carotenoid-porphyrin-C60 molecular triad in explicit tetrahydrofuran is presented. Based on the Meyer-Miller mapping Hamiltonian, various semiclassical and mixed quantum-classical dynamics are employed, including the linearized semiclassical, symmetrical quasiclassical, mean-field Ehrenfest, classical mapping model, and spin-mapping model approaches. The all-atom nonadiabatic dynamics were compared to multi-state harmonic models with a globally shared bath, and the models built using the ensemble averages on the initial electronic state could reproduce the all-atom results. The solvent effect was found to be critical for the photoinduced charge transfer, and the time-dependent solute-solvent radial distribution functions revealed that only the nonadiabatic dynamics started with the effective forces on the initial electronic state could capture the correct nuclear dynamics. The proposed strategy for modeling condensed-phase nonadiabatic dynamics with atomistic details is readily applied to complex condensed-phase systems.
Collapse
Affiliation(s)
- Zhubin Hu
- Division of Arts and Sciences, New York University Shanghai, 1555 Century Avenue, Shanghai 200122, China.,NYU-ECNU Center for Computational Chemistry, New York University Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China.,State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | - Xiang Sun
- Division of Arts and Sciences, New York University Shanghai, 1555 Century Avenue, Shanghai 200122, China.,NYU-ECNU Center for Computational Chemistry, New York University Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China.,State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
32
|
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
|
33
|
Sneyd A, Beljonne D, Rao A. A New Frontier in Exciton Transport: Transient Delocalization. J Phys Chem Lett 2022; 13:6820-6830. [PMID: 35857739 PMCID: PMC9340810 DOI: 10.1021/acs.jpclett.2c01133] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Accepted: 07/18/2022] [Indexed: 05/20/2023]
Abstract
Efficient exciton transport is crucial to the application of organic semiconductors (OSCs) in light-harvesting devices. While the physics of exciton transport in highly disordered media is well-explored, the description of transport in structurally and energetically ordered OSCs is less established, despite such materials being favorable for devices. In this Perspective we describe and highlight recent research pointing toward a highly efficient exciton transport mechanism which occurs in ordered OSCs, transient delocalization. Here, exciton-phonon couplings play a critical role in allowing localized exciton states to temporarily access higher-energy delocalized states whereupon they move large distances. The mechanism shows great promise for facilitating long-range exciton transport and may allow for improved device efficiencies and new device architectures. However, many fundamental questions on transient delocalization remain to be answered. These questions and suggested next steps are summarized.
Collapse
Affiliation(s)
- Alexander
J. Sneyd
- Department
of Physics, Cavendish Laboratory, University
of Cambridge, Cambridge CB3 0HE, United Kingdom
| | - David Beljonne
- Laboratory
for Chemistry of Novel Materials, University
of Mons, Mons 7000, Belgium
| | - Akshay Rao
- Department
of Physics, Cavendish Laboratory, University
of Cambridge, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
34
|
Pradhan CS, Jain A. Detailed Balance and Independent Electron Surface-Hopping Method: The Importance of Decoherence and Correct Calculation of Diabatic Populations. J Chem Theory Comput 2022; 18:4615-4626. [PMID: 35880817 DOI: 10.1021/acs.jctc.2c00320] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We benchmark and improve the independent electron surface-hopping (IESH) method developed by J. C. Tully's group for nonadiabatic simulations near metal surfaces. We have incorporated decoherence within the IESH method as well as implemented a scheme for the accurate calculation of diabatic populations. We benchmark the original IESH method with the above inclusions for a model system to calculate rate constants and long-time populations. The original IESH method fails to capture the detailed balance for some of the parameters, which is corrected with the inclusion of decoherence and accurate calculation of diabatic populations. Total rate constants are well captured both within the original IESH method as well as within our modified IESH.
Collapse
Affiliation(s)
- Chinmay S Pradhan
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
35
|
New Opportunities for Organic Semiconducting Polymers in Biomedical Applications. Polymers (Basel) 2022; 14:polym14142960. [PMID: 35890734 PMCID: PMC9318588 DOI: 10.3390/polym14142960] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 07/14/2022] [Accepted: 07/15/2022] [Indexed: 01/25/2023] Open
Abstract
The life expectancy of humans has been significantly elevated due to advancements in medical knowledge and skills over the past few decades. Although a lot of knowledge and skills are disseminated to the general public, electronic devices that quantitatively diagnose one’s own body condition still require specialized semiconductor devices which are huge and not portable. In this regard, semiconductor materials that are lightweight and have low power consumption and high performance should be developed with low cost for mass production. Organic semiconductors are one of the promising materials in biomedical applications due to their functionalities, solution-processability and excellent mechanical properties in terms of flexibility. In this review, we discuss organic semiconductor materials that are widely utilized in biomedical devices. Some advantageous and unique properties of organic semiconductors compared to inorganic semiconductors are reviewed. By critically assessing the fabrication process and device structures in organic-based biomedical devices, the potential merits and future aspects of the organic biomedical devices are pinpointed compared to inorganic devices.
Collapse
|
36
|
Wang B, Chu W, Wu Y, Casanova D, Saidi WA, Prezhdo OV. Electron-Volt Fluctuation of Defect Levels in Metal Halide Perovskites on a 100 ps Time Scale. J Phys Chem Lett 2022; 13:5946-5952. [PMID: 35732502 DOI: 10.1021/acs.jpclett.2c01452] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Metal halide perovskites (MHPs) have gained considerable attention due to their excellent optoelectronic performance, which is often attributed to unusual defect properties. We demonstrate that midgap defect levels can exhibit very large and slow energy fluctuations associated with anharmonic acoustic motions. Therefore, care should be taken classifying MHP defects as deep or shallow, since shallow defects may become deep and vice versa. As a consequence, charges from deep levels can escape into bands, and light absorption can be extended to longer wavelengths, improving material performance. The phenomenon, demonstrated with iodine vacancy in CH3NH3PbI3 using a machine learning force field, can be expected for a variety of defects and dopants in many MHPs and other soft inorganic semiconductors. Since large-scale anharmonic motions can be precursors to chemical decomposition, a known problem with MHPs, we propose that materials that are stiffer than MHPs but softer than traditional inorganic semiconductors, such as Si and TiO2, may simultaneously exhibit excellent performance and stability.
Collapse
Affiliation(s)
- Bipeng Wang
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, United States
| | - Weibin Chu
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Yifan Wu
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - David Casanova
- Donostia International Physics Center (DIPC), Donostia, 20018 Euskadi, Spain
- Basque Foundation for Science, IKERBASQUE, Bilbao, 48009 Euskadi, Spain
| | - Wissam A Saidi
- Department of Mechanical Engineering & Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Oleg V Prezhdo
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, United States
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
37
|
Shakiba M, Stippell E, Li W, Akimov AV. Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids. J Chem Theory Comput 2022; 18:5157-5180. [PMID: 35758936 DOI: 10.1021/acs.jctc.2c00297] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.
Collapse
Affiliation(s)
- Mohammad Shakiba
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Elizabeth Stippell
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, China
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
38
|
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
|
39
|
Ahart CS, Rosso KM, Blumberger J. Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package. J Chem Theory Comput 2022; 18:4438-4446. [PMID: 35700315 PMCID: PMC9281399 DOI: 10.1021/acs.jctc.2c00284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Constrained density
functional theory (CDFT) is a powerful tool
for the prediction of electron transfer parameters in condensed phase
simulations at a reasonable computational cost. In this work we present
an extension to CDFT in the popular mixed Gaussian/plane wave electronic
structure package CP2K, implementing the additional force terms arising
from a constraint based on Hirshfeld charge partitioning. This improves
upon the existing Becke partitioning scheme, which is prone to give
unphysical atomic charges. We verify this implementation for a variety
of systems: electron transfer in (H2O)2+ in a vacuum, electron tunnelling
between oxygen vacancy centers in solid MgO, and electron self-exchange
in aqueous Ru2+–Ru3+. We find good agreement
with previous plane-wave CDFT results for the same systems, but at
a significantly lower computational cost, and we discuss the general
reliability of condensed phase CDFT calculations.
Collapse
Affiliation(s)
- Christian S Ahart
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Kevin M Rosso
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
40
|
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
|