Matsuhiro Maedaa,
Kotoka Odaa,
Ichiro Hisakib and
Kazukuni Tahara*a
aDepartment of Applied Chemistry, School of Science and Technology, Meiji University, 1-1-1 Higashimita, Tama-ku, Kawasaki, 214-8571, Japan. E-mail: tahara@meiji.ac.jp
bDivision of Chemistry, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
First published on 10th October 2023
The effect of the core size on the structure and chirality of self-assembled molecular networks was investigated using two aromatic carboxylic acid derivatives with frameworks displaying C3h symmetry, triphenylene derivative H3TTCA and dehydrobenzo[12]annulene (DBA) derivative DBACOOH, each having three carboxy groups per molecule. Scanning tunneling microscopy observations at the 1-heptanoic acid/graphite interface revealed H3TTCA exclusively forming a chiral honeycomb structure, and DBACOOH forming three structures (type I, II, and III structures) depending on its concentration and whether the system is subjected to annealing treatment. Hydrogen bonding interaction patterns and chirality were carefully analyzed based on a modeling study using molecular mechanics simulations. Moreover, DBACOOH forms chiral honeycomb structures through the co-adsorption of guest molecules. Structural diversity observed for DBACOOH is attributed to its relatively large core size, with this feature modulating the balance between molecule–molecule and molecule–substrate interactions.
Dimeric pairs of aromatic carboxylic acids formed each by complementary hydrogen bonding interactions between the carboxy groups have been extensively used in the design of molecular packing structures in three-dimensional (3D) space because of their directionality and selectivity (Fig. 1(a)).15–17 This hydrogen bonding interaction can also be used in two-dimensional space to construct SAMNs on surfaces both in UHV and interfacial conditions.18 Many such investigations have been performed in particular under ambient interfacial conditions.18–20 Trigonal aromatic molecules having each three carboxy groups are ideal molecular building blocks for constructing porous SAMNs displaying hexagonal periodicity.21 Lackinger and coworkers reported systematic studies on structural control of SAMNs at the liquid/solid interfaces using trimesic acid (TMA) and 1,3,5-benzenetribenzoic acid (BTB), which are aromatic carboxylic acid derivatives having molecular frameworks displaying D3h symmetry (and have substitutions at the 1,3,5-positions of the central benzene rings, Fig. 1(b)).22–25 TMA is the smallest such molecule that forms honeycomb (chicken-wire) and flower structures at the liquid/graphite interfaces, depending on its concentration and the identity of the solvent.20,23 In the honeycomb structure, the TMA molecules are connected by complementary-type (R22(8)) hydrogen bonding interactions between the two carboxy groups, whereas in the flower structure, the TMA molecules are connected by both complementary-type (R22(8)) and cyclic-type (R33(12)) hydrogen bonding interactions, with the cyclic-type interaction formed by three carboxy groups (Fig. 1(a)). The size of the aromatic core has impacts on SAMN formation.20 BTB is a larger molecule having a framework displaying D3h symmetry, and in which the 1,3,5-position hydrogens of the central benzene ring are substituted with 4-carboxyphenyl groups (Fig. 1(b)); the molecule forms a porous honeycomb as well as oblique structures, and a non-porous close-packed structure at the liquid/graphite interfaces.20,25 In the oblique structure, a hydrogen bonding interaction between the carboxy group and a hydrogen atom of the benzene ring would be formed in addition to the complementary-type (R22(8)) hydrogen bonding interactions between the carboxy groups (Fig. 1(a)). Moreover, large molecules based on a TMA motif also produce various porous SAMNs.26–32 In some of these SAMNs, guest molecules fit into the pores are co-adsorbed.27,29 Moreover, a non-porous SAMN sometimes shows a structural rearrangement in response to the addition of guest molecules through co-adsorption.33 Recent work has focused on the adaptability of these SAMNs to external stimuli such as heat,34 to an external electric field,35–37 and to a combination of these stimuli and additives.33,38–41 Nevertheless, despite the many studies on TMA and its expanded motifs displaying D3h symmetry, only a few studies have been reported on SAMNs having molecular frameworks displaying C3h symmetry.42–44 The SAMN formation and guest co-adsorption of the simplest molecule having a molecular framework displaying C3h symmetry, namely triphenylene-2,6,10-tricarboxylic acid (H3TTCA), have been reported (Fig. 1(c)).43,44 However, to the best of our knowledge, there has to date never been reported a study on the effect of core size on the formation of SAMNs using a planar molecule having the molecular framework displaying C3h symmetry. Understanding the effect of symmetry group reduction and size modulation of the molecular building blocks on the formation of SAMNs is important for the establishment of 2D crystal engineering (Fig. 1(a) and (c)).
Adsorption of organic molecules onto a surface often generates single-point molecular chirality even for achiral molecules in 3D space because of symmetry group reduction.45–48 The aromatic carboxylic acid derivatives having molecular frameworks displaying D3h symmetry would remain achiral on the surface at the single-molecule level (Fig. 2(a) and (c)). Here the orientation of the carboxy groups on a surface is ignored because it does not influence complementary hydrogen bonding interactions between the carboxy groups (R22(8) and R33(12)) (Fig. 1(a) and 2). While the molecules displaying D3h symmetry are achiral at the single-molecule level, their SAMNs sometimes adopt one of the chiral plane groups (P1, P2, P3, P4, and P6).49 For instance, the domains of the honeycomb structures of TMA and BTB are chiral, and this supramolecular chirality is distinguished only by the domain orientations with respect to the graphite lattice underneath (Fig. 2(e)). Specifically, the supramolecular chirality is determined by the sign (positive or negative) of α angle between the unit cell vector of the honeycomb structure and the nearest normal to the main symmetry axes of graphite (Fig. 2(e)). On the other hand, the achiral aromatic carboxylic acid derivatives having molecular frameworks displaying C3h symmetry (Fig. 1(c) and 2(b)) show single-point molecular chirality on the surface. We defined this single-point chirality as an R-type or L-type (Fig. 2(d)). The chiral R- or L-type molecule appears with equal probability on the surface. Assembling such a racemic mixture on a surface can lead to two outcomes: the formation of separate domains of homochiral structures (2D-conglomerate) or the formation of an achiral structure composed of equal amounts of both chiral molecules in the domain (2D-racemate).50,51 The supramolecular chirality of the domains in 2D-conglomerates could be distinguished by the orientation of each molecule (R- or L-type) as well as the relative orientation of the domains with respect to the substrate lattice (Fig. 2(f)). These SAMNs display one of the chiral 2D plane groups. In the case of the 2D-racemates, the SAMNs can display the chiral as well as achiral plane groups.50,51 When the 2D-racemate displays one of the chiral plane groups, the domain handedness is defined by its relative orientation with respect to the substrate lattice. Even for the simplest case of a molecule having the molecular framework displaying C3h symmetry, a detailed analysis of the chirality of the SAMNs has, to the best of our knowledge, not yet been reported (Fig. 1(c)).42–44 Therefore, addressing supramolecular chirality aspects of these SAMNs is also an intriguing subject.49
In this context, we herein investigated the effect of the core size on the structure and chirality of SAMNs using two aromatic carboxylic acids having frameworks displaying C3h symmetry, namely H3TTCA and dehydrobenzo[12]annulene (DBA) derivative DBACOOH, each having three carboxy groups per molecule, at the 1-heptanoic acid (1-HA)/graphite interfaces. STM observations revealed an exclusive formation of a chiral honeycomb structure by H3TTCA, while DBACOOH shows structural polymorphism, with the structure formed depending on its solute concentration and whether the system is subjected to annealing treatment. DBACOOH was found to form three types of structures, denoted here as type I, II, and III structures. The structural diversity observed for DBACOOH is attributed to its relatively large core size, a feature modulating the balance between molecule–molecule and molecule–substrate interactions. We also investigated guest co-adsorption in the SAMNs of DBACOOH. Chiral honeycomb structures are formed upon the addition of the guest molecule, coronene or phthalocyanine (Fig. 3 and S1†). The guest co-adsorption is essential for the formation of the honeycomb structures by DBACOOH.
The H3TTCA molecules were observed to exclusively form a honeycomb structure at the 1-HA/graphite interface at all solute concentrations (Fig. 4 and S2†). Since π-conjugated cores were observed as bright features in STM images because of their higher tunneling efficiency,54,55 bright triangles composed of three small doughnut-like features were assigned to the π-conjugated cores of H3TTCA. The honeycomb structure consists of the H3TTCA molecules with two different orientations of 60° (Fig. S3†). In the STM images, short rod-like features were observed at the periphery of the hexagonal pore, and these were attributed to co-adsorbed 1-HA molecules. The unit cell parameters of the honeycomb structure were measured to be a = b = 2.3 ± 0.1 nm and γ = 60.1 ± 0.4° (Fig. 4(a)), values consistent with those in the previous reports.43,44 The honeycomb structure adopts the chiral P6 symmetry group. The distance between the centers of the benzene rings of adjacent H3TTCA molecules was measured to be 0.8 ± 0.1 nm. Based on the unit cell parameters and this intermolecular distance, the network models were optimized using molecular mechanics (MM) simulation (COMPASS force field). In the optimized model, shown in Fig. 4(b), the O⋯H distance in the dimeric pair of carboxy groups between adjacent H3TTCA molecules is 0.175 nm,56,57 supporting the formation of complementary-type hydrogen bonding interactions (R22(8)). The honeycomb structure consists of only one enantiomer (R or L-H3TTCA, Fig. 4(d) and (e)), thus indicating the honeycomb structure to be chiral (R or L). This chirality is identical throughout the domains. The two enantiomorphic R and L domains were identified by the negative and positive signs of the value of the α angle between the short unit cell vector a and the one nearest to normal to the main symmetry axes of graphite. The |α| average angle value was determined to be 12 ± 1° (Fig. 4(c) and S4†). Both enantiomorphous domains appear with equal probability, confirming formation of a 2D-conglomerate (Fig. S4†).
Next, we investigated the SAMNs of DBACOOH at the 1-HA/graphite interface. DBACOOH forms three different SAMNs depending on its concentration and whether annealing is carried out (Fig. S5†). At solute concentrations from 1.0 × 10−6 to 1.0 × 10−5 M, type I and II structures were observed (Fig. 5 and 6). After an annealing treatment, at the same concentration range, a type III structure is formed (Fig. 7). At a solute concentration of 1.5 × 10−4 M, and regardless of whether an annealing treatment is applied, DBACOOH forms a non-periodic structure (Fig. S6†) in which the DBACOOH molecules adopting various orientations are packed densely on the surface.
In the type I structure, hexagonal and rectangular pores were observed to be arranged in a periodic manner (Fig. 5(b) and (c)). In these pores, we observed dim rod-like features (Fig. 5(c)), attributed to co-adsorbed 1-HA molecules. The unit cell parameters for the type I structure were measured to be a = b = 5.12 ± 0.06 nm and γ = 59.9 ± 0.4°. And the chiral P6 group was identified to be the plane group of the type I structure, when two DBA molecules located between the cyclic hexamers are not regarded as asymmetry units. In this type I structure, eight structurally distinct DBACOOH molecules were identified in the unit cell: six of them arranged in a hexamer and two located between hexamers. The distance between the centers of nearest benzene rings of adjacent DBACOOH molecules located at the vertices of a hexagonal pore was measured to be 0.70 ± 0.04 nm (red arrow in Fig. 5(c)). And for each of the two non-hexamer DBACOOH molecules, the distance between either DBACOOH molecule and the nearest DBACOOH molecule in the hexamer—specifically the distance between the centers of the nearest benzene rings—happened to be about the same, namely 0.9 ± 0.1 nm, despite their geometries being structurally inequivalent (blue and green arrows in Fig. 5(c)). The difference in the intermolecular distances between that within the hexamer and that between the hexamer and outside the hexamer suggested two different hydrogen bonding interaction modes. Accordingly, a network model was carefully constructed using MM simulation based on the unit cell parameters and these intermolecular distances (Fig. 5(b)). The hexagonal pore consisting, as described above, of the molecules having two different orientations of 60° suggests the formation of a homochiral L or R hexamer (Fig. 2(f)). To assess the chirality of remaining DBACOOH molecules located in between the hexamers, we modeled three possible molecular arrangements: (1) one having the same chirality of the molecules in the hexamer (Fig. S7(a)†), and (2) the other two molecules having the opposite chirality (Fig. S7(b) and (c)†). The network model consisting of the same chiral orientation nicely matches the experimental results. While the orientation of the carboxy groups could not be determined from the STM images, we considered that, to maximize the number of hydrogen bonding interactions, the DBA molecules in the hexamer adopt a Cs symmetry group with a different orientation for one of the three carboxy groups. The MM simulation yields two mean intermolecular O⋯H distances, namely 0.172 and 0.188 nm, for the two different types (R22(8) and D types) of hydrogen bonding interactions between the carboxy groups (Fig. 5(d) and (e)).56,57 Moreover, in the D type interaction, there is a weak hydrogen bonding interaction between a hydrogen atom of the benzene ring and an oxygen atom of the carboxy group (0.249 nm) (Fig. 5(d) and (e)).56,58 Further details on the network modeling are described in ESI (Fig. S7†). The type I structure was observed to be chiral throughout the domains, and this chirality could be differentiated by the sign (positive or negative) of the α angle value (R or L). An |α| average angle value of 27.5 ± 0.6° was measured (Fig. S8†). Both chiral domains appear with equal probability and thus the type I structure is regarded as a 2D-conglomerate.
The type II structure was observed to contain a narrow rectangular pore, and to be composed of four independent DBACOOH molecules with different orientations. Our results indicated unit cell parameters of a = 2.8 ± 0.1 nm, b = 4.9 ± 0.1 nm, and γ = 54 ± 1° and a chiral P2 plane group for the type II structure. The distances between the benzene rings of the adjacent independent DBACOOH molecules were measured to be 0.8 ± 0.1 nm, 0.9 ± 0.1 nm, and 0.7 ± 0.1 nm, respectively (red, green, and blue arrows in Fig. 6(a)). A network model optimized by carrying out an MM simulation based on the unit cell parameters and these intermolecular distances are shown in Fig. 6(b).
Of the three possible network structures with a different combination of homo- or heterochiral molecular orientations—namely the homochiral (R,R,R,R or L,L,L,L) structure and heterochiral (R,R,L,L or L,L,R,R) and (R,R,R,L or L,L,L,R) structures (Fig. S9†)—the heterochiral (R,R,L,L or L,L,R,R) structure was found to be the one in good agreement with the π-core positions in the STM images (Fig. 6(d) and (e)). We assumed in the MM model the formation of a maximum number of intermolecular hydrogen bonding interactions between the carboxy groups (see ESI†). The four inequivalent DBACOOH molecules (R-C3h, R-Cs, L-C3h and L-Cs symmetry) are connected through the two types (R22(8) and D types, Fig. 6(e)) of hydrogen bonding interactions with the two mean O⋯H distances of 0.172 nm and 0.179 nm.56,57 Despite the type II structure containing equal numbers of the chiral R and L molecules, domains having two different orientations with respect to the graphite lattice underneath were identified. This handedness was again determined by the sign of the α angle (+ or − for R or L). The |α| average angle value was measured to be 0.2 ± 0.3° (Fig. 6 and S10†). Both enantiomorphous domains appear with equal probability.
The type III structure appears exclusively after annealing treatment using the sample solutions at concentrations from 3.0 × 10−6 to 1.0 × 10−5 M (Fig. 7(a)). The structural change indicates the type III structure to be a thermodynamically favored phase relative to the type I and II structures. In the type III structure, the DBACOOH molecules were observed to be densely arranged in a linear fashion. Our results showed unit cell parameters of a = 1.5 ± 0.1 nm, b = 4.5 ± 0.1 nm, and γ = 70 ± 1° and the chiral P2 symmetry group for this type III structure. The distances between the benzene rings of the adjacent independent DBACOOH molecules were measured to be 0.7 ± 0.1 nm, 0.9 ± 0.1 nm, and 0.8 ± 0.1 nm, respectively (red, green, and blue arrows in Fig. 7(a)). Two possible structures, a homochiral alignment consisting of the molecules with the R- or L-orientation and a heterochiral alignment consisting of the molecules with both R- and L-orientations, were modeled (Fig. S11†). The homochiral alignment optimized by performing an MM simulation was found to nicely match the positions of each triangular π-core of DBACOOH in the STM image (Fig. 7(b)). In the MM model, there is one type of hydrogen bonding interaction (D type) between the carboxy groups (Fig. 7(d) and (e)), with a mean O⋯H distance of 0.173 nm.56,57 The supramolecular chirality of the type III structure was also determined by the sign (positive or negative) of the α angle (R or L, Fig. 7(c)), with the|α| average angle value measured to be 13 ± 1° (Fig. S12†). As both chiral domains appear with equal probability, the type III structure is concluded to be a 2D-conglomerate.
Structural rearrangement of H3TTCA SAMNs upon addition of guest molecules such as phthalocyanine and coronene has been reported.43,44 In this context, we investigated the effect of these guest molecules on the SAMN formation of DBACOOH.
First phthalocyanine, with its large square π-core, was tested (Fig. 3(a) and S1†). We tested various pre-mixed solutions of DBACOOH and phthalocyanine with different guest/host (G/H) molar ratios but in which the concentration of DBACOOH was kept constant at 1.0 × 10−5 M. Extended domains of the honeycomb structure were observed at exclusively a G/H ratio of 5/1 (Fig. 8(a) and S13†). Inspection of the STM image of this SAMN showed fuzzy bright circular features at the hexagonal pores. We attributed these bright features to co-adsorbed phthalocyanine molecules (Fig. 8(c) and (d)). The fuzzy appearance is attributed to the lateral mobility of the guest molecules at the pore. Also, distorted hexagonal pores were observed at the domain borders, and a square bright feature was occasionally observed in these distorted pores (Fig. 8(c) and (d)). This square appearance is attributed to size complementarity between the guest molecule and distorted pore. Unit cell parameters for this honeycomb structure (P6 group) were measured to be a = b = 2.91 ± 0.07 nm and γ = 60.1 ± 0.4°. A corresponding molecular model optimized by performing an MM simulation based on experimental data is shown in Fig. 8(b). This optimization shows the DBACOOH molecules adopting two orientations with a difference of 60°, suggesting the formation of a homochiral hexamer. The complementary hydrogen bonding interactions of the R22(8) type between the carboxy groups sustain the cyclic hexamer (Fig. 8(b)). An O⋯H distance of 0.173 nm is in the MM model. The formation of L versus R supramolecular chirality was discerned by noting the sign (positive or negative) of the α angle. In the current work on DBACOOH with phthalocyanine, the |α| average angle value was measured to be 17 ± 1° (Fig. S14†). This result is indicative of the honeycomb structure being a 2D-conglomerate. The two-component honeycomb structure was transformed into a type III structure over the activation barrier for the structural rearrangement upon annealing (Fig. S13†). This result implied the two-component honeycomb structure to be a kinetically trapped phase.
Next, coronene with its D6h symmetry was tested as a guest molecule (Fig. 3(b) and S1†). Note the symmetry of this molecule fitting that of the pore, but its size being smaller than that of the pore. Extended domains of the honeycomb structure were observed at exclusively a G/H ratio of 10/1 (Fig. 9(a) and S15†). Fuzzy features with different levels of brightness were observed at the pore, indicating the co-adsorption of coronene molecules (Fig. 9(c) and (d)). The diameter of the brightest feature in the pore was measured to be 1.9 ± 0.1 nm, a value corresponding to twice the coronene diameter. Therefore, these features were attributed to two or three co-adsorbed coronene molecules. Lackinger and coworkers attributed the different levels of brightness of the co-adsorbed coronene molecules in the honeycomb structures of TMA and BTB to different numbers of co-adsorbed coronene molecules59 and/or to lateral mobility of the guest molecules in the pore.60 The coronene molecules with lateral dynamics in the pore would reflect their appearance in the STM images. Moreover, fuzzy oval features were observed in the distorted hexagonal pores at the domain boundaries (Fig. 9(c) and (d)). Lengths of 1.8 ± 0.1 nm and 1.1 ± 0.1 nm were measured for, respectively, the long and short sides of the fuzzy oval features, indicating these features to be two co-adsorbed coronene molecules. Our measurements indicated unit cell parameters of a = b = 2.91 ± 0.08 nm and γ = 60.1 ± 0.6° for this honeycomb structure. An MM-optimized network model based on the above result is shown in Fig. 9(b). Moreover, the two-component honeycomb structure was also found to be chiral, and its handedness to be determined by the sign (positive or negative) of the α angle (R or L). The |α| average angle value was measured to be 17 ± 1° (Fig. S16†). This two-component honeycomb structure also transformed into the type III structure through desorption of the guest molecules upon annealing (Fig. S15†). This two-component honeycomb structure was also concluded to be a kinetically trapped phase.
H3TTCA forms only the chiral honeycomb structure, regardless of the solute concentration, and does so as a 2D-conglomerate. In contrast, DBACOOH shows structural polymorphism, with the four structures depending on solute concentration, whether the system is subjected to annealing, and the presence of the guest molecules. The type I and III structures as well as the two-component honeycomb structures are 2D-conglomerates. The type II structure is a 2D-racemate in terms of the composition of the L and R-molecules within the domain, while the domains show two chiral orientations.
We considered the difference in absorbance abilities of H3TTCA and DBACOOH to be a possible explanation for the structural polymorphism observed for DBACOOH but not for H3TTCA. Table 1 summarizes the total energies per unit area (Etotal) estimated by MM calculations. These calculations indicated the Etotal value of the honeycomb structure of H3TTCA to be lower than that of DBACOOH. This result is consistent with the STM observation, which indicates that the H3TTCA molecules form the honeycomb structure while the DBACOOH molecules do not. This is due to a lower molecular density on the surface for DBACOOH. Next, the relative stabilities of the type I, II, and III structures of DBACOOH are compared. This comparison shows a small difference in the Etotal values between the type I and II structures. And on the other hand, the calculations showed a highly negative Etotal value for the type III structure. This difference between the type III structure and the other structures is attributed to the type III structure having the highest surface molecular density. Again, the highly negative Etotal value for the type III agrees with the experimental observations. The calculations showed roughly comparable contributions of Emol–sub and Emol–mol to Etotal for H3TTCA, but much greater contributions of the Emol–sub values than of the Emol–mol values to Etotal for DBACOOH. We considered that the enlargement of the core, upon switching from H3TTCA to DBACOOH, altered the balance of these interactions resulting in the formation of various, yet high-density phases for DBACOOH.
Compound | Structure | Na | Molecular density (nm−2) | Etotalb (kcal mol−1 nm−2) | Emol–molc (kcal mol−1 nm−2) | Emol–subd (kcal mol−1 nm−2) |
---|---|---|---|---|---|---|
a Number of the indicates molecules per unit cell.b Non-bonding interaction energy of the whole system including the substrate (bilayered graphene sheets).c Non-bonding interaction energy obtained from a single-point energy calculation of the molecules without the substrate.d Non-bonding interaction energy between the molecules and substrate, calculated using the equation Emol–sub = Etotal − Emol–mol.e A honeycomb structure of DBACOOH is a hypothetical model construed for a comparison. | ||||||
H3TTCA | Honeycomb | 2 | 0.38 | −30.8 | −12.1 | −18.6 |
DBACOOH | Honeycombe | 2 | 0.27 | −17.5 | −1.5 | −15.9 |
DBACOOH | Type I | 8 | 0.35 | −21.6 | −0.8 | −20.8 |
DBACOOH | Type II | 4 | 0.39 | −23.0 | 0.2 | −23.2 |
DBACOOH | Type III | 4 | 0.56 | −39.9 | −2.6 | −37.3 |
Next, we estimated the interaction energies for the host–guest systems of DBACOOH from the results of MM simulations (Table 2). The Etotal values of the honeycomb structures with one co-adsorbed phthalocyanine molecule per pore and three co-adsorbed coronene molecules per pore were calculated to be lower than those of the type I and II structures, consistent with the experimental observations. The calculations indicated a lower Etotal value for the type III structure than for the honeycomb structure with coronene, results agreeing with the favorable formation of the type III structure upon annealing. The Etotal values for the honeycomb structure with phthalocyanine and the type III structure are essentially identical. The formation of the type III structure upon annealing would relate to the factors that were not considered in the MM simulations, such as solvation and wetting energies.
Structure | Guest molecule | Etotal (kcal mol−1 nm−2) | Ehosta (kcal mol−1 nm−2) | Eguestb (kcal mol−1 nm−2) | Ehost–guestc (kcal mol−1 nm−2) |
---|---|---|---|---|---|
a Ehost values were obtained by performing single-point energy calculations for the host–guest systems in their optimized geometries without guest molecule(s).b Eguest values were obtained by performing single-point energy calculations for the host–guest systems in their optimized geometries without the host networks.c Calculated using the equation Ehost–guest = −(Etotal − Ehost − Eguest). | |||||
Honeycomb | One phthalocyanine | −40.0 | −17.4 | −22.3 | 0.3 |
Honeycomb | Three coronenes | −26.4 | −17.3 | −7.9 | 1.2 |
Footnote |
† Electronic supplementary information (ESI) available: Additional STM images, STM image analysis, details of MM calculations. See DOI: https://doi.org/10.1039/d3ra05762c |
This journal is © The Royal Society of Chemistry 2023 |