Thermal Diffusion Processes in Metal-Tip-Surface Interactions: Contact Formation and Adatom Mobility

Kinetic Monte Carlo modeling of oxide thin film growth

In spite of the increasing interest in and application of ultrathin film oxides in commercial devices, the understanding of the mechanisms that control the growth of these films at the atomic scale remains limited and scarce. This limited understanding prevents the rational design of novel solutions based on precise control of the structure and properties of ultrathin films. Such a limited understanding stems in no minor part from the fact that most of the available modeling methods are unable to access and robustly sample the nanosecond to second timescales required to simulate both atomic deposition and surface reorganization at ultrathin films. To contribute to this knowledge gap, here we have combined molecular dynamics and adaptive kinetic Monte Carlo simulations to study the deposition and growth of oxide materials over an extended timescale of up to ∼0.5 ms. In our pilot studies, we have examined the growth of binary oxide thin films on oxide substrates. We have investigated three scenarios: (i) the lattice parameter of both the substrate and thin film are identical, (ii) the lattice parameter of the thin film is smaller than the substrate, and (iii) the lattice parameter is greater than the substrate. Our calculations allow for the diffusion of ions between deposition events and the identification of growth mechanisms in oxide thin films. We make a detailed comparison with previous calculations. Our results are in good agreement with the available experimental results and demonstrate important limitations in former calculations, which fail to sample phase space correctly at the temperatures of interest (typically 300–1000 K) with self-evident limitations for the representative modeling of thin films growth. We believe that the present pilot study and proposed combined methodology open up for extended computational support in the understanding and design of ultrathin film growth conditions tailored to specific applications.

Zinc(II) tetraphenylporphyrin on Au(111) investigated by scanning tunnelling microscopy and photoemission spectroscopy measurements

Porphyrins are a versatile class of molecules, which have attracted attention over the years due to their electronic, optical and biological properties. Self-assembled monolayers of porphyrins were widely studied on metal surfaces in order to understand the supramolecular organization of these molecules, which is a crucial step towards the development of devices starting from the bottom-up approach. This perspective could lead to tailor the interfacial properties of the surface, depending on the specific interaction between the molecular assembly and the metal surface. In this study, we revisit the investigation of the assembly of zinc-tetraphenylporphyrins on Au(111) in order to explore the adsorption of the molecular network on the noble metal substrate. The combined analysis of scanning tunneling microscopy (STM) imaging and core levels photoemission spectroscopy measurements support a peculiar arrangement of the ZnTPP molecular network, with Zn atoms occupying the bridge sites of the Au surface atoms. Furthermore, we prove that, at few-layers coverage, the interaction between the deposited layers allows a relevant molecular mobility of the adlayer, as observed by STM and supported by core levels photoemission analysis.

An efficient algorithm for finding the minimum energy path for cation migration in ionic materials

The Nudged Elastic Band (NEB) is an established method for finding minimum-energy paths and energy barriers of ion migration in materials, but has been hampered in its general application by its significant computational expense when coupled with density functional theory (DFT) calculations. Typically, an NEB calculation is initialized from a linear interpolation of successive intermediate structures (also known as images) between known initial and final states. However, the linear interpolation introduces two problems: (1) slow convergence of the calculation, particularly in cases where the final path exhibits notable curvature; (2) divergence of the NEB calculations if any intermediate image comes too close to a non-diffusing species, causing instabilities in the ensuing calculation. In this work, we propose a new scheme to accelerate NEB calculations through an improved path initialization and associated energy estimation workflow. We demonstrate that for cation migration in an ionic framework, initializing the diffusion path as the minimum energy path through a static potential built upon the DFT charge density reproduces the true NEB path within a 0.2 Å deviation and yields up to a 25% improvement in typical NEB runtimes. Furthermore, we find that the locally relaxed energy barrier derived from this initialization yields a good approximation of the NEB barrier, with errors within 20 meV of the true NEB value, while reducing computational expense by up to a factor of 5. Finally, and of critical importance for the automation of migration path calculations in high-throughput studies, we find that the new approach significantly enhances the stability of the calculation by avoiding unphysical image initialization. Our algorithm promises to enable efficient calculations of diffusion pathways, resolving a long-standing obstacle to the computational screening of intercalation compounds for Li-ion and multivalent batteries.

Diffusive versus Displacive Contact Plasticity of Nanoscale Asperities: Temperature- and Velocity-Dependent Strongest Size

We predict a strongest size for the contact strength when asperity radii of curvature decrease below 10 nm. The reason for such strongest size is found to be correlated with the competition between the dislocation plasticity and surface diffusional plasticity. The essential role of temperature is calculated and illustrated in a comprehensive asperity size-strength-temperature map taking into account the effect of contact velocity. Such a map should be essential for various phenomena related to nanoscale contacts such as nanowire cold welding, self-assembly of nanoparticles and adhesive nanopillar arrays, as well as the electrical, thermal, and mechanical properties of macroscopic interfaces.

