AN-04 · CC0; linked components retain their own terms
Global time and the bicharacteristic relation
A compact characteristic tube supplies coordinates along one finite segment. A global propagation construction needs compatible geometry along every segment. The compact return condition in the preceding global solvability lesson gives precisely that geometry: a smooth space of trajectories, one section meeting each trajectory, and an invariant positive radial variable. We prove the global statements before using them in a parametrix.
All manifolds here are Hausdorff and second countable, smooth and without boundary. The full finite-coordinate flow proofs NF1–NF7 include existence, uniqueness, all parameter derivatives, maximal continuation, the flow law and coordinate gluing. The finite inverse and implicit proofs P2–P3 supply every local inverse used below. PS5 constructs compact exhaustions and smooth subordinate partitions on the actual Hausdorff second-countable manifolds. The choice and Baire proofs include the countable selection used in the dimension argument.
Compactness and finite calculus, integration with all compact parameters, real powers, logarithms and exponentials, and the full trigonometric proofs supply the scalar calculations. The Hamilton identities and normalized characteristic lifts give the exact symplectic convention and all-real-order reduction. We supply the global quotient, affine averaging, invariant-radius and relation arguments here; their proof map connects the exact earlier proofs.
Self-checked by the writing AI.
1. Three descriptions of a global flow
Let v be a smooth vector field on a manifold M, with maximal flow ϕt(y) on its open domain D⊂M×R. An integral curve is complete here when it is taken on its maximal interval. A stationary curve is included among periodic curves. Consider the following conditions.
(A) No complete integral curve is relatively compact. For each compact K⊂M, there is a compact K′⊂M containing every finite curve interval with both endpoints in K.
(B) There are no periodic integral curves, and the same-trajectory relation
R={(ϕt(y),y):(y,t)∈D}(FG1)
is a closed embedded smooth submanifold of M×M.
(C) There is a smooth manifold M0 and an open set Ω⊂M0×R containing M0×{0}, such that each fiber I(q)={t:(q,t)∈Ω} is an interval, and a diffeomorphism M≃Ω carries v to ∂t.
Theorem 1.1. Conditions (A), (B) and (C) are equivalent. Every nonempty component satisfying them has dimension d≥1; its trajectory space has dimension d−1, and its relation has dimension d+1. An empty manifold satisfies the conditions vacuously.
We will also use the stronger explicit consequence: under any of these conditions, the map
F:D⟶M×M,F(y,t)=(ϕt(y),y)(FG2)
is a proper injective immersion, and is a diffeomorphism onto R.
2. Compact return gives a proper relation
Under (A), v has no zero and no periodic orbit. Both would give a complete curve in a compact set. There is also no relatively compact positive or negative half-curve. To prove this, first use the compact ODE continuation argument to extend a trapped half-curve indefinitely in its outward time. For times tj→+∞, choose ϕtj(y)→z. For each finite T, the shifted intervals [tj−T,tj+T] lie eventually in the half-curve. Compact bounds and continuous dependence give the solution through z on [−T,T], still in the same compact set. Unique solutions for different T agree. Their union is a complete compact curve, contradicting (A). Use −v for the negative half.
The differential of (FG2) is injective. Its second component determines the variation in y, and the remaining time variation is a multiple of the nonzero vector v(ϕt(y)). The map itself is injective because equality of two times on one trajectory would give a periodic orbit.
We prove properness, including the maximal-domain issue. Suppose
yj⟶y,ϕtj(yj)⟶x.(FG3)
The endpoint sequences and their limits lie in one compact K. Every intervening interval lies in the compact K′ from (A). After taking a subsequence, tj has a limit in the extended real line. If that limit is +∞, compact ODE continuation and continuous dependence show that the positive half-curve through y exists for every finite positive time and remains in K′. This contradicts the preceding half-curve argument. The negative infinite limit is excluded in the same way.
Thus tj→T∈R. The entire bounded-time interval through y remains in K′; uniform short-time rectangles covering K′ continue it across any supposed finite maximal endpoint. Consequently (y,T)∈D, and continuity gives ϕT(y)=x. Every sequence in the inverse image of a compact set therefore has a convergent subsequence whose limit is still in D. The compact target does supply convergent endpoint subsequences: cover it by finitely many coordinate neighborhoods whose closures are compact inside larger charts, and use finite-dimensional sequential compactness on one closure containing infinitely many sequence terms. The same local argument applies on any of our manifolds. The sequence together with its limit is compact directly from the open-cover definition, so the endpoint sets used above really lie in a compact K. This proves compactness directly: second countability gives a countable subcover of any open cover of that inverse image; if no finite initial subfamily covered it, points chosen outside the first j members would have a convergent subsequence, contradicting eventual membership in a cover member containing its limit. This is properness.
A proper injective immersion is an embedding in this case without any extra quotient assumption. Around each domain point, an invertible derivative minor supplies local immersion coordinates. Properness prevents a sequence of other domain points from entering its image neighborhood while escaping those coordinates: a convergent image sequence has a convergent inverse subsequence, and injectivity identifies its only possible limit. Thus the inverse on the image is continuous, and the local immersion inverses make it smooth. The image is closed by the same compact subsequence argument. This proves (A) implies (B) and the stated property of F.
For completeness, (B) also gives this property. Under (B), the absence of periodic and stationary curves is an assumption. The injectivity and differential-rank calculation therefore applies independently of (A); the argument in this paragraph does not use the return bound it will prove. Write d=dimM. The smooth map F:D→R has rank d+1. The dimension of R is exactly d+1, rather than an unstated additional hypothesis in (B). If it were larger on an open component, select at each source point a square derivative minor of full source dimension. The inverse theorem makes those target components coordinates on a small source neighborhood; the remaining target components become smooth functions of them. Its image is therefore a graph of positive codimension in the larger-dimensional relation chart. Compact subpieces have compact images with empty interior, hence nowhere dense. Countably many such compact subpieces cover the relevant open subset of D, by its countable coordinate exhaustion. Their images cannot cover an open manifold. More precisely, restrict to one nonempty relation chart of the proposed larger dimension. Its inverse image under F is an open second-countable manifold. Refine its countable precompact coordinate cover so that each compact closure maps into one of the full-minor graph charts just constructed. The restriction of F is smooth into the embedded relation: in an ambient submanifold chart the normal coordinates of its image are identically zero, and its remaining coordinates are smooth. Each compact graph image is closed in the relation chart and has empty interior there; a transverse target coordinate can be varied off its graph in every neighborhood. These images give precisely the countable closed nowhere-dense cover excluded in the next argument. Here the needed Baire argument is elementary: in a coordinate ball, successively choose nonempty closed smaller balls avoiding the first, second and subsequent nowhere dense sets, with radii tending to zero and each contained in the preceding interior. The first ball is compact; the nested balls have a common point, which avoids every set. This proves the contradiction. The opposite dimension inequality already follows from the injective differential.
The inverse function theorem now makes the bijection F:D→R a local diffeomorphism, hence a diffeomorphism. Because R is closed in M×M, F is proper as a map to that product. A relatively compact complete curve would extend for all real time by compact continuation, and F(y,tj), with ∣tj∣→∞, would contradict properness. For compact K, put E=F−1(K×K), compact in D. The set
{(y,θt):(y,t)∈E, 0≤θ≤1}(FG4)
is compact, is contained in D by the interval property of maximal times, and its flow image is compact. That image contains every return interval in question. This proves (B) implies (A).
3. The quotient and one global section
Assume the equivalent conditions just proved. Give M0=M/R its quotient topology, with projection π. Saturation of an open subset of M is open: it is the union of its images under the local flow diffeomorphisms. Thus π is open. It follows that the images of a countable base in M form a countable base in M0.
The quotient is Hausdorff. If y and z belong to different trajectories, closedness of R gives neighborhoods U,V with (U×V)∩R=∅. The open sets π(U),π(V) are disjoint and separate their classes.
Choose a small smooth transversal S through y. The flow box follows from the smooth local flow and the inverse theorem: the derivative of (s,t)↦ϕt(s) at (y,0) has the transversal tangent space and v(y) as complementary columns. Shrink the box so that π∣S is injective. Indeed, the continuous inverse of (FG2) near (y,y) forces any relation pair in S×S to have a sufficiently small time, and in a flow box such a curve meets S only once. The image π(S) is open because it is the image of the open flow box. The same reasoning with smaller transversal subsets proves that π∣S is a homeomorphism onto its image.
These charts form a smooth atlas of dimension d−1. To check a transition, fix the finite time taking one representative to the other transversal. Apply that fixed-time smooth flow to nearby representatives, then use the second flow box's smooth transversal projection. It gives exactly the representative on the second transversal, by the uniqueness just proved. The reversed construction is its smooth inverse. The quotient is therefore a Hausdorff, second-countable smooth manifold, and π is locally the projection of a flow box. Each quotient chart supplies a smooth local section σi.
For the partition in the next paragraph, choose the small balls in PS5 inside the given section domains; its construction then makes each closed support lie inside one such domain. If M0 is zero dimensional, it is a countable discrete manifold. Choose one origin on each trajectory and use its singleton as a section domain; the singleton weights give the smooth locally finite partition directly. This includes the one-dimensional components of M.
Each trajectory has an intrinsic affine time coordinate: choosing any point as origin gives its open maximal interval, and changing the origin translates time. Differences of times are unique because there is no periodic orbit. Let {χi} be the smooth nonnegative locally finite partition of unity on M0, with supports inside the local section domains. Define the section σ(q) to be the affine weighted average of the finitely many σi(q) with nonzero weights.
This averaging is meaningful even for an incomplete trajectory. In any time coordinate on that trajectory all the relevant points lie in its open interval. Their convex combination lies in the same interval. Translation of the origin changes each time and its average by the same amount, so it gives the same point of the trajectory. Near a fixed q, use one local section as origin. The other section times are smooth by the smooth inverse of F, and only finitely many partition supports occur. Their weighted sum is smooth; terms extend by zero outside their support domains. Thus σ is a global smooth section, with πσ=id. At a point outside a section domain, the closed support of its partition weight has a neighborhood disjoint from that point. Its weighted term is identically zero there. This is why extension by zero remains smooth even if the unweighted time function has no limit at the boundary of its domain.
Put
Ω={(q,t):(σ(q),t)∈D},G(q,t)=ϕt(σ(q)).(FG5)
The domain is open, contains time zero, and has interval fibers. Every point lies on one selected trajectory, so G is onto. Unique times and unique trajectory classes make it injective. Locally its transversal columns and flow column are independent, so it is a local diffeomorphism. Its local inverses agree by injectivity, giving a global smooth inverse. It carries ∂t to v. This proves (A) and (B) imply (C).
Finally, (C) implies (A). A complete vertical trajectory cannot be relatively compact: an infinite time endpoint makes its time coordinate unbounded; a finite maximal endpoint approached inside a compact subset of Ω has a limit in that open set and therefore permits continuation. For the return condition, pairs of endpoints in a compact K⊂Ω with the same q form a closed subset of K×K, hence a compact set. Their interval interpolation (q,(1−θ)t1+θt2), 0≤θ≤1, has compact image contained in Ω by fiber convexity. It is the required K′. Theorem 1.1 is proved in full.
4. A radial coordinate constant along the flow
Let M now be a smooth conic manifold: the positive scaling action is free and has smooth local product coordinates, with Hausdorff second-countable quotient Ms. Let v commute with this action, and suppose its induced field on Ms satisfies Theorem 1.1. Choose a smooth positive degree-one function r on M. Such a function exists: local positive radial coordinates transform by positive smooth degree-zero factors; a subordinate smooth partition on Ms makes their locally finite positive weighted sum a global degree-one function. In a cotangent cone one can simply use a positive smooth cotangent norm.
The map y↦(πs(y),r(y)) is an equivariant diffeomorphism to Ms×R+: every positive scaling orbit meets r=1 exactly once, with smooth inverse in the local products. By Theorem 1.1 write Ms≃Ω⊂M0×R. In these coordinates scaling invariance gives
v=∂t+a(q,t)r∂r,b(q,t)=−∫0ta(q,s)ds,R=reb(q,t).(FG6)
The interval from zero to t stays in Ω. Integration with parameters makes b smooth. Direct differentiation gives vR=0, while R stays positive and has degree one. The inverse is r=Re−b(q,t), so this is an equivariant diffeomorphism, and in the resulting coordinates
M≃Ω×R+,v=∂t.(FG7)
No completeness in the original radial variable is assumed. The transformation is defined on the full actual interval fibers and retains their finite endpoints.
5. Transport on the whole conic manifold
Let c be smooth, complex and homogeneous of degree zero. In (FG7) it is independent of R. For any real μ, define Sμ(M) by the local estimates in conic product coordinates, on every compact set of (q,t) and for R≥1:
∣∂q,tα∂Rkf∣≤Cα,kRμ−k.(FG8)
These estimates are invariant under conic coordinate changes. To see every derivative, write the old coordinates as z=z(z′), R=h(z′)R′, where h>0 and z=(q,t). On a fixed compact set, h,h−1 and all derivatives needed in a given estimate are bounded. Each R′-derivative supplies one old radial derivative and a bounded factor h. A derivative in z′ either differentiates a bounded coefficient, supplies an old base derivative, or supplies an old radial derivative together with one factor R′. Repeated differentiation therefore gives finite terms of the form (R′)ℓA(z′)∂zβ∂Rk+ℓf. The old estimate bounds each by a constant times (R′)μ−k, since R and R′ are comparable. If R′≥1 corresponds to R<1, it lies in a bounded interval with a positive lower bound; smoothness on that compact interval gives the same inequality after enlarging the constant. This proves all orders of the coordinate invariance. We assert compact local symbol bounds, with no uniform bound over the whole time axis.
For f∈Sμ(M), put
C(q,t)=∫0tc(q,s)ds,u(q,t,R)=e−C(q,t)∫0teC(q,s)f(q,s,R)ds.(FG9)
For the complex integrating factor, use the actual definition
eC=eReC(cos(ImC)+isin(ImC)).(FGA1)
The earlier real exponential and trigonometric derivative identities give ∂eC=eC∂C for every real parameter derivative. Their product identities give eCe−C=1. Repeated differentiation is therefore a finite sum of products of derivatives of C and this same exponential; on each compact interpolation set all these factors are bounded. This justifies the complex factor without imposing a reality assumption on c.
Then u∈Sμ(M), u∣t=0=0, and (v+c)u=f. The last equation follows by the product rule and the fundamental theorem, with the displayed signs. For the estimates, a compact set E⊂Ω has a compact interpolation image {(q,θt):(q,t)∈E,0≤θ≤1}⊂Ω. On it all derivatives of c,C,e±C are bounded. Each differentiated integral is a finite sum of endpoint derivatives and integrals of derivatives, of bounded lengths. Each radial derivative falls only on f, and gives exactly Rμ−k. This proves every estimate in (FG8). If f is exactly homogeneous of degree μ, so is u.
The relation manifold in these coordinates consists of (q,t,s,R) with t,s∈I(q). The field (v,0) is ∂t. On this full conic relation its trajectory labels are (q,s,R). After quotienting by positive dilation, the radial coordinate disappears and the trajectory labels are (q,s). In particular its projected flow domain is
Ω={(q,s,t):(q,s)∈Ω, (q,t)∈Ω}.(FGA2)
This is open over the manifold of labels (q,s)∈Ω. Its t-fiber is the same interval I(q), which contains zero. It thus has exactly condition (C) of Theorem 1.1, and the full conic relation has the radial product above. Applying the already proved construction gives the relation-manifold straightening and transport assertion, with the radial variable retained in the full trajectory quotient.
6. The characteristic relation is closed and Lagrangian
Return to a proper scalar operator with real homogeneous principal symbol p of any real degree m. Assume global real principal type and the compact return condition of the preceding global solvability theorem. Following the traditional terminology, these two conditions say that the domain is pseudoconvex with respect to the operator. They concern characteristic geometry, rather than convexity of its base coordinates.
Let N={p=0}⊂T∗X∖0. If it is empty, the characteristic relation is empty and the assertions below are vacuous. Otherwise dp is nonzero at each point of N: a zero would make Hp vanish there and give a complete stationary characteristic in a compact base set, contrary to global real principal type. The written inverse and implicit theorem therefore makes N a smooth hypersurface. Choose a positive degree-one norm w. For the radial Euler field E=∑jξj∂ξj, homogeneity gives dp(E)=mp=0 on N, whereas dw(E)=w>0. Thus the differentials of p and w are independent on N∩{w=1}. The same implicit theorem makes that section a smooth manifold of dimension 2n−2, and scaling identifies it with the quotient Ns.
Set p=w1−mp. On N, its Hamilton field is w1−mHp, so it gives the same unparametrized trajectories. The full normalization, radial exclusion and homogeneous lifting argument in compact-set solvability, Section 1, applies. The projected field on Ns has no compact complete trajectory. A compact endpoint set in Ns projects into a compact base set; the assumed return bound confines all intervening base points to a compact set, whose cosphere bundle is compact. Thus that projected field satisfies (A).
If n=1, a nonempty characteristic cosphere is zero dimensional; its projected field would vanish and violate escape. Thus the nonempty case has n≥2. It follows that N has (FG7) coordinates, and its same-trajectory relation C has coordinates (q,t,s,R) as in Section 5. It is a closed embedded conic submanifold of N×N, of dimension 2n. Indeed dimN=2n−1, dimNs=2n−2 and dimM0=2n−3; adding two times and one radius gives 2n.
We prove closedness also in the larger T∗(X×X)∖0, where one covector could a priori tend to zero. Suppose relation pairs have a limit there. Their base endpoints lie in a compact set K; pseudoconvexity confines every connecting interval to a compact K′. On the characteristic unit cosphere over K′, the positive invariant radius satisfies
0<c≤R/w≤C<∞.(FG10)
It has the same value at the two endpoints. Consequently their cotangent norms have ratios between c/C and C/c. One cannot tend to zero while the other tends to a nonzero covector. Both limit covectors are therefore nonzero, lie in N, and closedness in N×N puts the limit in C. This proves the claimed ambient closedness with the zero-section issue included.
Use ω=∑jdξj∧dxj and the difference symplectic form ω⊕(−ω) on the two cotangent factors. For the convention Hp=∑j(∂ξjp∂xj−∂xjp∂ξj), the contraction is ιHpω=−dp. Here is a coordinate proof that Hamilton flow preserves ω. Order coordinates as (x,ξ), put J=(0−II0), and write S for the symmetric Hessian of p. The variational equation is Y˙=JSY; the matrix of ω is −J. Consequently the derivative of YT(−J)Z is zero, since (JS)T(−J)+(−J)JS=−S+S=0. This proves preservation in each cotangent chart. A finite flow interval is covered by finitely many such charts; composing their local flow maps proves preservation on the entire interval.
Parameterize a relation neighborhood by (y,τ)↦(ϕτ(y),y), y∈N. For two variations (Y,a),(Z,b), the equal-time terms cancel by symplectic invariance. The remaining mixed terms are −adp(DϕτZ)+bdp(DϕτY), zero because both transported variations are tangent to N. The two-time term is zero by skew symmetry. Thus the difference form vanishes on C. Its dimension is half the ambient dimension, so C is Lagrangian.
The kernel convention identifies (x,ξ;y,η) with (x,y;ξ,−η)∈T∗(X×X). It pulls the product cotangent symplectic form back to ω⊕(−ω), so the corresponding kernel relation is also Lagrangian, with the stated sign on the input covector. We have proved that the bicharacteristic relation is a homogeneous canonical relation, closed even where the product cotangent space allows only one nonzero component. Replacing p by a real smooth nonvanishing multiple changes trajectory parameters but preserves this relation; uniqueness on each characteristic curve proves that assertion, including a reversal by a negative multiple. On any fixed compact segment a real nonvanishing multiplier and its reciprocal are bounded, so integration of its reciprocal gives the local change of clock throughout that segment. Exhausting a maximal trajectory by its compact segments proves equality of the whole unparametrized trajectories, even if the new maximal time interval has different endpoints. These geometric results do not yet construct its Fourier-integral parametrix.
The construction also respects a symplectic diffeomorphism κ defined on the full characteristic domain. If q=p∘κ−1, preservation of the symplectic form gives ικ∗Hpω=−dq. Nondegeneracy makes κ∗Hp=Hq. Uniqueness of the flow then carries the full same-trajectory relation by κ×κ, with the new maximal time domain. If κ commutes with positive dilation, this also preserves the homogeneous structure. When a map is defined only on a smaller open set, the assertion applies to trajectory intervals contained in that set.
An exact model of the invariant radius

