Cyclohexadiene Revisited: A Time-Resolved ... - ACS Publications

Mar 28, 2016 - Cyclohexadiene Revisited: A Time-Resolved Photoelectron. Spectroscopy ...... (65) Hoover, W. Canonical Dynamics: Equilibrium Phase-Spac...
0 downloads 0 Views 4MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Cyclohexadiene Revisited – a Time Resolved Photoelectron Spectroscopy and Ab Initio Study Oliver Schalk, Ting Geng, Travis Thompson, Noel K Baluyot, Richard D. Thomas, Enrico Tapavicza, and Tony Hansson J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.5b10928 • Publication Date (Web): 28 Mar 2016 Downloaded from http://pubs.acs.org on March 28, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Cyclohexadiene Revisited – a Time Resolved Photoelectron Spectroscopy and Ab Initio Study Oliver Schalk,

1,*

1

Enrico Tapavicza, 1

2

2

Ting Geng, Travis Thompson, Noel Baluyot, Richard D. Thomas, 2,*

Tony Hansson,

1

1

Department of Chemical Physics, AlbaNova University Centre, Stockholm University, Roslagstullsbacken 21, 109 61 Stockholm, Sweden

2

Long Beach, Department of Chemistry and Biochemistry, California State University, Long Beach, CA, United States of America

* O. Schalk. Electronic address: [email protected]; E. Tapavicza. Electronic address: [email protected]

Abstract We have reinvestigated the excited state dynamics of cyclohexa-1,3-diene (CHD) with timeresolved photoelectron spectroscopy and fewest switches surface hopping molecular dynamics based on linear response time-dependent density functional theory after excitation to the lowest lying ππ* (1B)-state. The combination of both theory and experiment revealed several new results: First, the dynamics progress on one single excited state surface. After an incubation time of (35 ± 10) fs on the excited state, the dynamics proceed to the ground state in an additional (60 ± 10) fs, either via a conrotatory ring-opening to hexatriene or back to the CHD ground state. Moreover, ring-opening predominantly occurs when the wavepacket crosses the region of strong non-adiabatic coupling with a positive velocity in the bond alternation coordinate. After 100 fs, trajectories remaining in the excited state must return to the CHD ground state. This extra time-delay induces a revival of the photoelectron signal and is an experimental confirmation of the previously formulated model of two parallel reaction channels with distinct time constants. Finally, our simulations suggest that, after the initially formed cis-Z-cis HT rotamer, the trans-Z-trans isomer is formed, before the thermodynamical equilibrium of three possible rotamers is reached after 1 ps.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 42

I. INTRODUCTION Cyclohexa-1,3-diene (CHD) is one of the most thoroughly studied molecules in time resolved photoexcitation experiments1-16 and ab initio computations9, 15, 17-31 and is the subject of two recent reviews.32, 33 The large interest originally stems from CHD being the prototype for a pericyclic reaction used to explain the Woodward-Hoffmann rules.34 The earliest models of the photodynamics were extended later via the van der Lugt / Oosterhoff model35 and beyond.28 In addition, the photoinduced dynamics of many larger systems are based on the cyclohexadiene dynamophore36, fulgide,

38, 39

37

such as the photoswitches diarylethylene and indolyl-

and the ring-opening dynamics of provitamin D, which is an intermediate in the

synthesis of vitamin D.28-30 The general view of the photo-induced dynamics upon excitation to the lowest lying ππ*-state (1B-symmetry in the Franck-Condon region) is as follows: After excitation, the dynamics proceed to a first conical intersection (CoIn) with a state that is both dark to one photon absorption and has 2A-symmetry in the Franck-Condon region. In the gas phase, a time constant of (60 ± 15) fs for the initial excited state to reach this CoIn was extracted7, 10, 14, 40 (slightly varying between the experiments, the shortest value being 30 fs16 ). After switching to the 2A-surface, the seam with the ground state is reached within 60 - 80 fs after excitation7, 10, 14, 16, 40

where the molecules can either return to the ground state without bond breaking or

undergo a conrotatory ringopening reaction to give Z-hexatriene. Quantum yields for the ringopening process were calculated to 60 % 15 and 62 % 30 in the gas phase, while experimental studies on CHD solvated in pentane gave a value of 41 %.41 In the gas phase, early scattering experiments provided no evidence for return to the CHD ground state,8, 42 but later work using multiphoton probe provided a maximal hexatriene yield of 73%

43

while a recent study with

90 nm probe pulses showed the build-up and recovery traces of hexatriene and CHD, respectively,14 from which the authors extracted a hexatriene quantum yield as low as 30 %. On the ground state surface, the molecules rearrange from cis-Z-cis-hexatriene (cZc) to cis-Ztrans-hexatriene (cZt) and finally to trans-Z-trans-hexatriene (tZt) which was studied in the liquid phase and takes place within a few hundred of femtoseconds.4, 5, 6 In the gas phase, a recent study suggests that several rotamers can be formed immediately; however, only the first 200 fs of the dynamics were investigated.15 Despite the apparently well studied reaction path, there are still unresolved questions concerning the photoinduced dynamics of CHD. In a recent theoretical treatment, Tapavicza et al. found no evidence for two excited states in the dynamics of CHD.30, 44 Instead, they

