Exploring in-plane interactions beside an adsorbed molecule with lateral force microscopy.

Atomic force microscopy with a CO-functionalized tip can be used to directly image the internal structure of a planar molecule and to characterize chemical bonds. However, hydrogen atoms usually cannot be directly observed due to their small size. At the same time, these atoms are highly important, since they can direct on-surface chemical reactions. Measuring in-plane interactions at the sides of PTCDA (3,4,9,10-perylenetetracarboxylic dianhydride) molecules with lateral force microscopy allowed us to directly identify hydrogen atoms via their repulsive signature, which we confirmed with a model incorporating radially symmetric atomic interactions. Additional features were observed in the force data and could not be explained by H-bonding of the CO tip with the PTCDA sides. Instead, they are caused by electrostatic interaction of the large dipole of the metal apex, which we verified with density functional theory. This calculation allowed us to estimate the strength of the dipole at the metal tip apex. To further confirm that this dipole generally affects measurements on weakly polarized systems, we investigated the archetypical surface adsorbate of a single CO molecule. We determined the radially symmetric atomic interaction to be valid over a large solid angle of 5.4 sr, corresponding to 82°. We therefore find that in both the PTCDA and CO systems, the underlying interaction preventing direct observations of H-bonding and causing a collapse of the radially symmetric model is the dipole at the metal apex, which plays a significant role when approaching closer than standard imaging heights.

Autonomous Single-Molecule Manipulation Based on Reinforcement Learning.

Building nanostructures one-by-one requires precise control of single molecules over many manipulation steps. The ideal scenario for machine learning algorithms is complex, repetitive, and time-consuming. Here, we show a reinforcement learning algorithm that learns how to control a single dipolar molecule in the electric field of a scanning tunneling microscope. Using about 2250 iterations to train, the algorithm learned to manipulate the molecule toward specific positions on the surface. Simultaneously, it generates physical insights into the movement as well as orientation of the molecule, based on the position where the electric field is applied relative to the molecule. This reveals that molecular movement is strongly inhibited in some directions, and the torque is not symmetric around the dipole moment.

Note: Reproducibility of potential energy surfaces of organic/metal interfaces on the example of PTCDA on Ag(111).

We recently published a benchmark study of common local, semi-local, and non-local exchange correlation functionals in combination with various van der Waals (vdW) corrections, where we investigated the reproducibility of the potential energy surface of perylenetetracarboxylic dianhydride on Ag(111). This Note presents an additional benchmark of the recently developed non-local many body dispersion (MBD-NL) vdW correction, coupled with the Perdew-Burke-Ernzerhof (PBE) functional. We find that this computation method shows similar performance as the established approaches. Notably, it yields very similar results as PBE + MBD.

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice.

The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.

Reproducibility of potential energy surfaces of organic/metal interfaces on the example of PTCDA on Ag(111).

Molecular adsorption at organic/metal interfaces depends on a range of mechanisms: covalent bonds, charge transfer, Pauli repulsion, and van der Waals (vdW) interactions shape the potential energy surface (PES), making it key to understanding organic/metal interfaces. Describing such interfaces with density functional theory requires carefully selecting the exchange correlation (XC) functional and vdW correction scheme. To explore the reproducibility of the PES with respect to the choice of method, we present a benchmark of common local, semi-local, and non-local XC functionals in combination with various vdW corrections. We benchmark these methods using perylene-tetracarboxylic dianhydride on Ag(111), one of the most frequently studied organic/metal interfaces. For each method, we determine the PES using a Gaussian process regression algorithm, which requires only about 50 density functional theory calculations as input. This allows a detailed analysis of the PESs' features, such as the positions and energies of minima and saddle points. Comparing the results from different combinations of XC functionals and vdW corrections enables us to identify trends and differences between the approaches. PESs for different computation methods are in qualitative agreement but also display significant quantitative differences. In particular, the lateral positions of adsorption geometries agree well with experiment, while adsorption heights, energies, and barriers show larger discrepancies.

NIR Phosphorescent Intramolecularly Bridged Benzoporphyrins and Their Application in Oxygen-Compensated Glucose Optode.

A glucose optode measuring the internal oxygen gradient is presented. The multilayer biosensor is composed of (i) analyte-impermeable transparent support, (ii) first oxygen-sensing layer combined with an enzymatic layer, (iii) diffusion barrier, and (iv) second oxygen-sensing layer. To make this design suitable for measurement in subcutaneous tissue, a pair of NIR phosphorescent indicators with very different spectral properties is chosen. Combination of a conventional Pt(II) tetrabenzoporphyrin dye (absorption and emission maxima at 617 and 772 nm, respectively) used in the first layer and a new intramolecularly bridged Pt(II) complex (absorption and emission maxima at 673 and 872 nm, respectively) in the second layer enables efficient separation of both emission signals. This specially designed dye class is accessible via Scholl-reaction from tetraphenyltetrabenzoporphyrin complexes. For the first time, the new optode allows simultaneous glucose and oxygen measurement in a single spot and therefore accurate compensation of oxygen heterogeneities resulting from fluctuations in the tissue. The presented material covers the dynamic ranges from 0 to 150 hPa O2 and from 0 to 360 mg/dL (20 mM) glucose (at 37 °C).

