Photoionization Spectroscopic and Theoretical Study on the Molecular Structures of cis- and trans-3-Chlorothioanisole.

Resonance-enhanced two-photon ionization (R2PI) and mass-analyzed threshold ionization (MATI) spectra are measured for the cis- and trans-3-chlorothioanisole (3ClTA). The first electronic excitation energy (E 1) and the adiabatic ionization energy (IE) of the cis-rotamer are determined to be 33 959±3 and 65 326±5 cm-1, respectively, and those of the trans-rotamer are determined to be 34102±3 and 65 471±5 cm-1, respectively. Density functional theory (DFT) calculations confirm that both the cis- and trans-rotamers of 3ClTA are stable and coexist in their respective S0, S1, and D0 states. Both rotamers adopt planar structures with cis- being slightly more stable than trans- in the respective S0, S1, and D0 states. The conformation, substitution, and isotope effects on the molecular structure, active vibrations, and electronic transition and ionization energies of 3ClTA are analyzed.

Introduction

There have been numerous experimental and theoretical studies for investigating the stable structures and properties of anisole and thioanisole (TA) derivatives in their respective electronic ground (S0), first excited (S1), and cationic ground (D0) states.[1−9] Due to strong conjugation effect between the benzene ring and the substituent functional groups, most of those molecules are found to be planar. However, several exceptions have been observed. For example, a nonplanar stable structure was observed for trans-2-fluoroanisole (2-F anisole) in its S0 state with Fourier transform infrared (FTIR) spectroscopy[10] and a small structure tilt was observed for trans-2-fluorothioanisole (2FTA) in its S1 state through resonance-enhanced two-photon ionization (R2PI) spectroscopy.[11] In 2FTA, the enhanced steric effect and inductive effect caused by the 2-substituted (ortho-substituted) F atom play important roles.

For TA derivatives, there exist a lot of controversies, from previous calculations and spectroscopic results,[6,9,12−14] about their stable structures, particularly in their S1 states. Some theoretical calculations predict their molecular skeletons to be nonplanar in the S1 state,[12,15] in which the thiomethyl (−SCH3) group was nearly perpendicular to the benzene ring. The whole molecule has a vertical gauche-like structure. This contradicts the general rule that the benzene ring tends to keep a planar structure as favored by π-orbital delocalization.[16] The nonplanar structures obtained by theoretical calculations do not exactly agree with the spectroscopic experiments.[6,9,15] Nevertheless, most TA derivatives have now been determined to be planar or quasi-planar,[6,8,9,12−15] and the competition between the steric effect and the inductive effect caused by −SCH3 was used to explain the difference between calculations and spectroscopic measurements.[6,9,12] However, the reason for the structure tilt seen in the theoretical calculation is still unclear, and it is worthwhile to investigate the stable structures of halogen-substituted anisole and TA derivatives in various electronic states.

The deformation changes of benzene derivatives upon photoexcitation and photoionization are usually similar to each other for molecules with similar structures.[1,6,16−20] For example, benzene ring is usually enlarged in the S1 state and a quinoidal structure is often found in the D0 state. It is reasonable to assume that the conjugation effects of −OCH3 and −SCH3 groups with the benzene ring are similar since O and S atoms belong to the same group in the periodic table. Compared with the O atom, the S atom is larger in size, which could make the −SCH3 group more flexible than the −OCH3 group.[2,14] Investigation on how the −SCH3 group substitution affects the molecular properties is limited.

Besides the conjugation effect mentioned above, the steric effect and inductive effect caused by the halogen substitutions should also be considered. For example, the cis-rotamer is observed to be stable for 3-fluoroanisole (3-F anisole)[21] but not for 2-F anisole.[10,22] That is due to the fact that the steric effect between the F atom and the −OCH3 group is weaker in 3-substituted (meta-substituted) anisole than that in the 2-substituted one. The transition energies (E1) of Cl-substituted anisoles were found to be red-shifted compared to those of anisole,[23] while those of F-substituted anisoles were found to be blue-shifted.[17,24−26] This is mainly caused by the different inductive effects of F and Cl atoms.

Based on the above discussion, the 3-chlorothioanisole (3ClTA) molecule can be a suitable system for simultaneously comparing the substitution effects between −SCH3 and −OCH3 groups and between Cl and F atoms. To the best of our knowledge, the vibrational spectra of 3ClTA in both S1 and D0 states have not been reported, and its structural properties have not been investigated. In this work, we reported a combined theoretical and spectroscopic study on the cis- and trans-rotamers of 3ClTA. The vibrational spectra of 3ClTA in S1 and D0 states were measured by R2PI and mass-analyzed threshold ionization (MATI) spectroscopy. The stable molecular structures and the rotational barriers between the cis- and trans-rotamers of 3ClTA are calculated by comparing the spectroscopic measurements with the theoretical calculations. The structural changes induced by photoexcitation and photoionization are discussed in detail.

Results

Calculated Results

