# 5.3: Linear Combinations of Eigenfunctions

- Page ID
- 4500

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

It is not necessary that an electron be described by an eigenfunction of the Hamiltonian operator. Many problems encountered by quantum chemists and computational chemists lead to wavefunctions that are not eigenfunctions of the Hamiltonian operator. Science is like that; interesting problems are not simple to solve. They require adaptation of current techniques, creative energy, and a good set of skills developed by studying solutions to previously solved interesting problems.

Consider a free electron in one dimension that is described by the wavefunction

\[ \psi (x) = C_1\psi _1 (x) + C_2 \psi _2 (x) \label {5-21}\]

with

\[ \begin{align} \psi _1(x) &= \left ( \dfrac {1}{2L} \right )^{1/2} e^{ik_1x} \label {5-22} \\[4pt] \psi _2(x) &= \left ( \dfrac {1}{2L} \right )^{1/2} e^{ik_2x} \label {5-23} \end{align} \]

where \(k_1\) and \(k_2\) have different magnitudes. Although such a function is not an eigenfunction of the momentum operator or the Hamiltonian operator, we can calculate the average momentum and average energy of an electron in this state from the expectation value integral. (Note: "in-this-state" means "described-by-this-wavefunction".)

Show that the function \(\psi(x)\) defined by Equation \(\ref{5-21}\) is not an eigenfunction of the momentum operator or the Hamiltonian operator for a free electron in one dimension.

The function shown in Equation \(\ref{5-21}\) belongs to a class of functions known as **superposition functions**, which are linear combinations of eigenfunctions. A linear combination of functions is a sum of functions, each multiplied by a weighting coefficient, which is a constant. The adjective linear is used because the coefficients are constants. The constants, e.g. \(C_1\) and \(C_2\) in Equation \(\ref{5-21}\), give the weight of each component (\(\psi_1\) and \(\psi_2\)) in the total wavefunction. Notice from the discussion previously that each component in Equation \(\ref{5-21}\) is an eigenfunction of the momentum operator and the Hamiltonian operator although the linear combination function (i.e., \(\psi(x)\)) is **not**.

The expectation value, i.e. average value, of the momentum operator is found as follows. First, write the integral for the expectation value and then substitute into this integral the superposition function and its complex conjugate as shown below. Since we are considering a free particle in one dimension, the limits on the integration are \(–L\) and \(+L\) with \(L\) going to infinity.

\[ \begin{align} \left \langle p \right \rangle &= \int \psi ^* (x) \left ( -i\hbar \dfrac {d}{dx} \right ) \psi (x) dx \\[4pt]&= \dfrac {-i\hbar}{2L} \int \limits _{-L}^{+L} \left ( C_1^* e^{-ik_1x} + C_2^* e^{-ik_2x} \right ) \dfrac {d}{dx} \left ( C_1 e^{ik_1x} + C_2 e^{ik_2x} \right ) dx \\[4pt] &= \dfrac {-i\hbar}{2L} \int \limits _{-L}^{+L} \left ( C_1^* e^{-ik_1x} + C_2^* e^{-ik_2x} \right ) \left ( (ik_1)C_1 e^{ik_1x} + (ik_2)C_2 e^{ik_2x} \right ) dx \label {5-24} \end{align} \]

Cross-multiplying the two factors in parentheses yields four terms.

\[\left \langle p \right \rangle = I_1 + I_2 + I_3 + I_4 \nonumber \]

with

\[ \begin{align} I_1 &= \dfrac {\hbar k_1}{2L} C^*_1 C_1 \int \limits ^{+L} _{-L} dx = C^*_1 C_1 \hbar k_1 \\[4pt] I_2 &= \dfrac {\hbar k_2}{2L} C^*_2 C_2 \int \limits ^{+L} _{-L} dx = C^*_2 C_2 \hbar k_2 \\[4pt] I_3 &= \dfrac {\hbar k_1}{2L} C^*_1 C_2 \int \limits ^{+L} _{-L} e^{i(k_2 - k_1)x} dx \\[4pt] I_4 &= \dfrac {\hbar k_1}{2L} C^*_2 C_1 \int \limits ^{+L} _{-L} e^{i(k_1 - k_2)x} dx \label {5-25} \end{align} \]

An integral of two different functions, e.g. \(\int \psi _1^* \psi _2 dx\), is called an **overlap integral** or **orthogonality integral**. When such an integral equals zero, the functions are said to be orthogonal. The integrals in \(I_3\) and \(I_4\) are zero because the functions \(\psi_1\) and \(\psi_2\) are orthogonal. We know \(\psi_1\) and \(\psi_2\) are orthogonal because of the Orthogonality Theorem, described previously, that states that eigenfunctions of any Hermitian operator, such as the momentum operator or the Hamiltonian operator, with different eigenvalues, which is the case here, are orthogonal. Also, by using Euler's formula and following Example \(\PageIndex{1}\) below, you can see why these integrals are zero.

For the integral part of \(I_3\) obtain

\[ \int \cos [(k_2 - k_1 ) x ] dx + i \int \sin [(k_2 - k_1)x ] dx \nonumber\]

from Euler’s formula.

**Solution**

Here we have the integrals of a cosine and a sine function along the x-axis from minus infinity to plus infinity. Since these integrals are the area under the cosine and sine curves, they must be zero because the positive lobes are canceled by the negatives lobes when the integration is carried out from \(–∞\) to \(+∞\).

As a result of this orthogonality, \(\left \langle p \right \rangle\) is just \(I_1 + I_2\), which is

\[ \begin{align} \left \langle p \right \rangle &= C_1^* C_1 \hbar k_1 + C^*_2 C_2 \hbar k_2 \\[4pt] &= C_1^* C_1p_1 + C^*_2 C_2 p_2 \label {5-26} \end{align} \]

where \(\hbar k_1\) is the momentum \(p_1\) of state \(\psi_1\), and \(\hbar k_2\) is the momentum \(p_2\) of state \(\psi_2\). As explained in Chapter 3, an average value can be calculated by summing, over all possibilities, the possible values times the probability of each value. Equation \(\ref{5-26}\) has this form if we interpret \(C_1^*C_1\) and \(C_2^*C_2\) as the probability that the electron has momentum \(p_1\) and \(p_2\), respectively. These coefficients therefore are called probability amplitude coefficients, and their absolute value squared gives the probability that the electron is described by \(\psi_1\) and \(\psi_2\), respectively. This interpretation of these coefficients as probability amplitudes is very important.

Find the expectation value for the energy \(\left \langle E \right \rangle\) for the superposition wavefunction given by Equation \(\ref{5-23}\).

Explain why \(C_1^*C_1\) is the probability that the electron has energy \(\dfrac {\hbar ^2 k^2_1}{2m}\) and \(C_2^*C_2\) is the probability that the electron has energy \(\dfrac {\hbar ^2 k^2_2}{2m}\).

What is the expectation value for the energy when both components have equal weights in the superposition function, i.e. when \(C_1 = C_2 = 2^{-1/2}\)?