Doubling a boundary problem and computing its index

A first-order elliptic boundary problem can be turned into an operator on a closed doubled manifold. The construction has three parts that must agree exactly: the reflected half must carry the complementary stable modes, the boundary coupling must stay elliptic while separate boundary conditions become matching traces, and the resulting mixed operator must be approximated in operator norm before its continuous symbol can compute the index.

This lesson proves all three parts. It also resolves a possible factor-of-two ambiguity: the geometric double built here has a zero-index complementary half, so its index equals the boundary index. A one-half formula belongs to a different doubled datum whose complementary half has the same index as the original problem.

The named prerequisites are Stable modes and the algebra of boundary data, Fredholm boundary problems with first-order Calderón defects, Finite defects under perturbation, and Symbols, finite defects, and the index on a closed manifold. We use Dt=−i∂tD_t=-i\partial_t, inward collar coordinates, and Hermitian inner products linear in the first entry. Matrix factors retain their displayed order.

1. The split first-order model

Let XX be compact with boundary YY. Suppose

E|Y=E+⊕E−(DI1) E|_Y=E^+\oplus E^- \tag{DI1}

is orthogonal for a chosen Hermitian metric. Choose a product density and product metric in a collar Y×[0,δ)tY\times[0,\delta)_t. Let ϕ∈Cc∞([0,δ))\phi\in C_c^\infty([0,\delta)) be real, 0≤ϕ≤10\leq\phi\leq1, and ϕ=1\phi=1 near t=0t=0. Take Λ±∈Ψphg1(Y;E±)\Lambda^\pm\in\Psi_{\mathrm{phg}}^1(Y;E^\pm) with positive scalar real principal symbols λ±(y,η)IE±\lambda^\pm(y,\eta)I_{E^\pm}, and take an interior operator Λ∈Ψphg1(X∘;E)\Lambda\in\Psi_{\mathrm{phg}}^1(X^\circ;E) whose principal symbol λ(x,ξ)IE\lambda(x,\xi)I_E is positive for every ξ≠0\xi\ne0. Extend Λ±\Lambda^\pm constantly in tt where ϕ=1\phi=1, and put

Pbu=ϕ(Dt+iΛ+00−Dt+iΛ−)(u+u−),Piu=i(1−ϕ)Λ(1−ϕ)u,P=Pb+Pi,Bu=γ0u−.(DI2) \begin{aligned} P^bu&=\phi \begin{pmatrix} D_t+i\Lambda^+&0\\ 0&-D_t+i\Lambda^- \end{pmatrix} \binom{u^+}{u^-},\\ P^iu&=i(1-\phi)\Lambda(1-\phi)u,\qquad P=P^b+P^i,\qquad Bu=\gamma_0u^-. \end{aligned} \tag{DI2}

In the collar-transition region the principal symbol is

p(y,t,η,τ)=ϕ(t)(τ+iλ+(y,η)00−τ+iλ−(y,η))+i(1−ϕ(t))2λ(y,t,η,τ)IE.(DI3) p(y,t,\eta,\tau)= \phi(t) \begin{pmatrix} \tau+i\lambda^+(y,\eta)&0\\ 0&-\tau+i\lambda^-(y,\eta) \end{pmatrix} +i(1-\phi(t))^2\lambda(y,t,\eta,\tau)I_E. \tag{DI3}

If 0<ϕ<10<\phi<1, the imaginary part of each diagonal scalar is positive at every nonzero covector. Where ϕ=1\phi=1, a nonzero tangential covector gives a positive imaginary part and a nonzero pure normal covector gives the real entries ±τ\pm\tau. Where ϕ=0\phi=0, the last term is invertible. Thus PP is elliptic.

Freeze at (y,η)∈T*Y\0(y,\eta)\in T^*Y\setminus0. The two normal equations and their solutions are

(Dt+iλ+)v+=0,v+(t)=eλ+tv+(0),v+ bounded⇔v+(0)=0,(−Dt+iλ−)v−=0,v−(t)=e−λ−tv−(0),v− bounded for every v−(0).(DI4) \begin{array}{lll} (D_t+i\lambda^+)v^+=0, &v^+(t)=e^{\lambda^+t}v^+(0), &v^+\text{ bounded}\Longleftrightarrow v^+(0)=0,\\[1mm] (-D_t+i\lambda^-)v^-=0, &v^-(t)=e^{-\lambda^-t}v^-(0), &v^-\text{ bounded for every }v^-(0). \end{array} \tag{DI4}

Hence the stable Cauchy space is exactly Ey−E^-_y, and the principal boundary map v↦v−(0)v\mapsto v^-(0) is its identity. The generalized interior and boundary symbols are elliptic. Fredholm boundary problems with first-order Calderón defects therefore makes

(P,B)s:H‾s(X∘;E)→H‾s−1(X∘;E)⊕Hs−1/2(Y;E−)(DI5) (P,B)_s:\bar H_s(X^\circ;E)\longrightarrow \bar H_{s-1}(X^\circ;E)\oplus H^{s-1/2}(Y;E^-) \tag{DI5}

Fredholm for every real s≥1s\geq1.

2. The boundary signs and the shifted kernel

For a compactly supported smooth scalar or vector function ww in the collar, integration by parts with the stated inner-product convention gives the exact signs

2Im⁡(ϕDtw,w)X=∥γ0w∥Y2+∫Xϕ′(t)|w|2,2Im⁡(−ϕDtw,w)X=−∥γ0w∥Y2−∫Xϕ′(t)|w|2.(DI6) \begin{aligned} 2\operatorname{Im}(\phi D_tw,w)_X &=\|\gamma_0w\|_Y^2+\int_X\phi'(t)|w|^2,\\ 2\operatorname{Im}(-\phi D_tw,w)_X &=-\|\gamma_0w\|_Y^2-\int_X\phi'(t)|w|^2. \end{aligned} \tag{DI6}

Indeed, the first difference from its complex conjugate is −i∫ϕ∂t|w|2-i\int\phi\partial_t|w|^2, and the endpoint at t=0t=0 has the displayed positive sign. The second identity is its negative.

