On the Use of Perturbation Theory in Eigenvalue Problems

Show more

1. Introduction

The simplest eigenvalue problem in Algebra is that of finding the value of *λ* which makes the square matrix *A* diagonal, and this is solved by writing the characteristic equation which is

$\mathrm{det}\left(A-\lambda I\right)=0$ (1)

And solving for its roots, *I* is the unit matrix. However, the problem becomes more interesting when dealing with symmetric matrices since those are of importance in applied fields such as in sound in determining natural modes of vibrations or in elasticity. Such problems will stimulate using different computer software such as subroutine JACOBI [1].

Another kind of eigenvalue problems is that encountered in integral equations, e.g. to solve the eigenvalue problem [2]

${f}_{n}\left(x\right)={\lambda}_{n}{\displaystyle {\int}_{a}^{b}\text{d}y\text{\hspace{0.05em}}\text{\hspace{0.05em}}K\left(x,y\right){f}_{n}\left(y\right)}$ (2)

Most important are the eigenvalue problems, faced with, in quantum mechanics, *i.e.* when the time-dependent Schrodinger equation or the time-independent one are under consideration and the energy values are required for a certain system (e.g. the energy levels of the Hydrogen atom). The time-independent Schrodinger equation, as an eigenvalue problem, is written as

$H\phi =E\phi $ (3)

where
$H$ is the Hamiltonian, *E* is the energy eigenvalue, and *φ* is the eigenfunction. Note that
$H$ is given by

$H=-\frac{{\hslash}^{2}}{2m}{\nabla}^{2}+V\left(\stackrel{\to}{r}\right)$ (4)

${\nabla}^{2}$ is the Laplacian.

It is to be noted that dealing with eigenvalue problems, the operators encountered in the calculation are of the Hermitian type and that the Hamiltonian is Hermitian [3].

Keeping in mind the above small piece of information about the eigenvalue problems, we now turn to answer the important question: What is perturbation theory?

Perturbation theory is a theory built on certain mathematical methods used in order to get an approximate solution to some given problem with the condition that we start with an exact known solution for the related simpler problem. This means that the problem is partitioned into two parts: one is solvable and the other part represents the perturbation. Hence, perturbation theory is only used when the original problem cannot be solved exactly, but it can be formulated through the addition of a small perturbation term to the problem which can be solved exactly.

The ultimate solution to the original problem is obtained as a series called perturbation series. The first term of which is the exact solution of the related problem (without perturbation) and the other terms give the deviations in the exact solution; in fact, this gives us the approximate solution to the problem. It is clear that the more terms we take from the series, the better result we get.

In principle, the preliminary steps of performing perturbation techniques are as follows:

Let
$\u03f5$ be a small parameter and *S* is a certain solution to be approximated in the problem, then we arrive to a series of the type:

$S={S}_{0}+{\u03f5}^{1}{S}_{1}+{\u03f5}^{2}{S}_{2}+\cdots ={\displaystyle {\sum}_{i=0}^{\infty}{S}_{i}{\u03f5}^{i}}$ (5)

where ${S}_{0}$ is the solution to the unperturbed problem and the approximation in the original problem starts with the term ${\u03f5}^{1}{S}_{1}$. Hence, we get the approximate solution for the problem by truncating the above series; e.g. the approximate solution of the first order is

$S\cong {S}_{0}+{\u03f5}^{1}{S}_{1}$ (6)

In a way, perturbation theory is similar to the methods used in numerical analysis. The first early time, the theory was used was in space mechanics to evaluate the moon’s orbit; however, perturbation theory is now widely used in many fields, such as wave mechanics, statistical mechanics, quantum mechanics, and chemistry [4].

In quantum mechanics, perturbation theory is used to get approximate solutions for complex quantum systems in terms of simpler ones. The idea is to start with a simple system which is solvable then to add a weak perturbation Hamiltonian, and if the perturbation is small then the different physical phenomena (observables) which are connected with the perturbed system, such as energy levels and eingenfunctions, can be expressed as corrections to the corresponding simple system. Due to the smallness of these corrections they can be calculated via approximate methods in the form of converging series.

Perturbation theory is a good and important tool to describe real quantum systems noting that the exact solutions of Schrodinger equation for some Hamiltonians are difficult to obtain even if these Hamiltonians are not that complex [5].

In the next section, we give some examples on the use of perturbation theory in solving algebraic equations and in solving differential and integral equations. In Section 3, we direct our attention to its uses in solving eigenvalue problems in many applied fields. Finally we conclude shedding light on some recent uses of perturbation theory.

2. Algebraic Perturbation Theory

As mentioned, earlier, one of the important eigenvalue problems are of the type of Equation (2), *i.e.* it is an integral equation and where
$K\left(x,y\right)$ is the kernel; if this kernel is separable, which means that:

$K\left(x,y\right)=U\left(x\right)V\left(y\right)$ (7)

Then, one can get non-trivial solutions (eigenvalues and eigenfunctions).

Taking into account Equation (2) and Equation (7), we see that

${f}_{n}\left(x\right)={\lambda}_{n}AU\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}A\text{\hspace{0.17em}}\text{is}\text{\hspace{0.17em}}\text{a}\text{\hspace{0.17em}}\text{constant}$ (8)

Substituting, the eigenvalue will be given by

${\lambda}_{n}=\frac{1}{{\displaystyle {\int}_{a}^{b}U\left(y\right)V\left(y\right)\text{d}y}}$ (9)

Now if $U\left(y\right)=V\left(y\right)=\mathrm{sin}y$, then $f\left(x\right)=\lambda A\mathrm{sin}x$ and the eigenvalue is given by $\lambda =\frac{4}{\pi}$ [2].

2.1. Small Perturbation Parameter

If we consider the eigenvalue problem:

${\int}_{-\infty}^{\infty}{\text{e}}^{-\left|x\right|-\left|y\right|}U\left(y\right)\text{d}y}=\lambda U\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty <x<\infty $ (11)

And if we are looking for a value of *λ*giving a solution
$U\left(x\right)\ne 0$ for this equation fulfilling the condition that
${\int}_{-\infty}^{\infty}{\left|U\left(x\right)\right|}^{2}}<\infty $ ; then
$U\left(x\right)$ is an eigenfunction corresponding to the eigenvalue *λ*. In this case the spectrum of the eigenvalues is discrete.

For the eigenvalue *λ*, it is clear, with a straight forward manipulation that
$\lambda =1$ and the corresponding eigenfunction is
${\phi}_{0}={\text{e}}^{-\left|x\right|}$ [6].

Now, consider the following eigenvalue problem

${\int}_{-\infty}^{\infty}{\text{e}}^{-\left|x\right|-\left|y\right|}U\left(y\right)\text{d}y}+\u03f5xU\left(x\right)=\lambda U\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty <x<\infty $ (11)

where $\u03f5$ is a given real number. If $\u03f5$ is assumed to be of small value, then the problem in this case is to find the eigenvalue $\lambda =\lambda \left(\u03f5\right)$ and the corresponding eigenfunction $U=U\left(x,\u03f5\right)$. To solve the problem, we write the two following series

$\lambda ={\lambda}_{0}+\u03f5{\lambda}_{1}+{\u03f5}^{2}{\lambda}_{2}+\cdots ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\lambda}_{0}=1$ (12)

And

$U={\phi}_{0}+\u03f5{\phi}_{1}+{\u03f5}^{2}{\phi}_{2}+\cdots ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\phi}_{0}={\text{e}}^{-\left|x\right|}$ (13)

Moreover we put

$AU\equiv {\displaystyle {\int}_{-\infty}^{\infty}{\text{e}}^{-\left|x\right|-\left|y\right|}U\left(y\right)\text{d}y};\text{\hspace{0.17em}}\text{\hspace{0.17em}}BU\equiv xU\left(x\right)$ (14)

So as to get

$\left(A+\u03f5B\right)\left({\phi}_{0}+\u03f5{\phi}_{1}+{\u03f5}^{2}{\phi}_{2}+\cdots \right)=\left({\lambda}_{0}+\u03f5{\lambda}_{1}+{\u03f5}^{2}{\lambda}_{2}+\cdots \right)\left({\phi}_{0}+\u03f5{\phi}_{1}+{\u03f5}^{2}{\phi}_{2}+\cdots \right)$ (15)

And if we assume that the series ${\phi}_{0}+\u03f5{\phi}_{1}+{\u03f5}^{2}{\phi}_{2}+\cdots $ is convergent in order to be able to perform multiplication term by term, then on comparing coefficients of ${\u03f5}^{i}$ we get

