1
|
Toldo JM, Mattos RS, Pinheiro M, Mukherjee S, Barbatti M. Recommendations for Velocity Adjustment in Surface Hopping. J Chem Theory Comput 2024; 20:614-624. [PMID: 38207213 DOI: 10.1021/acs.jctc.3c01159] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.
Collapse
Affiliation(s)
- Josene M Toldo
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | - Rafael S Mattos
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | - Max Pinheiro
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | | | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| |
Collapse
|
2
|
Gómez S, Spinlove E, Worth G. Benchmarking non-adiabatic quantum dynamics using the molecular Tully models. Phys Chem Chem Phys 2024; 26:1829-1844. [PMID: 38170796 DOI: 10.1039/d3cp03964a] [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
On-the-fly non-adiabatic dynamics methods are becoming more important as tools to characterise the time evolution of a system after absorbing light. These methods, which calculate quantities such as state energies, gradients and interstate couplings at every time step, circumvent the requirement for pre-computed potential energy surfaces. There are a number of different algorithms used, the most common being Tully Surface Hopping (TSH), but all are approximate solutions to the time-dependent Schrödinger equation and benchmarking is required to understand their accuracy and performance. For this, a common set of systems and observables are required to compare them. In this work, we validate the on-the-fly direct dynamics variational multi-configuration Gaussian (DD-vMCG) method using three molecular systems recently suggested by Ibele and Curchod as molecular versions of the Tully model systems used to test one-dimensional non-adiabatic behaviour [Ibele et al., Phys. Chem. Chem. Phys. 2020, 22, 15183-15196]. Parametrised linear vibronic potential energy surfaces for each of the systems were also tested and compared to on-the-fly results. The molecules, which we term the Ibele-Curchod models, are ethene, DMABN and fulvene and the authors used them to test and compare several versions of the Ab Initio Multiple Spawning (AIMS) method alongside TSH. The three systems present different deactivation pathways after excitation to their ππ* bright states. When comparing DD-vMCG to AIMS and TSH, we obtain crucial differences in some cases, for which an explanation is provided by the classical nature and the chosen initial conditions of the TSH simulations.
Collapse
Affiliation(s)
- Sandra Gómez
- Departamento de Química Física, Universidad de Salamanca, 37008, Spain
| | - Eryn Spinlove
- Faculty of Science and Engineering, Theoretical Chemistry - Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
- Department of Chemistry, University College London, 20 Gordon St, London WC1H 0AJ, UK.
| | - Graham Worth
- Department of Chemistry, University College London, 20 Gordon St, London WC1H 0AJ, UK.
| |
Collapse
|
3
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 08/31/2023]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
4
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.1] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 03/12/2024]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
5
|
Swan E, Platts K, Blencowe A. An overview of the cycloaddition chemistry of fulvenes and emerging applications. Beilstein J Org Chem 2019; 15:2113-2132. [PMID: 31579091 PMCID: PMC6753682 DOI: 10.3762/bjoc.15.209] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 08/21/2019] [Indexed: 01/21/2023] Open
Abstract
The unusual electronic properties and unique reactivity of fulvenes have interested researchers for over a century. The propensity to form dipolar structures at relatively low temperatures and to participate as various components in cycloaddition reactions, often highly selectively, makes them ideal for the synthesis of complex polycyclic carbon scaffolds. As a result, fulvene cycloaddition chemistry has been employed extensively for the synthesis of natural products. More recently, fulvene cycloaddition chemistry has also found application to other areas including materials chemistry and dynamic combinatorial chemistry. This highlight article discusses the unusual properties of fulvenes and their varied cycloaddition chemistry, focussing on applications in organic and natural synthesis, dynamic combinatorial chemistry and materials chemistry, including dynamers, hydrogels and charge transfer complexes. Tables providing comprehensive directories of fulvene cycloaddition chemistry are provided, including fulvene intramolecular and intermolecular cycloadditions complete with reactant partners and their resulting cyclic adducts, which provide a useful reference source for synthetic chemists working with fulvenes and complex polycyclic scaffolds.
Collapse
Affiliation(s)
- Ellen Swan
- School of Pharmacy and Medical Sciences, The University of South Australia, Adelaide, South Australia 5000, Australia
| | - Kirsten Platts
- School of Pharmacy and Medical Sciences, The University of South Australia, Adelaide, South Australia 5000, Australia
| | - Anton Blencowe
- School of Pharmacy and Medical Sciences, The University of South Australia, Adelaide, South Australia 5000, Australia.,Future Industries Institute, The University of South Australia, Mawson Lakes, South Australia 5095, Australia
| |
Collapse
|
6
|
Tahchieva DN, Bakowies D, Ramakrishnan R, von Lilienfeld OA. Torsional Potentials of Glyoxal, Oxalyl Halides, and Their Thiocarbonyl Derivatives: Challenges for Popular Density Functional Approximations. J Chem Theory Comput 2018; 14:4806-4817. [PMID: 30011363 DOI: 10.1021/acs.jctc.8b00174] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides, and their thiocarbonyl derivatives. HF and 18 functionals of varying complexity, from local density to range-separated hybrid approximations and double-hybrid, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, most functionals fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on σ-hole based intermolecular binding in dimers and crystals.
Collapse
Affiliation(s)
- Diana N Tahchieva
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| | - Dirk Bakowies
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| | - Raghunathan Ramakrishnan
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| |
Collapse
|
7
|
Preethalayam P, Krishnan KS, Thulasi S, Chand SS, Joseph J, Nair V, Jaroschik F, Radhakrishnan KV. Recent Advances in the Chemistry of Pentafulvenes. Chem Rev 2017; 117:3930-3989. [DOI: 10.1021/acs.chemrev.6b00210] [Citation(s) in RCA: 87] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Preethanuj Preethalayam
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| | - K. Syam Krishnan
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| | - Sreeja Thulasi
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| | - S. Sarath Chand
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| | - Jomy Joseph
- Institut de Chimie Moléculaire de Reims CNRS (UMR 7312) and Université de Reims Champagne-Ardenne, Reims 51687 Cedex 2, France
| | - Vijay Nair
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| | - Florian Jaroschik
- Institut de Chimie Moléculaire de Reims CNRS (UMR 7312) and Université de Reims Champagne-Ardenne, Reims 51687 Cedex 2, France
| | - K. V. Radhakrishnan
- Organic
Chemistry Section, National Institute for Interdisciplinary Science and Technology (CSIR-NIIST), Trivandrum 695 019, Kerala, India
| |
Collapse
|
8
|
Alborzpour JP, Tew DP, Habershon S. Efficient and accurate evaluation of potential energy matrix elements for quantum dynamics using Gaussian process regression. J Chem Phys 2016; 145:174112. [DOI: 10.1063/1.4964902] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan P. Alborzpour
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
9
|
Tuna D, Sobolewski AL, Domcke W. Conical-Intersection Topographies Suggest That Ribose Exhibits Enhanced UV Photostability. J Phys Chem B 2016; 120:10729-10735. [DOI: 10.1021/acs.jpcb.6b09048] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Deniz Tuna
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | | | - Wolfgang Domcke
- Department
of Chemistry, Technische Universität München, 85747 Garching, Germany
| |
Collapse
|
10
|
Li H, Migani A, Blancafort L, Li Q, Li Z. Early events in the photochemistry of 5-diazo Meldrum's acid: formation of a product manifold in C–N bound and pre-dissociated intersection seam regions. Phys Chem Chem Phys 2016; 18:30785-30793. [DOI: 10.1039/c6cp06290c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Photolysis of 5-diazo Meldrum's acid yields three different primary products in a sub-ps scale. High-level computations show that this is due to the presence of an extended seam of intersection associated to different reaction paths.
Collapse
Affiliation(s)
- Huijing Li
- Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- 100081 Beijing
| | - Annapaola Migani
- Catalan Institute of Nanoscience and Nanotechnology (ICN2)
- CSIC and The Barcelona Institute of Science and Technology
- Campus UAB
- 08193 Barcelona
- Spain
| | - Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química
- Universitat de Girona
- Facultat de Ciències
- 17003 Girona
- Spain
| | - Quansong Li
- Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- 100081 Beijing
| | - Zesheng Li
- Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- 100081 Beijing
| |
Collapse
|
11
|
Blancafort L. Photochemistry and photophysics at extended seams of conical intersection. Chemphyschem 2014; 15:3166-81. [PMID: 25157686 DOI: 10.1002/cphc.201402359] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2014] [Indexed: 11/07/2022]
Abstract
The role of extended seams of conical intersection in excited-state mechanisms is reviewed. Seams are crossings of the potential energy surface in many dimensions where the decay from the excited to the ground state can occur, and the extended seam is composed of different segments lying along a reaction coordinate. Every segment is associated with a different primary photoproduct, which gives rise to competing pathways. This idea is first illustrated for fulvene and ethylene, and then it is used to explain more complex cases such as the dependence of the isomerisation of retinal chromophore isomers on the protein environment, the dependence of the efficiency of the azobenzene photochemical switch on the wavelength of irradiation and the direction of the isomerisation, and the coexistence of different mechanisms in the photo-induced Wolff rearrangement of diazonaphthoquinone. The role of extended seams in the photophysics of the DNA nucleobases and the relationship between two-state seams and three-state crossings is also discussed. As an outlook, the design of optical control strategies based on the passage of the excited molecule through the seam is considered, and it is shown how the excited-state lifetime of fulvene can be modulated by shaping the energy of the seam.
Collapse
Affiliation(s)
- Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona (Spain).
| |
Collapse
|
12
|
Jenkins S, Blancafort L, Kirk SR, Bearpark MJ. The response of the electronic structure to electronic excitation and double bond torsion in fulvene: a combined QTAIM, stress tensor and MO perspective. Phys Chem Chem Phys 2014; 16:7115-26. [PMID: 24618666 DOI: 10.1039/c4cp00003j] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
New insights into the double bond isomerization of fulvene in the ground and excited electronic states are provided by newly developed QTAIM and stress tensor tools. The S0 and S1 states follow the 'biradical' torsion model, but the double bond is stiffer in the S0 state; by contrast, the S2 state follows the 'zwitterionic' torsion. Differences are explained in terms of the ellipticity and bond critical point (BCP) stiffness for both QTAIM and the stress tensor. Overall, the wave-function based analysis is found to be in agreement with the work of Bonačić-Koutecký and Michl that the bond-twisted species can have biradical or zwitterionic character, depending on the state. Using QTAIM and the stress tensor a new understanding of bond torsion is revealed; the electronic charge density around the twisted bond is found not to rotate in concert with the nuclei of the rotated -CH2 methylene group. The ability to visualize how the bond stiffness varies between individual electronic states and how this correlates with the QTAIM and stress tensor bond stiffness is highlighted. In addition, the most and least preferred morphologies of bond-path torsion are visualized. Briefly we discuss the prospects for using this new QTAIM and stress tensor analysis for excited state chemistry.
Collapse
Affiliation(s)
- Samantha Jenkins
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China) and Key Laboratory of Resource Fine-Processing and Advanced Materials of Hunan Province, College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha Hunan 410081, China.
| | | | | | | |
Collapse
|
13
|
Joubert-Doriol L, Lasorne B, Lauvergnat D, Meyer HD, Gatti F. A generalised vibronic-coupling Hamiltonian model for benzopyran. J Chem Phys 2014; 140:044301. [DOI: 10.1063/1.4861226] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Sousa C, de Graaf C, Rudavskyi A, Broer R, Tatchen J, Etinski M, Marian CM. Ultrafast Deactivation Mechanism of the Excited Singlet in the Light‐Induced Spin Crossover of [Fe(2,2′‐bipyridine)
3
]
2+. Chemistry 2013; 19:17541-51. [DOI: 10.1002/chem.201302992] [Citation(s) in RCA: 126] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2013] [Indexed: 11/09/2022]
Affiliation(s)
- Carmen Sousa
- Departament de Química Física and Institut de Química, Teòrica i Computacional, Universitat de Barcelona, C/Martí i Franquès 1, 08028 Barcelona (Spain)
| | - Coen de Graaf
- Departament de Química Física i Inorgànica, Universitat Rovira i Virgili, Marcel‐lí Domingo s/n, 43007 Tarragona (Spain), Fax: (+34)‐977559563
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, 08010 Barcelona (Spain)
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen (The Netherlands)
| | - Andrii Rudavskyi
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen (The Netherlands)
| | - Ria Broer
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen (The Netherlands)
| | - Jörg Tatchen
- Institute of Theoretical and Computational Chemistry, Heinrich‐Heine University Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf (Germany)
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski Trg 12–16, 11158 Belgrade (Serbia)
| | - Christel M. Marian
- Institute of Theoretical and Computational Chemistry, Heinrich‐Heine University Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf (Germany)
| |
Collapse
|
15
|
Ruiz-Barragan S, Blancafort L. Photophysics of fulvene under the non-resonant stark effect. Shaping the conical intersection seam. Faraday Discuss 2013; 163:497-512; discussion 513-43. [PMID: 24020219 DOI: 10.1039/c3fd20155d] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We introduce a mechanistic strategy to control the excited state lifetime of fulvene based on shaping the topography of an extended seam of intersection with the non-resonant dynamic Stark effect. Fulvene has a very short excited state lifetime due to an energetically accessible seam of intersection which lies along the methylene torsion coordinate, and the initial decay occurs at the seam segment around the planar conical intersection structure. We have followed a three-step approach to simulate the control. First, we have calculated the effect of a non-resonant electric field on the potential energy surface at the ab initio level, including the field in a self-consistent way. The relative energy of the planar segment of the seam is increased by the non-resonant field. In the second step we simulate the control carrying out MCTDH quantum dynamics propagations under a static non-resonant field to derive the main control mechanisms. At moderately intense fields (epsilon < or = 0.03 a.u.) the decay is faster as compared to the field free case because the vibrational overlap between the excited and ground state vibrational functions is increased. However, at more intense fields (epsilon = 0.04 a.u.) the planar conical intersection is energetically inaccessible and the decay occurs at a slower time scale, at the segment of the seam with more twisted geometries. In the third step, the control over the dynamics is exerted with a non-resonant dynamic field. The acceleration of the decay due to the improved vibrational overlap does not occur, but the decay can be made slower with a dynamic field of 0.08 a.u. The results show the viability of our approach to control the photophysics shaping the topology of the conical intersection seam, and they prove that the extended nature of the seam is crucial for simulating and understanding the control.
Collapse
Affiliation(s)
- Sergi Ruiz-Barragan
- Institut de Quimica Computacional and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona, Spain
| | | |
Collapse
|
16
|
Zimmermann T, Vaníček J. Evaluation of the importance of spin-orbit couplings in the nonadiabatic quantum dynamics with quantum fidelity and with its efficient "on-the-fly" ab initio semiclassical approximation. J Chem Phys 2013; 137:22A516. [PMID: 23249053 DOI: 10.1063/1.4738878] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose to measure the importance of spin-orbit couplings (SOCs) in the nonadiabatic molecular quantum dynamics rigorously with quantum fidelity. To make the criterion practical, quantum fidelity is estimated efficiently with the multiple-surface dephasing representation (MSDR). The MSDR is a semiclassical method that includes nuclear quantum effects through interference of mixed quantum-classical trajectories without the need for the Hessian of potential energy surfaces. Two variants of the MSDR are studied, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics. The fidelity criterion and MSDR are first tested on one-dimensional model systems amenable to numerically exact quantum dynamics. Then, the MSDR is combined with "on-the-fly" computed electronic structure to measure the importance of SOCs and nonadiabatic couplings in the photoisomerization dynamics of CH(2)NH(2)(+) considering 20 electronic states and in the collision of F + H(2) considering six electronic states.
Collapse
Affiliation(s)
- Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
17
|
Mendive-Tapia D, Lasorne B, Worth GA, Robb MA, Bearpark MJ. Towards converging non-adiabatic direct dynamics calculations using frozen-width variational Gaussian product basis functions. J Chem Phys 2012; 137:22A548. [DOI: 10.1063/1.4765087] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Chartrand D, Castro Ruiz CA, Hanan GS. Diimine Triscarbonyl Re(I) of Isomeric Pyridyl-fulvene Ligands: an Electrochemical, Spectroscopic, and Computational Investigation. Inorg Chem 2012; 51:12738-47. [DOI: 10.1021/ic301559s] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniel Chartrand
- Department of Chemistry, Université de Montréal, Montréal,
Québec H3T 1J4, Canada
| | | | - Garry S. Hanan
- Department of Chemistry, Université de Montréal, Montréal,
Québec H3T 1J4, Canada
| |
Collapse
|
19
|
Sethi A, Keshavamurthy S. Driven coupled Morse oscillators: visualizing the phase space and characterizing the transport. Mol Phys 2012. [DOI: 10.1080/00268976.2012.667166] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
20
|
Li Q, Migani A, Blancafort L. Wave Packet Dynamics at an Extended Seam of Conical Intersection: Mechanism of the Light-Induced Wolff Rearrangement. J Phys Chem Lett 2012; 3:1056-1061. [PMID: 26286570 DOI: 10.1021/jz300235z] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Quantum dynamics calculations on a model surface based on CASPT2//CASSCF calculations are carried out to probe the traversal of a wave packet through an extended seam of conical intersection during the light-induced Wolff rearrangement of diazonaphtoquinone. The reaction is applied in the fabrication of integrated circuits. It consists of nitrogen elimination and ring rearrangement to yield a ketene. After excitation, the wave packet relaxes and reaches the extended seam. A fraction of the wave packet decays to the ground state at a region of the seam connected to a carbene intermediate, while the remaining part decays at a region leading to the ketene. The passage of the wave packet through the extended seam explains the competition between concerted ketene formation and a stepwise mechanism involving a carbene. The two primary photoproducts are formed in the first 100 fs of the simulation, in agreement with recent ultrafast spectroscopy measurements.
Collapse
Affiliation(s)
- Quansong Li
- Institut de Química Computacional and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona, Spain
| | - Annapaola Migani
- Institut de Química Computacional and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona, Spain
| | - Lluís Blancafort
- Institut de Química Computacional and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona, Spain
| |
Collapse
|