The dynamical effective potential felt by an electron tunneling in a one-dimensional model STM junction is considered. The electron is coupled inside the barrier to surface plasmons. The corresponding many body Schrödinger equation is solved exactly by means of a matrix method. Results for the electron effective potential and for tunneling times are presented. They are compared to calculations for the same model barrier described within the path integral formalism. As is well known, significant differences from the corresponding static image potential are obtained when tunneling times are shorter than the inverse surface plasmon frequency. However, our results show that path integral calculations underestimate the tunneling traversal time, leading to a larger effective electron potential relative to that obtained by the matrix method.The actual conductance of an STM junction, with a given geometry and tip-sample separation, depends crucially on the potential felt by a tunneling electron [1,2]. A substantial contribution to this potential depends on the extent to which the tunnel junction is polarized by the electron's field. As well known, if tunneling is sufficiently slow (in comparison to the inverse surface plasmon frequency), the electron will feel the "image force" potential. If the tunneling is very rapid, however, the redistribution of the charge around the junction cannot occur sufficiently rapidly to build up the image charge.The dynamical aspects of the image potential have been studied by several authors within different approaches and approximations. For example, the dynamical image potential has been treated (i) classically [3, 4], (ii) within the selfenergy approach by non self-consistent [5, 6], semiclassical self-consistent [7] or fully self-consistent [8] calculations and (iii) by using a path integral technique [9][10][11][12].In this paper, we present a general matrix approach to treat the electron-plasmon coupling inside a tunneling junction. The corresponding many body Schrödinger equation is solved on a real-space mesh. We apply the method to a onedimensional model case and compare the results with those obtained within the path integral formalism for the same model.
The modelsWe consider a model case where an electron tunnels through a static barrier V(r) between two planar metallic surfaces. Describing the metallic electrodes by a local scalar dielectric function, the electron is coupled, within the barrier, only to the surface plasmon (SP) modes. Each mode, of frequency ω ν and coupling matrix elements Γ ν (r), is characterized by a composite index ν = (q, α) where q is a vector parallel to the electrode surfaces and α = (±) represents the even and odd SP modes (i.e. in phase and out of phase oscillation of the charge on the two electrode surfaces).Due to the translation invariance parallel to the electrode surfaces, it is possible to recast the problem into a onedimensional system (only the motion of the electron perpendicular to the electrode surfaces is considered). The corresponding H...