# Spectral theory of ordinary differential equations

Last updated

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 TitchmarshKodaira 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 GelfandNaimark 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 D be the second order differential operator on (a,b) given by

${\displaystyle Df(x)=-p(x)f^{\prime \prime }(x)+r(x)f^{\prime }(x)+q(x)f(x),}$

where p is a strictly positive continuously differentiable function and q and r are continuous real-valued functions.

For x0 in (a, b), define the Liouville transformation ψ by

${\displaystyle \psi (x)=\int _{x_{0}}^{x}p(t)^{-1/2}\,dt}$

If

${\displaystyle U:L^{2}(a,b)\mapsto L^{2}(\psi (a),\psi (b))}$

is the unitary operator defined by

${\displaystyle (Uf)(\psi (x))=f(x)\times \left(\psi '(x)\right)^{-1/2},\ \ \forall x\in (a,b)}$

then

${\displaystyle U{\frac {\mathrm {d} }{\mathrm {d} x}}U^{-1}g=g'\psi '+{\frac {1}{2}}g{\frac {\psi ''}{\psi '}}}$

and

{\displaystyle {\begin{aligned}U{\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}U^{-1}g&=\left(U{\frac {\mathrm {d} }{\mathrm {d} x}}U^{-1}\right)\times \left(U{\frac {\mathrm {d} }{\mathrm {d} x}}U^{-1}\right)g\\&={\frac {\mathrm {d} }{\mathrm {d} \psi }}\left[g'\psi '+{\frac {1}{2}}g{\frac {\psi ''}{\psi '}}\right]\cdot \psi '+{\frac {1}{2}}\left[g'\psi '+{\frac {1}{2}}g{\frac {\psi ''}{\psi '}}\right]\times {\frac {\psi ''}{\psi '}}\\&=g''\psi '^{2}+2g'\psi ''+{\frac {1}{2}}g\times \left[{\frac {\psi '''}{\psi '}}+{\frac {\psi ''^{2}}{\psi '^{2}}}\right]\end{aligned}}}

Hence,

${\displaystyle UDU^{-1}g=-g^{\prime \prime }+Rg^{\prime }+Qg,}$

where

${\displaystyle R={\frac {p'+r}{p^{1/2}}}}$

and

${\displaystyle Q=q-{\frac {rp'}{4p}}+{\frac {p''}{4}}-{\frac {p'^{2}}{8p}}}$

The term in g' can be removed using an Euler integrating factor. If S' /S = R/2, then h = Sg satisfies

${\displaystyle (SUDU^{-1}S^{-1})h=-h^{\prime \prime }+Vh,}$

where the potential V is given by

${\displaystyle V=Q+{\frac {S''}{S}}}$

The differential operator can thus always be reduced to one of the form [1]

${\displaystyle Df=-f^{\prime \prime }+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 E. [2]

Let α, β be arbitrary elements of E, A a bounded operator on E and q a continuous function on [a,b].

Then, for c = a or b, the differential equation

Df = Af

has a unique solution f in C2([a,b],E) satisfying the initial conditions

f(c) = β, f '(c) = α.

In fact a solution of the differential equation with these initial conditions is equivalent to a solution of the integral equation

f = h + Tf

with T the bounded linear map on C([a,b], E) defined by

${\displaystyle Tf(x)=\int _{c}^{x}K(x,y)f(y)\,dy,}$

where K is the Volterra kernel

K(x,t)= (xt)(q(t) A)

and

h(x) = α(xc) + β.

Since ||Tk|| tends to 0, this integral equation has a unique solution given by the Neumann series

f = (IT)1h = h + Th + T2h + T3h + ···

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

### Fundamental eigenfunctions

If f is twice continuously differentiable (i.e. C2) on (a, b) satisfying Df = λf, then f is called an eigenfunction of L with eigenvalue λ.

• In the case of a compact interval [a, b] and q continuous on [a, b], the existence theorem implies that for c = a or b and every complex number λ there a unique C2 eigenfunction fλ on [a, b] with fλ(c) and f 'λ(c) prescribed. Moreover, for each x in [a, b], fλ(x) and f 'λ(x) are holomorphic functions of λ.
• For an arbitrary interval (a,b) and q continuous on (a, b), the existence theorem implies that for c in (a, b) and every complex number λ there a unique C2 eigenfunction fλ on (a, b) with fλ(c) and f 'λ(c) prescribed. Moreover, for each x in (a, b), fλ(x) and f 'λ(x) are holomorphic functions of λ.

### Green's formula

If f and g are C2 functions on (a, b), the Wronskian W(f, g) 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 x, y in (a, b)

${\displaystyle \int _{x}^{y}(Df)g-f(Dg)\,dt=W(f,g)(y)-W(f,g)(x).}$

When q is continuous and f, gC2 on the compact interval [a, b], this formula also holds for x = a or y = b.

When f and g are eigenfunctions for the same eigenvalue, then

${\displaystyle {d \over dx}W(f,g)=0,}$

so that W(f, g) is independent of x.

## Classical Sturm–Liouville theory

Let [a, b] be a finite closed interval, q a real-valued continuous function on [a, b] and let H0 be the space of C2 functions f on [a, b] satisfying the Robin boundary conditions

${\displaystyle \cos \alpha \,f(a)-\sin \alpha \,f^{\prime }(a)=0,\qquad \cos \beta \,f(b)-\sin \beta \,f^{\prime }(b)=0,}$

with inner product

${\displaystyle (f,g)=\int _{a}^{b}f(x){\overline {g(x)}}\,dx.}$

In practise usually one of the two standard boundary conditions:

is imposed at each endpoint c = a, b.

The differential operator D given by

${\displaystyle Df=-f^{\prime \prime }+qf}$

acts on H0. A function f in H0 is called an eigenfunction of D (for the above choice of boundary values) if Df = λ f for some complex number λ, the corresponding eigenvalue. By Green's formula, D is formally self-adjoint on H0, since the Wronskian W(f,g) vanishes if both f,g satisfy the boundary conditions:

(Df, g) = (f, Dg) for f, g in H0.

As a consequence, exactly as for a self-adjoint matrix in finite dimensions,

It turns out that the eigenvalues can be described by the maximum-minimum principle of RayleighRitz [3] (see below). In fact it is easy to see a priori that the eigenvalues are bounded below because the operator D is itself bounded below on H0:

• ${\displaystyle (Df,f)\geq M(f,f)}$ for some finite (possibly negative) constant ${\displaystyle M}$.

In fact integrating by parts

${\displaystyle (Df,f)=[-f^{\prime }{\overline {f}}]_{a}^{b}+\int |f^{\prime }|^{2}+\int q|f|^{2}.}$

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

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

"Given ε > 0, there is constant R >0 such that |f(x)|2 ≤ ε (f', f') + R (f, f) for all f in C1[a, b]."

In fact, since

|f(b) f(x)| ≤ (ba)½·||f '||2,

only an estimate for f(b) is needed and this follows by replacing f(x) in the above inequality by (xa)n·(ba)n·f(x) for n sufficiently large.

### Green's function (regular case)

From the theory of ordinary differential equations, there are unique fundamental eigenfunctions φλ(x), χλ(x) such that

• D φλ = λ φλ, φλ(a) = sin α, φλ'(a) = cos α
• D χλ = λ χλ, χλ(b) = sin β, χλ'(b) = cos β

which at each point, together with their first derivatives, depend holomorphically on λ. Let

ω(λ) = W(φλ, χλ),

This function ω(λ) plays the rôle of the characteristic polynomial of D. Indeed, the uniqueness of the fundamental eigenfunctions implies that its zeros are precisely the eigenvalues of D and that each non-zero eigenspace is one-dimensional. In particular there are at most countably many eigenvalues of D 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 D on H0, define the Green's function by

Gλ(x,y) = φλ (x) χλ(y) / ω(λ) for xy and χλ(x) φλ (y) / ω(λ) for yx.

This kernel defines an operator on the inner product space C[a,b] via

${\displaystyle (G_{\lambda }f)(x)=\int _{a}^{b}G_{\lambda }(x,y)f(y)\,dy.}$

Since Gλ(x,y) is continuous on [a, b] x [a, b], it defines a Hilbert–Schmidt operator on the Hilbert space completion H of C[a, b] = H1 (or equivalently of the dense subspace H0), taking values in H1. This operator carries H1 into H0. When λ is real, Gλ(x,y) = Gλ(y,x) is also real, so defines a self-adjoint operator on H. Moreover,

• Gλ (D λ) =I on H0
• Gλ carries H1 into H0, and (D λ) Gλ = I on H1.

Thus the operator Gλ can be identified with the resolvent (D λ)−1.

### Spectral theorem

Theorem.The eigenvalues of D are real of multiplicity one and form an increasing sequence λ1 < λ2 < ··· tending to infinity.

The corresponding normalised eigenfunctions form an orthonormal basis ofH0.

The kth eigenvalue of D is given by the minimax principle

${\displaystyle \lambda _{k}=\max _{{\rm {dim}}\,G=k-1}\,\min _{f\perp G}{(Df,f) \over (f,f)}.}$

In particular if q1 ≤ q2, then

${\displaystyle \lambda _{k}(D_{1})\leq \lambda _{k}(D_{2}).}$

In fact let T = Gλ for λ large and negative. Then T defines a compact self-adjoint operator on the Hilbert space H. By the spectral theorem for compact self-adjoint operators, H has an orthonormal basis consisting of eigenvectors ψn of T with Tψn = μn ψn, where μn tends to zero. The range of T contains H0 so is dense. Hence 0 is not an eigenvalue of T. The resolvent properties of T imply that ψn lies in H0 and that

D ψn = (λ + 1/μn) ψn

The minimax principle follows because if

${\displaystyle \lambda (G)=\min _{f\perp G}{(Df,f) \over (f,f)},}$

then λ(G)= λk for the linear span of the first k  1 eigenfunctions. For any other (k  1)-dimensional subspace G, some f in the linear span of the first k eigenvectors must be orthogonal to G. Hence λ(G) ≤ (Df,f)/(f,f) ≤ λk.

### Wronskian as a Fredholm determinant

For simplicity, suppose that mq(x) ≤ M on [0,π] with Dirichlet boundary conditions. The minimax principle shows that

${\displaystyle n^{2}+m\leq \lambda _{n}(D)\leq n^{2}+M.}$

It follows that the resolvent (D λ)−1 is a trace-class operator whenever λ is not an eigenvalue of D and hence that the Fredholm determinant det I μ(D λ)−1 is defined.

The Dirichlet boundary conditions imply that

ω(λ)= φλ(b).

Using Picard iteration, Titchmarsh showed that φλ(b), and hence ω(λ), is an entire function of finite order 1/2:

ω(λ) = O(e|λ|)

At a zero μ of ω(λ), φμ(b) = 0. Moreover,

${\displaystyle \psi (x)=\partial _{\lambda }\varphi _{\lambda }(x)|_{\lambda =\mu }}$

satisfies (D  μ)ψ = φμ. Thus

ω(λ) = (λ μ)ψ(b) + O( (λ μ)2).

This implies that [4]

• μ is a simple zero of ω(λ).

For otherwise ψ(b) = 0, so that ψ would have to lie in H0. But then

μ, φμ) = ((D μ)ψ, φμ) = (ψ, (D μ)φμ) = 0,

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 [5]

• ${\displaystyle \omega (\lambda )=C\prod (1-\lambda /\lambda _{n}),}$

for some non-zero constant C.

Hence

${\displaystyle {\rm {det}}\,(I-\mu (D-\lambda )^{-1})=\prod \left(1-{\mu \over \lambda _{n}-\lambda }\right)=\prod {1-(\lambda +\mu )/\lambda _{n} \over 1-\lambda /\lambda _{n}}={\omega (\lambda +\mu ) \over \omega (\lambda )}.}$

In particular if 0 is not an eigenvalue of D

${\displaystyle \omega (\mu )=\omega (0)\cdot {\rm {det}}\,(I-\mu D^{-1}).}$

## Tools from abstract spectral theory

### Functions of bounded variation

A function ρ(x) of bounded variation [6] on a closed interval [a, b] is a complex-valued function such that its total variation V(ρ), the supremum of the variations

${\displaystyle \sum _{r=0}^{k-1}|\rho (x_{r+1})-\rho (x_{r})|}$

over all dissections

${\displaystyle a=x_{0}

is finite. The real and imaginary parts of ρ are real-valued functions of bounded variation. If ρ is real-valued and normalised so that ρ(a)=0, it has a canonical decomposition as the difference of two bounded non-decreasing functions:

${\displaystyle \rho (x)=\rho _{+}(x)-\rho _{-}(x),}$

where ρ+(x) and ρ(x) are the total positive and negative variation of ρ over [a, x].

If f is a continuous function on [a, b] its Riemann–Stieltjes integral with respect to ρ

${\displaystyle \int _{a}^{b}f(x)\,d\rho (x)}$

is defined to be the limit of approximating sums

${\displaystyle \sum _{r=0}^{k-1}f(x_{r})(\rho (x_{r+1})-\rho (x_{r}))}$

as the mesh of the dissection, given by sup |xr+1 - xr|, tends to zero.

This integral satisfies

${\displaystyle \left|\int _{a}^{b}f(x)\,d\rho (x)\right|\leq V(\rho )\cdot \|f\|_{\infty }}$

and thus defines a bounded linear functional dρ on C[a, b] with norm ||dρ||=V(ρ).

Every bounded linear functional μ on C[a, b] has an absolute value |μ| defined for non-negative f by [7]

${\displaystyle |\mu |(f)=\sup _{0\leq |g|\leq f}|\mu (g)|.}$

The form |μ| extends linearly to a bounded linear form on C[a, b] with norm ||μ|| and satisfies the characterizing inequality

|μ(f)| ≤ |μ|(|f|)

for f in C[a, b]. If μ is real, i.e. is real-valued on real-valued functions, then

${\displaystyle \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 g by the formula [8]

${\displaystyle \mu (g)=\lim \mu (f_{n}),}$

where the non-negative continuous functions fn increase pointwise to g.

The same therefore applies to an arbitrary bounded linear form μ, so that a function ρ of bounded variation may be defined by [9]

${\displaystyle \rho (x)=\mu (\chi _{[a,x]}),}$

where χA denotes the characteristic function of a subset A of [a, b]. Thus μ = dρ and ||μ|| = ||dρ||. Moreover μ+ = dρ+ and μ = dρ.

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

The support of μ = dρ is the complement of all points x in [a,b] where ρ is constant on some neighborhood of x; by definition it is a closed subset A of [a,b]. Moreover, μ((1-χA)f) =0, so that μ(f) = 0 if f vanishes on A.

### Spectral measure

Let H be a Hilbert space and ${\displaystyle T}$ a self-adjoint bounded operator on H with ${\displaystyle 0\leq T\leq I}$, so that the spectrum ${\displaystyle \sigma (T)}$ of ${\displaystyle T}$ is contained in ${\displaystyle [0,1]}$. If ${\displaystyle p(t)}$ is a complex polynomial, then by the spectral mapping theorem

${\displaystyle \sigma (p(T))=p(\sigma (T))}$

and hence

${\displaystyle \|p(T)\|\leq \|p\|_{\infty }}$

where ${\displaystyle \|\,\|_{\infty }}$ denotes the uniform norm on ${\displaystyle C[0,1]}$. By the Weierstrass approximation theorem, polynomials are uniformly dense in ${\displaystyle C[0,1]}$. It follows that ${\displaystyle f(T)}$ can be defined ${\displaystyle \forall f\in C[0,1]}$, with

${\displaystyle \sigma (f(T))=f(\sigma (T))}$ and ${\displaystyle \|f(T)\|\leq \|f\|_{\infty }}$.

If ${\displaystyle 0\leq g\leq 1}$ is a lower semicontinuous function on ${\displaystyle [0,1]}$, for example the characteristic function ${\displaystyle \chi _{[0,\alpha ]}}$ of a subinterval of ${\displaystyle [0,1]}$, then ${\displaystyle g}$ is a pointwise increasing limit of non-negative ${\displaystyle f_{n}\in C[0,1]}$.

According to Szőkefalvi-Nagy, [10] if ${\displaystyle \xi }$ is a vector in H, then the vectors

${\displaystyle \eta _{n}=f_{n}(T)\xi }$

form a Cauchy sequence in H, since, for ${\displaystyle n\geq m}$,

${\displaystyle \|\eta _{n}-\eta _{m}\|^{2}\leq (\eta _{n},\xi )-(\eta _{m},\xi ),}$

and ${\displaystyle (\eta _{n},\xi )=(f_{n}(T)\xi ,\xi )}$ is bounded and increasing, so has a limit.

It follows that ${\displaystyle g(T)}$ can be defined by [11]

${\displaystyle g(T)\xi =\lim f_{n}(T)\xi }$.

If ${\displaystyle \xi }$ and ${\displaystyle \eta }$ are vectors in H, then

${\displaystyle \mu _{\xi ,\eta }(f)=(f(T)\xi ,\eta )}$

defines a bounded linear form ${\displaystyle \mu _{\xi ,\eta }}$ on H. By the Riesz representation theorem

${\displaystyle \mu _{\xi ,\eta }=d\rho _{\xi ,\eta }}$

for a unique normalised function ${\displaystyle \rho _{\xi ,\eta }}$ of bounded variation on ${\displaystyle [0,1]}$.

${\displaystyle d\rho _{\xi ,\eta }}$ (or sometimes slightly incorrectly ${\displaystyle \rho _{\xi ,\eta }}$ itself) is called the spectral measure determined by ${\displaystyle \xi }$ and ${\displaystyle \eta }$.

The operator ${\displaystyle g(T)}$ is accordingly uniquely characterised by the equation

${\displaystyle (g(T)\xi ,\eta )=\mu _{\xi ,\eta }(g)=\int _{0}^{1}g(\lambda )\,d\rho _{\xi ,\eta }(\lambda ).}$

The spectral projection ${\displaystyle E(\lambda )}$ is defined by

${\displaystyle E(\lambda )=\chi _{[0,\lambda ]}(T),}$

so that

${\displaystyle \rho _{\xi ,\eta }(\lambda )=(E(\lambda )\xi ,\eta ).}$

It follows that

${\displaystyle g(T)=\int _{0}^{1}g(\lambda )\,dE(\lambda ),}$

which is understood in the sense that for any vectors ${\displaystyle \xi }$ and ${\displaystyle \eta }$,

${\displaystyle (g(T)\xi ,\eta )=\int _{0}^{1}g(\lambda )\,d(E(\lambda )\xi ,\eta )=\int _{0}^{1}g(\lambda )\,d\rho _{\xi ,\eta }(\lambda ).}$

For a single vector ${\displaystyle \xi ,\,\mu _{\xi }=\mu _{\xi ,\xi }}$ is a positive form on ${\displaystyle [0,1]}$ (in other words proportional to a probability measure on ${\displaystyle [0,1]}$) and ${\displaystyle \rho _{\xi }=\rho _{\xi ,\xi }}$ is non-negative and non-decreasing. Polarisation shows that all the forms ${\displaystyle \mu _{\xi ,\eta }}$ can naturally be expressed in terms of such positive forms, since

${\displaystyle \mu _{\xi ,\eta }={\frac {1}{4}}{\bigg (}\mu _{\xi +\eta }+i\mu _{\xi +i\eta }-\mu _{\xi -\eta }-i\mu _{\xi -i\eta }{\bigg )}}$

If the vector ${\displaystyle \xi }$ is such that the linear span of the vectors ${\displaystyle (T^{n}\xi )}$ is dense in H, i.e. ${\displaystyle \xi }$ is a cyclic vector for ${\displaystyle T}$, then the map ${\displaystyle U}$ defined by

${\displaystyle U(f)=f(T)\xi ,\,C[0,1]\rightarrow H}$

satisfies

${\displaystyle (Uf_{1},Uf_{2})=\int _{0}^{1}f_{1}(\lambda ){\overline {f_{2}(\lambda )}}\,d\rho _{\xi }(\lambda ).}$

Let ${\displaystyle L_{2}([0,1],d\rho _{\xi })}$ denote the Hilbert space completion of ${\displaystyle C[0,1]}$ associated with the possibly degenerate inner product on the right hand side. [12] Thus ${\displaystyle U}$ extends to a unitary transformation of ${\displaystyle L_{2}([0,1],\rho _{\xi })}$ onto H. ${\displaystyle UTU^{\ast }}$ is then just multiplication by ${\displaystyle \lambda }$ on ${\displaystyle L_{2}([0,1],d\rho _{\xi })}$; and more generally ${\displaystyle Uf(T)U^{\ast }}$ is multiplication by ${\displaystyle f(\lambda )}$. In this case, the support of ${\displaystyle d\rho _{\xi }}$ is exactly ${\displaystyle \sigma (T)}$, so that

• the self-adjoint operator becomes a multiplication operator on the space of functions on its spectrum with inner product given by the spectral measure.

## Weyl–Titchmarsh–Kodaira theory

The eigenfunction expansion associated with singular differential operators of the form

${\displaystyle Df=-(pf^{\prime })^{\prime }+qf}$

on an open interval (a, b) requires an initial analysis of the behaviour of the fundamental eigenfunctions near the endpoints a and b to determine possible boundary conditions there. Unlike the regular Sturm–Liouville case, in some circumstances spectral values of D can have multiplicity 2. In the development outlined below standard assumptions will be imposed on p and q that guarantee that the spectrum of D 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 D, (D + R )−1 for R large and positive, is given by an operator T corresponding to a Green's function constructed from two fundamental eigenfunctions. In the classical case T was a compact self-adjoint operator; in this case T is just a self-adjoint bounded operator with 0 ≤ T ≤ I. The abstract theory of spectral measure can therefore be applied to T to give the eigenfunction expansion for D.

The central idea in the proof of Weyl and Kodaira can be explained informally as follows. Assume that the spectrum of D lies in [1,∞) and that T =D−1 and let

${\displaystyle E(\lambda )=\chi _{[\lambda ^{-1},1]}(T)}$

be the spectral projection of D corresponding to the interval [1,λ]. For an arbitrary function f define

${\displaystyle f(x,\lambda )=(E(\lambda )f)(x).}$

f(x,λ) may be regarded as a differentiable map into the space of functions of bounded variation ρ; or equivalently as a differentiable map

${\displaystyle x\mapsto (d_{\lambda }f)(x)}$

into the Banach space E of bounded linear functionals dρ on C[α,β] whenever [α,β] is a compact subinterval of [1, ∞).

Weyl's fundamental observation was that dλf satisfies a second order ordinary differential equation taking values in E:

${\displaystyle D(d_{\lambda }f)=\lambda \cdot d_{\lambda }f.}$

After imposing initial conditions on the first two derivatives at a fixed point c, this equation can be solved explicitly in terms of the two fundamental eigenfunctions and the "initial value" functionals

${\displaystyle (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: f(c,λ) and fx(c,λ) may be written as

${\displaystyle f(c,\lambda )=(f,\xi _{1}(\lambda )),\quad f_{x}(c,\lambda )=(f,\xi _{2}(\lambda )),}$

where ξ1(λ) and ξ2(λ) are given purely in terms of the fundamental eigenfunctions. The functions of bounded variation

${\displaystyle \sigma _{ij}(\lambda )=(\xi _{i}(\lambda ),\xi _{j}(\lambda ))}$

determine a spectral measure on the spectrum of D 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 q(x) be a continuous real-valued function on (0,∞) and let D be the second order differential operator

${\displaystyle Df=-f^{\prime \prime }+qf}$

on (0,∞). Fix a point c in (0,∞) and, for λ complex, let ${\displaystyle \varphi _{\lambda },\theta _{\lambda }}$ be the unique fundamental eigenfunctions of D on (0,∞) satisfying

${\displaystyle (D-\lambda )\varphi _{\lambda }=0,\quad (D-\lambda )\theta _{\lambda }=0}$

together with the initial conditions at c

${\displaystyle \varphi _{\lambda }(c)=1,\,\varphi _{\lambda }^{\prime }(c)=0,\,\theta _{\lambda }(c)=0,\,\theta _{\lambda }^{\prime }(c)=1.}$

Then their Wronskian satisfies

${\displaystyle W(\varphi _{\lambda },\theta _{\lambda })=\varphi _{\lambda }\theta _{\lambda }^{\prime }-\theta _{\lambda }\varphi _{\lambda }^{\prime }\equiv 1,}$

since it is constant and equal to 1 at c.

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

${\displaystyle {\rm {Im}}(\lambda )\int _{c}^{x}|\varphi +\mu \theta |^{2}={\rm {Im}}(\mu ).}$

Therefore, the set of ${\displaystyle \mu }$ satisfying this equation is not empty. This set is a circle in the complex ${\displaystyle \mu }$-plane. Points ${\displaystyle \mu }$ in its interior are characterized by

${\displaystyle \int _{c}^{x}|\varphi +\mu \theta |^{2}<{{\rm {Im}}(\mu ) \over {\rm {Im}}(\lambda )}}$

if x > c and by

${\displaystyle \int _{x}^{c}|\varphi +\mu \theta |^{2}<{{\rm {Im}}(\mu ) \over {\rm {Im}}(\lambda )}}$

if x < c.

Let Dx be the closed disc enclosed by the circle. By definition these closed discs are nested and decrease as x approaches 0 or ∞. So in the limit, the circles tend either to a limit circle or a limit point at each end. If ${\displaystyle \mu }$ is a limit point or a point on the limit circle at 0 or ∞, then ${\displaystyle f=\varphi +\mu \theta }$ is square integrable (L2) near 0 or ∞, since ${\displaystyle \mu }$ lies in Dx for all x>c (in the ∞ case) and so ${\displaystyle \int _{c}^{x}|\varphi +\mu \theta |^{2}<{{\rm {Im}}(\mu ) \over {\rm {Im}}(\lambda )}}$ is bounded independent of x. In particular: [13]

• there are always non-zero solutions of Df = λf which are square integrable near 0 resp. ∞;
• in the limit circle case all solutions of Df = λf are square integrable near 0 resp. ∞.

The radius of the disc Dx can be calculated to be

${\displaystyle \left|{1 \over {2{\rm {Im}}(\lambda )\int _{c}^{x}|\theta |^{2}}}\right|}$

and this implies that in the limit point case ${\displaystyle \theta }$ cannot be square integrable near 0 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 0 resp. ∞.

On the other hand, if Dg = λ' g for another value λ', then

${\displaystyle 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 Dh = λh, so that

${\displaystyle 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 (D-λ)g = (λ'-λ)g. Using this to estimate g, it follows that [13]

• the limit point/limit circle behaviour at 0 or ∞ is independent of the choice of λ.

More generally if Dg= (λ – r) g for some function r(x), then [14]

${\displaystyle 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 [14]

• if r is continuous at 0, then D + r is limit point or limit circle at 0 precisely when D is,

so that in particular [15]

• if q(x)- a/x2 is continuous at 0, then D is limit point at 0 if and only if a ≥ ¾.

Similarly

• if r has a finite limit at ∞, then D + r is limit point or limit circle at ∞ precisely when D is,

so that in particular [16]

• if q has a finite limit at ∞, then D 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

${\displaystyle D_{0}f=-(p_{0}f^{\prime })^{\prime }+q_{0}f}$

on (0,∞) with q0 positive and continuous on (0,∞) and p0 continuously differentiable in [0,∞), positive in (0,∞) and p0(0)=0.

Moreover, assume that after reduction to standard form D0 becomes the equivalent operator

${\displaystyle Df=-f^{\prime \prime }+qf}$

on (0,∞) where q has a finite limit at ∞. Thus

• D is limit point at ∞.

At 0, D may be either limit circle or limit point. In either case there is an eigenfunction Φ0 with DΦ0=0 and Φ0 square integrable near 0. In the limit circle case, Φ0 determines a boundary condition at 0:

${\displaystyle W(f,\Phi _{0})(0)=0.}$

For λ complex, let Φλ and Χλ satisfy

• (D – λ)Φλ = 0, (D – λ)Χλ = 0
• Χλ square integrable near infinity
• Φλ square integrable at 0 if 0 is limit point
• Φλ satisfies the boundary condition above if 0 is limit circle.

Let

${\displaystyle \omega (\lambda )=W(\Phi _{\lambda },\mathrm {X} _{\lambda }),}$

a constant which vanishes precisely when Φλ and Χλ are proportional, i.e. λ is an eigenvalue of D for these boundary conditions.

On the other hand, this cannot occur if Im λ ≠ 0 or if λ is negative. [13]

Indeed, if D f= λf with q0 – λ ≥ δ >0, then by Green's formula (Df,f) = (f,Df), since W(f,f*) is constant. So λ must be real. If f is taken to be real-valued in the D0 realization, then for 0 < x < y

${\displaystyle [p_{0}ff^{\prime }]_{x}^{y}=\int _{x}^{y}(q_{0}-\lambda )|f|^{2}+p_{0}(f^{\prime })^{2}.}$

Since p0(0) = 0 and f is integrable near 0, p0ff ' must vanish at 0. Setting x = 0, it follows that f(y) f '(y) >0, so that f2 is increasing, contradicting the square integrability of f near ∞.

Thus, adding a positive scalar to q, it may be assumed that

ω(λ) ≠ 0 when λ is not in [1,∞).

If ω(λ) ≠ 0, the Green's function Gλ(x,y) at λ is defined by

${\displaystyle G_{\lambda }(x,y)=\Phi _{\lambda }(x)\mathrm {X} _{\lambda }(y)/\omega (\lambda )\,\,(x\leq y),\,\,\,\,\mathrm {X} _{\lambda }(x)\Phi _{\lambda }(y)/\omega (\lambda )\,\,(x\geq y).}$

and is independent of the choice of λ and Χλ.

In the examples there will be a third "bad" eigenfunction Ψλ defined and holomorphic for λ not in [1, ∞) such that Ψλ satisfies the boundary conditions at neither 0 nor ∞. This means that for λ not in [1, ∞)

• Wλλ) is nowhere vanishing;
• Wλλ) is nowhere vanishing.

In this case Χλ is proportional to Φλ + m(λ) Ψλ, where

• m(λ) = – Wλλ) / Wλλ).

Let H1 be the space of square integrable continuous functions on (0,∞) and let H0 be

• the space of C2 functions f on (0,∞) of compact support if D is limit point at 0
• the space of C2 functions f on (0,∞) with W(f0)=0 at 0 and with f = 0 near ∞ if D is limit circle at 0.

Define T = G0 by

${\displaystyle (Tf)(x)=\int _{0}^{\infty }G_{0}(x,y)f(y)\,dy.}$

Then TD = I on H0, DT = I on H1 and the operator D is bounded below on H0:

${\displaystyle (Df,f)\geq (f,f).}$

Thus T is a self-adjoint bounded operator with 0 ≤ TI.

Formally T = D−1. The corresponding operators Gλ defined for λ not in [1,∞) can be formally identified with

${\displaystyle (D-\lambda )^{-1}=T(I-\lambda T)^{-1}}$

and satisfy Gλ (D – λ) = I on H0, (D – λ)Gλ = I on H1.

### Spectral theorem and Titchmarsh–Kodaira formula

Theorem. [13] [17] [18] For every real number λ let ρ(λ) be defined by theTitchmarsh–Kodaira formula:

${\displaystyle \rho (\lambda )=\lim _{\delta \downarrow 0}\lim _{\varepsilon \downarrow 0}{1 \over \pi }\int _{\delta }^{\lambda +\delta }{\rm {Im}}\,m(t+i\varepsilon )\,dt.}$

Then ρ(λ) is a lower semicontinuous non-decreasing function of λ and if

${\displaystyle (Uf)(\lambda )=\int _{0}^{\infty }f(x)\Phi (x,\lambda )\,dx,}$

then U defines a unitary transformation of L2(0,∞) onto L2([1,∞), dρ) such that UDU−1corresponds to multiplication by λ.

The inverse transformation U−1 is given by

${\displaystyle (U^{-1}g)(x)=\int _{1}^{\infty }g(\lambda )\Phi (x,\lambda )\,d\rho (\lambda ).}$

The spectrum of D equals the support of dρ.

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

In fact for T =D−1 with 0 ≤ TI, the spectral projection E(λ) of T is defined by

${\displaystyle E(\lambda )=\chi _{[\lambda ^{-1},1]}(T)}$

It is also the spectral projection of D corresponding to the interval [1,λ].

For f in H1 define

${\displaystyle f(x,\lambda )=(E(\lambda )f)(x).}$

f(x,λ) may be regarded as a differentiable map into the space of functions ρ of bounded variation; or equivalently as a differentiable map

${\displaystyle x\mapsto (d_{\lambda }f)(x)}$

into the Banach space E of bounded linear functionals dρ on C[α,β] for any compact subinterval [α,β] of [1, ∞).

The functionals (or measures) dλf(x) satisfies the following E-valued second order ordinary differential equation:

${\displaystyle D(d_{\lambda }f)=\lambda \cdot d_{\lambda }f,}$

with initial conditions at c in (0,∞)

${\displaystyle (d_{\lambda }f)(c)=d_{\lambda }f(c,\cdot )=\mu ^{(0)},\quad (d_{\lambda }f)^{\prime }(c)=d_{\lambda }f_{x}(c,\cdot )=\mu ^{(1)}.}$

If φλ and χλ are the special eigenfunctions adapted to c, then

${\displaystyle d_{\lambda }f(x)=\varphi _{\lambda }(x)\mu ^{(0)}+\chi _{\lambda }(x)\mu ^{(1)}.}$

Moreover,

${\displaystyle \mu ^{(k)}=d_{\lambda }(f,\xi _{\lambda }^{(k)}),}$

where

${\displaystyle \xi _{\lambda }^{(k)}=DE(\lambda )\eta ^{(k)},}$

with

${\displaystyle \eta _{z}^{(0)}(y)=G_{z}(c,y),\,\,\,\,\eta _{z}^{(1)}(x)=\partial _{x}G_{z}(c,y),\,\,\,\,(z\notin [1,\infty )).}$

(As the notation suggests, ξλ(0) and ξλ(1) do not depend on the choice of z.)

Setting

${\displaystyle \sigma _{ij}(\lambda )=(\xi _{\lambda }^{(i)},\xi _{\lambda }^{(j)}),}$

it follows that

${\displaystyle d_{\lambda }(E(\lambda )\eta _{z}^{(i)},\eta _{z}^{(j)})=|\lambda -z|^{-2}\cdot d_{\lambda }\sigma _{ij}(\lambda ).}$

On the other hand, there are holomorphic functions a(λ), b(λ) such that

• φλ + a(λ) χλ is proportional to Φλ;
• φλ + b(λ) χλ is proportional to Χλ.

Since Wλλ) = 1, the Green's function is given by

${\displaystyle G_{\lambda }(x,y)={(\varphi _{\lambda }(x)+a(\lambda )\chi _{\lambda }(x))(\varphi _{\lambda }(y)+b(\lambda )\chi _{\lambda }(y)) \over b(\lambda )-a(\lambda )}\,\,(x\leq y),\,\,\,\,{(\varphi _{\lambda }(x)+b(\lambda )\chi _{\lambda }(x))(\varphi _{\lambda }(y)+a(\lambda )\chi _{\lambda }(y)) \over b(\lambda )-a(\lambda )}\,\,(y\leq x).}$

Direct calculation [22] shows that

${\displaystyle (\eta _{z}^{(i)},\eta _{z}^{(j)})={\rm {Im}}\,M_{ij}(z)/{\rm {Im}}\,z,}$

where the so-called characteristic matrixMij(z) is given by

${\displaystyle M_{00}(z)={a(z)b(z) \over a(z)-b(z)},\,\,M_{01}(z)=M_{10}(z)={a(z)+b(z) \over 2(a(z)-b(z))},\,\,M_{11}(z)={1 \over a(z)-b(z)}.}$

Hence

${\displaystyle \int _{-\infty }^{\infty }({\rm {Im}}\,z)\cdot |\lambda -z|^{-2}\,d\sigma _{ij}(\lambda )={\rm {Im}}M_{ij}(z),}$

which immediately implies

${\displaystyle \sigma _{ij}(\lambda )=\lim _{\delta \downarrow 0}\lim _{\varepsilon \downarrow 0}\int _{\delta }^{\lambda +\delta }{\rm {Im}}\,M_{ij}(t+i\varepsilon )\,dt.}$

(This is a special case of the "Stieltjes inversion formula".)

Setting ψλ(0)λ and ψλ(1)λ, it follows that

${\displaystyle (E(\mu )f)(x)=\sum _{i.j}\int _{0}^{\mu }\int _{0}^{\infty }\psi _{\lambda }^{(i)}(x)\psi _{\lambda }^{(j)}(y)f(y)\,dy\,d\sigma _{ij}(\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 [23] [24] [25] concerns the eigenfunction expansion associated with the Legendre differential operator D

${\displaystyle Df=-((x^{2}-1)f^{\prime })^{\prime }=-(x^{2}-1)f^{\prime \prime }-2xf^{\prime }}$

on (1,∞). The eigenfunctions are the Legendre functions [26]

${\displaystyle P_{-1/2+i{\sqrt {\lambda }}}(\cosh r)={1 \over 2\pi }\int _{0}^{2\pi }\left({\sin \theta +ie^{-r}\cos \theta \over \cos \theta -ie^{-r}\sin \theta }\right)^{{1 \over 2}+i{\sqrt {\lambda }}}\,d\theta }$

with eigenvalue λ ≥ 0. The two Mehler–Fock transformations are [27]

${\displaystyle Uf(\lambda )=\int _{1}^{\infty }f(x)\,P_{-1/2+i{\sqrt {\lambda }}}(x)\,dx}$

and

${\displaystyle U^{-1}g(x)=\int _{0}^{\infty }g(\lambda )\,{1 \over 2}\tanh \pi {\sqrt {\lambda }}\,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, [28] consider the group G = SU(1,1) consisting of complex matrices of the form

${\displaystyle \left({\begin{matrix}\alpha &\beta \\{\overline {\beta }}&{\overline {\alpha }}\end{matrix}}\right)}$

with determinant |α|2 |β|2 = 1.

## Generalisations and alternative approaches

A Weyl function can be defined at a singular endpoint ${\displaystyle a}$ giving rise to a singular version of Weyl–Titchmarsh–Kodaira theory. [29] this applies for example to the case of radial Schrödinger operators

${\displaystyle Df=-f^{\prime \prime }+{\frac {l(l+1)}{x^{2}}}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. [30]

## Notes

1. Titchmarsh 1962 , p. 22
2. Dieudonné 1969, Chapter X.
3. Titchmarsh, E.C. (1939), Theory of Functions, Oxford University Press, §8.2.
4. Burkill, J.C. (1951), The Lebesgue Integral, Cambridge Tracts in Mathematics and Mathematical Physics, 40, Cambridge University Press, pp. 50–52, ISBN   978-0-521-04382-3
5. Loomis, Lynn H. (1953), An Introduction to Abstract Harmonic Analysis, van Nostrand, page 40.
6. Loomis 1953 , pp. 30–31
7. Kolmogorov, A.N.; Fomin, S.V. (1975), , Dover, pp.  374–376, ISBN   978-0-486-61226-3
8. Riesz & Szőkefalvi-Nagy 1990 , p. 263
9. This is a limit in the strong operator topology.
10. A bona fide inner product is defined on the quotient by the subspace of null functions ${\displaystyle f}$, i.e. those with ${\displaystyle \mu _{\xi }(|f|_{2})=0}$. Alternatively in this case the support of the measure is ${\displaystyle \sigma (T)}$, so the right hand side defines a (non-degenerate) inner product on ${\displaystyle C(\sigma (T))}$.
11. Weyl 1910 harvnb error: multiple targets (2×): CITEREFWeyl1910 (help)
12. Bellman 1969 , p. 116
13. Reed & Simon 1975 , p. 159
14. Reed & Simon 1975 , p. 154
15. Titchmarsh 1946, Chapter III.
16. Kodaira 1949 , pp. 935–936
17. Kodaira 1949 , pp. 929–932; for omitted details, see Kodaira 1950 , pp. 529–536
18. Stone 1932, Chapter X.
19. Kodaira 1950 , pp. 534–535
20. Mehler, F.G. (1881), "Ueber mit der Kugel- und Cylinderfunctionen verwandte Function und ihre Anwendung in der Theorie der Elektricitätsverteilung", Mathematische Annalen , 18 (2): 161–194, doi:10.1007/BF01445847
21. Fock, V.A. (1943), "On the representation of an arbitrary function by an integral involving Legendre's functions with a complex index", C. R. Acad. Sci. URSS, 39: 253–256
22. Terras, Audrey (1984), "Non-Euclidean harmonic analysis, the central limit theorem, and long transmission lines with random inhomogeneities", J. Multivariate Anal., 15 (2): 261–276, doi:
23. Lebedev, N.N. (1972), Special Functions and Their Applications, Dover, ISBN   978-0-486-60624-8
24. Vilenkin 1968, Chapter VI.
25. Kostenko, Aleksey; Sakhnovich, Alexander; Teschl, Gerald (2012), "Weyl–Titchmarsh Theory for Schrödinger Operators with Strongly Singular Potentials", Int Math Res Notices, 2012: 1699–1747, arXiv:, doi:10.1093/imrn/rnr065
26. Eckhardt, Jonathan; Teschl, Gerald (2013), "Sturm–Liouville operators with measure-valued coefficients", J. d'Analyse Math., 120: 151–224, arXiv:, doi:10.1007/s11854-013-0018-x

## Related Research Articles

In mathematics, a self-adjoint operator on a finite-dimensional complex vector space V with inner product is a linear map A that is its own adjoint: for all vectors v and w. If V is finite-dimensional with a given orthonormal basis, this is equivalent to the condition that the matrix of A is a Hermitian matrix, i.e., equal to its conjugate transpose A. By the finite-dimensional spectral theorem, V has an orthonormal basis such that the matrix of A relative to this basis is a diagonal matrix with entries in the real numbers. In this article, we consider generalizations of this concept to operators on Hilbert spaces of arbitrary dimension.

In mathematics, particularly in functional analysis, a projection-valued measure (PVM) is a function defined on certain subsets of a fixed set and whose values are self-adjoint projections on a fixed Hilbert space. Projection-valued measures are formally similar to real-valued measures, except that their values are self-adjoint projections rather than real numbers. As in the case of ordinary measures, it is possible to integrate complex-valued functions with respect to a PVM; the result of such an integration is a linear operator on the given Hilbert space.

In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

In the theory of general relativity, linearized gravity is the application of perturbation theory to the metric tensor that describes the geometry of spacetime. As a consequence, linearized gravity is an effective method for modeling the effects of gravity when the gravitational field is weak. The usage of linearized gravity is integral to the study of gravitational waves and weak-field gravitational lensing.

In physics, the Majorana equation is a relativistic wave equation. It is named after the Italian physicist Ettore Majorana, who proposed it in 1937 as a means of describing fermions that are their own antiparticle. Particles corresponding to this equation are termed Majorana particles, although that term now has a more expansive meaning, referring to any fermionic particle that is its own anti-particle.

In mathematics, specifically in symplectic geometry, the momentum map is a tool associated with a Hamiltonian action of a Lie group on a symplectic manifold, used to construct conserved quantities for the action. The momentum map generalizes the classical notions of linear and angular momentum. It is an essential ingredient in various constructions of symplectic manifolds, including symplectic (Marsden–Weinstein) quotients, discussed below, and symplectic cuts and sums.

In physics and fluid mechanics, a Blasius boundary layer describes the steady two-dimensional laminar boundary layer that forms on a semi-infinite plate which is held parallel to a constant unidirectional flow. Falkner and Skan later generalized Blasius' solution to wedge flow, i.e. flows in which the plate is not parallel to the flow.

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

Alternatives to general relativity are physical theories that attempt to describe the phenomenon of gravitation in competition to Einstein's theory of general relativity. There have been many different attempts at constructing an ideal theory of gravity.

In mathematics, the Prolate spheroidal wave functions (PSWF) are a set of orthogonal bandlimited functions. They are eigenfunctions of a timelimiting operation followed by a lowpassing operation. Let denote the time truncation operator, such that if and only if is timelimited within . Similarly, let denote an ideal low-pass filtering operator, such that if and only if is bandlimited within . The operator turns out to be linear, bounded and self-adjoint. For we denote with the n-th eigenfunction, defined as

In the Standard Model, using quantum field theory it is conventional to use the helicity basis to simplify calculations. In this basis, the spin is quantized along the axis in the direction of motion of the particle.

In mathematics, the Plancherel theorem for spherical functions is an important result in the representation theory of semisimple Lie groups, due in its final form to Harish-Chandra. It is a natural generalisation in non-commutative harmonic analysis of the Plancherel formula and Fourier inversion formula in the representation theory of the group of real numbers in classical harmonic analysis and has a similarly close interconnection with the theory of differential equations. It is the special case for zonal spherical functions of the general Plancherel theorem for semisimple Lie groups, also proved by Harish-Chandra. The Plancherel theorem gives the eigenfunction expansion of radial functions for the Laplacian operator on the associated symmetric space X; it also gives the direct integral decomposition into irreducible representations of the regular representation on L2(X). In the case of hyperbolic space, these expansions were known from prior results of Mehler, Weyl and Fock.

In fluid dynamics, a cnoidal wave is a nonlinear and exact periodic wave solution of the Korteweg–de Vries equation. These solutions are in terms of the Jacobi elliptic function cn, which is why they are coined cnoidal waves. They are used to describe surface gravity waves of fairly long wavelength, as compared to the water depth.

In mathematical physics, the Belinfante–Rosenfeld tensor is a modification of the energy–momentum tensor that is constructed from the canonical energy–momentum tensor and the spin current so as to be symmetric yet still conserved.

Symmetries in quantum mechanics describe features of spacetime and particles which are unchanged under some transformation, in the context of quantum mechanics, relativistic quantum mechanics and quantum field theory, and with applications in the mathematical formulation of the standard model and condensed matter physics. In general, symmetry in physics, invariance, and conservation laws, are fundamentally important constraints for formulating physical theories and models. In practice, they are powerful methods for solving problems and predicting what can happen. While conservation laws do not always give the answer to the problem directly, they form the correct constraints and the first steps to solving a multitude of problems.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

Functional principal component analysis (FPCA) is a statistical method for investigating the dominant modes of variation of functional data. Using this method, a random function is represented in the eigenbasis, which is an orthonormal basis of the Hilbert space L2 that consists of the eigenfunctions of the autocovariance operator. FPCA represents functional data in the most parsimonious way, in the sense that when using a fixed number of basis functions, the eigenfunction basis explains more variation than any other basis expansion. FPCA can be applied for representing random functions, or in functional regression and classification.

The generalized functional linear model (GFLM) is an extension of the generalized linear model (GLM) that allows one to regress univariate responses of various types on functional predictors, which are mostly random trajectories generated by a square-integrable stochastic processes. Similarly to GLM, a link function relates the expected value of the response variable to a linear predictor, which in case of GFLM is obtained by forming the scalar product of the random predictor function with a smooth parameter function . Functional Linear Regression, Functional Poisson Regression and Functional Binomial Regression, with the important Functional Logistic Regression included, are special cases of GFLM. Applications of GFLM include classification and discrimination of stochastic processes and functional data.

In the theory of Lie groups, Lie algebras and their representation theory, a Lie algebra extensione is an enlargement of a given Lie algebra g by another Lie algebra h. Extensions arise in several ways. There is the trivial extension obtained by taking a direct sum of two Lie algebras. Other types are the split extension and the central extension. Extensions may arise naturally, for instance, when forming a Lie algebra from projective group representations. Such a Lie algebra will contain central charges.

In astrophysics, Chandrasekhar's white dwarf equation is an initial value ordinary differential equation introduced by the Indian American astrophysicist Subrahmanyan Chandrasekhar, in his study of the gravitational potential of completely degenerate white dwarf stars. The equation reads as