The chemical structures and the atomic labelings of cis- and trans-3ClTA are shown in Figure . The optimized geometric parameters and atomic charges for the two rotamers in S0, S1, and D0 states are listed in Table , and their Cartesian coordinates are provided in the Supporting Information (Table S1). The calculated one-dimensional potential energy curves (PECs) of 3ClTA in S0, S1, and D0 states are presented in Figure . All PECs were scanned along the C2–C1–S–C8 dihedral angle from −10° to 190° with a 5° increment, and all other structural parameters were optimized without any constraints. The presented energies for the stationary points corresponding to minimum structures include zero-point energy (ZPE) corrections, but ZPE corrections are not included for other points on the PECs.

Figure 1

Chemical structures and atomic labelings of (left) cis- and (right) trans-3ClTA.

Table 1

Calculated Geometric Parameters and Atomic Charges of cis- and trans-Rotamers of 3ClTA

	cis-			trans-
	S0	S1	D0	S0	S1	D0
Bond Length (Å)
S–CH3	1.815	1.804	1.809	1.815	1.806	1.811
S–C1	1.774	1.754	1.717	1.775	1.754	1.716
C1–C2	1.394	1.395	1.406	1.399	1.412	1.417
C2–C3	1.389	1.434	1.377	1.383	1.431	1.375
C3–C4	1.385	1.401	1.415	1.390	1.389	1.409
C4–C5	1.392	1.390	1.393	1.387	1.402	1.400
C5–C6	1.384	1.440	1.378	1.391	1.433	1.375
C6–C1	1.400	1.403	1.427	1.395	1.389	1.419
C3–Cl	1.755	1.746	1.719	1.754	1.761	1.722
Bond Angle (o)
C8-S-C1	103.9	106.9	106.7	103.8	106.7	107.2
S-C1-C2	124.1	122.2	125.1	116.0	112.2	114.8
S-C1-C6	116.5	112.9	114.2	124.7	122.9	124.5
C2-C1-C6	119.4	124.9	120.7	119.4	124.8	120.7
C1-C2-C3	119.2	116.0	118.6	119.5	115.6	119.1
C2-C3-C4	122.1	120.7	121.0	121.8	120.9	120.3
C3-C4-C5	118.2	122.0	120.2	118.2	121.9	120.4
C4-C5-C6	121.0	119.0	119.9	121.2	118.7	120.5
C5-C6-C1	120.2	117.5	119.6	119.9	117.8	110.1
C2-C3-Cl	118.6	118.8	120.3	119.0	118.6	120.5
Atomic Charges (e)
H of −CH3 group	0.208	0.238	0.248	0.207	0.238	0.248
	0.208	0.238	0.248	0.207	0.235	0.248
	0.226	0.233	0.261	0.226	0.233	0.261
C8	–0.705	–0.736	–0.729	–0.704	–0.738	–0.729
S	0.270	0.638	0.668	0.270	0.646	0.681
C1	–0.137	–0.195	–0.170	–0.136	–0.200	–0.169
C2 (Ha)	–0.279	–0.348	–0.241	–0.257	–0.372	–0.199
	(0.218)	(0.208)	(0.240)	(0.222)	(0.212)	(0.248)
C3	0.017	–0.104	0.010	0.014	–0.074	–0.014
C4 (Ha)	–0.251	–0.153	–0.079	–0.250	–0.153	–0.072
	(0.218)	(0.208)	(0.243)	(0.218)	(0.207)	(0.243)
C5 (Ha)	–0.165	–0.281	–0.181	–0.162	–0.301	–0.162
	(0.209)	(0.204)	(0.246)	(0.207)	(0.202)	(0.244)
C6 (Ha)	–0.230	–0.320	–0.117	–0.253	–0.282	–0.161
	(0.210)	(0.203)	(0.239)	(0.205)	(0.192)	(0.229)
Cl	–0.018	–0.031	0.116	–0.015	–0.044	0.104

The H atom in parentheses is directly connected to the C atom before the parentheses.

Figure 2

Calculated one-dimensional potential energy curves (PECs) along the dihedral angle between the −SCH3 group and the benzene ring in S0, S1, and D0 states of 3ClTA. The presented energies for the stationary points corresponding to minimum structures include zero-point energy (ZPE) corrections, but ZPE corrections are not included for intermediate points on the PECs.

The calculations indicate that different chlorine isotopes have negligible effects on the investigated properties of 3ClTA, such as molecular geometry, excitation (E1) and ionization (IE) energies, and vibrational frequencies. Such finding is similar to previous studies on several aromatic molecules containing the chlorine substituent.[1,3,20,23,27,28] Hence, we will focus on the results of the 35Cl-substituted 3ClTA in the following sections for clarity.

The PECs presented in Figure indicate that the cis-rotamer of 3ClTA is slightly more stable than its trans-rotamer in S0 and D0 states, while the calculated energies including ZPE corrections for the stationary points in the S1 state show that the trans-rotamer is slightly more stable than the cis-rotamer. Nonetheless, these energy differences are relatively small, and a more confident assignment needs experimental support. In the S0 state, the calculated energy difference between the cis- and trans-rotamers is 68 cm–1, including ZPE corrections, and the isomerization energy barrier is about 500 cm–1 with respect to the cis-rotamer. A previous study showed that flexible organic molecules with interconversion barriers larger than 400 cm–1 would not relax significantly during the supersonic expansion.[25] Hence, both the cis- and trans-rotamers of 3ClTA could exist in the supersonic molecular beam with an estimated population ratio Ncis/Ntrans of about 1.4 according to the Maxwell–Boltzmann distribution under the assumption of thermal equilibrium.

