Journal of the Chilean Chemical Society
versión On-line ISSN 0717-9707
J. Chil. Chem. Soc. v.53 n.1 Concepción mar. 2008
J. Chil. Chem. Soc, 53, N° 1 (2008)
"CONFORMATIONAL PREFERENCE IN 4,6-DIMETHYL-l,3-THIOXANE"
J. GUILLERMO CONTRERAS* AND LORENA A. GERLI
Facultad de Ciencias Químicas, Universidad de Concepción, Casilla 160-C, Concepción, CHILE, e-mail: email@example.com
Trans-4,6-dimethyl-l,3-thioxane shows two conformations as a result of the heterocyclic ring inversion. The I and II conformers interconversion process has been studied by means of ab initio methods. The thermodynamics of this reaction calculated at the MP2/6-311+G**// HF/6-31G" and B3LYP/6-311+G"// B3LYP/6-31G" levels yield AG° values of ca. 1.5 kcal/mol in the gas phase, implying that II is the most important species in the gas phase and that its relative concentration is ca. 93%. In solution of low and medium polarity solvents, PCM method predicts an 85% concentration for this species, consistent with an 83% estimated from 'H-NMR spectra. The experimental NMR properties ( 8 and 'JC_H ) have not been reported since I and II separation could not be achieved. The results obtained for two long range "J2ax6e and 6J2e Je coupling constants of ca. 0.5 Hz are in good agreement with 0.41 Hz, calculated in the present work. The calculated chemical shifts for the C2 and C5 protons show that 8Heq < 8Hax implying that some anomalous effects are operating, though the 8 seems to be moderately insensitive C-H bond distances changes. The spin-spin coupling constants are much more sensitive to changes occurring in the ring. Thus, normal and reverse Perlin effects for 'J are predicted for C2 -H and C5-H are operating, respectively. The cC5 - Heq > (7*C - S interaction would be the responsible for the very small elongation of the C5 - Heq over the C5 - Hax bond distances.
Keywords: 4,6 - dimethyl - 1,3 - Thioxane; conformational preference, NMR properties, NBO analysis.
The study of the factors ruling the preferred conformation in sulfur-containing six-membered heterocycles has been an active area of research for a time already. Normally, cyclohexane has been taken as the model undistorted chair like structure to be compared with other substituted species such as 1,3-dioxane ,1,3-dithiane and 1,3-thioxane. Structure deviations and reactivity can be predicted for dialkyl derivatives of these heterocycles. In fact, the introduction of two heteroatoms produce important distortions in the "ideal" chair structure of cyclohexane. Thus, changes occurring in the aliphatic region of thioxane, show that this part of the molecule is less planar than dioxane, thought the former is more easily distorted upon hydrogen substitution by alkyl groups than the later. Conformational analyses on substituted derivatives of these heterocycles have provided additional information on the different factors governing the preferred conformation adopted by the ring and also allowed to understanding some anomalous structural and spectroscopic features1"5. In fact, particular interest have been paid to chemical shift and one-bond 13C-'H spin-spin coupling constants as they allow directly to study hyperconjugation. In the present work, we have studied the two conformers (I and II) equilibrium in trans-4,6-dimethyl-l,3-thioxane (DMT) and see whether the calculated model chemistry is consistent with the coexistence of both species as suggested by the 'H-NMR spectra6"7. We have also explored the effect of the solvent, by using PCM model8"11, on such equilibrium and approximately estimated the conformers relative concentrations. The transition state for the two conformers interconversion reaction have been found and the effect of the solvent on the activation energy also studied.
Conformers I and II initial structures were fully optimized at HF/6-31G", MP2/6-31G" and B3LYP/6-31G" levels using Gaussian 98 and 03 series of programs12. Frequency and IR intensities were predicted at the equilibrium geometries (B3LYP/6-31G") yielding all real frequencies so that all calculated structures are local minima. The calculated energies were corrected for zero point vibrational energies. Thermodynamic data ( 298 K , 1 atm) were derived from the vibrational modes within the harmonic oscillator, rigid rotor and ideal gas approximations. Energy calculations were carried out using a flexible basis set like 6-311+G" that includes both polarization and diffuse functions. Diffuse functions are important for molecules having lone pair of electrons (LP)13, and to into account the relative diffuseness of LP. On the other hand, polarization functions are useful in conformational studies where stereoelectronic effects play an important role14. Enthalpies and free energy changes of interconversion were obtained from the corresponding sums of electronic and thermal enthalpies and free energies, respectively. The solute-solvent interactions were calculated with the polarized continuum method (PCM )-SCRF method for which its capacity to reproduce solvation energies15 has been proved. In general PCM behaves well, failing just in few cases. In applying PCM the gas phase geometries were taken, since optimization in CS2 produce almost no change in the gas phase structural parameters of I and II (largest bond length deviation ca. 0.009 Á, whereas bond angles vary ca. 0.3 degrees) . Free energies in solution (G°SOLN )16Q were calculated from G°SOLN = AG°QAS + AAG°S where the later term is the PCM solvation energy that includes both electrostatic and non-electrostatic contributions. The concentration of each species was estimated from the corresponding AG°SOLN value.
In order to study the electronic structures of conformers I and II, Natural Bond Orbital (NBO)16b analyses were carried out on the B3LYP/6-31G" geometries and the energies of hyperconjugation interactions, the Fock (F.) matrix element corresponding to the orbital interactions and the energies of the donor and acceptor orbitals, derived. The hyperconjugation interaction energies were obtained from the second-order perturbation approach: E = -2<a I F | a*^2 /£ ,- £ , where the integral corresponds to the F.. matrix element between interacting orbitals i and j whose energies are £ and £ ,. 2 stands for the population of the donor orbital. There is a very good linear correlation between the second-order perturbation hyperconjugative energies and the deletion procedure23.
Theoretical :H shifts (8) and the one-bond 13C-'H coupling constants ('JC.H), were obtained at the B3LYP/6-31G" geometries. The TMS ,used as standard, geometry was also calculated at this level. Magnetic shielding tensors were calculated at B3LYP/6-311+G"// B3LYP/6-31G" level of theory using the GIAO (Gauge-Independent Atomic Orbital) method. Chemical shifts were derived with respect to the NMR isotropic magnetic shielding tensors (in ppm) from the corresponding standard tensor (TMS = 32.1 ppm for 1H)16. The transition state for the I <=> II interconversion reaction was located using the synchronous transit-guided quasi-Newton (STQN) method1718 characterized by a frequency calculation which yielded a single imaginary frequency. The animation of such a vibrational mode using Gauss View19 program, confirm that the imaginary frequency corresponds to a ring deformation mode that conducts either to the reactant ( conformer I) or the product (conformer II). The TS > I and TS > II paths were followed by an intrinsic reaction coordinate (IRC) procedure20. IRC calculations lead to the conclusion that the TS connects both reactant and product. These three independent methods to analyze this vibrational mode point towards a true TS
RESULTS AND DISCUSSION
A) Molecular Structures
For conformers I and II, selected bond distances and angles describing mainly the ring geometries are given in table 1, whereas figure 1 shows the atom numbering used in table 1 and in setting the Z-matrix. The chair structure of conformers I and II has been divided into two parts: the one formed by Sl-C2-03 set of atoms (set 1) and the hydrocarbon part formed by C4-C5 and C6 atoms (set 2). Thus in I, set 1 is down and set 2 is up; the opposite situation is found in II (see figure 1). Thioxane is not planar in the zone connecting sets 1 and 2. In fact, in dioxane, dithiane and thioxane, the dihedral S1C6C403 varies from 0.00 , 0.08 and 1.39° for dioxane, dithiane and thioxane, respectively. The larger value found in thioxane can be attributed to the C-S and C-O bond distances difference and to the dissimilar spatial requirements of the heteroat oms (S and O).
The non-regularity of thioxane is also reflected in the values of the torsion angles t and t with values of 50.7 and 64.8°, respectively. In 4,6-dimethyl derivatives, t (C2S1C6C5) are ca. 50.9 and 60.3° in I and II, respectively, whereas for t (C203C4C5) values of 47.4° and 65.4° were found for these conformers. Comparing these values with those of the parent compound it can be inferred that I is more distorted in the oxygen side, whereas II is distorted in S side. Conformer II would stabilize due to hyperconjugative interactions of the type LPS^ c* C5-C6, LPS^ c*C6-Me and LPS^c*C2-03 (see NBO analysis ).
The transition state for the I <=> II interconversion reaction looks like a large distorted species. In fact, as SI moves up by ca 12° , 03 remains in its position. The most drastic changes are experienced by C2 and C5 that move up and down, respectively. These combined movements leads to an almost planar ring geometry (C203C4C5C6 ~ Io) with the exception of S1. The TS methyl groups on C4 and C6 are in axial and equatorial, respectively.
B) Energetics and Thermodynamics
Energetic and thermodynamics of conformers I <=> II interconversion reaction are given in table 2. The MP2 and DFT methods include electron correlation that is important in conformational analysis13,21. In fact, at MP2/6-31G" and B3LYP/6-31G" levels the whole molecular system is stabilized by ca. 730 and 1701 kcal/mol (relative to de HF level), respectively. Basis sets with diffuse functions, such as 6-311+G" and 6-311++G", are also important to account the diffuse nature of lone pairs of electrons13, whereas the inclusion of polarization functions are useful when hyperconjugative stereoelectronic interactions and third row elements are present22"25. The thermodynamics of the I <=> II interconversion reaction were calculated at the MP2/6-311+G" // HF/6-31G" and B3LYP/6-311+G*7/B3LYP/6-31G" levels. Both levels of theory indicate that the use of a larger basis set than 6-31G" produces the system to gain ca. 99 (MP2) and 64 (DFT) kcal/mol. The addition of diffuse functions to the hydrogens leads to no further stabilization of the molecular system. An approximate relative concentration % of II ca. 93% points towards a very low concentration of conformer I in the gas phase. In fact, a AG° of ca.-1.566 kcal/mol indicates that the interconversion reaction is extensively displaced to the conformer II formation.
This conclusion is in good agreement with the experimental finding that estimates ca 85% of II concentration26. As usual in conformers interconversion reactions27,28, entropy remains constant (AS = 0.121 cal K"1 mol"1), whereas enthalpy change appears controlling the interconversion. On account of the almost complete planarity of the ring, TS does not resemble the reactant or the product. The animation of the imaginary normal mode shows that C5 and 03 move up and down with rather large amplitudes, The C2H2 group simultaneously performs the same movement in such a way that conformer I is clearly "recovered". Intrinsic reaction coordinate calculations (IRC) both in the forward and reverse directions ( both at HF and DFT/6-31G") lead to the reactant and product and thereby the TS is fully characterized. Figure 2 shows the energy profile which is consistent with a single saddle point at the TS. No other saddle points were "detected" most likely due to steric effects of sulfur atoms that would produce a larger repulsion with the bulky methyl group leading to rather rigid structures.
C) Solvent Effect
The solvent-solute interactions were calculated using PCM method. Solvents of low and medium high polarity (dielectric constants: 2.60,4.90) and 46.70) were employed. The free energy changes in solution were obtained, as usual, from AG°SOLI are given in table 3.
From table 3 it can be inferred that the solvent polarities do not influence the interconversion reaction since conformer II relative concentration remains constant throughout the solvent dielectric constant. It likely that the small effect one can derived by applying PCM method is due to the very small dipole moment difference (Au = 0.111 D). Conformer II is the main species in the gas phase ( ca. 93%) due to its intrinsic stability, but decreases to ca. 86% in solution in agreement with the experimental finding that in CS2 solution, this species is in ca. 85% concentration26.
D) NMR Spectra
The NMR properties were calculated al the B3LYP/6-3 lG(d,p) geometries and the GIAO magnetic shielding tensor obtained at the B3LYP/6-311+G(d,p) level of theory. To control the reliability of this level, the 1,3-thioxane chemical shifts and one-bond
13C-'H coupling constants were also derived and compared with the theoretical and experimental data3,4 (table 4). The good agreement between the calculated and experimental NMR data allows to predict confidently the spectral properties for DMT conformers I and II. The 'H-NMR spectra of a mixture of both conformers would show a typical two AB system pattern for the axial and equatorial protons attached to C2, i.e., two well separated four signal multiplets.
The two A parts of the I <=> II system ( Ha^II) and H|5ix(I)) would show up at lower field than the two B parts of the systems ( Ha (I) and H|5 (II) ) will appear at higher field in a similar fashion the 4-ethyl-6-methyl-l,3-dithiane behaves716. As usual, the AB system with the methyl group equatorial (conformer I) and next to the sulfur will presents the highest intensity. At room temperature, both multiplets would collapse into two doublets calculated at 4.82 and 4.61 ppm (A8=0.21 ppm) and 5.48 and 5.39 ppm (A8=0.09 ppm), respectively. In the parent 1,3-thioxane (see table 4) and conformers of DMT (see table 5) the calculations indicate that 'Jc.He > 1Jc_Stx implying a normal Perlin effect is working at C2, whereas at C5 the oposite situation is observed and thereby the reverse Perlin effect must be operating.
The calculated spin-spin coupling constants for conformers I and II qualitatively reflect the effect that larger C-H bond distances lead to smaller 'JCH. This observation agree well with the data given in tables 4 and 5. The elongation of C5-Heq can be explained in terms of a dominant cC5-Heq > (7*C-S two-electron/two-orbitals interctions. In fact, the second order perturbation energies for the above interactions are ca. 5.56 and 5.85 kcal/mol for I and II ,respectively. The reverse interactions (cC - S > (7*C5 - Heq are less favored with E ca. 3.20 and 1.32 kcal/mol) The sensitivity of coupling constants to the changes occurring in the heterocyclic ring upon substitution is a well stablished observation3"5M now. No relationship between the chemical shifts (8) and bond lengths have been found, most likely due to the rather insensitivity of 8 to geometric changes in the ring. However, it is worth noting that equatorial and axial protons are displaced to higher and lower field, respectively, due to rjC-S > (7*C-Heq hiperconjugative interactions.
The calculations carried out on the conformers I and II of 4,6-dimethyl-1,3-thioxane, allow to draw some interesting conclusions:
l.-The presence of two heteroatoms in the ring leads to important distortions that are reflected in the values of the torsional angles. The main deviations, with respect to dithiane, occur on the angles C2S1C4C5('U1) and C203C6C5('U2). T deviates by ca. 7° and 3°, whereas t departs by ca. 10° and 8°, in I and II respectively.
2.-Conformer II would be more stable than I due to hiperconjugative interactions of the type LPS^ a* C5-C6, LPS^ c*C6-Me and LPS^C7*C2-03. The sulfur electron lone pairs are delocalising onto antiperiplanar C5 - C6 and C2 - 03 empty antibonding orbitals. These delocalisation are supported by high bond populations and orders
3.-The gas phase calculations point to conformer II as the most stable species with a relative concentration of ca. 92%, whereas in solution conformer I increases its concentration to ca. 12%
4.-The calculated chemical shifts for the C2 and C5 protons show that SHeq < SHax implying that some anomalous effects are operating, though the 8 seems to be somewhat insensitive C-H bond distances.
5.-The spin-spin coupling constants are much more sensitive to changes occurring in the ring. Thus normal Perlin effect for ' J are predicted for C2 -H, whereas at C5 the reverse Perking effect is operating. The cC5 - Heq > rj*C - S interaction would be the responsible for the very small elongation of the C5 - Heq over the C5 - Hax bond distances.
1.- E. Juaristi, Ace. Chem. Rev. 22, 357 (1989) [ Links ]
2.- E. Juarist, G. Cuevas, A.Vela, J. Am. Chem. Soc. 116, 5796 (1994) [ Links ]
3.- G. Cuevas, E. Juaristi, A.Vela, J. Mol.Struct.(THEOCHEM) 418, 231 (1997) [ Links ]
4.- G. Cuevas, E. Juaristi, A. Vela, J. Phys. Chem. 103, 932 (1999) [ Links ]
5.- I.V. Alabugin , J. Org. Chem. 65, 3910 (2000) [ Links ]
6.- J. Gelan, M. Anteunis, Bull. Soc, Chim, Beiges 77, 423 (1968) [ Links ]
7.- J. Gelan , M. Anteunis, Bull. Soc. Chim. Beiges. 78, 599 (1969) [ Links ]
8.- J. Tomasi, M. Pérsico, Chem. Rev. 94, 2027 (1994) [ Links ]
9.- M.T. Canees, V. Mennucci, J. Tomasi, J. Chem. Phys. 107, 3032 (1997) [ Links ]
10.- V Barone, M. Cossi, J. Tomasi, J.Chem. Phys, 107, 3210(1997) [ Links ]
11.- V. Barone, M. Cossi, J. Tomasi, J. Comput. Chem. 19, 404 (1998) [ Links ]
12.- M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. Morokuma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, 0. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, 1. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, C. Gonzalez, and J. A. Pople, Gaussian 03, Revision C.02. Gaussian, Inc., Wallingford CT, 2004. [ Links ]
13.- J.B. Foresman, A. Frisch,./'Exploring Chemistry with Electronic Structure Methods" 2ntl edition , Gaussian, Inc; Pittsburgh, PA, 1996. [ Links ]
14.- F. Freeman, K.T. Le, J. Phys. Chem. 107, 2908 (2003) [ Links ]
15.- J.G. Contreras, S.T. Madariaga, J. Phys. Org. Chem. 16, 47 (2003) [ Links ]
16.- a) J.G. Contreras, L.A. Gerli, J. Chil. Chem. Soc. 52, 1271 (2007) [ Links ]
b).- Reed AE, Curtus LA, Weinhold F, Chem. Rev. 88, 899 (1988) [ Links ]
17.- C. Pen, P.Y. Ayala, H.B. Schlegel, M.J. Frisch , J. Comput. Chem., 17, 49(1996) [ Links ]
18.- C. Pen, H.B. Schlegel, Isr. J. Chem. 33, 449 (1994) [ Links ]
19.- GaussView, Revision 3.09, Gaussian Inc., Wallingford CT, 2004 [ Links ]
20.- C. Gonzalez, H.B. Schlegel, J. Phys. Chem. 94, 5523 (1990) [ Links ]
21.- D.M. Chapman, R.E. Hester. J. Phys. Chem. A 101, 3382(1997) [ Links ]
22.- G. Cuevas, E. Juaristi,. A. Vela, J. Phys. Chem. A 103, 932 (1999)
23.- F. Freeman, M.L Kasner, W.J. Hehre, J. Phys. Chem. A 105, 10123 (2001) [ Links ]
24.- F. Freeman, K.U. Do, J. Mol. Struct. (THEOCHEM), 577, 43 (2002) [ Links ]
25.- F. Freeman, J.A. Kasner, M.L. Kasner, W.J. Hehre J. Mol. Struct.(THEOCHEM) 496, 19 (2000) [ Links ]
26.- J. Gelan, M. Anteunis,. Bull. Soc. Chim. Beiges, 77, 447 (1968) [ Links ]
27.- J.G. Contreras, S.T. Madariaga, J. Mol. Struct. (THEOCHEM) 572, 235 (2001) [ Links ]
28.- J.G. Contreras, S.M. Hurtado, L.A. Gerli. S.T Madariaga J. Mol. Struct. (THEOCHEM) 713, 207 (2005) [ Links ]
(Received: 7August 2007 - Accepted: 11 January 2008)