Differential operator

From HandWiki
Short description: Typically linear operator defined in terms of differentiation of functions
A harmonic function defined on an annulus. Harmonic functions are exactly those functions which lie in the kernel of the Laplace operator, an important differential operator.

In mathematics, a differential operator is an operator defined as a function of the differentiation operator. It is helpful, as a matter of notation first, to consider differentiation as an abstract operation that accepts a function and returns another function (in the style of a higher-order function in computer science).

This article considers mainly linear differential operators, which are the most common type. However, non-linear differential operators also exist, such as the Schwarzian derivative.

Definition

Given a nonnegative integer m, an order-[math]\displaystyle{ m }[/math] linear differential operator is a map [math]\displaystyle{ P }[/math] from a function space [math]\displaystyle{ \mathcal{F}_1 }[/math] to another function space [math]\displaystyle{ \mathcal{F}_2 }[/math] that can be written as:

[math]\displaystyle{ P = \sum_{|\alpha|\le m}a_\alpha(x) D^\alpha\ , }[/math] where [math]\displaystyle{ \alpha = (\alpha_1,\alpha_2,\cdots,\alpha_n) }[/math] is a multi-index of non-negative integers, [math]\displaystyle{ |\alpha| = \alpha_1 + \alpha_2 + \cdots + \alpha_n }[/math], and for each [math]\displaystyle{ \alpha }[/math], [math]\displaystyle{ a_\alpha(x) }[/math] is a function on some open domain in n-dimensional space. The operator [math]\displaystyle{ D^\alpha }[/math] is interpreted as

[math]\displaystyle{ D^\alpha = \frac{\partial^{|\alpha|}}{\partial x_1^{\alpha_1}\partial x_2^{\alpha_2}\cdots\partial x_n^{\alpha_n}} }[/math]

Thus for a function [math]\displaystyle{ f \in \mathcal{F}_1 }[/math]:

[math]\displaystyle{ P f = \sum_{|\alpha|\le m}a_\alpha(x) \frac{\partial^{|\alpha|}f}{\partial x_1^{\alpha_1}\partial x_2^{\alpha_2}\cdots\partial x_n^{\alpha_n}} }[/math]

The notation [math]\displaystyle{ D^{\alpha} }[/math] is justified (i.e., independent of order of differentiation) because of the symmetry of second derivatives.

The polynomial p obtained by replacing D by variables [math]\displaystyle{ \xi }[/math] in P is called the total symbol of P; i.e., the total symbol of P above is: [math]\displaystyle{ p(x, \xi) = \sum_{|\alpha|\le m}a_\alpha(x) \xi^\alpha }[/math] where [math]\displaystyle{ \xi^\alpha = \xi_1^{\alpha_1} \cdots \xi_n^{\alpha_n}. }[/math] The highest homogeneous component of the symbol, namely,

[math]\displaystyle{ \sigma(x, \xi) = \sum_{|\alpha|= m}a_\alpha(x) \xi^\alpha }[/math]

is called the principal symbol of P. While the total symbol is not intrinsically defined, the principal symbol is intrinsically defined (i.e., it is a function on the cotangent bundle).[1]

More generally, let E and F be vector bundles over a manifold X. Then the linear operator

[math]\displaystyle{ P: C^\infty(E) \to C^\infty(F) }[/math]

is a differential operator of order [math]\displaystyle{ k }[/math] if, in local coordinates on X, we have

[math]\displaystyle{ Pu(x) = \sum_{|\alpha| = k} P^\alpha(x) \frac {\partial^\alpha u} {\partial x^{\alpha}} + \text{lower-order terms} }[/math]

where, for each multi-index α, [math]\displaystyle{ P^\alpha(x):E \to F }[/math] is a bundle map, symmetric on the indices α.

The kth order coefficients of P transform as a symmetric tensor

[math]\displaystyle{ \sigma_P: S^k (T^*X) \otimes E \to F }[/math]

whose domain is the tensor product of the kth symmetric power of the cotangent bundle of X with E, and whose codomain is F. This symmetric tensor is known as the principal symbol (or just the symbol) of P.