$A{\phi}_{0}-{\lambda}_{0}{\phi}_{0}=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}A{\phi}_{1}-{\lambda}_{0}{\phi}_{1}={\lambda}_{1}{\phi}_{0}-B{\phi}_{0}$ (16)

In general for the nth terms, one gets

$A{\phi}_{n}-{\lambda}_{0}{\phi}_{n}={\lambda}_{n}{\phi}_{0}+{\lambda}_{n-1}{\phi}_{1}+\cdots +{\lambda}_{1}{\phi}_{n-1}-B{\phi}_{n-1}$ (17)

From Equation (16) we see that $\lambda =1$ and ${\phi}_{0}={\text{e}}^{-\left|x\right|}$, as expected.

The last equation is written in the form

$A{\phi}_{n}-{\lambda}_{0}{\phi}_{n}=f$, (18)

and where

$f\equiv {\lambda}_{n}{\phi}_{0}+{\lambda}_{n-1}{\phi}_{1}+\cdots +{\lambda}_{1}{\phi}_{n-1}-B{\phi}_{n-1}$ (19)

Taking the inner product of Equation (18) with
${\phi}_{0}$ and noting that the kernel is symmetric in *x* and *y*, one gets

${\lambda}_{1}=\left({\phi}_{0},B{\phi}_{0}\right)$ (20)

${\lambda}_{1}$ is the first order correction to the eigenvalue. Equation (20) can be stated as “The approximation of first order in the eigenvalue is equal to the mean value of the perturbation operator *B* with respect to the normalized unperturbed eigenfunction
${\phi}_{0}$ ” [6] [7].

In this problem ${\lambda}_{1}={\displaystyle {\int}_{-\infty}^{\infty}x{\text{e}}^{-2\left|x\right|}\text{d}x}=0$. Moreover, the correction to the eigenfunction ${\phi}_{1}$ is given by ${\phi}_{1}=x{\text{e}}^{-\left|x\right|}$.

Calculating
${\lambda}_{2}$ and
${\phi}_{2}$ in the same manner, one obtains
$\lambda $ and *U* as

$\lambda =1+\frac{1}{2}{\u03f5}^{2}+\dots $ (21)

And

$U={\text{e}}^{-\left|x\right|}\left(1+x\u03f5+\left({x}^{2}-\frac{3}{4}\right){\u03f5}^{2}+\cdots \right)$ (22)

2.2. Algebraic Equations

To illustrate the use of perturbation theory in solving certain algebraic equations, we give the following examples:

Example 1

Consider the cubic equation ${x}^{3}-x+\u03f5=0$, where $\u03f5$ is a small perturbation number and assume that the solution is of the form

$x={x}_{0}+\u03f5{x}_{1}+{\u03f5}^{2}{x}_{2}+O(\; \u03f5\; 3\; )$

where $O\left({\u03f5}^{3}\right)$ is a term containing all terms with ${\u03f5}^{3}$ and higher. Substituting in the cubic equation and equating coefficients of ${\u03f5}^{i}$, we get

${x}_{0}^{3}-{x}_{0}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}3{x}_{0}^{2}{x}_{1}-{x}_{1}+1=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}3{x}_{0}{x}_{2}-{x}_{2}+3{x}_{0}{x}_{1}^{2}=0$

From which, we obtain ${x}_{0}=0,\pm 1$ and ${x}_{1}=\frac{1}{1-3{x}_{0}^{2}}$. Therefore, for ${x}_{0}=0$ we get ${x}_{1}=1$ and $x=\u03f5+O\left({\u03f5}^{2}\right)$ ; while for ${x}_{0}=\pm 1$ we get $x=\pm 1-\frac{1}{2}\u03f5+O\left({\u03f5}^{2}\right)$.

Continuing this process we can have the solution as convergent expansion in $\u03f5$. The truncation in the series is determined by the degree of precision we like to reach. Moreover, the same result can be obtained using the usual Taylor expansion [7].

Example 2

Again, we consider the cubic equation
$\u03f5{x}^{3}-x+1=0$ ; using this with the familiar expansion in
$\u03f5$ and on substitution, we get
$-{x}_{0}+1=0$,
$-{x}_{1}+{x}_{0}^{3}=0$,
$-{x}_{2}+3{x}_{1}{x}_{0}^{2}=0$ ; for
${x}_{0}={x}_{1}=1$, we have
$x\left(\u03f5\right)=1+\u03f5+O\left({\u03f5}^{2}\right)$, this gives the linear part of the solution. As to the two other solutions, we note that as
$\u03f5\to 0$,
$x\left(\u03f5\right)\to \infty $. This is actually a standard singular perturbation problem. To overcome this problem, some sort of scaling is needed; this is done by writing *x* as

$x\left(\u03f5\right)=\frac{1}{\rho \left(\u03f5\right)}y\left(\u03f5\right)$ so as $y=O\left(1\right)$ as $\u03f5\to 0$, and from the cubic equation one gets $\frac{\u03f5{y}^{3}}{{\rho}^{3}}-\frac{y}{\rho}+1=0$. $\rho $ needs to be determined in order to fulfill our

requirement, namely to get a definite answer. For that, the principle of dominant balance is to be used; this requires that the two leading terms in the equation, at

least, should have the same order of magnitude. Hence
$\frac{\u03f5}{{\rho}^{3}}=\frac{1}{\rho}$ giving
$\rho =\sqrt{\u03f5}$. Our cubic equation, now, becomes, in terms of the variable *y*,
${y}^{3}-y+\sqrt{\u03f5}=0$. Following the steps of the last example, we get
$y=\pm 1-\frac{1}{2}{\u03f5}^{\frac{1}{2}}+O\left(\u03f5\right)$ or
$x=\pm \frac{1}{\sqrt{\u03f5}}-\frac{1}{2}+O\left({\u03f5}^{1/2}\right)$ ; *which are the other two solutions of the problem* [7].

Example 3

If the quadratic equation $\left(1-\u03f5\right){x}^{2}-2x+1=0$ is to be solved, then making the usual expansion in $\u03f5$ and comparing coefficients of ${\u03f5}^{i}$, we get ${x}_{0}^{2}-2{x}_{0}+1=0$, $2\left({x}_{0}-1\right){x}_{1}={x}_{0}^{2}$ ; the first equation leads to the solution ${x}_{0}=1$ which is a repeated root, substituting in the second equation we get a contradiction (0 = 1?) which means that the series techniques do not work here. Hence, we calculate the roots of the equation using the simple general square root method to get

${x}_{\pm}=\frac{1\pm \sqrt{\u03f5}}{1-\u03f5}$ ; that gives us a hint for writing the expansion as

$x={x}_{0}+\sqrt{\u03f5}{x}_{1}+\u03f5{x}_{2}+O\left({\u03f5}^{\frac{3}{2}}\right)$ instead of the expansion $x={x}_{0}+\u03f5{x}_{1}+O\left({\u03f5}^{2}\right)$. Now

following the same techniques used before we get ${x}_{0}=1$, ${x}_{1}=\pm 1$ and hence the solutions are $x=1\pm \sqrt{\u03f5}+O\left(\u03f5\right)$ as expected [7].

Example 4

Given that $x{\text{e}}^{-x}=\u03f5$, then we see that there are two possible solutions for it for the case $\u03f5\to {0}^{+}$, namely

1) When $x\to 0$ and $x=\u03f5+{\u03f5}^{2}+O\left({\u03f5}^{3}\right)$.

2) When
${\text{e}}^{-x}\to 0$ as
$x\to \infty $, for this case *x* has to behave like
$\mathrm{log}\left(1/\u03f5\right)$ ; and in order to get a convergent solution, we write the problem as
${\text{e}}^{-x}=\u03f5/x$ which implies that
$x=\mathrm{log}x+\mathrm{log}\left(1/\u03f5\right)$, note that
${x}_{1}=\mathrm{log}\left(1/\u03f5\right)$ is the initial starting point for any recurrence techniques we choose.

Using the formula ${x}_{n+1}=\mathrm{log}{x}_{n}+\mathrm{log}\left(1/\u03f5\right),n\in N$ and defining ${L}_{1}=\mathrm{log}\left(1/\u03f5\right)$, ${L}_{2}=\mathrm{log}\left(\mathrm{log}\left(1/\u03f5\right)\right)$, we see that ${x}_{2}={L}_{1}+{L}_{2}$, ${x}_{3}={L}_{1}+\mathrm{log}\left({L}_{1}+{L}_{2}\right)$, etc.