The scalar positive principal symbols of Λ±\Lambda^\pm and Λ\Lambda give the matrix Gårding lower bounds for the real parts of their quadratic forms. Derivatives of ϕ\phi, nonsymmetric lower-order terms, and the bounded collar patching terms contribute at most C∥u∥X2C\|u\|_X^2. Combining these facts with (DI6) gives

2Im⁡(Pu,u)X≥∥γ0u+∥Y2−∥γ0u−∥Y2−C∥u∥X2.(DI7) 2\operatorname{Im}(Pu,u)_X \geq \|\gamma_0u^+\|_Y^2- \|\gamma_0u^-\|_Y^2-C\|u\|_X^2. \tag{DI7}

For T>C/2T>C/2, set PT=P+iTIEP_T=P+iT I_E. If PTu=0P_Tu=0 and Bu=0Bu=0, then

0=2Im⁡(PTu,u)X≥∥γ0u+∥Y2+(2T−C)∥u∥X2.(DI8) 0=2\operatorname{Im}(P_Tu,u)_X \geq\|\gamma_0u^+\|_Y^2+(2T-C)\|u\|_X^2. \tag{DI8}

Thus u=0u=0. Elliptic regularity from Fredholm boundary problems with first-order Calderón defects first makes every weak kernel element smooth, so the calculation applies to the full kernel of (DI5).

3. The adjoint boundary relation and the shifted cokernel

Take the base realization s=1s=1. A pair (v,h)(v,h) in the annihilator of the range of (PT,B)1(P_T,B)_1 is smooth by the dual regularity theorem in Fredholm boundary problems with first-order Calderón defects. The leading normal terms give Green’s identity

(Pu,v)X−(u,P*v)X=i(γ0u+,γ0v+)Y−i(γ0u−,γ0v−)Y.(DI9) (Pu,v)_X-(u,P^*v)_X =i(\gamma_0u^+,\gamma_0v^+)_Y -i(\gamma_0u^-,\gamma_0v^-)_Y. \tag{DI9}

The annihilation equation is therefore

0=(u,(P*−iT)v)X+i(γ0u+,γ0v+)Y−i(γ0u−,γ0v−)Y+(γ0u−,h)Y.(DI10) 0=(u,(P^*-iT)v)_X +i(\gamma_0u^+,\gamma_0v^+)_Y -i(\gamma_0u^-,\gamma_0v^-)_Y +(\gamma_0u^-,h)_Y. \tag{DI10}

Interior tests and arbitrary boundary traces give, with no change of the conjugate-linear slot,

(P*−iT)v=0,γ0v+=0,h=−iγ0v−.(DI11) (P^*-iT)v=0,\qquad \gamma_0v^+=0,\qquad h=-i\gamma_0v^-. \tag{DI11}

Put Q=−P*Q=-P^*. Exchange the names of the two orthogonal summands:

EQ+=E−,EQ−=E+.(DI12) E_Q^+=E^-,\qquad E_Q^-=E^+. \tag{DI12}

The full differentiated-cutoff terms are computed below in (DA3)–(DA5). The collar principal part of QQ, written in the order EQ+⊕EQ−E_Q^+\oplus E_Q^-, is

(Dt+i(Λ−)*00−Dt+i(Λ+)*),BQv=γ0v+=γ0vQ−.(DI13) \begin{pmatrix} D_t+i(\Lambda^-)^*&0\\ 0&-D_t+i(\Lambda^+)^* \end{pmatrix},\qquad B_Qv=\gamma_0v^+=\gamma_0v_Q^-. \tag{DI13}

Multiplying the first equation in (DI11) by −1-1 gives (Q+iT)v=0(Q+iT)v=0, while BQv=0B_Qv=0. The estimate (DI8), with the same sufficiently large TT after enlarging CC once, gives v=0v=0, and then h=0h=0. Thus (PT,B)1(P_T,B)_1 is bijective. The path

(P+irTIE,B),0≤r≤1,(DI14) (P+i rT I_E,B),\qquad 0\leq r\leq1, \tag{DI14}

keeps both principal symbols fixed and stays Fredholm. Its index is constant, so ind⁡(P,B)1=0\operatorname{ind}(P,B)_1=0. Fredholm boundary problems with first-order Calderón defects identifies the same smooth kernel and adjoint obstruction spaces at every s≥1s\geq1. Consequently

ind⁡(P,B)s=0(s≥1).(DI15) \boxed{\operatorname{ind}(P,B)_s=0\quad(s\geq1).} \tag{DI15}

The full positivity remainder and adjoint coefficients

Here is a direct proof of the lower bounds used in (DI7), with their actual lower-order operators retained. If AA is any of Λ+\Lambda^+, Λ−\Lambda^-, or the interior Λ\Lambda, its principal symbol is λAI\lambda_A I, with λA>0\lambda_A>0 off the zero section. Quantize the degree-one-half symbol λAI\sqrt{\lambda_A}I to an operator RAR_A. Use the full bundle composition and adjoint formulas in Symbols, operators and Sobolev scales and Detecting regularity without choosing coordinates. Their principal product is exactly λAI\lambda_A I. Hence

HA=A+A*2−RA*RA∈Ψ0,Re⁡(Aw,w)=∥RAw∥22+(HAw,w),CA=∥HA∥L2→L2.(DA1) H_A=\frac{A+A^*}{2}-R_A^*R_A\in\Psi^0, \qquad \operatorname{Re}(Aw,w)=\|R_Aw\|_2^2+(H_Aw,w), \qquad C_A=\|H_A\|_{L^2\to L^2}. \tag{DA1}

The principal symbol comparison proves the order-zero claim; the linked order-zero mapping theorem makes CAC_A finite. The formula retains the entire lower-order operator HAH_A, rather than replacing AA by its positive principal symbol. It gives Re⁡(Aw,w)≥−CA∥w∥22\operatorname{Re}(Aw,w)\geq-C_A\|w\|_2^2. For tangential families take the supremum of these constants on the compact support of ϕ\phi. Smooth family quantization and the finite seminorm bounds give finite suprema. For the interior operator the input (1−ϕ)u(1-\phi)u is supported away from the boundary; insert cutoffs equal to one on that support and use the same calculation on interior charts. All errors from these cutoffs are included in HΛH_\Lambda.

