1
|
Liu D, Wang B, Vasenko AS, Prezhdo OV. Decoherence ensures convergence of non-adiabatic molecular dynamics with number of states. J Chem Phys 2024; 161:064104. [PMID: 39120030 DOI: 10.1063/5.0222557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Accepted: 07/24/2024] [Indexed: 08/10/2024] Open
Abstract
Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems.
Collapse
Affiliation(s)
| | - Bipeng Wang
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Andrey S Vasenko
- HSE University, 101000 Moscow, Russia
- Donostia International Physics Center (DIPC), 20018 San Sebastián-Donostia, Euskadi, Spain
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
2
|
Schirato A, Sanders SK, Proietti Zaccaria R, Nordlander P, Della Valle G, Alabastri A. Quantifying Ultrafast Energy Transfer from Plasmonic Hot Carriers for Pulsed Photocatalysis on Nanostructures. ACS NANO 2024; 18:18933-18947. [PMID: 38990155 DOI: 10.1021/acsnano.4c01802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
Photocatalysis with plasmonic nanostructures has lately emerged as a transformative paradigm to drive and alter chemical reactions using light. At the surface of metallic nanoparticles, photoexcitation results in strong near fields, short-lived high-energy "hot" carriers, and light-induced heating, thus creating a local environment where reactions can occur with enhanced efficiencies. In this context, it is critical to understand how to manipulate the nonequilibrium processes triggered by light, as their ultrafast (femto- to picoseconds) relaxation dynamics compete with the process of energy transfer toward the reactants. Accurate predictions of the plasmon photocatalytic activity can lead to optimized nanophotonic architectures with enhanced selectivity and rates, operating beyond the intrinsic limitations of the steady state. Here, we report on an original modeling approach to quantify, with space, time, and energy resolution, the ultrafast energy exchange from plasmonic hot carriers (HCs) to molecular systems adsorbed on the metal nanoparticle surface while consistently accounting for photothermal bond activation. Our analysis, illustrated for a few typical cases, reveals that the most energetic nonequilibrium carriers (i.e., with energies well far from the Fermi level) may introduce a wavelength-dependence of the reaction rates, and it elucidates on the role of the carriers closer to the Fermi energy and the photothermally heated lattice, suggesting ways to enhance and optimize each contribution. We show that the overall reaction rates can benefit strongly from using pulsed illumination with the optimal pulse width determined by the properties of the system. Taken together, these results contribute to the rational design of nanoreactors for pulsed catalysis, which calls for predictive modeling of the ultrafast HC-hot adsorbate energy transfer.
Collapse
Affiliation(s)
- Andrea Schirato
- Department of Physics, Politecnico di Milano, Milano 20133, Italy
- Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, United States
- Istituto Italiano di Tecnologia, Genoa 16163, Italy
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| | - Stephen Keith Sanders
- Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, United States
| | | | - Peter Nordlander
- Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, United States
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
- Laboratory for Nanophotonics, Rice University, Houston, Texas 77005, United States
| | - Giuseppe Della Valle
- Department of Physics, Politecnico di Milano, Milano 20133, Italy
- Istituto di Fotonica e Nanotecnologie─Consiglio Nazionale delle Ricerche, Piazza Leonardo da Vinci 32, Milano I-20133, Italy
| | - Alessandro Alabastri
- Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, United States
- Laboratory for Nanophotonics, Rice University, Houston, Texas 77005, United States
- Smalley-Curl Institute, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
3
|
Huang L, Shi Z, Wang L. Detailed Complementary Consistency: Wave Function Tells Particle How to Hop, Particle Tells Wave Function How to Collapse. J Phys Chem Lett 2024; 15:6771-6781. [PMID: 38912973 DOI: 10.1021/acs.jpclett.4c01313] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
In mixed quantum-classical dynamics, the quantum subsystem can have both wave function and particle-like descriptions. However, they may yield inconsistent results for the expectation value of the same physical quantity. We here propose a novel detailed complementary consistency (DCC) method based on the principle of detailed internal consistency. Namely, the wave function along each trajectory tells the particle how to hop, while the particle tells the wave function how to collapse based on active states in the trajectory ensemble. As benchmarked in a diverse array of representative models with localized nonadiabatic couplings, DCC not only achieves fully consistent results (i.e., identical populations calculated based on wave functions and active states) but also closely reproduces the exact quantum results. Due to the high performance, our new DCC method has great potential to give a consistent and accurate mixed quantum-classical description of general nonadiabatic dynamics after further development.
Collapse
Affiliation(s)
- Lei Huang
- 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
|
4
|
Kim CW, Franco I. General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation. J Chem Phys 2024; 160:214111. [PMID: 38833366 DOI: 10.1063/5.0202860] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2024] [Accepted: 05/13/2024] [Indexed: 06/06/2024] Open
Abstract
We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima-Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Chonnam National University, Gwangju 61186, South Korea
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
- Department of Physics, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
5
|
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
|
6
|
Suchan J, Liang F, Durden AS, Levine BG. Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone. J Chem Phys 2024; 160:134310. [PMID: 38573851 DOI: 10.1063/5.0198333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/18/2024] [Indexed: 04/06/2024] Open
Abstract
Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?
Collapse
Affiliation(s)
- Jiří Suchan
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Yang KR, Kyro GW, Batista VS. The landscape of computational approaches for artificial photosynthesis. NATURE COMPUTATIONAL SCIENCE 2023; 3:504-513. [PMID: 38177419 DOI: 10.1038/s43588-023-00450-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 04/11/2023] [Indexed: 01/06/2024]
Abstract
Artificial photosynthesis is an attractive strategy for converting solar energy into fuels, largely because the Earth receives enough solar energy in one hour to meet humanity's energy needs for an entire year. However, developing devices for artificial photosynthesis remains difficult and requires computational approaches to guide and assist the interpretation of experiments. In this Perspective, we discuss current and future computational approaches, as well as the challenges of designing and characterizing molecular assemblies that absorb solar light, transfer electrons between interfaces, and catalyze water-splitting and fuel-forming reactions.
Collapse
Affiliation(s)
- Ke R Yang
- Department of Chemistry, Yale University, New Haven, CT, USA
- Energy Sciences Institute, Yale University, West Haven, CT, USA
| | - Gregory W Kyro
- Department of Chemistry, Yale University, New Haven, CT, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, CT, USA.
- Energy Sciences Institute, Yale University, West Haven, CT, USA.
| |
Collapse
|
10
|
Li W, Xue T, Mora-Perez C, Prezhdo OV. Ab initio quantum dynamics of plasmonic charge carriers. TRENDS IN CHEMISTRY 2023. [DOI: 10.1016/j.trechm.2023.02.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
|
11
|
Ding Y, Xue K, Zhang J, Yan L, Li Q, Yao Y, Zhou L. Two-Dimensional Octuple-Atomic-Layer M 2Si 2N 4 (M = Al, Ga and In) with Long Carrier Lifetime. MICROMACHINES 2023; 14:405. [PMID: 36838105 PMCID: PMC9966885 DOI: 10.3390/mi14020405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 02/03/2023] [Accepted: 02/03/2023] [Indexed: 06/18/2023]
Abstract
Bulk III-nitride materials MN (M = Al, Ga and In) and their alloys have been widely used in high-power electronic and optoelectronic devices, but stable two-dimensional (2D) III-nitride materials, except h-BN, have not been realized yet. A new kind of 2D III-nitride material M2Si2N4 (M = Al, Ga and In) is predicted by choosing Si as the appropriate passivation element. The stability, electronic and optical properties of 2D M2Si2N4 materials are studied systematically based on first-principles calculations. The results show that Al2Si2N4 and Ga2Si2N4 are found to be indirect bandgap semiconductors, while In2Si2N4 is a direct bandgap semiconductor. Moreover, Al2Si2N4 and In2Si2N4 have good absorption ability in the visible light region, while Ga2Si2N4 is an ultraviolet-light-absorbing material. Furthermore, the carrier lifetimes of Ga2Si2N4 and In2Si2N4 are as large as 157.89 and 103.99 ns, respectively. All these desirable properties of M2Si2N4 materials make them attractive for applications in electronics and photoelectronics.
Collapse
Affiliation(s)
- Yimin Ding
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
| | - Kui Xue
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Jing Zhang
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Luo Yan
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Qiaoqiao Li
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Yisen Yao
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
- Institute of Fundamental and Frontier Science, University of Electronic Science and Technology of China, Chengdu 611731, China
| | - Liujiang Zhou
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| |
Collapse
|
12
|
Xie RF, Zhang JB, Wu Y, Li L, Liu XY, Cui G. Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube. J Chem Phys 2023; 158:054108. [PMID: 36754819 DOI: 10.1063/5.0134353] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E11 excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E11 excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out-in+〉 and |out+in-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E22 state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.
Collapse
Affiliation(s)
- Rui-Fang Xie
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Jing-Bin Zhang
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Yang Wu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Laicai Li
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, 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
|
Mondelo-Martell M, Brey D, Burghardt I. Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models. J Chem Phys 2022; 157:094108. [DOI: 10.1063/5.0104729] [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 on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.
Collapse
Affiliation(s)
- Manel Mondelo-Martell
- Institut für Physikalische u. Theoretische Chemie, Goethe-Universitat Frankfurt am Main Institut fur Physikalische und Theoretische Chemie, Germany
| | | | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Germany
| |
Collapse
|
15
|
Mukherjee S, Pinheiro M, Demoulin B, Barbatti M. Simulations of molecular photodynamics in long timescales. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200382. [PMID: 35341303 PMCID: PMC8958277 DOI: 10.1098/rsta.2020.0382] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 08/12/2021] [Indexed: 05/04/2023]
Abstract
Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
| | - Max Pinheiro
- Aix Marseille University, CNRS, ICR, Marseille, France
| | | | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| |
Collapse
|
16
|
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
|
17
|
Modelling Ultrafast Dynamics at a Conical Intersection with Regularized Diabatic States: An Approach Based on Multiplicative Neural Networks. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
18
|
Xie BB, Jia PK, Wang KX, Chen WK, Liu XY, Cui G. Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems. J Phys Chem A 2022; 126:1789-1804. [PMID: 35266391 DOI: 10.1021/acs.jpca.1c10195] [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/30/2022]
Abstract
Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Pei-Ke Jia
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Ke-Xin Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, 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
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, Sichuan, 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
| |
Collapse
|
19
|
Zhang L, Chu W, Zheng Q, Zhao J. Effects of oxygen vacancies on the photoexcited carrier lifetime in rutile TiO 2. Phys Chem Chem Phys 2022; 24:4743-4750. [PMID: 35142307 DOI: 10.1039/d1cp04248c] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The photoexcited carrier lifetime in semiconductors plays a crucial role in solar energy conversion processes. The defects or impurities in semiconductors are usually proposed to introduce electron-hole (e-h) recombination centers and consequently reduce the photoexcited carrier lifetime. In this report, we investigate the effects of oxygen vacancies (OV) on the carrier lifetime in rutile TiO2, which has important applications in photocatalysis and photovoltaics. It is found that an OV introduces two excess electrons which form two defect states in the band gap. The lower state is localized on one Ti atom and behaves as a small polaron, and the higher one is a hybrid state contributed by three Ti atoms around the OV. Both the polaron and hybrid states exhibit strong electron-phonon (e-ph) coupling and their charge distributions become more and more delocalized when the temperature increases from 100 to 700 K. Such strong e-ph coupling and charge delocalization enhance the nonadibatic coupling between the electronic states along the hole relaxation path, where the defect states behave as intermediate states, leading to a distinct acceleration of e-h recombination. Our study provides valuable insights to understand the role of defects on photoexcited carrier lifetime in semiconductors.
Collapse
Affiliation(s)
- Lili Zhang
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China. .,Key Laboratory of Material Physics, Ministry of Education, School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450001, China
| | - Weibin Chu
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China. .,Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA
| | - Qijing Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
20
|
Gumber S, Agrawal S, Prezhdo OV. Excited State Dynamics in Dual-Defects Modified Graphitic Carbon Nitride. J Phys Chem Lett 2022; 13:1033-1041. [PMID: 35073096 DOI: 10.1021/acs.jpclett.1c04152] [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/14/2023]
Abstract
Significant efforts are focused on defect-engineering of metal-free graphitic carbon nitride (g-C3N4) to amplify its efficacy. A conceptually new multidefect-modified g-C3N4 having simultaneously two or more defects has attracted strong attention for its enhanced photocatalytic properties. We model and compare the excited state dynamics in g-C3N4 with (i) nitrogen defects (N vacancy and CN group) and (ii) dual defects (N vacancy, CN group, and O doping) and show that the nonradiative recombination of charge carriers in these systems follows the Shockley-Read-Hall mechanism. The nitrogen defects create three midgap states that trap charges and act as recombination centers. The dual-defect modified systems exhibit superior properties compared with pristine g-C3N4 because the defects facilitate rapid charge separation and extend the spectrum of absorbed light. The system doped with O shows better performance due to enhanced carrier lifetime and higher oxidation potential caused by a downshifted valence band. The study provides guidance for rational design of stable and efficient photocatalytic materials.
Collapse
Affiliation(s)
- Shriya Gumber
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Sraddha Agrawal
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
21
|
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
|
22
|
Akimov AV. Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain. J Phys Chem Lett 2021; 12:12119-12128. [PMID: 34913701 DOI: 10.1021/acs.jpclett.1c03823] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
23
|
Albareda G, Lively K, Sato SA, Kelly A, Rubio A. Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics. J Chem Theory Comput 2021; 17:7321-7340. [PMID: 34752108 PMCID: PMC8675140 DOI: 10.1021/acs.jctc.1c00772] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Indexed: 11/28/2022]
Abstract
We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.
Collapse
Affiliation(s)
- Guillermo Albareda
- Nano-Bio
Spectroscopy Group and European Theoretical Spectroscopy Facility
(ETSF), Universidad del País Vasco
(UPV/EHU), Av. Tolosa
72, 20018 San Sebastian, Spain
- Institute
of Theoretical and Computational Chemistry, University of Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Kevin Lively
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Shunsuke A. Sato
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center
for Computational Sciences, University of
Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Aaron Kelly
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Department
of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
| | - Angel Rubio
- Nano-Bio
Spectroscopy Group and European Theoretical Spectroscopy Facility
(ETSF), Universidad del País Vasco
(UPV/EHU), Av. Tolosa
72, 20018 San Sebastian, Spain
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Luruper Chaussee 149, 22761 Hamburg, Germany
- The
Hamburg Centre for Ultrafast Imaging, University
of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center
for Computational Quantum Physics (CCQ), Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United
States
| |
Collapse
|
24
|
Esch MP, Levine BG. An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics. J Chem Phys 2021; 155:214101. [PMID: 34879667 DOI: 10.1063/5.0070686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.
Collapse
Affiliation(s)
- Michael P Esch
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
25
|
|
26
|
Liu J, Lan Z, Yang J. An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics. Phys Chem Chem Phys 2021; 23:22313-22323. [PMID: 34591049 DOI: 10.1039/d1cp03477d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
27
|
Mangan SM, Zhou G, Chu W, Prezhdo OV. Dependence between Structural and Electronic Properties of CsPbI 3: Unsupervised Machine Learning of Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2021; 12:8672-8678. [PMID: 34472856 DOI: 10.1021/acs.jpclett.1c02361] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Using unsupervised machine learning on the trajectories from a nonadiabatic molecular dynamics simulation with time-dependent Kohn-Sham density functional theory, we elucidated the structural parameters with the largest influence on nonradiative recombination of charge carriers in CsPbI3, which forms the basis for solar energy and optoelectronic applications. The I-I-I angles between PbI6 octahedra, followed by the Cs-I distance, have the strongest impact on the bandgap and the nonadiabatic coupling. The importance of the Cs-I distance is unexpected, because Cs does not contribute to electron and hole wave functions. The nonadiabatic coupling is most influenced by static properties, which is also surprising, given its explicit dependence on atomic velocities.
Collapse
Affiliation(s)
- Spencer M Mangan
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, United States
| | - Guoqing Zhou
- Department of Physics and Astronomy, 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
| | - Oleg V Prezhdo
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
28
|
Xie XY, Yang JJ, Liu XY, Fang Q, Fang WH, Cui G. Interfacial photoinduced carrier dynamics tuned by polymerization of coronene molecules encapsulated in carbon nanotubes: bridging type-I and type-II heterojunctions. Phys Chem Chem Phys 2021; 23:13503-13511. [PMID: 34120157 DOI: 10.1039/d1cp01008e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Carbon nanomaterials play important roles in modern scientific research. Integrating different carbon-based building blocks into nano-hybrid architectures not only takes full advantage of each component, but also brings in novel interfacial properties. Herein, we have employed density functional theory (DFT) calculations to investigate the effects of polymerization degree of coronene molecules encapsulated in single-walled carbon nanotubes (SWNTs) (19,0) on their interfacial properties. The present results reveal that the interfacial properties of the formed heterojunctions are remarkably regulated by the polymerization degree. For example, monomer- and dimer-encapsulated SWNTs are type-I heterojunctions in which interfacial excitation energy transfer is preferred, whereas interfacial charge carrier transfer is favorable in trimer- and polymer-encapsulated SWNTs because they are type-II heterojunctions. On the other hand, we have employed the time-domain nonadiabatic dynamics simulation approach to explore the interfacial carrier dynamics in type-II polymer-encapsulated SWNT heterojunctions. It is found that the electron and hole transfer processes are asymmetric and occur in opposite directions and at different rates. The former takes place from polymers to SWNTs in an ultrafast way (ca. 370 fs), whereas the latter occurs slowly from SWNTs to polymers (ca. 24 ps). A closer analysis uncovers the fact that the different carrier transfer rates mainly originate from the different densities of the acceptor states, energy differences and inter-state couplings between the donor and acceptor states. Finally, the present work demonstrates that the polymerization degree could act as a new regulating strategy to tune the interfacial properties of molecule-encapsulated SWNT heterojunctions.
Collapse
Affiliation(s)
- Xiao-Ying Xie
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai, 264005, China
| | - Jia-Jia Yang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Qiu Fang
- 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
|
29
|
Li W, She Y, Vasenko AS, Prezhdo OV. Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites. NANOSCALE 2021; 13:10239-10265. [PMID: 34031683 DOI: 10.1039/d1nr01990b] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.
Collapse
Affiliation(s)
- Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China.
| | | | | | | |
Collapse
|
30
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
31
|
Popp W, Brey D, Binder R, Burghardt I. Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems. Annu Rev Phys Chem 2021; 72:591-616. [PMID: 33636997 DOI: 10.1146/annurev-physchem-090419-040306] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.
Collapse
Affiliation(s)
- Wjatscheslaw Popp
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| |
Collapse
|
32
|
Yamijala SSRKC, Huo P. Direct Nonadiabatic Simulations of the Photoinduced Charge Transfer Dynamics. J Phys Chem A 2021; 125:628-635. [DOI: 10.1021/acs.jpca.0c10151] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Sharma S. R. K. C. Yamijala
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department of Chemistry, Indian Institute of Technology-Madras, Chennai 600036, India
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
33
|
Di Maiolo F, Brey D, Binder R, Burghardt I. Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature. J Chem Phys 2020; 153:184107. [PMID: 33187420 DOI: 10.1063/5.0027588] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.
Collapse
Affiliation(s)
- Francesco Di Maiolo
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
34
|
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
|
35
|
Esch MP, Levine BG. Decoherence-corrected Ehrenfest molecular dynamics on many electronic states. J Chem Phys 2020; 153:114104. [DOI: 10.1063/5.0022529] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
36
|
Hegger R, Binder R, Burghardt I. First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature. J Chem Theory Comput 2020; 16:5441-5455. [PMID: 32786907 DOI: 10.1021/acs.jctc.0c00351] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σs ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10-2 cm2 s-1 to D ∼ 1 × 10-3 cm2 s-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.
Collapse
Affiliation(s)
- Rainer Hegger
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
37
|
Chu W, Saidi WA, Prezhdo OV. Long-Lived Hot Electron in a Metallic Particle for Plasmonics and Catalysis: Ab Initio Nonadiabatic Molecular Dynamics with Machine Learning. ACS NANO 2020; 14:10608-10615. [PMID: 32806073 DOI: 10.1021/acsnano.0c04736] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Multiple experiments provide evidence for photovoltaic, catalytic, optoelectronic, and plasmonic processes involving hot, i.e., high energy, electrons in nanoscale materials. However, the mechanisms of such processes remain elusive, because electrons rapidly lose energy by relaxation through dense manifolds of states. We demonstrate a long-lived hot electron state in a Pt nanocluster adsorbed on the MoS2 substrate. For this purpose, we develop a simulation technique, combining classical molecular dynamics based on machine learning potentials with ab initio nonadiabatic molecular dynamics and real-time time-dependent density functional theory. Choosing Pt20/MoS2 as a prototypical system, we find frequent shifting of a top atom in the Pt particle occurring on a 50 ps time scale. The distortion breaks particle symmetry and creates unsaturated chemical bonds. The lifetime of the localized state associated with the broken bonds is enhanced by a factor of 3. Hot electrons aggregate near the shifted atom and form a catalytic reaction center. Our findings prove that distortion of even a single atom can have important implications for nanoscale catalysis and plasmonics and provide insights for utilizing machine learning potentials to accelerate ab initio investigations of excited state dynamics in condensed matter systems.
Collapse
Affiliation(s)
- Weibin Chu
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| | - Wissam A Saidi
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Oleg V Prezhdo
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
38
|
Agrawal S, Lin W, Prezhdo OV, Trivedi DJ. Ab initio quantum dynamics of charge carriers in graphitic carbon nitride nanosheets. J Chem Phys 2020; 153:054701. [PMID: 32770911 DOI: 10.1063/5.0010628] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Graphitic carbon nitride (g-C3N4), a metal-free and visible light responsive photocatalyst, has garnered much attention due to its wide range of applications. In order to elucidate the role of dimensionality on the properties of photo-generated charge carriers, we apply nonadiabatic (NA) molecular dynamics combined with time-domain density functional theory to investigate nonradiative relaxation of hot electrons and holes, and electron-hole recombination in monolayer and bulk g-C3N4. The nonradiative charge recombination occurs on a nanosecond timescale and is faster in bulk than the nanosheet, in agreement with the experiment. The difference arises due to the smaller energy gap and participation of additional vibrations in the bulk system. The long carrier lifetimes are favored by small NA coupling and rapid phonon-induced loss of quantum coherence between the excited and ground electronic states. Decoherence is fast because g-C3N4 is soft and undergoes large scale vibrations. The NA coupling is small since electrons and holes are localized on different atoms, and the electron-hole overlap is relatively small. Phonon-driven relaxation of hot electrons and holes takes 100-200 fs and is slightly slower at higher initial energies due to participation of fewer vibrational modes. This feature of two-dimensional g-C3N4 contrasts traditional three-dimensional semiconductors, which exhibit faster relaxation at higher energies due to larger density of states, and can be used to extract hot carriers to perform useful functions. The ab initio quantum dynamics simulations present a comprehensive picture of the photo-induced charge carrier dynamics in g-C3N4, guiding design of photovoltaic and photocatalytic devices.
Collapse
Affiliation(s)
- Sraddha Agrawal
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Wei Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350108, People's Republic of China
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Dhara J Trivedi
- Department of Physics, Clarkson University, Potsdam, New York 13699, USA
| |
Collapse
|
39
|
He J, Fang WH, Long R, Prezhdo OV. Why Oxygen Increases Carrier Lifetimes but Accelerates Degradation of CH3NH3PbI3 under Light Irradiation: Time-Domain Ab Initio Analysis. J Am Chem Soc 2020; 142:14664-14673. [DOI: 10.1021/jacs.0c06769] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Jinlu He
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei-Hai Fang
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - Oleg V. Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
40
|
Wei Y, Tokina MV, Benderskii AV, Zhou Z, Long R, Prezhdo OV. Quantum dynamics origin of high photocatalytic activity of mixed-phase anatase/rutile TiO2. J Chem Phys 2020; 153:044706. [DOI: 10.1063/5.0014179] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yaqing Wei
- College of Chemistry, Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Marina V. Tokina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alexander V. Benderskii
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Zhaohui Zhou
- Chemical Engineering and Technology, School of Environmental Science and Engineering, Key Laboratory of Subsurface Hydrology and Ecological Effects in Arid Region, Ministry of Education, Chang’an University, Xi’an 710064, China
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Oleg V. Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
41
|
Esch MP, Levine BG. State-pairwise decoherence times for nonadiabatic dynamics on more than two electronic states. J Chem Phys 2020; 152:234105. [PMID: 32571062 DOI: 10.1063/5.0010081] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Independent trajectory (IT) nonadiabatic molecular dynamics simulation methods are powerful tools for modeling processes involving transitions between electronic states. Incorporation and refinement of decoherence corrections into popular IT methods, e.g., Ehrenfest dynamics and trajectory surface hopping, is an important means of improving their accuracies. In this work, we identify a new challenge in the development of such decoherence corrections; when a system exists in a coherent superposition of three or more electronic states, coherences may decay unphysically when the decoherence correction is based on decoherence times assigned on a state-wise basis. As a solution, we introduce decoherence corrected Ehrenfest schemes based on decoherence times assigned on a state-pairwise basis. By application of these methods to a set of very simple one-dimensional model problems, we show that one of these state-pairwise methods ("collapse to a block") correctly describes the loss of coherence between all pairs of states in our multistate model problems, whereas a method based on a state-wise description of coherence loss does not. The new one-dimensional models introduced here can serve as useful tests for other decoherence correction schemes.
Collapse
Affiliation(s)
- Michael P Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
42
|
He J, Casanova D, Fang WH, Long R, Prezhdo OV. MAI Termination Favors Efficient Hole Extraction and Slow Charge Recombination at the MAPbI 3/CuSCN Heterojunction. J Phys Chem Lett 2020; 11:4481-4489. [PMID: 32423207 DOI: 10.1021/acs.jpclett.0c01467] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Photoinduced charge separation is the key step determining the efficiency of photon-to-electron conversion in solar cells, while charge carrier lifetimes govern the overall solar cell performance. Experiments report that copper(I) thiocyanate (CuSCN) is a very promising hole extraction layer for perovskite solar cells. Using nonadiabatic molecular dynamics combined with ab initio time-domain density functional theory, we show that termination of CH3NH3PbI3 (MAPbI3) at MAPbI3/CuSCN heterojunctions has a strong influence on both charge separation and recombination. Both processes are favored by MAI termination, compared to PbI2 termination. Because the MAPbI3 valence band originates from iodine orbitals while the conduction band arises from Pb orbitals, MAI termination places holes close to CuSCN, favoring extraction, and creates an MAI barrier for recombination of electrons in MAPbI3 and holes in CuSCN. The opposite is true for PbI2 termination. The origin of these effects is attributed solely to the properties of the MAPbI3 surfaces, and therefore, the conclusions should apply to other hole-transporting materials and can be generalized to other perovskites. Importantly, the simulations show that the injected hole remains hot for several hundreds of femtoseconds, allowing it to escape the interfacial region and prevent formation of bound excitons. This study suggests that metal halide perovskites should be treated with an organic precursor, such as MAI, prior to the formation of their interfaces with hole-transporting materials. The reported results advance the fundamental understanding of the highly unusual properties of metal halide perovskites and provide specific guidelines for optimizing the performance of perovskite solar cells and other devices.
Collapse
Affiliation(s)
- Jinlu He
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - David Casanova
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain
- IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Euskadi, Spain
| | - Wei-Hai Fang
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, P. R. China
| | - Oleg V Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
43
|
Binder R, Bonfanti M, Lauvergnat D, Burghardt I. First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian. J Chem Phys 2020; 152:204119. [PMID: 32486686 DOI: 10.1063/5.0004510] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].
Collapse
Affiliation(s)
- Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Matteo Bonfanti
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR 8000, 91405 Orsay, France
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
44
|
Binder R, Burghardt I. First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. II. ML-MCTDH simulations of exciton dynamics at a torsional defect. J Chem Phys 2020; 152:204120. [DOI: 10.1063/5.0004511] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
45
|
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
|
46
|
Smith B, Akimov AV. Hot Electron Cooling in Silicon Nanoclusters via Landau-Zener Nonadiabatic Molecular Dynamics: Size Dependence and Role of Surface Termination. J Phys Chem Lett 2020; 11:1456-1465. [PMID: 31958367 DOI: 10.1021/acs.jpclett.9b03687] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We develop a new express methodology for modeling excited-state dynamics occurring in dense manifolds of electronic states in atomistic systems. The approach leverages a modified Landau-Zener formula, the neglect of a back-reaction approximation, and the highly efficient density functional tight-binding method. We study the hot electron dynamics in a series of H- and F-terminated silicon nanocrystals (NCs) containing up to several hundred atoms. We explain the slower electron cooling dynamics in F-terminated NCs by the larger energy gaps between the adjacent electronic states in these systems as well as their slower fluctuations. We conclude that both the mass and chemical identity of the surface termination groups equally influence the electron dynamics, on average. However, the mass effect becomes dominant for higher-energy excitations. We find that the electron decay dynamics in F-terminated NCs has a greater sensitivity to the mass of the surface ligands than do the H-terminated NCs and explain this observation by the details of the electron-phonon coupling in the systems. We find that in the H-terminated NCs, electronic transitions in the cooling process occur predominantly between the surface states, whereas in F-terminated Si NCs, both surface and NC core states are coupled to the nuclear vibrations. We find that electron energy relaxation is accelerated in larger NCs and attribute this effect to the higher densities of states and smaller energy gaps in these systems.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry , University at Buffalo, The State University of New York , Buffalo , New York 14260 , United States
| | - Alexey V Akimov
- Department of Chemistry , University at Buffalo, The State University of New York , Buffalo , New York 14260 , United States
| |
Collapse
|
47
|
Zhang Z, Qiao L, Mora-Perez C, Long R, Prezhdo OV. Pb dimerization greatly accelerates charge losses in MAPbI3: Time-domain ab initio analysis. J Chem Phys 2020; 152:064707. [DOI: 10.1063/1.5131342] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Zhaosheng Zhang
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Lu Qiao
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Carlos Mora-Perez
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Oleg V. Prezhdo
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
48
|
Shende P, Wakade VS. Biointerface: a nano-modulated way for biological transportation. J Drug Target 2020; 28:456-467. [DOI: 10.1080/1061186x.2020.1720218] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Pravin Shende
- Shobhaben Pratapbhai Patel School of Pharmacy and Technology Management, SVKM’s NMIMS, Mumbai, India
| | - Varun S. Wakade
- Shobhaben Pratapbhai Patel School of Pharmacy and Technology Management, SVKM’s NMIMS, Mumbai, India
| |
Collapse
|
49
|
Hsu CP. Reorganization energies and spectral densities for electron transfer problems in charge transport materials. Phys Chem Chem Phys 2020; 22:21630-21641. [DOI: 10.1039/d0cp02994g] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Various contributions to the outer reorganization energy of an electron transfer system and their theoretical and computational aspects have been discussed.
Collapse
Affiliation(s)
- Chao-Ping Hsu
- 128 Academia Road Section 2
- Institute of Chemistry
- Academia Sinica
- Taipei
- Taiwan
| |
Collapse
|
50
|
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]
|