Hongbing
Wang
ab,
Jinping
Hu
ab,
Zhaofeng
Liang
*c,
Huan
Zhang
ab,
Chaoqin
Huang
ab,
Lei
Xie
c,
Zheng
Jiang
abc,
Han
Huang
d and
Fei
Song
*abc
aShanghai Institute of Applied Physics, Zhangheng Road 239, Shanghai, 201000, China. E-mail: songfei@sinap.ac.cn
bUniversity of Chinese Academy of Sciences, Beijing, 100000, China
cShanghai Advanced Research Institute, Zhangheng Road 239, Shanghai, 201004, China. E-mail: liangzf@sari.ac.cn
dSchool of Physics and Electronics, Central South University, Changsha, 410073, China
First published on 26th January 2023
On-surface Ullmann coupling has been considered an appealing approach for the precise fabrication of carbon-based covalent nanostructures under solution-free conditions. However, chirality has seldom been discussed in Ullmann reactions. In this report, self-assembled two-dimensional chiral networks are initially constructed in a large area on Au(111) and Ag(111) after adsorption of the prochiral precursor, 6,12-dibromochrysene (DBCh). Self-assembled phases are then transformed into organometallic (OM) oligomers after debromination, preserving the chirality; in particular, the formation of scarcely reported OM species on Au(111) is discovered herein. With the aryl–aryl bonding induced after intensive annealing, covalent chains are fabricated via the cyclodehydrogenation between chrysene blocks, resulting in the formation of 8-armchair graphene nanoribbons with staggered valleys on both sides. Before chiral polymer chains are constructed by chrysene blocks, the high structural flexibility of OM intermediates on Ag(111) is also revealed during reactions, which is derived from the twofold coordination of Ag atoms and conformationally flexible metal–carbon bonding. Our report not only provides solid evidence of atomically precise fabrication of covalent nanostructures with a feasible bottom-up approach but also sheds insights into the comprehensive investigation of chirality variation from monomers to artificial architectures via surface coupling reactions.
With efforts devoted to the intrinsic prochirality of precursors, the structural transformation from prochiral reactants to advanced chiral polymers has been addressed through a sequence of single-molecule reactions,23–26 shedding new insights into on-surface chemistry.27,28 It is thus interesting to explore the crucial role of molecular prochirality and adsorbate chirality to on-surface Ullmann couplings with single-molecule-resolution imaging by monitoring variations of adsorbate chirality across polymerization reactions. On the other hand, it has been commonly recognized that metal surfaces possessing varying catalytic activities have significant influence on both the self-assembly and on-surface reactions.29–31 In this context, the chiral adsorption of precursors with rich enantiomers on different metal substrates is quite meaningful to unravel the variation of adsorption chirality across Ullmann couplings.
Within this context, 6,12-dibromochrysene (DBCh, Scheme 1) with the planar chirality, has been selected as the precursor to investigate the molecular chiral adsorption and chirality variation during surface Ullmann coupling on different metal substrates, namely Au(111) and Ag(111). Although a previous report has demonstrated the chiral adsorption of DBCh and the chiral arrangement of polymers on Au(111) and Cu(111),32 the process of chirality variation and preservation during reactions is still unclear due to the limited resolution of images. Recently, Cai et al. reported the formation of a chiral domain on Au(111) and one-dimensional metal–ligand chiral chains on Cu(111) with high-resolution scanning tunnelling microscopy (STM) and non-contact atomic force microscopy,33,34 and investigated the related electronic structures. However, the issue of chirality variation and structural variability of organometallic (OM) intermediates and polymers remains vague. As known, the activity of silver is between Au and Cu; Ag(111) is thus utilized, in addition to Au(111), to thoroughly explore the chiral molecular adsorbate and chirality variations during Ullmann coupling in this report. First, the chiral adsorption of DBCh with a hexagonal lattice is revealed on both Au(111) and Ag(111). Afterwards, transformation into chiral OM intermediates is revealed, including OM dimers on Au(111) and oligomers with flexible structures on Ag(111), which have not been properly addressed before and shall be attributed to the chiral adsorption of DBCh on the Au substrate and the flexible coordination of Ag, respectively. In the end, chiral nanoribbons (GNR) with alternating valleys on opposite sites are fabricated after overcoming the steric hindrance. Inspired by the high-resolution STM characterization and density functional theory (DFT) calculations, intermediate states with varied chirality during the reactions are clearly identified towards the fundamental understanding of chirality variation across Ullmann coupling.
As already mentioned, the ACW adsorption of DBCh enantiomers is rigid at the elbow site of reconstruction lines, as clearly confirmed in Fig. S1a in the ESI.† In contrast, the CW phase in Fig. 1d gets extended with flexibility at fcc sites of the herringbone reconstruction (see also in Fig. S1b in ESI†). Meanwhile, it can also be identified in Fig. 1d and e that the homochiral CW and ACW domains are constructed in order by intact DBCh building blocks with either the S or R adsorption chirality and are supposed to be stabilized mostly by the halogen bonding motifs in a triangular form as marked with dashed lines (the halogen bonding length is revealed to be about 0.33 ± 0.02 nm from statistics). It can thus be concluded that the as-formed nanostructure with the X3-synthon geometry36 is employed as the basic unit towards the fabrication of an extended porous network. Correspondingly, the tentative structural model is proposed in Fig. 1f for the CW domain (the model for the ACW phase is also shown in Fig. S1c in ESI†), highlighting the bonding motifs, which mediate the homochiral domains.
In fact, chains with limited length are also formed in connection to porous structures at the sub-monolayer coverage, as observed in Fig. 1b and c. Detailed information is shown in Fig. S1d and e, respectively, in ESI,† where precursors are closely packed with the intermolecular hydrogen and halogen bonding motifs. More intriguingly, R and S adsorbed enantiomers can also meet each other in the same chain, as presented in Fig. S1d.† Subsequently, once the coverage is increased to about 1 ML as presented in Fig. 2a, a two-dimensional (2D) porous network is ultimately established by expanding homochiral hexagonal domains, and regular nanopores are resolved within the network. Since the herringbone reconstruction can still be recognized at the full coverage, indicating the weak interaction between chiral adsorbates and substrate, the 2D hexagonal network is thus mediated by the delicate balance between the molecule–substrate (vdW force) and intermolecular (the halogen and hydrogen bonding) interactions. As inspired by the high-resolution STM in Fig. 2b, six DBCh molecules are tessellated into a hexagonal structure on the basis of the trimeric synthon, with a pore-to-pore distance of about 2.08 ± 0.02 nm (highlighted by the white rhomb). To verify the stability of homochiral domains, as an example, structural optimization of the CW arrangement has been performed with DFT calculations. As illustrated in Fig. 2c of the simulated STM hybridized with the model structure, the predicted adsorption configuration is, in general, consistent with the experimental observation, including topography (density of states) and lattice parameters. For instance, the periodicity of the CW domain is revealed to be 2.10 nm while the experimental value is 2.08 nm, and Br terminals are found to be interacting with neighbours in the form of X3 synthon geometry with the halogen bonding distance of 0.36 nm (0.33 nm from experiments). Nevertheless, the ACW homochiral phase is also discovered at the full monolayer coverage, which is presented in Fig. S1f in ESI.† In practice, the occurrence of CW and ACW phase is roughly equal from statistics, although they are locally separated from each other.
The new structure might be assigned to OM oligomers consisting of debrominated molecules with the same adsorption chirality (S) and Au adatoms in between,34 and a tentative model is overlaid on top of STM in Fig. 3a (yellow circle: Au atoms), which is thus expected to be chiral. Furthermore, the OM dimer is also discovered inside the self-assembled structure, as presented in Fig. S2a in ESI,† showing the structural transition from self-assembly to OM species. These observations suggest that OM oligomers with chiral adsorption can indeed be induced on Au(111), while OM intermediate states have been seldom reported in previous reports, suggesting the critical role of molecular chiral adsorbate towards the formation of OM species on Au(111) by capturing or confining Au adatoms in between.
Sequential annealing to 473 K induces the formation of linear oligomer chains with a large area, and alternative protrusions are identified on both sides of chains, as shown in Fig. 3b. The high-resolution STM in Fig. 3c undoubtedly shows that the protrusion actually consists of two parts with the opposite appearance: one part is imaged bright while the other one is relatively dim. Thus, one can conclude that while the covalent C–C bonding is formed between completely dehalogenated monomers, oligomer chains are not planar and apparently bear the chirality. As indicated by the structural model of covalent oligomers (both the side view and top view) overlaid atop in Fig. 3c, the steric hindrance between hydrogen atoms of neighbouring chrysene blocks results in the opposite tilting of consecutive units with respect to the substrate,34 which is consistent with the alternative bright and dim features in STM topography. Similar structures could also be found on the surface and are presented in Fig. S2b in ESI.† Subsequently, annealing to 573 K induces apparent changes to oligomer chains, as discovered in Fig. 3d. For instance, the previous feature with alternative brightness is missing, which might be assigned to the cyclodehydrogenation between neighbouring blocks by overcoming the steric hindrance and the corresponding formation of planar chains. The zoom-in view inset shows that the polymer chain is constructed by the alternating arrangement of S and R chrysene units with the armchair edge. Further annealing to 623 K makes the polymer chain clearly recognized with alternating valleys on the opposite site, namely, the N = 8 armchair graphene nanoribbon (8AGNR), as seen in Fig. 3e and S2c in ESI.† The DFT-optimized structure of the 1D chain is overlaid on top of STM in Fig. 3e, giving reasonable agreement with the experimental data (the distance between two neighbouring valleys on the same side is found to be 1.38 nm and 1.35 ± 0.02 nm from STM measurements). Further annealing to 673 K, no visible change is induced, as shown in Fig. S2d.† Interestingly, with intensive annealing to 723 K in Fig. 3f, the covalent chain is still robustly adsorbed on Au(111) with the inner structure clearly resolved. Consequently, it is reasonable to conclude that the self-assembly adsorbate homochirality on Au(111) is removed after Ullmann coupling, and OM chains with chirality can be induced as intermediate states.
In addition to the closely packed phase, DBCh enantiomers can also be anchored on Ag(111) with other assembled patterns. As displayed in Fig. 5a, homochiral porous domains, including both CW and ACW arrangements, are also constructed on Ag(111), which is held slightly below RT (275 K) during deposition, and features are analogous to the self-assembly of DBCh on Au(111). More intriguingly, CW and ACW tessellations are found on the same terrace with the domain boundary clearly resolved. High-resolution STM images of CW and ACW phases are illustrated in Fig. 5b and c, respectively (additional STM images are also provided in Fig. S5a and b in ESI†). Once again, the halogen bonding (bonding length: 0.34 ± 0.02 nm) is proposed between bromine atoms of neighbouring DBCh enantiomers, which is responsible for the formation of trimeric X3-synthons, and homochiral phases are thus constructed from R and S chiral adsorbate enantiomers. As an example, DFT predicted configuration for the CW phase is illustrated in Fig. 5d, which in general provides a satisfactory agreement in both topography and unit cell parameters: the rhombohedral lattice is measured to be 2.09 ± 0.02 nm from STM images, and 2.12 nm from DFT calculations. In fact, these self-assembled structures are quite similar to those on Au(111), for example, the coexistence of both CW and ACW homochiral architectures with a roughly equal ratio.
Intriguingly, the decorated porous structure is also observed connecting to normal pores, as shown in Fig. S4c† and elaborately illustrated in Fig. 6a. Apparently, extra trimeric X3-synthons (coloured dots) are constructed from the chiral adsorption of R or S enantiomers (the measured Br–Br distance: 0.37 ± 0.02 nm), embracing the original hexagonal patterns. STM images with the atomic resolution in Fig. 6b and c provide detailed information on the decorated hexagonal network, and the length of the rhombohedral unit cell is revealed to be 4.17 ± 0.02 nm. Obviously, homochiral hexamers are initially constructed by R and S chiral adsorbed enantiomers and surrounded by the homochiral X3-synthons with both the CW and ACW tessellations, as marked by the black or red circles in Fig. 6b and c. Surprisingly, regardless of its own homochirality, the inner hexamer is always surrounded by 6 homochiral X3-synthons with the regular alternation: 3 synthons are in the CW phase (marked with black dots) while another 3 are in the ACW phase (marked with blue dots). Nevertheless, quite a few exceptions are also found on the surface, where 4 and 5 synthons with the same adsorption chirality are positioned around the inner hexamer, as presented in Fig. S4d and e in ESI.† Notably, dim dots are resolved as well between DBCh precursors, which are assigned again to surface Ag adatoms, suggesting the existence of metal coordination bonding motifs. In contrast, the Ag coordination bonding is not evident in the regular porous phases in Fig. 5. Thus, one can conclude that the hexamers and trimeric clusters are stabilized together by the overall contribution from the halogen bonding, the intermolecular hydrogen bonding motifs as well as the metal coordinated bonds. To rationalize the multiple bonding motifs, DFT calculations have also been performed, and the optimized configuration corresponding to Fig. 6b is presented in Fig. 6d. The simulated STM image well corroborates the experimental findings where the dimension of the unit cell is calculated to be 4.22 nm, in good correspondence to the experimental value (4.17 ± 0.02 nm). It is also proposed by DFT calculations that CW and ACW trimeric synthons are adsorbed on different sites (top and hollow sites, respectively) on Ag(111). While the predominant self-assembled phase on Ag(111) at RT is the closely packed domain, the homochiral porous nanostructures are obtained with the substrate kept slightly below RT during deposition. Thus, it can be inferred that homochiral hexagonal phases with nanopores might be the preferred configuration on both Au(111) and Ag(111).
Notably, further annealing to 423 K results in the formation of substantially ordered OM chains on the surface, as shown in Fig. 7d. Apparently, detached Br atoms are still anchored in between oligomer chains. With the structural transition from tilted oligomers to linear chains, OM intermediates are still constituted by chrysene blocks bonded to Ag adatoms, as clearly visible in the zoom-in STM in Fig. 7e. While R or S chiral chrysene adsorbates are randomly arranged inside oligomers, a part with the only S chiral chrysene, for instance, is resolved as marked with the dashed quadrangle in Fig. 7e. Nevertheless, the local chirality can still be recognized regardless of the random connection of R or S chrysene blocks.
After further annealing to 473 K, OM chains are well preserved with an almost defect-free domain, as revealed in Fig. 7f. Notably, the characteristic of local chirality in OM chains becomes invisible since the majority of oligomer chains is now constructed by R or S chrysene without the specific order. Actually, OM species is revealed to be fairly stable up to 523 K annealing, and detached Br atoms are mediated by the metal coordination of Ag in combination with hydrogen bonding with the neighbouring chrysene blocks.37,38 Again, the C–Ag–C bond is still not linear at this stage with varied bonding angles. Thus, one can infer that the flexible coordination of Ag together with the rich conformations of oligomers, enable the high flexibility of OM chains on Ag(111). From the structural evolution illustrated in Fig. 7, one can also infer that the local chirality of OM oligomers declines with sequential annealing.
On the basis of flexible OM chains, more intensive annealing is applied to explore whether the structural transformation can occur from OM intermediate species to covalent polymers. Surprisingly, structure chaos is observed at the beginning after annealing to above 533 K (Fig. S6a in ESI†), and ordered nanostructures start to emerge only after annealing to 623 K, which prefer to accumulate along step edges as shown in Fig. 8a. Apparently, the total amount of species remaining on the surface gets significantly reduced due to desorption after intensive annealing. Moreover, a considerable number of bright dots are closely packed on both sides of polymers, which are attributed again to detached Br atoms as illustrated in detail in the zoom-in STM view in Fig. 8b. Meanwhile, such observation also agrees well with previous reports,31,37,38 where Br atoms falling off from precursors are stably adsorbed on the metal substrate.
Furthermore, several unique nanostructures are induced after the high-temperature annealing, as presented in Fig. 8b and c. First, the cyclodehydrogenation is induced after kicking off the coordinated Ag atom, and the C–C covalent dimer is thus fabricated by covalently connecting chrysene units,32,33 as highlighted with the molecular model in Fig. 8b. Besides, covalent dimers can also be linked together into a double-dimer architecture induced by the aryl–aryl bonding after further dehydrogenation, and the tentative model is also proposed atop. One should also bear in mind that the rotation around the single C–C bond between chrysene units shall make dimers tilted with respect to the surface and thus the partial brightness in STM. Besides, the enclosed-ring architecture is discovered as well in Fig. 8c. Judged from its appearance, such enclosed architecture is probably induced by the aryl–aryl connection between a homochiral dimer and a chrysene unit after the selective dehydrogenation, and the possible structural model is proposed atop the STM image. In fact, both the double-dimer structure and the enclosed ring architecture are induced by the selective dehydrogenation of prochiral chrysene blocks at the meta position, followed by the aryl–aryl coupling. Other similar nanostructures can also be found and are summarized in Fig. S6b and c in ESI.† Further annealing to 673 K in Fig. 8d makes the Ag(111) surface tidied up, and only the covalent polymer chains with staggered defects are formed via the alternative linkage of R and S chrysene units, which have the feature of 8AGNR but bear valley defects on the opposite site.32 However, it is also notable that the total amount of 8AGNR is small with the limited length, which is probably caused by the substantial desorption after the intensive annealing.
In a word, evidenced by the high-resolution STM, chiral tessellations of DBCh enantiomers are induced on both Au(111) and Ag(111) into hexagonal pores with self-assembly, and OM oligomers with local chirality are discovered on Au(111) and Ag(111) (bearing flexible structures) during Ullmann coupling. While the unique homochiral arrangement of self-assembled DBCh is induced on Au(111), multiple homochiral structures are fabricated on Ag(111), implying the significant templating effect from the substrate. For example, the coordination motifs via surface adatoms are fairly apparent on Ag(111), but instead, it is occasionally observed on Au(111) due to the shortage of surface adatoms. Notably, no debromination of DBCh is observed on Ag(111) at RT, while partial debromination of other halogenated precursors has been frequently observed on Ag(111) at RT.14,37 Such difference might be assigned to the relatively high activation energy of the C–Br bonding in DBCh, which is higher than the activation barrier of the C–Br bonding in Br2Py37 by about 0.4 eV as predicted by DFT calculations. In addition, DFT calculations have also been employed in understanding the intermolecular bonding motifs by optimizing adsorption configurations and providing reasonable structures to the experimental data. Importantly, the key role of intermolecular halogen bonding and metal coordinated interactions is witnessed not only in the chiral adsorption of DBCh enantiomers, but also in mediating the stability of partially chiral oligomers constructed from chrysene blocks, especially flexible polymer chains on Ag(111). While OM chains with local chirality are identified on Au(111) and Ag(111) at certain areas, the ultimate formation of 8AGNR completely removes the pristine adsorption chirality by alternatively linking S and R chrysene units, as such a configuration is the most stable structure from the energy point of view.32 it is also obvious that the quality of graphene nanoribbons is substantially influenced by metal substrates. Based on these observations, the gradual fading of adsorption chirality from the self-assembly to covalent polymers is identified across on-surface Ullmann coupling.
The Au(111) and Ag(111) (MaTecK, Germany) surfaces were cleaned by repeated Ar+ ion sputtering (alternating energy at 1.2 keV and 0.8 keV) and annealing cycles (623 K for Ag(111) and 673 K for Au(111), respectively). The DBCh precursor (Sigma-Aldrich, purity > 97%) was thoroughly degassed prior to deposition, and thermally evaporated onto the metal surface with a rate of 0.1 monolayer (ML) per min from the commercial evaporator (Omnivac, Germany). Herein, 1 ML is referred to the full coverage of DBCh molecules on the surface as judged from STM, and the evaporation rate was calibrated using a quartz crystal microbalance. Both the Au(111) and Ag(111) substrates were kept at RT (300 K) during deposition unless otherwise stated and were subsequently heated in sequence to initiate Ullmann coupling reactions.
DFT calculations were done by utilizing the Vienna ab initio simulation package (VASP).40,41 The PBE functional was adopted in combination with the third-generation van der Waals dispersion correction from Grimme (DFT-D3) and the projector-augmented wave (PAW) method for atomic cores.42–44 A plane-wave cut-off energy of 450 eV was employed, while the k-point sampling was chosen to be 3 × 3 × 1 to obtain a realistic and accurate picture of the energetics. Lattice parameters of Ag(111) and Au(111) after optimization were found to be 2.89 Å and 2.88 Å, respectively, in good agreement with previous reports.35,45 Optimized adsorption geometries for DBCh on the surface were obtained on the three-layer substrates based on the energy minimum, where the first layer and precursors are free to relax with the bottom two layers fixed until the forces on the active atoms dropped below 0.02 eV Å−1. The convergence criterion is 10−4 eV in total energy difference. STM simulations were performed in the framework of Tersoff–Hamann approximation with the gnuplot software by plotting partial charge densities, which were acquired by a static self-consistent calculation based on the DFT optimized configurations. In the end, the selection of a proper structure model is guided by the lattice parameters and STM features.
Footnote |
† Electronic supplementary information (ESI) available: Additional STM images of the self-assembly and the formation of various oligomers during Ullmann coupling. See DOI: https://doi.org/10.1039/d2na00789d |
This journal is © The Royal Society of Chemistry 2023 |