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
|
Filatov M, Mironov V, Kraka E. Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores. J Comput Chem 2024; 45:1033-1045. [PMID: 38216513 DOI: 10.1002/jcc.27304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 11/23/2023] [Accepted: 12/23/2023] [Indexed: 01/14/2024]
Abstract
The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S 1 /S 0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.
Collapse
Affiliation(s)
- Michael Filatov
- Center for Multidimensional Carbon Materials, Institute for Basic Science (IBS), Ulsan, Republic of Korea
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| | | | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| |
Collapse
|
3
|
Jaiswal VK, Montorsi F, Aleotti F, Segatta F, Keefer D, Mukamel S, Nenov A, Conti I, Garavelli M. Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level. J Chem Phys 2024; 160:164316. [PMID: 38686819 DOI: 10.1063/5.0203624] [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: 04/10/2024] [Indexed: 05/02/2024] Open
Abstract
We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.
Collapse
Affiliation(s)
- V K Jaiswal
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Montorsi
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Aleotti
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F Segatta
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Daniel Keefer
- Molecular Spectroscopy Department, Max Planck Institute for Polymer Research, Mainz 55128, Germany
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697, USA
| | - A Nenov
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - I Conti
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - M Garavelli
- Dipartimento di Chimica industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| |
Collapse
|
4
|
Vörös D, Mai S. Role of Ultrafast Internal Conversion and Intersystem Crossing in the Nonadiabatic Relaxation Dynamics of ortho-Nitrobenzaldehyde. J Phys Chem A 2023. [PMID: 37405967 PMCID: PMC10364085 DOI: 10.1021/acs.jpca.3c02899] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
ortho-Nitrobenzaldehyde (oNBA) is a well-known photoactivated acid and a prototypical photolabile nitro-aromatic compound. Despite extensive investigations, the ultrafast relaxation dynamics of oNBA is still not properly understood, especially concerning the role of the triplet states. In this work, we provide an in-depth picture of this dynamics by combining single- and multireference electronic structure methods with potential energy surface exploration and nonadiabatic dynamics simulations using the Surface Hopping including ARbitary Couplings (SHARC) approach. Our results reveal that the initial decay from the bright ππ* state to the S1 minimum is barrierless. It involves three changes in electronic structure from ππ* (ring) to nπ* (nitro group), to nπ* (aldehyde group), and then to another nπ* (nitro group). The decay of the ππ* takes 60-80 fs and can be tracked with time-resolved luminescence spectroscopy, where we predict for the first time a short-lived coherence of the luminescence energy with a 25 fs period. Intersystem crossing can occur already during the S4 → S1 deactivation cascade but also from S1, with a time constant of about 2.4 ps and such that first a triplet ππ* state localized on the nitro group is populated. The triplet population first evolves into an nπ* and then quickly undergoes hydrogen transfer to form a biradical intermediate, from where the ketene is eventually produced. The majority of the excited population decays from S1 through two conical intersections of equal utilization, a previously unreported one involving a scissoring motion of the nitro group that leads back to the oNBA ground state and the one involving hydrogen transfer that leads to the ketene intermediate.
Collapse
Affiliation(s)
- Dóra Vörös
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Vienna Doctoral School in Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| |
Collapse
|
5
|
Reiner M, Bachmair B, Tiefenbacher MX, Mai S, González L, Marquetand P, Dellago C. Nonadiabatic Forward Flux Sampling for Excited-State Rare Events. J Chem Theory Comput 2023; 19:1657-1671. [PMID: 36856706 PMCID: PMC10061683 DOI: 10.1021/acs.jctc.2c01088] [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: 11/01/2022] [Indexed: 03/02/2023]
Abstract
We present a rare event sampling scheme applicable to coupled electronic excited states. In particular, we extend the forward flux sampling (FFS) method for rare event sampling to a nonadiabatic version (NAFFS) that uses the trajectory surface hopping (TSH) method for nonadiabatic dynamics. NAFFS is applied to two dynamically relevant excited-state models that feature an avoided crossing and a conical intersection with tunable parameters. We investigate how nonadiabatic couplings, temperature, and reaction barriers affect transition rate constants in regimes that cannot be otherwise obtained with plain, traditional TSH. The comparison with reference brute-force TSH simulations for limiting cases of rareness shows that NAFFS can be several orders of magnitude cheaper than conventional TSH and thus represents a conceptually novel tool to extend excited-state dynamics to time scales that are able to capture rare nonadiabatic events.
Collapse
Affiliation(s)
- Madlen
Maria Reiner
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Maximilian Xaver Tiefenbacher
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Leticia González
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Christoph Dellago
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Faculty
of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
6
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, France
| |
Collapse
|
7
|
Merritt ICD, Jacquemin D, Vacher M. Nonadiabatic Coupling in Trajectory Surface Hopping: How Approximations Impact Excited-State Reaction Dynamics. J Chem Theory Comput 2023; 19:1827-1842. [PMID: 36897995 DOI: 10.1021/acs.jctc.2c00968] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
Abstract
Photochemical reactions are widely modeled using the popular trajectory surface hopping (TSH) method, an affordable mixed quantum-classical approximation to the full quantum dynamics of the system. TSH is able to account for nonadiabatic effects using an ensemble of trajectories, which are propagated on a single potential energy surface at a time and which can hop from one electronic state to another. The occurrences and locations of these hops are typically determined using the nonadiabatic coupling between electronic states, which can be assessed in a number of ways. In this work, we benchmark the impact of some approximations to the coupling term on the TSH dynamics for several typical isomerization and ring-opening reactions. We have identified that two of the schemes tested, the popular local diabatization scheme and a scheme based on biorthonormal wave function overlap implemented in the OpenMOLCAS code as part of this work, reproduce at a much reduced cost the dynamics obtained using the explicitly calculated nonadiabatic coupling vectors. The other two schemes tested can give different results, and in some cases, even entirely incorrect dynamics. Of these two, the scheme based on configuration interaction vectors gives unpredictable failures, while the other scheme based on the Baeck-An approximation systematically overestimates hopping to the ground state as compared to the reference approaches.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
8
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis. J Chem Theory Comput 2022; 18:5837-5855. [DOI: 10.1021/acs.jctc.2c00702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| |
Collapse
|
9
|
Park W, Filatov (Gulak) M, Sadiq S, Gerasimov I, Lee S, Joo T, Choi CH. A Plausible Mechanism of Uracil Photohydration Involves an Unusual Intermediate. J Phys Chem Lett 2022; 13:7072-7080. [PMID: 35900137 PMCID: PMC9358713 DOI: 10.1021/acs.jpclett.2c01694] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2022] [Accepted: 07/11/2022] [Indexed: 05/28/2023]
Abstract
It is well-known that photolysis of pyrimidine nucleobases, such as uracil, in an aqueous environment results in the formation of hydrate as one of the main products. Although several hypotheses regarding photohydration have been proposed in the past, e.g., the zwitterionic and "hot" ground-state mechanisms, its detailed mechanism remains elusive. Here, theoretical nonadiabatic simulations of the uracil photodynamics reveal the formation of a highly energetic but kinetically stable intermediate that features a half-chair puckered pyrimidine ring and a strongly twisted intracyclic double bond. The existence and the kinetic stability of the intermediate are confirmed by a variety of computational chemistry methods. According to the simulations, the unusual intermediate is mainly formed almost immediately (∼50-200 fs) upon photoabsorption and survives long enough to engage in a hydration reaction with a neighboring water. A plausible mechanism of uracil photohydration is proposed on the basis of the modeling of nucleophilic insertion of water into the twisted double bond of the intermediate.
Collapse
Affiliation(s)
- Woojin Park
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | | | - Saima Sadiq
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Igor Gerasimov
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Taiha Joo
- Department
of Chemistry, Pohang University of Science
and Technology (POSTECH), Pohang 37673, South Korea
| | - Cheol Ho Choi
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
10
|
Valverde D, Mai S, Canuto S, Borin AC, González L. Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water. JACS AU 2022; 2:1699-1711. [PMID: 35911449 PMCID: PMC9327080 DOI: 10.1021/jacsau.2c00250] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S2 state that converts internally to the dark S1 state, from which the population is transferred to the triplet T2 state via intersystem crossing and finally to the lowest T1 state. This S2 → S1 → T2 → T1 deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T1 lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T1/S0 crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ1 = 131 fs and τ2 = 191 fs) that are in excellent agreement with the experimentally reported value (τexp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.
Collapse
Affiliation(s)
- Danillo Valverde
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| | - Sylvio Canuto
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Antonio Carlos Borin
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| |
Collapse
|
11
|
Akimov AV. Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain. J Phys Chem Lett 2021; 12:12119-12128. [PMID: 34913701 DOI: 10.1021/acs.jpclett.1c03823] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.
Collapse
Affiliation(s)
- Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
12
|
Lin K, Peng J, Gu FL, Lan Z. Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network. J Phys Chem Lett 2021; 12:10225-10234. [PMID: 34647736 DOI: 10.1021/acs.jpclett.1c02672] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
13
|
M J C, Nag P, Vennapusa SR. Surface hopping dynamics reveal ultrafast triplet generation promoted by S 1-T 2-T 1 spin-vibronic coupling in 2-mercaptobenzothiazole. Phys Chem Chem Phys 2021; 23:20183-20192. [PMID: 34473155 DOI: 10.1039/d1cp02587b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the T1 formation upon populating the optically "bright" S2 in 2-mercaptobenzothiazole to interpret the underlying relaxation pathways associated with the experimental decay constants reported by D. Koyama and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 2016, 18, 26224-26235. Energetics, electronic populations and geometries of various stationary points of low-lying electronic states are computed using the semi-classical ab initio surface hopping dynamics simulations. Estimated decay constants of S2-S1 internal conversion (IC) and S1-T2 intersystem crossing (ISC) are in excellent agreement with the experiment. The observed ultrafast ISC is analyzed based on the S1-T2-T1 spin-vibronic coupling mechanism. In contrast to the previous assignment of 6 ps to the T2-T1 IC, our findings enable us to attribute this decay constant to the combined events of T2-T1 IC followed by relaxation of vibrationally hot T1.
Collapse
Affiliation(s)
- Chithra M J
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| | - Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Maruthamala PO, Vithura, Thiruvananthapuram 695551, Kerala, India.
| |
Collapse
|
14
|
Zobel JP, González L. The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes. JACS AU 2021; 1:1116-1140. [PMID: 34467353 PMCID: PMC8397362 DOI: 10.1021/jacsau.1c00252] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Indexed: 05/15/2023]
Abstract
This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
| |
Collapse
|
15
|
Zobel JP, Knoll T, González L. Ultrafast and long-time excited state kinetics of an NIR-emissive vanadium(iii) complex II. Elucidating triplet-to-singlet excited-state dynamics. Chem Sci 2021; 12:10791-10801. [PMID: 34476060 PMCID: PMC8372553 DOI: 10.1039/d1sc02149d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2021] [Accepted: 05/25/2021] [Indexed: 11/21/2022] Open
Abstract
We report the non-adiabatic dynamics of VIIICl3(ddpd), a complex based on the Earth-abundant first-row transition metal vanadium with a d2 electronic configuration which is able to emit phosphorescence in solution in the near-infrared spectral region. Trajectory surface-hopping dynamics based on linear vibronic coupling potentials obtained with CASSCF provide molecular-level insights into the intersystem crossing from triplet to singlet metal-centered states. While the majority of the singlet population undergoes back-intersystem crossing to the triplet manifold, 1-2% remains stable during the 10 ps simulation time, enabling the phosphorescence described in Dorn et al. Chem. Sci., 2021, DOI: 10.1039/D1SC02137K. Competing with intersystem crossing, two different relaxation channels via internal conversion through the triplet manifold occur. The nuclear motion that drives the dynamics through the different electronic states corresponds mainly to the increase of all metal-ligand bond distances as well as the decrease of the angles of trans-coordinated ligand atoms. Both motions lead to a decrease in the ligand-field splitting, which stabilizes the interconfigurational excited states populated during the dynamics. Analysis of the electronic character of the states reveals that increasing and stabilizing the singlet population, which in turn can result in enhanced phosphorescence, could be accomplished by further increasing the ligand-field strength.
Collapse
Affiliation(s)
- J Patrick Zobel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna Währingerstr. 19 1090 Vienna Austria
| | - Thomas Knoll
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna Währingerstr. 19 1090 Vienna Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna Währingerstr. 19 1090 Vienna Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna Währingerstr. 19 1090 Vienna Austria
| |
Collapse
|
16
|
Zhang RM, Xu X, Truhlar DG. Energy Dependence of Ensemble-Averaged Energy Transfer Moments and Its Effect on Competing Decomposition Reactions. J Phys Chem A 2021; 125:6303-6313. [PMID: 34232653 DOI: 10.1021/acs.jpca.1c03845] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We carried out a direct dynamics study on the internal-energy dependence of the ensemble-averaged energy transfer moments of the isobutyl radical in collisions with N2 bath gas. We find a linear dependence of the downward moment ⟨ΔEd⟩ and the root-mean-square moment ⟨ΔE2⟩ on the initial internal energy, but the upward moment ⟨ΔEu⟩ is found to be independent of the molecule's internal energy. We improved the exponential-down relaxation model by including a linear dependence of ⟨ΔEd⟩ on the initial energy, and we used the improved treatment in the 1D master equation for isobutyl radical decomposition reactions and for a model of competitive reactions with a larger difference in barrier heights. We calculated phenomenological rate constants and branching ratios from chemically significant eigenmodes of the master equation and showed that the energy dependence of ⟨ΔEd⟩ has a greater influence on channels with higher barriers in competitive reactions. Rate constants and branching ratios from master equation calculations indicate that for a given temperature and pressure, there is a constant ⟨ΔEd⟩ that can reproduce results obtained with an E-dependent ⟨ΔEd⟩. But a constant ⟨ΔEd⟩ cannot do this for all temperatures and pressures, with larger differences when the barriers for the competing channels differ more. We conclude that when the branching ratio of competitive reactions is sensitive to pressure, including the energy dependence of ⟨ΔEd⟩ in master equation simulations can make a significant difference in the results.
Collapse
Affiliation(s)
- Rui Ming Zhang
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Xuefei Xu
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
17
|
Park W, Lee S, Huix-Rotllant M, Filatov M, Choi CH. Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine. J Phys Chem Lett 2021; 12:4339-4346. [PMID: 33929858 DOI: 10.1021/acs.jpclett.1c00712] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S1-trapping, S2-trapping, and S1&S2-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S2-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S1-trapping mechanism with two lifetimes, τ1 = 30 ± 1 fs and τ2 = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S2 (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S2→S1 internal conversion and resides around the minimum on the S1 (nOπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S2→S1 internal conversion is mediated by fast bond length alternation distortion, the subsequent S1→S0 occurs through several conical intersections, involving a slow puckering motion.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
18
|
Milovanović B, Novak J, Etinski M, Domcke W, Došlić N. Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters. Phys Chem Chem Phys 2021; 23:2594-2604. [PMID: 33475644 DOI: 10.1039/d0cp05618a] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.
Collapse
Affiliation(s)
| | - Jurica Novak
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia. and Higher Medical and Biological School, South Ural State University, Chelyabinsk, Russia
| | - Mihajlo Etinski
- University of Belgrade, Faculty of Physical Chemistry, Belgrade, Serbia
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching, Germany
| | - Nađa Došlić
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia.
| |
Collapse
|
19
|
Filatov M, Lee S, Nakata H, Choi CH. Signatures of Conical Intersection Dynamics in the Time-Resolved Photoelectron Spectrum of Furan: Theoretical Modeling with an Ensemble Density Functional Theory Method. Int J Mol Sci 2021; 22:4276. [PMID: 33924097 PMCID: PMC8074317 DOI: 10.3390/ijms22084276] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 04/13/2021] [Accepted: 04/15/2021] [Indexed: 12/13/2022] Open
Abstract
The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck-Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA;
| | - Hiroya Nakata
- R & D Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan;
| | - Cheol-Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, Korea
| |
Collapse
|
20
|
Parker SM, Schiltz CJ. Surface hopping with cumulative probabilities: Even sampling and improved reproducibility. J Chem Phys 2020; 153:174109. [DOI: 10.1063/5.0024372] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Shane M. Parker
- Department of Chemistry, Case Western Reserve University, 10800 Euclid Ave., Cleveland, Ohio 44106, USA
| | - Colin J. Schiltz
- Department of Chemistry, Case Western Reserve University, 10800 Euclid Ave., Cleveland, Ohio 44106, USA
| |
Collapse
|
21
|
Talotta F, Boggio‐Pasqua M, González L. Early Relaxation Dynamics in the Photoswitchable Complex trans-[RuCl(NO)(py) 4 ] 2. Chemistry 2020; 26:11522-11528. [PMID: 32281169 PMCID: PMC7539916 DOI: 10.1002/chem.202000507] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 03/16/2020] [Indexed: 11/22/2022]
Abstract
The design of photoswitchable transition metal complexes with tailored properties is one of the most important challenges in chemistry. Studies explaining the underlying mechanisms are, however, scarce. Herein, the early relaxation dynamics towards NO photoisomerization in trans-[RuCl(NO)(py)4 ]2+ is elucidated by means of non-adiabatic dynamics, which provided time-resolved information and branching ratios. Three deactivation mechanisms (I, II, III) in the ratio 3:2:4 were identified. Pathways I and III involve ultrafast intersystem crossing and internal conversion, whereas pathway II involves only internal conversion.
Collapse
Affiliation(s)
- Francesco Talotta
- lnstitute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
- Laboratoire de Chimie et Physique Quantiques, UMR 5626, lRSAMCCNRS et Université Toulouse 3118 route de Narbonne31062ToulouseFrance
- Present Address: Laboratoire de Chimie PhysiqueUMR 8000 CNRS/University Paris-SudUniversity Paris-Saclay, 91403 Orsay, andlnstitut de Sciences Moleculaires d'Orsay, UMR8214 CNRS/University Paris-Sud, University Paris-Saclay91403OrsayFrance
| | - Martial Boggio‐Pasqua
- Laboratoire de Chimie et Physique Quantiques, UMR 5626, lRSAMCCNRS et Université Toulouse 3118 route de Narbonne31062ToulouseFrance
| | - Leticia González
- lnstitute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
- Vienna Research Platform on Accelerating Photoreaction DiscoveryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
22
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
23
|
Curchod BFE, Glover WJ, Martínez TJ. SSAIMS-Stochastic-Selection Ab Initio Multiple Spawning for Efficient Nonadiabatic Molecular Dynamics. J Phys Chem A 2020; 124:6133-6143. [PMID: 32580552 DOI: 10.1021/acs.jpca.0c04113] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Ab initio multiple spawning provides a powerful and accurate way of describing the excited-state dynamics of molecular systems, whose strength resides in the proper description of coherence effects during nonadiabatic processes thanks to the coupling of trajectory basis functions. However, the simultaneous propagation of a large number of trajectory basis functions can be numerically inconvenient. We propose here an elegant and simple solution to this issue, which consists of (i) detecting uncoupled groups of coupled trajectory basis functions and (ii) selecting stochastically one of these groups to continue the ab initio multiple spawning dynamics. We show that this procedure can reproduce the results of full ab initio multiple spawning dynamics in cases where the uncoupled groups of trajectory basis functions stay uncoupled throughout the dynamics (which is often the case in high-dimensional problems). We present and discuss the aforementioned idea in detail and provide simple numerical applications on indole, ethylene, and protonated formaldimine, highlighting the potential of stochastic-selection ab initio multiple spawning.
Collapse
Affiliation(s)
- Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - William J Glover
- NYU Shanghai, 1555 Century Ave., Shanghai 200122, 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
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
24
|
Gong Q, Li Y, Liu F. Trajectory surface-hopping dynamics of light-driven imine switches. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
25
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
26
|
Zobel JP, González L. Nonadiabatic Dynamics Simulation Predict Intersystem Crossing in Nitroaromatic Molecules on a Picosecond Time Scale. CHEMPHOTOCHEM 2019; 3:833-845. [PMID: 31681833 PMCID: PMC6813632 DOI: 10.1002/cptc.201900108] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2019] [Revised: 05/06/2019] [Indexed: 12/19/2022]
Abstract
Previous time-resolved spectroscopic experiments and static quantum-chemical calculations attributed nitronaphthalene derivatives one of the fastest time scales for intersystem crossing within organic molecules, reaching the 100 fs mark. Nonadiabatic dynamics simulations on three nitronaphthalene derivatives challenge this view, showing that the experimentally observed ∼100 fs process corresponds to internal conversion in the singlet manifolds. Intersystem crossing, instead, takes place on a longer time scale of ∼1 ps. The dynamics simulations further reveal that the spin transitions occur via two distinct pathways with different contribution for the three systems, which are determined by electronic factors and the torsion of the nitro group. This study, therefore, indicates that the existence of sub-picosecond intersystem crossing in other nitroaromatic molecules should be questioned.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Division of Theoretical Chemistry, KemicentrumLund UniversityP.O. Box 124SE-221 00LundSweden
| | - Leticia González
- Institute of Theoretical ChemistryUniversity of ViennaWähringer Straße 17A-1090ViennaAustria
| |
Collapse
|
27
|
Wang H, Wen K, You X, Mao Q, Luo KH, Pilling MJ, Robertson SH. Energy transfer in intermolecular collisions of polycyclic aromatic hydrocarbons with bath gases He and Ar. J Chem Phys 2019; 151:044301. [PMID: 31370521 DOI: 10.1063/1.5094104] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Classical trajectory simulations of intermolecular collisions were performed for a series of polycyclic aromatic hydrocarbons interacting with the bath gases helium and argon for bath gas temperature from 300 to 2500 K. The phase-space average energy transferred per deactivating collision, ⟨∆Edown⟩, was obtained. The Buckingham pairwise intermolecular potentials were validated against high-level quantum chemistry calculations and used in the simulations. The reactive force-field was used to describe intramolecular potentials. The dependence of ⟨∆Edown⟩ on initial vibrational energy is discussed. A canonical sampling method was compared with a microcanonical sampling method for selecting initial vibrational energy at high bath gas temperatures. Uncertainties introduced by the initial angular momentum distribution were identified. The dependence of the collisional energy transfer parameters on the type of bath gas and the molecular structure of polycyclic aromatic hydrocarbons was examined.
Collapse
Affiliation(s)
- Hongmiao Wang
- Center for Combustion Energy, Tsinghua University, Beijing 100084, China
| | - Kaicheng Wen
- Center for Combustion Energy, Tsinghua University, Beijing 100084, China
| | - Xiaoqing You
- Center for Combustion Energy, Tsinghua University, Beijing 100084, China
| | - Qian Mao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, China
| | - Kai Hong Luo
- Department of Mechanical Engineering, University College London, Torrington Place, London WC1E 7JE, United Kingdom
| | - Michael J Pilling
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Struan H Robertson
- Dassault Systèmes, BIOVIA, 334, Cambridge Science Park, Cambridge CB4 0WN, United Kingdom
| |
Collapse
|
28
|
Plasser F, Mai S, Fumanal M, Gindensperger E, Daniel C, González L. Strong Influence of Decoherence Corrections and Momentum Rescaling in Surface Hopping Dynamics of Transition Metal Complexes. J Chem Theory Comput 2019; 15:5031-5045. [DOI: 10.1021/acs.jctc.9b00525] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Felix Plasser
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K
| | - Sebastian Mai
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| | - Maria Fumanal
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Etienne Gindensperger
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Chantal Daniel
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Leticia González
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|
29
|
Polyak I, Hutton L, Crespo-Otero R, Barbatti M, Knowles PJ. Ultrafast Photoinduced Dynamics of 1,3-Cyclohexadiene Using XMS-CASPT2 Surface Hopping. J Chem Theory Comput 2019; 15:3929-3940. [DOI: 10.1021/acs.jctc.9b00396] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Iakov Polyak
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, United Kingdom
| | - Lewis Hutton
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, United Kingdom
| | - Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | | - Peter J. Knowles
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
30
|
Mai S, Wolf AP, González L. Curious Case of 2-Selenouracil: Efficient Population of Triplet States and Yet Photostable. J Chem Theory Comput 2019; 15:3730-3742. [DOI: 10.1021/acs.jctc.9b00208] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Anna-Patricia Wolf
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| |
Collapse
|
31
|
Parker SM, Roy S, Furche F. Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation. Phys Chem Chem Phys 2019; 21:18999-19010. [DOI: 10.1039/c9cp03127h] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.
Collapse
Affiliation(s)
- Shane M. Parker
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Saswata Roy
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Filipp Furche
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| |
Collapse
|
32
|
Filatov M, Paolino M, Min SK, Choi CH. Design and photoisomerization dynamics of a new family of synthetic 2-stroke light driven molecular rotary motors. Chem Commun (Camb) 2019; 55:5247-5250. [DOI: 10.1039/c9cc01955c] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Synthetic 2-stroke light driven molecular rotary motors with ultrafast function and high quantum efficiency.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry
- Kyungpook National University
- Daegu 702-701
- South Korea
| | - Marco Paolino
- Dipartimento di Biotecnologie
- Chimica e Farmacia (Dipartimento di Eccellenza 2018-2022)
- Università di Siena
- 53100 Siena
- Italy
| | - Seung Kyu Min
- Department of Chemistry
- School of Natural Sciences
- Ulsan National Institute of Science and Technology
- Ulsan 44919
- South Korea
| | - Cheol Ho Choi
- Department of Chemistry
- Kyungpook National University
- Daegu 702-701
- South Korea
| |
Collapse
|
33
|
Zobel JP, Nogueira JJ, González L. Finite-temperature Wigner phase-space sampling and temperature effects on the excited-state dynamics of 2-nitronaphthalene. Phys Chem Chem Phys 2019; 21:13906-13915. [DOI: 10.1039/c8cp03273d] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The concept of finite temperature Wigner phase-space sampling allowing the population of vibrationally excited states is introduced and employed to study temperature effects on the absorption spectrum of 2-nitronaphtalene (2NN) and its relaxation dynamics.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute of Theoretical Chemistry
- University of Vienna
- A-1090 Vienna
- Austria
| | - Juan J. Nogueira
- Institute of Theoretical Chemistry
- University of Vienna
- A-1090 Vienna
- Austria
| | - Leticia González
- Institute of Theoretical Chemistry
- University of Vienna
- A-1090 Vienna
- Austria
| |
Collapse
|
34
|
Plasser F, Gómez S, Menger MFSJ, Mai S, González L. Highly efficient surface hopping dynamics using a linear vibronic coupling model. Phys Chem Chem Phys 2018; 21:57-69. [PMID: 30306987 DOI: 10.1039/c8cp05662e] [Citation(s) in RCA: 72] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.
Collapse
Affiliation(s)
- Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, UK.
| | | | | | | | | |
Collapse
|
35
|
Mai S, Marquetand P, González L. Nonadiabatic dynamics: The SHARC approach. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2018; 8:e1370. [PMID: 30450129 PMCID: PMC6220962 DOI: 10.1002/wcms.1370] [Citation(s) in RCA: 224] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2018] [Revised: 04/05/2018] [Accepted: 04/06/2018] [Indexed: 12/12/2022]
Abstract
We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.
Collapse
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry University of Vienna Vienna Austria
| |
Collapse
|
36
|
Mai S, Gattuso H, Monari A, González L. Novel Molecular-Dynamics-Based Protocols for Phase Space Sampling in Complex Systems. Front Chem 2018; 6:495. [PMID: 30386775 PMCID: PMC6199692 DOI: 10.3389/fchem.2018.00495] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Accepted: 09/27/2018] [Indexed: 11/13/2022] Open
Abstract
The adequate exploration of the phase space of a chromophore is a fundamental necessity for the simulation of their optical and photophysical properties, taking into account the effects of vibrational motion and, most importantly, the coupling with a (non-homogeneous) molecular environment. A representative set of conformational snapshots around the Franck-Condon region is also required to perform non-adiabatic molecular dynamics, for instance in the framework of surface hopping. Indeed, in the latter case one needs to prepare a set of initial conditions providing a meaningful and complete statistical base for the subsequent trajectory propagation. In this contribution, we propose two new protocols for molecular dynamics-based phase space sampling, called "local temperature adjustment" and "individual QM/MM-based relaxation." These protocols are intended for situations in which the popular Wigner distribution sampling procedure is not applicable-as it is the case when anharmonic or nonlinear vibrations are present-and where regular molecular dynamics sampling might suffer from an inaccurate distribution of internal energy or from inaccurate force fields. The new protocols are applied to the case of phase space sampling of [Re(CO)3(Im)(Phen)]+ (im, imidazole; phen, phenanthroline) in aqueous solution, showing the advantages and limitations of regular Wigner and molecular dynamics sampling as well as the strengths of the new protocols.
Collapse
Affiliation(s)
- Sebastian Mai
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Hugo Gattuso
- Université de Lorraine and CNRS, LPTC UMR 7019, Nancy, France
| | - Antonio Monari
- Université de Lorraine and CNRS, LPTC UMR 7019, Nancy, France
| | - Leticia González
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| |
Collapse
|
37
|
Filatov M, Min SK, Kim KS. Non-adiabatic dynamics of ring opening in cyclohexa-1,3-diene described by an ensemble density-functional theory method. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1519200] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan, Korea
| | - Seung Kyu Min
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan, Korea
| | - Kwang S. Kim
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan, Korea
| |
Collapse
|
38
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
39
|
Zobel JP, Nogueira JJ, González L. Mechanism of Ultrafast Intersystem Crossing in 2-Nitronaphthalene. Chemistry 2018; 24:5379-5387. [PMID: 29377370 PMCID: PMC5947663 DOI: 10.1002/chem.201705854] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2017] [Indexed: 01/10/2023]
Abstract
Nitronaphthalene derivatives efficiently populate their electronically excited triplet states upon photoexcitation through ultrafast intersystem crossing (ISC). Despite having been studied extensively by time-resolved spectroscopy, the reasons behind their ultrafast ISC remain unknown. Herein, we present the first ab initio nonadiabatic molecular dynamics study of a nitronaphthalene derivative, 2-nitronaphthalene, including singlet and triplet states. We find that there are two distinct ISC reaction pathways involving different electronic states at distinct nuclear configurations. The high ISC efficiency is explained by the very small electronic and nuclear alterations that the chromophore needs to undergo during the singlet-triplet transition in the dominating ISC pathway after initial dynamics in the singlet manifold. The insights gained in this work are expected to shed new light on the photochemistry of other nitro polycyclic aromatic hydrocarbons that exhibit ultrafast intersystem crossing.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute of Theoretical Chemistry, Faculty of ChemistryUniversity of ViennaWähringer Straße 171090ViennaAustria
| | - Juan J. Nogueira
- Institute of Theoretical Chemistry, Faculty of ChemistryUniversity of ViennaWähringer Straße 171090ViennaAustria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of ChemistryUniversity of ViennaWähringer Straße 171090ViennaAustria
| |
Collapse
|
40
|
Atkins AJ, González L. Trajectory Surface-Hopping Dynamics Including Intersystem Crossing in [Ru(bpy) 3] 2. J Phys Chem Lett 2017; 8:3840-3845. [PMID: 28766339 DOI: 10.1021/acs.jpclett.7b01479] [Citation(s) in RCA: 90] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Surface-hopping dynamics coupled to linear response TDDFT and explicit nonadiabatic and spin-orbit couplings have been used to model the ultrafast intersystem crossing (ISC) dynamics in [Ru(bpy)3]2+. Simulations using an ensemble of trajectories starting from the singlet metal-to-ligand charge transfer (1MLCT) band show that the manifold of 3MLCT triplet states is first populated from high-lying singlet states within 26 ± 3 fs. ISC competes with an intricate internal conversion relaxation process within the singlet manifold to the lowest singlet state. Normal-mode analysis and principal component analysis, combined with further dynamical simulations where the nuclei are frozen, unequivocally demonstrate that it is not only the high density of states and the large spin-orbit couplings of the system that promote ISC. Instead, geometrical relaxation involving the nitrogen atoms is required to allow for state mixing and efficient triplet population transfer.
Collapse
Affiliation(s)
- Andrew J Atkins
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna , Währinger Straße 17, A-1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna , Währinger Straße 17, A-1090 Vienna, Austria
| |
Collapse
|
41
|
Peccati F, Mai S, González L. Insights into the deactivation of 5-bromouracil after ultraviolet excitation. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:20160202. [PMID: 28320905 PMCID: PMC5360901 DOI: 10.1098/rsta.2016.0202] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 11/03/2016] [Indexed: 06/06/2023]
Abstract
5-Bromouracil is a nucleobase analogue that can replace thymine in DNA strands and acts as a strong radiosensitizer, with potential applications in molecular biology and cancer therapy. Here, the deactivation of 5-bromouracil after ultraviolet irradiation is investigated in the singlet and triplet manifold by accurate quantum chemistry calculations and non-adiabatic dynamics simulations. It is found that, after irradiation to the bright ππ* state, three main relaxation pathways are, in principle, possible: relaxation back to the ground state, intersystem crossing (ISC) and C-Br photodissociation. Based on accurate MS-CASPT2 optimizations, we propose that ground-state relaxation should be the predominant deactivation pathway in the gas phase. We then employ different electronic structure methods to assess their suitability to carry out excited-state dynamics simulations. MRCIS (multi-reference configuration interaction including single excitations) was used in surface hopping simulations to compute the ultrafast ISC dynamics, which mostly involves the 1nOπ* and 3ππ* states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.
Collapse
Affiliation(s)
- Francesca Peccati
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| |
Collapse
|
42
|
Carlos Borin A, Mai S, Marquetand P, González L. Ab initio molecular dynamics relaxation and intersystem crossing mechanisms of 5-azacytosine. Phys Chem Chem Phys 2017; 19:5888-5894. [DOI: 10.1039/c6cp07919a] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Nonadiabatic SHARC dynamics simulations reveal the molecular deformations involved in the photodeactivation pathways of 5-azacytosine.
Collapse
Affiliation(s)
- Antonio Carlos Borin
- Department of Fundamental Chemistry
- Institute of Chemistry
- University of São Paulo
- NAP-PhotoTech the USP Consortium for Photochemical Technology
- São Paulo
| | - Sebastian Mai
- Institute of Theoretical Chemistry
- Faculty of Chemistry
- University of Vienna
- 1090 Vienna
- Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry
- Faculty of Chemistry
- University of Vienna
- 1090 Vienna
- Austria
| | - Leticia González
- Institute of Theoretical Chemistry
- Faculty of Chemistry
- University of Vienna
- 1090 Vienna
- Austria
| |
Collapse
|
43
|
|
44
|
Mai S, Marquetand P, González L. Intersystem Crossing Pathways in the Noncanonical Nucleobase 2-Thiouracil: A Time-Dependent Picture. J Phys Chem Lett 2016; 7:1978-83. [PMID: 27167106 PMCID: PMC4893732 DOI: 10.1021/acs.jpclett.6b00616] [Citation(s) in RCA: 109] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
The deactivation mechanism after ultraviolet irradiation of 2-thiouracil has been investigated using nonadiabatic dynamics simulations at the MS-CASPT2 level of theory. It is found that after excitation the S2 quickly relaxes to S1, and from there intersystem crossing takes place to both T2 and T1 with a time constant of 400 fs and a triplet yield above 80%, in very good agreement with recent femtosecond experiments in solution. Both indirect S1 → T2 → T1 and direct S1 → T1 pathways contribute to intersystem crossing, with the former being predominant. The results contribute to the understanding of how some noncanonical nucleobases respond to harmful ultraviolet light, which could be relevant for prospective photochemotherapeutic applications.
Collapse
|
45
|
Nangia S, Jasper AW, Miller TF, Truhlar DG. Erratum: “Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method” [J. Chem. Phys. 120, 3586 (2004)]. J Chem Phys 2016; 144:139901. [DOI: 10.1063/1.4944694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Shikha Nangia
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Ahren W. Jasper
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Thomas F. Miller
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
46
|
Lingerfelt DB, Williams-Young DB, Petrone A, Li X. Direct ab Initio (Meta-)Surface-Hopping Dynamics. J Chem Theory Comput 2016; 12:935-45. [DOI: 10.1021/acs.jctc.5b00697] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- David B. Lingerfelt
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | | - Alessio Petrone
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
47
|
Soloviov M, Meuwly M. CO-dynamics in the active site of cytochrome c oxidase. J Chem Phys 2015; 140:145101. [PMID: 24735320 DOI: 10.1063/1.4870264] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The transfer of CO from heme a3 to the Cu(B) site in Cytochrome c oxidase (CcO) after photolysis is studied using molecular dynamics simulations using an explicitly reactive, parametrized potential energy surface based on density functional theory calculations. After photodissociation from the heme-Fe, the CO ligand rebinds to the Cu(B) site on the sub-picosecond time scale. Depending on the simulation protocol the characteristic time ranges from 260 fs to 380 fs which compares with an estimated 450 fs from experiment based on the analysis of the spectral changes as a function of time delay after the photodissociating pulse. Following photoexcitation ≈90% of the ligands are found to rebind to either the Cu(B) (major component, 85%) or the heme-Fe (minor component, 2%) whereas about 10% remain in an unbound state. The infrared spectra of unbound CO in the active site is broad and featureless and no appreciable shift relative to gas-phase CO is found, which is in contrast to the situation in myoglobin. These observations explain why experimentally, unbound CO in the binuclear site of CcO has not been found as yet.
Collapse
Affiliation(s)
- Maksym Soloviov
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
48
|
DeSalvo SC, Liu Y, Choudhary GS, Ren D, Nangia S, Sureshkumar R. Signaling factor interactions with polysaccharide aggregates of bacterial biofilms. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:1958-1966. [PMID: 25607640 DOI: 10.1021/la504721b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Biofilms are surface-attached colonies of bacteria embedded in an extracellular polymeric substance (EPS). Inside the eukaryotic hosts, bacterial biofilms interact with the host cells through signaling factors (SFs). These signaling processes play important roles in the interaction between bacteria and host cells and the outcome of infections and symbiosis. However, how host immune factors diffuse through biofilms is not well understood. Here, we describe synergistic molecular dynamics and experimental approaches for studying the translocation of signaling factors through polysaccharide chain aggregates present in the extracellular matrix of bacterial biofilms. The effect of polysaccharide chain degradation on the energetics of SF-EPS interactions was examined by simulating an EPS consisting of various polysaccharide chain lengths. It is shown that the SF stabilization energy, defined as the average potential of mean force difference between the environments outside and within the matrix, increases linearly with decreasing chain length. This effect has been explained based on the changes in the polysaccharide configurations around the SF. Specifically, shorter chains are packed tightly around the SF, promoting favorable SF-EPS interactions, while longer chains are packed loosely resulting in screening of interactions with neighboring chains. We further investigated the translocation of SFs through the host cell membrane using molecular dynamics simulations. Further, simulations predict the existence of energy barriers greater than 1000 kJ mol(-1) associated with the translocation of the signaling factors necrosis factor-alpha (TNF-α) and granulocyte macrophage colony stimulating factor (GM-CSF) across the lipid bilayer. The agreement of computational and experimental findings motivates future computational studies using a more detailed description of the EPS aimed at understanding the role of the extracellular matrix on biofilm drug resistance.
Collapse
Affiliation(s)
- Stephen C DeSalvo
- Department of Biomedical and Chemical Engineering, §Department of Civil and Environmental Engineering, ∥Department of Biology, ⊥Department of Physics, ‡Syracuse Biomaterials Institute, Syracuse University , Syracuse, New York 13244, United States
| | | | | | | | | | | |
Collapse
|
49
|
Malhado JP, Bearpark MJ, Hynes JT. Non-adiabatic dynamics close to conical intersections and the surface hopping perspective. Front Chem 2014; 2:97. [PMID: 25485263 PMCID: PMC4240175 DOI: 10.3389/fchem.2014.00097] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2014] [Accepted: 10/16/2014] [Indexed: 11/13/2022] Open
Abstract
Conical intersections play a major role in the current understanding of electronic de-excitation in polyatomic molecules, and thus in the description of photochemistry and photophysics of molecular systems. This article reviews aspects of the basic theory underlying the description of non-adiabatic transitions at conical intersections, with particular emphasis on the important case when the dynamics of the nuclei are treated classically. Within this classical nuclear motion framework, the main aspects of the surface hopping methodology in the conical intersection context are presented. The emerging picture from this treatment is that of electronic transitions around conical intersections dominated by the interplay of the nuclear velocity and the derivative non-adiabatic coupling vector field.
Collapse
Affiliation(s)
| | | | - James T Hynes
- Department of Chemistry and Biochemistry, University of Colorado, Boulder Boulder, CO, USA ; Département de Chimie, École Normale Supérieur, UMR ENS-CNRS-UPMC 8640 Paris, France
| |
Collapse
|
50
|
Jasper AW, Kamarchik E, Miller JA, Klippenstein SJ. First-principles binary diffusion coefficients for H, H2, and four normal alkanes + N2. J Chem Phys 2014; 141:124313. [DOI: 10.1063/1.4896368] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ahren W. Jasper
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, USA
| | - Eugene Kamarchik
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, USA
| | - James A. Miller
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Stephen J. Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|