Tristan J.
Keller§
a,
Christopher
Sterzenbach§
a,
Joshua
Bahr
a,
Taria L.
Schneiders
a,
Markus
Bursch
b,
Julia
Kohn
b,
Theresa
Eder
c,
John M.
Lupton
*c,
Stefan
Grimme
*b,
Sigurd
Höger
*a and
Stefan-S.
Jester
*a
aKekulé-Institut für Organische Chemie und Biochemie, Rheinische Friedrich-Wilhelms-Universität Bonn, Gerhard-Domagk-Str. 1, 53121 Bonn, Germany. E-mail: hoeger@uni-bonn.de; stefan.jester@uni-bonn.de
bMulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Bering str. 4, 53115 Bonn, Germany. E-mail: grimme@thch.uni-bonn.de
cInstitut für Experimentelle und Angewandte Physik, Universität Regensburg, Universitätsstr. 31, 93053 Regensburg, Germany. E-mail: john.lupton@physik.uni-regensburg.de
First published on 9th June 2021
Molecular spoked wheels with D3h and Cs symmetry are synthesized by Vollhardt trimerization of C2v-symmetric dumbbell structures with central acetylene units and subsequent intramolecular ring closure. Scanning tunneling microscopy of the D3h-symmetric species at the solid/liquid interface on graphite reveals triporous chiral honeycomb nanopatterns in which the alkoxy side chains dominate the packing over the carboxylic acid groups, which remain unpaired. In contrast, Cs-symmetric isomers partially allow for pairing of the carboxylic acids, which therefore act as a probe for the reduced alkoxy chain nanopattern stabilization. This observation also reflects the adsorbate substrate symmetry mismatch, which leads to an increase of nanopattern complexity and unexpected templating of alkoxy side chains along the graphite armchair directions. State-of-the-art GFN-FF calculations confirm the overall structure of this packing and attribute the unusual side-chain orientation to a steric constraint in a confined environment. These calculations go far beyond conventional simple space-filling models and are therefore particularly suitable for this special case of molecular packing.
When one of the acetylene ends contains the aforementioned 3,4,5-trihexadecyloxybenzyloxy unit (bulky corner in Scheme 1), the processability of the MSWs is maintained even when the remaining corner positions of the final product are widely varied. In case of carboxylic groups, giant size homologues of trimesic acid11a,b and also trimellitic acid (benzene 1,2,4-tricarboxylic acid)18 are available and allow a comparison between the packing behavior of the respective esters and acids. The latter generally dimerize intermolecularly at the surface and stabilize the adsorbate packing. However, it is not predictable if the bulky flexible side groups allow this stabilization also in the present case or if instead the alkyl packing dominates the pattern formation.
Scheme 1 Schematic synthesis of molecular spoked wheels (MSWs) with reduced symmetry via Vollhardt trimerization of asymmetrically substituted acetylenes. |
All compounds, as well as the aggregation of 1a in dichloromethane, were fully characterized in solution by NMR and PL lifetime measurements (for details see ESI‡).19 Since the wheel precursors as well as the wheels 1a and 1b show identical NMR spectra, the assignment of the respective isomers is thus far based only on their different polarities (Rf values). However, a final proof of the structure of the wheels 1a/b and 2a/b is obtained by scanning tunneling microscopy (STM, Fig. 2). All four compounds form self-assembled monolayers (SAMs) at the solid/liquid interface of highly oriented pyrolytic graphite (HOPG) and a dilute (1–3 × 10−6 M) solution using octanoic acid (OA) as a solvent. The choice of solvent was optimized in thorough previous work on different MSWs.20 Initial studies on 1a in OA (see ESI‡) showed that thermal treatment of the adlayers at sufficiently high temperatures (of 100–110 °C) is required for the formation of 2D nanopatterns visible by STM at ambient temperature. This requirement is not uncommon for such large molecules and is related to the high compound affinity to the substrate and the thermal activation barrier of 2D diffusion,10e as well as kinetic trapping of molecules in amorphous film regions, where the molecules retain some residual mobility and therefore reduced visibility in the STM. The monolayers were investigated with the STM tip immersed into the solution phase (for details see ESI‡). Bright and dark image regions indicate low and high tunneling resistivities, respectively. For all MSWs dark regions in the central hub and at the corners are observed, probably due to a sterically driven electronic decoupling of these molecule parts from the substrate.21
1a forms domains (of about 100 × 100 nm2 size, see ESI‡) with a chiral 2D crystalline honeycomb packing (Fig. 2a, left part) that is a union of two triangular sublattices (cf. bright and dark blue colored molecules, Fig. 2e and i). The structure is triporous, i.e. a nanopattern containing pores of three different shapes or sizes, while it is biporous with respect to the intermolecular pores, and each molecule inherits another six smaller intramolecular triangular pores. The domains are separated by uncovered regions (see ESI‡) or nonporous domain boundaries (Fig. 2a, right part).22 A unit cell with parameters a = b = (8.2 ± 0.2) nm and γ(a,b) = (60 ± 2)° containing two molecules with identical conformations but different (relative) orientations cA and cB with γ(cA,cB) = 180° (or 60°/120°, Fig. 2j) is indexed to the 2D lattice. The phenylene units carrying the three hexadecyloxy side chains appear as bright dots (cf. white circles in Fig. 2a and p), and prove the overall threefold symmetry of the backbone. Its side chains are oriented along two different HOPG zig–zag directions (with an angle of 120° relative to each other). Hexagonal nanopores (arrows 1, 3, and 5 in Fig. 2a, e, and i) of diameter D1 = 2.8 nm are surrounded by six intermolecular triples of hexadecyloxy side chains of ABB ordering (cf.Fig. 2i) with the innermost chains aligned in clockwise (−) (Fig. 2a, e, i and p) or counterclockwise (+) orientation (see ESI‡).23 While some contrast variation is observed in these pores (Fig. 2a), we cannot assign a specific arrangement of intercalated OA molecules. We relate this to rapid dynamics on the timescale of STM and/or a mismatch of the solvent cluster size relative to the dimensions of the nanopore (cf. ESI‡). Smaller intermolecular nanopores (arrows 2, 4, and 6) between pairs of molecules are surrounded by two side-chain triples of ABB ordering with the innermost chains aligned in either clockwise (−) or counterclockwise (+) orientation, and two MSW backbone segments, containing the methyl carboxylate units in a distance of D2 = 1.5 nm (Fig. 2p).24 In these pores, four intercalated hydrogen-bonded OA dimers are clearly resolved by STM (arrows 7 to 10 in Fig. 2p, S4c and d in the ESI‡). 1b (Fig. 2b) forms a packing identical to 1a despite the presence of three carboxylic acid functions per molecule. As observed for the esters, the carboxylic acids point towards the smaller nanopore and point past each other. They do not form intermolecular hydrogen-bonded (homo)dimers, in contrast to trimesic acid11a,b and its larger homologues.11c,d The functional group distances (D3 = 1.5 nm, Fig. 2q) are unaltered compared to 1a. This similarity demonstrates that the alkyl interdigitation pattern and the interaction of the 1,3,5-isomer with the HOPG substrate is so stable that it overcomes the driving force of carboxylic acids to form hydrogen-bonded dimers. However, OA as a non-inert solvent can coordinate to the carboxylic acid functions and thus prevent dimerization.25 The absence of any OA dimers in the packing of 1b (which is in contrast to 1a, cf. arrows 7 to 10 in Fig. 2p and ESI‡) provides further evidence for an altered solvent intercalation scenario with less dense packing and therefore dynamics faster than the timescale of STM (cf. Fig. S7 in the ESI‡). Moreover, in contrast to the acceptably resolved STM images of 2D crystals of 1a after adsorption from a solution of 1a in 1-phenyloctane (PHO; cf. Fig. S6 in the ESI‡), and the well resolved 2D crystalline packings of 1b using OA (Fig. 2b, q, and S9 in the ESI‡), we observed a highly disordered monolayer for 1b using PHO (cf. Fig. S10 in the ESI‡). More precisely, 1b showed some tendency to form a honeycomb-like assembly at least in some surface regions, however, with irregular distances and alignments. In addition, some contrast variation, probably due to solvent intercalation, is found in the large hexagonal nanopores of 1b using OA as solvent, as it was already discussed for 1a.
2a forms an oblique packing (Fig. 2c, domain sizes of about 40 × 40 nm2, see ESI‡) to which a unit cell of a = (9.3 ± 0.2) nm, b = (5.7 ± 0.2) nm, γ(a,b) = (88 ± 2)° containing two molecules (with backbones rotated by 180°) is indexed. Of each pair of molecules, 14 of the 18 hexadecyloxy side chains (shown in bright and dark blue color in Fig. 2g, l, and m) are oriented along the zig–zag directions of the HOPG substrate, and therefore define the overall packing. As a result of steric constraint, the remaining four side chains of each pair of 2a (shown in purple color in Fig. 2g, l, and m and marked by dashed boxes in Fig. 2c, g, l, and r) cannot adsorb along one of the HOPG zig–zag directions. We interpret the medium bright (Fig. 2c and ESI‡) to bright (Fig. 2r) lines as hexadecyloxy side chains that are aligned along one of the three HOPG armchair directions, indicated as dashed lines in Fig. 2c, g, and r. Being energetically less favorable,13 domains of (cyclo-) alkanes oriented along the armchair direction in a constrained environment have been observed previously,14 although this arrangement is rather rare. In addition, in the STM images medium-bright features similar in shape and length to those of the alkoxy side chains are clearly visible, which cannot be attributed to MSW entities. As in the case of earlier results, OA (i.e. solvent) molecules form dimers on the surface, which fill the otherwise uncovered space on the surface, align in parallel to the (MSW) alkoxy chains, and therefore stabilize the packing. As such, we suspect that the alkoxy chains aligned along the less favorable armchair direction are stabilized by two intercalating solvent dimers, which are clearly visible in the STM images (see ESI‡).
To scrutinize this assumption, for this highly specific question we performed geometry optimizations of four MSWs 2a with four intercalated OA dimers (shown in pale purple and green) on a graphene cutout (C5644H190) containing the relevant monolayer regions (Fig. 3, with 3652 atoms). Our calculations were conducted with a new general force field (called GFN-FF). In brief, this method outperforms other FFs in terms of generality and accuracy, and may reach even more elaborate quantum-mechanical methods in terms of performance. The method therefore appears suitable to verify the robustness of the rather unexpected structural model found in Fig. 2, helping us to clearly distinguish the assignment of STM image contrasts to adsorbate film structures from possible misleading image artefacts. The calculations require only starting coordinates and elemental composition as input, basically treating any chemical structure consisting of elements up to radon automatically, and allows access to systems composed of roughly 9500 atoms arranged in, e.g., highly unsaturated components (such as arylenes and graphite or graphene), which would be prohibitively expensive using even fast semi-empirical quantum mechanical tight-binding calculations (e.g. GFN2-xTB) due to the number of atoms involved. Additional details of the calculations are found in the ESI‡ and ref. 26 and 27. First, the results (Fig. 3) agree perfectly with the models (cf.Fig. 2g, l, and m) fitted to the data obtained in STM imaging (Fig. 2c and r).28 More precisely, after geometry optimization of the aforementioned starting geometries by GFN-FF, the alkoxy side chains in Fig. 3 (bright and dark blue, purple) and the OA dimers (pale purple, green) retain their general orientations, even if they are orientated along the unfavorable HOPG armchair directions.13a,14 Moreover, hydrogen bond dimerization of the OA carboxylic acid groups is also retained during optimization. Both these results resemble the high degree of descriptiveness of non-covalent interactions by the GFN-FF method. A more detailed discussion of this issue is given in the ESI.‡
As for 1a, the ester groups in 2a do not interact, and are D4 = 2.3 nm apart (Fig. 2r). However, while both ester 1a and acid 1b MSWs form the same pattern on HOPG, this does not hold for the ester 2a and acid 2b MSWs. The 1,2,4-substitution pattern in 2b allows a packing of the molecules such that the carboxylic acids drive the formation of 1D supramolecular chains via hydrogen bonds of the acids in the 1- and 4-positions – despite using OA as a solvent. To this 1D packing, a unit vector of a = (4.4 ± 0.2) nm is indexed, however, neglecting different backbone orientations. According to this packing distance along a, we assume that only 10 of the 12 side chains along arrow 12 (Fig. 2h) or arrow 14 (Fig. 2n) can adsorb along the HOPG main axis directions (assuming the parallel conformer),13a while the residual chains point towards the solution phase. In terms of their packing distances, the alkoxy chains have to adapt to the distances given by backbone sizes and cyclic carboxylic acid dimer packing distances. In addition, only 2 of 3 chains per molecule can adopt the conformation marked by arrow 11 (Fig. 2h) or arrow 13 (Fig. 2n), while the third chain may share the remaining space with intercalated solvent molecules.29 Regarding the 1D packing, 2b can be viewed as a larger homologue of terephthalic acid.11b The carboxylic acids in the 2-position in turn do not interact with other carboxylic acids and are spaced apart D5 = 2.4 nm in the nearest possible orientations. The two possible orientations of the molecules (cA and c′A in Fig. 2) within a chain occur randomly, introducing some disorder. The reduced STM image resolution of 1b (compared to the images of 1a, 1b, and 2a) is attributed to the reduced order owing to the aforementioned factors. However, each molecule determines the orientation of the other molecule (cB or c′B, respectively) by side-chain interactions (arrows 11 and 13), whereas the orientation of the molecules in the subsequent row is not influenced (arrows 12 and 14). Occasionally, hydrogen-bonded trimers are observed (Fig. 2s–u), which result from the interaction of carboxylic acids in the 1- and 2-positions, which is sterically unfavorable for phthalic acid.11b In addition, trapezoid superstructures containing 11 and 12, as well as both enantiomers of a rhomboid superstructure containing 10 molecules were observed under alike imaging conditions (see Fig. S17 in ESI‡). These moieties include carboxylic acid dimers at the 1-/4-positions as well as at the 1-/2-positions. The superstructures formed are interpreted as a consequence of two energetically similar packaging motives as a result of the increased size as compared to benzene dicarboxylic acids.11b
Footnotes |
† Dedicated to Fred Wudl on the occasion of his 80th birthday. |
‡ Electronic supplementary information (ESI) available: Description of the synthesis, aggregation and the photophysics as well as computational and STM details. See DOI: 10.1039/d1sc01381e |
§ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2021 |