3.6: Wavefunctions Must Be Normalized
- Page ID
- 210798
\( \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}\)Since wavefunctions can in general be complex functions, the physical significance cannot be found from the function itself because the \(\sqrt {-1}\) is not a property of the physical world. Rather, the physical significance is found in the product of the wavefunction and its complex conjugate, i.e. the absolute square of the wavefunction, which also is called the square of the modulus.
\[ \Psi^*(r , t ) \psi (r , t ) = {|\Psi (r , t)|}^2 \label {3.6.1}\]
where \(r\) is a vector (x, y, z) specifying a point in three-dimensional space. The square is used, rather than the modulus itself, just like the intensity of a light wave depends on the square of the electric field. Remember that the Born interpretation is that \(\Psi^*(r_i)\Psi(r_i)\, d\tau\) is the probability that the electron is in the volume dτ located at \(r_i\). The Born interpretation therefore calls the wavefunction the probability amplitude, the absolute square of the wavefunction is called the probability density, and the probability density times a volume element in three-dimensional space (\(d\tau\)) is the probability.
Show that the square of the modulus of \(\Psi(r,t) = \psi(r) e^{-iωt}\) is time independent. What insight regarding stationary states do you gain from this proof?
Solution
The square of the modulus of a wavefunction is \(\Psi(r,t)^*\Psi(r,t)\) so for wavefunctions of this form, the square of the modulus is
\[ \begin{align*} \Psi(r,t)^*\Psi(r,t) &= \psi(r) \cancel{e^{+iωt}} \psi(r) \cancel{e^{-iωt}} \\[4pt] &=\psi(r)^2 \end{align*}\]
Hence, there is no time dependence to the modulus of wavefunctions of this work, which from the probability interpretation of the wavefunction means that the probability density is time-independent.
Normalization of the Wavefunction
A probability is a real number between 0 and 1, inclusive. An outcome of a measurement which has a probability 0 is an impossible outcome, whereas an outcome which has a probability 1 is a certain outcome. According to Equation \(\ref{3.6.1}\), the probability of a measurement of \(x\) yielding a result between \(-\infty\) and \(+\infty\) is
\[P_{x \in -\infty:\infty}(t) = \int_{-\infty}^{\infty}\vert\psi(x,t)\vert^{ 2} dx. \label{3.6.2}\]
However, a measurement of \(x\) must yield a value between \(-\infty\) and \(+\infty\), since the particle has to be located somewhere. It follows that \(P_{x \in -\infty:\infty}(t) =1\), or
\[\int_{-\infty}^{\infty}\vert\psi(x,t)\vert^{ 2} dx = 1, \label{3.6.3}\]
which is generally known as the normalization condition for the wavefunction.
Normalize the wavefunction of a Gaussian wave packet, centered on \(x=x_o\) with characteristic width \(\sigma\):
\[\psi(x) = \psi_0 {\rm e}^{-(x-x_0)^{ 2}/(4 \sigma^2)}.\label{3.6.4}\]
Solution
To determine the normalization constant \(\psi_0\), we simply substitute Equation \(\ref{3.6.4}\) into Equation \(\ref{3.6.3}\), to obtain
\[\vert\psi_0\vert^{ 2}\int_{-\infty}^{\infty}{\rm e}^{-(x-x_0)^{ 2}/(2 \sigma^2)} dx = 1. \nonumber\]
Changing the variable of integration to \(y=(x-x_0)/(\sqrt{2} \sigma)\), we get
\[\vert\psi_0\vert^{ 2}\sqrt{2} \sigma \int_{-\infty}^{\infty}{\rm e}^{-y^2} dy=1. \nonumber\]
However, from an integral table we know
\[\int_{-\infty}^{\infty}{\rm e}^{-y^2} dy = \sqrt{\pi},\nonumber\]
which implies that
\[\vert\psi_0\vert^{ 2} = \dfrac{1}{(2\pi \sigma^2)^{1/2}}. \nonumber\]
Hence, a general normalized Gaussian wavefunction takes the form
\[\psi(x) = \dfrac{e^{\rm{i} \phi}}{(2\pi \sigma^2)^{1/4}}e^{-(x-x_0)^2/(4 \sigma^2)} \nonumber\]
where \(\phi\) is an arbitrary real phase-angle.
Calculate the probability of finding an electron at \(L/2\) in a box of infinite height within an interval ranging from \(\dfrac {L}{2} - \dfrac {L}{200}\) to \(\dfrac {L}{2} + \dfrac {L}{200}\) for the \(n = 1\) and \(n = 2\) states. Since the length of the interval, \(L/100\), is small compared to \(L\), you can get an approximate answer without integrating.
Solution
The \(n = 1\) state
\[\psi(x)=\sqrt{\frac{2}{L}} \sin \left(\frac{n \pi x}{L}\right) \nonumber\]
For \(n=1\), \(\psi_{n=1}=\sqrt{\frac{2}{2}} \sin \left(\frac{\pi x}{2}\right)\)
\[\left|\psi_{n=1}\right|^{2}=\frac{2}{L} \int_{\frac{99 L}{200}}^{\frac{101L}{200}} \sin ^{2}\left(\frac{\pi x}{L}\right) dx \label{Ex3}\]
We can solve this, but we can also recognize that Equation \ref{Ex3} is just calculating an area that can be approximated as the area of a rectangle with a height (\(\frac{2}{L} \sin ^{2}\left(\frac{\pi x}{L}\right)\) at \(x=L/2\)) and width (\(\Delta x = L/100\)):
\[ \begin{align*} \left|\psi_{n=1}\right|^{2} &\approx \frac{2}{L} \cancelto{1}{\sin ^{2}\left(\frac{\pi (L/2)}{L}\right)} \Delta x \\[4pt] & \approx \left(\frac{2}{L} \right) (L/100) \\[4pt] & \approx 1/50 = 0.02\end{align*}\]
Given that the wavefunction is sinusoidal, the actual probability of finding an electron within the given interval at \(\frac{L}{2}\) should be slightly less because of the behavior of the sinusoid at \(\frac{L}{2}\) is at its peak of the wavefunction.
The \(n = 2\) state
For the \(n=2\) state, \(\psi_{n=2}=\sqrt{\frac{2}{2}} \sin \left(\frac{2 \pi x}{2}\right)\)
\[\left|\psi_{n=2}\right|^{2}=\frac{2}{L} \int_{\frac{99 L}{200}}^{\frac{101L}{200}} \sin ^{2}\left(\frac{2 \pi x}{L}\right) dx \label{Ex4}\]
We can use the same approximation as above
\[ \begin{align*} \left|\psi_{n=1}\right|^{2} &\approx \frac{2}{L} \cancelto{0}{\sin ^{2}\left(\frac{ 2\pi (L/2)}{L}\right)} \Delta x \\[4pt] & \approx 0\end{align*}\]
The probability of finding an electron in a box at \(\frac{L}{2}\) for \(n=2\) is approximately zero.
Time Dependence to the Wavefunction
Now, it is important to demonstrate that if a wavefunction is initially normalized then it stays normalized as it evolves in time according to the time-dependent Schrödinger's equation. If this is not the case then the probability interpretation of the wavefunction is untenable, since it does not make sense for the probability that a measurement of \(x\) yields any possible outcome (which is, manifestly, unity) to change in time. Hence, we require that
\[\dfrac{d}{dt} \int_{-\infty}^{\infty}\vert \psi(x,t)\vert^{ 2} dx = 0 \label{3.6.10}\]
for wavefunctions satisfying the time-dependent Schrödinger's equation (this results from the time-dependent Schrödinger's equation and Equation \(\ref{3.6.3}\)). The above equation gives
\[\dfrac{d}{dt} \int_{-\infty}^{\infty}\psi^* \psi dx=\int_{-\infty}^{\infty} \left(\dfrac{\partial\psi^*}{\partial t} \psi+\psi^* \dfrac{\partial\psi}{\partial t}\right) dx=0. \label{3.6.11}\]
Now, multiplying Schrödinger's equation by \(\psi^{*}/({\rm i} \hbar)\), we obtain
\[\psi^* \dfrac{\partial\psi}{\partial t}= \dfrac{\rm{i}}{2m} \psi^* \dfrac{\partial^2 \psi}{\partial x^2} - \dfrac{\rm{i}}{\hbar} V \vert \psi \vert^2 \label{3.6.12}\]
The complex conjugate of this expression yields
\[\psi \dfrac{\partial\psi^*}{\partial t}= -\dfrac{\rm{i}}{2m} \psi \dfrac{\partial^2 \psi^*}{\partial x^2} + \dfrac{\rm{i}}{\hbar} V \vert \psi \vert^2 \label{3.6.13}\]
since
- \((A B)^* = A^* B^*\),
- \(A^{* *}=A\), and
- \(i^*= -i\).
Summing Equation \ref{3.6.12} and \ref{3.6.13} results in
\[\dfrac{\partial\psi^*}{\partial t} \psi + \psi^{*} \dfrac{\partial\psi}{\partial t} = \dfrac{\rm{i}}{2m} \left(\psi^* \dfrac{\partial^2 \psi}{\partial x^2} - \psi \dfrac{\partial^2\psi^*}{\partial x^2}\right)=\dfrac{\rm{i}}{2m} \dfrac{\partial}{\partial x} \left(\psi^* \dfrac{\partial \psi}{\partial x} - \psi \dfrac{\partial\psi^*}{\partial x}\right). \label{3.6.14}\]
Equations \(\ref{3.6.11}\) and \(\ref{3.6.14}\) can be combined to produce
\[\dfrac{d}{dt}\int_{-\infty}^{\infty} \vert\psi\vert^2 dx= \dfrac{\rm{i}}{2m} \left[\psi^*\dfrac{\partial\psi}{\partial x}- \psi\dfrac{\partial\psi^*}{\partial x} \right]_{-\infty}^{\infty} = 0. \label{3.6.15}\]
The above equation is satisfied provided the wavefunction converges
\[\lim_{\vert x\vert\rightarrow\infty} \vert\psi\vert = 0 \label{3.6.16}\]
However, this is a necessary condition for the integral on the left-hand side of Equation \(\ref{3.6.3}\) to converge. Hence, we conclude that all wavefunctions which are square-integrable [i.e., are such that the integral in Equation \ref{3.6.3} converges] have the property that if the normalization condition Equation \(\ref{3.6.3}\) is satisfied at one instant in time then it is satisfied at all subsequent times.
Not all wavefunctions can be normalized according to the scheme set out in Equation \(\ref{3.6.3}\). For instance, a planewave wavefunction for a quantum free particle
\[\psi(x,t) = \psi_0 {\rm e}^{ {\rm i} (k x-\omega t)} \nonumber \]
is not square-integrable, and, thus, cannot be normalized. For such wavefunctions, the best we can say is that
\[P_{x \in a:b}(t) \propto \int_{a}^{b}\vert\psi(x,t)\vert^{ 2} dx.\nonumber \]
In the following, all wavefunctions are assumed to be square-integrable and normalized, unless otherwise stated.
Contributors
David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski ("Quantum States of Atoms and Molecules")
Richard Fitzpatrick (Professor of Physics, The University of Texas at Austin)