ACS Paragon Plus Environment

Page 3 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

report that the character of the excited state changed along the reaction coordinate, resulting in a decrease of the oscillator strength.44 The absence of a presumably present ‘dark’ valence state in the reaction path was also earlier discussed with respect to the structurally similar cyclopenta-1,3-diene (CPD) molecule. Early experiments on CPD seemed to confirm the presence of two excited states, i.e. analogous to CHD.45, 46 However, later calculations and further experiments revealed that the dynamics take place not only on one potential energy surface but that rearrangements first take place in the plane of the C5-ring before out-of-plane motions break the C2-symmetry of the system and drive the molecules toward a CoIn with the ground state.47,

48

This step-wise mechanism thus was misinterpreted as dynamics on two

different states. A similar symmetry break was also observed for CHD28, and although the dynamics localize at different dynamophors (the CPD-CoIn is similar to ethylene while CHD shows delocalized dynamics), the similarities between the molecules might be large enough to rationalize the absence of the ‘dark’ 2A-state from the reaction coordinate. Further questions related to CHD concern the branching ratio between an unreactive return to the ground state and the ring opening reaction. As described earlier, gas phase experiments found no evidence for an unreactive return to the ground state8, 42 although theory predicts large quantum yields.30 Several mechanisms were proposed. The CoIn with the ground state was found to be near a ridge where the wavepacket could split, to proceed to either the ground state or to hexatriene.17 Later, several passes over the CoIn were discussed such that momentum toward the ground state structure is gained.17,

24

Alternatively, calculations

suggested that the wavepacket split much earlier in the reaction process, which then could lead to the two different products.28 Finally, there used to be the issue of the absence of a prompt response in photoelectron pumpprobe experiments where only a delayed signal via ionization through Rydberg states was observed after excitation to the bright ππ*-state and a subsequent probe at 400 nm.10 In a very recent TRPES-study, Pemberton et al. resolved this issue and found a prompt response from CHD which they assigned to the 1B-state.16 However, in that study, the authors claim a static photoelectron band for the 1B-state on one hand, while, on the other, they analyze the shift of the photoelectron band along the 2A-state. It is evident that more research along that line is required. In the present study, we used time-resolved photoelectron spectroscopy (TRPES) in combination with the linear response time-dependent density functional theory fewest switching surface hopping (TDDFT-SH) method30,

44, 49

and tried to answer the above

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

mentioned questions. In the following sections, we will first describe our experimental and theoretical methods (section II) before we present and analyze our results (section III). Thereafter, we discuss the impact of our findings and compare them to previous results (section IV) before we summarize and give an outlook to future experiments (section V).

II. METHODS A. Experimental methods CHD was purchased from Sigma-Aldrich with a nominal purity of 97% and used without further purification. Absorption spectra were taken in a 1 cm quartz cuvette (Hellma) under saturated vapor pressure using a Cary 5e photospectrometer (Varian). Our time-resolved photoelectron spectrometer experiment was partly described before.50 Briefly, the magnetic bottle type spectrometer consists of a 5.5 m long flight tube which allows for an E/∆E resolution of better than 100 which means that the energy resolution in the low energy part of the spectrometer is largely defined by the time resolution of our experiment. Femtosecond laser pulses were obtained from a Ti:Sapphire regenerative amplifier (Coherent Legend USP-HE) with an output energy of 4 mJ/pulse at a repetition rate of 1 kHz. Pump and probe pulses at 267, and 400 nm were generated by frequency doubling and consecutive mixing with the fundamental laser beam in beta-barium-borate crystals. In our experiments, the laser beams were attenuated to 1 µJ (267 nm pump) and 20 µJ (400 nm probe) for pump and probe pulses, respectively, and focused weakly into the interaction region by a concave spherical aluminum mirror (f/150 for the pump and f/125 for the probe pulse, ensuring a smaller spot size of the probe pulse in the interaction region). The cross-correlation between pump and probe pulses was measured in xenon and was (125 ± 15) fs. The spectral bandwidth of each of the pump −1

and the probe pulses was around 200 cm . The time delay between the two pulses was controlled by a motorized linear translation stage. Perpendicular to the incoming laser pulses, the sample was inserted in the interaction region of the magnetic bottle spectrometer via a gas needle. At each time delay, the measured pump-probe signal was corrected by subtracting the background signals due to the pump and probe laser pulses alone.

B. Computational methods All calculations were performed with the density functional (DFT)51 and time-dependent density functional theory (TDDFT)52 modules of the TURBOMOLE quantum chemistry

ACS Paragon Plus Environment

Page 4 of 42

Page 5 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

program.53 Non-adiabatic molecular dynamics (MD) were carried out with the TDDFT-SH implementation.30, 44 All calculations employed the def2-SVP54 basis set and the PBE055, 56 hybrid exchange correlation functional. For CHD, usage of the larger basis sets def2-TZVP54 and aug-cc-pVTZ57, 58 leads to excitation energies that are 0.14 and 0.22 eV higher than the def2-SVP results. This leads to a def2-SVP S1 ionization potential that is less than 0.2 eV lower than the ones obtained by def2-TZVP and aug-cc-pVTZ calculations (Table S1 in the supporting information). Since this error is below the accuracy expected for the DFT/TDDFT approach,

