In this blog post I will just record some things I’ve been trying to learn about lately, largely just so I can have a place to collect my thoughts. Most of this is in Hörmander’s monograph on differential operators, and is motivated by trying to understand Vasy’s method and Atiyah-Singer index theory.

**Pseudodifferential operators on manifolds.**

Let us recall that a symbol on an open subset X of is by definition a smooth function on the cotangent bundle of X (for which certain seminorms are finite). This was curious to me — you can motivate it by saying that a symbol is an observable and the cotangent bundle is “phase space” in the sense that a point consists of a position x and a momentum , but why should the momentum live in a cotangent space and not the fiber of some other vector bundle? When we quantize a symbol a, defining an operator a(D) by formally substituting the differential operator in place of the momentum, we by definition obtain a pseudodifferential operator. Now let be a diffeomorphism, and introduce the pushforward symbol . This is the “right” definition in the sense that .

If a is a symbol of order m, then modulo symbols of order m – 1. But is invariantly defined as an isomorphism of tangent bundles , so its transpose should be an isomorphism of the dual bundle. This only makes sense if is a covector at y.

The above paragraphs are totally obvious, and yet puzzled me for the past three years, until last week when I sat down and decided to work out the details for myself.

The consequence is that we cannot define the symbol of a pseudodifferential operator invariantly. Rather, we declare that a pseudodifferential operator A has the property that for every chart and every pair of cutoffs on Y, then the operator is a pseudodifferential operator on Y (in the sense that it is the quantization of a symbol on Y; here the pushforward is defined to be the inverse of the pullback ). Since Y is an open subset of this makes sense.

Previously we have discussed pseudodifferential operators on manifolds M. These can be viewed more abstractly as acting on sections of the trivial line bundle . However, in geometry one frequently has to deal with sections of more general vector bundles over M. For example, a 1-form is a section of the cotangent bundle. If E, F are vector bundles over M of rank r, s respectively, one may define the Hom-bundle Hom(E, F), which locally is isomorphic to the matrix bundle . Then a pseudodifferential operator from sections of E to sections of F is nothing more than a linear map which, after trivialization of E and F, looks like a $s \times r$ matrix of pseudodifferential operators on M. The principal symbol of such an operator sends the cotangent bundle of M into the Hom-bundle Hom(E, F).

**Wavefront sets.**

In this section we will impose that all pseudodifferential operators have Schwartz kernels K such that the projections of supp K are both proper maps. Modulo the space of pseudodifferential operators of order , this assumption is no loss of generality. Under this assumption, the top-order term of a symbol — that is, the principal symbol — satisfies the pushforward formula , so the principal symbol is well-defined as an element of (here is the th symbol class). The principal symbol encodes important information about the nature of the operator; for example we have:

*Definition*. An *elliptic pseudodifferential operator* of order m is one whose principal symbol is near infinity of each cotangent space.

The important property is that if A is an elliptic pseudodifferential operator, then A is also invertible modulo the quantization of . For example the Laplace-Beltrami operator is elliptic on Riemannian manifolds since its symbol is ; since the quadratic form induced by a Lorentzian metric is not positive-definite, it follows that on Lorentzian manifolds, the Laplace-Beltrami operator is not elliptic. Since a Lorentzian Laplace-Beltrami operator is really just the d’Alembertian, whose symbol is , this should be no surprise.

Recall that a conic set in a vector space is a set which is closed under multiplication by conic scalars. A conic set in a vector bundle, then, is one which is conic in every fiber.

*Definition*. Let a be the principal symbol of a pseudodifferential operator A of order m. We say that A is *noncharacteristic* near if there is a conic neighborhood of wherein near infinity. Otherwise, we say that is a *characteristic point*. The set of characteristic points is denoted Char A and the set of noncharacteristic points is denoted Ell A.

Thus a pseudodifferential operator A is noncharacteristic at if in a neighborhood of x, A is elliptic when restricted to the direction . By definition, Char A is closed, so we may make the following definition.

*Definition*. Let u be a distribution. The *wavefront set* WF(u) is the intersection of all sets Char A, where A ranges over pseudodifferential operators such that .

Then WF(u) is a closed conic subset of the cotangent bundle , and its projection to M is exactly the singular support ss(u). Indeed, iff for every pseudodifferential operator A in a sufficiently small neighborhood of x, ; in other words no matter how hard we try, we cannot force u to become singular without differentiating it away from x. The wavefront set also remembers the direction in which this singularity happens; by elliptic invertibility, it will not happen in a direction that A is noncharacteristic.

