We will fix a manifold N, which will be thought of as ∂M, and consider the manifold N×R≥0. We will write y∈N and u∈R≥0. We let E be a Hermitian vector bundle on N.
Let A be a first-order self-adjoint elliptic operator on E, and consider the differential operator
D=∂u∂+A.
This is the differential operator of interest on N×R≥0. While the definition resembles that of a heat operator, it may not be the best idea to think of it as a heat operator. We will later introduce a further time variable t and consider e−tD∗D, and confusion might arise.
The basic problem we want to solve is the equation
Df=g.
For now, we assume g∈Γc(E×R≥0), and later extend to more general functions after we establish the right bounds.
As before, since A is self-adjoint and elliptic, we can find an eigenbasis {ψλ} of A. We write
f(y,u)=λ∑fλ(u)ψλ(y),g(y,t)=λ∑gλ(u)ψλ(y).
The we see that for the equation to hold, we must have
(dud+λ)fλ(u)=gλ(u).
We have previously decided that the boundary condition we want is fλ(0)=0 for λ≥0. There is another way we can think about this boundary condition. Observe that the solution to this differential equation is well defined up to adding Ce−λu. If λ<0, then there is at most one solution with reasonable growth as u→∞. However, if λ≥0, then we can add any multiples of e−λu and have a sensible solution (λ=0 is somewhat of an edge case). Thus, fixing fλ(0)=0 gives us a unique “sensible” solution.
For convenience, we write P for the composition
H1(E×R≥0)→rH0(E×{0})→H0(E×{0}),
where r is restriction and the last map is the spectral projection onto the subspace spanned by the eigenvectors with non-negative eigenvalues. Then our boundary condition is Pf=0. By an abuse of notation, we write 1−P for r−P, and the boundary condition for D∗ will then by (1−P)f=0. We shall write Γ(E×R≥0;P) for the space of smooth sections f such that Pf=0.
Imposing this boundary condition, we can write down explicit solutions for fλ:
We check that for large u, the function fλ(u) is either exponentially decreasing (if λ>0), constant (if λ=0) or identically zero (if λ<0). The constant term wrecks our hope that this takes values in Hs, but we can still hope it takes values in Hlocs.
Note also that for λ<0, the value of fλ(u) depends on the future values of gλ, which ties in with our previous point — at large u, all future values of gλ are zero, hence fλ(u)=0.
We let Q:Γc(E×R≥0)→Γ(E×R≥0;P) be the function that sends gλ to fλ above, which is clearly linear. The next proposition will, in particular, show that this map is well-defined, i.e. Qf is actually smooth:
Proposition6.1
Q extends to a continuous map Hs→Hlocs+1.
Proof
□
Recall that since A is elliptic, we have
∥f∥Hs+1≤C(∥f∥Hs+∥Af∥Hs).
So we can equivalently define the H1 norm by
∥f∥H12=∥f∥L22+∥∂uf∥L22+∥Af∥L22.
Using that the eigenspace decomposition is orthogonal in L2 and the defining equations
Aψλ=λψλ,∂u∂fλ=gλ−λfλ,
we get an inequality
∥f∥H12≤C⋅λ∑(1+λ2)∥fλ∥L22+∥gλ∥L22.
To bound ∥fλ∥L2 in terms of ∥gλ∥L2, we use the (rotated) Laplace transform
f^λ(ξ)=∫0∞e−iuξgλ(u)du.
One checks easily that up to a constant which we omit, we have Parseval's identity:
∥f^λ∥L2=∥fλ∥L2.
We can explicitly compute
f^λ(ξ)=λ+iξg^λ(ξ)+fλ(0).
Using Parseval's identity, and the fact that
fλ(0)=−∫0∞eλugλ(u)du if λ<0,
we obtain bounds
∣λ∣∥fλ∥L2≤2∥gλ∥L2.
So if we ignore the λ=0 term, Q would send H0 into H1. The f0 term is eventually constant, given by the integral of the compactly supported function ∫0∞g0(v)dv. So Q maps H0 into Hloc1.
Using the equation
dusds(dud+λ)fλ=dusdsgλ
and calculating as above, we see that Hs gets mapped into Hlocs+1.
Proof
□
The main takeaway from the section is then the following theorem:
Theorem6.2
There is a linear operator
Q:Γc(E×R≥0)→Γ(E×R≥0;P)
such that
DQgQDf=g for all g∈Γc(E×R≥0),=f for all f∈Γc(E×R≥0;P).
Moreover, Q extends to a continuous map Hs→Hlocs+1 for all integral s. It is given by convolution with a kernel Qu(y,z), where u is now allowed to take negative values. The kernel is smooth away from u=0.