1
|
Salazar EX, Menger MFSJ, Faraji S. Ultrafast Photoinduced Dynamics in 1,3-Cyclohexadiene: A Comparison of Trajectory Surface Hopping Schemes†. J Chem Theory Comput 2024. [PMID: 38949625 DOI: 10.1021/acs.jctc.4c00012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/02/2024]
Abstract
Photoinduced nonadiabatic processes play a crucial role in a wide range of disciplines, from fundamental steps in biology to modern applications in advanced materials science. A theoretical understanding of these processes is highly desirable, and trajectory surface hopping (TSH) has proven to be a well-suited framework for a wide range of systems. In this work, we present a comprehensive comparison between two TSH algorithms, the conventional Tully's fewest switches surface hopping (FSSH) scheme and the Landau-Zener surface hopping (LZSH), to study the photoinduced ring-opening of 1,3-cyclohexadiene (CHD) to 1,3,5-hexatriene at the spin-flip time-dependent density functional theory (SF-TDDFT) level of theory. Additionally, we compare our results with a literature study at the extended multistate complete active space second-order perturbation theory method (XMS-CASPT2) level of theory. Our results show that the average population and lifetimes estimated with LZSH using SF-TDDFT are closer to the literature (using multireference methods) than those estimated with FSSH using SF-TDDFT. The latter speaks in favor of applying LZSH in combination with the SF-TDDFT method to study larger and more complex systems such as molecular photoswitches where the CHD molecule acts as a backbone. In addition, we present an implementation of Tully's FSSH algorithm as an extension to the PySurf software package.
Collapse
Affiliation(s)
- Edison X Salazar
- Instituut-Lorentz, Universiteit Leiden, 2300 RA Leiden, The Netherlands
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Maximilian F S J Menger
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, 69120 Heidelberg, Germany
| | - Shirin Faraji
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
- Institute of Theoretical and Computational Chemistry, Faculty of Mathematics and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
2
|
Zhang Q, Shao X, Li W, Mi W, Pavanello M, Akimov AV. Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:385901. [PMID: 38866023 DOI: 10.1088/1361-648x/ad577d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 06/12/2024] [Indexed: 06/14/2024]
Abstract
In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.
Collapse
Affiliation(s)
- Qingxin Zhang
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| | - Xuecheng Shao
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China
| | - Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, People's Republic of China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People's Republic of China
| | - Michele Pavanello
- Department of Physics, Rutgers University, The State University of New Jersey, Newark, NJ 07102, United States of America
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260, United States of America
| |
Collapse
|
3
|
Shakiba M, Akimov AV. Machine-Learned Kohn-Sham Hamiltonian Mapping for Nonadiabatic Molecular Dynamics. J Chem Theory Comput 2024; 20:2992-3007. [PMID: 38581699 DOI: 10.1021/acs.jctc.4c00008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2024]
Abstract
In this work, we report a simple, efficient, and scalable machine-learning (ML) approach for mapping non-self-consistent Kohn-Sham Hamiltonians constructed with one kind of density functional to the nearly self-consistent Hamiltonians constructed with another kind of density functional. This approach is designed as a fast surrogate Hamiltonian calculator for use in long nonadiabatic dynamics simulations of large atomistic systems. In this approach, the input and output features are Hamiltonian matrices computed from different levels of theory. We demonstrate that the developed ML-based Hamiltonian mapping method (1) speeds up the calculations by several orders of magnitude, (2) is conceptually simpler than alternative ML approaches, (3) is applicable to different systems and sizes and can be used for mapping Hamiltonians constructed with arbitrary density functionals, (4) requires a modest training data, learns fast, and generates molecular orbitals and their energies with the accuracy nearly matching that of conventional calculations, and (5) when applied to nonadiabatic dynamics simulation of excitation energy relaxation in large systems yields the corresponding time scales within the margin of error of the conventional calculations. Using this approach, we explore the excitation energy relaxation in C60 fullerene and Si75H64 quantum dot structures and derive qualitative and quantitative insights into dynamics in these systems.
Collapse
Affiliation(s)
- Mohammad Shakiba
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
4
|
Lee IS, Filatov M, Min SK. Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states. J Chem Phys 2024; 160:154103. [PMID: 38624116 DOI: 10.1063/5.0202095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 03/31/2024] [Indexed: 04/17/2024] Open
Abstract
A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.
Collapse
Affiliation(s)
- In Seong Lee
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Michael Filatov
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Center for Multidimensional Carbon Materials (CMCM), Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, Republic of Korea
| |
Collapse
|
5
|
Guo X, Li G, Shi Z, Wang L. Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching. J Phys Chem Lett 2024:3345-3353. [PMID: 38498301 DOI: 10.1021/acs.jpclett.4c00437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.
Collapse
Affiliation(s)
- Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
6
|
Akimov AV. Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:11673-11683. [PMID: 38109379 DOI: 10.1021/acs.jpclett.3c03029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
7
|
Stanton R, Trivedi DJ. Charge Carrier Dynamics at the Interface of 2D Metal-Organic Frameworks and Hybrid Perovskites for Solar Energy Harvesting. NANO LETTERS 2023; 23:11932-11939. [PMID: 38100376 DOI: 10.1021/acs.nanolett.3c04054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
Interfacing perovskites with two-dimensional materials such as metal-organic frameworks (MOFs) for improved stability and electron or hole extraction has emerged as a promising path forward for the generation of highly efficient and stable solar cells. In this work, we examine the structural properties and excitation dynamics of two MOF-perovskite systems: UMCM309-a@MAPbI3 and ZrL3@MAPbI3. We find that precise band alignment and electronegativity of the MOF-linkers are necessary to facilitate the capture of excited charge carriers. Furthermore, we demonstrate that intraband relaxation of hot electrons to the MOF subsystem results in optically disallowed transitions across the band gap, suppressing radiative recombination. Furthermore, we elucidate the key mechanisms associated with improved structural stability afforded to the perovskites by the two-dimensional MOFs, highlighting the necessity of broad surface coverage and strong MOF-perovskite interaction.
Collapse
Affiliation(s)
- Robert Stanton
- Department of Physics, Clarkson University, Potsdam, New York 13699, United States
| | - Dhara J Trivedi
- Department of Physics, Clarkson University, Potsdam, New York 13699, United States
| |
Collapse
|
8
|
Huang DM, Green AT, Martens CC. A first principles derivation of energy-conserving momentum jumps in surface hopping simulations. J Chem Phys 2023; 159:214108. [PMID: 38047505 DOI: 10.1063/5.0178534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Accepted: 11/06/2023] [Indexed: 12/05/2023] Open
Abstract
The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]-along with its many later variations-forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH.
Collapse
Affiliation(s)
| | - Austin T Green
- University of California, Irvine, California 92697-2025, USA
| | - Craig C Martens
- University of California, Irvine, California 92697-2025, USA
| |
Collapse
|
9
|
Liu XY, Chen WK, Fang WH, Cui G. Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications. J Chem Theory Comput 2023. [PMID: 37984502 DOI: 10.1021/acs.jctc.3c00960] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Hefei National Laboratory, Hefei 230088, P. R. China
| |
Collapse
|
10
|
Freixas VM, Oldani N, Tretiak S, Fernandez-Alberti S. Twisting Aromaticity and Photoinduced Dynamics in Hexapole Helicenes. J Phys Chem Lett 2023; 14:10145-10150. [PMID: 37924328 DOI: 10.1021/acs.jpclett.3c02628] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2023]
Abstract
Curved aromatic molecules are attractive electronic materials, where an additional internal strain uniquely modifies their structure, aromaticity, dynamics, and optical properties. Helicenes are examples of such twisted conjugated systems. Herein, we analyze the photoinduced dynamics in different stereoisomers of a hexapole helicene by using nonadiabatic excited-state molecular dynamics simulations. We explore how changes in symmetry and structural distortion modulate the intramolecular energy redistribution. We find that distinct helical assembly leads to different rigid distorted structures that in turn impact the nonradiative energy relaxation and ultimately formation of the self-trapped exciton. Subsequently, the value of the twisting angles relative to the central triphenylene core structure controls the global molecular aromaticity and electronic localization during the internal conversion process. Our work sheds light on how the future synthesis of novel curved aromatic compounds can be directed to attain specific desired electronic properties through the modulation of their twisted aromaticity.
Collapse
Affiliation(s)
- Victor M Freixas
- Department of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | - Nicolas Oldani
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | | |
Collapse
|
11
|
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
|
12
|
Shao C, Shi Z, Xu J, Wang L. Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics. J Phys Chem Lett 2023; 14:7680-7689. [PMID: 37606199 DOI: 10.1021/acs.jpclett.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.
Collapse
Affiliation(s)
- Cancan Shao
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Jiabo Xu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
13
|
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
|
14
|
Gumber S, Prezhdo OV. Zeno and Anti-Zeno Effects in Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2023; 14:7274-7282. [PMID: 37556319 PMCID: PMC10440816 DOI: 10.1021/acs.jpclett.3c01831] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Accepted: 08/01/2023] [Indexed: 08/11/2023]
Abstract
Decoherence plays an important role in nonadiabatic (NA) molecular dynamics (MD) simulations because it provides a physical mechanism for trajectory hopping and can alter transition rates by orders of magnitude. Generally, decoherence effects slow quantum transitions, as exemplified by the quantum Zeno effect: in the limit of infinitely fast decoherence, the transitions stop. If the measurements are not sufficiently frequent, an opposite quantum anti-Zeno effect occurs, in which the transitions are accelerated with faster decoherence. Using two common NA-MD approaches, fewest switches surface hopping and decoherence-induced surface hopping, combined with analytic examination, we demonstrate that including decoherence into NA-MD slows down NA transitions; however, many realistic systems operate in the anti-Zeno regime. Therefore, it is important that NA-MD methods describe both Zeno and anti-Zeno effects. Numerical simulations of charge trapping and relaxation in graphitic carbon nitride suggest that time-dependent NA Hamiltonians encountered in realistic systems produce robust results with respect to errors in the decoherence time, a favorable feature for NA-MD simulations.
Collapse
Affiliation(s)
- Shriya Gumber
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department
of Physics and Astronomy, University of
Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
15
|
Negrin-Yuvero H, Freixas VM, Ondarse-Alvarez D, Alfonso-Hernandez L, Rojas-Lorenzo G, Bastida A, Tretiak S, Fernandez-Alberti S. Vibrational Funnels for Energy Transfer in Organic Chromophores. J Phys Chem Lett 2023; 14:4673-4681. [PMID: 37167537 DOI: 10.1021/acs.jpclett.3c00748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Photoinduced intramolecular energy transfers in multichromophoric molecules involve nonadiabatic vibronic channels that act as energy transfer funnels. They commonly take place through specific directions of motion dictated by the nonadiabatic coupling vectors. Vibrational funnels may support persistent coherences between electronic states and sometimes delineate the presence of minor alternative energy transfer pathways. The ultimate confirmation of their role on the interchromophoric energy transfer can be achieved by performing nonadiabatic excited-state molecular dynamics simulations by selectively freezing the nuclear motions in question. Our results point out this strategy as a useful tool to identify and evaluate the impact of these vibrational funnels on the energy transfer processes and guide the in silico design of materials with tunable properties and enhanced functionalities. Our work encourages applications of this methodology to different chemical and biochemical processes such as reactive scattering and protein conformational changes, to name a few.
Collapse
Affiliation(s)
- Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Victor Manuel Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Dianelys Ondarse-Alvarez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Laura Alfonso-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - German Rojas-Lorenzo
- Departamento de Física Atómica y Molecular, Instituto Superior de Tecnologías y Ciencias Aplicadas, Universidad de La Habana, La Habana 10400, Cuba
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos 87545, New Mexico, USA
| | | |
Collapse
|
16
|
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
|
17
|
Weight BM, Sifain AE, Gifford BJ, Htoon H, Tretiak S. On-the-Fly Nonadiabatic Dynamics Simulations of Single-Walled Carbon Nanotubes with Covalent Defects. ACS NANO 2023; 17:6208-6219. [PMID: 36972076 DOI: 10.1021/acsnano.2c08579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Single-walled carbon nanotubes (SWCNTs) with covalent surface defects have been explored recently due to their promise for use in single-photon telecommunication emission and in spintronic applications. The all-atom dynamic evolution of electrostatically bound excitons (the primary electronic excitations) in these systems has only been loosely explored from a theoretical perspective due to the size limitations of these large systems (>500 atoms). In this work, we present computational modeling of nonradiative relaxation in a variety of SWCNT chiralities with single-defect functionalizations. Our excited-state dynamics modeling uses a trajectory surface hopping algorithm accounting for excitonic effects with a configuration interaction approach. We find a strong chirality and defect-composition dependence on the population relaxation (varying over 50-500 fs) between the primary nanotube band gap excitation E11 and the defect-associated, single-photon-emitting E11* state. These simulations give direct insight into the relaxation between the band-edge states and the localized excitonic state, in competition with dynamic trapping/detrapping processes observed in experiment. Engineering fast population decay into the quasi-two-level subsystem with weak coupling to higher-energy states increases the effectiveness and controllability of these quantum light emitters.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, United States
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E Sifain
- Department of Chemistry, Princeton University, Princeton, New Jersey 08540 United States
| | - Brendan J Gifford
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Han Htoon
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
18
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, France
| |
Collapse
|
19
|
Jain A, Sindhu A. Pedagogical Overview of the Fewest Switches Surface Hopping Method. ACS OMEGA 2022; 7:45810-45824. [PMID: 36570264 PMCID: PMC9773185 DOI: 10.1021/acsomega.2c04843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Accepted: 11/18/2022] [Indexed: 06/17/2023]
Abstract
The fewest switches surface hopping method continues to grow in popularity to capture electronic nonadiabaticity and quantum nuclear effects due to its simplicity and accuracy. Knowing the basics of the method is essential for the correct implementation and interpretation of results. This review covers the fundamentals of the fewest switches surface hopping method with a detailed discussion of the nuances such as decoherence schemes and frustrated hops and the correct approach to calculating populations. The consequences of incorrect implementation are further discussed toward calculating kinetic and thermodynamic properties. Some tips for practitioners and a step-by-step algorithm for developers are provided. Finally, some of the finer technicalities of the fewest switches surface hopping method that are buried deep in the literature are pointed out to help graduate students better appreciate this method.
Collapse
|
20
|
Tang D, Jia L, Shen L, Fang WH. Fewest-Switches Surface Hopping with Long Short-Term Memory Networks. J Phys Chem Lett 2022; 13:10377-10387. [PMID: 36317657 DOI: 10.1021/acs.jpclett.2c02299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The mixed quantum-classical dynamical simulation is essential for studying nonadiabatic phenomena in photophysics and photochemistry. In recent years, many machine learning models have been developed to accelerate the time evolution of the nuclear subsystem. Herein, we implement long short-term memory (LSTM) networks as a propagator to accelerate the time evolution of the electronic subsystem during the fewest-switches surface hopping (FSSH) simulations. A small number of reference trajectories are generated using the original FSSH method, and then the LSTM networks can be built, accompanied by careful examination of typical LSTM-FSSH trajectories that employ the same initial condition and random numbers as the corresponding reference. The constructed network is applied to FSSH to further produce a trajectory ensemble to reveal the mechanism of nonadiabatic processes. Taking Tully's three models as test systems, we qualitatively reproduced the collective results. This work demonstrates that LSTM can be applied to the most popular surface hopping simulations.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Luyang Jia
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| |
Collapse
|
21
|
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
|
22
|
Alfonso Hernandez L, Freixas VM, Rodriguez-Hernandez B, Tretiak S, Fernandez-Alberti S, Oldani N. Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring. Phys Chem Chem Phys 2022; 24:24095-24104. [PMID: 36178044 DOI: 10.1039/d2cp03162k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure-dynamics-optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.
Collapse
Affiliation(s)
- Laura Alfonso Hernandez
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| | - Victor M Freixas
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| | | | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Nicolas Oldani
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| |
Collapse
|
23
|
Freixas VM, Tretiak S, Fernandez-Alberti S. Infinitene: Computational Insights from Nonadiabatic Excited State Dynamics. J Phys Chem Lett 2022; 13:8495-8501. [PMID: 36066077 DOI: 10.1021/acs.jpclett.2c02296] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Progress in organic synthesis opens exploration of a rich diversity of molecules with interesting new structural topologies. This is the case of a recently synthesized helically twisted figure-eight molecule coined infinitene. The molecule belongs to a numerous family of looped polyarenes, where the degree of π-conjugation is controlled by high strain energies and steric hindrances. A particular balance of these ingredients leads to unusual optoelectronic properties potentially suitable for a range of applications in nanoelectronics and photonics. Due to its recent discovery, the photophysical properties of infinitene remain unexplored. In this Letter, atomistic nonadiabatic excited state molecular dynamics modeling unveils unique features of intramolecular electronic and vibrational energy relaxation and redistribution that take place after molecular photoexcitation. Our results detail relationships between optical and electronic properties providing useful knowledge for future molecular designs related to infinitene.
Collapse
Affiliation(s)
- Victor Manuel Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - 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
|
24
|
Kroonblawd MP, Yoshimura A, Goldman N, Maiti A, Lewicki JP, Saab AP. Multiscale Strategy for Predicting Radiation Chemistry in Polymers. J Chem Theory Comput 2022; 18:5117-5124. [PMID: 35960960 DOI: 10.1021/acs.jctc.2c00582] [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
A primary mode for radiation damage in polymers arises from ballistic electrons that induce electronic excitations, yet subsequent chemical mechanisms are poorly understood. We develop a multiscale strategy to predict this chemistry starting from subatomic scattering calculations. Nonadiabatic molecular dynamics simulations sample initial bond-breaking events following the most likely excitations, which feed into semiempirical simulations that approach chemical equilibrium. Application to polyethylene reveals a mechanism explaining the low propensity to cross-link in crystalline samples.
Collapse
Affiliation(s)
- Matthew P Kroonblawd
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Anthony Yoshimura
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Nir Goldman
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States.,Department of Chemical Engineering, University of California, Davis, California 95616, United States
| | - Amitesh Maiti
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - James P Lewicki
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Andrew P Saab
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
25
|
Wu X, Wen S, Song H, Frauenheim T, Tretiak S, Yam C, Zhang Y. Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method. J Chem Phys 2022; 157:084114. [DOI: 10.1063/5.0100339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.
Collapse
Affiliation(s)
- Xiaoyan Wu
- Shenzhen JL Computational Science and Applied Research Institute, Shenzhen JL Computational Science and Applied Research Institute, China
| | | | - Huajing Song
- Los Alamos National Laboratory, United States of America
| | | | - Sergei Tretiak
- Theoretical Division, T-1, Los Alamos National Laboratory, United States of America
| | - ChiYung Yam
- Beijing Computational Science Research Center, Beijing Computational Science Research Center, China
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, United States of America
| |
Collapse
|
26
|
Shakiba M, Stippell E, Li W, Akimov AV. Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids. J Chem Theory Comput 2022; 18:5157-5180. [PMID: 35758936 DOI: 10.1021/acs.jctc.2c00297] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.
Collapse
Affiliation(s)
- Mohammad Shakiba
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Elizabeth Stippell
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, China
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
27
|
Ha JK, Min SK. Independent Trajectory Mixed Quantum-Classical Approaches Based on the Exact Factorization. J Chem Phys 2022; 156:174109. [DOI: 10.1063/5.0084493] [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
Mixed quantum-classical dynamics based on the exact factorization exploits the "derived" electron-nuclear correlation (ENC) term aiming for the description of quantum coherences. The ENC contains interactions between the phase of electronic states and nuclear quantum momenta which depend on the spatial shape of the nuclear density.The original surface hopping based on the exact factorization (SHXF) [\textit{J. Phys. Chem. Lett.} \textbf{2018}, \textit{9}, 1097] exploits frozen Gaussian functions to construct the nuclear density in the ENC term while the phase of electronic states is approximated as a fictitious nuclear momentum change.However, in reality, the width of nuclear wave packets varies in time depending on the shape of potential energy surfaces.In this work, we present a modified SHXF approach and a newly-developed Ehrenfest dynamics based on the exact factorization (EhXF) with time-dependent Gaussian functions and phases by enforcing total energy conservation.We perform numerical tests for various one-dimensional two-state model Hamiltonians.Overall, the time-dependent width of Gaussian functions and the energy conserving phase show a reliable decoherence compared to the original frozen Gaussian-based SHXF and the exact quantum mechanical calculation.Especially, the energy conserving phase is crucial for EhXF to reproduce the correct quantum dynamics.
Collapse
Affiliation(s)
- Jong-Kwon Ha
- Chemistry, Ulsan National Institute of Science and Technology, Korea, Republic of (South Korea)
| | - Seung Kyu Min
- Ulsan National Institute of Science and Technology, Korea, Republic of (South Korea)
| |
Collapse
|
28
|
Qiu J, Lu Y, Wang L. Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States. J Chem Theory Comput 2022; 18:2803-2815. [PMID: 35380833 DOI: 10.1021/acs.jctc.2c00130] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We present a multilayer subsystem surface hopping (MSSH) method to deal with nonadiabatic dynamics in large-scale systems. A small subsystem instead of the full system is adopted for surface hopping and is updated on-the-fly to achieve a reliable description of important adiabatic states and the wave function evolution. Additional subsystems for molecular dynamics and statistical description are introduced to further improve the simulation reliability. The global flux hopping probabilities with optimal state assignments are utilized to treat the complex surface crossings. As demonstrated in a series of one- and two-dimensional Holstein models with up to hundreds of thousands of states, MSSH shows weak parameter dependence in all investigated systems. Especially, the computational costs are reduced by 2-6 orders of magnitude compared to traditional surface hopping simulations in full systems, and size-independent results are achieved with a large time-step size of 2-5 fs. The new method is compatible with different decoherence correction strategies and achieves a much better balance between efficiency and reliability, thus promising for applications in general charge and exciton dynamics simulations.
Collapse
Affiliation(s)
- Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Yao Lu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
29
|
Xie BB, Jia PK, Wang KX, Chen WK, Liu XY, Cui G. Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems. J Phys Chem A 2022; 126:1789-1804. [PMID: 35266391 DOI: 10.1021/acs.jpca.1c10195] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Pei-Ke Jia
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Ke-Xin Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, Sichuan, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
30
|
Li B, Xu J, Li G, Shi Z, Wang L. A Mixed Deterministic-Stochastic Algorithm of the Branching Corrected Mean Field Method for Nonadiabatic Dynamics. J Chem Phys 2022; 156:114116. [DOI: 10.1063/5.0084013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new algorithm of the branching corrected mean field (BCMF) method for nonadiabatic dynamics [J. Xu and L. Wang, J. Phys. Chem. Lett. 11, 8283 (2020)], which combines the key advantages of the two existed algorithms, i.e., the deterministic BCMF algorithm based on weights of trajectory branches (BCMF-w) and the stochastic BCMF algorithm with random collapse of the electronic wavefunction (BCMF-s). The resulting mixed deterministic-stochastic BCMF algorithm (BCMF-ws) is benchmarked in a series of standard scattering problems with potential wells on the excited-state surfaces, which are common in realistic systems. In all investigated cases, BCMF-ws holds the same high accuracy while the computational time is reduced about two orders of magnitude compared to the original BCMF-w and BCMF-s algorithms, thus promising for nonadiabatic dynamics simulations of general systems.
Collapse
Affiliation(s)
| | | | | | | | - Linjun Wang
- Department of Chemistry, Zhejiang University, China
| |
Collapse
|
31
|
Negrin-Yuvero H, Mukazhanova A, Freixas VM, Tretiak S, Sharifzadeh S, Fernandez-Alberti S. Vibronic Photoexcitation Dynamics of Perylene Diimide: Computational Insights. J Phys Chem A 2022; 126:733-741. [PMID: 35084863 DOI: 10.1021/acs.jpca.1c09484] [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
Perylene diimide (PDI) represents a prototype material for organic optoelectronic devices because of its strong optical absorbance, chemical stability, efficient energy transfer, and optical and chemical tunability. Herein, we analyze in detail the vibronic relaxation of its photoexcitation using nonadiabatic excited-state molecular dynamics simulations. We find that after the absorption of a photon, which excites the electron to the second excited state, S2, induced vibronic dynamics features persistent modulations in the spatial localization of electronic and vibrational excitations. These energy exchanges are dictated by strong vibronic couplings that overcome structural disorders and thermal fluctuations. Specifically, the electronic wavefunction periodically swaps between localizations on the right and left sides of the molecule. Within 1 ps of such dynamics, a nonradiative transition to the lowest electronic state, S1, takes place, resulting in a complete delocalization of the wavefunction. The observed vibronic dynamics emerges following the electronic energy deposition in the direction that excites a combination of two dominant vibrational normal modes. This behavior is maintained even with a chemical substitution that breaks the symmetry of the molecule. We believe that our findings elucidate the nature of the complex dynamics of the optically excited states and, therefore, contribute to the development of tunable functionalities of PDIs and their derivatives.
Collapse
Affiliation(s)
- Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Aliya Mukazhanova
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sahar Sharifzadeh
- Division of Materials Science and Engineering, Boston University, Boston, Massachusetts 02215, United States.,Department of Electrical and Computer Engineering, Boston University, Boston, Massachusetts 02215, United States
| | | |
Collapse
|
32
|
Coupled- and Independent-Trajectory Approaches Based on the Exact Factorization Using the PyUNIxMD Package. Top Curr Chem (Cham) 2022; 380:8. [PMID: 35083549 DOI: 10.1007/s41061-021-00361-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Accepted: 11/11/2021] [Indexed: 10/19/2022]
Abstract
We present mixed quantum-classical approaches based on the exact factorization framework. The electron-nuclear correlation term in the exact factorization enables us to deal with quantum coherences by accounting for electronic and nuclear nonadiabatic couplings effectively within classical nuclei approximation. We compare coupled- and independent-trajectory approximations with each other to understand algorithms in description of the bifurcation of nuclear wave packets and the correct spatial distribution of electronic wave functions along with nuclear trajectories. Finally, we show numerical results for comparisons of coupled- and independent-trajectory approaches for the photoisomerization of a protonated Schiff base from excited state molecular dynamics (ESMD) simulations with the recently developed Python-based ESMD code, namely, the PyUNIxMD program.
Collapse
|
33
|
Akimov AV. Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain. J Phys Chem Lett 2021; 12:12119-12128. [PMID: 34913701 DOI: 10.1021/acs.jpclett.1c03823] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
34
|
Esch MP, Levine BG. An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics. J Chem Phys 2021; 155:214101. [PMID: 34879667 DOI: 10.1063/5.0070686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.
Collapse
Affiliation(s)
- Michael P Esch
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
35
|
Galindo JF, Freixas VM, Tretiak S, Fernandez-Alberti S. Back-and-Forth Energy Transfer during Electronic Relaxation in a Chlorin-Perylene Dyad. J Phys Chem Lett 2021; 12:10394-10401. [PMID: 34669398 DOI: 10.1021/acs.jpclett.1c03034] [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/13/2023]
Abstract
Donor-acceptor dyads represent a practical approach to tuning the photophysical properties of linear conjugated polymers in materials chemistry. Depending on the absorption wavelength, the acceptor and donor roles can be interchanged, and as such, the directionality of the energy transfer can be controlled. Herein, nonadiabatic excited state molecular dynamics simulations have been performed in an arylethylene-linked perylene-chlorin dyad. After an initial photoexcitation at the Soret band of chlorin, we observe an ultrafast sequential electronic relaxation to the lowest excited state. This process is accomplished through an efficient round-trip chlorin-to-perylene-to-chlorin energy transfer. It is characterized by successive intermittent localized and delocalized vibronic dynamics. Nonradiative relaxation takes place mainly through energy transfer events with perylene acting as a "heat sink" through which the nonradiative relaxation is efficiently funneled, and the excess energy is dispersed in a larger space of vibrational degrees of freedom. Thus, our findings suggest the use of donor-acceptor dyads as a useful strategy when one needs to deactivate an electronic excitation.
Collapse
Affiliation(s)
- Johan F Galindo
- Department of Chemistry, Universidad Nacional de Colombia, Bogotá 111321, Colombia
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - 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
|
36
|
Akimov AV. Excited state dynamics in monolayer black phosphorus revisited: Accounting for many-body effects. J Chem Phys 2021; 155:134106. [PMID: 34624981 DOI: 10.1063/5.0065606] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dynamics of electron-hole recombination in pristine and defect-containing monolayer black phosphorus (ML-BP) has been studied computationally by several groups relying on the one-particle description of electronic excited states. Our recent developments enabled a more sophisticated and accurate treatment of excited states dynamics in systems with pronounced excitonic effects, including 2D materials such as ML-BP. In this work, I present a comprehensive characterization of optoelectronic properties and nonadiabatic dynamics of the ground state recovery in pristine and divacancy-containing ML-BP, relying on the linear-response time-dependent density functional theory description of excited states combined with several trajectory surface hopping methodologies and decoherence correction schemes. This work presents a revision and new implementation of the decoherence-induced surface hopping methodology. Several popular algorithms for nonadiabatic dynamics algorithms are assessed. The kinetics of nonradiative relaxation of lower-lying excited states in ML-BP systems is revised considering the new methodological developments. A general mechanism that explains the sensitivity of the nonradiative dynamics to the presence of divacancy defect in ML-BP is proposed. According to this mechanism, the excited states' relaxation may be inhibited by the presence of energetically close higher-energy states if electronic decoherence is present in the system.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, USA
| |
Collapse
|
37
|
Freixas VM, Wilhelm P, Nelson T, Hinderer F, Höger S, Tretiak S, Lupton JM, Fernandez-Alberti S. Excitation Energy Transfer between bodipy Dyes in a Symmetric Molecular Excitonic Seesaw. J Phys Chem A 2021; 125:8404-8416. [PMID: 34542292 DOI: 10.1021/acs.jpca.1c06332] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We examine the redistribution of energy between electronic and vibrational degrees of freedom that takes place between a π-conjugated oligomer, a phenylene-butadiynylene, and two identical boron-dipyrromethene (bodipy) end-caps using femtosecond transient absorption spectroscopy, single-molecule spectroscopy, and nonadiabatic excited-state molecular dynamics (NEXMD) modeling techniques. The molecular structure represents an excitonic seesaw in that the excitation energy on the oligomer backbone can migrate to either one end-cap or the other, but not to both. The NEXMD simulations closely reproduce the characteristic time scale for redistribution of electronic and vibrational energy of 2.2 ps and uncover the vibrational modes contributing to the intramolecular relaxation. The calculations indicate that the dihedral angle between the bodipy dye and the oligomer change upon excitation of the oligomer. Single-molecule experiments reveal a difference in photoluminescence lifetime of the bodipy dyes depending on whether they are excited by direct absorption or by redistribution of energy from the backbone. This difference in lifetime may be attributed to the difference in dihedral angle. The simulations also suggest that a strong coupling can occur between the two end-caps, giving rise to a reversible shuttling of excitation energy between them. Strong coupling should lead to a pronounced loss in polarization memory of the fluorescence since the oligomer backbone tends to be slightly distorted and the two bodipy transition dipoles have different orientations. A sensitive single-molecule technique is presented to test for such coupling. However, although redistribution of electronic and vibrational energy between the end-caps can occur, it appears to be unidirectional and irreversible, suggesting that an additional localization mechanism is at play which is, as yet, not fully accounted for in the simulations.
Collapse
Affiliation(s)
- Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Philipp Wilhelm
- Institut für Angewandte und Experimentelle Physik, Universität Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | - Tammie Nelson
- Theoretical Division and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Florian Hinderer
- Kekulé-Institut für Organische Chemie und Biochemie, Universität Bonn, Gerhard-Domagk-Straße 1, 53121 Bonn, Germany
| | - Sigurd Höger
- Kekulé-Institut für Organische Chemie und Biochemie, Universität Bonn, Gerhard-Domagk-Straße 1, 53121 Bonn, Germany
| | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies (CINT), 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
| | | |
Collapse
|
38
|
Lee IS, Ha JK, Han D, Kim TI, Moon SW, Min SK. PyUNIxMD: A Python-based excited state molecular dynamics package. J Comput Chem 2021; 42:1755-1766. [PMID: 34197646 PMCID: PMC8362049 DOI: 10.1002/jcc.26711] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 06/10/2021] [Accepted: 06/16/2021] [Indexed: 01/17/2023]
Abstract
Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Jong-Kwon Ha
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Daeho Han
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Tae In Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Sung Wook Moon
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| |
Collapse
|
39
|
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
|
40
|
Alfonso-Hernandez L, Oldani N, Athanasopoulos S, Lupton JM, Tretiak S, Fernandez-Alberti S. Photoinduced Energy Transfer in Linear Guest-Host Chromophores: A Computational Study. J Phys Chem A 2021; 125:5303-5313. [PMID: 34106721 DOI: 10.1021/acs.jpca.1c02644] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Polymer-based guest-host systems represent a promising class of materials for efficient light-emitting diodes. The energy transfer from the polymer host to the guest is the key process in light generation. Therefore, microscopic descriptions of the different mechanisms involved in the energy transfer can contribute to enlighten the basis of the highly efficient light harvesting observed in this kind of materials. Herein, the nature of intramolecular energy transfer in a dye-end-capped conjugated polymer is explored by using atomistic nonadiabatic excited-state molecular dynamics. Linear perylene end-capped (PEC) polyindenofluorenes (PIF), consisting of n (n = 2, 4, and 6) repeat units, i.e., PEC-PIFn oligomers, are considered as model systems. After photoexcitation at the oligomer absorption maximum, an initial exciton becomes self-trapped on one of the monomer units (donors). Thereafter, an efficient ultrafast through-space energy transfer from this unit to the perylene acceptor takes place. We observe that this energy transfer occurs equally well from any monomer unit on the chain. Effective specific vibronic couplings between each monomer and the acceptor are identified. These oligomer → end-cap energy transfer steps do not match with the rates predicted by Förster-type energy transfer. The through-space and through-bond mechanisms are two distinct channels of energy transfer. The former dominates the overall process, whereas the through-bond energy transfer between indenofluorene monomer units along the oligomer backbone only makes a minor contribution.
Collapse
Affiliation(s)
- L Alfonso-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - N Oldani
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - S Athanasopoulos
- Departamento de Física, Universidad Carlos III de Madrid, Avenida Universidad 30, 28911 Leganés, Madrid, Spain
| | - J M Lupton
- Institut für Angewandte und Experimentelle Physik, Universität Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | - 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
| | - S Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| |
Collapse
|
41
|
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
|
42
|
Smith B, Shakiba M, Akimov AV. Crystal Symmetry and Static Electron Correlation Greatly Accelerate Nonradiative Dynamics in Lead Halide Perovskites. J Phys Chem Lett 2021; 12:2444-2453. [PMID: 33661640 DOI: 10.1021/acs.jpclett.0c03799] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Using a recently developed many-body nonadiabatic molecular dynamics (NA-MD) framework for large condensed matter systems, we study the phonon-driven nonradiative relaxation of excess electronic excitation energy in cubic and tetragonal phases of the lead halide perovskite CsPbI3. We find that the many-body treatment of the electronic excited states significantly changes the structure of the excited states' coupling, promotes a stronger nonadiabatic coupling of states, and ultimately accelerates the relaxation dynamics relative to the single-particle description of excited states. The acceleration of the nonadiabatic dynamics correlates with the degree of configurational mixing, which is controlled by the crystal symmetry. The higher-symmetry cubic phase of CsPbI3 exhibits stronger configuration mixing than does the tetragonal phase and subsequently yields faster nonradiative dynamics. Overall, using a many-body treatment of excited states and accounting for decoherence dynamics are important for closing the gap between the computationally derived and experimentally measured nonradiative excitation energy relaxation rates.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Mohammad Shakiba
- Department of Materials Science and Engineering, Shahid Bahonar University of Kerman, Kerman, Iran
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
43
|
Sindhu A, Jain A. Benchmarking the Surface Hopping Method to Include Nuclear Quantum Effects. J Chem Theory Comput 2021; 17:655-665. [PMID: 33432812 DOI: 10.1021/acs.jctc.0c01065] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have benchmarked the surface hopping method to capture nuclear quantum effects in the spin-Boson model in the deep tunneling regime. The thermal populations and the rate constants calculated using the surface hopping method are compared with those calculated using Boltzmann theory and Fermi's golden rule, respectively. Additionally, we have proposed a simple kinetic model that partially includes nuclear quantum effects within Marcus theory, and the results of the surface hopping method are analyzed under the framework of this simple kinetic model. A broad range of parameters are investigated to identify the regimes for the successes and failures of the surface hopping method. This work shows that with the accurate treatment of decoherence and velocity reversal, surface hopping can generally capture the nuclear quantum effects in the deep tunneling and weak diabatic coupling regime.
Collapse
Affiliation(s)
- Aarti Sindhu
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Amber Jain
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
44
|
Smith B, Shakiba M, Akimov AV. Nonadiabatic Dynamics in Si and CdSe Nanoclusters: Many-Body vs Single-Particle Treatment of Excited States. J Chem Theory Comput 2021; 17:678-693. [DOI: 10.1021/acs.jctc.0c01009] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| | - Mohammad Shakiba
- Department of Materials Science and Engineering, Shahid Bahonar University of Kerman, Kerman 76169-14111, Iran
| | - Alexey V. Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260 United States
| |
Collapse
|
45
|
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
|
46
|
De Sio A, Sommer E, Nguyen XT, Groß L, Popović D, Nebgen BT, Fernandez-Alberti S, Pittalis S, Rozzi CA, Molinari E, Mena-Osteritz E, Bäuerle P, Frauenheim T, Tretiak S, Lienau C. Intermolecular conical intersections in molecular aggregates. NATURE NANOTECHNOLOGY 2021; 16:63-68. [PMID: 33199882 DOI: 10.1038/s41565-020-00791-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 10/07/2020] [Indexed: 06/11/2023]
Abstract
Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.
Collapse
Affiliation(s)
- Antonietta De Sio
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany.
| | - Ephraim Sommer
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| | - Xuan Trung Nguyen
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| | - Lynn Groß
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - Duško Popović
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | | | - Sebastian Fernandez-Alberti
- National University of Quilmes/CONICET, Department of Science and Technology, Bernal (B1876BXD), Buenos Aires Province, Argentina
| | | | | | - Elisa Molinari
- Istituto Nanoscienze-CNR, Modena, Italy
- Università di Modena e Reggio Emilia, Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Modena, Italy
| | - Elena Mena-Osteritz
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | - Peter Bäuerle
- Institut für Organische Chemie II und Neue Materialien, Universität Ulm, Ulm, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
- Computational Science Research Center, Beijing and Computational Science and Applied Research Institute Shenzhen, Shenzhen, China
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Christoph Lienau
- Institut für Physik and Center of Interface Science, Carl von Ossietzky Universität, Oldenburg, Germany
| |
Collapse
|
47
|
Tang D, Shen L, Fang WH. Evaluation of mixed quantum-classical molecular dynamics on cis-azobenzene photoisomerization. Phys Chem Chem Phys 2021; 23:13951-13964. [PMID: 34142685 DOI: 10.1039/d1cp01374b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantitative prediction of nonadiabatic transitions between different electronic states is important to understand ultrafast processes in photochemistry. A variety of mixed quantum-classical molecular dynamics methods such as surface hopping and Ehrenfest mean-field have been developed. However, how to choose an appropriate one from a wide diversity of dynamics algorithms to study a realistic photochemical process is still unclear. In this work, we implemented 30 combinations of different mixed quantum-classical dynamics methods, including 24 surface hopping models with 8 decoherence corrections and 3 momentum rescaling strategies as well as 6 mean-field models. Then we performed numerical investigations by simulating the photoisomerization of cis-azobenzene combined with on-the-fly electronic structure calculations. Predictions of the S1 lifetime and the quantum yield of the photoproduct using different models are distinct. Surface hopping is more robust than mean-field in our test system. Moreover, the choice of momentum rescaling methods in surface hopping brings more significant changes than decoherence corrections, while a large discrepancy between simulation results with different mean-field algorithms has been observed.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| |
Collapse
|
48
|
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
|
49
|
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
|
50
|
Negrin-Yuvero H, Freixas VM, Rodriguez-Hernandez B, Rojas-Lorenzo G, Tretiak S, Bastida A, Fernandez-Alberti S. Photoinduced Dynamics with Constrained Vibrational Motion: FrozeNM Algorithm. J Chem Theory Comput 2020; 16:7289-7298. [PMID: 33201709 DOI: 10.1021/acs.jctc.0c00930] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. Herein, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.
Collapse
Affiliation(s)
- H Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - V M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - B Rodriguez-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - G Rojas-Lorenzo
- Departamento de Física Atómica y Molecular, Instituto Superior de Tecnologías y Ciencias Aplicadas, Universidad de La Habana , La Habana, Cuba
| | - 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
| | - A Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - S Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| |
Collapse
|