Two-dimensional short-range disordered crystalline networks from flexible molecular modules

Anisotropic functionalized platelets: percolation, porosity and network properties

Anisotropic functionalized platelets are able to model the assembly behaviour of molecular systems in two dimensions thanks to the unique combination of steric and bonding constraints. The assembly scenarios can vary from open to close-packed crystals, finite clusters and chains, according to the features of the imposed constraints. In this work, we focus on the assembly of equilibrium networks. These networks can be seen as disordered, porous monolayers and can be of interest for instance in nano-filtration and optical applications. We investigate the formation and properties of two dimensional networks from shape anisotropic colloids functionalized with four patches. We characterize the connectivity properties, the typical local bonding motives, as well as the geometric features of the emerging networks for a large variety of different systems. Our results show that networks of shape anisotropic colloids assemble into highly versatile network topologies, that may be utilized for applications at the nanoscale.

On-Surface Synthesis of Multiple Cu Atom-Bridged Organometallic Oligomers

A metal-metal bond between coordination complexes has the nature of a covalent bond in hydrocarbons. While bimetallic and trimetallic compounds usually have three-dimensional structures in solution, the high directionality and robustness of the bond can be applied for on-surface syntheses. Here, we present a systematic formation of complex organometallic oligomers on Cu(111) through sequential ring opening of 11,11,12,12-tetraphenyl-1,4,5,8-tetraazaanthraquinodimethane and bonding of phenanthroline derivatives by multiple Cu atoms. A detailed characterization with a combination of scanning tunneling microscopy and density functional theory calculations revealed the role of the Cu adatoms in both enantiomers of the chiral oligomers. Furthermore, we found sufficient strength of the bonds against sliding friction by manipulating the oligomers up to a hexamer. This finding may help to increase the variety of organometallic nanostructures on surfaces.

Order-Disorder Transition of Two-Dimensional Molecular Networks through a Stoichiometric Design

Materials with disordered structures may exhibit interesting properties. Metal-organic frameworks (MOFs) are a class of hybrid materials composed of metal nodes and coordinating organic linkers. Recently, there has been growing interest in MOFs with structural disorder and the investigations of amorphous structures on surfaces. Herein, we demonstrate a bottom-up method to construct disordered molecular networks on metal surfaces by selecting two organic molecule linkers with the same symmetry but different sizes for preparing two-component samples with different stoichiometric ratios. The amorphous networks are directly imaged by scanning tunneling microscopy under ultrahigh vacuum with a submolecular resolution, allowing us to quantify its degree of disorder and other structural properties. Furthermore, we resort to molecular dynamics simulations to understand the formation of the amorphous metal-organic networks. The results may advance our understanding of the mechanism of formation of monolayer molecular networks with structural disorders, facilitating the design and exploration of amorphous MOF materials with intriguing properties.

Predicting Molecular Self-Assembly on Metal Surfaces Using Graph Neural Networks Based on Experimental Data Sets

The application of supramolecular chemistry on solid surfaces has received extensive attention in the past few decades. To date, combining experiments with quantum mechanical or molecular dynamic methods represents the key strategy to explore the molecular self-assembled structures, which is, however, often laborious. Recently, machine learning (ML) has become one of the most exciting tools in material research, allowing for both efficiency and accuracy in predicting molecular properties. In this work, we constructed a graph neural network to predict the self-assembly of functional polycyclic aromatic hydrocarbons (PAHs) on metal surfaces. Using scanning tunneling microscopy (STM), we characterized the self-assembled nanostructures of a homologous series of PAH molecules on different metal surfaces to construct an experimental data set for model training. Compared with traditional ML algorithms, our model exhibits better predictive performance. Finally, the generalization of the model is further verified by comparing the ML predictions and experimental results of different functionalized molecule. Our results demonstrate training experimental data sets to produce a predictive ML model of molecular self-assembly with generalization performance, which allows for the predictive design of nanostructures with functional molecules.

High-Throughput Preparation of Supramolecular Nanostructures on Metal Surfaces

