Nonadiabatic Effect in Photodissociation Dynamics of Thiophenol via

Jun 1, 2018 - Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque , New Mexico 87131 , United States. J. Phys. Chem...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of Winnipeg Library

A: Spectroscopy, Molecular Structure, and Quantum Chemistry

Nonadiabatic Effect in Photodissociation Dynamics of Thiophenol via the ##* State 1

Guang-Shuang-Mu Lin, Changjian Xie, and Daiqian Xie J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.8b03460 • Publication Date (Web): 01 Jun 2018 Downloaded from http://pubs.acs.org on June 1, 2018

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 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 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.

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 35 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

Nonadiabatic Effect in Photodissociation Dynamics of Thiophenol via the 1ππ* State Guang-Shuang-Mu Lin,† Changjian Xie,*,¶ and Daiqian Xie*,†



Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China



Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM 87131, USA

1

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

Abstract Nonadiabatic photodissociation dynamics of thiophenol (PhSH) and deuterated thiophenol (PhSD) via the 1ππ* state was investigated by a reduced three-dimensional (3D) quantum model based on the associated 3D diabatic potential energy surfaces constructed at the explicitly correlated multireference configuration interaction (MRCI-F12) level with the cc-pVTZ-F12 basis. The lifetimes of the low-lying vibronic S1 states for PhSH and PhSD were calculated using a low-storage filter diagonalization method and in reasonably good agreement with the available experimental results. The nonadiabatic effect was further examined in the photodissociation process by comparing the results in diabatic and adiabatic models. It was found that the adiabatic lifetimes are about 2~4 times shorter than the exact ones in diabatic model for both PhSH and PhSD. More importantly, the exact ground wavefunction including geometric phase (GP) possesses a node along the C-C-S-H/C-C-S-D torsional coordinate, while the node is absent when GP is not included in adiabatic model. The node structure of the wavefunction is a hallmark of GP, which sheds light on the nonadiabatic photodissociation dynamics facilitated by the conical intersections.

2

ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35 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

1. Introduction Nonadiabatic dynamics facilitated by the conical intersections (CIs) has received increasing attentions, because of its important role for understanding the radiationless relaxation mechanisms in many heteroaromatic molecules, such as phenols, pyrroles, and indoles.1-13 The ultraviolet (UV) absorptions of the heteroaromatic molecules dominated by the S1←S0 transitions are the key in photoprotection for humans on the earth.14-15 It is well known that the first and second excited states (S1 and S2) of phenol are coupled in the Franck-Condon (FC) region by the S1/S2 CI,16-17 which is responsible for the 1πσ*-mediated photodissociation ascribed to the nature that the 1

ππ* state has the characteristic of the bound state and the 1πσ* is repulsive along the

O-H dissociation coordinate. The nonadiabatically predissociative nature of the 1ππ* state allows to determine its lifetime, which differs from the case of resonances in complex-forming reactions.18-19 Compared to phenol system, the characteristics of the excited electronic states of thiophenol are quite similar due to S and O atoms with identical number of valence electrons. There are multi-states coupled in the dissociation pathways to the products of H and thiophenoxyl radicals. Apart from the S1/S2 CI in FC region, there is another CI formed between the 1πσ* state and the ground state 1ππ at longer S-H distance (see Fig. 1). It is clarified throughout this paper that 1ππ, 1ππ*, and 1πσ* denote the diabatic states, and S0, S1, and S2 are for the adiabatic labels. There have been many experimental studies on the photodissociation dynamics of thiophenol (PhSH) and deuterated thiophenol (PhSD), including the low bands of the 3

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

