ScalingStacks

3.2. Green’s currents for Riemannian submanifolds [04ZJ]

Original official author HTML, exact retained edition. Historical TeX conversion verdicts remain unchanged. Cited-edition alignment and mathematical self-containment are not assessed.

Complete original source context · Original author HTML

3.2. Green’s currents for Riemannian submanifolds

First we recall and introduce the basic terminology. Let (Q,g)(Q,g) be an oriented Riemannian manifold of dimension mm. Denote by Ω0l​(Q)\Omega^{l}_{0}(Q) the space of differential ll-forms with compact supports in QQ.

Definition 3.5 (kk-current).

A kk-current on QQ is a linear functional T:Ω0m−k​(Q)→ℝT:\Omega^{m-k}_{0}(Q)\to\mathbb{R} which is continuous in the sense of distributions, i.e. suppose χj∈Ω0m−k​(Q)\chi_{j}\in\Omega_{0}^{m-k}(Q) is a sequence of differential forms with all derivatives uniformly converging to 00 as j→∞j\to\infty, then limj→∞(T,χj)=0\lim\limits_{j\to\infty}(T,\chi_{j})=0.

The notion of currents unifies the notion of differential forms and submanifolds. In particular, a locally integrable kk-form β\beta can be naturally viewed as a kk-current via the pairing

(3.36) (β,χ)≡∫Qβ∧χ,χ∈Ω0m−k​(Q),(\beta,\chi)\equiv\int_{Q}\beta\wedge\chi,\ \ \ \ \chi\in\Omega^{m-k}_{0}(Q),

and an oriented submanifold PP of co-dimension kk also defines a kk-current via

(3.37) (δP,χ)≡∫Pχ,χ∈Ω0m−k​(Q).(\delta_{P},\chi)\equiv\int_{P}\chi,\ \ \ \ \chi\in\Omega^{m-k}_{0}(Q).

The usual exterior differential dd and the Hodge star operator ∗* on differential forms then naturally extend to currents. Given a kk-current TT and χ∈Ω0m−k​(Q)\chi\in\Omega_{0}^{m-k}(Q), then we define

(3.38) (d​T,χ)≡(−1)k+1​(T,d​χ),\displaystyle(dT,\chi)\equiv(-1)^{k+1}(T,d\chi),
(3.39) (∗T,χ)≡(−1)k⁡(m−k)(T,∗χ).\displaystyle(*T,\chi)\equiv(-1)^{k(m-k)}(T,*\chi).

Let d∗d^{*} be the codifferential operator and denote by Δ≡d​d∗+d∗​d\Delta\equiv dd^{*}+d^{*}d the Hodge Laplacian, then it follows that for every kk-current TT and χ∈Ω0m−k​(Q)\chi\in\Omega_{0}^{m-k}(Q),

(3.40) (d∗​T,χ)=(−1)k​(T,d∗​χ),\displaystyle(d^{*}T,\chi)=(-1)^{k}(T,d^{*}\chi),
(3.41) (Δ​T,χ)=(T,Δ​χ).\displaystyle(\Delta T,\chi)=(T,\Delta\chi).

A kk-current TT is called harmonic if Δ​T=0\Delta T=0. It follows from the standard elliptic regularity theory that a harmonic kk-current can be represented by a smooth harmonic kk-form.

Now let P⊂QP\subset Q be a (not necessarily closed) embedded oriented submanifold. Although the following discussion applies to more general setting, for our purpose in the following we will only consider the case when PP is of co-dimension 33 in QQ. The importance of the co-dimension 33 case in our setting is related to the fact that there is a Hopf fibration ℝ4→ℝ3\mathbb{R}^{4}\rightarrow\mathbb{R}^{3} which is a singular S1S^{1} fibration with a smooth total space and co-dimension discriminant locus on the base. The co-dimension 3 condition also appears in other geometric settings, for example, Hitchin’s theory of Gerbes [Hit01].

Definition 3.6 (Green’s current).

Suppose PP is of co-dimension 3 in QQ. A Green’s current GPG_{P} for PP in QQ is a locally integrable 33-form which solves the following current equation on QQ

(3.42) Δ​GP=2​π⋅δP.\Delta G_{P}=2\pi\cdot\delta_{P}.
Example 3.7.

The above normalization constant is chosen such that in the case Q≡ℝ3Q\equiv\mathbb{R}^{3} and P≡03∈ℝ3P\equiv 0^{3}\in\mathbb{R}^{3}, then

(3.43) GP=12​|y|​d​y1∧d​y2∧d​y3G_{P}=\frac{1}{2|y|}dy_{1}\wedge dy_{2}\wedge dy_{3}

solves the current equation Δ0​GP=2​π⋅δP\Delta_{0}G_{P}=2\pi\cdot\delta_{P} for the standard Hodge Laplacian Δ0\Delta_{0} on ℝ3\mathbb{R}^{3}.

In particular GPG_{P} is harmonic outside PP hence is smooth. Notice a Green’s current GPG_{P} for PP is not unique, but it is unique up to the addition of a harmonic 33-form, so the singular behavior near PP does not depend on the particular choice of GPG_{P}. Also it is clear that if Q′⊂QQ^{\prime}\subset Q is an open submanifold, then the restriction of GPG_{P} to Q′Q^{\prime} is a Green’s current for P′=P∩Q′P^{\prime}=P\cap Q^{\prime} in Q′Q^{\prime}, so that we can study the regularity problem locally. Our goal in this subsection is to understand the local existence and regularity of GPG_{P} via approximation by the standard model, which is the product space ℝ3×ℝn−3\mathbb{R}^{3}\times\mathbb{R}^{n-3}.

To begin with, we have the following simple regularity result for d⁡(GP)d(G_{P}).

Proposition 3.8.

Given a Green’s current GPG_{P}, its differential d⁡(GP)d(G_{P}) extends to a smooth 44-form across PP.

Proof.

This is a local result so we can work with the geodesic ball Br​(p)B_{r}(p) for any p∈Pp\in P such that Br​(p)¯⊂⊂Q\overline{B_{r}(p)}\subset\subset Q and Br​(p)¯∩P⊂⊂P\overline{B_{r}(p)}\cap P\subset\subset P. We will show that the 44-current d​GPdG_{P} is a harmonic in Br​(p)B_{r}(p) in the distributional sense. In fact, for any test form χ∈Ω0m−4​(Br​(p))\chi\in\Omega^{m-4}_{0}(B_{r}(p)) we have

(3.44) (Δ⁡(d⁡(GP)),χ)=(d​Δ​GP,χ)=(Δ​GP,𝑑χ)=(2​π​δP,𝑑χ)=2​π​∫P𝑑χ=2​π​∫∂Br​(p)∩Pχ=0.(\Delta(d(G_{P})),\chi)=(d\Delta G_{P},\chi)=(\Delta G_{P},d\chi)=(2\pi\delta_{P},d\chi)=2\pi\int_{P}d\chi=2\pi\int_{\partial B_{r}(p)\cap P}\chi=0.

Therefore, d⁡(GP)d(G_{P}) is a harmonic 44-current in Br​(p)B_{r}(p) and hence it is smooth in Br​(p)B_{r}(p). ∎

In the rest of this section, we will frequently use the following notation.

Notation 3.9.

Given k∈ℤ+k\in\mathbb{Z}_{+}, a capital Greek letter with index kk, Πp(k)\Pi_{p}^{(k)}, always denotes a general local pp-form with 0≤p≤30\leq p\leq 3, which is of the form

(3.45) Πp(k)≡Pi,j1,…,jp−1(k)​(x,y)​(d​yj1∧…∧d​yjp−1)∧d​xi+Qj1,…,jp(k)​(x,y)​d​yj1∧…∧d​yjp,\Pi_{p}^{(k)}\equiv P_{i,j_{1},\ldots,j_{p-1}}^{(k)}(x,y)(dy_{j_{1}}\wedge\ldots\wedge dy_{j_{p-1}})\wedge dx_{i}+Q_{j_{1},\ldots,j_{p}}^{(k)}(x,y)dy_{j_{1}}\wedge\ldots\wedge dy_{j_{p}},

where Pi,j1,…,jp−1(k)​(x,y)P_{i,j_{1},\ldots,j_{p-1}}^{(k)}(x,y) and Qj1,…,jp(k)​(x,y)Q_{j_{1},\ldots,j_{p}}^{(k)}(x,y) are homogeneous polynomial functions in yy of degree kk whoses coefficient functions are smooth on U⊂PU\subset P.

Remark 3.9.1.

Notice this expression depends on the choice of local coordinates, but under a change of coordinates, a pp-form Πp(k)\Pi_{p}^{(k)} will still have such an expression, modulo a term which is of order O~​(rk+1)\widetilde{O}(r^{k+1}).

Now we are ready to state the first main theorem of this section, which gives a local existence for Green’s current, and its leading singular behavior.

Theorem 3.10.

If P⊂QP\subset Q is an embedded submanifold of co-dimension 33, then for any p∈Pp\in P, there is a neighborhood 𝒰\mathcal{U} of pp in QQ, and a Green’s current GUG_{U} for U=𝒰∩PU=\mathcal{U}\cap P in 𝒰\mathcal{U} satisfying

(3.46) Δ​GU=2​π⋅δUin𝒰,\displaystyle\Delta G_{U}=2\pi\cdot\delta_{U}\ \ \ \ \text{in}\ \ \ \ \ \mathcal{U},

and the local expansion

GU=\displaystyle G_{U}= 12​r​(1−Hα​yα2)​d​y1∧d​y2∧d​y3+12​r​yβ​Ai​α​β​d​xi∧d​yα^−14​Ai​j​α​β​r⋅d​yα​β^∧d​xi∧d​xj\displaystyle\frac{1}{2r}(1-\frac{H^{\alpha}y_{\alpha}}{2})dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}y_{\beta}A_{i\alpha\beta}dx_{i}\wedge dy_{\widehat{\alpha}}-\frac{1}{4}A_{ij\alpha\beta}r\cdot dy_{\widehat{\alpha\beta}}\wedge dx_{i}\wedge dx_{j}
+\displaystyle+ 316​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​d​(r​yα)∧d​xi∧d​xj+r−3​Π3(4)+O′​(r2),\displaystyle\frac{3}{16}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})d(ry_{\alpha})\wedge dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{(4)}+O^{\prime}(r^{2}),

where H→=Hα​eα\overrightarrow{H}=H^{\alpha}e_{\alpha} is the mean curvature of P⊂QP\subset Q and the 33-form Π3(4)\Pi_{3}^{(4)} is defined in (3.45).

Here and in the following we use the notation that for α∈{1,2,3}\alpha\in\{1,2,3\}, α^=(α+1,α+2)\widehat{\alpha}=(\alpha+1,\alpha+2) (with the convention 3+1=13+1=1), d​yα​β≡d​yα∧d​yβdy_{\alpha\beta}\equiv dy_{\alpha}\wedge dy_{\beta}, and that

(3.48) d​yα​β^≡{d​yα+2,β=α+1−d​yα+1,β=α+20,β=α.\displaystyle dy_{\widehat{\alpha\beta}}\equiv\begin{cases}dy_{\alpha+2},&\beta=\alpha+1\\ -dy_{\alpha+1},&\beta=\alpha+2\\ 0,&\beta=\alpha.\end{cases}
Remark 3.10.1.

By the above discussion, if GPG_{P} is any Green’s current for PP in QQ, then locally near p∈Pp\in P, GPG_{P} will also have an expansion of the form (3.10).

Remark 3.10.2.

At a given point p∈Up\in U, we can always choose a special frame {eα}\{e_{\alpha}\} such that Ai​α​β=0A_{i\alpha\beta}=0 at pp. On the other hand, since the singular behavior of GUG_{U} does not depend on the choice of coordinates, one sees that in general we need the second term in the expansion.

Remark 3.10.3.

By Proposition 3.8, d​GUdG_{U} is smooth. This is compatible with the above expansion. For example, from the expansion we see the leading term involving forms of type d​xi∧d​xj∧d​yβ∧d​yγdx_{i}\wedge dx_{j}\wedge dy_{\beta}\wedge dy_{\gamma} is given by

(3.49) (12​r​yβ​Ai​j​α​β​d​yα^−14​Ai​j​α​β​d​r∧d​yα​β^)∧d​xi∧d​xj.(\frac{1}{2r}y_{\beta}A_{ij\alpha\beta}dy_{\hat{\alpha}}-\frac{1}{4}A_{ij\alpha\beta}dr\wedge dy_{\widehat{\alpha\beta}})\wedge dx_{i}\wedge dx_{j}.

Elementary calculation shows this vanishes.

Remark 3.10.4.

In the proof we shall not keep track of the explicit form of Π3(4)\Pi_{3}^{(4)} because it is not needed in our applications. However, it is possible to obtain the precise expression with more work. Given the above expansion, there are also some constraint for Π3(4)\Pi_{3}^{(4)} following from the fact that d⁡(GU)d(G_{U}) is smooth by Proposition 3.8.

Before starting the proof of Theorem 3.10, we need some preparations. For the convenience of our calculations, we introduce three 11-forms

(3.50) ηα≡d​yα+pi​α​d​xi\eta_{\alpha}\equiv dy_{\alpha}+p_{i\alpha}dx_{i}

such that

(3.51) ⟨ηα,d​xj⟩=0\langle\eta_{\alpha},dx_{j}\rangle=0

for all jj and α\alpha at all points of 𝒰\mathcal{U}. Then the linear span of the ηα\eta_{\alpha}’s is orthogonal to the linear span of the d​xidx_{i}’s.

The lemma below is s crucial in the proof of Theorem 3.10.

Lemma 3.11.

For any 1≤k≤m−31\leq k\leq m-3 and 1≤α≤31\leq\alpha\leq 3,

(3.52) pk​α=gk​α+O~​(r3).p_{k\alpha}=g_{k\alpha}+\widetilde{O}(r^{3}).
Proof.

We write the full matrix expression of the metric gg as

(3.53) g=[gα​β00gi​j]+[0SSt0],g=\left[{\begin{array}[]{cc}g_{\alpha\beta}&0\\ 0&g_{ij}\\ \end{array}}\right]+\left[{\begin{array}[]{cc}0&S\\ S^{t}&0\\ \end{array}}\right],

where S=(gα​i)=O~​(r)S=(g_{\alpha i})=\widetilde{O}(r). We denote by (hα​β)(h_{\alpha\beta}) and (hi​j)(h_{ij}) the inverse matrix of (gα​β)(g_{\alpha\beta}) and (gi​j)(g_{ij}) respectively. Then by elementary consideration

(3.74) g−1\displaystyle g^{-1} =\displaystyle= [hα​β00hi​j]−[hα​β00hi​j]​[0SSt0]​[hα​β00hi​j]\displaystyle\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]-\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]\left[{\begin{array}[]{cc}0&S\\ S^{t}&0\\ \end{array}}\right]\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]
+\displaystyle+ [hα​β00hi​j]​[0SSt0]​[hα​β00hi​j]​[0SSt0]​[hα​β00hi​j]\displaystyle\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]\left[{\begin{array}[]{cc}0&S\\ S^{t}&0\\ \end{array}}\right]\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]\left[{\begin{array}[]{cc}0&S\\ S^{t}&0\\ \end{array}}\right]\left[{\begin{array}[]{cc}h_{\alpha\beta}&0\\ 0&h_{ij}\\ \end{array}}\right]
+\displaystyle+ O~​(r3).\displaystyle\widetilde{O}(r^{3}).

Notice the third term does not have off-diagonal contributions, so the inverse matrix g−1=(gI​J)g^{-1}=(g^{IJ}) satisfies

