Since STS relies on
tunneling phenomena and measurement of the tunneling current or its
derivative, understanding the expressions for the tunneling current is very important. Using the modified Bardeen transfer Hamiltonian method, which treats tunneling as a
perturbation, the tunneling current (I) is found to be {{NumBlk|| I = \frac{4\pi e}{\hbar}\int_{-\infty}^{\infty}\left[f\left(E_F-eV+\varepsilon\right)-f\left(E_F+\varepsilon\right)\right] \rho_S\left(E_F-eV+\varepsilon\right) \rho_T\left(E_F+\varepsilon\right) \left|M_{\mu\nu}\right|^2 \, d\varepsilon \ ,|}} where f\left(E\right) is the
Fermi distribution function, \rho_s and \rho_T are the
density of states (DOS) in the sample and tip, respectively, and M_{\mu\nu} is the tunneling matrix element between the modified wavefunctions of the tip and the sample surface. The tunneling matrix element, {{NumBlk|| M_{\mu\nu} = -\frac{\hbar^2}{2m} \int_\Sigma\left( \chi_\nu^*\nabla\psi_\mu - \psi_\mu \nabla \chi_\nu^*\right) \cdot \,d{\mathbf {S}}\ ,|}} describes the energy lowering due to the interaction between the two states. Here \psi and \chi are the sample wavefunction modified by the tip potential, and the tip wavefunction modified by sample potential, respectively. For low temperatures and a constant tunneling matrix element, the tunneling current reduces to {{NumBlk|| I \propto\int_0^{eV} \rho_S\left(E_F-eV+\varepsilon\right) \rho_T\left(E_F+\varepsilon\right) \, d\varepsilon\ ,|}} which is a convolution of the DOS of the tip and the sample. Generally, STS experiments attempt to probe the sample DOS, but equation (3) shows that the tip DOS must be known for the measurement to have meaning. Equation (3) implies that {{NumBlk|| \frac{dI}{dV} \propto\rho_S\left(E_F-eV\right)\ ,|}} under the gross assumption that the tip DOS is constant. For these ideal assumptions, the tunneling conductance is directly proportional to the sample DOS. For higher bias voltages, the predictions of simple planar tunneling models using the Wentzel-Kramers Brillouin (WKB) approximation are useful. In the WKB theory, the tunneling current is predicted to be {{NumBlk|| I \propto \int_0^{eV}\rho_S\left(r,E\right)\rho_T\left(r,E-eV\right) T\left(E,eV,r\right)\,dE\ ,|}} where \rho_s and \rho_t are the density of states (DOS) in the sample and tip, respectively. The energy- and bias-dependent electron tunneling transition probability, T, is given by {{NumBlk||T=\exp\left(- \frac{2Z\sqrt{2m}}{\hbar}\sqrt{\frac{\phi_s+\phi_t}{2}+\frac{eV}{2}-E}\right)\ ,|}} where \phi_s and \phi_t are the respective
work functions of the sample and tip and Z is the distance from the sample to the tip. The tip is often regarded to be a single molecule, essentially neglecting further shapes induced effects. This approximation is the Tersoff-Hamann approximation, which suggests the tip to be a single ball-shaped molecule of certain radius. The tunneling current therefore becomes proportional to the local density of states (LDOS). ==Experimental methods==