Editing Scanning tunneling microscope (section)

== Principle of operation ==
Quantum tunneling of electrons is a functioning concept of STM that arises from [[quantum mechanics]]. Classically, a particle hitting an impenetrable barrier will not pass through. If the barrier is described by a potential acting along ''z'' direction, in which an electron of mass ''m''<sub>e</sub> acquires the potential energy ''U''(''z''), the electron's trajectory will be deterministic and such that the sum ''E'' of its kinetic and potential energies is at all times conserved:
: <math>E = \frac{p^2}{2m_\text{e}} + U(z).</math>

The electron will have a defined, non-zero momentum ''p'' only in regions where the initial energy ''E'' is greater than ''U''(''z''). In quantum physics, however, the electron can pass through classically forbidden regions. This is referred to as [[Quantum tunnelling|tunneling]].<ref name="Chen" />

=== Rectangular barrier model ===
[[File:Scanning tunneling microscope - rectangular potential barrier model.svg|thumb|300x300px|The real and imaginary parts of the wave function in a rectangular potential barrier model of the scanning tunneling microscope]]
The simplest model of tunneling between the sample and the tip of a scanning tunneling microscope is that of a [[rectangular potential barrier]].<ref name="Lounis">{{cite arXiv |last=Lounis |first=Samir |name-list-style = vanc |date=2014-04-03 |title=Theory of Scanning Tunneling Microscopy |class=cond-mat.mes-hall |eprint=1404.0961}}</ref><ref name="Chen" /> An electron of energy ''E'' is incident upon an energy barrier of height ''U'', in the region of space of width ''w''. An electron's behavior in the presence of a potential ''U''(''z''), assuming one-dimensional case, is described by [[wave function]]s <math>\psi(z)</math> that satisfy [[Schrödinger's equation]]
: <math>-\frac{\hbar^2}{2m_\text{e}} \frac{\partial^2\psi(z)}{\partial z^2} + U(z)\,\psi(z) = E\,\psi(z),</math>

where ''ħ'' is the [[reduced Planck constant]], ''z'' is the position, and ''m''<sub>e</sub> is the [[electron mass]]. In the zero-potential regions on two sides of the barrier, the wave function takes on the forms

: <math>\psi_L(z) = e^{ikz} + r\,e^{-ikz}</math> for ''z'' < 0,
: <math>\psi_R(z) = t\,e^{ikz}</math> for ''z'' > ''w'',

where <math>k = \tfrac{1}{\hbar} \sqrt{2m_\text{e}E}</math>. Inside the barrier, where ''E'' < ''U'', the wave function is a superposition of two terms, each decaying from one side of the barrier:
: <math>\psi_B (z) = \xi e^{-\kappa z} + \zeta e^{\kappa z}</math> for 0 < ''z'' < ''w'',
where <math>\kappa = \tfrac{1}{\hbar} \sqrt{2m_\text{e}(U - E)}</math>.

The coefficients ''r'' and ''t'' provide measure of how much of the incident electron's wave is reflected or transmitted through the barrier. Namely, of the whole impinging particle current <math>j_i = \hbar k/m_\text{e}</math> only <math>j_t = |t|^2\, j_i</math> is transmitted, as can be seen from the [[probability current]] expression
: <math>j_t = -i \frac{\hbar}{2m_\text{e}} \left\{\psi_R^* \frac{\partial}{\partial z}\psi_R - \psi_R \frac{\partial}{\partial z}\psi_R^*\right\},</math>