(3.75) gi​α\displaystyle g^{i\alpha} =−hi​j​gj​β​hβ​α+O~​(r3)\displaystyle=-h_{ij}g_{j\beta}h_{\beta\alpha}+\widetilde{O}(r^{3})
(3.76) gi​j\displaystyle g^{ij} =hi​j+O~​(r2),\displaystyle=h_{ij}+\widetilde{O}(r^{2}),
(3.77) gα​β\displaystyle g^{\alpha\beta} =hα​β+O~​(r2)=δα​β+O~​(r2),\displaystyle=h_{\alpha\beta}+\widetilde{O}(r^{2})=\delta_{\alpha\beta}+\widetilde{O}(r^{2}),

where we used Lemma 3.4. The definition of ηα\eta_{\alpha} requires ⟨ηα,d​xj⟩=0\langle\eta_{\alpha},dx_{j}\rangle=0, which implies that

(3.78) pi​α​gi​j+gj​α=0.p_{i\alpha}g^{ij}+g^{j\alpha}=0.

Let (g^i​j)(\hat{g}_{ij}) be the inverse of the matrix (gi​j)(g^{ij}) with 1≤i,j≤m−31\leq i,j\leq m-3 such that gi​j​g^j​k=δi​kg^{ij}\hat{g}_{jk}=\delta_{ik}. Multiplying by g^j​k\hat{g}_{jk}, we have

(3.79) pi​α​gi​j​g^j​k+gj​α​g^j​k=0,p_{i\alpha}g^{ij}\hat{g}_{jk}+g^{j\alpha}\hat{g}_{jk}=0,

and hence

(3.80) pk​α=−gj​α​g^j​k.p_{k\alpha}=-g^{j\alpha}\hat{g}_{jk}.

We claim that for any 1≤i,j≤m−31\leq i,j\leq m-3,

(3.81) g^i​j−gi​j=O~​(r2).\hat{g}_{ij}-g_{ij}=\widetilde{O}(r^{2}).

In fact, since

(3.82) gi​j​gj​k+gi​α​gα​k=δi​k.g^{ij}g_{jk}+g^{i\alpha}g_{\alpha k}=\delta_{ik}.

Multipling by the inverse of the submatrix (gi​j)(g^{ij}),

(3.83) g^l​i​(gi​j​gj​k+gi​α​gα​k)=g^l​k.\hat{g}_{li}(g^{ij}g_{jk}+g^{i\alpha}g_{\alpha k})=\hat{g}_{lk}.

So this implies that

(3.84) g^l​k−gl​k=g^l​i​gi​α​gα​k=O~​(r2).\hat{g}_{lk}-g_{lk}=\hat{g}_{li}g^{i\alpha}g_{\alpha k}=\widetilde{O}(r^{2}).

Therefore, combining (3.75),(3.77), (3.80) and (3.84), we obtain

(3.85) pk​α\displaystyle p_{k\alpha} =\displaystyle= −(gk​j+O~​(r2))​gj​α\displaystyle-(g_{kj}+\widetilde{O}(r^{2}))g^{j\alpha}
=\displaystyle= −gk​j​gj​α+O~​(r3)\displaystyle-g_{kj}g^{j\alpha}+\widetilde{O}(r^{3})
=\displaystyle= hα​β​gk​β+O~​(r3)\displaystyle h_{\alpha\beta}g_{k\beta}+\widetilde{O}(r^{3})
=\displaystyle= gk​α+O~​(r3).\displaystyle g_{k\alpha}+\widetilde{O}(r^{3}).

∎

The following symmetry property of pi​αp_{i\alpha} will be frequently used in our later calculations. By Lemma 3.4 and Lemma 3.11, we may write

(3.86) pi​α=Ai​α​β​yβ+12​Bi​α​β​γ​yβ​yγ+O~​(r3).p_{i\alpha}=A_{i\alpha\beta}y_{\beta}+\frac{1}{2}B_{i\alpha\beta\gamma}y_{\beta}y_{\gamma}+\widetilde{O}(r^{3}).

Here the connection term Ai​α​β≡⟨∇∂xi∂yβ,yα⟩|(x,0)A_{i\alpha\beta}\equiv\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\beta}},y_{\alpha}\rangle|_{(x,0)} is skew-symmetric in α\alpha, β\beta, and the curvature term Bi​α​β​γ≡−43​Ri​β​γ​αB_{i\alpha\beta\gamma}\equiv-\frac{4}{3}R_{i\beta\gamma\alpha} is symmetric in β\beta, γ\gamma.

Lemma 3.12.

For every 1≤α,β≤31\leq\alpha,\beta\leq 3,

(3.87) ⟨ηα,ηβ⟩=δα​β+O~​(r2).\langle\eta_{\alpha},\eta_{\beta}\rangle=\delta_{\alpha\beta}+\widetilde{O}(r^{2}).
Proof.

By definition

(3.88) ⟨ηα,ηβ⟩=⟨d​yα+pi​α​d​xi,d​yβ+pj​β​d​xj⟩.\langle\eta_{\alpha},\eta_{\beta}\rangle=\langle dy_{\alpha}+p_{i\alpha}dx_{i},dy_{\beta}+p_{j\beta}dx_{j}\rangle.

By (3.77) we get

(3.89) ⟨d​yα,d​yβ⟩=δα​β+O~​(r2).\langle dy_{\alpha},dy_{\beta}\rangle=\delta_{\alpha\beta}+\widetilde{O}(r^{2}).

Also we have pi​α=O~​(r)p_{i\alpha}=\widetilde{O}(r) and ⟨d​xi,d​yβ⟩=O~​(r)\langle dx_{i},dy_{\beta}\rangle=\widetilde{O}(r) for all ii and α\alpha. The conclusion then follows. ∎

Using the above differential forms ηα\eta_{\alpha}’s, we can decompose the volume form dvolg\dvol_{g} in the horizontal and vertical directions, which will substantially simplify the computations regarding the Hodge Laplacian. The volume form of gg is given by

(3.90) dvolg=det(g)⋅d​y1∧d​y2∧d​y3∧d​x1∧⋯∧d​xm−3,\dvol_{g}=\sqrt{\det(g)}\cdot dy_{1}\wedge dy_{2}\wedge dy_{3}\wedge dx_{1}\wedge\cdots\wedge dx_{m-3},

where we have used the orientation fixed above. We define the normal and tangential volume forms by

(3.91) {dvolN≡η1∧η2∧η3dvolT≡det(gi​jP)⋅d​x1∧⋯∧d​xm−3.\displaystyle\begin{cases}{\dvol}_{N}\equiv\eta_{1}\wedge\eta_{2}\wedge\eta_{3}\\ {\dvol}_{T}\equiv\sqrt{\det(g_{ij}^{P})}\cdot dx_{1}\wedge\cdots\wedge dx_{m-3}.\end{cases}

By the expansion formula (3.86), the normal volume form dvolN\dvol_{N} has the following expansion,

dvolN=d​y1∧d​y2∧d​y3+(Ai​α​β​yβ+12​Bi​α​β​γ​yβ​yγ)​d​xi∧d​yα^\displaystyle{\dvol}_{N}=dy_{1}\wedge dy_{2}\wedge dy_{3}+(A_{i\alpha\beta}y_{\beta}+\frac{1}{2}B_{i\alpha\beta\gamma}y_{\beta}y_{\gamma})dx_{i}\wedge dy_{\widehat{\alpha}}
(3.92) +Ai,α+1,β​Aj,α+2,γ​yβ​yγ​d​yα∧d​xi∧d​xj+O~​(r3).\displaystyle+A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}y_{\beta}y_{\gamma}dy_{\alpha}\wedge dx_{i}\wedge dx_{j}+\widetilde{O}(r^{3}).

In addition, by the definition of ηα\eta_{\alpha}’s, it holds that for each 1≤α≤31\leq\alpha\leq 3,

(3.93) ηα∧dvolT=d​yα∧dvolT,\eta_{\alpha}\wedge\dvol_{T}=dy_{\alpha}\wedge\dvol_{T},

and hence

(3.94) dvolN∧dvolT=d​y1∧d​y2∧d​y3∧dvolT=det(gi​jP)det(g)⋅dvolg.\dvol_{N}\wedge\dvol_{T}=dy_{1}\wedge dy_{2}\wedge dy_{3}\wedge\dvol_{T}=\frac{\sqrt{\det(g_{ij}^{P})}}{\sqrt{\det(g)}}\cdot\dvol_{g}.
Lemma 3.13.

Denote by ∗:Ωk​(Q)→Ωm−k​(Q)*:\Omega^{k}(Q)\to\Omega^{m-k}(Q) the Hodge ∗* operator, then we have the following:

  1. (1)
    (3.95) ∗dvolT=(−1)m+1​dvolN⋅(1−Hα​yα+O~​(r2)),*\dvol_{T}=(-1)^{m+1}\dvol_{N}\cdot(1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2})),
  2. (2)

    For any α∈{1,2,3}\alpha\in\{1,2,3\},

    (3.96) ∗(ηα∧dvolT)=λ1​ηα^+λ2​ηα+1^+λ3​ηα+2^,*(\eta_{\alpha}\wedge\dvol_{T})=\lambda_{1}\eta_{\widehat{\alpha}}+\lambda_{2}\eta_{\widehat{\alpha+1}}+\lambda_{3}\eta_{\widehat{\alpha+2}},

    where λ1=1−Hα​yα+O~​(r2)\lambda_{1}=1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2}), λ2=O~​(r2)\lambda_{2}=\widetilde{O}(r^{2}), λ3=O~​(r2)\lambda_{3}=\widetilde{O}(r^{2}).

Proof.

First, we prove Item (1). By (3.51) we have

(3.97) (−1)m+1∗dvolT=λ⋅dvolN(-1)^{m+1}*\dvol_{T}=\lambda\cdot\dvol_{N}

for a function λ>0\lambda>0. The function λ\lambda is given by

(3.98) λ=|dvolT|2​det(g)det(gi​jP)=det(g)​det(gi​j)​det(gi​jP).\lambda=\frac{|\dvol_{T}|^{2}\sqrt{\det(g)}}{\sqrt{\det(g_{ij}^{P})}}=\sqrt{\det(g)}\det(g^{ij})\sqrt{\det(g_{ij}^{P})}.

Now we compute the expansion of λ\lambda. Applying the expansions of gi​jg_{ij}, gα​βg_{\alpha\beta} and gi​αg_{i\alpha} in Lemma 3.4, one can directly obtain the following,

(3.99) det(g)\displaystyle\det(g) =det(gα​β)⋅det(gi​j)+O~​(r2),\displaystyle=\det(g_{\alpha\beta})\cdot\det(g_{ij})+\widetilde{O}(r^{2}),
(3.100) det(gα​β)\displaystyle\det(g_{\alpha\beta}) =1+O~​(r2),\displaystyle=1+\widetilde{O}(r^{2}),
(3.101) det(gi​j)\displaystyle\det(g_{ij}) =det(gi​jP)⋅(1+2​Hα​yα)+O~​(r2).\displaystyle=\det(g_{ij}^{P})\cdot(1+2H^{\alpha}y_{\alpha})+\widetilde{O}(r^{2}).

Plugging (3.100) and (3.101) into (3.99),

(3.102) det(g)=det(gi​jP)⋅(1+2​Hα​yα)+O~​(r2).\det(g)=\det(g_{ij}^{P})\cdot(1+2H^{\alpha}y_{\alpha})+\widetilde{O}(r^{2}).

Let (hi​j)(h_{ij}) be the inverse of the matrix (gi​j)(g_{ij}). Since gi​j=hi​j+O~​(r2)g^{ij}=h_{ij}+\widetilde{O}(r^{2}) by (3.76), so it follows that

(3.103) det(gi​j)=det(hi​j)+O~​(r2)=(det(gi​j))−1+O~​(r2).\det(g^{ij})=\det(h_{ij})+\widetilde{O}(r^{2})=(\det(g_{ij}))^{-1}+\widetilde{O}(r^{2}).

Plugging (3.101) into the above,

(3.104) det(gi​j)=det(gi​jP)−1⋅(1−2​Hα​yα)+O~​(r2).\det(g^{ij})=\det(g_{ij}^{P})^{-1}\cdot(1-2H^{\alpha}y_{\alpha})+\widetilde{O}(r^{2}).

Therefore, substituting (3.102) and (3.104) into (3.98),

(3.105) λ=1−Hα​yα+O~​(r2),\lambda=1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2}),

which completes the proof of Item (1).

Now we prove Item (2). For each α∈{1,2,3}\alpha\in\{1,2,3\}, we can write

(3.106) ∗(ηα∧dvolT)=λ1⋅ηα^+λ2⋅ηα+1^+λ3⋅ηα+2^.*(\eta_{\alpha}\wedge\dvol_{T})=\lambda_{1}\cdot\eta_{\widehat{\alpha}}+\lambda_{2}\cdot\eta_{\widehat{\alpha+1}}+\lambda_{3}\cdot\eta_{\widehat{\alpha+2}}.

Taking point-wise wedge product with ηα∧dvolT\eta_{\alpha}\wedge\dvol_{T}, and noticing ηα+1^∧ηα\eta_{\widehat{\alpha+1}}\wedge\eta_{\alpha}, ηα+2^∧ηα\eta_{\widehat{\alpha+2}}\wedge\eta_{\alpha} are both zero, then we obtain

(3.107) λ1⋅ηα^∧ηα∧dvolT=(ηα∧dvolT)∧∗(ηα∧dvolT).\lambda_{1}\cdot\eta_{\widehat{\alpha}}\wedge\eta_{\alpha}\wedge\dvol_{T}=(\eta_{\alpha}\wedge\dvol_{T})\wedge*(\eta_{\alpha}\wedge\dvol_{T}).
(3.108) λ1⋅dvolN∧dvolT=|ηα∧dvolT|2​dvolg\displaystyle\lambda_{1}\cdot\dvol_{N}\wedge\dvol_{T}=|\eta_{\alpha}\wedge\dvol_{T}|^{2}\dvol_{g}

Therefore, by (3.51) and (3.87) we get

(3.109) λ1=|ηα∧dvolT|2​det(g)det(gi​jP)=1−Hα​yα+O~​(r2),\lambda_{1}=\frac{|\eta_{\alpha}\wedge\dvol_{T}|^{2}\sqrt{\det(g)}}{\sqrt{\det(g_{ij}^{P})}}=1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2}),

Similarly taking wedge product with ηα+1∧dvolT\eta_{\alpha+1}\wedge\dvol_{T} and ηα+2∧dvolT\eta_{\alpha+2}\wedge\dvol_{T} respectively, and again by (3.87) we obtain that

(3.110) λ2=O~​(r2),λ3=O~​(r2).\lambda_{2}=\widetilde{O}(r^{2}),\lambda_{3}=\widetilde{O}(r^{2}).

These imply that

(3.111) ∗(ηα∧dvolT)=λ1⋅ηα^+λ2⋅ηα+1^+λ3⋅ηα+2^,*(\eta_{\alpha}\wedge\dvol_{T})=\lambda_{1}\cdot\eta_{\widehat{\alpha}}+\lambda_{2}\cdot\eta_{\widehat{\alpha+1}}+\lambda_{3}\cdot\eta_{\widehat{\alpha+2}},

where λ1=1+Hα​yα+O~​(r2)\lambda_{1}=1+H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2}), λ2=O~​(r2)\lambda_{2}=\widetilde{O}(r^{2}) and λ3=O~​(r2)\lambda_{3}=\widetilde{O}(r^{2}). ∎

Now we proceed to prove Theorem 3.10. This will be done in several steps.

Step 1. We start by defining a 3-form

(3.112) ϕ1≡(−1)m+1​12​r∗dvolT.\phi_{1}\equiv(-1)^{m+1}\frac{1}{2r}*\dvol_{T}.

By Item (1) of Lemma 3.13, immediately we have

(3.113) ϕ1=dvolN2​r⋅(1−Hα​yα+O~​(r2)).\displaystyle\phi_{1}=\frac{\dvol_{N}}{2r}\cdot(1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2})).

Then applying the expansion of dvolN\dvol_{N} in (3.92), ϕ1\phi_{1} has a further expansion,