Note that the successive terms of the series converge slowly as $\u03f5\to 0$ [7].

2.3. Examples from Different Mathematics Fields

In this subsection, we give applications of perturbation theory in other different fields of mathematics such as differential equations and geometry. These are illustrated by few examples.

Example 5

Given the ordinary differential equation $\frac{{\text{d}}^{2}y}{\text{d}{x}^{2}}=-\u03f5\frac{\text{d}y}{\text{d}x}-1$, $y\left(0\right)=0$, $\stackrel{\dot{}}{y}\left(0\right)=1$,

which represents a projectile vertical motion in air with air resistance force taken into account.
$\u03f5=kv/mg$, where *v* is the initial velocity, *m* is the projectile mass, *k* is the friction constant, and *g* is the gravitational constant. If we want now to

solve for *y*, we get the exact solution as
$y\left(x\right)=\frac{1+\u03f5}{{\u03f5}^{2}}\left(1-{\text{e}}^{-\u03f5x}\right)-x/\u03f5$ ; and if

$0<\u03f5\le 1$ which corresponds to the case when the air resistance is weak, then we use Taylor’s series to obtain the first few terms of the expansion as

$y\left(x\right)=x-\frac{{x}^{2}}{2}+\u03f5\left(-\frac{{x}^{2}}{2}+\frac{{x}^{3}}{6}\right)+{\u03f5}^{2}\left(\frac{{x}^{3}}{6}-\frac{{x}^{4}}{24}\right)+O\left({\u03f5}^{3}\right)$. The expansion is valid for $0<x\ll 1/\u03f5$.

Now, assuming that we are not aware of the exact solution, then we use perturbation theory to write $y\left(x\right)={y}_{0}\left(x\right)+\u03f5{y}_{1}\left(x\right)+{\u03f5}^{2}{y}_{3}\left(x\right)+O\left({\u03f5}^{3}\right)$ ; and where ${y}_{i}\left(x\right),i=0,1,2,\cdots $ are to be determined; as before, we substitute this in the original equation, following the same old techniques, keeping in mind that the initial conditions become ${y}_{0}\left(0\right)+\u03f5{y}_{1}\left(0\right)+{\u03f5}^{2}{y}_{2}\left(0\right)+O\left({\u03f5}^{3}\right)=0$, ${\stackrel{\dot{}}{y}}_{0}\left(0\right)-1+\u03f5{\stackrel{\dot{}}{y}}_{1}+{\u03f5}^{2}{\stackrel{\dot{}}{y}}_{2}+O\left({\u03f5}^{3}\right)=0$ ; then we get the system of differential equations ${\stackrel{\xa8}{y}}_{0}+1=0,{y}_{0}\left(0\right)=0,{\stackrel{\dot{}}{y}}_{0}\left(0\right)=1$ ; ${\stackrel{\xa8}{y}}_{1}+{\stackrel{\dot{}}{y}}_{0}=0,{y}_{1}\left(0\right)=0,{\stackrel{\dot{}}{y}}_{1}=0$ ;

${\stackrel{\xa8}{y}}_{2}+{\stackrel{\dot{}}{y}}_{1}=0,{y}_{2}\left(0\right)=0,{\stackrel{\dot{}}{y}}_{2}\left(0\right)=0$, on equating the coefficients of ${\u03f5}^{i}$. Solving these

equations, we get ${y}_{0}\left(x\right)=x-\frac{{x}^{2}}{2}$, ${y}_{1}\left(x\right)=-\frac{{x}^{2}}{2}+\frac{{x}^{3}}{6}$, ${y}_{2}\left(x\right)=\frac{{x}^{3}}{6}-\frac{{x}^{4}}{24}$, in conformation with the results obtained using Taylor’s series expansion [8].

Example 6

Consider a right angle triangle with sides 1, 0.1, and $1.0049\cong \sqrt{{1}^{2}+{0.1}^{2}}$, it is difficult to differentiate between this right angle triangle and its being an isosceles one; since there is not that difference between the two sides 1 and 1.0049. Let us make a simulation of the situation using perturbation theory by considering the side as 1, $\u03f5$ and $\sqrt{1+{\u03f5}^{2}}$. Using binomial theorem we have

$\sqrt{1+{\u03f5}^{2}}=1+{\u03f5}^{2}+O\left({\u03f5}^{4}\right)$. Hence, the difference between the side 1 and the hypotenuse is of order ${\u03f5}^{2}\left(=0.01\right)$ which is a small number. Note that the angle between

them is given by $\theta ={\mathrm{tan}}^{-1}\left(\frac{\u03f5}{1}\right)\cong \u03f5=5.71\u02da$ which is a sensible angle and can be

observed better than the difference between the side and the hypotenuse. Actually, if we make the value of $\u03f5$ smaller, we will be misled to the conclusion that the triangle is segment of a straight line [9].

2.4. The Smallness of the Perturbation Parameter and Its Relation with the Unperturbed Operator

If we consider the following eigenvalue problem

$-\frac{{\text{d}}^{2}u}{\text{d}{x}^{2}}+{a}^{2}{x}^{2}u=\lambda u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty <x<\infty $ (23)

where *a* is a given positive constant, then the discrete eigenvalues are given by

${\lambda}_{n}=\left(2n+1\right)a,\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,3,\cdots $ (24)

And the corresponding eigenfunctions are

${u}_{n}={\text{e}}^{-a{x}^{2}/2}{H}_{n}\left(\sqrt{a}x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,3,\cdots $ (25)

${H}_{n}\left(x\right)={\left(-1\right)}^{n}{\text{e}}^{{x}^{2}}\frac{{\text{d}}^{2}}{\text{d}{x}^{2}}\left({\text{e}}^{{x}^{2}}\right)$ is a Hermite polynomial.

This problem describes the eigenvalue problem for the simple harmonic oscillator [10].

Now, adding a perturbation term of the form $\u03f5xu$, we get

$-\frac{{\text{d}}^{2}u}{\text{d}{x}^{2}}+{a}^{2}{x}^{2}u+\u03f5xu=\lambda u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\infty <x<\infty $ (26)

We may look at the term
$\u03f5xu$ as not small. This is because the variable value *x* takes values from −∞ to ∞; however if we compare this perturbation parameter with the term
${a}^{2}{x}^{2}u$ which is the non-perturbed one, the perturbation parameter (for small
$\u03f5$ ) is then considered as a small perturbation. One concludes then that: In perturbation problems, the perturbing parameter has to be small in comparison with the unperturbed operator [6].

Now the eigenvalues and eigenfunctions for the perturbed problem are given by

${\lambda}_{n}=\left(2n+1\right)a-\frac{{\u03f5}^{2}}{4{a}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,3,\cdots $ (27a)

and

${u}_{n}={\text{e}}^{-a/2{\left(x+\frac{\u03f5}{2{a}^{2}}\right)}^{2}}{H}_{n}\left(\sqrt{a}\left[x+\frac{\u03f5}{2{a}^{2}}\right]\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,3,\cdots $ (27b)

Note that the functions ${u}_{n}$ are regular functions in $\u03f5$.

At the end of this subsection, we note that the main goal of perturbation theory is the precise definition of smallness of the perturbing parameter with respect to the unperturbed operator [6].

Before we move to the applications of perturbation theory in different applied fields, we make a little digression, in the next section, on spectral perturbation theory.

2.5. Spectral Perturbation Theory

The theory deals with the way the perturbation affects the spectrum of a certain operator when it is perturbed. This compels us to reformulate the problem in a more precise and different manner [6].

Definition 1

If $\mathbb{H}$ is a Hilbert space with inner product $\langle .,.\rangle $ and that ${A}^{\u03f5}$ is a linear operator in $\mathbb{H}$ such that $D\left({A}^{\u03f5}\right)\subset \mathbb{H}\to \mathbb{H}$ and $D\left({A}^{\u03f5}\right)$ is the domain of ${A}^{\u03f5}$ which depends on a small parameter $\u03f5$ ; assume that

1)
${A}^{\u03f5}$ Is a self-adjoint operator, *i.e.*

$\langle x,{A}^{\u03f5}y\rangle =\langle {A}^{\u03f5}x,y\rangle ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall x,y\in D\left({A}^{\u03f5}\right)$ (28)

2) ${A}^{\u03f5}$ has a smooth branch of simple eigenvalues ${\lambda}^{\u03f5}\in R$ with eigenvectors ${x}^{\u03f5}\in \mathbb{H}$ ; namely

${A}^{\u03f5}{x}^{\u03f5}={\lambda}^{\u03f5}{x}^{\u03f5}$ (29)

Then, the job of perturbation theory is the calculation of the eigenvalues with respect to their values at $\u03f5=0$ when $\u03f5$ is small but not equal to zero [6].

To clarify this idea, we take as an example the perturbation in the eigenvalues for a symmetric matrix; in this case $\mathbb{H}={R}^{n}$, with the Euclidean inner product

$\langle x,y\rangle ={x}^{\text{T}}y$ (30)

And ${A}^{\u03f5}:{R}^{n}\to {R}^{n}$ is a linear transformation for a symmetric matrix of the order $n\times n$ ( ${a}_{ij}^{\u03f5}$ ). [note that we are concentrating on symmetric matrices because these are what we face in physics]

Moreover, the perturbation in the eigenvalues for a Hermitian matrix is of our concern, in this case $\mathbb{H}={C}^{n}$ with the inner product

$\langle x,y\rangle ={\stackrel{\xaf}{x}}^{\text{T}}y$ (31)

Problems of this sort are found in many applications especially in quantum mechanics and where Hilbert spaces of infinite dimensions are dealt with. In this case expansions of the form

${A}^{\u03f5}={A}_{0}+\u03f5{A}_{1}+\cdots +{\u03f5}^{n}{A}_{n}+\cdots $, (32)

${x}^{\u03f5}=x+\u03f5{x}_{1}+\cdots +{\u03f5}^{n}{x}_{n}+\cdots $, (33)

${\lambda}^{\u03f5}={\lambda}_{0}+\u03f5{\lambda}_{1}+\cdots +{\u03f5}^{n}{\lambda}_{n}+\cdots $ (34)

Are used.

Now following the same techniques we have used in solving previous examples and in particular equating coefficients of ${\u03f5}^{n}$, we get the following equations

$\left({A}_{0}-{\lambda}_{0}I\right){x}_{0}=0$ (35)

$\left({A}_{0}-{\lambda}_{0}I\right){x}_{1}=-{A}_{1}{x}_{0}+{\lambda}_{1}{x}_{0}$ (36)

$\cdots $

$\left({A}_{0}-{\lambda}_{0}I\right){x}_{n}={\displaystyle {\sum}_{i=1}^{n}\left(-{A}_{i}{x}_{n-i}+{\lambda}_{i}{x}_{n-i}\right)}$ (37)

[*I* is the
$n\times n$ unit matrix].

Now, suppose that ${x}_{0}\ne 0$, then from Equation (35) ${\lambda}_{0}$ is the eigenvalue of the unperturbed operator ${A}_{0}$ ; and the corresponding eigenvector is ${x}_{0}$ [7].

As to the second correction ${\lambda}_{1}$ and ${x}_{1}$, we note that Equation (36) is a singular equation and in order to be able to solve it, we need a certain hypothesis and the so-called Fredholm alternative [7].

3. Applications

In this section, we give some representative well-known problems in different branches of physics; we start first with the ones from classical physics, then we move to the ones from modern physics (or more specific in quantum mechanics).

3.1. Problems from Classical Physics

1) String Problem in the Theory of Sound