Structure Prediction for Surface-Induced Phases of Organic Monolayers: Overcoming the Combinatorial Bottleneck.

Structure determination and prediction pose a major challenge to computational material science, demanding efficient global structure search techniques tailored to identify promising and relevant candidates. A major bottleneck is the fact that due to the many combinatorial possibilities, there are too many possible geometries to be sampled exhaustively. Here, an innovative computational approach to overcome this problem is presented that explores the potential energy landscape of commensurate organic/inorganic interfaces where the orientation and conformation of the molecules in the tightly packed layer is close to a favorable geometry adopted by isolated molecules on the surface. It is specifically designed to sample the energetically lowest lying structures, including the thermodynamic minimum, in order to survey the particularly rich and intricate polymorphism in such systems. The approach combines a systematic discretization of the configuration space, which leads to a huge reduction of the combinatorial possibilities with an efficient exploration of the potential energy surface inspired by the Basin-Hopping method. Interfacing the algorithm with first-principles calculations, the power and efficiency of this approach is demonstrated for the example of the organic molecule TCNE (tetracyanoethylene) on Au(111). For the pristine metal surface, the global minimum structure is found to be at variance with the geometry found by scanning tunneling microscopy. Rather, our results suggest the presence of surface adatoms or vacancies that are not imaged in the experiment.

Exploring the driving forces behind the structural assembly of biphenylthiolates on Au(111).

In this contribution, we use dispersion-corrected density functional theory to study inter- and intramolecular interactions in a prototypical self-assembled monolayer (SAM) consisting of biphenylthiolates bonded to Au(111) via thiolate groups. The goal is to identify the nature of the interactions that drive the monolayer into a specific conformation. Particular focus is laid on sampling realistic structures rather than high symmetry model configurations. This is achieved by studying conceptually different local minimum structures of the SAM that are obtained via exploring the potential energy surface from systematically varied starting geometries. The six obtained packing motifs differ in the relative arrangement of the two molecules in the unit cell (co-planar versus herringbone) and in the intramolecular configuration (twisted versus planar rings). We find that van der Waals interactions within the organic adsorbate and between the adsorbate and substrate are the main reason that these molecular assemblies can form stable structures at all. The van der Waals interactions are, however, very similar for all observed motifs; by analyzing various types of interactions in the course of three notional SAM-formation steps, we find that the main driving force stabilizing the actual global minimum structure originates from electrostatic interactions between the molecules.

Partial restoration of aromaticity of pentacene-5,7,12,14-tetrone on Cu(111).

The π-conjugation of organic molecules can be strongly influenced when functional groups are added to a molecule, for example when pentacene is converted into pentacene-5,7,12,14-tetrone (P4O) by substitution of four H-atoms with four O-atoms, leading to four CO double bonds. In fact, although free P4O resembles the parent hydrocarbon pentacene structurally at a first glance, its electronic properties differ drastically and can be more accurately described by three benzene units connected via four carbonyl groups. If P4O is deposited onto Cu(111), the electronic interaction across the interface has previously been reported to fully restore the π-conjugation through a weakening of the CO double bonds and a redistribution of electrons, both of which have been explained with the model of surface-induced aromatic stabilization. Here, we observe for the case of P4O on Cu(111) that the molecule does not exhibit full π-conjugation upon interaction with the surface, likely because of the special electronic nature of the hybridized P4O on Cu(111). Our results are derived from CO-functionalized noncontact atomic force microscopy measurements in combination with dispersion-corrected density functional theory calculations yielding bond lengths and molecular geometries. To characterize the aromaticity, we apply the harmonic oscillator model of aromaticity.

Kinetic Trapping of Charge-Transfer Molecules at Metal Interfaces.

Despite the common expectation that conjugated organic molecules on metals adsorb in a flat-lying layer, several recent studies have found coverage-dependent transitions to upright-standing phases, which exhibit notably different physical properties. In this work, we argue that from an energetic perspective, thermodynamically stable upright-standing phases may be more common than hitherto thought. However, for kinetic reasons, this phase may often not be observed experimentally. Using first-principles kinetic Monte Carlo simulations, we find that the structure with lower molecular density is (almost) always formed first, reminiscent of Ostwald's rule of stages. The phase transitions to the upright-standing phase are likely to be kinetically hindered under the conditions typically used in surface science. The simulation results are experimentally confirmed for the adsorption of tetracyanoethylene on Cu(111) using infrared and X-ray photoemission spectroscopy. Investigating both the role of the growth conditions and the energetics of the interface, we find that the time for the phase transition is determined mostly by the deposition rate and, thus, is mostly independent of the nature of the molecule.

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.