ϕ1\displaystyle\phi_{1} =1−Hα​yα2​r​d​y1∧d​y2∧d​y3+12​r​Ai​α​β​yβ​d​xi∧d​yα^\displaystyle=\frac{1-H^{\alpha}y_{\alpha}}{2r}dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}A_{i\alpha\beta}y_{\beta}dx_{i}\wedge dy_{\widehat{\alpha}}
+12​r​Ai,α+1,β​Aj,α+2,γ​yβ​yγ​d​yα∧d​xi∧d​xj+r−1​Π3(2)+O′​(r2)\displaystyle+\frac{1}{2r}A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}y_{\beta}y_{\gamma}dy_{\alpha}\wedge dx_{i}\wedge dx_{j}+r^{-1}\Pi_{3}^{(2)}+O^{\prime}(r^{2})
=1−Hα​yα2​r​d​y1∧d​y2∧d​y3+12​r​Ai​α​β​yβ​d​xi∧d​yα^\displaystyle=\frac{1-H^{\alpha}y_{\alpha}}{2r}dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}A_{i\alpha\beta}y_{\beta}dx_{i}\wedge dy_{\hat{\alpha}}
(3.114) +14(Ai​α,α+1Aj​α,α+2−Ai​α,α+2Aj​α,α+1)yα⋅dr∧dxi∧dxj+r−1Π3(2)+O′(r2),\displaystyle+\frac{1}{4}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})y_{\alpha}\cdot dr\wedge dx_{i}\wedge dx_{j}+r^{-1}\Pi_{3}^{(2)}+O^{\prime}(r^{2}),

where Π3(2)\Pi_{3}^{(2)} is the 33-form introduced in Notation 3.9 and the last step can be achieved by applying the following lemma:

Lemma 3.14 (Rearrangement Lemma).
(3.115) Ai,α+1,β​Aj,α+2,γ​yβ​yγ​d​yα∧d​xi∧d​xj=12​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​yα⋅r​d​r∧d​xi∧d​xj.A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}y_{\beta}y_{\gamma}dy_{\alpha}\wedge dx_{i}\wedge dx_{j}=\frac{1}{2}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})y_{\alpha}\cdot rdr\wedge dx_{i}\wedge dx_{j}.
Proof.

First by writing out the terms and re-arranging the subscripts and using the skew symmetry of Ai​α​βA_{i\alpha\beta} we get

(3.116) Ai,α+1,β​Aj,α+2,γ​yβ​yγ​d​yα\displaystyle A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}y_{\beta}y_{\gamma}dy_{\alpha}
=\displaystyle= (Ai,α+1,α​Aj,α+2,α​yα2+Ai,α+1,α​Aj,α+2,α+1​yα​yα+1CLOSE\displaystyle\Big(A_{i,\alpha+1,\alpha}A_{j,\alpha+2,\alpha}y_{\alpha}^{2}+A_{i,\alpha+1,\alpha}A_{j,\alpha+2,\alpha+1}y_{\alpha}y_{\alpha+1}
OPEN+Ai,α+1,α+2​Aj,α+2,α​yα​yα+2+Ai,α+1,α+2​Aj,α+2,α+1​yα+1​yα+2)​d​yα\displaystyle+A_{i,\alpha+1,\alpha+2}A_{j,\alpha+2,\alpha}y_{\alpha}y_{\alpha+2}+A_{i,\alpha+1,\alpha+2}A_{j,\alpha+2,\alpha+1}y_{\alpha+1}y_{\alpha+2}\Big)dy_{\alpha}
=\displaystyle= Ai,α,α+1​Aj,α,α+2​yα2​d​yα−Ai​α,α+2​Aj​α,α+1​yα​yα+2​d​yα+2\displaystyle A_{i,\alpha,\alpha+1}A_{j,\alpha,\alpha+2}y_{\alpha}^{2}dy_{\alpha}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1}y_{\alpha}y_{\alpha+2}dy_{\alpha+2}
−Ai,α,α+2​Aj,α,α+1​yα​yα+1​d​yα+1−Ai,α,α+1​Aj,α,α+1​yα​yα+1​d​yα+2\displaystyle-A_{i,\alpha,\alpha+2}A_{j,\alpha,\alpha+1}y_{\alpha}y_{\alpha+1}dy_{\alpha+1}-A_{i,\alpha,\alpha+1}A_{j,\alpha,\alpha+1}y_{\alpha}y_{\alpha+1}dy_{\alpha+2}

Now we can skew-symmetrize with respect to ii and jj

(3.117) Ai,α+1,β​Aj,α+2,γ​yβ​yγ​d​yα∧d​xi∧d​xj=12​(Ai,α+1,β​Aj,α+2,γ−Aj,α+1,β​Ai,α+2,γ)​yβ​yγ​d​yα∧d​xi∧d​xjA_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}y_{\beta}y_{\gamma}dy_{\alpha}\wedge dx_{i}\wedge dx_{j}=\frac{1}{2}(A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}-A_{j,\alpha+1,\beta}A_{i,\alpha+2,\gamma})y_{\beta}y_{\gamma}dy_{\alpha}\wedge dx_{i}\wedge dx_{j}

Correspondingly by skew-symmetrizing each term of (3.116) with respect to ii and jj, we get

(3.118) 12​(Ai,α+1,β​Aj,α+2,γ−Aj,α+1,β​Ai,α+2,γ)​yβ​yγ​d​yα\displaystyle\frac{1}{2}(A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}-A_{j,\alpha+1,\beta}A_{i,\alpha+2,\gamma})y_{\beta}y_{\gamma}dy_{\alpha}
=\displaystyle= 12​(Ai,α,α+1​Aj,α,α+2−Ai​α,α+2​Aj​α,α+1)​yα​(yα​d​yα+yα+1​d​yα+1+yα+2​d​yα+2)\displaystyle\frac{1}{2}(A_{i,\alpha,\alpha+1}A_{j,\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})y_{\alpha}(y_{\alpha}dy_{\alpha}+y_{\alpha+1}dy_{\alpha+1}+y_{\alpha+2}dy_{\alpha+2})
=\displaystyle= 12​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​yα⋅r​d​r\displaystyle\frac{1}{2}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})y_{\alpha}\cdot rdr

∎

Step 2. In this step will explicitly compute the singular (unbounded) terms of Δ​ϕ1\Delta\phi_{1}. Mainly, we will prove the following proposition.

Proposition 3.15.

Let ϕ1\phi_{1} be the 33-form defined in (3.112), then Δ​ϕ1\Delta\phi_{1} has the following expansion,

Δ​ϕ1=\displaystyle\Delta\phi_{1}= −Hα​yα2​r3​d​y1∧d​y2∧d​y3−Ωi​j​α​β​(14​r​d​yα​β^+14​r2​yα​β^​d​r)∧d​xi∧d​xj\displaystyle-\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}-\Omega_{ij\alpha\beta}(\frac{1}{4r}dy_{\widehat{\alpha\beta}}+\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr)\wedge dx_{i}\wedge dx_{j}
(3.119) +Ai​j​α​β​(14​r2​yα​β^​d​r−14​r​d​yα​β^)∧d​xi∧d​xj+r−5​Π3(4)+O′​(1).\displaystyle+A_{ij\alpha\beta}(\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr-\frac{1}{4r}dy_{\widehat{\alpha\beta}})\wedge dx_{i}\wedge dx_{j}+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).
Proof.

The proof consists of two steps.

The first step focuses on the computation for d∗​d​ϕ1d^{*}d\phi_{1}. Starting with the expansion of ϕ1\phi_{1} in (3.113), we have

(3.120) d​ϕ1=−1+Hα​yα+O~​(r2)2​r3⋅r​d​r∧dvolN+12​r​d​((1−Hα​yα+O~​(r2))​dvolN).d\phi_{1}=\frac{-1+H_{\alpha}y_{\alpha}+\widetilde{O}(r^{2})}{2r^{3}}\cdot rdr\wedge\dvol_{N}+\frac{1}{2r}d\Big((1-H^{\alpha}y_{\alpha}+{\widetilde{O}}(r^{2}))\dvol_{N}\Big).

To deal with the first term, we use Lemma 3.11 and (3.13) in Lemma 3.2, then

(3.121) yα​ηα\displaystyle y_{\alpha}\eta_{\alpha} =\displaystyle= yα​d​yα−yα​pi​α​d​xi\displaystyle y_{\alpha}dy_{\alpha}-y_{\alpha}p_{i\alpha}dx_{i}
=\displaystyle= r​d​r−yα​gi​α​d​xi+O~​(r4)\displaystyle rdr-y_{\alpha}g_{i\alpha}dx_{i}+\widetilde{O}(r^{4})
=\displaystyle= r​d​r+O~​(r4),\displaystyle rdr+\widetilde{O}(r^{4}),

which yields

(3.122) r​d​r∧dvolN=(yα​ηα)∧dvolN+O~​(r4)=O~​(r4).rdr\wedge\dvol_{N}=(y_{\alpha}\eta_{\alpha})\wedge\dvol_{N}+\widetilde{O}(r^{4})=\widetilde{O}(r^{4}).

So it follows that

(3.123) d​ϕ1=12​r​d​((1−Hα​yα+O~​(r2))​dvolN)+O′​(r).d\phi_{1}=\frac{1}{2r}d\Big((1-H^{\alpha}y_{\alpha}+\widetilde{O}(r^{2}))\dvol_{N}\Big)+O^{\prime}(r).

It is easy to see that

(3.124) d⁡(O~​(r2)​dvolN)=O~​(r2).d(\widetilde{O}(r^{2})\dvol_{N})=\widetilde{O}(r^{2}).

So we obtain

(3.125) d​ϕ1=12​r​d​((1−Hα​yα)​dvolN)+O′​(r).d\phi_{1}=\frac{1}{2r}d\Big((1-H^{\alpha}y_{\alpha})\dvol_{N}\Big)+O^{\prime}(r).

Next we will compute the expansion for d⁡(dvolN)d(\dvol_{N}). By definition,

(3.126) d⁡(dvolN)=d⁡(η1∧η2∧η3)=d​ηα∧ηα^.d(\dvol_{N})=d(\eta_{1}\wedge\eta_{2}\wedge\eta_{3})=d\eta_{\alpha}\wedge\eta_{\widehat{\alpha}}.

By (3.86),

(3.127) d​ηα\displaystyle d\eta_{\alpha} =\displaystyle= d⁡(pi​α)∧d​xi\displaystyle d(p_{i\alpha})\wedge dx_{i}
=\displaystyle= Ai​α​β​d​yβ∧d​xi+Aj​i​α​β​yβ​d​xj∧d​xi+Bi​α​β​γ​yβ​d​yγ∧d​xi+O~​(r2).\displaystyle A_{i\alpha\beta}dy_{\beta}\wedge dx_{i}+A_{ji\alpha\beta}y_{\beta}dx_{j}\wedge dx_{i}+B_{i\alpha\beta\gamma}y_{\beta}dy_{\gamma}\wedge dx_{i}+\widetilde{O}(r^{2}).

So we have

(3.128) d⁡(dvolN)=Ai​α​β​d​yβ∧d​xi∧ηα^+yβ​(Aj​i​α​β​d​xj∧d​xi+Bi​α​β​γ​d​yγ∧d​xi)∧d​yα^+O~​(r2).d(\dvol_{N})=A_{i\alpha\beta}dy_{\beta}\wedge dx_{i}\wedge\eta_{\widehat{\alpha}}+y_{\beta}(A_{ji\alpha\beta}dx_{j}\wedge dx_{i}+B_{i\alpha\beta\gamma}dy_{\gamma}\wedge dx_{i})\wedge dy_{\widehat{\alpha}}+\widetilde{O}(r^{2}).

Now we need to rearrange the above expansion. Since Ai​α​βA_{i\alpha\beta} is skew symmetric in α\alpha and β\beta, we have for α∈{1,2,3}\alpha\in\{1,2,3\},

(3.129) Ai​α​α=0,A_{i\alpha\alpha}=0,

so the leading order in the first term vanishes, hence

(3.130) Ai​α​β​d​yβ∧d​xi∧ηα^\displaystyle A_{i\alpha\beta}dy_{\beta}\wedge dx_{i}\wedge\eta_{\widehat{\alpha}} =\displaystyle= Ai​α​β​Aj​μ​γ​yγ​d​yβ∧d​xi∧d​xj∧d​yα​μ^+O~​(r2)\displaystyle A_{i\alpha\beta}A_{j\mu\gamma}y_{\gamma}dy_{\beta}\wedge dx_{i}\wedge dx_{j}\wedge dy_{\widehat{\alpha\mu}}+\widetilde{O}(r^{2})
=\displaystyle= Ai​α​β​Aj​β​γ​yγ​d​yα^∧d​xi∧d​xj+O~​(r2)\displaystyle A_{i\alpha\beta}A_{j\beta\gamma}y_{\gamma}dy_{\widehat{\alpha}}\wedge dx_{i}\wedge dx_{j}+\widetilde{O}(r^{2})
=\displaystyle= 12​(Ai​α​β​Aj​β​γ−Ai​γ​β​Aj​β​α)​yγ​d​yα^∧d​xi∧d​xj+O~​(r2).\displaystyle\frac{1}{2}(A_{i\alpha\beta}A_{j\beta\gamma}-A_{i\gamma\beta}A_{j\beta\alpha})y_{\gamma}dy_{\widehat{\alpha}}\wedge dx_{i}\wedge dx_{j}+\widetilde{O}(r^{2}).

Therefore,

(3.131) d⁡(dvolN)=Ωi​j​α​β⋅yβ⋅d​yα^∧d​xi∧d​xj+Bi​α​β​α⋅yβ⋅d​y1∧d​y2∧d​y3∧d​xi+O~​(r2).d(\dvol_{N})=\Omega_{ij\alpha\beta}\cdot y_{\beta}\cdot dy_{\widehat{\alpha}}\wedge dx_{i}\wedge dx_{j}+B_{i\alpha\beta\alpha}\cdot y_{\beta}\cdot dy_{1}\wedge dy_{2}\wedge dy_{3}\wedge dx_{i}+\widetilde{O}(r^{2}).

By (3.92) we have

(3.132) d⁡(Hα​yα)∧dvolN=(Hα​Ai​α​β−∂i(Hβ))⋅yβ⋅d​y1∧d​y2∧d​y3∧d​xi+O~​(r2).d(H^{\alpha}y_{\alpha})\wedge\dvol_{N}=(H^{\alpha}A_{i\alpha\beta}-\partial_{i}(H^{\beta}))\cdot y_{\beta}\cdot dy_{1}\wedge dy_{2}\wedge dy_{3}\wedge dx_{i}+\widetilde{O}(r^{2}).

Now substituting (3.131) and (3.132) into (3.125),

