To describe vibrationally mediated configuration changes of adsorbates on surfaces we have developed a theory to calculate both reaction rates and pathways. The method uses the T-matrix to describe excitations of vibrational states by the electrons of the substrate, adsorbate, and tunneling electrons from a scanning tunneling probe. In addition to reaction rates, the theory also provides the reaction pathways by going beyond the harmonic approximation and using the full potential energy surface of the adsorbate which contains local minima corresponding to the adsorbates different configurations. To describe the theory, we reproduce the experimental results in [T. Kumagai et al., Phys. Rev. B 79, 035423 (2009)], where the hydrogen/deuterium atom of an adsorbed hydroxyl (OH/OD) exhibits back and forth flipping between two equivalent configurations on a Cu(110) surface at T = 6 K. We estimate the potential energy surface and the reaction barrier, similar to 160 meV, from DFT calculations. The calculated flipping processes arise from (i) at low bias, tunneling of the hydrogen through the barrier, (ii) intermediate bias, tunneling electrons excite the vibrations increasing the reaction rate although over the barrier processes are rare, and (iii) higher bias, overtone excitations increase the reaction rate further.

We present a method capable of calculating elastic scanning tunneling microscopy (STM) currents from localized atomic orbital density functional theory (DFT). To overcome the poor accuracy of the localized orbital description of the wave functions far away from the atoms, we propagate the wave functions, using the total DFT potential. From the propagated wave functions, the Bardeen's perturbative approach provides the tunneling current. To illustrate the method we investigate carbon monoxide adsorbed on a Cu(111) surface and recover the depression/protrusion observed experimentally with normal/CO-functionalized STM tips. The theory furthermore allows us to discuss the significance of s- and p-wave tips.

Achieving a high intensity in inelastic scanning tunneling spectroscopy (IETS) is important for precise measurements. The intensity of the IETS signal can vary by up to a factor of 3 for various tips without an apparent reason accessible by scanning tunneling microscopy (STM) alone. Here, we show that combining STM and IETS with atomic force microscopy enables carbon monoxide front-atom identification, revealing that high IETS intensities for CO/Cu(111) are obtained for single-atom tips, while the intensity drops sharply for multiatom tips. Adsorption of the CO molecule on a Cu adatom [CO/Cu/Cu(111)] such that the molecule is elevated over the substrate strongly diminishes the tip dependence of IETS intensity, showing that an elevated position channels most of the tunneling current through the CO molecule even for multiatom tips, while a large fraction of the tunneling current bypasses the CO molecule in the case of CO/Cu(111).