Electrostatic potentials of atomic nanostructures at metal surfaces quantified by scanning quantum dot microscopy.

The discrete and charge-separated nature of matter - electrons and nuclei - results in local electrostatic fields that are ubiquitous in nanoscale structures and relevant in catalysis, nanoelectronics and quantum nanoscience. Surface-averaging techniques provide only limited experimental access to these potentials, which are determined by the shape, material, and environment of the nanostructure. Here, we image the potential over adatoms, chains, and clusters of Ag and Au atoms assembled on Ag(111) and quantify their surface dipole moments. By focusing on the total charge density, these data establish a benchmark for theory. Our density functional theory calculations show a very good agreement with experiment and allow a deeper analysis of the dipole formation mechanisms, their dependence on fundamental atomic properties and on the shape of the nanostructures. We formulate an intuitive picture of the basic mechanisms behind dipole formation, allowing better design choices for future nanoscale systems such as single-atom catalysts.

Accurate quantification of the stability of the perylene-tetracarboxylic dianhydride on Au(111) molecule-surface interface.

Studying inorganic/organic hybrid systems is a stepping stone towards the design of increasingly complex interfaces. A predictive understanding requires robust experimental and theoretical tools to foster trust in the obtained results. The adsorption energy is particularly challenging in this respect, since experimental methods are scarce and the results have large uncertainties even for the most widely studied systems. Here we combine temperature-programmed desorption (TPD), single-molecule atomic force microscopy (AFM), and nonlocal density-functional theory (DFT) calculations, to accurately characterize the stability of a widely studied interface consisting of perylene-tetracarboxylic dianhydride (PTCDA) molecules on Au(111). This network of methods lets us firmly establish the adsorption energy of PTCDA/Au(111) via TPD (1.74 ± 0.10 eV) and single-molecule AFM (2.00 ± 0.25 eV) experiments which agree within error bars, exemplifying how implicit replicability in a research design can benefit the investigation of complex materials properties.

Design Principles for Metastable Standing Molecules.

Molecular nanofabrication with a scanning probe microscope (SPM) is a promising route toward the prototyping of metastable functional molecular structures and devices which do not form spontaneously. The aspect of mechanical stability is crucial for such structures, especially if they extend into the third dimension vertical to the surface. A prominent example is freestanding molecules fabricated on a metal which can function as field emitters or electric field sensors. Improving the stability of such molecular configurations is an optimization task involving many degrees of freedom and therefore best tackled by computational nanostructure design. Here, we use density functional theory to study 3,4,9,10-perylene-tetracarboxylic dianhydride (PTCDA) standing on the Ag(111) surface as well as on the tip of a scanning probe microscope. We cast our results into a simple set of design principles for such metastable structures, the validity of which we subsequently demonstrate in two computational case studies. Our work proves the capabilities of computational nanostructure design in the field of metastable molecular structures and offers the intuition needed to fabricate new devices without tedious trial and error.

The stabilization potential of a standing molecule.

The part-by-part assembly of functional nanoscale machinery is a central goal of nanotechnology. With the recent fabrication of an isolated standing molecule with a scanning probe microscope, the third dimension perpendicular to the surface will soon become accessible to molecule-based construction. Beyond the flatlands of the surface, a wealth of structures and functionalities is waiting for exploration, but issues of stability are becoming more critical. Here, we combine scanning probe experiments with ab initio potential energy calculations to investigate the thermal stability of a prototypical standing molecule. We reveal its generic stabilization mechanism, a fine balance between covalent and van der Waals interactions including the latter's long-range screening by many-body effects, and find a remarkable agreement between measured and calculated stabilizing potentials. Beyond their relevance for the design and construction of three-dimensional molecular devices at surfaces, our results also indicate that standing molecules may serve as tunable mechanical gigahertz oscillators.

Communication: Self-assembly of a model supramolecular polymer studied by replica exchange with solute tempering.

