1
|
Forde A, Maity S, Freixas VM, Fernandez-Alberti S, Neukirch AJ, Kleinekathöfer U, Tretiak S. Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers. J Phys Chem Lett 2024; 15:4142-4150. [PMID: 38593451 DOI: 10.1021/acs.jpclett.4c00362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.
Collapse
Affiliation(s)
- Aaron Forde
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sayan Maity
- School of Science, Constructor University, Campus Ring 1, 28759 Bremen, Germany
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologiia, Univresidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | | | - Amanda J Neukirch
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
2
|
Sun Z, Li S, Meng Y, An Z. Surface hopping simulations on charge photogeneration in conjugated polymers. Phys Chem Chem Phys 2024; 26:10982-10988. [PMID: 38526436 DOI: 10.1039/d3cp06013f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2024]
Abstract
The mechanism of charge photogeneration in neat conjugated polymers has long been controversial and a unified explanation has not been achieved so far. In this paper, we use a surface hopping method to simulate the excited-state dynamics of a system composed of five π-stacked polymer chains. In this system, polaron pairs (PPs) with large electron-hole distance can be seen as free charges. The surface hopping method is based on the Pariser-Parr-Pople (PPP) Hamiltonian and the excited states of the system are calculated with configuration interaction singles (CIS) formalism. During the simulations, the yields of PPs and free charges are calculated using a statistical method. By comparison, it is found that impurity and excess energy have significant effects on the yields of PPs and free charges. Free charges are difficult to be generated in neat systems with small excess energy. Free charges come from direct dissociation of high-energy excitons.
Collapse
Affiliation(s)
- Zhen Sun
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China.
- Zhejiang Institute of Photoelectronics & Zhejiang Institute for Advanced Light Source, Zhejiang Normal University, Jinhua 321004, China
| | - Sheng Li
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China.
- Zhejiang Institute of Photoelectronics & Zhejiang Institute for Advanced Light Source, Zhejiang Normal University, Jinhua 321004, China
| | - Yan Meng
- Department of Physics, Xingtai University, Xingtai 054001, China
| | - Zhong An
- College of Physics, Hebei Normal University, Shijiazhuang 050024, China.
| |
Collapse
|
3
|
Tao Z, Qiu T, Bhati M, Bian X, Duston T, Rawlinson J, Littlejohn RG, Subotnik JE. Practical phase-space electronic Hamiltonians for ab initio dynamics. J Chem Phys 2024; 160:124101. [PMID: 38526114 DOI: 10.1063/5.0192084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 02/14/2024] [Indexed: 03/26/2024] Open
Abstract
Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Jonathan Rawlinson
- Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
4
|
Vandaele E, Mališ M, Luber S. A Local Diabatisation Method for Two-State Adiabatic Conical Intersections. J Chem Theory Comput 2024; 20:856-872. [PMID: 38174710 DOI: 10.1021/acs.jctc.3c01008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
5
|
Götze JP, Lokstein H. Excitation Energy Transfer between Higher Excited States of Photosynthetic Pigments: 2. Chlorophyll b is a B Band Excitation Trap. ACS OMEGA 2023; 8:40015-40023. [PMID: 37929150 PMCID: PMC10620878 DOI: 10.1021/acsomega.3c05896] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 09/21/2023] [Indexed: 11/07/2023]
Abstract
Chlorophylls (Chls) are known for fast, subpicosecond internal conversion (IC) from ultraviolet/blue absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates such as photosynthetic pigment-protein complexes. We show that according to a Förster resonance energy transport (FRET) scheme, unmodulated B-B EET has an unexpectedly high range. Unsuppressed, it could pose an existential threat-the damage potential of blue light for photochemical reaction centers (RCs) is well-known. This insight reveals so-far undescribed roles for carotenoids (Crts, cf. previous article in this series) and Chl b (this article) of possibly vital importance. Our model system is the photosynthetic antenna pigment-protein complex (CP29). The focus of the study is on the role of Chl b for EET in the Q and B bands. Further, the initial excited pigment distribution in the B band is computed for relevant solar irradiation and wavelength-centered laser pulses. It is found that both accessory pigment classes compete efficiently with Chl a absorption in the B band, leaving only 40% of B band excitations for Chl a. B state population is preferentially relocated to Chl b after excitation of any Chls, due to a near-perfect match of Chl b B band absorption with Chl a B state emission spectra. This results in an efficient depletion of the Chl a population (0.66 per IC/EET step, as compared to 0.21 in a Chl a-only system). Since Chl b only occurs in the peripheral antenna complexes of plants and algae, and RCs contain only Chl a, this would automatically trap potentially dangerous B state population in the antennae, preventing forwarding to the RCs.
Collapse
Affiliation(s)
- Jan P. Götze
- Institut
für Chemie und Biochemie, Fachbereich Biologie Chemie Pharmazie, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Heiko Lokstein
- Department
of Chemical Physics and Optics, Charles
University, Ke Karlovu
3, 121 16 Prague
2, Czech Republic
| |
Collapse
|
6
|
Mejia-Rodriguez D, Aprà E, Autschbach J, Bauman NP, Bylaska EJ, Govind N, Hammond JR, Kowalski K, Kunitsa A, Panyala A, Peng B, Rehr JJ, Song H, Tretiak S, Valiev M, Vila FD. NWChem: Recent and Ongoing Developments. J Chem Theory Comput 2023; 19:7077-7096. [PMID: 37458314 DOI: 10.1021/acs.jctc.3c00421] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff R Hammond
- Accelerated Computing, NVIDIA Helsinki Oy, Porkkalankatu 1, 00180 Helsinki, Finland
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marat Valiev
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
Sun Z, Li S, Xie S, Meng Y, An Z. Surface hopping simulations on charge separation in an organic donor-acceptor system. Phys Chem Chem Phys 2023; 25:26203-26210. [PMID: 37740356 DOI: 10.1039/d3cp02164e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/24/2023]
Abstract
Charge separation in organic solar cells is a long-lasting and heavily debated issue. Here, we use the surface hopping method based on the Pariser-Parr-Pople (PPP) Hamiltonian and configuration interaction singles (CIS) approximation to simulate the charge separation process in an organic donor-acceptor system. The system is composed of one donor polymer chain and four acceptor polymer chains, and they are all stacked face-to-face. We let the system to relax from a photoexcited state, and then we observed that the electron is transferred from the donor chain to different acceptor chains and the hole is left on the donor chain, forming polaron pairs with different electron-hole distances. By performing statistical analysis on a number of trajectories, we found that the electron and the hole are fully separated before the system relaxes to its lowest excited state. The yield of free charges shows a significant dependence on the donor-acceptor band offset which provides the driving force for charge separation, while showing negligible dependence on the photoexcitation energy. The external electric field has a remarkable effect on the yield of free charges.
Collapse
Affiliation(s)
- Zhen Sun
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China.
- Zhejiang Institute of Photoelectronics & Zhejiang Institute for Advanced Light Source, Zhejiang Normal University, Jinhua 321004, China
| | - Sheng Li
- Department of Physics, Zhejiang Normal University, Jinhua 321004, China.
- Zhejiang Institute of Photoelectronics & Zhejiang Institute for Advanced Light Source, Zhejiang Normal University, Jinhua 321004, China
| | - Shijie Xie
- School of Physics, State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, China
| | - Yan Meng
- Department of Physics, Xingtai University, Xingtai 054001, China
| | - Zhong An
- College of Physics, Hebei Normal University, Shijiazhuang 050024, China.
| |
Collapse
|
8
|
Freixas VM, Malone W, Li X, Song H, Negrin-Yuvero H, Pérez-Castillo R, White A, Gibson TR, Makhov DV, Shalashilin DV, Zhang Y, Fedik N, Kulichenko M, Messerly R, Mohanam LN, Sharifzadeh S, Bastida A, Mukamel S, Fernandez-Alberti S, Tretiak S. NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2023; 19:5356-5368. [PMID: 37506288 DOI: 10.1021/acs.jctc.3c00583] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2023]
Abstract
We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.
Collapse
Affiliation(s)
- Victor M Freixas
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | - Walter Malone
- Department of Physics, Tuskegee University, Tuskegee, Alabama 36088, United States
| | - Xinyang Li
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Royle Pérez-Castillo
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Alexander White
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie R Gibson
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | | - Yu Zhang
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nikita Fedik
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Maksim Kulichenko
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Richard Messerly
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Luke Nambi Mohanam
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Sahar Sharifzadeh
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - Shaul Mukamel
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | | | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
9
|
Bonilla V, Freixas VM, Fernandez-Alberti S, Galindo JF. Impact of the core on the inter-branch exciton exchange in dendrimers. Phys Chem Chem Phys 2023; 25:12097-12106. [PMID: 37133823 DOI: 10.1039/d2cp06009d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Organic dendrimers with π conjugated systems are capable of capturing solar energy as a renewable source for human use. Nonetheless, further study regarding the relationship between the structure and the energy transfer mechanism in these types of molecules is still necessary. In this work, nonadiabatic excited state molecular dynamics (NEXMD) were carried out to study the intra- and inter-branch exciton migration in two tetra-branched dendrimers, C(dSSB)4 and Ad(BuSSB)4, which differ in their respective carbon and adamantane core. Both systems undergo a ladder decay mechanism between excited states, with back-and-forth transitions between S1 and S2. Despite presenting very similar absorption-emission spectra, differences in the photoinduced energy relaxation are observed. The size of the core impacts the inter-branch energy exchange and transient exciton localization/delocalization, which ultimately condition the relative energy relaxation rates, being faster in Ad(BuSSB)4 with respect to C(dSSB)4. Nevertheless, the photoinduced processes lead to a progressive final exciton-self-trapping in one of the branches of both dendrimers, which is a desirable feature in organic photovoltaic applications. Our results can inspire the design of more efficient dendrimers with the desired magnitude of inter-branch exciton exchange and localization/delocalization according to changes in their core.
Collapse
Affiliation(s)
- Valeria Bonilla
- Departamento de Química, Universidad Nacional de Colombia, Sede Bogotá, 111321, Bogotá, Colombia.
| | - Victor M Freixas
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | | | - Johan Fabian Galindo
- Departamento de Química, Universidad Nacional de Colombia, Sede Bogotá, 111321, Bogotá, Colombia.
| |
Collapse
|
10
|
Niehaus TA. Exact non-adiabatic coupling vectors for the time-dependent density functional based tight-binding method. J Chem Phys 2023; 158:054103. [PMID: 36754796 DOI: 10.1063/5.0136838] [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/15/2023] Open
Abstract
We report on non-adiabatic coupling vectors between electronic excited states for the time-dependent-density functional theory based tight-binding (TD-DFTB) method. The implementation includes orbital relaxation effects that have been previously neglected and covers also the case of range-separated exchange-correlation functionals. Benchmark calculations with respect to first principles TD-DFT highlight the large dependence of non-adiabatic couplings on the functional. Closer investigations of the topology around a conical intersection between excited states show that TD-DFTB delivers near-exact values of the Berry phase, which paves the way for consistent non-adiabatic molecular dynamics simulations for large systems.
Collapse
Affiliation(s)
- Thomas A Niehaus
- University Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| |
Collapse
|
11
|
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
|
12
|
Nelson T, Huestis PL, Manner VW. Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes. J Phys Chem A 2022; 126:7094-7101. [PMID: 36196028 PMCID: PMC9574918 DOI: 10.1021/acs.jpca.2c03404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The photolytic stability of explosives and energetic
functional
groups is of importance for those who regularly handle or are exposed
to explosives in typical environmental conditions. This study models
the photolytic degradation of dodecane substituted with various energetic
functional groups: azide, nitro, nitrate ester, and nitramine. For
the studied molecules, it was found that excitons localize on the
energetic functional group, no matter where they were initially formed,
and thus, the predominant degradation pathway involves the degradation
of the energetic functional group. The relative trends for both 4
and 8 eV excitation energies followed with what is expected from the
relative stability of the energetic functional groups to thermal and
sub-shock degradation. The one notable exception was the azide functional
group; more work should be done to further understand the photolytic
effects on the azide functional group.
Collapse
Affiliation(s)
- Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Patricia L Huestis
- High Explosives Science & Technology, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Virginia W Manner
- High Explosives Science & Technology, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
13
|
Forde A, Freixas VM, Fernandez-Alberti S, Neukirch AJ, Tretiak S. Charge-Transfer Luminescence in a Molecular Donor-Acceptor Complex: Computational Insights. J Phys Chem Lett 2022; 13:8755-8760. [PMID: 36099248 DOI: 10.1021/acs.jpclett.2c02479] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Donor-acceptor molecular complexes are a popular class of materials utilizing charge-transfer states for practical applications. A recent class of donor-acceptor dyads based on the fluorescent BODIPY functionalized with triphenylamine (TPA) shows the peculiar property of dual fluorescence. It is hypothesized that instead of the sensitized charge-transfer state being optically dark, it provides an additional bright radiative pathway. Here we use time-dependent density functional theory to characterize the energetic alignment of excitonic and charge-transfer states in a BODIPY-TPA molecular complex. We observe that using a long-range exchange corrected functional in combination with state-specific solvation scheme gives a qualitatively correct alignment of the exciton and charge-transfer states and an enhancement in oscillator strength for the equilibrium solvated charge-transfer state, in agreement with experiment. This work provides rationalization of charge-transfer state emission and provides a foundation to explore charge-transfer using ab initio excited-state nonadiabatic dynamics.
Collapse
Affiliation(s)
- Aaron Forde
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologiia, Univresidad Nacional de Quilmes/CONICET, B1876BXD, Bernal, Argentina
| | | | - Amanda J Neukirch
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
14
|
Tracy DA, Fernandez-Alberti S, Tretiak S, Roitberg AE. Adiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation. J Chem Theory Comput 2022; 18:5213-5220. [PMID: 36044726 DOI: 10.1021/acs.jctc.2c00561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a method to link the Nonadiabatic EXcited-state Molecular Dynamics (NEXMD) package to the SANDER package supplied by AMBERTOOLS to provide excited-state adiabatic quantum mechanics/molecular mechanics (QM/MM) simulations. NEXMD is a computational package particularly developed to perform simulations of the photoexcitation and subsequent nonadiabatic electronic and vibrational energy relaxation in large multichromophoric conjugated molecules involving several coupled electronic excited states. The NEXMD-SANDER exchange has been optimized in order to achieve excited-state adiabatic dynamics simulations of large conjugated materials in a QM/MM environment, such as an explicit solvent. Dynamics of a substituted polyphenylene vinylene oligomer (PPV3-NO2) in vacuum and different explicit solvents has been used as a test case by performing comparative analysis of changes in its optical spectrum, state-dependent conformational changes, and quantum bond orderings. The method has been tested and compared with respect to previous implicit solvent implementations. Also, the impact on the expansion of the QM region by including a variable number of solvent molecules has been analyzed. Altogether, these results encourage future implementations of NEXMD simulations using the same combination of methods.
Collapse
Affiliation(s)
- Dustin A Tracy
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | | | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS) and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Adrian E Roitberg
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
15
|
Zheng Y, Han Y, Weight BM, Lin Z, Gifford BJ, Zheng M, Kilin D, Kilina S, Doorn SK, Htoon H, Tretiak S. Photochemical spin-state control of binding configuration for tailoring organic color center emission in carbon nanotubes. Nat Commun 2022; 13:4439. [PMID: 35915090 PMCID: PMC9343348 DOI: 10.1038/s41467-022-31921-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 07/04/2022] [Indexed: 12/04/2022] Open
Abstract
Incorporating fluorescent quantum defects in the sidewalls of semiconducting single-wall carbon nanotubes (SWCNTs) through chemical reaction is an emerging route to predictably modify nanotube electronic structures and develop advanced photonic functionality. Applications such as room-temperature single-photon emission and high-contrast bio-imaging have been advanced through aryl-functionalized SWCNTs, in which the binding configurations of the aryl group define the energies of the emitting states. However, the chemistry of binding with atomic precision at the single-bond level and tunable control over the binding configurations are yet to be achieved. Here, we explore recently reported photosynthetic protocol and find that it can control chemical binding configurations of quantum defects, which are often referred to as organic color centers, through the spin multiplicity of photoexcited intermediates. Specifically, photoexcited aromatics react with SWCNT sidewalls to undergo a singlet-state pathway in the presence of dissolved oxygen, leading to ortho binding configurations of the aryl group on the nanotube. In contrast, the oxygen-free photoreaction activates previously inaccessible para configurations through a triplet-state mechanism. These experimental results are corroborated by first principles simulations. Such spin-selective photochemistry diversifies SWCNT emission tunability by controlling the morphology of the emitting sites. Chemical functionalization of the sidewalls of single-wall carbon nanotubes (SWCNTs) is an emerging route to introduce fluorescent quantum defects and tailor the emission properties. Here, authors demonstrate that spin-selective photochemistry diversifies SWCNT emission tunability by controlling the morphology of the emitting sites.
Collapse
Affiliation(s)
- Yu Zheng
- Center for Integrated Nanotechnologies, Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Yulun Han
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND, 58102, USA
| | - Braden M Weight
- Center for Integrated Nanotechnologies, Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.,Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND, 58102, USA.,Department of Physics, North Dakota State University, Fargo, ND, 58102, USA.,Department of Physics and Astronomy, University of Rochester, Rochester, NY, 14627, USA
| | - Zhiwei Lin
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD, 20899, USA
| | - Brendan J Gifford
- Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Ming Zheng
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD, 20899, USA
| | - Dmitri Kilin
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND, 58102, USA
| | - Svetlana Kilina
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND, 58102, USA
| | - Stephen K Doorn
- Center for Integrated Nanotechnologies, Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Han Htoon
- Center for Integrated Nanotechnologies, Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Sergei Tretiak
- Center for Integrated Nanotechnologies, Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA. .,Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| |
Collapse
|
16
|
Ten Brink M, Gräber S, Hopjan M, Jansen D, Stolpp J, Heidrich-Meisner F, Blöchl PE. Real-time non-adiabatic dynamics in the one-dimensional Holstein model: Trajectory-based vs exact methods. J Chem Phys 2022; 156:234109. [PMID: 35732530 DOI: 10.1063/5.0092063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We benchmark a set of quantum-chemistry methods, including multitrajectory Ehrenfest, fewest-switches surface-hopping, and multiconfigurational-Ehrenfest dynamics, against exact quantum-many-body techniques by studying real-time dynamics in the Holstein model. This is a paradigmatic model in condensed matter theory incorporating a local coupling of electrons to Einstein phonons. For the two-site and three-site Holstein model, we discuss the exact and quantum-chemistry methods in terms of the Born-Huang formalism, covering different initial states, which either start on a single Born-Oppenheimer surface, or with the electron localized to a single site. For extended systems with up to 51 sites, we address both the physics of single Holstein polarons and the dynamics of charge-density waves at finite electron densities. For these extended systems, we compare the quantum-chemistry methods to exact dynamics obtained from time-dependent density matrix renormalization group calculations with local basis optimization (DMRG-LBO). We observe that the multitrajectory Ehrenfest method, in general, only captures the ultrashort time dynamics accurately. In contrast, the surface-hopping method with suitable corrections provides a much better description of the long-time behavior but struggles with the short-time description of coherences between different Born-Oppenheimer states. We show that the multiconfigurational Ehrenfest method yields a significant improvement over the multitrajectory Ehrenfest method and can be converged to the exact results in small systems with moderate computational efforts. We further observe that for extended systems, this convergence is slower with respect to the number of configurations. Our benchmark study demonstrates that DMRG-LBO is a useful tool for assessing the quality of the quantum-chemistry methods.
Collapse
Affiliation(s)
- M Ten Brink
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - S Gräber
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M Hopjan
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - D Jansen
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - J Stolpp
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - F Heidrich-Meisner
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - P E Blöchl
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
17
|
Bian X, Wu Y, Teh HH, Subotnik JE. Incorporating Berry Force Effects into the Fewest Switches Surface-Hopping Algorithm: Intersystem Crossing and the Case of Electronic Degeneracy. J Chem Theory Comput 2022; 18:2075-2090. [PMID: 35263116 DOI: 10.1021/acs.jctc.1c01103] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a preliminary surface-hopping approach for modeling intersystem crossing (ISC) dynamics between four electronic states: one singlet and one (triply degenerate) triplet. In order to incorporate all Berry force effects, the algorithm requires that, when moving along an adiabatic surface associated with the triplet manifold, one must also keep track of a quasi-diabatic index (akin to a "ms" quantum number) for each trajectory. For a simple model problem, we find that a great deal of new physics can be captured by our algorithm, setting the stage for larger, more realistic (or perhaps even ab initio) simulations in the future.
Collapse
Affiliation(s)
- Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Yanze Wu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Hung-Hsuan Teh
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
18
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
19
|
Gima K, Inerbaev TM, Kilin DS. Excited state dynamics in a sodium and iodine co-doped lead telluride nanowire. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1874557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Kevin Gima
- Department of Physics, North Dakota State University, Fargo, ND, USA
| | - Talgat M. Inerbaev
- L. N. Gumilyov Eurasian National University, Nur-Sultan, Kazakhstan
- Sobolev Institute of Geology and Mineralogy SB RAS, Novosibirsk, Russia
- National University of Science and Technology MISIS, Moscow, Russian Federation
| | - D. S. Kilin
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND, USA
| |
Collapse
|
20
|
Song H, Nam Y, Keefer D, Garavelli M, Mukamel S, Tretiak S. Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene. J Phys Chem Lett 2021; 12:5716-5722. [PMID: 34128675 DOI: 10.1021/acs.jpclett.1c01397] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S2/S3(D2d) → S1(D8h) → S0(D4h) → S0(D2d) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. The obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| | - Yeonsig Nam
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Daniel Keefer
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Marco Garavelli
- Department of Industrial Chemistry, "T. Montanari", Università degli Studi di Bologna, Viale del Risorgimento, 4, Bologna 40136, Italy
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine 92697, California, United States
- Department of Physics and Astronomy, University of California, Irvine 92697, California, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| |
Collapse
|
21
|
Song H, Freixas VM, Fernandez-Alberti S, White AJ, Zhang Y, Mukamel S, Govind N, Tretiak S. An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2021; 17:3629-3643. [PMID: 34014085 DOI: 10.1021/acs.jctc.1c00131] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD, Bernal, Argentina
| | | | - Alexander J White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Shaul Mukamel
- Departments of Chemistry, Physics, and Astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
22
|
Brian D, Liu Z, Dunietz BD, Geva E, Sun X. Three-state harmonic models for photoinduced charge transfer. J Chem Phys 2021; 154:174105. [PMID: 34241055 DOI: 10.1063/5.0050289] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensed-phase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi's golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid-porphyrin-C60 (CPC60) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor-acceptor, donor-ground, and acceptor-ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC60/THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC60/THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.
Collapse
Affiliation(s)
- Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Zengkui Liu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
23
|
Freixas VM, White AJ, Nelson T, Song H, Makhov DV, Shalashilin D, Fernandez-Alberti S, Tretiak S. Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence. J Phys Chem Lett 2021; 12:2970-2982. [PMID: 33730495 DOI: 10.1021/acs.jpclett.1c00266] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum-classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes, Roque Saénz Peña 352, B1876BXD Bernal, Argentina
| | - Alexander J White
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, U.K
- School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K
| | | | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
24
|
Erdmann E, Aguirre NF, Indrajith S, Chiarinelli J, Domaracka A, Rousseau P, Huber BA, Bolognesi P, Richter R, Avaldi L, Díaz-Tendero S, Alcamí M, Łabuda M. A general approach to study molecular fragmentation and energy redistribution after an ionizing event. Phys Chem Chem Phys 2021; 23:1859-1867. [PMID: 33439170 DOI: 10.1039/d0cp04890a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We propose to combine quantum chemical calculations, statistical mechanical methods, and photoionization and particle collision experiments to unravel the redistribution of internal energy of the furan cation and its dissociation pathways. This approach successfully reproduces the relative intensity of the different fragments as a function of the internal energy of the system in photoelectron-photoion coincidence experiments and the different mass spectra obtained when ions ranging from Ar+ to Xe25+ or electrons are used in collision experiments. It provides deep insights into the redistribution of the internal energy in the ionized molecule and its influence on the dissociation pathways and resulting charged fragments. The present pilot study demonstrates the efficiency of a statistical exchange of excitation energy among various degrees of freedom of the molecule and proves that the proposed approach is mature to be extended to more complex systems.
Collapse
Affiliation(s)
- Ewa Erdmann
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland.
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
25
|
Rodríguez-Hernández B, Nelson T, Oldani N, Martínez-Mesa A, Uranga-Piña L, Segawa Y, Tretiak S, Itami K, Fernandez-Alberti S. Exciton Spatial Dynamics and Self-Trapping in Carbon Nanocages. J Phys Chem Lett 2021; 12:224-231. [PMID: 33326240 DOI: 10.1021/acs.jpclett.0c03364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Three-dimensional cage-shaped molecules formed from chainlike structures hold potential as unique optoelectronic materials and host compounds. Their optical, structural, and dynamical features are tunable by changes in shape and size. We perform a comparison of these properties for three sizes of strained conjugated [n.n.n]carbon nanocages composed of three paraphenylene chains (bridges) of length n = 4, 5, or 6. The exciton intramolecular redistribution occurring during nonradiative relaxation has been explored using nonadiabatic excited-state molecular dynamics. Our results provide atomistic insight into the conformational features associated with the observed red- and blue-shift trends in the absorption and fluorescence spectra, respectively, with increasing nanocage size. Their internal conversion processes involve intramolecular energy transfer that leads to exciton self-trapping on a few phenylene units at the center of a single bridge. The dependence of these dynamical features on the size of the nanocage can be used to tune their host-guest chemical properties and their use for organic electronics and catenane-like applications.
Collapse
Affiliation(s)
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nicolas Oldani
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Aliezer Martínez-Mesa
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
- DynAMoS (Dynamical processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana, San Lázaro y L, La Habana 10400, Cuba
| | - Llinersy Uranga-Piña
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
- DynAMoS (Dynamical processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana, San Lázaro y L, La Habana 10400, Cuba
| | - Yasutomo Segawa
- Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- JST, ERATO, Itami Molecular Nanocarbon Project, Nagoya University, Nagoya 464-8602, Japan
- Institute for Molecular Science, Myodaiji, Okazaki 444-8787, Japan
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Myodaiji, Okazaki 444-8787, Japan
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Kenichiro Itami
- Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
- JST, ERATO, Itami Molecular Nanocarbon Project, Nagoya University, Nagoya 464-8602, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Nagoya 464-8602, Japan
| | | |
Collapse
|
26
|
Mukazhanova A, Malone W, Negrin-Yuvero H, Fernandez-Alberti S, Tretiak S, Sharifzadeh S. Photoexcitation dynamics in perylene diimide dimers. J Chem Phys 2020; 153:244117. [PMID: 33380092 DOI: 10.1063/5.0031485] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We utilize first-principles theory to investigate photo-induced excited-state dynamics of functionalized perylene diimide. This class of materials is highly suitable for solar energy conversion because of the strong optical absorbance, efficient energy transfer, and chemical tunability. We couple time-dependent density functional theory to a recently developed time-resolved non-adiabatic dynamics approach based on a semi-empirical description. By studying the monomer and dimer, we focus on the role stacking plays on the time-scales associated with excited-state non-radiative relaxation from a high excitonic state to the lowest energy exciton. We predict that the time-scale for energy conversion in the dimer is significantly faster than that in the monomer when equivalent excited states are accounted for. Additionally, for the dimer, the decay from the second to the nearly degenerate lowest energy excited-state involves two time-scales: a rapid decay on the order of ∼10 fs followed by a slower decay of ∼100 fs. Analysis of the spatial localization of the electronic transition density during the internal conversion process points out the existence of localized states on individual monomers, indicating that the strength of thermal fluctuations exceeds electronic couplings between the states such that the exciton hops between localized states throughout the simulation.
Collapse
Affiliation(s)
- Aliya Mukazhanova
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, USA
| | - Walter Malone
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Hassiel Negrin-Yuvero
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Sahar Sharifzadeh
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, USA
| |
Collapse
|
27
|
Freixas VM, Nelson T, Ondarse-Alvarez D, Nijjar P, Mikhailovsky A, Zhou C, Fernandez-Alberti S, Bazan GC, Tretiak S. Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful? J Chem Phys 2020; 153:244114. [PMID: 33380074 DOI: 10.1063/5.0028126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born-Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor-acceptor-donor (D-A-D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are "sufficiently good" to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Parmeet Nijjar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| | - Alexander Mikhailovsky
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Cheng Zhou
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
| | | | - Guillermo C Bazan
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
28
|
Song H, Fischer SA, Zhang Y, Cramer CJ, Mukamel S, Govind N, Tretiak S. First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2020; 16:6418-6427. [DOI: 10.1021/acs.jctc.0c00295] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Sean A. Fischer
- Chemistry Division, U.S. Naval Research Laboratory, Washington, District of Columbia 20375, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Christopher J. Cramer
- Department of Chemistry, Supercomputing Institute and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Shaul Mukamel
- Departments of Chemistry, and physics and astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
29
|
Malone W, Nebgen B, White A, Zhang Y, Song H, Bjorgaard JA, Sifain AE, Rodriguez-Hernandez B, Freixas VM, Fernandez-Alberti S, Roitberg AE, Nelson TR, Tretiak S. NEXMD Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5771-5783. [DOI: 10.1021/acs.jctc.0c00248] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Walter Malone
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Benjamin Nebgen
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Alexander White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Josiah A. Bjorgaard
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E. Sifain
- U.S. Army Research Laboratory, Aberdeen Proving Ground, Aberdeen, Maryland 21005, United States
| | | | - Victor M. Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Adrian E. Roitberg
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Tammie R. Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
30
|
Rodríguez-Hernández B, Oldani N, Martínez-Mesa A, Uranga-Piña L, Tretiak S, Fernandez-Alberti S. Photoexcited energy relaxation and vibronic couplings in π-conjugated carbon nanorings. Phys Chem Chem Phys 2020; 22:15321-15332. [PMID: 32628225 DOI: 10.1039/d0cp01452d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Conjugated carbon nanorings exhibit unique photophysical properties that, combined with their tunable sizes and conformations, make them suitable for a variety of practical applications. These properties are intimately associated to their strained, bent and sterically hindered cyclic structures. Herein we perform a comparative analysis of the photoinduced dynamics in carbon nanorings composed of nine phenyl units([9]CPP) and nine naphthyl units ([9]CN) respectively. The sterically demanding naphthyl units lead to large dihedral angles between neighboring units. Nevertheless, the ultrafast electronic and vibrational energy relaxation and redistribution is found to be similar for both systems. We observe that vibronic couplings, introduced by nonadiabatic energy transfer between electronic excited states, ensure the intramolecular vibrational energy redistribution through specific vibrational modes. The comparative impact of the internal conversion process on the exciton spatial localization and intra-ring migration indicates that naphthyl units in [9]CN achieve more efficient but less dynamical self-trapping compared to that of phenyl units in [9]CPP. That is, during the photoinduced process, the exciton in [9]CN is more static and localized than the exciton in [9]CPP. The internal conversion processes take place through a specific set of middle- to high-frequency normal modes, which directly influence the spatial exciton redistribution during the internal conversion, self-trapping and intra-ring migration.
Collapse
Affiliation(s)
- B Rodríguez-Hernández
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| | | | | | | | | | | |
Collapse
|
31
|
Straub S, Stubbe J, Lindner J, Sarkar B, Vöhringer P. Vibrational Relaxation Dynamics of an Azido–Cobalt(II) Complex from Femtosecond UV-Pump/MIR-Probe Spectroscopy and Model Simulations with Ab Initio Anharmonic Couplings. Inorg Chem 2020; 59:14629-14642. [DOI: 10.1021/acs.inorgchem.0c00553] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Steffen Straub
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| | - Jessica Stubbe
- Institut für Chemie und Biochemie, Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/34, 14195 Berlin, Germany
| | - Jörg Lindner
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| | - Biprajit Sarkar
- Institut für Chemie und Biochemie, Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/34, 14195 Berlin, Germany
- Lehrstuhl für Anorganische Koordinationschemie, Institut für Anorganische Chemie Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter Vöhringer
- Lehrstuhl für Molekulare Physikalische Chemie, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms-Universität, Wegelerstraße 12, 53115 Bonn, Germany
| |
Collapse
|
32
|
Freixas VM, Tretiak S, Makhov DV, Shalashilin DV, Fernandez-Alberti S. Vibronic Quantum Beating between Electronic Excited States in a Heterodimer. J Phys Chem B 2020; 124:3992-4001. [DOI: 10.1021/acs.jpcb.0c01685] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- V. M. Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - S. Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - D. V. Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, U.K
- School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K
| | | | - S. Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| |
Collapse
|
33
|
Zhang Y, Li L, Tretiak S, Nelson T. Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems. J Chem Theory Comput 2020; 16:2053-2064. [DOI: 10.1021/acs.jctc.9b00928] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Linqiu Li
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
34
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
35
|
Alfonso-Hernandez L, Athanasopoulos S, Tretiak S, Miguel B, Bastida A, Fernandez-Alberti S. Vibrational energy redistribution during donor–acceptor electronic energy transfer: criteria to identify subsets of active normal modes. Phys Chem Chem Phys 2020; 22:18454-18466. [DOI: 10.1039/d0cp03102j] [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/21/2022]
Abstract
Photoinduced electronic energy transfer in conjugated donor–acceptor systems is naturally accompanied by intramolecular vibrational energy redistributions accepting an excess of electronic energy.
Collapse
Affiliation(s)
- L. Alfonso-Hernandez
- CONICET – Centro Integral de Medicina Nuclear y Radioterapia
- Centro Atómico Bariloche
- Río Negro
- Argentina
| | - S. Athanasopoulos
- Departamento de Física, Universidad Carlos III de Madrid
- 28911 Leganés
- Spain
| | - S. Tretiak
- Theoretical Division
- Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT)
- Los Alamos National Laboratory
- Los Alamos
- USA
| | - B. Miguel
- Departamento de Ingeniería Química y Ambiental
- Universidad Politécnica de Cartagena
- 30203 Cartagena
- Spain
| | - A. Bastida
- Departamento de Química Física
- Universidad de Murcia
- 30100 Murcia
- Spain
| | - S. Fernandez-Alberti
- Departamento de Ciencia y Tecnologia
- Universidad Nacional de Quilmes/CONICET
- B1876BXD Bernal
- Argentina
| |
Collapse
|
36
|
Carof A, Giannini S, Blumberger J. How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm. Phys Chem Chem Phys 2019; 21:26368-26386. [PMID: 31793569 DOI: 10.1039/c9cp04770k] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.
Collapse
Affiliation(s)
- Antoine Carof
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Samuele Giannini
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK. and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
37
|
Zhang Y, Nelson T, Tretiak S. Non-adiabatic molecular dynamics of molecules in the presence of strong light-matter interactions. J Chem Phys 2019; 151:154109. [PMID: 31640366 DOI: 10.1063/1.5116550] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
38
|
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
|
39
|
Miao G, Subotnik J. Revisiting the Recoherence Problem in the Fewest Switches Surface Hopping Algorithm. J Phys Chem A 2019; 123:5428-5435. [DOI: 10.1021/acs.jpca.9b03188] [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)
- Gaohan Miao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
40
|
Sifain AE, Wang L, Tretiak S, Prezhdo OV. Numerical tests of coherence-corrected surface hopping methods using a donor-bridge-acceptor model system. J Chem Phys 2019; 150:194104. [DOI: 10.1063/1.5092999] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew E. Sifain
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089-0485, USA
- Theoretical Division, Center for Nonlinear Studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Angeles, New Mexico 87545, USA
| | - Linjun Wang
- Center for Chemistry of Novel and High-Performance Materials and Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Angeles, New Mexico 87545, USA
- Skolkovo Institute of Science and Technology, Moscow 143026, Russia
| | - Oleg V. Prezhdo
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089-0485, USA
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| |
Collapse
|
41
|
Miao G, Bellonzi N, Subotnik J. An extension of the fewest switches surface hopping algorithm to complex Hamiltonians and photophysics in magnetic fields: Berry curvature and “magnetic” forces. J Chem Phys 2019; 150:124101. [DOI: 10.1063/1.5088770] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gaohan Miao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
42
|
Parker SM, Roy S, Furche F. Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation. Phys Chem Chem Phys 2019; 21:18999-19010. [DOI: 10.1039/c9cp03127h] [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/21/2022]
Abstract
We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.
Collapse
Affiliation(s)
- Shane M. Parker
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Saswata Roy
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Filipp Furche
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| |
Collapse
|
43
|
Ondarse-Alvarez D, Nelson T, Lupton JM, Tretiak S, Fernandez-Alberti S. Let Digons be Bygones: The Fate of Excitons in Curved π-Systems. J Phys Chem Lett 2018; 9:7123-7129. [PMID: 30508376 DOI: 10.1021/acs.jpclett.8b03160] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We explore the diverse origins of unpolarized absorption and emission of molecular polygons consisting of π-conjugated oligomer chains held in a bent geometry by strain controlled at the vertex units. For this purpose, we make use of atomistic nonadiabatic excited-state molecular dynamics simulations of a bichromophore molecular polygon (digon) with bent chromophore chains. Both structural and photoexcited dynamics were found to affect polarization features. Bending strain induces exciton localization on individual chromophore units of the conjugated chains. The latter display different transition dipole moment orientations, a feature not present in the linear oligomer counterparts. In addition, bending makes exciton localization very sensitive to molecular distortions induced by thermal fluctuations. The excited-state dynamics reveals an ultrafast intramolecular energy redistribution that spreads the exciton equally among spatially separated chromophore fragments within the molecular system. As a result, digons become virtually unpolarized absorbers and emitters, in agreement with recent experimental studies on the single-molecule level.
Collapse
Affiliation(s)
| | - Tammie Nelson
- Theoretical Division, Physics and Chemistry of Materials (T-1) , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - John M Lupton
- Institut für Angewandte und Experimentelle Physik , Universität Regensburg , Universitätsstrasse 31 , 93053 Regensburg , Germany
| | - Sergei Tretiak
- Theoretical Division, Physics and Chemistry of Materials (T-1) , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | |
Collapse
|
44
|
Mandal A, Hunt KLC. Nonadiabatic transition probabilities in a time-dependent Gaussian pulse or plateau pulse: Toward experimental tests of the differences from Dirac’s transition probabilities. J Chem Phys 2018; 149:204110. [DOI: 10.1063/1.5054313] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anirban Mandal
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Katharine L. C. Hunt
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
45
|
Cozza C, Bonomi M, Pietropaolo A. A Versatile Computational Strategy To Characterize the Free-Energy Landscape of Excited States in Oligofluorenes. J Chem Theory Comput 2018; 14:5441-5445. [PMID: 30365313 DOI: 10.1021/acs.jctc.8b00949] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The π-rich rings of conjugated polymers and molecular rotors shape their typical properties, allowing a variety of chemical and photoresponsive phenomena. Herein, we present a computational method in the framework of classical simulations to estimate the free-energy gap between ground and excited states of oligofluorenes, which is a class of molecular rotors widely used in optoelectronic devices, because of the inner torsional rotation triggered by light irradiation. We devised multiple sets of free-energy simulations in combination with free-energy perturbation theory to predict the free-energy gap between the ground state and the first excited state. The computed excitation energies show good agreement with experiments. The approach presented herein allows one to achieve comprehensive sampling of the conformational landscape and accurate estimates of the excited state free-energy landscapes at the same time.
Collapse
Affiliation(s)
- Concetta Cozza
- Dipartimento di Scienze della Salute , Università di Catanzaro , Viale Europa , 88100 Catanzaro , Italy
| | - Massimiliano Bonomi
- Department of Chemistry , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - Adriana Pietropaolo
- Dipartimento di Scienze della Salute , Università di Catanzaro , Viale Europa , 88100 Catanzaro , Italy
| |
Collapse
|
46
|
Sifain AE, Gifford BJ, Gao DW, Lystrom L, Nelson TR, Tretiak S. NEXMD Modeling of Photoisomerization Dynamics of 4-Styrylquinoline. J Phys Chem A 2018; 122:9403-9411. [DOI: 10.1021/acs.jpca.8b09103] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrew E. Sifain
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089-0485, United States
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Brendan J. Gifford
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - David W. Gao
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Los Alamos High School, Los Alamos, New Mexico 87544, United States
| | - Levi Lystrom
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, North Dakota 58108, United States
| | - Tammie R. Nelson
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
47
|
Wang R, Chernyak VY. Dynamical consequences of time-reversal symmetry for systems with odd number of electrons: Conical intersections, semiclassical dynamics, and topology. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.09.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
48
|
Lystrom L, Zhang Y, Tretiak S, Nelson T. Site-Specific Photodecomposition in Conjugated Energetic Materials. J Phys Chem A 2018; 122:6055-6061. [DOI: 10.1021/acs.jpca.8b04381] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Levi Lystrom
- Theoretical Division, Physics and Chemistry of Materials (T-1), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, North Dakota 58108-6050, United States
| | - Yu Zhang
- Theoretical Division, Physics and Chemistry of Materials (T-1), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Theoretical Division, Physics and Chemistry of Materials (T-1), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Theoretical Division, Physics and Chemistry of Materials (T-1), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
49
|
Nelson TR, Ondarse-Alvarez D, Oldani N, Rodriguez-Hernandez B, Alfonso-Hernandez L, Galindo JF, Kleiman VD, Fernandez-Alberti S, Roitberg AE, Tretiak S. Coherent exciton-vibrational dynamics and energy transfer in conjugated organics. Nat Commun 2018; 9:2316. [PMID: 29899334 PMCID: PMC5998141 DOI: 10.1038/s41467-018-04694-8] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2018] [Accepted: 04/27/2018] [Indexed: 11/09/2022] Open
Abstract
Coherence, signifying concurrent electron-vibrational dynamics in complex natural and man-made systems, is currently a subject of intense study. Understanding this phenomenon is important when designing carrier transport in optoelectronic materials. Here, excited state dynamics simulations reveal a ubiquitous pattern in the evolution of photoexcitations for a broad range of molecular systems. Symmetries of the wavefunctions define a specific form of the non-adiabatic coupling that drives quantum transitions between excited states, leading to a collective asymmetric vibrational excitation coupled to the electronic system. This promotes periodic oscillatory evolution of the wavefunctions, preserving specific phase and amplitude relations across the ensemble of trajectories. The simple model proposed here explains the appearance of coherent exciton-vibrational dynamics due to non-adiabatic transitions, which is universal across multiple molecular systems. The observed relationships between electronic wavefunctions and the resulting functionalities allows us to understand, and potentially manipulate, excited state dynamics and energy transfer in molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division, Center for Nonlinear studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, NM, 81545, USA
| | | | - Nicolas Oldani
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD, Bernal, Argentina
| | | | | | - Johan F Galindo
- Department of Chemistry, Universidad Nacional de Colombia, Av. Cra 30 # 45-03, Bogotá, 111321, Colombia
| | - Valeria D Kleiman
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA
| | | | - Adrian E Roitberg
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear studies and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, NM, 81545, USA.
| |
Collapse
|
50
|
Carof A, Giannini S, Blumberger J. Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping. J Chem Phys 2018; 147:214113. [PMID: 29221382 DOI: 10.1063/1.5003820] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.
Collapse
Affiliation(s)
- Antoine Carof
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Samuele Giannini
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|