absorption spectra,20-21 lifetimes of the S1 vibronic states,20-22 and the product branching ratios (X/A) between the ground and first electronic excited states of the C6H5S radical.21 Particularly, the kinetic energy distributions for the ground and excited states of the products were measured at several wavelengths as well,17, 20, 23 and the dominant product vibrational states in the distributions were found to be an odd number of quanta in aʹʹ C-C twist mode v16a.17 This mode specificity in the product state distributions were attributed to the geometric phase (GP) effect in the nonadiabatic photodissociation facilitated by the CIs.24 The GP is associated with the sign change (eiπ = -1) of the electronic wavefunction when the system encircles the CI in a close loop, which makes the wavefunction double-valued.25-28 The recent theoretical study demonstrated that the size of the nonadiabatic effect in the tunneling-facilitated photodissociation of phenol (S1) was extremely large, the adiabatic predicted lifetimes were about two orders of magnitude shorter than the exact ones.29 More importantly, the adiabatic wavefunctions without GP have been shown to be significantly different from the exact wavefunctions (implicitly including GP),29-32 which suggested that the neglect of GP would lead to incorrect results in describing the dynamics near the CIs. Compared to phenol, the lifetimes of the S1(1ππ*) state of thiophenol were found to be shorter and in femtosecond by several experimental groups.20-22 As discussed above, the nonadiabatic effect had been found to have significant impact on the lifetimes of phenol (S1), but for the nonadiabatic photodissociation of thiophenol S1 state, the nonadiabatic effect is still unclear. The size of the nonadiabatic effect in the photodissociation dynamics can be taken into 4

ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35 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

account by comparing the results in diabatic (exact) and adiabatic (on a single surface) models. The adiabatic treatment on a single surface ignores the couplings with other states and the derivative couplings (including GP and diagonal Born-Oppenheimer corrections33). It should be noted that the GP effect is in quantum community and thus cannot be treated by classical approach. The diabatic potential energy surfaces (PESs) are the cornerstone in nonadiabatic quantum dynamical calculations, a deep understanding of the photodissociation dynamics can be achieved by the state-to-state quantum dynamical calculations on the available diabatic PES, such as in the cases of H2O and NH3.34-47 It is still a challenging to construct the full dimensional diabatic PESs for polyatomic systems. In the past few years, the full dimensional diabatic PESs of phenol have been constructed based on the ab initio energy gradients and derivative couplings.48-50 While for the thiophenol, there were only two-dimensional (2D) diabatic PESs available

at

low

computational

level.51-52

Recently,

we

constructed

a

three-dimensional(3D) diabatic PESs with multi-reference configuration interaction method with explicitly correction (MRCI-F12) based on the linear-vibronic coupling method.53 The vibronic energy levels of S1 states of PhSH and PhSD calculated by a reduced 3D quantum dynamical model were in good agreement with the available experimental results, which validates the PESs and reduce dynamical model are reasonably accurate.53 In this work, the lifetimes of PhSH and PhSD were calculated in a reduced 3D model, and the nonadiabatic effect in the photodissociation of PhSH and PhSD was 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

Page 6 of 35

further investigated. The theoretical results are in reasonably good agreement with the available experimental values. This paper is organized as follows, section 2 presents the description of PESs and the wave packet propagation method. Results and discussions are given in section 3, and a brief summary is presented in section 4. 2. Theory In order to investigate the photodissociation dynamics of thiophenol, the adiabatic PESs are taken from our recent work,53 and the construction of the diabatic PESs were achieved by the linear-vibronic coupling method.54-55 The reduced 3D PESs include three degrees of freedom, S-H stretch, C-S-H bend, and C-C-S-H torsion, in which the CIs are one-dimensional (1D). As shown in Fig. 2, the energy of the CI seam is sensitively dependent on the C-S-H bending coordinate. At the planar geometry (φ=0), the 1ππ, 1ππ*, and 1πσ* states are of Aʹ, Aʹ, and Aʹʹ symmetries, respectively, whereas at φ=90 deg, the states have Aʹ, Aʹʹ, and Aʹ symmetries, respectively. In this work, in order to constrain the correct symmetry changes for the C-C-S-H torsional angle varying from 0 to 90 deg, we instead used the expressions for the non-diagonal elements,

W12 = λ sin ϕ ,

(1)

and

W23 =

λ 2

sin 2ϕ ,

(2) in which the symbols 1, 2, and 3 denote the 1ππ, 1πσ*, and 1ππ* states, respectively. It is clear that W12 and W23 in Eqs. (1) and (2) satisfy the linear-vibronic coupling 6

ACS Paragon Plus Environment

Page 7 of 35 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

