1
|
Wang Y, Li C, Zheng X. Markov State Models Reveal How Folding Kinetics Influence Absorption Spectra of Foldamers. J Chem Theory Comput 2024. [PMID: 38900275 DOI: 10.1021/acs.jctc.4c00202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
Self-assembly of platinum(II) complex foldamers is an essential approach to fabricate advanced luminescent materials. However, a comprehensive understanding of folding kinetics and their absorption spectra remains elusive. By constructing Markov state models (MSMs) from large-scale molecular dynamics simulations, we reveal that two largely similar dinuclear alknylplatinum(II) terpyridine foldamers, Pt-PEG and Pt-PE with slightly different bridges, exhibit distinctive folding kinetics. Particularly, Pt-PEG bears bridge-dominant, plane-dominant, and cooperative pathways, while Pt-PE only prefers the plane-dominant pathway. Such preference originates from their difference in intrabridge electrostatic interactions, leading to contrastive distributions of metastable states. We also found that the bridge-dominant pathway for Pt-PEG becomes more favorable when lowering the temperature. Interestingly, based on the comprehensive conformation ensembles from our MSMs, we reveal the conformation-dependent absorption spectra of Pt-PEG and Pt-PE. Our theoretical spectra not only align with experimental results but also reveal the contributions of diverse conformations to the overall absorption bands explicitly, facilitating the rational design of stimuli-responsive smart luminescent molecules.
Collapse
Affiliation(s)
- Yijia Wang
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, Key Laboratory of Medicinal Molecule Science and Pharmaceutics Engineering of Ministry of Industry and Information Technology, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China
| | - Chu Li
- Department of Chemistry, The Hong Kong University of Science and Technology, Hong Kong, China
| | - Xiaoyan Zheng
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, Key Laboratory of Medicinal Molecule Science and Pharmaceutics Engineering of Ministry of Industry and Information Technology, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China
| |
Collapse
|
2
|
Zhang L, Pios SV, Martyka M, Ge F, Hou YF, Chen Y, Chen L, Jankowska J, Barbatti M, Dral PO. MLatom Software Ecosystem for Surface Hopping Dynamics in Python with Quantum Mechanical and Machine Learning Methods. J Chem Theory Comput 2024. [PMID: 38836623 DOI: 10.1021/acs.jctc.4c00468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.
Collapse
Affiliation(s)
- Lina Zhang
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Sebastian V Pios
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Mikołaj Martyka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Fuchun Ge
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lipeng Chen
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| | - Pavlo O Dral
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| |
Collapse
|
3
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
4
|
Xu J, Shi Z, Wang L. Consistent Construction of the Density Matrix from Surface Hopping Trajectories. J Chem Theory Comput 2024; 20:2349-2361. [PMID: 38490993 DOI: 10.1021/acs.jctc.4c00137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2024]
Abstract
Proper construction of the density matrix based on surface hopping trajectories remains a difficult problem. Due to the well-known overcoherence in traditional surface hopping simulations, the electronic wave function cannot be used directly. In this work, we propose a consistent density matrix construction method, which takes the advantage of occupation of active states to rescale the coherence calculated by wave functions and ensures the intrinsic consistency of the density matrix. This new trajectory analysis method can be used for both Tully's fewest switches surface hopping (FSSH) and our recently proposed branching corrected surface hopping (BCSH). As benchmarked in both one- and two-dimensional standard scattering models, the new approach combined with BCSH trajectories achieves highly accurate time-dependent spatial distributions of adiabatic populations and coherence compared to exact quantum results.
Collapse
Affiliation(s)
- Jiabo Xu
- 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
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
5
|
Guo X, Li G, Shi Z, Wang L. Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching. J Phys Chem Lett 2024:3345-3353. [PMID: 38498301 DOI: 10.1021/acs.jpclett.4c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.
Collapse
Affiliation(s)
- Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Guijie Li
- 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
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
6
|
Cao Y, Halls MD, Friesner RA. Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method. J Chem Phys 2024; 160:084106. [PMID: 38385510 DOI: 10.1063/5.0188277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 01/31/2024] [Indexed: 02/23/2024] Open
Abstract
A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc., 1540 Broadway, 24th Floor, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger Inc., 9868 Scranton, Suite 3200, San Diego, California 92121, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
7
|
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
|
8
|
Liu XY, Chen WK, Fang WH, Cui G. Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications. J Chem Theory Comput 2023. [PMID: 37984502 DOI: 10.1021/acs.jctc.3c00960] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| |
Collapse
|
9
|
Li G, Shi Z, Guo X, Wang L. What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching. J Phys Chem Lett 2023; 14:9855-9863. [PMID: 37890155 DOI: 10.1021/acs.jpclett.3c02690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/29/2023]
Abstract
When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.
Collapse
Affiliation(s)
- Guijie Li
- 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
| | - Xin Guo
- 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
|
10
|
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
|
11
|
Zhao X, Merritt ICD, Lei R, Shu Y, Jacquemin D, Zhang L, Xu X, Vacher M, Truhlar DG. Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation. J Chem Theory Comput 2023; 19:6577-6588. [PMID: 37772732 DOI: 10.1021/acs.jctc.3c00813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/30/2023]
Abstract
Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.
Collapse
Affiliation(s)
- Xiaorui Zhao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, P. R. China
| | | | - Ruiqing Lei
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris 75005, France
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xuefei Xu
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, P. R. China
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
12
|
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
|
13
|
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
|
14
|
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: 19] [Impact Index Per Article: 9.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
|
15
|
Li H, Gong X, Ni H, Lu P, Luo X, Wen J, Yang Y, Qian X, Sun Z, Wu J. Light-Induced Ultrafast Molecular Dynamics: From Photochemistry to Optochemistry. J Phys Chem Lett 2022; 13:5881-5893. [PMID: 35730581 PMCID: PMC9251772 DOI: 10.1021/acs.jpclett.2c01119] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
By precisely controlling the waveform of ultrashort laser fields, electronic and nuclear motions in molecules can be steered on extremely short time scales, even in the attosecond regime. This new research field, termed "optochemistry", presents the light field in the time-frequency domain and opens new avenues for tailoring molecular reactions beyond photochemistry. This Perspective summarizes the ultrafast laser techniques employed in recent years for manipulating the molecular reactions based on waveform control of intense ultrashort laser pulses, where the chemical reactions can take place in isolated molecules, clusters, and various nanosystems. The underlying mechanisms for the coherent control of molecular dynamics are explicitly explored. Challenges and opportunities coexist in the field of optochemistry. Advanced technologies and theoretical modeling are still being pursued, with great prospects for controlling chemical reactions with unprecedented spatiotemporal precision.
Collapse
Affiliation(s)
- Hui Li
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Xiaochun Gong
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Hongcheng Ni
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Peifen Lu
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Xiao Luo
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200241, China
| | - Jin Wen
- State
Key Laboratory for Modification of Chemical Fibers and Polymer Materials,
College of Materials Science and Engineering, Donghua University, Shanghai 201620, China
| | - Youjun Yang
- State
Key Laboratory of Bioreactor Engineering, Shanghai Key Laboratory
of Chemical Biology, School of Pharmacy, East China University of Science and Technology, Shanghai 200237, China
| | - Xuhong Qian
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200241, China
| | - Zhenrong Sun
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| | - Jian Wu
- State
Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| |
Collapse
|
16
|
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
|
17
|
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
|
18
|
Buren B, Chen M. Stereodynamics-Controlled Product Branching in the Nonadiabatic H + NaD → Na(3s, 3p) + HD Reaction at Low Temperatures. J Phys Chem A 2022; 126:2453-2462. [PMID: 35434992 DOI: 10.1021/acs.jpca.2c00114] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nonadiabatic processes play an important role at energies near or higher than conical intersection of adiabatic potential energy surfaces in chemical reactions. In this work, dynamics of the nonadiabatic H + NaD reaction at low temperatures are studied by using the quantum wave packet method based on an improved L-shaped grid. The nonadiabatic H + NaD reaction has two exothermic reaction channels: Na(3s) + HD and Na(3p) + HD; the latter can only occur via nonadiabatic transition. The dynamics results show that the product branching of the H + NaD reaction at collision energies ranging from 20 to 80 cm-1 is controlled by stereodynamics. The Na(3s) and Na(3p) reaction channels occur through collinear collision and side-on collision, respectively. When the collision energy is lower than 20 cm-1, the resonance-mediated reaction mechanism is dominant in both the Na(3s) and Na(3p) reaction channels.
Collapse
Affiliation(s)
- Bayaer Buren
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| |
Collapse
|
19
|
Xie H, Xu X, Wang L, Zhuang W. Surface hopping dynamics in periodic solid-state materials with a linear vibronic coupling model. J Chem Phys 2022; 156:154116. [PMID: 35459287 DOI: 10.1063/5.0085759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a surface hopping approach in which the implemented linear vibronic coupling Hamiltonian is constructed and the electronic wavefunction is propagated in the reciprocal space. The parameters of the linear vibronic coupling model, including onsite energies, phonon frequencies, and electron-phonon couplings, are calculated with density-functional theory and density-functional perturbation theory and interpolated in fine sampling points of the Brillouin zone with maximally localized Wannier functions. Using this approach, we studied the relaxation dynamics of the photo-excited hot carrier in a one-dimensional periodic carbon chain. The results show that the completeness of the number of Hilbert space k points and the number of phonon q points plays an important role in the hot carrier relaxation processes. By calculating the relaxation times of hot carriers under different reciprocal space sampling and extrapolating with the stretched-compressed exponential function, the relaxation times of hot electrons and holes in the quasi-continuous energy band are obtained. By considering the feedback effect in the hopping processes and analyzing the time-dependent phonon energy in different normal modes, we found that the long-wave longitudinal optical phonons play a major role in the relaxation dynamics of hot electrons and holes. We, therefore, provided herein an efficient and accurate approach for modeling the photophysical processes in periodic solid-state material systems.
Collapse
Affiliation(s)
- Hua Xie
- Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, School of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaoliang Xu
- Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, School of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, Fujian 350002, China
| |
Collapse
|
20
|
Qiu J, Lu Y, Wang L. Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States. J Chem Theory Comput 2022; 18:2803-2815. [PMID: 35380833 DOI: 10.1021/acs.jctc.2c00130] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We present a multilayer subsystem surface hopping (MSSH) method to deal with nonadiabatic dynamics in large-scale systems. A small subsystem instead of the full system is adopted for surface hopping and is updated on-the-fly to achieve a reliable description of important adiabatic states and the wave function evolution. Additional subsystems for molecular dynamics and statistical description are introduced to further improve the simulation reliability. The global flux hopping probabilities with optimal state assignments are utilized to treat the complex surface crossings. As demonstrated in a series of one- and two-dimensional Holstein models with up to hundreds of thousands of states, MSSH shows weak parameter dependence in all investigated systems. Especially, the computational costs are reduced by 2-6 orders of magnitude compared to traditional surface hopping simulations in full systems, and size-independent results are achieved with a large time-step size of 2-5 fs. The new method is compatible with different decoherence correction strategies and achieves a much better balance between efficiency and reliability, thus promising for applications in general charge and exciton dynamics simulations.
Collapse
Affiliation(s)
- Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Yao Lu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
21
|
Ren J, Li W, Jiang T, Wang Y, Shuai Z. Time‐dependent density matrix renormalization group method for quantum dynamics in complex systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| |
Collapse
|
22
|
Liang R, Bakhtiiari A. Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug. J Phys Chem B 2022; 126:2382-2393. [PMID: 35297246 DOI: 10.1021/acs.jpcb.1c10819] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| | - Amirhossein Bakhtiiari
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| |
Collapse
|
23
|
Giannini S, Blumberger J. Charge Transport in Organic Semiconductors: The Perspective from Nonadiabatic Molecular Dynamics. Acc Chem Res 2022; 55:819-830. [PMID: 35196456 PMCID: PMC8928466 DOI: 10.1021/acs.accounts.1c00675] [Citation(s) in RCA: 32] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
Organic semiconductors (OSs) are an exciting
class of materials
that have enabled disruptive technologies in this century including
large-area electronics, flexible displays, and inexpensive solar cells.
All of these technologies rely on the motion of electrical charges
within the material and the diffusivity of these charges critically
determines their performance. In this respect, it is remarkable that
the nature of the charge transport in these materials has puzzled
the community for so many years, even for apparently simple systems
such as molecular single crystals: some experiments would better fit
an interpretation in terms of a localized particle picture, akin to
molecular or biological electron transfer, while others are in better
agreement with a wave-like interpretation, more akin to band transport
in metals. Exciting recent progress in the theory and simulation
of charge
carrier transport in OSs has now led to a unified understanding of
these disparate findings, and this Account will review one of these
tools developed in our laboratory in some detail: direct charge carrier
propagation by quantum-classical nonadiabatic molecular dynamics.
One finds that even in defect-free crystals the charge carrier can
either localize on a single molecule or substantially delocalize over
a large number of molecules depending on the relative strength of
electronic couplings between the molecules, reorganization, or charge
trapping energy of the molecule and thermal fluctuations of electronic
couplings and site energies, also known as electron–phonon
couplings. Our simulations predict that in molecular OSs exhibiting
some of
the highest measured charge mobilities to date, the charge carrier
forms “flickering” polarons, objects that are delocalized
over 10–20 molecules on average and that constantly change
their shape and extension under the influence of thermal disorder.
The flickering polarons propagate through the OS by short (≈10
fs long) bursts of the wave function that lead to an expansion of
the polaron to about twice its size, resulting in spatial displacement,
carrier diffusion, charge mobility, and electrical conductivity. Arguably
best termed “transient delocalization”, this mechanistic
scenario is very similar to the one assumed in transient localization
theory and supports its assertions. We also review recent applications
of our methodology to charge transport in disordered and nanocrystalline
samples, which allows us to understand the influence of defects and
grain boundaries on the charge propagation. Unfortunately, the
energetically favorable packing structures of
typical OSs, whether molecular or polymeric, places fundamental constraints
on charge mobilities/electronic conductivity compared to inorganic
semiconductors, which limits their range of applications. In this
Account, we review the design rules that could pave the way for new
very high-mobility OS materials and we argue that 2D covalent organic
frameworks are one of the most promising candidates to satisfy them. We conclude that our nonadiabatic dynamics method is a powerful
approach for predicting charge carrier transport in crystalline and
disordered materials. We close with a brief outlook on extensions
of the method to exciton transport, dissociation, and recombination.
This will bring us a step closer to an understanding of the birth,
survival, and annihiliation of charges at interfaces of optoelectronic
devices.
Collapse
Affiliation(s)
- Samuele Giannini
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
24
|
Li B, Xu J, Li G, Shi Z, Wang L. A Mixed Deterministic-Stochastic Algorithm of the Branching Corrected Mean Field Method for Nonadiabatic Dynamics. J Chem Phys 2022; 156:114116. [DOI: 10.1063/5.0084013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new algorithm of the branching corrected mean field (BCMF) method for nonadiabatic dynamics [J. Xu and L. Wang, J. Phys. Chem. Lett. 11, 8283 (2020)], which combines the key advantages of the two existed algorithms, i.e., the deterministic BCMF algorithm based on weights of trajectory branches (BCMF-w) and the stochastic BCMF algorithm with random collapse of the electronic wavefunction (BCMF-s). The resulting mixed deterministic-stochastic BCMF algorithm (BCMF-ws) is benchmarked in a series of standard scattering problems with potential wells on the excited-state surfaces, which are common in realistic systems. In all investigated cases, BCMF-ws holds the same high accuracy while the computational time is reduced about two orders of magnitude compared to the original BCMF-w and BCMF-s algorithms, thus promising for nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
| | | | | | | | - Linjun Wang
- Department of Chemistry, Zhejiang University, China
| |
Collapse
|
25
|
Dilmurat R, Prodhan S, Wang L, Beljonne D. Thermally activated intra-chain charge transport in high charge-carrier mobility copolymers. J Chem Phys 2022; 156:084115. [DOI: 10.1063/5.0082569] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Disordered or even seemingly amorphous, donor–acceptor type, conjugated copolymers with high charge-carrier mobility have emerged as a new class of functional materials, where transport along the conjugated backbone is key. Here, we report on non-adiabatic molecular dynamics simulations of charge-carrier transport along chains of poly (indacenodithiophene-co-benzothiadiazole), within a model Hamiltonian parameterized against first-principles calculations. We predict thermally activated charge transport associated with a slightly twisted ground-state conformation, on par with experimental results. Our results also demonstrate that the energy mismatch between the hole on the donor vs the acceptor units of the copolymer drives localization of the charge carriers and limits the intra-chain charge-carrier mobility. We predict that room-temperature mobility values in excess of 10 cm2 V−1 s−1 can be achieved through proper chemical tuning of the component monomer units.
Collapse
Affiliation(s)
- Rishat Dilmurat
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| | - Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Place du Parc, 20, 7000 Mons, Belgium
| |
Collapse
|
26
|
Wang Z, Dong J, Qiu J, Wang L. All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning. ACS APPLIED MATERIALS & INTERFACES 2022; 14:22929-22940. [PMID: 35100503 DOI: 10.1021/acsami.1c22181] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.
Collapse
Affiliation(s)
- Zedong Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiawei Dong
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
27
|
Kaiser S, Neumann T, Symalla F, Schlöder T, Fediai A, Friederich P, Wenzel W. De Novo Calculation of the Charge Carrier Mobility in Amorphous Small Molecule Organic Semiconductors. Front Chem 2022; 9:801589. [PMID: 35004618 PMCID: PMC8738089 DOI: 10.3389/fchem.2021.801589] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Accepted: 11/18/2021] [Indexed: 11/13/2022] Open
Abstract
Organic semiconductors (OSC) are key components in applications such as organic photovoltaics, organic sensors, transistors and organic light emitting diodes (OLED). OSC devices, especially OLEDs, often consist of multiple layers comprising one or more species of organic molecules. The unique properties of each molecular species and their interaction determine charge transport in OSCs-a key factor for device performance. The small charge carrier mobility of OSCs compared to inorganic semiconductors remains a major limitation of OSC device performance. Virtual design can support experimental R&D towards accelerated R&D of OSC compounds with improved charge transport. Here we benchmark a de novo multiscale workflow to compute the charge carrier mobility solely on the basis of the molecular structure: We generate virtual models of OSC thin films with atomistic resolution, compute the electronic structure of molecules in the thin films using a quantum embedding procedure and simulate charge transport with kinetic Monte-Carlo protocol. We show that for 15 common amorphous OSC the computed zero-field and field-dependent mobility are in good agreement with experimental data, proving this approach to be an effective virtual design tool for OSC materials and devices.
Collapse
Affiliation(s)
- Simon Kaiser
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | | | | | - Tobias Schlöder
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Artem Fediai
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Pascal Friederich
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,Institute of Theoretical Informatics, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
28
|
Bai X, Guo X, Wang L. Machine Learning Approach to Calculate Electronic Couplings between Quasi-diabatic Molecular Orbitals: The Case of DNA. J Phys Chem Lett 2021; 12:10457-10464. [PMID: 34672582 DOI: 10.1021/acs.jpclett.1c03053] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Diabatization of one-electron states in flexible molecular aggregates is a great challenge due to the presence of surface crossings between molecular orbital (MO) levels and the complex interaction between MOs of neighboring molecules. In this work, we present an efficient machine learning approach to calculate electronic couplings between quasi-diabatic MOs without the need of nonadiabatic coupling calculations. Using MOs of rigid molecules as references, the MOs that can be directly regarded to be quasi-diabatic in molecular dynamics are selected out, state tracked, and phase corrected. On the basis of this information, artificial neural networks are trained to characterize the structure-dependent onsite energies of quasi-diabatic MOs and the intermolecular electronic couplings. A representative sequence of DNA is systematically studied as an illustration. Smooth time evolution of electronic couplings in all base pairs is obtained with quasi-diabatic MOs. In particular, our method can calculate electronic couplings between different quasi-diabatic MOs independently, and thus, this possesses unique advantages in many applications.
Collapse
Affiliation(s)
- Xin Bai
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
29
|
Prodhan S, Giannini S, Wang L, Beljonne D. Long-Range Interactions Boost Singlet Exciton Diffusion in Nanofibers of π-Extended Polymer Chains. J Phys Chem Lett 2021; 12:8188-8193. [PMID: 34415752 DOI: 10.1021/acs.jpclett.1c02275] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Raising the distance covered by singlet excitons during their lifetimes to values maximizing light absorption (a few hundred nm) would solve the exciton diffusion bottleneck issue and lift the constraint for fine (∼10 nm) phase segregation in bulk heterojunction organic solar cells. In that context, the recent report of highly ordered conjugated polymer nanofibers featuring singlet exciton diffusion length, LD, in excess of 300 nm is both appealing and intriguing [Jin, X.; et al. Science 2018, 360 (6391), 897-900]. Here, on the basis of nonadiabatic molecular dynamics simulations, we demonstrate that singlet exciton diffusion in poly(3-hexylthiophene) (P3HT) fibers is highly sensitive to the interplay between delocalization along the polymer chains and long-range interactions along the stacks. Remarkably, the diffusion coefficient is predicted to rocket by 3 orders of magnitude when going beyond nearest-neighbor intermolecular interactions in fibers of extended (30-mer) polymer chains and to be resilient to interchain energetic and positional disorders.
Collapse
Affiliation(s)
- Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| | - Samuele Giannini
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| |
Collapse
|
30
|
Sneyd AJ, Fukui T, Paleček D, Prodhan S, Wagner I, Zhang Y, Sung J, Collins SM, Slater TJA, Andaji-Garmaroudi Z, MacFarlane LR, Garcia-Hernandez JD, Wang L, Whittell GR, Hodgkiss JM, Chen K, Beljonne D, Manners I, Friend RH, Rao A. Efficient energy transport in an organic semiconductor mediated by transient exciton delocalization. SCIENCE ADVANCES 2021; 7:7/32/eabh4232. [PMID: 34348902 PMCID: PMC8336960 DOI: 10.1126/sciadv.abh4232] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Accepted: 06/14/2021] [Indexed: 05/12/2023]
Abstract
Efficient energy transport is desirable in organic semiconductor (OSC) devices. However, photogenerated excitons in OSC films mostly occupy highly localized states, limiting exciton diffusion coefficients to below ~10-2 cm2/s and diffusion lengths below ~50 nm. We use ultrafast optical microscopy and nonadiabatic molecular dynamics simulations to study well-ordered poly(3-hexylthiophene) nanofiber films prepared using living crystallization-driven self-assembly, and reveal a highly efficient energy transport regime: transient exciton delocalization, where energy exchange with vibrational modes allows excitons to temporarily re-access spatially extended states under equilibrium conditions. We show that this enables exciton diffusion constants up to 1.1 ± 0.1 cm2/s and diffusion lengths of 300 ± 50 nm. Our results reveal the dynamic interplay between localized and delocalized exciton configurations at equilibrium conditions, calling for a re-evaluation of exciton dynamics and suggesting design rules to engineer efficient energy transport in OSC device architectures not based on restrictive bulk heterojunctions.
Collapse
Affiliation(s)
- Alexander J Sneyd
- Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK
| | - Tomoya Fukui
- Department of Chemistry, University of Victoria, Victoria, BC V8P 5C2, Canada
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - David Paleček
- Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK
| | - Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| | - Isabella Wagner
- MacDiarmid Institute for Advanced Materials and Nanotechnology and School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6010, New Zealand
| | - Yifan Zhang
- Department of Chemistry, University of Victoria, Victoria, BC V8P 5C2, Canada
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Jooyoung Sung
- Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK
| | - Sean M Collins
- School of Chemical and Process Engineering and School of Chemistry, University of Leeds, Leeds LS2 9JT, UK
| | - Thomas J A Slater
- Electron Physical Science Imaging Centre, Diamond Light Source Ltd., Oxfordshire OX11 0DE, UK
| | | | - Liam R MacFarlane
- Department of Chemistry, University of Victoria, Victoria, BC V8P 5C2, Canada
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - J Diego Garcia-Hernandez
- Department of Chemistry, University of Victoria, Victoria, BC V8P 5C2, Canada
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | | | - Justin M Hodgkiss
- MacDiarmid Institute for Advanced Materials and Nanotechnology and School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6010, New Zealand
| | - Kai Chen
- MacDiarmid Institute for Advanced Materials and Nanotechnology and School of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6010, New Zealand
- Robinson Research Institute, Faculty of Engineering, Victoria University of Wellington, Wellington 6012, New Zealand
- The Dodd-Walls Centre for Photonic and Quantum Technologies, Dunedin 9016, New Zealand
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium.
| | - Ian Manners
- Department of Chemistry, University of Victoria, Victoria, BC V8P 5C2, Canada.
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Richard H Friend
- Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK
| | - Akshay Rao
- Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK.
| |
Collapse
|
31
|
Shao C, Xu J, Wang L. Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems. J Chem Phys 2021; 154:234109. [PMID: 34241240 DOI: 10.1063/5.0056224] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
32
|
Wang Y, Ren J, Shuai Z. Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG. J Chem Phys 2021; 154:214109. [PMID: 34240969 DOI: 10.1063/5.0052804] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%-40% higher than the rates under HA.
Collapse
Affiliation(s)
- Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
33
|
Mandal A, Hunt KLC. Quantum transition probabilities due to overlapping electromagnetic pulses: Persistent differences between Dirac's form and nonadiabatic perturbation theory. J Chem Phys 2021; 154:024116. [PMID: 33445917 DOI: 10.1063/5.0020169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The probability of transition to an excited state of a quantum system in a time-dependent electromagnetic field determines the energy uptake from the field. The standard expression for the transition probability has been given by Dirac. Landau and Lifshitz suggested, instead, that the adiabatic effects of a perturbation should be excluded from the transition probability, leaving an expression in terms of the nonadiabatic response. In our previous work, we have found that these two approaches yield different results while a perturbing field is acting on the system. Here, we prove, for the first time, that differences between the two approaches may persist after the perturbing fields have been completely turned off. We have designed a pair of overlapping pulses in order to establish the possibility of lasting differences, in a case with dephasing. Our work goes beyond the analysis presented by Landau and Lifshitz, since they considered only linear response and required that a constant perturbation must remain as t → ∞. First, a "plateau" pulse populates an excited rotational state and produces coherences between the ground and excited states. Then, an infrared pulse acts while the electric field of the first pulse is constant, but after dephasing has occurred. The nonadiabatic perturbation theory permits dephasing, but dephasing of the perturbed part of the wave function cannot occur within Dirac's method. When the frequencies in both pulses are on resonance, the lasting differences in the calculated transition probabilities may exceed 35%. The predicted differences are larger for off-resonant perturbations.
Collapse
Affiliation(s)
- Anirban Mandal
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Katharine L C Hunt
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
34
|
Chu W, Zheng Q, Akimov AV, Zhao J, Saidi WA, Prezhdo OV. Accurate Computation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials. J Phys Chem Lett 2020; 11:10073-10080. [PMID: 33179939 DOI: 10.1021/acs.jpclett.0c03080] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Synergy of nonadiabatic molecular dynamics with real-time time-dependent density functional theory has led to significant progress in modeling excited-state dynamics in nanoscale and condensed matter systems over the past decade. Nonadiabatic coupling (NAC) is the central quantity in such simulations, and its accurate and efficient evaluation is an enduring challenge in time-dependent Kohn-Sham theory, particularly in conjunction with planewave basis sets and projector augmented-wave (PAW) pseudopotentials because of the complexity of the PAW "all-electron" wave function. We report a method for rigorous evaluation of the NAC with PAW wave functions and demonstrate an efficient approximation to the rigorous NAC that gives comparable accuracy. As a validation, we intensely examine the NAC matrix elements calculated using both pseudo- and all-electron wave functions under the PAW formalism in six representative systems. The approximate NAC obtained with pseudowave functions is close to the exact all-electron NAC, with the largest deviations observed when subshell d-electrons are involved in the transitions. The developed approach provides a rigorous and convenient methodology for the numerical computation of NAC in the Kohn-Sham theory framework.
Collapse
Affiliation(s)
- Weibin Chu
- Department of Chemistry and Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| | - Qijing Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- ICQD/Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Strongly Coupled Quantum Matter Physics, and Synergetic Innovation Center of Quantum Information & Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh Pennsylvania 15260, United States
| | - Wissam A Saidi
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Oleg V Prezhdo
- Department of Chemistry and Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
35
|
Abstract
When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
- Jiabo Xu
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials, and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
36
|
Shuai Z, Li W, Ren J, Jiang Y, Geng H. Applying Marcus theory to describe the carrier transports in organic semiconductors: Limitations and beyond. J Chem Phys 2020; 153:080902. [PMID: 32872875 DOI: 10.1063/5.0018312] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.
Collapse
Affiliation(s)
- Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Yuqian Jiang
- Laboratory for Nanosystem and Hierarchy Fabrication, National Center for Nanoscience and Technology, Chinese Academy of Sciences, 100084 Beijing, People's Republic of China
| | - Hua Geng
- Department of Chemistry, Capital Normal University, 100048 Beijing, People's Republic of China
| |
Collapse
|
37
|
Prodhan S, Qiu J, Ricci M, Roscioni OM, Wang L, Beljonne D. Design Rules to Maximize Charge-Carrier Mobility along Conjugated Polymer Chains. J Phys Chem Lett 2020; 11:6519-6525. [PMID: 32692920 DOI: 10.1021/acs.jpclett.0c01793] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The emergence of polymeric materials displaying high charge-carrier mobility values despite poor interchain structural order has spawned a renewal of interest in the identification of structure-property relationships pertaining to the transport of charges along conjugated polymer chains and the subsequent design of optimized architectures. Here, we present the results of intrachain charge transport simulations obtained by applying a robust surface hopping algorithm to a phenomenological Hamiltonian parametrized against first-principles simulations. Conformational effects are shown to provide a clear signature in the temperature-dependent charge-carrier mobility that complies with recent experimental observations. We further contrast against molecular crystals the evolution with electronic bandwidth and electron-phonon interactions of the room-temperature mobility in polymers, showing that intrachain charge-carrier mobility values in excess of 100 cm2/(V s) can be achieved through a proper chemical engineering of the backbones.
Collapse
Affiliation(s)
- Suryoday Prodhan
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| | - Jing Qiu
- Center for Chemistry of Novel & High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | | | | | - Linjun Wang
- Center for Chemistry of Novel & High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons 7000, Belgium
| |
Collapse
|
38
|
Uratani H, Nakai H. Non-adiabatic molecular dynamics with divide-and-conquer type large-scale excited-state calculations. J Chem Phys 2020; 152:224109. [DOI: 10.1063/5.0006831] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
39
|
Nematiaram T, Troisi A. Modeling charge transport in high-mobility molecular semiconductors: Balancing electronic structure and quantum dynamics methods with the help of experiments. J Chem Phys 2020; 152:190902. [DOI: 10.1063/5.0008357] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Tahereh Nematiaram
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Alessandro Troisi
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| |
Collapse
|
40
|
Chu W, Zheng Q, Prezhdo OV, Zhao J. CO2 Photoreduction on Metal Oxide Surface Is Driven by Transient Capture of Hot Electrons: Ab Initio Quantum Dynamics Simulation. J Am Chem Soc 2020; 142:3214-3221. [DOI: 10.1021/jacs.9b13280] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Weibin Chu
- ICQD/Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, and Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemistry and Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| | - Qijing Zheng
- ICQD/Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, and Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Oleg V. Prezhdo
- Department of Chemistry and Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| | - Jin Zhao
- ICQD/Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, and Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
- Synergetic Innovation Center of Quantum Information & Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
41
|
Zhong YJ, Lan CF, Lin BC, Hu CD, Cheng YC, Hsu CP. The anisotropy and temperature dependence in the mobility of rubrene. ADVANCES IN QUANTUM CHEMISTRY 2020. [DOI: 10.1016/bs.aiq.2020.04.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|