Conformational Structure, Infrared Spectra and Light-Induced Transformations of Thymol Isolated in Noble Gas Cryomatrices

: The conformational space of the natural product thymol (2-isopropyl-5-methylphenol) was investigated using quantum chemical calculations at the B3LYP and MP2 levels, which revealed the existence of four types of conformers differing in the orientation of the isopropyl and hydroxyl groups. Thymol monomers were isolated in noble gas (Ar and Xe) matrices (at 15 K) and characterized by IR spectroscopy. With the support of B3LYP harmonic vibrational calculations, the two most stable trans-OH-conformers, differing in the isopropyl orientation, were identiﬁed in the cryomatrices. The two less stable cis-OH conformers were not detected as they shall undergo fast tunneling to the most stable ones. Annealing experiments in a Xe matrix up to 75 K did not lead to any conversion between the two isolated conformers, which is in accordance with the signiﬁcative energy barrier computed for rotamerization of the bulky isopropyl group (~24 kJ mol − 1 ). Vibrational excitation promoted by broadband or by narrowband irradiation, at the 2 ν (OH) frequencies of the isolated conformers, did not lead to any conversion either, which was interpreted in terms of a more efﬁcient energy transfer to the hydroxyl rotamerization (associated with a lower energy barrier and a light H-atom) than to the isopropyl rotamerization coordinate. Broadband UV irradiation experiments ( λ > 200 nm) led to a prompt transformation of matrix isolated thymol, with spectroscopic evidence suggesting the formation of isomeric alkyl-substituted cyclohexadienones, Dewar isomers and open-chain conjugated ketenes. The photochemical mechanism interpretation concords with that reported for analogous phenol derivatives.

Despite the wide range of bioactivities, the number of studies on the structure of thymol available in literature is relatively scarce [9][10][11].X-ray data show that in the crystalline state, the thymol molecules are held together through O-H•••O hydrogen bonds forming hexamers, which in turn are linked to each other by van der Waals contacts, while the crystal conformation is characterized by having the O-H group pointing away from the isopropyl fragment and the tertiary C-H bond of this fragment pointing towards the hydroxyl oxygen [9].In gas phase, combination of microwave spectroscopy with theoretical calculations led to the experimental identification of three conformers for thymol, with one of them being much more populated than the other two [10].
Despite the wide range of bioactivities, the number of studies on the structure of thymol available in literature is relatively scarce [9][10][11].X-ray data show that in the crystalline state, the thymol molecules are held together through O−H•••O hydrogen bonds forming hexamers, which in turn are linked to each other by van der Waals contacts, while the crystal conformation is characterized by having the O−H group pointing away from the isopropyl fragment and the tertiary C−H bond of this fragment pointing towards the hydroxyl oxygen [9].In gas phase, combination of microwave spectroscopy with theoretical calculations led to the experimental identification of three conformers for thymol, with one of them being much more populated than the other two [10].
Some infrared (IR) spectroscopic studies have been performed for thymol in the solid state [12,13] and in solution [14].However, as far as we are aware, IR spectra of the monomeric compound have never been reported, nor has its light-induced transformations.Therefore, one of the objectives of the present paper is to fill this gap and to present, for the first time, a detailed analysis of the IR spectra of thymol isolated in lowtemperature Ar and Xe matrices, and of the transformations triggered by exposing the matrix-isolated compound to IR-and UV-radiation.
Thymol is a positional isomer of carvacrol (5-isopropyl-2-methylphenol, see Scheme 1).The two molecules only differ in the relative position of the OH and isopropyl groups.In thymol, they are in ortho-position, while in carvacrol, they are in meta-position.Very recently, we have studied the conformational, vibrational and photochemical properties of carvacrol isolated in a low-temperature Ar matrix [15].In that study, it has been shown that the molecule exists as two conformers in the cryogenic environment, differing from each other by the orientation of the isopropyl fragment.Exposition of the matrix-isolated carvacrol to broadband IR radiation was found to be effective in changing the relative population of the two conformers, by internal rotation of the bulky isopropyl fragment.A more specific goal of the current work is to verify whether this channel of vibrational energy dissipation is still effective in thymol, where the OH and isopropyl groups are closer together, making the barrier for internal rotation of the isopropyl group higher than in carvacrol.

Experimental Methods
Solid thymol (melting point of 49-51 °C) was supplied by Fluka with a purity degree of 99.5%.The matrix host gases (Ar N60 or Xe N48) were provided by Air Liquide.The low-temperature matrices were prepared by placing a small amount of thymol crystals into a glass tube, which was then connected to a cryostat (APD Cryogenics, Allentown, PA, US, closed-cycle He refrigerator, with a DE-202A expander) through a needle valve, kept at ~298 K.The vapor of the sublimating compound was co-deposited with an excess of the matrix gas onto a CsI window cooled to 15 K.The temperature of the CsI window was measured directly at the sample holder by using a silicon diode temperature sensor connected to a digital controller (Scientific Instruments, Allentown, PA, US, model 9650-1), providing the stabilization of temperature with an accuracy of 0.1 K. IR spectra of the Some infrared (IR) spectroscopic studies have been performed for thymol in the solid state [12,13] and in solution [14].However, as far as we are aware, IR spectra of the monomeric compound have never been reported, nor has its light-induced transformations.Therefore, one of the objectives of the present paper is to fill this gap and to present, for the first time, a detailed analysis of the IR spectra of thymol isolated in low-temperature Ar and Xe matrices, and of the transformations triggered by exposing the matrix-isolated compound to IR-and UV-radiation.
Thymol is a positional isomer of carvacrol (5-isopropyl-2-methylphenol, see Scheme 1).The two molecules only differ in the relative position of the OH and isopropyl groups.In thymol, they are in ortho-position, while in carvacrol, they are in meta-position.Very recently, we have studied the conformational, vibrational and photochemical properties of carvacrol isolated in a low-temperature Ar matrix [15].In that study, it has been shown that the molecule exists as two conformers in the cryogenic environment, differing from each other by the orientation of the isopropyl fragment.Exposition of the matrix-isolated carvacrol to broadband IR radiation was found to be effective in changing the relative population of the two conformers, by internal rotation of the bulky isopropyl fragment.A more specific goal of the current work is to verify whether this channel of vibrational energy dissipation is still effective in thymol, where the OH and isopropyl groups are closer together, making the barrier for internal rotation of the isopropyl group higher than in carvacrol.

