ScalingStacks

3. Green’s currents [04ZC]

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. Green’s currents

In this section, we study in detail some existence and regularity theory of Green’s currents. Our main motivation for studying these arises from Section 2, where we see the Green’s currents appear as Dirac type singular solutions to the linearization of dimension reduced Calabi-Yau equation by the S1S^{1}-symmetry. It is possible that this study will also have applications to other geometric problems, especially to those concerning adiabatic limits.

This Section is organized as follows. In Section 3.1 we recall the generalized geodesic normal coordinates for an an embedded submanifold. In Section 3.2 we discuss the definition, local existence and regularity properties of Green’s currents in the general Riemannian setting. In Section 3.3, we refine these results in the special case related to Kähler geometry. In Section 3.4, we will prove a global existence result which will be immediately used in Section 4.

3.1. Normal coordinates for an embedded submanifold

We start our discussion by introducing the basic notions of the normal exponential map and normal coordinates with respect to an embedded submanifold. This part seems to be standard in Riemannian geometry. For the consistence of the notations and the completeness of the paper, here we include detailed discussions and proofs.

Let (Q,g)(Q,g) be an oriented Riemannian manifold of dimension mm and let P⊂QP\subset Q be a closed embedded oriented submanifold of codimension k0k_{0} in QQ. In our later applications, we only need the case k0=3k_{0}=3. Denote by N{N} the normal bundle of PP in QQ, equipped with the induced fiberwise Riemannian inner products. For any p∈Pp\in P, N⁡(p)N(p) denotes the fiber of pp in NN. The normal exponential map of PP in QQ is defined as

(3.1) ExpP:N→Q,(p,v)↦Expp⁡(v),v∈N⁡(p),\Exp_{P}:{N}\rightarrow Q,\ (p,v)\mapsto\Exp_{p}(v),\ v\in N(p),

where Expp:Tp​Q→Q\Exp_{p}:T_{p}Q\to Q is the standard exponential map at p∈Qp\in Q. By standard implicit function theorem, it is straightforward that ExpP:N→Q\Exp_{P}:{N}\to Q is a diffeomorphism from some neighborhood of the zero section of N{N} to some tubular neighborhood of PP in QQ.

Now we define the normal coordinates. Fix p∈Pp\in P. First we choose local coordinates {x1′,…,xm−k0′}\{x_{1}^{\prime},\ldots,x_{m-k_{0}}^{\prime}\} in a small neighborhood U⊂PU\subset P of pp such that

(3.2) ⟨∂xi′,∂xj′⟩=δi​j, 1≤i,j≤m−k0,\langle\partial_{x_{i}^{\prime}},\partial_{x_{j}^{\prime}}\rangle=\delta_{ij},\ 1\leq i,j\leq m-k_{0},

at pp. We also assume that ∂x1′∧⋯∧∂xm−k0′\partial_{x_{1}^{\prime}}\wedge\cdots\wedge\partial_{x_{m-k_{0}}^{\prime}} is compatible with the orientation on UU. Next, we pick local orthonormal sections {e1,…,ek0}\{e_{1},\ldots,e_{k_{0}}\} of the normal bundle NN such that

(3.3) ⟨eα,eβ⟩=δα​β, 1≤α,β≤k0,\langle e_{\alpha},e_{\beta}\rangle=\delta_{\alpha\beta},\ 1\leq\alpha,\beta\leq k_{0},

on UU. Again we assume that e1∧⋯∧ek0e_{1}\wedge\cdots\wedge e_{k_{0}} is compatible with the orientation on NN, i.e. e1∧⋯∧ek0∧∂x1′∧⋯∂xk−m0′e_{1}\wedge\cdots\wedge e_{k_{0}}\wedge\partial_{x_{1}^{\prime}}\wedge\cdots\partial_{x_{k-m_{0}}^{\prime}} is compatible with the orientation on QQ. Then we can find ϵ>0\epsilon>0 such that Expp\Exp_{p} restricts to a diffeomorphism from

(3.4) 𝔖0={(q,v)∈N|q∈U,|v|<ϵ}\mathfrak{S}_{0}=\{(q,v)\in N|q\in U,\ |v|<\epsilon\}

to a neighborhood 𝒰\mathcal{U} of pp. In particular, U=𝒰∩PU=\mathcal{U}\cap P.

Definition 3.1 (Normal coordinates).

For any (q,v)∈𝔖0(q,v)\in\mathfrak{S}_{0} with v=∑α=1k0vα​eαv=\sum\limits_{\alpha=1}^{k_{0}}v_{\alpha}e_{\alpha}, the local normal coordinates are defined as follows

(3.5) {xj​(ExpP⁡(q,v))≡xj′​(q),1≤j≤m−k0,yβ​(ExpP⁡(q,v))≡vβ,1≤β≤k0.\displaystyle\begin{cases}x_{j}(\Exp_{P}(q,v))\equiv x_{j}^{\prime}(q),&1\leq j\leq m-k_{0},\\ y_{\beta}(\Exp_{P}(q,v))\equiv v_{\beta},&1\leq\beta\leq k_{0}.\end{cases}

By definition for each fixed point (y1,…,yk0,x1,…,xm−k0)(y_{1},\ldots,y_{k_{0}},x_{1},\ldots,x_{m-k_{0}}) in 𝒰\mathcal{U}, the curve

(3.6) ϑ⁡(t)≡(t​y1,…,t​yk0,x1,…,xm−k0),\vartheta(t)\equiv(ty_{1},\ldots,ty_{k_{0}},x_{1},\ldots,x_{m-k_{0}}),

is a normal geodesic which is orthogonal to U⊂PU\subset P. Let

(3.7) ∂y1,…,∂yk0,∂x1,…,∂xm−k0\partial_{y_{1}},\ldots,\partial_{y_{k_{0}}},\partial_{x_{1}},\ldots,\partial_{x_{m-k_{0}}}

be the induced coordinate vector fields. Then ∂yα|y=0=eα\partial_{y_{\alpha}}|_{y=0}=e_{\alpha} when both are viewed as sections of NN over U⊂PU\subset P. For 1≤i,j≤m−k01\leq i,j\leq m-k_{0} and 1≤α≤k01\leq\alpha\leq k_{0}, we denote

(3.8) gi​j≡⟨∂xi,∂xj⟩,gi​α≡⟨∂xi,∂yα⟩,gα​β≡⟨∂yα,∂yβ⟩.g_{ij}\equiv\langle\partial_{x_{i}},\partial_{x_{j}}\rangle,\ g_{i\alpha}\equiv\langle\partial_{x_{i}},\partial_{y_{\alpha}}\rangle,\ g_{\alpha\beta}\equiv\langle\partial_{y_{\alpha}},\partial_{y_{\beta}}\rangle.

By definition, we have

(3.9) gi​j​(0,0)=δi​j,gα​β​(x,0)=δα​β,gi​α​(x,0)=0,g_{ij}(0,0)=\delta_{ij},\ g_{\alpha\beta}(x,0)=\delta_{\alpha\beta},\ g_{i\alpha}(x,0)=0,

for all x∈U⊂Px\in U\subset P. The second fundamental form of the embedding U↪𝒰U\hookrightarrow\mathcal{U} can be written as II=IIi​jα∂yα⊗(dxi⊗dxj)\IIs=\IIs_{ij}^{\alpha}\partial_{y_{\alpha}}\otimes(dx_{i}\otimes dx_{j}), where

(3.10) IIi​jα≡⟨∇∂xieα,∂xj⟩\IIs^{\alpha}_{ij}\equiv\langle\nabla_{\partial_{x_{i}}}e_{\alpha},\partial_{x_{j}}\rangle

Denote by H→=Hα​eα\overrightarrow{H}=H^{\alpha}e_{\alpha} the mean curvature vector, then

(3.11) Hα≡gi​j​IIi​jα.H^{\alpha}\equiv g^{ij}\IIs^{\alpha}_{ij}.

In the above coordinates, we define the normal distance function

(3.12) r≡|y|=(∑α=1k0yα2)12.r\equiv|y|=\Big(\sum\limits_{\alpha=1}^{k_{0}}y_{\alpha}^{2}\Big)^{\frac{1}{2}}.

A straightforward extension of the usual Gauss Lemma gives the following and we omit the proof.

Lemma 3.2 (Generalized Gauss Lemma).

For any p∈Pp\in P, there is some sufficiently small neighborhood U⊂PU\subset P of pp and a tubular neighborhood 𝒰⊂Q\mathcal{U}\subset Q with U=𝒰∩PU=\mathcal{U}\cap P such that the function rr defined by (3.12) satisfies the following properties:

  1. (1)

    ∇r=∂r\nabla r=\partial_{r} holds in 𝒰∖U\mathcal{U}\setminus U. In particular, rr is the normal distance function in 𝒰\mathcal{U}, i.e., r⁡(q)=d⁡(q,P)r(q)=d(q,P) for all q∈𝒰q\in\mathcal{U}.

  2. (2)

    ∂r\partial_{r} is orthogonal to ∂xi\partial_{x_{i}}’s in 𝒰⊂Q\mathcal{U}\subset Q, and hence

    (3.13) ⟨∂xi,r∂r⟩=⟨∂xi,∑α=1k0yα∂yα⟩=∑α=1k0yαgi​α=0, 1≤i≤m−k0.\langle\partial_{x_{i}},r\partial_{r}\rangle=\langle\partial_{x_{i}},\sum_{\alpha=1}^{k_{0}}y_{\alpha}\partial_{y_{\alpha}}\rangle=\sum_{\alpha=1}^{k_{0}}y_{\alpha}g_{i\alpha}=0,\ 1\leq i\leq m-k_{0}.

For the convenience of later discussion, we introduce several notations concerning the normal regularity order near the submanifold PP. It will be frequently used throughout the paper.

Definition 3.3 (Normal regularity order).

Let T⁡(x,y)T(x,y) be a tensor locally defined in 𝒰\mathcal{U} which is C∞C^{\infty} on 𝒰∖U\mathcal{U}\setminus U, then for a non-negative integer kk we say as r→0r\rightarrow 0

  1. (1)

    T⁡(x,y)=O′​(rk)T(x,y)=O^{\prime}(r^{k}) if for each ϵ>0\epsilon>0

    (3.14) |∂xI∂yJT⁡(x,y)|={O⁡(r−ϵ),|J|≤k,O⁡(rk−|J|−ϵ),|J|>k,\displaystyle\Big|\partial_{x}^{I}\partial_{y}^{J}T(x,y)\Big|=\begin{cases}O(r^{-\epsilon}),&|J|\leq k,\\ O(r^{k-|J|-\epsilon}),&|J|>k,\end{cases}

    for all multi-indices II and JJ. In particular, if T∈C∞​(𝒰)T\in C^{\infty}(\mathcal{U}), then T=O′​(rk)T=O^{\prime}(r^{k}) for all k∈ℤk\in\mathbb{Z}.

  2. (2)

    T⁡(x,y)=rk​O′​(1)T(x,y)=r^{k}O^{\prime}(1) if r−k​T​(x,y)=O′​(1)r^{-k}T(x,y)=O^{\prime}(1).

  3. (3)

    T​(x,y)=O~​(rk)T(x,y)=\widetilde{O}(r^{k}) if T⁡(x,y)∈C∞​(𝒰)T(x,y)\in C^{\infty}(\mathcal{U}) and

    (3.15) T⁡(x,y)=rk​O′​(1).T(x,y)=r^{k}O^{\prime}(1).

    In other words, T⁡(x,y)T(x,y) is smooth in 𝒰\mathcal{U} and has vanishing normal derivatives along UU up to order k−1k-1.

Notice that the defining condition does not depend on the choice of the local coordinates, since if we have another coordinate system {y~1,⋯,y~k0,x~1,⋯,x~m−k0}\{\tilde{y}_{1},\cdots,\tilde{y}_{k_{0}},\tilde{x}_{1},\cdots,\tilde{x}_{m-k_{0}}\}, then we have

(3.16) ∂y~α=∂yβ∂y~α⋅∂yβ+∂xl∂y~α⋅∂xl,\displaystyle\partial_{\tilde{y}_{\alpha}}=\frac{\partial y_{\beta}}{\partial\tilde{y}_{\alpha}}\cdot\partial_{y_{\beta}}+\frac{\partial x_{l}}{\partial\tilde{y}_{\alpha}}\cdot\partial_{x_{l}},
(3.17) ∂x~j=∂xl∂x~j⋅∂xl+∂yα∂x~j⋅∂yα,\displaystyle\partial_{\tilde{x}_{j}}=\frac{\partial x_{l}}{\partial\tilde{x}_{j}}\cdot\partial_{x_{l}}+\frac{\partial y_{\alpha}}{\partial\tilde{x}_{j}}\cdot\partial_{y_{\alpha}},

and ∂yα∂x~j∈r​O′​(1)\frac{\partial y_{\alpha}}{\partial\tilde{x}_{j}}\in rO^{\prime}(1). Similarly, we can also use any local coordinate system {yα,xi}\{y_{\alpha},x_{i}\} such that yα=0y_{\alpha}=0 along PP for 1≤α≤k01\leq\alpha\leq k_{0}.

Lemma 3.4.

In the above normal coordinates, we have the following expansions of the metric tensor gg of QQ along the normal directions,

(3.18) gα​β|(x,y)\displaystyle g_{\alpha\beta}|_{(x,y)} =δα​β−13​Rmα​γ​ξ​β|(x,0)​yγ​yξ+O~​(r3),\displaystyle=\delta_{\alpha\beta}-\frac{1}{3}\Rm_{\alpha\gamma\xi\beta}\Big|_{(x,0)}y_{\gamma}y_{\xi}+\widetilde{O}(r^{3}),
(3.19) gi​j|(x,y)\displaystyle g_{ij}|_{(x,y)} =gi​jP(x)+2IIi​jα|(x,0)yα−(Rmi​γ​ξ​j+⟨∇∂xi∂yγ,∇∂xj∂yξ⟩)|(x,0)yγyξ+O~(r3),\displaystyle=g^{P}_{ij}(x)+2\IIs^{\alpha}_{ij}\Big|_{(x,0)}y_{\alpha}-(\Rm_{i\gamma\xi j}+\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\gamma}},\nabla_{\partial_{x_{j}}}\partial_{y_{\xi}}\rangle)\Big|_{(x,0)}y_{\gamma}y_{\xi}+\widetilde{O}(r^{3}),
(3.20) gi​α|(x,y)\displaystyle g_{i\alpha}|_{(x,y)} =⟨∇∂xi∂yγ,∂yα⟩|(x,0)yγ−23Rmi​γ​ξ​α|(x,0)yγyξ+O~(r3),\displaystyle=\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\gamma}},\partial_{y_{\alpha}}\rangle\Big|_{(x,0)}y_{\gamma}-\frac{2}{3}\Rm_{i\gamma\xi\alpha}\Big|_{(x,0)}y_{\gamma}y_{\xi}+\widetilde{O}(r^{3}),