59

we use the computationally more efficient def2-SVP basis set for our calculations.

Along the entire ring-opening pathway, differences between the calculations using these basis sets never exceed 0.22 eV. To minimize problems associated with singlet instabilities at conical intersections,60,61 we used the Tamm-Dancoff approximation62 (TDA) to TDDFT. The ground state Boltzmann ensemble of CHD at 300 K was generated by ab initio MD60 with a time step of 50 au (~1.1 fs). A Nosé-Hoover63,

64

thermostate with characteristic

response time of 500 au was used. 200 structures were taken from the MD trajectory to calculate the electronic absorption spectrum and Gaussian broadening with a FWHM of 0.4 eV was applied to each individual spectrum as described in ref 63. The broadened spectra were averaged to obtain the spectrum of the Boltzmann ensemble at 300K. For TDDFT-SH we followed the procedure previously described in refs 30 and 44. 119 structures and velocities were taken from the ground state Boltzmann ensemble and used as starting structures and starting velocities for the excited state TDDFT-SH simulations. To simulate the time-resolved photoelectron ionization spectrum (TRPES), we calculated the total ground state energy of a singly positively charged molecule for each molecular structure from each of the 119 TDDFT-SH trajectories. If the molecule is in the ground state, the ionization potential EIP is calculated by the difference of the ground state energy of the neutral molecule E0(0) and the total energy of the positively charged molecule E0(+1): (1) In the case that the trajectory is currently in the first excited state, the ionization potential is calculated by the ground state ionization potential plus the excitation energyω1. (2) The photoelectron kinetic energy is then obtained by subtracting the energy of the two-photon

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

probe pulse energy of 6.18 eV plus a correction of 1.19 eV from the calculated ionization potential. This correction accounts for the systematic underestimation of the ionization potential of our calculation compared to the experimental value (see Table S1). The energy was chosen to obtain the same vertical ionization potential as in the experiment. For each trajectory the electron kinetic energy was broadened by 0.1 eV in the energy dimension. Afterward, the ionization energies of all trajectories were averaged to obtain the macroscopic ionization potential. The obtained signal was then temporally broadened using a Gaussian line shape with a FWHM of 120 fs.

III. RESULTS A. UV-Vis absorption The UV-absorption spectrum of CHD is shown in Figure 1 and agrees with the spectra from literature.67, 68 The first absorption maximum lies at 251 nm (5.06 eV) and can be assigned to the lowest lying ππ*-state. This means that our pump laser pulse excites the sample on the red edge of the absorption band with little excess energy. The higher lying states around 200 nm were assigned to the π3px-Rydberg state and include a short vibrational progression.67

B. Time-resolved photoelectron spectroscopy The time resolved photoelectron spectrum of CHD excited at 267 nm and probed with 400 nm is shown in Figure 2. The highest lying peak of photoelectron kinetic energies lies at 2.48 eV while the onset of the band is located at 2.65 eV. The peak maximum is slightly below the energy cut-off of one pump (267 nm = 4.65 eV) and two probe photons (400 nm = 3.10 eV), which is 4.65 eV + 2 x 3.10 eV – 8.25 eV = 2.60 eV, where 8.25 eV is the vertical ionization potential (IP) of CHD.69 The differences of the two values will be discussed below. In Figure 3, we present the time zero spectra of CHD (open circles) which we fitted with a set of Gaussians (blue fits). As can be seen, the spectrum consists of a series of narrow banded peaks (FWHM = 30-100 meV) which are superimposed on a broad-banded background. These peaks originate from intermediate Rydberg states which are passed during the two photon probe process. Their assignment is challenging.10, 16 More information on this topic can be found in the supporting information. In addition to the Rydberg peaks, we observe a signal from a broad background. We tentatively fitted the spectrum with two Gaussians (green curves in Figure 3). This background must be assigned either to ionization via a valence state or to direct two photon

ACS Paragon Plus Environment

Page 6 of 42

Page 7 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ionization. More details on this band are found in the upcoming sections. The time-resolved photoelectron spectrum of CHD can be fitted with a monoexponential decay of 60 ± 10 fs with a tiny residual ( 100 fs), the BAC-velocity seems to decrease due to coupling with other modes (dephasing) and is not sufficient to overcome a possible barrier for breaking the C5-C6 bond so that the molecule is forced to return to the ring-closed ground state. This picture emphasizes that not only the momentum but also the geometry where the seam of a conical intersection is crossed plays an important role for product formation. The slow component found in the experimental TRPES, confirms the existence of a slow unreactive channel besides a fast parallel reactive channel.30, 33 For the hexatriene ground state, the simulations allowed us to follow the dynamics for 5 ps and we found that, after the Woodward-Hoffmann allowed conrotatory ring-opening, the tZt conformer is formed within 200 fs while the other conformers are formed afterwards, in slight disagreement with earlier work based on Ehrenfest dynamics simulations.15 Therefore, more work will be necessary to disentangle the ground state rearrangements on a 100 femtosecond time scale. A temperature dependent study of the reaction pathway might give further information.