Experimental Methods
Solid thymol (melting point of 49-51 • C) was supplied by Fluka with a purity degree of 99.5%.The matrix host gases (Ar N60 or Xe N48) were provided by Air Liquide.The low-temperature matrices were prepared by placing a small amount of thymol crystals into a glass tube, which was then connected to a cryostat (APD Cryogenics, Allentown, PA, USA, closed-cycle He refrigerator, with a DE-202A expander) through a needle valve, kept at ~298 K.The vapor of the sublimating compound was co-deposited with an excess of the matrix gas onto a CsI window cooled to 15 K.The temperature of the CsI window was measured directly at the sample holder by using a silicon diode temperature sensor connected to a digital controller (Scientific Instruments, Allentown, PA, USA, model 9650-1), providing the stabilization of temperature with an accuracy of 0.1 K. IR spectra of the matrixisolated compound were recorded in the mid-IR range (4000-400 cm −1 ) with 0.5 cm −1 resolution using a Nicolet 6700 Fourier-transform infrared (FTIR) spectrometer (Porto, Portugal), equipped with a deuterated triglycine sulphate (DTGS) detector (Porto, Portugal) and a KBr beam splitter (Porto, Portugal).For the experiments carried out in xenon, spectra in the near-IR (NIR) range (7500-4000 cm −1 , 1 cm −1 resolution) were also recorded using the same spectrometer equipped with an Indium Gallium Arsenide (InGaAs) detector and a CaF 2 beamsplitter.
The thymol/Xe matrix was irradiated through the outer quartz window of the cryostat by using narrowband frequency-tunable NIR light provided by the idler beam of a Quanta-Ray MOPO-SL optical parametric oscillator (Mountain View, CA, USA) (pulse energy 10 mJ, duration 10 ns, repetition rate 10 Hz) pumped with a pulsed Nd:YAG laser.Alternatively, the sample was irradiated with broadband IR light provided by a kanthal wire electrically heated to a reddish-orange glow.Broadband UV-irradiations were also performed on the thymol/Ar matrix by using a 200 W output from a medium pressure Xe/Hg lamp (Oriel, Newport, UK).

Computational Methods
The thymol molecule has two conformational degrees of freedom, corresponding to the H-C7-C2-C1 (α) and H-O11-C1-C2 (β) dihedral angles (see Scheme 1).For identification of the conformers exhibited by this molecule, a two-dimensional (2D) potential energy surface (PES) was computed by increasing both dihedral angles in steps of 20 • over the 0-360 • range, in all combinations, while all the remaining internal coordinates were optimized by using the B3LYP [16][17][18] hybrid DFT method in conjunction with the Pople-type 6-311++G(d,p) basis set.Then, the equilibrium geometries of the thymol conformers, as well as of the putative photoproducts resulting from the UV-photolysis of the matrix-isolated compound, were fully optimized followed by harmonic vibrational calculations at the same DFT level.For the thymol conformers, calculations were also carried out by combining the second-order Møller-Plesset (MP2) perturbation method [19] with the aug-cc-pVDZ Dunning's correlation consistent basis set.Single-point energy calculations were further carried out on the MP2 optimized geometries using the QCISD method [20] combined with the same aug-cc-pVDZ basis set.The Cartesian coordinates of the B3LYP and MP2 optimized geometries of the thymol conformers are provided in Tables S1 and S2.
All quantum-mechanical computations reported in this work were performed with the Gaussian 09 software package [23], while the vibrational assignments were carried out with the vibAnalysis software [24].