d​ϕ1\displaystyle d\phi_{1} =12​r​Ωi​j​α​β⋅yβ⋅d​yα^∧d​xi∧d​xj\displaystyle=\frac{1}{2r}\Omega_{ij\alpha\beta}\cdot y_{\beta}\cdot dy_{\widehat{\alpha}}\wedge dx_{i}\wedge dx_{j}
(3.133) +12​r​(Bi​α​β​α−(Hα​Ai​α​β−∂i(Hβ))⋅yβ⋅d​y1∧d​y2∧d​y3∧d​xi+O′​(r)CLOSE.\displaystyle+\frac{1}{2r}\Big(B_{i\alpha\beta\alpha}-(H^{\alpha}A_{i\alpha\beta}-\partial_{i}(H^{\beta})\Big)\cdot y_{\beta}\cdot dy_{1}\wedge dy_{2}\wedge dy_{3}\wedge dx_{i}+O^{\prime}(r).

Now we need to take d∗d^{*} of this. Notice that the leading order of d∗​d​ϕ1d^{*}d\phi_{1} can be computed by using the operators in the Euclidean case, so we obtain

(3.134) d∗​d​ϕ1\displaystyle d^{*}d\phi_{1} =\displaystyle= Ωi​j​α​β​(12​r​d​yβ​α^−12​r3​yμ​yβ​d​yμ​α^)​d​xi∧d​xj+r−3​Π3(2)+O′​(1).\displaystyle\Omega_{ij\alpha\beta}(\frac{1}{2r}dy_{\widehat{\beta\alpha}}-\frac{1}{2r^{3}}y_{\mu}y_{\beta}dy_{\widehat{\mu\alpha}})dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{(2)}+O^{\prime}(1).

In our next step, we will compute d​d∗​ϕ1dd^{*}\phi_{1}. First,

(3.135) ∗ϕ1=12​r​dvolT.*\phi_{1}=\frac{1}{2r}\dvol_{T}.

Notice that d⁡(dvolT)=0d(\dvol_{T})=0, so

(3.136) d∗ϕ1=−12​r3⋅rdr∧dvolT.d*\phi_{1}=-\frac{1}{2r^{3}}\cdot rdr\wedge\dvol_{T}.

By (3.121), r​d​r=yα​ηα+O~​(r4)rdr=y_{\alpha}\eta_{\alpha}+{\widetilde{O}}(r^{4}), then

(3.137) d∗ϕ1=−yα2​r3​ηα∧dvolT+O′​(r).d*\phi_{1}=-\frac{y_{\alpha}}{2r^{3}}{}\eta_{\alpha}\wedge\dvol_{T}+O^{\prime}(r).

Applying Item (2) of Lemma 3.13,

(3.138) ∗d∗ϕ1=−yα2​r3​(1−Hβ​yβ+O~​(r2))​ηα^+O′​(r).*d*\phi_{1}=-\frac{y_{\alpha}}{2r^{3}}{}(1-H^{\beta}y_{\beta}+\widetilde{O}(r^{2}))\eta_{\widehat{\alpha}}+O^{\prime}(r).

So it follows that

(3.139) d∗​ϕ1\displaystyle d^{*}\phi_{1} =\displaystyle= −∗d∗ϕ1=yα2​r3(1−Hβyβ)ηα^+O~​(r2)r3yαdyα^+O′(r).\displaystyle-*d*\phi_{1}=\frac{y_{\alpha}}{2r^{3}}{}(1-H^{\beta}y_{\beta})\eta_{\widehat{\alpha}}{}+\frac{\widetilde{O}(r^{2})}{r^{3}}y_{\alpha}dy_{\widehat{\alpha}}+O^{\prime}(r).

Taking dd and applying Lemma 3.2,

d​d∗​ϕ1=\displaystyle dd^{*}\phi_{1}= (1−Hβ​yβ)​(−3​yα2​r5​r​d​r∧ηα^+12​r3​d​(yα​ηα^))−Hβ​yα2​r3​ηα^∧d​yβ+r−5​Π3(4)+O′​(1)\displaystyle(1-H^{\beta}y_{\beta})\Big(-\frac{3y_{\alpha}}{2r^{5}}{}rdr\wedge\eta_{\widehat{\alpha}}+\frac{1}{2r^{3}}d(y_{\alpha}\eta_{\widehat{\alpha}})\Big){}-\frac{H^{\beta}y_{\alpha}}{2r^{3}}\eta_{\widehat{\alpha}}\wedge dy_{\beta}+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1)
=\displaystyle= (1−Hβ​yβ)​(−3​yα2​r5​r​d​r∧ηα^+12​r3​d​(yα​ηα^))−Hα​yα2​r3​d​y1∧d​y2∧d​y3\displaystyle(1-H^{\beta}y_{\beta})\Big(-\frac{3y_{\alpha}}{2r^{5}}{}rdr\wedge\eta_{\widehat{\alpha}}+\frac{1}{2r^{3}}d(y_{\alpha}\eta_{\widehat{\alpha}})\Big)-\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}
(3.140) +\displaystyle+ r−5​Π3(4)+O′​(1).\displaystyle r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

Now we simplify this expression. By (3.121),

(3.141) −3​yα2​r5​r​d​r∧ηα^=−3​yα​yβ2​r5​ηβ∧ηα^+O′​(1)=−32​r3​dvolN+O′​(1).-\frac{3y_{\alpha}}{2r^{5}}{}rdr\wedge\eta_{\widehat{\alpha}}=-\frac{3y_{\alpha}y_{\beta}}{2r^{5}}{}\eta_{\beta}\wedge\eta_{\widehat{\alpha}}+O^{\prime}(1)=-\frac{3}{2r^{3}}{}\dvol_{N}+{O}^{\prime}(1).

Also

(3.142) 12​r3​d​(yα​ηα^)\displaystyle\frac{1}{2r^{3}}d(y_{\alpha}\eta_{\widehat{\alpha}}) =\displaystyle= 12​r3​d​yα∧ηα^+12​r3​yα​d​ηα^\displaystyle\frac{1}{2r^{3}}dy_{\alpha}\wedge\eta_{\widehat{\alpha}}+\frac{1}{2r^{3}}y_{\alpha}d\eta_{\widehat{\alpha}}
=\displaystyle= 12​r3​(ηα−pi​α​d​xi)∧ηα^+12​r3​yα​(d​ηα+1∧ηα+2−ηα+1∧d​ηα+2)\displaystyle\frac{1}{2r^{3}}(\eta_{{\alpha}}-p_{i\alpha}dx_{i})\wedge\eta_{\widehat{\alpha}}+\frac{1}{2r^{3}}y_{\alpha}(d{\eta_{\alpha+1}}\wedge{\eta_{\alpha+2}}-{\eta_{\alpha+1}}\wedge d{\eta_{\alpha+2}})
=\displaystyle= 32​r3​dvolN−12​r3​(Ai​α​β​yβ+12​Bi​α​β​γ​yβ​yγ)​d​xi∧ηα^\displaystyle\frac{3}{2r^{3}}\dvol_{N}-\frac{1}{2r^{3}}(A_{i\alpha\beta}y_{\beta}+\frac{1}{2}B_{i\alpha\beta\gamma}y_{\beta}y_{\gamma})dx_{i}\wedge\eta_{\widehat{\alpha}}
+\displaystyle+ 12​r3​yα​(d​ηα+1∧ηα+2−ηα+1∧d​ηα+2)+O′​(1).\displaystyle\frac{1}{2r^{3}}y_{\alpha}(d{\eta_{\alpha+1}}\wedge{\eta_{\alpha+2}}-{\eta_{\alpha+1}}\wedge d{\eta_{\alpha+2}})+{O^{\prime}}(1).

So it follows that

d​d∗​ϕ1=\displaystyle dd^{*}\phi_{1}= −12​r3​(Ai​α​β⋅yβ⋅d​xi∧ηα^−yα​(d​ηα+1∧ηα+2−ηα+1∧d​ηα+2))\displaystyle-\frac{1}{2r^{3}}\Big(A_{i\alpha\beta}\cdot y_{\beta}\cdot dx_{i}\wedge\eta_{\widehat{\alpha}}-y_{\alpha}(d{\eta_{\alpha+1}}\wedge{\eta_{\alpha+2}}-{\eta_{\alpha+1}}\wedge d{\eta_{\alpha+2}})\Big)
(3.143) −\displaystyle- Hα​yα2​r3​d​y1∧d​y2∧d​y3+r−5​Π3(4)+O′​(1).\displaystyle\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

Next, we will show a crucial cancellation for the first term of the above d​d∗​ϕ1dd^{*}\phi_{1}, which gives a further order improvement.

Lemma 3.16 (Cancellation Lemma).
Ai​α​β⋅yβ⋅d​xi∧ηα^−yα​(d​ηα+1∧ηα+2−ηα+1∧d​ηα+2)\displaystyle A_{i\alpha\beta}\cdot y_{\beta}\cdot dx_{i}\wedge\eta_{\widehat{\alpha}}-y_{\alpha}(d{\eta_{\alpha+1}}\wedge{\eta_{\alpha+2}}-{\eta_{\alpha+1}}\wedge d{\eta_{\alpha+2}})
(3.144) =\displaystyle= −Ai​j​α​β​yβ​yμ​d​yμ​α^∧d​xi∧d​xj+Π3(2)+O~​(r3).\displaystyle-A_{ij\alpha\beta}y_{\beta}y_{\mu}dy_{\widehat{\mu\alpha}}\wedge dx_{i}\wedge dx_{j}+\Pi_{3}^{(2)}+{\widetilde{O}}(r^{3}).
Proof.

Directly applying the definition of ηα\eta_{\alpha}, then we have

Ai​α​β⋅yβ⋅d​xi∧ηα^\displaystyle A_{i\alpha\beta}\cdot y_{\beta}\cdot dx_{i}\wedge\eta_{\widehat{\alpha}}
(3.145) =\displaystyle= Ai​α​β​yβ​d​xi∧d​yα^+Ai​α​β​yβ​yγ​(Aj,α+1,γ​d​yα+2−Aj,α+2,γ​d​yα+1)∧d​xi∧d​xj+O~​(r3).\displaystyle A_{i\alpha\beta}y_{\beta}dx_{i}\wedge dy_{\widehat{\alpha}}+A_{i\alpha\beta}y_{\beta}y_{\gamma}(A_{j,\alpha+1,\gamma}dy_{\alpha+2}-A_{j,\alpha+2,\gamma}dy_{\alpha+1})\wedge dx_{i}\wedge dx_{j}+\widetilde{O}(r^{3}).

By (3.127), we get

yα​(d​ηα+1∧ηα+2−ηα+1∧d​ηα+2)\displaystyle y_{\alpha}(d{\eta_{\alpha+1}}\wedge{\eta_{\alpha+2}}-{\eta_{\alpha+1}}\wedge d{\eta_{\alpha+2}})
=\displaystyle= yα​(Ai,α+1,β​d​yβ∧d​xi∧d​yα+2−Ai,α+2,β​d​yβ∧d​xi∧d​yα+1)\displaystyle y_{\alpha}(A_{i,\alpha+1,\beta}dy_{\beta}\wedge dx_{i}\wedge dy_{\alpha+2}-A_{i,\alpha+2,\beta}dy_{\beta}\wedge dx_{i}\wedge dy_{\alpha+1})
+\displaystyle+ yα​yγ​(Ai,α+1,β​Aj,α+2,γ−Ai,α+2,β​Aj,α+1,γ)​d​yβ∧d​xi∧d​xj\displaystyle y_{\alpha}y_{\gamma}(A_{i,\alpha+1,\beta}A_{j,\alpha+2,\gamma}-A_{i,\alpha+2,\beta}A_{j,\alpha+1,\gamma})dy_{\beta}\wedge dx_{i}\wedge dx_{j}
+\displaystyle+ yα​yβ​(Ai​j,α+1,β​d​yα+2−Ai​j,α+2,β​d​yα+1)∧d​xi∧d​xj\displaystyle y_{\alpha}y_{\beta}(A_{ij,\alpha+1,\beta}dy_{\alpha+2}-A_{ij,\alpha+2,\beta}dy_{\alpha+1})\wedge dx_{i}\wedge dx_{j}
(3.146) +\displaystyle+ Π3(2)+O~​(r3).\displaystyle\Pi_{3}^{(2)}+\widetilde{O}(r^{3}).

Rearranging the subscripts of the first groups of terms in (3.146),

(3.147) yα​(Ai,α+1,β​d​yβ∧d​xi∧d​yα+2−Ai,α+2,β​d​yβ∧d​xi∧d​yα+1)\displaystyle y_{\alpha}(A_{i,\alpha+1,\beta}dy_{\beta}\wedge dx_{i}\wedge dy_{\alpha+2}-A_{i,\alpha+2,\beta}dy_{\beta}\wedge dx_{i}\wedge dy_{\alpha+1})
(3.148) =\displaystyle= yα​(Ai,α+1,α​d​yα∧d​xi∧d​yα+2−Ai,α+2,α​d​yα∧d​xi∧d​yα+1)\displaystyle y_{\alpha}(A_{i,\alpha+1,\alpha}dy_{\alpha}\wedge dx_{i}\wedge dy_{\alpha+2}-A_{i,\alpha+2,\alpha}dy_{\alpha}\wedge dx_{i}\wedge dy_{\alpha+1})
(3.149) =\displaystyle= yα+2​Ai,α,α+2​d​yα+2∧d​xi∧d​yα+1−yα+1​Ai,α,α+1​d​yα+1∧d​xi∧d​yα+2\displaystyle y_{\alpha+2}A_{i,\alpha,\alpha+2}dy_{\alpha+2}\wedge dx_{i}\wedge dy_{\alpha+1}-y_{\alpha+1}A_{i,\alpha,\alpha+1}dy_{\alpha+1}\wedge dx_{i}\wedge dy_{\alpha+2}
(3.150) =\displaystyle= Ai​α​β​yβ​d​xi∧d​yα^,\displaystyle A_{i\alpha\beta}y_{\beta}dx_{i}\wedge dy_{\widehat{\alpha}},

which matches the first term of (3.145). As in the proof of Lemma 3.14, one can see that the second groups of terms in (3.145) and (3.146) are both equal to

(3.151) (Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​yα⋅r​d​r∧d​xi∧d​xj.(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})y_{\alpha}\cdot rdr\wedge dx_{i}\wedge dx_{j}.

Next, the third group of terms in (3.146) can be rewritten as follows,

(3.152) yα​yβ​(Ai​j,α+1,β​d​yα+2−Ai​j,α+2,β​d​yα+1)∧d​xi∧d​xj\displaystyle y_{\alpha}y_{\beta}(A_{ij,\alpha+1,\beta}dy_{\alpha+2}-A_{ij,\alpha+2,\beta}dy_{\alpha+1})\wedge dx_{i}\wedge dx_{j}
=\displaystyle= Ai​j​α​β​yβ​(yα+2​d​yα+1−yα+1​d​yα+2)∧d​xi∧d​xj\displaystyle A_{ij\alpha\beta}y_{\beta}(y_{\alpha+2}dy_{\alpha+1}-y_{\alpha+1}dy_{\alpha+2})\wedge dx_{i}\wedge dx_{j}
=\displaystyle= Ai​j​α​β​yβ​yμ​d​yμ​α^∧d​xi∧d​xj.\displaystyle A_{ij\alpha\beta}y_{\beta}y_{\mu}dy_{\widehat{\mu\alpha}}\wedge dx_{i}\wedge dx_{j}.

The conclusion just follows.

∎

Now we return to the expansion of d​d∗​ϕ1dd^{*}\phi_{1} given by (3.143). Applying Lemma 3.16, finally we obtain

d​d∗​ϕ1\displaystyle dd^{*}\phi_{1} =−Hα​yα2​r3​d​y1∧d​y2∧d​y3+12​r3​Ai​j​α​β​yβ​yμ​d​yμ​α^∧d​xi∧d​xj\displaystyle=-\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r^{3}}A_{ij\alpha\beta}y_{\beta}y_{\mu}dy_{\widehat{\mu\alpha}}\wedge dx_{i}\wedge dx_{j}
(3.153) +r−5​Π3(4)+O′​(1).\displaystyle+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

In the last step of the proof, we will further simplify d∗​d​ϕ1d^{*}d\phi_{1} and d​d∗​ϕ1dd^{*}\phi_{1}. For this purpose, we need the following lemma.