In the near future, we plan further investigations on cyclic dienes before focusing more on the interplay of TRPES with TDDFT-SH trajectory calculations on larger molecular systems.

Acknowledgment We acknowledge the generous support and technical advice of Prof. Dr. Raimund Feifel and the outstanding technical help of Dr. John Alexander, especially with the development of our

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

measuring software. This work was supported by the Swedish Research Council (VR), the Knut and Alice Wallenberg Foundation, the Lars Hiertas Minne Foundation, and the WennerGren Foundation. In addition, ET thanks the California State University, Long Beach for a start-up fund as well as the CSUPERB biotechnology fund.

ACS Paragon Plus Environment

Page 16 of 42

Page 17 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

References [1] Trulson, M. O.; Dollinger, G. D.; Mathies R. A. Femtosecond Photochemical Ring Opening Dynamics of 1,3-Cyclohexadiene from Resonance Raman Intensities. J. Am. Chem. Soc. 1987,109, 586–587. [2] Trulson, M. O.; Dollinger G. D.; Mathies R. A. Excited State Structure and Femtosecond Ring-Opening Dynamics of 1,3-Cyclohexadiene from Absolute Resonance Raman Intensities. J. Chem. Phys. 1989, 90, 4274–4281. [3] Lawless, M. K.; Wickham S. D.; Mathies, R. A. Resonance Raman view of Pericyclic Photochemical Ring-Opening Reactions: Beyond the Woodward-Hoffmann Rules. Acc. Chem. Res. 1995, 28, 493–502. [4] Pullen, S. H.; Walker, L. A.; Donovan, B.; Sension, R. J. Femtosecond Tansient Absorption Study of the Ring Opening Reaction of 1,3-Cyclohexadiene. Chem. Phys. Lett. 1995, 242, 415–20. [5] Pullen, S. H.; Anderson, N. A.; Walker, L. A. II; Sension, R. J. The Ultrafast Photochemical Ring-Opening Reaction of 1,3-Cyclohexadiene in Cyclohexane. J. Chem. Phys. 1998, 108, 556–563. [6] Lochbrunner, S.; Fuß, W.; Schmid, W. E.; Kompa, K. L. Electronic Relaxation and Ground-State Dynamics of 1,3-Cyclohexadiene and cis-Hexatriene in Ethanol. J. Phys. Chem. A. 1998, 102, 9334–9344. [7] Fuß, W.; Schmid, W. E.; Trushin, S. A. Time-Resolved Dissociative Intense-Laser Field Ionization for Probing Dynamics: Femtosecond Photochemical Ring Opening of 1,3Cyclohexadiene. J. Chem. Phys. 2000, 112, 8347-8362. [8] Dudek, R. C.; Weber, P. M. Ultrafast Diffraction Imaging of the Electrocyclic RingOpening Reaction of 1,3-Cyclohexadiene. J. Phys. Chem. A. 2001, 17, 4167–4171. [9] Garavelli, M.; Page, C. S.; Celani, P.; Olivucci, M.; Schmid, W. E.; Trushin, S. A.; Fuß, W. Reaction Path of a Sub-200 fs Photochemical Electrocyclic Reaction. J. Phys. Chem.

2001, 105, 4458–4469. [10] Kuthirummal, N.; Rudakov, F. M.; Evans, C. L.; Weber, P. M. Spectroscopy and Femtosecond Dynamics of the Ring Opening Reaction of 1,3-Cyclohexadiene. J. Chem. Phys. 2006, 125, 133307.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[11] Rudakov, F. M.; Weber, P. M. Ground State Recovery and Molecular Structure upon Ultrafast Transition through Conical Intersections in Cyclic Dienes. Chem. Phys. Lett.

2009, 470, 187–190. [12] Bühler, C. C.; Minitti, M. P.; Deb, S.; Bao, J.; Weber, P. M.; Ultrafast Dynamics of 1,3 Cyclohexadiene in Highly Excited States. J. At. Mol. Opt. Phys. 2011, 2011, Article ID 637593. [13] Petrovic, V. S.; Schorb, S.; Kim, J.; White, J.; Cryan, J. P.; Glownia, J. M.; Zipp, L.; Broege, D.; Miyabe, S.; Tao, H. et al. Enhancement of Strong-Field Multiple Ionization in the Vicinity of the Conical Intersection in 1,3-Cyclohexadiene Ring Opening. J. Chem. Phys. 2013, 139, 184309. [14] Adachi, S.; Sato, M.; Suzuki, T. Direct Observation of Ground-State Product Formation in a 1,3-Cyclohexadiene Ring-Opening Reaction. J. Phys. Chem. Lett. 2015, 6, 343−346. [15] Minitti, M. P.; Budarz, J. M.; Kirrander, A.; Robinson, J. S.; Ratner, D.; Lane, T. J.; Zhu, D.; Glownia, J. M.; Kozina, M.; Lemke, H. T. et al. Imaging Molecular Motion: Femtosecond X-Ray Scattering of an Electrocyclic Chemical Reaction. Phys. Rev. Lett.