Conformers and Barriers to Internal Interconversion
The B3LYP/6-311++G(d,p) 2D-PES calculated for thymol, as a function of the H-C7-C2-C1 (α) and H-O11-C1-C2 (β) dihedrals, is represented in Figure 1 as a contour map.Two of the four minima found on the PES exist as pairs of mirror-image structures lying within relatively broad and shallow potential energy valleys along the α coordinate, extending from −40 to +40 • .One of those doubly degenerate minima is found at (α, β) = (±35 • , 0 • ) and is designated as gc (g + c/g − c), while the other is located at (α, β) = (±35 • , 180 • ) and is labeled as gt (g + t/g − t), where g ± , c and t are abbreviations of gauche plus/minus, cis and trans orientations, respectively.The first and second letters refer to the orientation around the α and β dihedrals, respectively.The g + c/g − c and g + t/g − t pairs of mirror-image structures are separated by first-order saddle points characterized by α = 0 • .Two additional unique minima are found on the PES at (α, β) = (180 • , 180 • ) and (180 • , 0 • ) and are respectively called tt and tc.In Figure 2 are represented the geometries of the four identified thymol conformers, while in Table 1 are given their electronic (∆E el ), zero-point corrected (∆E 0 ) and Gibbs energies (∆G) at 298.15 K (temperature of the gaseous compound before the matrix deposition), calculated with the B3LYP/6-311++G(d,p) and MP2/aug-cc-pVDZ model chemistries on geometries fully optimized at the respective levels.In Table 1 are included the QCISD/aug-cc-pVDZ electronic energies calculated for the MP2/aug-cc-pVDZ optimized geometries.The results of our computations are concordant with those published earlier by Schmitz et al. [10].
structures are separated by first-order saddle points characterized by α = 0°.Two additional unique minima are found on the PES at (α, β) = (180°, 180°) and (180°, 0°) and are respectively called tt and tc.In Figure 2 are represented the geometries of the four identified thymol conformers, while in Table 1 are given their electronic (ΔEel), zero-point corrected (ΔE0) and Gibbs energies (ΔG) at 298.15 K (temperature of the gaseous compound before the matrix deposition), calculated with the B3LYP/6-311++G(d,p) and MP2/aug-cc-pVDZ model chemistries on geometries fully optimized at the respective levels.In Table 1 are included the QCISD/aug-cc-pVDZ electronic energies calculated for the MP2/aug-cc-pVDZ optimized geometries.The results of our computations are concordant with those published earlier by Schmitz et al. [10].Both angles were incrementally fixed in steps of 20°, in all combinations, while all the remaining parameters were optimized at the B3LYP/6-311++G(d,p) level of theory.The minima are indicated in the plot as gt (g + t/g − t), gc (g + c/g − c), tt and tc (geometries are shown in Figure 2).The color bar on the top represents the relative energy scale (zero is assigned to the energy of the global minimum, gt).The isoenergy lines are traced with steps of 1 kJ mol -1 ., in all combinations, while all the remaining parameters were optimized at the B3LYP/6-311++G(d,p) level of theory.The minima are indicated in the plot as gt (g + t/g − t), gc (g + c/g − c), tt and tc (geometries are shown in Figure 2).The color bar on the top represents the relative energy scale (zero is assigned to the energy of the global minimum, gt).The isoenergy lines are traced with steps of 1 kJ mol -1 .
Photochem 2022, 2, FOR PEER REVIEW 5 Figure 2. Geometries of the thymol conformers, including some geometrical parameters calculated at the MP2/aug-cc-pVDZ.For the doubly generated gt and gc conformers, only one of the mirrorimage structures (g + t and g + c) is shown.At all levels of theory, gt is the lowest-energy conformer.It has the hydroxyl group pointing away from of the isopropyl fragment, while the C7-H bond is ± gauche with respect to the C2-C1 bond (see Figure 2), with α taking values of ±35.6 • (B3LYP) or ±39.0 • (MP2).The deviation of the C7-H bond from the aromatic ring plane should be attributed to the presence of the hydroxyl group in ortho-position relative to the isopropyl fragment, a feature also observed in a structural microwave spectroscopic study of the 2,6-diisopropylphenol analog [25].Note that in carvacrol or p-isopropylphenol, where the OH and isopropyl substituents occupy metaor para-positions, respectively, all conformers identified for these molecules were found to have the tertiary C-H bond of the isopropyl group lying in the aromatic ring plane [15,26,27].The geometry of gt is not too different from that adopted by the thymol molecules in the crystal, although in the latter case the C7-H bond is less deviated from the ring plane (α = ±13.3• ) [9].The second most stable conformer, tt, also has the hydroxyl group pointing away from the isopropyl fragment, but the C7-H and C2-C1 bonds are trans to each other, bringing the two methyl groups of the isopropyl close to the oxygen atom (Figure 2).This conformational arrangement results in an increase of the electronic, zero-point corrected, or Gibbs energies by 1.1-2.5 kJ mol −1 relative to the global minimum.
As shown in Figure 2, the (C7 in tt are smaller than the sum of the van der Waals radii of the oxygen and hydrogen atoms (2.72 Å).This suggests the formation of C-H•••O interactions in both trans-OH-conformers, as has been claimed in previous studies carried out for thymol [10] and 2-isopropylphenol [28].These interactions, however, should be very weak because of the significant deviation of the C-H•••O bonding angles (96-116 • ) from the linearity: it is well known that the strongest nonbonded interactions are observed for the linear arrangements (optimal value of 180 • ) [29,30].In the gc and tc conformers, the hydroxyl group is directed towards the isopropyl fragment.Such an orientation does not allow the formation of stabilizing C-H•••O interactions.Instead, it originates one (gc) or two (tc) repulsive CH•••HO close contacts (see Figure 2).Both factors explain the lower stability of these cis-OH-conformers, and in particular of tc, relative to the gt and tt ones.
To characterize the conformational composition of the gaseous sample of thymol immediately before the matrix deposition, the populations of the thymol conformers in the gaseous phase were estimated from the values of ∆G at 298.15 K, through the application of the Boltzmann statistics, and the obtained values have been included in Table 1.According to these estimates, gt is by far the most abundant conformer in the gaseous phase, with a population of 57-67%.Conformer tt is also predicted to be present in the gaseous phase, but in a smaller amount (15-20%), while the population of conformer gc is of only 9% when it is estimated from the B3LYP Gibbs energies, increasing to 24% when obtained from the MP2 Gibbs energies.Nevertheless, as will be shown below, the difference of population predicted for this conformer by the two theoretical methods is not relevant from the experimental viewpoint.Finally, a very small the gaseous phase abundance (3-4%) is predicted for the less-stable tc conformer.
For the interpretation of the results of the matrix-isolation experiments, it is also relevant to have an idea of the energy barriers corresponding to the conversion of the higher into the lower energy conformers [31][32][33].In Figure 3 are represented the OH-torsional barriers interconnecting the pair of conformers (gc and gt), and the pair of conformers (tc and tt).The heights of these barriers in the direction of conformational relaxation were calculated to be 10-12 kJ mol −1 .These values are of the same order as those computed for the conversion of higher into lower energy OH-rotamers in various compounds containing the phenol moiety, such as 2-cyanophenol (15 kJ mol −1 ) [34], 2-isocyanophenol (12 kJ mol −1 ) [35], hydroquinone (11 kJ mol −1 ) [36], tetrachlorohydroquinone (19 kJ mol −1 ) [37], and carvacrol (16 kJ mol −1 ) [21].Although these energy barriers are too high to be surmounted in a low-temperature matrix environment, for all these previously studied derivatives of phenol compounds, the higher energy OH-rotamers escaped from spectroscopic detection in matrixes of inert gases because of their fast depopulation in favor of the corresponding lower energy OH-rotamers by means of quantum tunneling of hydrogen atom [38,39].The same behavior has been reported in various matrix-isolation studies carried out for other compounds containing conformers differing from each other by 180 • rotation of an OH group, namely carboxylic acids [40][41][42][43][44] and aminoacids [45].Therefore, it is expected that in the case of thymol, the gc and tc cis-OH-conformers shall also undergo fast tunneling OH rotamerization to the respective gt and tt trans-OH-forms in the low temperature Ar and Xe matrices, thus preventing experimental identification of the cis-OH-forms.
To theoretically address this issue, we carried out computations of reaction paths for the flip of the OH group in the (gc, gt) and (tc, tt) pairs of conformers.Initially, the firstorder transition states for OH-torsion were optimized, and the respective force constants were computed analytically.The intrinsic reaction paths were then followed from the transition states in both directions.The intrinsic reaction coordinate (IRC) was set using the "IRC = Cartesian" option of the Gaussian software and is expressed in units of Bohr.The computed reaction paths for the OH torsion in the (gc, gt) and (tc, tt) pairs are presented in Figure 3.
The probability of tunneling through a parabolic barrier was estimated by using the equation where a particle with mass m tunnels through a barrier with height E and width w, and h is the Planck constant [38,39].Using the calculated barrier heights E of 10.32 and 8.57 kJ mol −1 and widths w of 2.915 and 2.715 Bohr for the gc → gt and tc → tt conversions, respectively, one can estimate the tunneling probability (transmission coefficient) of the light H-atom (m = 1) in gc as 2.97 × 10 −8 , and in tc as 4.08 × 10 −7 .The tunneling rate is a product of the transmission coefficient and the frequency of attempts.Assuming that the OH group is vibrating at an OH torsional frequency of 317 cm −1 in gc and 275 cm −1 in tc (computed values), tunneling rates of 2.82 × 10 5 s −1 (half-life time of ~2.5 microseconds) and 3.36 × 10 6 s −1 (half-life time of ~0.2 microseconds) were estimated for the gc → gt and tc → tt rotamerizations, respectively.According to these results, there is no doubt that the gc and tc conformers are unlikely to survive long enough in low temperature Ar and Xe matrices to allow their spectroscopic detection.
Photochem 2022, 2, FOR PEER REVIEW 7 in the case of thymol, the gc and tc cis-OH-conformers shall also undergo fast tunneling OH rotamerization to the respective gt and tt trans-OH-forms in the low temperature Ar and Xe matrices, thus preventing experimental identification of the cis-OH-forms.To theoretically address this issue, we carried out computations of reaction paths for the flip of the OH group in the (gc, gt) and (tc, tt) pairs of conformers.Initially, the firstorder transition states for OH-torsion were optimized, and the respective force constants were computed analytically.The intrinsic reaction paths were then followed from the transition states in both directions.The intrinsic reaction coordinate (IRC) was set using the "IRC = Cartesian" option of the Gaussian software and is expressed in units of Bohr.The computed reaction paths for the OH torsion in the (gc, gt) and (tc, tt) pairs are presented in Figure 3.
The probability of tunneling through a parabolic barrier was estimated by using the equation where a particle with mass m tunnels through a barrier with height E and width w, and h is the Planck constant [38,39].Using the calculated barrier heights E of 10.32 and 8.57 kJ mol −1 and widths w of 2.915 and 2.715 Bohr for the gc → gt and tc → tt conversions, Conversion of the tt into the gt conformer involves the internal rotation of a bulky isopropyl group and should only occur through an over-the-barrier mechanism.The torsional barrier for this isomerization was computed to be ~24 kJ mol −1 (see Figures 1 and S1).This value is too high to allow the occurrence of a tt → gt relaxation in a matrix environment at a temperature as low as 15 K [46,47].Accordingly, both conformers are expected to be trapped in the as-deposited Ar and Xe matrixes with populations of 76-81% (gt) and 24-19% (tt).These values were obtained by considering the conversion of the estimated gas phase populations of gc to gt and of tc to tt (see Table 1).
It is interesting to compare the potential energy profile for the rotation of the isopropyl group in thymol with that computed for the closely related carvacrol and p-cymene molecules (see Figure S1).The heights of the barriers separating the conformers of these last two molecules differing from each other by internal rotation of the isopropyl group are ~10 kJ mol −1 , which is 2.5 times lower than the one computed for thymol.For these three molecules, the first-order transition state separating the conformers is characterized by having the tertiary C-H bond of the isopropyl group nearly perpendicular to the aromatic ring plane.In p-cymene and carvacrol, this gives rise to steric repulsions between the two methyl groups of the isopropyl fragment and their neighboring aromatic C-H bonds, whereas in thymol, one of the methyl groups stays too close to the OH group, leading to stronger steric repulsions, which is responsible for the significant increase in the barrier height associated with the isopropyl group rotamerization.As it will be shown below, this result has important implications for the outcomes of the experiments of thermal and IR excitations carried out for matrix-isolated thymol.

Infrared Spectra of Matrix-Isolated Thymol, Annealing and IR Irradiations
Monomers of thymol were trapped in Ar and Xe matrixes at 15 K, as described in Section 2.2, and the mid-IR spectra recorded for the freshly deposited samples are shown in Figure 4a,b (3675-3580 and 1650-700 cm −1 regions), while the spectral positions and relative intensities (expressed qualitatively) of the observed bands are collected in Table 2.The two experimental spectra are well reproduced by the population-weighted theoretical spectrum of the matrix-isolated compound, which is displayed in Figure 4c.This spectrum was obtained by summing the individual spectra calculated for the conformers present in the matrix, shown in Figure 4d, with the computed IR intensities weighted by their expected populations in the noble gas matrices: 0.79 for gt and 0.21 for tt (mean values obtained from the B3LYP and MP2 computations, see Table 1).From the comparison of the experimental and theoretical spectra, an assignment of the experimental bands was carried out and has been included in Table 2.     a Absorptions falling in the 3100-2800 cm −1 range, corresponding to the CH and CH 3 stretching vibrations, were not analyzed.Experimental intensities (Int.) are expressed qualitatively: vs = very strong; s = strong; m = medium; w = weak; vw = very weak.The most intense component of each multiplet band is highlighted in bold.b Calculated harmonic wavenumbers are scaled by 0.980 and 0.950, for the regions below 2000 cm −1 and above 3200 cm −1 , respectively.Some weak calculated absorptions with no correspondence in the experimental spectrum are not included (a full list of the vibrations calculated for the two conformers is provided in Table S3).c Based on the results provided by the "vibAnalysis" software [24] supported by ChemCraft animation of the vibrations.Atom numbering is presented in Scheme 1. Abbreviations: ν, stretching; δ, in-plane deformation; γ, out-of-plane deformation; τ, torsion; ρ, rocking; as, antisymmetric; s, symmetric; iso, isopropyl group, sh, shoulder.Signs "+" and "−" designate combinations of vibrations in the same phase and in the opposite phase, respectively.
The intense band observed in the experimental spectra at wavenumbers above 3600 cm −1 is assigned to the stretching vibration of the OH group (νOH).In the spectrum of thymol recorded in Ar matrix, this band appears as a relatively broad and split feature with two main components at 3643 and 3638 cm −1 and a shoulder at ~3635 cm −1 .OH stretching bands at practically the same positions have been identified in the IR spectra of phenol [48] (also having two main components) and of some phenol derivatives such as tyramine [49] and carvacrol [15], isolated in solid Ar.In the spectrum of thymol trapped in Xe, the νOH band is sharper and appears as a doublet at 3617/3611 cm −1 , shifted to lower frequencies, as compared with thymol in Ar.Such a shift of the OH-stretching frequency is typical for matrix-isolated compounds [32,33].In this range, the spectrum of thymol in a xenon matrix exhibits a very similar profile to that of the population-weighted theoretical spectrum.This similarity suggests that the higher and lower frequency band components in the experimental spectrum are assigned, respectively, to the gt and tt conformers, which agrees with the order of wavenumbers predicted for the νOH absorptions of the two rotamers: 3645.5 cm −1 (A th = 56.6 km mol −1 ) for gt and 3641.5 cm −1 (A th = 52.8km mol −1 ) for tt.
In the fingerprint region, a detailed comparison of the experimental and theoretical spectra allows the identification of seven bands that can be assigned to the most populated gt conformer.These bands appear between 1300 and 1000 cm −1   2), with the respective calculated absorptions predicted at 1278.9, 1218.4,1123.1 and 1047.4 cm −1 .
As we have mentioned in the previous section, the gc and tc conformers are not expected to be present in the as-deposited matrixes because of their fast depopulation (on the microsecond time scale) in favor of the gt and tt forms by quantum mechanical tunneling during deposition of the matrixes.This is experimentally supported by a thorough comparison between the experimental and calculated spectra shown in Figure S2, where intense absorptions calculated for the two higher-energy forms have no correspondence in the experimental spectra.
The experimental populations of the two matrix-populated conformers were estimated from the integrated intensities of selected pairs of band components (Ar: 1294/1278, 1219/1228, 1149/1123; Xe: 3617/3611, 1291/1274, 1216/1228, 1145/1120), normalized by the calculated absolute intensities of the corresponding vibration modes.The obtained values were 75 ± 6% for gt and 25 ± 6% for tt, which are in good concordance with the theoretical predictions.This gives an approximate gt:tt experimental ratio of 3:1, which will be discussed later on in this paper.
One attempt to promote conversion between the two rotamers in the matrix, in order to allow their unequivocal identification, consisted of performing annealing experiments.Taking into consideration the magnitude of the energy barrier corresponding to the conversion of the tt into the gt conformer (~24 kJ mol −1 ), this conformational relaxation may only be observed in a matrix with stronger relaxant properties and with the possibility of heating to the highest possible temperature.Xe matrix was therefore the obvious choice for performing these experiments.The sample was annealed up to 75 K, which is the upper limit of thermal stability of solid xenon.However, no changes in the band intensities were observed that could be attributed to the conversion of the higher energy conformer into the most stable form.Apparently, the high energy barrier separating the two matrix-isolated conformers, combined with necessity of internal torsion of a bulky isopropyl group, precludes the occurrence of a tt → gt relaxation in the Xe matrix, even upon heating the sample at a temperature as high as 75 K.According to the correlation proposed by Barnes [50], a temperature above 80 K would be needed to overcome an energy barrier of ~24 kJ mol −1 .Hence, the non-observation of a thermally induced tt → gt relaxation is not surprising under the present experimental conditions.
Another attempt to induce changes in the population of the two matrix populated conformers consisted of the excitation of monomers of thymol isolated in solid Xe with narrowband or broadband IR radiation.The narrowband IR excitations were undertaken by tuning the optical parametric oscillator at the first overtone of the OH stretching vibration (2νOH).Before performing the irradiations, a near-IR spectrum was recorded to identify the position of the 2νOH band.This band appears as a doublet with a main component centered at 7061 cm −1 , and a shoulder at 7047 cm −1 (see Figure S3), which is very close to its position observed in the near-IR spectrum of thymol dissolved in CCl 4 solution (7056 cm −1 ) [14].The shape of the 2νOH band is identical to that of the fundamental νOH feature observed in the mid-IR spectrum, meaning that the higher and lower frequency components should be ascribed, respectively, to the 2νOH absorptions of the gt and tt conformers.Absorption of a near-IR photon with wavenumber in the 7061−7047 cm -1 range introduces in the molecule an energy of ~84 kJ mol -1 .This amount of excitation energy is 3.5 times larger than the barrier height computed for the interconversion between the conformers in both directions, which, in principle, is more than sufficient to overcome the barrier for conformational isomerization for a free monomer in the gas phase by internal rotation of the isopropyl group.However, exposure of thymol isolated in a Xe matrix to laser light tuned to 7061 or 7047 cm −1 for ~30 min did not result in any spectral modifications that could indicate that a conformational transformation had occurred.We also attempted to promote changes in the population of the two matrix-isolated conformers by letting the matrix be exposed to broadband IR-radiation emitted from a kanthal wire electrically heated to a reddish-orange glow, a methodology that has been successfully applied by us to induce changes in the relative abundances of the two conformers of the structurally similar carvacrol isolated in an Ar matrix [15].However, in the present case, we did not observe any spectral indications of the occurrence of conformational changes.
A possible explanation for the lack of infrared-induced conformational isomerization of the isopropyl group in thymol tt or gt is that the energy deposited in the 2νOH mode is more efficiently transferred to the lower energy OH-rotamerization coordinate (Figure 5), which will result in the isomerization of tt into tc or gt into gc.However, because tc and gc are high-energy cis-OH conformers that should decay very fast by H-tunneling back to the most stable tt and gt trans-OH-forms, respectively, no spectral changes would be observed in the course of the infrared irradiation experiments monitored using steady-state spectroscopy.Interestingly, as shown in Figure 5, in comparison with thymol, the carvacrol molecule has the C 3 H 7 and OH competitive rotamerizations with an opposite barrier height trend.That feature is likely to provide for carvacrol a more efficient energy transfer to its lower energy C 3 H 7 rotamerization coordinate, even when the energy is deposited in the 2ν(OH) or ν(OH) mode.This explains why the conformational isomerization of the isopropyl group in carvacrol was previously observed upon broadband IR irradiation [15].

