1
|
Brady RP, Drury C, Yurchenko SN, Tennyson J. Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules. J Chem Theory Comput 2024; 20:2127-2139. [PMID: 38171539 PMCID: PMC10938500 DOI: 10.1021/acs.jctc.3c01150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 12/19/2023] [Accepted: 12/20/2023] [Indexed: 01/05/2024]
Abstract
The (time-independent) Schrödinger equation for atomistic systems is solved by using the adiabatic potential energy curves (PECs) and the associated adiabatic approximation. In cases where interactions between electronic states become important, the associated nonadiabatic effects are taken into account via derivative couplings (DDRs), also known as nonadiabatic couplings (NACs). For diatomic molecules, the corresponding PECs in the adiabatic representation are characterized by avoided crossings. The alternative to the adiabatic approach is the diabatic representation obtained via a unitary transformation of the adiabatic states by minimizing the DDRs. For diatomics, the diabatic representation has zero DDR and nondiagonal diabatic couplings ensue. The two representations are fully equivalent and so should be the rovibronic energies and wave functions, which result from the solution of the corresponding Schrödinger equations. We demonstrate (for the first time) the numerical equivalence between the adiabatic and diabatic rovibronic calculations of diatomic molecules using the ab initio curves of yttrium oxide (YO) and carbon monohydride (CH) as examples of two-state systems, where YO is characterized by a strong NAC, while CH has a strong diabatic coupling. Rovibronic energies and wave functions are computed using a new diabatic module implemented in the variational rovibronic code Duo. We show that it is important to include both the diagonal Born-Oppenheimer correction and nondiagonal DDRs. We also show that the convergence of the vibronic energy calculations can strongly depend on the representation of nuclear motion used and that no one representation is best in all cases.
Collapse
Affiliation(s)
- Ryan P. Brady
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Charlie Drury
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Sergei N. Yurchenko
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| |
Collapse
|
2
|
Hetherington CV, Mohan T M N, Tilluck RW, Beck WF, Levine BG. Origin of Vibronic Coherences During Carrier Cooling in Colloidal Quantum Dots. J Phys Chem Lett 2023; 14:11651-11658. [PMID: 38109055 DOI: 10.1021/acs.jpclett.3c02384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Recent two-dimensional electronic spectroscopy experiments [Tilluck et al. J. Phys. Chem. Lett. 2021, 12 (39), 9677-9683] indicate the creation of coherent vibronic wavepackets in the first femtoseconds of hot carrier cooling in hexadecylamine-passivated CdSe quantum dots. Here we present a quantum chemical study of the origin of these coherences in a CdSe nanocrystal. We find that coherent wavepacket motions along vibrational coordinates with alkylamine character promote nonradiative relaxation through conical intersections between the exciton states of the inorganic core. Electronic excitations in the core are found to pass energy to the vibrations of the ligands via two distinct mechanisms: excitation of core phonon modes that are coupled to the ligand vibrations and direct excitation of ligand vibrations by delocalization of the exciton onto the ligands, both of which naturally arise within a photochemical framework based on many-electron potential energy surfaces. If these findings are demonstrated to be general, vibronic coherences may be leveraged to control photophysical outcomes in colloidal quantum dots.
Collapse
Affiliation(s)
- Caitlin V Hetherington
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University Stony Brook, New York 11733 United States
| | - Nila Mohan T M
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824 United States
| | - Ryan W Tilluck
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824 United States
| | - Warren F Beck
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824 United States
| | - Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University Stony Brook, New York 11733 United States
| |
Collapse
|
3
|
Akher FB, Shu Y, Varga Z, Truhlar DG. Semiclassical Multistate Dynamics for Six Coupled 5A' States of O + O 2. J Chem Theory Comput 2023. [PMID: 37441750 DOI: 10.1021/acs.jctc.3c00517] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Dynamics simulations of high-energy O2-O collisions play an important role in simulating thermal energy content and heat flux in flows around hypersonic vehicles. To carry out such dynamics simulations efficiently requires accurate global potential energy surfaces and (in most algorithms) state couplings for many energetically accessible electronic states. The ability to treat collisions involving many coupled electronic states has been a challenge for decades. Very recently, a new diabatization method, the parametrically managed diabatization by deep neural network (PM-DDNN), has been developed. The PM-DDNN method uses a deep neural network architecture with an activation function parametrically dependent on input data to discover and fit the diabatic potential energy matrix (DPEM) as a function of geometry, and the adiabatic potential energy surfaces are obtained by diagonalization of a small matrix with analytic matrix elements. Here, we applied the PM-DDNN method to the six lowest-energy potential energy surfaces in the 5A' manifold of O3 to perform simultaneous diabatization and fitting; the data are obtained by extended multistate complete-active-space second-order perturbation theory. We then used the adiabatic surfaces for dynamics calculations with three methods: coherent switching with decay of mixing (CSDM), curvature-driven CSDM (κCSDM), and electronically curvature-driven CSDM (eκCSDM). The κCSDM calculations require only adiabatic potential energies and gradients. The three dynamical methods are in good agreement. We then calculated electronically nonadiabatic, electronically inelastic, and dissociative cross sections for seven initial collision energies, five initial vibrational levels, and four initial rotational levels. Trends in the electronically inelastic cross sections as functions of the initial collision energy and vibrational level were rationalized in terms of the coordinate ranges where the gaps between the second and third potential energy surfaces are small.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
4
|
Shu Y, Zhang L, Wu D, Chen X, Sun S, Truhlar DG. New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States. J Chem Theory Comput 2023; 19:2419-2429. [PMID: 37079755 DOI: 10.1021/acs.jctc.2c01173] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Dihua Wu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Xiye Chen
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
5
|
Zhao X, Shu Y, Zhang L, Xu X, Truhlar DG. Direct Nonadiabatic Dynamics of Ammonia with Curvature-Driven Coherent Switching with Decay of Mixing and with Fewest Switches with Time Uncertainty: An Illustration of Population Leaking in Trajectory Surface Hopping Due to Frustrated Hops. J Chem Theory Comput 2023; 19:1672-1685. [PMID: 36877830 DOI: 10.1021/acs.jctc.2c01260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2023]
Abstract
Mixed quantum-classical nonadiabatic dynamics is a widely used approach to simulate molecular dynamics involving multiple electronic states. There are two main categories of mixed quantum-classical nonadiabatic dynamics algorithms, namely, trajectory surface hopping (TSH) in which the trajectory propagates on a single potential energy surface, interrupted by hops, and self-consistent-potential (SCP) methods, such as semiclassical Ehrenfest, in which propagation occurs on a mean-field surface without hops. In this work, we will illustrate an example of severe population leaking in TSH. We emphasize that such leaking is a combined effect of frustrated hops and long-time simulations that drive the final excited-state population toward zero as a function of time. We further show that such leaking can be alleviated-but not eliminated-by the fewest switches with time uncertainty TSH algorithm (here implemented in the SHARC program); the time uncertainty algorithm slows down the leaking process by a factor of 4.1. The population leaking is not present in coherent switching with decay of mixing (CSDM), which is an SCP method with non-Markovian decoherence included. Another result in this paper is that we find very similar results with the original CSDM algorithm, with time-derivative CSDM (tCSDM), and with curvature-driven CSDM (κCSDM). Not only do we find good agreement for electronically nonadiabatic transition probabilities but also we find good agreement of the norms of the effective nonadiabatic couplings (NACs) that are derived from the curvature-driven time-derivative couplings as implemented in κCSDM with the time-dependent norms of the nonadiabatic coupling vectors computed by state-averaged complete-active-space self-consistent field theory.
Collapse
Affiliation(s)
- Xiaorui Zhao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China.,School of Aerospace Engineering, Tsinghua University, Beijing 100084, P. R. China
| | - Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xuefei Xu
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
6
|
Boeije Y, Olivucci M. From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions. Chem Soc Rev 2023; 52:2643-2687. [PMID: 36970950 DOI: 10.1039/d2cs00719c] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
Abstract
This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.
Collapse
Affiliation(s)
- Yorrick Boeije
- Van 't Hoff Institute for Molecular Sciences (HIMS), University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Massimo Olivucci
- Chemistry Department, University of Siena, Via Aldo Moro n. 2, 53100 Siena, Italy
- Chemistry Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, USA
| |
Collapse
|
7
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
8
|
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: 4] [Impact Index Per Article: 2.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
|
9
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
10
|
Lingerfelt DB, Yu T, Yoshimura A, Ganesh P, Jakowski J, Sumpter BG. Nonadiabatic Effects on Defect Diffusion in Silicon-Doped Nanographenes. NANO LETTERS 2021; 21:236-242. [PMID: 33337886 DOI: 10.1021/acs.nanolett.0c03587] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Single atom impurities in graphene, substitutional silicon defects in particular, have been observed to diffuse under electron beam irradiation. However, the relative importance of elastic and inelastic scattering in facilitating their mobility remains unclear. Here, we employ excited-state electronic structure calculations to explore potential inelastic effects, and find an electronically nonadiabatic excited-state silicon diffusion pathway involving "softened" Si-C bonding that presents an ∼2 eV lower diffusion barrier than the ground-state pathway. Beam-induced transition rates to this state indicate that the excited-state pathway is accessible through irradiation of the defect site. However, even in the limit of fully elastic scattering, upward nonadiabatic transitions are also possible along the diffusion coordinate, increasing the diffusion barrier and further demonstrating the potential for electronic nonadiabaticity to influence beam-induced atomic transformations in materials. We also propose some experimentally testable signatures of such excited-state pathways.
Collapse
Affiliation(s)
- David B Lingerfelt
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Tao Yu
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, United States
| | - Anthony Yoshimura
- Department of Physics, Applied Physics, and Astronomy, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Panchapakesan Ganesh
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Jacek Jakowski
- Computational Sciences and Engineering Division Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Bobby G Sumpter
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
11
|
Schaumann EN, Tian B. Biological Interfaces, Modulation, and Sensing with Inorganic Nano-Bioelectronic Materials. SMALL METHODS 2020; 4:1900868. [PMID: 34295965 PMCID: PMC8294120 DOI: 10.1002/smtd.201900868] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Accepted: 02/16/2020] [Indexed: 05/30/2023]
Abstract
The last several years have seen a large and increasing interest in scientific developments that combine methods and materials from nanotechnology with questions and applications in bioelectronics. This follows with a number of broader trends: the rapid increase in functionality for materials at the nanoscale; a growing recognition of the importance of electric fields in diverse physiological processes; and continuous improvements in technologies that are naturally complementary with bioelectronics, such as optogenetics. Here, a progress report is provided on several of the most exciting recent developments in this field. The three critical functions of biointerface formation, biological modulation, and biological sensing using newly developed nanoscale materials are considered.
Collapse
Affiliation(s)
- Erik N Schaumann
- Department of Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Bozhi Tian
- Department of Chemistry, The University of Chicago, Chicago, IL 60637, USA
| |
Collapse
|
12
|
Zhang L, Chu W, Zheng Q, Benderskii AV, Prezhdo OV, Zhao J. Suppression of Electron-Hole Recombination by Intrinsic Defects in 2D Monoelemental Material. J Phys Chem Lett 2019; 10:6151-6158. [PMID: 31553184 DOI: 10.1021/acs.jpclett.9b02620] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
The Shockley-Read-Hall (SRH) model, in which the deep trap defect states in the band gap are proposed as nonradiative electron-hole (e-h) recombination centers, has been widely used to describe the nonradiative e-h recombination through the defects in semiconductor. By using the ab initio nonadiabatic molecular dynamics method, we find that the SRH model fails to describe the e-h recombination behavior for defects in 2D monoelemental material such as monolayer black phosphorus (BP). Through the investigation of three intrinsic defects with shallow and deep defect states in monolayer BP, it is found that, surprisingly, none of these defects significantly accelerates the e-h recombination. Further analysis shows that because monolayer BP is a monoelemental material, the distinct impurity phonon, which often induces fast e-h recombination, is not formed. Moreover, because of the flexibility of 2D material, the defects scatter the phonons present in pristine BP, generating multiple modes with lower frequencies compared with the pristine BP, which further suppresses the e-h recombination. We propose that the conclusion can be extended to other monoelemental 2D materials, which is important guidance for the future design of functional semiconductors.
Collapse
Affiliation(s)
- Lili Zhang
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
- Department of Chemistry , University of Southern California , Los Angeles , California 90089 , United States
| | - Weibin Chu
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
- Department of Chemistry , University of Southern California , Los Angeles , California 90089 , United States
| | - Qijing Zheng
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
- Synergetic Innovation Center of Quantum Information & Quantum Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
| | - Alexander V Benderskii
- 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
| | - Jin Zhao
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
- Synergetic Innovation Center of Quantum Information & Quantum Physics , University of Science and Technology of China , Hefei , Anhui 230026 , China
| |
Collapse
|
13
|
Shu Y, Kryven J, Sampaio de Oliveira-Filho AG, Zhang L, Song GL, Li SL, Meana-Pañeda R, Fu B, Bowman JM, Truhlar DG. Direct diabatization and analytic representation of coupled potential energy surfaces and couplings for the reactive quenching of the excited 2Σ+ state of OH by molecular hydrogen. J Chem Phys 2019; 151:104311. [DOI: 10.1063/1.5111547] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Joanna Kryven
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Antonio Gustavo Sampaio de Oliveira-Filho
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
- Departamento de Química, Laboratório Computacional de Espectroscopia e Cinética, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, 14040-901 Ribeirão Preto-SP, Brazil
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Guo-Liang Song
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Shaohong L. Li
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Rubén Meana-Pañeda
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Bina Fu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
14
|
Levine BG, Esch MP, Fales BS, Hardwick DT, Peng WT, Shu Y. Conical Intersections at the Nanoscale: Molecular Ideas for Materials. Annu Rev Phys Chem 2019; 70:21-43. [DOI: 10.1146/annurev-physchem-042018-052425] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ability to predict and describe nonradiative processes in molecules via the identification and characterization of conical intersections is one of the greatest recent successes of theoretical chemistry. Only recently, however, has this concept been extended to materials science, where nonradiative recombination limits the efficiencies of materials for various optoelectronic applications. In this review, we present recent advances in the theoretical study of conical intersections in semiconductor nanomaterials. After briefly introducing conical intersections, we argue that specific defects in materials can induce conical intersections between the ground and first excited electronic states, thus introducing pathways for nonradiative recombination. We present recent developments in theoretical methods, computational tools, and chemical intuition for the prediction of such defect-induced conical intersections. Through examples in various nanomaterials, we illustrate the significance of conical intersections for nanoscience. We also discuss challenges facing research in this area and opportunities for progress.
Collapse
Affiliation(s)
- Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - B. Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Dylan T. Hardwick
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
15
|
Levine BG, Peng WT, Esch MP. Locality of conical intersections in semiconductor nanomaterials. Phys Chem Chem Phys 2019; 21:10870-10878. [PMID: 31106323 DOI: 10.1039/c9cp01584a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A predictive theory connecting atomic structure to the rate of recombination would enable the rational design of semiconductor nanomaterials for optoelectronic applications. Recently our group has demonstrated that the theoretical study of conical intersections can serve this purpose. Here we review recent work in this area, focusing on the thesis that low-energy conical intersections in nanomaterials share a common feature: locality. We define a conical intersection as local if (a) the intersecting states differ by the excitation of an electron between spatially local orbitals, and (b) the intersection is accessed when the energies of these orbitals are tuned by local distortions of the geometry. After illustrating the locality of the conical intersection responsible for recombination at dangling bond defects in silicon, we demonstrate the locality of low-energy conical intersections in cases where locality may be a surprise. First, we demonstrate the locality of low-energy self-trapped conical intersections in a pristine silicon nanocrystal, which has no defects that one would expect to serve as the center of a local intersection. Second, we demonstrate that the lowest energy intersection in a silicon system with two neighboring dangling bond defects localizes to a single defect site. We discuss the profound implications of locality for predicting the rate of recombination and suggest that the locality of intersections could be exploited in the experimental study of recombination, where spectroscopic studies of molecular models of defects could provide new insights.
Collapse
Affiliation(s)
- Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Wei-Tao Peng
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| | - Michael P Esch
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
16
|
Esch MP, Shu Y, Levine BG. A Conical Intersection Perspective on the Low Nonradiative Recombination Rate in Lead Halide Perovskites. J Phys Chem A 2019; 123:2661-2673. [DOI: 10.1021/acs.jpca.9b00952] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
17
|
Fedorov DA, Levine BG. A discontinuous basis enables numerically exact solution of the Schrödinger equation around conical intersections in the adiabatic representation. J Chem Phys 2019; 150:054102. [PMID: 30736673 DOI: 10.1063/1.5058268] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Solving the vibrational Schrödinger equation in the neighborhood of conical intersections in the adiabatic representation is a challenge. At the intersection point, first- and second-derivative nonadiabatic coupling matrix elements become singular, with the singularity in the second-derivative coupling (diagonal Born-Oppenheimer correction) being non-integrable. These singularities result from discontinuities in the vibronic functions associated with the individual adiabatic states, and our group has recently argued that these divergent matrix elements cancel when discontinuous adiabatic vibronic functions sum to a continuous total nonadiabatic wave function. Here we describe the realization of this concept: a novel scheme for the numerically exact solution of the Schrödinger equation in the adiabatic representation. Our approach is based on a basis containing functions that are discontinuous at the intersection point. We demonstrate that the individual adiabatic nuclear wave functions are themselves discontinuous at the intersection point. This proves that discontinuous basis functions are essential to any tractable method that solves the Schrödinger equation around conical intersections in the adiabatic representation with high numerical precision. We establish that our method provides numerically exact results by comparison to reference calculations performed in the diabatic representation. In addition, we quantify the energetic error associated with constraining the density to be zero at the intersection point, a natural approximation. Prospects for extending the present treatment of a two-dimensional model to systems of higher dimensionality are discussed.
Collapse
Affiliation(s)
- Dmitry A Fedorov
- 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
|
18
|
Zhang X, Carter EA. Subspace Density Matrix Functional Embedding Theory: Theory, Implementation, and Applications to Molecular Systems. J Chem Theory Comput 2018; 15:949-960. [DOI: 10.1021/acs.jctc.8b00990] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Zhang
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - Emily A. Carter
- School of Engineering and Applied Science, Princeton University, Princeton, New Jersey 08544-5263, United States
| |
Collapse
|
19
|
Improved potential energy surfaces of thioanisole and the effect of upper surface variations on the product distribution upon photodissociation. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.06.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
20
|
Han Y, Anderson K, Hobbie EK, Boudjouk P, Kilin DS. Unraveling Photodimerization of Cyclohexasilane from Molecular Dynamics Studies. J Phys Chem Lett 2018; 9:4349-4354. [PMID: 30004709 DOI: 10.1021/acs.jpclett.8b01691] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Photoinduced reactions of a pair of cyclohexasilane (CHS) monomers are explored by time-dependent excited-state molecular dynamics (TDESMD) calculations. In TDESMD trajectories, one observes vivid reaction events including dimerization and fragmentation. A general reaction pathway is identified as (i) ring-opening formation of a dimer, (ii) rearrangement induced by bond breaking, and (iii) decomposition through the elimination of small fragments. The identified pathway supports the chemistry proposed for the fabrication of silicon-based materials using CHS as a precursor. In addition, we find dimers have smaller HOMO-LUMO gaps and exhibit a red shift and line-width broadening in the computed photoluminescence spectra compared with a pair of CHS monomers.
Collapse
Affiliation(s)
- Yulun Han
- Department of Chemistry and Biochemistry , North Dakota State University , Fargo , North Dakota 58102 , United States
| | - Kenneth Anderson
- Department of Chemistry and Biochemistry , North Dakota State University , Fargo , North Dakota 58102 , United States
| | - Erik K Hobbie
- Department of Chemistry and Biochemistry , North Dakota State University , Fargo , North Dakota 58102 , United States
| | - Philip Boudjouk
- Department of Chemistry and Biochemistry , North Dakota State University , Fargo , North Dakota 58102 , United States
| | - Dmitri S Kilin
- Department of Chemistry and Biochemistry , North Dakota State University , Fargo , North Dakota 58102 , United States
| |
Collapse
|
21
|
Peng WT, Fales BS, Shu Y, Levine BG. Dynamics of recombination via conical intersection in a semiconductor nanocrystal. Chem Sci 2018; 9:681-687. [PMID: 29629136 PMCID: PMC5869574 DOI: 10.1039/c7sc04221c] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Accepted: 11/13/2017] [Indexed: 12/17/2022] Open
Abstract
Conical intersections are well known to introduce nonradiative decay pathways in molecules, but have only recently been implicated in nonradiative recombination processes in materials. Here we apply excited state ab initio molecular dynamics simulations based on a multireference description of the electronic structure to defective silicon nanocrystals up to 1.7 nm in diameter to search for accessible nonradiative recombination pathways. Dangling bond defects are found to induce conical intersections between the ground and first excited electronic states of five systems of various sizes. These defect-induced conical intersections are accessible at energies that are in the visible range (2.4-2.7 eV) and very weakly dependent on particle size. The dynamic simulations suggest that these intersections are accessed 40-60 fs after creation of a defect-localized excitation. This ultrafast recombination is attributed to the fact that Jahn-Teller distortion on the first excited state drives the defect directly towards a conical intersection with the ground electronic state.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department of Chemistry , Michigan State University , East Lansing , MI 48824 , USA .
| | - B Scott Fales
- Department of Chemistry , The PULSE Institute , Stanford University , Stanford , CA 94305 , USA.,SLAC National Accelerator Laboratory , Menlo Park , CA 94025 , USA
| | - Yinan Shu
- Department of Chemistry , University of Minnesota , Minneapolis , MN 55455 , USA
| | - Benjamin G Levine
- Department of Chemistry , Michigan State University , East Lansing , MI 48824 , USA .
| |
Collapse
|
22
|
Shu Y, Parker KA, Truhlar DG. Dual-Functional Tamm–Dancoff Approximation with Self-Interaction-Free Orbitals: Vertical Excitation Energies and Potential Energy Surfaces near an Intersection Seam. J Phys Chem A 2017; 121:9728-9735. [DOI: 10.1021/acs.jpca.7b11400] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Kelsey A. Parker
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|