At each fixed tt, the scalar ϕ(t)\phi(t) commutes with the tangential operators. It is nonnegative, so multiplying (DA1) by it and integrating in tt gives the tangential lower bound without taking a square root of the cutoff. The exact interior quadratic form is Im⁡(i(1−ϕ)Λ(1−ϕ)u,u)=Re⁡(Λ(1−ϕ)u,(1−ϕ)u)\operatorname{Im}(i(1-\phi)\Lambda(1-\phi)u,u) =\operatorname{Re}(\Lambda(1-\phi)u,(1-\phi)u). Combining these statements with both full identities in (DI6) proves (DI7) with the concrete choice

Cb=max⁡(CΛ+,CΛ−),Ci=CΛ,CP=∥ϕ′∥∞+2Cb+2Ci.(DA2) C_b=\max(C_{\Lambda^+},C_{\Lambda^-}),\qquad C_i=C_\Lambda,\qquad C_P=\|\phi'\|_\infty+2C_b+2C_i. \tag{DA2}

No self-adjointness of AA was assumed. Its geometric adjoint has the same real quadratic form, so the same constants work for Λ*\Lambda^* and (Λ±)*(\Lambda^\pm)^*.

To compute Q=−P*Q=-P^*, use the original product metric and density. Multiplication by the real ϕ\phi is self-adjoint and Dtϕ=ϕDt−iϕ′D_t\phi=\phi D_t-i\phi'. The full collar adjoint, before exchanging summands, is

(Pb)*=(ϕDt−iϕ′−iϕ(Λ+)*00−ϕDt+iϕ′−iϕ(Λ−)*),(Pi)*=−i(1−ϕ)Λ*(1−ϕ).(DA3) (P^b)^*= \begin{pmatrix} \phi D_t-i\phi'-i\phi(\Lambda^+)^*&0\\ 0&-\phi D_t+i\phi'-i\phi(\Lambda^-)^* \end{pmatrix},\qquad (P^i)^*=-i(1-\phi)\Lambda^*(1-\phi). \tag{DA3}

Tangential operators commute with ϕ(t)\phi(t), including when their coefficients depend on tt. Let S(v+,v−)=(v−,v+)S(v^+,v^-)=(v^-,v^+) be the unitary change of order on the collar. In precisely the order (DI12) the full operator is

SQbS−1=ϕ(Dt+i(Λ−)*00−Dt+i(Λ+)*)+(−iϕ′00iϕ′),SQiS−1=i(1−ϕ)SΛ*S−1(1−ϕ),BQv=γ0v+=γ0(Sv)−.(DA4) \begin{aligned} SQ^bS^{-1} &=\phi\begin{pmatrix} D_t+i(\Lambda^-)^*&0\\ 0&-D_t+i(\Lambda^+)^* \end{pmatrix} +\begin{pmatrix}-i\phi'&0\\0&i\phi'\end{pmatrix},\\ SQ^iS^{-1}&=i(1-\phi)S\Lambda^*S^{-1}(1-\phi),\\ B_Qv&=\gamma_0v^+=\gamma_0(Sv)^-. \end{aligned} \tag{DA4}

The second formula is a collar expression only; it does not assert that the splitting or SS extends over the whole interior bundle. On the original bundle the full interior operator is i(1−ϕ)Λ*(1−ϕ)i(1-\phi)\Lambda^*(1-\phi). These descriptions give the same quadratic form on their respective charts. The displayed multiplication remainder contributes exactly

2Im⁡((−iϕ′00iϕ′)vQ,vQ)=−2∫ϕ′|vQ+|2+2∫ϕ′|vQ−|2.(DA5) 2\operatorname{Im} \left(\begin{pmatrix}-i\phi'&0\\0&i\phi'\end{pmatrix} v_Q,v_Q\right) =-2\int\phi'|v_Q^+|^2+2\int\phi'|v_Q^-|^2. \tag{DA5}

Its absolute value is at most 2∥ϕ′∥∞∥v∥222\|\phi'\|_\infty\|v\|_2^2. Therefore (DI7) holds for QQ with CQ=CP+2∥ϕ′∥∞C_Q=C_P+2\|\phi'\|_\infty, with its plus and minus boundary spaces exactly as in (DI12). A single T>max⁡(CP,CQ)/2T>\max(C_P,C_Q)/2 supplies both estimates used in (DI8) and (DI11). This proves the shifted bijectivity without discarding either differentiated-cutoff term.

4. A surjective boundary constraint preserves the index

Let X0,Y0,Z0X_0,Y_0,Z_0 be Hilbert spaces, let A=(P,B):X0→Y0⊕Z0A=(P,B):X_0\to Y_0\oplus Z_0 be Fredholm, and suppose B:X0→Z0B:X_0\to Z_0 is surjective. Its restriction to (ker⁡B)⟂(\ker B)^\perp is a bounded bijection onto Z0Z_0, so the Banach inverse theorem supplies a bounded right inverse R:Z0→X0R:Z_0\to X_0. Set XB=ker⁡BX_B=\ker B. Define

J:XB⊕Z0→X0,J(x,z)=x+Rz,J−1u=(u−RBu,Bu),S:Y0⊕Z0→Y0⊕Z0,S(y,z)=(y−PRz,z).(DI16) \begin{aligned} J:X_B\oplus Z_0&\longrightarrow X_0, &J(x,z)&=x+Rz,\\ J^{-1}u&=(u-RBu,Bu),\\ S:Y_0\oplus Z_0&\longrightarrow Y_0\oplus Z_0, &S(y,z)&=(y-PRz,z). \end{aligned} \tag{DI16}

Both JJ and SS are bounded isomorphisms, and direct substitution gives the exact triangular reduction

SAJ(x,z)=(Px,z),SAJ=(P|XB)⊕IZ0.(DI17) SAJ(x,z)=(Px,z),\qquad SAJ=(P|_{X_B})\oplus I_{Z_0}. \tag{DI17}