For example, the only way that can be made smooth is by cutting off u to away from , which can be done by pseudodifferential operators of order 0 which are elliptic in the x-direction, but not possibly in the y-direction, along the x-axis.

**Pseudotransport equations.**

Hyperbolic operators are meant to generalize the transport equation . Let us therefore begin by studying the “pseudotransport” equation .

We assume that is uniformly bounded in and continuous in , and the real part of a is uniformly bounded from below. Then we have the energy estimate

valid for any and large enough depending on s. Applying the Hanh-Banach theorem we conclude that for every initial data in we can find which solves the pseudotransport equation. In particular, given Schwartz initial data, it follows that u is smooth.

Now fix initial data and assume that the principal symbol exists and is imaginary. (This forces the transport operator to be real and of order 1.) Let q be a symbol of order 0 on space, with principal symbol . If in fact Q(D) is a pseudodifferential operator on spacetime such that such at time 0, Q(0) = q, and Q(t, D) commutes with then Qu solves the pseudotransport equation. (Actually, we will find Q so that is a pseudodifferential operator of order ; this is good enough.) In particular if then WF(u) is contained in Char Q, and WF(u) should be the intersection of all such sets Char Q.

To compute WF(u), let be the principal symbol of a(D) and suppose that , where is principal, is given. Then the principal symbol of is the Poisson bracket

where is the Hamilton vector field of a symbol p. By inducting on j, we can use this computation to compute and conclude that modulo an error term of order , we can choose Q to be invariant along the Hamiltonian flow given by the Hamiltonian . That is, if , then . This result is a sort of “propagation of singularities” for the pseudotransport equation, which generalizes the fact that the transport equation acts on Dirac masses by transporting them, as expected.

**Solving the hyperbolic Cauchy problem.**

Let X be a manifold that represents “spacetime”. A priori we may not have a Lorentzian metric to work with, so instead we fix a function that is a “time coordinate”. The level surfaces of can be viewed as “spacelike hypersurfaces” in X.

Throughout we will let and denote the present and future, respectively.

*Definition*. A *hyperbolic operator* is a differential operator P of principal symbol p and order m such that and for every such that is not in the span of , there are m distinct such that .

Since P is a differential operator, p(x) is a homogeneous polynomial of order m. To make sense of the condition, let me restrict to the case that with its usual Riemannian metric and is the projection onto the t-axis. Then after rotating the first coordinate so that is a covector dual to the x-axis, the condition says that given we can find exactly m real numbers such that . In the case of the d’Alembertian, we have , and indeed given we can set .

To state the initial-value problem with initial data in the “initial-time slice” , let v be a vector field such that , so v points “forward in time”. The action of v is “differentiating with respect to time”. Note that this hypothesis prevents from degenerating.

*Theorem* (solving the hyperbolic Cauchy problem). Let P be a hyperbolic operator of order m with smooth coefficients, Y a precompact open submanifold of X, and . Assume we are given an inhomogeneous term satisfying and initial data , j < m. Then there is supported in such that Pu = f in and in .

The proof is in Chapter 23.2 of Hörmander. The idea is to first prove uniqueness of solutions. By compactness, we may cover Y with finitely many charts U which are isomorphic to open subsets of Minkowski spacetime in which level sets of are spacelike hypersurfaces and orbits of v are worldlines. Since Minkowski spacetime has an honest-to-god time coordinate, the hyperbolicity hypothesis allows us to factor the principal symbol p into first-order factors, and hence factor P into pseudotransport operators on U, at least modulo a lower-order error. We may then apply the solution of the Cauchy problem for pseudotransport operators to solve the Cauchy problem for Pu = f in each chart U, and since there were only finitely many, uniqueness allows us to stitch the local solutions together into a global solution.

The proof outlined in the above paragraph is motivated by the special case when P is the d’Alembertian, which already appears in Chapter 2 of Evans. In that proof, one first observes that the Cauchy problem for the transport equation has an explicit solution. Then one reduces to the case that spacetime is two-dimensional, in which case there is an explicit factorization of P into transport operators, namely .

**Propagation of singularities, part I.**