Lemma 3.17.
(3.154) Ωi​j​α​β​(12​r​d​yβ​α^−12​r3​yμ​yβ​d​yμ​α^)\displaystyle\Omega_{ij\alpha\beta}(\frac{1}{2r}dy_{\widehat{\beta\alpha}}-\frac{1}{2r^{3}}y_{\mu}y_{\beta}dy_{\widehat{\mu\alpha}}) =−Ωi​j​α​β​(14​r​d​yα​β^+14​r2​yα​β^​d​r),\displaystyle=-\Omega_{ij\alpha\beta}(\frac{1}{4r}dy_{\widehat{\alpha\beta}}+\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr),
(3.155) 12​r3​Ai​j​α​β​yβ​yμ​d​yμ​α^\displaystyle\frac{1}{2r^{3}}A_{ij\alpha\beta}y_{\beta}y_{\mu}dy_{\widehat{\mu\alpha}} =Ai​j​α​β​(14​r2​yα​β^​d​r−14​r​d​yα​β^).\displaystyle=A_{ij\alpha\beta}(\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr-\frac{1}{4r}dy_{\widehat{\alpha\beta}}).
Proof.

We only prove (3.154) because the other equality follows from the same computations. Using the fact that Ωi​j​α​β=−Ωi​j​β​α\Omega_{ij\alpha\beta}=-\Omega_{ij\beta\alpha}, we can write out the left hand side as

Ωi​j​α​β​(12​r​d​yβ​α^−12​r3​yμ​yβ​d​yμ​α^)\displaystyle\Omega_{ij\alpha\beta}(\frac{1}{2r}dy_{\widehat{\beta\alpha}}-\frac{1}{2r^{3}}y_{\mu}y_{\beta}dy_{\widehat{\mu\alpha}})
=\displaystyle= Ωi​j​α,α+1​(−1r​d​yα+2−12​r3​(yα+2​yα+1​d​yyα+1−yα+12​d​yα+2)+12​r3​(yα2​d​yα+2−yα​yα+2​d​yα))\displaystyle\Omega_{ij\alpha,\alpha+1}\Big(-\frac{1}{r}dy_{\alpha+2}-\frac{1}{2r^{3}}(y_{\alpha+2}y_{\alpha+1}dy_{y_{\alpha+1}}-y_{\alpha+1}^{2}dy_{\alpha+2})+\frac{1}{2r^{3}}(y_{\alpha}^{2}dy_{\alpha+2}-y_{\alpha}y_{\alpha+2}dy_{\alpha})\Big)
=\displaystyle= Ωi​j​α,α+1​(−12​r​d​yα+2−12​r2​yα+2​d​r)\displaystyle\Omega_{ij\alpha,\alpha+1}\Big(-\frac{1}{2r}dy_{\alpha+2}-\frac{1}{2r^{2}}y_{\alpha+2}dr\Big)
(3.156) =\displaystyle= −12​Ωi​j​α,β​(12​r​d​yα​β^+12​r2​yα​β^​d​r).\displaystyle-\frac{1}{2}\Omega_{ij\alpha,\beta}\Big(\frac{1}{2r}dy_{\widehat{\alpha\beta}}+\frac{1}{2r^{2}}y_{\widehat{\alpha\beta}}dr\Big).

∎

Applying the above lemma, now (3.134) and (3.153) can be simplified as follows,

(3.157) d∗​d​ϕ1\displaystyle d^{*}d\phi_{1} =−Ωi​j​α​β​(14​r​d​yα​β^+14​r2​yα​β^​d​r)∧d​xi∧d​xj+r−3​Π3−2+O′​(1),\displaystyle=-\Omega_{ij\alpha\beta}(\frac{1}{4r}dy_{\widehat{\alpha\beta}}+\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr)\wedge dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{-2}+O^{\prime}(1),
d​d∗​ϕ1\displaystyle dd^{*}\phi_{1} =−Hα​yα2​r3​d​y1∧d​y2∧d​y3+Ai​j​α​β​(14​r2​yα​β^​d​r−14​r​d​yα​β^)∧d​xi∧d​xj\displaystyle=-\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}+A_{ij\alpha\beta}(\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr-\frac{1}{4r}dy_{\widehat{\alpha\beta}})\wedge dx_{i}\wedge dx_{j}
(3.158) +r−5​Π3(4)+O′​(1).\displaystyle+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

Therefore,

Δ​ϕ1=\displaystyle\Delta\phi_{1}= (d∗​d+d​d∗)​ϕ1\displaystyle(d^{*}d+dd^{*})\phi_{1}
=\displaystyle= −Hα​yα2​r3​d​y1∧d​y2∧d​y3−Ωi​j​α​β​(14​r​d​yα​β^+14​r2​yα​β^​d​r)∧d​xi∧d​xj\displaystyle-\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}-\Omega_{ij\alpha\beta}(\frac{1}{4r}dy_{\widehat{\alpha\beta}}+\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr)\wedge dx_{i}\wedge dx_{j}
(3.159) +Ai​j​α​β​(14​r2​yα​β^​d​r−14​r​d​yα​β^)∧d​xi∧d​xj+r−5​Π3(4)+O′​(1).\displaystyle+A_{ij\alpha\beta}(\frac{1}{4r^{2}}y_{\widehat{\alpha\beta}}dr-\frac{1}{4r}dy_{\widehat{\alpha\beta}})\wedge dx_{i}\wedge dx_{j}+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

The proof is done.

∎

Step 3. In this step we modify ϕ1\phi_{1} to kill the unbounded terms on the right hand side of (3.119). We first we recall some elementary computations involving the standard Euclidean Hodge Laplacian.

Lemma 3.18.

