Adatoms underneath single porphyrin molecules on Au(111).

The adsorption of porphyrin derivatives on a Au(111) surface was studied by scanning tunneling microscopy and spectroscopy at low temperatures in combination with density functional theory calculations. Different molecular appearances were found and could be assigned to the presence of single gold adatoms bonded by a coordination bond underneath the molecular monolayer, causing a characteristic change of the electronic structure of the molecules. Moreover, this interpretation could be confirmed by manipulation experiments of individual molecules on and off a single gold atom. This study provides a detailed understanding of the role of metal adatoms in surface-molecule bonding and anchoring and of the appearance of single molecules, and it should prove relevant for the imaging of related molecule-metal systems.

Towards design rules for covalent nanostructures on metal surfaces.

The covalent molecular assembly on metal surfaces is explored, outlining the different types of applicable reactions. Density functional calculations for on-surface reactions are shown to yield valuable insights into specific reaction mechanisms and trends across the periodic table. Finally, it is shown how design rules could be derived for nanostructures on metal surfaces.

Mechanisms of halogen-based covalent self-assembly on metal surfaces.

We computationally study the reaction mechanisms of halogen-based covalent self-assembly, a major route for synthesizing molecular nanostructures and nanographenes on surfaces. Focusing on biphenyl as a small model system, we describe the dehalogenation, recombination, and diffusion processes. The kinetics of the different processes are also investigated, in particular how diffusion and coupling barriers affect recombination rates. Trends across the periodic table are derived from three commonly used close-packed (111) surfaces (Cu, Ag, and Au) and two halogens (Br and I). We show that the halogen atoms can poison the surface, thus hindering long-range ordering of the self-assembled structures. Finally, we present core-level shifts of the relevant carbon and halogen atoms, to provide reference data for reliably detecting self-assembly without the need for atomic-resolution scanning tunneling microscopy.

Thermally and vibrationally induced tautomerization of single porphycene molecules on a Cu(110) surface.

We report the direct observation of intramolecular hydrogen atom transfer reactions (tautomerization) within a single porphycene molecule on a Cu(110) surface by scanning tunneling microscopy. It is found that the tautomerization can be induced via inelastic electron tunneling at 5 K. By measuring the bias-dependent tautomerization rate of isotope-substituted molecules, we can assign the scanning tunneling microscopy-induced tautomerization to the excitation of specific molecular vibrations. Furthermore, these vibrations appear as characteristic features in the dI/dV spectra measured over individual molecules. The vibrational modes that are associated with the tautomerization are identified by density functional theory calculations. At higher temperatures above ∼75 K, tautomerization is induced thermally and an activation barrier of about 168 meV is determined from an Arrhenius plot.

A Workflow for Identifying Viable Crystal Structures with Partially Occupied Sites Applied to the Solid Electrolyte Cubic Li7La3Zr2O12.

To date, experimental and theoretical works have been unable to uncover the ground-state configuration of the solid electrolyte cubic Li7La3Zr2O12 (c-LLZO). Computational studies rely on an initial low-energy structure as a reference point. Here, we present a methodology for identifying energetically favorable configurations of c-LLZO for a crystallographically predicted structure. We begin by eliminating structures that involve overlapping Li atoms based on nearest neighbor counts. We further reduce the configuration space by eliminating symmetry images from all remaining structures. Then, we perform a machine learning-based energetic ordering of all remaining structures. By considering the geometrical constraints that emerge from this methodology, we determine that a large portion of previously reported structures may not be feasible or stable. The method developed here could be extended to other ion conductors. We provide a database containing all of the generated structures with the aim of improving accuracy and reproducibility in future c-LLZO research.

Kinetically Corrected Monte Carlo-Molecular Dynamics Simulations of Solid Electrolyte Interphase Growth.

We present a kinetic approach to the Monte Carlo-molecular dynamics (MC-MD) method for simulating reactive liquids using nonreactive force fields. A graphical reaction representation allows definition of reactions of arbitrary complexity, including their local solvation environment. Reaction probabilities and molecular dynamics (MD) simulation times are derived from ab initio calculations. Detailed validation is followed by studying the development of the solid electrolyte interphase (SEI) in lithium-ion batteries. We reproduce the experimentally observed two-layered structure on graphite, with an inorganic layer close to the anode and an outer organic layer. This structure develops via a near-shore aggregation mechanism.

Zipping up: cooperativity drives the synthesis of graphene nanoribbons.

