Sturm–Liouville theory

Sturm–Liouville theory is a part of mathematics that studies special kinds of equations called Sturm–Liouville problems. These problems involve solving a type of second-order linear ordinary differential equation with certain conditions at the ends of the interval being studied. The main goals are to find special numbers, called eigenvalues, for which there are non-trivial solutions, and then to find the matching solutions, called eigenfunctions.

This theory is very important in applied mathematics because Sturm–Liouville problems appear often, especially when solving separable linear partial differential equations. For example, in quantum mechanics, the one-dimensional time-independent Schrödinger equation is a Sturm–Liouville problem.

The theory is named after Jacques Charles François Sturm and Joseph Liouville, who developed it. For regular Sturm–Liouville problems, there are infinitely many eigenvalues, each with a unique eigenfunction, and these eigenfunctions can be used to build a special space of functions called a Hilbert space.

Main results

Sturm–Liouville theory studies special kinds of equations used in mathematics. These equations help us understand how certain functions behave over a range of values.

The main focus is on solving these equations under specific conditions. The goals are to find special numbers, called eigenvalues, and the functions that go with them, called eigenfunctions. These solutions help us describe many natural patterns and phenomena.

For these problems to be "regular," certain conditions must be met. The functions involved need to be smooth and positive over the interval being studied. There are also specific rules about how the solutions behave at the ends of the interval.

The theory tells us important facts about these solutions, such as whether they cross zero and how often they do so. These properties are useful in many areas of science and engineering.

d d x [ p ( x ) d y d x ] + q ( x ) y = − λ w ( x ) y {\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\left[p(x){\frac {\mathrm {d} y}{\mathrm {d} x}}\right]+q(x)y=-\lambda \,w(x)y}

{ α 1 y ( a ) + α 2 y ′ ( a ) = 0 , α 1 , α 2 not both 0 , β 1 y ( b ) + β 2 y ′ ( b ) = 0 , β 1 , β 2 not both 0. {\displaystyle {\begin{cases}\alpha _{1}y(a)+\alpha _{2}y'(a)&=0,\qquad \alpha _{1},\alpha _{2}{\text{ not both }}0,\\\beta _{1}y(b)+\beta _{2}y'(b)&=0,\qquad \beta _{1},\beta _{2}{\text{ not both }}0.\end{cases}}}

Reduction to Sturm–Liouville form

The equation in the form given is called Sturm–Liouville form or self-adjoint form. All second-order linear homogenous ordinary differential equations can be rewritten into this form by multiplying both sides by a special kind of number called an integrating factor. However, this is not possible for second-order partial differential equations or when y is a vector.

Some examples are shown below.

Bessel equation

This equation can be rewritten into Sturm–Liouville form by dividing and grouping terms.

Legendre equation

This equation can also be rewritten into Sturm–Liouville form by noticing a pattern in its terms.

Example using an integrating factor

An example shows how multiplying by an integrating factor helps rewrite the equation.

Integrating factor for general second-order homogeneous equation

A general method uses an integrating factor to turn any second-order homogenous equation into Sturm–Liouville form.

Sturm–Liouville equations as self-adjoint differential operators

This section talks about a special kind of math problem. It shows how a certain math rule, called L, works with functions. We can think of L as a tool that changes one function into another.

In math, we often look for special numbers called eigenvalues and special functions called eigenfunctions. For this problem, we want to find numbers λ and functions u that follow the rule L u = λ u. This means when we use the tool L on the function u, we get the same function back, but multiplied by the number λ.

This math rule L has a special property: it is self-adjoint. This means that when we use L on two functions and then compare the results in a certain way, it’s the same as comparing the original functions using L in the opposite order. This helps mathematicians understand the problem better.

The study of these problems is important in areas like quantum mechanics, where similar math rules describe how particles behave.

Application to inhomogeneous second-order boundary value problems

When solving certain types of math problems, we sometimes need to find a special solution that fits two different points. This can be tricky, but Sturm–Liouville theory helps by using special functions called eigenfunctions. These functions can be added together in a series to match the needed shape at the two points.

For example, if we want a solution that fits the shape of a straight line between two points, we can use a series of sine waves. Even though these waves might look messy when added together, they still give the right answer when we need it. This method works for many different shapes and conditions, making it a powerful tool in solving complex math problems.

Application to partial differential equations

Certain big math ideas can be solved using Sturm–Liouville theory. For example, think about a thin piece of material stretched tight like a drumhead, held in a square shape. We can study how it wobbles and shakes using special math steps.

One way to understand the wobbling is by breaking the motion into simple patterns that repeat over time, called "normal modes". Each pattern vibrates at its own special pace, and together they make up all possible ways the material can move.

Representation of solutions and numerical calculation

The Sturm–Liouville equation with boundary conditions can be solved in different ways. One way is to use special math methods like the Rayleigh–Ritz method or the matrix-variational method.

There are also numerical methods to solve these problems. Some of these methods include shooting methods, the finite difference method, and the spectral parameter power series method. These methods help find the special numbers, called eigenvalues, that solve the equation correctly.

X ( n ) ( t ) = { − ∫ a x X ( n − 1 ) ( t ) p ( t ) − 1 y 0 ( t ) − 2 d t n odd , ∫ a x X ( n − 1 ) ( t ) y 0 ( t ) 2 w ( t ) d t n even {\displaystyle X^{(n)}(t)={\begin{cases}\displaystyle -\int _{a}^{x}X^{(n-1)}(t)p(t)^{-1}y_{0}(t)^{-2}\,dt&n{\text{ odd}},\\[6pt]\displaystyle \quad \int _{a}^{x}X^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt&n{\text{ even}}\end{cases}}}

X ~ ( n ) ( t ) = { ∫ a x X ~ ( n − 1 ) ( t ) y 0 ( t ) 2 w ( t ) d t n odd , − ∫ a x X ~ ( n − 1 ) ( t ) p ( t ) − 1 y 0 ( t ) − 2 d t n even. {\displaystyle {\tilde {X}}^{(n)}(t)={\begin{cases}\displaystyle \quad \int _{a}^{x}{\tilde {X}}^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt&n{\text{ odd}},\\[6pt]\displaystyle -\int _{a}^{x}{\tilde {X}}^{(n-1)}(t)p(t)^{-1}y_{0}(t)^{-2}\,dt&n{\text{ even.}}\end{cases}}}