Energy for strict complex oblique boundary conditions
The boundary operator can contain an imaginary tangential derivative. A real multiplier alone does not give the required boundary sign. We incorporate that derivative in a Hermitian divergence form of the actual wave operator, then choose a future multiplier with a small outward normal component.
This lesson proves a weighted estimate for smooth functions, with the volume energy, terminal energy, full boundary gradient and inhomogeneous boundary data all retained. It includes variable principal coefficients and arbitrary smooth complex lower terms. The Dirichlet lesson provides the separate existence theorem for Dirichlet data.
1. Geometry, source conventions and the theorem
H0. Fix the actual operators. Let X be a smooth manifold with boundary, g a tangent Lorentz metric of signature (1,n−1), and G=g−1. Fix a smooth positive density μ. In local coordinates the scalar operator is
L=Gαβ∂α∂β+dβ∂β+d0.(OB1)
The coefficients Gαβ are real and symmetric; dβ,d0 may be complex. A real clock ϕ is proper and G(dϕ,dϕ)>0. Its positive dual fixes the future cone. The boundary is timelike: its nonzero conormals have negative G-square. Thus its tangent metric has one positive direction.
On the boundary let N be the inward unit spacelike normal, and V,W real tangent vector fields. The normalized first-order operator is
B=N−V−iW+b,b∈C∞(∂X;C).(OB2)
Here vector fields act by ordinary derivatives. Write
f2e2=g(V,V)>0,c=g(W,V)/f,W⊥=W−(c/f)V,=−g(W⊥,W⊥),δ=1−e2>0,(OB3)
and assume that V is future timelike. Nonnegativity of e2, the smooth Lorentz frames used below, and the positive time cone follow from the Lorentz algebra in DC:T002, DC:T008–T010. The notation e is its nonnegative square root; the proof never requires a smooth choice of W⊥/e at a zero.
Choose a compact time interval [a,T] inside a slightly larger interval on which these hypotheses hold. Properness gives a compact slab. All norms below use fixed smooth positive volume and hypersurface densities and a fixed positive definite norm on covectors. Changing these choices on this slab changes constants only. Let
Ω={a<ϕ<T},ΣT={ϕ=T},Y=∂X∩{a<ϕ<T}.
The temporal levels meet the boundary transversely by DC:T003. We take u smooth up to the boundary and zero on {ϕ≤a}. A local coordinate version allows compact spatial support instead of properness.
S0. Weighted energy theorem. There are C,λ0>0, depending on this fixed slab and its coefficients but not on u or λ, such that for λ≥λ0,
λ∫Ωe−λϕ(∣du∣2+∣u∣2)dμ+e−λT∫ΣT(∣du∣2+∣u∣2)dST+∫Ye−λϕ(∣du∣2+∣u∣2)dSY≤λC∫Ωe−λϕ∣Lu∣2dμ+C∫Ye−λϕ∣Bu∣2dSY.(OB4)
The derivative on Y is the full gradient of u, including its normal derivative. Sections 2–6 prove this statement. The constants need not be uniform over unbounded time intervals.
A0. Match the source's strict inequalities. In Hörmander II, Example 12.9.14(b), write the tangential coefficient vector as F+iE=(b2,…,bn−1,−bn), with the last entry temporal. The Lorentz pairing is positive on that last direction. In our notation these are the metric-dual covectors of V,W. Consequently
[F,F]=f2,[E,F]=fc,[E,E]=c2−e2,[F,F](1+[E,E])−[E,F]2=f2δ.(OB5)
Thus the strict conditions are precisely that F is future timelike and δ>0; then 1+[E,E]=c2+δ>0 too. To check signs, with N=∂r, D=−i∂, multiplying B(D)=Dr+∑bjDyj+bnDt and its datum by i gives (OB2): Vt=−Rebn, Vyj=−Rebj, and the same formulas with imaginary parts give W. In the real-principal convention of the Dirichlet lesson one also uses L=−P; every forcing is multiplied by the same scalar. Formula (OB5) establishes the exact strict part of the source condition used here.
A1. Change no second-order operator. Extend N,V,W smoothly into a boundary collar for the moment. Set
A=N∧W,M=G+iA,(N∧W)αβ=NαWβ−WαNβ.(OB6)
The matrix M is Hermitian. Since Aαβuαβ=0, the actual operator (OB1) is exactly
L=divμ(Mdu)+Cβuβ+d0u,Cβ=dβ−μ−1∂α(μMαβ).(OB7)
All derivatives of A, the density and the principal coefficients are retained in C.
Let Z be a real vector field and m>0 a constant. With QM(du)=Mβγuγuβ∈R, define
Jα=2Re(MαβuβZu)−ZαQM(du)+mZα∣u∣2.(OB8)
For every real U, adding the bivector Z∧U to A leaves J unchanged. Indeed, put z=(Uu)Zu. The first term changes by −2ZαImz, while the negative quadratic term changes by +2ZαImz. They cancel identically. We may therefore discard these “electric” bivectors when testing positivity.
A2. Retain the complete divergence remainder. Differentiating (OB8), the second-derivative cross terms cancel by Hermitian symmetry and commutation of scalar coordinate derivatives. The result is
divμJ=2Re(LuZu)+Rm(u),∣Rm(u)∣≤Cm(∣du∣2+∣u∣2).(OB9)
Here is the full local remainder, which also verifies that no unestimated second derivative occurs:
Rm=−2Re((Cβuβ+d0u)Zu)+2Re(Mαβuβ(∂αZγ)uγ)−(divμZ)QM−Zα(∂αMβγ)uγuβ+m(divμZ)∣u∣2+2mRe(Zuu).(OB10)
For example the two apparent Hessian terms are
2Re(MαβuβZγuγα) and
−2Re(ZαMβγuγαuβ).
Conjugate the first, use Mαβ=Mβα, and rename the dummy indices to obtain their cancellation. Every remaining factor in (OB10) has at most one derivative of u, and bounded coefficients on the fixed compact slab give (OB9). These formulas are computed in charts; their left side and the current are intrinsic, so the local bounds combine over a finite cover.
3. Positive clock and causal flux
F0. A multiplier close to the tangent future field. On the boundary take
Z=V−κN,0<κ≤min{2fδ,f2+1+c2f2δ}.(OB11)
A single positive κ satisfies these bounds on a compact boundary slab. This multiplier points outward in its normal component and remains future timelike.
The component (c/f)N∧V of A is electric relative to Z. Put
NZ=N−f2−κ2κZ.
Direct scalar products give
g(NZ,Z)=0,g(NZ,NZ)=−f2−κ2f2,g(NZ,W⊥)=g(Z,W⊥)=0.(OB12)
Modulo electric bivectors, A=NZ∧W⊥. In a Lorentz frame with unit temporal vector Z/g(Z,Z), its spatial skew block has norm
eZ=f2−κ2fe<1.(OB13)
Indeed κ2/f2≤δ/4 and e2=1−δ. Choosing an orthonormal basis of the plane of NZ,W⊥ gives a single 2-by-2 skew block; if W⊥=0, that block is zero. This is pointwise linear algebra and requires no division by a vanishing field.
In this frame let H=I−iAsp. Its two possibly nonunit eigenvalues are 1−eZ,1+eZ: on the skew plane the eigenvectors are (1,i) and (1,−i), with their order depending on the orientation. All other eigenvalues are 1. Thus H is positive Hermitian, uniformly so on the compact slab.
F1. Keep the physical characteristic cone. Let ℓ=(ℓ0,ℓsp) be a real future conormal, so ℓ0≥∣ℓsp∣. For the derivative part of the current,
g(Z,Z)ℓαJderivα=ℓ0(∣u0∣2+usp∗Husp)−2Re(u0ℓsp∗Husp)≥(ℓ0−∣ℓsp∣)(∣u0∣2+usp∗Husp).(OB14)
For the inequality, ℓsp is real, so
ℓsp∗Hℓsp=∣ℓsp∣2.
The Hermitian Cauchy–Schwarz inequality follows by minimizing
(x−zy)∗H(x−zy)≥0 over z∈C; it bounds the cross term by
2∣u0∣∣ℓsp∣(usp∗Husp)1/2.
Now use 2ab≤a2+b2.
The mass term is mℓ(Z)∣u∣2, positive for every nonzero future causal ℓ. In particular dϕ⋅J and the future terminal flux control ∣du∣2+m∣u∣2 uniformly. Future null faces have nonnegative flux. All cones here are the original metric cones.
4. The full physical-boundary flux
B0. Compute the exact boundary matrix. The outward unit conormal is ν=g(N,⋅); it evaluates to −1 on the inward N. Its natural conormal derivative is
ναMαβuβ=Nu−iWu=Vu+Bu−bu.(OB15)
This follows directly from ναGαβ=Nβ,
ν(N)=−1 and ν(W)=0, including in coordinates with mixed coefficients.
At one point choose a Lorentz frame N=∂r, V=f∂t, W=c∂t−e∂y1. On the kernel of the principal part B0=N−V−iW,
ur=(f+ic)ut−ieuy1.(OB16)
Substitution in Jderiv⋅ν gives the form in (uy1,ut)
Kκ=(κ(1−e2)κe(c+if)κe(c−if)2f2−κ(f2+1+c2)),detKκ=2κf2δ−κ2(f2+1+c2−e2).(OB17)
Every other tangential spatial derivative contributes κ∣uyj∣2.
For an explicit verification of the signs, in the ordered frame (t,r,y1),
M=1ic0−ic−1ie0−ie−1,Zu=fut−κur,Jderiv⋅ν=−2Re((Mdu)rZu)−κQM.(OB18)
Use (OB16), expand, and collect the coefficients of ∣uy1∣2, ∣ut∣2 and uy1ut; they are precisely (OB17).
The second bound in (OB11) gives
(Kκ)tt≥f2,detKκ≥κf2δ>0.(OB19)
For the determinant, replace f2+1+c2−e2 by the larger f2+1+c2. The first diagonal entry is positive, and completing the complex square proves positivity of the entire 2-by-2 form. Relation (OB16) also bounds the normal derivative by the tangential ones. Compactness therefore gives uniform control of the full gradient on this principal boundary kernel.
When e=0, choose any spatial frame and retain the same formula with zero cross term. If there is no tangential spatial direction, namely ambient dimension two, only the positive tt entry is needed. This handles that case directly.
B1. Keep the boundary datum and zeroth term. For an arbitrary gradient put h=B0u. Subtract the vector with normal component h and all tangential components zero. The remainder lies in the kernel just considered. Expanding the quadratic form and using 2ab≤ϵa2+ϵ−1b2 gives
Jderiv⋅ν≥c1∣du∣2−C1∣h∣2.
All coefficients are uniformly bounded, and h=Bu−bu; thus, for constants independent of m,
J⋅ν≥c1∣du∣2−C2∣Bu∣2+(mκ−Cb)∣u∣2.(OB20)
The splitting argument also controls the normal component of the arbitrary gradient, since its modulus is bounded by that of the kernel component plus ∣h∣. Choose m once so that mκ−Cb≥1. No derivative of the boundary datum occurs in (OB20).
5. Construct one current on a compact slab
G0. Extend in two separate stages. The boundary fields and the collar needed above can be constructed on a neighborhood of the compact boundary slab. Choose finitely many boundary charts with inward defining functions rj, and bumps whose positive sets cover this compact set. Divide the bumps by their positive sum to form a partition. The function r=∑χjrj has, on the boundary, a positive weighted sum of inward conormals; all these lie on the same positive ray. Thus dr=0 and G(dr,dr)<0 there. Both properties persist in a smaller neighborhood. The field
N=−G♯dr/−G(dr,dr) is smooth and inward unit normal on each level of r.
The smooth flow and compact continuation proved in DC:T005–T006, followed by the inverse theorem, give a common short collar. Extend V,W smoothly along it and project them to N⊥. Their boundary values are unchanged. The strict margins persist after shrinking. Choose the common κ>0 smaller if necessary. Equations (OB12)–(OB14) then hold throughout this collar.
Keep Z=V−κN fixed while multiplying A=N∧W by a scalar cutoff 0≤χ≤1, equal to one near the boundary and zero before the collar ends. Its spatial skew norm relative to Z is χeZ<1; equivalently its current is the convex combination of the currents with A and with A=0. Positivity persists.
Only where A=0 blend Z with the global future field G♯dϕ. The future timelike cone is convex: in a Lorentz frame its vectors satisfy v0>∣vsp∣, and the triangle inequality preserves that strict inequality under positive convex combinations. Hence the blended field stays future timelike. Extend it as G♯dϕ in the remaining interior. Finite charts give smooth cutoffs and uniform margins. This constructs a single smooth A,Z on a neighborhood of the compact slab with the precise prescribed boundary values.
The construction uses no partition of u or its equation. It therefore introduces no boundary commutator. If the boundary is empty, use A=0 and Z=G♯dϕ throughout. This also proves the interior version in dimension one.
6. Integrate the actual identity
E0. Prove the weighted estimate. The exact identity is
λe−λϕdϕ⋅J+divμ(e−λϕJ)=e−λϕ{2Re(LuZu)+Rm}.(OB21)
Use the divergence formula DC:D1 over Ω. It follows from iterated integration on boundary graphs and extends to the present transverse corners by a finite partition and shrinking strips; no normal derivative of a rough function is being assumed here. The past flux is zero because u is smooth and zero for ϕ≤a.
On the physical boundary the divergence density equals the flux in (OB20) times a smooth positive factor: in a defining chart ν is a negative multiple of dr, so this is the ordinary outward divergence sign. The same observation on ΣT, with a positive multiple of dϕ, gives its future sign. All density factors have positive upper and lower bounds on the compact faces.
Equations (OB14) and (OB20), with m now fixed, put positive multiples of the three left-hand terms of (OB4) on the left, except for C2∥Bu∥L2(Y,e−λϕ)2, which is moved to the right. The remainder satisfies (OB9). Boundedness of Z and Young's inequality give
2∣Lu∣∣Zu∣≤4c0λ∣du∣2+λC∣Lu∣2,(OB22)
where c0>0 is a clock coercivity constant. Choose λ0 so that Cm(∣du∣2+∣u∣2) consumes at most another quarter of the positive volume term. Absorb these terms. Norm equivalence on the fixed slab gives exactly (OB4).
Additional artificial faces with future causal outward conormal contribute nonnegative flux by (OB14) and can be discarded. Compactly supported local functions give the half-space version by the same proof. Spatial cutoff limits are available when every term converges, as in the Dirichlet lesson. This completes the smooth weighted theorem. □
7. Exact models and a reproducible diagram
X0. A strict numerical model with exact constants. Take
L=∂t2−∂r2−∂y2,B=∂r−∂t+2i∂y,f=1, e=21, c=0, κ=41.
Then
K=(3/16i/8−i/83/2),detK=17/64,trK=27/16,K≥10817I.(OB23)
The last bound follows because the smaller positive eigenvalue is the determinant divided by the larger one, and the latter is no greater than the trace. The quadratic formula gives those two real eigenvalues, or completing the square verifies the displayed bound directly. On the principal boundary kernel,
∣ur∣2=∣ut−2iuy∣2≤2∣ut∣2+21∣uy∣2.
Thus ∣du∣2≤3(∣ut∣2+∣uy∣2), and the boundary derivative flux is at least 17∣du∣2/324. The shifted spatial block in the unit-Z frame has eZ=2/15, not 1/2.
V0. Read the planes separately. The first panel is physical (r,t)-space with r≥0, g=dt2−dr2, N=(1,0), V=(0,1) and Z=(−1/4,1). Its Z arrow has negative normal component and g(Z,Z)=15/16. Null directions are drawn at slope one. The second panel is a real two-dimensional slice of the complex derivative space of the spatial Hermitian block in the unit-Z frame. If a is the NZ derivative and the W⊥ derivative is ib, its form is a2+b2+2(2/15)ab; the two principal values are 1±2/15. It is not a spatial section of X. The third panel plots absolute values of the exact modes below at k=1, before cutoffs; the neutral and growing time factors are distinguished. All curves, arrows and constants are generated from the formulas in the editable figure script.
X1. Exact boundary-layer modes. On a periodic y-coordinate set
Be=∂r−∂t+ie∂y,wk=ek(σt−αr)e−iky,σ=2ee2−1,α=2ee2+1,e>0.(OB24)
The identities α2−σ2=1 and α+σ=e show, by direct differentiation, that Lwk=Bewk=0. An inactive periodic coordinate places this example in the source's ambient dimension four. A finite normal interval, with the function zero near its far endpoint, makes the clock proper.
X2. Why strictness matters. Fix T>0, a smooth θ zero for t≤0 and one for t≥T/2, and a normal cutoff ρ equal to one for 0≤r≤r0 and zero for r≥2r0, inside that finite interval.
At e=1, take uk=k−1/2θ(t)ρ(r)e−kr−iky. On a fixed later time interval, integrating ∣∂yuk∣2 for 0<r<1/k gives a positive lower bound independent of k. But
∥Luk∥L2(Ω)=O(k−1)+O(kAe−kr0),∥B1uk∥L2(Y)=O(k−1/2)(OB25)
for a fixed integer A. Indeed the time term is k−1/2θ′′ρe−kr−iky, whose normal integral contributes another factor k−1/2. The normal cutoff terms have support r≥r0 and are exponentially small. At r=0 only −k−1/2θ′e−iky remains in the boundary datum. For fixed λ, weights on the finite slab are equivalent to one. Therefore (OB4) cannot be extended to this nonstrict edge.
For e>1, put uk=e−σkTθρwk. The time-cutoff residual is bounded by a polynomial in k times e−σkT/2; the normal-cutoff residual is bounded by a polynomial times e−αkr0. The boundary datum is again a time-cutoff term and is exponentially small. Any fixed finite number of integer coordinate derivatives costs only another polynomial in k, so all these forcing and boundary norms tend to zero even after such derivatives.
On T−1/k<t<T, 0<r<1/k, however, ∣e−σkTwk∣≥e−σ−α. Its y-derivative has modulus at least ke−σ−α, while the (r,t) region has area k−2. The fixed periodic factors have positive measure. Thus ∥duk∥L2(Ω) has a positive lower bound. An estimate by any fixed finite collection of those differentiated forcing and boundary L2 norms also fails. Future cutoffs may be placed after T; they do not affect any norm just used. These are precise failures of these energy estimates, rather than a classification of every weaker boundary problem.
8. Exercises with complete solutions
Exercise 1. Normalize a complex boundary operator
Let B(D)=Dr−21Dy−Dt+i(41Dy+2Dt). Determine V,W,f,c,e,δ in the flat metric and decide whether the strict theorem applies.
Solution. Multiplication by i gives the ordinary operator ∂r−21∂y−∂t+i(41∂y+2∂t). Thus V=∂t+21∂y, W=−2∂t−41∂y. The metric gives
f2=3/4,g(W,V)=−15/8,c=−4315,g(W,W)=63/16.
Hence e2=c2−g(W,W)=75/16−63/16=3/4, and δ=1/4>0. Also V is future. Directly the source determinant is (3/4)(1+63/16)−(15/8)2=3/16=f2δ. The strict theorem applies on every compact slab. The datum must also be multiplied by i; no conjugation of the boundary operator is involved.
Exercise 2. Keep a large imaginary temporal coefficient
Take f=1,e=1/2,c=2. Give one permitted κ, compute its boundary matrix and a quantitative positive lower bound. Explain what changes if the boundary has no spatial direction.
Solution. Now δ=3/4, and (OB11) permits κ=1/64, since 1/64<3/4 and 1/64<1/8. Formula (OB17) gives
K=(3/256(2+i)/128(2−i)/12861/32),detK=16384361,trK=256491.(OB26)
The determinant/trace argument from (OB23) yields K≥361I/31424. The principal relation bounds
∣ur∣2=∣(1+2i)ut−2iuy∣2≤10∣ut∣2+21∣uy∣2,
so ∣du∣2≤11(∣ut∣2+∣uy∣2). A full-gradient flux bound is consequently 361∣du∣2/345664. If there is no spatial direction then W=c∂t, e=0, and the time entry is 2f2−κ(f2+1+c2)>0. There is no y block to discard or divide by.
Exercise 3. Check the unstable mode and its normalization
For e=3/2, compute σ,α. Verify both operator identities in (OB24), and prove that the normalized cutoff family has a positive gradient norm while every prescribed finite collection of differentiated residual norms tends to zero.
Solution. The values are σ=5/12,α=13/12. Therefore
σ2−α2+1=25/144−169/144+1=0, the multiplier of k2wk in Lwk; and −α−σ+e=−13/12−5/12+3/2=0, the multiplier of kwk in Bewk.
For uk=e−5kT/12θρwk, the product rule gives
Luk=e−5kT/12{(θ′′+2σkθ′)ρ−θ(ρ′′−2αkρ′)}wk.(OB27)
The first bracket is supported where t≤T/2; the second where r≥r0. Each derivative of total order at most a prescribed integer j is bounded by
Cjkj+2(e−5kT/24+e−13kr0/12) on the finite slab. Its L2 norm has the same bound up to its fixed volume. On the physical boundary only −e−5kT/12θ′wk remains, with the analogous exponential bound for each boundary derivative. The family is zero near the other normal endpoint. In the last 1/k units of time and first 1/k units of normal distance, θ=ρ=1 for large k, so the lower bound from Section 7 is e−3/2 for the amplitude and ke−3/2 for its y-derivative. Integrating over that region proves the positive lower bound. Every exponential upper bound beats any of the finitely many derivative powers. This proves the asserted finite-order failure with all cutoffs retained.
9. Sources and prerequisite proofs
The primary antecedents are Lars Hörmander, The Analysis of Linear Partial Differential Operators II, Springer 2005 eBook, Example 12.9.14(b), printed 178–179, and III, Springer 2007 eBook, Section 24.1, especially its final energy-extension paragraph on printed 423. The source's Lorentz form is positive in the temporal direction (II, Section 12.1). The argument above is independent exposition; the source example's full classification and its illustration are not reproduced.
The proof map identifies the exact earlier programme proofs. In particular DC:T002–T006 provide Lorentz and compact-slab geometry, DC:T008–T010 provide the real flux algebra and frames, DC:T005 gives the full smooth-flow construction, DC:D1 supplies the divergence formula including corners, and U001:P3 and U001:P14.3 provide the inverse theorem and smooth cutoffs. The external human foundations of those earlier proofs are explicitly recorded in the map. Their presence outside this export is retained in its dependency metadata.
The new exposition, proofs, exercises and original figure are by GPT-6.1 Sol (OpenAI). Original text and figure sources: public domain under CC0-1.0. Linked programme components retain their own authorship and terms.