Let Δ0\Delta_{0} be the standard Hodge Laplacian on the Euclidean space ℝ3\mathbb{R}^{3}, then the following holds:

  1. (1)

    Let {y1,y2,y3}\{y_{1},y_{2},y_{3}\} be the Cartesian coordinates of ℝ3\mathbb{R}^{3}, then

    (3.160) {Δ0​r=−2r,Δ0​(yα​yβr)=4​yα​yβr3,α≠β,Δ0​((yα2r−r))=4​yα2r3,Δ0​(yαr)=2​yαr3.\displaystyle\begin{cases}\Delta_{0}r=-\frac{2}{r},\\ \Delta_{0}(\frac{y_{\alpha}y_{\beta}}{r})=\frac{4y_{\alpha}y_{\beta}}{r^{3}},&\alpha\neq\beta,\\ \Delta_{0}((\frac{y_{\alpha}^{2}}{r}-r))=\frac{4y_{\alpha}^{2}}{r^{3}},\\ \Delta_{0}(\frac{y_{\alpha}}{r})=\frac{2y_{\alpha}}{r^{3}}.\end{cases}
  2. (2)

    Denote by 𝒫4\mathcal{P}_{4} the space of all homogeneous degree 4 polynomials on ℝ3\mathbb{R}^{3}, then the operator

    (3.161) □:𝒫4→𝒫4;f↦r5​Δ0​(r−3​f)\square:\mathcal{P}_{4}\rightarrow\mathcal{P}_{4};f\mapsto r^{5}\Delta_{0}(r^{-3}f)

    is an isomorphism.

Proof.

The first item is a direct calculation. An convenient way to see this is to use the following two facts

  1. (1)

    A homogeneous polynomial degree kk polynomial restricts to an eigenfunction of the Hodge-Laplacian ΔS2\Delta_{S^{2}} on the unit sphere, with eigenvalue k⁡(k+1)k(k+1).

  2. (2)

    Given an eigenfunction hh of ΔS2\Delta_{S^{2}} on the unit sphere with eigenvalue kk, for any ll, we can extend hh to a homogeneous function hlh_{l} on ℝ3∖{0}\mathbb{R}^{3}\setminus\{0\} of degree ll, and

    (3.162) Δ0​hl=r−2​(k−l⁡(l+1))​hl\Delta_{0}h_{l}=r^{-2}(k-l(l+1))h_{l}

For the second item it is possible to write down an explicit inverse to Δ0\Delta_{0}. Here we provide a quick abstract proof. First we notice □:𝒫4→𝒫4\square:\mathcal{P}_{4}\to\mathcal{P}_{4} is a well-defined linear map. This follows from the standard computations

r5​Δ0​(r−3​f)\displaystyle r^{5}\Delta_{0}(r^{-3}f) =\displaystyle= r5Δ0(r−3)⋅f−2r5∇(r−3)⋅∇f+r2Δ0f\displaystyle r^{5}\Delta_{0}(r^{-3})\cdot f-2r^{5}\nabla(r^{-3})\cdot\nabla f+r^{2}\Delta_{0}f
=\displaystyle= −6f−3∇(r2)⋅∇f+r2Δ0f.\displaystyle-6f-3\nabla(r^{2})\cdot\nabla f+r^{2}\Delta_{0}f.

Since each term in the above formula is a polynomial in 𝒫4\mathcal{P}_{4}, so □​f∈𝒫4\square f\in\mathcal{P}_{4}.

Now to prove □\square is an isomorphism it suffices to prove it has a trivial kernel in 𝒫4\mathcal{P}_{4}. Let u≡r−3​fu\equiv r^{-3}f, then u=O⁡(r)u=O(r) for both r→0r\to 0 and r→∞r\to\infty. If Δ0​(u)=0\Delta_{0}(u)=0, then uu is harmonic on ℝ3∖{0}\mathbb{R}^{3}\setminus\{0\}. The removable singularity theorem implies that uu extends smoothly on ℝ3\mathbb{R}^{3}. Since u=O⁡(r)u=O(r) as r→∞r\to\infty, applying the standard derivative estimate for harmonic functions, we conclude ∇2u≡0\nabla^{2}u\equiv 0. Therefore, uu must be a linear function. Noticing f∈𝒫4f\in\mathcal{P}_{4}, we conclude f≡0f\equiv 0. The proof is done.

∎

Next, we want to find a bounded correction 33-form 𝔅0=O′​(1)\mathfrak{B}_{0}=O^{\prime}(1) such that Δ​ϕ1\Delta\phi_{1} is corrected to a bounded term on 𝒰∖P\mathcal{U}\setminus P, i.e.,

(3.163) Δ⁡(ϕ1+𝔅0)=O′​(1)​on​𝒰∖P.\Delta(\phi_{1}+\mathfrak{B}_{0})=O^{\prime}(1)\ \text{on}\ \mathcal{U}\setminus P.

Now the main part is to eliminate the unbounded terms in Δ​ϕ1\Delta\phi_{1} which relies on the following explicit calculations for Δ​𝔅0\Delta\mathfrak{B}_{0}. In fact, the leading terms of Δ​𝔅0\Delta\mathfrak{B}_{0} are exactly given by the Euclidean Laplacian Δ0\Delta_{0} acting on the normal components such that the explicit computations in Lemma 3.18 can be effectively used in our context. Precisely, we have the following lemma.

Lemma 3.19.

Let Δ0\Delta_{0} be the Hodge Laplacian on ℝ3\mathbb{R}^{3}, then the following holds:

  1. (1)

    Denote by νy\nu_{y} one of the following differential forms d​yαdy_{\alpha}, d​yα∧d​yβdy_{\alpha}\wedge dy_{\beta} or d​y1∧d​y2∧d​y3dy_{1}\wedge dy_{2}\wedge dy_{3}. Similarly, let τx\tau_{x} be a tangential pp-form given by τx≡d​x1∧…∧d​xαp\tau_{x}\equiv dx_{1}\wedge\ldots\wedge dx_{\alpha_{p}} with 0≤p≤m−30\leq p\leq m-3. Let

    (3.164) ω≡f⁡(x)​h​(y)​νy∧τx,\omega\equiv f(x)h(y)\nu_{y}\wedge\tau_{x},

    where f⁡(x)f(x) is a smooth function defined on UU and h⁡(y)=O′​(|y|k)h(y)=O^{\prime}(|y|^{k}) for some k∈ℤ+k\in\mathbb{Z}_{+}, then

    (3.165) Δ​ω−f⁡(x)⋅Δ0​(h⁡(y))⋅νy∧τx=O′​(|y|k−1).\Delta\omega-f(x)\cdot\Delta_{0}(h(y))\cdot\nu_{y}\wedge\tau_{x}=O^{\prime}(|y|^{k-1}).
  2. (2)

    Let ff be a smooth function defined on U⊂PU\subset P and let

    (3.166) ω≡f⁡(x)⋅yαr​d​y1∧d​y2∧d​y3,\omega\equiv f(x)\cdot\frac{y_{\alpha}}{r}dy_{1}\wedge dy_{2}\wedge dy_{3},

    then

    (3.167) Δ​ω=2​r−2​ω+r−5​Γ3(4)+O′​(1),\displaystyle\Delta\omega=2r^{-2}\omega+r^{-5}\Gamma_{3}^{(4)}+O^{\prime}(1),

    where the definition of the 33-form Γ3(4)\Gamma_{3}^{(4)} is in (3.45) of Notation 3.9.

Proof.

First, we prove Item (1). By definition, Δ=d​d∗+d∗​d\Delta=dd^{*}+d^{*}d. We only prove the case νy=d​yα\nu_{y}=dy_{\alpha} for 1≤α≤31\leq\alpha\leq 3 and 1≤p≤m−31\leq p\leq m-3. The proof of the remaining cases is identical.

First, we compute d∗​d​ωd^{*}d\omega.

(3.168) d​ω\displaystyle d\omega =f⁡(x)⋅∂h⁡(y)∂yβ​d​yβ∧d​yα∧τx+∂f⁡(x)∂xj⋅h⁡(y)⋅d​xj∧νy∧τx,\displaystyle=f(x)\cdot\frac{\partial h(y)}{\partial y_{\beta}}dy_{\beta}\wedge dy_{\alpha}\wedge\tau_{x}+\frac{\partial f(x)}{\partial x_{j}}\cdot h(y)\cdot dx_{j}\wedge\nu_{y}\wedge\tau_{x},

which implies that

∗d​ω\displaystyle*d\omega =(−1)pf(x)⋅∂h⁡(y)∂yβdyβ​α^∧∗T(τx)+O′(|y|k)\displaystyle=(-1)^{p}f(x)\cdot\frac{\partial h(y)}{\partial y_{\beta}}dy_{\widehat{\beta\alpha}}\wedge*_{T}(\tau_{x})+O^{\prime}(|y|^{k})
(3.169) =(−1)pf(x)(∂h⁡(y)∂yα−1dyα+1−∂h⁡(y)∂yα+1dyα−1)∧∗T(τx)+O′(|y|k).\displaystyle=(-1)^{p}f(x)\Big(\frac{\partial h(y)}{\partial y_{\alpha-1}}dy_{\alpha+1}-\frac{\partial h(y)}{\partial y_{\alpha+1}}dy_{\alpha-1}\Big)\wedge*_{T}(\tau_{x})+O^{\prime}(|y|^{k}).

Differentiating the above equality,

d∗d​ω\displaystyle d*d\omega =(−1)p​f⋅(∂2h∂yα​∂yα−1​d​yα∧d​yα+1−∂2h∂yα​∂yα+1​d​yα∧d​yα−1CLOSE\displaystyle=(-1)^{p}f\cdot\Big(\frac{\partial^{2}h}{\partial y_{\alpha}\partial y_{\alpha-1}}dy_{\alpha}\wedge dy_{\alpha+1}-\frac{\partial^{2}h}{\partial y_{\alpha}\partial y_{\alpha+1}}dy_{\alpha}\wedge dy_{\alpha-1}
(3.170) +(∂2h∂yα−12+∂2h∂yα+12)dyα−1∧dyα+1)∧∗T(τx)+O′(|y|k−1).\displaystyle+\Big(\frac{\partial^{2}h}{\partial y_{\alpha-1}^{2}}+\frac{\partial^{2}h}{\partial y_{\alpha+1}^{2}}\Big)dy_{\alpha-1}\wedge dy_{\alpha+1}\Big)\wedge*_{T}(\tau_{x})+O^{\prime}(|y|^{k-1}).

Then it follows that

d∗​d​ω\displaystyle d^{*}d\omega =(−1)m​p+m+1∗d∗d​ω\displaystyle=(-1)^{mp+m+1}*d*d\omega
=f⋅(∂2h∂yα​∂yα−1​d​yα−1+∂2h∂yα​∂yα+1​d​yα+1−(∂2h∂yα−12+∂2h∂yα+12)​d​yα)∧τx\displaystyle=f\cdot\Big(\frac{\partial^{2}h}{\partial y_{\alpha}\partial y_{\alpha-1}}dy_{\alpha-1}+\frac{\partial^{2}h}{\partial y_{\alpha}\partial y_{\alpha+1}}dy_{\alpha+1}-\Big(\frac{\partial^{2}h}{\partial y_{\alpha-1}^{2}}+\frac{\partial^{2}h}{\partial y_{\alpha+1}^{2}}\Big)dy_{\alpha}\Big)\wedge\tau_{x}
(3.171) +O′​(|y|k−1).\displaystyle+O^{\prime}(|y|^{k-1}).

On the other hand,

(3.172) ∗ω=f(x)h(y)dyα^∧∗T(τx),\displaystyle*\omega=f(x)h(y)dy_{\widehat{\alpha}}\wedge*_{T}(\tau_{x}),

which implies

(3.173) d∗ω=f(x)∂h∂yαdvolN∧∗T(τx)+O′(|y|k).d*\omega=f(x)\frac{\partial h}{\partial y_{\alpha}}\dvol_{N}\wedge*_{T}(\tau_{x})+O^{\prime}(|y|^{k}).

So it follows that

(3.174) d∗ω=(−1)m​p+1∗d∗ω=−f⋅∂h∂yα⋅τx+O′(|y|k),d^{*}\omega=(-1)^{mp+1}*d*\omega=-f\cdot\frac{\partial h}{\partial y_{\alpha}}\cdot\tau_{x}+O^{\prime}(|y|^{k}),

and hence

(3.175) dd∗ω=−f⋅(∂2h∂yα−1​∂yαdyα−1+∂2h∂yα2dyα+∂2h∂yα+1​∂yαdyα+1)⋅dyα∧τx+O′(|y|k−1).dd^{*}\omega=-f\cdot\Big(\frac{\partial^{2}h}{\partial y_{\alpha-1}\partial y_{\alpha}}dy_{\alpha-1}+\frac{\partial^{2}h}{\partial y_{\alpha}^{2}}dy_{\alpha}+\frac{\partial^{2}h}{\partial y_{\alpha+1}\partial y_{\alpha}}dy_{\alpha+1}\Big)\cdot dy_{\alpha}\wedge\tau_{x}+O^{\prime}(|y|^{k-1}).

Therefore, combining (3.171) and (3.175),

(3.176) Δω=(d∗d+dd∗)ω=−f⋅(Δ0h(y))dyα∧τx+O′(|y|k−1),\Delta\omega=(d^{*}d+dd^{*})\omega=-f\cdot(\Delta_{0}h(y))dy_{\alpha}\wedge\tau_{x}+O^{\prime}(|y|^{k-1}),

where Δ0​(h⁡(y))=−∂2h∂y12−∂2h∂y22−∂2h∂y32\Delta_{0}(h(y))=-\frac{\partial^{2}h}{\partial y_{1}^{2}}-\frac{\partial^{2}h}{\partial y_{2}^{2}}-\frac{\partial^{2}h}{\partial y_{3}^{2}}. The proof of (1) is done.

Now we prove Item (2). Let ω=f⋅yαr​d​y1∧d​y2∧d​y3\omega=f\cdot\frac{y_{\alpha}}{r}dy_{1}\wedge dy_{2}\wedge dy_{3} and the first step is to compute the term d∗​d​ωd^{*}d\omega. By Lemma 3.2, d​r=yγ​d​yγrdr=\frac{y_{\gamma}dy_{\gamma}}{r}, then

(3.177) d⁡(yαr)∧d​y1∧d​y2∧d​y3=0.d(\frac{y_{\alpha}}{r})\wedge dy_{1}\wedge dy_{2}\wedge dy_{3}=0.

This implies that

(3.178) d​ω=∂f∂xi⋅yαr​d​xi∧d​y1∧d​y2∧d​y3,d\omega=\frac{\partial f}{\partial x_{i}}\cdot\frac{y_{\alpha}}{r}dx_{i}\wedge dy_{1}\wedge dy_{2}\wedge dy_{3},

and hence

(3.179) ∗dω=−∂f∂xi⋅yαr∗T(dxi)+O′(r).*d\omega=-\frac{\partial f}{\partial x_{i}}\cdot\frac{y_{\alpha}}{r}*_{T}(dx_{i})+O^{\prime}(r).

Differentiating the above equality and applying Lemma 3.2 again,

(3.180) d∗dω=−∂f∂xi(d​yαr−yα​yβ⋅d​yβr3)∗(dxi)+O′(1).\displaystyle d*d\omega=-\frac{\partial f}{\partial x_{i}}\Big(\frac{dy_{\alpha}}{r}-\frac{y_{\alpha}y_{\beta}\cdot dy_{\beta}}{r^{3}}\Big)*(dx_{i})+O^{\prime}(1).

It follows that

(3.181) ∗d∗d​ω=(−1)m+1⋅∂f∂xi⋅d​yα^∧d​xir+(−1)m​∂f∂xi⋅yα​yβr3⋅d​yβ^∧d​xi+O′​(1)\displaystyle*d*d\omega=(-1)^{m+1}\cdot\frac{\partial f}{\partial x_{i}}\cdot\frac{dy_{\widehat{\alpha}}\wedge dx_{i}}{r}+(-1)^{m}\frac{\partial f}{\partial x_{i}}\cdot\frac{y_{\alpha}y_{\beta}}{r^{3}}\cdot dy_{\widehat{\beta}}\wedge dx_{i}+O^{\prime}(1)

Therefore,

d∗​d​ω\displaystyle d^{*}d\omega =(−1)m+1∗d∗d​ω\displaystyle=(-1)^{m+1}*d*d\omega
(3.182) =∂f∂xi⋅d​yα^∧d​xir−∂f∂xi⋅yα​yβr3⋅d​yβ^∧d​xi+O′​(1).\displaystyle=\frac{\partial f}{\partial x_{i}}\cdot\frac{dy_{\widehat{\alpha}}\wedge dx_{i}}{r}-\frac{\partial f}{\partial x_{i}}\cdot\frac{y_{\alpha}y_{\beta}}{r^{3}}\cdot dy_{\widehat{\beta}}\wedge dx_{i}+O^{\prime}(1).

Now we compute d​d∗​ωdd^{*}\omega. By Lemma 3.13 and the expansion of dvolN\dvol_{N} in (3.92),

∗(d​y1∧d​y2∧d​y3)\displaystyle*(dy_{1}\wedge dy_{2}\wedge dy_{3}) =∗(dvolN+Ai​γ​βyγdyγ^∧dxi)+O~(r2)\displaystyle=*(\dvol_{N}+A_{i\gamma\beta}y_{\gamma}dy_{\widehat{\gamma}}\wedge dx_{i})+\widetilde{O}(r^{2})
(3.183) =(1+Hβyβ)dvolT+Ai​γ​βyγdyγ∧∗T(dxi)+O~(r2),\displaystyle=(1+H^{\beta}y_{\beta})\dvol_{T}+A_{i\gamma\beta}y_{\gamma}dy_{\gamma}\wedge*_{T}(dx_{i})+\widetilde{O}(r^{2}),

so we have

∗ω\displaystyle*\omega =f⋅yαr⋅dvolT+f⋅Hβ⋅yα​yβr⋅dvolT+f⋅Ai​γ​βyα​yγrdyγ∧∗T(dxi)+O′(r2)\displaystyle=f\cdot\frac{y_{\alpha}}{r}\cdot\dvol_{T}+f\cdot H^{\beta}\cdot\frac{y_{\alpha}y_{\beta}}{r}\cdot\dvol_{T}+f\cdot A_{i\gamma\beta}\frac{y_{\alpha}y_{\gamma}}{r}dy_{\gamma}\wedge*_{T}(dx_{i})+O^{\prime}(r^{2})
(3.184) ≡𝔗1+𝔗2+𝔗3+O′​(r2).\displaystyle\equiv\FT_{1}+\FT_{2}+\FT_{3}+O^{\prime}(r^{2}).

By collecting the leading terms, it is easy to compute the leading term in the above equality,

(3.185) d∗d⁡(𝔗1)\displaystyle d*d(\FT_{1}) =∂f∂xi​(d​xi∧d​yα^r−yα​yβ​d​xi∧d​yβ^r3)−2​ω+O′​(1)\displaystyle=\frac{\partial f}{\partial x_{i}}\Big(\frac{dx_{i}\wedge dy_{\widehat{\alpha}}}{r}-\frac{y_{\alpha}y_{\beta}dx_{i}\wedge dy_{\widehat{\beta}}}{r^{3}}\Big)-2\omega+O^{\prime}(1)
(3.186) d∗d⁡(𝔗2)\displaystyle d*d(\FT_{2}) =r−5​Π3(4),d∗d⁡(𝔗3)=r−5​Π3(4).\displaystyle=r^{-5}\Pi_{3}^{(4)},\ d*d(\FT_{3})=r^{-5}\Pi_{3}^{(4)}.

Therefore,

(3.187) dd∗ω=−∂f∂xi⋅d​xi∧d​yα^r+∂f∂xi⋅yα​yβr3⋅dxi∧dyβ^+2ω+r−5Π3(4)+O′(1).dd^{*}\omega=-\frac{\partial f}{\partial x_{i}}\cdot\frac{dx_{i}\wedge dy_{\widehat{\alpha}}}{r}+\frac{\partial f}{\partial x_{i}}\cdot\frac{y_{\alpha}y_{\beta}}{r^{3}}\cdot dx_{i}\wedge dy_{\widehat{\beta}}+2\omega+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

By (3.182) and (3.187) we obtain the expansion

(3.188) Δ​ω=(d∗​d+d​d∗)​ω=2​ω+r−5​Π3(4)+O′​(1).\displaystyle\Delta\omega=(d^{*}d+dd^{*})\omega=2\omega+r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

So the proof is done.

∎

Now we finish Step 2 by proving the following

Proposition 3.20.

There is some 33-form Λ3(4)\Lambda_{3}^{(4)} (given in Notation 3.9) such that if we choose

𝔅0\displaystyle\mathfrak{B}_{0} ≡Hα​yα4​r​d​y1∧d​y2∧d​y3\displaystyle\equiv\frac{H^{\alpha}y_{\alpha}}{4r}dy_{1}\wedge dy_{2}\wedge dy_{3}
(3.189) +(Ωi​j​α​β​(116​yα​β^​d​r−316​r​d​yα​β^)−116​Ai​j​α​β​(yα​β^​d​r+r​d​yα​β^))∧d​xi∧d​xj+r−3​Λ3(4),\displaystyle+\Big(\Omega_{ij\alpha\beta}(\frac{1}{16}y_{\widehat{\alpha\beta}}dr-\frac{3}{16}rdy_{\widehat{\alpha\beta}})-\frac{1}{16}A_{ij\alpha\beta}(y_{\widehat{\alpha\beta}}dr+rdy_{\widehat{\alpha\beta}})\Big)\wedge dx_{i}\wedge dx_{j}+r^{-3}\Lambda_{3}^{(4)},

then the corrected 33-form of ϕ1\phi_{1},

(3.190) ϕ2≡\displaystyle\phi_{2}\equiv ϕ1+𝔅0\displaystyle\phi_{1}+\mathfrak{B}_{0}

satisfies

(3.191) Δ​ϕ2=O′​(1)​on​𝒰∖P,\Delta\phi_{2}=O^{\prime}(1)\ \text{on}\ \mathcal{U}\setminus P,

and has the expansion

ϕ2=\displaystyle\phi_{2}= 12​r​(1−Hα​yα2)​d​y1∧d​y2∧d​y3+12​r​yβ​Ai​α​β​d​xi∧d​yα^−14​Ai​j​α​β​r⋅d​yα​β^∧d​xi∧d​xj\displaystyle\frac{1}{2r}(1-\frac{H^{\alpha}y_{\alpha}}{2})dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}y_{\beta}A_{i\alpha\beta}dx_{i}\wedge dy_{\widehat{\alpha}}-\frac{1}{4}A_{ij\alpha\beta}r\cdot dy_{\widehat{\alpha\beta}}\wedge dx_{i}\wedge dx_{j}
(3.192) +316​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​d​(r​yα)∧d​xi∧d​xj+r−3​Π3(4)+O′​(r2).\displaystyle+\frac{3}{16}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})d(ry_{\alpha})\wedge dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{(4)}+O^{\prime}(r^{2}).
Proof.

Let 𝔟0≡Hα​yα4​r​d​y1∧d​y2∧d​y3\mathfrak{b}_{0}\equiv\frac{H^{\alpha}y_{\alpha}}{4r}dy_{1}\wedge dy_{2}\wedge dy_{3}, then Item (2) of Lemma 3.19 tells us that

(3.193) Δ​𝔟0=Hα​yα2​r3​d​y1∧d​y2∧d​y3+r−5​Γ3(4)+O′​(1).\Delta\mathfrak{b}_{0}=\frac{H^{\alpha}y_{\alpha}}{2r^{3}}dy_{1}\wedge dy_{2}\wedge dy_{3}+r^{-5}\Gamma_{3}^{(4)}+O^{\prime}(1).

Let Π(4)\Pi^{(4)} be the 33-form in the expansion of Δ​ϕ1\Delta\phi_{1} given by (3.119) in Proposition 3.15.

Next, Lemma 3.18 and Lemma 3.19 tell us that there are 33-forms Γ^3(4)\widehat{\Gamma}_{3}^{(4)} and Π^3(4)\widehat{\Pi}_{3}^{(4)} which are also of the form as in (3.45) such that

(3.194) Δ⁡(r−3​Γ^3(4))=−r−5​Γ3(4)+O′​(1),\displaystyle\Delta(r^{-3}\widehat{\Gamma}_{3}^{(4)})=-r^{-5}\Gamma_{3}^{(4)}+O^{\prime}(1),
(3.195) Δ⁡(r−3​Π^3(4))=−r−5​Π3(4)+O′​(1).\displaystyle\Delta(r^{-3}\widehat{\Pi}_{3}^{(4)})=-r^{-5}\Pi_{3}^{(4)}+O^{\prime}(1).

Now let

(3.196) 𝔟1≡r−3​Γ^3(4)+r−3​Π^3(4),\mathfrak{b}_{1}\equiv r^{-3}\widehat{\Gamma}_{3}^{(4)}+r^{-3}\widehat{\Pi}_{3}^{(4)},

then the correction term 𝔟0+𝔟1\mathfrak{b}_{0}+\mathfrak{b}_{1} is chosen as the above such that Δ⁡(𝔟0+𝔟1)\Delta(\mathfrak{b}_{0}+\mathfrak{b}_{1}) in fact eliminates the O′​(r−2)O^{\prime}(r^{-2})-term and implicit O′​(r−1)O^{\prime}(r^{-1})-terms in the expansion of Δ​ϕ1\Delta\phi_{1} (see Proposition 3.15).

In the following, we will make a further correction such that those explicit O′​(r−1)O^{\prime}(r^{-1})-terms will be cancelled out as well. In fact, we define

(3.197) 𝔟2≡(Ωi​j​α​β​(116​yα​β^​d​r−316​r​d​yα​β^)−116​Ai​j​α​β​(yα​β^​d​r+r​d​yα​β^))∧d​xi∧d​xj,\mathfrak{b}_{2}\equiv\Big(\Omega_{ij\alpha\beta}(\frac{1}{16}y_{\widehat{\alpha\beta}}dr-\frac{3}{16}rdy_{\widehat{\alpha\beta}})-\frac{1}{16}A_{ij\alpha\beta}(y_{\widehat{\alpha\beta}}dr+rdy_{\widehat{\alpha\beta}})\Big)\wedge dx_{i}\wedge dx_{j},

applying Lemma 3.18 and Lemma 3.19 again, then

