Editing Helmholtz decomposition (section)

=== Tensor approach ===

In a <math>d</math>-dimensional vector space with <math>d\neq 3</math>, <math display="inline">-\frac{1}{4\pi\left|\mathbf{r}-\mathbf{r}'\right|}</math> can be replaced by the appropriate [[Green's function#Green's functions for the Laplacian|Green's function for the Laplacian]], defined by
<math display="block">
\nabla^2 G(\mathbf{r},\mathbf{r}') = \frac{\partial}{\partial r_\mu}\frac{\partial}{\partial r_\mu}G(\mathbf{r},\mathbf{r}') = \delta^d(\mathbf{r}-\mathbf{r}')
</math>
where [[Einstein notation|Einstein summation convention]] is used for the index <math>\mu</math>. For example, <math display="inline">G(\mathbf{r},\mathbf{r}')=\frac{1}{2\pi}\ln\left|\mathbf{r}-\mathbf{r}'\right|</math> in 2D.

Following the same steps as above, we can write
<math display="block">
F_\mu(\mathbf{r}) = \int_V F_\mu(\mathbf{r}') \frac{\partial}{\partial r_\mu}\frac{\partial}{\partial r_\mu}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
 = \delta_{\mu\nu}\delta_{\rho\sigma}\int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\rho}\frac{\partial}{\partial r_\sigma}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
</math>
where <math>\delta_{\mu\nu}</math> is the [[Kronecker delta]] (and the summation convention is again used). In place of the definition of the vector Laplacian used above, we now make use of an identity for the [[Levi-Civita symbol]] <math>\varepsilon</math>,
<math display="block">
\varepsilon_{\alpha\mu\rho}\varepsilon_{\alpha\nu\sigma} = (d-2)!(\delta_{\mu\nu}\delta_{\rho\sigma} - \delta_{\mu\sigma}\delta_{\nu\rho})
</math>
which is valid in <math>d\ge 2</math> dimensions, where <math>\alpha</math> is a <math>(d-2)</math>-component [[Multi-index notation|multi-index]]. This gives
<math display="block">
F_\mu(\mathbf{r}) = \delta_{\mu\sigma}\delta_{\nu\rho}\int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\rho}\frac{\partial}{\partial r_\sigma}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
+ \frac{1}{(d-2)!}\varepsilon_{\alpha\mu\rho}\varepsilon_{\alpha\nu\sigma} \int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\rho}\frac{\partial}{\partial r_\sigma}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
</math>

We can therefore write
<math display="block">
F_\mu(\mathbf{r}) = -\frac{\partial}{\partial r_\mu} \Phi(\mathbf{r}) + \varepsilon_{\mu\rho\alpha}\frac{\partial}{\partial r_\rho} A_{\alpha}(\mathbf{r})
</math>
where
<math display="block">
\begin{aligned}
\Phi(\mathbf{r}) &= -\int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\nu}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'\\
A_{\alpha} &= \frac{1}{(d-2)!}\varepsilon_{\alpha\nu\sigma} \int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\sigma}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
\end{aligned}
</math>
Note that the vector potential is replaced by a rank-<math>(d-2)</math> tensor in <math>d</math> dimensions.

Because <math>G(\mathbf{r},\mathbf{r}')</math> is a function of only <math>\mathbf{r}-\mathbf{r}'</math>, one can replace <math>\frac{\partial}{\partial r_\mu}\rightarrow - \frac{\partial}{\partial r'_\mu}</math>, giving
<math display="block">
\begin{aligned}
\Phi(\mathbf{r}) &= \int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r'_\nu}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'\\
A_{\alpha} &= -\frac{1}{(d-2)!}\varepsilon_{\alpha\nu\sigma} \int_V F_\nu(\mathbf{r}') \frac{\partial}{\partial r_\sigma'}G(\mathbf{r},\mathbf{r}') \,\mathrm{d}^d \mathbf{r}'
\end{aligned}
</math>
[[Integration_by_parts#Higher_dimensions|Integration by parts]] can then be used to give
<math display="block">
\begin{aligned}
\Phi(\mathbf{r}) &= -\int_V G(\mathbf{r},\mathbf{r}')\frac{\partial}{\partial r'_\nu}F_\nu(\mathbf{r}') \,\mathrm{d}^d \mathbf{r}' + \oint_{S} G(\mathbf{r},\mathbf{r}') F_\nu(\mathbf{r}') \hat{n}'_\nu \,\mathrm{d}^{d-1} \mathbf{r}'\\
A_{\alpha} &= \frac{1}{(d-2)!}\varepsilon_{\alpha\nu\sigma} \int_V G(\mathbf{r},\mathbf{r}') \frac{\partial}{\partial r_\sigma'}F_\nu(\mathbf{r}')  \,\mathrm{d}^d \mathbf{r}'- \frac{1}{(d-2)!}\varepsilon_{\alpha\nu\sigma} \oint_{S} G(\mathbf{r},\mathbf{r}') F_\nu(\mathbf{r}')  \hat{n}'_\sigma \,\mathrm{d}^{d-1} \mathbf{r}'
\end{aligned}
</math>
where <math>S=\partial V</math> is the boundary of <math>V</math>. These expressions are analogous to those given above for [[#Three-dimensional_space|three-dimensional space]].

For a further generalization to manifolds, see the discussion of [[Hodge decomposition]] [[#Differential forms|below]].