We report on a joint experimental and theoretical study of the ordered structures of melamine molecules formed on the Au( 111)-( 22× 3) surface. Scanning tunneling microscopy (STM) images taken under UHV conditions reveal two distinct monolayers one of which has never been reported before on gold. We also find that one of the structures may serve as a transition region ("domain wall") between islands formed by the other arrangement. Using state-of-the-art density functional calculations in conjunction with a systematic gas-phase analysis based on considering all planar structures melamine molecules can form with each other, we propose atomistic models for both structures and the transition region.

The self-assembly of flat organic molecules on metal surfaces is controlled, apart from the kinetic factors, by the interplay between the molecule-molecule and molecule-surface interactions. These are typically calculated using standard density functional theory within the generalized gradient approximation, which significantly underestimates nonlocal correlations, i.e. van der Waals (vdW) contributions, and thus affects interactions between molecules and the metal surface in the junction. In this paper we address this question systematically for the Au(111) surface and a number of popular flat organic molecules which form directional hydrogen bonds with each other. This is done using the recently developed first-principles vdW-DF method which takes into account the nonlocal nature of electron correlation [M. Dion et al., Phys. Rev. Lett. 2004, 92, 246401]. We report here a systematic study of such systems involving completely self-consistent vdW-DF calculations with full geometry relaxation. We find that the hydrogen bonding between the molecules is only insignificantly affected by the vdW contribution, both in the gas phase and on the gold surface. However, the adsorption energies of these molecules on the surface increase dramatically as compared with the ordinary density functional (within the generalized gradient approximation, GGA) calculations, in agreement with available experimental data and previous calculations performed within approximate or semiempirical models, and this is entirely due to the vdW contribution which provides the main binding mechanism. We also stress the importance of self-consistency in calculating the binding energy by the vdW-DF method since the results of non-self-consistent calculations in some cases may be off by up to 20%. Our calculations still support the usually made assumption of the molecule-surface interaction changing little laterally suggesting that single molecules and their small clusters should be quite mobile at room temperature on the surface. These findings support a gas-phase modeling for some flat metal surfaces, such as Au(111), and flat molecules, at least as a first approximation.

Two molecular phases of the DNA base adenine (A) on a Au(111) surface are observed by using STM under ultrahigh-vacuum conditions. One of these phases is reported for the first time. A systematic approach that considers all possible gas-phase two-dimensional arrangements of A molecules connected by double hydrogen bonds with each other and subsequent ab initio DFT calculations are used to characterize and identify the two phases. The influence of the gold surface on the structure of A assemblies is also discussed. DFT is found to predict a smooth corrugation potential of the gold surface that will enable A molecules to move freely across the surface at room temperature. This conclusion remains unchanged if van der Waals interaction between A and gold is also approximately taken into account. DFT calculations of the A pairs on the Au(111) surface show its negligible effect on the hydrogen bonding between the molecules. These results justify the gas-phase analysis of possible assemblies on flat metal surfaces. Nevertheless, the fact that it is not the most stable gas-phase monolayer that is actually observed on the gold surface indicates that the surface still plays a subtle role, which needs to be properly addressed.

Using a combination of scanning tunneling microscopy (STM) in ultrahigh vacuum (UHV) and a systematic theoretical method based on considering all possible hydrogen bond connections between molecules with subsequent density functional theory (DFT) calculations, we studied supramolecular assemblies of highly symmetrical rectangular PTCDI molecules on the Au(111) surface. We show, using a systematic prediction procedure followed by ab initio density functional calculations, that just over 10 monolayer structures are possible assuming two molecules in the primitive cell, some of which would appear indistinguishable in the STM images. By breaking down these structures into distinct assemblies, we predict six possible phases. Two of these had been observed previously: a canted phase seen on a number of surfaces including Au(111) and a brick wall phase seen so far only on the NaCl(001) surface. Using STM imaging of PTCDI molecules on the Au(111) surface in ultrahigh vacuum, we discovered a completely new “domino” phase, also predicted by our theory, in which molecules attach to each other rather like dominoes, to form squares repeated periodically across the surface. The interaction of the molecules with the gold surface seems to influence the orientation of the phases but not necessarily their stability.