The coordinate system xi permits a local trivialization of the cotangent bundle by the coordinate differentials dxi, which determine fiber coordinates ξi. In terms of a basis of frames eμ, fν of E and F, respectively, the differential operator P decomposes into components

[math]\displaystyle{ (Pu)_\nu = \sum_\mu P_{\nu\mu}u_\mu }[/math]

on each section u of E. Here Pνμ is the scalar differential operator defined by

[math]\displaystyle{ P_{\nu\mu} = \sum_{\alpha} P_{\nu\mu}^\alpha\frac{\partial}{\partial x^\alpha}. }[/math]

With this trivialization, the principal symbol can now be written

[math]\displaystyle{ (\sigma_P(\xi)u)_\nu = \sum_{|\alpha|=k} \sum_{\mu}P_{\nu\mu}^\alpha(x)\xi_\alpha u_\mu. }[/math]

In the cotangent space over a fixed point x of X, the symbol [math]\displaystyle{ \sigma_P }[/math] defines a homogeneous polynomial of degree k in [math]\displaystyle{ T^*_x X }[/math] with values in [math]\displaystyle{ \operatorname{Hom}(E_x, F_x) }[/math].

Fourier interpretation

A differential operator P and its symbol appear naturally in connection with the Fourier transform as follows. Let ƒ be a Schwartz function. Then by the inverse Fourier transform,

[math]\displaystyle{ Pf(x) = \frac{1}{(2\pi)^{\frac{d}{2}}} \int\limits_{\mathbf{R}^d} e^{ ix\cdot\xi} p(x,i\xi)\hat{f}(\xi)\, d\xi. }[/math]

This exhibits P as a Fourier multiplier. A more general class of functions p(x,ξ) which satisfy at most polynomial growth conditions in ξ under which this integral is well-behaved comprises the pseudo-differential operators.

Examples

  • The differential operator [math]\displaystyle{ P }[/math] is elliptic if its symbol is invertible; that is for each nonzero [math]\displaystyle{ \theta \in T^*X }[/math] the bundle map [math]\displaystyle{ \sigma_P (\theta, \dots, \theta) }[/math] is invertible. On a compact manifold, it follows from the elliptic theory that P is a Fredholm operator: it has finite-dimensional kernel and cokernel.
  • In the study of hyperbolic and parabolic partial differential equations, zeros of the principal symbol correspond to the characteristics of the partial differential equation.
  • In applications to the physical sciences, operators such as the Laplace operator play a major role in setting up and solving partial differential equations.
  • In differential topology, the exterior derivative and Lie derivative operators have intrinsic meaning.
  • In abstract algebra, the concept of a derivation allows for generalizations of differential operators, which do not require the use of calculus. Frequently such generalizations are employed in algebraic geometry and commutative algebra. See also Jet (mathematics).
  • In the development of holomorphic functions of a complex variable z = x + iy, sometimes a complex function is considered to be a function of two real variables x and y. Use is made of the Wirtinger derivatives, which are partial differential operators: [math]\displaystyle{ \frac{\partial}{\partial z} = \frac{1}{2} \left( \frac{\partial}{\partial x} - i \frac{\partial}{\partial y} \right) \ ,\quad \frac{\partial}{\partial\bar{z}}= \frac{1}{2} \left( \frac{\partial}{\partial x} + i \frac{\partial}{\partial y} \right) \ . }[/math] This approach is also used to study functions of several complex variables and functions of a motor variable.
  • The differential operator del, also called nabla, is an important vector differential operator. It appears frequently in physics in places like the differential form of Maxwell's equations. In three-dimensional Cartesian coordinates, del is defined as
[math]\displaystyle{ \nabla = \mathbf{\hat{x}} {\partial \over \partial x} + \mathbf{\hat{y}} {\partial \over \partial y} + \mathbf{\hat{z}} {\partial \over \partial z}. }[/math]
Del defines the gradient, and is used to calculate the curl, divergence, and Laplacian of various objects.

History

