Jumping, Rotating, and Flapping: The Atomic-Scale Motion of Thiophene on Cu(111)

Molecular motion of a nanoscopic moonlander via translations and rotations of triphenylphosphine on graphite

Mass transport at surfaces determines the kinetics of processes such as heterogeneous catalysis and thin-film growth, with the diffusivity being controlled by excitation across a translational barrier. Here, we use neutron spectroscopy to follow the nanoscopic motion of triphenylphosphine (P(C6H5)3 or PPh3) adsorbed on exfoliated graphite. Together with force-field molecular dynamics simulations, we show that the motion is similar to that of a molecular motor, i.e. PPh3 rolls over the surface with an almost negligible activation energy for rotations and motion of the phenyl groups and a comparably small activation energy for translation. While rotations and intramolecular motion dominate up to about 300 K, the molecules follow an additional translational jump-motion across the surface from 350-500 K. The unique behaviour of PPh3 is due to its three-point binding with the surface: Along with van der Waals corrected density functional theory calculations, we illustrate that the adsorption energy of PPh3 increases considerably compared to molecules with flat adsorption geometry, yet the effective diffusion barrier for translational motion increases only slightly. We rationalise these results in terms of molecular symmetry, structure and contact angle, illustrating that the molecular degrees of freedom in larger molecules are intimately connected with the diffusivity.

Chiral and Catalytic Effects of Site-Specific Molecular Adsorption

The changes of properties and preferential interactions based on subtle energetic differences are important characteristics of organic molecules, particularly for their functionalities in biological systems. Only slightly energetically favored interactions are important for the molecular adsorption and bonding to surfaces, which define their properties for further technological applications. Here, prochiral tetracenothiophene molecules are adsorbed on the Cu(111) surface. The chiral adsorption configurations are determined by Scanning Tunneling Microscopy studies and confirmed by first-principles calculations. Remarkably, the selection of the adsorption sites by chemically different moieties of the molecules is dictated by the arrangement of the atoms in the first and second surface layers. Furthermore, we have investigated the thermal effects on the direct desulfurization reaction that occurs under the catalytic activity of the Cu substrate. This reaction leads to a product that is covalently bound to the surface in chiral configurations.

Material properties particularly suited to be measured with helium scattering: selected examples from 2D materials, van der Waals heterostructures, glassy materials, catalytic substrates, topological insulators and superconducting radio frequency materials

Helium Atom Scattering (HAS) and Helium Spin-Echo scattering (HeSE), together helium scattering, are well established, but non-commercial surface science techniques. They are characterised by the beam inertness and very low beam energy (<0.1 eV) which allows essentially all materials and adsorbates, including fragile and/or insulating materials and light adsorbates such as hydrogen to be investigated on the atomic scale. At present there only exist an estimated less than 15 helium and helium spin-echo scattering instruments in total, spread across the world. This means that up till now the techniques have not been readily available for a broad scientific community. Efforts are ongoing to change this by establishing a central helium scattering facility, possibly in connection with a neutron or synchrotron facility. In this context it is important to clarify what information can be obtained from helium scattering that cannot be obtained with other surface science techniques. Here we present a non-exclusive overview of a range of material properties particularly suited to be measured with helium scattering: (i) high precision, direct measurements of bending rigidity and substrate coupling strength of a range of 2D materials and van der Waals heterostructures as a function of temperature, (ii) direct measurements of the electron-phonon coupling constant λ exclusively in the low energy range (<0.1 eV, tuneable) for 2D materials and van der Waals heterostructures (iii) direct measurements of the surface boson peak in glassy materials, (iv) aspects of polymer chain surface dynamics under nano-confinement (v) certain aspects of nanoscale surface topography, (vi) central properties of surface dynamics and surface diffusion of adsorbates (HeSE) and (vii) two specific science case examples - topological insulators and superconducting radio frequency materials, illustrating how combined HAS and HeSE are necessary to understand the properties of quantum materials. The paper finishes with (viii) examples of molecular surface scattering experiments and other atom surface scattering experiments which can be performed using HAS and HeSE instruments.

Quantifying entropic barriers in single-molecule surface diffusion

The quantitative role of entropy in the surface diffusion of molecules with many degrees of freedom is still not well understood. Here, we quantify entropic diffusion barriers as well as attempt frequencies by performing a systematic decomposition of the Arrhenius equation for single oligophenyl molecules of various lengths (two to six phenyl rings and benzene as the reference) on an amorphous silica surface using extensive molecular dynamics simulations. Attempt frequencies evaluated from velocity auto-correlation functions are found close to k_BT/h, the frequency factor of transition state theory. Importantly, we find large positive entropy contributions to the free energy barrier of diffusion up to 55%, increasing with molecular length with 4.1 kJ/mol/phenyl ring. The entropic barrier is about 40%-60% of the entropy of the molecule surface adsorption free energy, revealing that at the transition states, the molecules can liberate a major part of their conformational states, increasing with length. The substantial role of the internal degrees of freedom for the diffusive dynamics is explicitly demonstrated by studying internally constrained, "rigid" version of the molecules. Finally, we discuss also rotational diffusion and the role of surface vibrations. Our results affirm that it is essential for quantitative studies and interpretation of surface diffusion of complex molecules to consider internal entropic effects.