It is noted that in the calculated PEC for the S1 state (Figure ), the benzene ring is distorted and the Cl atom is optimized to be out-of-plane for most of the data points, except at the points where the C2–C1–S–C8 dihedral angles are 0 and 180°. Similar behavior has been found previously for the S1 state of TA.[9,12] Therefore, caution should be taken for quantitative evaluation of the calculated PEC and rotational barrier of 3ClTA in the S1 state. Nonetheless, at the C2–C1–S–C8 dihedral angles of 0 and 180°, the benzene ring was all optimized to be planar with the Cl atom in the same plane for S0, S1, and D0 states.

A discontinuity at 90° was noted in the PEC of the D0 state. It is mainly caused by the rotation of the −CH3 group, and such discontinuity does not affect the optimized results for both cis- and trans-rotamers at 0 and 180° in the D0 state. Similar discontinuities were observed previously in the PECs of 2-N-methylaminopyridine.[26] In this regard, the calculated IEs and relative stabilities of the two rotamers in the D0 state are assumed to be not affected by the discontinuity at 90°.

1C-R2PI Spectra

A typical time-of-flight (TOF) mass spectrum of 3ClTA recorded at the UV wavelength of ∼285 nm is presented in Figure . The mass peaks at 159 and 161 amu are assigned, respectively, to the ion signals of 35Cl-3ClTA and 37Cl-3ClTA isotopomers, showing that the mass resolution of the setup is enough for collecting the individual spectra of the two isotopomers of 3ClTA. The 1C-R2PI spectra of the 35Cl and 37Cl isotopomers of 3ClTA in the range of 0–1200 cm–1 are shown in Figure . Tentative assignments of the observed vibrational features are listed in Table , and the Varsányi’s labeling system is adopted to approximately describe the observed benzene-like vibrational modes of 3ClTA, which should be classified as an meta-di-heavy substituted benzene derivative.[29]

Figure 3

Typical time-of-flight (TOF) mass spectrum of 3ClTA.

Figure 4

1C-R2PI spectra of 35Cl-3ClTA (top) and 37Cl-3ClTA (bottom). Letters “c” and “t” in parentheses represent the cis-rotamer and trans-rotamer, respectively.

Table 2

Observed Vibrational Frequencies (in cm–1) in the 1C-R2PI Spectra of cis- and trans-3ClTA and Their Tentative Assignments

cis-3ClTA				trans-3ClTA
exptl.	vib.	theo.a	assignmentb	exptl.	vib.	theo.a	assignmentb
33 959	0		band origin	34 102	0		band origin
34 005	46	48	τ01,torsion	34 126	24	33	τ01,torsion
34 056	97	96c	τ02,torsion	34 143	41	66c	τ02,torsion
34 168	209	213	1501	35 065	963	954	101, breathing
34 308	349	360	6a01, β(C–C–C)
34 368	409	408c	6a01τ01
34 576	617	621	νs1, stretching
34 774	815	790	1201, β(C–C–C)
34 912	953	960	101, breathing

The theoretical values have been scaled by a scaling factor of 0.967.

The torsion vibration and stretching vibration are denoted by τ and ν, respectively. The bending mode is designated by β.

The calculated frequency for these overtone and combination modes are taken from the sum of the calculated individual frequency of each mode.

As shown in Figure , the first electronic excitation energies of cis- and trans-35Cl-3ClTA are identical to those of 37Cl-3ClTA, confirming that the Cl isotope substitution has negligible effects on the first electronic excitation energies of 3ClTA. Most of the vibrational frequencies are close to each other for the two 3ClTA isotopomers within 3 cm–1. For example, the 6a01, νs1, and 101 bands locate at 349, 617, and 953 cm–1 for 35Cl-3ClTA, and 351, 615, and 955 cm–1 for 37Cl-3ClTA, respectively. The relative intensities of several vibrational bands are different between the two isotopomers, which is presumably attributed to the edge eliminating effect in choosing the mass range to obtain the reduced spectra during experiments. Those experimentally observed similarities between the two 3ClTA isotopomers verified our predictions derived from theoretical calculations, allowing us to discuss only the spectrum of 35Cl-3ClTA in the following sections.

3ClTA has 42 normal modes, and the intensity of the vibronic bands is proportional to the Franck–Condon factors (FCFs), and thus not all vibrational modes, for instance, in the cases of very small FCFs, could be observed in the 1C-R2PI process. The assignments of the observed bands to a specific rotamer have been confirmed by measuring the respective IEs using the MATI spectroscopic method.

