I'm trying to derive the incompressible Navier-Stokes equation in Fourier space, given in e.g. Kraichnan & Montgomery (1980), Rep. Prog. Phys. 43 547 (PDF) (but probably countless other places), as $$ (\partial_t + \nu k^2)u_i(\mathbf{k}) = -\mathrm{i}\left(\delta_{ij}- \frac{k_i k_j}{k^2}\right) k_m \sum_{\mathbf{k}=\mathbf{k}_1+\mathbf{k}_2} u_j(\mathbf{k}_1) u_m(\mathbf{k}_2) \tag{0}\label{0} $$ where roman indices are tensor components, repeated indices are summed, $\mathbf{k}_1$ and $\mathbf{k}_2$ are dummy wavevectors in the sum over resonant wave triads, and I have omitted hats on the Fourier variables (i.e. $u_i(\mathbf{k}) = \hat{u}_i(\mathbf{k})= \int u_i(\mathbf{x})e^{-\mathrm{i}\mathbf{k}\cdot\mathbf{x}}d^dx$, modulo some normalising constants).
The left-hand side of $\eqref{0}$ is easy enough but I'm having trouble making the right-hand side nice and symmetric:
Navier-Stokes (with density normalised so that $\rho=1$) is $$ \partial_t\mathbf{u} + (\mathbf{u}\cdot\nabla)\mathbf{u}=-\nabla p + \nu\nabla^2\mathbf{u} \tag{1}\label{1} $$ and incompressibility ($\nabla\cdot\mathbf{u}=0$) gives for the pressure $$ \nabla^2p = -\nabla\cdot[(\mathbf{u}\cdot\nabla)\mathbf{u}] .\tag{2}\label{2} $$
I put $\eqref{2}$ in index notation and write $p, \mathbf{u}$ in Fourier series, e.g. $u_i(\mathbf{x})=\sum_\mathbf{k'}u_i(\mathbf{k'})e^{\mathrm{i}\mathbf{k'}\cdot\mathbf{x}}$. I then multiply by $e^{-\mathrm{i}\mathbf{k}\cdot\mathbf{x}}$, integrate over space and use $\int e^{\mathrm{i}(\mathbf{k'}-\mathbf{k})\cdot\mathbf{x}}d^dx=\delta_{\mathbf{k'}\mathbf{k}}$(modulo constants) to get $$ p(\mathbf{k})=-\frac{1}{k^2}k_j \sum_{\mathbf{k}=\mathbf{k}_1+\mathbf{k}_2} k_{1m} u_j(\mathbf{k}_1)u_m(\mathbf{k}_2). \tag{3}\label{3} $$ The $k_j$ on the RHS of $\eqref{3}$ comes from the outer del $\nabla\cdot[\ldots]$ on the RHS of $\eqref{2}$, and from the $j$-th component of the resonance condition $\mathbf{k}=\mathbf{k}_1+\mathbf{k}_2$.
Now I write $\eqref{1}$ in index notation and put the variables in Fourier series, using $\eqref{3}$ to eliminate the pressure. Using all the same tricks I get $$ (\partial_t + \nu k^2)u_i(\mathbf{k}) = -\mathrm{i}\left(\delta_{ij}- \frac{k_i k_j}{k^2}\right) \sum_{\mathbf{k}=\mathbf{k}_1+\mathbf{k}_2} k_{1m} u_j(\mathbf{k}_1) u_m(\mathbf{k}_2) \tag{4}\label{4}. $$ (Again, roman indices are components, numeric indices represent dummy wavevectors in the sum over wave triads).
The RHS of my eq. $\eqref{4}$ is similar to, but not the same as the desired RHS $\eqref{0}$. Specifically I just need to show that $$ P_{ij} \sum_{\mathbf{k}=\mathbf{k}_1+\mathbf{k}_2} k_{1m} u_j(\mathbf{k}_1) u_m(\mathbf{k}_2) = P_{ij}k_m\sum_{\mathbf{k} = \mathbf{k}_1+\mathbf{k}_2} u_j(\mathbf{k}_1) u_m(\mathbf{k}_2) $$ where $P_{ij}=\left(\delta_{ij}- \frac{k_i k_j}{k^2}\right)$. So I need to turn the $k_{1m}$ inside the sum into a $k_m=k_{1m}+k_{2m}$ that can be brought outside the sum. Now I believe the trick is to do some index manipulation with the dummy indices $j,m$ and the dummy wavevectors $\mathbf{k}_1,\mathbf{k}_2$. But I don't see it.
Any advice?