(3.198) Δ​𝔟2=Ωi​j​α​β​(14​r​d​yβ​α^+yα​β^4​r2​d​r)∧d​xi∧d​xj−Ai​j​α​β​(yα​β^4​r2​d​r−14​r​d​yα​β^)∧d​xi∧d​xj,\Delta\mathfrak{b}_{2}=\Omega_{ij\alpha\beta}\Big(\frac{1}{4r}dy_{\widehat{\beta\alpha}}+\frac{y_{\widehat{\alpha\beta}}}{4r^{2}}dr\Big)\wedge dx_{i}\wedge dx_{j}-A_{ij\alpha\beta}\Big(\frac{y_{\widehat{\alpha\beta}}}{4r^{2}}dr-\frac{1}{4r}dy_{\widehat{\alpha\beta}}\Big)\wedge dx_{i}\wedge dx_{j},

and hence

(3.199) Δ⁡(ϕ1+𝔟0+𝔟1+𝔟2)=O′​(1).\Delta(\phi_{1}+\mathfrak{b}_{0}+\mathfrak{b}_{1}+\mathfrak{b}_{2})=O^{\prime}(1).

Therefore, it suffices to choose the correction term

(3.200) 𝔅0≡𝔟0+𝔟1+𝔟2,\mathfrak{B}_{0}\equiv\mathfrak{b}_{0}+\mathfrak{b}_{1}+\mathfrak{b}_{2},

which gives Δ⁡(ϕ1+𝔅0)=O′​(1)\Delta(\phi_{1}+\mathfrak{B}_{0})=O^{\prime}(1).

Notice that, 𝔟2\mathfrak{b}_{2} has a further cancellation,

𝔟2=\displaystyle\mathfrak{b}_{2}= (Ωi​j​α​β​(116​yα​β^​d​r−316​r​d​yα​β^)−116​Ai​j​α​β​(yα​β^​d​r+r​d​yα​β^))∧d​xi∧d​xj,\displaystyle\Big(\Omega_{ij\alpha\beta}(\frac{1}{16}y_{\widehat{\alpha\beta}}dr-\frac{3}{16}rdy_{\widehat{\alpha\beta}})-\frac{1}{16}A_{ij\alpha\beta}(y_{\widehat{\alpha\beta}}dr+rdy_{\widehat{\alpha\beta}})\Big)\wedge dx_{i}\wedge dx_{j},
=\displaystyle= −14​Ai​j​α​β​r​d​yα​β^∧d​xi∧d​xj−116​(Ai,α,α+1​Aj,α,α+2−Ai,α,α+2​Aj,α,α+1)​yα​d​r∧d​xi∧d​xj\displaystyle-\frac{1}{4}A_{ij\alpha\beta}rdy_{\widehat{\alpha\beta}}\wedge dx_{i}\wedge dx_{j}-\frac{1}{16}(A_{i,\alpha,\alpha+1}A_{j,\alpha,\alpha+2}-A_{i,\alpha,\alpha+2}A_{j,\alpha,\alpha+1})y_{\alpha}dr\wedge dx_{i}\wedge dx_{j}
(3.201) +316​(Ai,α,α+1​Aj,α,α+2−Ai,α,α+2​Aj,α,α+1)​r​d​yα∧d​xi∧d​xj.\displaystyle+\frac{3}{16}(A_{i,\alpha,\alpha+1}A_{j,\alpha,\alpha+2}-A_{i,\alpha,\alpha+2}A_{j,\alpha,\alpha+1})rdy_{\alpha}\wedge dx_{i}\wedge dx_{j}.

Therefore,

ϕ2=\displaystyle\phi_{2}= ϕ1+𝔅0\displaystyle\phi_{1}+\mathfrak{B}_{0}
=\displaystyle= 12​r​(1−Hα​yα2)​d​y1∧d​y2∧d​y3+12​r​yβ​Ai​α​β​d​xi∧d​yα^−14​Ai​j​α​β​r⋅d​yα​β^∧d​xi∧d​xj\displaystyle\frac{1}{2r}(1-\frac{H^{\alpha}y_{\alpha}}{2})dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}y_{\beta}A_{i\alpha\beta}dx_{i}\wedge dy_{\widehat{\alpha}}-\frac{1}{4}A_{ij\alpha\beta}r\cdot dy_{\widehat{\alpha\beta}}\wedge dx_{i}\wedge dx_{j}
(3.202) +316​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​d​(r​yα)∧d​xi∧d​xj+r−3​Π3(4)+O′​(r2),\displaystyle+\frac{3}{16}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})d(ry_{\alpha})\wedge dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{(4)}+O^{\prime}(r^{2}),

and

(3.203) Δ​ϕ2=O′​(1).\Delta\phi_{2}=O^{\prime}(1).

∎

Step 4. In this step we compute Δ​ϕ2\Delta\phi_{2} as a current on 𝒰\mathcal{U}.

Lemma 3.21.

In 𝒰\mathcal{U}, we have

(3.204) Δ​ϕ2=2​π​δU+O′​(1).\Delta\phi_{2}=2\pi\delta_{U}+O^{\prime}(1).
Proof.

Suppose we are given a compactly supported test form χ∈Ω0m−3​(𝒰)\chi\in\Omega_{0}^{m-3}(\mathcal{U}), then we apply integration by parts once and we have

(3.205) (ϕ2,Δ​χ)=∫𝒰ϕ2∧(d​d∗+d∗​d)​χ=∫𝒰d​ϕ2∧d∗​χ−∫𝒰d∗​ϕ2∧𝑑χ.(\phi_{2},\Delta\chi)=\int_{\mathcal{U}}\phi_{2}\wedge(dd^{*}+d^{*}d)\chi=\int_{\mathcal{U}}d\phi_{2}\wedge d^{*}\chi-\int_{\mathcal{U}}d^{*}\phi_{2}\wedge d\chi.

Here there is no boundary term because ϕ2=O⁡(r−1)\phi_{2}=O(r^{-1}). Notice that (3.133) and (3.190) implies d​ϕ2=O′​(1)d\phi_{2}=O^{\prime}(1), so

(3.206) ∫𝒰d​ϕ2∧d∗​χ=∫𝒰d∗​d​ϕ2∧χ.\int_{\mathcal{U}}d\phi_{2}\wedge d^{*}\chi=\int_{\mathcal{U}}d^{*}d\phi_{2}\wedge\chi.

On the other hand, by (3.139),

(3.207) d∗​ϕ2=−yα2​r3​d​yα^+ζ,d^{*}\phi_{2}=-\frac{y_{\alpha}}{2r^{3}}{}dy_{\hat{\alpha}}+\zeta,

where ζ\zeta is a 22-form satisfying ζ=O′​(r−1)\zeta=O^{\prime}(r^{-1}). Denote by Sϵ2S_{\epsilon}^{2} the normal geodesic sphere bundle {r=ϵ}\{r=\epsilon\}, then we get that

(3.208) −∫𝒰d∗ϕ2∧dχ\displaystyle-\int_{\mathcal{U}}d^{*}\phi_{2}\wedge d\chi =\displaystyle= ∫𝒰d​d∗​ϕ2∧χ+∫Sϵ2d∗​ϕ2∧χ\displaystyle\int_{\mathcal{U}}dd^{*}\phi_{2}\wedge\chi+\int_{S_{\epsilon}^{2}}d^{*}\phi_{2}\wedge\chi
=\displaystyle= ∫𝒰d​d∗​ϕ2∧χ+limϵ→012​ϵ3​∫Sϵ(yα​d​yα^+ϵ2​ζ)∧χ.\displaystyle\int_{\mathcal{U}}dd^{*}\phi_{2}\wedge\chi+\lim_{\epsilon\rightarrow 0}\frac{1}{2\epsilon^{3}}{}\int_{S_{\epsilon}}(y_{\alpha}dy_{\hat{\alpha}}+\epsilon^{2}\zeta)\wedge\chi.

By direct calculation of the last term on the right hand side we obtain

(3.209) −∫𝒰d∗ϕ2∧dχ=∫𝒰dd∗ϕ2∧χ+2π∫Pχ.-\int_{\mathcal{U}}d^{*}\phi_{2}\wedge d\chi=\int_{\mathcal{U}}dd^{*}\phi_{2}\wedge\chi+2\pi{}\int_{P}\chi.

This concludes the proof. ∎

Step 5. Now we solve the Laplace equation with right hand side in O′​(1).O^{\prime}(1).

Lemma 3.22.

Given a local 3-form vv defined on a neighborhood 𝒰\mathcal{U} of pp in QQ with v=O′​(1)v=O^{\prime}(1), then there is some smaller neighborhood 𝒱⊂⊂𝒰\mathcal{V}\subset\subset\mathcal{U} such that there exists a local solution TT to the equation

(3.210) Δ​T=v,\Delta T=v,

with T=O′​(r2)T=O^{\prime}(r^{2}). Here rr is the distance to the submanifold PP.

Proof.

We just need to establish the following:

  1. (1)

    (General derivatives estimate) For each k∈ℕk\in\mathbb{N} and ϵ>0\epsilon>0, it holds that

    (3.211) |∇k+2T|=O⁡(r−(k+ϵ)).|\nabla^{k+2}T|=O(r^{-(k+\epsilon)}).
  2. (2)

    (Mixed derivatives estimate) For each k∈ℕk\in\mathbb{N}, ℓ∈ℕ\ell\in\mathbb{N} and ϵ>0\epsilon>0, it holds that

    (3.212) |(∇t)k​∇ℓ+2T|=O⁡(r−(ℓ+ϵ)),|(\nabla^{t})^{k}\nabla^{\ell+2}T|=O(r^{-(\ell+\epsilon)}),

    where ∇t\nabla^{t} denotes the tangential derivative.

The above estimates will be proved by induction.

First, we will prove the following order estimate for ∇2T\nabla^{2}T in a smaller neighborhood 𝒰′⊂⊂𝒰\mathcal{U}^{\prime}\subset\subset\mathcal{U}

(3.213) |∇2T|=O⁡(r−ϵ).|\nabla^{2}T|=O(r^{-\epsilon}).

This can be viewed as the base step for carrying out the inductive argument.

To begin with, by definition, for any p>0p>0, we have |v|Lp​(𝒰)≤Cp|v|_{L^{p}(\mathcal{U})}\leq C_{p}. Applying the standard elliptic W2,pW^{2,p}-estimate, for each p>0p>0, there is some constant Cp>0C_{p}>0 such that in a smaller neighborhood 𝒰1⊂⊂𝒰\mathcal{U}_{1}\subset\subset\mathcal{U} such that

(3.214) ‖∇2T‖Lp​(𝒰1)≤Cp.\|\nabla^{2}T\|_{L^{p}(\mathcal{U}_{1})}\leq C_{p}.

Then Sobolev embedding theorem tells us that

(3.215) T∈W2,p​(𝒰1)∩C1,α​(𝒰1)T\in W^{2,p}(\mathcal{U}_{1})\cap C^{1,\alpha}(\mathcal{U}_{1})

for any p>1p>1 and 0<α<10<\alpha<1.

To prove (3.213), we need to differentiate the equation, which schematically yields that

(3.216) Δ​∇2T+Q1∗∇3T+Q2∗∇2T+Q3∗∇T=w,\Delta\nabla^{2}T+Q_{1}*\nabla^{3}T+Q_{2}*\nabla^{2}T+Q_{3}*\nabla T=w,

where |w|=O′​(r−2)|w|=O^{\prime}(r^{-2}) and QiQ_{i}’s are smooth terms arising from differentiating the coefficients of Δ\Delta. The above equation can be viewed as an elliptic system in terms of the Hessian of TT. Let Φ≡∇2T\Phi\equiv\nabla^{2}T, noticing ∇T∈Cα​(𝒰1)\nabla T\in C^{\alpha}(\mathcal{U}_{1}), so the terms involving ∇T\nabla T can be absorbed to the right hand side of the equation. Then Φ\Phi can be treated as vector valued functions, once we fix a local frame. So it follows that

(3.217) Δ​Φ+Q1∗∇Φ+Q2∗Φ=η,\Delta\Phi+Q_{1}*\nabla\Phi+Q_{2}*\Phi=\eta,

where |η|=O′​(r−2)|\eta|=O^{\prime}(r^{-2}). Since η\eta has unbounded LpL^{p}-norm for large pp, the standard W2,pW^{2,p}-estimate for Φ\Phi does not directly apply.

For improving the regularity of ∇2T\nabla^{2}T, we will rescale the metric gg. For each xx in an even smaller neighborhood 𝒰2\mathcal{U}_{2} with r⁡(x)=rxr(x)=r_{x}, we rescale the metric gg in Brx​(x)B_{r_{x}}(x) by letting

(3.218) g~=(rx)−2​g,\tilde{g}=(r_{x})^{-2}g,

then the following equation holds in the rescaled geodesic ball B1g~​(x)B_{1}^{\tilde{g}}(x),

(3.219) Δ~​Φ~+rx​Q~1∗∇~​Φ~+(rx)2​Q~2∗Φ~=η~,|η~|=O′​(1),\widetilde{\Delta}\widetilde{\Phi}+r_{x}\widetilde{Q}_{1}*\widetilde{\nabla}\widetilde{\Phi}+(r_{x})^{2}\widetilde{Q}_{2}*\widetilde{\Phi}=\tilde{\eta},\quad|\tilde{\eta}|=O^{\prime}(1),

where Φ~​(y)=Φ⁡(rx⋅y)\widetilde{\Phi}(y)=\Phi(r_{x}\cdot y) for each y∈B1g~​(x)y\in B_{1}^{\tilde{g}}(x). In the above equation, all the coefficients are uniformly bounded independent of xx. Since we have shown |Φ|Lp≤Cp|\Phi|_{L^{p}}\leq C_{p} in (3.214), so simple rescaling gives rise to the following estimate for any p>0p>0,

(3.220) |Φ~|Lp​(B1g~​(x))≤Cp⋅(rx)−np.|\widetilde{\Phi}|_{L^{p}(B_{1}^{\tilde{g}}(x))}\leq C_{p}\cdot(r_{x})^{-\frac{n}{p}}.

Now applying the W2,pW^{2,p}-estimate for Φ~\widetilde{\Phi}, then for each p>0p>0

(3.221) |Φ~|W2,p​(B1/2g~​(x))≤Cp⋅(rx)−np|\widetilde{\Phi}|_{W^{2,p}(B_{1/2}^{\tilde{g}}(x))}\leq C_{p}\cdot(r_{x})^{-\frac{n}{p}}

with Cp>0C_{p}>0 independent of xx. By the Sobolev embedding

(3.222) |Φ~|C1,α​(B1/4g~​(x))≤Cα,p⋅(rx)−np|\widetilde{\Phi}|_{C^{1,\alpha}(B_{1/4}^{\tilde{g}}(x))}\leq C_{\alpha,p}\cdot(r_{x})^{-\frac{n}{p}}

with Cα,p>0C_{\alpha,p}>0 independent of xx. Scale back to the original metric, for any p>0p>0, there is some Cα,p>0C_{\alpha,p}>0 independent of the base point x∈𝒰2x\in\mathcal{U}_{2} such that

(3.223) |Φ|L∞​(Brx/4​(x))+|rx∇Φ|L∞​(Brx/4​(x))≤Cα,p⋅(rx)−np.|\Phi|_{L^{\infty}(B_{r_{x}/4}(x))}+|r_{x}\nabla\Phi|_{L^{\infty}(B_{r_{x}/4}(x))}\leq C_{\alpha,p}\cdot(r_{x})^{-\frac{n}{p}}.

This completes the proof of (3.213).

Now we will finish the proof of Item (1) by using the induction. Based on (3.223), the key induction step is to prove the following: Given any ℓ∈ℤ+\ell\in\mathbb{Z}_{+}, if for each ϵ>0\epsilon>0 and 0≤k≤ℓ−10\leq k\leq\ell-1,