Space-charge transfer in hybrid inorganic-organic systems.

We discuss density functional theory calculations of hybrid inorganic-organic systems that explicitly include the global effects of doping (i.e., position of the Fermi level) and the formation of a space-charge layer. For the example of tetrafluoro-tetracyanoquinodimethane on the ZnO(0001[over ¯]) surface we show that the adsorption energy and electron transfer depend strongly on the ZnO doping. The associated work function changes are large, for which the formation of space-charge layers is the main driving force. The prominent doping effects are expected to be quite general for charge-transfer interfaces in hybrid inorganic-organic systems and important for device design.

Large work function reduction by adsorption of a molecule with a negative electron affinity: pyridine on ZnO(1010).

Using thermal desorption and photoelectron spectroscopy to study the adsorption of pyridine on ZnO(1010), we find that the work function is significantly reduced from 4.5 eV for the bare ZnO surface to 1.6 eV for one monolayer of adsorbed pyridine. Further insight into the interface morphology and binding mechanism is obtained using density functional theory. Although semilocal density functional theory provides unsatisfactory total work functions, excellent agreement of the work function changes is achieved for all coverages. In a closed monolayer, pyridine is found to bind to every second surface Zn atom. The strong polarity of the Zn-pyridine bond and the molecular dipole moment act cooperatively, leading to the observed strong work function reduction. Based on simple alignment considerations, we illustrate that even larger work function modifications should be achievable using molecules with negative electron affinity. We expect the application of such molecules to significantly reduce the electron injection barriers at ZnO/organic heterostructures.

Impact of Static Distortion Waves on Superlubricity.

Friction is a major source of energy loss in mechanical devices. This energy loss may be minimized by creating interfaces with extremely reduced friction, i.e., superlubricity. Conventional wisdom holds that incommensurate interface structures facilitate superlubricity. Accurately describing friction necessitates the precise modeling of the interface structure. This, in turn, requires the use of accurate first-principles electronic structure methods, especially when studying organic/metal interfaces, which are highly relevant due to their tunability and propensity to form incommensurate structures. However, the system size required to calculate incommensurate structures renders such calculations intractable. As a result, studies of incommensurate interfaces have been limited to very simple model systems or strongly simplified methodology. We overcome this limitation by developing a machine-learned interatomic potential that is able to determine energies and forces for structures containing thousands to tens of thousands of atoms with an accuracy comparable to conventional first-principles methods but at a fraction of the cost. Using this approach, we quantify the breakdown of superlubricity in incommensurate structures due to the formation of static distortion waves. Moreover, we extract design principles to engineer incommensurate interface systems where the formation of static distortion waves is suppressed, which facilitates low friction coefficients.

Polymorphism mediated by electric fields: a first principles study on organic/inorganic interfaces.

Organic/inorganic interfaces are known to exhibit rich polymorphism, where different polymorphs often possess significantly different properties. Which polymorph forms during an experiment depends strongly on environmental parameters such as deposition temperature and partial pressure of the molecule to be adsorbed. To prepare desired polymorphs these parameters are varied. However, many polymorphs are difficult to access within the experimentally available temperature-pressure ranges. In this contribution, we investigate how electric fields can be used as an additional lever to make certain structures more readily accessible. On the example of tetracyanoethylene (TCNE) on Cu(111), we analyze how electric fields change the energy landscape of interface systems. TCNE on Cu(111) can form either lying or standing polymorphs, which exhibit significantly different work functions. We combine first-principles calculations with a machine-learning based structure search algorithm and ab initio thermodynamics to demonstrate that electric fields can be exploited to shift the temperature of the phase transition between standing and lying polymorphs by up to 100 K.

Interplay of Adsorption Geometry and Work Function Evolution at the TCNE/Cu(111) Interface.