conditions as well. It should be noted that the new treatment for the non-diagonal elements in diabatic PESs would not qualitatively change the photodissociation dynamics near the FC region, as discussed below. The nonadiabatic coupling parameters λ were determined in the high-symmetry Cs planar subspace with numerical differentiation of the following equation54-55 12

 1 ∂ 2 (V1 − V2 ) 2  λ= ,  ∂ϕ 2 8  R ,θ ,ϕ =0 (3) where Vi are the adiabatic potential energies. The final diabatic PESs were fitted by the cubic spline method.56-57 Fig. 3 shows the 3D plots of the diabatic and adiabatic PESs. Following our recent work,53 the reduced 3D quantum dynamical model was used to study the photodissociation dynamics of thiophenol. The nuclear Hamiltonian in the diabatic representation is given as follows

 1 0 0  V11 V12 V13    Hˆ = Tˆ  0 1 0  + V21 V22 V23  ,      0 0 1 V31 V32 V33  d

(4) in which Tˆ is the kinetic operator,

ˆj 2 h2 ∂2 , Tˆ = − + 2 µ R ∂R 2 2 µ R R 2 (5)

where

µR

is reduced mass µ R =

m H m C 6 H 5S m H + mC 6 H 5S

, j is the angular momentum operator of

7

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 8 of 35

thiophenoxyl, and the total angular moment is set to zero. The total wavefunction is expressed as a vector,

ψ 1    Ψ = ψ 2  , ψ 3    (6) in which three components are the wavefunctions for the 1ππ, 1πσ*, and 1ππ* states, respectively. Each nuclear wavefunction is expanded in a direct-product basis in a mixed representation:34

ψ i = ∑ Cαi jm α jm , α jm

(7) i in which Cα jm are the expansion coefficients, α denotes the DVR (discrete variable

representation58) grid index for radial coordinates R. Specifically, the Fourier bases eimφ (m=0, ±1, …) were used to represent the wavefunctions of the whole nonadiabatic system, which allows both the odd and even reflection symmetries.

jm

imϕ

is the spherical harmonics Ne Pj (θ ) , where Pj (θ ) is the associated m

m

Legendre polynomial, and N is the normalized factor. The photodissociation on the S1 state was simulated assuming a Condon excitation, in which the initial wave packet is a vibrational eigenfunction of the S0 state of thiophenol placed vertically on the first singlet excited 1ππ* (S1) state. The wave packet is propagated using the Chebyshev propagator59

Ψk = 2DHsΨk−1 − D2Ψk −2 , k ≥ 2 ,

(8) 8

ACS Paragon Plus Environment

Page 9 of 35 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

with

Ψ1 = DHsΨ0

(9)

and

0  Ψ0 =  0  Ψ   i

,

(10)

where

Ψi is the S0 vibrational eigenfunction, which is obtained by the iterative

Lanczos algorithm.60 The Hamiltonian is scaled to the spectral range of (-1,1) via

Ηs = (Η − H+I)/ H− , (11) +

in which the spectral medium ( H = (Hmax + Hmin )/ 2 ) and half width −

( H = (Hmax − Hmin )/ 2 ) were determined by the spectral extrema,

Hmax and Hmin ,

which can be readily estimated. To avoid reflection, the damping functions (D) were used at the edge of the radial grids. For the S1 vibronic resonances, we have also determined their complex energies (E-iГ/2) using a low-storage filter diagonalization method.60 The Chebyshev correlation function was used to build an energy-localized Hamiltonian matrix, from which the complex energy of the resonance is obtained by diagonalization. The numerical parameters used in wave packet calculations are listed in Table 1.

9

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 10 of 35

3. Results and discussions Table 2 lists the energy levels of S1 vibronic states of thiophenol (PhSH) and deuterated thiophenol (PhSD) calculated on the new diabatic PESs. In particular, the calculated vibrational frequency for the C-S-D bend state is 41015.8-40418.5=597.3 cm-1, which is close to the experimental value (608 cm-1)21 and that on our previous PES.53 It can be readily seen that the energies of PhSD vibronic states are lower than those of the PhSH ones, due to the isotope effect. The comparison between diabatic and adiabatic energies is discussed below. Due to the nonadiabatic coupling between the bound state 1ππ* and repulsive state πσ*, the S1 state is non-stable and has a lifetime in femtosecond range.20-22 This lifetime is shorter than that of phenol, due apparently to the energy barrier caused by the S1/S2 CI much lower than that of the phenol PESs. Table 3 lists the full-width at half maximum (FWHM) and lifetimes of three low-lying vibronic states of S1 state for PhSD and PhSH along with the available experimental data21 for comparison. It is clear that the lifetimes are mode specific and isotope dependent, which is consistent with previous quantum dynamical calculations in phenol.