According to the calculations, the cis-rotamer is slightly more stable than the trans-rotamer in the S0 state and should have a higher population in the supersonic molecular beam and thus corresponds to a stronger band origin peak in the 1C-R2PI spectrum. Hence, the observed vibrational band at 33959 cm–1 is tentatively assigned to the cis-rotamer’s band origin (S1 ← S0; 000 band). This band origin corresponds to E1. The band at 34102 cm–1 is tentatively assigned to the trans-rotamer’s 000 band. The experimentally measured 000 band of the trans-rotamer is blue-shifted by 143 cm–1 with respect to that of the cis-rotamer. These assignments are also consistent with the calculated results that the cis-rotamer is more stable than the trans-rotamer in the D0 state. However, these results seem not to be consistent with the calculations in the S1 state. Such discrepancy could be attributed to the expected difficulty of TDDFT in accurately predicting the electronically excited states.

Besides the 000 bands, the bands at 46 and 97 cm–1 are assigned as the methyl torsion vibrational mode τ01 and its overtone τ02 of the cis-rotamer, in good agreement with the calculated methyl torsion vibration of 48 cm–1 and its overtone of 96 cm–1 (Table ; also see Figure S1 in the Supporting Information for the vibrational vectors); the band at 209 cm–1 is assigned as the cis-1501 mode, and the corresponding calculated value is 213 cm–1; the bands at 349, 617, 815, and 953 cm–1 are assigned to the in-plane ring deformation modes 6a01, νs1, 1201, and 101 of the cis-rotamer, also agreeing with the calculated values of 360, 621, 790, and 960 cm–1; and a combination band of 6a01 and τ01 is observed at 409 cm–1 (Figure ; the calculated 6a01 + τ01 sum frequency is 408 cm–1). As for the trans-rotamer, the observed bands at 24 and 41 cm–1 are assigned to the methyl torsion mode τ01 and its overtone τ02, and the band at 963 cm–1 is assigned to the in-plane ring deformation mode 101. Those measured values are well reproduced by the calculations, i.e., 33 and 954 cm–1 for τ01 and 101, respectively. The above mode assignments are summarized in Table , and in Figure S1 in the Supporting Information, their vibrational vectors are provided.

MATI Spectra

The MATI spectra and the assignments of the vibrational bands observed for the cis- and trans-rotamers of 35Cl- and 37Cl-3ClTA are presented in this section. For the cis-rotamer, MATI spectra were recorded via the vibrational peaks at 0, 46, 349, and 409 cm–1, which show relatively high intensities in the 1C-R2PI spectra; for the trans-rotamer, MATI spectra via the band origin 000 at 0 cm–1 and the τ01 mode at 24 cm–1 were collected. Except for several vibrational peaks that were not observed due to the low intensity, the MATI spectra of 37Cl-3ClTA are nearly identical to those of 35Cl-3ClTA for both rotamers (Figure ). All vibrational frequencies of the two isotopomers were observed within 7 cm–1, close to the experimental uncertainty (±5 cm–1). Thus, only the MATI spectra of 35Cl-3ClTA are discussed here.

Figure 5

MATI spectra of cis-35Cl-3ClTA (top) and 37Cl-3ClTA (bottom) recorded via the 000 vibrational mode in the S1 state.

As shown in Figures and 6, a strong prominent band origin in the D0 state was observed for both cis- and trans-rotamers via the 000 vibrational mode in the S1 state, suggesting that the molecular shapes of both rotamers do not change significantly upon photoionization from the S1 to D0 state. The structural change has been quantitatively evaluated by calculating the root-mean-square deviation (RMSD) for distances, and it can be seen that the RMSD values between S1 and D0 are 0.029 and 0.123 Å for cis- and trans-rotamers, respectively (Figure S2). The band origins were measured to be 65 323 and 65 468 cm–1 for the cis- and trans-rotamers, respectively. After correcting the direct-current Stark shift, the IEs of cis- and trans-rotamers were determined, respectively, to be 65 326±5 and 65 471±5 cm–1, giving a measured difference of IE values to be 145 cm–1. These values agree reasonably with the calculated respective values of 63 205, 63 456, and 251 cm–1 including ZPE corrections (Figure ).

Figure 6

MATI spectra of trans-35Cl-3ClTA (top) and 37Cl-3ClTA (bottom) recorded via the 000 vibrational mode in the S1 state.

The assignments of all of the vibrational modes observed in the MATI spectra are listed in Table . Most of the observed active vibrations are found to be related to the ring deformation and substituent-sensitive modes. For the cis-35Cl-3ClTA, the bands at 88 and 177 cm–1 are assigned as the methyl torsion vibrational mode τ1 and its overtone τ2, respectively. Note that the calculated methyl torsion vibrational frequency is 83 cm–1. The bands at 223, 652, and 980 cm–1 are assigned as the vibrational modes of 151, νs1, and 11, respectively, agreeing with the calculated respective values of 270, 629, and 978 cm–1; the bands at 289 and 446 cm–1 are assigned as combination modes. The assignment of the prominent band at 398 cm–1 in the MATI spectra via the 000 mode (Figures and 6) is uncertain. It could be due to either the 6a1 or 16b1 mode. The 6a1 mode is observed to be prominent in the 1C-R2PI spectra (Figure ) with a vibrational frequency of 349 cm–1. This observation seems to support the assignment of the calculated 381 cm–1 band to the 6a1 mode. However, the vibrational frequency of the 6a1 mode is measured to be ∼366 cm–1 instead of 398 cm–1 in the MATI spectra when the τ01 and 6a01 modes in the S1 state were used as the intermediate levels (Figure a,b). Similar phenomena have been observed previously for 2FTA.[11] For example, the 6a1 band in the D0 state of 2FTA was observed at 411 cm–1via the 000 mode in the S1 state, while it was measured to be 400 cm–1 when the τ026a01 mode was used as the intermediate level. For the trans-35Cl-3ClTA, the frequencies of most benzene-ring-related vibrational modes are observed to be close to those of the cis-35Cl-3ClTA. The most prominent peak in the MATI spectrum of trans-35Cl-3ClTA via the methyl torsion mode τ01 at 24 cm–1 (Figure e) is tentatively assigned as β(C-SCH3) + β(C–Cl), in reasonable agreement with the calculated value of 154 cm–1 (Table ; see also Figure S1 in the Supporting Information for the vibrational vectors). The vibrational frequency for overtone τ3 should be also close to ∼120 cm–1, while activation of τ3 is usually highly unfavored.

