Solution (source code)

= Solution

Use real <differential forms> and the usual convention that a <smooth manifold> has no boundary unless one is specified. On a noncompact manifold the formal-adjoint identity below uses compactly supported forms. On a compact manifold it has no boundary contribution; if a boundary were allowed, extra boundary conditions would be necessary.

The <Riemannian metric> gives the dual <inner product> on $T_x^*M$ and hence the <inner product on exterior powers of a cotangent space>:
$$
\boxed{\langle\xi_1\wedge\cdots\wedge\xi_p,\zeta_1\wedge\cdots\wedge\zeta_p\rangle=\det\bigl(\langle\xi_i,\zeta_j\rangle\bigr).}
$$
Extend this bilinearly. Equivalently, the increasing wedge products of an orthonormal coframe are an orthonormal basis of $\Lambda^pT_x^*M$. The <Hodge star operator> is the unique map $*:\Lambda^pT_x^*M\to\Lambda^{n-p}T_x^*M$ satisfying $\alpha\wedge*\beta=\langle\alpha,\beta\rangle\operatorname{vol}_g$. If $e^1,\ldots,e^n$ is a positively oriented orthonormal coframe, $*e^I=\varepsilon(I,I^c)e^{I^c}$. Moving the $p$ factors of $I$ across the $n-p$ factors of its complement changes the sign by $(-1)^{p(n-p)}$, so
$$
\boxed{*^2\alpha=(-1)^{p(n-p)}\alpha.}
$$
The basis formula also proves that $*$ is an isometry.

For $\beta\in\Omega^{p-1}(M)$ and $\alpha\in\Omega^p(M)$ with compactly supported product, the <Stokes theorem> and the <Leibniz rule> for the <exterior derivative> give
$$
0=\int_Md(\beta\wedge*\alpha)=\int_Md\beta\wedge*\alpha+(-1)^{p-1}\int_M\beta\wedge d*\alpha.
$$
Thus $(d\beta,\alpha)_{L^2}=(-1)^p\int_M\beta\wedge d*\alpha$. The formal <adjoint operator> $\delta$ satisfies $*\delta\alpha=(-1)^p d*\alpha$. Applying the inverse star and its square on degree $p-1$ gives the <codifferential>
$$
\boxed{\delta|_{\Omega^p}=(-1)^{n(p+1)+1}*d*.}
$$
The exponent is congruent to $p+(p-1)(n-p+1)$ modulo two, so this convention is consistent in every degree. Define the <Hodge Laplacian> by $\Delta=d\delta+\delta d$ and a <harmonic differential form> by $\Delta\alpha=0$. On a compact manifold,
$$
(\Delta\alpha,\alpha)_{L^2}=\|d\alpha\|_{L^2}^2+\|\delta\alpha\|_{L^2}^2.
$$
Both nonnegative terms vanish for a harmonic form, giving \b[$d\alpha=\delta\alpha=0$]. Conversely these two equations imply harmonicity.

Now let $n=4$. The printed eigenvalue relation omits the subscript on its right-hand side; the intended equation is $*\alpha_\pm=\pm\alpha_\pm$. On two-forms $*^2=1$, so
$$
\boxed{\alpha_+=\tfrac12(\alpha+*\alpha),\qquad\alpha_-=\tfrac12(\alpha-*\alpha).}
$$
These are respectively <self-dual two-forms> and <anti-self-dual two-forms>, and sum to $\alpha$. Their uniqueness follows because the $+1$ and $-1$ eigenspaces have zero intersection. Star is an orthogonal involution here, hence is self-adjoint and its opposite eigenspaces are orthogonal.

For closed two-forms define $Q([\alpha],[\beta])=\int_N\alpha\wedge\beta$. Adding $d\xi$ to $\alpha$ changes the integral by $\int_Nd(\xi\wedge\beta)=0$; adding an exact form to $\beta$ has the same effect. This is the <real de Rham intersection form in dimension four>. The <wedge product of differential forms> in degrees two and two is symmetric because $(-1)^{2\cdot2}=1$; hence $Q$ is a symmetric bilinear form on $H_{\mathrm{dR}}^2(N;\mathbb R)$.

In four dimensions $\delta=-*d*$ in every degree. If $h$ is a harmonic two-form, $\delta h=0$ gives $d*h=0$, and $dh=0$ gives $\delta*h=-*dh=0$. Thus $*h$ is harmonic, and $\mathcal H^2=\mathcal H^+\oplus\mathcal H^-$ is the harmonic star-eigenspace decomposition. Use the permitted unique harmonic representative for each <de Rham cohomology> class. For harmonic $h,k$,
$$
Q(h,k)=(h,*k)_{L^2}= (h_+,k_+)_{L^2}-(h_-,k_-)_{L^2}.
$$
The mixed terms vanish by orthogonality. The first restriction is positive definite and the second negative definite. In particular a nonzero harmonic $h$ has $Q(h,*h)=\|h\|_{L^2}^2>0$, proving nondegeneracy. Hence the <signature of the intersection form from harmonic duality> is
$$
\boxed{\operatorname{signature}Q=(b^+(N),b^-(N)),\qquad b^\pm(N)=\dim\mathcal H^\pm(N).}
$$
The finite dimensions and the harmonic description also follow from the <Hodge decomposition theorem>.

For completeness, the <Hodge decomposition theorem> on a compact oriented <Riemannian manifold> without boundary states the finite-dimensional harmonic space and the $L^2$-orthogonal direct sum
$$
\boxed{\Omega^p=\mathcal H^p\oplus d\Omega^{p-1}\oplus\delta\Omega^{p+1}.}
$$
The three summand components are unique, although their potentials need not be. Given an exact three-form $\beta=d\omega$ on $N$, decompose its smooth two-form potential as $\omega=h+d\xi+\delta\psi$, with $\psi\in\Omega^3$. Since $dh=0$ and $d^2=0$, $\beta=d\delta\psi$. Put $a=\delta\psi$. In dimension four, $a=-*d*\psi$, and $*^2=1$ on two-forms gives $*a=-d*\psi$, so $d*a=0$. Consequently
$$
\boxed{\sigma=a+*a=\delta\psi+*\delta\psi,\qquad *\sigma=\sigma,\qquad d\sigma=\beta.}
$$
This constructs a <self-dual primitive of an exact three-form>. The sum contains no factor $1/2$: it is twice the self-dual projection of the coexact potential. Using the unscaled projection alone would give $d\sigma=\beta/2$. This proves the final unheaded request.