2015, 114, 255501. [16] Pemberton, C. C.; Zhang, Y.; Saita, K.; Kirrander, A.; Weber. P. M. From the (1B) Spectroscopic State to the Photochemical Product of the Ultrafast Ring-Opening of 1,3Cyclohexadiene: A Spectral Observation of the Complete Reaction Path. J. Phys. Chem. A 2015, 119, 8832−8845. [17] Celani, P.; Ottani, S.; Olivucci, M.; Bernardi, F.; Robb, M. A. What Happens During the Picoseberecond Lifetime of 2A1 Cyclohexa-1,3-diene? A CAS-SCF Study of the Cyclohexadiene/Hexatriene Photochemical Interconversion. J. Am. Chem. Soc. 1994, 116, 10141–10151. [18] Celani, P.; Bernardi, F.; Robb, M. A; Olivucci, M. J. Do Photochemical Ring-Openings 1

Occur in the Spectroscopic State? B2 Pathways for the Cyclohexadiene/Hexatriene Photochemical Interconversion. Phys. Chem. 1996, 100, 19364-19366. [19] Bernardi, F.; Olivucci, M.; Robb, M. A. Potential Energy Surface Crossings in Organic Photochemistry. Chem. Soc. Rev. 1996, 25,321-328. [20] Garavelli, M.; Celani, P.; Fato, M.; Bearpark, M. J.; Smith, B. R.; Olivucci, M. J.; Robb, M. A. Relaxation Paths from a Conical Intersection: The Mechanism of Product

ACS Paragon Plus Environment

Page 18 of 42

Page 19 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Formation in the Cyclohexadiene/Hexatriene Photochemical Interconversion. J. Phys. Chem. A, 1997, 101, 2023-2032. [21] Hofmann, A.; Kurtz, L.; de Vivie-Riedle, R. Interaction of Electronic Structure and Nuclear Dynamics on the S1 Reaction Surface for the Ring Opening of Cyclohexadiene. Appl. Phys. B 2000, 71, 391–396. [22] Nonnenberg, C.; Grimm, S; Frank, I. Restricted Open-Shell Kohn--Sham Theory for pi-pi* Transitions. II. Simulation of Photochemical Reactions. J. Chem. Phys. 2003, 119, 11585. [23] Tamura, H.; Nanbu, S.; Nakamura, H.; Ishada, T. A Theoretical Study of Cyclohexadiene/Hexatriene

Photochemical

Interconversion:

Multireference

Configuration Interaction Potential Energy Surfaces and Transition Probabilities for the Radiationless Decays. Chem. Phys. Lett. 2005, 401, 487 – 491. [24] Tamura, H.; Nanbu, S.; Ishada, T.; Nakamura, H. Ab initio Nonadiabatic Quantum Dynamics of Cyclohexadiene/Hexatriene Ultrafast Photoisomerization. J. Chem. Phys.

2006, 124,084313/1-13. [25] Li, A.; Yuan, S.; Dou, Y,; Wang, Y.; Wen, Z, Semiclassical Dynamics Simulation of Photon Induced Ring-Opening of Cyclohexadiene to Hexatriene, Chem. Phys. Lett.,

2009, 478, 28. [26] Li, A.; Yuan, S.; Dou, Y.; Wang, Y.; Wen, Z. Semiclassical Dynamic Simulation of Photon Induced Ring-Opening. Chem. Phys. Lett. 2009, 478, 28–32. [27] Schöborn, J. B.; Sielk, J.; Hartke, B. Photochemical Ring-Opening of Cyclohexadiene: Quantum Wavepacket Dynamics on a Global Ab Initio Potential Energy Surface. J. Phys. Chem. A 2010, 114, 4036-4044. [28] Nenov, A.; Kölle P.; Robb, M. A.; de Vivie-Riedle, R. Beyond the van der Lugt/Oosterhoff Model: When the Conical Intersection Seam and the S1 Minimum Energy Path Do Not Cross. J. Org. Chem. 2010, 75, 123-129. [29] Nenov, A.; de Vivie-Riedle, R. Geometrical and Substituent Effects in Conical Intersections: Linking Chemical Structure and Photoreactivity in Polyenes. J. Chem. Phys. 2011, 135, 034304. [30] Tapavicza, E.; Meyer, A. M.; Furche, F. Unravelling the Details of Vitamin D Photosynthesis by Non-Adiabatic Molecular Dynamics Simulations. Phys. Chem. Chem.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Phys. 2011, 13, 20986–20998. [31] Kim, J.; Tao, H.; Martinez, T. J.; Bucksbaum, P. Ab initio multiple spawning on laserdressed states: a study of 1,3-cyclohexadiene photoisomerization via light-induced conical intersections. J. Phys. B: At. Mol. Opt. Phys. 2015, 48, 164003. [32] Deb, S.; Weber, P. M. The Ultrafast Pathway of Photon-Induced Electrocyclic RingOpening Reactions: The Case of 1,3-Cyclohexadiene. Annu. Rev. Phys. Chem. 2011, 62, 19-39. [33] Arruda, B. C.; Sension, R. J. Ultrafast Polyene Dynamics: the Ring Opening of 1,3Cyclohexadiene Derivatives. Phys. Chem. Chem. Phys. 2014, 16, 4439-4455. [34] Woodward, R. B.; Hoffmann, R. The Conservation of Orbital Symmetry. Angew. Chem. Internat. Edit. 1969, 8, 781–853. [35] van der Lugt, W. T. A. M.; Oosterhoff, L. J. Symmetry Control and Photoinduced Reactions. J. Am. Chem. Soc. 1969, 91, 6042-6049. [36] Schalk, O.; Boguslavskiy, A. E.; Stolow, A.; Schuurman, M., Through-Bond Interactions and the Localization of Excited-State Dynamics. J. Am. Chem. Soc. 2011, 133, 1645116458. [37] Schalk, O.; Boguslavskiy, A. E.; Schuurman, M. S.; Stolow, A. The Dynamophore – Localization of Excited State Dynamics Studied by Time-Resolved Photoelectron Spectroscopy. EPJ Web of Conferences 2013, 41, 02037. [38] Feringa, B. L. (ed.), Molecular Switches, Wiley-VCH, Weinheim, 2001. [39] Griesbeck, A. G.; Oelgemöller, M.; Ghetti, F. (ed.), CRC Handbook of Organic Photochemistry and Photobiology CRC Press, 2012. [40] Kosma, K.; Trushin, S. A.; Fuß, W.; Schmid, W. E. Cyclohexadiene Ring Opening Observed with 13 fs Resolution: Coherent Oscillations Confirm the Reaction Path. Phys. Chem. Chem. Phys. 2009, 11, 172–181. [41] Minnaard N. G., Havinga E. Some Aspects of the Solution Photochemistry of 1,3Cyclohexadiene, (Z)- and (E)- 1,3,5-Hexatriene. Recl. Trav. Chim. Pays. Bas. 1973, 92, 1315–1320. [42] Ruan, C.-Y.; Lobastov, V. A.; Srinivasan, R.; Goodson, B. M.; Ihee, H.; Zewail, A. H. Ultrafast Diffraction and Structural Dynamics: The Nature of Complex Molecules far