One of the contemporary challenges in materials science lies in the rapid materials screening and discovery. Experimental sample libraries can be generated by high-throughput parallel synthesis to map the composition space for rapid material discoveries. Molecular self-assembly on surfaces has proved a useful way to construct nanostructures with interesting topologies or properties. Despite the strong dependence of molecular stoichiometry on the structures, high-throughput preparations of supramolecular surface nanostructures have been far less explored. Here, by integrating a physical mask into the standard ultra-high-vacuum (UHV) molecular preparation system we show a high-throughput approach for preparing supramolecular nanostructures of continuous composition spreads on metal surfaces. The spatially addressable sample libraries of supramolecular self-assemblies are characterized by high-resolution scanning probe microscopy. We could explore different binary nanostructures of varying molecular ratios on one single substrate. Moreover, we use the minimum spanning tree approach to qualitatively and quantitatively study the structural properties of the formed nanostructures. This high-throughput approach may accelerate the screening and exploration of surface-supported, low-dimensional nanostructures not limited to supramolecular interactions.

Monte Carlo Simulations of the Metal-Directed Self-Assembly of Y-Shaped Positional Isomers

The rational fabrication of low-dimensional materials with a well-defined topology and functions is an incredibly important aspect of nanotechnology. In particular, the on-surface synthesis (OSS) methods based on the bottom-up approach enable a facile construction of sophisticated molecular architectures unattainable by traditional methods of wet chemistry. Among such supramolecular constructs, especially interesting are the surface-supported metal–organic networks (SMONs), composed of low-coordinated metal atoms and π-aromatic bridging linkers. In this work, the lattice Monte Carlo (MC) simulation technique was used to extract the chemical information encoded in a family of Y-shaped positional isomers co-adsorbed with trivalent metal atoms on a flat metallic surface with (111) geometry. Depending on the intramolecular distribution of active centers (within the simulated molecular bricks, we observed a metal-directed self-assembly of two-dimensional (2D) openwork patterns, aperiodic mosaics, and metal–organic ladders. The obtained theoretical findings could be especially relevant for the scanning tunneling microscopy (STM) experimentalists interested in a surface-assisted construction of complex nanomaterials stabilized by directional coordination bonds.

Flexible Alkali-Halogen Bonding in Two Dimensional Alkali-Metal Organic Frameworks

On-surface fabrication of two-dimensional (2D) metal organic frameworks (MOFs) has been continuously attracting attentions for years. However, the synthesis of 2D MOFs with large-amplitude flexibility was rarely carried out since the bonding configurations in the coordination nodes are typically highly directional. Here we demonstrate that single alkali ions, which are fully isotropic in ionic bonding, can act as pivot joints for constructing tunable 2D MOFs by bonding to dihalogen groups in organic molecules. We take 2,3,6,7,10,11-hexabromotriphenylene, a 3-fold polycyclic molecule with three ortho-dibromo groups, and sodium (Na) atoms as a model system and successfully construct Na-based MOFs on Au(111) surface. The deflection angle of the Na coordination nodes is variable in an unprecedentedly large range between ±36° that allows the construction of multiple 2D MOF architectures. Such a flexible alkali-halogen bonding may provide a unique toolbox for designing and constructing more tunable MOFs by choosing various alkali atoms and halogen moieties.

Rotation in an Enantiospecific Self-Assembled Array of Molecular Raffle Wheels

Tailored nano-spaces can control enantioselective adsorption and molecular motion. We report on the spontaneous assembly of a dynamic system-a rigid kagome network with each pore occupied by a guest molecule-employing solely 2,6-bis(1H-pyrazol-1-yl)pyridine-4-carboxylic acid on Ag(111). The network cavity snugly hosts the chemically modified guest, bestows enantiomorphic adsorption and allows selective rotational motions. Temperature-dependent scanning tunnelling microscopy studies revealed distinct anchoring orientations of the guest unit switching with a 0.95 eV thermal barrier. H-bonding between the guest and the host transiently stabilises the rotating guest, as the flapper on a raffle wheel. Density functional theory investigations unravel the detailed molecular pirouette of the guest and how the energy landscape is determined by H-bond formation and breakage. The origin of the guest's enantiodirected, dynamic anchoring lies in the specific interplay of the kagome network and the silver surface.

Creating supramolecular semiregular Archimedean tilings via gas-mediated deprotonation of a terminal alkyne derivative

Combining surface-confined reactions with supramolecular self-assembly allows the chemical transformation of simple molecular precursors into higher-level tectons to generate complex tessellations with unique structural and functional properties.

Multiscale Analysis of Phase Transformations in Self-Assembled Layers of 4,4'-Biphenyl Dicarboxylic Acid on the Ag(001) Surface

Understanding the nucleation and growth kinetics of thin films is a prerequisite for their large-scale utilization in devices. For self-assembled molecular phases near thermodynamic equilibrium the nucleation-growth kinetic models are still not developed. Here, we employ real-time low-energy electron microscopy (LEEM) to visualize a phase transformation induced by the carboxylation of 4,4'-biphenyl dicarboxylic acid on Ag(001) under ultra-high-vacuum conditions. The initial (α) and transformed (β) molecular phases are characterized in detail by X-ray photoemission spectroscopy, single-domain low-energy electron diffraction, room-temperature scanning tunneling microscopy, noncontact atomic force microscopy, and density functional theory calculations. The phase transformation is shown to exhibit a rich variety of phenomena, including Ostwald ripening of the α domains, burst nucleation of the β domains outside the α phase, remote dissolution of the α domains by nearby β domains, and a structural change from disorder to order. We show that all phenomena are well described by a general growth-conversion-growth (GCG) model. Here, the two-dimensional gas of admolecules has a dual role: it mediates mass transport between the molecular islands and hosts a slow deprotonation reaction. Further, we conclude that burst nucleation is consistent with a combination of rather weak intermolecular bonding and the onset of an additional weak many-body attractive interaction when a molecule is surrounded by its nearest neighbors. In addition, we conclude that Ostwald ripening and remote dissolution are essentially the same phenomenon, where a more stable structure grows at the expense of a kinetically formed, less stable entity via transport through the 2D gas. The proposed GCG model is validated through kinetic Monte Carlo (kMC) simulations.

Coordination-Controlled C-C Coupling Products via ortho-Site C-H Activation

The coordination-restricted ortho-site C-H bond activation and dehydrogenative homocoupling of 4,4'-(1,3-phenylene)dipyridine (1,3-BPyB) and 4,4'-(1,4-phenylene)dipyridine (1,4-BPyB) on different metal surfaces were studied by a combination of scanning tunneling microscopy, noncontact atomic force microscopy, and density functional theory calculations. The coupling products on Cu(111) exhibited certain configurations subject to the spatial restriction of robust two-fold Cu-N coordination bonds. Compared to the V-shaped 1,3-BPyB, the straight backbone of 1,4-BPyB helped to further reduce the variety of reactive products. By utilizing the three-fold coordination of Fe atoms with 1,4-BPyB molecules on Au(111), a large-scale network containing single products was constructed. Our results offer a promising protocol for controllable on-surface synthesis with the aid of robust coordination interactions.

Polymorphic self-assembly of pyrazine-based tectons at the solution-solid interface

Exploring the surface self-assembly of small molecules that act as building blocks (tectons) for complex supramolecular structures is crucial for realizing surface-supported functional molecular devices. Here, we report on the synthesis and surface self-assembly of a new pyrazine-derived molecule with pyridine pendants. Ambient scanning tunneling microscopy investigation at the solution-solid interface reveals polymorphic self-assembly of these molecules on a HOPG substrate. Two different molecular packing structures with equal distribution are observed. Detailed analysis of the STM images emphasizes the crucial role of weak intermolecular hydrogen bonding, and molecule-substrate interactions in the formation of the observed polymorphs. Such weak hydrogen bonding interactions are highly desirable for the formation of modular supramolecular architectures since they can provide sufficiently robust molecular structures and also facilitate error correction.

An experimental and theoretical study of metallorganic coordination networks of tetrahydroxyquinone on Cu(111)

DFT modeling of STM and XAS evidences investigated the adsorption of THQ@Cu(111) that generates different ordered configurations at different temperatures.

Concentration-dependent supramolecular patterns of C3 and C2 symmetric molecules at the solid/liquid interface

Here we report on a scanning tunnelling microscopy (STM) investigation on the self-assembly of C₃- and C₂-symmetric molecules at the solution/graphite interface. 1,3,5-tris((E)-2-(pyridin-4-yl)vinyl)benzene and 1,1,2,2-tetrakis(4-(pyridin-4-yl)phenyl)ethane are used as model systems. These molecules displayed a concentration dependent self-assembly behaviour on graphite, resulting in highly ordered supramolecular structures, which are stabilized jointly by van der Waals substrate-adsorbate interactions and in-plane intermolecular H-bonding. Denser packing is obtained when applying a relatively high concentration solution to the basal plane of the surface whereas a less dense porous network is observed upon lowering the concentration. We show that the molecular conformation does not influence the stability of the self-assembly and a twisted molecule can pack into dense and porous architectures under the concentration effect.

Two-dimensional crystal engineering using halogen and hydrogen bonds: towards structural landscapes

Two-dimensional (2D) crystallization on solid surfaces is governed by a subtle balance of supramolecular and interfacial interactions. However, these subtle interactions often make the prediction of supramolecular structure from the molecular structure impossible. As a consequence, surface-based 2D crystallization has often been studied on a case-by-case basis, which hinders the identification of structure-determining relationships between different self-assembling systems. Here we begin the discussion on such structure-determining relationships by comparing the 2D crystallization of two identical building blocks based on a 1,3,5-tris(pyridine-4-ylethynyl)benzene unit at the solution-solid interface. The concepts of supramolecular synthons and structural landscapes are introduced in the context of 2D crystallization on surfaces to identify common structural elements. The systems are characterized using scanning tunneling microscopy (STM). This strategy involves carrying out minor structural modifications on the parent compound to access supramolecular patterns that are otherwise not obtained. We demonstrate that this chemical perturbation strategy translates equally well for 2D co-crystallization experiments with halogen bond donors yielding porous bi-component networks. The holistic approach described here represents a stepping stone towards gaining predictive power over the 2D crystallization of molecules on solid surfaces.

Turning off the majority-rules effect in two-dimensional hierarchical chiral assembly by introducing a chiral mismatch

Understanding the mechanism in chiral transmission from a single molecule to a supramolecular level is fundamentally important to decipher the nonlinear amplification effect in the two-dimensional (2D) chiral assembly process. In this contribution, we report on the dramatically different nonlinear amplification effect in the chiral co-adsorber induced homochiral assemblies constructed by a series of homologous achiral building blocks on the graphite surface under control of the majority-rules principle. Homologous hexagonal networks are formed for 5-(benzyloxy)-isophthalic acid (BIC) derivatives with different alkyl lengths. While globally homochiral monolayers of BIC-C6 or BIC-C16 can be obtained by using a mixture of chiral co-adsorber 2-octanol with a small enantiomeric excess, such majority-rules principle based nonlinear chiral amplification is inoperative for the assembly of BIC-C10. Molecular mechanistic analysis indicates that BIC-C10 assembly can accommodate a chiral mismatched motif to form long-range ordered but short-range disordered crystalline networks, leading to the co-adsorption of enantiomers without enantioselectivity. The present results shed important insights into the significance of chirality mismatch during chiral transmission and benefit the understanding of chiral communication in a surface monolayer.

Structure and role of metal clusters in a metal-organic coordination network determined by density functional theory

We present a comprehensive theoretical investigation of the structures formed by self-assembly of tetrahydroxybenzene (THB)-derivatives on Cu(111). The THB molecule is known to dehydrogenate completely during annealing, forming a reactive radical which assembles into a close-packed structure or a porous metal-coordinated network depending on the coverage of the system. Here, we present details on how the structures are determined by density functional theory calculations, using scanning tunneling microscopy-derived information on the periodicity. The porous network is based on adatom trimers. By analysing the charge distribution of the structure, it is found that this unusual coordination motif is preferred because it simultaneously provides a good coordination of all oxygen atoms and allows for the formation of a two-dimensional network on the surface.

Influence of the supramolecular order on the electrical properties of 1D coordination polymers based materials

The generation, under self-assembly conditions, of coordination polymers on surface based combinations of a terpyridine-antracene-pyridine based tecton and Co(II) or Pd(II) cations is primarily governed by the coordination geometry of the metal center (octahedral and square planar respectively). While the octahedral Co(II) based polymer self-assembles in insulating films exhibiting randomly oriented crystalline domains, the planarity of Pd(II) based polymers leads to the formation of conductive π-π stacked fibrillar structures exhibiting anisotropically oriented domains. In the latter case, the favorable Pd-Pd and anthracene-anthracene wavefunction overlaps along the fiber direction are responsible for the large electronic couplings between adjacent chains, whereas small electronic couplings are instead found along individual polymer chains. These results provide important guidelines for the design of conductive metal coordination polymers, highlighting the fundamental role of both intra- as well as inter-chain interactions, thus opening up new perspectives towards their application in functional devices.

Atomically Precise Prediction of 2D Self-Assembly of Weakly Bonded Nanostructures: STM Insight into Concentration-Dependent Architectures

A joint experimental and computational study is reported on the concentration-dependant self-assembly of a flat C3 -symmetric molecule on a graphite surface. As a model system a tripodal molecule, 1,3,5-tris(pyridin-3-ylethynyl)benzene, has been chosen, which can adopt either C3h or Cs symmetry when planar, as a result of pyridyl rotation along the alkynyl spacers. Density functional theory (DFT) simulations of 2D nanopatterns with different surface coverage reveal that the molecule can generate different types of self-assembled motifs. The stability of fourteen 2D patterns and the influence of concentration are analyzed. It is found that ordered, densely packed monolayers and 2D porous networks are obtained at high and low concentrations, respectively. A concentration-dependent scanning tunneling microscopy (STM) investigation of this molecular self-assembly system at a solution/graphite interface reveals four supramolecular motifs, which are in perfect agreement with those predicted by simulations. Therefore, this DFT method represents a key step forward toward the atomically precise prediction of molecular self-assembly on surfaces and at interfaces.

Dysprosium-carboxylate nanomeshes with tunable cavity size and assembly motif through ionic interactions

Tunability of cavity size and assembly motif of dysprosium-directed nanomeshes through ionic interactions.

Surface-Supported Robust 2D Lanthanide-Carboxylate Coordination Networks

Lanthanide-based metal-organic compounds and architectures are promising systems for sensing, heterogeneous catalysis, photoluminescence, and magnetism. Herein, the fabrication of interfacial 2D lanthanide-carboxylate networks is introduced. This study combines low- and variable-temperature scanning tunneling microscopy (STM) and X-ray photoemission spectroscopy (XPS) experiments, and density functional theory (DFT) calculations addressing their design and electronic properties. The bonding of ditopic linear linkers to Gd centers on a Cu(111) surface gives rise to extended nanoporous grids, comprising mononuclear nodes featuring eightfold lateral coordination. XPS and DFT elucidate the nature of the bond, indicating ionic characteristics, which is also manifest in appreciable thermal stability. This study introduces a new generation of robust low-dimensional metallosupramolecular systems incorporating the functionalities of the f-block elements.

Lattice-Directed Formation of Covalent and Organometallic Molecular Wires by Terminal Alkynes on Ag Surfaces

Surface reactions of 2,5-diethynyl-1,4-bis(phenylethynyl)benzene on Ag(111), Ag(110), and Ag(100) were systematically explored and scrutinized by scanning tunneling microscopy, molecular mechanics simulations, and density functional theory calculations. On Ag(111), Glaser coupling reaction became dominant, yielding one-dimensional molecular wires formed by covalent bonds. On Ag(110) and Ag(100), however, the terminal alkynes reacted with surface metal atoms, leading to one-dimensional organometallic nanostructures. Detailed experimental and theoretical analyses revealed that such a lattice dependence of the terminal alkyne reaction at surfaces originated from the matching degree between the periodicities of the produced molecular wires and the substrate lattice structures.

Orthogonal insertion of lanthanide and transition-metal atoms in metal-organic networks on surfaces

The orthogonal coordinative properties of tetrapyrrole macrocycles and nitrile ligands have been used in a multistep procedure towards interfacial d-f hetero-bimetallic nanoarchitectures based on a free-base porphyrin derivative functionalized with meso-cyanobiphenylene substituents. Molecular-level scanning tunneling microscopy studies reveal that the porphyrin module alone self-assembles on Ag(111) in a close-packed layer with a square unit cell. Upon co-deposition of Gd atoms, a square-planar motif is formed that reflects the fourfold coordination of CN ligands to the rare-earth centers. The resulting nanoporous network morphology is retained following exposure to a beam of Co atoms, which induces selective porphyrin metalation and ultimately yields a gridlike 2D metallosupramolecular architecture.

Quantum confinement in self-assembled two-dimensional nanoporous honeycomb networks at close-packed metal surfaces

Quantum confinement of a two-dimensional electron gas by supramolecular nanoporous networks is investigated using the boundary elements method based on Green's functions for finite geometries and electron plane wave expansion for periodic systems. The "particle in a box" picture was analyzed for cases with selected symmetries that model previously reported architectures constructed from organic and metal-organic scattering centers confining surface state electrons of Ag(111) and Cu(111). First, by analyzing a series of cases with systematically defined parameters (scattering geometry, potentials, and effective broadening), we demonstrate how the scattering processes affect the properties of the confined electrons. For the features of the local density of states reported by scanning tunneling spectroscopy (STS), we disentangle the contributions of lifetime broadening and splitting of quantum well states due to coupling of neighboring quantum dots. For each system, we analyze the local electron density distribution and relate it to the corresponding band structure as calculated within the plane-wave expansion framework. Then, we address two experimental investigations, where in one case only STS data and in the other case mainly angle-resolved photoemission spectroscopy (ARPES) data were reported. In both cases, the experimental findings can be successfully simulated. Furthermore, the missing information can be complemented because our approach allows to correlate the information obtained by STS with that of ARPES. The combined analysis of several observations suggests that the scattering potentials created by the network originate primarily from the adsorbate-induced changes of the local surface dipole barrier.

Two-dimensional soft supramolecular networks

Soft networks are self-assembled at the solid/liquid interface and characterized by local disorder arising from multivalent flexible intermolecular interactions.

Bestowing structure upon the pores of a supramolecular network

Trigonal molecules compartmentalise the pores of a honeycomb network of 3,4:9,10-tetracarboxylic diimide (PTCDI) and 1,3,5-triazine-2,4,6-triamine (melamine). Extending the 1,3,5-tri(phenylene-ethynylene)benzene core by a phenyl group allows for a well-defined accommodation of the molecule into two symmetry equivalent positions in the pore. The corresponding styryl or phenylene-ethynylene derivatives exceed the pore size and, thus, impede pore modification.

Triangular-shaped molecular random tiling and molecular rotation in two-dimensional glassy networks

Macrocycle-1 molecules can self-assemble into glassy state networks via van der Waals force and form many triangular nanopores in networks. The nanopores can be expressed by triangular tilings, which lead to a particularly rich range of arrangements. Moreover an interesting molecular rotation phenomenon was observed in the glassy networks.

A tale of tails: alkyl chain directed formation of 2D porous networks reveals odd-even effects and unexpected bicomponent phase behavior

Supramolecular self-assembly of suitably functionalized building blocks on surfaces can serve as an excellent test-bed to gain understanding and control over multicomponent self-assembly in more complex matter. Here we employ a powerful combination of scanning tunnelling microscopy (STM) and molecular modeling to uncover two-dimensional (2D) crystallization and mixing behavior of a series of alkylated building blocks based on dehydrobenzo[12]annulene, forming arrays of nanowells. Thorough STM investigation employing high-resolution spatial imaging, use of specially designed marker molecules, statistical analysis and thermal stability measurements revealed rich and complex supramolecular chemistry, highlighting the impact of odd-even effects on the phase behavior. The methodology and analysis presented in this work can be easily adapted to the self-assembly of other alkylated building blocks.

How surface bonding and repulsive interactions cause phase transformations: ordering of a prototype macrocyclic compound on Ag(111)

We investigated the surface bonding and ordering of free-base porphine (2H-P), the parent compound of all porphyrins, on a smooth noble metal support. Our multitechnique investigation reveals a surprisingly rich and complex behavior, including intramolecular proton switching, repulsive intermolecular interactions, and density-driven phase transformations. For small concentrations, molecular-level observations using low-temperature scanning tunneling microscopy clearly show the operation of repulsive interactions between 2H-P molecules in direct contact with the employed Ag(111) surface, preventing the formation of islands. An increase of the molecular coverage results in a continuous decrease of the average intermolecular distance, correlated with multiple phase transformations: the system evolves from an isotropic, gas-like configuration via a fluid-like phase to a crystalline structure, which finally gives way to a disordered layer. Herein, considerable site-specific molecule-substrate interactions, favoring an exclusive adsorption on bridge positions of the Ag(111) lattice, play an important role. Accordingly, the 2D assembly of 2H-P/Ag(111) layers is dictated by the balance between adsorption energy maximization while retaining a single adsorption site counteracted by the repulsive molecule-molecule interactions. The long-range repulsion is associated with a charge redistribution at the 2H-P/Ag(111) interface comprising a partial filling of the lowest unoccupied molecular orbital, resulting in long-range electrostatic interactions between the adsorbates. Indeed, 2H-P molecules in the second layer that are electronically only weakly coupled to the Ag substrate show no repulsive behavior, but form dense-packed islands.