Jun Miyazaki*ab,
Hiroshi Takiyamab and
Munetaka Nakata*c
aFaculty of Pharmaceutical Sciences, Hokuriku University, Ho-3, Kanagawa-machi, Kanazawa, Ishikawa 920-1181, Japan. E-mail: j-miyazaki@hokuriku-u.ac.jp
bDepartment of Chemical Engineering, Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
cGraduate School of BASE (Bio-Applications and Systems Engineering), Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan. E-mail: necom816@cc.tuat.ac.jp
First published on 17th January 2017
UV-induced photoreactions of thiazole isolated in low-temperature argon matrices have been investigated by a joint use of infrared spectroscopy and density-functional-theory calculations. Photoproducts have been identified by comparison of the observed infrared spectra with the corresponding calculated spectral patterns, leading to the conclusion that undetected open-chain molecules, syn-2-isocyanoethenethiol (CN–CHCH–SH) and 2-isocyanothiirane, are initially produced by cleavage of the CS–CN bond with hydrogen-atom migration, when the matrix samples are exposed to UV radiation coming from a super high-pressure mercury lamp for 3 min. In the secondary photolysis, syn-2-isocyanoethenethiol and 2-isocyanothiirane change to another unknown molecule, 2-isocyanoethanethial (CN–CH2–CHS), by hydrogen-atom migration with generation of the CS double bond. These photoreaction pathways are supported by kinetic analysis of the absorbance changes of IR bands against irradiation time. We have also found that HCN and the ˙CHCH–S˙ biradical are photodecomposed from thiazole by cleavage of the CN–CC bond following the cleavage of the CS–CN bond, where the hydrogen atom on the center carbon atom of ˙CHCH–S˙ immediately migrates to the end carbon atom to form CH2CS or to the sulfur atom to form HCC–SH. In addition, weak bands of the species of interest in astrophysics and astrochemistry such as HCCH, NC–SH, HNCS, HCNS, and the ˙CN radical are detected, but the photoconversion from thiazole to isothiazole or Dewar thiazole is not found. The ring-opening photoreaction, photoisomerization and photodecomposition pathways of thiazole isolated in low-temperature argon matrices are discussed comprehensively.
In astrochemical and astrophysical studies, UV photochemistry of heterocyclic molecules, which have yet to be detected in the interstellar medium (ISM), gives important information on small photofragments composed of H, C, N, O, and S because they are exposed to the strong UV radiation. In addition, the heterocyclic compounds are interested in astrobiology as the origins of S and N involved in biological materials.4 Indeed, 2-aminothiazole5 and 2-aminooxazole6 have been recently considered as important nucleotide precursors, implying that azole derivatives have potentials of prebiotic molecules. Since they have attractive electronically excited states due to calcogens (O, S, and Se) and N, a various photochemical reaction pathways are possible, implying that the photochemical reaction mechanism of azoles derivatives is more complicated than that of aromatic hydrocarbons. Although photochemical behavior of some five-membered ring containing two hetero atoms has been investigated,7 the whole reaction pathways have not been completed yet. In the present study, we focus our attention on the photoreaction pathways of thiazole.
Several research groups have investigated the photochemistry of thiazole (1) and isothiazole (1′) experimentally and theoretically. For example, Cateau et al. reported that the photoconversion from (1′) to (1) proceeds but not the reverse photoconversion from (1) to (1′) in any solvents.8 By contrast, the final photodecomposition products of thiazole (1) in the gas phase were found to be cyano radical (˙CN) and thiocyanato radical (˙NCS) using a flash photolysis technique,9 although the photodecomposition pathways have not been elucidated yet.
The low-temperature rare-gas matrix-isolation experiment is one of the most powerful techniques to study photoreaction pathways, because unstable or noble species are immediately frozen and trapped in low-temperature solid, which is chemically inert and no absorption in UV and visible-light regions, and stabilized in the environment without interactions between species and solvents. For example, the structure of a thiazole/carbon suboxide (C3O2) complex was investigated in argon matrices,10 while the phosphorescence of thiazole in Ne, Ar, Xe, SF6, and CH4 matrices at 4 K was measured.11 Most recently, the conformational changes of a thiazole derivative (thiazole-2-carboxylic acid) isolated in Ar and N2 matrices have been reported by Halasa et al., using narrowband excitation with near-infrared (IR) and UV light.12 However, the photoreaction pathways of matrix-isolated thiazole, the parent molecule with no substitution groups, have never been reported yet because of their complication, even though thiazole is one of the most fundamental heterocyclic compounds having two hetero atoms.
In the present study, the photochemical reactivity and stability of thiazole isolated in solid argon matrices are investigated by a joint use of IR spectroscopy and density-functional-theory (DFT) calculations. Our purposes are to identify intermediates and final products in the photolysis of thiazole as many as possible, to examine absorbance changes of IR bands against irradiation time, to elucidate when and where and why ring-opening reaction, photoisomerization, and photodecomposition occur, and finally to propose the whole UV-induced photoreaction pathways of thiazole.
IR spectra were measured with an FT-IR spectrometer (JEOL, JIR-WINSPEC50) equipped with an MCT detector cooled by liquid N2. All the spectra were measured at 0.5 cm−1 resolution and averaged over 64 scans. A super high-pressure mercury lamp (SHPML) (USHIO, BA-H500, λ > 200 nm) was used to induce photoreactions. The IR region in the light source was filtered out using a water filter, the wavelength of the light source was controlled with/without optical glass filters.
DFT calculations were performed to obtain the optimized geometry, relative energies and IR spectral patterns of reactants, intermediates, and products using the GAUSSIAN03 program15 at the UB3LYP/aug-cc-pVTZ level among some calculation levels, where their electronic ground state is singlet. Other calculation conditions are default values in the program. Scaling factors of 0.96, 0.97 and 0.98 were applied to the regions over 2800 cm−1, between 2800 and 1900 cm−1, and below 1900 cm−1, respectively, in the present study so as to reproduce the observed wavenumbers for bands of the reactant, thiazole. TD-DFT calculations were performed to obtain the vertical transition energy and oscillator strength at the same level.
A difference IR spectrum between the spectra measured before and after irradiation for only 3 min is shown in Fig. 3a, where the negative and positive bands indicate the reactant and the product, respectively. For example, the intensity of the 1484 and 1043 cm−1 bands of thiazole (1) decreased, and new bands appeared at 2113, 1591, and 1013 cm−1. The intense 2113 cm−1 band is characteristic of species containing an isocyano (CN–) or cyano (NC–) group.20 There are a few candidates of C3H3NS species including the triple bond, which have the same chemical formula as the parent molecule, thiazole (1). 2-Isocyanoethenethiol (CN–CHCH–SH) (2) is one of the candidates for the photoproducts of thiazole (1), which is produced by cleavage of the S1–C2 bond with the hydrogen-atom migration from C2 to S1 (see Scheme 1). The existence of (2) was predicted previously by the DFT calculation (B3LYP/6-31G*),21 but neither spectral parameters nor experimental evidences have yet been reported. Thus we have performed the DFT calculation for the geometrical optimization and the vibrational analysis of (2) to assign the bands newly appearing in Fig. 3.
Two conformations are possible around both the C–S single bond and the CC double bond of (2), resulting in the four planar conformers shown in Fig. 4. We call anti and syn for the former conformation and (E) and (Z) for the latter conformation. The optimized geometrical parameters and the relative energies are shown in Fig. 4. Syn-(Z)-2 is the most stable one among the four conformers, and anti-(Z)-2 is less stable than syn-(Z)-2 by 4.6 kJ mol−1, while both syn-(E)-2 and anti-(E)-2 are less stable than syn-(Z)-2 by ∼8 kJ mol−1. The calculated wavenumbers and IR intensities of the bands for the four conformers are listed in Table S2 of the ESI.†
The calculated spectral pattern of the most stable syn-(Z)-2 is shown in Fig. 3b to compare with the observed spectrum of the photoproduct. The most intense band at 2113 cm−1 is assignable to the CN– stretching mode of syn-(Z)-2, while the 1591 and 1013 cm−1 bands are assignable to the CC stretching and S–H bending modes, respectively. Since the wavenumbers predicted by the DFT calculation are consistent with the corresponding observed values within ∼10 cm−1, as compared in Table 1, we conclude that syn-(Z)-2 is the initial photoproduct of thiazole (1) in argon matrices.
Molecules | Observed | Calculateda | Assignment | ||
---|---|---|---|---|---|
ν/cm−1 | Intensityb | ν/cm−1 | Intensityb | ||
a Calculated at the DFT/UB3LYP/aug-cc-pVTZ level. Scaling factors of 0.96, 0.97, and 0.98 are applied to the regions over 2800 cm−1, between 2800 to 1900 cm−1, and below 1900 cm−1, respectively.b Relative intensity is normalized to the most intense band.c Overlapped with a band of (4).d Overlapped with a band of ethenethione.e Overlapped with a band of syn-3. | |||||
syn-(Z)-2-Isocyaonoethenethiol(syn-(Z)-2) | 1013 | 29.4 | 1003.99 | 23.70 | S–H bending |
1591 | 7.7 | 1599.32 | 15.69 | CC stretching | |
2113 | 100.0 | 2115.16 | 100.00 | CN– stretching | |
syn-2-Isocyanoethanethial(syn-3) | 627c | 624.94 | 6.47 | C–C–N bending | |
700d | 704.46 | 3.93 | CH2 rocking + C–H bending | ||
804 | 2.4 | 782.45 | 2.95 | C–C stretching | |
970 | 1.4 | 971.92 | 6.89 | C–H bending | |
998 | 9.4 | 992.26 | 5.77 | C–N stretching | |
1139 | 8.1 | 1134.31 | 19.49 | CS stretching | |
1140 | 5.4 | ||||
1302 | 3.2 | 1297.50 | 11.01 | C–H bending | |
1303 | 11.0 | ||||
1304 | 7.2 | ||||
1368 | 12.6 | 1375.91 | 18.99 | C–H bending | |
1425 | 4.1 | 1425.35 | 11.45 | CH2 scissoring | |
2157 | 100.0 | 2163.30 | 100.00 | CN– stretching | |
2915 | 3.7 | 2889.28 | 5.20 | C–H stretching | |
2971 | 2.2 | 2952.58 | 12.37 | C–H stretching | |
2972 | 2.7 | ||||
anti-2-Isocyanoethanethial(anti-3) | 700d | 704.37 | 4.95 | CH2 rocking + C–H bending | |
936 | 4.8 | 942.31 | 11.72 | C–N stretching | |
954 | 24.8 | 947.40 | 7.12 | C–H bending | |
1135 | 10.2 | 1117.97 | 22.74 | CS stretching | |
1285 | 10.6 | 1287.05 | 26.39 | C–H bending | |
2151 | 100.0 | 2155.71 | 100.00 | CN– stretching | |
2975 | 1.3 | 2976.14 | 2.64 | C–H stretching | |
2-Isocyanothiirane (4) | 627e | 625.79 | 17.72 | Ring stretching | |
855 | 5.4 | 840.54 | 4.49 | CH2 rocking | |
987 | 18.2 | 977.62 | 8.72 | C–H bending | |
1348 | 16.3 | 1347.44 | 12.73 | C–H bending | |
2132 | 100.0 | 2138.63 | 100.00 | CN– stretching | |
2133 | 68.9 |
IR bands of syn-3 and anti-3 in other spectral regions are very weak after UV irradiation for 45 min. Thus we tried to find the bands of syn-3 and anti-3 in the difference spectrum between the spectra measured after 180 min minus after 45 min UV irradiation (see Fig. 6). It is noted that the bands due to syn-(Z)-2 initially produced from thiazole (1) are negligible in this difference spectrum because their intensities are almost unchanged during this prolonged irradiation period, as described in Section 3.2.4. Several weak bands of syn-3 and anti-3 are detectable in Fig. 6 and marked with “syn-3” and “anti-3”, respectively. The observed wavenumbers and relative intensities of the bands of syn-3 (12 vibrational modes) and anti-3 (7 vibrational modes) are consistent with the corresponding calculated values, as shown in Table 1. Therefore, we conclude that both syn-3 and anti-3 are produced from thiazole (1) upon prolonged UV irradiation. The photoproduction pathways for syn-3 and anti-3 are discussed in Section 3.2.5 using the kinetic analysis of the IR absorbance changes against the irradiation time.
Initially, we considered that 2-isocyanothiirane (4) was produced from thiazole (1) via Dewar thiazole, because Dewar structures of five-membered ring molecules such as furan,22 thiophene,22 and cyclopentadiene24 were previously detected as intermediates in photolysis using the matrix-isolation technique. In addition, Dewar thiazole derivatives were theoretically proposed as candidates of the photochemical intermediates for 2-phenylthiazole and 2-acetylthiazole.25 However, no IR bands of Dewar thiazole are found in our IR spectra around the wavenumbers predicted by our DFT calculations (see Table S4, ESI†). For example, the most intense band due to the C–H in-plane bending mode is predicted at 1262.75 cm−1, but undetectable. The relative energy of Dewar thiazole is estimated to be higher than those of 2-isocyanothiirane (4) and thiazole (1) by 94.5 and 258.8 kJ mol−1, respectively, implying that Dewar thiazole is so unstable under our experimental condition to be detected. The photoproduction pathway shown in Scheme 2, where (4) is directly produced from thiazole (1), is supported by the kinetic analysis of the absorbance change against the irradiation time, as explained in the following section.
The photochemical behavior of thiazole derivatives such as methylthiazole and phenylthiazole in solvents were investigated theoretically25,26 and experimentally.27,28 It is reported that the photoconversion from the thiazole derivatives to the corresponding isothiazole derivatives occurs, where Dewar structures play a role of intermediates. If the photoconversion from thiazole (1) to isothiazole (1′) is accelerated in low-temperature argon matrices, we could detect the IR bands of isothiazole (1′). However, no corresponding bands near the 726, 819, 1239, and 1391 cm−1 bands reported in the gas phase29 were detected in the difference spectrum shown in Fig. 6. This fact also supports our conclusion that Dewar thiazole is not produced from thiazole (1) under our experimental condition.
The absorbance of syn-(Z)-2 slightly decreased upon prolonged irradiation over 45 min, as shown in Fig. 7, indicating that syn-(Z)-2 changes to other species by the secondary photo-induced reactions. On the other hand, the absorbance of syn-3 and anti-3 increases more slowly than that of syn-(Z)-2 and clearly reveals induction periods at the early irradiation stage. Thus we assume that syn-3 and anti-3 are photoproduced from syn-(Z)-2 by hydrogen-atom migration from S1 to C4 with generation of the CS double bond (see Scheme 3). If the SH group in the photoexcited state of syn-(Z)-2 is perpendicular to the other part (CN–CH), it is reasonable to assume that both syn-3 and anti-3 are produced from syn-(Z)-2.
In contrast to the fact that the absorbance of thiazole (1) continues to decrease upon prolonged irradiation over 45 min, the absorbance of (4) is nearly constant as shown in Fig. 7, indicating that a partial amount of (4) changes to other species by the secondary photo-induced reactions like syn-(Z)-2. It is noted that the absorbance of the band for syn-3 increases more strongly over 45 min irradiation than that for anti-3. Thus we assume that syn-3 is also photoproduced from (4) by hydrogen-atom migration from C5 to C4 with generation of the CS double bond (see Scheme 4). To explain this selective photoproduction of syn-3 but not anti-3, we assume that the isocyano group in (4) moves to the opposite side against the hydrogen atom migrating from the methylene group. In this case, the attraction between the sulfur atom and the nitrogen atom in the photoexcited state of (4) may play an important role for the selectivity, which is an interesting problem especially for theoretical chemists. We support the selective photoreaction pathway by the kinetic analysis of the absorbance changes described in the following section.
(1) → (syn-(Z)-2) + (4), rate constant: k1 | (1) |
(syn-(Z)-2) → syn-3 + anti-3, rate constant: k2 | (2) |
(4) → syn-3, rate constant: k3 | (3) |
Five rate equations are derived from the above reaction pathways of (1)–(3) as follows,
d[(1)]/dt = −k1[(1)], | (4) |
d[syn-(Z)-2]/dt = k1[(1)] − k2[syn-(Z)-2], | (5) |
d[(4)]/dt = k1[(1)] − k3[(4)], | (6) |
d[syn-3]/dt = k2[syn-(Z)-2] + k3[(4)], | (7) |
d[anti-3]/dt = k2[(syn-(Z)-2)], | (8) |
A(1) = A0exp(–k1t), | (9) |
A(syn-(Z)-2) = (ε(syn-(Z)-2)/ε(1))A0{αexp(–k1t) − αexp(–k2t)}, | (10) |
A(4) = (ε(4)/ε(1))A0{βexp(–k1t) − βexp(–k3t)}, | (11) |
A(syn-3) = (ε(syn-3)/ε(1))A0{2–(k2/k1)αexp(–k1t) − (k3/k1)βexp(–k1t) + αexp(–k2t) + βexp(–k3t)}, | (12) |
A(anti-3) = (ε(anti-3)/ε(1))A0[1 − (k2/k1)αexp(–k1t) + αexp(–k2t)], | (13) |
We have performed the least-squares fitting for the absorbance changes using the eqn (9)–(13) to determine three rate constants of k1, k2, and k3 and four ratios of absorption coefficients such as ε(4)/ε(1). The obtained rate constants are 0.00236 ± 0.00001, 0.112 ± 0.010, and 0.0228 ± 0.0040 min−1 for k1, k2, and k3, respectively, and the obtained ratios of absorption coefficients are 2.23 ± 0.18, 0.190 ± 0.028, 0.506 ± 0.013, and 0.536 ± 0.007 for ε(syn-(Z)-2)/ε(1), ε(4)/ε(1), ε(syn-3)/ε(1), and ε(anti-3)/ε(1), respectively, where the uncertainty represents one standard deviation. The rate constant k2 is five times larger than k3, meaning that the reactivity of syn-(Z)-2 is higher than that of (4). The ratio of absorption coefficients for syn-3 and anti-3 is calculated to be 0.506/0.536 = 0.944, which is almost consistent with the corresponding value obtained by the DFT calculation, 174.5/183.6 = 0.950 (see Table S3, ESI†). The calculated absorbance for each species is drawn by solid lines in Fig. 7, which reproduces the observed values satisfactorily, strongly supporting our proposed reaction pathways shown in Schemes 1–4.
Vertical transition energy/nm | Oscillator strength | Obsd. wavelengthb/nm | |
---|---|---|---|
a Calculated at the TD-DFT/UB3LYP/aug-cc-pVTZ level. Confromation of biradicals is assumed to be unchanged from thiazole (1). Transition energy with the largest value of oscillator coefficient is listed.b In solvents. Ref. 14. | |||
Thiazole | 221.54 | 0.0715 | ∼235 |
199.38 | 0.0362 | ∼205 | |
(S1–C2 cleavage) ˙CHN–CHCH–S˙ | 293.15 | 0.1669 | |
(C5–S1 cleavage) ˙CHCH–NCH–S˙ | 251.48 | 0.1600 | |
(N3–C4 cleavage) ˙NCH–S–CHCH˙ | 195.14 | 0.1314 |
Several papers reported the bands of HCN in the gas phase30 and in the low-temperature matrices23,31–34 to be ∼3303 and ∼720 cm−1 for the C–H stretching and bending modes, respectively. The 3304 cm−1 and 715 cm−1 bands observed in Fig. 6 are consistent with the values of HCN (3306 and 721 cm−1) produced by the photolysis of s-tetrazine and s-triazine in argon matrices.35 Thus we conclude that HCN is produced by the cleavage of the N3–C4 following the cleavage of the S1–C2 bond according to Scheme 5.
There is an intense band at 1752 cm−1 marked with “TK” in Fig. 6. This band is characteristic of the CCS stretching mode of thioketene derivatives. For example, the bands of trimethylenethioketene and cyclopentylidenthioketene were observed at 1793 and 1790 cm−1 in argon matrices at 10 K,36,37 respectively, while the corresponding bands of methylthioketene (CH3–CHCS) and vinylthioketene (CH2CH–CHCS) were observed at 1777 and 1740 cm−1 in argon matrices at 12 K.38,39 Since the wavenumber of the 1752 cm−1 band observed in Fig. 6 is close to that of the CCS stretching mode of CH2CS in argon matrices,40 1755 cm−1, we conclude that CH2CS is produced in the photolysis of thiazole (1). The C–H stretching band of CH2CS, which is reported to be 3010 cm−1 in argon matrices, is also detected at 3011 cm−1. Six bands marked with “TK” in Fig. 6 are assignable to CH2CS. The observed and reported wavenumbers and the relative intensities are summarized in Table 3.
Obsd. ν/cm−1 | Intensitya | Ar matrixb | Assignment | |
---|---|---|---|---|
a Relative intensity is normalized to the most intense band.b Ref. 40. Symbols of vs, s, m, w, and vw denote very strong, strong, medium, weak, and very weak, respectively.c Overlapped with bands of syn-3 and anti-3.d Undetectable.e Ref. 32.f Ref. 35. | ||||
Ethenethione (CH2CS) | ||||
700c | 692 | s | CH2 wagging | |
1330 | 11.5 | 1322 | m | |
1410 | 2.1 | 1410 | w | |
1752 | 100.0 | 1755 | vs | CCS stretching |
1843 | 3.4 | 1840 | m | |
3011 | 20.9 | 3010 | m | C–H stretching |
Ethynethiol (HCC–SH) | ||||
d | 558 | w | C–H bending | |
966 | 10.3 | 959 | w | S–H bending |
d | 1112 | m | ||
2062 | 19.5 | 2065 | w | CC stretching |
2532 | 10.9 | 2575 | vw | S–H stretching |
3313 | 100.0 | 3315 | vs | C–H stretching |
HCC–SH was previously identified as a photoproduct of 1,2,3-thiadiazole,40 where the characteristic IR bands due to the C–H, S–H, and CC stretching modes were observed at 3315, 2575, and 2065 cm−1, respectively. The 3313 cm−1 band marked with “ET” in Fig. 6 is consistent with the reported 3315 cm−1 band. In addition, the bands due to the CC stretching and S–H bending modes are detected at 2062 and 966 cm−1, respectively (see Table 3), while the S–H stretching band is detected at 2532 cm−1. Thus we conclude that HCC–SH is co-produced from thiazole (1) with HCN by the cleavage of the N3–C4 bond following the cleavage of the S1–C2 bond like CH2CS. The detection of the IR bands for HCN, CH2CS and HCC–SH supports the photodecomposition pathway shown in Scheme 5.
It is known that HCC–SH is converted to CH2CS by hydrogen-atom migration upon UV irradiation.40 However, the absorbance of the bands for both CH2CS and HCC–SH increase linearly with no induction periods at the early irradiation state, as shown in Fig. 8. Thus we assume that CH2CS and HCC–SH are independently produced from thiazole (1) by cleavage of the N3–C4 bond following the cleavage of the S1–C2 bond with the hydrogen-atom migration according to Scheme 5. One may claim that HCN, CH2CS, and HCC–SH are produced by cleavage of the S1–C2 bond following the cleavage of the N3–C4 bond. However, the vertical transition energy of ˙NCH–S–CHCH˙, which is produced by cleavage of the N3–C4 bond, is calculated to be 195.14 nm with the oscillator strength of 0.1314, as listed in Table 2. Since the intensity of the light source of our SHPML with shorter wavelength than 200 nm is extremely weak, the secondary photolysis of ˙NCH–S–CHCH˙ is difficult in contrast to the hydrogen-atom migration in the primary photolysis described in Section 3.4.3.
Fig. 8 Absorbance changes against irradiation time for HCN, CH2CS, and HCCSH. The absorbance is normalized using the absorbance of the band at 862 cm−1 for thiazole (1) before UV irradiation in Fig. 7. |
Fig. 9 Enlarged difference spectra (3450–3200, 1790–1720, and 1290–1250 cm−1) of Fig. 6. (a) Observed difference spectra. (b) Calculated spectral patterns of anti-2-isocyanoethanethial (anti-3), (methyleneamino)-ethenethione (anti-5 and syn-5), and syn-N-ethynylthioformamide (syn-6) at the DFT/UB3LYP/aug-cc-pVTZ level. Scaling factors of 0.96 and 0.98 are applied to regions over 2800 cm−1 and below 1900 cm−1, respectively. |
This species (5) has two conformations around the N3–C4 single bond (syn-5 and anti-5). The optimized geometrical parameters and the relative energies are shown in Fig. 4. All the calculated wavenumbers and IR intensities of both anti-5 and syn-5 are listed in Table S5 of the ESI.† The conformation of (5) photoproduced initially could be syn-5 but not anti-5, because the molecular shape of syn-5 is more similar to that of the reactant, thiazole (1), than that of anti-5. However, the DFT calculation reveals that syn-5 is less stable than anti-5 by 10.04 kJ mol−1, probably because of the strong repulsion between π electrons on the CN and CCS double bonds like 1,3-butadiene.41,42 Note that the optimized structure of syn-5 is non-planar and the dihedral angle around CH2N–CHC is calculated to be 21.7° (see Fig. 4). The weak doublet band detected at 1760 and 1763 cm−1 in Fig. 9 is found to be consistent with the calculated value of the CCS stretching mode of anti-5, 1761.49 cm−1, but not with that of syn-5, 1725.47 cm−1. In addition, the C–H bending mode of anti-5 is detected at 1278 cm−1 in Fig. 9. Thus we conclude that anti-5 is produced from thiazole (1) via syn-5 initially photoproduced in an argon-matrix cage according to Scheme 6. The observed and calculated wavenumbers and the relative intensities of anti-5 are compared in Table 4.
Molecules | Observed | Calculateda | Assignment | ||
---|---|---|---|---|---|
ν/cm−1 | Intensityb | ν/cm−1 | Intensityb | ||
a Calculated at the DFT/UB3LYP/aug-cc-pVTZ level. Scaling factors of 0.96, 0.97, and 0.98 are applied to the regions over 2800 cm−1, between 2800 and 1900 cm−1, and below 1900 cm−1, respectively.b Relative intensity is normalized to the most intense band.c Overlapped with bands of the reactant. | |||||
anti-(Methyleneamino)-ethenethione (anti-5) | 1278 | 90.5 | 1276.71 | 14.77 | C–H bending |
1760 | 100.0 | 1761.49 | 100.00 | CC stretching | |
1763 | 59.0 | ||||
syn-N-Ethynylthioformamide (syn-6) | c | 1042.22 | 16.27 | C–H bending | |
1252 | 35.0 | 1255.41 | 77.36 | C–H bending | |
c | 1389.64 | 100.00 | C–H bending | ||
2167 | 19.8 | 2178.52 | 19.28 | CC stretching | |
2169 | 73.1 | ||||
2989 | 0.9 | 2976.37 | 10.17 | C–H stretching | |
3326 | 36.7 | 3335.53 | 42.28 | CC–H stretching | |
3329 | 60.0 | ||||
3431 | 100.0 | 3418.55 | 24.45 | N–H stretching | |
3436 | 51.1 |
Fig. 10 Enlarged difference spectra (2300–2160 and 2100–1900 cm−1) of Fig. 6. The spectral region between 2160 and 2100 cm−1 is omitted because the intensity of the IR bands of syn-3 and anti-3 is too strong. Bands marked with “a”, “b”, “c”, “d”, “e”, and “f” are tentatively assigned to 2-cyanoethenethiol (NC–CHCH–SH), 2-cyanothiirane, NC–SH, HNCS, ˙CN radical, and HCNS, respectively. |
We could not identify the other conformer, anti-6, although the relative energy of anti-6 estimated to be lower than that of syn-6 by 9.20 kJ mol−1. For example, neither the most intense band around 1250 cm−1 nor the second intense band around 1500 cm−1 is detectable in Fig. 6. We have performed the one-step optimization between anti-6 and syn-6, resulting in the barrier height of ∼80 kJ mol−1 at the dihedral angles of 105 and 285° (see Fig. S1, ESI†). Thus we assume that the barrier height is so high that the conversion from syn-6 initially produced from thiazole (1) to the more stable conformer, anti-6, is unable to occur in an argon-matrix cage.
It is noted that the absorbance of the bands for anti-5 and syn-6 is much weaker than that of anti-3, syn-3, HCN, CH2CS, and HCC–SH, as shown in Fig. 6. Thus we conclude that the photoreaction pathways to form anti-5 by the cleavage of the S1–C2 bond shown in Scheme 6 and to form syn-6 by the cleavage of the C5–S1 bond shown in Scheme 7 are minor ring-opening reactions.
2-Cyanothiirane could be caused by detachment and recombination of the CN group in 2-isocianothiirane (4) like the photoconversion from CN–CHCH–SH to NC–CHCH–SH. The calculated wavenumbers and IR intensities of 2-cyanothiirane are shown in Table S4 of the ESI.† The NC– stretching band is predicted at 2276 cm−1, and the corresponding band is detected at 2275 cm−1, marked with “b” in Fig. 10. The most intense band due to the C–S symmetric stretching mode, 629 cm−1, and the thirdly intense band due to the CH2 waging mode, 1065 cm−1, are also detected at 658 and 1079 cm−1 in Fig. 6, respectively. Thus we conclude that a small amount of syn-(Z)-2 and (4), which are initially produced from thiazole (1), convert to the corresponding cyano compounds upon prolonged UV irradiation.
A lot of papers have reported the wavenumbers of the vibrational modes for HCCH in low-temperature matrices46–50 and in the gas phase.51 We detected the 3238, 3234 and 3231 cm−1 triplet band in Fig. 6, which is in good agreement with the reported value for HCCH, 3240 cm−1. The NC– stretching mode of NC–SH52 and the NCS stretching mode of HNCS53 are reported to be 2182 cm−1 and 1981.8 cm−1 in low-temperature matrices, respectively. We detected the corresponding bands at 2186 cm−1 marked with “c” and at 1987 cm−1 marked with “d” in Fig. 10, respectively. The detection of HCCH, NC–SH, and HNCS supports the photodecomposition pathway shown in Scheme 9, although it is a minor reaction pathway.
If the hydrogen atom of ˙NCH–S˙ detaches but migrates neither to the sulfur atom nor to the nitrogen atom, ˙NCS radical could be co-produced with HCCH. The band of ˙NCS radical is reported to be 1942.2 cm−1 in free jet experiment.54,55 We tried to detect the corresponding band, but there is no band in the NCS stretching region between 1970 and 1900 cm−1 in Fig. 10. Instead of the band due to ˙NCS radical, we found a band at 2046 cm−1 marked with “e” in Fig. 10, which is consistent with the band of ˙CN radical produced from HCN by vacuum UV light in argon matrices at 14 K, 2046 cm−1.31
Venkatasubramanian and Krishnamachari claimed that ˙NCS radical changes to ˙CN radical in a flash photolysis experiment of thiazole (1).9 Note that our difference spectra were obtained by experiments without any short-cutoff glass filters, meaning that all the radiation coming from the SHPML is used to induce photoreactions. Thus one possibility for the detection of ˙CN radical is that the photodecomposition from ˙NCS radical to S and ˙CN radical immediately occurs even when ˙NCS radical is produced by detachment of the hydrogen atom from ˙NCH–S˙.
However, we propose another possibility that ˙CN radical is detached from the isocyano compounds, syn-(Z)-2 and (4). As explained in Section 3.5, a small amount of the cyano compounds are photoproduced from the corresponding isocyano compounds by detachment and recombination of ˙CN radical. If ˙CN radical detached from the isocyano compounds stays in a low-temperature argon-matrix cage without recombination, the IR band of ˙CN radical is detectable.
On the other hand, another minor ring-opening photoreaction is caused by the cleavage of the S1–C2 bond upon prolonged UV irradiation. The hydrogen-atom migration from C5 to C2 in ˙CHN–CHCH–S˙ forms anti-5 via syn-5 initially produced in an argon-matrix cage, as shown in the green area of Scheme 10. Other minor ring-opening photoreaction is caused by the cleavage of the C5–S1 bond upon prolonged UV irradiation to produce ˙CHCH–NCH–S˙ biradical, from which syn-6 is produced by the hydrogen-atom migration from C4 to N3, where the syn conformation is kept in an argon-matrix cage without photoconversion to anti-6. The photoproducts of (2), (3), (4), (5) and (6) are identified for the first time in the present study by a joint use of IR spectroscopy and DFT calculations.
When the N3–C4 bond cleaves following the cleavage of the C5–S1 bond, HCCH and ˙NCH–S˙ biradical are photodecomposed from ˙CHCH–NCH–S˙. The hydrogen atom of ˙NCH–S˙ immediately migrates to the nitrogen atom to form HNCS or to the sulfur atom to form NC–SH, as shown in the blue area of Scheme 10. In addition, several weak bands in the difference spectra measured after prolonged UV irradiation are tentatively assigned to cyano compounds caused by photoconversion from isocyano compounds, ˙CN radical caused by detachment from isocyano compounds, and HCNS caused by photoisomerization from photofragments such as NC–SH.
In the present study, we investigated the photochemical reactivity and stability of thiazole isolated in solid argon matrices by a joint use of IR spectroscopy and DFT calculations. By examination of absorbance changes of IR bands against irradiation time, we identified intermediates and final products. Some of them are unknown isocyano compounds, and some of them are species interested in astrochemistry and astrophysics. We expect that the present results will contribute to develop these research fields and help to detect other simple but unknown species.
Footnote |
† Electronic supplementary information (ESI) available: Tables listing observed wavenumbers and IR intensities of thiazole (1) isolated in solid argon matrices with references, and calculated wavenumbers and IR intensities of 2-isocyanoethenethiol (2), 2-isocyanoethenethial (3), 2-isocyanothiirane (4), (methyleneamino)-ethenethione (5), N-ethynylthiformamide (6), Dewar thiazole, and 2-cyanothiirane obtained at the DFT/UB3LYP/aug-cc-pVTZ level, and a figure showing calculated potential energy around the C–NH–CHS dihedral angle of N-ethynylthioformamide (6)obtained by the one-step optimization at an interval of 15°. See DOI: 10.1039/c6ra27496j |
This journal is © The Royal Society of Chemistry 2017 |