The vibrations of a string, with a small stiffness can be described the following eigenvalue problem

$-\frac{{\text{d}}^{2}u}{\text{d}{x}^{2}}+\u03f5\frac{{\text{d}}^{4}u}{\text{d}{x}^{4}}=\lambda u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<x<1$ (38)

where
$\u03f5$ is a small positive number and is proportional to the stiffness, *i.e.*
$\u03f5$ represents the perturbation. To solve this problem boundary conditions are needed. For instance to impose the conditions at the end points;
$x=1$ and
$x=0$ [6].

If the two ends of the string are to be supported, then the boundary conditions will be $u\left(0\right)=\stackrel{\xa8}{u}\left(0\right)=0$, $u\left(1\right)=\stackrel{\xa8}{u}\left(1\right)=0$ ; with these conditions in mind, one can get the eigenfunctions and eigenvalues as

${u}_{n}\left(x\right)=\mathrm{sin}n\pi x\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\lambda}_{n}={n}^{2}{\pi}^{2}+\u03f5{n}^{4}{\pi}^{4},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=1,2,3,\cdots $ (39)

The set ${\left\{\mathrm{sin}n\pi x\right\}}_{n=1}^{\infty}$ is complete in the interval $\left(0,1\right)$ ; this means that ${\left\{{\lambda}_{n}\right\}}_{n=1}^{\infty}$ represents the full spectrum of the problem. For $\u03f5\ge 0$, the spectrum is bounded below and the smallest value of the eigenvalues is ${\lambda}_{1}={\pi}^{2}+\u03f5{\pi}^{4}$. However, for $\u03f5<0$ the spectrum is not bounded below and this is not allowed because $\u03f5$ measure the stiffness which is a physical quantity and is not negative.

2) The Simple Pendulum

The equation of motion for the simple pendulum is

$mL\stackrel{\xa8}{\theta}+mg\mathrm{sin}\theta =0$ (40)

where
$\theta $ is the angle between the light string (of length *L*) and the vertical equilibrium position, *m* is the mass, and *g* is the gravitational acceleration. The equation can be written as

$\stackrel{\xa8}{\theta}+{k}^{2}\mathrm{sin}\theta =0$ (41)

where ${k}^{2}=g/L$.

Now, taking into account the initial conditions, the problem, for a small angle, then becomes

$\stackrel{\xa8}{\theta}+{k}^{2}\mathrm{sin}\theta =0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\theta \left(0\right)=\u03f5,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\dot{}}{\theta}\left(0\right)=0$ (42)

And $0<\u03f5\ll 1$, $\theta =O\left(\u03f5\right)$.

Putting $\theta =\u03f5\Psi $ and dividing by $\u03f5$, the last equation can be written as

$\stackrel{\xa8}{\Psi}+{k}^{2}\left(\Psi -\frac{1}{6}{\u03f5}^{2}{\Psi}^{3}\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\Psi \left(0\right)=1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\dot{}}{\Psi}=0$ (43)

If the solution is required to be of second order, then the higher order terms can be neglected. Making the transformation
$\tau =\omega t$, where *t* is the time, Equation (43) is then written as

${\omega}^{2}\stackrel{\xa8}{\Psi}+{k}^{2}\left(\Psi -\frac{1}{6}{\u03f5}^{2}{\Psi}^{3}\right)=0$ (44)

And where the differentiation is with respect to $\tau $. Now, expanding in terms of powers of ${\u03f5}^{2}$, we get

$\omega ={\omega}_{0}+{\u03f5}^{2}{\omega}_{1}$ (45)

And

$\Psi ={\Psi}_{0}+{\u03f5}^{2}{\Psi}_{1}$ (46)

Substituting with Equations (45)-(46) in Equation (44) and equating coefficients of ${\u03f5}^{2}$, one obtains

${\omega}_{0}{\stackrel{\xa8}{\Psi}}_{0}+{k}^{2}{\Psi}_{0}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Psi}_{0}=1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\dot{}}{\Psi}}_{0}=0$ (47)

And

${\omega}_{0}{\stackrel{\xa8}{\Psi}}_{1}+{k}^{2}{\Psi}_{1}=-2{\omega}_{0}{\omega}_{1}{\stackrel{\xa8}{\Psi}}_{0}+\frac{1}{6}{k}^{2}{\Psi}_{0}^{3};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Psi}_{1}\left(0\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\dot{}}{\Psi}}_{1}\left(0\right)=0$ (48)

Choosing ${\omega}_{0}=k$, the solution ${\Psi}_{0}$ is given by

${\Psi}_{0}=\mathrm{cos}\tau ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\omega}_{0}=k$, (50)

while ${\Psi}_{1}$ satisfies the equation

${\stackrel{\xa8}{\Psi}}_{1}+{\Psi}_{1}=\frac{2}{k}{\omega}_{1}\mathrm{cos}\tau +\frac{1}{8}\mathrm{cos}\tau +\frac{1}{24}\mathrm{cos}3\tau $ (51)

If we choose ${\omega}_{1}=-\frac{k}{16}$ the solution ${\Psi}_{1}$ will then be

${\Psi}_{1}=\frac{1}{192}\left(\mathrm{cos}\tau -\mathrm{cos}3\tau \right)$ (52)

and the final solution is

$\theta \left(t\right)=\u03f5\mathrm{cos}\tau +\frac{1}{192}{\u03f5}^{3}\left(\mathrm{cos}\tau -\mathrm{cos}3\tau \right)+O\left({\u03f5}^{5}\right)$ (53)