ACS Paragon Plus Environment

Page 20 of 42

Page 21 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

from Equilibrium. Proc. Natl. Acad. Sci. 2001, 98, 7117–7122. [43] Kotur, M.; Weinacht, T; Pearson, B. J.; Matsika, S. Closed-Loop Learning Control of Isomerization using Shaped Ultrafast Laser Pulses in the Deep Ultraviolet. J. Chem. Phys., 2009, 130, 134311. [44] Tapavicza, E.; Bellchambers, G. D.; Vincent, J. C.; Furche, F. Ab Initio Non-Adiabatic Molecular Dynamics Simulations. Phys. Chem. Chem. Phys. 2013, 15, 18336–18348. [45] Fuß, W.; Schmid, W. E.; Trushin, S. A. Ultrafast Photochemistry of Cyclopentadiene: Competing Hydrogen Migration and Electrocyclic Ring Closure. Chem. Phys. 2005, 316, 225–234. [46] Schalk, O.; Boguslavskiy, A. E.; Stolow, A., Substituent Effects on Dynamics at Conical Intersections: Cyclopentadienes. J. Phys. Chem. A 2010, 114, 4058-4064. [47] Kuhlman, T. S.; Glover, W. J.; Mori, T.; Møller, K. B.; Martínez, T. J. Between Ethylene and Polyenes - the Non-Adiabatic Dynamics of cis-Dienes. Faraday Discuss. 2012, 157, 193-212. [48] Wolf, T. J. A.; Kuhlman, T. S.; Schalk, O.; Martinez, T. J.; Møller, K. B.; Stolow, A.; Unterreiner, A.-N. Hexamethylcyclopentadiene: A Test Case for the Combination of Time Resolved Photoelectron Spectroscopy and ab initio Multiple Spawning Simulations, Phys. Chem. Chem. Phys. 2014, 16, 11770-11779. [49] Tapavicza, E.; Tavernelli, I.; Rothlisberger, U. Trajectory Surface Hopping Within Linear Response Time-Dependent Density-Functional Theory. Phys. Rev. Lett. 2007, 98, 023001. [50] Kłoda, T.; Matsuda, A.; Karlsson, H. O.; Elshakre, M.; Linusson, P.; Eland, J. H. D.; Feifel, R.; Hansson, T. Strong-Field Photoionization of O2 at Intermediate Light Intensity. Phys. Rev. A 2010, 82, 033431. [51] Häser, M.; Ahlrichs, R. Improvements on the Direct SCF Method. J. Comput. Chem.

