Martin
Dračínský
*a,
Jan
Storch
b,
Vladimír
Církva
b,
Ivana
Císařová
c and
Jan
Sýkora
*b
aInstitute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo 2, 16610 Prague, Czech Republic. E-mail: dracinsky@uochb.cas.cz
bInstitute of Chemical Process Fundamentals, Czech Academy of Sciences, Rozvojová 135, 16502 Prague, Czech Republic. E-mail: sykora@icpf.cas.cz
cFaculty of Science, Charles University in Prague, Hlavova 2030, 12840 Prague, Czech Republic
First published on 3rd January 2017
The conformational behaviour of a prototype helical molecule, [6]helicene, and its derivatives was studied in solution and the solid state. Available crystal structures of [6]helicene revealed surprisingly large flexibility of this molecule and variable-temperature NMR experiments provided unusual temperature dependence of chemical shifts of hydrogen, carbon and fluorine atoms in the peripheral aromatic rings of [6]helicene and tetrafluoro[6]helicene. These chemical shift changes were interpreted as a consequence of the helicene ‘pitch’ opening with elevated temperature, and the experimental data were corroborated by DFT calculations of the chemical shift dependence on the helicene conformation and by variable-temperature DFT molecular dynamics simulations.
The highly delocalised large π-electron system of fully aromatic helicenes along with the previously mentioned inherent chirality provide them with unique optical and electronic properties. The helicenes achieve e.g. exceptionally high values of specific rotation (thousands of degrees per mole).2 They also provide strong circular dichroism,3 electroluminescence,4,5 circularly polarised luminescence,6,7 and non-linear optical properties.8–11 Furthermore, helicenes exhibit semi-conductive behaviour.12,13 All of these properties predetermine their utilisation in a wide range of applications in organic electronic devices.14–16 The chirality of helicenes has also been utilised in the chiral discotic liquid-crystalline materials, or these molecules have served as ligands in asymmetric catalysis.17–19 In addition, they can be used as building blocks in functional polymers.20 Undoubtedly, these valuable properties are closely connected with the overall structure of the molecule.
Unsubstituted helicenes usually crystallise in non-centrosymmetric space groups (P21, P212121) despite the fact that the helicenes are racemic.21–25 This phenomenon involves so-called ‘lamellar twinning’, where the layers of individual enantiomers alternate within a single crystal, forming homo-chiral domains.26 Some of the derivatives have a tendency towards co-facial stacking in bulk state forming macroscopic fibres consisting of columnar aggregates. Similar self-assembly motifs can be also found in Langmuir–Blodgett films prepared by dip or spin coating.27,28 Due to the perpendicular arrangement of the aromatic rings with respect to the surface, the deposited films are usually anisotropic.29
Starting from [6]helicene, the twisted molecule of helicenes is highly constrained due to the steric hindrance of the peripheral rings, which gives the impression of rigidity. However, a comparison of accessible crystal structures of [6]helicene reveals surprising flexibility within a certain range of the molecular conformation. On the one hand, there is a crystal structure of a charge-transfer complex of [6]helicene with trimesic acid, which has provided the most compressed conformation of the [6]helicene molecule with the distance between C2 and C2′ atoms of the peripheral rings of 3.94 Å (Fig. 1).30 On the other hand, the helicene conformation can be stretched by the insertion of another molecule or atom into the helicene pitch, such as in the crystal structure of plain [6]helicene (the C2–C2′ distance of 4.58 Å in the original structure from 1973),21 which consists of pairs of interlocked helicene molecules, or a [6]helicene complex with silver salts, where the silver ion is bound by η1 bond to C2 and C2′ atoms, giving the C2–C2′ distance of 4.61 Å.31 The adopted conformation is a consequence of intra- and intermolecular interactions in the crystal lattice and induces deformations of individual aromatic rings.
Fig. 1 The crystal structure of (a) [6]helicene, (b) [6]helicene-trimesic acid (benzene-1,3,5-tricarboxylic acid), (c) Ag–[6]helicene complex, (d) [6]helicene–perfluorobenzene complex and (e) tetrafluoro[6]helicene. Complete atom numbering is shown in the ESI.† |
This paper presents an original study of internal dynamics in [6]helicenes both in solution and the solid state with special attention to the helicene pitch opening. It analyses the temperature dependence of the structural parameters determined by X-ray diffraction and of NMR chemical shifts. The experimental results are subsequently confronted with DFT calculations and molecular dynamics simulations. It is shown that a combination of experimental and theoretical approaches allows a precise description of the conformational behaviour of helicenes.
Diffraction data were collected at 150 K on a Nonius Kappa CCD diffractometer (Enraf-Nonius) with graphite-monochromated Mo-Kα radiation. A Cryostream Cooler (Oxford Cryosystem) was used for low-temperature measurements. The structures were solved by direct methods (SIR92,33 SHELXL9734) and refined by full-matrix least-squares on F2 values (CRYSTALS35). All heavy atoms were refined anisotropically. Hydrogen atoms were localised from the expected geometry and difference electron density maps and were refined isotropically. PARST9736 was used for the calculation of ring-puckering parameters.37 The details of atom numbering and the complete list of puckering parameters are listed in the ESI.† ORTEP-338 was used for structure presentation.
NMR spectra were measured at room temperature on a Bruker Avance III spectrometer operating at 600.1 MHz for 1H and at 150.9 MHz for 13C in DMSO-d6 and aceton-d6. A combination of 1D and 2D experiments (COSY, HSQC, HMBC) was used for the assignment of all 1H and 13C resonances. Variable-temperature NMR experiments were performed on a Bruker Avance II spectrometer operating at 499.9 MHz for 1H, 470.3 MHz for 19F and at 125.7 MHz for 13C. Temperature was calibrated with methanol and ethyleneglycol; a 200–400 K temperature range with 50 K steps was used for the experiments. To minimise the solvent effects and temperature dependence of bulk magnetic susceptibility on the reported variable-temperature NMR data,39 relative chemical shift dependences on temperature with respect to a selected signal are reported, i.e. the chemical shift of zero ppm was attributed to the selected signal at every temperature in both solvents. Signals of hydrogen H8, carbon C8 and fluorine F4 were selected for the referencing because the signals of these atoms are the least affected by the through-space ring-current effects of peripheral aromatic rings.
Two different approaches were used for DFT computations. First, geometry optimisation of isolated molecules with B3LYP40,41 and PBE42 functionals and standard Pople-type basis sets were performed with the Gaussian09 program package.43 The NMR parameters were calculated using the GIAO method and B3LYP functional. Second, geometry optimisation and molecular dynamics (MD) simulations of fully periodic systems and of isolated molecules were run in the CASTEP program,44 which exploits periodic plane waves as the basis-set functions. DFT molecular dynamics (MD) simulation of crystalline [6]helicene under periodic boundary conditions was done using an NVT ensemble maintained at a constant temperature of 300 K using a Langevin thermostat, a 0.5 fs integration time step, total simulation length of 20 ps, ultrasoft pseudopotentials,45 a planewave cutoff energy of 500 eV, and with integrals taken over the Brillouin zone using a Monkhorst–Pack46 grid of a minimum k-point sampling of 0.1 Å−1. Electron-correlation effects were modelled using the generalised gradient approximation of Perdew, Burke and Ernzerhof.42 The atomic positions were optimised at the same computational level prior to the MD runs, while lattice parameters were fixed to the experimental values. No symmetry constraints were applied during the runs as these are only relevant to the time-averaged structure. The same procedure was used for MD simulations of isolated [6]helicene at 200 and 300 K modelled as an isolated molecule in a cubic periodic box of 20 × 20 × 20 Å3. The effect of empirical corrections for dispersion47 was also examined both in the simulations of isolated molecules and crystals.
As mentioned above, the overall conformation of the helicene molecule is reflected in the distance between the peripheral rings. Extreme values were found in the complex of [6]helicene with trimesic acid (C2–C2′ distance of 3.94 Å) and, surprisingly, in the structure of plain [6]helicene (4.64 Å at 100 K). Both values describe conformations severely influenced by intermolecular interactions: compression produced by the grid of the crystallisation partner (trimesic acid) on the one hand and the insertion of the second helicene molecule into the pitch on the other. The insertion of a silver cation induced a similar opening of the pitch to 4.62 Å. The two extreme values gave a surprisingly large conformational range of 0.7 Å. The relaxed helicene molecule in the structure of the charge-transfer complex with perfluorobenzene provided a C2–C2′ distance of 4.31 Å, which is close to the middle of the observed conformational range (4.29 Å). On the other hand, the C2–C2′ distance found in the relaxed structure of tetrafluoro[6]helicene is shorter (4.14 Å); the π–π interaction of the electron-deficient with the electron-rich aromatic ring leads to a significant shortening of the distance between the peripheral rings.
A detailed conformational analysis of individual aromatic rings based on a comparison of standard ring-puckering parameters37 has revealed the central aromatic rings as the most distorted from planarity, while the peripheral rings remain almost planar. The total puckering amplitude Q ranges from 0.01 to 0.07 Å for the peripheral rings 1 and 6, from 0.11 to 0.16 Å for the rings 2 and 5, and from 0.18 to 0.21 Å for the central rings 3 and 4. According to standard visualisation of the ring-puckering parameters, the conformations of all aromatic rings are located near the ‘equator’, adopting relatively narrow range of twisted or boat conformations 2T4–B4,1–3T1 for right-handed helix (Fig. 2) and 1T3–1,4B–4T2 for the left-handed helix. For a detailed explanation and the list of puckering parameters, see the ESI.† The range of the individual benzene ring conformations is similar to that found in previous computational studies of benzene and other aromatic compounds.48–52
Fig. 2 The visualisation of ring-puckering parameters calculated for the right-handed helix. Only the most relevant part of the conformational space is depicted. For full details see the ESI.† |
As the main focus of this paper is on the internal dynamics in these helical molecules, a variable-temperature X-ray diffraction study of [6]helicene in the temperature range of 100–300 K has also been performed. The analysis has shown that the higher temperature increases only the oscillations of individual atoms manifested by larger thermal ellipsoids without any significant effect on molecular conformation. The C2–C2′ distance even slightly decreases with the elevated temperature, which may look surprising at first sight; however, this is only a consequence of higher amplitudes of librational motion, which leads to shorter distances between average atomic positions.53 A similar shortening of the C2–C2′ distance has also been observed in the structure of the [6]helicene–trimesic acid complex measured at 105 K and room temperature.30 Nevertheless, the restricted motion of the [6]helicene molecule fixed in the crystal lattice is in contrast to the prolongation of the cell parameters. Especially, the cell parameters b and c were extended by 0.21 and 0.14 Å, respectively, which corresponds to 1.6% and 0.8% of their value, respectively. The elongation of the parameter a was significantly smaller (0.04 Å, 0.5%). A detailed analysis of the crystal packing revealed a noticeable change only in the intermolecular distance between the two interlocked molecules, which increased by 0.14 Å (measured as the distance between molecular centroids) in the temperature range studied. This effect takes place mainly in the direction [012] and therefore strongly affects solely the cell parameters b and c in the corresponding ratio. The increased thermal atomic motion results in the ousting of the second molecule from the helicene pitch rather than in any observable conformational change (Fig. 3).
Fig. 3 The ousting of the second molecule from the helicene pitch in the structure of [6]helicene. The structure measured at 100 K in black, at 300 K in grey. |
Chemical shifts of hydrogen and carbon atoms in an unsubstituted peripheral ring of [6]helicene and tetrafluoro[6]helicene molecules are shown in ESI.† The largest 1H chemical-shift differences in the studied compounds were found for the H1 atom, which is found at 7.43 ppm in the spectrum of unsubstituted [6]helicene and at 7.01 ppm in the spectrum of tetrafluoro[6]helicene. The hydrogen atom H1 is the closest atom to the opposite peripheral aromatic ring in helicenes, as a result of which it is the most affected by the ring-current effects. The fact that the H1 in tetrafluoro[6]helicene is more shielded indicates that its helix is more ‘squeezed’ and the H1 is closer to the opposite fluorinated aromatic ring. The described chemical-shift difference can mostly be attributed to through-space effects as the fluorine atoms are separated by at least eight covalent bonds. On the other hand, the chemical-shift differences of the hydrogen atoms H2, H3 and H4 in [6]helicene and tetrafluoro[6]helicene are rather small. The slightly higher values of these chemical shifts in tetrafluoro[6]helicene can be explained by the slightly lower electron density in this ring caused by electron-density transfer to electron-withdrawing fluorine substituents on the other side of the molecule. Similar chemical-shift changes were also found in 13C NMR spectra. The carbon atoms C1, C1a and C2 in [6]helicene are significantly less shielded than the corresponding atoms in tetrafluoro[6]helicene.
The flexibility of helicene molecules found in various crystal structures might be reflected in large conformational fluctuations of the isolated molecules. Therefore, variable-temperature NMR measurements in a broad temperature range (200–400 K) have been performed to gain insight into the dynamics of the studied helicenes. In the absence of any structural changes, chemical shifts usually exhibit only weak temperature dependence caused by vibrational averaging.54–56 The temperature dependence of 1H and 13C chemical shifts in [6]helicene is shown in Fig. 4. The atoms H1 and C1 in this molecule have revealed significant temperature dependence leading to higher chemical shifts at higher temperature. This behaviour may be explained by increased molecular dynamics at higher temperatures, leading to a larger average distance between the opposite peripheral aromatic rings and hence to a decrease of the effects of the ring currents. The effect of temperature on H2, H3 and H4 chemical shifts is significantly lower than in the case of H1, and the rest of the hydrogen chemical shifts are almost temperature-independent. Similar behaviour was also observed for 13C resonances, where C1 shielding is significantly more temperature-dependent than the shielding of other carbons. The 1H, 13C and 19F chemical shifts in tetrafluoro[6]helicene have revealed similar temperature dependence (ESI†). The atoms that are close to the opposite peripheral ring are significantly deshielded at a higher temperature, while the temperature effect on other chemical shifts is much less pronounced. For comparison, the temperature dependence of chemical shifts of phenanthrene (a reference planar aromatic compound without through-space shielding by ring-current effects) is also shown in the ESI.†
Fig. 4 The temperature dependence of relative 1H chemical shifts with respect to H8 (left) and of 13C chemical shifts with respect to C8 (right) in [6]helicene measured in acetone and DMSO. For atom numbering see Fig. 1 and ESI.† |
Geometry optimisation under fully periodic conditions performed in the CASTEP provided excellent agreement of the optimised geometries with X-ray structures. The inclusion of empirical corrections for dispersion did not change the resulting geometry much, probably because the crystal packing itself forces the molecules to adopt the preferred conformation and neglecting dispersion leads to error cancelation of intramolecular and intermolecular weak interactions.
On the other hand, when geometry optimisation is performed at the same computational level but for isolated molecules, the resulting conformation is significantly different from the X-ray structure and the inclusion of dispersion significantly affects the conformation (Table 1). The structure of tetrafluoro[6]helicene optimised in vacuum with dispersion correction is almost identical to that found in the crystalline state, while [6]helicene in the crystal is ca. 0.4 Å more ‘open’ than in vacuum, which can be explained by the intermolecular interactions discussed above (the interlocking of helicene molecules). The influence of the functional and basis set used for the geometry optimisation was also examined (see details in the ESI†).
Compound | Exp. | Calc. | |||
---|---|---|---|---|---|
X-ray | Periodic | Periodic + D | Isolated | Isolated + D | |
[6]Helicene | 4.622 | 4.627 | 4.622 | 4.460 | 4.198 |
F4[6]helicene | 4.141 | 4.152 | 4.159 | 4.254 | 4.153 |
The intramolecular potential along the C2–C2′ distance was obtained by geometry optimisations of both compounds with the distance fixed at selected values (Fig. 5). These calculations confirm the results found in fully relaxed structures, i.e. the dispersion correction leads to a shortening of the C2–C2′ distance and [6]helicene is more ‘open’ than tetrafluoro[6]helicene. The potential is, however, relatively flat in both compounds (an increase of the distance by 0.4 Å leads to an increase of potential energy by only 0.5 kcal mol−1). At room temperature, Boltzmann thermal energy quantum kBT corresponds to ca. 0.6 kcal mol−1, hence the distance may significantly fluctuate around the equilibrium value.54 These small energy penalties for the deformation of the helicene molecule may be easily paid by intermolecular interactions in crystals, which explains the conformational variability of helicenes in the solid state.
Fig. 5 The calculated relative energy dependence on the C2–C2′ distance in isolated [6]helicene. For tetrafluoro[6]helicene, see the ESI.† |
Chemical shifts of [6]helicene were also calculated for partially optimised structures with a fixed C2–C2′ distance (used above for the construction of the intramolecular potential). The distance dependence of hydrogen chemical shifts is shown in Fig. 6, which clearly resembles the experimental temperature dependence of relative hydrogen chemical shifts in [6]helicene (Fig. 4 left). The chemical shift of H1 is the most distance dependent with higher shifts observed for larger distances. This observation thus confirms the influence of ring-current effects on nuclear shielding, and the similarity with the experimental temperature dependence of chemical shifts confirms the expected trend in the helix opening with increasing temperature. At higher temperatures, the helix is more open, the effect of ring currents is weaker, and the chemical shift of H1 increases. Through a comparison of the calculated chemical-shift dependence on the C2–C2′ distance with the experimental temperature dependence, the opening of the helix in an isolated [6]helicene molecule can be estimated to be ca. 0.2 Å per 100 K.
The temperature dependence of experimental chemical shifts in [6]helicene and tetrafluoro[6]helicene indicates significant internal dynamics in these helical molecules. In order to explore this internal dynamics in detail, molecular dynamics simulations of a fully periodic crystal structure of [6]helicene at 300 K and of isolated [6]helicene at 200 and 300 K were performed.
The distance found with the highest probability in the periodic simulation (4.6 Å) is very close to the experimental C2–C2′ distance found in the crystal (Fig. 7). The crystal packing restricts the flexibility of individual [6]helicene molecules, which leads to a relatively sharp and symmetric distribution of the distance probability around the most probable value.
On the other hand, the helix is less open in the isolated molecule of [6]helicene than in the crystalline state; this is reflected in a shift of the probability function in the case of the isolated molecule to a smaller distance than in periodic simulations, both performed at 300 K. The probability distribution function is also slightly broader, which can be explained by the lack of crystal-packing forces. The distance probability distribution calculated for isolated molecules is significantly affected by the temperature of the simulation. With increasing temperature, the distribution function broadens and moves towards higher values, which reflects larger amplitudes of motion and a larger average C2–C2′ distance at higher temperatures. The most probable C2–C2′ distance increases by ca. 0.2 Å when going from 200 to 300 K, which is in excellent agreement with the estimate made from the comparison of experimental and calculated chemical-shift changes. Average C2–C2′ distance calculated throughout the whole MD simulations increased by 0.1 Å.
DFT calculations have confirmed high flexibility of [6]helicene conformation; the intramolecular ‘opening’ potential is surprisingly flat with the minimum at 4.2 Å, and the 0.4 Å stretching/squeezing penalty being well below 1 kcal mol−1, which can be easily overcome by intermolecular interactions in solid materials. Furthermore, DFT calculations have confirmed the trend of increasing chemical shifts of peripheral atoms in conformations with a larger opening of the helicene molecule. A comparison of the experimental chemical-shift dependence on temperature with the calculated chemical-shift dependence on the helicene opening has made it possible to estimate that the helicene opening in isolated molecules increases by approximately 0.2 Å per 100 K. A similar result was provided by molecular dynamics simulations of isolated molecules, where the most populated C2–C2′ distance increased by 0.2 Å between the simulations performed at 200 and 300 K.
The experimental and theoretical investigation performed within this work describes the unexpectedly large conformation flexibility and dynamic behaviour of [6]helicene molecules. These properties should be taken into account in the design of helicene-like structures for future applications.
Footnote |
† Electronic supplementary information (ESI) available: Additional experimental and computational data. CCDC 1514731–1514735. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6cp07552e |
This journal is © the Owner Societies 2017 |