where gP=(gi​jP)g^{P}=(g^{P}_{ij}) denotes the restriction of the metric gg to U⊂PU\subset P, Rm\Rm denotes the Riemann curvature tensor of gg.

Proof.

The above expansions can be proved using the Jacobi fields. Fix a point q=(0k0,x1,…,xm−k0)∈U⊂Pq=(0^{k_{0}},x_{1},\ldots,x_{m-k_{0}})\in U\subset P, we we choose a unit vector v=∑α=1m−k0vα∂yα∈N(q)≅ℝk0v=\sum\limits_{\alpha=1}^{m-k_{0}}v_{\alpha}\partial_{y_{\alpha}}\in N(q)\cong\mathbb{R}^{k_{0}} with |v|=1|v|=1. Let ϑ\vartheta be the following radial geodesic in 𝒰\mathcal{U},

(3.21) ϑ⁡(t)=ExpP⁡(q,t​v)≡Expq⁡(t​v)\vartheta(t)=\Exp_{P}(q,tv)\equiv\Exp_{q}(tv)

such that ϑ′​(0)=v\vartheta^{\prime}(0)=v. In the normal coordinates, the geodesic ϑ\vartheta can represented as ϑ⁡(t)=(t​v1,…,t​vk0,x1,…,xm−k0)\vartheta(t)=(tv_{1},\ldots,tv_{k_{0}},x_{1},\ldots,x_{m-k_{0}}).

For each 1≤α≤k01\leq\alpha\leq k_{0} and 1≤i≤m−k01\leq i\leq m-k_{0}, we define the geodesic variations

(3.22) σα​(t,s)\displaystyle\sigma_{\alpha}(t,s) ≡((t​v1,…,t⁡(vα+s),…,t​vk0,x1,…,xm−k0)CLOSE,\displaystyle\equiv((tv_{1},\ldots,t(v_{\alpha}+s),\ldots,tv_{k_{0}},x_{1},\ldots,x_{m-k_{0}}),
(3.23) σi​(t,s)\displaystyle\sigma_{i}(t,s) ≡(t​v1,…,t​vk0,x1,…,xi+s,…​xm−k0).\displaystyle\equiv(tv_{1},\ldots,tv_{k_{0}},x_{1},\ldots,x_{i}+s,\ldots x_{m-k_{0}}).

Then variation fields of σα​(t,s)\sigma_{\alpha}(t,s) and σi​(t,s)\sigma_{i}(t,s) give the following Jacobi fields along the radial geodesic ϑ⁡(t)\vartheta(t) respectively:

(3.24) {Jα(t)=t⋅∂yα,1≤α≤k0Ji(t)=∂xi,1≤i≤m−k0.\displaystyle\begin{cases}J_{\alpha}(t)=t\cdot\partial_{y_{\alpha}},&1\leq\alpha\leq k_{0}\\ J_{i}(t)=\partial_{x_{i}},&1\leq i\leq m-k_{0}.\end{cases}

By definition,

(3.25) Jα(0)=0,Ji(0)=∂xi.J_{\alpha}(0)=0,\ J_{i}(0)=\partial_{x_{i}}.

Taking first derivatives at t=0t=0,

(3.26) Jα′(0)=∂yα,Ji′(0)=vα∇∂xi∂yα.J_{\alpha}^{\prime}(0)=\partial_{y_{\alpha}},\ J_{i}^{\prime}(0)=v_{\alpha}\nabla_{\partial_{x_{i}}}\partial_{y_{\alpha}}.

Then applying the Jacobi equation along the geodesic ϑ\vartheta,

(3.27) {Jα′′+Rm⁡(Jα,ϑ′)​ϑ′=0,Ji′′+Rm⁡(Ji,ϑ′)​ϑ′=0,\displaystyle\begin{cases}J_{\alpha}^{\prime\prime}+\Rm(J_{\alpha},\vartheta^{\prime})\vartheta^{\prime}=0,\\ J_{i}^{\prime\prime}+\Rm(J_{i},\vartheta^{\prime})\vartheta^{\prime}=0,\\ \end{cases}

where Rm⁡(X,Y)​Z≡∇X∇Y​Z−∇Y∇X​Z−∇[X,Y]Z\Rm(X,Y)Z\equiv\nabla_{X}\nabla_{Y}Z-\nabla_{Y}\nabla_{X}Z-\nabla_{[X,Y]}Z denotes the Riemann curvature tensor of gg, so it follows that

(3.28) Jα′′(0)=0,Jα′′′(0)=−vγvξRm(∂yα,∂yγ)∂yξ,\displaystyle J^{\prime\prime}_{\alpha}(0)=0,\ J^{\prime\prime\prime}_{\alpha}(0)=-v_{\gamma}v_{\xi}\Rm(\partial_{y_{\alpha}},\partial_{y_{\gamma}})\partial_{y_{\xi}},
(3.29) Ji′′(0)=−vγvξRm(∂xi,∂yγ)∂yξ.\displaystyle J^{\prime\prime}_{i}(0)=-v_{\gamma}v_{\xi}\Rm(\partial_{x_{i}},\partial_{y_{\gamma}})\partial_{y_{\xi}}.

Therefore,

(3.30) gα​β\displaystyle g_{\alpha\beta} =t−2​g​(Jα,Jβ)=δα​β−13​Rmα​γ​ξ​β​vγ​vξ​t2+O~​(t3),\displaystyle=t^{-2}g(J_{\alpha},J_{\beta})=\delta_{\alpha\beta}-\frac{1}{3}\Rm_{\alpha\gamma\xi\beta}v_{\gamma}v_{\xi}t^{2}+\widetilde{O}(t^{3}),
(3.31) gi​j\displaystyle g_{ij} =g(Ji,Jj)=gi​jP+2vαIIi​jαt−(Rmi​γ​ξ​j+⟨∇∂xi∂yγ,∇∂xj∂yξ⟩)vγvξt2+O~(t3),\displaystyle=g(J_{i},J_{j})=g^{P}_{ij}+2v_{\alpha}\IIs_{ij}^{\alpha}t-\Big(\Rm_{i\gamma\xi j}+\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\gamma}},\nabla_{\partial_{x_{j}}}\partial_{y_{\xi}}\rangle\Big)v_{\gamma}v_{\xi}t^{2}+\widetilde{O}(t^{3}),
(3.32) gi​α\displaystyle g_{i\alpha} =t−1g(Ji,Jα)=⟨∇∂xi∂yγ,∂yα⟩vγt−23Rmi​γ​ξ​αvγvξt2+O~(t3).\displaystyle=t^{-1}g(J_{i},J_{\alpha})=\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\gamma}},\partial_{y_{\alpha}}\rangle v_{\gamma}t-\frac{2}{3}\Rm_{i\gamma\xi\alpha}v_{\gamma}v_{\xi}t^{2}+\widetilde{O}(t^{3}).

Let yα=t​vαy_{\alpha}=tv_{\alpha}, then we obtain the desired expansions. ∎

As a digression we briefly discuss the intrinsic meaning of the above expansion. The point is that locally the Riemannian metric gg is approximated by a Riemannian metric gNg_{N} on the normal bundle NN up to the first order. Notice we have the natural projection π:N→P\pi:N\rightarrow P and NN is a Riemannian vector bundle together with an induced “normal” connection, given by the normal component of the Levi-Civita connection of gg. The latter hence gives rise to a distribution of horizontal subspaces at each point of NN, which in our coordinates is spanned by ∂xi−Ai​α​βyβ∂yα\partial_{x_{i}}-A_{i\alpha\beta}y_{\beta}\partial_{y_{\alpha}}, where

(3.33) Ai​α​β(x)≡⟨∇∂xi∂yβ,∂yα⟩|y=0A_{i\alpha\beta}(x)\equiv\langle\nabla_{\partial_{x_{i}}}\partial_{y_{\beta}},\partial_{y_{\alpha}}\rangle|_{y=0}

is a smooth function on UU. We define gNg_{N} so that at each point of NN, the vertical and horizontal subspaces are orthogonal and on the vertical part is given by the bundle metric on NN, and on the horizontal part is given by the perturbation of the base metric gPg_{P} using the second fundamental form. In this way we get a coordinate free description of the above expansion up to the first order.

We also define

(3.34) Ai​j​α​β≡12​(∂xiAj​α​β−∂xjAi​α​β).A_{ij\alpha\beta}\equiv\frac{1}{2}(\partial_{x_{i}}A_{j\alpha\beta}-\partial_{x_{j}}A_{i\alpha\beta}).

Then the curvature of the normal connection is given by

