research33 min studyLesson 829 in Mathematics

Propagation of singularities

Follow microlocal nonnegligibility along real Hamilton trajectories and distinguish residual wavefront propagation from weaker measure transport.

On this page

What you’ll understand

  • State finite-phase-space real-principal-type propagation away from a source.
  • Separate elliptic regularity, characteristic flow and the strength of residual assumptions.
  • Compute a phase graph and a translated wavefront under an exact transport equation.

Microlocal irregularity does not move arbitrarily through phase space. For an equation with real principal symbol, the Hamiltonian flow organizes where a nonnegligible oscillation can persist. Elliptic regions exclude such oscillations unless the source supplies them, while characteristic trajectories carry them between microlocal neighbourhoods. The theorem is sensitive to the residual scale: a statement about all powers of h is stronger than a limiting statement about quadratic mass.

Separate elliptic points from characteristic trajectories

Let P_h be a self-adjoint semiclassical pseudodifferential operator with smooth real principal symbol p, in a controlled local symbol calculus. Fix a real energy E and an h-tempered family u_h. Where p−E is nonzero, a microlocal parametrix shows that u_h is residual if (P_h−E)u_h is residual there. Thus away from the source wavefront, WF_h(u_h) can lie only on the characteristic set p=E. This is elliptic regularity and does not yet describe transport along that set.

At real-principal-type points, dp is nonzero on p=E, so H_p is a nonvanishing characteristic vector field. On an open phase-space region where the source is O(h^∞), the wavefront set of u_h is a union of segments of H_p trajectories as long as they stay in that region. Nonnegligibility cannot simply stop in the middle of a source-free segment. Fixed cutoffs, polynomial h bounds and a suitable pseudodifferential calculus are part of the hypothesis, rather than consequences of a ray drawing.

(Ph−E)uh=O(h∞) microlocally,dp≠0 on p=E⟹WF⁡h(uh)⊂{p=E} and propagates along Hp(P_h-E)u_h=O(h^\infty)\ \text{microlocally},\quad dp\ne0\text{ on }p=E\quad\Longrightarrow\quad\operatorname{WF}_h(u_h)\subset\{p=E\}\text{ and propagates along }H_p
The statement is local at finite phase-space points and away from source wavefront. The operator has real principal type and the family is h-tempered.↔ Scroll wide equations horizontally.

The residual scale determines the conclusion

The positive-commutator method compares a localized norm before and after motion along a characteristic. Its proof uses the sign of H_p acting on a carefully selected test symbol, then estimates lower-order and source terms. This is not an assertion that the waveform itself remains unchanged along a trajectory. Amplitudes, phases and transverse spreading can evolve while microlocal nonnegligibility persists. At turning or radial configurations, the chosen local description and hypotheses must be checked rather than replacing the theorem by an unrestricted propagation slogan.

An O(h) residual is not O(h^∞). For an L²-bounded normalized quasimode, a stronger o(h) residual under self-adjointness and suitable localization gives invariance of a semiclassical measure through the limiting commutator identity. That measure-level conclusion does not automatically determine every point of the wavefront set. A polynomially weak component may disappear from the measure and remain in wavefront support. It is therefore necessary to state whether one is transporting a residual regularity estimate, a leading quadratic measure, or only a formal phase.

∥(Ph−E)uh∥=o(h)⟹∫Hpa dμ=0(a∈Cc∞)\|(P_h-E)u_h\|=o(h)\quad\Longrightarrow\quad\int H_pa\,d\mu=0\qquad(a\in C_c^\infty)
This measure-invariance implication uses a self-adjoint realization, an L²-bounded family, the commutator calculus and appropriate localization. It is weaker than full O(h^∞) wavefront propagation.↔ Scroll wide equations horizontally.

Cutoffs create legitimate sources

A smooth spatial cutoff can make an exact oscillatory solution square integrable, but derivatives of that cutoff create a source. Those terms are often only O(h), not residual. If an oscillation disappears at the edge of the cutoff, this does not contradict source-free propagation because the edge is precisely where the cutoff commutator acts. Inside a region where the cutoff equals one, the exact equation and its characteristic interpretation remain valid. Keeping these regions separate avoids a common false counterexample.

For a first-order transport operator with linear momentum symbol, the solution translates exactly. Its support shifts at the group velocity, while its momentum stays fixed. A finite sine-wave lab can illustrate this speed only for its own dispersion relation. The displayed wave sin(2x−t) has speed one half, and should not be silently identified with a different Hamiltonian whose flow speed is one. The operator symbol determines the characteristic velocity, not a generic animation clock.

A phase graph preserved by a real characteristic flow

Locally consider P_h=hDₓ+x at E=0 and u_h=exp(−ix²/(2h)). Verify the equation, locate its phase graph, and check Hamiltonian invariance. Explain what a compact cutoff changes.

  1. The derivative is u′=−ixu/h, so hDₓu=−xu and (hDₓ+x)u=0 exactly. This is a local h-tempered solution, though its modulus one is not globally square integrable on the line.

  2. Its phase derivative is −x, giving the characteristic graph ξ=−x. The symbol p=ξ+x has Hamilton equations xdot=1 and ξdot=−1, so x(t)+ξ(t) stays zero and trajectories remain on the graph.

  3. Multiplying by a compact smooth χ gives residual −ihχ′u. This source is supported where the cutoff changes and is generally only O(h). Where χ equals one, the original exact equation and source-free propagation still apply.

The result

The exact phase graph follows the characteristic flow. A cutoff's endpoint behaviour is explained by its source term rather than a failure of propagation.

Explore the model

Interactive lab

Interactive lab · waves

The model's own characteristic speed

Pause and Reset, then set phase difference = 0° and ratio = 0. After two Steps, t=1 and the first wave is sin(2x−1), with crests shifted right by 1/2. Reset and set both phase difference = 0° and ratio = 1, then take two Steps: the sum is twice that wave and has the same shift. These phases match the transport symbol p=ξ/2 and speed one half. The instrument does not implement a general pseudodifferential propagation theorem or a semiclassical residual estimate.

The model's own characteristic speedTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forChanging amplitude does not change this transport velocity; the symbol and phase relation fix the speed.

Check your understanding

A cutoff exact solution has residual −ihχ′u near its endpoints. What follows?

Try it yourself

Solve ih∂ₜu=(hDₓ/2)u for initial state χ(x)exp(3ix/h). Track the finite wavefront at time t=2 if χ is concentrated near x=1. State the symbol and speed.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The equation is ∂ₜu=−(1/2)∂ₓu. Translate the entire initial function, including its phase, not just the envelope.

Compare with the solution

The exact solution is u(t,x)=χ(x−t/2)exp⁡(3i(x−t/2)/h)u(t,x)=\chi(x-t/2)\exp(3i(x-t/2)/h). The symbol is p=ξ/2, so xdot=1/2 and ξdot=0. At t=2, the envelope is concentrated near x=2 and the momentum remains three; the finite wavefront is (support χ+1) times {3}. The extra time-dependent scalar phase does not change this support. This is an exact transport example with known group velocity, rather than a claim that the same translation formula solves a general Schrödinger operator with nonlinear dispersion.

Take it with you

  • Elliptic regularity removes wavefront away from the characteristic set and source.
  • Real-principal-type propagation follows Hamilton trajectories under residual microlocal assumptions.
  • Cutoff sources, measure invariance and full wavefront propagation require distinct estimates.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.