Chem. Mater. 1998, 10, 1459-1467
1459
Structural and Conformational Studies of Quaterthiophenes Having Extra Bulky Terminal Groups Hitoshi Muguruma,*,†,| Kazuya Kobiro,‡,| and Shu Hotta*,§ Research Center for Advanced Science and Technology, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153, Japan; Niihama National College of Technology, 7-1 Yagumocho, Niihama, Ehime 792, Japan; and Advanced Materials Research Laboratory, Matsushita Research Institute Tokyo, Inc., 3-10-1 Higashimita, Tama-ku, Kawasaki 214, Japan Received January 23, 1998
An integrated set of experiments including crystallographic analysis and spectroscopy have been carried out to investigate structural and conformational features of tetramer oligothiophene compounds (quaterthiophenes). These quaterthiophenes are characterized by the presence of extra bulky terminal groups N,N-bis(trimethylsilyl)aminomethyl or (2,2,5,5tetramethyl-1-aza-2,5-disila-1-cyclopentyl)alkyl, where the alkyl is methyl or ethyl. These bulky groups produce variation in molecular conformation and morphology. The X-ray analysis shows that two types of conformers are presently stable in the crystals: one is all-anti and the other is syn-anti-syn with respect to the molecular long axis. The molecular conformations have been further studied in solutions by spectroscopic methods, especially temperature-dependent ones. Of these, the UV-vis spectroscopy gives important information about the relative stability between the conformers. A clearly resolved isosbestic point means that two interconvertible conformers are present as the major species in the solutions. Numerical analysis of the spectra indicates that the free energy difference between these two is ∼0.8-1.2 kcal/mol and that an absorption edge of the less stable conformer is significantly blue-shifted relative to the other. The former conformer is thought to be generated on account of the destruction of π-conjugation caused by distortion around the bonds connecting the thiophene rings. This associates the conformers with rotamers. The interconversion rate processes between these rotamers have been studied by the dynamic 1H NMR spectroscopy. The line-width analysis of the coalesced spectral lines that are observed for the coupled protons implies that the free energy of activation for the interconversion is around 20 kcal/mol. These spectroscopic data are compared with the results of the quantum chemical calculations.
Introduction Crystal structures and conformations of conjugated polymers and oligomers are currently being studied extensively. Among these materials, polythiophenes and oligothiophenes continue drawing particularly great attention. A good deal of crystallographic and spectroscopic information has been accumulated to interpret their structural and conformational features. In a rigid lattice of the polythiophene, polymer chains assume the fully stretched S-anti conformation1 and are packed in a well-known herringbone structure2 where the S-anti planar macromolecules somewhat turn aside from the face-to-face arrangement; this structure is quite often observed for many organic crystals. Mean†
University of Tokyo. Niihama National College of Technology. Matsushita Research Institute Tokyo, Inc. | Present address: Department of Environmental Systems Engineering, Kochi University of Technology, Tosayamada, Kochi 782-8502, Japan. Fax: +81-8875-7-2520. E-mail:
[email protected]. (1) Bru¨ckner, S.; Porzio, W. Makromol. Chem. 1988, 189, 961. (2) (a) Gavezzotti, A.; Filippini, G. Synth. Met. 1991, 40, 257. (b) Bernstein, J.; Sarma, J. A. R. P.; Gavezzotti, A. Chem. Phys. Lett. 1990, 174, 361. (c) Desiraju, G. R.; Gavezzotti, A. Acta Crystallogr. 1989, B45, 473. (d) Porzio, W.; Destri, S.; Mascherpa, M.; Bru¨ckner, S. Acta Polym. 1993, 44, 266. ‡ §
while, the polythiophenes appropriately substituted with various chemical groups such as alkyl3,4 show interesting structural characteristics on account of these substituents. In the alkyl-substituted polymers, for instance, the polythiophene backbone planes that have the S-anti conformation as well are forced to take the rigorous face-to-face polymer chain packing. Other interesting aspects of these polymers can be observed in solutions as chromisms.5 These chromisms result from the fact that a variety of conformations allowed in the solutions are coupled to the electronic states of the π-conjugated system. The conformational change accompanied by a color change is believed to be associated with the rotational defects (rotamers) introduced upon dissolution of the polymers, these defects breaking up the π-conjugation and causing segmentation or localization of the electronic wave functions.6 To study the variation in the crystal structures and molecular conformations more precisely, the oligoth(3) Tashiro, K.; Ono, K.; Minagawa, Y.; Kobayashi, M.; Kawai, T.; Yoshino, K. J. Polym. Sci., Polym. Phys. Ed. 1991, 29, 1223. (4) Winokur, M. J.; Wamsley, P.; Moulton, J.; Smith, P.; Heeger, A. J. Macromolecules 1991, 24, 3812. (5) Rughooputh, S. D. D. V.; Hotta, S.; Heeger, A. J.; Wudl, F. J. Polym. Sci., Polym. Phys. Ed. 1987, 25, 1071.
S0897-4756(98)00041-6 CCC: $15.00 © 1998 American Chemical Society Published on Web 04/21/1998
1460 Chem. Mater., Vol. 10, No. 5, 1998
iophenes are particularly suited since their chemical structure can be tailored as desired and the crystals of high quality for the structure analysis are readily accessible. In particular, the oligothiophenes having long or bulky side/end groups on various positions of the backbone give an interesting variation in the crystal structures and molecular conformations. For instance, in crystals of R,ω-bis(triisopropylsilyl)sexithiophene7 and related compounds of pentamer8 and hexamer8,9 various conformers (or rotamers) can be visualized as a consequence of distortion around the inter-ring σ-bonding. Having these circumstances as a background, we have studied the structural and conformational characteristics of a specially selected set of oligothiophenes in both the crystals and solutions. These compounds are featured by extra bulky terminal groups10 and the following aspects ensue from their presence: (i) The bulky terminal groups introduce relatively large intermolecular steric hindrance between them. At the same time, the presence of these bulky groups results in loose packing of the oligothiophene backbones, which may lead to their conformational variation in the solid state. (ii) In the solutions, in turn, the bulky terminal groups are expected to influence the dynamics of the whole molecule because of, e.g., the large moment of inertia of these bulky groups. The compounds investigated in the present studies are 5,5′′′-bis[(2,2,5,5-tetramethyl-1-aza-2,5-disila-1-cyclopentyl)methyl]-2,2′:5′,2′′:5′′,2′′′-quaterthiophene (1), 5,5′′′-bis[(2,2,5,5-tetramethyl-1-aza-2,5-disila-1-cyclopentyl)ethyl]-2,2′:5′,2′′:5′′,2′′′-quaterthiophene (2), and 5,5′′′-bis[N,N-bis(trimethylsilyl)aminomethyl]-2,2′:5′,2′′: 5′′,2′′′-quaterthiophene (3); see Figure 1 for the chemical structure. Although these compounds have been originally obtained as intermediates for aminoalkyl-terminated oligothiophenes,10,11 the intermediates turned out to be quite stable in the absence of acids. In such an environment the materials crystallize in needles or plates of high quality and can be even vacuum-processed into thin films without causing chemical degradation or alteration.12 We present in this article the results obtained from a coordinated set of experiments. The results of the quantum chemical calculations are also presented. Experimental Section Synthesis of 5,5′′′-Bis[N,N-bis(trimethylsilyl)aminomethyl]-2,2′:5′,2′′:5′′,2′′′-quaterthiophene (3). Syntheses and crystal growths of 1 and 2 are described in the previous report.10 Compound 3 was synthesized as follows: Ethylmagnesium bromide ether solution (3 M, 78.7 mL, 0.236 mol) was added to a solution of 2-(aminomethyl)thiophene (12.4 g, 0.11 mol, Lancaster) in 80 mL of ether and the mixture was (6) Hotta, S.; Rughooputh, S. D. D. V.; Heeger, A. J.; Wudl, F. Macromolecules 1987, 20, 212. (7) Yassar, A.; Garnier, F.; Deloffre, F.; Horowitz, G.; Ricard, L. Adv. Mater. 1994, 6, 660. (8) Liao, J. H.; Benz, M.; LeGoff, E.; Kanatzidis, M. G. Adv. Mater. 1994, 6, 135. (9) Herrema, J. K.; Wildeman, J.; van Bolhuis, F.; Hadziioannou, G. Synth. Met. 1993, 60, 239. (10) Muguruma, H.; Saito, T.; Sasaki, S.; Hotta, S.; Karube, I. J. Heterocycl. Chem. 1996, 33, 173. (11) Muguruma, H.; Saito, T.; Hiratsuka, A.; Karube, I.; Hotta, S. Langmuir 1996, 12, 5451. (12) Muguruma, H.; Hotta, S., to be published.
Muguruma et al.
Figure 1. Chemical structure and conformation of 1, 2, and 3. The backbone conformation is either all-anti (for 1 and 3) or syn-anti-syn (for 2). Ha, Hb, Hc, and Hd indicate the assignments of the protons in the 1H NMR spectra (see Figure 5 and the relevant text). The carbon numbering is supposed to start from the leftmost ring. refluxed for 2 h. Chlorotrimethylsilane (25.7 g, 0.236 mol) in 50 mL of ether was then added, and the solution was refluxed further for 4 h. The resulting solution was washed, dried over K2CO3, and evaporated to yield a pale yellow liquid. This liquid was distilled at 85 °C under reduced pressure 2-3 Torr, giving a viscous colorless liquid of 2-[N,N-bis(trimethylsilyl)aminomethyl]thiophene (4), yield 15.4 g, 55%. 1H NMR (500 MHz, CDCl3) δ (ppm) 7.12 (dd, 1H, 3J5,4 ) 4.9 Hz, 4J5,3 ) 1.2 Hz, 5-H of thiophene); 7.00 (dd, 1H, 3J4,5 ) 4.9 Hz, 3J4,3 ) 3.7 Hz, 4-H of thiophene); 6.81 (dd, 1H, 3J3,4 ) 3.7 Hz, 4J3,5 ) 1.2 Hz, 3-H of thiophene); 4.22 (s, 2H, CH2N); 0.11 (s, 18H, SiCH3). A solution of 4 (5.1 g, 20 mmol) in 80 mL of ether was cooled to 0 °C under dry nitrogen. n-BuLi/hexane (1.6 M, 13.75 mL, 22 mmol, Aldrich) was then added dropwise, causing the solution to turn dark purple. After 30 min, MgBr2‚OEt2 (5.68 g, 22 mmol) was added and allowed to warm to room temperature. After a further 30 min, Ni(dppp)Cl2 (0.16 g, 0.033 mmol) was fed, followed by the immediate addition of 5,5′-dibromo2,2′-bithiophene (2.6 g, 8 mmol), which was prepared beforehand by mixing 2,2′-bithiophene (Aldrich) and N-bromosuccinimide (NBS) in methanol and subsequently recrystallizing the generated white precipitates from acetone.13 The solution was stirred overnight and refluxed for 4 h. The resulting deposits were filtered and recrystallized from hexane to give yellow thin blade-shaped crystals (3), yield 1.9 g, 36%. 1H NMR (500 MHz, CDCl3) δ (ppm) 7.03 (d, 2H, 3J ) 3.7 Hz, ring); 7.01 (d, 2H, 3J ) 3.7 Hz, ring); 6.98 (d, 2H, 3J ) 3.7 Hz, ring); 6.72 (d, 2H, 3J ) 3.7 Hz, ring); 4.19 (s, 4H, CH2N); 0.14 (s, 36H, SiCH3). 13C NMR (125 MHz, CDCl3) δ (ppm) 1.85, 44.45, 123.21, 123.23, 123.56, 123.97, 134.67, 135.36, 136.84, 150.45. UV-vis in CHCl3, λmax 404 nm at room temperature. Anal. Calcd for C30H48N2S4Si4: C, 53.20; H, 7.14; N, 4.14; S, 18.94. Found: C, 53.50; H, 7.09; N, 4.25; S, 18.61. X-ray Crystallography. X-ray data of the single crystals 1-3 were collected using a Mac Science MXC18 four-circle diffractometer with Cu KR radiation at room temperature under ambient environment. Neither decay corrections nor absorption corrections were made. Data analysis was carried out with the CRYSTAN program. The structure was solved by the direct method on a MULTAN and refined by the fullmatrix least-squares method using reflections whose intensities were more than 3 times greater than their standard (13) Kellogg, R. M.; Schaap, A. P.; Wynberg, H. J. Org. Chem. 1969, 34, 343.
Quaterthiophenes Having Extra Bulky Terminal Groups
Chem. Mater., Vol. 10, No. 5, 1998 1461
Table 1. Crystallographic Data for 1, 2, and 3
a
compound
1
2
3
formula formula weight crystal color, habit crystal size (mm) crystal system space group a (Å) b (Å) c (Å) β (deg) V (Å3) Z F(000) Dcalc (g/cm3) λ(Å) µ(Cu KR) (cm) 2θ range (deg) index range h, k, l total reflns unique reflns reflns used in refinement no. of variables Ra Rwb GOF
C30H44N2S4Si4 673.3 yellow, needle 0.45 × 0.45 × 0.15 monoclinic P21/c 16.609(4) 12.408(2) 9.073(1) 91.88(2) 1868.7(6) 2 716 1.20 1.54178 35.80 3-130 0-19, 0-14, (10 3495 3107 2878 248 0.060 0.093 1.27
C32H48N2S4Si4 701.4 yellow, plate 0.35 × 0.30 × 0.15 monoclinic P21/c 19.809(4) 6.158(1) 15.749(2) 97.20(1) 1906.0(5) 2 748 1.22 1.54178 35.26 3-130 0-23, -7-0, (18 3731 3170 2924 263 0.044 0.046 2.47
C30H48N2S4Si4 677.3 yellow, plate 0.35 × 0.30 × 0.20 monoclinic P21/c 17.459(3) 6.3377(8) 17.001(3) 94.62(1) 1875.0(5) 2 724 1.20 1.54178 35.68 3-130 (20, 0-7, 0-19 3686 3127 2887 254 0.055 0.063 1.95
R ) ∑||Fo| - |Fc||/∑|Fo|. b Rw ) [∑w(|Fo| - |Fc|)2/∑w|Fo|2]1/2.
deviations. All the non-hydrogen atoms were refined anisotropically. Spectroscopic Measurements. UV-vis spectra were recorded in a 20 µM toluene solution of the oligothiophene at various temperatures covering -10 to 80 °C with a Ubest-50 (Japan Spectroscopic Co., Ltd.) equipped with a temperature controller using a Peltier effect. The toluene solutions were sealed in a quartz cuvette to ensure against loss of volatile solvent during the course of the measurements. The measurement temperatures may involve errors within (1 °C. 1H NMR spectra were measured at 500 MHz using a JEOL JNM-A500 spectrometer equipped with standard variabletemperature devices. The measurements were carried out either in deuteriochloroform at room temperature or in deuteriotoluene to cover a temperature range -30 to 110 °C. The measurement temperatures may involve errors within (3 °C. 1H COSY spectra were also collected with the same apparatus at room temperature. The sample concentrations were 1540 mM, and the chemical shifts were referred to tetramethylsilane (TMS) for both the 1H NMR and 1H COSY experiments.
Quantum Chemical Calculations. Semiempirical calculations were performed using CACheMOPAC 94 software on a Macintosh Centris 650 equipped with a CXP 40/40 coprocessor. Preliminarily, we tried the PM314 and AM115 methods and compared them by optimizing the molecular geometry, starting from that determined crystallographically. The optimized geometry reached that includes the bond lengths and angles as well as the overall molecular shape is related to the starting geometry with the PM3 method better than with the AM1, and so we chose the PM3 for subsequent calculations.
Results A. Crystallographic Data. Crystal structures: Crystallographic results of 1-3 are collected in Table 1. The atom numbering scheme is shown in Figure 2. (14) Stewart, J. J. P. J. Comput. Chem. 1989, 10, 209. (15) Dewar, M. J. S.; Zoebisch, E. G.; Healy, E. F.; Stewart, J. J. P. J. Am. Chem. Soc. 1985, 107, 3902.
Figure 2. Molecular structure and atomic numbering for 1 (a), 2 (b), and 3 (c). C(10′), S(2′), etc., are given by the inversion of the corresponding C(10), S(2), etc. relative to the molecule center.
The crystals of all the compounds belong to the space group P21/c, and the molecules have the center of symmetry accordingly. For instance, in 1 the coordinates of C(10′), S(2′), etc., are given by the inversion of the corresponding C(10), S(2), etc., relative to the molecular center [i.e., the center of the C(10)-C(10′) bond] which is positioned at the center of the crystal axis of the unit cell (see Figure 3). All the thiophene rings are nearly coplanar, and the methylene carbons directly bonded to the outer thiophene rings are posi-
1462 Chem. Mater., Vol. 10, No. 5, 1998
Muguruma et al.
Figure 3. Stereoviews of the molecular packing for the crystals 1 (a), 2 (b), and 3 (c).
tioned very close to the coplanar backbone plane comprising the four thiophene rings. Figure 3 shows stereoviews of crystals 1-3. The overall crystallographic profile is characterized by a layered structure in which two-dimensionally spreading molecular “sheets” are stacked up along the a axis. The molecular long axis meets with the a axis at an angle of 62.4°, 53.2°, and 59.1° for 1, 2, and 3, respectively. These crystallographic features accompanied by the P21/c (or P21/a) space group are very often observed for crystals of both aromatic2b and aliphatic16 hydrocarbons whose molecules have a well-defined molecular long axis. The bulky terminal groups largely turn aside from the backbone plane. As regards 1 and 2, the leastsquares plane of the five-membered ring of azadisilacyclopentyl in the terminal group and that of the outer thiophene ring cross at an angle 100.3° and 121.8°, respectively. For 3 the plane consisting of the nitrogen and two silicon atoms of the terminal group meets at an angle 107.3° with the least-squares plane of the outer thiophene ring. Owing to the presence of these terminal groups, the molecules glide over one another. As a result, portions covering merely two thiophene rings overlap each other within a pair of adjacent molecules. This is similarly the case with a sexithiophene derivative having bulky triisopropylsilyl groups.7 In 1, furthermore, the bulky terminal groups force the crystals to lack the herringbone structure and the overlapping molecules to be related by the translation.17 The dihedral angle between the two least-squares planes of the overlapping thiophene rings is 3.6° in this arrangement. Another specific feature of 1 is that within the terminal groups 2,2,5,5-tetramethyl-1-aza2,5-disila-1-cyclopentyl, a short intermolecular contact of 3.731(8) Å is observed between two methyl carbons C(18) and C(19). Notice that these two methyl carbons face almost right in front of each other. This contact is significantly smaller than two times van der Waals radius of a methyl group (∼2 Å),18 meaning the presence of strong intermolecular interaction (or steric hindrance) between the terminal groups. The shortest atomic (16) Shearer, H. M. M.; Vand, V. Acta Crystallogr. 1956, 9, 379. (17) Barbarella, G.; Zambianchi, M.; Bongini, A.; Antolini, L. Adv. Mater. 1992, 4, 282. (18) Encyclopaedia Chimica; Kyoritsu Shuppan: Tokyo, 1964; Vol. 7, p 639 (in Japanese).
contact between the molecules is 3.566(3) Å for S(2)‚‚‚C(7). This is pretty close to the sum of the van der Waals radius of each atom (∼1.7 Å for carbon19 and ∼1.8 Å for sulfur19). A bithiophene compound20 with the same terminal groups as 1 shows a virtually identical molecular packing scheme. Molecular structures: The most striking feature associated with the molecular structure is that two distinct nearly coplanar conformations can be noticed (Figures 1 and 2). One is all-anti (for 1 and 3) and the other is syn-anti-syn (for 2). For both conformations, the centrosymmetry in the crystals necessarily requires the two centermost thiophene rings to take the anti conformation relative to the molecular long axis (a line linking the methylene carbons bonded to the thiophene rings). Consequently, the coplanarity encompassing these two rings is rigorously achieved; the deviation of the individual atoms on these rings from the leastsquares plane defined by those atoms is at most 0.007 Å for all the compounds. The outer ring is almost completely planar as well; the deviation of the atoms on the outer ring from its least-squares plane is also at most 0.007 Å. The two least-squares planes meet at an angle of 3.6, 5.0, and 12.1° for 1, 2, and 3, respectively. From the symmetry requirement the said angle for 1 is exactly the same as the dihedral angle between the thiophene rings overlapping in the adjacent pair of molecules. Note that rather unusual S-syn forms are found as disorder or as crystallographically independent components in the compounds also having relatively large substituents.8,21 A “noncoplanar” syn conformer is also observed for a β-methyl-substituted quaterthiophene.22 The selected bond lengths and angles for 1, 2, and 3 are collected in Tables 2, 3, and 4, respectively. The S-C bonds of 1, 2, and 3 exhibit closely related lengths in the ranges 1.723(3)-1.727(3), 1.721(3)-1.732(3), and 1.721(3)-1.735(3) Å, respectively, being contrasted with 1.717(7)-1.750(6) Å for 5,5′′′-dimethyl-2,2′:5′,2′′:5′′,2′′′(19) Bondi, A. J. Phys. Chem. 1964, 68, 441. (20) Muguruma, H.; Hotta, S.; Karube, I. Chem. Lett. 1996, 501. (21) Armes, S. P.; Chaloner, P. A.; Hitchcock, P. B.; Simmons, M. R. Acta Crystallogr. 1994, C50, 1945. (22) Barbarella, G.; Zambianchi, M.; Bongini, A.; Antolini, L. Adv. Mater. 1993, 5, 834.
Quaterthiophenes Having Extra Bulky Terminal Groups
Chem. Mater., Vol. 10, No. 5, 1998 1463
Table 2. Selected Bond Lengths (Å) and Bond Angles (deg) for 1a S1-C14 S1-C8 S2-C10 S2-C12 C6-C8 C6-C7 C7-C14 C14-S1-C8 C10-S2-C12 C8-C6-C7 C14-C7-C6 C6-C8-C12 C6-C8-S1 C12-C8-S1 C10-C9-C13 C9-C10-S2
1.725(3) 1.727(3) 1.723(3) 1.726(3) 1.360(5) 1.425(5) 1.347(4) 92.4(1) 92.4(1) 113.3(3) 113.1(3) 127.8(3) 110.2(2) 122.0(2) 113.5(3) 110.2(2)
C8-C12 C9-C10 C9-C13 C11-C14 C12-C13 C10-C10′ C13-C12-C8 C13-C12-S2 C8-C12-S2 C12-C13-C9 C7-C14-C11 C7-C14-S1 C11-C14-S1 C9-C10-C10′ S2-C10-C10′
1.453(4) 1.363(5) 1.407(5) 1.503(4) 1.354(4) 1.453(6) 129.9(3) 110.5(2) 119.6(2) 113.4(3) 127.4(3) 111.0(2) 121.6(2) 129.9(4) 119.9(3)
a The data are only for the atoms included in the quaterthiophene backbone. This also applies to Tables 3 and 4.
Table 3. Selected Bond Lengths (Å) and Bond Angles (deg) for 2 S1-C13 S1-C7 S2-C9 S2-C12 C7-C8 C7-C9 C8-C10 C13-S1-C7 C9-S2-C12 C8-C7-C9 C8-C7-S1 C9-C7-S1 C7-C8-C10 C11-C9-C7 C11-C9-S2 C7-C9-S2
1.721(3) 1.726(3) 1.723(3) 1.732(3) 1.364(4) 1.455(4) 1.435(5) 92.6(2) 92.7(1) 128.0(3) 110.9(2) 121.0(2) 111.9(3) 127.7(3) 110.3(2) 122.0(2)
C9-C11 C10-C13 C11-C14 C12-C14 C13-C21 C12-C12′ C13-C10-C8 C9-C11-C14 C14-C12-S2 C10-C13-C21 C10-C13-S1 C21-C13-S1 C12-C14-C11 C14-C12-C12′ S2-C12-C12′
1.362(4) 1.350(5) 1.408(5) 1.358(4) 1.505(5) 1.450(5) 114.0(3) 113.2(3) 109.7(2) 129.6(3) 110.5(2) 119.9(3) 114.1(3) 129.7(3) 120.6(3)
Table 4. Selected Bond Lengths (Å) and Bond Angles (deg) for 3 S1-C11 S1-C7 S2-C8 S2-C6 C6-C9 C6-C7 C7-C10 C11-S1-C7 C8-S2-C6 C9-C6-C7 C9-C6-S2 C7-C6-S2 C10-C7-C6 C10-C7-S1 C6-C7-S1 C14-C8-C12
1.725(3) 1.735(3) 1.721(3) 1.732(3) 1.361(5) 1.442(4) 1.364(5) 92.5 92.7(2) 129.9(3) 109.7(2) 120.2(2) 130.0(3) 110.2(2) 119.7(2) 129.5(3)
C8-C14 C8-C12 C9-C14 C10-C13 C11-C13 C11-C11′ C14-C8-S2 C12-C8-S2 C6-C9-C14 C7-C10-C13 C13-C11-S1 C11-C13-C10 C8-C14-C9 C13-C11-C11′ S1-C11-C11′
1.355(5) 1.507(5) 1.409(5) 1.406(5) 1.368(5) 1.451(6) 110.3(2) 120.2(2) 113.6(3) 113.4(3) 110.2(2) 113.8(3) 113.7(3) 129.3(4) 120.5(3)
quaterthiophene (DMQtT)23a in which the bond lengths range considerably widely. Moreover, the bond angles C-S-C of 1, 2, and 3 also show related values 92.4(1), 92.6(2)-92.7(1), and 92.5(2)-92.7(2)°, respectively. These angles are also compared with the corresponding angles 91.4(3)-92.6(3)° for DMQtT.23 These data verify that the ring deformation is smaller in the quaterthiophenes of the present studies than in DMQtT. Seeing that the ring deformation in DMQtT results from the tight packing of the molecules in the crystal lattice,23 the (23) (a) Hotta, S.; Waragai, K. J. Mater. Chem. 1991, 1, 835. (b) Hotta, S.; Waragai, K. Adv. Mater. 1993, 5, 896.
Figure 4. UV-vis spectra of 2 in toluene measured at -10, 0, 10, 20, 40, 60, and 80 °C. The direction of the arrows follows the increasing measuring temperatures.
small ring deformation in 1-3 is ascribable to the loose packing of the oligothiophene backbones. This packing scheme is likely to result from the incomplete overlap between the molecules and to allow various conformers to be stabilized according to a change in the crystallographic environment that is brought about by the subtle modification of the bulky terminal groups. This, at the same time, suggests that the energy level of those conformers is close (vide infra). B. Spectroscopic Data in Solutions. UV-Vis spectra: A typical example of the temperature-dependent UV-vis spectra in toluene is shown for 2 in Figure 4. The spectral profiles of all the compounds are closely related to one another, reflecting the electronic structure of the quaterthiophene backbone. The absorption maxima of those spectra are located around 400 nm. Distinct isosbestic points are observed at 376, 374, and 376 nm for 1, 2, and 3, respectively, within (1 nm accuracy. This demonstrates that two interconvertible conformers are present as the major species under thermodynamic equilibrium. In this case we can readily determine on the basis of the data of Figure 4 the free energy difference (∆G) between the two conformers using a standard process.24,25 We have carried out the numerical computation examining various trial combinations A and B (where A and B are molar absorption coefficients of the conformers A and B, respectively, with A more stable than B) at an appropriately fixed wavelength. For each combination of A and B a standard deviation (∆) of ∆G can be defined as follows:
∑i [(∆Gi - ∆Gh )/∆Gh ]2}1/2
∆){
(1)
where ∆G h is an average of ∆Gi (the subscript i denotes the individual temperature measured). If minimum ∆ (24) Kittel, C.; Kroemer, H. Thermal Physics, 2nd ed.; Freeman: New York, 1980; pp 58-64. (25) Harris, D. C.; Bertolucci, M. D. Symmetry and Spectroscopy; Dover Publications: New York, 1989; pp 307-343, 416-418.
1464 Chem. Mater., Vol. 10, No. 5, 1998
Muguruma et al.
∂T of Hc is larger than that of Hd implies that during the conformers interconversion the chemical environment in the vicinity of Hc changes larger than that for Hd. This may be the origin of the coalescence. Although the coalescence is reported on various organic compounds,29 it is for the first time to the best of our knowledge that the coalescence has been clearly observed regarding the oligothiophene system. We remark in this respect that Blanchard and co-workers27 did report the occurrence of the coalescence for an oligothiophene in their first report but afterward withdrew it as apparent. For the coupled system the interconversion rate constant (k) can be expressed as29
k ≈ (π/x2)(∆νc2 + 6J2)1/2
Figure 5. Temperature-dependent 1H NMR spectra of the aromatic protons of 2 measured in deuteriotoluene. The assignments of these protons are indicated in Figure 1.
is reached after the trial calculations, the ∆G h determined correspondingly to that minimum most likely represents the free energy difference between the two conformers in the temperature range of the measurements. Thus, a reasonable value of ∆G for all the compounds turned out to range from 0.8 to 1.2 kcal/mol, for which the smallest deviations ∆ within a few percent were reached. During the course of the numerical computations, we have found out that at a wavelength longer than ∼390 nm (corresponding to ∼3.18 eV) B/A ≈ 0 gives the most probable solution. This implies that the less stable conformer B has an absorption edge at ∼3.18 eV, even though the absorption edge in the real spectra occurs at 2.58 eV (480 nm) owing to the more stable species A (Figure 4). The less stable species have an increasing contribution in the spectra with increasing temperatures, causing a blue-shift in the absorption maxima. 1H NMR spectra: The temperature-dependent spectra of the aromatic protons in deuteriotoluene are represented for 2 in Figure 5. The assignment of the peaks has been done by COSY spectra and on the basis of the pervious assignment for analogous compounds.26 Figure 5 shows that the coupled quartet peaks due to Hc and Hd cause a downfield shift with increasing temperatures with the difference in the chemical shift (∆ν) between these protons smaller. These peaks finally coalesce at 65 °C. Although the corresponding peaks of 3 also coalesce at 107 °C, those peaks of 1 do not coalesce within the experimental range of the temperatures (up to 110 °C). Since the chemical shift is strongly affected by the vicinal chemical environment,27,28 the fact that ∂(∆ν)/ (26) Albers, W. M.; Canters, G. W.; Reedijk, J. Tetrahedron 1995, 51, 3895.
(2)
where ∆νc is the chemical shift difference between the coupled protons at an instant of the coalescence and is approximated by the fwhm (full width at half-maximum) of the coalesced line; J is their coupling constant. The applicability of eq 2, however, is valid only if ∆νc exceeds the coupling constant.29d In cases where ∆νc is well below the coupling constant, the complete lineshape analysis29d shows that the rate constant is rather approximated by a simpler form:29a
k ≈ (π/x2)∆νc
(3)
this expression being usually used for the uncoupled protons system. Since ∆νc of 2 and 3 are 1.33 and 2.52 Hz, respectively, and J is 3.66 Hz, we presume that the rate constant is between those evaluated from eqs 2 and 3 (see Figure 2 of ref 29d). It is well-established that the free energy of activation for the conformers interconversion (∆Gq) is expressed by30
∆Gq ) 2.303 × 10-3RTc log(kBoltzmannTc/hk)
(4)
where Tc is the coalescence temperature, kBoltzmann the Boltzmann constant, and h the Planck constant. We assume the transmission coefficient to be unity.30,31 Since in our case the rate constant for the process A to B (kA) and that for B to A (kB) will not be identical, the individual constants may be defined using the averaged constant k obtained either in eq 2 or in eq 3 as32 kA ) 2PBk and kB ) 2PAk where PA and PB are populations (27) (a) Horne, J. C.; Blanchard, G. J.; LeGoff, E. J. Am. Chem. Soc. 1995, 117, 9551. (b) DeWitt, L.; Blanchard, G. J.; LeGoff, E.; Benz, M. E.; Liao, J. H.; Kanatzidis, M. G. J. Am. Chem. Soc. 1993, 115, 12158. (28) Jackman, J. M.; Sternhell, S. Applications of Nuclear Magnetic Resonance Spectroscopy in Organic Chemistry, 2nd ed.; Pergamon Press: Oxford, 1969; p 312. (29) (a) Oki, M. Applications of Dynamic NMR Spectroscopy to Organic Chemistry; VCH Publishers: Deerfield Beach, FL, 1985; pp 1-12. (b) Oki, M.; Iwamura, H.; Hayakawa, N. Bull. Chem. Soc. Jpn. 1964, 37, 1865. (c) Kurland, R. J.; Rubin, M. B.; Wise, W. B. J. Chem. Phys. 1964, 40, 2426. (d) Kost, D.; Carlson, E, H.; Raban, M. J. Chem. Soc., Chem. Commun. 1971, 656. (e) Kost, D.; Zeichner, A. Tetrahedron Lett. 1974, 4533. (30) Sandstro¨m, J. Dynamic NMR Spectroscopy; Academic Press: London, 1982; pp 93-123. (31) Glasstone, S.; Laidler, K. J.; Eyring, H. The Theory of Rate Processes; McGraw-Hill: New York, 1941; p 13. (32) (a) Kaplan, F.; Meloy, G. K. Tetrahedron Lett. 1964, 2427. (b) Kaplan, F.; Meloy, G. K. J. Am. Chem. Soc. 1966, 88, 950. (c) ShananAtidi, H.; Bar-Eli, K. H. J. Phys. Chem. 1970, 74, 961.
Quaterthiophenes Having Extra Bulky Terminal Groups
Chem. Mater., Vol. 10, No. 5, 1998 1465 Table 5. Enthalpy Data (kcal/mol) and the Point Group for Various Rotamers ∆Hq (from ∆Hq (from no. of axis of point syn)c anti)d compound axes distortion groupa ∆Hb 1 1 2 3 1 2 3 1
1 1 1 1 2 2 2 3
5′,2′′ 2,2′ 5′,2′′ 5′,2′′ 2,2′; 5′′,2′′′ 2,2′-; 5′′,2′′′ 2,2′-; 5′′,2′′′ all bonds
C2 C1 C2 C2 Ci Ci Ci C2
0.5 0.4 0.4 0.5 0.9 1.1 1.0 1.3
0.7 0.6 0.7 0.7 1.3 1.4 1.3 1.9
0.2 0.2 0.3 0.2 0.4 0.3 0.3 0.6
a The point group is referred only to the quaterthiophene backbone. b Enthalpy difference between the two minima. c Barrier enthalpy starting from the syn form. d Barrier enthalpy starting from the anti form.
Figure 6. Enthalpy of the oligothiophenes as a function of the torsional angle θ around the 2,2′ and 5′′,2′′′ axes. The centermost rings are held at the anti conformation and the backbones retain the Ci symmetry after the inter-ring distortion.
of A and B, respectively. The constant k in eq 4 then must be replaced with either kA or kB. The populations can be determined on the basis of the electronic spectra (Figure 4). Thus, we obtained as the free energy of activation ∼17.3-19.0 (18.2-20.2) kcal/mol for the species B (A) of 2 and ∼19.6-21.0 (20.5-22.1) kcal/mol for the species B (A) of 3. For this estimation both the variation in the free energy difference (0.8-1.2 kcal/mol) and an error in the measuring temperatures ((3 °C) were taken into account. C. Quantum Chemical Calculations. There exist many rotamers according to torsional angles around the axes interconnecting the thiophene rings. The energy (E) of these rotamers can be expressed as follows as a function of these angles:
E ) E(θ1,θ2,θ3)
(5)
where θ1, θ2, and θ3 represent the torsional angles around the 2,2′, 5′,2′′, and 5′′,2′′′ axes, respectively (see Figure 1 where the carbon numbering is supposed to start from the leftmost ring). After optimizing the molecular geometry, we allowed only the inter-ring rotational freedom within the quaterthiophene backbone. The torsional angle is defined as a dihedral angle between the two planes determined by the thick lines drawn in Figure 6 such that this angle is zero at the syn position. We carried out the calculations, permitting various molecular symmetries and changing the number of axes that participate in the rotamers interconversion. The molecular symmetry was referred only to the quaterthiophene backbone.
In all the above cases, the energy (represented by enthalpy in the figure) of the rotamers shows a double minima as a function of the torsional angle θ, these minima being found at θ ) 30° (syn-periplanar or syngauche) and 150° (anti-periplanar or anti-gauche) with the former more stable than the latter by ∼0.4-1.3 kcal/ mol. The rotamers are destabilized most when each torsional angle is 90°. An example is shown in Figure 6, and the related results are summarized in Table 5. The relative stability between the rotamers and the barrier energies depend principally upon the number of the axes that participate in the inter-ring distortion regardless of the molecular symmetry. This is because the distortion around the individual bond gives rise to a similar level of energy instability33 of ∼0.4-0.5 kcal/ mol per bond participating in the ring distortion. Therefore, if two (three) such bonds are involved in the distortion, the barrier energies and their difference are expected to be doubled (triplicated).33 The energy difference between the two minima is expected to represent the free energy difference between the conformers; the barrier energies correspond to the free energy of activation for the conformers interconversion. Note for the latter case that the energy of 1.9 kcal/mol (Table 5) sets an upper limit to the barrier energy, because during the course of the interconversion all the three angles θ1, θ2, and θ3 do not need to assume 90° all at once. Discussion On the basis of the structural and spectroscopic results collected in both the solids and solutions, we discuss their physicochemical implications referring to the results of the quantum chemical studies. Although the vibronic structure reflecting the presence of the well-defined molecular symmetry23,25,34,35 is always inconspicuous in the solution spectra, the absorption edge observed at 2.58 eV for the more stable conformer is pretty close to that for the thin films (∼2.5 eV)12 where the molecules assume a planar conformation. This implies that this conformer in the solution (33) (a) Belleteˆte, M.; Ce´sare, N. D.; Leclerc, M.; Durocher, G. Chem. Phys. Lett. 1996, 250, 31. (b) dos Santos, D. A.; Galva˜o, D. S.; Laks, B.; dos Santos, M. C. Chem. Phys. Lett. 1991, 184, 579. (34) Herzberg, G. Molecular Spectra and Molecular Structure, III. Electronic Spectra and Electronic Structure of Polyatomic Molecules; D. Van Nostrand: Princeton, 1967; pp 128-157. (35) Wilson, E. B., Jr.; Decius, J. C.; Cross, P. C. Molecular Vibrations; Dover Publications: New York, 1980; pp 146-156.
1466 Chem. Mater., Vol. 10, No. 5, 1998
is quasiplanar despite causing symmetry lowering. The occurrence of an absorption edge at a considerably short wavelength in the other conformer is due to the interring distortion that breaks the π-conjugation.5,6 In this context, the conformers we are dealing with are rotamers. We conclude from this that the more stable rotamer is more extended and that the less stable one is more twisted or coiled. Consequently, the conformers interconversion may be characterized as the quasiplanar-coiled one where the symmetry of both the rotamers is lowered to a greater or lesser extent. It should be emphasized that only two rotamers are present as the major species, as evidenced by the UV-vis spectra, despite the possibility of the presence of various conformers that correspond to a variety of combinations of θ1, θ2, and θ3 giving the energy minima. Our conclusion is consistent with the results of the IR measurements by Herna´ndez et al.36 carried out in solution. They treat the conformers of the oligothiophenes in terms of a mean conjugation length (MCL), and they conclude that two conformers exist in the solution regarding a related material, DMQtT. One has MCL near 4 and the other has MCL between 2 and 3. Although the quantum chemical calculations reproduce the free energy difference between the conformers pretty well, the barrier heights separating the conformers obtained from these calculations are admittedly too low to account for the free energy of activation for their interconversion. The calculated barrier height of at most 1.9 kcal/mol is contrasted with the free energy of activation estimated to be ∼20 kcal/mol from the 1H NMR. This made us assume the presence of the origins that the quantum chemical calculations ignored. This may be legitimated, e.g., by the inference by Ortı´ et al.37 that the experimentally determined barrier height for bithiophene38 is considerably larger than the theoretical estimates because of impediment from surrounding liquid-crystalline media. One such extrinsic origin is the rotational diffusion of the molecules in solution that is based upon the interaction between those molecules as the solute and solvent molecules. This model deals with reorientation in liquid media of nonspherical molecules that are considered as rigid rotors, in terms of the rotational diffusion tensor.39,40 The tensor is associated with the moment of inertia. More precisely, after diagonalized the diffusion tensor components with respect to the principal axes are approximately inversely proportional to the three principal moments of inertia.39 Although the free energy of activation for the rotational diffusion of the molecules was not expressed explicitly as a function of the tensor components, studies on cyclohexene by Pajak et al.40 showed a trend to which a larger rotational diffusion tensor is associated with a smaller free energy of activation for the molecular reorientation. If this is generally the case, the rotational diffusion around the principal axis relevant to the smallest (36) Herna´ndez, V.; Casado, J.; Ramı´rez, F. J.; Alemany, L. J.; Hotta, S.; Lo´pez Navarrete, J. T. J. Phys. Chem. 1996, 100, 289. (37) Ortı´, E.; Viruela, P. M.; Sa´nchez-Marı´n, J.; Toma´s, F. J. Phys. Chem. 1995, 99, 4955. (38) Bucci, P.; Longeri, M.; Veracini, C. A.; Lunazzi, L. J. Am. Chem. Soc. 1974, 96, 1305. (39) Assink, R. A.; Jonas, J. J. Phys. Chem. 1969, 73, 2445. (40) Pajak, Z.; Latanowicz, L.; Jurga, K. Ber. Bunsen-Ges. Phys. Chem. 1980, 84, 769.
Muguruma et al. Table 6. Moment of Inertia of the Oligothiophenes and Free Energy of Activation for the Interconversion of the Rotamersa
compound 1 2 3 2,2′-bithiopheneb
2,2′ axis 5′,2′′ axis (×10-46 kgm2) (×10-46 kgm2) 413 357 165 16
274 1497 680
free energy of activation (kcal/mol) 17.3-20.2 19.6-22.1 ∼5b
a The moment of inertia is only for the left part relative to the axis (see text). b Reference 38.
principal moment of inertia will be first activated with the smallest free energy of activation. Following this consideration, we infer that the above applies with the rotamers of the oligothiophenes. By intuition, the moment of inertia around the molecular long axis of the oligothiophenes is thought to be by far smaller than the others; think of the small separations of the individual atoms from this axis. Therefore, the rotation about this axis (more exactly, the 2,2′ or 5′,2′′ axis) is expected to be most easily activated. Note, however, that the interconversion of the rotamers differs from the conventional rotational diffusion in that during this process the molecules can be regarded not as the rigid rotors but as a flexible system having inner rotational freedom around the chemical bonds. Here let us tentatively deal with the molecular motion assuming that the left portion relative to the 2,2′ axis (or 5′,2′′ axis) causes the rotational diffusion, while the remaining right portion is permitted to remain stationary as a whole relative to a surrounding liquid medium. In this situation, we must consider both the energy given to activate the rotational diffusion and that for breaking the π-conjugation. To make the treatment easier, we calculated the moments of inertia around the individual chemical bonds connecting the thiophene rings, limiting the freedom only to the rotation around either the 2,2′ or 5′,2′′ axis and borrowing the relevant geometrical data from the crystals. The results are shown in Table 6 along with the free energies of activation for the rotamer interconversion determined from the coalescence data. The moments of inertia are only for the left portion relative to the axis. Again, Table 6 shows the trend to which the larger moment of inertia is apparently associated with the larger free energy of activation. Furthermore, the large difference between the moments of inertia around the 2,2′ axis and the 5′,2′′ axis for 2 and 3 is likely to be responsible for the coalescence. This is because the rotational diffusion around the 5′,2′′ axis is hindered because of the large moment of inertia around this axis, leading to smaller change in the chemical environment in the vicinity of Hd compared to that for Hc. Conclusion The structures and molecular conformations of the oligothiophenes having the extra bulky terminal groups have been investigated by X-ray analysis and UV-vis and 1H NMR spectroscopic methods. These bulky terminal groups make the structural and conformational features of the oligothiophenes stand out among the related compounds. The crystals are characterized by
Quaterthiophenes Having Extra Bulky Terminal Groups
the molecular layered structure. In these crystals, the molecules of compounds 1 and 3 have the all-anti conformation, whereas those of the compound 2 have the syn-anti-syn form. The conformational variation in the crystals produced by the subtle modification in the terminal groups suggests that the relevant conformers have a close energy level. We have studied the relative stability between them through the UV-vis spectroscopy carried out in the solutions. The results show that two interconvertible conformers are present under thermodynamic equilibrium with the free energy difference between the two ∼0.8-1.2 kcal/mol. The detailed numerical analysis of the spectral data implies that of the two conformers one exhibits an absorption edge at ∼480 nm (2.58 eV), while that of the other is located around 390 nm (3.18 eV). The latter conformer has an increasing contribution in the spectra with raised temperatures, causing the blueshift in the absorption maximum. These spectroscopic results indicate that the two conformers are characterized as the rotamers. One is more extended, and the other is more distorted. This is consistent with the results of the IR measurements done in solution.36 These rotamers can be best described quantitatively in terms of the torsional angles around the bonds connecting the thiophene rings. Within this framework, we carried out the quantum chemical calculations based on the PM3 method. These calculations show that the rotamers exhibit energy minima when they assume the syn-periplanar and anti-periplanar conformations, the former being more stable than the latter. The energy difference between the two was estimated to be about 0.4-0.5 kcal/mol per bond participating in the ring distortion. This reproduces appreciably well the level
Chem. Mater., Vol. 10, No. 5, 1998 1467
of the free energy difference between the two conformers estimated from the numerical analysis of the UV-vis spectra in the solutions. The dynamic methods of the 1H NMR spectroscopy offer important information about the rate processes of the interconversion between the conformers. Analyzing widths of the lines at the instant of the coalescence that come from the coupled protons Hc and Hd, we estimated that the free energy of activation for the interconversion ranges from roughly 17 to 22 kcal/mol. This is significantly larger than the energies obtained by the quantum chemical calculations. The discrepancy is indicative of the presence of the origin that was disregarded in the quantum chemical calculations. We sought this origin from the rotational diffusion of the molecules in the solutions and found out that this model interprets the large free energy of activation pretty satisfactorily. Acknowledgment. We would like to thank Prof. Isao Karube for his encouragement to our studies. We would like to thank Dr. Kazushi Kinbara for his technical assistance. The present studies were partly supported by grants to the Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists (H.M.). Supporting Information Available: Tables of positional parameters, equivalent isotropic thermal parameters, bond lengths and angles, symmetry operations, and the relevant crystallographic information for 1-3 and temperature-dependent 1H NMR spectra of aromatic protons in deuteriotoluene and temperature-dependent UV-vis spectra in toluene for compounds 1 and 3 (34 pages). See any current masthead page for ordering and Internet access instructions. CM980041I