spectral theory of ordinary differential equations

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh– Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

Introduction

Spectral theory for second order ordinary differential equations on a compact interval was developed by Jacques Charles François Sturm and Joseph Liouville in the nineteenth century and is now known as Sturm–Liouville theory. In modern language, it is an application of the spectral theorem for compact operators due to David Hilbert. In his dissertation, published in 1910, Hermann Weyl extended this theory to second order ordinary differential equations with singularities at the endpoints of the interval, now allowed to be infinite or semi-infinite. He simultaneously developed a spectral theory adapted to these special operators and introduced boundary conditions in terms of his celebrated dichotomy between ''limit points'' and ''limit circles''.

In the 1920s, John von Neumann established a general spectral theorem for unbounded self-adjoint operators, which Kunihiko Kodaira used to streamline Weyl's method. Kodaira also generalised Weyl's method to singular ordinary differential equations of even order and obtained a simple formula for the spectral measure. The same formula had also been obtained independently by E. C. Titchmarsh in 1946 (scientific communication between Japan and the United Kingdom had been interrupted by World War II). Titchmarsh had followed the method of the German mathematician Emil Hilb, who derived the eigenfunction expansions using complex function theory instead of operator theory. Other methods avoiding the spectral theorem were later developed independently by Levitan, Levinson and Yoshida, who used the fact that the resolvent of the singular differential operator could be approximated by compact resolvents corresponding to Sturm–Liouville problems for proper subintervals. Another method was found by Mark Grigoryevich Krein; his use of ''direction functionals'' was subsequently generalised by Izrail Glazman to arbitrary ordinary differential equations of even order.

Weyl applied his theory to Carl Friedrich Gauss's hypergeometric differential equation, thus obtaining a far-reaching generalisation of the transform formula of Gustav Ferdinand Mehler (1881) for the Legendre differential equation, rediscovered by the Russian physicist Vladimir Fock in 1943, and usually called the Mehler–Fock transform. The corresponding ordinary differential operator is the radial part of the Laplacian operator on 2-dimensional hyperbolic space. More generally, the Plancherel theorem for SL(2,R) of Harish Chandra and Gelfand– Naimark can be deduced from Weyl's theory for the hypergeometric equation, as can the theory of spherical functions for the isometry groups of higher dimensional hyperbolic spaces. Harish Chandra's later development of the Plancherel theorem for general real semisimple Lie groups was strongly influenced by the methods Weyl developed for eigenfunction expansions associated with singular ordinary differential equations. Equally importantly the theory also laid the mathematical foundations for the analysis of the Schrödinger equation and scattering matrix in quantum mechanics.

Solutions of ordinary differential equations

Reduction to standard form

Let be the second order differential operator on given by
$Df(x) = -p(x) f''(x) + r(x) f'(x) + q(x) f(x),$
where is a strictly positive continuously differentiable function and and are continuous real-valued functions.

For in , define the Liouville transformation by
$\psi(x) = \int_^x p(t)^\, dt$

