Supplementary Information
Impact of heterocirculene molecular symmetry upon two-dimensional crystallization
W.D. Xiao1,2*, Y.Y. Zhang2, L. Tao2, K. Aït-Mansour1,4, K.Y. Chernichenko5,6, V.G. Nenajdenko6, P. Ruffieux1, S.X. Du2, H.-J. Gao2, and R. Fasel1,3*
1 Empa, Swiss Federal Laboratories for Materials Science and Technology, 8600 Dübendorf, Switzerland
2 Institute of Physics, Chinese Academy of Sciences, 100190 Beijing, P.R. China
3 Department of Chemistry and Biochemistry, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
4 ICMP, Ecole Polytechnique Fédérale de Lausanne, Station 3, 1015 Lausanne, Switzerland
5 Department of Chemistry, University of Helsinki, P.O. Box 55, 00014, Finland
6 Department of Chemistry, Moscow State University, Leninskiye Gory, 119992 Moscow, Russia
Contents
- Interaction between heterocirculenesand Au(111)
- Formation of a second layer of 1on Au(111) at a coverage of ~1.5 ML
- Intermolecular interaction between heterocirculenes
*To whom correspondence should be addressed. Email: ; .
1.Interaction between heterocirculenesand Au(111)
The interactions between 1 and Au(111)arestudied using the vdW-DF method.1-8Periodic boundary conditions were appliedto model the adsorption of 1on Au(111).The calculation cells consisted of repeated p(5×5) four-layer Au(111) slabs separated by 18Åof vacuum. A single molecule of1was placed in each slab. The structures were relaxed until residual forces were less than 0.02 eV. Supplementary Figure 1 shows the optimized configuration of 1adsorbed on Au(111).
According to our vdW-DF calculations, the binding energy between 1 and Au(111)is 2.2 eV/molecule.This binding energy is comparable with that of 2.3 eV/molecule between coronene and Au(111) reported by Medeiros et al..6 The vertical distance between 1 and the toplayer of the substrate is 3.5 Å, also comparable to that of 3.3 Åbetween coronene and Au(111)reported by Medeiros et al..6Similar results are obtained for 2 adsorbed on Au(111). Thus, for the adsorption of both 1 and 2 on Au(111) we conclude that the molecule-substrate interaction mainly derives from aweak coupling between the molecular π-electrons and the metal, akin to the one of planar polycyclic aromatic hydrocarbon molecules on Au(111), and in contrast to the strong S-Au bonding widely reported for thiols adsorbed on gold substrates.9-11
Supplementary Figure 1 | Optimized configuration of 1 adsorbedon Au(111).
2.Formation of a second layer of 1 on Au(111) at a coverage of ~1.5 ML
Supplementary Figure 2 | STM images of 1.5 ML of 1 on Au(111). (a) Overview showing the formation of asecond layer on top of a completed first layer of 1. The first and second layers are labelled as I and II, respectively. (b) Close-up of the second layer exhibitingarectangular superlattice with a= 1.05±0.05nm and b= 1.89±0.05nm, respectively. Note that the molecules adopt tilted configurations, similar to thosewithin the (001) plane in the bulk crystal of 1.12,13
3.Intermolecular interaction between heterocirculenes
In plots of the electron localization function (ELF), which is often used to determine formation of bonds, we find no evidence for strong S-S interaction between adjacent molecules. Supplementary Figure 3 shows the iso-surface of ELF = 0.8 of 1 in the centered rectangular lattice. The iso-surface locates only at the molecule. We analyzed the electron density difference, which reflects the charge transfer between adjacent molecules (Supplementary Figure 3(b) and (c)). To get a significant iso-surface, iso-values have to be chosen as small as ±10-6 Å-3. This is clearly also within the regime of vdW interactions. For comparison, this value is around 10-3 Å-3 for FePc molecules adsorbed on Au(111).14
Supplementary Figure 3| Iso-surfaces of ELF and electron density difference of 1 in the centered rectangular lattice.(a) Iso-surface ofELF=0.8. (b, c) Iso-surface of electron density difference of Δρ=±10-6 Å-3 for positive (b) and for negative (c) values of Δρ.
Supplementary References:
1.Kresse, G. Hafner,J. Ab initio molecular dynamics for liquid metals.Phys. Rev. B47, 558-561 (1993).
2.Kresse, G. Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.Phys. Rev. B54, 11169-11186 (1996).
3.Perdew, J. P., Burke, K. Ernzerhof, M. Generalized Gradient Approximation Made Simple.Phys. Rev. Lett.77, 3865-3868 (1996).
4.Blöchl, P. E.Projector augmented-wave method.Phys. Rev. B50, 17953-17979 (1994).
5.Kresse, G. Joubert, D.From ultrasoft pseudopotentials to the projector augmented-wave method.Phys. Rev. B59, 1758-1775 (1999).
6.Medeiros, P. V. C., Gueorguiev, G. K. Stafström, S. Benzene, coronene, and circumcoronene adsorbed on gold, and a gold cluster adsorbed on graphene: Structural and electronic properties. Phys. Rev. B85, 205423 (2012).
7.Langreth, D. C.et al.Van der Waals density functional theory with applications. Int.J.QuantumChem.101, 599-610 (2005).
8.Dion, M., Rydberg, H., Schröder, E., Langreth, D. C. Lundqvist, B. I.Van der Waals Density Functional for General Geometries. Phys. Rev. Lett.92, 246401 (2004).
9.Nion, A., Jiang, P., Popoff, A. & Fichou, D. Rectangular Nanostructuring of Au(111) Surfaces by Self-Assembly of Size-Selected Thiacrown Ether Macrocycles. J. Am. Chem. Soc.129, 2450-2451 (2007).
10.Bashir, A. et al. Selenium as a Key Element for Highly Ordered Aromatic Self-Assembled Monolayers.Angew. Chem. Int. Ed.47, 5250 (2008).
11.Hohman, J. N. et al.Dynamic Double Lattice of 1-Adamantaneselenolate Self-Assembled Monolayers on Au{111}. J. Am. Chem. Soc.133, 19422-19431 (2011).
12.Chernichenko, K. Y. et al. “Sulflower”: A New Form of Carbon Sulfide Angew. Chem. Int. Ed.45, 7367-7370 (2006).
13.Fujimoto, T., Suizu, R., Yoshikawa, H. & Awaga, K. Molecular, Crystal, and Thin-Film Structures of Octathio[8]circulene: Release of Antiaromatic Molecular Distortion and Lamellar Structure of Self-Assembling Thin Films. Chem. Eur. J.14, 6053-6056 (2008).
14.Zhang, Y.Y., Du, S.X. & Gao, H.-J. Binding configuration, electronic structure, and magnetic properties of metal phthalocyanines on a Au(111) surface studied with ab initio calculations. Phys. Rev. B84, 125446 (2011)
1