Moreover, $\tau =k\left[1-\frac{1}{16}{\u03f5}^{2}+O\left({\u03f5}^{4}\right)\right]$ [11].

3) Steady State Diffusion Through a Corrugated Slab

In this application, the simplest type of heat transfer is dealt with, which is the steady state type, *i.e.* when
$\frac{\partial T}{\partial t}=0$. The flow of heat is through a corrugated

slab of thickness *H*, with the assumption that the lower side (
$z=0$ ) is at temperature
${T}_{0}$ ; and the upper side (
$z=H$ ) is at temperature
${T}_{0}+\Delta T$. Then, for steady state flow we have

$k{T}_{zz}=0$ (54)

(this is because the temperature does not depend on time explicitly $\frac{\partial T}{\partial t}=0$ ).

From the boundary conditions, the solution is given by

$T={T}_{0}+\frac{\Delta T}{H}z$ (55)

And the heat flux through the slab is determined by

$\stackrel{\to}{F}=-\frac{k\Delta T}{H}\frac{\stackrel{\to}{z}}{z}$ (56)

Now, assuming that the slab occupies the region $-\infty <x<\infty $, $h\left(x\right)<z<H$, and suppose that the slab is not homogenous with a boundary condition on the lower part at $z=h\left(x\right)$ given by

$h\left(x\right)={h}_{\mathrm{max}}\mathrm{cos}kx$ (57)

where ${h}_{\mathrm{max}}<H$ ; in this case Laplace equation for the steady state case is given by

${T}_{xx}+{T}_{zz}=0$ (58)

with Dirichlet boundary conditions given as

$T\left(x,h\left(x\right)\right)={T}_{0};\text{\hspace{0.17em}}\text{\hspace{0.17em}}T\left(x,H\right)={T}_{0}+\Delta T$ (59)

Introducing the dimensionless variables

$\stackrel{\xaf}{T}=\frac{T-{T}_{0}}{\Delta T},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\xaf}{x}=kx,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\xaf}{z}=kz$ (60)

The problem becomes dimensionless and the boundary conditions will be as

$T\left(x,\u03f5\mathrm{cos}x\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}T\left(x,\beta \right)=1$ (61)

where $\u03f5=k{h}_{\mathrm{max}},\beta =kH$. It is clear that ${h}_{\mathrm{max}}<H$ and hence $\u03f5<\beta $.

For small amplitudes, $\u03f5\ll 1$, and this is fulfilled when ${h}_{\mathrm{max}}$ is made to be very small ( $\u03f5\to 0$ with $\beta $ kept as a constant) [12].

Now, the boundary condition $T\left(x,\u03f5h\right)=0$ can be written as a Taylor series as

$T\left(x,0\right)+\u03f5h{T}_{z}\left(x,0\right)+\frac{1}{2}{\left(\u03f5h\right)}^{2}{T}_{zz}\left(x,0\right)+\cdots $ (62)

Also, using perturbation theory, *T* is written as

$T={T}_{0}+\u03f5{T}_{1}+{\u03f5}^{2}{T}_{2}+\cdots $ (63)

Using the last two equations, following the same old procedure, and equating coefficients of ${\u03f5}^{n}$, we get

${T}_{0}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{T}_{1}+h{T}_{0z}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{T}_{2}+h{T}_{1z}+\frac{1}{2}{h}^{2}{T}_{0zz}=0$ (64)

This is for the first boundary condition; as for the second one $T\left(x,\beta \right)=1$, it tends to ${T}_{0}=1$ and ${T}_{n\ge 1}=0$, with all ${T}_{n}\left(x,z\right)$ satisfying Laplace equation

$\left(\frac{{\partial}^{2}}{\partial {x}^{2}}+\frac{{\partial}^{2}}{\partial {z}^{2}}\right){T}_{n}=0$. putting, now $h=\mathrm{cos}x$ ; we get the solution as

$\begin{array}{c}T=\frac{z}{\beta}-\u03f5\mathrm{cos}x\frac{\mathrm{sinh}\left(\beta -z\right)}{\beta \mathrm{sinh}h\beta}\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}-\frac{1}{2}{\u03f5}^{2}\frac{\mathrm{coth}\beta}{\beta}\left[\frac{\beta -z}{\beta}+\mathrm{cos}2x\left(\frac{\mathrm{sinh}2\left(\beta -z\right)}{\mathrm{sinh}2\beta}\right)\right]+O\left({\u03f5}^{3}\right)\end{array}$ (65)

4) Incompressible Flow Past a Circle

Here, the stream function is given by

$\text{d}\Psi =u\text{d}y-v\text{d}x$ (66)

$\left(u,v\right)$ is the velocity; and if the motion is inviscid, the function satisfies the equation

${\nabla}^{2}\Psi =-\omega \left(\Psi \right)$ (67)

where $\omega $ is the vorticity, which is equal to zero upstream if the flow is uniform; hence Equation (66) becomes, in cylindrical coordinates, as

${\Psi}_{rr}+\frac{1}{r}{\Psi}_{r}+\frac{1}{{r}^{2}}{\Psi}_{\theta \theta}=0$ (68)

Also, the boundary conditions are given by

$\Psi \left(r,\theta \right)\to ur\mathrm{sin}\theta $ as $r\to \infty $ (upstream) (69)

and

$\Psi \left(a,\theta \right)=0$ (on the surface) (70)

where *a* is the radius of the circle. Moreover, it is assumed that there is no rotational motion, this means that the flow should be symmetrical about the line
$\theta =0$ and this means that
$\Psi \left(r,\theta \right)=-\Psi \left(r,-\theta \right)$ ; hence the solution, in this case, will be

${\Psi}_{0}=u\left(r-\frac{{a}^{2}}{r}\right)\mathrm{sin}\theta $ (71)

Now, in the presence of the circle and if the vorticity number $\u03f5$ is small, then we expect that the solution will not depart from the one we got in Equation (70). So, we write the solution in the presence of vorticity as

$\Psi \left(r,\theta ;\u03f5\right)={\Psi}_{0}\left(r,\theta \right)+\u03f5{\Psi}_{1}\left(r,\theta \right)$ (72)

${\Psi}_{0}\left(r,\theta \right)$ is the unperturbed solution given by Equation (70).

On substituting Equation (71) in Equation (67), and equating coefficients of ${\u03f5}^{n}$, we get

$\frac{{\partial}^{2}{\Psi}_{1}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial {\Psi}_{1}}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial}^{2}{\Psi}_{1}}{\partial {\theta}^{2}}=\frac{u}{a}$ (73)

Also, we get the boundary conditions as

${\Psi}_{1}\to \frac{1}{4}\frac{u}{a}{r}^{2}\left(1-\mathrm{cos}2\theta \right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}r\to \infty ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Psi}_{1}\left(a,\theta \right)=0$ (74)

Therefore, the solution is written as

${\Psi}_{1}=\frac{1}{4}\frac{u}{a}{r}^{2}\left(1-\mathrm{cos}2\theta \right)+\phi \left(r,\theta \right)$ (75)

And the equation in $\phi \left(r,\theta \right)$ is given by

$\frac{{\partial}^{2}{\phi}_{1}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial {\phi}_{1}}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial}^{2}{\phi}_{1}}{\partial {\theta}^{2}}=0$ (76)

With the conditions

${\phi}_{1}\to \text{constant}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{as}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}r\to \infty ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\phi}_{1}\left(a,\theta \right)=-\frac{1}{4}ua\left(1-\mathrm{cos}2\theta \right)$ (77)

Using separation of variables, the solution in final form is

$\Psi =u\left(r-\frac{{a}^{2}}{r}\right)\mathrm{sin}\theta +\frac{1}{4}\u03f5uu\left[\frac{{r}^{2}}{a}\left(1-\mathrm{cos}2\theta \right)+\frac{{a}^{3}}{{r}^{2}}\mathrm{cos}2\theta -a\right]$ (78)

Equation (77) gives the solution in first order in $\u03f5$ [12].

3.2. Problems in Modern Physics

In some physics subjects, one can solve for the eigenvalues and eigenfunctions of the Hamiltonian of the problem exactly and precisely such as the simple harmonic oscillator and the Hydrogen atom in atomic physics. However, many physical systems in nature cannot be solved exactly and that requires the need for the development of appropriate procedures and tools to deal with such situations. As we stated before, perturbation theory was successful in treating these cases in which the distortion in the system, that can be solved exactly, is small.

Mathematically, the Hamiltonian of the system $H$ is written as

