[step:Equip forms with the $L^2$ inner product and record the Laplacian identity]
Let $n=\dim M$, and let $\operatorname{vol}_g$ denote the Riemannian volume measure on $M$. We use the endpoint convention that $\Omega^r(M)=\{0\}$ for every integer $r<0$ or $r>n$, and that $d$ and $\delta_g$ are the corresponding zero maps whenever their source or target lies outside degrees $0,\dots,n$. For each integer $r$ with $0 \leq r \leq n$, let
\begin{align*}
(\cdot,\cdot)_{L^2}: \Omega^r(M) \times \Omega^r(M) &\to \mathbb{R}
\end{align*}
be the $L^2$ inner product defined by
\begin{align*}
(\alpha,\beta)_{L^2}
=
\int_M \langle \alpha(p),\beta(p)\rangle_{\Lambda^r T_p^*M,g}\,d\operatorname{vol}_g(p),
\end{align*}
where $\langle \cdot,\cdot\rangle_{\Lambda^r T_p^*M,g}$ is the inner product induced by $g_p$ on $\Lambda^r T_p^*M$.
The codifferential $\delta_g: \Omega^{r+1}(M) \to \Omega^r(M)$ is the formal $L^2$ adjoint of $d: \Omega^r(M) \to \Omega^{r+1}(M)$, so for all $\alpha \in \Omega^r(M)$ and $\beta \in \Omega^{r+1}(M)$,
\begin{align*}
(d\alpha,\beta)_{L^2} = (\alpha,\delta_g\beta)_{L^2}.
\end{align*}
Since $M$ is closed, there are no boundary terms in this adjoint relation.
For any $\omega \in \Omega^k(M)$, using $\Delta_g=d\delta_g+\delta_gd$ and the adjoint relation gives
\begin{align*}
(\Delta_g\omega,\omega)_{L^2}
&=
(d\delta_g\omega,\omega)_{L^2}+(\delta_gd\omega,\omega)_{L^2} \\
&=
(\delta_g\omega,\delta_g\omega)_{L^2}+(d\omega,d\omega)_{L^2} \\
&=
\|\delta_g\omega\|_{L^2}^2+\|d\omega\|_{L^2}^2.
\end{align*}
[/step]