Optimal shape and regularity of magnetic needle domains

We study the nonlocal energy $$\begin{array}{r}\mathcal{E}(\mathrm{\Omega })=\mathcal{P}(\mathrm{\Omega })+{\int }_{{\mathbb{R}}^n}{\left|\mathrm{\nabla }{\mathrm{\Phi }}_{\mathrm{\Omega }}\right|}^2\text{rm d}x\end{array}$$ where $\mathcal{P}(\mathrm{\Omega })$ denotes the perimeter of a set $\mathrm{\Omega }\subset {\mathbb{R}}^n$ with prescribed volume $|\mathrm{\Omega }|=V>0$ and ${\mathrm{\Phi }}_{\mathrm{\Omega }}\in {\dot{H}}^1({\mathbb{R}}^n)$ is the weak solution of $\mathrm{\Delta }{\mathrm{\Phi }}_{\mathrm{\Omega }}={\partial }_1{\chi }_{\mathrm{\Omega }}$. This energy is related to the study of long and slender "needle" domains in uniaxial ferromagnetic materials where the second term models dipole interactions. We prove existence of minimizers and show that they are (${\mathcal{L}}^n$ equivalent to) bounded, connected sets with smooth boundary. In particular we show that $\mathrm{\nabla }{\mathrm{\Phi }}_{\mathrm{\Omega }}\in L^{\mathrm{\infty }}$ for local minimizers. For the physically important case $n=3$, we furthermore establish a scaling law for the minimal energy in terms of the prescribed mass $V$.