(3.35) Ωi​j​α​β≡Ai​j​α​β+12​(Ai​α​γ​Aj​γ​β−Ai​β​γ​Aj​γ​α).\Omega_{ij\alpha\beta}\equiv A_{ij\alpha\beta}+\frac{1}{2}(A_{i\alpha\gamma}A_{j\gamma\beta}-A_{i\beta\gamma}A_{j\gamma\alpha}).

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}).

∎

3.3. Green’s currents on a cylinder

In this subsection we assume Q≡D×ℝQ\equiv D\times\mathbb{R} is a Riemannian product of a Kähler manifold (D,ωD,JD)(D,\omega_{D},J_{D}) of complex dimension n−1n-1, and the real line ℝ\mathbb{R} with coordinate zz. Given a smooth divisor H⊂DH\subset D, let P≡H×{0}⊂D×{0}P\equiv H\times\{0\}\subset D\times\{0\}. In our discussion sometimes we also naturally identify HH with PP. The results of this subsection will be purely local so DD and HH are not necessarily compact.

The splitting of a line ℝ\mathbb{R} allows us to study the normal exponential map in QQ in terms of the normal exponential map in DD. Notice the normal bundle of PP in QQ is naturally a Riemannian direct sum

(3.253) N=N0⊕ℝz,{N}=N_{0}\oplus\mathbb{R}_{z},

where N0N_{0} is the normal bundle of HH in DD given as the orthogonal complement (T​H)⟂(TH)^{\perp} of T​HTH in T​D|HTD|_{H} (with respect to ωD\omega_{D}). So N0N_{0} is naturally a hermitian line bundle. We also naturally identify N0N_{0} with the holomorphic normal bundle T​D|H/T​HTD|_{H}/TH, as complex line bundles. Therefore, N0N_{0} can be viewed as a holomorphic hermitian line bundle. The normal exponential map of HH in DD is defined by

(3.254) ExpH:N0→D,(p,v)↦Expp⁡(v),\Exp_{H}:N_{0}\to D,\ (p,v)\mapsto\Exp_{p}(v),

which gives a local diffeomorphism from a neighborhood of the zero section in N0N_{0} to a tubular neighborhood of HH in DD. Immediately,

(3.255) d​ExpH:T​N0|H⟶T​D|Hd\text{Exp}_{H}:TN_{0}|_{H}\longrightarrow TD|_{H}

is the identity map under the natural isomorphisms T​N0|H≅N0⊕T​HTN_{0}|_{H}\cong N_{0}\oplus TH and T​D|H≅N0⊕T​HTD|_{H}\cong N_{0}\oplus TH.

Given any point p∈Hp\in H, we may choose local holomorphic coordinates {wi}i=1n−1\{w_{i}\}_{i=1}^{n-1} on DD, centered at pp, such that HH is locally defined by w1=0w_{1}=0. Then d​w1,w2,…,wn−1dw_{1},w_{2},\ldots,w_{n-1} induces local holomorphic coordinates on N0N_{0}, which we denote by {ζ,w2′,…,wn−1′}\{\zeta,w_{2}^{\prime},\ldots,w_{n-1}^{\prime}\}. Given any (p,v)∈N0(p,v)\in N_{0}, its coordinates are by definition given as