Table 3

Observed Vibrational Frequencies (in cm–1) in the MATI Spectra of the cis- and trans-Rotamers of 3ClTA and Their Tentative Assignments

cis-							trans-
exptl.							exptl.
via 000	via 46 cm–1	via 349 cm–1	via 409 cm–1	via 617 cm–1	theo.a	assignmentb	via 000	via 24 cm–1	theo.a	assignmentb
0	0	0	0	0			0	0
88					83	τ1	34	33	70	τ1
177	178				166c	τ2		81		τ2
								127	154	β(C-SCH3) + β(C–Cl) or τ3
223	225				270	151	177	186	182	151
289	310				353c	151τ1	223	210	252c	151τ1
398	367	366			381	6a1 or 16b1	399		400	6a1 or 16b1
446	446		441		464c	6a1τ1
652				653	629	νs1	655		627	νs1
980					978	11	981		968	11

The theoretical values have been scaled by a scaling factor of 0.967.

The denotations are according to the 1C-R2PI spectra in Table , and the torsion vibration and stretching vibration are denoted by τ and ν, respectively. The bending mode is designated by β.

The calculated frequency for these overtone and combination modes are taken from the sum of the calculated individual frequency of each mode.

Figure 7

MATI spectra of cis- and trans-35Cl-3ClTA recorded via different vibrational modes in the S1 state.

In Figure , MATI spectra via several other vibrational bands are presented. It can be seen that when the MATI spectra were recorded via the in-plane ring deformation modes 6a01, νs1, and the combination mode 6a01τ01 as the intermediate levels, the most intensive vibrational peaks observed in the MATI spectra correspond to the same modes of the intermediate levels (Figure b,d). This indicates that the molecular geometry and the corresponding vibrational coordinates of the 6a01, νs1, and 6a01τ01 modes do not change much upon photoionization from the S1 state to the D0 state, as verified by the RMSD analysis (Figure S2). Similar phenomena via the benzene-ring-related vibrational levels were previously observed for anisole and its derivatives.[1,17,21,23,30] However, when the methyl torsion mode τ01 was used as the intermediate level, the MATI spectra showed dramatically different characteristics as those via the benzene-ring-related vibrational levels. For the cis-rotamer, the most prominent peak was due to the 0+ band, and many other vibrational modes were also significantly excited (Figure a). The MATI spectrum of the trans-rotamer via the methyl torsion mode τ01 at 24 cm–1 is completely different from that of the cis-rotamer, the 0+ band was much weaker, and several low-frequency vibrational modes were strongly excited (Figure e). These observations indicate that the methyl torsion mode vibrations are strongly coupled to many other vibrational modes for both the cis- and trans-rotamers during the photoionization process, but their detailed coupling mechanisms should be quite different due to the different orientations of the −SCH3 groups in the two rotamers (see the torsion mode vibration vectors in Figure S1).

Discussion

Molecular Structures of cis- and trans-3ClTA in S0, S1, and D0 States

Molecular Structures in S0 State

Both cis- and trans-3ClTA are stable and adopt planar structures possessing a Cs point group in the S0 state. The cis-rotamer is slightly more stable than the trans-rotamer, suggesting that the inductive effect of the Cl atom dominates over the steric effect between the two substituent groups in determining the relative stability of the two rotamers. As seen in Table , the S–C8 bond in the −SCH3 group is 1.815 Å in the cis-rotamer and trans-rotamer, and ∠SC1C2 is 124.1° in the cis-rotamer and 116.0° in the trans-rotamer. The slightly larger ∠SC1C2 in the cis-rotamer may arise from the larger steric repelling effect between the Cl atom and the −SCH3 group in the cis-rotamer than that in the trans-rotamer. This is similar to the larger ∠NC1C2 in cis-3-Cl-N-methylaniline than that in the trans-rotamer.[16]

As for the benzene ring, all of the calculated bond length differences between the two rotamers are less than 0.007 Å and all the bond angle differences are less than 0.4°, indicating that the influences on the structure of the benzene ring caused by the different orientations of the substituent groups in the two rotamers are small.