It follows that P|XB:XB→Y0P|_{X_B}:X_B\to Y_0 is Fredholm and that its kernel, cokernel, and index agree with those of AA:

ker⁡(P|XB)=ker⁡A,coker⁡(P|XB)≃coker⁡A,ind⁡(P|XB)=ind⁡A.(DI18) \ker(P|_{X_B})=\ker A,\qquad \operatorname{coker}(P|_{X_B})\simeq\operatorname{coker}A,\qquad \operatorname{ind}(P|_{X_B})=\operatorname{ind}A. \tag{DI18}

This also proves the cokernel isomorphism, rather than only an index count.

5. The reflected complementary half

Now let P1:E1→F1P_1:E_1\to F_1 be any elliptic first-order mixed operator on a copy X1=XX_1=X whose collar part, after a fixed collar identification κ:F1≃E1\kappa:F_1\simeq E_1, is the split expression in (DI2). Its compactly supported interior part need not be the positive model used to prove (DI15). Give a second copy X2=XX_2=X the inward coordinate t2≥0t_2\geq0, copy the split bundle E2+⊕E2−E_2^+\oplus E_2^-, and set near its boundary

P2b=ϕ(−Dt2+iΛ+00Dt2+iΛ−),B2u2=γ0u2+.(DI19) P_2^b=\phi \begin{pmatrix} -D_{t_2}+i\Lambda^+&0\\ 0&D_{t_2}+i\Lambda^- \end{pmatrix},\qquad B_2u_2=\gamma_0u_2^+. \tag{DI19}

Complete P2bP_2^b away from the boundary by the reflected version of the positive interior term in (DI2): choose Λ2∈Ψphg1(X2∘;E2)\Lambda_2\in\Psi_{\mathrm{phg}}^1(X_2^\circ;E_2) with positive scalar principal symbol and set P2i=i(1−ϕ)Λ2(1−ϕ)P_2^i=i(1-\phi)\Lambda_2(1-\phi), P2=P2b+P2iP_2=P_2^b+P_2^i. Its frozen solutions are

u2+(t2)=e−λ+t2u2+(0),u2−(t2)=eλ−t2u2−(0).(DI20) u_2^+(t_2)=e^{-\lambda^+t_2}u_2^+(0),\qquad u_2^-(t_2)=e^{\lambda^-t_2}u_2^-(0). \tag{DI20}

Thus its stable trace is exactly E+E^+, and B2B_2 is the identity there. After ordering the summands as E2−⊕E2+E_2^-\oplus E_2^+, this is the model of Section 1. Hence

ind⁡(P2,B2)=0.(DI21) \operatorname{ind}(P_2,B_2)=0. \tag{DI21}

6. Smooth gluing of the doubled bundles and operator

Glue X1X_1 and X2X_2 along YY. In the resulting two-sided collar use the signed coordinate

r=t1on X1,r=−t2on X2.(DI22) r=t_1\quad\hbox{on }X_1,\qquad r=-t_2\quad\hbox{on }X_2. \tag{DI22}

Glue E1E_1 to E2E_2 by the identity on E+⊕E−E^+\oplus E^-, producing Ê\widehat E. Glue F1F_1 to the target E2E_2 by the collar identification κ\kappa, producing F̂\widehat F. Since Dt2=−DrD_{t_2}=-D_r, both collar formulas become

P̂b=(Dr+iΛ+00−Dr+iΛ−)(DI23) \widehat P^b= \begin{pmatrix} D_r+i\Lambda^+&0\\ 0&-D_r+i\Lambda^- \end{pmatrix} \tag{DI23}

on a full two-sided neighborhood of the seam. The interior terms vanish on that neighborhood because 1−ϕ=01-\phi=0 there. Thus all coefficient jets and the bundle maps match, and the piecewise operator defines

P̂:Hs(X̂;Ê)→Hs−1(X̂;F̂).(DI24) \widehat P:H^s(\widehat X;\widehat E) \longrightarrow H^{s-1}(\widehat X;\widehat F). \tag{DI24}

This is a smooth mixed collar operator. It need not be an ordinary global pseudodifferential operator, because Λ±\Lambda^\pm act only in the tangential variables near the seam. This proves both the construction and the distinction between a mixed collar operator and an ordinary global pseudodifferential operator.

7. The coupling path and the H1H^1 gluing domain

At the base level put

𝒳=H‾1(X1;E1)⊕H‾1(X2;E2),𝒴=L2(X1;F1)⊕L2(X2;E2),𝒵=H1/2(Y;E+)⊕H1/2(Y;E−),Cτ(u1,u2)=(γ0u2+−τγ0u1+,γ0u1−−τγ0u2−),0≤τ≤1,𝒜τ=(P1⊕P2,Cτ):𝒳→𝒴⊕𝒵.(DI25) \begin{aligned} \mathcal X&=\bar H_1(X_1;E_1)\oplus\bar H_1(X_2;E_2),\\ \mathcal Y&=L^2(X_1;F_1)\oplus L^2(X_2;E_2),\\ \mathcal Z&=H^{1/2}(Y;E^+)\oplus H^{1/2}(Y;E^-),\\ C_\tau(u_1,u_2)&= \bigl(\gamma_0u_2^+-\tau\gamma_0u_1^+, \gamma_0u_1^--\tau\gamma_0u_2^-\bigr), \qquad0\leq\tau\leq1,\\ \mathcal A_\tau&=(P_1\oplus P_2,C_\tau): \mathcal X\longrightarrow\mathcal Y\oplus\mathcal Z. \end{aligned} \tag{DI25}

The stable traces of P1⊕P2P_1\oplus P_2 have the form (0,a−;a+,0)(0,a^-;a^+,0). On them,

Cτ(0,a−;a+,0)=(a+,a−),(DI26) C_\tau(0,a^-;a^+,0)=(a^+,a^-), \tag{DI26}

independently of τ\tau. The principal boundary map is therefore an isomorphism for every parameter. The generalized Fredholm theorem and homotopy stability give