UV-Induced Transformations
Monomers of thymol 1 isolated in a low temperature Ar matrix were subjected to irradiation with broadband UV-light at λ > 200 nm emitted directly from a Xe/Hg lamp.The experimental conditions for the UV-irradiations were established taking into consideration the UV-spectrum of thymol in methanol solution, which exhibits absorption maxima at λ = 221 and 280 nm [13].The times of exposition of the sample to the UV-light varied from 2 to 120 min.After each irradiation, an IR spectrum was collected to monitor the transformations occurring in the matrix-isolated compound.
Spectral indication of occurrence of photoreactions started to be observed shortly after the first four minutes of irradiation.The most notorious spectral change was the appearance of a relatively broad feature within the 2125-2090 cm −1 range, centered at ~2106 cm −1 (see Figure 6a).This band is very characteristic of the antisymmetric stretching vibration of the ketene group [ν as ( C=C=O)] [21,48,[51][52][53][54], thus revealing the photogeneration of compounds bearing this fragment.According to the results of previous photochemical studies carried out for matrix-isolated compounds containing the phenol [15,21,48,49,54,55] moiety, the photoreaction pathway leading to the formation of ketenes is initiated by detachment of the H-atom from the O-H bond in 1.Then, recombination of the released Hatom with the highly reactive phenoxyl-type radical 2 at the orthoor pararing carbon atoms with respect to the CO group, gives rise to alkyl-substituted cyclohexadienones, which may exist in three isomeric forms, labeled as 3a, 3b and 3c (see Scheme 2).The preference of the H-atom to reconnect at these ring positions is theoretically supported by the results of an NBO analysis carried out at the B3LYP/6-311++G(d,p) level for 2. In fact, besides the oxygen atom, atoms C2, C4 and C6 are those exhibiting, by far, the largest natural spin densities (see Figure S4).The presence of 3 in the UV-irradiated thymol/Ar matrix is spectroscopically confirmed by the appearance of two bands at 1677 and 1659 cm −1 (Figure 6a), which are typical positions of the stretching vibration of a C=O group in cyclohexadienones [15,21,48,56].Finally, cleavage of the weaker C-C bond at α-position with respect to the C=O group in 3a or 3b, yields open-chain conjugated ketenes 5a or 5b, respectively (Scheme 2).

UV-Induced Transformations
Monomers of thymol 1 isolated in a low temperature Ar matrix were subjected to irradiation with broadband UV-light at λ > 200 nm emitted directly from a Xe/Hg lamp.The experimental conditions for the UV-irradiations were established taking into consideration the UV-spectrum of thymol in methanol solution, which exhibits absorption maxima at λ = 221 and 280 nm [13].The times of exposition of the sample to the UV-light varied from 2 to 120 min.After each irradiation, an IR spectrum was collected to monitor the transformations occurring in the matrix-isolated compound.
Spectral indication of occurrence of photoreactions started to be observed shortly after the first four minutes of irradiation.The most notorious spectral change was the appearance of a relatively broad feature within the 2125−2090 cm −1 range, centered at ~2106 cm −1 (see Figure 6a).This band is very characteristic of the antisymmetric stretching vibration of the ketene group [νas( C=C=O)] [21,48,51-54], thus revealing the To theoretically support the proposed assignments, we have carried out full geometry optimizations and harmonic vibrational calculations for the relevant isomeric forms adopted by 3 and 5, which are characterized in Tables S4-S6.For species 5, only the structures adopting Z-configuration around the central C=C bond have been considered in the calculations.This assumption was based on the results of a recent study on the photochemistry of 1,3-cyclohexadiene isolated in solid parahydrogen, where only the Z-type isomers of the open-ring photoproducts were found to appear during the initial stages of the UV-irradiations [57].The range of wavenumbers predicted for the ν as (C=C=O ) mode of 5 was 2106-2084 cm -1 (A th = 637-1485 km mol -1 ) and for the ν(C=O ) mode of 3 was 1679-1654 cm -1 (A th = 165-318 km mol -1 ).These values are very close to the positions of the new observed absorptions of photoproducts, emerging upon UV-irradiations of thymol (Figure 6).We also tried to identify a spectral signature of 2. However, the few strong infrared absorptions predicted for this species could not be observed experimentally, which is not surprising, due to its high reactivity.
the ortho-or para-ring carbon atoms with respect to the CO group, gives rise to alkyl-substituted cyclohexadienones, which may exist in three isomeric forms, labeled as 3a, 3b and 3c (see Scheme 2).The preference of the H-atom to reconnect at these ring positions is theoretically supported by the results of an NBO analysis carried out at the B3LYP/6-311++G(d,p) level for 2. In fact, besides the oxygen atom, atoms C2, C4 and C6 are those exhibiting, by far, the largest natural spin densities (see Figure S4).The presence of 3 in the UV-irradiated thymol/Ar matrix is spectroscopically confirmed by the appearance of two bands at 1677 and 1659 cm −1 (Figure 6a), which are typical positions of the stretching vibration of a C=O group in cyclohexadienones [15,21,48,56].Finally, cleavage of the weaker C−C bond at α-position with respect to the C=O group in 3a or 3b, yields open-chain conjugated ketenes 5a or 5b, respectively (Scheme 2).To theoretically support the proposed assignments, we have carried out full geometry optimizations and harmonic vibrational calculations for the relevant isomeric forms adopted by 3 and 5, which are characterized in Tables S4-S6.For species 5, only the structures adopting Z-configuration around the central C=C bond have been considered Alongside with the IR bands assigned to 3 and 5, a new band centered at ~1773 cm −1 was also observed in the infrared spectra recorded in the early stages of the UV-irradiations (Figure 6a).A band at nearly the same position has been previously detected in the photochemical studies performed for matrix-isolated phenol (1789 cm -1 ) [48] and carvacrol (1775 cm -1 ) [15] and has been ascribed to the stretching vibration of a C=O bond of a Dewar valence isomer.Therefore, an analogous species, labeled as 4 in Scheme 2, must be also formed during the UV-irradiations of 1.It is most likely produced with the intermediacy of 3a or 3b, by analogy with the identical process observed for the electronically equivalent α-pyrone [52, 58,59] and methyl-substituted α-pyrones [51,60].The wavenumbers predicted for the possible structures adopted by 4 (see Table S7) fall between 1801 and 1794 cm −1 (A th = 281-340 km mol −1 ), see Figure 6b.
Increasing the time of exposition of the sample to the UV-radiation led to the occurrence of photodecomposition reactions, as revealed by the appearance of a narrow band at ~2137 cm -1 which is a spectral indication of the formation of carbon monoxide (CO) (Figure 6a and Figure S5).This band was found to continuously grow up as the time of irradiation increases, whereas the opposite behavior was observed for the bands assigned to cyclohexadienones 3 (and the open-chain ketenes 5), see Figure S5.By analogy with our recent discussion on the photochemical behavior of matrix-isolated carvacrol [15], the most probable species resulting from the photodecomposition of thymol and its photoproducts are alkyl-substituted cyclopentadienes.These photoproducts are most likely formed by decarbonylation of the alkyl-substituted cyclohexadienones 3, which exist in photoequilibrium with the open-chain ketenes 5 (see Scheme 2) and as long as isomers 3 are consumed, isomers 5 are also consumed (indirectly) [15], see Figure S5.These species are constituted exclusively by carbon and hydrogen atoms, thus having intrinsically low intensities in the IR spectra.Furthermore, such decarbonylated products share common groups (methyl, isopropyl) with their precursors and their new bands may overlap with already existing bands of the precursors.Therefore, it is difficult to experimentally confirm their presence in the UV-irradiated matrix.

Concluding Discussion
In this work, we investigated the conformational structure of monomeric thymol by a combined experimental and computational approach.Geometry optimizations at the B3LYP/6-311++G(d,p) and MP2/aug-cc-pVDZ levels of theory predicted that thymol may adopt four unique conformations defined by the mutual orientation of the vicinal isopropyl and hydroxyl groups.Two of these conformations, with the gauche-isopropyl group, are doubly degenerate and represent mirror-equivalent structures with broad and shallow double-minima on the potential energy surface of the isopropyl torsion.The two conformations with the trans-isopropyl group represent narrow and deep single minima.The two cis-OH-conformers are destabilized due to the internal repulsions between the hydrogen atoms of the hydroxyl and isopropyl groups.The two trans-OH-forms (gt and tt) are the most stable and are expected to be the most populated forms in the gas phase of thymol and (consequently) also for thymol isolated in the cryogenic matrix.
Infrared spectroscopy was used to characterize the conformational composition of thymol experimentally.In agreement with the computations, the infrared spectrum of the as-prepared matrix revealed that only the two most stable thymol conformers (trans-OH-forms gt and tt, with the OH group pointing away from the vicinal isopropyl group) co-exist in the matrix.The experimental integrated intensities of the marker bands due to the gt and tt forms, reduced by their respective computed absolute infrared intensities, suggest an approximate 3:1 ratio in favor of gt.The two higher energy cis-OH forms, which may have residual populations in the gas phase prior to deposition, were shown to convert into their more stable trans-OH-counterparts by fast quantum mechanical tunneling during deposition of the matrix.
By analogy with the previously studied carvacrol (differing from thymol only by position of the OH group), we attempted to promote a partial conformational isomerization between the gt and tt conformers.For this purpose, we irradiated the matrix-isolated thymol either with broadband IR light (provided by electrically heated kanthal wire) or with narrowband near-IR light (provided by an optical parametric oscillator tuned to the frequency of the first overtone of the OH-stretching mode).Despite all of these attempts, the initial populations of the gt and tt conformers did not change.The annealing of thymol isolated in a xenon matrix up to the temperatures over 75 K (i.e., up to the limit of the thermal stability of solid xenon), aiming at relaxation of the higher-energy tt form into the global energy minimum gt form, did not result in the conformational interconversion either.From these experiments, one could safely postulate that conformational populations of the gt and tt forms in a freshly deposited matrix, in the 3:1 ratio, correspond to the ratio of conformational populations of these forms in the gas phase (at room temperature) during the deposition of the sample.Taking into consideration the computed relative Gibbs free energies of the gt and tt forms at room temperature, such a conformational composition is only possible when the gt and gc forms contribute to the Boltzmann equilibrium with the weighting factors of two and the tt and tc forms with the weighting factors of one.This represents an experimental proof of the fact that the two mirror-equivalent g + t and g − t forms of thymol, separated by a torsional energy barrier not exceeding 0.7 kJ mol −1 , are indeed two individual minima, and not a single symmetry-averaged ct minimum.
The lack of conformational relaxation of matrix-isolated thymol was explained with the aid of computations of potential energy surfaces for the torsion of the hydroxyl and isopropyl groups in thymol and in the related molecules.Furthermore, the torsional barrier for isopropyl group in thymol is nearly 2.5 times higher than in its isomeric carvacrol.On the other hand, the torsional barrier for the hydroxyl group in carvacrol and thymol are approximately equal, but in thymol it is about a half as low compared to the torsional barrier of its isopropyl group.In our mechanistic interpretation, the energy provided by vibrational excitation (infrared or near-IR) of thymol is dissipated much more efficiently via the hydroxyl torsion, and therefore the isopropyl torsion does not occur.
The matrix-isolated thymol was subsequently irradiated with broadband UV-light (λ > 200 nm), and the resulting transformations were monitored by IR spectroscopy.As is typical of phenol derivatives, the photoexcitation of thymol resulted in production of corresponding alkyl-substituted cyclohexadienones.Mechanistically, such a photochemical behavior was explained in terms of the OH group cleavage (with generation of an alkylsubstituted phenoxyl radical), followed by recombination of the released H-atom (also a radical) at the orthoor parapositions of the ring (with respect to the CO group).The natural spin densities of the phenoxyl-type radical, extracted from the NBO calculations, indicate that precisely the orthoor paracarbon atoms of the ring are the most favorable positions for the above-mentioned recombination of the radical pair.The cyclohexadienone photoproducts were found to undergo further isomerizations, yielding Dewar isomeric species and open-chain conjugated ketenes.Decarbonylation of the photoproducts was also observed for longer irradiation times.
On the whole, the present investigation fills a gap in the knowledge about the structural, vibrational, and photochemical properties of monomeric thymol.This fundamental knowledge constitutes a basis for the understanding of the bioactivity of this and related compounds.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/photochem2020028/s1. Tables S1 and S2, Cartesian coordinates for the conformers of thymol optimized at B3LYP/6-311++G(d,p) and MP2/aug-cc-pVDZ levels of theory; Table S3, Scaled wavenumbers (cm −1 ) and IR intensities (km mol −1 ) calculated for the thymol conformers by means of B3LYP harmonic vibrational calculations; Tables S4-S7, Geometries and selected energetic and vibrational parameters calculated at the B3LYP level for different photoproducts resulting from the UV-isomerizations of matrix-isolated thymol; Figure S1, B3LYP potential energy scans for the intramolecular torsion of the isopropyl group computed for thymol, carvacrol, and p-cymene; Figure S2, Detailed comparison of the experimental spectra of thymol isolated in Ar and Xe matrices with the spectra calculated for the four conformers of thymol; Figure S3, Fragment (7100-7020 cm −1 ) of the near-IR spectrum of thymol isolated in Xe showing the band assigned to the first OH-stretching overtone; Figure S4, Natural spin densities computed for the heavy atoms of the two conformers of the phenoxyl-type radical of thymol by means of a Natural Bond Orbital (NBO) analysis carried out at the UB3LYP/6-311++G(d,p) level of theory; Figure S5, Evolution of the bands due to photoproducts emerging in the 2160-2080 cm −1 and 1690-1650 cm −1 regions, as function of the time of UV-irradiation.

Scheme 1 .
Scheme 1.Molecular structures of thymol (a) and carvacrol (b), including numbering of the heavy atoms and identification of the two conformationally relevant dihedral angles for the thymol structure.The H atoms connected to the aromatic ring are implicit.

Figure 1 .
Figure 1.Contour plot of the potential energy surface computed for the thymol molecule as a function of the H−C7−C2−C1 (α) and H−O11−C1−C2 (β) dihedral angles (specified in Scheme 1).Both angles were incrementally fixed in steps of 20°, in all combinations, while all the remaining parameters were optimized at the B3LYP/6-311++G(d,p) level of theory.The minima are indicated in the plot as gt (g + t/g − t), gc (g + c/g − c), tt and tc (geometries are shown in Figure2).The color bar on the top represents the relative energy scale (zero is assigned to the energy of the global minimum, gt).The isoenergy lines are traced with steps of 1 kJ mol -1 .

Figure 1 .
Figure 1.Contour plot of the potential energy surface computed for the thymol molecule as a function of the H−C7−C2−C1 (α) and H−O11−C1−C2 (β) dihedral angles (specified in Scheme 1).Both angles were incrementally fixed in steps of 20• , in all combinations, while all the remaining parameters were optimized at the B3LYP/6-311++G(d,p) level of theory.The minima are indicated in the plot as gt (g + t/g − t), gc (g + c/g − c), tt and tc (geometries are shown in Figure2).The color bar on the top represents the relative energy scale (zero is assigned to the energy of the global minimum, gt).The isoenergy lines are traced with steps of 1 kJ mol -1 .

Figure 3 .
Figure 3. IRC profiles for the flip of OH-group in (gc, gt) (a) and (tc, tt) (b) pairs of thymol conformers, computed at the B3LYP/6-311++G(d,p) level of theory in Cartesian (non-massweighted) coordinates.The vertical arrows (E = 10.32 or 8.57 kJ mol −1 ) designate the computed ZPEcorrected energies of the transition states with respect to the reactants (gc or tc, left in each frame).The horizontal arrows designate the estimated barrier widths (w = 2.915 or 2.715 Bohr) of the barrier at the ZPE level of the reactants.The relative zero energies correspond to the products (gt or tt, right in each frame).

Figure 3 .
Figure 3. IRC profiles for the flip of OH-group in (gc, gt) (a) and (tc, tt) (b) pairs of thymol conformers, computed at the B3LYP/6-311++G(d,p) level of theory in Cartesian (non-mass-weighted) coordinates.The vertical arrows (E = 10.32 or 8.57 kJ mol −1 ) designate the computed ZPE-corrected energies of the transition states with respect to the reactants (gc or tc, left in each frame).The horizontal arrows designate the estimated barrier widths (w = 2.915 or 2.715 Bohr) of the barrier at the ZPE level of the reactants.The relative zero energies correspond to the products (gt or tt, right in each frame).

Figure 4 .
Figure 4. Experimental IR spectra of thymol isolated in (a) Ar and (b) Xe matrices at 15 K; (c) Theoretical IR spectrum of thymol, obtained as a sum of the B3LYP/6-311++G(d,p) spectra calculated for the two experimentally relevant conformers with their intensities scaled by the respective expected populations in the noble gas matrices.This spectrum was simulated by applying Lorentzian functions (fwhm = 2 cm −1 ) centered at the scaled wavenumbers; (d) Scaled wavenumbers and unscaled IR intensities extracted from the harmonic vibrational calculations carried out at the B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Figure 4 .
Figure 4. Experimental IR spectra of thymol isolated in (a) Ar and (b) Xe matrices at 15 K; (c) Theoretical IR spectrum of thymol, obtained as a sum of the B3LYP/6-311++G(d,p) spectra calculated for the two experimentally relevant conformers with their intensities scaled by the respective expected populations in the noble gas matrices.This spectrum was simulated by applying Lorentzian functions (fwhm = 2 cm −1 ) centered at the scaled wavenumbers; (d) Scaled wavenumbers and unscaled IR intensities extracted from the harmonic vibrational calculations carried out at the B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Figure 5 .
Figure 5. B3LYP/6-311++G(d,p) relaxed potential energy scans calculated for (a) carvacrol and (b) thymol regarding the rotamerization of the hydroxy (OH) and isopropyl (C 3 H 7 ) groups.The bottom and top x-axis (red and blue) refer to the corresponding dihedral angle changes of the first and second most stable conformations, respectively, which are illustrated for each species at the right side of the panel.The conformers of carvacrol are designated according to the scheme adopted in this work for thymol: the first letter designates the orientation of the tertiary C-H bond of the isopropyl fragment relative to the C-O bond, while the second refers to the orientation of the O-H group relative to the vicinal alkyl group.

Figure 6 .
Figure 6.(a) Fragment of the IR spectrum of thymol 1 (Ar, 15 K) recorded after 4 min of UVirradiation (black trace) and before any UV-irradiation (red trace).(b) B3LYP/6-311++G(d,p) calculated wavenumbers (scaled) and IR intensities corresponding to the antisymmetric C=C=O stretching modes of all surveyed isomeric forms of the open-ring ketenes 5 (green), C=O stretching modes of all Dewar valence isomers 4 (blue), and C=O stretching modes of all isomeric forms of the alkyl-substituted cyclohexadienones 3 (orange).For the structures of 3, 4, 5, see Scheme 2.

Table 1 .
Relative energies and Boltzmann populations calculated for the thymol conformers at different levels of theory.a Conformer Figure 2. Geometries of the thymol conformers, including some geometrical parameters calculated at the MP2/aug-cc-pVDZ.For the doubly generated gt and gc conformers, only one of the mirror-image structures (g

+ t and g + c) is shown.Table 1 .
Relative energies and Boltzmann populations calculated for the thymol conformers at different levels of theory a .

/gc or tt/tc, respectively
).Values in parentheses correspond to the expected populations of conformers gt and tt with the contributions of the populations of conformers gc and tc, respectively.b Single-point calculations performed on the MP2/aug-cc-pVDZ optimized geometries.