We investigate the cooperative effects controlling the synthesis of a graphene nanoribbon on the Au(111) surface starting from an anthracene polymer using density functional calculations including van der Waals interactions. We focus on the high-temperature cyclodehydrogenation step of the reaction and find that the reaction proceeds by simultaneously transferring two H-atoms from the anthracene units to the Au surface, leaving behind a C-C bond in the process. This step is significantly more favorable than the three other potential reaction paths. Moreover, we find that successive dehydrogenations proceed from one end of the polyanthracene and propagate step-by-step through the polymer in a domino-like fashion.

Clean coupling of unfunctionalized porphyrins at surfaces to give highly oriented organometallic oligomers.

The direct coupling of complex, functional organic molecules at a surface is one of the outstanding challenges in the road map to future molecular devices. Equally demanding is to meet this challenge without recourse to additional functionalization of the molecular building blocks and via clean surface reactions that leave no surface contamination. Here, we demonstrate the directional coupling of unfunctionalized porphyrin molecules--large aromatic multifunctional building blocks--on a single crystal copper surface, which generates highly oriented one-dimensional organometallic macromolecular nanostructures (wires) in a reaction which generates gaseous hydrogen as the only byproduct. In situ scanning tunneling microscopy and temperature programmed desorption, supported by theoretical modeling, reveal that the process is driven by C-H bond scission and the incorporation of copper atoms in between the organic components to form a very stable organocopper oligomer comprising organometallic edge-to-edge porphyrin-Cu-porphyrin connections on the surface that are unprecedented in solution chemistry. The hydrogen generated during the reaction leaves the surface and, therefore, produces no surface contamination. A remarkable feature of the wires is their stability at high temperatures (up to 670 K) and their preference for 1D growth along a prescribed crystallographic direction of the surface. The on-surface formation of directional organometallic wires that link highly functional porphyrin cores via direct C-Cu-C bonds in a single-step synthesis is a new development in surface-based molecular systems and provides a versatile approach to create functional organic nanostructures at surfaces.

Sensitivity analysis and uncertainty calculation for dispersion corrected density functional theory.

The precision of binding energies and distances computed with dispersion-corrected density functional theory (DFT-D) is investigated by propagation of uncertainties, yielding relative uncertainties of several percent. Sensitivity analysis is used to calculate the geometry-dependent relative importance of each input parameter for the dispersion correction. While DFT-Ds are exact at asymptotically large distances, their damping functions are shown to play a significant role in binding geometries. This is demonstrated in detail for the interlayer binding of graphite. The techniques presented allow practitioners to quickly compute error bars and to get an a posteriori estimate about the transferability of their results. They can also aid the development of future dispersion corrections.

Rigid dimers formed through strong interdigitated H-bonds yield compact 1D supramolecular helical polymers.

Hierarchical self-assembly of small abiotic molecular modules interacting through noncovalent forces is increasingly being used to generate functional structures and materials for electronic, catalytic, and biomedical applications. The greatest control over the geometry in H-bond supramolecular architectures, especially in H-bonded supramolecular polymers, can be achieved by using conformationally rigid molecular modules undergoing self-assembly through strong H-bonds. Their binding strength depends on the multiplicity of the H-bonds, the nature of donor/acceptor pairs and their secondary attractive/repulsive interactions. Here a functionalized molecular module is described, which is capable of self-associating through self-complementary H-bonding patterns comprising four strong and two medium-strength H-bonds to form dimers. The self-association of these phenylpyrimidine-based dimers through directional H-bonding between two lateral pyridin-2(1H)-one units of neighboring molecules allows the formation of highly compact 1D supramolecular polymers by self-assembly on graphite. A concentration-dependent study by scanning tunneling microscopy at the solid-liquid interface, corroborated by dispersion-corrected density functional studies, reveals the controlled generation of either linear supramolecular 2D arrays, or long helical supramolecular polymers with a high shape persistence.

High-resolution molecular orbital imaging using a p-wave STM tip.

Individual pentacene and naphthalocyanine molecules adsorbed on a bilayer of NaCl grown on Cu(111) were investigated by means of scanning tunneling microscopy using CO-functionalized tips. The images of the frontier molecular orbitals show an increased lateral resolution compared with those of the bare tip and reflect the modulus squared of the lateral gradient of the wave functions. The contrast is explained by tunneling through the p-wave orbitals of the CO molecule. Comparison with calculations using a Tersoff-Hamann approach, including s- and p-wave tip states, demonstrates the significant contribution of p-wave tip states.

Dynamics of single-molecule force-ramp experiments: The role of fluctuations.

