Jürgen
Rotzler
a,
Heiko
Gsellinger
a,
Markus
Neuburger
a,
David
Vonlanthen
a,
Daniel
Häussinger
*a and
Marcel
Mayor
*ab
aDepartment of Chemistry, University of Basel, St. Johannsring 19, 4056 Basel, Switzerland. E-mail: marcel.mayor@unibas.ch; Fax: +41-61-267-1016; Tel: +41-61-267-1006
bInstitute of Nanotechnology, Karlsruhe Institute of Technology, P. O. Box 3640, 76021 Karlsruhe, Germany
First published on 18th November 2010
The thermodynamics of the atropisomerisation of torsion angle restricted, axial chiral biphenyl-based push-pull cyclophanes were studied. Using 1H NMR coalescence measurements the rotation barrier around the central C–C bond was determined to be 50 kJ mol−1 for the propyl-bridged biphenyl derivative 1b, displaying only a negligible solvent dependence. By protonation of the piperidinyl nitrogen as electron donor, the free energy ΔG‡(T) of the rotation barrier increased, indicating that the tendency of the push-pull system to planarise may be considered as a driving force for the atropisomerisation. For the more restricted butyl-bridged cyclophane 1c a rotation barrier of ΔG‡(T) = 90 kJ mol−1 was measured using dynamic chromatography. The difference in the free energy of rotation around the central C–C bond probably reflects the crowdedness of the transition states.
Recently, studies about torsion angle restricted biphenyl-based push-pull systems were published where the influence of the chromophore’s delocalization on the hyperpolarizability was investigated by EFISH measurements.23 The investigated 2 and 2′ position alkyl-bridged biphenyl-based cyclophanes 1a–1d are axial chiral and therefore consist of a racemic mixture of two atropisomers (Fig. 1). In crystallography, it is well known that optically pure enantiomers tend to crystallise in noncentrosymmetric point groups, whereas racemic mixtures especially of donor-π-acceptor molecules crystallise in a centrosymmetric fashion since their dipoles force such an arrangement (Fig. 2). Therefore the described push-pull systems 1a–1d are potential candidates to realise such a transformation of the noncentrosymmetry from a molecular level to the macroscopic scale by enantiomeric resolution and subsequent crystal growth.
Fig. 1 (A) Sketch of the concept to restrict the interphenyl torsion angle by an additional interlinking alkyl chain of various lengths; (B) The two atropisomers of 1c (structures calculated with a MM2 basic set) and the rotation around the central C–C bond showing the atropisomerisation process. |
Fig. 2 Centrosymmetric pairwise arrangement of the P (right) and M (left) atropisomers of the push-pull derivative 1c in the solid state. |
To gain further insight into the conformational stability of the non-linear active cyclophanes 1a–1d and to check if it is possible to separate the two atropisomers, the rotation barrier around the central C–C bond moved into the focus of interest.
For such thermodynamic studies mainly three different experimental methods are applicable. Coalescence measurements by 1H NMR,24 dynamic chromatography25–29 and monitoring the change in the optical rotation of a chiral non-racemic mixture. For the latter, separated enantiomers are required, whereas the first two methods can easily be performed using racemic mixtures of the molecules under investigation. Fast dynamic molecular processes can be investigated by NMR spectroscopy, whereas dynamic chromatography is more suitable for isomerisation studies of more conformationally stable compounds. Herein, the enantiomerisation dynamics of torsion angle restricted biphenyl push-pull cyclophanes 1a–1d are reported to obtain further insight into their dynamic behaviour and the interaction of the two phenyls under these conditions.
Scheme 1 (a) Benzophenone imine, Pd2(dba)3·CHCl3, BINAP, NaOtBu, toluene, 80 °C, 4 h, then 3% aq. HCl, THF, rt, 2 h, 80%–quant.; (b) benzophenone imine, Pd(OAc)2, BINAP, Cs2CO3, THF, 65 °C, 17 h, then 3% aq. HCl, THF, rt, 2 h, 80%–quant.; c) 1,5-dibromopentane, K2CO3, toluene–EtOH: 1/1, MW, 150 °C, 40 min, 34–44%; (d) NaBO3·4H2O, H3PW12O40, CTAB (hexadecyltrimethylammonium bromide) (10 cmc in water), 55–60 °C, 16 h, 45–64%. |
Fig. 3 Solid state structure of the P enantiomer of 1c. Thermal ellipsoids are set at the 50% probability level. |
ΔG = 0.0191Tc(9.97 + log(Tc/Δν)) | (1) |
The coalescence temperature (Tc) of the push-pull system shown in Fig. 4 (1b) was measured in different solvents to obtain changes in the activation energy due to solvent effects. As the chemical shifts for 1b can be solvent dependent, the coalescence temperature must also show a variation with the solvent, even in the case of unchanged activation parameters. It is, however, possible that differences in the solvation of the sterically demanding transition state of the racemisation lead to differences in the activation energy. The results of the measurements of 1b in different solvents are shown in Table 1.
Fig. 4 Stacked plot of the low temperature 1H NMR spectra of 1b measured in CDCl3. Starting from fast exchange at 298 K, reaching the coalescence temperature at 252 K, ending in the slow exchange at 233 K where the signals for the four protons are separated. |
It was possible to record NMR spectra below and above as well as exactly at the coalescence temperature for most solvents used. Only the melting point of benzene was too high, so that further cooling was not feasible. For toluene, two Tc values for two different sets of resonances were observed. The two different data sets lead to the very same activation energy, as expected. As seen in Table 1, the influence of different solvents on the activation energy is not very pronounced. The polarity of the solvent does not seem to significantly influence the activation energy for the interconversion of the enantiomers. Only in the case of strongly protic conditions (trifluoroacetic acid), a significantly higher activation energy was obtained. This is most likely due to the protonation of the amine group, leading to a pull-pull, rather than a push-pull substitution pattern. In a pull-pull system, the tendency to form a conjugated and, therefore, nearly planar arrangement will be less pronounced compared to the push-pull case, where this resonance stabilisation factor determines the racemisation energy.34 Additionally, two electron-withdrawing groups reduce the π-electron density at the 1 and 1′ positions, which makes an out-of-plane bending of the axis bond more difficult and therefore becomes the dominating contribution to the energy barrier.35,36 In turn, a higher activation is necessary to force the cyclophane to racemisation.
In the case of the ethyl-bridged biphenyl 1a, it was not possible to lower the temperature far enough (<230 K) to reach the coalescence in any solvent. In marked contrast, for compounds 1c and 1d there was no fast exchange regime obtained, even when heating was continued to 353 K in deuterated tetrachloroethane, consistent with findings of earlier studies on unsubstituted biphenyl cyclophanes.37,38 In addition, 2D-EXSY experiments at the elevated temperature indicated that the rate constant of the atropisomerisation has to be below 10−2s−1.39
As mentioned above, it was not possible to determine the rotation barrier of the C–C bond of the butyl- and pentyl-bridged biphenyl derivatives 1c and 1d by coalescence measurements using 1H NMR spectroscopy. Careful analysis of the 1H NMR spectra clearly showed a fixed torsion angle for both push-pull systems. Another useful method to estimate isomerisation energies in cases where NMR spectroscopy fails is dynamic chromatography.25–29 As previously described for such investigations only minute sample amounts are required and stereoisomers do not have to be isolated prior to the analysis. A chromatogram with two separated peaks is expected in the case of a fixed torsion angle of the compounds under investigation when they are separated on a chiral stationary phase. When racemisation occurs on the chromatographic time scale a temperature dependent plateau formation between the two peaks is expected. A third case that can occur in dynamic HPLC is peak coalescence. This occurs when the enantioconversion is significantly faster than the chromatographic separation process.40 Hence, we tried to separate the two atropisomers of 1c and 1d by chiral HPLC. Screening of several cellulose- and amylose-derived HPLC columns revealed that the enantiomers of 1c showed a tendency to separate on Chiralpak AD-H. At 20 °C the two peaks could be separated by approximately 2 min when the compound was eluted with a mixture of n-hexane–iPrOH 97:3 at a flow rate of 0.5 mL min−1. Furthermore, plateau formation between the two peaks was observed pointing at the expected interconversion process of atropisomers (Fig. 6). The UV/vis spectra of the two distinctive peaks in the chromatogram and the plateau were the same (Fig. 5), indicating the presence of only one compound and therefore verifying the presence of two atropisomers. The rate constant of enantioconversions can be directly calculated from these typical chromatograms by iterative comparison of experimental and simulated chromatograms.
Fig. 5 Representative chromatogram of push-pull cyclophane 1c, separated by a Chiralpak AD-H column eluting with a mixture of n-hexane–iPrOH 97:3 and a flow rate of 0.5 mL min−1 at 25 °C. In the third dimension the UV/vis spectra of the separated parts are shown. |
Fig. 6 Selected experimental enantiomerisation profiles of push-pull cyclophane 1c at different temperatures on Chiralpak AD-H using a mixture of n-hexane–iPrOH 97:3 as eluent. |
Thanks to the remarkable work of Schurig27 and Trapp,33,41 powerful computer simulation tools are now available. The simulation programs profit from the so-called “theoretical plate model”,26 its statistical description namely the “stochastic model”41 or more recently the “unified equation”32 to determine the rate constants of enantiomerisation. Temperature dependent dynamic chromatography gives access to thermodynamic data by applying the Eyring equation (eqn (2)).
k1 = (kBT)/h exp(−ΔG/RT) | (2) |
Thus, temperature dependent dynamic HPLC experiments were performed for compound 1c. Repeated series of measurements were performed between 10 °C and 25 °C in 5 °C steps. To overcome reproducibility problems the column was preconditioned for 10 h operating at the desired conditions. After two series the conditioning was repeated and the measurements were then performed again with a different sample of the same compound. A statistical evaluation was then performed. Peak coalescence or baseline separation were not observed in this temperature range. By computer simulations, using the DCXplorer software,33 of the experimentally obtained elution profiles (Fig. 6) the enantiomerisation rate constants at different temperatures were obtained (Table 2). Since enantiomerisations are degenerate processes k1 and k−1 were assumed to be equal.
With the unified equation obtained, enantiomerisation rate constants of 1c were plotted according to the Eyring equation (ln(k1/T) versus 1/T). By linear regression analysis (Fig. 7), ΔH‡ was found to be ΔH‡ = 54.5 ± 3.9 kJ mol−1 and ΔS‡ = −120 ± 14 J K−1 mol−1. From this, ΔG‡ was calculated to be ΔG298 K‡ = 90.3 ± 0.2 kJ mol−1, which is in the region of 2-phenyl-2′-isopropylbiphenyl and 2-cyclohexyl-2′-phenylbiphenyl (both approx. 91 kJ mol−1) but above the rotational energy barrier of 2,2′-dimethylbiphenyl (approx. 78 kJ mol−1) and below that of 2,2′-diisopropylbiphenyl (approx. 110 kJ mol−1) reported by C. Wolf.42
Fig. 7 Eyring plot for the atropisomerisation of 1c: temperature dependent dynamic HPLC measurements analysed by linear regression. The upper and lower curves represent the error bands of the linear regression with a level of confidence of 95%. For the linear regression, four measurements for each temperature were considered. |
Unfortunately, such an investigation for the pentyl-bridged biphenyl-based push-pull system 1d was not possible by dynamic chromatography as coalescence was observed even at low temperatures. Thus, an enantioconversion faster than the chromatographic separation process can be assumed. Since 1H NMR coalescence measurements showed a fixed conformation of the two phenyls of compound 1d, it is obvious that the free energy of rotation is between 50 kJ mol−1 obtained for the 7-membered ring derivative 1b and 90 kJ mol−1 obtained for the 8-membered cyclophane 1c.
By comparison of the rotational energy barriers of the torsion angle restricted biphenyl-based push-pull cyclophanes 1a–1d it is obvious that the butyl-bridged system 1c is the conformationally most stable compound of the series followed by the pentyl-bridged derivative 1d. Compound 1c atropisomerises with a half life time τ1/2 of 12.5 min which means that this compound can be assumed to have a fixed torsion angle in most physical investigations such as EFISH measurements. For the propyl-bridged push-pull system 1b this is not the case since the energy barrier for atropisomerisation is relatively low and therefore, depending on the time frame of the experiment, isomerisation can take place during measurements which would lead to a distribution of different torsion angles and therefore to an averaged torsion angle over the course of the experiment. The same is true for the ethyl-bridged derivative 1a which is even more flexible than 1b.
To explain this rather unexpected grading in enantiomerisation energies, we hypothesise a sterically more crowded transition through a planar, linear conformation in the isomerisation process. Linearity of the biphenyl core can be assumed since partial double-bond character of the central C–C bond was observed in the crystal structure of the 2 and 2′ position unsubstituted push-pull system. This means that the longer the alkyl-bridge, the closer the methylene groups have to pass by each other in the planar conformation. The pentyl-bridge in system 1d seems to have too many degrees of freedom and hence it can arrange in a sterically more favourable fashion to allow isomerisation in an energetically more convenient way.
Unfortunately, the separation of the two atropisomers failed for all non-linear optically active compounds 1a–1d, and thus it was not possible to grow noncentrosymmetric crystals. Nevertheless, we obtained insight in the conformational stability of our systems 1a–1d. For compound 1a, the rotation around the central C–C bond is rather fast and so the system has to be considered to be flexible. Also, the rotation around the central axis of the propyl-bridged push-pull cyclophane 1b is quite fast at room temperature. Depending on the physical experiment, compound 1c and 1d can be assumed as torsion angle fixed systems when the experimental timescale is in the region of seconds. Nevertheless, the maximal possible torsion angle was defined by bridging the biphenyl core with alkyl chains. Motivated by these results, it will now be of interest to determine the rotational energy barrier of 2,2′-alkyl-bridged and 4 and 4′ position donor-substituted biphenyl systems and 4,4′-diacceptor-substituted biphenyl systems which are available in our group and have already been investigated in single molecule conductance measurements.30,43,44
Footnote |
† Crystal data and structure refinement for cyclophane 1c: The crystal was measured on a Bruker Kappa Apex2 diffractometer at 123 K using graphite-monochromated Mo-Kα-radiation with λ = 0.71073 Å, Θmax = 27.655°. Minimal/maximal transmission 0.99/0.99, μ = 0.083 mm−1. 1c formula C21H24N2O2, M = 336.43 g mol−1, F(000) = 1440, yellow block, size 0.070 0.130 0.290 mm3, monoclinic, space groupC2/c, Z = 8, a = 25.437(6) Å, b = 10.510(2) Å, c = 17.615(5) Å, α = 90°, β = 132.55(2)°, γ = 90°, V = 3469.1(18) Å3, Dc = 1.288 Mg m−3. The Apex2 suite has been used for data collection and integration. From a total of 21978 reflections, 4017 were independent (merging r = 0.111). From these, 2392 were considered as observed (I > 2.0σ(I)) and were used to refine 226 parameters. The structure was solved by direct methods using the program SIR92. Least-squares refinement against F was carried out on all non-hydrogen atoms using the program CRYSTALS. R = 0.0533 (observed data), wR = 0.1511 (all data), GOF = 1.0519. Minimal/maximal residual electron density = −0.31/0.33 e Å−3. Chebychev polynomial weights were used to complete the refinement. Plots were produced using Mercury 2.3 and ORTEP3 for Windows. |
This journal is © The Royal Society of Chemistry 2011 |