A theory of action spectroscopy, i.e., a reaction rate or yield as a function of bias voltage, is presented for single-molecule reactions induced by the inelastic tunneling current with a scanning tunneling microscope. A formula for the reaction yield is derived using the adsorbate resonance model, which provides a versatile tool to analyze vibrationally mediated reactions of single adsorbates on conductive surfaces. This allows us to determine the energy quantum of the excited vibrational mode, the effective broadening of the vibrational density of states (as described by Gaussian or Lorentzian functions), and a prefactor characterizing the elementary process behind the reaction. The underlying approximations are critically discussed. We point out that observation of reaction yields at both bias voltage polarities can provide additional insight into the adsorbate density of states near the Fermi level. As an example, we apply the theory to the case of flip motion of a hydroxyl dimer (OD)(2) on Cu(110) which was experimentally observed by Kumagai et al. [Phys. Rev. B 79, 035423 (2009)]. In combination with density functional theory calculations for the vibrational modes, the vibrational damping due to electron-hole pair generation, and the potential energy landscape for the flip motion, a detailed microscopic picture for the switching process is established. This picture reveals that the predominant mechanism is excitation of the OD stretch modes which couple anharmonically to the low-energy frustrated rotation mode.
Electronic transport mechanisms in molecular junctions are investigated by a combination of first-principles calculations and current-voltage measurements of several well-characterized structures. We study self-assembled layers of alkanethiols grown on Au(111) and form tunnel junctions by contacting the molecular layers with the tip of a conductive force microscope. Measurements done under low-load conditions permit us to obtain reliable tilt-angle and molecular length dependencies of the low-bias conductance through the alkanethiol layers. The observed dependence on tilt-angle is stronger for the longer molecular chains. Our calculations confirm the observed trends and explain them as a result of two mechanisms, namely, a previously proposed intermolecular tunneling enhancement as well as a hitherto overlooked tilt-dependent molecular gate effect.
We describe a first-principles method to calculate scanning tunneling microscopy (STM) images, and compare the results to well-characterized experiments combining STM with atomic force microscopy (AFM). The theory is based on density functional theory with a localized basis set, where the wave functions in the vacuum gap are computed by propagating the localized-basis wave functions into the gap using a real-space grid. Constant-height STM images are computed using Bardeen's approximation method, including averaging over the reciprocal space. We consider copper adatoms and single CO molecules adsorbed on Cu(111), scanned with a single-atom copper tip with and without CO functionalization. The calculated images agree with state-of-the-art experiments, where the atomic structure of the tip apex is determined by AFM. The comparison further allows for detailed interpretation of the STM images.
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.
We present a method used to intuitively interpret the scanning tunneling microscopy (STM) contrast by investigating individual wave functions originating from the substrate and tip side. We use localized basis orbital density functional theory, and propagate the wave functions into the vacuum region at a real-space grid, including averaging over the lateral reciprocal space. Optimization by means of the method of Lagrange multipliers is implemented to perform a unitary transformation of the wave functions in the middle of the vacuum region. The method enables (i) reduction of the number of contributing tip-substrate wave function combinations used in the corresponding transmission matrix, and (ii) to bundle up wave functions with similar symmetry in the lateral plane, so that (iii) an intuitive understanding of the STM contrast can be achieved. The theory is applied to a CO dimer adsorbed on a Cu(1 1 1) surface scanned by a single-atom Cu tip, whose STM image is discussed in detail by the outlined method.
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.
Electron transport through nanotubes is studied theoretically using the Landauer formalism. The studies are carried out for finite metallic nanotubes that bridge two contacts pads. The current is observed to increase stepwise with the applied voltage. Each step corresponds to resonance tunneling including one single-particle eigenstate of the nanotube. Moderate bending of the nanotube results in a shift of the single-particle levels but the overall current remains essentially unaffected. For large bending, however, the pi electron system becomes more disturbed, which introduces backscattering and a marked decrease in the conductivity along the tube. A single carbon vacancy in the nanotube is shown to have very small effect on the conductivity in the center of the metallic band whereas, by increasing the defect concentration the conductivity decreases in the same way as for the strongly bent tubes.
Theory of heat transfer between adsorbate vibrational degrees of freedom and ultrafast laser heated hot electrons including vibrational intermode coupling is applied to calculate two-pulse correlation, laser fluence dependence and time dependence of lateral hopping of CO molecules from a step to terrace site on a stepped Pt (111) surface. The intermode coupling is a key ingredient to describe vibrational heating of the frustrated translation mode responsible for the CO hopping. The calculated results are in good agreement with the experimental results, especially if we scale down the experimentally determined absorbed fluence. It is found that CO hopping is induced by indirect heating of the FT mode by the FR mode with a strong frictional coupling to hot electrons.
For the rational design of single-molecular electronic devices, it is essential to understand environmental effects on the electronic properties of a working molecule. Here we investigate the impact of molecular interactions on the single-molecule conductance by accurately positioning individual molecules on the electrode. To achieve reproducible and precise conductivity measurements, we utilize relatively weak pi-bonding between a phenoxy molecule and a STM-tip to form and cleave one contact to the molecule. The anchoring to the other electrode is kept stable using a chalcogen atom with strong bonding to a Cu(110) substrate. These non-destructive measurements permit us to investigate the variation in single-molecule conductance under different but controlled environmental conditions. Combined with density functional theory calculations, we clarify the role of the electrostatic field in the environmental effect that influences the molecular level alignment.
Mechanical methods for single-molecule control have potential for wide application in nanodevices and machines. Here we demonstrate the operation of a single-molecule switch made functional by the motion of a phenyl ring, analogous to the lever in a conventional toggle switch. The switch can be actuated by dual triggers, either by a voltage pulse or by displacement of the electrode, and electronic manipulation of the ring by chemical substitution enables rational control of the on-state conductance. Owing to its simple mechanics, structural robustness, and chemical accessibility, we propose that phenyl rings are promising components in mechanical molecular devices.
We present first-principles calculations of the dI/dV characteristics of an H-2 molecule sandwiched between Au and Pt electrodes in the presence of electron-phonon interactions. The conductance is found to decrease by a few percentages at threshold voltages corresponding to the excitation energy of longitudinal vibrations of the H-2 molecule. In the case of Pt electrodes, the transverse vibrations can mediate transport through otherwise nontransmitting Pt d channels leading to an increase in the differential conductance even though the hydrogen junction is characterized predominately by a single almost fully open transport channel. In the case of Au, the transverse modes do not affect the dI/dV because the Au d states are too far below the Fermi level. A simple explanation of the first-principles results is given using scattering theory. Finally, we compare and discuss our results in relation to experimental data.
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).
Recent investigations of inelastic electron tunneling spectroscopy (IETS) for alkanethiol self-assembled monolayers (SAMs) are reviewed. Alkanethiol SAMs are usually prepared by immersing a gold substrate into a solution of alkanethiol molecules, and they are very stable, even under ambient conditions. Thus, alkanethiol SAMs have been used as typical molecules for research into molecular electronics. Infrared spectroscopy and electron energy loss spectroscopy (EELS) have frequently been employed to characterize SAMs on the macroscopic scale. For characterization of alkanethiol SAMs on the nanometer scale region, or for single alkanethiol molecules through which electrons actually tunnel, IETS has proven to be an effective method. However, IETS experiments for alkanethiol SAMs employing different methods have shown large differences, i.e., there is a lack of standard data for alkanethiol SAMs with which to understand the IET process or to satisfactorily compare with theoretical investigations. An effective means of acquiring standard data is the formation of a tunneling junction with scanning tunneling microscopy (STM). After explanation of the STM experimental techniques, standard IETS data are presented whereby a contact condition between the tip and SAM is tuned. We have found that many vibrational modes are detected by STM-IETS, as is also the case for EELS. These results are compared with LET spectra measured with different tunneling junctions. In order to precisely investigate which vibrational modes are active in IETS, isotope labeling of alkanethiols with specifically synthesized isotopically substituted molecule has been examined. This method provides unambiguous assignments of IET spectra peaks and site selectivity for alkanethiol SAMs such that all parts of the alkanethiol molecules almost equally contribute to the IET process. The LET process is also discussed based on density functional theory and nonequilibrium Green's function calculations. These results quantitatively reproduce many the experimentally observed features, whereas Fermi's golden rule for JETS qualitatively explains the propensity rule and site selectivity observed in the experiments. However, comparison between experiment and theory reveals a large difference in JETS intensity for the C H stretching mode that originates from the side chains of the alkanethiol molecules. In order to explain this difference, we discuss the importance of an intermolecular tunneling process in the SAM. Application of STM-IETS to identify a hydrogenated alkanethiol molecule inserted into a deuterated alkanethiol SAM matrix is also demonstrated. (C) 2012 Elsevier Ltd. All rights reserved.
The vibration of a molecule adsorbed on a surface contains essential information on the molecule–surface bond, which is important to understand the surface reactions that occur, e.g., in catalytic reactions. Accessing the vibrational energies of a single molecule is possible by combining scanning tunneling microscopy with inelastic electron spectroscopy. However, the tip of a microscope exerts a force on a nearby molecule, and possibly even induces slight structural changes. To study this problem, we have further incorporated atomic force microscopy. The relationship between the exerted forces and vibrational energies is well reproduced by a classical mechanical model. This combined technique opens the possibility to study the atomic-scale interaction of a molecule on a surface with unprecedented precision.The oscillation frequencies of a molecule on a surface are determined by the mass distribution in the molecule and the restoring forces that occur when the molecule bends. The restoring force originates from the atomic-scale interaction within the molecule and with the surface, which plays an essential role in the dynamics and reactivity of the molecule. In 1998, a combination of scanning tunneling microscopy with inelastic tunneling spectroscopy revealed the vibrational frequencies of single molecules adsorbed on a surface. However, the probe tip itself exerts forces on the molecule, changing its oscillation frequencies. Here, we combine atomic force microscopy with inelastic tunneling spectroscopy and measure the influence of the forces exerted by the tip on the lateral vibrational modes of a carbon monoxide molecule on a copper surface. Comparing the experimental data to a mechanical model of the vibrating molecule shows that the bonds within the molecule and with the surface are weakened by the proximity of the tip. This combination of techniques can be applied to analyze complex molecular vibrations and the mechanics of forming and loosening chemical bonds, as well as to study the mechanics of bond breaking in chemical reactions and atomic manipulation.
Vibrationally induced configurational change and nonlinear current-voltage (I-V) characteristics are investigated within the scanning tunneling microscope (STM) junction, including hydroxyl dimers on a Cu(110) surface. H-bonded hydroxyl dimers composed of OH and/or OD have a unique inclined geometry that can be switched back and forth by vibrational excitations via the inelastic electron tunneling process of the STM. The relative occupation change between the high-and low-conductance states as a function of bias voltage critically depends on the isotopic compositions, and thus the I-V characteristics can be modified to exhibit negative differential resistance by H/D substitution. The experimental results of the occupation change and I-V curves are nicely reproduced using a recently proposed analytical model combined with comprehensive density functional calculations for the input parameters (vibrational modes and their emission rates by tunneling electrons, conductance, and relative occupation change of high-and low-conductance states), and they underlines the different roles played by the free and shared O-H(D) stretch modes of the hydroxyl dimers on a Cu(110) surface.
Computationally inexpensive approximations describing electron-phonon scattering in molecular-scale conductors are derived from the nonequilibrium Green's function method. The accuracy is demonstrated with a first-principles calculation on an atomic gold wire. Quantitative agreement between the full nonequilibrium Green's function calculation and the newly derived expressions is obtained while simplifying the computational burden by several orders of magnitude. In addition, analytical models provide intuitive understanding of the conductance including nonequilibrium heating and provide a convenient way of parameterizing the physics. This is exemplified by fitting the expressions to the experimentally observed conductances through both an atomic gold wire and a hydrogen molecule.
We present a method to analyze the results of first-principles based calculations of electronic currents including inelastic electron-phonon effects. This method allows us to determine the electronic and vibrational symmetries in play, and hence to obtain the so-called propensity rules for the studied systems. We show that only a few scattering states—namely those belonging to the most transmitting eigenchannels—need to be considered for a complete description of the electron transport. We apply the method on first-principles calculations of four different systems and obtain the propensity rules in each case.
We study formation and conductance of alkanedithiol junctions using density functional based molecular dynamics. The formation involves straightening of the molecule, migration of thiol end-groups, and pulling out Au atoms. Plateaus are found in the low-bias conductance traces which decrease by 1 order of magnitude when gauche defects are present. We further show that the inelastic electron tunneling spectra depend on the junction geometry. In particular, our simulations suggest ways to identify gauche defects.
A comprehensive analysis of the elementary processes behind the scanning tunneling microscope controlled rotation of C2H2 and C2D2, isotopologues of a single acetylene molecule adsorbed on the Cu(001) surface, is given, with a focus on the isotope effects. With the help of density-functional theory we calculate the vibrational modes of C2H2 and C2D2 on Cu(001) and estimate the anharmonic couplings between them, using a simple strings-on-rods model. The probability of the elementary processes, nonlinear and combination band, is estimated using the Keldysh diagram technique. This allows us to clarify the main peculiarities and the isotope effects of the C2H2 and C2D2 on Cu(001) rotation, discovered in the pioneering work [B. C. Stipe et al., Phys. Rev. Lett. 81, 1263 (1998)], which have not been previously understood.
We study the elementary processes behind one of the pioneering works on scanning tunneling microscope controlled reactions of single molecules [Stipe et al., Phys. Rev. Lett. 81, 1263 (1998)]. Using the Keldysh-Green function approach for the vibrational generation rate in combination with density functional theory calculations to obtain realistic parameters we reproduce the experimental rotation rate of an acetylene molecule on a Cu(100) surface as a function of bias voltage and tunneling current. This combined approach allows us to identify the reaction coordinate mode of the acetylene rotation and its anharmonic coupling with the C-H stretch mode. We show that three different elementary processes, the excitation of C-H stretch, the overtone ladder climbing of the hindered rotational mode, and the combination band excitation together explain the rotation of the acetylene molecule on Cu(100).