(3.256) {ζ=(d​w1)p​(v),wj′=wj​(p),j≥2.\displaystyle\begin{cases}\zeta=(dw_{1})_{p}(v),\\ w_{j}^{\prime}=w_{j}(p),&j\geq 2.\end{cases}

Under the normal exponential map ExpH\Exp_{H}, these coordinates can also be viewed as local (non-holomorphic) coordinates on DD, and when restricted to HH we have wj′=wj​(j≥2)w_{j}^{\prime}=w_{j}(j\geq 2) and d​ζ=d​w1d\zeta=dw_{1}. In particular, {w2′,⋯,wn−1′}\{w_{2}^{\prime},\cdots,w_{n-1}^{\prime}\} still gives holomorphic coordinates on HH.

Similarly using ExpH\Exp_{H}, the coordinate vector field ∂ζ\partial_{\zeta}, originally defined on the normal bundle N0N_{0}, can also be viewed as a local (non-holomorphic) vector field on DD. When restricted to HH, the vector field ∂ζ\partial_{\zeta} can hence be identified with the local section σ\sigma of (T​H)⟂⊂T​D|H(TH)^{\perp}\subset TD|_{H} given by the orthogonal projection of the holomorphic vector field ∂w1\partial_{w_{1}}. Then we obtain a local unitary frame ee of (T​H)⟂(TH)^{\perp} given by

(3.257) e≡σ/|σ|.e\equiv\sigma/|\sigma|.

These generate fiber coordinates y,y¯y,\bar{y} on N0N_{0} such that

(3.258) y=|σ|⋅ζ.y=|\sigma|\cdot\zeta.

In this way we obtain local coordinates {y,y¯,y3=z,w2′,w¯2′,…,wn−1′,w¯n−1′}\{y,\bar{y},y_{3}=z,w_{2}^{\prime},\bar{w}_{2}^{\prime},\ldots,w_{n-1}^{\prime},\bar{w}_{n-1}^{\prime}\} in a neighborhood of pp in QQ. To match with the notation in the previous subsection, with respect to the local orthonormal basis {σ+σ¯,−1(σ−σ¯),∂z}\{\sigma+\bar{\sigma},\sqrt{-1}(\sigma-\bar{\sigma}),\partial_{z}\}, the normal geodesic coordinates are given by {y1=R​e​(y),y2=I​m​(y),y3=z}\{y_{1}=Re(y),y_{2}=Im(y),y_{3}=z\}, and

(3.259) r2=|y|2+z2.r^{2}=|y|^{2}+z^{2}.

Also, the convention for orientation is given such that

(3.260) 2−(n−1)​−1​d​y∧d​y¯∧d​z∧∏j=2n−1(−1​d​wj′∧d​w¯j′)2^{-(n-1)}\sqrt{-1}dy\wedge d\bar{y}\wedge dz\wedge\prod\limits_{j=2}^{n-1}(\sqrt{-1}dw_{j}^{\prime}\wedge d\bar{w}_{j}^{\prime})

defines a positive volume form. In the following, we will also use ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle to denote the hermitian inner product on (1,0)(1,0)-type vectors. The relation with the Riemannian inner product is seen as

(3.261) ⟨ξ,ξ⟩=2​⟨R​e​(ξ),R​e​(ξ)⟩.\langle\xi,\xi\rangle=2\langle Re(\xi),Re(\xi)\rangle.

What the notation means will be clear in the context.

By making PP and QQ smaller we get local existence of Green’s current GPG_{P} for PP in QQ, by Theorem 3.10, with the expansion given there. In our case the formula can be written in terms of the above complex coordinates

Proposition 3.24.

let GPG_{P} be a Green’s current for P≡H×{0}P\equiv H\times\{0\} in QQ, then locally

(3.262) GP=ψ∧d​z+ℛ,G_{P}=\psi\wedge dz+\mathcal{R},

where ℛ∈C∞\mathcal{R}\in C^{\infty} and ψ\psi is family of real-valued (1,1)(1,1)-forms on DD parametrized by zz, satisfying

(3.263) ψ⁡(−z)−ψ⁡(z)∈C∞.\psi(-z)-\psi(z)\in C^{\infty}.

Moreover, in terms of the above local coordinates we can write

(3.264) ψ=−14​r​d​y∧d​y¯+12​r​(y​d​y¯+y¯​d​y)∧Γ+r⋅d​Γ+r−3​Π2(4)+O′​(r2),\psi=\frac{\sqrt{-1}}{4r}dy\wedge d\bar{y}+\frac{1}{2r}(yd\bar{y}+\bar{y}dy)\wedge\Gamma+r\cdot d\Gamma+r^{-3}\Pi_{2}^{(4)}+O^{\prime}(r^{2}),

where Γ\Gamma is a smooth real-valued 11-form locally defined on HH given by

(3.265) Γ(v)≡−−12⟨∇v∂y,∂y⟩|H,v∈TpH,\Gamma(v)\equiv-\frac{\sqrt{-1}}{2}\langle\nabla_{v}\partial_{y},\partial_{y}\rangle\Big|_{H},\ v\in T_{p}H,

and Π2(4)\Pi_{2}^{(4)} is the 22-form given by Notation 3.9 such that it contains at least one of the d​ydy or d​y¯d\bar{y}.

Proof.

This essentially follows from the fact that PP is located on the slice {z=0}\{z=0\} and HH is a complex submanifold of DD. Indeed, we can decompose

(3.266) GP=ψ1∧d​z+ℛ1,G_{P}=\psi_{1}\wedge dz+\mathcal{R}_{1},

where ℛ1\mathcal{R}_{1} does not involve d​zdz. Given any compactly supported test form χ∈Ω02​n−4​(Q)\chi\in\Omega_{0}^{2n-4}(Q), we can write

(3.267) χ=β∧d​z+γ,\chi=\beta\wedge dz+\gamma,

where γ\gamma does not involve d​zdz. Immediately,

(3.268) (ℛ1,Δ​γ)=∫Qℛ1∧Δ​γ=0(\mathcal{R}_{1},\Delta\gamma)=\int_{Q}\mathcal{R}_{1}\wedge\Delta\gamma=0

and

(3.269) (ψ1∧𝑑z,Δ⁡(β∧𝑑z))=∫Qψ1∧𝑑z∧Δ⁡(β∧𝑑z)=0.(\psi_{1}\wedge dz,\Delta(\beta\wedge dz))=\int_{Q}\psi_{1}\wedge dz\wedge\Delta(\beta\wedge dz)=0.

So it follows that

(3.270) (ℛ1,Δ​χ)=(ℛ1,Δ⁡(β∧𝑑z))=(GP,Δ⁡(β∧𝑑z))=2​π​∫P(β∧𝑑z)=0.(\mathcal{R}_{1},\Delta\chi)=(\mathcal{R}_{1},\Delta(\beta\wedge dz))=(G_{P},\Delta(\beta\wedge dz))=2\pi\int_{P}(\beta\wedge dz)=0.

This implies that Δ​ℛ1=0\Delta\mathcal{R}_{1}=0 in the distributional sense. By the standard elliptic regularity, we have ℛ1∈C∞\mathcal{R}_{1}\in C^{\infty}.

Now write

(3.271) ψ1=ψ+ψ2,\psi_{1}=\psi+\psi_{2},

where ψ\psi is JDJ_{D}-invariant, i.e. of type (1,1)(1,1) in DD, and ψ2\psi_{2} is anti-JDJ_{D}-invariant. Since HH is a complex submanifold of DD, the Dirac current δP\delta_{P} is JDJ_{D}-invariant, hence JD​(GP)J_{D}(G_{P}) is also a Green’s current for PP, so we see that ψ2=12​(GP−J⁡(GP))\psi_{2}=\frac{1}{2}(G_{P}-J(G_{P})) is smooth. Then we have

(3.272) GP=ψ∧d​z+ℛ,G_{P}=\psi\wedge dz+\mathcal{R},

where ℛ=ℛ1+ψ2∧d​z\mathcal{R}=\mathcal{R}_{1}+\psi_{2}\wedge dz is smooth. Similarly since the δP\delta_{P} is invariant under z↦−zz\mapsto-z, the difference ψ⁡(z)−ψ⁡(−z)\psi(z)-\psi(-z) is smooth.

To see the expansion of ψ\psi, we notice that HH is a Kähler, in particular minimal, submanifold of DD. So the mean curvature of HH in DD vanishes. Also notice ∂z\partial_{z} is parallel on QQ so Ai​α​β=0A_{i\alpha\beta}=0 if either α=3\alpha=3 or β=3\beta=3. This then implies that

(3.273) ψ=−14​r​d​y∧d​y¯+12​r​(y​d​y¯+y¯​d​y)∧Γ+r⋅𝔸+r−3​Π2(4)+O~​(r2),\psi=\frac{\sqrt{-1}}{4r}dy\wedge d\bar{y}+\frac{1}{2r}(yd\bar{y}+\bar{y}dy)\wedge\Gamma+r\cdot\mathbb{A}+r^{-3}\Pi_{2}^{(4)}+\widetilde{O}(r^{2}),

where

(3.274) Γ\displaystyle\Gamma =−12​Ai​12​d​xi,\displaystyle=-\frac{1}{2}A_{i12}dx_{i},
(3.275) 𝔸\displaystyle\mathbb{A} =−14​Ai​j​α​β​d​xi∧d​xj.\displaystyle=-\frac{1}{4}A_{ij\alpha\beta}dx_{i}\wedge dx_{j}.

In particular, 𝔸=d​Γ\mathbb{A}=d\Gamma. Re-writing

(3.276) Ai​12=⟨∇∂xi∂y2,∂y1⟩A_{i12}=\langle\nabla_{\partial_{x_{i}}}\partial_{y_{2}},\partial_{y_{1}}\rangle

in terms of the complex coordinates y,y¯y,\bar{y} and bearing in mind (3.261) we obtain the desired formula for Γ\Gamma. ∎

Proposition 3.24 has a quick corollary which will be used in our later calculations.

Corollary 3.24.1.

For any positive integer k≥2k\geq 2, we have

(3.277) ψk=O′​(rk−1).\psi^{k}=O^{\prime}(r^{k-1}).
Proof.

By Proposition 3.24, we write

(3.278) ψ\displaystyle\psi =𝔗1+𝔗2+𝔗3,\displaystyle=\FT_{1}+\FT_{2}+\FT_{3},
(3.279) 𝔗1\displaystyle\FT_{1} ≡−14​r​d​y∧d​y¯,\displaystyle\equiv\frac{\sqrt{-1}}{4r}dy\wedge d\bar{y},
(3.280) 𝔗2\displaystyle\FT_{2} ≡12​r​(y​d​y¯+y¯​d​y)∧Γ=O′​(1),\displaystyle\equiv\frac{1}{2r}(yd\bar{y}+\bar{y}dy)\wedge\Gamma=O^{\prime}(1),
(3.281) 𝔗3\displaystyle\FT_{3} ≡r⋅d​Γ+r−3​Π2(4)+O′​(r2)=O′​(r).\displaystyle\equiv r\cdot d\Gamma+r^{-3}\Pi_{2}^{(4)}+O^{\prime}(r^{2})=O^{\prime}(r).

Immediately we have (𝔗1)2=(𝔗2)2=𝔗1∧𝔗2=0(\FT_{1})^{2}=(\FT_{2})^{2}=\FT_{1}\wedge\FT_{2}=0 and for all k≥2k\geq 2, 𝔗2∧(𝔗3)k=(𝔗3)k=O′​(rk)\FT_{2}\wedge(\FT_{3})^{k}=(\FT_{3})^{k}=O^{\prime}(r^{k}). Moreover,

(3.282) 𝔗1∧𝔗3=−14​dy∧d​y¯∧d​Γ+−14​r4​dy∧d​y¯∧Π2(4)+O′​(r).\FT_{1}\wedge\FT_{3}=\frac{\sqrt{-1}}{4}dy\wedge d\bar{y}\wedge d\Gamma+\frac{\sqrt{-1}}{4r^{4}}dy\wedge d\bar{y}\wedge\Pi_{2}^{(4)}+O^{\prime}(r).

Notice that −14​d​y∧d​y¯∧d​Γ\frac{\sqrt{-1}}{4}dy\wedge d\bar{y}\wedge d\Gamma is a smooth term and by definition d​y∧d​y¯∧Π2(4)=0dy\wedge d\bar{y}\wedge\Pi_{2}^{(4)}=0, then

(3.283) 𝔗1∧𝔗3=O′​(r).\FT_{1}\wedge\FT_{3}=O^{\prime}(r).

So for all k≥1k\geq 1,

(3.284) 𝔗1∧(𝔗3)k=O′​(rk).\FT_{1}\wedge(\FT_{3})^{k}=O^{\prime}(r^{k}).

Now by direct calculation,

(3.285) ψk=k⁡((𝔗1)∧(𝔗3)k−1+(𝔗2)∧(𝔗3)k−1)+(𝔗3)k.\displaystyle\psi^{k}=k\Big((\FT_{1})\wedge(\FT_{3})^{k-1}+(\FT_{2})\wedge(\FT_{3})^{k-1}\Big)+(\FT_{3})^{k}.

The conclusion then follows. ∎

Notice that the above local coordinates {y,y¯}\{y,\bar{y}\} are not canonical, and depend on the initial choice of the local coordinates {wi}i=1n−1\{w_{i}\}_{i=1}^{n-1} on DD. However, a different choice of local holomorphic coordinates on DD will induce the coordinates y~,y~¯\tilde{y},\bar{\tilde{y}} on fibers of N0N_{0} such that

(3.286) y=e−1​ϕ⋅y~y=e^{\sqrt{-1}\phi}\cdot\tilde{y}

for some real function ϕ\phi on HH. In particular, we have the transformation

(3.287) d​y∧d​y¯=d​y~∧d​y~¯−−1​d​|y|2∧d​ϕdy\wedge d\bar{y}=d\tilde{y}\wedge d\bar{\tilde{y}}-\sqrt{-1}d|y|^{2}\wedge d\phi

and

(3.288) Γ=Γ~−12​d​ϕ,\Gamma=\widetilde{\Gamma}-\frac{1}{2}d\phi,

This suggests viewing Γ\Gamma as a connection 1-form on the normal bundle. Indeed this is exactly the case.

Lemma 3.25.

2​−1​Γ2\sqrt{-1}\Gamma is the Chern connection 1-form of the normal bundle N0N_{0} with respect to the above hermitian holomorphic structure, in the local holomorphic frame σ\sigma. In other words,

(3.289) Γ=12​dHc​log⁡|σ|.\Gamma=\frac{1}{2}d^{c}_{H}\log|\sigma|.
Proof.

By definition

(3.290) σ=f∂y=∂w1−∑j≥2μj∂wj,\sigma=f\partial_{y}=\partial_{w_{1}}-\sum_{j\geq 2}\mu_{j}\partial_{w_{j}},

where f=|σ|>0f=|\sigma|>0 is local real valued function on HH, and μ2,⋯,μn−1\mu_{2},\cdots,\mu_{n-1} are local complex valued function on HH. The key property we will use is that along HH, ∇∂w¯kσ\nabla_{\partial_{\bar{w}_{k}}}\sigma is tangential to HH for k≥2k\geq 2. In fact, the Kähler condition implies ∇∂¯wk∂wj=0\nabla_{\bar{\partial}_{w_{k}}}\partial_{w_{j}}=0 for all jj, and hence

(3.291) ∇∂w¯kσ=∇∂w¯k(∂w1−∑j=2n−1μj∂wj)=−∑j=2n−1∂w¯k(μj)∂wj.\nabla_{\partial_{\bar{w}_{k}}}\sigma=\nabla_{\partial_{\bar{w}_{k}}}\Big(\partial_{w_{1}}-\sum_{j=2}^{n-1}\mu_{j}\partial_{w_{j}}\Big)=-\sum_{j=2}^{n-1}\partial_{\bar{w}_{k}}(\mu_{j})\partial_{w_{j}}.

Therefore,

(3.292) ∂wkf=∂wk⟨∂y,σ⟩=⟨∇∂wk∂y,f∂y⟩+⟨∂y,∇∂w¯kσ⟩=f⟨∇∂wk∂y,∂y⟩,\partial_{w_{k}}f=\partial_{w_{k}}\langle\partial_{y},\sigma\rangle=\langle\nabla_{\partial_{w_{k}}}\partial_{y},f\partial_{y}\rangle+\langle\partial_{y},\nabla_{\partial_{\bar{w}_{k}}}\sigma\rangle=f\langle\nabla_{\partial_{w_{k}}}\partial_{y},\partial_{y}\rangle,

and hence

(3.293) ⟨∇∂wk∂y,∂y⟩=f−1∂wkf=∂wk(logf).\langle\nabla_{\partial_{w_{k}}}\partial_{y},\partial_{y}\rangle=f^{-1}\partial_{w_{k}}f=\partial_{w_{k}}(\log f).

Differentiating |∂y|2=1|\partial_{y}|^{2}=1, we get

(3.294) ⟨∇∂wk∂y,∂y⟩+⟨∂y,∇∂w¯k∂y⟩=0,\langle\nabla_{\partial_{w_{k}}}\partial_{y},\partial_{y}\rangle+\langle\partial_{y},\nabla_{\partial_{\bar{w}_{k}}}\partial_{y}\rangle=0,

which implies

(3.295) ⟨∇∂w¯k∂y,∂y⟩=−∂w¯klogf.\langle\nabla_{\partial_{\bar{w}_{k}}}\partial_{y},\partial_{y}\rangle=-\partial_{\bar{w}_{k}}\log f.

Therefore,

Γ\displaystyle\Gamma =−−12(∑k≥2⟨∇∂wk∂y,∂y⟩dwk+∑k≥2⟨∇∂w¯k∂y,∂y⟩dw¯k)\displaystyle=-\frac{\sqrt{-1}}{2}(\sum_{k\geq 2}\langle\nabla_{\partial_{w_{k}}}\partial_{y},\partial_{y}\rangle dw_{k}+\sum_{k\geq 2}\langle\nabla_{\partial_{\bar{w}_{k}}}\partial_{y},\partial_{y}\rangle d\bar{w}_{k})
=−−12​(∑k≥2∂wk(log⁡f)​d​wk−∑k≥2∂w¯k(log⁡f)​d​w¯k)\displaystyle=-\frac{\sqrt{-1}}{2}(\sum_{k\geq 2}\partial_{w_{k}}(\log f)dw_{k}-\sum_{k\geq 2}\partial_{\bar{w}_{k}}(\log f)d\bar{w}_{k})
=−−12​(∂Hlog⁡f−∂¯H​log⁡f)\displaystyle=-\frac{\sqrt{-1}}{2}(\partial_{H}\log f-\bar{\partial}_{H}\log f)
(3.296) =12​dHc​log⁡f.\displaystyle=\frac{1}{2}d^{c}_{H}\log f.

∎

For later applications we will need a few more local expansion results. We will also use the notation O′O^{\prime} and O~\widetilde{O} in Definition 3.3. The meaning is similar, but here we work on a neighborhood of HH in DD, and the distance function is locally given by |y||y|. Notice the following expansions are given in the local (non-holomorphic) coordinates {y,y¯,w2′,w¯2′,⋯,wn−1,w¯n−1′}\{y,\bar{y},w_{2}^{\prime},\bar{w}_{2}^{\prime},\cdots,w_{n-1},\bar{w}_{n-1}^{\prime}\}, and by definition we have wj′|H=wj|Hw_{j}^{\prime}|_{H}=w_{j}|_{H} for j≥2j\geq 2.

Proposition 3.26.

The following holds locally near the point p∈Hp\in H,

(3.297) dDc​|y|2=−1​(y​d​y¯−y¯​d​y)+4​|y|2⋅Γ+O~​(|y|3).d^{c}_{D}|y|^{2}=\sqrt{-1}(yd\bar{y}-\bar{y}dy)+4|y|^{2}\cdot\Gamma+\widetilde{O}(|y|^{3}).

The proof relies on the following expansions of the holomorphic coordinate functions wjw_{j}.

Lemma 3.27.

We have the expansion

(3.298) {w1=a1​y+a2​y2+O~​(|y|3)wj=wj′+cj​y+dj​y2+O~​(|y|3),j≥2,\displaystyle\begin{cases}w_{1}=a_{1}y+a_{2}y^{2}+\widetilde{O}(|y|^{3})\\ w_{j}=w_{j}^{\prime}+c_{j}y+d_{j}y^{2}+\widetilde{O}(|y|^{3}),&j\geq 2,\end{cases}

where a1=|σ|−1>0a_{1}=|\sigma|^{-1}>0, a2a_{2}, cjc_{j}, djd_{j} are local smooth functions on HH.

Proof.

By definition, σ\sigma is the orthogonal projection of ∂w1\partial_{w_{1}} onto (T​H)⟂(TH)^{\perp}, so we have along HH,

(3.299) a1−1∂y=∂w1+∑j=2n−1bj∂wj,a_{1}^{-1}\partial_{y}=\partial_{w_{1}}+\sum_{j=2}^{n-1}b_{j}\partial_{w_{j}},

where a1=|σ|−1>0a_{1}=|\sigma|^{-1}>0 and bjb_{j} are smooth functions on HH. Now write

(3.300) ∂y=∑j=1n−1∂wj∂y∂wj+∑j=1n−1∂w¯j∂y∂w¯j.\partial_{y}=\sum_{j=1}^{n-1}\frac{\partial w_{j}}{\partial y}\partial_{w_{j}}+\sum_{j=1}^{n-1}\frac{\partial\bar{w}_{j}}{\partial{y}}\partial_{\bar{w}_{j}}.

then we get that along HH,

(3.301) ∂w¯j∂y=0,j≥1,\displaystyle\frac{\partial\bar{w}_{j}}{\partial y}=0,\ j\geq 1,

which in particular implies

(3.302) ∂wj∂y¯=∂w¯j∂y¯=0,j≥1.\frac{\partial w_{j}}{\partial\bar{y}}=\overline{\frac{\partial\bar{w}_{j}}{\partial y}}=0,\ j\geq 1.

Now by the definition of the normal exponential map, we have at pp,

(3.303) ∇∂y∂y=∇∂y¯∂y=∇∂y¯∂y¯=0.\nabla_{\partial_{y}}\partial_{y}=\nabla_{\partial_{\bar{y}}}\partial_{y}=\nabla_{\partial_{\bar{y}}}\partial_{\bar{y}}=0.

Using the Kähler condition we have

(3.304) ∇∂wj∂w¯k=∇∂w¯j∂wk=0,j,k≥1.\nabla_{\partial_{w_{j}}}{\partial_{\bar{w}_{k}}}=\nabla_{\partial_{\bar{w}_{j}}}\partial_{w_{k}}=0,\ \ j,k\geq 1.

Then by (3.300) we get

(3.305) ∂2wj∂y​∂y¯=∂2wj∂y¯2=0,j≥1.\frac{\partial^{2}w_{j}}{\partial y\partial\bar{y}}=\frac{\partial^{2}w_{j}}{\partial\bar{y}^{2}}=0,\ \ \ j\geq 1.

Therefore, the conclusion follows.

∎

Proof of Proposition 3.26.

Given the above Lemma we first obtain that

(3.306) d​w¯1=a1​d​y¯+y¯​(d​a1+2​a¯2​d​y¯)+O~​(|y|2),d\bar{w}_{1}=a_{1}d\bar{y}+\bar{y}(da_{1}+2\bar{a}_{2}d\bar{y})+\widetilde{O}(|y|^{2}),

then

(3.307) w1​d​w¯1=a12​y​d​y¯+|y|2​a1​d​a1+a1​y​(2​a¯2​y¯+a2​y)​d​y¯+O~​(|y|3).w_{1}d\bar{w}_{1}=a_{1}^{2}yd\bar{y}+|y|^{2}a_{1}da_{1}+a_{1}y(2\bar{a}_{2}\bar{y}+a_{2}y)d\bar{y}+\widetilde{O}(|y|^{3}).

Hence

(3.308) dDc​|w1|2=−1​a12​(y​d​y¯−y¯​d​y)+−1​a1​a¯2​y¯​(2​y​d​y¯−y¯​d​y)−−1​a1​a2​y​(2​y¯​d​y−y​d​y¯)+O~​(|y|3).d_{D}^{c}|w_{1}|^{2}=\sqrt{-1}a_{1}^{2}(yd\bar{y}-\bar{y}dy)+\sqrt{-1}a_{1}\bar{a}_{2}\bar{y}(2yd\bar{y}-\bar{y}dy)-\sqrt{-1}a_{1}a_{2}y(2\bar{y}dy-yd\bar{y})+\widetilde{O}(|y|^{3}).

On the other hand, we have

(3.309) |w1|2=a12​|y|2+a1​(a2​y+a¯2​y¯)​|y|2+O~​(|y|4).|w_{1}|^{2}=a_{1}^{2}|y|^{2}+a_{1}(a_{2}y+\bar{a}_{2}\bar{y})|y|^{2}+\widetilde{O}(|y|^{4}).

So

(3.310) dDc​|w1|2=a12​dDc​|y|2+|y|2​dDc​a12+dDc​(a1​(a2​y+a¯2​y¯)​|y|2)+O~​(|y|3).d_{D}^{c}|w_{1}|^{2}=a_{1}^{2}d_{D}^{c}|y|^{2}+|y|^{2}d_{D}^{c}a_{1}^{2}+d_{D}^{c}(a_{1}(a_{2}y+\bar{a}_{2}\bar{y})|y|^{2})+\widetilde{O}(|y|^{3}).

Now by Lemma 3.27,

(3.311) dDc​(a1​y)=dDc​w1+O~​(|y|)=−−1​d​w1+O~​(|y|),d_{D}^{c}(a_{1}y)=d_{D}^{c}w_{1}+\widetilde{O}(|y|)=-\sqrt{-1}dw_{1}+\widetilde{O}(|y|),

so

(3.312) dDc​y=−−1​d​y+O~​(|y|).d_{D}^{c}y=-\sqrt{-1}dy+\widetilde{O}(|y|).

Similarly, dDc​y¯=−1​d​y¯+O~​(|y|)d_{D}^{c}\bar{y}=\sqrt{-1}d\bar{y}+\widetilde{O}(|y|). Plugging these into (3.310), and compare with (3.308) we obtain

(3.313) dDc​|y|2=−1​(y​d​y¯−y¯​d​y)−2​|y|2​dDc​log⁡a1+O~​(|y|3).\displaystyle d_{D}^{c}|y|^{2}=\sqrt{-1}(yd\bar{y}-\bar{y}dy)-2|y|^{2}d_{D}^{c}\log a_{1}+\widetilde{O}(|y|^{3}).

Thanks to Lemma 3.27, a1=|σ|−1a_{1}=|\sigma|^{-1} which is a smooth function on HH, so

(3.314) dDc​|y|2=−1​(y​d​y¯−y¯​d​y)+2​|y|2​dHc​log⁡|σ|+O~​(|y|3).\displaystyle d_{D}^{c}|y|^{2}=\sqrt{-1}(yd\bar{y}-\bar{y}dy)+2|y|^{2}d_{H}^{c}\log|\sigma|+\widetilde{O}(|y|^{3}).

By Lemma 3.25, Γ=12​dHc​log⁡|σ|\Gamma=\frac{1}{2}d_{H}^{c}\log|\sigma|, so we conclude

(3.315) dDc​|y|2=−1​(y​d​y¯−y¯​d​y)+4​|y|2​Γ+O~​(|y|3).d_{D}^{c}|y|^{2}=\sqrt{-1}(yd\bar{y}-\bar{y}dy)+4|y|^{2}\Gamma+\widetilde{O}(|y|^{3}).

∎

Now we prove an expansion result for the trace of ψ\psi.

Proposition 3.28.

Let ψ\psi be the 22-form on QQ given as in (3.264), then we have the following expansion near PP

(3.316) TrωD⁡ψ=12​r+O′​(r).\Tr_{\omega_{D}}\psi=\frac{1}{2r}+O^{\prime}(r).

Using (3.264) it is easy to see TrωD⁡ψ\Tr_{\omega_{D}}\psi admits an expansion of the form

(3.317) TrωD⁡ψ=A0r+A1​y+A¯1​y¯r+O′​(r).\Tr_{\omega_{D}}\psi=\frac{A_{0}}{r}+\frac{A_{1}y+\bar{A}_{1}\bar{y}}{r}+O^{\prime}(r).

for local functions A0,A1A_{0},A_{1} defined on HH. It suffices to show A0≡1A_{0}\equiv 1 and A1≡0A_{1}\equiv 0. Since the left hand side is independent of the choice of local holomorphic coordinates, it suffices to we only need to work on the slice z=0z=0 with special local holomorphic coordinates in a neighborhood of p∈Hp\in H, and it suffices to understand the Taylor expansion along the fiber N0​(p)N_{0}(p) of N0N_{0} over the fixed point pp.

Lemma 3.29.

We may choose the above holomorphic coordinates {wi}i=1n−1\{w_{i}\}_{i=1}^{n-1} centered at pp, so that HH is given by w1=0w_{1}=0 and

(3.318) ωD=−12​gi​j¯​d​wi∧d​w¯j,\omega_{D}=\frac{\sqrt{-1}}{2}g_{i\bar{j}}dw_{i}\wedge d\bar{w}_{j},

where

(3.319) {gi​j¯​(0)=δi​j,1≤i,j≤n−1,∂w1gi​j¯​(0)=0,1≤i,j≤n−1,∂wkg1​1¯​(0)=∂wkgi​j¯​(0)=0,2≤i,j,k≤n−1.\displaystyle\begin{cases}g_{i\bar{j}}(0)=\delta_{ij},&1\leq i,j\leq n-1,\\ \partial_{w_{1}}g_{i\bar{j}}(0)=0,&1\leq i,j\leq n-1,\\ \partial_{w_{k}}g_{1\bar{1}}(0)=\partial_{w_{k}}g_{i\bar{j}}(0)=0,&2\leq i,j,k\leq n-1.\end{cases}
Remark 3.29.1.

In fact, the only non-trivial Christoffel symbols at pp are

(3.320) Γi​j1​(0)=∂igj​1¯​(0),Γi¯​j¯1¯=∂i¯g1​j¯​(0)\Gamma_{ij}^{1}(0)=\partial_{i}g_{j\bar{1}}(0),\ \ \Gamma_{\bar{i}\bar{j}}^{\bar{1}}=\partial_{\bar{i}}g_{1\bar{j}}(0)

for i,j≥2i,j\geq 2. This is due to the constraint that the equation w1=0w_{1}=0 defines HH, which prevents us from using substitutions like

(3.321) w1=z1+∑i,j=2n−1C1​i​j​zi​zj.w_{1}=z_{1}+\sum\limits_{i,j=2}^{n-1}C_{1ij}z_{i}z_{j}.

Intrinsically, {Γi​j1}i,j≥2\{\Gamma^{1}_{ij}\}_{i,j\geq 2} captures the second fundamental form of the complex hypersurface HH at pp.

Proof of Lemma 3.29.

This follows from elementary manipulation. First, the holomorphic coordinates {wi}i=1n−1\{w_{i}\}_{i=1}^{n-1} can be chosen such that gi​j¯​(0)=δi​jg_{i\bar{j}}(0)=\delta_{ij} for all 1≤i,j≤n−11\leq i,j\leq n-1. By the substitution of the form

(3.322) {wi=zi+12∑j,k=2n−1Ci​j​kzjzk+∑j=2n−1Di​jz1zj+Eiz12, 2≤i≤n−1,w1=z1+∑j=1n−1Fj​z1​zj,\begin{cases}w_{i}=z_{i}+\frac{1}{2}\sum\limits_{j,k=2}^{n-1}C_{ijk}z_{j}z_{k}+\sum\limits_{j=2}^{n-1}D_{ij}z_{1}z_{j}+E_{i}z_{1}^{2},\ \ 2\leq i\leq n-1,\\ w_{1}=z_{1}+\sum\limits_{j=1}^{n-1}F_{j}z_{1}z_{j},\end{cases}

with suitable choices of coefficients, where Ci​j​k=Ci​k​jC_{ijk}=C_{ikj} for 2≤i,j,k≤n−12\leq i,j,k\leq n-1. One can plug both the Taylor expansion of gi​j¯g_{i\bar{j}} along zkz_{k}’s and (3.322) into ωD\omega_{D}. Comparing the coefficients, then it follows that,

(3.323) {Ci​j​k=−∂wkgj​i¯(0),Di​j=−∂wjg1​i¯(0),Ei=−12∂w1g1​j¯(0),Fi=−∂wjg1​1¯(0),F1=−12∂w1g1​1¯(0),\displaystyle\begin{cases}C_{ijk}=-\partial_{w_{k}}g_{j\bar{i}}(0),\\ D_{ij}=-\partial_{w_{j}}g_{1\bar{i}}(0),\\ E_{i}=-\frac{1}{2}\partial_{w_{1}}g_{1\bar{j}}(0),\\ F_{i}=-\partial_{w_{j}}g_{1\bar{1}}(0),\\ F_{1}=-\frac{1}{2}\partial_{w_{1}}g_{1\bar{1}}(0),\end{cases}

where 2≤i,j,k≤n−12\leq i,j,k\leq n-1. Then we can achieve (3.319) with {wi}i=1n−1\{w_{i}\}_{i=1}^{n-1} replaced by {zi}i=1n−1\{z_{i}\}_{i=1}^{n-1}.

∎

Now we prove Proposition 3.28.

Proof of Proposition 3.28.

The goal is to show A0=1A_{0}=1 and A1=0A_{1}=0 in the expansion (3.317). We work in the above special coordinates centered at p∈Hp\in H.

The first step is to show that the O′​(1)O^{\prime}(1)-term in the expansion of ψ\psi given by Proposition 3.24 in fact vanishes along N0​(p)N_{0}(p). To this end, notice that ∂y=σ=∂w1\partial_{y}=\sigma=\partial_{w_{1}} at p∈Hp\in H and hence by Lemma 3.27,

(3.324) w1=y+a2​y2+O~​(|y|3).w_{1}=y+a_{2}y^{2}+\widetilde{O}(|y|^{3}).

Since the only non-trivial Christofell symsbols at pp are Γi​j1\Gamma_{ij}^{1} and Γi¯​j¯1¯\Gamma_{\bar{i}\bar{j}}^{\bar{1}} for i,j≥2i,j\geq 2, it easily follows that

(3.325) d​|σ|​(p)=0.d|\sigma|(p)=0.

Combining (3.325) and Lemma 3.25,

(3.326) Γ⁡(p)=12​(dHc​log⁡|σ|)​(p)=0,\Gamma(p)=\frac{1}{2}(d_{H}^{c}\log|\sigma|)(p)=0,

for each p∈Hp\in H. Therefore, along the fiber N0​(p)N_{0}(p) of the normal bundle N0​(p)N_{0}(p), the expansion of ψ\psi in Proposition 3.24 becomes

(3.327) ψ=−14​|y|​d​y∧d​y¯+O⁡(|y|).\psi=\frac{\sqrt{-1}}{4|y|}dy\wedge d\bar{y}+O(|y|).

Next, we will compute the coefficients A0​(p)A_{0}(p) and A1​(p)A_{1}(p) in (3.317). As in the proof of Lemma 3.27, we obtain that

(3.328) ∂wj∂y​(p)=∂wj∂y¯=0,j≥2,\frac{\partial w_{j}}{\partial y}(p)=\frac{\partial w_{j}}{\partial\bar{y}}=0,\ \ j\geq 2,

and

(3.329) ∂2wj∂y2​(p)=∂2wj∂y​∂y¯​(p)=∂2wj∂y¯2​(p)=0,j≥1.\frac{\partial^{2}w_{j}}{\partial y^{2}}(p)=\frac{\partial^{2}w_{j}}{\partial y\partial\bar{y}}(p)=\frac{\partial^{2}w_{j}}{\partial\bar{y}^{2}}(p)=0,\ \ j\geq 1.

This particularly implies that a2​(p)=0a_{2}(p)=0 and along the fiber N0​(p)N_{0}(p),

(3.330) wj=O⁡(|y|3),j≥2.w_{j}=O(|y|^{3}),\ \ j\geq 2.

By Lemma 3.29, ∂w1gi​j¯​(p)=0\partial_{w_{1}}g_{i\bar{j}}(p)=0 for all 1≤i,j≤n−11\leq i,j\leq n-1, then the expansion of ωD\omega_{D} along the fiber N0​(p)N_{0}(p) is at least quadratic in the w1w_{1}-direction, i.e.

(3.331) ωD\displaystyle\omega_{D} =\displaystyle= −12​(d​w1∧d​w¯1+∑j=2n−1d​wj∧d​w¯j)+O⁡(∑j=2n−1|wj|2)+O⁡(|w1|2)\displaystyle\frac{\sqrt{-1}}{2}\Big(dw_{1}\wedge d\bar{w}_{1}+\sum_{j=2}^{n-1}dw_{j}\wedge d\bar{w}_{j}\Big)+O\Big(\sqrt{\sum_{j=2}^{n-1}|w_{j}|^{2}}\Big)+O(|w_{1}|^{2})
=\displaystyle= −12​(d​w1∧d​w¯1+∑j=2n−1d​wj∧d​w¯j)+O⁡(|y|2).\displaystyle\frac{\sqrt{-1}}{2}(dw_{1}\wedge d\bar{w}_{1}+\sum_{j=2}^{n-1}dw_{j}\wedge d\bar{w}_{j})+O(|y|^{2}).

By (3.324) and (3.330), along the fiber N0​(p)N_{0}(p), we have

(3.332) d​w1=d​y+O⁡(|y|2).dw_{1}=dy+O(|y|^{2}).

and

(3.333) d​wj=d​wj′+O⁡(|y|2),j≥2.dw_{j}=dw_{j}^{\prime}+O(|y|^{2}),\ \ j\geq 2.

So we get

(3.334) ωD=−12​(d​y∧d​y¯+∑j=2n−1d​wj′∧d​w¯j′)+O⁡(|y|2)\omega_{D}=\frac{\sqrt{-1}}{2}(dy\wedge d\bar{y}+\sum_{j=2}^{n-1}dw_{j}^{\prime}\wedge d\bar{w}_{j}^{\prime})+O(|y|^{2})

Since by definition,

(3.335) (TrωD⁡ψ)⋅ωDn−1(n−1)!=ψ∧ωDn−2(n−2)!.\Big(\Tr_{\omega_{D}}\psi\Big)\cdot\frac{\omega_{D}^{n-1}}{(n-1)!}=\psi\wedge\frac{\omega_{D}^{n-2}}{(n-2)!}.

by elementary manipulations we get that A0​(p)=1A_{0}(p)=1 and A1​(p)=0A_{1}(p)=0. ∎

We close this subsection by proving an expansion of a local holomorphic volume form on DD. Given the choice of local holomorphic coordinates on DD as before, let ΩD\Omega_{D} be a local holomorphic volume form in a neighborhood of pp, then we can always write

(3.336) ΩD=f⋅d​w1∧d​w2∧⋯∧d​wn−1,\Omega_{D}=f\cdot dw_{1}\wedge dw_{2}\cdots\wedge dw_{n-1},

for a local nowhere vanishing holomorphic function ff. Denote the local holomorphic volume form on HH

(3.337) ΩH≡dw2′∧⋯dwn−1′=(dw2∧⋯∧dwn−1)|H\Omega_{H}\equiv dw_{2}^{\prime}\wedge\cdots dw_{n-1}^{\prime}=(dw_{2}\wedge\cdots\wedge dw_{n-1})|_{H}

Then ΩH\Omega_{H} can be naturally viewed as a complex (n−2)(n-2)-form in some neighborhood of pp in DD, in the coordinate system given by {y,y¯,w2′,w¯2′,⋯,wn−1′,w¯n−1′}\{y,\bar{y},w_{2}^{\prime},\bar{w}_{2}^{\prime},\cdots,w_{n-1}^{\prime},\bar{w}_{n-1}^{\prime}\}.

Proposition 3.30.

We have the following expansion

(3.338) ΩD=F⁡(d​y+2​−1​y​Γ)∧ΩH+O~​(|y|)​d​y+O~​(|y|2)\Omega_{D}=F(dy+2\sqrt{-1}y\Gamma)\wedge\Omega_{H}+\widetilde{O}(|y|)dy+\widetilde{O}(|y|^{2})

for some local smooth function FF on HH.

Proof.

We need to calculate the expansion for d​w1∧…∧d​wn−1dw_{1}\wedge\ldots\wedge dw_{n-1}. First, by Lemma 3.27,

(3.339) d​w1=a1​(d​y+y​dH​log⁡a1)+O~​(|y|)​d​y+O~​(|y|2),dw_{1}=a_{1}(dy+yd_{H}\log a_{1})+\widetilde{O}(|y|)dy+\widetilde{O}(|y|^{2}),

where a1=|σ|−1a_{1}=|\sigma|^{-1}. Notice that

(3.340) dH​log⁡a1=2​∂Hlog⁡a1−−1​dHc​log⁡a1.d_{H}\log a_{1}=2\partial_{H}\log a_{1}-\sqrt{-1}d_{H}^{c}\log a_{1}.

Applying Lemma 3.25,

(3.341) dH​log⁡a1=2​(∂Hlog⁡a1+−1​Γ).d_{H}\log a_{1}=2(\partial_{H}\log a_{1}+\sqrt{-1}\Gamma).

Next, applying Lemma 3.27 to wjw_{j}’s for j≥2j\geq 2,

(3.342) d​wj=d​wj′+cj​d​y+y​d​cj+O~​(|y|)​d​y+O~​(|y|2).dw_{j}=dw_{j}^{\prime}+c_{j}dy+ydc_{j}+\widetilde{O}(|y|)dy+\widetilde{O}(|y|^{2}).

Since it holds that

(3.343) ∂Hlog⁡a1∧ΩH≡0,\partial_{H}\log a_{1}\wedge\Omega_{H}\equiv 0,

then taking the wedge product,

(3.344) d​w1∧⋯∧d​wn−1=a1​(d​y+2​−1​y​Γ)∧ΩH+O~​(|y|)​d​y+O~​(|y|2).dw_{1}\wedge\cdots\wedge dw_{n-1}=a_{1}(dy+2\sqrt{-1}y\Gamma)\wedge\Omega_{H}+\widetilde{O}(|y|)dy+\widetilde{O}(|y|^{2}).

On the other hand, we have the expansion of ff,

(3.345) f=f|H+∂f∂y|H⋅y+∂f∂y¯|H⋅y¯+O~​(|y|2).f=f|_{H}+\frac{\partial f}{\partial y}|_{H}\cdot y+\frac{\partial f}{\partial\bar{y}}|_{H}\cdot\bar{y}+\widetilde{O}(|y|^{2}).

Therefore,

(3.346) ΩD=f|H​a1​(d​y+2​−1​y​Γ)∧ΩH+O~​(|y|)​d​y+O~​(|y|2).\Omega_{D}=f|_{H}a_{1}(dy+2\sqrt{-1}y\Gamma)\wedge\Omega_{H}+\widetilde{O}(|y|)dy+\widetilde{O}(|y|^{2}).

So we obtain the conclusion by taking F=f|H⋅a1F=f|_{H}\cdot a_{1}.

∎

3.4. A global existence result

In this subsection, we will prove a global existence result for Green’s currents. Although the results hold in general Riemannian settings, for our application we shall only state the result in a special setting. We assume now (D,ωD,JD)(D,\omega_{D},J_{D}) is a compact Kähler manifold, HH is a smooth divisor Poincaré dual to k2​π​[ωD]\frac{k}{2\pi}[\omega_{D}] for some positive k∈ℤ+k\in\mathbb{Z}_{+}.

Proposition 3.31.

In the above context, given any constants k−,k+∈ℝk_{-},k_{+}\in\mathbb{R} with

(3.347) k−−k+=k,k_{-}-k_{+}=k,

there exists a unique global Green’s current GPG_{P} for PP in QQ such that the following properties hold:

  1. (1)

    GPG_{P} is of the form

    (3.348) GP=ψ⁡(z)∧d​z.G_{P}=\psi(z)\wedge dz.

    Moreover, for each z∈ℝz\in\mathbb{R}, ψ⁡(z)\psi(z) is a closed real (1,1)(1,1)-current on DD.

  2. (2)

    For any nonnegative integer k∈ℕk\in\mathbb{N} and for any δ∈(0,10−2)\delta\in(0,10^{-2}),

    (3.349) {|∇k(ψ⁡(z)−(k−​z)⋅ωD)|=O⁡(e(1−δ)​λ1​z),z→−∞,|∇k(ψ⁡(z)−(k+​z)⋅ωD)|=O⁡(e−(1−δ)​λ1​z),z→∞,\displaystyle\begin{cases}|\nabla^{k}(\psi(z)-(k_{-}z)\cdot\omega_{D})|=O(e^{(1-\delta)\sqrt{\lambda_{1}}z}),&z\rightarrow-\infty,\\ |\nabla^{k}(\psi(z)-(k_{+}z)\cdot\omega_{D})|=O(e^{-(1-\delta)\sqrt{\lambda_{1}}z}),&z\rightarrow\infty,\end{cases}

    where λ1>0\lambda_{1}>0 is the first eigenvalue of the Hodge Laplacian acting on closed real (1,1)(1,1)-forms on DD.

Remark 3.31.1.

This proposition can be seen as a generalization of theorem 2.6 in [HSVZ18] whose proof uses the general existence result of Green’s function on 𝕋2×ℝ\mathbb{T}^{2}\times\mathbb{R}. We thank Lorenzo Foscolo for discussions concerning the following proof via Fourier expansion, which is more constructive.

The proof of the proposition relies on the spectral analysis of the Hodge Laplacian. In particular, we need the following CkC^{k}-estimate of the eigenforms in terms of the eigenvalues. This lemma will be also used in Section 5. The proof follows from standard W2,pW^{2,p}-elliptic regularity and the Sobolev embedding theorems, so we omit it.

Lemma 3.32.

Let (Mm,g)(M^{m},g) be a closed Riemannian manifold of dimension m≥2m\geq 2. For any p∈ℕp\in\mathbb{N}, denote by Λ(p)≡{λj}j=0∞\Lambda^{(p)}\equiv\{\lambda_{j}\}_{j=0}^{\infty} with λ0=0\lambda_{0}=0 the spectrum of the Hodge Laplacian Δ\Delta acting on the pp-forms. For any k∈ℕk\in\mathbb{N}, there is some constant C>0C>0 depending only on (M,g)(M,g) and kk, pp such that for all ϕj∈Ωp​(Mm)\phi_{j}\in\Omega^{p}(M^{m}) satisfying

(3.350) {Δ​ϕj=λj​ϕj,‖ϕj‖L2​(Mm)=1,\displaystyle\begin{cases}\Delta\phi_{j}=\lambda_{j}\phi_{j},\\ \|\phi_{j}\|_{L^{2}(M^{m})}=1,\end{cases}

we have

(3.351) ‖∇kϕj‖Ck​(Mm)≤C⋅(λj)12​[m2]+k+12.\|\nabla^{k}\phi_{j}\|_{C^{k}(M^{m})}\leq C\cdot(\lambda_{j})^{\frac{1}{2}[\frac{m}{2}]+\frac{k+1}{2}}.

Now we are ready to prove Proposition 3.31.

Proof of Proposition 3.31.

The uniqueness follows from the fact that the difference of any two Green’s currents for PP differ by a harmonic form, which must vanish by the asymptotic condition (3.349).

Now we focus on the proof of the global existence of GPG_{P} on QQ. Let {ϕj}j=0∞\{\phi_{j}\}_{j=0}^{\infty} be a complete orthonormal basis of eigenvectors for the Hodge Laplacian acting on real-valued 22-forms on DD, and let ΛD≡{λj}j=0∞\Lambda_{D}\equiv\{\lambda_{j}\}_{j=0}^{\infty} be the corresponding spectrum. Our basic strategy is to first obtain a formal series expression of GPG_{P} and then prove the convergence of this series.

To begin with, the Dirac 33-current δP\delta_{P} of P⊂QP\subset Q has a formal expansion along the DD direction

(3.352) δP=∑j=0∞fj​(z)​ϕj∧d​z,\delta_{P}=\sum_{j=0}^{\infty}f_{j}(z)\phi_{j}\wedge dz,

where fj​(z)f_{j}(z) is a 00-current on ℝ\mathbb{R} and given by

(3.353) fj(z)≡2π(∫H∗Dϕj)δ0(z),f_{j}(z)\equiv 2\pi\Big(\int_{H}*_{D}\phi_{j}\Big)\delta_{0}(z),

where δ0​(z)\delta_{0}(z) is the standard Dirac 00-current acting on functions on ℝ\mathbb{R}, supported at the slice {z=0}\{z=0\}. For λj>0\lambda_{j}>0, by Hodge theory, fjf_{j} is non-zero only when ϕj\phi_{j} is a closed real (1,1)(1,1)-form because HH is a closed complex submanifold in DD. So we only restrict to the subset of such λj\lambda_{j}’s. Furthermore, if ϕj\phi_{j} is harmonic, then

(3.354) ∫H∗Dϕj=∫Dd2​πωD∧∗Dϕj=d2​π⟨ωD,ϕj⟩L2​(D).\int_{H}*_{D}\phi_{j}=\int_{D}\frac{d}{2\pi}\omega_{D}\wedge*_{D}\phi_{j}=\frac{d}{2\pi}\langle\omega_{D},\phi_{j}\rangle_{L^{2}(D)}.

It follows that there is exactly one jj, which we may assume to be 00, such that λj=0\lambda_{j}=0 and fjf_{j} is non-zero. The corresponding eigenform is normalized to be

(3.355) ϕ0≡1(∫DωDn/n!)1/2⋅ωD.\phi_{0}\equiv\frac{1}{(\int_{D}\omega_{D}^{n}/n!)^{1/2}}\cdot\omega_{D}.

Now let GPG_{P} be the formal series

(3.356) GP=∑j=0∞hj⋅(ϕj∧d​z),G_{P}=\sum\limits_{j=0}^{\infty}h_{j}\cdot(\phi_{j}\wedge dz),

where hjh_{j} satisfies

(3.357) −hj′′​(z)+λj⋅hj​(z)=fj​(z).-h_{j}^{\prime\prime}(z)+\lambda_{j}\cdot h_{j}(z)=f_{j}(z).

For each j∈ℤ+j\in\mathbb{Z}_{+}, we can write a formal solution

hj​(z)\displaystyle h_{j}(z) =π−λj(e−λj⋅z∫−∞zeλj⋅ufj(u)du+eλj​z∫z∞e−λj⋅ufj(u)du)\displaystyle=\frac{\pi}{-\sqrt{\lambda_{j}}}\Big(e^{-\sqrt{\lambda_{j}}\cdot z}\int_{-\infty}^{z}e^{\sqrt{\lambda_{j}}\cdot u}f_{j}(u)du+e^{\sqrt{\lambda_{j}}z}\int_{z}^{\infty}e^{-\sqrt{\lambda_{j}}\cdot u}f_{j}(u)du\Big)
(3.358) ={π−λj⋅e−λj⋅z⋅∫H∗D(ϕj),z>0,π−λj⋅eλj⋅z⋅∫H∗D(ϕj),z≤0.\displaystyle=\begin{cases}\frac{\pi}{-\sqrt{\lambda_{j}}}\cdot e^{-\sqrt{\lambda_{j}}\cdot z}\cdot\int_{H}*_{D}(\phi_{j}),&z>0,\\ \frac{\pi}{-\sqrt{\lambda_{j}}}\cdot e^{\sqrt{\lambda_{j}}\cdot z}\cdot\int_{H}*_{D}(\phi_{j}),&z\leq 0.\end{cases}

For j=0j=0, a solution is given by a piecewise linear function

(3.359) h0​(z)={k+​z,z≥0,k−​z,z≤0.\displaystyle h_{0}(z)=\begin{cases}k_{+}z,&z\geq 0,\\ k_{-}z,&z\leq 0.\end{cases}

Notice that the formal solution h0​(z)h_{0}(z) is unique up to the addition of a linear function in zz. Fixing a choice of h0h_{0} we then obtain a formal solution GPG_{P}.

Next we show that the above formal series GPG_{P} is well-defined by showing the formal solution indeed converges in the weak sense and has some exponential decaying rate as |z||z| large, which consists of two steps.

In the first step, we claim that globally the formal expansion

(3.360) GP≡∑j=0∞hj⋅(ϕj∧d​z)G_{P}\equiv\sum\limits_{j=0}^{\infty}h_{j}\cdot(\phi_{j}\wedge dz)

in fact gives a well-defined 33-current on QQ and the series converges in the following sense: for any test form χ∈Ω0m−3​(Q)\chi\in\Omega_{0}^{m-3}(Q),

(3.361) ∑j=0N(hj​(z)​ϕj∧d​z,χ)⟶(GP,χ)​as​N→∞.\displaystyle\sum\limits_{j=0}^{N}\Big(h_{j}(z)\phi_{j}\wedge dz,\chi\Big)\longrightarrow(G_{P},\chi)\ \text{as}\ N\to\infty.

It suffices to show that for any smooth test form χ∈Ω0m−3​(Q)\chi\in\Omega_{0}^{m-3}(Q) and for any N∈ℕN\in\mathbb{N},

(3.362) ∑j=0N|(hj​(z)​ϕj∧d​z,χ)|≤C0,\sum\limits_{j=0}^{N}\Big|\Big(h_{j}(z)\phi_{j}\wedge dz,\chi\Big)\Big|\leq C_{0},

where C0>0C_{0}>0 is independent of NN. To see this, for each jj, we write

(3.363) (hj​(z)​ϕj∧d​z,χ)\displaystyle\Big(h_{j}(z)\phi_{j}\wedge dz,\chi\Big) =\displaystyle= ∫Qhj​(z)​ϕj∧𝑑z∧χ\displaystyle\int_{Q}h_{j}(z)\phi_{j}\wedge dz\wedge\chi
=\displaystyle= ∫ℝ(hj(z)⋅∫D⟨χ,∗D(ϕj)⟩)dz.\displaystyle\int_{\mathbb{R}}\Big(h_{j}(z)\cdot\int_{D}\langle\chi,*_{D}(\phi_{j})\rangle\Big)dz.

The estimate (3.363) can be accomplished in the following manner. To begin with, we will show that the integral ∫D⟨χ,∗D(ϕj)⟩\int_{D}\langle\chi,*_{D}(\phi_{j})\rangle has an uniform bound which is independent of zz. In fact, notice that (Δ)k​ϕj=(λj)k​ϕj(\Delta)^{k}\phi_{j}=(\lambda_{j})^{k}\phi_{j} holds for any k∈ℤ+k\in\mathbb{Z}_{+}, then

(3.364) ∫D⟨χ(z),∗D(ϕj)⟩dvolg\displaystyle\int_{D}\langle\chi(z),*_{D}(\phi_{j})\rangle\dvol_{g} =\displaystyle= 1(λj)k∫D⟨χ(z),∗D(Δ)k(ϕj)⟩dvolg\displaystyle\frac{1}{(\lambda_{j})^{k}}\int_{D}\langle\chi(z),*_{D}(\Delta)^{k}(\phi_{j})\rangle\dvol_{g}
=\displaystyle= 1(λj)k​∫D⟨χ⁡(z),(Δ)k∗D(ϕj)⟩​dvolg\displaystyle\frac{1}{(\lambda_{j})^{k}}\int_{D}\langle\chi(z),(\Delta)^{k}*_{D}(\phi_{j})\rangle\dvol_{g}
=\displaystyle= 1(λj)k∫D⟨(Δ)kχ(z),∗D(ϕj)⟩.\displaystyle\frac{1}{(\lambda_{j})^{k}}\int_{D}\langle(\Delta)^{k}\chi(z),*_{D}(\phi_{j})\rangle.

Lemma 3.32 implies

(3.365) ‖ϕj‖C0​(D)≤C⋅(λj)n2,\|\phi_{j}\|_{C^{0}(D)}\leq C\cdot(\lambda_{j})^{\frac{n}{2}},

where CC depends only on nn and the metric gg. So it follows that

(3.366) |∫D⟨χ(z),∗D(ϕj)⟩dvolg|≤C⋅∥χ∥C2​k​(Q)⋅1(λj)k−n2.\Big|\int_{D}\langle\chi(z),*_{D}(\phi_{j})\rangle\dvol_{g}\Big|\leq C\cdot\|\chi\|_{C^{2k}(Q)}\cdot\frac{1}{(\lambda_{j})^{k-\frac{n}{2}}}.

Next, we will estimate the integral ∫ℝhj​(z)​𝑑z\int_{\mathbb{R}}h_{j}(z)dz. To this end, for each j∈ℤ+j\in\mathbb{Z}_{+}, let ϕj\phi_{j} satsify

(3.367) {Δ​ϕj=λj⋅ϕj‖ϕj‖L2​(D)=1.\displaystyle\begin{cases}\Delta\phi_{j}=\lambda_{j}\cdot\phi_{j}\\ \|\phi_{j}\|_{L^{2}(D)}=1.\end{cases}

By Lemma 3.32, for each k∈ℕk\in\mathbb{N},

(3.368) |∇ωDkϕj|≤C​(λj)k+n2.|\nabla^{k}_{\omega_{D}}\phi_{j}|\leq C(\lambda_{j})^{\frac{k+n}{2}}.

For fixed constant z0>102z_{0}>10^{2}, applying (3.358) and (3.368),

(3.369) |∫ℝhj​(z)​𝑑z|\displaystyle\Big|\int_{\mathbb{R}}h_{j}(z)dz\Big| ≤\displaystyle\leq ∫−∞−z0|hj​(z)|​𝑑z+∫−z0z0|hj​(z)|​𝑑z+∫z0+∞|hj​(z)|​𝑑z\displaystyle\int_{-\infty}^{-z_{0}}|h_{j}(z)|dz+\int_{-z_{0}}^{z_{0}}|h_{j}(z)|dz+\int_{z_{0}}^{+\infty}|h_{j}(z)|dz
≤\displaystyle\leq C⁡(1+λjn2−1).\displaystyle C(1+\lambda_{j}^{\frac{n}{2}-1}).

Combining the above estimates, we have

(3.370) |(hj(z)ϕj∧dz,χ)|=|∫ℝ(hj(z)⋅∫D⟨χ,∗D(ϕj)⟩)dz|≤C⋅∥χ∥C2​k​(Q)⋅C⁡(1+λjn2−1)(λj)k−n2.|(h_{j}(z)\phi_{j}\wedge dz,\chi)|=\Big|\int_{\mathbb{R}}\Big(h_{j}(z)\cdot\int_{D}\langle\chi,*_{D}(\phi_{j})\rangle\Big)dz\Big|\leq C\cdot\|\chi\|_{C^{2k}(Q)}\cdot\frac{C(1+\lambda_{j}^{\frac{n}{2}-1})}{(\lambda_{j})^{k-\frac{n}{2}}}.

Then applying Weyl’s law, if kk is sufficiently large, then the above series converges as stated in (3.362), which completes the proof of the claim.

At our next stage, we will study the exponential decaying behavior of the current GPG_{P} defined in (3.360). For any j∈ℤ+j\in\mathbb{Z}_{+} and for any number z>10n2+k2z>10^{n^{2}+k^{2}}, we have

(3.371) |∇Qk(hj​(z)⋅ϕj)|≤C​(λj)2​n+k−12​e−λj​z.|\nabla^{k}_{Q}(h_{j}(z)\cdot\phi_{j})|\leq C(\lambda_{j})^{\frac{2n+k-1}{2}}e^{-\sqrt{\lambda_{j}}z}.

Notice that by elementary computations, for each δ∈(0,10−2)\delta\in(0,10^{-2}), there is some z0>10n2+k2z_{0}>10^{n^{2}+k^{2}} such that for all z∈(z0,+∞)z\in(z_{0},+\infty) and j∈ℤ+j\in\mathbb{Z}_{+},

(3.372) (λj)2​n+k−12⋅e−δ⋅λjz≤(λj)−6​n.(\lambda_{j})^{\frac{2n+k-1}{2}}\cdot e^{-\delta\cdot\sqrt{\lambda_{j}}z}\leq(\lambda_{j})^{-6n}.

This implies that

(3.373) ∑j=0∞|∇Qk(hj​(z)⋅ϕj)|≤C​e−(1−δ)​λ1​z⋅∑j=0∞(λj)−6​n.\sum\limits_{j=0}^{\infty}|\nabla^{k}_{Q}(h_{j}(z)\cdot\phi_{j})|\leq Ce^{-(1-\delta)\sqrt{\lambda_{1}}z}\cdot\sum\limits_{j=0}^{\infty}(\lambda_{j})^{-6n}.

By Weyl’s law implies that the above numerical series converges, and hence for each k∈ℕk\in\mathbb{N} ∇Qk(GP−h0​(z))\nabla_{Q}^{k}(G_{P}-h_{0}(z)) has an exponential decaying rate as z→+∞z\rightarrow+\infty. The argument is identical for z<0z<0.

The only remaining part is to show that the series GPG_{P} defined by (3.360) satisfies the current equation

(3.374) Δ​GP=2​π​δP\Delta G_{P}=2\pi\delta_{P}

in the distributional sense, i.e., for any χ∈Ω0m−3​(Q)\chi\in\Omega_{0}^{m-3}(Q),

(3.375) (GP,Δ​χ)=2​π​∫Pχ.(G_{P},\Delta\chi)=2\pi\int_{P}\chi.

Applying the definition of fjf_{j}, hjh_{j} and integration by parts, it is straightforward that for each j∈ℕj\in\mathbb{N},

(3.376) (hj​(z)​ϕj∧d​z,Δ​χ)=(fj​(z)​ϕj∧d​z,χ).(h_{j}(z)\phi_{j}\wedge dz,\Delta\chi)=(f_{j}(z)\phi_{j}\wedge dz,\chi).

Since χ∈Ω0m−3​(Q)\chi\in\Omega_{0}^{m-3}(Q), the smooth 22-form ∗D(χ)*_{D}(\chi) has the following L2L^{2}-expansion on the slice D×{0}D\times\{0\},

(3.377) ∗D(χ(0))=∑j=0∞⟨∗D(χ(0)),ϕj⟩D⋅ϕj=∑j=0∞⟨χ(0),∗D(ϕj)⟩D⋅ϕj*_{D}(\chi(0))=\sum\limits_{j=0}^{\infty}\langle*_{D}(\chi(0)),\phi_{j}\rangle_{D}\cdot\phi_{j}=\sum\limits_{j=0}^{\infty}\langle\chi(0),*_{D}(\phi_{j})\rangle_{D}\cdot\phi_{j}

and hence

(3.378) χ⁡(0)=∑j=0∞(∫Dϕj∧χ⁡(0))∗D(ϕj).\chi(0)=\sum\limits_{j=0}^{\infty}\Big(\int_{D}\phi_{j}\wedge\chi(0)\Big)*_{D}(\phi_{j}).

This implies that

(3.379) (∑j=0∞fj⋅ϕj∧dz,χ)=2π∑j=0∞(∫P∗D(ϕj))⋅∫Dϕj∧χ(0)=∫Pχ(0)=2π∫Pχ.\Big(\sum\limits_{j=0}^{\infty}f_{j}\cdot\phi_{j}\wedge dz,\chi\Big)=2\pi\sum\limits_{j=0}^{\infty}\Big(\int_{P}*_{D}(\phi_{j})\Big)\cdot\int_{D}\phi_{j}\wedge\chi(0)=\int_{P}\chi(0)=2\pi\int_{P}\chi.

Therefore,

(3.380) (GP,Δ​χ)=(∑j=0∞hj⋅ϕj∧𝑑z,Δ​χ)=(∑j=0∞fj⋅ϕj∧𝑑z,χ)=2​π​∫Pχ=2​π​δP​(χ),(G_{P},\Delta\chi)=\Big(\sum\limits_{j=0}^{\infty}h_{j}\cdot\phi_{j}\wedge dz,\Delta\chi\Big)=\Big(\sum\limits_{j=0}^{\infty}f_{j}\cdot\phi_{j}\wedge dz,\chi\Big)=2\pi\int_{P}\chi=2\pi\delta_{P}(\chi),

which completes the proof. ∎

The constants k−k_{-} and k+k_{+} determines some information of the above ψ\psi.

Lemma 3.33.

Let ψ\psi be the (1,1)(1,1)-current in Proposition 3.31, then the following holds:

  1. (1)

    The cohomology class [∂zψ⁡(z)]∈H2​(D,ℝ)[\partial_{z}\psi(z)]\in H^{2}(D;\mathbb{R}) is given by k−​[ωD]k_{-}[\omega_{D}] and k+​[ωD]k_{+}[\omega_{D}] for z<0z<0 and z>0z>0 respectively.

  2. (2)

    At z=0z=0, we have

    (3.381) ∂zψ⁡(0)=12​(k−+k+)​ωD.\partial_{z}\psi(0)=\frac{1}{2}(k_{-}+k_{+})\omega_{D}.

    In particular, it extends smoothly across PP.

Proof.

First, we prove Item (1). Since QQ is a Riemannian product, we have for z≠0z\neq 0,

(3.382) d2d​z2​ψ​(z)=ΔD​ψ​(z)=d​d∗​ψ​(z)\frac{d^{2}}{dz^{2}}\psi(z)=\Delta_{D}\psi(z)=dd^{*}\psi(z)

is exact, which implies that the cohomology class [ψ⁡(z)]∈H2​(D,ℝ)[\psi(z)]\in H^{2}(D;\mathbb{R}) is locally constant for z∈ℝ∖{0}z\in\mathbb{R}\setminus\{0\}. On the other hand, by the exponential decay property in (3.349) we see that

(3.383) limz→±∞[ψ⁡(z)]=limz→∞k±​[ωD].\lim_{z\rightarrow\pm\infty}[\psi(z)]=\lim_{z\rightarrow\infty}k_{\pm}[\omega_{D}].

For Item (2), denote

(3.384) ψ~​(z)≡ψ⁡(−z)+(k−+k+)​z​ωD.\tilde{\psi}(z)\equiv\psi(-z)+(k_{-}+k_{+})z\omega_{D}.

Then ψ~∧d​z\tilde{\psi}\wedge dz is also a Green current for PP and it is also asymptotic to k±​z⋅ωDk_{\pm}z\cdot\omega_{D} as z→±∞z\rightarrow\pm\infty. Therefore by uniqueness, ψ~​(z)=ψ​(z)\tilde{\psi}(z)=\psi(z). Taking the zz-derivative at z=0z=0 we get the conclusion. ∎

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