7.7: Quantum Tunneling of Particles through Potential Barriers
- Page ID
- 4532
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)By the end of this section, you will be able to:
- Describe how a quantum particle may tunnel across a potential barrier
- Identify important physical parameters that affect the tunneling probability
- Identify the physical phenomena where quantum tunneling is observed
- Explain how quantum tunneling is utilized in modern technologies
Quantum tunneling is a phenomenon in which particles penetrate a potential energy barrier with a height greater than the total energy of the particles. The phenomenon is interesting and important because it violates the principles of classical mechanics. Quantum tunneling is important in models of the Sun and has a wide range of applications, such as the scanning tunneling microscope and the tunnel diode.
Tunneling and Potential Energy
To illustrate quantum tunneling, consider a ball rolling along a surface with a kinetic energy of 100 J. As the ball rolls, it encounters a hill. The potential energy of the ball placed atop the hill is 10 J. Therefore, the ball (with 100 J of kinetic energy) easily rolls over the hill and continues on. In classical mechanics, the probability that the ball passes over the hill is exactly 1—it makes it over every time. If, however, the height of the hill is increased—a ball placed atop the hill has a potential energy of 200 J—the ball proceeds only part of the way up the hill, stops, and returns in the direction it came. The total energy of the ball is converted entirely into potential energy before it can reach the top of the hill. We do not expect, even after repeated attempts, for the 100-J ball to ever be found beyond the hill. Therefore, the probability that the ball passes over the hill is exactly 0, and probability it is turned back or “reflected” by the hill is exactly 1. The ball never makes it over the hill. The existence of the ball beyond the hill is an impossibility or “energetically forbidden.”
However, according to quantum mechanics, the ball has a wave function and this function is defined over all space. The wave function may be highly localized, but there is always a chance that as the ball encounters the hill, the ball will suddenly be found beyond it. Indeed, this probability is appreciable if the “wave packet” of the ball is wider than the barrier.
View this interactive simulation for a simulation of tunneling.
In the language of quantum mechanics, the hill is characterized by a potential barrier. A finite-height square barrier is described by the following potential-energy function:
\[U(x) =
\begin{cases}
0, & \mbox{when } x < 0 \\[4pt] U_0, & \mbox{when } 0 \leq x \leq L \\[4pt] 0, & \mbox{when } x > L
\end{cases} \label{PIBPotential} \]
The potential barrier is illustrated in Figure \(\PageIndex{1}\). When the height \(U_0\) of the barrier is infinite, the wave packet representing an incident quantum particle is unable to penetrate it, and the quantum particle bounces back from the barrier boundary, just like a classical particle. When the width \(L\) of the barrier is infinite and its height is finite, a part of the wave packet representing an incident quantum particle can filter through the barrier boundary and eventually perish after traveling some distance inside the barrier.
When both the width \(L\) and the height \(U_0\) are finite, a part of the quantum wave packet incident on one side of the barrier can penetrate the barrier boundary and continue its motion inside the barrier, where it is gradually attenuated on its way to the other side. A part of the incident quantum wave packet eventually emerges on the other side of the barrier in the form of the transmitted wave packet that tunneled through the barrier. How much of the incident wave can tunnel through a barrier depends on the barrier width \(L\) and its height \(U_0\), and on the energy \(E\) of the quantum particle incident on the barrier. This is the physics of tunneling.
Barrier penetration by quantum wave functions was first analyzed theoretically by Friedrich Hund in 1927, shortly after Schrӧdinger published the equation that bears his name. A year later, George Gamow used the formalism of quantum mechanics to explain the radioactive \(\alpha\)-decay of atomic nuclei as a quantum-tunneling phenomenon. The invention of the tunnel diode in 1957 made it clear that quantum tunneling is important to the semiconductor industry. In modern nanotechnologies, individual atoms are manipulated using a knowledge of quantum tunneling.
Tunneling and the Wavfunction
Suppose a uniform and time-independent beam of electrons or other quantum particles with energy \(E\) traveling along the x-axis (in the positive direction to the right) encounters a potential barrier described by Equation \ref{PIBPotential}. The question is: What is the probability that an individual particle in the beam will tunnel through the potential barrier? The answer can be found by solving the boundary-value problem for the time-independent Schrӧdinger equation for a particle in the beam. The general form of this equation is given by Equation \ref{TISE}, which we reproduce here:
\[-\frac{\hbar^2}{2m} \frac{d^2\psi(x)}{dx^2} + U(x)\psi(x) = E\psi(x), \label{TISE} \]
where \(-\infty < x < +\infty\).
The potential function \(U(x)\) in Equation \ref{TISE} is defined by Equation \ref{PIBPotential}. We assume that the given energy \(E\) of the incoming particle is smaller than the height \(U_0\) of the potential barrier, \(E < U_0\), because this is the interesting physical case. Knowing the energy \(E\) of the incoming particle, our task is to solve Equation \ref{TISE} for a function \(\psi(x)\) that is continuous and has continuous first derivatives for all x. In other words, we are looking for a “smooth-looking” solution (because this is how wave functions look) that can be given a probabilistic interpretation so that \(|\psi(x)|^2 = \psi^*(x)\psi(x)\) is the probability density.
We divide the real axis into three regions with the boundaries defined by the potential function in Equation \ref{PIBPotential} (illustrated in Figure \(\PageIndex{1}\)) and transcribe Equation \ref{TISE} for each region. Denoting by \(\psi_I(x)\) the solution in region \(I\) for \(x < 0\), by \(\psi_{II}(x)\) the solution in region \(II\) for \(0 \leq x \leq L\), and by \(\psi_{III}(x)\) the solution in region \(III\) for \(x > L\), the stationary Schrӧdinger equation has the following forms in these three regions:
\[-\dfrac{\hbar^2}{2m} \frac{d^2\psi_I(x)}{dx^2} = E\psi_I(x), \label{R1} \]
in region \(I\): \(-\infty < x < 0,\)
\[-\dfrac{\hbar^2}{2m} \frac{d^2\psi_{II}(x)}{dx^2} + U_0\psi_{II}(x) = E\psi_{II}(x) \label{R2} \]
in region \(II\): \(0 < x < L,\)
\[-\dfrac{\hbar^2}{2m} \frac{d^2\psi_{III}(x)}{dx^2} = E\psi_{III}(x) \label{R3} \]
in region \(III\): \(L < x < +\infty,\)
The continuity condition at region boundaries requires that:
\[\psi_I(0) = \psi_{II}(0) \label{Bound1} \]
at the boundary between regions \(I\) and \(II\)
and
\[\psi_{II}(L) = \psi_{III}(L) \label{Bound2} \]
at the boundary between regions \(II\) and \(III\).
The “smoothness” condition requires the first derivative of the solution be continuous at region boundaries:
\[ \left.\dfrac{d\psi_I(x)}{dx}\right|_{x = 0} = \left.\dfrac{d\psi_{II}(x)}{dx}\right|_{x = 0} \label{smooth1} \]
at the boundary between regions \(I\) and \(II\)
and
\[ \left.\dfrac{d\psi_{II}(x)}{dx}\right|_{x = L} = \left.\dfrac{d\psi_{III}(x)}{dx}\right|_{x = L} \label{smooth2} \]
at the boundary between regions \(II\) and \(III\).
In what follows, we find the functions \(\psi_I(x), \, \psi_{II}(x)\), and \(\psi_{III}(x)\).
We can easily verify (by substituting into the original equation and differentiating) that in regions \(I\) and \(III\), the solutions must be in the following general forms:
\[\psi_I(x) = Ae^{+ikx} + Be^{-ikx} \label{Eq5A} \]
\[\psi_{III}(x) = Fe^{+ikx} + Ge^{-ikx} \label{Eq5B} \]
where \(k = \sqrt{2mE}/\hbar\) is a wave number and the complex exponent denotes oscillations,
\[e^{\pm ikx} = \cos \, kx \pm i \, \sin \, kx. \nonumber \]
The constants \(A\), \(B\), \(F\), and \(G\) in Equations \ref{Eq5A} and \ref{Eq5B} may be complex. These solutions are illustrated in Figure \(\PageIndex{2}\). In region I, there are two waves—one is incident (moving to the right) and one is reflected (moving to the left)—so none of the constants \(A\) and \(B\) in Equation \ref{Eq5A} may vanish. In region III, there is only one wave (moving to the right), which is the transmitted wave, so the constant \(G\) must be zero in Equation \ref{Eq5B}, \(G = 0\). We can write explicitly that the incident wave is \(\psi_{in}(x) = Ae^{+ikx}\) and that the reflected wave is \(\psi_{ref}(x) = Be^{-ikx}\), and that the transmitted wave is \(\psi_{tra}(x) = Fe^{+ikx}\). The amplitude of the incident wave is
\[ \begin{align*}|\psi_{in}(x)|^2 &= \psi_{in}^*(x)\psi_{in}(x) \\[4pt] &= (Ae^{+ikx})^*Ae^{+ikx} \\[4pt] &= A^*e^{-ikx}Ae^{+ikx} \\[4pt] &= A^*A = |A|^2. \end{align*} \nonumber \]
Similarly, the amplitude of the reflected wave is \(|\psi_{ref}(x)|^2 = |B|^2\) and the amplitude of the transmitted wave is \(|\psi_{tra}(x)|^2 = |F|^2\). We know from the theory of waves that the square of the wave amplitude is directly proportional to the wave intensity. If we want to know how much of the incident wave tunnels through the barrier, we need to compute the square of the amplitude of the transmitted wave. The transmission probability or tunneling probability is the ratio of the transmitted intensity (\(|F|^2\)) to the incident intensity (\(|A|^2\)), written as
\[ \begin{align} T(L, E) &= \frac{|\psi_{tra}(x)|^2}{|\psi_{in}(x)|^2} \\[4pt] &= \frac{|F|^2}{|A|^2} \\[4pt] &= \left|\frac{F}{A}\right|^2 \label{trans} \end{align} \]
where \(L\) is the width of the barrier and \(E\) is the total energy of the particle. This is the probability an individual particle in the incident beam will tunnel through the potential barrier. Intuitively, we understand that this probability must depend on the barrier height \(U_0\).
In region II, the terms in equation Equation \ref{R2} can be rearranged to
\[\frac{d^2\psi_{II}(x)}{dx^2} = \beta^2 \psi_{II}(x) \label{eq10} \]
where \(\beta^2\) is positive because \(U_0 > E\) and the parameter \(\beta\) is a real number,
\[\beta^2 = \frac{2m}{\hbar^2}(U_0 - E). \label{eq20} \]
The general solution to Equation \ref{eq10} is not oscillatory (unlike in the other regions) and is in the form of exponentials that describe a gradual attenuation of \(\psi_{II}(x)\),
\[\psi_{II}(x) = Ce^{-\beta x} + De^{+\beta x}. \label{eq25B} \]
The two types of solutions in the three regions are illustrated in Figure \(\PageIndex{2}\).
Now we use the boundary conditions to find equations for the unknown constants. Equations \ref{Eq5A} and \ref{eq25B} are substituted into Equation \ref{Bound1} to give
\[A + B = C + D. \nonumber \]
Equations \ref{eq25B} and \ref{Eq5B} are substituted into Equation \ref{Bound2} to give
\[Ce^{-\beta L} + De^{+\beta L} = Fe^{+ik L}. \nonumber \]
Similarly, we substitute Equations \ref{Eq5A} and \ref{eq25B} into Equation \ref{smooth1}, differentiate, and obtain
\[-ik(A - B) = \beta(D - C). \nonumber \]
Similarly, the boundary condition Equation \ref{smooth2} reads explicitly
\[\beta (De^{+\beta L} - Ce^{-\beta L}) = +ikFe^{+ikL}. \nonumber \]
We now have four equations for five unknown constants. However, because the quantity we are after is the transmission coefficient (\(T\)), defined in Equation \ref{trans} by the fraction \(F/A\), the number of equations is exactly right because when we divide each of the above equations by \(A\), we end up having only four unknown fractions: \(B/A\), \(C/A\), \(D/A\), and \(F/A\), three of which can be eliminated to find \(F/A\). The actual algebra that leads to expression for \(F/A\) is pretty lengthy, but it can be done either by hand or with a help of computer software. The end result is
\[\frac{F}{A} = \frac{e^{-ikL}}{\cosh \, (\beta L) + i (\gamma /2) \, \sinh \, (\beta L)}. \label{eq40} \]
In deriving Equation \ref{eq40}, to avoid the clutter, we use the substitutions \(\gamma \equiv \beta/k - k/\beta\), and the definition of hyperbolic functions:
\[\cosh \, y = \frac{e^y + e^{-y}}{2} \nonumber \]
and
\[\sinh \, y = \frac{e^y - e^{-y}}{2}. \nonumber \]
We substitute Equation \ref{eq40} into Equation \ref{trans} and obtain the exact expression for the transmission coefficient for the barrier,
\[T(L,E) = \left(\frac{F}{A}\right)^*\frac{F}{A} = \frac{e^{+ikL}}{\cosh \, (\beta L) - i (\gamma /2) \, \sinh \, (\beta L)} \cdot \frac{e^{-ikL}}{\cosh \, (\beta L) + i (\gamma /2) \, \sinh \, (\beta L)}. \nonumber \]
or
\[T(L,E) = \frac{1}{\cosh^2 \, (\beta L) + (\gamma /2)^2 \, \sinh^2 \, (\beta L)}. \label{eqTrans1} \]
where
\[\left(\frac{\gamma}{2}\right)^2 = \frac{1}{4}\left( \frac{1 - E/U_0}{E/U_0} + \frac{E/U_0}{1 - E/U_0} - 2\right). \nonumber \]
For a wide and high barrier that transmits poorly, Equation \ref{eqTrans1} can be approximated by
\[T(L,E) \approx 16 \frac{E}{U_0}\left(1 - \dfrac{E}{U_0}\right) e^{-2\beta L}. \label{eqTrans2} \]
Whether it is the exact expression (Equation \ref{eqTrans1}) or the approximate expression (Equation \ref{eqTrans2}), we see that the tunneling effect very strongly depends on the width \(L\) of the potential barrier. In the laboratory, we can adjust both the potential height \(U_0\) and the width \(L\) to design nano-devices with desirable transmission coefficients.
Two copper nanowires are insulated by a copper oxide nano-layer that provides a 10.0-eV potential barrier. Estimate the tunneling probability between the nanowires by 7.00-eV electrons through a 5.00-nm thick oxide layer. What if the thickness of the layer were reduced to just 1.00 nm? What if the energy of electrons were increased to 9.00 eV?
Strategy
Treating the insulating oxide layer as a finite-height potential barrier, we use Equation \ref{eqTrans2}. We identify \(U_0 = 10.0 \, eV, \, E_1 = 7.00 \, eV, \, E_2 = 9.00 \, eV, \, L_1 = 5.00 \, nm\), and \(L_2 = 1.00 \, nm\). We use Equation \ref{eq20} to compute the exponent. Also, we need the rest mass of the electron \(m = 511 \, keV/c^2\) and Planck’s constant \(\hbar = 0.1973 \, keV \cdot nm/c\). It is typical for this type of estimate to deal with very small quantities that are often not suitable for handheld calculators. To make correct estimates of orders, we make the conversion \(e^y = 10^{y/ln \, 10}\).
Solution
Constants:
\[\frac{2m}{\hbar^2} = \frac{2(511 \, keV/c^2)}{(0.1973 \, keV \cdot nm/c^2)^2} = 26,254 \frac{1}{keV \cdot (nm)^2}, \nonumber \]
\[\beta = \sqrt{\frac{2m}{\hbar^2}(U_0 - E)} = \sqrt{26,254\frac{(10.0 \, eV - E)}{keV \cdot (nm)^2}} = \sqrt{26.254(10.00 - E)/eV} \frac{1}{nm}.\nonumber \]
For a lower-energy electron with \(E_1 = 7.00 \, eV\):
\[\beta_1 = \sqrt{26.254(10.00 \, eV - E_1)/eV}\frac{1}{nm} = \sqrt{26.254(10.00 - 7.00)}\frac{1}{nm} = \frac{8.875}{nm},\nonumber \]
\[T(L,E) = 16\frac{E_1}{U_0}\left(1 - \frac{E_1}{U_0}\right)e^{-2\beta_1L} = 16\frac{7}{10}\left(1 - \frac{7}{10}\right)e^{-17.75 \, L/nm} = 3.36 \, e^{-17.75 \, L/nm}\nonumber \]
For a higher-energy electron with \(E_2 = 9.00 \, eV\):
\[\beta_2 = \sqrt{26.254(10.00 \, eV - E_2)/eV}\frac{1}{nm} = \sqrt{26.254 (10.00 - 9.00)}\frac{1}{nm} = \frac{5.124}{nm},\nonumber \]
\[T(L,E_2) = 16\frac{E_2}{U_0}\left( 1 - \frac{E_2}{U_0}\right) e^{-2\beta_2 L} = 16 \frac{9}{10} \left(1 - \frac{9}{10}\right) e^{-10.25 \, L/nm} = 1.44 \, e^{-10.25 \, L/nm}\nonumber \]
For a broad barrier with \(L_1 = 5.00 \, nm\):
\[T(L_1,E_1) = 3.36 e^{-17.75 \, L_1/nm} = 3.36 e^{-17.75 \cdot 5.00 nm/nm} = 3.36e^{-88} = 3.36(6.2 \times 10^{-39}) = 2.1\% \times 10^{-36}\nonumber \]
\[T(L_1,E_2) = 1.44 e^{-10.25 \, L_1/nm} = 1.44 e^{-10.25 \cdot 5.00 \, nm/nm} = 1.44^{-51.2} = 1.44(5.81 \times 10^{-12}) = 8.36\% \times 10^{-25}\nonumber \]
For a narrower barrier with \(L_2 = 1.00 \, nm\):
\[T(L_2,E_1) = 3.36e^{−17.75 \, L_2/nm} = 3.36 e^{−17.75 \cdot 1.00 \, nm/nm} = 3.36 e^{−17.75} = 3.36(5.1 \times 10^{−7}) = 1.7\% \times 10^{−4},\nonumber \]
\[T(L_2,E_2) = 1.44e^{−10.25 \, L_2/nm} = 1.44e^{-10.25 \cdot 1.00 \, nm/nm} = 1.44e^{−10.25} = 1.44(3.53 \times 10^{−5}) = 5.09\% \times 10^{-7}.\nonumber \]
Significance
We see from these estimates that the probability of tunneling is affected more by the width of the potential barrier than by the energy of an incident particle. In today’s technologies, we can manipulate individual atoms on metal surfaces to create potential barriers that are fractions of a nanometer, giving rise to measurable tunneling currents. One of many applications of this technology is the scanning tunneling microscope (STM), which we discuss later in this section.
A proton with kinetic energy 1.00 eV is incident on a square potential barrier with height 10.00 eV. If the proton is to have the same transmission probability as an electron of the same energy, what must the width of the barrier be relative to the barrier width encountered by an electron?
- Answer
-
\[L_{proton}/L_{electron} = \sqrt{m_e/m_p} = 2.3\% \nonumber \]
Radioactive Decay
In 1928, Gamow identified quantum tunneling as the mechanism responsible for the radioactive decay of atomic nuclei. He observed that some isotopes of thorium, uranium, and bismuth disintegrate by emitting α-particles (which are doubly ionized helium atoms or, simply speaking, helium nuclei). In the process of emitting an α-particle, the original nucleus is transformed into a new nucleus that has two fewer neutrons and two fewer protons than the original nucleus. The α-particles emitted by one isotope have approximately the same kinetic energies. When we look at variations of these energies among isotopes of various elements, the lowest kinetic energy is about 4 MeV and the highest is about 9 MeV, so these energies are of the same order of magnitude. This is about where the similarities between various isotopes end.
When we inspect half-lives (a half-life is the time in which a radioactive sample loses half of its nuclei due to decay), different isotopes differ widely. For example, the half-life of polonium-214 is 160 µs and the half-life of uranium is 4.5 billion years. Gamow explained this variation by considering a ‘spherical-box’ model of the nucleus, where α-particles can bounce back and forth between the walls as free particles. The confinement is provided by a strong nuclear potential at a spherical wall of the box. The thickness of this wall, however, is not infinite but finite, so in principle, a nuclear particle has a chance to escape this nuclear confinement. On the inside wall of the confining barrier is a high nuclear potential that keeps the α-particle in a small confinement. But when an α-particle gets out to the other side of this wall, it is subject to electrostatic Coulomb repulsion and moves away from the nucleus. This idea is illustrated in Figure \(\PageIndex{3}\). The width \(L\) of the potential barrier that separates an α-particle from the outside world depends on the particle’s kinetic energy \(E\). This width is the distance between the point marked by the nuclear radius \(R\) and the point \(R_0\) where an α-particle emerges on the other side of the barrier, \(L = R_0 - R\). At the distance \(R_0\), its kinetic energy must at least match the electrostatic energy of repulsion, \(E = (4\pi\epsilon_0)^{-1}Ze^2/R_0\) (where +Ze is the charge of the nucleus). In this way we can estimate the width of the nuclear barrier,
\[L = \frac{e^2}{4\pi\epsilon_0} \frac{Z}{E} - R. \nonumber \]
We see from this estimate that the higher the energy of α-particle, the narrower the width of the barrier that it is to tunnel through. We also know that the width of the potential barrier is the most important parameter in tunneling probability. Thus, highly energetic α-particles have a good chance to escape the nucleus, and, for such nuclei, the nuclear disintegration half-life is short. Notice that this process is highly nonlinear, meaning a small increase in the α-particle energy has a disproportionately large enhancing effect on the tunneling probability and, consequently, on shortening the half-life. This explains why the half-life of polonium that emits 8-MeV α-particles is only hundreds of milliseconds and the half-life of uranium that emits 4-MeV α-particles is billions of years.
Field Emission
Field emission is a process of emitting electrons from conducting surfaces due to a strong external electric field that is applied in the direction normal to the surface (Figure \(\PageIndex{4}\)). As we know from our study of electric fields in earlier chapters, an applied external electric field causes the electrons in a conductor to move to its surface and stay there as long as the present external field is not excessively strong. In this situation, we have a constant electric potential throughout the inside of the conductor, including its surface. In the language of potential energy, we say that an electron inside the conductor has a constant potential energy \(U(x) - -U_0\) (here, the x means inside the conductor).
In the situation represented in Figure \(\PageIndex{4}\), where the external electric field is uniform and has magnitude \(E_g\), if an electron happens to be outside the conductor at a distance x away from its surface, its potential energy would have to be \(U(x) = -eE_gx\) (here, x denotes distance to the surface). Taking the origin at the surface, so that \(x = 0\) is the location of the surface, we can represent the potential energy of conduction electrons in a metal as the potential energy barrier shown in Figure \(\PageIndex{5}\). In the absence of the external field, the potential energy becomes a step barrier defined by \(U(x \leq 0) = -U_0\) and by \(U(x > 0) = 0\).
When an external electric field is strong, conduction electrons at the surface may get detached from it and accelerate along electric field lines in a direction antiparallel to the external field, away from the surface. In short, conduction electrons may escape from the surface. The field emission can be understood as the quantum tunneling of conduction electrons through the potential barrier at the conductor’s surface. The physical principle at work here is very similar to the mechanism of α-emission from a radioactive nucleus.
Suppose a conduction electron has a kinetic energy E (the average kinetic energy of an electron in a metal is the work function ϕ for the metal and can be measured, as discussed for the photoelectric effect in Photons and Matter Waves), and an external electric field can be locally approximated by a uniform electric field of strength \(E_g\). The width L of the potential barrier that the electron must cross is the distance from the conductor’s surface to the point outside the surface where its kinetic energy matches the value of its potential energy in the external field. In Figure \(\PageIndex{5}\), this distance is measured along the dashed horizontal line \(U(x) = E\) from \(x = 0\) to the intercept with \(U(x) = -eE_gx\), so the barrier width is
\[L = \frac{e^{-1}E}{E_g} = \frac{e^{-1}\phi}{E_g} \nonumber \]
We see that L is inversely proportional to the strength \(E_g\) of an external field. When we increase the strength of the external field, the potential barrier outside the conductor becomes steeper and its width decreases for an electron with a given kinetic energy. In turn, the probability that an electron will tunnel across the barrier (conductor surface) becomes exponentially larger. The electrons that emerge on the other side of this barrier form a current (tunneling-electron current) that can be detected above the surface. The tunneling-electron current is proportional to the tunneling probability. The tunneling probability depends nonlinearly on the barrier width L, and L can be changed by adjusting \(E_g\). Therefore, the tunneling-electron current can be tuned by adjusting the strength of an external electric field at the surface. When the strength of an external electric field is constant, the tunneling-electron current has different values at different elevations L above the surface.
The quantum tunneling phenomenon at metallic surfaces, which we have just described, is the physical principle behind the operation of the scanning tunneling microscope (STM), invented in 1981 by Gerd Binnig and Heinrich Rohrer. The STM device consists of a scanning tip (a needle, usually made of tungsten, platinum-iridium, or gold); a piezoelectric device that controls the tip’s elevation in a typical range of 0.4 to 0.7 nm above the surface to be scanned; some device that controls the motion of the tip along the surface; and a computer to display images. While the sample is kept at a suitable voltage bias, the scanning tip moves along the surface (Figure \(\PageIndex{6}\)) and the tunneling-electron current between the tip and the surface is registered at each position.
The amount of the current depends on the probability of electron tunneling from the surface to the tip, which, in turn, depends on the elevation of the tip above the surface. Hence, at each tip position, the distance from the tip to the surface is measured by measuring how many electrons tunnel out from the surface to the tip. This method can give an unprecedented resolution of about 0.001 nm, which is about 1% of the average diameter of an atom. In this way, we can see individual atoms on the surface, as in the image of a carbon nanotube in Figure \(\PageIndex{7}\).
Resonant Quantum Tunneling
Quantum tunneling has numerous applications in semiconductor devices such as electronic circuit components or integrated circuits that are designed at nanoscales; hence, the term ‘nanotechnology.’ For example, a diode (an electric-circuit element that causes an electron current in one direction to be different from the current in the opposite direction, when the polarity of the bias voltage is reversed) can be realized by a tunneling junction between two different types of semiconducting materials. In such a tunnel diode, electrons tunnel through a single potential barrier at a contact between two different semiconductors. At the junction, tunneling-electron current changes nonlinearly with the applied potential difference across the junction and may rapidly decrease as the bias voltage is increased. This is unlike the Ohm’s law behavior that we are familiar with in household circuits. This kind of rapid behavior (caused by quantum tunneling) is desirable in high-speed electronic devices.
Another kind of electronic nano-device utilizes resonant tunneling of electrons through potential barriers that occur in quantum dots. A quantum dot is a small region of a semiconductor nanocrystal that is grown, for example, in a silicon or aluminum arsenide crystal. Figure \(\PageIndex{8a}\) shows a quantum dot of gallium arsenide embedded in an aluminum arsenide wafer. The quantum-dot region acts as a potential well of a finite height (Figure \(\PageIndex{8b}\)) that has two finite-height potential barriers at dot boundaries. Similarly, as for a quantum particle in a box (that is, an infinite potential well), lower-lying energies of a quantum particle trapped in a finite-height potential well are quantized. The difference between the box and the well potentials is that a quantum particle in a box has an infinite number of quantized energies and is trapped in the box indefinitely, whereas a quantum particle trapped in a potential well has a finite number of quantized energy levels and can tunnel through potential barriers at well boundaries to the outside of the well. Thus, a quantum dot of gallium arsenide sitting in aluminum arsenide is a potential well where low-lying energies of an electron are quantized, indicated as \(E_{dot}\) in part (b) in the figure. When the energy \(E_{electron}\) of an electron in the outside region of the dot does not match its energy \(E_{dot}\) that it would have in the dot, the electron does not tunnel through the region of the dot and there is no current through such a circuit element, even if it were kept at an electric voltage difference (bias). However, when this voltage bias is changed in such a way that one of the barriers is lowered, so that \(E_{dot}\) and \(E_{electron}\) become aligned, as seen in part (c) of the figure, an electron current flows through the dot. When the voltage bias is now increased, this alignment is lost and the current stops flowing. When the voltage bias is increased further, the electron tunneling becomes improbable until the bias voltage reaches a value for which the outside electron energy matches the next electron energy level in the dot. The word ‘resonance’ in the device name means that the tunneling-electron current occurs only when a selected energy level is matched by tuning an applied voltage bias, such as in the operation mechanism of the resonant-tunneling diode just described. Resonant-tunneling diodes are used as super-fast nano-switches.