ind⁡𝒜τ=ind⁡𝒜0(0≤τ≤1).(DI27) \operatorname{ind}\mathcal A_\tau =\operatorname{ind}\mathcal A_0 \qquad(0\leq\tau\leq1). \tag{DI27}

The coupling map itself is surjective. For data (a+,b−)(a^+,b^-), choose fixed Sobolev trace extensions with

γ0u2+=a+,γ0u1−=b−,γ0u1+=0,γ0u2−=0.(DI28) \gamma_0u_2^+=a^+,\quad \gamma_0u_1^-=b^-,\quad \gamma_0u_1^+=0,\quad \gamma_0u_2^-=0. \tag{DI28}

This gives one bounded right inverse R:𝒵→𝒳R:\mathcal Z\to\mathcal X for every CτC_\tau. Apply (DI18):

ind⁡((P1⊕P2)|ker⁡Cτ)=ind⁡𝒜τ.(DI29) \operatorname{ind}\bigl((P_1\oplus P_2)|_{\ker C_\tau}\bigr) =\operatorname{ind}\mathcal A_\tau. \tag{DI29}

At τ=0\tau=0, the two constraints are γ0u2+=0\gamma_0u_2^+=0 and γ0u1−=0\gamma_0u_1^-=0. After a fixed permutation of the target factors, the constrained operator is the direct sum of P1|ker⁡B1P_1|_{\ker B_1} and P2|ker⁡B2P_2|_{\ker B_2}. Equation (DI18) identifies their indices with those of (P1,B1)(P_1,B_1) and (P2,B2)(P_2,B_2). Equations (DI21), (DI27), and (DI29) give

ind⁡((P1⊕P2)|ker⁡C1)=ind⁡(P1,B1).(DI30) \operatorname{ind}\bigl((P_1\oplus P_2)|_{\ker C_1}\bigr) =\operatorname{ind}(P_1,B_1). \tag{DI30}

The varying kernels can also be trivialized explicitly. Put

K(u1,u2)=(γ0u1+,γ0u2−),Cτ=C0−τK,KR=0,C0R=I𝒵.(DI31) K(u_1,u_2)=(\gamma_0u_1^+,\gamma_0u_2^-),\qquad C_\tau=C_0-\tau K,\qquad KR=0,\qquad C_0R=I_{\mathcal Z}. \tag{DI31}

For x∈ker⁡C0x\in\ker C_0, define

Jτx=x+τRKx.(DI32) J_\tau x=x+\tau RKx. \tag{DI32}

Then CτJτx=0C_\tau J_\tau x=0. Its inverse on ker⁡Cτ\ker C_\tau is y↦y−RC0yy\mapsto y-RC_0y, because C0y=τKyC_0y=\tau Ky and KR=0KR=0. Hence Jτ:ker⁡C0→ker⁡CτJ_\tau:\ker C_0\to\ker C_\tau is a continuous family of bounded isomorphisms. This makes the Fredholm deformation in (DI30) an actual fixed-domain path, not merely an index comparison between unnamed spaces.

At τ=1\tau=1, the constraint is equality of the two full boundary traces. The exact H1H^1 gluing identity is

H1(X̂;Ê)≃{(u1,u2)∈𝒳:γ0u1+=γ0u2+,γ0u1−=γ0u2−}=ker⁡C1.(DI33) H^1(\widehat X;\widehat E) \simeq\{(u_1,u_2)\in\mathcal X: \gamma_0u_1^+=\gamma_0u_2^+, \ \gamma_0u_1^-=\gamma_0u_2^-\} =\ker C_1. \tag{DI33}

One direction follows by restriction. Conversely, the distributional first derivative of a piecewise H1H^1 section has a seam delta whose coefficient is its trace jump; equality of traces removes that delta, leaving every weak first derivative in L2L^2. Under (DI33), the constrained piecewise operator is exactly (DI24). Thus

ind⁡P̂=ind⁡(P1,B1)(DI34) \boxed{\operatorname{ind}\widehat P =\operatorname{ind}(P_1,B_1)} \tag{DI34}

at H1→L2H^1\to L^2. The statement at this step is deliberately made at s=1s=1: equality of value traces alone is not the complete gluing condition for arbitrary higher Sobolev order.

The stable summands on the reflected halves, the coupling path, and the two index strata

The left panel keeps both inward coordinates and the signed collar coordinate visible. The middle panel shows why the boundary symbol stays invertible throughout the coupling. The right panel records the complementary index that distinguishes equality from a one-half formula.

The complete seam distribution

We verify the asserted gluing in the actual product collar. Use its signed coordinate r=t1=−t2r=t_1=-t_2, product density dμYdrd\mu_Y\,dr, and the domain-bundle identification from Section 6. In any product frame, write