In the force-ramp mode of the atomic force microscope, the force with which a macromolecule is stretched is increased linearly in time by properly controlling the motion of the cantilever through a feedback loop. Using a master equation approach for the coupled cantilever-macromolecule system, we mimic such a feedback loop, to study nonequilibrium effects in the measurements of force-extension curves and fluctuations. In particular, it is shown that the fluctuations are the same for force-ramp experiments and for the more commonly used constant velocity experiments. Thus the exact same statistics suffice for the explanation of either experiment. Specific results are presented for the stretching of Dextran.

The simplest supramolecular helix.

Diethylamine is the smallest and simplest molecule that features a supramolecular helix as its lowest energy aggregate. Structural studies and large scale sampling simulations show that the helical arrangement is more stable than cyclic structures, which are the dominant species for other small hydrogen bonding molecules.

Breaking bonds in the atomic force microscope: theory and analysis.

A theoretical framework is developed to analyze molecular bond breaking in dynamic force spectroscopy using atomic force microscopy. An analytic expression of the observed bond breaking probability as a function of force is obtained in terms of the relevant physical parameters. The force-ramp mode is discussed in detail, which gives the best framework to extract the relevant physical parameters such as the potential depth and its width, if a set of widely different force-loading rates are used. We also show that the commonly used Ritchie-Evans model is incomplete and that it is only applicable for forces well below the maximum permitted by the potential. Statistical complications arising from the use of constant velocity experiments are discussed in detail.

The solid state structures of the high and low temperature phases of dimethylcadmium.

The solid state structure of dimethylcadmium, a classic organometallic compound with a long history, has remained elusive for almost a century. X-ray crystallography and density functional theory reveal similar phase behaviour as in dimethylzinc. The high temperature tetragonal phase, α-Me2Cd, exhibits two-dimensional disorder, while the low temperature monoclinic phase, ß-Me2Cd, is ordered. Both phases contain linearly coordinated cadmium atoms. While the methyl groups are staggered in the α-phase, they are eclipsed in the ß-phase.

Nonequilibrium theory of polymer stretching based on the master equation.

We present a model for fast polymer-stretching experiments. We use the master equation and argue that the end-to-end extension of a polymer molecule can be used as a stochastic variable after appropriate coarse graining. The main effect of increasing pulling speed or force loading rate is a marked hysteresis in the force-extension curve as well as an overall shift of the curve to higher forces when compared to the equilibrium curve. This can be understood in terms of the moments of the transition probability in the master equation. An analysis of the fluctuations and relaxation times is also given in the framework of our theory.

Controlling intramolecular hydrogen transfer in a porphycene molecule with single atoms or molecules located nearby.

Although the local environment of a molecule can play an important role in its chemistry, rarely has it been examined experimentally at the level of individual molecules. Here we report the precise control of intramolecular hydrogen-transfer (tautomerization) reactions in single molecules using scanning tunnelling microscopy. By placing, with atomic precision, a copper adatom close to a porphycene molecule, we found that the tautomerization rates could be tuned up and down in a controlled fashion, surprisingly also at rather large separations. Furthermore, we extended our study to molecular assemblies in which even the arrangement of the pyrrolic hydrogen atoms in the neighbouring molecule influences the tautomerization reaction in a given porphycene, with positive and negative cooperativity effects. Our results highlight the importance of controlling the environment of molecules with atomic precision and demonstrate the potential to regulate processes that occur in a single molecule.

Versatile bottom-up construction of diverse macromolecules on a surface observed by scanning tunneling microscopy.

The heterocoupling of organic building blocks to give complex multicomponent macromolecules directly at a surface holds the key to creating advanced molecular devices. While "on-surface" synthesis with prefunctionalized molecules has recently led to specific one- and two- component products, a central challenge is to discover universal connection strategies that are applicable to a wide range of molecules. Here, we show that direct activation of C-H bonds intrinsic to π-functional molecules is a highly generic route for connecting different building blocks on a copper surface. Scanning tunneling microscopy (STM) reveals that covalent π-functional macromolecular heterostructures, displaying diverse compositions, structures and topologies, are created with ease from seven distinct building blocks (including porphyrins, pentacene and perylene). By exploiting differences in C-H bond reactivity in the deposition and heating protocols we also demonstrate controlled synthesis of specific products, such as block copolymers. Further, the symmetry and geometry of the molecules and the surface also play a critical role in determining the outcome of the covalent bond forming reactions. Our "pick-mix-and-link" strategy opens up the capability to generate libraries of multivariate macromolecules directly at a surface, which in conjunction with nanoscale probing techniques could accelerate the discovery of functional interfaces.