$H={H}_{0}+\u03f5V$ (79)

${H}_{0}$ is the unperturbed part of the Hamiltonian, $\u03f5$ is a small parameter, and $V$ is the potential, related to the perturbation, and supposed to be independent of time.

We assume that we can solve the time-independent Schrodinger equation for the unperturbed Hamiltonian *i.e.* we can get the eigenvalues
${E}_{n}$ and the eigenfunctions
${\Psi}_{n}\left(x\right)$, such that

${H}_{0}{\Psi}_{n}\left(x\right)={E}_{n}{\Psi}_{n}\left(x\right)$ (80)

Such solutions may give unperturbed eigenvalues ${E}_{n}$, nondegenerate or degenerate.

3.3. Nondegenerate Solutions

Let the time-independent Schrodinger equation for the total Hamiltonian be

$H\Psi \left(x\right)=E\Psi \left(x\right)$ (81)

Since $\u03f5$ is a small parameter, we can use Taylor series to expand the wave function $\Psi \left(x\right)$ in terms of $\u03f5$ as

$\Psi \left(x\right)={\Psi}_{0}\left(x\right)+\u03f5{\eta}_{1}\left(x\right)+{\u03f5}^{2}{\eta}_{2}\left(x\right)+\cdots $ (82)

The same thing can be done with the state energy to get

$E={E}_{0}+\u03f5{\mathcal{E}}_{1}+{\u03f5}^{2}{\mathcal{E}}_{2}+\cdots $ (83)

Substituting with these two equations in Equation (81) one gets

$\begin{array}{l}\left({H}_{0}+\u03f5V\right)\left({\Psi}_{0}\left(x\right)+\u03f5{\eta}_{1}\left(x\right)+{\u03f5}^{2}{\eta}_{2}\left(x\right)+\cdots \right)\\ =\left({E}_{0}+\u03f5{\mathcal{E}}_{1}+{\u03f5}^{2}{\mathcal{E}}_{2}+\cdots \right)\left({\Psi}_{0}\left(x\right)+\u03f5{\eta}_{1}\left(x\right)+{\u03f5}^{2}{\eta}_{2}\left(x\right)+\cdots \right)\end{array}$ (84)

[
$\left({\mathcal{E}}_{1},{\eta}_{1}\right),\left({\mathcal{E}}_{2},{\eta}_{2}\right),\cdots $ are the corrections in *E* and
$\Psi $ ].

Comparing terms of ${\u03f5}^{n}$, one gets

For the order ${\u03f5}^{0}$

$\left({H}_{0}-{E}_{0}\right){\Psi}_{0}=0$ (85)

This equation actually shows that ${\Psi}_{0}$ is the unperturbed eigenfunction corresponding to the eigenvalue ${E}_{0}$.

For the order ${\u03f5}^{1}$

$\left({H}_{0}-{E}_{0}\right){\eta}_{1}+\left(V-{\mathcal{E}}_{1}\right){\Psi}_{0}=0$ (86)

Taking the scalar product of this equation with ${\Psi}_{0}$, one gets

$\langle {\Psi}_{0},{H}_{0},{\eta}_{1}\rangle +\langle {\Psi}_{0},V{\Psi}_{0}\rangle ={E}_{0}\langle {\Psi}_{0},{\eta}_{1}\rangle +{\mathcal{E}}_{1}\langle {\Psi}_{0},{\Psi}_{0}\rangle $ (87)

Using the hermiticity of ${H}_{0}$ and that ${E}_{0}$ is the eigenvalue of ${H}_{0}$, one gets

${\mathcal{E}}_{1}=\frac{\langle {\Psi}_{0},V{\Psi}_{0}\rangle}{\langle {\Psi}_{0},{\Psi}_{0}\rangle}$ (88)

${\mathcal{E}}_{1}$ represents the shift in *E* to first order; which means that

$E={E}_{0}+\u03f5\frac{\langle {\Psi}_{0},V{\Psi}_{0}\rangle}{\langle {\Psi}_{0},{\Psi}_{0}\rangle}+O\left({\u03f5}^{2}\right)$ (89)

This is a very important result which can be stated as

{The perturbed energy levels can be calculated from those which are not perturbed} [13].

Now, taking the scalar product of Equation (86) with ${\Psi}_{m}$ and using the orthogonality of the eigenfunctions, and the hermiticty of the unperturbed Hamiltonian, one obtains

$\langle {\Psi}_{m},{\eta}_{1}\rangle =\frac{\langle {\Psi}_{m},V{\Psi}_{0}\rangle}{{E}_{0}-{E}_{m}}$ (90)

Hence and in first order, the perturbed eigenfunction is given by

$\Psi \left(x\right)={\Psi}_{0}\left(x\right)+\u03f5{\displaystyle {\sum}_{m\ne 0}\frac{\langle {\Psi}_{m},V{\Psi}_{0}\rangle}{{E}_{0}-{E}_{m}}{\Psi}_{m}\left(x\right)}+O\left({\u03f5}^{2}\right)$ (91)

This result is for the ground state; a similar result can be obtained for the other energy levels ${E}_{n}$ and

$\Psi \left(x\right)={\Psi}_{n}\left(x\right)+\u03f5{\displaystyle {\sum}_{m\ne n}\frac{\langle {\Psi}_{m},V{\Psi}_{n}\rangle}{{E}_{n}-{E}_{m}}{\Psi}_{m}\left(x\right)}+O\left({\u03f5}^{2}\right)$ (92)

Note that $\Psi \left(x\right)$ is normalized to first order in $\u03f5$, namely $\langle \Psi |\Psi \rangle =1+O\left({\u03f5}^{2}\right)$ ; where Dirac notation is used here.

Example 7. Infinite Square Well

Consider a particle moving under the potential described as