u(y,r)={u1(y,r),r>0,u2(y,−r),r<0,j(y)=γ0u1(y)−γ0u2(y).(DG1) u(y,r)=\begin{cases} u_1(y,r),&r>0,\\ u_2(y,-r),&r<0, \end{cases} \qquad j(y)=\gamma_0u_1(y)-\gamma_0u_2(y). \tag{DG1}

Let HH be the Heaviside function; its value at zero does not affect the distribution. Integration by parts on each half, with the common boundary measure, gives every first derivative:

Dru=H(r)Dt1u1(y,r)−H(−r)Dt2u2(y,−r)−ij(y)δ(r),Dyau=H(r)Dyau1(y,r)+H(−r)Dyau2(y,−r),1≤a≤dim⁡Y.(DG2) \begin{aligned} D_ru&=H(r)D_{t_1}u_1(y,r) -H(-r)D_{t_2}u_2(y,-r)-i j(y)\delta(r),\\ D_{y_a}u&=H(r)D_{y_a}u_1(y,r) +H(-r)D_{y_a}u_2(y,-r), \qquad 1\leq a\leq\dim Y. \end{aligned} \tag{DG2}

For smooth half-sections this follows directly from ∂rH(r)=δ(r)\partial_rH(r)=\delta(r) and ∂rH(−r)=−δ(r)\partial_rH(-r)=-\delta(r). For half-sections in H1H^1, approximate on each half by smooth sections. The continuous H1→H1/2H^1\to H^{1/2} trace map from Mixed symbols on every real two-parameter Sobolev scale makes the seam coefficients converge, while the half-space functions and weak derivatives converge in L2L^2. Pairing with each compactly supported smooth test section passes the identities to the limit. Thus no differentiability of the traces is assumed.

All bulk terms in (DG2) are in L2L^2. If j=0j=0, every weak first derivative is in L2L^2, proving the converse in (DI33). If j≠0j\ne0, choose a smooth compactly supported tangential test section bb with ⟨j,b⟩≠0\langle j,b\rangle\ne0. Such a test exists because a nonzero distribution cannot vanish on every test section. Choose h∈Cc∞(ℝ)h\in C_c^\infty(\mathbb R) with h(0)=1h(0)=1, and put hε(r)=h(r/ε)h_\varepsilon(r)=h(r/\varepsilon). Then

⟨jδ(r),b(y)hε(r)⟩=⟨j,b⟩,∥b(y)hε(r)∥L2(dμYdr)=ε1/2∥b∥L2(dμY)∥h∥L2(dr).(DG3) \begin{aligned} \langle j\delta(r),b(y)h_\varepsilon(r)\rangle &=\langle j,b\rangle,\\ \|b(y)h_\varepsilon(r)\|_{L^2(d\mu_Y\,dr)} &=\varepsilon^{1/2}\|b\|_{L^2(d\mu_Y)}\|h\|_{L^2(dr)}. \end{aligned} \tag{DG3}

An L2L^2 distribution has pairings bounded by its L2L^2 norm times the second line, which tends to zero. The first line is a fixed nonzero number. Therefore the seam delta is not in L2L^2. In (DG2) it cannot cancel an L2L^2 bulk term. This proves necessity of trace equality as well.

The restriction map and the piecewise inverse are bounded for the finite-chart H1H^1 norms: (DG2) with zero jump bounds every coordinate first derivative by the two half-norms, and restriction bounds each half-norm by the global norm. Smooth transition matrices contribute their actual bounded first derivatives; the fixed product metric and density compare these finitely many chart norms. This proves the bounded isomorphism (DI33), including its full domain. It makes no claim that matching values alone glues higher Sobolev orders, which would require the corresponding normal-jet conditions.

8. Approximation by ordinary operators and the symbol index

In the signed collar let (η,ρ)(\eta,\rho) be the covariables dual to (y,r)(y,r). The principal symbol of (DI23), extended by the ordinary interior symbols, is continuous and degree one on T*X̂\0T^*\widehat X\setminus0. At η=0\eta=0, positive homogeneity gives λ±(y,0)=0\lambda^\pm(y,0)=0, so its two collar blocks equal ρ\rho and −ρ-\rho. Thus it is still invertible there. Denote this continuous elliptic symbol by

p̂:T*X̂\0→Hom⁡(π*Ê,π*F̂).(DI35) \widehat p:T^*\widehat X\setminus0 \longrightarrow \operatorname{Hom}(\pi^*\widehat E,\pi^*\widehat F). \tag{DI35}

We now verify the operator-norm approximation required by Symbols, finite defects, and the index on a closed manifold rather than inferring it from pointwise smoothing. Choose the scalar cutoff from Symbols, finite defects, and the index on a closed manifold, Section 10,

χ(η,ρ)=1if |ρ|≤max⁡(1,|η|),χ(η,ρ)=0if |ρ|≥2max⁡(1,|η|),(DI36) \chi(\eta,\rho)=1\ \text{if }|\rho|\leq\max(1,|\eta|), \qquad \chi(\eta,\rho)=0\ \text{if }|\rho|\geq2\max(1,|\eta|), \tag{DI36}

and for 0<ε≤10<\varepsilon\leq1 replace every partial tangential symbol a(y,r,η)a(y,r,\eta) of order one by

aε(y,r,η,ρ)=a(y,r,η)χ(η,ερ).(DI37) a_\varepsilon(y,r,\eta,\rho) =a(y,r,\eta)\chi(\eta,\varepsilon\rho). \tag{DI37}

For fixed ε>0\varepsilon>0, this is a classical full-variable symbol. On the total unit cosphere its leading term differs from the continuous zero extension of the partial leading symbol only where |η|≤2ε|ρ||\eta|\leq2\varepsilon|\rho|. Since that leading term has size at most C|η|C|\eta|,

sup|η|2+ρ2=1∥aε,1−a1∥≤Cε.(DI38) \sup_{|\eta|^2+\rho^2=1} \|a_{\varepsilon,1}-a_1\|\leq C\varepsilon. \tag{DI38}

The exact partial-operator estimate in Section 10 of Symbols, finite defects, and the index on a closed manifold gives, for every real ss,

∥Op⁡(aε)−a(y,r,Dy)∥Hs→Hs−1≤Csε.(DI39) \|\operatorname{Op}(a_\varepsilon)-a(y,r,D_y)\| _{H^s\to H^{s-1}}\leq C_s\varepsilon. \tag{DI39}

Apply this construction in a finite signed-collar atlas. Keep DrD_r and every ordinary interior operator unchanged, and patch with the same fixed bundle charts and cutoffs. The compactness of X̂\widehat X turns the local estimates into ordinary classical operators P̂ε\widehat P_\varepsilon satisfying

P̂ε→P̂in ℒ(Hs,Hs−1) for every real s,p̂ε→p̂ uniformly on S*X̂.(DI40) \widehat P_\varepsilon\longrightarrow\widehat P \quad\text{in }\mathcal L(H^s,H^{s-1}) \text{ for every real }s,\qquad \widehat p_\varepsilon\longrightarrow\widehat p \text{ uniformly on }S^*\widehat X. \tag{DI40}

Uniform invertibility of p̂\widehat p on the compact cosphere makes p̂ε\widehat p_\varepsilon elliptic for all sufficiently small ε\varepsilon. The norm-limit theorem and its continuous-symbol identity, proved in Sections 5 and 8 of that lesson, now give

ind⁡P̂=sind⁡(p̂).(DI41) \operatorname{ind}\widehat P =\operatorname{sind}(\widehat p). \tag{DI41}

Because (DI40) is compatible at every Sobolev order, the compatible norm-limit conclusion in that lesson identifies the same smooth kernel and adjoint obstruction spaces at all levels. Hence (DI34) and (DI41) hold for every s≥1s\geq1:

ind⁡(P1,B1)s=ind⁡(P̂:Hs→Hs−1)=sind⁡(p̂).(DI42) \boxed{\operatorname{ind}(P_1,B_1)_s =\operatorname{ind}(\widehat P:H^s\to H^{s-1}) =\operatorname{sind}(\widehat p).} \tag{DI42}

This proves the doubled symbol-index formula with the precise approximation hypotheses.

9. The two exact kinds of doubled index datum

Two published index wordings become exact once the complementary half is made part of the datum. Let A:XA→YAA:X_A\to Y_A be Fredholm. A deformation double datum for AA is a Fredholm operator D:XD→YDD:X_D\to Y_D, a Fredholm complement C:XC→YCC:X_C\to Y_C, fixed domain and target isomorphisms at the two endpoints, and a continuous fixed-space Fredholm path joining DD to A⊕CA\oplus C. Homotopy invariance and direct-sum additivity give the exact morphism from this space of data to the integers:

ℑ(D,C)=ind⁡D=ind⁡A+ind⁡C.(DI43) \mathfrak I(D,C)=\operatorname{ind}D =\operatorname{ind}A+\operatorname{ind}C. \tag{DI43}

There are two mathematically different strata:

𝔇0(A)={(D,C):ind⁡C=0}⇒ind⁡D=ind⁡A,𝔇=(A)={(D,C):ind⁡C=ind⁡A}⇒ind⁡D=2ind⁡A.(DI44) \begin{array}{lll} \mathfrak D_0(A)=\{(D,C):\operatorname{ind}C=0\} &\Longrightarrow&\operatorname{ind}D=\operatorname{ind}A,\\[1mm] \mathfrak D_{=}(A)=\{(D,C):\operatorname{ind}C=\operatorname{ind}A\} &\Longrightarrow&\operatorname{ind}D=2\operatorname{ind}A. \end{array} \tag{DI44}

Equations (DI31)–(DI34) supply all the data for the geometric double: the fixed-domain path is (P1⊕P2)Jτ(P_1\oplus P_2)J_\tau, its endpoint complement is (P2,B2)(P_2,B_2), and (DI21) gives its index zero. Therefore the concrete operator P̂\widehat P belongs to 𝔇0(P1,B1)\mathfrak D_0(P_1,B_1), and its formula is the equality in (DI42).

A one-half formula is valid only after a separately specified double datum is proved to lie in 𝔇=(A)\mathfrak D_{=}(A):

ind⁡A=12ind⁡D⇔ind⁡C=ind⁡Awithin a deformation double datum.(DI45) \operatorname{ind}A=\frac12\operatorname{ind}D \quad\Longleftrightarrow\quad \operatorname{ind}C=\operatorname{ind}A \quad\text{within a deformation double datum.} \tag{DI45}

Thus the introductory one-half wording cannot be applied to the operator P̂\widehat P constructed here. Proposition 20.3.2’s zero-index complementary half gives equality. A summary using one half must name and prove a different same-index complement and its deformation data. Without that information the one-half statement is underdetermined; it is not a replacement for (DI42). This is the exact conceptual correction. It keeps the equality theorem and the conditional one-half formula attached to their respective doubled constructions.

10. Scalar sign check

For E+=E−=ℂE^+=E^-=\mathbb C and Λ+=Λ−=|Dy|\Lambda^+=\Lambda^-=|D_y|, the first-half stable mode is (0,e−|η|t1a−)(0,e^{-|\eta|t_1}a^-), while the second-half stable mode is (e−|η|t2a+,0)(e^{-|\eta|t_2}a^+,0). The coupling gives

Cτ(0,e−|η|t1a−;e−|η|t2a+,0)=(a+,a−).(DI46) C_\tau(0,e^{-|\eta|t_1}a^-; e^{-|\eta|t_2}a^+,0)=(a^+,a^-). \tag{DI46}

In the signed coordinate, both halves have symbol diag⁡(ρ+i|η|,−ρ+i|η|)\operatorname{diag}(\rho+i|\eta|,-\rho+i|\eta|). This checks the reflection sign, the stable summands, the order of the coupling components, and the seam compatibility in one calculation.

11. Three concrete models

A finite-dimensional boundary constraint

Take X0=ℂ3X_0=\mathbb C^3, Y0=ℂ2Y_0=\mathbb C^2, Z0=ℂZ_0=\mathbb C, and

P(x1,x2,x3)=(x1,x2+x3),B(x1,x2,x3)=x3.(DI47) P(x_1,x_2,x_3)=(x_1,x_2+x_3),\qquad B(x_1,x_2,x_3)=x_3. \tag{DI47}

The combined map (P,B)(P,B) is invertible. On ker⁡B\ker B, the restricted map is (x1,x2,0)↦(x1,x2)(x_1,x_2,0)\mapsto(x_1,x_2), also invertible. The right inverse Rz=(0,0,z)Rz=(0,0,z) makes (DI16)–(DI17) an ordinary triangular matrix factorization. This finite model displays the exact cokernel isomorphism used for Sobolev boundary spaces.

Why matching values is the right H1H^1 condition

Let u1∈H1([0,1])u_1\in H^1([0,1]) and u2∈H1([−1,0])u_2\in H^1([-1,0]), and join them as a piecewise function uu. For every compactly supported smooth test function φ\varphi, integration by parts on the two halves gives

⟨∂ru,φ⟩=∫−10u2′φdr+∫01u1′φdr+(u1(0)−u2(0))φ(0).(DI48) \langle \partial_ru,\varphi\rangle =\int_{-1}^0u_2'\varphi\,dr+\int_0^1u_1'\varphi\,dr +(u_1(0)-u_2(0))\varphi(0). \tag{DI48}

The distributional derivative lies in L2L^2 exactly when the trace jump vanishes. This is the one-dimensional form of (DI33). Higher Sobolev gluing needs further matching derivatives, which is why the index identification is first made at H1→L2H^1\to L^2.

The angular cap estimate

For the scalar partial symbol a(η)=|η|a(\eta)=|\eta|, the leading-symbol difference in (DI37) is supported where |η|≤2ε|ρ||\eta|\leq2\varepsilon|\rho|. On |η|2+ρ2=1|\eta|^2+\rho^2=1,

|aε,1−a1|≤|η|≤2ε.(DI49) |a_{\varepsilon,1}-a_1|\leq|\eta|\leq2\varepsilon. \tag{DI49}

This is the exact small factor that turns the tangential operator into an operator-norm limit of ordinary full-variable pseudodifferential operators.

12. Exercises with complete solutions

1. Recover the two endpoint signs. Starting from Dt=−i∂tD_t=-i\partial_t, prove (DI6) without assuming that ww is scalar.

Solution. The fiber metric gives ∂t|w|2=(∂tw,w)+(w,∂tw)\partial_t|w|^2=(\partial_tw,w)+(w,\partial_tw). Hence

2iIm⁡(ϕDtw,w)=−i∫0∞ϕ∂t|w|2dt=i|w(0)|2+i∫0∞ϕ′|w|2dt.(DI50) \begin{aligned} 2i\operatorname{Im}(\phi D_tw,w) &=-i\int_0^\infty\phi\,\partial_t|w|^2dt\\ &=i|w(0)|^2+i\int_0^\infty\phi'|w|^2dt. \end{aligned} \tag{DI50}

Divide by ii. Replacing DtD_t by −Dt-D_t reverses both terms. The calculation is component-free, so it holds for bundle-valued ww.

2. Check the cokernel multiplier. In (DI10), use arbitrary minus traces to recover the exact formula for hh.

Solution. For a=γ0u−a=\gamma_0u^-, the minus terms are −i(a,γ0v−)+(a,h)-i(a,\gamma_0v^-)+(a,h). Since the inner product is linear in the first entry, (a,iγ0v−)=−i(a,γ0v−)(a,i\gamma_0v^-)=-i(a,\gamma_0v^-). Vanishing for every aa gives

h+iγ0v−=0,h=−iγ0v−.(DI51) h+i\gamma_0v^-=0,\qquad h=-i\gamma_0v^-. \tag{DI51}

3. Verify the triangular reduction. Prove directly that the maps JJ and SS in (DI16) are invertible and calculate SAJSAJ.

Solution. Because BR=IBR=I, u−RBu∈ker⁡Bu-RBu\in\ker B, so the displayed formula for J−1J^{-1} is defined. Direct substitution gives J−1J(x,z)=(x,z)J^{-1}J(x,z)=(x,z) and JJ−1u=uJJ^{-1}u=u. The inverse of SS is S−1(y,z)=(y+PRz,z)S^{-1}(y,z)=(y+PRz,z). Finally,

SAJ(x,z)=S(Px+PRz,z)=(Px,z).(DI52) SAJ(x,z)=S(Px+PRz,z)=(Px,z). \tag{DI52}

4. Trivialize the coupling kernels. Prove that JτJ_\tau in (DI32) maps ker⁡C0\ker C_0 bijectively onto ker⁡Cτ\ker C_\tau.

Solution. For x∈ker⁡C0x\in\ker C_0, use Cτ=C0−τKC_\tau=C_0-\tau K, C0R=IC_0R=I, and KR=0KR=0:

Cτ(x+τRKx)=0+τKx−τKx−τ2KRKx=0.(DI53) C_\tau(x+\tau RKx)=0+\tau Kx-\tau Kx-\tau^2KRKx=0. \tag{DI53}

If y∈ker⁡Cτy\in\ker C_\tau, then C0y=τKyC_0y=\tau Ky. Put x=y−RC0yx=y-RC_0y. This lies in ker⁡C0\ker C_0, Kx=KyKx=Ky, and x+τRKx=yx+\tau RKx=y. Thus y↦y−RC0yy\mapsto y-RC_0y is the inverse.

5. Test the two index strata. Let ind⁡A=−3\operatorname{ind}A=-3. Compute the doubled index for a zero-index complement and for a same-index complement.

Solution. Equation (DI43) gives

ind⁡D={−3,ind⁡C=0,−6,ind⁡C=−3.(DI54) \operatorname{ind}D= \begin{cases} -3,&\operatorname{ind}C=0,\\ -6,&\operatorname{ind}C=-3. \end{cases} \tag{DI54}

The first construction gives equality. The second gives ind⁡A=12ind⁡D=−3\operatorname{ind}A=\frac12\operatorname{ind}D=-3. Moving the factor 1/21/2 between these two data would give the wrong integer.

6. Why must the approximation be in operator norm? Explain why pointwise convergence of the full symbols in (DI37) would not by itself prove the index identity.

Solution. Fredholm index stability applies to bounded operators in their operator-norm topology. Pointwise symbol convergence gives no uniform control on the maps between Sobolev spaces and does not exclude a loss concentrated in a shrinking angular cap. The partial-operator estimate supplies the needed bound:

∥P̂ε−P̂∥Hs→Hs−1≤Csε.(DI55) \|\widehat P_\varepsilon-\widehat P\|_{H^s\to H^{s-1}} \leq C_s\varepsilon. \tag{DI55}

Together with uniform principal-symbol convergence and ellipticity, this is exactly the hypothesis of the norm-limit index theorem.

13. Reading notes and references

Equations (DI43)–(DI45) identify the missing datum and give the exact correction for the actual chosen complementary half.

Stable modes and the algebra of boundary data proves the stable-space algebra used in Sections 1 and 5. Fredholm boundary problems with first-order Calderón defects supplies the generalized Fredholm and regularity theorem. Symbols, finite defects, and the index on a closed manifold proves the partial-operator approximation and continuous symbol-index theorem used in Section 8.

The arguments, examples, figure, and exercises in this lesson are independently written.

Written and dedicated to the public domain by Codex under CC0 1.0.