To study the propagation of singularities we need to recall some symplectic geometry. Let Q be a pseudodifferential operator on X and q its principal symbol. Then the Hamilton vector field induces a flow on which preserves q.

*Definition*. The *bicharacteristic flow* of a pseudodifferential operator Q of principal symbol q is the flow of on . A *bicharacteristic* of Q is an orbit of the bicharacteristic flow.

The intuition for the bicharacteristic flow is that its projection to X is “lightlike”, at least if Q is the d’Alembertian.

*Theorem* (Hörmander’s propagation of singularities). Let P be a pseudodifferential operator of order m such that the Schwartz kernel of P has proper support, and the principal symbol of P is real. Then for every distribution u, WF(u) – WF(f) is invariant under the bicharacteristic flow of P.

By definition of the wavefront set, for every distribution u, WF(u) – WF(Qu) is contained in Char Q. But if Q is a differential operator, then Char Q is exactly the “characteristic variety” , which is exactly the variety where the bicharacteristic flow of Q is defined. Therefore we can ask that WF(u) – WF(Qu) be invariant under the bicharacteristic flow.

If P is a hyperbolic operator of principal symbol p, then the solutions of the equation are all real and distinct, and modulo lower-order terms this can be used to enforce that the coefficients of p are real. We phrase this more simply by saying that the principal symbol of every hyperbolic operator is real.

A partial converse to the reality of principal symbols of hyperbolic operators holds. If Q is a differential operator, then its principal symbol q is a homogeneous polynomial on each cotangent space. Fixing a particular cotangent space, we can write where ranges over all multiindices of order m and . In order that the characteristic variety of Q have more than one real point, there must be some positive and some negative. But this is exactly the situation of the d’Alembertian, whose principal symbol is .

Thus, while the propagation of singularities theorem only assumes that the principal symbol is real, if the operator P is (for example) elliptic or parabolic, then the conclusion of the theorem is degenerate in the sense that the characteristic variety only has a single real point, so that WF(u) – WF(f) is invariant under EVERY group action on the characteristic variety, not just the bicharacteristic flow.

The interpretation of the propagation of singularities theorem is that P is something like the d’Alembertian, in which case p is something like a Lorentzian metric. The bicharacteristic flow is a flow on the characteristic bundle, which is the space whose points consist of a position x and a lightlike momentum . Therefore the projection of any bicharacteristic to X consists of a worldline. Thus, if the initial data is something like a Dirac mass at x, then the Dirac mass travels along the worldline containing x.

To prove the propagation of singularities theorem, we need a propagation estimate. Recall that if A is a pseudodifferential operator, then WF(A) denotes the microsupport of A; that is, the complement of the largest conic set on which A has order .

*Theorem* (propagation estimate). Let U be an open conic set, and let . Let P be a pseudodifferential operator of real principal symbol p and order m.

For every N > 0 and there is C > 0 such that for every distribution u and every inhomogeneous term f with Pu = f,

given that the following criteria are met:

- The projection of U is precompact in X.
- For every , if , then and the radial vector field are linearly independent at .
- WF(A) and WF(B) are contained in U, while .
- For every trajectory of with , there is T < 0 such that for every , and .

The term is an error term created by the use of pseudodifferential operators and is not interesting. The operator is a cutoff which microlocalizes the problem to a neighborhood to the conic set U. We are interested in WF(u) – WF(f), so we want and . Actually, since we only care about the complement of WF(f), we might as well take f Schwartz, in which case we can take and simplify the propagation estimate to

The interesting point here is the relationship between the operators A and B. We can optimize the propagation estimate by assuming that WF(B) = Ell B. This is because we really desperately want B to be elliptic on its microsupport, so that it does not introduce any new singularities. Under the assumption WF(B) = Ell B, B is a microlocalization to WF(B), and if , then got to WF(A) after passing through WF(B). The point is that if u has a singularity at , then (if the regularity exponent s is taken large enough) , but we assumed f Schwartz, so this implies , so that if we traveled back along the bicharacteristic flow from for long enough, we would see that u already had a singularity at some time with T < 0.

Moreover, the propagation estimate is time-reversible in the sense we can replace T < 0 with -T > 0. Thus the bicharacteristic flow neither creates nor destroys singularities in the distribution u. This readily implies the propagation of singularities theorem.

The proof of the propagation estimate is quite technical and this post is meant as a more of a conceptual discussion so I will omit it.