1989, 10, 104. [52] Furche, F.; Ahlrichs, R. Adiabatic Time-Dependent Density Functional Methods for Excited State Properties. J. Chem. Phys. 2002, 117, 7433-7447. [53] TURBOMOLE V6.1, TURBOMOLE GmbH, Karlsruhe, 2009; available from http://www.turbomole.com.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[54] Schäfer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian Basis Sets for Atoms Li to Kr. J. Chem. Phys. 1992, 97, 2571-2577. [55] Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett., 1996, 77, 3865-3868. [56] Perdew, J. P.; Ernzerhof, M.; Burke, K. Rationale for Mixing Exact Exchange with Density Functional Approximations. J. Chem. Phys. 1996, 105, 9982-9985. [57] Dunning, T. H. J. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen, J. Chem. Phys. 1989, 90, 1007 [58] Kendall, R. A.; Dunning, T. H. J.; Harrison, R. J. Electron Affinities of the First‐Row Atoms Revisited. Systematic Basis Sets and Wave Functions, J. Chem. Phys. 1992, 96, 6796. [59] Send, R., Kühn, M.; Furche F. Assessing Excited State Methods by Adiabatic Excitation Energies, J. Chem. Theory Comput. 2011, 7, 2376-2386. [60] Levine, B. G.; Ko, C.; Quenneville, J.; Martinez, T. J. Conical Intersections and Double Excitations in Density Functional Theory. Mol. Phys., 2006, 104, 1039-1051. [61] Tapavicza, E.; Tavernelli, I.; Rothlisberger, U.; Filippi, C.; Casida, M. E. Mixed TimeDependent Density-Functional Theory/Classical Trajectory Surface Hopping Study of Oxirane Photochemistry. J. Chem. Phys. 2008, 129, 124108. [62] Hirata, S.; Head-Gordon, M. Time-Dependent Density Functional Theory within the Tamm-Dancoff Approximation. Chem. Phys. Lett. 1999, 314, 291. [63] Elliott, S. D.; Ahlrichs, R.; Hampe, O.; Kappes, M. M. Auto-Ionised Products from the Reaction of Sodium Clusters with Dioxygen: Theory and Experiment. Phys. Chem. Chem. Phys. 2000, 2, 3415-3424. [64] Nosé, S. A. A Molecular Dynamics Method for Simulations in the Canonical Ensemble, Mol. Phys. 1984, 52, 2, 255 -268. [65] Hoover, W. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A,

1985, 31, 1695-1697. [66] Epstein, S.A.; Tapavicza, E.; Furche, F.; Nizkorodov, S. A. Direct Photolysis of Carbonyl Compounds Dissolved in Cloud and Fog Droplets. Atmos. Chem. Phys. 2013, 13, 9461–9477.

ACS Paragon Plus Environment

Page 22 of 42

Page 23 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[67] Merchan, M.; Serrano-Andres, L.; Slater, L. S.; Roos, B. O.; McDiarmid, R.; Xing, X. Electronic Spectra of 1,4-Cyclohexadiene and 1,3-Cyclohexadiene: A Combined Experimental and Theoretical Investigation. Phys. Chem. A 1999, 103, 5468-5476. [68] McDiarmid, R.; Sabljic, A.; Doering, J. P.. Valence Transitions in 1,3-Cyclopentadiene, 1,3-Cyclohexadiene, and 1,3-Cycloheptadiene. J. Chem. Phys. 1985, 83, 2147–2152. [69] Bischof, P.; Heilbronner, E. Photoelektron-Spektren von Cycloalkenen und Cycloalkadienen. Helv. Chim. Acta 1970, 53, 1677-1682 . [70] Schalk, O.; Schuurman, M. S.; Wu, G.; Lang, P.; Mucke, M.; Feifel, R.; Stolow, A. Internal Conversion versus Intersystem Crossing: What Drives the Gas Phase Dynamics of cyclic α,β-enones? J. Phys. Chem. A 2014, 118, 2279-2287. [71] Blanchet, V.; Zgierski, M. Z.; Seideman, T.; Stolow, A. Discerning Vibronic Molecular Dynamics Using Time-Resolved Photoelectron Spectroscopy. Nature 1999, 401, 52-54. [72] Levine, B. G.; Todd J. Martınez, T. J. Ab Initio Multiple Spawning Dynamics of Excited Butadiene: Role of Charge Transfer. J. Phys. Chem. A 2009, 113, 12815–12824. [73] Kim, J.; Tao, H.; White, J. L.; Petrovic, V. S.; Martinez, T. J.; Bucksbaum, P. H. Control of 1,3-Cyclohexadiene Photoisomerization Using Light-Induced Conical Intersections. J. Phys. Chem. A 2012, 116, 2758–2763. [74] Lei, Y.; Wua, H.; Zhenga, X.; Zhaia, G.; Zhub, C. Photo-Induced 1,3-Cyclohexadiene Ring Opening Reaction: Ab initio On-the-Fly Nonadiabatic Molecular Dynamics Simulation. J. Photochem. Photobiol. A: Chemistry 2016, 317, 39–49.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figures

Figure 1. Experimental (black) and calculated (red) normalized absorption spectra of cyclohexa-1,3-diene. The pump pulse of our experiment at 267 nm is indicated with a blue arrow.

ACS Paragon Plus Environment

Page 24 of 42

Page 25 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Time-resolved photoelectron spectrum of cyclohexa-1,3-diene excited at 267 nm and probed at 400 nm.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Time zero spectrum of cyclohexa-1,3-diene (filled circles). The spectrum was fitted with a set of Gaussians (blue color for the Rydberg states and green color for the valence states). Superimposed on the spectrum is the time zero shift of the spectrum (open circles). See text for more details.

ACS Paragon Plus Environment

Page 26 of 42

