 Bogomolny square completion for an Abelian Higgs vortex

With $D_i={\partial }_i-i{A}_i$ and $B={\partial }_1{A}_2-{\partial }_2{A}_1$, use the critical-coupling energy $E=\int [B^2/2+|D_i\varphi {|}^2/2+(1-|\varphi {|}^2{)}^2/8]d^{2}x$. For positive vortex number, integration by parts givesThe boundary divergence vanishes for the usual decaying vortex fields, and magnetic flux is $2\pi N$. The two squares vanish exactly at the Bogomolny vortex equations. Reversing both signs treats negative $N$. Energy coefficients and the numerical bound change together under other normalizations.