Table 2 .
Bands observed in the experimental IR spectra of thymol isolated in Ar and Xe matrices, compared with the harmonic wavenumbers (/cm −1 ) and absolute IR intensities (A th /km mol −1 ) calculated for the gt and tt conformers at the B3LYP/6-311++G(d,p) level, and vibrational assignment.
a Xe, 15 K a Calc.gt b Calc.tt b Assignment c ῦ Int.

Table 2 .
Bands observed in the experimental IR spectra of thymol isolated in Ar and Xe matrices, compared with the harmonic wavenumbers ( and unscaled IR intensities extracted from the harmonic vibratio B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Table 2 .
Bands observed in the experimental IR spectra of thymo compared with the harmonic wavenumbers (/cm −1 ) and absol calculated for the gt and tt conformers at the B3LYP/6-311 assignment./cm−1 ) and absolute IR intensities (A th /km mol −1 ) calculated for the gt and tt conformers at the B3LYP/6-311++G(d,p) level, and vibrational assignment.Figure 4. Experimental IR spectra of thymol isolated in (a) Ar and (b) Xe matrices at 15 K; (c) Theoretical IR spectrum of thymol, obtained as a sum of the B3LYP/6-311++G(d,p) spectra calculated for the two experimentally relevant conformers with their intensities scaled by the respective expected populations in the noble gas matrices.This spectrum was simulated by applying Lorentzian functions (fwhm = 2 cm −1) centered at the scaled wavenumbers; (d) Scaled wavenumbers and unscaled IR intensities extracted from the harmonic vibrational calculations carried out at the B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Table 2 .
Bands observed in the experimental IR spectra of thymol isolated in Ar and Xe matrices, compared with the harmonic wavenumbers (/cm −1) and absolute IR intensities (A th /km mol −1 ) calculated for the gt and tt conformers at the B3LYP/6-311++G(d,p) level, and vibrational assignment.Experimental IR spectra of thymol isolated in (a) Ar and (b) Xe matrices at 15 K; (c) Theoretical IR spectrum of thymol, obtained as a sum of the B3LYP/6-311++G(d,p) spectra calculated for the two experimentally relevant conformers with their intensities scaled by the respective expected populations in the noble gas matrices.This spectrum was simulated by applying Lorentzian functions (fwhm = 2 cm −1 ) centered at the scaled wavenumbers; (d) Scaled wavenumbers and unscaled IR intensities extracted from the harmonic vibrational calculations carried out at the B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Table 2 .
Bands observed in the experimental IR spectra of thymol isolated in Ar and Xe matrices, compared with the harmonic wavenumbers (/cm −1) and absolute IR intensities (A th /km mol −1 ) calculated for the gt and tt conformers at the B3LYP/6-311++G(d,p) level, and vibrational assignment.Experimental IR spectra of thymol isolated in (a) Ar and (b) Xe matrices Theoretical IR spectrum of thymol, obtained as a sum of the B3LYP/6-311++G(d,p) spec for the two experimentally relevant conformers with their intensities scaled by th expected populations in the noble gas matrices.This spectrum was simulated Lorentzian functions (fwhm = 2 cm −1 ) centered at the scaled wavenumbers; (d) Scaled w and unscaled IR intensities extracted from the harmonic vibrational calculations carri B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Table 2 .
Bands observed in the experimental IR spectra of thymol isolated in Ar and compared with the harmonic wavenumbers (/cm −1) and absolute IR intensities (A calculated for the gt and tt conformers at the B3LYP/6-311++G(d,p) level, and assignment.Experimental IR spectra of thymol isolated in (a Theoretical IR spectrum of thymol, obtained as a sum of the B for the two experimentally relevant conformers with thei expected populations in the noble gas matrices.This sp Lorentzian functions (fwhm = 2 cm −1 ) centered at the scaled w and unscaled IR intensities extracted from the harmonic vib B3LYP/6-311++G(d,p) level for the gt and tt conformers.

Table 2 .
Bands observed in the experimental IR spectra of t compared with the harmonic wavenumbers (/cm −1) and calculated for the gt and tt conformers at the B3LYP/ assignment.