Conventional molecular-dynamics (cMD) simulation has a well-known limitation in accessible time and length scales, and thus various enhanced sampling techniques have been proposed to alleviate the problem. In this paper, we explore the utility of replica exchange with solute tempering (REST) (i.e., a variant of Hamiltonian replica exchange methods) to simulate the self-assembly of a supramolecular polymer in explicit solvent and compare the performance with temperature-based replica exchange MD (T-REMD) as well as cMD. As a test system, we consider a relatively simple all-atom model of supramolecular polymerization (namely, benzene-1,3,5-tricarboxamides in methylcyclohexane solvent). Our results show that both REST and T-REMD are able to predict highly ordered polymer structures with helical H-bonding patterns, in contrast to cMD which completely fails to obtain such a structure for the present model. At the same time, we have also experienced some technical challenge (i.e., aggregation-dispersion transition and the resulting bottleneck for replica traversal), which is illustrated numerically. Since the computational cost of REST scales more moderately than T-REMD, we expect that REST will be useful for studying the self-assembly of larger systems in solution with enhanced rearrangement of monomers.

Binary functionalization of H:Si(111) surfaces by alkyl monolayers with different linker atoms enhances monolayer stability and packing.

Alkyl monolayer modified Si forms a class of inorganic-organic hybrid materials with applications across many technologies such as thin-films, fuel/solar-cells and biosensors. Previous studies have shown that the linker atom, through which the monolayer binds to the Si substrate, and any tail group in the alkyl chain, can tune the monolayer stability and electronic properties. In this paper we study the H:Si(111) surface functionalized with binary SAMs: these are composed of alkyl chains that are linked to the surface by two different linker groups. Aiming to enhance SAM stability and increase coverage over singly functionalized Si, we examine with density functional theory simulations that incorporate vdW interactions, a range of linker groups which we denote as -X-(alkyl) with X = CH2, O(H), S(H) or NH(2) (alkyl = C6 and C12 chains). We show how the stability of the SAM can be enhanced by adsorbing alkyl chains with two different linkers, e.g. Si-[C, NH]-alkyl, through which the adsorption energy is increased compared to functionalization with the individual -X-alkyl chains. Our results show that it is possible to improve stability and optimum coverage of alkyl functionalized SAMs linked through a direct Si-C bond by incorporating alkyl chains linked to Si through a different linker group, while preserving the interface electronic structure that determines key electronic properties. This is important since any enhancement in stability and coverage to give more densely packed monolayers will result in fewer defects. We also show that the work function can be tuned within the interval of 3.65-4.94 eV (4.55 eV for bare H:Si(111)).

Density functional theory with van der waals corrections study of the adsorption of alkyl, alkylthiol, alkoxyl, and amino-alkyl chains on the H:Si(111) surface.

Surface modification of silicon with organic monolayers tethered to the surface by different linkers is an important process in realizing future miniaturized electronic and sensor devices. Understanding the roles played by the nature of the linking group and the chain length on the adsorption structures and stabilities of these assemblies is vital to advance this technology. This paper presents a density functional theory (DFT) study of the hydrogen passivated Si(111) surface modified with alkyl chains of the general formula H:Si-(CH2)n-CH2 and H:Si-X-(CH2)n-CH3, where X = NH, O, S and n = (0, 1, 3, 5, 7, 9, 11), at half coverage. For (X)-hexane and (X)-dodecane functionalization, we also examined various coverages up to full monolayer grafting in order to validate the result of half covered surface and the linker effect on the coverage. We find that it is necessary to take into account the van der Waals interaction between the alkyl chains. The strongest binding is for the oxygen linker, followed by S, N, and C, irrespective of chain length. The result revealed that the sequence of the stability is independent of coverage; however, linkers other than carbon can shift the optimum coverage considerably and allow further packing density. For all linkers apart from sulfur, structural properties, in particular, surface-linker-chain angles, saturate to a single value once n > 3. For sulfur, we identify three regimes, namely, n = 0-3, n = 5-7, and n = 9-11, each with its own characteristic adsorption structures. Where possible, our computational results are shown to be consistent with the available experimental data and show how the fundamental structural properties of modified Si surfaces can be controlled by the choice of linking group and chain length.