From the perspective of intramolecular charge distributions, the two rotamers are close to each other. Table shows that their calculated atomic charges in the S0 state are almost identical, except for the C2 and C6 atoms. The C2 atom is more negative in the cis-rotamer (−0.279 e for the cis-rotamer and −0.257 e for the trans-rotamer), while the C6 atom is more negative in the trans-rotamer (−0.230 e for the cis-rotamer and −0.253 e for the trans-rotamer). This could be caused by the electron-donating effect of the −CH3 group, and the C atom closer to the −CH3 group gets more negative charge.

Structural Changes during the S1←S0 Transition

The intense 000 bands observed in the 1C-R2PI spectra in Figure suggest that there are only minor geometric changes upon excitation from the S0 state to the S1 state, as verified by the RMSD analysis (Figure S2). According to the calculated results shown in Table , the C1–C2, C2–C3, C3–C4, C5–C6, and C6–C1 bonds in the benzene ring become a little longer upon excitation while the C4–C5 bond stays about the same for the cis-rotamer; for the trans-rotamer, the C1–C2, C2–C3, C4–C5, and C5–C6 bonds become longer during the S1 ← S0 transition, while the C3–C4 and C6–C1 bonds are contracted by 0.001 and 0.006 Å, respectively. The above observations are consistent with the previous studies, which showed that the S1 ← S0 transitions of benzene derivatives are mainly subject to the π* ← π electronic excitation and a benzene ring distortion usually occurs, making it deviate from a perfect hexagon.[31,32] Such transition scheme is confirmed by analyzing the molecular orbitals involved in the first electronic transition (Figure S3).

Besides the benzene ring distortion, the C3–Cl bond connecting the benzene ring and the Cl atom is shortened by 0.009 Å for the cis-rotamer during the S1 ← S0 transition, while it is elongated slightly by 0.007 Å for the trans-rotamer. This may indicate that the steric effect between the Cl atom and the −SCH3 group becomes weaker and/or the interaction between the Cl atom and the benzene ring becomes stronger in the S1 state for the cis-rotamer. The −SCH3 group is also distorted during the S1 ← S0 transition. The S–C1 bond connecting the −SCH3 group with the benzene ring is shortened by 0.020 Å for the cis-rotamer and 0.021 Å for the trans-rotamer. The S–C8 bond connecting the S atom with the −CH3 group is shortened by 0.011 Å for the cis-rotamer and 0.009 Å for the trans-rotamer. This is consistent with the change of the intramolecular charge distribution upon photoexcitation. As shown in Table , all H atoms in the −CH3 group and the S atom become more positive upon the excitation, while all H atoms directly connected to the benzene ring, and all of the C atoms of the benzene ring become more negative in the S1 state, except for C4, which becomes less negative and is furthest from the −SCH3 group. The shortened bond lengths and the intramolecular charge distribution changes in the S1 state enhance the p−π conjugation interaction between the −SCH3 group and the benzene ring (Figure S3). Some anisole derivatives, such as 3-F anisole[21] and 3-Cl anisole,[23] undergo similar ring deformation upon excitation from the S0 state to the S1 state.

Structural Changes during Photoionization

Both the measured spectra and theoretical calculations suggest that both of the 3ClTA rotamers adopt a stable planar structure in the D0 state. The calculated results show that for both rotamers, the benzene ring of 3ClTA has longer C1–C2, C3–C4, C4–C5, and C6–C1 bonds and shorter C2–C3 and C5–C6 bonds in the D0 state than those in the S0 state, showing a quinoidal structure. The chlorine atom and all the C atoms of the benzene ring except for C1 that connects with the −SCH3 group become more positive (or less negative) in the D0 state than those in the S1 state upon photoionization; while the changes of the charges on C1 and the −SCH3 group are comparably smaller. This implies that the ejected photoelectron is mainly contributed by the nonbonding p orbitals of the S and Cl atoms and the π orbital of the benzene ring (see the plotted highest molecular orbitals in Figure S3).

For the −SCH3 group, the S–C1 bond shows a very strong double bond character in the D0 state, and it is significantly shorter (1.718 Å for the cis-rotamer and 1.716 Å for the trans-rotamer) than that in S0 (1.774 Å for the cis-rotamer and 1.775 Å for the trans-rotamer) and S1 states (1.754 Å for the cis-rotamer and 1.754 Å for the trans-rotamer). This is consistent with the spectroscopic measurements that the frequencies of the methyl torsion mode τ1 for both rotamers become significantly larger in the D0 state (88 cm–1 for the cis-rotamer and 34 cm–1 for the trans-rotamer) than those in the S1 state (46 cm–1 for the cis-rotamer and 24 cm–1 for the trans-rotamer). The above observations indicate that the conjugation interaction between the p orbital occupied by the remaining nonbonding electron on the S atom and the π orbital on the benzene ring is enhanced in the photoionization process.