(3.224) |∇kΦ|=O⁡(r−(k+ϵ)),|∇~k​Φ~|W2,p​(B1/2g~​(x))≤Ck,p,ϵ⋅(rx)−ϵ,|\nabla^{k}\Phi|=O(r^{-(k+\epsilon)}),\quad|\widetilde{\nabla}^{k}\widetilde{\Phi}|_{W^{2,p}(B_{1/2}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}\cdot(r_{x})^{-\epsilon},

then for each ϵ>0\epsilon>0, we have

(3.225) |∇ℓΦ|=O⁡(r−(ℓ+ϵ)),|∇~ℓ​Φ~|W2,p​(B1/4g~​(x))≤Cℓ,p,ϵ⋅(rx)−ϵ.|\nabla^{\ell}\Phi|=O(r^{-(\ell+\epsilon)}),\quad|\widetilde{\nabla}^{\ell}\widetilde{\Phi}|_{W^{2,p}(B_{1/4}^{\tilde{g}}(x))}\leq C_{\ell,p,\epsilon}\cdot(r_{x})^{-\epsilon}.

Indeed, then differentiating (3.216) by ∇k\nabla^{k},

(3.226) Δ⁡(∇ℓΦ)+∑j=1ℓ+1Qj∗∇jΦ=wℓ,\Delta(\nabla^{\ell}\Phi)+\sum\limits_{j=1}^{\ell+1}Q_{j}*\nabla^{j}\Phi=w_{\ell},

where |wℓ|=O′​(r−(ℓ+2))|w_{\ell}|=O^{\prime}(r^{-(\ell+2)}). As before, we rescale the metric gg by taking g~=(rx)−2​g\tilde{g}=(r_{x})^{-2}g, then

(3.227) Δ~​(∇~ℓ​Φ~)+∑j=1ℓ+1(rx)ℓ−j+2⋅Qj∗∇jΦ~=w~ℓ,\widetilde{\Delta}(\widetilde{\nabla}^{\ell}\widetilde{\Phi})+\sum\limits_{j=1}^{\ell+1}(r_{x})^{\ell-j+2}\cdot Q_{j}*\nabla^{j}\widetilde{\Phi}=\tilde{w}_{\ell},

where |w~ℓ|=O′​(1)|\tilde{w}_{\ell}|=O^{\prime}(1). Let k=ℓ−1k=\ell-1, applying the induction hypothesis (3.224) and Sobolev embedding, we have

(3.228) |∇~ℓ​Φ~|L∞​(B1/2g~​(x))=O⁡((rx)−ϵ).|\widetilde{\nabla}^{\ell}\widetilde{\Phi}|_{L^{\infty}(B_{1/2}^{\tilde{g}}(x))}=O((r_{x})^{-\epsilon}).

The above enables us to apply the W2,pW^{2,p}-elliptic estimate, so we obtain the following estimate for each ϵ>0\epsilon>0,

(3.229) |∇~ℓ​Φ~|W2,p​(B1/4g~​(x))≤Cℓ,p,ϵ⋅(rx)−ϵ,|\widetilde{\nabla}^{\ell}\widetilde{\Phi}|_{W^{2,p}(B_{1/4}^{\tilde{g}}(x))}\leq C_{\ell,p,\epsilon}\cdot(r_{x})^{-\epsilon},

where Cℓ,p,ϵ>0C_{\ell,p,\epsilon}>0 is independent of the base point xx. Applying the Sobolev embedding W2,p⊂C1,αW^{2,p}\subset C^{1,\alpha} and scaling back to the original metric gg,

(3.230) |∇ℓΦ|L∞​(Brx/8​(x))≤Cℓ,ϵ⋅(rx)−(ℓ+ϵ)|\nabla^{\ell}\Phi|_{L^{\infty}(B_{r_{x}/8}(x))}\leq C_{\ell,\epsilon}\cdot(r_{x})^{-(\ell+\epsilon)}

for each ϵ>0\epsilon>0. So we complete the proof of Item (1).

Now we are ready to finish the proof of Item (2). We only focus on the case ℓ=0\ell=0 and the case for ℓ>0\ell>0 can be directly achieved by applying the above rescaling arguments. To this end, we need the following claim for the tangential derivatives estimate.

Claim. Let |Φ~|W2,p​(B1/2g~​(x))≤Cp,ϵ⋅(rx)−ϵ|\widetilde{\Phi}|_{W^{2,p}(B_{1/2}^{\tilde{g}}(x))}\leq C_{p,\epsilon}\cdot(r_{x})^{-\epsilon} for any ϵ>0\epsilon>0 and for any x∈𝒰x\in\mathcal{U}. Assume that Φ~\widetilde{\Phi} solves the elliptic equation

(3.231) Δ​Φ~+Q1∗∇Φ~+Q2∗Φ~=ζ,\Delta\widetilde{\Phi}+Q_{1}*\nabla\widetilde{\Phi}+Q_{2}*\widetilde{\Phi}=\zeta,

where QiQ_{i}’s are smooth coefficients, |ζ|∈O′​(1)|\zeta|\in O^{\prime}(1). Then for any x∈𝒰x\in\mathcal{U}, the estimate

(3.232) |(∇t)k​Φ~|W2,p​(B1/8g~​(x))≤Ck,p,ϵ⋅(rx)−ϵ,|(\nabla^{t})^{k}\widetilde{\Phi}|_{W^{2,p}(B_{1/8}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}\cdot(r_{x})^{-\epsilon},

holds for all k≥1k\geq 1, p>1p>1 and ϵ>0\epsilon>0.

Taking the first tangential derivative ∇t\nabla^{t} for Δ​Φ~\Delta\widetilde{\Phi},

(3.233) (∇t)​(Δ​Φ~)=Δ​∇tΦ~+Q1∗∇2Φ~+Q2∗∇Φ~,(\nabla^{t})(\Delta\widetilde{\Phi})=\Delta\nabla^{t}\widetilde{\Phi}+Q_{1}*\nabla^{2}\widetilde{\Phi}+Q_{2}*\nabla\widetilde{\Phi},

where QiQ_{i}’s are smooth functions. Hence differentiating (3.231) once by the tangential derivative ∇t\nabla^{t}, we have

(3.234) Δ⁡(∇tΦ~)+Q1∗∇2Φ~+Q2∗∇Φ~=w1,\Delta(\nabla^{t}\widetilde{\Phi})+Q_{1}*\nabla^{2}\widetilde{\Phi}+Q_{2}*\nabla\widetilde{\Phi}=w_{1},

where QiQ_{i}’s are smooth functions, w1≡∇tζw_{1}\equiv\nabla^{t}\zeta and |w1|∈O′​(1)|w_{1}|\in O^{\prime}(1). Since we have already assumed |∇2Φ~|Lp​(B1/2g~​(x))≤Cp,ϵ⋅(rx)−ϵ|\nabla^{2}\widetilde{\Phi}|_{L^{p}(B_{1/2}^{\tilde{g}}(x))}\leq C_{p,\epsilon}\cdot(r_{x})^{-\epsilon} for all ϵ>0\epsilon>0, applying the standard W2,pW^{2,p}-estimate, then for any p>1p>1 and ϵ>0\epsilon>0,

(3.235) |∇tΦ~|W2,p​(B1/3g~​(x))≤Cp,ϵ⋅(rx)−ϵ.|\nabla^{t}\widetilde{\Phi}|_{W^{2,p}(B_{1/3}^{\tilde{g}}(x))}\leq C_{p,\epsilon}\cdot(r_{x})^{-\epsilon}.

Now we prove the higher order mixed derivatives estimate by induction. Repeat taking the tangential derivatives and let Ψ(k)≡(∇t)k​Φ~\Psi^{(k)}\equiv(\nabla^{t})^{k}\widetilde{\Phi} for all k>1k>1. Assume that |Ψ(j)|W2,p​(B1/4g~​(x))≤Ck,p,ϵ​(rx)−ϵ|\Psi^{(j)}|_{W^{2,p}(B_{1/4}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}(r_{x})^{-\epsilon} holds for all j≤k−1j\leq k-1 and p>1p>1, then

(3.236) Δ⁡(Ψ(k))+∑μ=1,21≤μ+ν≤kQμ​ν∗∇μ(Ψ(ν))=wk,\Delta(\Psi^{(k)})+\sum_{\begin{subarray}{c}\mu=1,2\\ 1\leq\mu+\nu\leq k\end{subarray}}Q_{\mu\nu}*\nabla^{\mu}(\Psi^{(\nu)})=w_{k},

where wk≡(∇t)k​ζw_{k}\equiv(\nabla^{t})^{k}\zeta and |wk|=O′​(1)|w_{k}|=O^{\prime}(1). Applying the induction hypothesis, it follows that for any k≥1k\geq 1,

(3.237) |Δ⁡(Ψ(k))|Lp​(B1/4g~​(x))≤Ck,p,ϵ⋅(rx)−ϵ|\Delta(\Psi^{(k)})|_{L^{p}(B_{1/4}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}\cdot(r_{x})^{-\epsilon}

Therefore, for any k≥1k\geq 1, p>1p>1 and ϵ>0\epsilon>0, there is some constant Ck,p,ϵ>0C_{k,p,\epsilon}>0 such that

(3.238) |(∇t)k​Φ~|W2,p​(B1/8g~​(x))=|Ψ(k)|W2,p​(B1/8g~​(x))≤Ck,p,ϵ⋅(rx)−ϵ.|(\nabla^{t})^{k}\widetilde{\Phi}|_{W^{2,p}(B_{1/8}^{\tilde{g}}(x))}=|\Psi^{(k)}|_{W^{2,p}(B_{1/8}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}\cdot(r_{x})^{-\epsilon}.

This completes the proof of the claim.

Now we are in a position to finish the proof of the lemma by completing the induction arguments for Item (2). As before, for any xx, under the rescaled metric g~=(rx)−2​g\tilde{g}=(r_{x})^{-2}g, we start with the equation for Φ~\widetilde{\Phi} in the rescaled geodesic ball B1g~​(x)B_{1}^{\tilde{g}}(x),

(3.239) Δ~​Φ~+rx​Q~1∗∇~​Φ~+(rx)2​Q~2∗Φ~=η~, 0<rx<1,\widetilde{\Delta}\widetilde{\Phi}+r_{x}\widetilde{Q}_{1}*\widetilde{\nabla}\widetilde{\Phi}+(r_{x})^{2}\widetilde{Q}_{2}*\widetilde{\Phi}=\tilde{\eta},\ 0<r_{x}<1,

where Q~i\widetilde{Q}_{i}’s are smooth functions and |η~|=O′​(1)|\tilde{\eta}|=O^{\prime}(1). The above claim tells us that for any k≥1k\geq 1, p>1p>1 and ϵ>0\epsilon>0,

(3.240) |(∇t)k​Φ~|W2,p​(B1/2g~​(x))≤Ck,p,ϵ⋅(rx)−ϵ,|(\nabla^{t})^{k}\widetilde{\Phi}|_{W^{2,p}(B_{1/2}^{\tilde{g}}(x))}\leq C_{k,p,\epsilon}\cdot(r_{x})^{-\epsilon},

where Ck,p,ϵ>0C_{k,p,\epsilon}>0 is independent of xx. Applying the Sobolev embedding, then for any 0<α<10<\alpha<1,

(3.241) |(∇t)k​Φ~|C1,α​(B1/4g~​(x))≤Ck,α,ϵ⋅(rx)−ϵ.|(\nabla^{t})^{k}\widetilde{\Phi}|_{C^{1,\alpha}(B_{1/4}^{\tilde{g}}(x))}\leq C_{k,\alpha,\epsilon}\cdot(r_{x})^{-\epsilon}.

In particular, |(∇t)k​Φ~|=O⁡(r−ϵ)|(\nabla^{t})^{k}\widetilde{\Phi}|=O(r^{-\epsilon}) for all ϵ>0\epsilon>0. Notice that the above estimate is independent of the choice of xx. Rescaling back to the original metric, then for each ϵ>0\epsilon>0,

(3.242) |(∇t)k​∇2T|=O⁡(r−ϵ).|(\nabla^{t})^{k}\nabla^{2}T|=O(r^{-\epsilon}).

The proof of the lemma is done.

∎

With all the above preparations, now we are ready to finish the proof of Theorem 3.10.

Proof of Theorem 3.10.

Let ϕ2\phi_{2} be the 33-current defined in Lemma 3.21 such that

(3.243) Δ​ϕ2=2​π​δU+𝔅1,\Delta\phi_{2}=2\pi\delta_{U}+\mathfrak{B}_{1},

with 𝔅1=O′​(1)\mathfrak{B}_{1}=O^{\prime}(1). So Lemma 3.22 implies that there is some 33-current ℜ=O′​(r2)\mathfrak{R}=O^{\prime}(r^{2}) such that

(3.244) Δ​ℜ=𝔅1\Delta\mathfrak{R}=\mathfrak{B}_{1}

and hence the 33-current

(3.245) GU≡ϕ2+ℜG_{U}\equiv\phi_{2}+\mathfrak{R}

satisfies the equation

(3.246) Δ​GU=2​π​δU.\Delta G_{U}=2\pi\delta_{U}.

Moreover, by Lemma 3.21, GUG_{U} has the expansion

GU=\displaystyle G_{U}= 12​r​(1−Hα​yα2)​d​y1∧d​y2∧d​y3+12​r​yβ​Ai​α​β​d​xi∧d​yα^−14​Ai​j​α​β​r⋅d​yα​β^∧d​xi∧d​xj\displaystyle\frac{1}{2r}(1-\frac{H^{\alpha}y_{\alpha}}{2})dy_{1}\wedge dy_{2}\wedge dy_{3}+\frac{1}{2r}y_{\beta}A_{i\alpha\beta}dx_{i}\wedge dy_{\widehat{\alpha}}-\frac{1}{4}A_{ij\alpha\beta}r\cdot dy_{\widehat{\alpha\beta}}\wedge dx_{i}\wedge dx_{j}
(3.247) +316​(Ai​α,α+1​Aj​α,α+2−Ai​α,α+2​Aj​α,α+1)​d​(r​yα)∧d​xi∧d​xj+r−3​Π3(4)+O′​(r2).\displaystyle+\frac{3}{16}(A_{i\alpha,\alpha+1}A_{j\alpha,\alpha+2}-A_{i\alpha,\alpha+2}A_{j\alpha,\alpha+1})d(ry_{\alpha})\wedge dx_{i}\wedge dx_{j}+r^{-3}\Pi_{3}^{(4)}+O^{\prime}(r^{2}).

The proof of Theorem 3.10 is done.

∎

For our purpose later, we also need the following lemma.

Lemma 3.23.

Let Δ\Delta denote the Hodge Laplacian, then

(3.248) Δ⁡(r2)=−6+O~​(r2).\Delta(r^{2})=-6+\widetilde{O}(r^{2}).
Proof.

This follows from similar, and simpler arguments as above. First,

(3.249) d​r2=2​r​d​r=2​yα​ηα+O~​(r4),dr^{2}=2rdr=2y_{\alpha}\eta_{\alpha}+\widetilde{O}(r^{4}),

so it follows that

(3.250) ∗d​r2=2​yα​ηα^∧dvolT+O~​(r3)*dr^{2}=2y_{\alpha}\eta_{\widehat{\alpha}}\wedge\dvol_{T}+\widetilde{O}(r^{3})

and

(3.251) d∗d​r2=6​dvolN∧dvolT+O~​(r2)=6​dvolg+O~​(r2).d*dr^{2}=6\dvol_{N}\wedge\dvol_{T}+\widetilde{O}(r^{2})=6\dvol_{g}+\widetilde{O}(r^{2}).

Hence

(3.252) Δ(r2)=d∗dr2=−∗d∗dr2=−6+O~(r2).\Delta(r^{2})=d^{*}dr^{2}=-*d*dr^{2}=-6+\widetilde{O}(r^{2}).

∎

Original mathematics by the credited authors. Source-backed reader collection; mathematical self-containment is not assessed.