If
$U: L^2(a, b) \mapsto L^2(\psi(a), \psi(b))$
is the unitary operator defined by
$(Uf)(\psi(x)) = f(x) \times \left(\psi'(x)\right)^,\ \ \forall x \in (a, b)$
then
$U \frac U^ g = g' \psi' + \frac 1 2 g \frac$
and
\begin
U \frac U^ g & = \left( U \frac U^ \right) \times \left( U \frac U^ \right) g \\ ex& = \frac \left g' \psi' + \frac 12 g \frac \right\cdot \psi' + \frac 1 2 \left g' \psi' + \frac 12 g \frac \right\cdot \frac \\ ex& = g'' \psi'^2 + 2 g' \psi'' + \frac 12 g \cdot \left \frac - \frac 1 2 \frac \right\end

Hence,
$UDU^ g= -g'' + R g' + Q g,$
where
$R = \frac$
and
$Q = q - \frac + \frac 4 - \frac$

The term in can be removed using an Euler integrating factor. If , then satisfies
$(S U D U^ S^) h = -h'' + V h,$
where the potential ''V'' is given by
$V = Q + \frac$

The differential operator can thus always be reduced to one of the form
$Df = - f'' + qf.$

Existence theorem

The following is a version of the classical Picard existence theorem for second order differential equations with values in a Banach space .

Let , be arbitrary elements of , a bounded operator on and a continuous function on .

Then, for or , the differential equation
$D f = A f$
has a unique solution in satisfying the initial conditions
$f(c) = \beta \, , \; f'(c) = \alpha.$

In fact a solution of the differential equation with these initial conditions is equivalent to a solution
of the integral equation
$f = h + T f$
with the bounded linear map on defined by
$Tf(x) = \int_c^x K(x,y)f(y) \, dy,$
where is the Volterra kernel
$K(x, t) = (x - t) ( q(t) - A)$
and
$h(x) = \alpha (x - c) + \beta.$

Since tends to 0, this integral equation has a unique solution given by the Neumann series
$f = (I - T)^ h = h + T h + T^2 h + T^3 h + \cdots$

This iterative scheme is often called ''Picard iteration'' after the French mathematician Charles Émile Picard.

Fundamental eigenfunctions

If is twice continuously differentiable (i.e. ) on satisfying , then is called an eigenfunction of with eigenvalue .

* In the case of a compact interval and continuous on , the existence theorem implies that for or and every complex number there a unique eigenfunction on with and prescribed. Moreover, for each in , and are holomorphic functions of .
* For an arbitrary interval and continuous on , the existence theorem implies that for in and every complex number there a unique eigenfunction on with and prescribed. Moreover, for each in , and are holomorphic functions of .

Green's formula

If and are functions on , the Wronskian is defined by
$W(f, g) (x) = f(x) g'(x) - f'(x) g(x).$

Green's formula - which in this one-dimensional case is a simple integration by parts - states that for , in
$\int_x^y (Df) g - f (Dg) \, dt = W(f,g)(y) - W(f,g)(x).$

When is continuous and , are on the compact interval , this formula also holds for or .

When and are eigenfunctions for the same eigenvalue, then
$\frac W(f,g) = 0,$
so that is independent of .

Classical Sturm–Liouville theory

Let be a finite closed interval, a real-valued continuous function on and let be the space of functions on satisfying the Robin boundary conditions
\begin
\cos \alpha \,f(a) - \sin \alpha \,f'(a)=0, \\ .5ex \cos \beta \,f(b) - \sin \beta \, f'(b)=0,
\end
with inner product
$(f,g) = \int_a^b f(x) \overline \, dx.$

In practice usually one of the two standard boundary conditions:

*Dirichlet boundary condition
*Neumann boundary condition

is imposed at each endpoint .

The differential operator given by
$Df = -f'' + qf$
acts on . A function in is called an eigenfunction of (for the above choice of boundary values) if for some complex number , the corresponding eigenvalue. By Green's formula, is formally self-adjoint on , since the Wronskian vanishes if both , satisfy the boundary conditions:
$(Df, g) = (f, Dg), \quad \text f, g \in H_0.$

As a consequence, exactly as for a self-adjoint matrix in finite dimensions,
*the eigenvalues of are real;
*the eigenspaces for distinct eigenvalues are orthogonal.

It turns out that the eigenvalues can be described by the maximum-minimum principle of Rayleigh– Ritz (see below). In fact it is easy to see ''a priori'' that the eigenvalues are bounded below because the operator is itself ''bounded below'' on :

In fact, integrating by parts,
(Df,f) = \left f' \overline\righta^b + \int , f', ^2 + \int q , f, ^2.

For Dirichlet or Neumann boundary conditions, the first term vanishes and the inequality holds with .

For general Robin boundary conditions the first term can be estimated using an elementary ''Peter-Paul'' version of Sobolev's inequality:

"Given , there is constant such that for all in ."

In fact, since
$, f(b) - f(x), \leq (b-a)^ \cdot \, f' \, _2,$
only an estimate for is needed and this follows by replacing in the above inequality by for sufficiently large.

Green's function (regular case)

From the theory of ordinary differential equations, there are unique fundamental eigenfunctions , such that

* , ,
* , ,

which at each point, together with their first derivatives, depend holomorphically on . Let
$\omega(\lambda) = W(\phi_\lambda, \chi_\lambda),$
be an entire holomorphic function.

This function plays the role of the characteristic polynomial of . Indeed, the uniqueness of the fundamental eigenfunctions implies that its zeros are precisely the eigenvalues of and that each non-zero eigenspace is one-dimensional. In particular there are at most countably many eigenvalues of and, if there are infinitely many, they must tend to infinity. It turns out that the zeros of also have mutilplicity one (see below).

If is not an eigenvalue of on , define the Green's function by
G_\lambda(x, y) = \begin
\phi_\lambda(x) \chi_\lambda(y) / \omega(\lambda) & \text x \geq y \\ ex\chi_\lambda(x) \phi_\lambda(y) / \omega(\lambda) & \text y \geq x.
\end

This kernel defines an operator on the inner product space via
$(G_\lambda f)(x) =\int_a^b G_\lambda(x,y) f(y)\, dy.$

Since is continuous on , it defines a Hilbert–Schmidt operator on the Hilbert space completion of (or equivalently of the dense subspace ), taking values in . This operator carries into . When is real, is also real, so defines a self-adjoint operator on . Moreover,

* on
* carries into , and on .

Thus the operator can be identified with the resolvent .

Spectral theorem

In fact let for large and negative. Then defines a compact self-adjoint operator on the Hilbert space .
By the spectral theorem for compact self-adjoint operators, has an orthonormal basis consisting of eigenvectors of with , where tends to zero. The range of contains so is dense. Hence 0 is not an eigenvalue of . The resolvent properties of imply that lies in and that
$D \psi_n = \left(\lambda + \frac\right) \psi_n$

The minimax principle follows because if
$\lambda(G) = \min_ \frac,$
then for the linear span of the first eigenfunctions. For any other -dimensional subspace , some in the linear span of the first eigenvectors must be orthogonal to . Hence .

Wronskian as a Fredholm determinant

For simplicity, suppose that on with Dirichlet boundary conditions. The minimax principle shows that
$n^2 + m \le \lambda_n(D) \le n^2 + M.$

It follows that the resolvent is a trace-class operator whenever is not an eigenvalue of and hence that the Fredholm determinant is defined.

The Dirichlet boundary conditions imply that
$\omega(\lambda) = \phi_\lambda(b).$

Using Picard iteration, Titchmarsh showed that , and hence , is an entire function of finite order :
$\omega(\lambda) = \mathcal O\left(e^\right)$

At a zero of , . Moreover,
$\psi(x)=\partial_\lambda \varphi_\lambda(x), _$
satisfies . Thus
$\omega(\lambda) = (\lambda - \mu) \psi(b) + \mathcal O( (\lambda - \mu)^2)$

This implies that

For otherwise , so that would have to lie in . But then
$(\phi_\mu, \phi_\mu) = ((D - \mu)\psi, \phi_\mu) = (\psi, (D - \mu) \phi_\mu) = 0,$
a contradiction.

On the other hand, the distribution of the zeros of the entire function ω(λ) is already known from the minimax principle.

By the Hadamard factorization theorem, it follows that
$\omega(\lambda) = C \prod (1 -\lambda/\lambda_n),$
for some non-zero constant .

Hence
$\det ( I - \mu(D - \lambda)^) = \prod \left( 1 - \right) = \prod = .$

In particular if 0 is not an eigenvalue of
$\omega(\mu) = \omega(0) \cdot \det ( I - \mu D^) .$

Tools from abstract spectral theory

Functions of bounded variation

A function of bounded variation on a closed interval is a complex-valued function such that its total variation , the supremum of the variations
$\sum_^ , \rho(x_) - \rho(x_r),$
over all dissections
$a= x_0 < x_1 < \dots < x_k =b$
is finite. The real and imaginary parts of are real-valued functions of bounded variation. If is real-valued and normalised so that , it has a canonical decomposition as the difference of two bounded non-decreasing functions:
$\rho(x) = \rho_+(x) - \rho_-(x),$
where and are the total positive and negative variation of over .

If is a continuous function on its Riemann–Stieltjes integral with respect to
$\int_a^b f(x)\, d\rho(x)$
is defined to be the limit of approximating sums
$\sum_^ f(x_r)(\rho(x_)-\rho(x_r))$
as the mesh of the dissection, given by , tends to zero.

This integral satisfies
$\left, \int_a^b f(x)\, d\rho(x)\\le V(\rho)\cdot \, f\, _\infty$

and thus defines a bounded linear functional on with norm .

Every bounded linear functional on has an absolute value , μ, defined for non-negative by
$, \mu, (f) = \sup_ , \mu(g), .$

The form extends linearly to a bounded linear form on with norm and satisfies the characterizing inequality
$, \mu(f), \leq , \mu, (, f, )$
for in . If is ''real'', i.e. is real-valued on real-valued functions, then
$\mu = , \mu, - (, \mu, -\mu)\equiv \mu_+-\mu_-$
gives a canonical decomposition as a difference of ''positive'' forms, i.e. forms that are non-negative on non-negative functions.

Every positive form extends uniquely to the linear span of non-negative bounded lower semicontinuous functions by the formula
$\mu(g) = \lim \mu(f_n),$
where the non-negative continuous functions increase pointwise to .

The same therefore applies to an arbitrary bounded linear form , so that a function of bounded variation may be defined by
$\rho(x) = \mu(\chi_),$
where denotes the characteristic function of a subset of . Thus and .
Moreover and .

This correspondence between functions of bounded variation and bounded linear forms is a special case of the Riesz representation theorem.

The support of is the complement of all points in where is constant on some neighborhood of ; by definition it is a closed subset of . Moreover, , so that if vanishes on .

Spectral measure

Let be a Hilbert space and $T$ a self-adjoint bounded operator on with $0 \leq T \leq I$, so that the spectrum $\sigma(T)$ of $T$ is contained in ,1/math>. If $p(t)$ is a complex polynomial, then by the spectral mapping theorem
$\sigma (p(T)) = p (\sigma(T))$
and hence
$\, p(T)\, \leq \, p\, _\infty$
where $\, \cdot \, _\infty$ denotes the uniform norm on . By the Weierstrass approximation theorem, polynomials are uniformly dense in . It follows that $f(T)$ can be defined \forall f \in C ,1/math>, with
$\sigma (f(T)) = f(\sigma(T))$ and $\, f(T) \, \leq \, f\, _\infty.$

If $0 \leq g \leq 1$ is a lower semicontinuous function on , for example the characteristic function $\chi_$ of a subinterval of , then
$g$ is a pointwise increasing limit of non-negative f_n \in C ,1/math>.

If $\xi$ is a vector in , then the vectors
$\eta_n = f_n(T)\xi$
form a Cauchy sequence in , since, for $n \geq m$,
$\, \eta_n-\eta_m\, ^2 \le (\eta_n,\xi) - (\eta_m,\xi),$
and $( \eta_n, \xi) = (f_n(T)\xi, \xi)$ is bounded and increasing, so has a limit.

It follows that $g(T)$ can be defined by
$g(T)\xi = \lim f_n(T)\xi.$

If $\xi$ and are vectors in , then
$\mu_(f) = (f(T) \xi,\eta)$
defines a bounded linear form $\mu_$ on . By the Riesz representation theorem
$\mu_ = d\rho_$
for a unique normalised function $\rho_$ of bounded variation on .

$d\rho_$ (or sometimes slightly incorrectly $\rho_$ itself) is called the spectral measure determined by $\xi$ and .

The operator $g(T)$ is accordingly uniquely characterised by the equation
$(g(T)\xi,\eta) = \mu_(g) = \int_0^1 g(\lambda) \, d\rho_(\lambda).$

The spectral projection $E(\lambda)$ is defined by
$E(\lambda) = \chi_(T),$
so that
$\rho_(\lambda) = (E(\lambda)\xi,\eta).$

It follows that
$g(T) = \int_0^1 g(\lambda) \, dE(\lambda),$
which is understood in the sense that for any vectors $\xi$ and $\eta$,
$(g(T)\xi,\eta) = \int_0^1 g(\lambda)\, d(E(\lambda)\xi,\eta) = \int_0^1 g(\lambda)\, d\rho_ (\lambda).$

For a single vector $\xi, \, \mu_ = \mu_$ is a positive form on
(in other words proportional to a probability measure on ) and $\rho_ = \rho_$ is non-negative and non-decreasing. Polarisation shows that all the forms $\mu_$ can naturally be expressed in terms of such positive forms, since
$\mu_ = \frac\left(\mu_+i\mu_-\mu_-i\mu_\right)$

If the vector $\xi$ is such that the linear span of the vectors $(T^n\xi)$ is dense in , i.e. $\xi$ is a ''cyclic vector'' for $T$, then the map $U$ defined by
U(f) = f(T)\xi, \, C ,1\rightarrow H
satisfies
$(Uf_1,Uf_2)= \int_0^1 f_1(\lambda) \overline \, d\rho_\xi(\lambda).$

Let L_2( ,1 d\rho_\xi) denote the Hilbert space completion of C ,1/math> associated with the possibly degenerate inner product on the right hand side.
Thus $U$ extends to a unitary transformation of L_2( ,1 \rho_\xi) onto . $U T U^\ast$ is then just multiplication by $\lambda$ on L_2( ,1 d\rho_\xi); and more generally $U f(T) U^\ast$ is multiplication by $f(\lambda)$. In this case, the support of $d\rho_\xi$ is exactly $\sigma(T)$, so that

Weyl–Titchmarsh–Kodaira theory

The eigenfunction expansion associated with singular differential operators of the form
$Df = -(pf')' + qf$
on an open interval requires an initial analysis of the behaviour of the fundamental eigenfunctions near the endpoints and to determine possible boundary conditions there. Unlike the regular Sturm–Liouville case, in some circumstances spectral values of can have multiplicity 2. In the development outlined below standard assumptions will be imposed on and that guarantee that the spectrum of has multiplicity one everywhere and is bounded below. This includes almost all important applications; modifications required for the more general case will be discussed later.

Having chosen the boundary conditions, as in the classical theory the resolvent of , for large and positive, is given by an operator corresponding to a Green's function constructed from two fundamental eigenfunctions. In the classical case was a compact self-adjoint operator; in this case is just a self-adjoint bounded operator with . The abstract theory of spectral measure can therefore be applied to to give the eigenfunction expansion for .

The central idea in the proof of Weyl and Kodaira can be explained informally as follows. Assume that the spectrum of lies in and that and let
$E(\lambda) = \chi_(T)$
be the spectral projection of corresponding to the interval . For an arbitrary function define
$f(x,\lambda) = (E(\lambda)f)(x).$
may be regarded as a differentiable map into the space of functions of bounded variation ρ; or equivalently as a differentiable map
$x\mapsto (d_\lambda f)(x)$
into the Banach space of bounded linear functionals on whenever is a compact subinterval of .

Weyl's fundamental observation was that satisfies a second order ordinary differential equation taking values in :
$D (d_\lambda f) = \lambda \cdot d_\lambda f.$

After imposing initial conditions on the first two derivatives at a fixed point , this equation can be solved explicitly in terms of the two fundamental eigenfunctions and the "initial value" functionals
$(d_\lambda f)(c) = d_\lambda f(c,\cdot), \quad (d_\lambda f)^\prime(c) = d_\lambda f_x(c,\cdot).$

This point of view may now be turned on its head: and may be written as
$f(c,\lambda) = (f,\xi_1(\lambda)), \quad f_x(c,\lambda) = (f,\xi_2(\lambda)),$
where and are given purely in terms of the fundamental eigenfunctions.
The functions of bounded variation
$\sigma_(\lambda) = (\xi_i(\lambda),\xi_j(\lambda))$
determine a spectral measure on the spectrum of and can be computed explicitly from the behaviour of the fundamental eigenfunctions (the Titchmarsh–Kodaira formula).

Limit circle and limit point for singular equations

Let be a continuous real-valued function on and let be the second order differential operator
$Df = -f'' + qf$
on . Fix a point in and, for complex , let $\varphi_\lambda, \theta_\lambda$ be the unique fundamental eigenfunctions of on satisfying
$(D-\lambda)\varphi_\lambda = 0, \quad (D-\lambda)\theta_\lambda =0$
together with the initial conditions at
$\varphi_\lambda(c) = 1,\, \varphi_\lambda'(c) = 0, \, \theta_\lambda(c) = 0, \, \theta_\lambda'(c) = 1.$

Then their Wronskian satisfies
$W(\varphi_\lambda,\theta_\lambda) = \varphi_\lambda\theta_\lambda'- \theta_\lambda \varphi_\lambda' \equiv 1,$

since it is constant and equal to 1 at .

Let be non-real and . If the complex number $\mu$ is such that $f = \varphi + \mu \theta$ satisfies the boundary condition $\cos\beta\, f(x) - \sin\beta\, f'(x) = 0$ for some $\beta$ (or, equivalently, $f'(x) / f(x)$ is real) then, using integration by parts, one obtains
$\operatorname(\lambda) \int_c^x , \varphi +\mu \theta, ^2 = \operatorname(\mu).$

Therefore, the set of satisfying this equation is not empty. This set is a circle in the complex -plane. Points in its interior are characterized by
$\int_c^x , \varphi + \mu \theta, ^2 <$
if and by
$\int_x^c , \varphi + \mu \theta, ^2 <$
if .

Let be the closed disc enclosed by the circle. By definition these closed discs are nested and decrease as approaches or . So in the limit, the circles tend either to a limit circle or a limit point at each end. If $\mu$ is a limit point or a point on the limit circle at or , then $f = \varphi + \mu\theta$ is square integrable () near or , since $\mu$ lies in for all (in the ∞ case) and so $\int_c^x , \varphi +\mu \theta, ^2 <$ is bounded independent of . In particular:.

* there are always non-zero solutions of which are square integrable near resp. ;
* in the limit circle case all solutions of are square integrable near resp. .

The radius of the disc can be calculated to be
$\left, \$
and this implies that in the limit point case $\theta$ cannot be square integrable near resp. . Therefore, we have a converse to the second statement above:

* in the limit point case there is exactly one non-zero solution (up to scalar multiples) of Df = λf which is square integrable near resp. .

On the other hand, if for another value , then
$h(x) = g(x) -(\lambda^\prime-\lambda) \int_c^x (\varphi_\lambda(x) \theta_\lambda(y) - \theta_\lambda(x)\varphi_\lambda(y))g(y)\, dy$
satisfies , so that
$g(x )= c_1 \varphi_\lambda + c_2 \theta_\lambda + (\lambda^\prime-\lambda) \int_c^x (\varphi_\lambda(x) \theta_\lambda(y) - \theta_\lambda(x)\varphi_\lambda(y))g(y)\, dy.$

This formula may also be obtained directly by the variation of constant method from .
Using this to estimate , it follows that

* the limit point/limit circle behaviour at or is independent of the choice of .

More generally if for some function , then
$g(x)=c_1 \varphi_\lambda + c_2 \theta_\lambda - \int_c^x (\varphi_\lambda(x) \theta_\lambda(y) - \theta_\lambda(x) \varphi_\lambda(y)) r(y) g(y)\, dy.$

From this it follows that

* if is continuous at , then is limit point or limit circle at precisely when is,

so that in particular

* if is continuous at , then is limit point at if and only if .

Similarly

* if has a finite limit at , then is limit point or limit circle at precisely when is,

so that in particular

* if has a finite limit at , then is limit point at .

Many more elaborate criteria to be limit point or limit circle can be found in the mathematical literature.

Green's function (singular case)

Consider the differential operator
$D_0 f = -(p_0 f')' + q_0 f$
on with positive and continuous on and continuously differentiable in , positive in and .

Moreover, assume that after reduction to standard form becomes the equivalent operator
$Df= -f'' + qf$
on where has a finite limit at . Thus
* is limit point at .

At 0, may be either limit circle or limit point. In either case there is an eigenfunction with and square integrable near . In the limit circle case, determines a boundary condition at :
$W(f,\Phi_0)(0) = 0.$

For complex , let and satisfy

* ,
* square integrable near infinity
* square integrable at if is ''limit point''
* satisfies the boundary condition above if is ''limit circle''.

Let
$\omega(\lambda) = W(\Phi_\lambda,\Chi_\lambda),$
a constant which vanishes precisely when and are proportional, i.e. is an eigenvalue of for these boundary conditions.

On the other hand, this cannot occur if or if is negative.

Indeed, if with , then by Green's formula , since is constant. So must be real. If is taken to be real-valued in the realization, then for
_0 f f'x^y = \int_x^y (q_0 -\lambda), f, ^2 + p_0 (f')^2 .

Since and is integrable near , must vanish at . Setting , it follows that , so that is increasing, contradicting the square integrability of near .

Thus, adding a positive scalar to , it may be assumed that
$\omega(\lambda) \neq 0 ~~ \text \lambda \notin [1, \infty).$

If , the Green's function at is defined by
G_\lambda(x,y) = \begin
\Phi_\lambda(x) \Chi_\lambda(y)/\omega(\lambda) & (x\le y), \\ ex\Chi_\lambda(x) \Phi_\lambda(y)/\omega(\lambda) & (x\ge y).
\end
and is independent of the choice of and .

In the examples there will be a third "bad" eigenfunction defined and holomorphic for not in such that satisfies the boundary conditions at neither nor . This means that for not in

* is nowhere vanishing;
* is nowhere vanishing.

In this case is proportional to , where
$m(\lambda) = - W(\Phi_\lambda, \Chi_\lambda) / W(\Psi_\lambda, \Chi_\lambda).$

Let be the space of square integrable continuous functions on and let be

* the space of functions on of compact support if is limit point at
* the space of functions on with at and with near if is limit circle at .

Define by
$(Tf)(x) =\int_0^\infty G_0(x,y)f(y) \, dy.$

Then on , on and the operator is bounded below on :
$(Df,f) \ge (f,f).$

Thus is a self-adjoint bounded operator with .

Formally . The corresponding operators defined for not in can be formally identified with
$(D-\lambda)^ = T(I-\lambda T)^$
and satisfy on , on .

Spectral theorem and Titchmarsh–Kodaira formula

Kodaira gave a streamlined version of Weyl's original proof. ( M.H. Stone had previously shown how part of Weyl's work could be simplified using von Neumann's spectral theorem.)

In fact for with , the spectral projection of is defined by
$E(\lambda) =\chi_(T)$

It is also the spectral projection of corresponding to the interval .

For in define
$f(x,\lambda) = (E(\lambda)f)(x).$

may be regarded as a differentiable map into the space of functions of bounded variation; or equivalently as a differentiable map
$x \mapsto (d_\lambda f)(x)$
into the Banach space of bounded linear functionals on for any compact subinterval of .

The functionals (or measures) satisfies the following -valued second order ordinary differential equation:
$D (d_\lambda f) = \lambda \cdot d_\lambda f,$
with initial conditions at in
$(d_\lambda f)(c)= d_\lambda f(c,\cdot)=\mu^, \quad (d_\lambda f)^\prime(c)= d_\lambda f_x(c,\cdot)=\mu^.$

If and are the special eigenfunctions adapted to , then
$d_\lambda f (x) = \varphi_\lambda(x) \mu^ + \chi_\lambda(x) \mu^.$

Moreover,
$\mu^= d_\lambda (f,\xi^_\lambda),$
where
$\xi^_\lambda = D E(\lambda) \eta^,$
with
$\eta_z^(y) = G_z(c,y), \,\,\,\, \eta_z^(x)=\partial_x G_z(c,y), \,\,\,\, (z \notin [1,\infty)).$
(As the notation suggests, and do not depend on the choice of .)

Setting
$\sigma_(\lambda) = (\xi^_\lambda, \xi^_\lambda),$
it follows that
$d_\lambda (E(\lambda)\eta_z^,\eta_z^) = , \lambda - z, ^ \cdot d_\lambda \sigma_(\lambda).$

On the other hand, there are holomorphic functions , such that

* is proportional to ;
* is proportional to .

Since , the Green's function is given by
G_\lambda(x,y) = \begin
\dfrac & (x\le y), \\ ex\dfrac & (y\le x).
\end

Direct calculation shows that
$(\eta_z^,\eta_z^) = \operatorname M_(z)/ \operatorname z,$
where the so-called is given by
$M_(z)= \frac,\,\, M_(z) = M_(z) = \frac, \,\, M_(z) = \frac.$

Hence
$\int_^\infty (\operatorname z)\cdot, \lambda-z, ^\, d\sigma_(\lambda) = \operatorname M_(z),$
which immediately implies
$\sigma_(\lambda) = \lim_ \lim_ \int_\delta^ \operatorname M_(t +i\varepsilon)\, dt.$
(This is a special case of the Stieltjes transformation, "Stieltjes inversion formula".)

Setting and , it follows that
$(E(\mu)f)(x) = \sum_\int_0^\mu \int_0^\infty\psi^_\lambda(x)\psi^_\lambda(y) f(y) \, dy \, d\sigma_(\lambda) = \int_0^\mu \int_0^\infty\Phi_\lambda(x) \Phi_\lambda(y) f(y)\, dy \, d\rho(\lambda).$

This identity is equivalent to the spectral theorem and Titchmarsh–Kodaira formula.

Application to the hypergeometric equation

The Mehler–Fock transform concerns the eigenfunction expansion associated with the Legendre differential operator
$Df = -((x^2-1) f')' =-(x^2-1)f'' -2x f'$
on . The eigenfunctions are the Legendre functions
$P_(\cosh r) = \int_0^ \left( \right)^\, d\theta$
with eigenvalue . The two Mehler–Fock transformations are
$Uf(\lambda) = \int_1^\infty f(x)\, P_(x) \, dx$
and
$U^g(x)=\int_0^\infty g(\lambda) \, \tanh \pi \sqrt\,d\lambda.$

(Often this is written in terms of the variable .)

Mehler and Fock studied this differential operator because it arose as the radial component of the Laplacian on 2-dimensional hyperbolic space.
More generally, consider the group consisting of complex matrices of the form
$\begin \alpha & \beta \\ \overline & \overline \end$

with determinant .

Application to the hydrogen atom

Generalisations and alternative approaches

A Weyl function can be defined at a singular endpoint giving rise to a singular version of Weyl–Titchmarsh–Kodaira theory. this applies for example to the case of radial Schrödinger operators
$Df = -f'' + \frac f+ V(x) f, \qquad x\in(0,\infty)$

The whole theory can also be extended to the case where the coefficients are allowed to be measures.

Gelfand–Levitan theory

Notes

References

Citations

Bibliography

*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*

{{SpectralTheory

Ordinary differential equations
Operator theory
Spectral theory