Page 27 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4: a) Time evolution of the C5-C6 distance (density plot) upon excitation to the S1state of cyclohexa-1,3-diene. The red, solid line represents the population in the first excited state; the green curve represents the excited state population of the successful trajectories and the yellow curve the excited state population of the unsuccessful trajectories. b) Time evolution of the extended bond alternation coordinate BAC (BAC + d(C5-C6) - d(C2-C3), see section III.D for details). The red curve shows the absolute value of the non-adiabatic coupling between the S1- and S0-state in atomic units.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Time evolution of the bond distance between C5 and C6 for all trajectories. The average values of the trajectories that lead to a certain rotational isomer are indicated by the colored lines.

ACS Paragon Plus Environment

Page 28 of 42

Page 29 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6: Time-evolution of the dihedral angles C6-C1-C2-C3 (φ1) and C5-C4-C3-C2 (φ2) for different time intervals of the excited state dynamics. After ring opening, all successful trajectories (non-black) first form the trans-trans conformer (0.15-0.25 ps), characterized by the corners of the plot, i.e. φ1/φ2 combinations -180/-180, -180/180, 180/-180, and 180/180. At later times, during hot ground state dynamics, the trans-Z-trans (tZt, red) cis-Z-trans, trans-Z-cis (tZc or cZt, yellow), and cis-Z-cis (cZc, green) conformers are formed. In the last frame, the different reaction products correspond to the different C5C6 bond distances in Figure 5.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Calculated time-resolved photoelectron spectrum of cyclohexa-1,3-diene upon excitation to the S1-state. The probe pulse was assumed to be 6.2 eV (two photons of 400 nm). A shift of 1.9 eV was applied to correct for the systematic overestimation of ionization potentials.

ACS Paragon Plus Environment

Page 30 of 42

Page 31 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8: Sketch of the potential energy surface and the reaction pathways of cyclohexa-1,3diene after excitation to the lowest lying excited state.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 42

Tables Assignment

3p

3d

4s

4p

4d

5s

5p

6s

Energetic position / eV

0.72 / 1 0.83

1.31 / 1 1.41

1.62

1.81

2.03

2.15

2.29

2.46

Time delay / fs

44

32

32

31

27

19

17

10

Table 1. Energetic positions and time delays for the intermediate Rydberg states observed in 1

the time-resolved photoelectron spectrum of cyclohexa-1,3-diene. Two peaks are assigned for ionization through the 3p- and 3d-Rydberg states.

ACS Paragon Plus Environment

Page 33 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure for TOC

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Experimental (black) and calculated (red) normalized absorption spectra of cyclohexa-1,3-diene. The pump pulse of our experiment at 267 nm is indicated with a blue arrow.

ACS Paragon Plus Environment

Page 34 of 42

Page 35 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Time-resolved photoelectron spectrum of cyclohexa-1,3-diene excited at 267 nm and probed at 400 nm.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Time zero spectrum of cyclohexa-1,3-diene (filled circles). The spectrum was fitted with a set of Gaussians (blue color for the Rydberg states and green color for the valence states). Superimposed on the spectrum is the time zero shift of the spectrum (open circles). See text for more details.

ACS Paragon Plus Environment

Page 36 of 42

Page 37 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a) Time evolution of the C5-C6 distance (density plot) upon excitation to the S1-state of cyclohexa-1,3diene. The red, solid line represents the population in the first excited state; the green curve represents the excited state population of the successful trajectories and the yellow curve the excited state population of the unsuccessful trajectories. b) Time evolution of the extended bond alternation coordinate BAC (BAC + d(C5-C6) + d(C2-C3)). The red curve shows the absolute value of the non-adiabatic coupling between the S1- and S0-state in atomic units. 202x172mm (96 x 96 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Time evolution of the bond distance between C5 and C6 for all trajectories. The average values of the trajectories that lead to a certain rotational isomer are indicated by the colored lines.

ACS Paragon Plus Environment

Page 38 of 42

Page 39 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Time-evolution of the dihedral angles C6-C1-C2-C3 (φ1) and C5-C4-C3-C2 (φ2) for different time intervals of the excited state dynamics. After ring opening, all successful trajectories (non-black) first form the transtrans conformer (0.15-0.25 ps), characterized by the corners of the plot, i.e. φ1/φ2 combinations -180/180, -180/180, 180/-180, and 180/180. At later times, during hot ground state dynamics, the trans-Z-trans (tZt, red) cis-Z-trans, trans-Z-cis (tZc or cZt, yellow), and cis-Z-cis (cZc, green) conformers are formed. In the last frame, the different reaction products correspond to the different C5-C6 bond distances in Figure 5.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Calculated time-resolved photoelectron spectrum of cyclohexa-1,3-diene upon excitation to the S1-state. The probe pulse was assumed to be 6.2 eV (two photons of 400 nm). A shift of 1.9 eV was applied to correct for the systematic overestimation of ionization potentials.

ACS Paragon Plus Environment

Page 40 of 42

Page 41 of 42

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Sketch of the potential energy surface and the reaction pathways of cyclohexa-1,3-diene after excitation to the lowest lying excited state.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Figure

ACS Paragon Plus Environment

Page 42 of 42