Implementation of atomically defined field ion microscopy tips in scanning probe microscopy

The field ion microscope (FIM) can be used to characterize the atomic configuration of the apices of sharp tips. These tips are well suited for scanning probe microscope (SPM) use since they predetermine the SPM resolution and the electronic structure for spectroscopy. A protocol is proposed for preserving the atomic structure of the tip apex from etching due to gas impurities during the period of transfer from the FIM to the SPM, and estimations are made regarding the time limitations of such an experiment due to contamination with ultra-high vacuum rest gases. While avoiding any current setpoint overshoot to preserve the tip integrity, we present results from approaches of atomically defined tungsten tips to the tunneling regime with Au(111), HOPG (highly oriented pyrolytic graphite) and Si(111) surfaces at room temperature. We conclude from these experiments that adatom mobility and physisorbed gas on the sample surface limit the choice of surfaces for which the tip integrity is preserved in tunneling experiments at room temperature. The atomic structure of FIM tip apices is unchanged only after tunneling to the highly reactive Si(111) surface.

Single-layer metal-on-metal islands driven by strong time-dependent forces

Nonlinear transport properties of single-layer metal-on-metal islands driven with strong static and time-dependent forces are studied. We apply a semiempirical lattice model and use master-equation and kinetic Monte Carlo simulation methods to compute observables such as the velocity and the diffusion coefficient. Two types of time-dependent driving are considered: a pulsed rotated field and an alternating field with a zero net force (electrophoretic ratchet). Small islands up to 12 atoms were studied in detail with the master-equation method and larger ones with simulations. Results are presented mainly for a parametrization of Cu on Cu(001) surface, which has been the main system of interest in several previous studies. The main results are that the pulsed field can increase the current in both diagonal and axis direction when compared to static field, and there exists a current inversion in the electrophoretic ratchet. Both of these phenomena are a consequence of the coupling of the internal dynamics of the island with its transport. In addition to the previously discovered "magic size"effect for islands in equilibrium, a strong odd-even effect was found for islands driven far out of equilibrium. Master-equation computations revealed nonmonotonous behavior for the leading relaxation constant and effective Arrhenius parameters. Using cycle optimization methods, typical island transport mechanisms are identified for small islands.

Fabrication of highly stable configurable metal quantum point contacts

Metal quantum point contacts (MQPCs), with dimensions comparable to the de Broglie wavelength of conducting electrons, reveal ballistic transport of electrons and quantized conductance in units of G0=2e(2)/h. While these contacts hold great promise for applications such as coherent controlled devices and atomic switches, their realization is mainly based on the scanning tunneling microscope (STM) and mechanically controlled break junction (MCBJ), which cannot be integrated into electronic circuits. MQPCs produced by these techniques have also limited stability at room temperature. Here we report on a new method to form MQPCs with quantized conductance values in the range of 1-4G0. The contacts appear to be stable at room temperature for hours and can be deterministically switched between conductance values, or reform in case they break, using voltage pulses. The method enables us to integrate MQPCs within nanoscale circuits to fully harness their unique advantages.

A flexible nudged elastic band program for optimization of minimum energy pathways using ab initio electronic structure methods

A driver program for carrying out nudged elastic band optimizations of minimum energy reaction pathways is described. This approach allows for the determination of minimum energy pathways using only energies and gradient information. The driver code has been interfaced with the GAUSSIAN 98 program. Applications to two isomerization reactions and to a cluster model for H(2) desorption from the Si(100)-2 x 1 surface are presented.

Study on tip-substrate interactions by STM and APFIM

Processes occurring at the interface of two materials coming in contact, separating or moving with respect to each other have been studied with the scanning tunnelling microscope (STM) and atom-probe (AP) field ion microscopy (APFIM). STM probe tips have been first characterised by field ion microscopy (FIM), brought into well-defined contact in the STM and afterwards inspected by time-of-flight AP. The results from mechanical contact and indentation experiments, showing material transfer and neck formation, are in reasonable good agreement with computer-based simulations on metal tip-surface interactions.

Theoretical studies of atomic-scale processes relevant to crystal growth

The study of adsorption, diffusion, island formation, and interlayer transport of atoms on a growing surface has been an active field in the past decade, because of both experimental and theoretical advances. Experiments can give detailed images of patterns formed on growing surfaces. An important challenge to the theoretical studies is the identification of dynamical processes controlling the pattern formation and overall surface morphology. This can be achieved by accurate modeling of the atomic interactions, a thorough search for active atomic-scale processes, and simulation of the growth on the experimental timescale to allow for detailed comparison with the experimental measurements. An overview of some of the theoretical methodology used in these studies and results obtained for one of the most extensively studied systems, Pt(111), is given here. A remarkable richness of phenomena has emerged from these studies, where apparently small effects can shift the balance between competing molecular processes and thereby change the morphology of a growing surface. The article concludes with a discussion of possible future directions in this research area.