For the Cl atom, the C3–Cl bond is also significantly shorter in the D0 state (1.719 Å for the cis-rotamer and 1.722 Å for the trans-rotamer) than that in S0 (1.755 Å for the cis-rotamer and 1.754 Å for the trans-rotamer) and S1 states (1.746 Å for the cis-rotamer and 1.761 Å for the trans-rotamer). The Cl atom can interact with the benzene ring through the p-π conjugative effect, which arises from the donation of a lone pair electron from the Cl atom to the benzene ring, thus the conjugative effect between the Cl atom and the benzene ring may also be enhanced due to the shortened C3–Cl bond in the D0 state and the rigidity of the benzene ring could be increased. This is in accordance with the spectroscopic observations. For example, the frequencies of the in-plane ring deformation mode 11 in the D0 state (980 cm–1 for the cis-rotamer and 981 cm–1 for the trans-rotamer) are slightly larger than those in the S1 state (953 cm–1 for the cis-rotamer and 963 cm–1 for the trans-rotamer), indicating a more rigid benzene ring in the D0 state.

Active Vibrational Modes and Their Frequencies

The present 1C-R2PI and MATI spectroscopic measurements show that the in-plane ring deformation modes 6a01, νs1, and 101 of the cis-rotamer have frequencies of 349, 617, and 953 cm–1 in the S1 state, and, respectively, become 366 (or 398), 652, and 980 cm–1 in the D0 state; the in-plane ring deformation mode 101 of the trans-rotamer has a frequency of 963 cm–1 in the S1 state, and becomes 981 cm–1 in the D0 state. Thus, almost all of the benzene ring-related vibrational modes have higher frequencies for both the cis- and trans-rotamers during the photoionization process, indicating that the conjugation effects between the benzene ring and the substituent groups are enhanced in the D0 state. This increases the rigidity of the benzene ring. Similar vibrational frequency enhancement by photoionization was observed for the substituent-sensitive modes, for example, the methyl torsion mode τ01, as presented above.

The current study along with several previous studies on other anisole derivatives[21,23,27,30] and TA derivatives[6,11] show that the frequencies of the benzene-ring-related modes are not very sensitive to the orientations of the substituents. In this study, the frequencies of the substituent-sensitive modes, e.g., the methyl torsion mode τ01, are measured for both the cis- and trans-3ClTA. Its frequency for the cis-rotamer is 46 cm–1 in the S1 state, and 88 cm–1 in the D0 state, while it is only 24 and 34 cm–1 in the S1 and D0 states of the trans-rotamer. The MATI spectra via the methyl torsion mode τ01 for the two rotamers as shown in Figure a,e show that the coupling scheme of the methyl torsion mode τ01 to other vibrational modes is strongly rotamer-dependent. This may be caused by the different through-space interactions between the two substituents Cl and −SCH3 in the two rotamers. The two substituents are closer to each other, and thus should have stronger through-space interaction in the cis-rotamer than that in the trans-rotamer. That makes the methyl group in the cis-rotamer harder to move, resulting in a larger vibrational frequency.

It should also be noted that the methyl torsion mode τ01 has been observed to be active in TA,[6,9,12]trans-2FTA,[11] and cis- and trans-3ClTA, while it was not observed for stable anisole derivatives, such as 2-F anisole,[22] 3-Cl anisole,[23] and 3-F anisole.[21] This indicates that the methyl group on the more flexible −SCH3 group is easier to be activated.

Substitution Effect on the Electronic Excitation and Ionization Energies

The measured electronic excitation transition energies (E1) and ionization energies (IE) of several anisole and TA derivatives are listed in Table . It can be seen that the E1 and IE values depend on both the nature and location of the substituent halogen atom.

Table 4

Transition and Ionization Energies (in cm–1) of Anisole and TA Derivatives

	E1	ΔE1	IE	ΔIE		E1	ΔE1	IE	ΔIE
thioanisole (TA)a	34 506		63 906		anisole	36 383d	0	66 399e
trans-o-fluorothioanisoleb	34 974	468	65 114	1208	trans-o-fluoroanisolef	36 611	228	67 354	955
cis-m-fluorothioanisolec	34 820	314	65 468	1562	cis-m-fluoroanisoleg	36 662	279	67 867	1468
trans-m-fluorothioanisolec	35 047	541	64 644	1738	trans-m-fluoroanisoleg	36 819	436	68 304	1905
cis-m-chlorothioanisole	33 959	–547	65 326	1420	cis-m-chloroanisolei	35 822	–561	67 645	1246
	(33 580)h		(63 205)h
trans-m-chlorothioanisole	34 102	–404	65 471	1565	trans-m-chloroanisolei	35 868	–515	68 008	1609
	(33 434)h		(63 456)h

From ref (6).

From ref (11).

From our unpublished data.

From ref (2)

From ref (44).

From ref (17).

From ref (21).

Values in parentheses represent theoretical values calculated in this work.

From ref (23).[23]

For E1 values shown in Table , the E1 values of the Cl-substituted TA derivatives, cis- and trans-3ClTA, are red-shifted by 547 and 404 cm–1, respectively, relative to that of TA,[6,9,12] while the E1 values of the F-substituted TA derivatives, the trans-2FTA,[11] and cis- and trans-3FTA, are all blue-shifted compared to that of TA.[6,9,12] Similar phenomena have been observed between anisole and its F- and Cl-substituted derivatives.[2,17,21,23] A previous study showed that the substitution of a functional group on the benzene ring can lower its zero-point energy (ZPE) level.[21] The red-shifted E1 values of 3ClTA imply that the interaction of Cl with the benzene ring is stronger in the S1 state than that in the S0 state, resulting in a larger ZPE lowering of the S1 state than that of the S0 state.