The dynamics of benzene on Cu(111): a combined helium spin echo and dispersion-corrected DFT study into the diffusion of physisorbed aromatics on metal surfaces

We use helium spin-echo spectroscopy (HeSE) to investigate the dynamics of the diffusion of benzene adsorbed on Cu(111). The results of these measurements show that benzene moves on the surface through an activated jump-diffusion process between the adsorption sites on a Bravais lattice. Density Functional Theory (DFT) calculations with van der Waals (vdW) corrections help us understand that the molecule diffuses by jumping through non-degenerate hollow sites. The results of the calculations shed light on the nature of the binding interaction between this prototypical aromatic molecule and the metallic surface. The highly accurate HeSE experimental data provide a quantitatively stringent benchmark for the vdW correction schemes applied to the DFT calculations and we compare the performances of several dispersion interaction schemes.

Mass Transport in Surface Diffusion of van der Waals Bonded Systems: Boosted by Rotations?

Mass transport at a surface is a key factor in heterogeneous catalysis. The rate is determined by excitation across a translational barrier and depends on the energy landscape and the coupling to the thermal bath of the surface. Here we use helium spin-echo spectroscopy to track the microscopic motion of benzene adsorbed on Cu(001) at low coverage (θ ∼ 0.07 ML). Specifically, our combined experimental and computational data determine both the absolute rate and mechanism of the molecular motion. The observed rate is significantly higher by a factor of 3.0 ± 0.1 than is possible in a conventional, point-particle model and can be understood only by including additional molecular (rotational) coordinates. We argue that the effect can be described as an entropic contribution that enhances the population of molecules in the transition state. The process is generally relevant to molecular systems and illustrates the importance of the pre-exponential factor alongside the activation barrier in studies of surface kinetics.

Coupling between diffusion and orientation of pentacene molecules on an organic surface

The realization of efficient organic electronic devices requires the controlled preparation of molecular thin films and heterostructures. As top-down structuring methods such as lithography cannot be applied to van der Waals bound materials, surface diffusion becomes a structure-determining factor that requires microscopic understanding. Scanning probe techniques provide atomic resolution, but are limited to observations of slow movements, and therefore constrained to low temperatures. In contrast, the helium-3 spin-echo (HeSE) technique achieves spatial and time resolution on the nm and ps scale, respectively, thus enabling measurements at elevated temperatures. Here we use HeSE to unveil the intricate motion of pentacene admolecules diffusing on a chemisorbed monolayer of pentacene on Cu(110) that serves as a stable, well-ordered organic model surface. We find that pentacene moves along rails parallel and perpendicular to the surface molecules. The experimental data are explained by admolecule rotation that enables a switching between diffusion directions, which extends our molecular level understanding of diffusion in complex organic systems.

Vibrational lifetimes and friction in adsorbate motion determined from quasi-elastic scattering

The vibrational excitation of molecules adsorbed on a surface is typically probed by spectroscopic techniques such as infrared or Raman spectroscopy. In the present article we demonstrate an alternative method to determine vibrational lifetimes of adsorbate molecules using quasi-elastic helium atom scattering (QHAS). As a probe of diffusive motion of molecules on surfaces QHAS is well established. Here, we demonstrate that QHAS can also be used to probe the vibrational lifetime of a molecule in its adsorption well. Measurements of cyclopentadienyl, C5H5, on Cu(111) allow us to distinguish two substrate phonon modes as well as two molecular vibrational modes, perpendicular and parallel to the surface. We further find that the dephasing of the vibrational motion corresponds to the friction determined in previous diffusion measurements.

Combined diffraction and density functional theory calculations of halogen-bonded cocrystal monolayers

This work describes the combined use of synchrotron X-ray diffraction and density functional theory (DFT) calculations to understand the cocrystal formation or phase separation in 2D monolayers capable of halogen bonding. The solid monolayer structure of 1,4-diiodobenzene (DIB) has been determined by X-ray synchrotron diffraction. The mixing behavior of DIB with 4,4'-bipyridyl (BPY) has also been studied and interestingly is found to phase-separate rather than form a cocrystal, as observed in the bulk. DFT calculations are used to establish the underlying origin of this interesting behavior. The DFT calculations are demonstrated to agree well with the recently proposed monolayer structure for the cocrystal of BPY and 1,4-diiodotetrafluorobenzene (DITFB) (the perfluorinated analogue of DIB), where halogen bonding has also been identified by diffraction. Here we have calculated an estimate of the halogen bond strength by DFT calculations for the DITFB/BPY cocrystal monolayer, which is found to be ∼20 kJ/mol. Computationally, we find that the nonfluorinated DIB and BPY are not expected to form a halogen-bonded cocrystal in a 2D layer; for this pair of species, phase separation of the components is calculated to be lower energy, in good agreement with the diffraction results.