# 4.2: Square Potential Barrier

- Page ID
- 15744

\( \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}}} \)

Consider a particle of mass \(m\) and energy \(E>0\) interacting with the simple square potential barrier

\begin{equation}V(x)=\left\{\begin{array}{ll}

V_{0} & \text { for } 0 \leq x \leq a \\

0 & \text { otherwise }

\end{array}\right.\end{equation}

where \(V_0>0\). In the regions to the left and to the right of the barrier, \(\psi(x)\) satisfies \[\label{e5.15} \frac{d^{\,2} \psi}{d x^{\,2}} = - k^{\,2}\,\psi,\] where \(k\) is given by Equation ([e5.6]).

Let us adopt the following solution of the previous equation to the left of the barrier (i.e., \(x<0\)): \[\psi(x) = {\rm e}^{\,{\rm i}\,k\,x} + R\,{\rm e}^{-{\rm i}\,k\,x}.\] This solution consists of a plane-wave of unit amplitude traveling to the right [because the time-dependent wavefunction is multiplied by \(\exp(-{\rm i}\,\omega\,t)\), where \(\omega=E/\hbar>0\)], and a plane wave of complex amplitude \(R\) traveling to the left. We interpret the first plane wave as an incoming particle (or, rather, a stream of incoming particles), and the second as a particle (or stream of particles) reflected by the potential barrier. Hence, \(|R|^{\,2}\) is the probability of reflection. This can be seen by calculating the probability current ([eprobc]) in the region \(x<0\), which takes the form \[j_l=v\,(1-|R|^{\,2}),\] where \(v = p/m=\hbar\,k/m\) is the classical particle velocity.

Let us adopt the following solution to Equation ([e5.15]) to the right of the barrier (i.e. \(x>a\)): \[\psi(x) = T\,{\rm e}^{\,{\rm i}\,k\,x}.\] This solution consists of a plane-wave of complex amplitude \(T\) traveling to the right. We interpret this as a particle (or stream of particles) transmitted through the barrier. Hence, \(|T|^{\,2}\) is the probability of transmission. The probability current in the region \(x>a\) takes the form \[j_r = v\, |T|^{\,2}.\] Now, according to Equation ([ediffp]), in a stationary state (i.e., \(\partial |\psi|^{\,2}/\partial t = 0\)), the probability current is a spatial constant (i.e., \(\partial j/\partial x = 0\)). Hence, we must have \(j_l=j_r\), or \[\label{e5.20} |R|^{\,2} + |T|^{\,2} = 1.\] In other words, the probabilities of reflection and transmission sum to unity, as must be the case, because reflection and transmission are the only possible outcomes for a particle incident on the barrier.

Inside the barrier (i.e., \(0\leq x \leq a\)), \(\psi(x)\) satisfies \[\label{e5.21} \frac{d^{\,2} \psi}{d x^{\,2}} = - q^{\,2}\,\psi,\] where \[q^{\,2} = \frac{2\,m\,(E-V_0)}{\hbar^{\,2}}.\]

Let us, first of all, consider the case where \(E> V_0\). In this case, the general solution to Equation ([e5.21]) inside the barrier takes the form \[\psi(x) = A\,{\rm e}^{\,{\rm i}\,q\,x} +B\,{\rm e}^{-{\rm i}\,q\,x},\] where \(q=\sqrt{2\,m\,(E-V_0)/\hbar^{\,2}}\).

Now, the boundary conditions at the edges of the barrier (i.e., at \(x=0\) and \(x=a\)) are that \(\psi\) and \(d\psi/d x\) are both continuous. These boundary conditions ensure that the probability current ([eprobc]) remains finite and continuous across the edges of the boundary, as must be the case if it is to be a spatial constant.

Continuity of \(\psi\) and \(d\psi/d x\) at the left edge of the barrier (i.e., \(x=0\)) yields \[\begin{aligned} 1 + R &= A+B,\\[0.5ex] k\,(1-R) &= q\,(A-B).\end{aligned}\] Likewise, continuity of \(\psi\) and \(d\psi/d x\) at the right edge of the barrier (i.e., \(x=a\)) gives \[\begin{aligned} A\, {\rm e}^{\,{\rm i}\,q\,a}+ B \,{\rm e}^{-{\rm i}\,q\,a} &= T\,{\rm e}^{\,{\rm i}\,k\,a},\\[0.5ex] q\left(A\,{\rm e}^{\,{\rm i}\,q\,a} -B \,{\rm e}^{-{\rm i}\,q\,a}\right) &=k\,T\,{\rm e}^{\,{\rm i}\,k\,a}.\end{aligned}\] After considerable algebra, the previous four equations yield \[\label{e5.28} |R|^{\,2} = \frac{(k^{\,2}-q^{\,2})^{\,2}\,\sin^2(q\,a)}{4\,k^{\,2}\,q^{\,2} + (k^{\,2}-q^{\,2})^{\,2}\,\sin^2(q\,a)},\] and \[\label{e5.29} |T|^{\,2} = \frac{4\,k^{\,2}\,q^{\,2}}{4\,k^{\,2}\,q^{\,2} + (k^{\,2}-q^{\,2})^{\,2}\,\sin^2(q\,a)}.\] Note that the previous two expression satisfy the constraint ([e5.20]).

It is instructive to compare the quantum mechanical probabilities of reflection and transmission—([e5.28]) and ([e5.29]), respectively—with those derived from classical physics. Now, according to classical physics, if a particle of energy \(E\) is incident on a potential barrier of height \(V_0<E\) then the particle slows down as it passes through the barrier, but is otherwise unaffected. In other words, the classical probability of reflection is zero, and the classical probability of transmission is unity.

**Figure 10:** *Transmission (solid-curve) and reflection (dashed-curve) probabilities for a square potential barrier of width *\(\begin{equation}a=1.25 \lambda\end{equation}\)*, where *\(\lambda\) *the free-space de Broglie wavelength, as a function of the ratio of the height of the barrier, *\(\begin{equation}V_{0}\end{equation}\)*, to the energy, *\(\begin{equation}E\end{equation}\)*, of the incident particle.*

**Figure 11:** *Transmission (solid-curve) and reflection (dashed-curve) probabilities for a particle of energy incident on a square potential barrier of height *\(\begin{equation}V_{0}=0.75 E\end{equation}\)*, as a function of the ratio of the width of the barrier, *\(\begin{equation}a\end{equation}\), *to the free-space de Broglie wavelength, *\(\lambda\)

The reflection and transmission probabilities obtained from Equations ([e5.28]) and ([e5.29]), respectively, are plotted in Figures [fb1] and [fb2]. It can be seen, from Figure [fb1], that the classical result, \(|R|^{\,2}=0\) and \(|T|^{\,2}=1\), is obtained in the limit where the height of the barrier is relatively small (i.e., \(V_0\ll E\)). However, when \(V_0\) is of order \(E\), there is a substantial probability that the incident particle will be reflected by the barrier. According to classical physics, reflection is impossible when \(V_0 < E\).

It can also be seen, from Figure [fb2], that at certain barrier widths the probability of reflection goes to zero. It turns out that this is true irrespective of the energy of the incident particle. It is evident, from Equation ([e5.28]), that these special barrier widths correspond to \[q\,a = n\,\pi,\] where \(n=1,2,3,\cdots\). In other words, the special barriers widths are integer multiples of half the de Broglie wavelength of the particle inside the barrier. There is no reflection at the special barrier widths because, at these widths, the backward traveling wave reflected from the left edge of the barrier interferes destructively with the similar wave reflected from the right edge of the barrier to give zero net reflected wave.

Let us, now, consider the case \(E< V_0\). In this case, the general solution to Equation ([e5.21]) inside the barrier takes the form \[\psi(x) = A\,{\rm e}^{\,q\,x} + B\,{\rm e}^{-q\,x},\] where \(q=\sqrt{2\,m\,(V_0-E)/\hbar^{\,2}}\). Continuity of \(\psi\) and \(d \psi/d x\) at the left edge of the barrier (i.e., \(x=0\)) yields \[\begin{aligned} 1 + R &=A+B,\\[0.5ex] {\rm i}\,k\,(1-R) &= q\,(A-B).\end{aligned}\] Likewise, continuity of \(\psi\) and \(d\psi/d x\) at the right edge of the barrier (i.e., \(x=a\)) gives \[\begin{aligned} A\, {\rm e}^{\,q\,a}+ B \,{\rm e}^{-q\,a} &= T\,{\rm e}^{\,{\rm i}\,k\,a},\\[0.5ex] q\left(A\, {\rm e}^{\,q\,a}-B \,{\rm e}^{-q\,a}\right) &= {\rm i}\,k\,T\,{\rm e}^{\,{\rm i}\,k\,a}.\end{aligned}\] After considerable algebra, the previous four equations yield \[\label{e5.36} |R|^{\,2} = \frac{(k^{\,2}+q^{\,2})^{\,2}\,\sinh^2(q\,a)}{4\,k^{\,2}\,q^{\,2} + (k^{\,2}+q^{\,2})^{\,2}\,\sinh^2(q\,a)},\] and \[\label{e5.37} |T|^{\,2} = \frac{4\,k^{\,2}\,q^{\,2}}{4\,k^{\,2}\,q^{\,2} + (k^{\,2}+q^{\,2})^{\,2}\,\sinh^2(q\,a)}.\] These expressions can also be obtained from Equations ([e5.28]) and ([e5.29]) by making the substitution \(q\rightarrow -{\rm i}\,q\). Note that Equations ([e5.36]) and ([e5.37]) satisfy the constraint ([e5.20]).

It is again instructive to compare the quantum mechanical probabilities of reflection and trans/-mission—([e5.36]) and ([e5.37]), respectively—with those derived from classical physics. Now, according to classical physics, if a particle of energy \(E\) is incident on a potential barrier of height \(V_0>E\) then the particle is reflected. In other words, the classical probability of reflection is unity, and the classical probability of transmission is zero.

**Figure 12:** *Transmission (solid-curve) and reflection (dashed-curve) probabilities for a square potential barrier of width *\(\begin{equation}a=0.5 \lambda \text { , where } \lambda \text { is }\end{equation}\)*the free-space de Broglie wavelength, as a function of the ratio of the energy, , of the incoming particle to the height, *\(\begin{equation}V_{0}\end{equation}\),*of the barrier.*

**Figure 13:** *Transmission (solid-curve) and reflection (dashed-curve) probabilities for a particle of energy incident on a square potential barrier of height *\(\begin{equation}V_{0}=(4 / 3) E\end{equation}\)*, as a function of the ratio of the width of the barrier, , to the free-space de Broglie wavelength, *\(\lambda\)

The reflection and transmission probabilities obtained from Equations ([e5.36]) and ([e5.37]), respectively, are plotted in Figures [fb3] and [fb4]. It can be seen, from Figure [fb3], that the classical result, \(|R|^{\,2}=1\) and \(|T|^{\,2}=0\), is obtained for relatively thin barriers (i.e., \(q\,a\sim 1\)) in the limit where the height of the barrier is relatively large (i.e., \(V_0\gg E\)). However, when \(V_0\) is of order \(E\), there is a substantial probability that the incident particle will be transmitted by the barrier. According to classical physics, transmission is impossible when \(V_0 > E\).

It can also be seen, from Figure [fb4], that the transmission probability decays exponentially as the width of the barrier increases. Nevertheless, even for very wide barriers (i.e., \(q\,a\gg 1\)), there is a small but finite probability that a particle incident on the barrier will be transmitted. This phenomenon, which is inexplicable within the context of classical physics, is called *tunneling*.

## Contributors and Attributions

Richard Fitzpatrick (Professor of Physics, The University of Texas at Austin)

\( \newcommand {\ltapp} {\stackrel {_{\normalsize<}}{_{\normalsize \sim}}}\) \(\newcommand {\gtapp} {\stackrel {_{\normalsize>}}{_{\normalsize \sim}}}\) \(\newcommand {\btau}{\mbox{\boldmath$\tau$}}\) \(\newcommand {\bmu}{\mbox{\boldmath$\mu$}}\) \(\newcommand {\bsigma}{\mbox{\boldmath$\sigma$}}\) \(\newcommand {\bOmega}{\mbox{\boldmath$\Omega$}}\) \(\newcommand {\bomega}{\mbox{\boldmath$\omega$}}\) \(\newcommand {\bepsilon}{\mbox{\boldmath$\epsilon$}}\)