Take the conic manifold M=(−1,1)×R+, with coordinates (t,r), dilation (t,r)↦(t,λr), and v=∂t+(log2)r∂r. The projected field is ∂t on the open interval; its trajectory space is a single point. Formula (FG6) gives R=r2−t, so its curves are r(t)=R2t. The upper panel draws the three actual curves with R=1/2,1,2. The lower panel draws their images, where R is constant and the field is ∂t. The section t=0 meets every curve exactly once. Arrowheads specify increasing time, not the magnitude of a velocity vector.
The open circles at t=−1,1 show excluded endpoints; they are outside M. Their radii in the upper panel are exactly R/2 and 2R. A maximal curve approaches these missing boundary points, so its image is not relatively compact in M. The dotted lines mark those boundaries and the zero-time section; the drawing only truncates the otherwise positive radial axis. This is a coordinate example of (FG6)–(FG7), not an identification of these two coordinates with a cotangent chart. The real power and derivative identities used here are the earlier complete scalar proofs.
7. Exercises with complete solutions
1. A quotient that fails to be Hausdorff. In the slit plane X=R2∖{(0,y):y≥0}, take v=∂x. Explain why no complete trajectory is relatively compact, but the trajectory quotient is not Hausdorff. Identify the missing hypothesis in Theorem 1.1.
Solution. Below height zero, a trajectory is the full horizontal line. At nonnegative height it is one of the two open half-lines. Each maximal trajectory has an unbounded base direction, so none is relatively compact. Let the two quotient points be the negative and positive half-lines at height zero. Any open saturated neighborhoods of these trajectories contain points at sufficiently small negative heights near (−1,0) and (1,0), respectively. For a common sufficiently small negative height these two points lie on the same full horizontal trajectory. The quotient neighborhoods therefore intersect, and the two classes cannot be separated. The compact return condition fails: the intervals from (−1,−1/j) to (1,−1/j) have endpoints in the two compact vertical segments, but their midpoints approach the excluded origin and escape every compact subset of X. This is exactly the missing second clause of (A), rather than a periodic orbit or a zero of v.
2. Radial drift and resonant scalar transport. On R×R+, let v=∂t+αr∂r, with real constant α, and c=β∈C. Find the invariant radius and solve (v+β)u=rμ, u(0,r)=0, for real μ, including μα+β=0.
Solution. The invariant radius is R=re−αt. In its coordinates, the forcing is Rμeμαt, and (FG9) gives
u(t,r)={rμ(1−e−(μα+β)t)/(μα+β),trμ,μα+β=0,μα+β=0.(FG11)
For a function rμh(t), the left side is rμ[h′(t)+(μα+β)h(t)]. Both displayed branches satisfy this scalar equation with right side one and zero initial value. The solution is homogeneous of degree μ; every radial derivative has order μ−k on compact time sets. Exponential growth on an unbounded time axis does not contradict these local symbol bounds. The resonant branch also follows from the limit (1−e−δt)/δ→t.
3. The free transport canonical relation. For n≥2, use x=(t,z) on Rn and P=Dt. Write its characteristic relation, check its dimension and symplectic form, and verify ambient closedness.
Solution. The characteristic covectors are (0,ζ), with ζ=0, and trajectories change only t. The relation is
C={((t,z;0,ζ),(s,z;0,ζ)):t,s∈R, z∈Rn−1, ζ=0}.(FG12)
Its coordinates have dimension 2+(n−1)+(n−1)=2n. The two transverse symplectic terms cancel because both factors have the same z,ζ, and the longitudinal terms vanish since both longitudinal covectors are zero. Thus it is Lagrangian for the difference form. Its endpoint covector norms are equal; in a product cotangent limit with at least one nonzero component, both retain the same nonzero ζ. The relation remains closed. For compact endpoints, the transverse coordinates stay in a compact projection and the longitudinal coordinate lies between two bounded endpoint values, giving the compact return bound. When n=1, the characteristic set for Dt is empty, so the displayed nonzero transverse-covector model has no points.
References and current scope
Hörmander IV, Definition 26.1.10, Lemmas 26.1.11–26.1.12 and their two remarks, and Proposition 26.1.13, printed 67–69/PDF 78–80, supply the mathematical antecedents. The quotient dimension, proper embedding, smooth affine section, all symbol derivatives and ambient zero-section exclusion are developed here in full, using the exact written programme foundations. The exposition and exercises are independently written; no source error or new theorem is claimed. The cited AN-03 proofs are licensed under GFDL 1.2 only, with no invariant sections, front-cover texts or back-cover texts. Their expression and assets are not reproduced here. Original expression is eligible for the course's CC0 dedication.
The next analytic step is to construct a Fourier-integral parametrix on this relation and solve its successive symbol transport equations. The geometric and scalar transport results proved here supply the global coordinates for that construction.
The new exact-model drawing and its reproducible source are original programme work under CC0. Its embedded font retains the DejaVu licence.