29, 61-64

The lifetimes of

these low-lying vibronic states are ranging from 10-40 fs, which is consistent with the observed lifetime in femtosecond range for S1 state in experiment. 20-22 In particular, it can be readily seen that the predicted widths of the ground vibronic state 00 of PhSD and PhSH are 143.9 and 186.7 cm-1, respectively, which are in good agreement with the experimental data (190 and 195 cm-1). On the other hand, the lifetimes of PhSD vibronic states were found to slightly (3.5 with   R − Rend   =1, otherwise Chebyshev steps: 10,000

16

ACS Paragon Plus Environment

Rend = 6.0 ;

Page 16 of 35

Page 17 of 35 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

Table 2. Energies (in cm-1) of S1 state relative to ground state minimum on PES for PhSH/PhSD in adiabatic and diabatic models. m and n in the label mn denote the vibrations (2 and 3 represent the CSH/CSD bending and CCSH/CCSD torsion modes, respectively) and quantum numbers. Diabatic 0

0

2

1

31

Adiabatic

PhSH

PhSD

PhSH

PhSD

40957.2 41889.9 40999.7

40418.5 41015.8 40430.2

40854.9 41740.9 40872.6

40357.2 40936.3 40358.7

17

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 18 of 35

Table 3. Full-width at half maximum (FWHM, in cm-1) and lifetime (in fs) of the low-lying vibronic S1 state PhSD and PhSH. Experimental results were taken from Refs 20 and 21.

FWHM 00 Lifetime FWHM 31 Lifetime FWHM 21 Lifetime

Adiabatic

Diabatic

Exp.

PhSD

428.7

143.9

190

PhSH

834.7

186.7

195

PhSD

12.3

36.6

27.7

PhSH

6.3

28.2

27.0

PhSD

474.5

181.9

PhSH

967.1

284.2

PhSD

11.1

29.0

PhSH

5.4

18.5

PhSD

564.7

285.5

PhSH

1108

413.2

PhSD

9.3

18.4

PhSH

4.8

12.7

18

ACS Paragon Plus Environment

Page 19 of 35 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.

Blancafort, L.; Cohen, B.; Hare, P. M.; Kohler, B.; Robb, M. A. Singlet excited-state dynamics of

5-fluorocytosine and cytosine: An experimental and computational study. J. Phys. Chem. A 2005, 109, 4431-4436. 2.

Tseng, C.-M.; Lee, Y. T.; Lin, M.-F.; Ni, C.-K.; Liu, S.-Y.; Lee, Y.-P.; Xu, Z. F.; Lin, M. C.

Photodissociation dynamics of phenol. J. Phys. Chem. A 2007, 111, 9463-9470. 3.

Tseng, C.-M.; Lee, Y. T.; Ni, C.-K.; Chang, J.-L. Photodissociation dynamics of the chromophores of

the amino acid tyrosine: p-methylphenol, p-ethylphenol, and p-(2-Aminoethyl)phenol. J. Phys. Chem. A 2007, 111, 6674-6678. 4.

Iqbal, A.; Pegg, L.-J.; Stavros, V. G. Direct versus indirect H atom elimination from photoexcited

phenol molecules. J. Phys. Chem. A 2008, 112, 9531-9534. 5.

Iqbal, A.; Cheung, M. S. Y.; Nix, M. G. D.; Stavros, V. G. Exploring the time-scales of H-atom

detachment from photoexcited phenol-h6 and phenol-d5: Statistical vs nonstatistical decay. J. Phys. Chem. A 2009, 113, 8157-8163. 6.

King, G. A.; Oliver, T. A. A.; Nix, M. G. D.; Ashfold, M. N. R. High resolution photofragment

translational spectroscopy studies of the ultraviolet photolysis of phenol-d5. J. Phys. Chem. A 2009, 113, 7984-7993. 7.

Tseng, C.-M.; Lee, Y. T.; Ni, C.-K. Photodissociation dynamics of n-methylindole, n-methylpyrrole,

and anisole. J. Phys. Chem. A 2009, 113, 3881-3885. 8.

Iqbal, A.; Stavros, V. G. Exploring the time scales of H-atom elimination from photoexcited indole.

J. Phys. Chem. A 2010, 114, 68-72. 9.

Wang, F.; Neville, S. P.; Wang, R.; Worth, G. A. Quantum dynamics study of photoexcited aniline. J.

Phys. Chem. A 2013, 117, 7298-7307. 10. Benda, Z.; Szalay, P. G. Details of the excited-state potential energy surfaces of adenine by coupled cluster techniques. J. Phys. Chem. A 2014, 118, 6197-6207. 11. Capello, M. C.; Broquier, M.; Ishiuchi, S.-I.; Sohn, W. Y.; Fujii, M.; Dedonder-Lardeux, C.; Jouvet, C.; Pino, G. A. Fast nonradiative decay in o-aminophenol. J. Phys. Chem. A 2014, 118, 2056-2062. 12. Han, S.; You, H. S.; Kim, S.-Y.; Kim, S. K. Dynamic role of the intramolecular hydrogen bonding in nonadiabatic chemistry revealed in the UV photodissociation reactions of 2-fluorothiophenol and 2-chlorothiophenol. J. Phys. Chem. A 2014, 118, 6940-6949. 13. Staniforth, M.; Young, J. D.; Cole, D. R.; Karsili, T. N. V.; Ashfold, M. N. R.; Stavros, V. G. Ultrafast excited-state dynamics of 2,4-dimethylpyrrole. J. Phys. Chem. A 2014, 118, 10909-10918. 14. Sobolewski, A. L.; Domcke, W.; Dedonder-Lardeux, C.; Jouvet, C. Excited-state hydrogen detachment and hydrogen transfer driven by repulsive 1πσ* states: A new paradigm for nonradiative decay in aromatic biomolecules. Phys. Chem. Chem. Phys. 2002, 4, 1093-1100. 15. Baker, L. A.; Marchetti, B.; Karsili, T. N. V.; Stavros, V. G.; Ashfold, M. N. R. Photoprotection: extending lessons learned from studying natural sunscreens to the design of artificial sunscreen constituents. Chem. Soc. Rev. 2017, 46, 3770-3791. 16. Ashfold, M. N. R.; Cronin, B.; Devine, A. L.; Dixon, R. N.; Nix, M. G. D. The role of πσ* excited states in the photodissociation of heteroaromatic molecules. Science 2006, 312, 1637-1640. 17. Ashfold, M. N. R.; Devine, A. L.; Dixon, R. N.; King, G. A.; Nix, M. G. D.; Oliver, T. A. A. Exploring nuclear motion through conical intersections in the UV photodissociation of phenols and thiophenol. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 12701-12706. 19

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

18. Guo, H. Quantum dynamics of complex-forming bimolecular reactions. Int. Rev. Phys. Chem. 2012, 31, 1-68. 19. Zuo, J.; Hu, X.; Xie, D. Quantum dynamics of oxyhydrogen complex-forming reactions for the HO2 and HO3 systems. Chin. J. Chem. Phys. 2018, 31, 123-134. 20. Devine, A. L.; Nix, M. G. D.; Dixon, R. N.; Ashfold, M. N. R. Near-ultraviolet photodissociation of thiophenol. J. Phys. Chem. A 2008, 112, 9563-9574. 21. You, H. S.; Han, S.; Lim, J. S.; Kim, S. K. (ππ*/πσ*) Conical intersection seam experimentally observed in the S-D bond dissociation reaction of thiophenol-d1. J. Phys. Chem. Lett. 2015, 6, 3202-3208. 22. Ovejas, V.; Fernández-Fernández, M.; Montero, R.; Longarte, A. On the ultrashort lifetime of electronically excited thiophenol. Chem. Phys. Lett. 2016, 661, 206-209. 23. Lim, I. S.; Lim, J. S.; Lee, Y. S.; Kim, S. K. Experimental and theoretical study of the photodissociation reaction of thiophenol at 243 nm: Intramolecular orbital alignment of the phenylthiyl radical. J. Chem. Phys. 2007, 126, 034306. 24. Nix, M. G. D.; Devine, A. L.; Dixon, R. N.; Ashfold, M. N. R. Observation of geometric phase effect induced photodissociation dynamics in phenol. Chem. Phys. Lett. 2008, 463, 305-308. 25. Herzberg, G.; Longuet-Higgins, H. C. Intersection of potential energy surfaces in polyatomic molecules. Discuss. Faraday Soc. 1963, 35, 77-82. 26. Mead, C. A.; Truhlar, D. G. On the determination of Born–Oppenheimer nuclear motion wave functions including complications due to conical intersections and identical nuclei. J. Chem. Phys. 1979, 70, 2284-2296. 27. Berry, M. V. Quantal phase factors accompanying adiabatic changes. Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences 1984, 392, 45-57. 28. Mead, C. A. The geometric phase in molecular systems. Rev. Mod. Phys. 1992, 64, 51-85. 29. Xie, C.; Ma, J.; Zhu, X.; Yarkony, D. R.; Xie, D.; Guo, H. Nonadiabatic tunneling in photodissociation of phenol. J. Am. Chem. Soc. 2016, 138, 7828-7831. 30. Xie, C.; Yarkony, D. R.; Guo, H. Nonadiabatic tunneling via conical intersections and the role of the geometric phase. Phys. Rev. A 2017, 95, 022104. 31. Xie, C.; Kendrick, B. K.; Yarkony, D. R.; Guo, H. Constructive and destructive interference in nonadiabatic tunneling via conical intersections. J. Chem.Theory Comput. 2017, 13, 1902-1910. 32. Ryabinkin, I. G.; Izmaylov, A. F. Geometric phase effects in dynamics near conical intersections: Symmetry breaking and spatial localization. Phys. Rev. Lett. 2013, 111, 220406. 33. Baer, M. Beyond Born-Oppenheimer: Electronic Nonadiabatic Coupling Terms and Conical Intersections. John Wiley & Sons, Inc.: Hoboken, New Jersey, 2006. 34. Hu, X.; Zhou, L.; Xie, D. State-to-state photodissociation dynamics of the water molecule. WIREs: Comput. Mol. Sci. 2017, e1350. 35. Zhou, L.; Jiang, B.; Xie, D.; Guo, H. State-to-state photodissociation dynamics of H2O in the B-band: Competition between two coexisting nonadiabatic pathways. J. Phys. Chem. A 2013, 117, 6940-6947. 36. Zhou, L.; Lin, G.; Xie, D. State to state photodissociation dynamics of D2O in the B band. J. Chem. Phys. 2013, 139, 114303. 37. Ma, J.; Xie, C.; Zhu, X.; Yarkony, D. R.; Xie, D.; Guo, H. Full-dimensional quantum dynamics of vibrationally mediated photodissociation of NH3 and ND3 on coupled ab initio potential energy surfaces: Absorption spectra and NH2(A2A1)/NH2(X2B1) branching ratios. J. Phys. Chem. A 2014, 118, 20

ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35 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

11926-11934. 38. Xie, C.; Ma, J.; Zhu, X.; Zhang, D.; Yarkony, D. R.; Xie, D.; Guo, H. Full-dimensional quantum state-to-state nonadiabatic dynamics for photodissociation of ammonia in its A-band. J. Phys. Chem. Lett. 2014, 5, 1055-1060. 39. Zhou, L.; Xie, D.; Sun, Z.; Guo, H. Product fine-structure resolved photodissociation dynamics: The A band of H2O. J. Chem. Phys. 2014, 140, 024310. 40. Han, S.; Zhou, L.; Xie, D. State to state photodissociation dynamics of vibrationally excited D2O in B band. Chin. J. Chem. Phys. 2015, 28, 396-402. 41. Zhou, L.; Xie, D. Full-dimensional quantum dynamics of vibrational mediated photodissociation of HOD in its B band. J. Phys. Chem. A 2015, 119, 12062-12072. 42. Zhou, L.; Xie, D.; Guo, H. Signatures of non-adiabatic dynamics in the fine-structure state distributions of the OH(X/A) products in the B-band photodissociation of H2O. J. Chem. Phys. 2015, 142, 124317. 43. Jiang, B.; Xie, D.; Guo, H. Communication: State-to-state differential cross sections for H2O(B) photodissociation. J. Chem. Phys. 2011, 134, 231103. 44. Jiang, B.; Xie, D. Theoretical studies for photodissociation dynamics of small molecules. Progress in Chemistry 2012, 24, 1120-1128. 45. Jiang, B.; Xie, D.; Guo, H. State-to-state photodissociation dynamics of triatomic molecules: H2O in the B band. J. Chem. Phys. 2012, 136, 034302. 46. Lai, W.; Lin, S. Y.; Xie, D.; Guo, H. Nonadiabatic dynamics of A-state photodissociation of ammonia: A four-dimensional wave packet study. J. Phys. Chem. A 2010, 114, 3121-3126. 47. Lin, G.; Zhou, L.; Xie, D. Theoretical study of the state-to-state photodissociation dynamics of the vibrationally excited water molecule in the B band. J. Phys. Chem. A 2014, 118, 9220-9227. 48. Zhu, X.; Yarkony, D. R. Fitting coupled potential energy surfaces for large systems: Method and construction of a 3-state representation for phenol photodissociation in the full 33 internal degrees of freedom using multireference configuration interaction determined data. J. Chem. Phys. 2014, 140, 024112. 1

49. Zhu, X.; Malbon, C. L.; Yarkony, D. R. An improved quasi-diabatic representation of the 1, 2, 3 A coupled adiabatic potential energy surfaces of phenol in the full 33 internal coordinates. J. Chem. Phys. 2016, 144, 124312. 50. Yang, K. R.; Xu, X.; Zheng, J.; Truhlar, D. G. Full-dimensional potentials and state couplings and multidimensional tunneling calculations for the photodissociation of phenol. Chem. Sci. 2014, 5, 4661-4680. 51. Venkatesan, T. S.; Ramesh, S. G.; Lan, Z.; Domcke, W. Theoretical analysis of photoinduced H-atom elimination in thiophenol. J. Chem. Phys. 2012, 136, 174312. 52. An, H.; Choi, H.; Lee, Y. S.; Baeck, K. K. Factors affecting the branching ratio of photodissociation: Thiophenol studied through quantum wavepacket dynamics. Chemphyschem 2015, 16, 1529-1534. 53. Lin, G.; Xie, C.; Xie, D. Three-dimensional diabatic potential energy surfaces for the photodissociation of thiophenol. J. Phys. Chem. A 2017, 121, 8432-8439. 54. Köppel, H.; Gronki, J.; Mahapatra, S. Construction scheme for regularized diabatic states. J. Chem. Phys. 2001, 115, 2377-2388. 55. Köppel, H.; Schubert, B. The concept of regularized diabatic states for a general conical intersection. Mol. Phys. 2006, 104, 1069-1079. 56. Sathyamurthy, N.; Raff, L. M. Quasiclassical trajectory studies using 3D spline interpolation of ab 21

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 22 of 35

initio surfaces. J. Chem. Phys. 1975, 63, 464-473. 2

57. Xu, C.; Xie, D.; Zhang, D.; Lin, S. Y.; Guo, H. A new ab initio potential-energy surface of HO2(X A″) and quantum studies of HO2 vibrational spectrum and rate constants for the H+O2↔O+OH reactions. J. Chem. Phys. 2005, 122, 244305. 58. Light, J. C.; Carrington, T. Discrete-Variable Representations and their Utilization. In Adv. Chem. Phys., John Wiley & Sons, Inc.: 2007; pp 263-310. 59. Chen, R.; Guo, H. The Chebyshev propagator for quantum systems. Comput. Phys. Commun. 1999, 119, 19-31. 60. Guo, H. Recursive Solutions to Large Eigenproblems in Molecular Spectroscopy and Reaction Dynamics. In Rev. Comput. Chem., John Wiley & Sons, Inc.: 2007; pp 285-348. 61. Lan, Z.; Domcke, W.; Vallet, V.; Sobolewski, A. L.; Mahapatra, S. Time-dependent quantum 1

wave-packet description of the πσ* photochemistry of phenol. J. Chem. Phys. 2005, 122, 224315. 62. Vieuxmaire, O. P. J.; Lan, Z.; Sobolewski, A. L.; Domcke, W. Ab initio characterization of the conical intersections involved in the photochemistry of phenol. J. Chem. Phys. 2008, 129, 224307. 63. An, H.; Baeck, K. K. Quantum wave packet propagation study of the photochemistry of phenol: Isotope effects (Ph-OD) and the direct excitation to the 1πσ* state. J. Phys. Chem. A 2011, 115, 13309-13315. 64. Dixon, R. N.; Oliver, T. A. A.; Ashfold, M. N. R. Tunnelling under a conical intersection: Application to the product vibrational state distributions in the UV photodissociation of phenols. J. Chem. Phys. 2011, 134, 194303. 65. Xie, C.; Malbon, C. L.; Yarkony, D. R.; Guo, H. Dynamic mapping of conical intersection seams: A general method for incorporating the geometric phase in adiabatic dynamics in polyatomic systems. J. Chem. Phys. 2017, 147, 044109. 66. Xie, C.; Malbon, C. L.; Yarkony, D. R.; Xie, D.; Guo, H. Signatures of a conical intersection in adiabatic dissociation on the ground electronic state. J. Am. Chem. Soc. 2018, 140, 1986-1989. 67. Barclay, V. J.; Dateo, C. E.; Hamilton, I. P.; Kendrick, B.; Pack, R. T.; Schwenke, D. W. Anomalous symmetries of the rovibrational states of HO2: Consequences of a conical intersection. J. Chem. Phys. 1995, 103, 3864-3867. 68. Joubert-Doriol, L.; Ryabinkin, I. G.; Izmaylov, A. F. Geometric phase effects in low-energy dynamics near conical intersections: A study of the multidimensional linear vibronic coupling model. J. Chem. Phys. 2013, 139, 234103. 69. Ryabinkin, I. G.; Joubert-Doriol, L.; Izmaylov, A. F. Geometric phase effects in nonadiabatic dynamics near conical intersections. Acc. Chem. Res. 2017, 50, 1785-1793. 70. Juanes-Marcos, J. C.; Althorpe, S. C.; Wrede, E. Theoretical study of geometric phase effects in the hydrogen-exchange reaction. Science 2005, 309, 1227-1230. 71. Althorpe, S. C. General explanation of geometric phase effects in reactive systems: Unwinding the nuclear wave function using simple topology. J. Chem. Phys. 2006, 124, 084105. 72. Bouakline, F. Investigation of geometric phase effects in photodissociation dynamics at a conical intersection. Chem. Phys. 2014, 442, 31-40.

22

ACS Paragon Plus Environment

Page 23 of 35 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 captions: Figure 1. The 1D adiabatic and diabatic PES as a function of S-H stretch at several C-C-S-H torsion angles (the C-S-H bend angle was fixed at the ground state equilibrium geometry). Figure 2. The energy of CI seam as function of C-S-H bend angle. Figure 3. The 3D adiabatic (left panel) and diabatic (right panel) PESs as functions of S-H stretch and C-C-S-H torsion (the C-S-H bend angle was fixed at the ground state equilibrium geometry). Figure 4. Moduli of the wavefunction of 00 for the S1 PhSD (upper panel) and PhSH (lower panel) from the diabatic(left) and adiabatic(right) calculations. The angle θ fixed at 97.4 and 94.5 deg for PhSH and PhSD, respectively. Red dot denotes the CI. Figure 5. Moduli of the wavefunction of 31 for the S1 PhSD (upper panel) and PhSH (lower panel) from the diabatic(left) and adiabatic(right) calculations. The angle θ fixed at 97.4 and 94.5 deg for PhSH and PhSD, respectively. Red dot denotes the CI.

23

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

Fig. 1

24

ACS Paragon Plus Environment

Page 24 of 35

Page 25 of 35 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

Fig. 2

25

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

Fig. 3

26

ACS Paragon Plus Environment

Page 26 of 35

Page 27 of 35 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

Fig. 4

27

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

Fig. 5

28

ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35 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

TOC Graphic

29

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

219x151mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35 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

148x118mm (300 x 300 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

176x91mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35 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

621x644mm (150 x 150 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

611x644mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 34 of 35

Page 35 of 35 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

611x329mm (150 x 150 DPI)

ACS Paragon Plus Environment