1
|
Miller ER, Hoehn SJ, Kumar A, Jiang D, Parker SM. Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory. J Chem Phys 2024; 161:034105. [PMID: 39007373 DOI: 10.1063/5.0203679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 06/26/2024] [Indexed: 07/16/2024] Open
Abstract
We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96 ps (S2 → S1) and 498 fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548 ps decay and (ii) a 1.69 ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.
Collapse
Affiliation(s)
- Ericka Roy Miller
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Sean J Hoehn
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Abhijith Kumar
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Dehua Jiang
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Ave., Cleveland, Ohio 44106, USA
| |
Collapse
|
2
|
Vindel-Zandbergen P, González-Vázquez J. Non-adiabatic dynamics of photoexcited cyclobutanone: Predicting structural measurements from trajectory surface hopping with XMS-CASPT2 simulations. J Chem Phys 2024; 161:024104. [PMID: 38984954 DOI: 10.1063/5.0203722] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 06/12/2024] [Indexed: 07/11/2024] Open
Abstract
Over the years, theoretical calculations and scalable computer simulations have complemented ultrafast experiments, as they offer the advantage of overcoming experimental restrictions and having access to the whole dynamics. This synergy between theory and experiment promises to yield a deeper understanding of photochemical processes, offering valuable insights into the behavior of complex systems at the molecular level. However, the ability of theoretical models to predict ultrafast experimental outcomes has remained largely unexplored. In this work, we aim to predict the electron diffraction signals of an upcoming ultrafast photochemical experiment using high-level electronic structure calculations and non-adiabatic dynamics simulations. In particular, we perform trajectory surface hopping with extended multi-state complete active space with second order perturbation simulations for understanding the photodissociation of cyclobutanone (CB) upon excitation at 200 nm. Spin-orbit couplings are considered for investigating the role of triplet states. Our simulations capture the bond cleavage after ultrafast relaxation from the 3s Rydberg state, leading to the formation of the previously observed primary photoproducts: CO + cyclopropane/propene (C3 products), ketene, and ethene (C2 products). The ratio of the C3:C2 products is found to be about 1:1. Within 700 fs, the majority of trajectories transition to their electronic ground state, with a small fraction conserving the initial cyclobutanone ring structure. We found a minimal influence of triplet states during the early stages of the dynamics, with their significance increasing at later times. We simulate MeV-ultrafast electron diffraction (UED) patterns from our trajectory results, linking the observed features with specific photoproducts and the underlying structural dynamics. Our analysis reveals highly intense features in the UED signals corresponding to the photochemical processes of CB. These features offer valuable insights into the experimental monitoring of ring opening dynamics and the formation of C3 and C2 photoproducts.
Collapse
Affiliation(s)
| | - Jesús González-Vázquez
- Departamento de Química, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
- Institute of Advanced Chemistry (IADChem), Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
3
|
Yao Q, Ji Q, Li X, Zhang Y, Chen X, Ju MG, Liu J, Wang J. Machine Learning Accelerates Precise Excited-State Potential Energy Surface Calculations on a Quantum Computer. J Phys Chem Lett 2024; 15:7061-7068. [PMID: 38950102 DOI: 10.1021/acs.jpclett.4c01445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2024]
Abstract
Electronically excited-state problems represent a crucial research field in quantum chemistry, closely related to numerous practical applications in photophysics and photochemistry. The emerging of quantum computing provides a promising computational paradigm to solve the Schrödinger equation for predicting potential energy surfaces (PESs). Here, we present a deep neural network model to predict parameters of the quantum circuits within the framework of variational quantum deflation and subspace search variational quantum eigensolver, which are two popular excited-state algorithms to implement on a quantum computer. The new machine learning-assisted algorithm is employed to study the excited-state PESs of small molecules, achieving highly accurate predictions. We then apply this algorithm to study the excited-state properties of the ArF system, which is essential to a gas laser. Through this study, we believe that with future advancements in hardware capabilities, quantum computing could be harnessed to solve excited-state problems for a broad range of systems.
Collapse
Affiliation(s)
- Qianjun Yao
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Qun Ji
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Xiaopeng Li
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei 230026, China
| | - Yehui Zhang
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Xinyu Chen
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Ming-Gang Ju
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Jie Liu
- Hefei National Laboratory, Hefei 230088, China
| | - Jinlan Wang
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
- Suzhou Laboratory, Suzhou 215009, China
| |
Collapse
|
4
|
Gandon A, Baiardi A, Ollitrault P, Tavernelli I. Nonadiabatic Molecular Dynamics with Fermionic Subspace-Expansion Algorithms on Quantum Computers. J Chem Theory Comput 2024. [PMID: 38967621 DOI: 10.1021/acs.jctc.4c00233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/06/2024]
Abstract
We introduce a novel computational framework for excited-state molecular quantum dynamics simulations driven by quantum-computing-based electronic-structure calculations. This framework leverages the fewest-switches surface-hopping method for simulating the nuclear dynamics and calculates the required excited-state transition properties with different flavors of the quantum subspace expansion and quantum equation-of-motion algorithms. We apply our method to simulate the collision reaction between a hydrogen atom and a hydrogen molecule. For this system, we critically compare the accuracy and efficiency of different quantum subspace expansion and equation-of-motion algorithms and show that only methods that can capture both weak and strong electron correlation effects can properly describe the nonadiabatic effects that tune the reactive event.
Collapse
Affiliation(s)
- Anthony Gandon
- IBM Quantum, IBM Research - Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
- Institute for Theoretical Physics, ETH Zürich, Wolfgang-Pauli-Str. 27, 8093 Zürich, Switzerland
| | - Alberto Baiardi
- IBM Quantum, IBM Research - Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | | | - Ivano Tavernelli
- IBM Quantum, IBM Research - Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
5
|
Zhang L, Pios SV, Martyka M, Ge F, Hou YF, Chen Y, Chen L, Jankowska J, Barbatti M, Dral PO. MLatom Software Ecosystem for Surface Hopping Dynamics in Python with Quantum Mechanical and Machine Learning Methods. J Chem Theory Comput 2024; 20:5043-5057. [PMID: 38836623 DOI: 10.1021/acs.jctc.4c00468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.
Collapse
Affiliation(s)
- Lina Zhang
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Sebastian V Pios
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Mikołaj Martyka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Fuchun Ge
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lipeng Chen
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| | - Pavlo O Dral
- College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| |
Collapse
|
6
|
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
|
7
|
Polonius S, Lehrner D, González L, Mai S. Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations. J Chem Theory Comput 2024; 20:4738-4750. [PMID: 38768386 PMCID: PMC11171268 DOI: 10.1021/acs.jctc.4c00169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/22/2024] [Accepted: 05/07/2024] [Indexed: 05/22/2024]
Abstract
Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - David Lehrner
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
8
|
Shu Y, Truhlar DG. Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients. J Chem Theory Comput 2024; 20:4396-4426. [PMID: 38819014 DOI: 10.1021/acs.jctc.4c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
9
|
Bustamante CM, Todorov T, Gadea ED, Tarasi F, Stella L, Horsfield A, Scherlis DA. Modeling the electroluminescence of atomic wires from quantum dynamics simulations. J Chem Phys 2024; 160:214102. [PMID: 38828814 DOI: 10.1063/5.0201447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
Static and time-dependent quantum-mechanical approaches have been employed in the literature to characterize the physics of light-emitting molecules and nanostructures. However, the electromagnetic emission induced by an input current has remained beyond the realm of molecular simulations. This is the challenge addressed here with the help of an equation of motion for the density matrix coupled to a photon bath based on a Redfield formulation. This equation is evolved within the framework of the driven-Liouville von Neumann approach, which incorporates open boundaries by introducing an applied bias and a circulating current. The dissipated electromagnetic power can be computed in this context from the time derivative of the energy. This scheme is applied in combination with a self-consistent tight-binding Hamiltonian to investigate the effects of bias and molecular size on the electroluminescence of metallic and semiconducting chains. For the latter, a complex interplay between bias and molecular length is observed: there is an optimal number of atoms that maximizes the emitted power at high voltages but not at low ones. This unanticipated behavior can be understood in terms of the band bending produced along the semiconducting chain, a phenomenon that is captured by the self-consistency of the method. A simple analytical model is proposed that explains the main features revealed by the simulations. The methodology, applied here at a self-consistent tight-binding level but extendable to more sophisticated Hamiltonians such as density functional tight binding and time dependent density functional theory, promises to be helpful for quantifying the power and quantum efficiency of nanoscale electroluminescent devices.
Collapse
Affiliation(s)
- Carlos M Bustamante
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg 22761, Germany
| | - Tchavdar Todorov
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, United Kingdom
| | - Esteban D Gadea
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Facundo Tarasi
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Lorenzo Stella
- Centre for Light-Matter Interactions, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, United Kingdom
| | - Andrew Horsfield
- Department of Materials, Thomas Young Centre, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Damián A Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| |
Collapse
|
10
|
Qiu T, Bian X, Tao Z, Subotnik JE. A fast and smooth one-electron approach for investigating charge transfer states and D1-D0 crossings for systems with odd numbers of electrons. J Chem Phys 2024; 160:214115. [PMID: 38832731 DOI: 10.1063/5.0209036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Accepted: 05/13/2024] [Indexed: 06/05/2024] Open
Abstract
We propose an efficient version of ensemble Hartree-Fock/density functional theory to calculate a set of two charge-transfer states for systems with odd-numbers of electrons. The approach does require definitions of donor/acceptor fragments; however, the approach is not very sensitive to such definitions-even in the limit of very strong electronic coupling. The key ansatz is that, by mandating that the vector space spanned by the active orbitals projects equally onto the donor and acceptor fragments, such a constraint eliminates all intra-molecular local excitations and makes it far easier to generate potential energy surfaces that are smooth over a wide region of configuration space. The method is fast, working with only two electron configurations, and should be useful for ab initio non-adiabatic dynamics in the near future.
Collapse
Affiliation(s)
- Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
11
|
Montero-Cabrera LA, Montero-Alejo AL, Aspuru-Guzik A, García de la Vega JM, Piris M, Díaz-Fernández LA, Pérez-Badell Y, Guerra-Barroso A, Alfonso-Ramos JE, Rodríguez J, Fuentes ME, de Armas CM. Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties. J Chem Phys 2024; 160:214108. [PMID: 38828812 DOI: 10.1063/5.0208809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb-exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.
Collapse
Affiliation(s)
- Luis A Montero-Cabrera
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
- Donostia International Physics Center (DIPC), 20018 Donostia - San Sebastián, Basque Country, Spain
| | - Ana L Montero-Alejo
- Departamento de Física, Facultad de Ciencias Naturales, Matemática y del Medio Ambiente (FCNMM), Universidad Tecnológica Metropolitana; Ñuñoa, Santiago 7800002, Chile
| | - Alan Aspuru-Guzik
- Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | | | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia - San Sebastián, Basque Country, Spain
| | - Lourdes A Díaz-Fernández
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Yoana Pérez-Badell
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Alberto Guerra-Barroso
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Javier E Alfonso-Ramos
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Javier Rodríguez
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - María E Fuentes
- Facultad de Ciencias Químicas, Universidad Autónoma de Chihuahua, Chihuahua, 31100 Chihuahua, Mexico
| | - Carlos M de Armas
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| |
Collapse
|
12
|
Araujo L, Lasser C, Schmidt B. FSSH-2: Fewest Switches Surface Hopping with Robust Switching Probability. J Chem Theory Comput 2024; 20:3413-3419. [PMID: 38696709 DOI: 10.1021/acs.jctc.4c00089] [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/2024]
Abstract
This study introduces the FSSH-2 scheme, a redefined and numerically stable adiabatic Fewest Switches Surface Hopping (FSSH) method for mixed quantum-classical dynamics. It reformulates the standard FSSH hopping probability without using nonadiabatic coupling vectors and allows for numerical time integration with larger step sizes. The advantages of FSSH-2 are demonstrated by numerical experiments for five different model systems in one and two spatial dimensions with up to three electronic states.
Collapse
Affiliation(s)
- Leonardo Araujo
- Department of Mathematics, TUM School of Computation, Information and Technology, Technische Universität München, Boltzmannstr. 3, 85748 Garching bei München, Germany
| | - Caroline Lasser
- Department of Mathematics, TUM School of Computation, Information and Technology, Technische Universität München, Boltzmannstr. 3, 85748 Garching bei München, Germany
| | - Burkhard Schmidt
- Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany
| |
Collapse
|
13
|
Liu Z, Song Z, Sun X. All-Atom Photoinduced Charge Transfer Dynamics in Condensed Phase via Multistate Nonlinear-Response Instantaneous Marcus Theory. J Chem Theory Comput 2024; 20:3993-4006. [PMID: 38657208 PMCID: PMC11099976 DOI: 10.1021/acs.jctc.4c00010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/30/2024] [Accepted: 04/08/2024] [Indexed: 04/26/2024]
Abstract
Photoinduced charge transfer (CT) in the condensed phase is an essential component in solar energy conversion, but it is challenging to simulate such a process on the all-atom level. The traditional Marcus theory has been utilized for obtaining CT rate constants between pairs of electronic states but cannot account for the nonequilibrium effects due to the initial nuclear preparation. The recently proposed instantaneous Marcus theory (IMT) and its nonlinear-response formulation allow for incorporating the nonequilibrium nuclear relaxation to electronic transition between two states after the photoexcitation from the equilibrium ground state and provide the time-dependent rate coefficient. In this work, we extend the nonlinear-response IMT method for treating photoinduced CT among general multiple electronic states and demonstrate it in the organic photovoltaic carotenoid-porphyrin-fullerene triad dissolved in explicit tetrahydrofuran solvent. All-atom molecular dynamics simulations were employed to obtain the time correlation functions of energy gaps, which were used to generate the IMT-required time-dependent averages and variances of the relevant energy gaps. Our calculations show that the multistate IMT could capture the significant nonequilibrium effects due to the initial nuclear state preparation, and this is corroborated by the substantial differences between the population dynamics predicted by the multistate IMT and the Marcus theory, where the Marcus theory underestimates the population transfer. The population dynamics by multistate IMT is also shown to have a better agreement with the all-atom nonadiabatic mapping dynamics than the Marcus theory does. Because the multistate nonlinear-response IMT is straightforward and cost-effective in implementation and accounts for the nonequilibrium nuclear effects, we believe this method offers a practical strategy for studying charge transfer dynamics in complex condensed-phase systems.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Zailing Song
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
14
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
15
|
Wang K, Xu Y, Xie X, Ma H. Theoretical investigation of distal charge separation in a perylenediimide trimer. J Chem Phys 2024; 160:164303. [PMID: 38647303 DOI: 10.1063/5.0205671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 04/03/2024] [Indexed: 04/25/2024] Open
Abstract
An exciton-phonon (ex-ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex-ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore's reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.
Collapse
Affiliation(s)
- Ke Wang
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xiaoyu Xie
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
16
|
Vandaele E, Mališ M, Luber S. The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols. J Chem Theory Comput 2024; 20:3258-3272. [PMID: 38606908 PMCID: PMC11044273 DOI: 10.1021/acs.jctc.3c01400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 03/27/2024] [Accepted: 03/29/2024] [Indexed: 04/13/2024]
Abstract
The photochemistry of nitrophenols is a source of smog as nitrous acid is formed from their photolysis. Nevertheless, computational studies of the photochemistry of these widespread toxic molecules are scarce. In this work, the initial photodeactivation of ortho-nitrophenol and para-nitrophenol is modeled, both in gas phase and in aqueous solution to simulate atmospheric and aerosol environments. A large number of excited states, six for ortho-nitrophenol and 11 for para-nitrophenol, have been included and were all populated during the decay. Moreover, periodic time-dependent density functional theory (TDDFT) is used for both the explicitly included solvent and the solute. A comparison to periodic QM/MM (TDDFT/MM), with electrostatic embedding, is made, showing notable differences between the decays of solvated nitrophenols simulated with QM/MM and full (TD)DFT. A reduced intersystem crossing in aqueous solution could be observed thanks to the surface hopping approach using explicit, periodic TDDFT solvation including spin-orbit couplings.
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
17
|
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
|
18
|
Mukherjee S, Mattos RS, Toldo JM, Lischka H, Barbatti M. Prediction Challenge: Simulating Rydberg photoexcited cyclobutanone with surface hopping dynamics based on different electronic structure methods. J Chem Phys 2024; 160:154306. [PMID: 38624122 DOI: 10.1063/5.0203636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 03/28/2024] [Indexed: 04/17/2024] Open
Abstract
This research examines the nonadiabatic dynamics of cyclobutanone after excitation into the n → 3s Rydberg S2 state. It stems from our contribution to the Special Topic of the Journal of Chemical Physics to test the predictive capability of computational chemistry against unseen experimental data. Decoherence-corrected fewest-switches surface hopping was used to simulate nonadiabatic dynamics with full and approximated nonadiabatic couplings. Several simulation sets were computed with different electronic structure methods, including a multiconfigurational wavefunction [multiconfigurational self-consistent field (MCSCF)] specially built to describe dissociative channels, multireference semiempirical approach, time-dependent density functional theory, algebraic diagrammatic construction, and coupled cluster. MCSCF dynamics predicts a slow deactivation of the S2 state (10 ps), followed by an ultrafast population transfer from S1 to S0 (<100 fs). CO elimination (C3 channel) dominates over C2H4 formation (C2 channel). These findings radically differ from the other methods, which predicted S2 lifetimes 10-250 times shorter and C2 channel predominance. These results suggest that routine electronic structure methods may hold low predictive power for the outcome of nonadiabatic dynamics.
Collapse
Affiliation(s)
| | - Rafael S Mattos
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
| | - Josene M Toldo
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| |
Collapse
|
19
|
Makhov DV, Armstrong G, Chuang HH, Ambalampitiya H, Lemishko K, Mohr S, Nelson A, Tennyson J, Shalashilin D. Dissociation of Hydrofluorocarbon Molecules after Electron Impact in Plasma. J Phys Chem Lett 2024; 15:3404-3411. [PMID: 38502942 PMCID: PMC10983062 DOI: 10.1021/acs.jpclett.4c00348] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 03/08/2024] [Accepted: 03/12/2024] [Indexed: 03/21/2024]
Abstract
The process of dissociation for two hydrofluorocarbon molecules in low triplet states excited by electron impact in plasma is investigated by ab initio molecular dynamics (AIMD). The interest in the dissociation of hydrofluorocarbons in plasma is motivated by their role in plasma etching in microelectronic technologies. Dissociation of triplet states is very fast, and the reaction products can be predicted. In this work, it was found that higher triplet states relax into the lowest triplet state within a few femtoseconds due to nonadiabatic dynamics, such that the simplest ab initio MD on the lowest triplet state seems to give a reasonable estimate of the reaction channels branching ratios. We provide evidence of the existence of simple rules for the dissociation of hydrofluorocarbon molecules in triplet states. For molecules with a double bond, the bonds adjacent to the double bond dissociate faster than the other bonds.
Collapse
Affiliation(s)
- Dmitry V. Makhov
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School
of Mathematics, University of Bristol, Fry Building, Woodland Road, Bristol BS8 1UG, United Kingdom
| | - Gregory Armstrong
- Quantemol
Ltd., 320 City Road,
The Angel, London EC1V
2NZ, United Kingdom
| | - Hsiao-Han Chuang
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School
of Mathematics, University of Bristol, Fry Building, Woodland Road, Bristol BS8 1UG, United Kingdom
| | | | - Kateryna Lemishko
- Quantemol
Ltd., 320 City Road,
The Angel, London EC1V
2NZ, United Kingdom
| | - Sebastian Mohr
- Quantemol
Ltd., 320 City Road,
The Angel, London EC1V
2NZ, United Kingdom
| | - Anna Nelson
- Quantemol
Ltd., 320 City Road,
The Angel, London EC1V
2NZ, United Kingdom
| | - Jonathan Tennyson
- Quantemol
Ltd., 320 City Road,
The Angel, London EC1V
2NZ, United Kingdom
- Department
of Physics and Astronomy, University College
London, London WC1E 6BT, United
Kingdom
| | | |
Collapse
|
20
|
Pios SV, Gelin MF, Ullah A, Dral PO, Chen L. Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra. J Phys Chem Lett 2024; 15:2325-2331. [PMID: 38386692 DOI: 10.1021/acs.jpclett.4c00107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024]
Abstract
Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.
Collapse
Affiliation(s)
- Sebastian V Pios
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou, Zhejiang 310018, People's Republic of China
| | - Arif Ullah
- School of Physics and Optoelectronic Engineering, Anhui University, Hefei, Anhui 230601, People's Republic of China
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Lipeng Chen
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| |
Collapse
|
21
|
Perpendicular state of stilbene unveiled with ultrafast Raman spectroscopy. Nat Chem 2024; 16:308-309. [PMID: 38191618 DOI: 10.1038/s41557-023-01398-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
|
22
|
Cao Y, Halls MD, Friesner RA. Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method. J Chem Phys 2024; 160:084106. [PMID: 38385510 DOI: 10.1063/5.0188277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 01/31/2024] [Indexed: 02/23/2024] Open
Abstract
A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc., 1540 Broadway, 24th Floor, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger Inc., 9868 Scranton, Suite 3200, San Diego, California 92121, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
23
|
Brey D, Burghardt I. Coherent Transient Localization Mechanism of Interchain Exciton Transport in Regioregular P3HT: A Quantum-Dynamical Study. J Phys Chem Lett 2024; 15:1836-1845. [PMID: 38334949 DOI: 10.1021/acs.jpclett.3c03290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Transient localization has been proposed as a transport mechanism in organic materials, for both charge carriers and excitons. Here, we characterize a quantum coherent transient localization mechanism using full quantum simulations of an H-aggregated model system representative of regioregular poly(3-hexylthiophene) (rrP3HT). A Frenkel-Holstein Hamiltonian parametrized from first principles is considered, including local high-frequency modes and anharmonic, site-correlated interchain modes. Quantum-dynamical calculations are carried out using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a 13-site system with 195 vibrational modes, under periodic boundary conditions. It is shown that temporary localization of exciton polarons alternates with resonant transfer driven by interchain modes. While the transport process is mainly determined by exciton-polarons at the low-energy band edge, persistent coupling with the excitonic manifold is observed, giving rise to a nonadiabatic excitonic flux. This elementary transport mechanism remains preserved for limited static disorder and gives way to Anderson localization when the static disorder becomes dominant.
Collapse
Affiliation(s)
- Dominik Brey
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| |
Collapse
|
24
|
Titov E. Visible Light Induced Exciton Dynamics and trans-to- cis Isomerization in Azobenzene Aggregates: Insights from Surface Hopping/Semiempirical Configuration Interaction Molecular Dynamics Simulations. ACS OMEGA 2024; 9:8520-8532. [PMID: 38405525 PMCID: PMC10882624 DOI: 10.1021/acsomega.3c09900] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 01/19/2024] [Accepted: 01/24/2024] [Indexed: 02/27/2024]
Abstract
Assemblies of photochromic molecules feature exciton states, which govern photochemical and photophysical processes in multichromophoric systems. Understanding the photoinduced dynamics of the assemblies requires nonadiabatic treatment involving multiple exciton states and numerous nuclear degrees of freedom, thus posing a challenge for simulations. In this work, we address this challenge for aggregates of azobenzene, a prototypical molecular switch, performing on-the-fly surface hopping calculations combined with semiempirical configuration interaction electronic structure and augmented with transition density matrix analysis to characterize exciton evolution. Specifically, we consider excitation of azobenzene tetramers in the nπ* absorption band located in the visible (blue) part of the electromagnetic spectrum, thus extending our recent work on dynamics after ππ* excitation corresponding to the ultraviolet region [Titov, J. Phys. Chem. C2023, 127, 13678-13688]. We find that the nπ* excitons, which are initially strongly localized by ground-state conformational disorder, undergo further (very strong) localization during short-time photodynamics. This excited-state localization process is extremely ultrafast, occurring within the first 10 fs of photodynamics. We observe virtually no exciton transfer of the localized excitons in the nπ* manifold. However, the transfer may occur via secondary pathways involving ππ* states or the ground state. Moreover, we find that the nπ* quantum yields of the trans-to-cis isomerization are reduced in the aggregated state.
Collapse
Affiliation(s)
- Evgenii Titov
- Institute of Chemistry, Theoretical
Chemistry, University of Potsdam, Karl-Liebknecht-Straße 24-25, 14476 Potsdam, Germany
| |
Collapse
|
25
|
Shi Y, Zhang G, Xiang C, Liu C, Hu J, Wang J, Ge R, Ma H, Niu Y, Xu Y. Defect-Engineering-Mediated Long-Lived Charge-Transfer Excited-State in Fe-Gallate Complex Improves Iron Cycle and Enables Sustainable Fenton-Like Reaction. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2305162. [PMID: 37708316 DOI: 10.1002/adma.202305162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 08/18/2023] [Indexed: 09/16/2023]
Abstract
Fenton reactions are inefficient because the Fe(II) catalyst cannot be recycled in time due to the lack of a rapid electron transport pathway. This results in huge H2 O2 wastage in industrial applications. Here, it is shown that a sustainable heterogeneous Fenton system is attainable by enhancing the ligand-to-metal charge-transfer (LMCT) excited-state lifetime in Fe-gallate complex. By engineering oxygen defects in the complex, the lifetime is improved from 10-90 ps. The lengthened lifetime ensures sufficient concentrations of excited-states for an efficient Fe cycle, realizing previously unattainable H2 O2 activation kinetics and hydroxyl radical (• OH) productivity. Spectroscopic and electrochemical studies show the cyclic reaction mechanism involves in situ Fe(II) regeneration and synchronous supply of oxygen atoms from water to recover dissociated Fe─O bonds. Trace amounts of this catalyst effectively destroy two drug-resistant bacteria even after eight reaction cycles. This work reveals the link among LMCT excited-state lifetime, Fe cycle, and catalytic activity and stability, with implications for de novo design of efficient and sustainable Fenton-like processes.
Collapse
Affiliation(s)
- Yanfeng Shi
- Institute of Biomedical Engineering, College of Life Sciences, Qingdao University, Qingdao, 266071, China
| | - Gong Zhang
- Center for Water and Ecology, State Key Joint Laboratory of Environment Simulation and Pollution Control School of Environment, Tsinghua University, Beijing, 100084, China
| | - Chao Xiang
- Center for Water and Ecology, State Key Joint Laboratory of Environment Simulation and Pollution Control School of Environment, Tsinghua University, Beijing, 100084, China
| | - Chengzhen Liu
- Institute of Biomedical Engineering, College of Life Sciences, Qingdao University, Qingdao, 266071, China
| | - Jun Hu
- School of Chemical Engineering, Northwest University, Xi' an, 710069, China
| | - Junhu Wang
- Center for Advanced Mössbauer Spectroscopy, Mössbauer Effect Data Center, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116000, China
| | - Rile Ge
- Center for Advanced Mössbauer Spectroscopy, Mössbauer Effect Data Center, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116000, China
| | - Haixia Ma
- School of Chemical Engineering, Northwest University, Xi' an, 710069, China
| | - Yusheng Niu
- Institute of Biomedical Engineering, College of Life Sciences, Qingdao University, Qingdao, 266071, China
| | - Yuanhong Xu
- Institute of Biomedical Engineering, College of Life Sciences, Qingdao University, Qingdao, 266071, China
| |
Collapse
|
26
|
Li X, Lubbers N, Tretiak S, Barros K, Zhang Y. Machine Learning Framework for Modeling Exciton Polaritons in Molecular Materials. J Chem Theory Comput 2024; 20:891-901. [PMID: 38168674 DOI: 10.1021/acs.jctc.3c01068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon, and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupled to a cavity mode is theoretically and computationally challenging. Recent advances in machine learning (ML) techniques have shown promising capabilities in modeling ground-state chemical systems. This work presents a general protocol to predict excited-state properties, such as energies, transition dipoles, and nonadiabatic coupling vectors with the hierarchically interacting particle neural network. ML predictions are then applied to compute the potential energy surfaces and electronic spectra of a prototype azomethane molecule in the collective coupling scenario. These computational tools provide a much-needed framework to model and understand many molecules' emerging excited-state polariton chemistry.
Collapse
Affiliation(s)
- Xinyang Li
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nicholas Lubbers
- Information Sciences, Computer, Computational, and Statistical Sciences 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
| | - Kipton Barros
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
27
|
LeCroy G, Ghosh R, Untilova V, Guio L, Stone KH, Brinkmann M, Luscombe C, Spano FC, Salleo A. Polaron absorption in aligned conjugated polymer films: breakdown of adiabatic treatments and going beyond the conventional mid-gap state model. MATERIALS HORIZONS 2024; 11:545-553. [PMID: 37982315 DOI: 10.1039/d3mh01278f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
This study provides the first experimental polarized intermolecular and intramolecular optical absorption components of field-induced polarons in regioregular poly(3-hexylthiophene-2,5-diyl), rr-P3HT, a polymer semiconductor. Highly aligned rr-P3HT thin films were prepared by a high temperature shear-alignment process that orients polymer backbones along the shearing direction. rr-P3HT in-plane molecular orientation was measured by electron diffraction, and out-of-plane orientation was measured through series of synchrotron X-ray scattering techniques. Then, with molecular orientation quantified, polarized charge modulation spectroscopy was used to probe mid-IR polaron absorption in the ℏω = 0.075 - 0.75 eV range and unambiguously assign intermolecular and intramolecular optical absorption components of hole polarons in rr-P3HT. This data represents the first experimental quantification of these polarized components and allowed long-standing theoretical predictions to be compared to experimental results. The experimental data is discrepant with predictions of polaron absorption based on an adiabatic framework that works under the Born-Oppenheimer approximation, but the data is entirely consistent with a more recent nonadiabatic treatment of absorption based on a modified Holstein Hamiltonian. This nonadiabatic treatment was used to show that both intermolecular and intramolecular polaron coherence break down at length scales significantly smaller than estimated structural coherence in either direction. This strongly suggests that polaron delocalization is fundamentally limited by energetic disorder in rr-P3HT.
Collapse
Affiliation(s)
- Garrett LeCroy
- Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305, USA.
| | - Raja Ghosh
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92037, USA
| | | | - Lorenzo Guio
- Department of Materials Science and Engineering, University of Washington, Seattle, WA, USA
| | - Kevin H Stone
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Martin Brinkmann
- Université de Strasbourg, CNRS, ICS UPR 22, F-67000 Strasbourg, France
| | - Christine Luscombe
- pi-Conjugated Polymers Unit, Okinawa Institute of Science and Technology, Onna, Okinawa 904-0495, Japan
| | - Frank C Spano
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Alberto Salleo
- Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
28
|
Dines A, Ellis M, Blumberger J. Stabilized coupled trajectory mixed quantum-classical algorithm with improved energy conservation: CTMQC-EDI. J Chem Phys 2023; 159:234118. [PMID: 38117021 DOI: 10.1063/5.0183589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 11/26/2023] [Indexed: 12/21/2023] Open
Abstract
Coupled trajectory mixed quantum-classical (CTMQC) dynamics is a rigorous approach to trajectory-based non-adiabatic dynamics, which has recently seen an improvement to energy conservation via the introduction of the CTMQC-E algorithm. Despite this, the method's two key quantities distinguishing it from Ehrenfest dynamics, the modified Born-Oppenheimer momentum and the quantum momentum, require regularization procedures in certain circumstances. Such procedures in the latter can cause instabilities, leading to undesirable effects, such as energy drift and spurious population transfer, which is expected to become increasingly prevalent when the system gets larger as such events would happen more frequently. We propose a further modification to CTMQC-E, which includes a redefinition of the quantum momentum, CTMQC-EDI (double intercept), such that it has no formal divergences. We then show for Tully models I-III and the double arch model that the algorithm has greatly improved total energy conservation and negligible spurious population transfer at all times, in particular in regions of strong non-adiabatic coupling. CTMQC-EDI, therefore, shows promise as a numerically robust non-adiabatic dynamics technique that accounts for decoherence from first principles and that is scalable to large molecular systems and materials.
Collapse
Affiliation(s)
- Aaron Dines
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Matthew Ellis
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
29
|
Mondal S, Chowdhury U, Dey S, Habib M, Mora Perez C, Frauenheim T, Sarkar R, Pal S, Prezhdo OV. Controlling Charge Carrier Dynamics in Porphyrin Nanorings by Optically Active Templates. J Phys Chem Lett 2023; 14:11384-11392. [PMID: 38078872 PMCID: PMC10749466 DOI: 10.1021/acs.jpclett.3c03304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 12/01/2023] [Accepted: 12/01/2023] [Indexed: 12/22/2023]
Abstract
Understanding the dynamics of photogenerated charge carriers is essential for enhancing the performance of solar and optoelectronic devices. Using atomistic quantum dynamics simulations, we demonstrate that a short π-conjugated optically active template can be used to control hot carrier relaxation, charge carrier separation, and carrier recombination in light-harvesting porphyrin nanorings. Relaxation of hot holes is slowed by 60% with an optically active template compared to that with an analogous optically inactive template. Both systems exhibit subpicosecond electron transfer from the photoactive core to the templates. Notably, charge recombination is suppressed 6-fold by the optically active template. The atomistic time-domain simulations rationalize these effects by the extent of electron and hole localization, modification of the density of states, participation of distinct vibrational motions, and changes in quantum coherence. Extension of the hot carrier lifetime and reduction of charge carrier recombination, without hampering charge separation, demonstrate a strategy for enhancing efficiencies of energy materials with optically active templates.
Collapse
Affiliation(s)
- Shrabanti Mondal
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
| | - Uttam Chowdhury
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
| | - Subhajit Dey
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
| | - Md Habib
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
- Department
of Chemistry, Sripat Singh College, Jiaganj 742122, India
| | - Carlos Mora Perez
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Thomas Frauenheim
- Bremen
Center
for Computational Materials Science, Universität
Bremen, Bremen 28359, Germany
- Beijing
Computational Science Research Center, Beijing 100193, China
- Shenzhen
JL Computational Science and Applied Research Institute, Shenzhen 518109, China
| | - Ritabrata Sarkar
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
- Bremen
Center
for Computational Materials Science, Universität
Bremen, Bremen 28359, Germany
| | - Sougata Pal
- Department
of Chemistry, University of Gour Banga, Malda 732103, India
| | - 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
|
30
|
Buttarazzi E, Perrella F, Rega N, Petrone A. Watching the Interplay between Photoinduced Ultrafast Charge Dynamics and Nuclear Vibrations. J Chem Theory Comput 2023; 19:8751-8766. [PMID: 37991892 PMCID: PMC10720350 DOI: 10.1021/acs.jctc.3c00855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 11/02/2023] [Accepted: 11/03/2023] [Indexed: 11/24/2023]
Abstract
Here is presented the ultrafast hole-electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4- (dcbpy = 4,4'-dicarboxy-2,2'-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS- ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.
Collapse
Affiliation(s)
- Edoardo Buttarazzi
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
| | - Fulvio Perrella
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Nadia Rega
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| | - Alessio Petrone
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| |
Collapse
|
31
|
do Casal MT, Veys K, Bousquet MHE, Escudero D, Jacquemin D. First-Principles Calculations of Excited-State Decay Rate Constants in Organic Fluorophores. J Phys Chem A 2023; 127:10033-10053. [PMID: 37988002 DOI: 10.1021/acs.jpca.3c06191] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
In this Perspective, we discuss recent advances made to evaluate from first-principles the excited-state decay rate constants of organic fluorophores, focusing on the so-called static strategy. In this strategy, one essentially takes advantage of Fermi's golden rule (FGR) to evaluate rate constants at key points of the potential energy surfaces, a procedure that can be refined in a variety of ways. In this way, the radiative rate constant can be straightforwardly obtained by integrating the fluorescence line shape, itself determined from vibronic calculations. Likewise, FGR allows for a consistent calculation of the internal conversion (related to the non-adiabatic couplings) in the weak-coupling regime and intersystem crossing rates, therefore giving access to estimates of the emission yields when no complex photophysical phenomenon is at play. Beyond outlining the underlying theories, we summarize here the results of benchmarks performed for various types of rates, highlighting that both the quality of the vibronic calculations and the accuracy of the relative energies are crucial to reaching semiquantitative estimates. Finally, we illustrate the successes and challenges in determining the fluorescence quantum yields using a series of organic fluorophores.
Collapse
Affiliation(s)
- Mariana T do Casal
- Department of Chemistry, Physical Chemistry and Quantum Chemistry Division, KU Leuven, 3001 Leuven, Belgium
| | - Koen Veys
- Department of Chemistry, Physical Chemistry and Quantum Chemistry Division, KU Leuven, 3001 Leuven, Belgium
| | | | - Daniel Escudero
- Department of Chemistry, Physical Chemistry and Quantum Chemistry Division, KU Leuven, 3001 Leuven, Belgium
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), FR-75005 Paris, France
| |
Collapse
|
32
|
Liu L, Li S, Zhang N, Shi Q, Liu K, Liu T, Huang Z, Ding L, Fang Y. Comparative Observation of Distinct Dynamic Stokes Shifts in Diaryl BODIPY Triads with Broadband Two-Photon Absorption. J Phys Chem B 2023; 127:10171-10178. [PMID: 37967951 DOI: 10.1021/acs.jpcb.3c06757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2023]
Abstract
Time-resolved evolution of excited states in the twist-conjugated chromophores is of great fundamental interest for photoluminescent applications. The four diaryl BODIPY triads modified with diverse end-cappers at 2,6-positions were investigated properly, and considerable two-photon absorption capabilities in the first biological spectral window were obtained. Fast relaxations from the initially twisted conformation to the planarized conformation in the excited state were resolved spectrally and kinetically, accompanied by the discernible phenomenon of the fluorescence dynamic Stokes shift (DSS). Along with increasing electron donating capabilities and solvent polarities, the characteristics of structural rearrangement and intramolecular charge transfer have been estimated by enhanced DSS behaviors. Especially, the blue-shifted DSS was rationalized as the sequence conversion between the planarized state and the twisted charge transfer state. A molecular-level picture for relaxation pathways in different polarities was depicted and supported by the theoretical simulations. Significant and fast structural motions in this work contribute to the excited-state dynamics and rational development of versatile BODIPY chromophores.
Collapse
Affiliation(s)
- Lu Liu
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Sheng Li
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Nan Zhang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Qiyuan Shi
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Ke Liu
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Taihong Liu
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Zhiyan Huang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Liping Ding
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| | - Yu Fang
- Key Laboratory of Applied Surface and Colloid Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, People's Republic of China
| |
Collapse
|
33
|
Weight BM, Li X, Zhang Y. Theory and modeling of light-matter interactions in chemistry: current and future. Phys Chem Chem Phys 2023; 25:31554-31577. [PMID: 37842818 DOI: 10.1039/d3cp01415k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2023]
Abstract
Light-matter interaction not only plays an instrumental role in characterizing materials' properties via various spectroscopic techniques but also provides a general strategy to manipulate material properties via the design of novel nanostructures. This perspective summarizes recent theoretical advances in modeling light-matter interactions in chemistry, mainly focusing on plasmon and polariton chemistry. The former utilizes the highly localized photon, plasmonic hot electrons, and local heat to drive chemical reactions. In contrast, polariton chemistry modifies the potential energy curvatures of bare electronic systems, and hence their chemistry, via forming light-matter hybrid states, so-called polaritons. The perspective starts with the basic background of light-matter interactions, molecular quantum electrodynamics theory, and the challenges of modeling light-matter interactions in chemistry. Then, the recent advances in modeling plasmon and polariton chemistry are described, and future directions toward multiscale simulations of light-matter interaction-mediated chemistry are discussed.
Collapse
Affiliation(s)
- Braden M Weight
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
- Department of Physics and Astronomy, University of Rochester, Rochester, NY, 14627, USA
| | - Xinyang Li
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| |
Collapse
|
34
|
Wasowicz TJ, Dąbkowska I, Kivimäki A, Richter R. Two-body dissociation of isoxazole following double photoionization - an experimental PEPIPICO and theoretical DFT and MP2 study. Phys Chem Chem Phys 2023; 25:31655-31666. [PMID: 37964643 DOI: 10.1039/d3cp03760f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
The dissociative double photoionization of isoxazole molecules has been investigated experimentally and theoretically. The experiment has been carried out in the 27.5-36 eV photon energy range using vacuum ultraviolet (VUV) synchrotron radiation excitation combined with ion time-of-flight (TOF) spectrometry and photoelectron-photoion-photoion coincidence (PEPIPICO) technique. Five well-resolved two-body dissociation channels have been identified in the isoxazole's coincidence maps, and their appearance energies have been determined. The coincidence yield curves of these dissociation channels have been obtained in the photon energy ranges from their appearance energies up to 36 eV. The double photoionization of isoxazole produces a C3H3NO2+ transient dication, which decomposes into fragments differing from previously reported photofragmentation products of isoxazole. We have found no evidence of pathways leading to the C3H2NO+, HCN+, C2H2O+, C3HN+, or C2H2+ fragments or their neutral counterparts that have been observed in previous neutral photodissociation and single photoionization studies. Instead, the dissociation of isoxazole after the ejection of two electrons is bond-selective and is governed by two reactions, HCO+ + H2CCN+ and H2CO+ + HCCN+, whose appearance energies are 28.6 (±0.3) and 29.4 (±0.3) eV, respectively. A third dissociation channel turns out to be a variant of the most intense channel (HCO+ + H2CCN+), where one of the fragment ions contains a heavy isotope. Two minor dissociation channels occurring at higher energies, CO+ + CH3CN+ and CN+ + H3CCO+, are also identified. The density functional and ab initio quantum chemical calculations have been performed to elucidate the dissociative charge-separating mechanisms and determine the energies of the observed photoproducts. The present work unravels hitherto unexplored photodissociation mechanisms of isoxazole and thus provides deeper insight into the photophysics of five-membered heterocyclic molecules containing two heteroatoms.
Collapse
Affiliation(s)
- Tomasz J Wasowicz
- Division of Complex Systems Spectroscopy, Institute of Physics and Computer Science, Gdansk University of Technology, ul. G. Narutowicza 11/12, 80-233 Gdansk, Poland.
| | - Iwona Dąbkowska
- Department of Analytical Chemistry, Faculty of Chemistry, University of Gdansk, ul. Wita Stwosza 63, 80-308 Gdansk, Poland
| | | | | |
Collapse
|
35
|
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
|
36
|
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
|
37
|
Dillon AD, Gieseking RLM. Convergence of Time-Derivative Nonadiabatic Couplings in Plane-Wave DFT Calculations. J Phys Chem A 2023; 127:9612-9620. [PMID: 37924298 DOI: 10.1021/acs.jpca.3c04858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2023]
Abstract
Accurate prediction of charge carrier relaxation rates is essential to design molecules and materials with the desired photochemical properties for applications like photocatalysis and solar energy conversion. Nonadiabatic molecular dynamics allows one to simulate the relaxation process of excited charge carriers. Plane-wave density functional theory (DFT) calculations make the time-derivative nonadiabatic couplings (TNACs) simple to compute because the basis is independent of the atomic positions. However, the effect of the kinetic energy cutoff for the plane-wave basis on the accuracy of the dynamics has not been studied. Here, we examine the effect of the kinetic energy cutoff on the TNACs and decay time scales for the prototypical model system of tetracene. These calculations show that the choice of kinetic energy cutoff can change the relaxation time by up to 30%. The relaxation times of states that have small TNACs to other states or are far from degenerate are more sensitive to the kinetic energy cutoff than those of states with large TNACs or near degeneracies. A kinetic energy cutoff of 60 Ry is sufficient for all states to reach semiquantitative agreement (absolute error <10%) with the decay times of our 110 Ry reference data, and a cutoff of 80 Ry is required for all states to reach quantitative agreement (absolute error <2%).
Collapse
Affiliation(s)
- Alva D Dillon
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| | - Rebecca L M Gieseking
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| |
Collapse
|
38
|
Wu GY, Zhu HM, Li H, Zhang K, Zhang X, Yan D, Zhang XD, Lin L, Lu Z. The impact of aggregation of AIE and ACQ moiety-integrating material on the excited state dynamics. RSC Adv 2023; 13:33911-33917. [PMID: 38020029 PMCID: PMC10658659 DOI: 10.1039/d3ra06359c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 10/21/2023] [Indexed: 12/01/2023] Open
Abstract
The investigation of the properties of aggregate materials is highly interesting because the process of aggregation can result in the disappearance of original properties and the emergence of new ones. Here, a novel fluorescent material (TPEIP), which synergistically combines aggregation-induced emission (AIE) and aggregation caused quenching (ACQ) moieties, was first synthesized by the cyclization reaction of 2,3-diamino-phenazine with 4-tetraphenylenthenealdehyde. We controlled the degree of aggregation of TPEIP to shed light on the impact of the aggregation on the excited state dynamics. TPEIP aggregation realized control over the Intersystem Crossing (ISC) rates and, in turn, the suppression of triplet excited states in MeOH, EtOH or via the simple addition of water to TPEIP solutions in DMSO. From global target analysis, the time scale was 966.2 ps for ISC for TPEIP in DMSO, but it was 860 ps in the case of TPEIP solutions featuring 5% water. The dynamics of TPEIP excited states undergo significant changes as the degree of aggregation increases. Notably, the lifetime of singlet excited states decreases, and there was a gradual diminishment in triplet states.
Collapse
Affiliation(s)
- Gui-Yuan Wu
- Anhui Province Key Laboratory of Optoelectronic Material Science and Technology, School of Physics and Electronic Information, Anhui Normal University Wuhu 241002 China
| | - Hui-Min Zhu
- Anhui Province Key Laboratory of Optoelectronic Material Science and Technology, School of Physics and Electronic Information, Anhui Normal University Wuhu 241002 China
| | - Hao Li
- Anhui Province Key Laboratory of Optoelectronic Material Science and Technology, School of Physics and Electronic Information, Anhui Normal University Wuhu 241002 China
| | - Kai Zhang
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University Jinan 250014 China
| | - Xianyi Zhang
- Anhui Province Key Laboratory of Optoelectronic Material Science and Technology, School of Physics and Electronic Information, Anhui Normal University Wuhu 241002 China
| | - Dong Yan
- Key Laboratory of Functional Molecular Solids, Ministry of Education, College of Chemistry and Materials Science, Anhui Normal University Wuhu 241002 China
| | - Xiu-Du Zhang
- Key Laboratory of Functional Molecular Solids, Ministry of Education, College of Chemistry and Materials Science, Anhui Normal University Wuhu 241002 China
| | - Lili Lin
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University Jinan 250014 China
| | - Zhou Lu
- Anhui Province Key Laboratory of Optoelectronic Material Science and Technology, School of Physics and Electronic Information, Anhui Normal University Wuhu 241002 China
| |
Collapse
|
39
|
Li G, Shi Z, Guo X, Wang L. What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching. J Phys Chem Lett 2023; 14:9855-9863. [PMID: 37890155 DOI: 10.1021/acs.jpclett.3c02690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/29/2023]
Abstract
When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.
Collapse
Affiliation(s)
- Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
40
|
Giannini S, Di Virgilio L, Bardini M, Hausch J, Geuchies JJ, Zheng W, Volpi M, Elsner J, Broch K, Geerts YH, Schreiber F, Schweicher G, Wang HI, Blumberger J, Bonn M, Beljonne D. Transiently delocalized states enhance hole mobility in organic molecular semiconductors. NATURE MATERIALS 2023; 22:1361-1369. [PMID: 37709929 DOI: 10.1038/s41563-023-01664-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 08/14/2023] [Indexed: 09/16/2023]
Abstract
Evidence shows that charge carriers in organic semiconductors self-localize because of dynamic disorder. Nevertheless, some organic semiconductors feature reduced mobility at increasing temperature, a hallmark for delocalized band transport. Here we present the temperature-dependent mobility in two record-mobility organic semiconductors: dinaphtho[2,3-b:2',3'-f]thieno[3,2-b]-thiophene (DNTT) and its alkylated derivative, C8-DNTT-C8. By combining terahertz photoconductivity measurements with atomistic non-adiabatic molecular dynamics simulations, we show that while both crystals display a power-law decrease of the mobility (μ) with temperature (T) following μ ∝ T -n, the exponent n differs substantially. Modelling reveals that the differences between the two chemically similar semiconductors can be traced to the delocalization of the different states that are thermally accessible by charge carriers, which in turn depends on their specific electronic band structure. The emerging picture is that of holes surfing on a dynamic manifold of vibrationally dressed extended states with a temperature-dependent mobility that provides a sensitive fingerprint for the underlying density of states.
Collapse
Affiliation(s)
- Samuele Giannini
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium.
| | | | - Marco Bardini
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium
| | - Julian Hausch
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | | | - Wenhao Zheng
- Max Planck Institute for Polymer Research, Mainz, Germany
| | - Martina Volpi
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Jan Elsner
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, UK
| | - Katharina Broch
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | - Yves H Geerts
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
- International Solvay Institutes for Physics and Chemistry, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Frank Schreiber
- Institut für Angewandte Physik, Universität Tübingen, Tübingen, Germany
| | - Guillaume Schweicher
- Laboratoire de Chimie des Polymères, Faculté des Sciences, Université Libre de Bruxelles (ULB), Bruxelles, Belgium
| | - Hai I Wang
- Max Planck Institute for Polymer Research, Mainz, Germany.
- Nanophotonics, Debye Institute for Nanomaterials Science, Utrecht University, Utrecht, The Netherlands.
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London, UK
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Mainz, Germany.
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, University of Mons, Mons, Belgium.
| |
Collapse
|
41
|
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
|
42
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
43
|
Polonius S, Zhuravel O, Bachmair B, Mai S. LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution. J Chem Theory Comput 2023; 19:7171-7186. [PMID: 37788824 PMCID: PMC10601485 DOI: 10.1021/acs.jctc.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Indexed: 10/05/2023]
Abstract
We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - Oleksandra Zhuravel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
44
|
Menger MFSJ, Ou Q, Shao Y, Faraji S, Subotnik JE, Cofer-Shabica DV. Nature of Hops, Coordinates, and Detailed Balance for Nonadiabatic Simulations in the Condensed Phase. J Phys Chem A 2023; 127:8427-8436. [PMID: 37782887 DOI: 10.1021/acs.jpca.3c03546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Photoinduced processes play a crucial role in a multitude of important molecular phenomena. Accurately modeling these processes in an environment other than a vacuum requires a detailed description of the electronic states involved as well as how energy flows are coupled to the surroundings. Nonadiabatic effects must also be included in order to describe the exchange of energy between electronic and nuclear degrees of freedom correctly. In this work, we revisit the ring-opening reaction 1,3-cylohexadiene (CHD) in a solvent environment. Using our newly developed Interface for Non-Adiabatic Quantum mechanics/molecular mechanics in Solvent (INAQS) we trace the evolution of the reaction via hybrid quantum mechanics/molecular mechanics (QM/MM) surface hopping with a focus on the solvent's participation in the nonadiabatic relaxation process and the long-time approach to equilibrium. We explicitly include the MM solvent contribution to the nonadiabatic coupling vector─enabling an accurate approach to equilibrium at long times─and find that in highly multidimensional systems gradients can have little or nothing to do with the nonadiabatic couplings.
Collapse
Affiliation(s)
- Maximilian F S J Menger
- Theoretische Chemie, Physikalisch-Chemisches Institut, University Heidelberg, INF 229, 69120 Heidelberg, Germany
- Zernike Institute for Advanced Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Qi Ou
- AI for Science Institute, Beijing 100084, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Joseph E Subotnik
- Chemistry, University of Pennsylvania, 231 S. 34 Street, Cret Wing 141D, Philadelphia, Pennsylvania 19104-6243, United States
| | - D Vale Cofer-Shabica
- Chemistry, University of Pennsylvania, 231 S. 34 Street, Cret Wing 141D, Philadelphia, Pennsylvania 19104-6243, United States
| |
Collapse
|
45
|
Zhao X, Merritt ICD, Lei R, Shu Y, Jacquemin D, Zhang L, Xu X, Vacher M, Truhlar DG. Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation. J Chem Theory Comput 2023; 19:6577-6588. [PMID: 37772732 DOI: 10.1021/acs.jctc.3c00813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/30/2023]
Abstract
Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.
Collapse
Affiliation(s)
- Xiaorui Zhao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, P. R. China
| | | | - Ruiqing Lei
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris 75005, France
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xuefei Xu
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, P. R. China
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
46
|
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
|
47
|
Mandal A, Taylor MA, Weight BM, Koessler ER, Li X, Huo P. Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics. Chem Rev 2023; 123:9786-9879. [PMID: 37552606 PMCID: PMC10450711 DOI: 10.1021/acs.chemrev.2c00855] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 08/10/2023]
Abstract
When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | - Michael A.D. Taylor
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| | - Braden M. Weight
- Department
of Physics and Astronomy, University of
Rochester, Rochester, New York 14627, United
States
| | - Eric R. Koessler
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Xinyang Li
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
48
|
Khaldi O, Ferhi H, Larbi T, Jomni F, Ben Younes R. Dielectric response of high- κ hafnium oxide under finite electric field: nonlinearities from ab initio and experimental points of view. Phys Chem Chem Phys 2023; 25:22310-22318. [PMID: 37578192 DOI: 10.1039/d3cp01552a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Herein, we report on the dielectric-voltage nonlinearities under a constant electric field in metal insulator metal (MIM) capacitor-based hafnium oxide (HfO2) with respect to the frequency range. Via the Schottky emission mechanism obtained from the current-voltage characteristic (I-V), we calculated the optical dielectric constant εr,opt for different external DC bias values. The extracted εr,opt revealed a quadratic dependence on the applied external field. This confirmed that such dependence is a common feature of high-κ oxides in the low and high frequency ranges. The results were correlated with the ab initio calculations using the finite field (FF) method as implemented in the CRYSTAL 17 code. Good agreement between the results from the FF method, I-V plots, as well as the UV-visible spectrometry is observed. To assess any change in the dielectric constant upon the application of an external electric field, several parameters such as exchange-correlation functional, basis sets (BSs), as well as supercell expansion factor (N) were tested. The corresponding parameters have a great influence on the macroscopic electron density and voltage along the field direction and thus on the optical response. For N > 2 and rich basis sets, the hybrid functional B3LYP revealed good agreement with the experimental results as compared to other Hamiltonian's forms such as LDA, PW-GGA and HF.
Collapse
Affiliation(s)
- Othmen Khaldi
- University of Tunis El Manar, Laboratory of Materials, Organization and Properties, LMOP(LR99ES17), 2092, Tunis, Tunisia.
| | - Hanen Ferhi
- University of Gafsa, Laboratory of Technology, Energy and Innovative Materials, Department of Physics, Faculty of Sciences of Gafsa, 2112, Gafsa, Tunisia
| | - Tarek Larbi
- University of Tunis El Manar, Laboratory of Nanomaterials, Nanotechnology and Energy, 2092, Tunis, Tunisia
| | - Fethi Jomni
- University of Tunis El Manar, Laboratory of Materials, Organization and Properties, LMOP(LR99ES17), 2092, Tunis, Tunisia.
| | - Rached Ben Younes
- University of Gafsa, Laboratory of Technology, Energy and Innovative Materials, Department of Physics, Faculty of Sciences of Gafsa, 2112, Gafsa, Tunisia
| |
Collapse
|
49
|
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
|
50
|
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
|