$V\left(x\right)=\{\begin{array}{l}{V}_{0}\mathrm{cos}\left(\frac{\pi x}{2a}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left|x\right|<a\\ \infty ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left|x\right|\ge a\end{array}$, and assume that it is required to compute the

ground state energy to first order using perturbation theory. Then, to do that we need to seek solutions of the corresponding infinite square well problem. The

potential in this case is given by ${V}_{p}\left(x\right)=\{\begin{array}{l}0,\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\left|x\right|<a\\ \infty ,\text{\hspace{0.17em}}\left|x\right|\ge a\end{array}$ ; this means that the particle

(*m*) is moving inside a box (
$\left|x\right|<a$ ) freely, while there is a very strong repulsion force at the walls (
$x=\pm a$ ). Hence, the time-independent Schrodinger equation

is given by $Eu=-\frac{{\hslash}^{2}}{2m}\frac{{\text{d}}^{2}u}{\text{d}{x}^{2}}$, its solution is then ${u}_{n}\left(x\right)=\frac{1}{\sqrt{a}}\{\begin{array}{l}\mathrm{cos}\left(\frac{n\pi x}{2a}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{odd}\text{\hspace{0.17em}}n\\ \mathrm{sin}\left(\frac{n\pi x}{2a}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{even}\text{\hspace{0.17em}}n\end{array}$ ; moreover, the energy eigenvalues are given by ${E}_{n}=\frac{{n}^{2}{\pi}^{2}{\hslash}^{2}}{8m{a}^{2}}$ ; and the least value energy ( ground state) is ${E}_{1}=\frac{{\pi}^{2}{\hslash}^{2}}{8m{a}^{2}}$ [13].

Now, going back to the original problem, we see that the ground state energy needs to be calculated to first order where the perturbation is given by

${H}^{\prime}={V}_{0}\mathrm{cos}\left(\frac{\pi x}{2a}\right)$ ; this is small compared with the term ${E}_{2}-{E}_{1}$.

Or ${V}_{0}\ll {E}_{2}-{E}_{1}$. Using perturbation theory techniques, one obtains the shift in ground state energy as

$\Delta E={\mathcal{E}}_{1}={{H}^{\prime}}_{11}={\displaystyle {\int}_{-\infty}^{\infty}{u}_{1}^{*}{H}^{\prime}{u}_{1}\text{d}x}=\frac{{V}_{0}}{a}{\displaystyle {\int}_{-a}^{a}{\mathrm{cos}}^{3}\left(\frac{\pi x}{2a}\right)\text{d}x}=\frac{8{V}_{0}}{3\pi}=0.85{V}_{0}$

(it is to be noted that *n* = 1 for the ground state energy) [13] [14].

Example 8. Simple Harmonic Oscillator under a Linear Perturbation

The unperturbed Hamiltonian for the simple harmonic oscillator in one dimension is given by

${H}_{0}=-\frac{{\hslash}^{2}}{2m}\frac{{\text{d}}^{2}}{\text{d}{x}^{2}}+\left(\frac{m{\omega}_{0}^{2}}{2}\right){x}^{2}=\frac{1}{2}\left(-\frac{{\text{d}}^{2}}{\text{d}{\rho}^{2}}+{\rho}^{2}\right)\hslash {\omega}_{0}$,

where $\rho =\sqrt{\frac{m{\omega}_{0}}{\hslash}}x$. Putting $\hslash {\omega}_{0}=1$, one gets for the unperturbed states energy ${E}_{n}^{\left(0\right)}=n+\frac{1}{2}$, $n=0,1,2,\cdots $ ; the superscript 0 means that the energy is for the unperturbed case. Also, the unperturbed normalized eigenfunctions are given by ${\Psi}_{n}^{\left(0\right)}\left(\rho \right)={N}_{n}{H}_{n}\left(\rho \right){\text{e}}^{-\frac{{\rho}^{2}}{2}}$ ; ${N}_{n}={\pi}^{-1/4}{\left({2}^{n}n!\right)}^{2}$ ; ${H}_{n}\left(\rho \right)$ are Hermite polynomials.

Now, in the presence of the linear perturbation $\beta \rho $, we need to compute the matrix elements of $\rho $ ; which is $\langle n|\rho |m\rangle ={\displaystyle {\int}_{-\infty}^{\infty}{\Psi}_{n}^{\left(0\right)}\left(\rho \right)\rho \text{\hspace{0.05em}}{\Psi}_{m}^{\left(0\right)}\left(\rho \right)\text{d}\rho}$ ; again Dirac notation is used here for the eigenfunctions. Using the recurrence relation for Hermite polynomials [10], ${H}_{n+1}\left(\rho \right)=2\rho {H}_{n}\left(\rho \right)-2n{H}_{n-1}\left(\rho \right)$, and the form

of ${\Psi}_{n}^{\left(0\right)}\left(\rho \right)$, one obtains that $\rho {\Psi}_{n}^{\left(0\right)}\left(\rho \right)=\sqrt{\frac{n+1}{2}}{\Psi}_{n+1}^{\left(0\right)}\left(\rho \right)+\sqrt{\frac{n}{2}}{\Psi}_{n-1}^{\left(0\right)}\left(\rho \right)$. Taking the product of this equation with ${\Psi}_{m}^{\left(0\right)}\left(\rho \right)$ and using the orthogonality of these functions, we get $\langle m|\rho |n\rangle =\langle n|\rho |m\rangle ={\delta}_{m,n+1}\sqrt{\frac{n+1}{2}}+{\delta}_{m,n-1}\sqrt{\frac{n}{2}}$ ; and hence, for a linear perturbation of the form ${H}^{\prime}=\beta \rho $, the total Hamiltonian is $H={H}_{0}+{H}^{\prime}=\frac{1}{2}\left(-\frac{{\text{d}}^{2}}{\text{d}{\rho}^{2}}+{\rho}^{2}\right)+\beta \rho =\frac{1}{2}\left(-\frac{{\text{d}}^{2}}{\text{d}{\mu}^{2}}+{\mu}^{2}\right)-\frac{1}{2}{\beta}^{2}$ ; where $\mu =\rho +\beta $. Therefore, the total energy, including the perturbation term, is given by ${E}_{n}={E}_{n}^{\left(0\right)}-\frac{1}{2}{\beta}^{2}=n+\frac{1}{2}-\frac{1}{2}{\beta}^{2}$ [14].

3.4. Degenerate States

The correction in the eigenfunctions, in first order of perturbation, is small only when

$\frac{\langle {\Psi}_{m}|V|{\Psi}_{n}\rangle}{{E}_{n}-{E}_{m}}\cong 1$ (93)

If the difference between the unperturbed energy eigenvalues is small compared with the matrix element in the numerator ${V}_{mn}$, then the perturbation will be large and the approximation fails. In particular, if there exist levels of the degenerate type (levels with the same energy), the denominator will be equal to zero which is a singular situation and the solution obtained for the nondegenerate states is not acceptable.

Now, if unperturbed state is *l*-fold degenerate, then we use the projection
$P$ on the state and the projection
$Q$ perpendicular to it. Acting with
$P$ on the space of these degenerate states, one gets

$P\left(V-{\mathcal{E}}_{1}\right){\Psi}_{0}=0$ (94)

Choosing an appropriate basis for the space of the degenerate states ${\left\{{\varnothing}_{i}\right\}}_{i=1}^{l}$ ; ${\Psi}_{0}$ can be written as

${\Psi}_{0}={\displaystyle {\sum}_{i=1}^{l}{c}_{i}{\varnothing}_{i}}$ (95)

So, one gets

$\langle {\varnothing}_{i}|V|{\varnothing}_{j}\rangle {c}_{j}={\mathcal{E}}_{1}{c}_{i}$ (96)

Hence,
${\mathcal{E}}_{1}$ is an eigenvalue of the matrix
$\left({V}_{ij}\right)$. this will lead to *l* roots. Equation (96) is a generalization of Equation (89) for the nondegenerate case. Moreover, if the states are degenerate, one can, partially, remove the degeneracy as will be shown in the following example.

Example 9. Stark Effect

The Schrodinger equation for the Hydrogen atom is given by [10]

$-\frac{{\hslash}^{2}}{2m}{\nabla}^{2}u-\frac{{e}^{2}}{r}u=Eu$ (97)

$m,e$ are the mass and charge of the electron which is orbiting around the proton.

Equation (96) can be written as
$Hu=Eu$, where
$H$ is the Hamiltonian and *u* is the eigenfunction which describes the electron motion under the attractive

force
$-\frac{{e}^{2}}{r}$ ; while *E* is the eigenvalue (the energy). Using separation of variables in spherical coordinates, the eigenvalues are found to be [10]

${E}_{n}=-\frac{m{e}^{4}}{2{n}^{2}{\hslash}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=1,2,3,\cdots $ (98)

*n* is the principal quantum number and for each *n* there exist an orbital angular quantum number *l* such that
$l<n$ and
$l=0,1,2,\cdots ,n-1$ ; also, there exist a magnetic quantum number *m* such that
$-l\le m\le l$. This means that for each orbital quantum number *l* the degeneracy in energy states equals
$2l+1$.

The state with
$l=0$ is labelled as *s*-state, while the one with
$l=1$ is called the *p*-state, … and so on. Moreover, the ground eigenstate for Hydrogen atom is

given ${u}_{100}\left(\stackrel{\to}{r}\right)=\frac{1}{\sqrt{\pi}}{\left(\frac{1}{{a}_{0}}\right)}^{\frac{3}{2}}\mathrm{exp}\left(-\frac{r}{a}\right)$ [10]. Note that, in this case,

$n=1\to l=0=m$. In terms of the fine structure formula, the wave function in this situation is expressed as $|1s\rangle $. If $n=2$, then $l=0,1$ and there is a 4-fold degeneracy with the same energy and the eigenstates are $|2s\rangle $, $|2p,0\rangle $, $|2p,+1\rangle $, $|2p,-1\rangle $,

Now, if the Hydrogen atom is acted upon by a weak electric field $\stackrel{\to}{E}$ in the z-direction, and $\stackrel{\to}{E}=\mathbb{E}\stackrel{^}{z}$ ; $\mathbb{E}$ is the magnitude of the electric field. The potential, in this case, is $V=-ez\mathbb{E}$.

To compute the matrix elements of
${V}_{ij}$, for the unperturbed states for the Hydrogen atom with energy
${E}_{2}$, we note that we are computing a 4 × 4 Hermitian matrix keeping in mind that the perturbing potential is odd under parity operation; and hence the only non-zero elements are those with different parity (e.g.
${V}_{01}\ne 0$, and
${V}_{10}\ne 0$, this because
$|2s\rangle $ is of even parity, while
$|2p,0\rangle $ is of odd parity). Note, also that the Hydrogen eigenstates are eigenstates of the total angular momentum
${L}^{2}$ and the *z*-component of the angular momentum
${L}_{z}$. Therefore,
${V}_{ij}$ between states *p* and *s* may not be zeros; moreover *V* commutes with
${L}_{z}$. Taking previous information into account, the only non-zero elements are
$\langle 2s|V|2p,0\rangle $ and its Hermitian conjugate. *i.e.*

$\left({V}_{ij}\right)=\left(\begin{array}{cccc}0& 3e\mathbb{E}{a}_{0}& 0& 0\\ 3e\mathbb{E}{a}_{0}& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right)$. Note that to reach these results, the forms of

the Hydrogen eigenstates were used; ${a}_{0}$ is Bohr radius. From these results we conclude that the electric field was able to remove the degeneracy between the two states $|2p,0\rangle $ and $|2s\rangle $ (it made what was hidden apparent); while the other two states $|2p,\pm 1\rangle $ are still without any change in their status. What happened for the states $|2s\rangle $ and $|2p,0\rangle $ was a discretization in the energy, so that the energy of one of them becomes ${E}_{0}+3e\mathbb{E}{a}_{0}$ and the other is ${E}_{0}-3e\mathbb{E}{a}_{0}$. This is the so-called Stark effect [13].

Example 10. The Approximation in Calculating the Ground State Energy for the Helium Atom

In the He atom, the nucleus is composed of two protons and two neutrons with two electrons moving around the nucleus. The unperturbed total Hamiltonian is given by ${H}_{0}={H}_{1}+{H}_{2}$ where ${H}_{i},i=1,2$ are the Hamiltonians for the first and second electrons; the total, ground state, eigenfunction of the total unperturbed Hamiltonian is given by ${\Psi}_{0}={u}_{100}\left({\stackrel{\to}{r}}_{1}\right){u}_{100}\left({\stackrel{\to}{r}}_{2}\right){\psi}_{0,0}$ ; ${\stackrel{\to}{r}}_{1},{\stackrel{\to}{r}}_{2}$ are the position vectors for electron 1 and electron 2 respectively. While, ${\psi}_{0,0}$ has to do with the spin of the ground state.

For the total Hamiltonian $H$, when the repulsion between the two electrons is taken into account ( ${H}^{\prime}$ ) and which is considered as the perturbation in the

problem, is given by $H={H}_{0}+{H}^{\prime}$ and ${H}^{\prime}=\frac{{e}^{2}}{4\pi {\u03f5}_{0}{r}_{12}}$, ${r}_{12}=\left|{\stackrel{\to}{r}}_{1}-{\stackrel{\to}{r}}_{2}\right|$.

Now, to compute the first order correction in the energy, we compute the expectation value for ${H}^{\prime}$, and where ${u}_{100}\left(\stackrel{\to}{r}\right)=\frac{1}{\sqrt{\pi}}{\left(\frac{Z}{{a}_{0}}\right)}^{\frac{3}{2}}\mathrm{exp}\left(-\frac{Zr}{{a}_{0}}\right)$ ; hence $\Delta {E}_{1}=\frac{{a}_{0}}{4\pi {\u03f5}_{0}}{\left(\frac{{Z}^{3}}{\pi {a}_{0}^{3}}\right)}^{2}{{\displaystyle \int}}^{\text{}}\frac{1}{{r}_{12}}\mathrm{exp}\left(-\frac{2Z\left({r}_{1}+{r}_{2}\right)}{{a}_{0}}\right)\text{d}{\tau}_{1}\text{d}{\tau}_{2}=\frac{5}{4}Z{R}_{y}=34\text{\hspace{0.17em}}\text{eV}$ ; ( $Z=2$ is the atomic number for Helium and ${R}_{y}\left[=13.606\text{\hspace{0.17em}}\text{eV}\right]$ is Rydberg constant) [13].

Therefore, ${E}_{1}=-108.8\text{\hspace{0.17em}}\text{eV}+34\text{\hspace{0.17em}}\text{eV}=-74.8\text{\hspace{0.17em}}\text{eV}$, and this represents the corrected value of the energy, in first order, for the Helium ground state [13] [15].

Example 11. Zeeman Effect

In this problem, an electron (−*e*) moves under the effect of a field resulting from a nucleus with positive charge *Ze* positioned at the center of the coordinates system
$x,y,z$ ; moreover, there is a uniform magnetic field in the z-direction

so that the wave equation is given by $-{\nabla}^{2}u-\frac{2Z}{r}u+\frac{\u03f5}{i}\left(x{u}_{y}-y{u}_{x}\right)=\lambda u$, $\u03f5$ is a

small parameter representing the perturbation and $\lambda $ is related to the energy (it is the eigenvalue of the equation).

Now, if $\u03f5=0$ (no perturbation), the eigenvalues, ${\lambda}_{n}$, are then given by

${\lambda}_{n}=-\frac{{Z}^{2}}{{n}^{2}}$, $n=1,2,3,\cdots $ ; and the corresponding eigentates are ${u}_{n;i,m}\left(r,\theta ,\phi \right)={f}_{n;l}\left(r\right){P}_{l}^{m}\left(\mathrm{cos}\theta \right){\text{e}}^{im\phi}$, $l=0,1,\cdots ,n-1$ ; $m=0,\pm 1,\cdots ,\pm l$. Bu, $\frac{1}{i}\left(x{u}_{y}-y{u}_{x}\right)=\frac{1}{i}\frac{\partial u}{\partial \phi}$. The effect of this last operator is clear, since $\frac{1}{i}\frac{\partial {u}_{n;i,m}}{\partial \phi}=m{u}_{n;l,m}$. Hence, the eigenvalue, in the presence of perturbation, is given by ${\lambda}_{n}\left(\u03f5\right)=-\frac{{Z}^{2}}{{n}^{2}}+\u03f5m$. While the eigenfunctions do not change with the perturbation [6].

Thence, the magnetic field removed the degeneracy in eigenvalues, e.g. if $l=1$

then $m=0,\pm 1$ ; and the energy levels become $-\frac{{Z}^{2}}{{n}^{2}}-\u03f5,-\frac{{Z}^{2}}{{n}^{2}},-\frac{{Z}^{2}}{{n}^{2}}+\u03f5$, respectively.

4. Concluding Discussion

In this brief review article, we have shown the importance of perturbation theory in evaluating the eigenvalues and eigenfunctions for physical systems in many applied fields, in classical mechanics, fluid mechanics, heat flow in certain bodies, and most prominent in quantum mechanics where we gave sample examples of well-known problems in the literature and cited in the references given below.

Perturbation theory will keep to be a vital subject leading to many applications in many fields; as an example, recently a PhD dissertation is written on perturbation theory and its role in the evaluation of Feyman Integrals using integration by parts identities [16].

References

[1] Awin, A.M. (2010) Numerical Methods. 1st Edition, Misurata University Publishing, Misurata.

[2] Awin, A.M. (2003) Lectures on Mathematical Methods. 1st Edition, Dar Alkitaab Aljadeed, Beirut.

[3] Benmousa, K.A. (2014) The Numerical Solution of Schrodinger Equation. MSc Thesis, University of Tripoli, Tripoli.

[4] Robertson, E. (2019) Perturbation Theory.

http://en.wikpedia.org/wiki/perturbationtheory

[5] Jenann, I. (2019) Theory-Quantum Mechanics.

http://en.wikpedia.org/wiki/perturbation

[6] Rellich, F. (1969) Perturbation Theory of Eigenvalue Problems. 1st Edition, Gordon and Breach Science Publishers, New York.

[7] Hunter, J. (2017) Asymptotic Analysis and Singular Perturbation Theory. Chapter 1, 1-17.

http://www.math.ucdavis.edu

[8] Dushman, S. (2018) Element of the Quantum Theory.

https://www.scribd.com/document/378218401/Perturbation-theory

[9] Young, W.R. (2017) Part A: Perturbation Theory. University of California at San Diego, La Jolla, 1-9.

[10] Ziock, K. (1969) Basic Quantum Mechanics. 1st Edition. Willey and Sons, New York.

[11] Rienstra, S. (2018) JMBC Course.

https://www.win.tue.nl/~sjoerdr/Q/Onderwijs/JMBC-Cursus-2018.pdf

[12] Van Dyke, M. (1975) Perturbation Methods in Fluid Mechanics. Annotated Edition, Academic Press, New York.

[13] Debbio, L. (2013) Perturbation Theory.

https://www2.ph.ed.ac.uk/~ldeldebb/docs/QM/lect17.pdf

[14] Kerley, G.I. (2008) New Developments in Time-Independent Quantum-Mechanical Perturbation Theory.

https://arxiv.org/ftp/arxiv/papers/1306/1306.6912.pdf

[15] Bransden, B.H. and Joachain, C.J. (2000) Quantum Mechanics. 2nd Edition, Pearcon Education Ltd., Essex.

[16] Georgoudis, A. (2019) Topics in Perturbation Theory. PhD Thesis, Acta Universitatis, Upsaliensis.