which evaluates to <math>j_t = \tfrac{\hbar k}{m_\text{e}} \vert t \vert^2</math>. The transmission coefficient is obtained from the continuity condition on the three parts of the wave function and their derivatives at ''z'' = 0 and ''z'' = ''w'' (detailed derivation is in the article [[Rectangular potential barrier#Analysis of the obtained expressions|Rectangular potential barrier]]). This gives <math>|t|^2 = \big[1 + \tfrac{1}{4} \varepsilon^{-1}(1 - \varepsilon)^{-1} \sinh^2\kappa w\big]^{-1},</math> where <math>\varepsilon = E/U</math>. The expression can be further simplified, as follows:

In STM experiments, typical barrier height is of the order of the material's surface [[work function]] ''W'', which for most metals has a value between 4 and 6&nbsp;eV.<ref name="Lounis" /> The [[work function]] is the minimum energy needed to bring an electron from an occupied level, the highest of which is the [[Fermi level]] (for metals at ''T'' = 0&nbsp;K), to [[vacuum level]]. The electrons can tunnel between two metals only from occupied states on one side into the unoccupied states of the other side of the barrier. Without bias, Fermi energies are flush, and there is no tunneling. Bias shifts electron energies in one of the electrodes higher, and those electrons that have no match at the same energy on the other side will tunnel. In experiments, bias voltages of a fraction of 1&nbsp;V are used, so <math>\kappa</math> is of the order of 10 to 12&nbsp;nm<sup>−1</sup>, while ''w'' is a few tenths of a nanometre. The barrier is strongly attenuating. The expression for the transmission probability reduces to <math>|t|^2 = 16\,\varepsilon(1 - \varepsilon)\,e^{-2\kappa w}.</math> The tunneling current from a single level is therefore<ref name="Lounis" />
: <math>j_t = \left[\frac{4k\kappa}{k^2 + \kappa^2}\right]^2 \, \frac{\hbar k}{m_\text{e}}\,e^{-2\kappa w},</math>

where both wave vectors depend on the level's energy ''E'', <math>k = \tfrac{1}{\hbar} \sqrt{2m_\text{e}E},</math> and <math>\kappa = \tfrac{1}{\hbar}\sqrt{2m_\text{e}(U - E)}.</math>

Tunneling current is exponentially dependent on the separation of the sample and the tip, typically reducing by an order of magnitude when the separation is increased by 1&nbsp;Å (0.1&nbsp;nm).<ref name="Chen" /> Because of this, even when tunneling occurs from a non-ideally sharp tip, the dominant contribution to the current is from its most protruding atom or orbital.<ref name="Lounis" />

=== Tunneling between two conductors ===
[[File:Scanning tunneling microscope - tunneling - Density of states.svg|thumb|300x300px|Negative sample bias ''V'' raises its electronic levels by ''e⋅V''. Only electrons that populate states between the Fermi levels of the sample and the tip are allowed to tunnel.]]
As a result of the restriction that the tunneling from an occupied energy level on one side of the barrier requires an empty level of the same energy on the other side of the barrier, tunneling occurs mainly with electrons near the Fermi level. The tunneling current can be related to the density of available or filled states in the sample. The current due to an applied voltage ''V'' (assume tunneling occurs from the sample to the tip) depends on two factors: 1) the number of electrons between the Fermi level ''E''<sub>F</sub> and ''E''<sub>F</sub>&nbsp;−&nbsp;''eV'' in the sample, and 2) the number among them which have corresponding free states to tunnel into on the other side of the barrier at the tip.<ref name="Chen" /> The higher the density of available states in the tunneling region the greater the tunneling current. By convention, a positive ''V'' means that electrons in the tip tunnel into empty states in the sample; for a negative bias, electrons tunnel out of occupied states in the sample into the tip.<ref name="Chen" />

For small biases and temperatures near absolute zero, the number of electrons in a given volume (the electron concentration) that are available for tunneling is the product of the density of the electronic states ''ρ''(''E''<sub>F</sub>) and the energy interval between the two Fermi levels, ''eV''.<ref name="Chen" /> Half of these electrons will be travelling away from the barrier. The other half will represent the [[Electric current#Drift speed|electric current]] impinging on the barrier, which is given by the product of the electron concentration, charge, and velocity ''v'' (''I''<sub>i</sub>&nbsp;=&nbsp;''nev''),<ref name="Chen" />
: <math>I_i = \tfrac{1}{2} e^2 v\,\rho(E_\text{F})\,V.</math>

The tunneling electric current will be a small fraction of the impinging current. The proportion is determined by the transmission probability ''T'',<ref name="Chen" /> so\
: <math>I_t = \tfrac{1}{2} e^2v\,\rho(E_\text{F})\,V\,T.</math>

In the simplest model of a rectangular potential barrier the transmission probability coefficient ''T'' equals |''t''|<sup>2</sup>.

=== Bardeen's formalism ===
[[File:Scanning tunneling microscope - tip, barrier and sample wave functions.svg|thumb|300x300px|Tip, barrier and sample wave functions in a model of the scanning tunneling microscope. Barrier width is ''w''. Tip bias is ''V''. Surface work functions are ''ϕ''.]]

A model that is based on more realistic wave functions for the two electrodes was devised by [[John Bardeen]] in a study of the [[metal–insulator–metal]] junction.<ref name="Bardeen">{{cite journal |vauthors = Bardeen J |year=1961 |title=Tunneling from a many particle point of view |journal=Phys. Rev. Lett. |volume=6 |issue=2| pages=57–59 |bibcode=1961PhRvL...6...57B |doi=10.1103/PhysRevLett.6.57}}</ref> His model takes two separate orthonormal sets of wave functions for the two electrodes and examines their time evolution as the systems are put close together.<ref name="Chen" /><ref name="Lounis" /> Bardeen's novel method, ingenious in itself,<ref name="Chen" /> solves a time-dependent perturbative problem in which the perturbation emerges from the interaction of the two subsystems rather than an external potential of the standard [[Perturbation theory (quantum mechanics)|Rayleigh–Schrödinger perturbation theory]].

Each of the wave functions for the electrons of the sample (S) and the tip (T) decay into the vacuum after hitting the surface potential barrier, roughly of the size of the surface work function. The wave functions are the solutions of two separate Schrödinger's equations for electrons in potentials ''U''<sub>S</sub> and ''U''<sub>T</sub>. When the time dependence of the states of known energies <math>E^\text{S}_\mu</math> and <math>E^\text{T}_\nu</math> is factored out, the wave functions have the following general form
: <math>\psi^\text{S}_\mu(t) = \psi^\text{S}_\mu \exp\left(-\frac{i}{\hbar} E^\text{S}_\mu t\right),</math>
: <math>\psi^\text{T}_\nu(t) = \psi^\text{T}_\nu \exp\left(-\frac{i}{\hbar} E^\text{T}_\nu t\right).</math>

If the two systems are put closer together, but are still separated by a thin vacuum region, the potential acting on an electron in the combined system is ''U''<sub>T</sub> + ''U''<sub>S</sub>. Here, each of the potentials is spatially limited to its own side of the barrier. Only because the tail of a wave function of one electrode is in the range of the potential of the other, there is a finite probability for any state to evolve over time into the states of the other electrode.<ref name="Chen" /> The future of the sample's state ''μ'' can be written as a linear combination with time-dependent coefficients of <math>\psi^\text{S}_\mu(t)</math> and all <math>\psi^\text{T}_\nu(t)</math>:
: <math id="ansatz">\psi(t) = \psi^\text{S}_\mu(t) + \sum_\nu c_\nu(t) \psi^\text{T}_\nu(t)</math>
with the initial condition <math>c_\nu(0) = 0</math>.<ref name="Chen" /> When the new wave function is inserted into the Schrödinger's equation for the potential ''U''<sub>T</sub> + ''U''<sub>S</sub>, the obtained equation is projected onto each separate <math>\psi^\text{T}_\nu </math> (that is, the equation is multiplied by a <math>{\psi^\text{T}_\nu}^* </math> and integrated over the whole volume) to single out the coefficients <math>c_\nu.</math> All <math>\psi^\text{S}_\mu</math> are taken to be ''nearly orthogonal'' to all <math>\psi^\text{T}_\nu</math> (their overlap is a small fraction of the total wave functions), and only first-order quantities retained. Consequently, the time evolution of the coefficients is given by
: <math>\frac{\mathrm{d}}{\mathrm{d}t} c_\nu(t) = -\frac{i}{\hbar} \int \psi^\text{S}_\mu \,U_\text{T}\, {\psi^\text{T}_\nu}^* \,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z\,\exp\left[-\frac{i}{\hbar} (E^\text{S}_\mu - E^\text{T}_\nu) t\right].</math>

Because the potential ''U''<sub>T</sub> is zero at the distance of a few atomic diameters away from the surface of the electrode, the integration over ''z'' can be done from a point ''z''<sub>0</sub> somewhere inside the barrier and into the volume of the tip (''z''&nbsp;>&nbsp;''z''<sub>0</sub>).

If the tunneling matrix element is defined as
: <math>M_{\mu\nu} = \int_{z > z_0} \psi^\text{S}_\mu \,U_\text{T}\, {\psi^\text{T}_\nu}^* \,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z,</math>
the probability of the sample's state ''μ'' evolving in time ''t'' into the state of the tip ''ν'' is
: <math>|c_\nu (t)|^2 = |M_{\mu\nu}|^2 \frac{4 \sin^2\big[\tfrac{1}{2\hbar}(E^\text{S}_\mu - E^\text{T}_\nu )t\big]}{(E^\text{S}_\mu - E^\text{T}_\nu)^2}.</math>

In a system with many electrons impinging on the barrier, this probability will give the proportion of those that successfully tunnel. If at a time ''t'' this fraction was <math>|c_\nu (t)|^2,</math> at a later time ''t''&nbsp;+&nbsp;d''t'' the total fraction of <math>|c_\nu(t + \mathrm{d}t)|^2</math> would have tunneled. The ''current'' of tunneling electrons at each instance is therefore proportional to <math>|c_\nu(t + \mathrm{d}t)|^2 - |c_\nu(t)|^2</math> divided by <math>\mathrm{d}t,</math> which is the time derivative of <math>|c_\nu(t)|^2,</math><ref name="Lounis" />
: <math>\Gamma_{\mu \to \nu}\ \overset{\text{def}}{=}\ \frac{\mathrm{d}}{\mathrm{d}t} |c_\nu(t)|^2 = \frac{2\pi}{\hbar} |M_{\mu\nu}|^2\frac{\sin\big[(E^\text{S}_\mu - E^\text{T}_\nu) \tfrac{t}{\hbar}\big]}{\pi(E^\text{S}_\mu - E^\text{T}_\nu)}.</math>

The time scale of the measurement in STM is many orders of magnitude larger than the typical [[femtosecond]] time scale of electron processes in materials, and <math>t/\hbar</math> is large. The fraction part of the formula is a fast-oscillating function of <math>(E^\text{S}_\mu - E^\text{T}_\nu)</math> that rapidly decays away from the central peak, where <math>E^\text{S}_\mu = E^\text{T}_\nu</math>. In other words, the most probable tunneling process, by far, is the elastic one, in which the electron's energy is conserved. The fraction, as written above, is a representation of the [[Dirac delta function#Oscillatory integrals|delta function]], so
: <math>\Gamma_{\mu \to \nu} = \frac{2\pi}{\hbar} |M_{\mu\nu}|^2 \delta(E^\text{S}_\mu - E^\text{T}_\nu).</math>

Solid-state systems are commonly described in terms of continuous rather than discrete energy levels. The term <math>\delta(E^\text{S}_\mu - E^\text{T}_\nu)</math> can be thought of as the [[Density of states#Definition|density of states]] of the tip at energy <math>E^\text{S}_\mu,</math> giving
: <math>\Gamma_{\mu \to \nu} = \frac{2\pi}{\hbar} |M_{\mu\nu}|^2 \rho_\text{T}(E^\text{S}_\mu).</math>

The number of energy levels in the sample between the energies <math>\varepsilon</math> and <math>\varepsilon + \mathrm{d}\varepsilon</math> is <math>\rho_\text{S}(\varepsilon)\,\mathrm{d}\varepsilon.</math> When occupied, these levels are spin-degenerate (except in a few special classes of materials) and contain charge <math>2e \cdot \rho_\text{S}(\varepsilon)\,\mathrm{d}\varepsilon</math> of either spin. With the sample biased to voltage <math>V,</math> tunneling can occur only between states whose occupancies, given for each electrode by the [[Fermi–Dirac distribution]] <math>f</math>, are not the same, that is, when either one or the other is occupied, but not both. That will be for all energies <math>\varepsilon</math> for which <math>f(E_\text{F} - eV + \varepsilon) - f(E_\text{F} + \varepsilon)</math> is not zero. For example, an electron will tunnel from energy level <math>E_\text{F} - eV</math> in the sample into energy level <math>E_\text{F}</math> in the tip (<math>\varepsilon = 0</math>), an electron at <math>E_\text{F}</math> in the sample will find unoccupied states in the tip at <math>E_\text{F} + eV</math> (<math>\varepsilon = eV</math>), and so will be for all energies in between. The tunneling current is therefore the sum of little contributions over all these energies of the product of three factors: <math>2e \cdot \rho_\text{S}(E_\text{F} - eV + \varepsilon)\,\mathrm{d}\varepsilon</math> representing available electrons, <math>f(E_\text{F} - eV + \varepsilon) - f(E_\text{F} + \varepsilon)</math> for those that are allowed to tunnel, and the probability factor <math>\Gamma</math> for those that will actually tunnel:
: <math>I_t = \frac{4 \pi e}{\hbar} \int_{-\infty}^{+\infty} [f(E_\text{F} - eV + \varepsilon) - f(E_\text{F} + \varepsilon)] \, \rho_\text{S}(E_\text{F} - eV + \varepsilon) \, \rho_\text{T}(E_\text{F} + \varepsilon) \, |M|^2 \, d \varepsilon.</math>

Typical experiments are run at a liquid-helium temperature (around 4&nbsp;K), at which the Fermi-level cut-off of the electron population is less than a millielectronvolt wide. The allowed energies are only those between the two step-like Fermi levels, and the integral becomes
: <math>I_t = \frac{4 \pi e}{\hbar} \int_0^{eV} \rho_\text{S}(E_\text{F} - eV + \varepsilon) \, \rho_\text{T}(E_\text{F} + \varepsilon) \, |M|^2 \, d \varepsilon.</math>

When the bias is small, it is reasonable to assume that the electron wave functions and, consequently, the tunneling matrix element do not change significantly in the narrow range of energies. Then the tunneling current is simply the convolution of the densities of states of the sample surface and the tip:
: <math>I_t \propto \int_0^{eV} \rho_\text{S}(E_\text{F} - eV + \varepsilon) \, \rho_\text{T}(E_\text{F} + \varepsilon) \, d \varepsilon.</math>

How the tunneling current depends on distance between the two electrodes is contained in the tunneling matrix element
: <math>M_{\mu\nu} = \int_{z > z_0} \psi^\text{S}_\mu \,U_\text{T}\, {\psi^\text{T}_\nu}^* \,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z.</math>

This formula can be transformed so that no explicit dependence on the potential remains. First, the <math>U_\text{T}\, {\psi^\text{T}_\nu}^*</math> part is taken out from the Schrödinger equation for the tip, and the elastic tunneling condition is used so that
: <math>M_{\mu\nu} = \int_{z > z_0} \left({\psi^\text{T}_\nu}^* E_\mu \psi^\text{S}_\mu + \psi^\text{S}_\mu \frac{\hbar^2}{2m} \frac{\partial^2}{\partial z^2}{\psi^\text{T}_\nu}^*\right) \,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z.</math>

Now <math>E_\mu\, {\psi^\text{S}_\mu}</math> is present in the Schrödinger equation for the sample and equals the kinetic plus the potential operator acting on <math>\psi^\text{S}_\mu.</math> However, the potential part containing ''U''<sub>S</sub> is on the tip side of the barrier nearly zero. What remains,
: <math>M_{\mu\nu} = -\frac{\hbar^2}{2m} \int_{z > z_0} \left({\psi^\text{T}_\nu}^* \frac{\partial^2}{\partial z^2}{\psi^\text{S}_\mu} - {\psi^\text{S}_\mu} \frac{\partial^2}{\partial z^2}{\psi^\text{T}_\nu}^*\right) \,\mathrm{d}x\,\mathrm{d}y\,\mathrm{d}z,</math>
can be integrated over ''z'' because the integrand in the parentheses equals <math>\partial_z\left({\psi^\text{T}_\nu}^* \, \partial_z \psi^\text{S}_\mu - {\psi^\text{S}_\mu} \, \partial_z {\psi^\text{T}_\nu}^*\right).</math>

Bardeen's tunneling matrix element is an integral of the wave functions and their gradients over a surface separating the two planar electrodes:
: <math>M_{\mu\nu} = \frac{\hbar^2}{2m} \int_{z = z_0} \left(
 {\psi^\text{S}_\mu} \frac{\partial}{\partial z}{\psi^\text{T}_\nu}^* -
 {\psi^\text{T}_\nu}^* \frac{\partial}{\partial z}{\psi^\text{S}_\mu}
\right) \,\mathrm{d}x\,\mathrm{d}y.</math>

The exponential dependence of the tunneling current on the separation of the electrodes comes from the very wave functions that ''leak'' through the potential step at the surface and exhibit exponential decay into the classically forbidden region outside of the material.

The tunneling matrix elements show appreciable energy dependence, which is such that tunneling from the upper end of the ''eV'' interval is nearly an order of magnitude more likely than tunneling from the states at its bottom. When the sample is biased positively, its unoccupied levels are probed as if the density of states of the tip is concentrated at its Fermi level. Conversely, when the sample is biased negatively, its occupied electronic states are probed, but the spectrum of the electronic states of the tip dominates. In this case it is important that the density of states of the tip is as flat as possible.<ref name="Chen" />

The results identical to Bardeen's can be obtained by considering adiabatic approach of the two electrodes and using the standard time-dependent perturbation theory.<ref name="Lounis" /> This leads to [[Fermi's golden rule]] for the transition probability <math>\Gamma_{\mu \to \nu}</math> in the form given above.

Bardeen's model is for tunneling between two planar electrodes and does not explain scanning tunneling microscope's lateral resolution. Tersoff and Hamann<ref>{{cite journal| vauthors = Tersoff J, Hamann DR |date=1983-06-20 |title=Theory and Application for the Scanning Tunneling Microscope |journal=Physical Review Letters |volume=50 |issue=25 |pages=1998–2001 |doi=10.1103/PhysRevLett.50.1998 |bibcode=1983PhRvL..50.1998T |doi-access=free}}</ref><ref>{{cite journal | vauthors = Tersoff J, Hamann DR | title = Theory of the scanning tunneling microscope | journal = Physical Review B | volume = 31 | issue = 2 | pages = 805–813 | date = January 1985 | pmid = 9935822 | doi = 10.1103/PhysRevB.31.805 | bibcode = 1985PhRvB..31..805T | url = https://link.aps.org/doi/10.1103/PhysRevB.31.805 }}</ref><ref>{{cite journal |last1=Hansma |first1=Paul K. |last2=Tersoff |first2=Jerry |name-list-style = vanc |date=1987-01-15 |title=Scanning tunneling microscopy |url=https://aip.scitation.org/doi/10.1063/1.338189 |journal=Journal of Applied Physics |volume=61 |issue=2 |pages=R1–R24 |doi=10.1063/1.338189 |bibcode=1987JAP....61R...1H |issn=0021-8979}}</ref> used Bardeen's theory and modeled the tip as a structureless geometric point.<ref name="Chen" /> This helped them disentangle the properties of the tip—which are hard to model—from the properties of the sample surface. The main result was that the tunneling current is proportional to the local density of states of the sample at the Fermi level taken at the position of the center of curvature of a spherically symmetric tip (''s''-wave tip model). With such a simplification, their model proved valuable for interpreting images of surface features bigger than a nanometre, even though it predicted atomic-scale corrugations of less than a picometre. These are well below the microscope's detection limit and below the values actually observed in experiments.

In sub-nanometre-resolution experiments, the convolution of the tip and sample surface states will always be important, to the extent of the apparent inversion of the atomic corrugations that may be observed within the same scan. Such effects can only be explained by modeling of the surface and tip electronic states and the ways the two electrodes interact from [[Ab initio quantum chemistry methods|first principles]].