Key words
sequence-defined macromolecules - shape-defined macromolecules - Sonogashira coupling - iterative exponential growth - density functional theory - π-conjugation
Introduction
Poly(p-phenylene ethynylenes) (PPEs)[2] are special in that they belong (together with other poly(p-phenylene)s)[3] to a class of conjugated, non-ladder polymers, whose backbone structures are fully shape-defined. This shape persistence arises from the fact that, with fully linear triple bonds, rotation around any of the single/triple bonds present in the backbone of a PPE does not change a PPE's overall end-to-end distance. With other macromolecules, a similar degree of shape persistence can only be obtained by introducing rings into the backbones, for example, with chirality-assisted synthesis.[4]
Due to their linearity, high degree of π-conjugation, and associated electronic communication between the different phenylene units, π-conjugated macromolecules[5] (including PPEs) have found applications in the fields of sensing,[3b]
[3e]
[3j]
[6] organic electronics,[3h]
[7] and biological imaging.[2d]
[3j]
[8] However, open questions still remain as to how the geometrical and photophysical properties of PPEs are affected by electron-withdrawing substituents[3g]
[9] like ester groups. Initial studies in this regard have focused on ester-functionalized PPE systems with a distribution of different lengths.[10] However, a size distribution in chain lengths can make it difficult to correlate the detailed photophysical properties with chain length, since the spectra are naturally broadened due to the inherent length distributions present in each sample. While more challenging to access,[11] the study of unimolecular macromolecules offers valuable additional information, in particular as to how absorption linewidths are affected by conformational disorder.[12] Yet, the prior literature investigating[13] unimolecular models of PPEs has been focused primarily on unsubstituted and/or alkoxy-substituted oligo(p-phenylene ethynylenes) (OPEs), which behave quite differently from ester-functionalized OPE systems, as detailed in the Results and Discussions section. Here we now synthesized some of the largest, unimolecular, ester-functionalized OPEs with iterative exponential growth[14] (IEG).
Results and Discussion
We started this work with density functional theory (DFT) calculations[15] to predict the exact geometries of the triple bonds in OPEs with various substituents. To account for dispersion interactions, we utilized the B3LYP-MM functional[16] with the cc-pVDZ++ basis set (for single-point calculations) and the LACVP* basis set (for geometry optimizations). The B3LYP-MM parameters were carefully optimized with a large dataset of non-covalent interaction energies to accurately reproduce dispersion interactions, even in the presence of basis set superposition error.[16e]
As a simple model for ester-functionalized OPEs, we utilized tetramethyl 2,2'-(ethyne-1,2-diyl)diterephthalate (1) for our DFT analysis. To our surprise, we noted that in the optimized structure of 1 (lowest energy structure in vacuum, see [Figure 3B] for an alternative low-energy conformation), the triple bonds are bent ([Figure 1A]), with C–C≡C angles of 171.7°. This finding is in stark contrast to the larger C–C≡C angles of 180.0° and 177.9°, which we observed at the same level of theory for the corresponding unsubstituted (1,2-diphenylethyne) as well as for methoxyl-substituted (1,2-bis(2,5-dimethoxyphenyl)ethyne) systems (see the Supporting Information for the optimized structures). We then calculated the critical points of the electron density and used them to visualize ([Figure 1A]) the non-covalent interactions that are primarily responsible for the bending of the triple bonds with the NCI code[17] implemented in the Jaguar[18] software package. The NCI critical points, which were calculated from the electron density ([Figure 1B]), clearly demonstrate the presence of attractive, supramolecular interactions between the carbonyl groups of the ester groups, and the triple bonds. This result is consistent with triple-bond bending, driven by carbonyl-to-alkyne electron donation. Further experimental evidence for these interactions stems from a published[19] crystal structure ([Figure 1C]) of a model compound (dimethyl 2,2'-(ethyne-1,2-diyl)bis(3-(2-((t-butoxycarbonyl)amino)-propanamido)benzoate), which also shows the bent triple bonds arising from the carbonyl-to-alkyne interactions).
Figure 1
A. DFT-optimized tetramethyl 2,2'-(ethyne-1,2-diyl)diterephthalate (1) as a model for an OPE repeat unit. The DFT-optimized structure (lowest energy conformation in vacuum, see [Figure 3B] for an alternate low-energy conformation) illustrates how the triple bonds in the OPEs bend due to carbonyl-to-alkyne electron donation effects. NCI critical points, calculated with the Jaguar software package from the electron density (see Panel B), are illustrated with blue spheres. As has been established by Johnson et al. see Ref. ([17] these NCI critical points represent attractive supramolecular interactions (NCI interaction strength = 9.0 kcal mol−1 in vacuum and 9.1 kcal mol−1 in CHCl3 with a PBF solvent model). B. DFT-calculated electron density (isosurface at 0.014 a.u.) of tetramethyl 2,2'-(ethyne-1,2-diyl)diterephthalate (1). Arrows indicate the enhanced sections of the electron density, which correspond to the attractive supramolecular interactions between the carbonyl groups and the alkyne units of 1 identified by the NCI analysis shown in Panel A. C. Single-crystal X-ray structure of a model compound (dimethyl 2,2'-(ethyne-1,2-diyl)bis(3-(2-((t-butoxycarbonyl)amino)-propanamido)benzoate), reported in Ref. [19], CCDC 915930), which clearly shows the bent triple bonds arising due to carbonyl-to-alkyne interactions.
Next, we discovered that the carbonyl-to-alkyne interactions also significantly alter the barriers for rotation around the triple bonds in the OPEs. Notably, we found ([Figure 2]) the barrier for rotation around the triple bond in the ester-functionalized model system tetramethyl-2,2'-(ethyne-1,2-diyl)diterephthalate (1) to be nearly twice as high as in the unsubstituted model system, 1,2-diphenylethyne (2). This finding is explained by the carbonyl-to-alkyne interactions, which desymmetrize ([Figure 2C]) the two orthogonal π-bonds of the alkynes in the backbone of the OPEs.
Figure 2
A. Carbonyl-to-alkyne electron donation effects on the torsional profiles (triple bond rotation) of 2,2'-(ethyne-1,2-diyl)diterephthalate (1) and 1,2-diphenylethyne. Torsional profiles were calculated at the B3LYP-MM/cc-pVDZ + +//B3LYP-MM/LACVP* level of theory in vacuum. B. Definition of the dihedral angles ω used as the abscissa for the torsional plots. C. Carbonyl-to-alkyne electron donation effects lead to angled triple bonds with non-degenerate π-bonds.
Based on these computational results, which demonstrate the unique geometrical and conformational properties of ester-functionalized OPEs, we next set out to synthesize such macromolecules in a unimolecular fashion. As shown in [Scheme 1], the synthesis of up to ∼10 nm long, unimolecular OPEs was accomplished with Sonogashira coupling-based IEG synthesis.[20]
Scheme 1 Iterative convergent/divergent (a.k.a. iterative-exponential growth – IEG) synthesis of unimolecular oligo(p-phenylene ethynylenes) (OPEs), substituted with up to 32 ester functional groups. Conditions for activation (see the Supporting Information for details): i) diazotization: NaNO2, HCl, H2O, CH3CN, toluene (used for the synthesis of the longer oligomers 4, 5, and 6 to enhance solubility), 0° C. ii) Iodination: KI, 0° C. Conditions for triisopropylsilyl (TIPS) deprotection (see the Supporting Information for details): iii) tetrabutylammonium fluoride (TBAF), CH2Cl2, room temperature. Conditions for Sonogashira couplings (see the Supporting Information for details): iv) Pd(PPh3)4 (3 mol%), CuI (6 mol%), NEt3, DMF, 70 °C. Notably, all intermediates, including the unprotected acetylenes with free amino groups, are air stable. Furthermore, the carbonyl groups of the ester substituents assist the oxidative addition step of the Sonogashira couplings.
The starting material, bis(2-ethylhexyl)-2-amino-5-((triisopropylsilyl)ethynyl)terephthalate (3), was synthesized as detailed in the Supporting Information. Briefly, 2-amino-5-iodo-1,4-benzenedicarboxylic acid (synthesized as described previously in the literature)[21] was deprotonated with potassium carbonate, and the resulting bis(carboxylate) derivative alkylated with 3-(bromomethyl)heptane. Sonogashira coupling with (triisopropylsilyl)acetylene then afforded bis(2-ethylhexyl)-2-amino-5-((triisopropylsilyl)ethynyl)terephthalate (3) as the starting point for IEG growth.
As detailed in [Scheme 1], IEG growth of the OPEs then consisted of three simple steps, which were applied iteratively.[22] (i) About half of the triisopropylsilyl (TIPS)-protected sample at each growth stage is deprotected with tetrabutylammonium fluoride (TBAF) to afford the terminal alkyne derivative (which can directly engage as the alkyne donor in a Sonogashira cross-coupling). (ii) The other half of the sample is then activated to become the alkyne acceptor for the Sonogashira coupling step by converting the terminal aniline group into an aryl iodide, via a one-pot diazotization/iodination[23] reaction sequence. (iii) Finally, the aryl iodide component is linked to the component containing the free acetylene group in a Sonogashira coupling step to double the chain length. Notably, the presence of the electron-withdrawing ester groups along the OPE backbone renders all the intermediates (including the unprotected acetylenes) fully air-stable. Related intermediates for alkoxy-substituted OPE derivatives can show air-sensitivity for increased polymer lengths, a challenge[24] which seems to be completely avoided by our ester-functionalized backbones.
With this IEG approach we were able to isolate the OPEs 4–7.[25]
[26]
[27]
[28]
[29]
7 represents, to the best of our knowledge, the longest fully ester-functionalized OPE synthesized to date. With 7 in hand, we set out to investigate the effects of the ester substituents on the photo-physics of the unimolecular OPEs. The UV/Vis absorption spectrum of 7 displayed ([Figure 3A]) a similar absorption maximum (at 405 nm) as the previously reported[30] heptadecameric OPE 8, which contains both unsubstituted phenylenes and alkoxy-substituted phenylene units. Interestingly, however, the UV/Vis absorption spectrum of 7 was clearly broadened ([Figure 3A]), compared to the UV/Vis absorption spectrum of 8. Significant line-broadening was also observed (see [Supplementary Figures S15] and [S16]) for the shorter oligomers 5 and 6.
Figure 3 Enhanced conformational disorder resulting from carbonyl-to-alkyne-derived alkyne-bending contributes to broadening of the UV/Vis absorption spectrum of the hexadecamer 7. A. Comparison of the UV/Vis absorption spectrum (CHCl3) of 7 to the UV/Vis absorption spectrum (CHCl3) of the mixed unsubstituted/alkoxy-substituted OPE 8. The UV/Vis absorption data for 8 were extracted with the WebPlotDigitizer from Ref. [30]. Broadening of the UV/Vis absorption spectrum for the OPE 7 with the ester groups is observed, compared to the OPE 8, which lacks the ability to engage in carbonyl-to-alkyne electron donation. B. Alternative low-energy conformation of the OPE model compound 1. Relative energy (relative to the conformation shown in [Figure 1A]) in vacuum, E
rel, vacuum = 0.05 kcal mol−1. Relative energy with the CHCl3 Poisson–Boltzmann finite element (PBF) solvent model implemented in Jaguar, E
rel, CHCl3 = –0.38 kcal mol−1. We hypothesize that this secondary low-energy conformation (which is accessible to each of the 15 internal triple bonds of the ester-functionalized OPE 7) contributes to the observed ([Figure 3A]) line broadening of the UV/Vis absorption spectrum of 7. The relative energies were calculated at the B3LYP-MM/cc-pVDZ + +//B3LYP-MM/LACVP* level of theory. Similar line broadening is observed for the shorter oligomers 5 and 6 (see [Supplementary Figures S15] and [S16]).
DFT calculations (performed like before at the B3LYP-MM/cc-pVDZ++//B3LYP-MM/LACVP* level[16] of theory to account for dispersion interactions as well as basis set superposition error) were able to explain the observed broadening of the UV/Vis absorption spectrum caused by the ester substituents on the OPE 7. Specifically, the DFT calculations showed that, due to the bent triple bonds, an alternate low-energy conformation exists for each triple bond. This secondary low-energy conformation is only 0.05 kcal mol−1 higher in energy than the most stable conformation (shown in [Figure 1]) in vacuum, while it becomes slightly favored in energy when applying a Poisson–Boltzmann finite element (PBF) solvent model in CHCl3. In general, solvation favors the alternate low-energy conformation shown in [Figure 3B], as it possesses a larger dipole moment (3.4 debye) than the conformation shown in [Figure 1A] (0.0 debye). The alternate conformation of the triple bond, in which both ortho-ester groups are located on the same side of the triple bond, is stabilized by [C–H…O]-hydrogen bonds,[31] which are shown as critical points of the electron density (blue spheres) in [Figure 3B]. Given the small energetic differences between these two very distinct conformations, we conclude that both of these low-energy conformations very likely coexist in solution for each of the 15 triple bonds in 7, which induces significant conformational disorder, and associated conformational line broadening of the UV/Vis absorption spectrum of 7.[32]
Conclusions
We investigated the effects of electron-withdrawing ester functions on the geometry and photophysical properties of OPEs, for the first time with unimolecular, ester-functionalized OPEs up to ∼10 nm in length (with up to 16 repeating units). We demonstrated that – in contrast to unsubstituted and alkoxy-substituted OPEs – the triple bonds in the ester-functionalized OPEs have a tendency to bend, departing from their idealized, fully linear conformations. The observed bending of the triple bonds in the ester-functionalized OPEs is driven by carbonyl-to-alkyne electron donation effects, which also increase the rotational barriers around the triple bonds, and lead to enhanced conformational disorder and broadening of the UV/Vis absorption spectra. Our results advance the fundamental understanding of how the geometrical and associated photophysical properties of unimolecular, π-conjugated oligomers and polymers can be tuned with electron-withdrawing ester substituents. We are currently utilizing our new unimolecular, ester-functionalized macromolecules as templates for polymer replication, as well as for sensing applications.