The conceptual step of writing a differential operator as something free-standing is attributed to Louis François Antoine Arbogast in 1800.[2]

Notations

The most common differential operator is the action of taking the derivative. Common notations for taking the first derivative with respect to a variable x include:

[math]\displaystyle{ {d \over dx} }[/math], [math]\displaystyle{ D }[/math], [math]\displaystyle{ D_x, }[/math] and [math]\displaystyle{ \partial_x }[/math].

When taking higher, nth order derivatives, the operator may be written:

[math]\displaystyle{ {d^n \over dx^n} }[/math], [math]\displaystyle{ D^n }[/math], [math]\displaystyle{ D^n_x }[/math], or [math]\displaystyle{ \partial_x^n }[/math].

The derivative of a function f of an argument x is sometimes given as either of the following:

[math]\displaystyle{ [f(x)]' }[/math]
[math]\displaystyle{ f'(x). }[/math]

The D notation's use and creation is credited to Oliver Heaviside, who considered differential operators of the form

[math]\displaystyle{ \sum_{k=0}^n c_k D^k }[/math]

in his study of differential equations.

One of the most frequently seen differential operators is the Laplacian operator, defined by

[math]\displaystyle{ \Delta = \nabla^2 = \sum_{k=1}^n \frac{\partial^2}{\partial x_k^2}. }[/math]

Another differential operator is the Θ operator, or theta operator, defined by[3]

[math]\displaystyle{ \Theta = z {d \over dz}. }[/math]

This is sometimes also called the homogeneity operator, because its eigenfunctions are the monomials in z: [math]\displaystyle{ \Theta (z^k) = k z^k,\quad k=0,1,2,\dots }[/math]

In n variables the homogeneity operator is given by [math]\displaystyle{ \Theta = \sum_{k=1}^n x_k \frac{\partial}{\partial x_k}. }[/math]

As in one variable, the eigenspaces of Θ are the spaces of homogeneous functions. (Euler's homogeneous function theorem)

In writing, following common mathematical convention, the argument of a differential operator is usually placed on the right side of the operator itself. Sometimes an alternative notation is used: The result of applying the operator to the function on the left side of the operator and on the right side of the operator, and the difference obtained when applying the differential operator to the functions on both sides, are denoted by arrows as follows:

[math]\displaystyle{ f \overleftarrow{\partial_x} g = g \cdot \partial_x f }[/math]
[math]\displaystyle{ f \overrightarrow{\partial_x} g = f \cdot \partial_x g }[/math]
[math]\displaystyle{ f \overleftrightarrow{\partial_x} g = f \cdot \partial_x g - g \cdot \partial_x f. }[/math]

Such a bidirectional-arrow notation is frequently used for describing the probability current of quantum mechanics.

Adjoint of an operator

Given a linear differential operator [math]\displaystyle{ T }[/math] [math]\displaystyle{ Tu = \sum_{k=0}^n a_k(x) D^k u }[/math] the adjoint of this operator is defined as the operator [math]\displaystyle{ T^* }[/math] such that [math]\displaystyle{ \langle Tu,v \rangle = \langle u, T^*v \rangle }[/math] where the notation [math]\displaystyle{ \langle\cdot,\cdot\rangle }[/math] is used for the scalar product or inner product. This definition therefore depends on the definition of the scalar product (or inner product).

Formal adjoint in one variable

In the functional space of square-integrable functions on a real interval (a, b), the scalar product is defined by [math]\displaystyle{ \langle f, g \rangle = \int_a^b \overline{f(x)} \,g(x) \,dx , }[/math]

where the line over f(x) denotes the complex conjugate of f(x). If one moreover adds the condition that f or g vanishes as [math]\displaystyle{ x \to a }[/math] and [math]\displaystyle{ x \to b }[/math], one can also define the adjoint of T by [math]\displaystyle{ T^*u = \sum_{k=0}^n (-1)^k D^k \left[ \overline{a_k(x)} u \right]. }[/math]

This formula does not explicitly depend on the definition of the scalar product. It is therefore sometimes chosen as a definition of the adjoint operator. When [math]\displaystyle{ T^* }[/math] is defined according to this formula, it is called the formal adjoint of T.

A (formally) self-adjoint operator is an operator equal to its own (formal) adjoint.

Several variables

If Ω is a domain in Rn, and P a differential operator on Ω, then the adjoint of P is defined in L2(Ω) by duality in the analogous manner:

[math]\displaystyle{ \langle f, P^* g\rangle_{L^2(\Omega)} = \langle P f, g\rangle_{L^2(\Omega)} }[/math]

for all smooth L2 functions f, g. Since smooth functions are dense in L2, this defines the adjoint on a dense subset of L2: P* is a densely defined operator.

Example

The Sturm–Liouville operator is a well-known example of a formal self-adjoint operator. This second-order linear differential operator L can be written in the form

[math]\displaystyle{ Lu = -(pu')'+qu=-(pu''+p'u')+qu=-pu''-p'u'+qu=(-p) D^2 u +(-p') D u + (q)u. }[/math]

This property can be proven using the formal adjoint definition above.[4]

This operator is central to Sturm–Liouville theory where the eigenfunctions (analogues to eigenvectors) of this operator are considered.

Properties of differential operators

Differentiation is linear, i.e.

[math]\displaystyle{ D(f+g) = (Df)+(Dg), }[/math]
[math]\displaystyle{ D(af) = a(Df), }[/math]

where f and g are functions, and a is a constant.

Any polynomial in D with function coefficients is also a differential operator. We may also compose differential operators by the rule

[math]\displaystyle{ (D_1 \circ D_2)(f) = D_1(D_2(f)). }[/math]

Some care is then required: firstly any function coefficients in the operator D2 must be differentiable as many times as the application of D1 requires. To get a ring of such operators we must assume derivatives of all orders of the coefficients used. Secondly, this ring will not be commutative: an operator gD isn't the same in general as Dg. For example we have the relation basic in quantum mechanics:

[math]\displaystyle{ Dx - xD = 1. }[/math]

The subring of operators that are polynomials in D with constant coefficients is, by contrast, commutative. It can be characterised another way: it consists of the translation-invariant operators.

The differential operators also obey the shift theorem.

Ring of polynomial differential operators

Ring of univariate polynomial differential operators

Main page: Weyl algebra

If R is a ring, let [math]\displaystyle{ R\langle D,X \rangle }[/math] be the non-commutative polynomial ring over R in the variables D and X, and I the two-sided ideal generated by DXXD − 1. Then the ring of univariate polynomial differential operators over R is the quotient ring [math]\displaystyle{ R\langle D,X\rangle/I }[/math]. This is a non-commutative simple ring. Every element can be written in a unique way as a R-linear combination of monomials of the form [math]\displaystyle{ X^a D^b \text{ mod } I }[/math]. It supports an analogue of Euclidean division of polynomials.

Differential modules[clarification needed] over [math]\displaystyle{ R[X] }[/math] (for the standard derivation) can be identified with modules over [math]\displaystyle{ R\langle D,X\rangle/I }[/math].

Ring of multivariate polynomial differential operators

If R is a ring, let [math]\displaystyle{ R\langle D_1,\ldots,D_n,X_1,\ldots,X_n\rangle }[/math] be the non-commutative polynomial ring over R in the variables [math]\displaystyle{ D_1,\ldots,D_n,X_1,\ldots,X_n }[/math], and I the two-sided ideal generated by the elements

[math]\displaystyle{ (D_i X_j-X_j D_i)-\delta_{i,j},\ \ \ D_i D_j -D_j D_i,\ \ \ X_i X_j - X_j X_i }[/math]

for all [math]\displaystyle{ 1 \le i,j \le n, }[/math] where [math]\displaystyle{ \delta }[/math] is Kronecker delta. Then the ring of multivariate polynomial differential operators over R is the quotient ring [math]\displaystyle{ R\langle D_1,\ldots,D_n,X_1,\ldots,X_n\rangle/I }[/math].

This is a non-commutative simple ring. Every element can be written in a unique way as a R-linear combination of monomials of the form [math]\displaystyle{ X_1^{a_1} \ldots X_n^{a_n} D_1^{b_1} \ldots D_n^{b_n} }[/math].

Coordinate-independent description

In differential geometry and algebraic geometry it is often convenient to have a coordinate-independent description of differential operators between two vector bundles. Let E and F be two vector bundles over a differentiable manifold M. An R-linear mapping of sections P : Γ(E) → Γ(F) is said to be a kth-order linear differential operator if it factors through the jet bundle Jk(E). In other words, there exists a linear mapping of vector bundles

[math]\displaystyle{ i_P: J^k(E) \to F }[/math]

such that

[math]\displaystyle{ P = i_P\circ j^k }[/math]

where jk: Γ(E) → Γ(Jk(E)) is the prolongation that associates to any section of E its k-jet.

This just means that for a given section s of E, the value of P(s) at a point x ∈ M is fully determined by the kth-order infinitesimal behavior of s in x. In particular this implies that P(s)(x) is determined by the germ of s in x, which is expressed by saying that differential operators are local. A foundational result is the Peetre theorem showing that the converse is also true: any (linear) local operator is differential.

Relation to commutative algebra

An equivalent, but purely algebraic description of linear differential operators is as follows: an R-linear map P is a kth-order linear differential operator, if for any k + 1 smooth functions [math]\displaystyle{ f_0,\ldots,f_k \in C^\infty(M) }[/math] we have

[math]\displaystyle{ [f_k,[f_{k-1},[\cdots[f_0,P]\cdots]]=0. }[/math]

Here the bracket [math]\displaystyle{ [f,P]:\Gamma(E)\to \Gamma(F) }[/math] is defined as the commutator

[math]\displaystyle{ [f,P](s)=P(f\cdot s)-f\cdot P(s). }[/math]

This characterization of linear differential operators shows that they are particular mappings between modules over a commutative algebra, allowing the concept to be seen as a part of commutative algebra.

Variants

A differential operator of infinite order

A differential operator of infinite order is (roughly) a differential operator whose total symbol is a power series instead of a polynomial.

Bidifferential operator

A differential operator acting on two functions [math]\displaystyle{ D(g,f) }[/math] is called a bidifferential operator. The notion appears, for instance, in an associative algebra structure on a deformation quantization of a Poisson algebra.[5]

Microdifferential operator

A microdifferential operator is a type of operator on an open subset of a cotangent bundle, as opposed to an open subset of a manifold. It is obtained by extending the notion of a differential operator to the cotangent bundle.[6]

See also


References

  1. Schapira 1985, 1.1.7
  2. James Gasser (editor), A Boole Anthology: Recent and classical studies in the logic of George Boole (2000), p. 169; Google Books.
  3. E. W. Weisstein. "Theta Operator". http://mathworld.wolfram.com/ThetaOperator.html. 
  4. [math]\displaystyle{ \begin{align} L^*u & {} = (-1)^2 D^2 [(-p)u] + (-1)^1 D [(-p')u] + (-1)^0 (qu) \\ & {} = -D^2(pu) + D(p'u)+qu \\ & {} = -(pu)''+(p'u)'+qu \\ & {} = -p''u-2p'u'-pu''+p''u+p'u'+qu \\ & {} = -p'u'-pu''+qu \\ & {} = -(pu')'+qu \\ & {} = Lu \end{align} }[/math]
  5. Omori, Hideki; Maeda, Y.; Yoshioka, A. (1992). "Deformation quantization of Poisson algebras" (in en). Proceedings of the Japan Academy, Series A, Mathematical Sciences 68 (5). doi:10.3792/PJAA.68.97. 
  6. Schapira 1985, § 1.2. § 1.3.

Further reading

  • Fedosov, Boris; Schulze, Bert-Wolfgang; Tarkhanov, Nikolai (2002). "Analytic index formulas for elliptic corner operators" (in en). Annales de l'Institut Fourier 52 (3): 899–982. doi:10.5802/aif.1906. ISSN 1777-5310. 

External links