[proof]
Set $w(t) := (1-\Delta)^{s/2}u(t) \in L^2(\mathbb{R}^n; \mathbb{R}^n)$. By the Plancherel identity derived above, $\|u(t)\|_{H^s}^2 = \|w(t)\|_{L^2}^2$. Since $(1-\Delta)^{s/2}$ acts only in the spatial variable, it commutes with $\partial_t$, giving $\partial_t w = (1-\Delta)^{s/2}(\partial_t u)$. The chain rule $\frac{d}{dt}\|w\|_{L^2}^2 = 2\int w \cdot \partial_t w\,d\mathcal{L}^n$ therefore gives
\begin{align*}
\frac{1}{2}\frac{d}{dt}\|u\|_{H^s}^2 &= \int_{\mathbb{R}^n} (1-\Delta)^{s/2} u \cdot (1-\Delta)^{s/2}(\partial_t u)\,d\mathcal{L}^n.
\end{align*}
Substituting $\partial_t u = -\mathbb{P}[(u \cdot \nabla)u]$ from Step 1:
\begin{align*}
&= -\int_{\mathbb{R}^n} (1-\Delta)^{s/2} u \cdot (1-\Delta)^{s/2}\bigl(\mathbb{P}[(u \cdot \nabla)u]\bigr)\,d\mathcal{L}^n.
\end{align*}
The Leray projector $\mathbb{P}$ can now be removed. Three facts are needed: (i) $\mathbb{P}$ is self-adjoint on $L^2(\mathbb{R}^n; \mathbb{R}^n)$ (from the [Helmholtz Decomposition](/theorems/654)); (ii) $\mathbb{P}$ commutes with $(1-\Delta)^{s/2}$ — both are Fourier multipliers, one acting component-wise as multiplication by $\delta_{ij} - \xi_i\xi_j/|\xi|^2$ and the other by the scalar $(1+|\xi|^2)^{s/2}$, and scalar multiplication commutes with matrix multiplication; (iii) $\mathbb{P}((1-\Delta)^{s/2}u) = (1-\Delta)^{s/2}u$ — since $u$ is divergence-free ($\xi \cdot \hat{u}_{\mathrm{rep}}(\xi) = 0$ for $\mathcal{L}^n$-a.e. $\xi$), and multiplying each component by the scalar $(1+|\xi|^2)^{s/2}$ preserves this condition, $(1-\Delta)^{s/2}u$ is also divergence-free, so $\mathbb{P}$ acts as the identity. Using self-adjointness of $\mathbb{P}$ to move it from the right factor to the left, then applying (iii):
\begin{align*}
&= -\int_{\mathbb{R}^n} (1-\Delta)^{s/2} u \cdot (1-\Delta)^{s/2}\bigl[(u \cdot \nabla)u\bigr]\,d\mathcal{L}^n.
\end{align*}
Now expand $(1-\Delta)^{s/2}[(u \cdot \nabla)u]$ using the commutator identity $AB = BA + [A,B]$ with $A = (1-\Delta)^{s/2}$ and $B = u \cdot \nabla$:
\begin{align*}
(1-\Delta)^{s/2}\bigl[(u \cdot \nabla)u\bigr] &= (u \cdot \nabla)(1-\Delta)^{s/2}u + \bigl[(1-\Delta)^{s/2},\, u \cdot \nabla\bigr]u.
\end{align*}
Substituting produces two integrals:
\begin{align*}
\frac{1}{2}\frac{d}{dt}\|u\|_{H^s}^2 &= -\underbrace{\int_{\mathbb{R}^n} (1-\Delta)^{s/2} u \cdot (u \cdot \nabla)(1-\Delta)^{s/2}u\,d\mathcal{L}^n}_{(I)} - \underbrace{\int_{\mathbb{R}^n} (1-\Delta)^{s/2} u \cdot \bigl[(1-\Delta)^{s/2}, u \cdot \nabla\bigr]u\,d\mathcal{L}^n}_{(II)}.
\end{align*}
It remains to show $(I) = 0$. Recalling $w = (1-\Delta)^{s/2}u$ and expanding component by component:
\begin{align*}
(I) &= \sum_{i=1}^n \sum_{j=1}^n \int_{\mathbb{R}^n} w_i\, u_j\, \partial_j w_i\,d\mathcal{L}^n.
\end{align*}
Using $w_i\,\partial_j w_i = \frac{1}{2}\partial_j(w_i^2)$ and summing over $i$:
\begin{align*}
(I) &= \frac{1}{2}\sum_{j=1}^n \int_{\mathbb{R}^n} u_j\,\partial_j(|w|^2)\,d\mathcal{L}^n = \frac{1}{2}\int_{\mathbb{R}^n} u \cdot \nabla(|w|^2)\,d\mathcal{L}^n.
\end{align*}
Integration by parts moves $\nabla$ from $|w|^2$ to $u$ (justified since $u, w \in H^s \hookrightarrow L^\infty$ for $s > n/2$ by the [Sobolev Embedding Theorem](/theorems/903), so all terms are in $L^1(\mathbb{R}^n)$):
\begin{align*}
(I) &= -\frac{1}{2}\int_{\mathbb{R}^n} |w|^2\,(\nabla \cdot u)\,d\mathcal{L}^n = 0,
\end{align*}
since $\nabla \cdot u = 0$. Therefore only $(II)$ survives, giving the claimed identity.
[/proof]