The adsorption of organic electron acceptors on metal surfaces is a powerful way to change the effective work function of the substrate through the formation of charge-transfer-induced dipoles. The work function of the interfaces is hence controlled by the redistribution of charges upon adsorption of the organic layer, which depends not only on the electron affinity of the organic material but also on the adsorption geometry. As shown in this work, the latter dependence controls the work function also in the case of adsorbate layers exhibiting a mixture of various adsorption geometries. Based on a combined experimental (core-level and infrared spectroscopy) and theoretical (density functional theory) study for tetracyanoethylene (TCNE) on Cu(111), we find that TCNE adsorbs in at least three different orientations, depending on TCNE coverage. At low coverage, flat lying TCNE dominates, as it possesses the highest adsorption energy. At a higher coverage, additionally, two different standing orientations are found. This is accompanied by a large increase in the work function of almost 3 eV at full monolayer coverage. Our results suggest that the large increase in work function is mainly due to the surface dipole of the free CN groups of the standing molecules and less dependent on the charge-transfer dipole of the differently oriented and charged molecules. This, in turn, opens new opportunities to control the work function of interfaces, e.g., by synthetic modification of the adsorbates, which may allow one to alter the adsorption geometries of the molecules as well as their contributions to the interface dipoles and, hence, the work function.

Dimensionality effects in the electronic structure of organic semiconductors consisting of polar repeat units.

In conjugated organic molecules, excitation gaps typically decrease reciprocally with increasing the number of repeat units, n. This usually holds for individual molecules as well as for the corresponding bulk materials. Here, we show using density-functional theory calculations that a qualitatively different evolution is found for layers built from molecules consisting of polar repeat units. Whereas a 1/n-dependence is still observed in the case of isolated polar molecules, the global gap decreases essentially linearly with n in the corresponding 2D-periodic systems and vanishes beyond a certain molecular length, with the frontier states being localized at opposite ends of the layer. The latter is accompanied by a saturation of the dipole moment per molecule, an effect not observed in the isolated polar molecules. Interestingly, in both cases the limit of the gap for long (but finite) molecules differs qualitatively from that of infinite length obtained in 1D-periodic and 3D-periodic calculations, the latter serving as models for polymers and the bulk. We rationalize these dimensionality effects as a consequence of the potential gradient within the finite-length layers. They arise from the collective action of intra-molecular dipoles in the 2D periodic layers and can be traced back to surface effects.

From a bistable adsorbate to a switchable interface: tetrachloropyrazine on Pt(111).

Virtually all organic (opto)electronic devices rely on organic/inorganic interfaces with specific properties. These properties are, in turn, inextricably linked to the interface structure. Therefore, a change in structure can introduce a shift in function. If this change is reversible, it would allow constructing a switchable interface. We accomplish this with tetrachloropyrazine on Pt(111), which exhibits a double-well potential with a chemisorbed and a physisorbed minimum. These minima have significantly different adsorption geometries allowing the formation of switchable interface structures. Importantly, these structures facilitate different work function changes and coherent fractions (as would be obtained from X-ray standing wave measurements), which are ideal properties to read out the interface state. We perform surface structure search using a modified version of the SAMPLE approach and account for thermodynamic conditions using ab initio thermodynamics. This allows investigating millions of commensurate as well as higher-order commensurate interface structures. We identify three different classes of structures exhibiting different work function changes and coherent fractions. Using temperature and pressure as handles, we demonstrate the possibility of reversible switching between those different classes, creating a dynamic interface for potential applications in organic electronics.

Role of Adatoms for the Adsorption of F4TCNQ on Au(111).

Organic adlayers on inorganic substrates often contain adatoms, which can be incorporated within the adsorbed molecular species, forming two-dimensional metal-organic frameworks at the substrate surface. The interplay between native adatoms and adsorbed molecules significantly changes various adlayer properties such as the adsorption geometry, the bond strength between the substrate and the adsorbed species, or the work function at the interface. Here, we use dispersion-corrected density functional theory to gain insight into the energetics that drive the incorporation of native adatoms within molecular adlayers based on the prototypical, experimentally well-characterized system of F4TCNQ on Au(111). We explain the adatom-induced modifications in the adsorption geometry and the adsorption energy based on the electronic structure and charge transfer at the interface.

Toward Targeted Kinetic Trapping of Organic-Inorganic Interfaces: A Computational Case Study.

Properties of inorganic-organic interfaces, such as their interface dipole, strongly depend on the structural arrangements of the organic molecules. A prime example is tetracyanoethylene (TCNE) on Cu(111), which shows two different phases with significantly different work functions. However, the thermodynamically preferred phase is not always the one that is best suited for a given application. Rather, it may be desirable to selectively grow a kinetically trapped structure. In this work, we employ density functional theory and transition state theory to discuss under which conditions such a kinetic trapping might be possible for the model system of TCNE on Cu. Specifically, we want to trap the molecules in the first layer in a flat-lying orientation. This requires temperatures that are sufficiently low to suppress the reorientation of the molecules, which is thermodynamically more favorable for high dosages, but still high enough to enable ordered growth through diffusion of molecules. On the basis of the temperature-dependent diffusion and reorientation rates, we propose a temperature range at which the reorientation can be successfully suppressed.