Conclusions

The 1C-R2PI and MATI spectra of cis- and trans-3ClTA are obtained, implying that both the cis- and trans-rotamers coexist in the experiment and have planar structures in their respective S0, S1, and D0 states. The cis-rotamer is slightly more stable than the trans-rotamer with an isomerization barrier of ∼500 cm–1. The first electronic excitation and ionization energies of cis-3ClTA are determined to be 33 959±3 and 65 326±5 cm–1, and those of trans-3ClTA to be 34 102±3 and 65 471±5 cm–1, respectively. The two Cl isotopes are found to have negligible effects on the molecular properties of 3ClTA. The E1 and IE values are red-shifted by the Cl substitution for both rotamers, similar to the anisole derivatives.

Most of the active vibrations are found to be the in-plane ring deformation and substituent-sensitive modes. Their vibrational frequencies are generally measured to be higher in the D0 state than those in the S1 state, indicating that the conjugation interaction between the benzene ring and the substituted functional groups and the through-space interaction between the two substituents are enhanced during the photoionization process. This is confirmed by the molecular structural changes due to photoionization as revealed by theoretical calculations. The frequencies of the in-plane ring deformation modes are found to be not sensitive to the relative orientation between the two substituents, while the frequencies of the substituent-sensitive modes and their coupling schemes with other vibrational modes are found to be strongly rotamer-dependent.

Experimental and Computational Methods

Experimental Methods

The photoionization time-of-flight (TOF) mass spectrometer used in the current study has been described previously,[5,20] and only a brief description is given here. The 3ClTA sample with a purity of 99% was purchased from Alfa Aesar and used without further purification. It was seeded in the Ar carrier gas (∼2 atm) and then expanded into a vacuum chamber through a pulsed General valve with a nozzle diameter of 0.25 mm. The molecular beam passed through a skimmer with a diameter of 1 mm located at 20 mm downstream from the nozzle and then entered the photoionization region of the mass spectrometer, which is 70 mm from the nozzle. In the photoionization region, the molecular beam was crossed by one or two ultraviolet (UV) laser beams perpendicularly and photoionized for detection.

In the one color resonant two-photon ionization (1C-R2PI) spectroscopic measurement, only one UV laser beam was used. The 3ClTA molecule absorbed a single UV photon to be excited to the S1 state first, and then the 3ClTA molecule in the S1 state absorbed a second UV photon to go above the photoionization threshold. The ions thus produced by the UV laser were immediately accelerated by two direct-current electric fields of 200 and 2100 V/cm, respectively. In the MATI experiment, two counterpropagating UV laser beams of different wavelengths were used. The first UV laser beam was the same as that used in the 1C-R2PI experiment and was used to excite the 3ClTA molecule to its S1 state; the second UV beam excited the molecules from the S1 state to the high-n-Rydberg states, which are slightly below the ionization thresholds. The photoionization region was field-free during the photoexcitation process, and after about 200 ns, a pulsed electric field of −0.5 V/cm was switched on to reject the prompt ions, which were usually generated together with the high-n-Rydberg neutrals. After a delay of 11 μs, two pulsed electric fields of +200 and +2000 V/cm were switched on synchronously to field ionize the high-n-Rydberg neutrals and accelerate the produced ions to the detector. The ions produced in the above two processes were focused by an Einzel lens and flew through a 1.0 m long field-free tube toward a dual-stacked microchannel plate (MCP) detector. After passing through a preamplifier (SR445, Stanford Research System), the TOF signal was amplified by 25 times and collected by a multichannel scaler (SR430, Stanford Research System).

The UV lasers used in this experiment were generated by doubling the outputs of two tunable dye lasers (Sirah-CSTR), which were pumped by the second harmonic output of an Nd:YAG laser (Quanta-Ray, Pro-230-10E) at a repetition rate of 10 Hz. Three laser dyes, Pyrromethene 597, Pyrromethene 580, and DCM, were used in this work. Two pulse delay generators (DG535, Stanford Research System) were used to synchronize the whole system. Typical pressures of ∼3.0 × 10–3 and ∼3.0 × 10–5 Pa were maintained in the source and ionization chambers during operation.

Computational Methods

All calculations in this study were carried out with the Gaussian09 program package.[33] The stable structures and harmonic vibrational frequencies of 3ClTA in S0, S1, and D0 states were calculated with density functional theory (DFT), in which time-dependent DFT (TDDFT)[34,35] was used for the excited state S1. The B3LYP functional[36,37] and cc-pVTZ basis set[38−40] were employed for all the calculations, and no imaginary frequencies were found for the optimized stationary points. The relaxed potential energy curve (PEC) in the S0 state was scanned to illustrate the isomerization barrier between the cis- and trans-3ClTA. The calculated vibrational frequencies of 3ClTA in S0, S1, and D0 states were scaled by 0.967[41] to approximately account for the vibrational anharmonicity effect. The charge distributions based on natural population analysis (NPA)[42] were conducted by Natural Bond Orbital (NBO) version 3.1 program implemented in the Gaussian09 program package.[43]