Vortices, Flux Quantization, and Core Scales

In the Abelian Higgs model, scalar winding and a compensating gauge field produce a smooth vortex with quantized magnetic flux and finite tension. The gauge field cancels the long-range phase gradient, unlike a global vortex whose tension grows logarithmically with the transverse system size. Two independent inverse masses set the scalar and magnetic core scales, and their ratio controls whether vortices attract, repel, or satisfy first-order equations.

Required background. Finite-energy boundary data supplies the transverse-circle classification, while gauge fields and observable content fixes the gauge redundancy and covariant derivative. Helpful background. Quantized topological terms and global consistency helps distinguish flux quantization from a merely local field equation.

Abelian Higgs data and the vortex ansatz

Work in $2+1$ dimensions, or per unit length of a straight static string in $3+1$ dimensions. With $D_\mu=\partial_\mu-ieA_\mu$, take

$$\mathcal L =-\frac14F_{\mu\nu}F^{\mu\nu} +(D_\mu\phi)^*D^\mu\phi -\frac{\lambda}{2}\left(|\phi|^2-v^2\right)^2 ,$$

where $e,\lambda,v>0$ and the minimally charged field has charge $e$. For a rotationally symmetric configuration of nonzero winding $n\in\mathbb Z\setminus\{0\}$,

$$\phi(r,\theta)=v f(r)e^{in\theta}, \qquad A=\frac{n}{e}a(r)\,\mathrm d\theta .$$

regularity at the origin and finite energy at infinity require

$$f(0)=0,\qquad a(0)=0, \qquad f(\infty)=1,\qquad a(\infty)=1 .$$

The $n=0$ vacuum sector is different: the scalar need not vanish at the origin, and $a(r)$ drops out of this particular parametrization because $A=0$.

The angular covariant derivative is

$$D_\theta\phi =in\big(1-a(r)\big)\phi,$$

so the gauge field cancels the scalar phase gradient asymptotically. The tension is

$$\begin{aligned} T=2\pi\int_0^\infty r\,\mathrm dr\, \Bigg[ &v^2(f')^2 +\frac{n^2v^2}{r^2}(1-a)^2f^2\\ &+\frac{n^2}{2e^2r^2}(a')^2 +\frac{\lambda v^4}{2}(f^2-1)^2 \Bigg]. \end{aligned}$$

Every term is finite with the stated boundary conditions. The radial equations obtained by varying this functional determine $f$ and $a$; the boundary data alone do not.

Flux quantization and orientation

Orient the plane by $\mathrm dx^1\wedge\mathrm dx^2=r\,\mathrm dr\wedge\mathrm d\theta>0$. Then

$$B=F_{12}=\frac{F_{r\theta}}{r} =\frac{n}{er}a'(r).$$

Thus an increasing profile $a:0\to1$ has positive magnetic field and flux when $n>0$. The magnetic two-form is $F=\mathrm dA$, and Stokes’ theorem gives

$$\Phi_B =\int_{\mathbb R^2}F =\oint_{S^1_\infty}A =\frac{2\pi n}{e}.$$

Equivalently, $D_\theta\phi\to0$ requires $eA_\theta\to\partial_\theta\arg\phi=n$. The integer belongs to the phase map on the circle at infinity, while the physical flux quantum also depends on the minimal electric charge and the global gauge group. A field redefinition that moves $e$ from $D_\mu$ into the gauge kinetic term changes the displayed formula but not the Aharonov–Bohm phase factor $\exp(iq\Phi_B)$ of an allowed probe, or equivalently the angle $q\Phi_B$ modulo $2\pi$.

The original relativistic construction is Nielsen and Olesen 1973, §§ 2–3, pp. 47–54. Modern treatments of the topology, profiles, forces, and moduli include Manton and Sutcliffe 2004, ch. 7, pp. 158–240 and Tong 2005, lecture 3.

Scalar and magnetic core scales

Linearizing around $|\phi|=v$ in this normalization gives

$$m_{\mathrm H}^2=2\lambda v^2, \qquad m_{\mathrm A}^2=2e^2v^2 .$$

Thus the scalar modulus approaches its vacuum over $\xi_{\mathrm H}\sim m_{\mathrm H}^{-1}$, while magnetic flux spreads over $\xi_{\mathrm A}\sim m_{\mathrm A}^{-1}$. The dimensionless ratio

$$\beta=\frac{m_{\mathrm H}^2}{m_{\mathrm A}^2} =\frac{\lambda}{e^2}$$

compares the cores. At $\beta=1$ the energy admits the critical Bogomolny completion and

$$T=2\pi v^2|n|$$

for a saturated solution. For $\beta>1$, the same rearrangement leaves the additional nonnegative term $(\lambda-e^2)(|\phi|^2-v^2)^2/2$, so $T\geq2\pi v^2|n|$ still holds but is strict for a nontrivial vortex. For $\beta<1$, this square completion does not provide that lower bound.

For two well-separated vortices of the same sign in this single-complex-scalar model, the longer-ranged scalar tail gives attraction for $\beta<1$, the longer-ranged vector tail gives repulsion for $\beta>1$, and the forces cancel at $\beta=1$ Jacobs and Rebbi 1979, pp. 4486–4494. This is the relativistic type-I/type-II distinction in the normalization used here. In the common Ginzburg–Landau convention,

$$\kappa_{\mathrm{GL}}=\sqrt{\frac{\beta}{2}},$$

so $\kappa_{\mathrm{GL}}=1/\sqrt2$ is the same critical point as $\beta=1$. The force statement does not describe a vortex–antivortex pair. It also assumes infinite flat transverse space; additional charged fields, Chern–Simons terms, boundaries, or non-Abelian structure can change both moduli and forces.

The equations have useful near-core checks. Regularity gives

$$f(r)\sim c_f r^{|n|}, \qquad a(r)\sim c_a r^2 \quad (r\to0),$$

while the far tails solve massive linear equations and decay exponentially with the two masses. A numerical profile that uses one fitted length for both tails away from $\beta=1$ has lost physical information.

Why a global vortex is different

Set $e=0$ while retaining a spontaneously broken global $U(1)$. Far outside the scalar core, $\phi\simeq ve^{in\theta}$ and

$$|\boldsymbol\nabla\phi|^2 \simeq\frac{n^2v^2}{r^2}.$$

The tension between a core radius $\xi$ and an infrared radius $R$ is therefore

$$T_{\mathrm{global}} \simeq 2\pi n^2v^2\log\!\frac{R}{\xi} +T_{\mathrm{core}}.$$

It is finite in a finite container but diverges logarithmically as $R\to\infty$. Calling it a finite-tension local vortex suppresses the order of limits and the massless Goldstone tail. In the gauged theory, $A_\theta$ cancels the phase gradient and the remaining fields are massive, so the infinite-plane tension is finite.

Boundary map, existence, and stability

The shared boundary-family map emphasizes the transverse circle and its relation to other codimensions: for the Abelian Higgs vortex, finite energy ties scalar winding on $S^1_\infty$ to magnetic flux. The soliton boundary and stability comparison keeps that classification separate from profile existence, coupling-dependent forces, positional moduli, and quantum scope.

The integer flux sector prevents a smooth finite-energy path to the vacuum when the allowed boundary conditions and charge lattice are held fixed. It does not ensure that an axially symmetric $|n|>1$ profile is the energy minimum: away from critical coupling it may prefer separated unit vortices or a bound multivortex. Energetic stability must be stated for the given $n$, coupling ratio, geometry, and allowed perturbations.

At critical coupling, the first-order vortex equations imply the second-order field equations and yield a $2|n|$-dimensional real moduli space, interpreted as the transverse positions of $|n|$ vortices Weinberg 1979, pp. 3008–3012. Away from critical coupling, those static moduli are generally lifted by forces. Exact non-Abelian vortex moduli and their supersymmetric dynamics require additional field content and belong to the supersymmetry treatment rather than this Abelian model.

Checks and failure modes

Flux. Compute $\Phi_B$ both from the area integral of $B$ and from the asymptotic line integral. Disagreement diagnoses a polar-coordinate or gauge-patch error.

Dimensions. In $2+1$ dimensions, $[v]=[e]=1/2$ and $[\lambda]=1$, so $\beta$ is dimensionless and $2\pi v^2|n|$ has the dimension of energy. In $3+1$ dimensions, $[v]=1$ while $e$ and $\lambda$ are dimensionless, so the same expression is a string tension of mass dimension two.

Tail hierarchy. Fit scalar and magnetic tails separately. The fitted masses should agree with the vacuum spectrum before a core-scale claim is trusted.

Infrared order of limits. For a global vortex, state $R$ and take the infinite-volume limit explicitly. For a local vortex, verify that covariant rather than ordinary gradients vanish.

Common pitfalls

Writing $A_\theta\to0$ in the winding gauge. In the one-form convention used here, $A_\theta\to n/e$. The gauge-invariant requirement is $D_\theta\phi\to0$; another gauge may move the winding and the asymptotic potential together.

Equating winding with flux without declaring the charge lattice. The integer winding fixes $e\Phi_B/2\pi$ for the minimally charged Higgs field in this model. A different global form or minimally allowed charge changes which fluxes are distinct.

Calling the global vortex tension finite. Its Goldstone gradient produces a logarithm in the transverse infrared. A finite numerical box regulates rather than removes it.

Exercises

1. Starting from the ansatz, compute $F$ and verify $\Phi_B=2\pi n/e$.

Solution

Since $A=(n/e)a(r)\mathrm d\theta$,

$$F=\frac{n}{e}a'(r)\,\mathrm dr\wedge\mathrm d\theta.$$

Therefore

$$\Phi_B =\frac{n}{e}\int_0^\infty\mathrm dr\,a'(r) \int_0^{2\pi}\mathrm d\theta =\frac{2\pi n}{e}[a(\infty)-a(0)] =\frac{2\pi n}{e}.$$

2. Show that the global-vortex angular gradient gives a logarithmic tension and identify the ultraviolet and infrared regulators.

Solution

Outside the core, the angular energy is $n^2v^2/r^2$. Integrating with $\mathrm d^2x=2\pi r\,\mathrm dr$ gives

$$2\pi n^2v^2\int_\xi^R\frac{\mathrm dr}{r} =2\pi n^2v^2\log(R/\xi).$$

The core size $\xi$ regulates the short-distance approximation, and the system size or inter-vortex separation $R$ regulates the infrared Goldstone field.

3. At critical coupling and positive winding, derive the radial first-order equations from

$$(D_1+iD_2)\phi=0, \qquad B=e(v^2-|\phi|^2).$$

Solution

With the chosen orientation,

$$D_1+iD_2 =e^{i\theta}\left(D_r+\frac{i}{r}D_\theta\right), \qquad B=\frac{n}{er}a'.$$

Substituting the ansatz gives

$$f'=\frac{n}{r}(1-a)f, \qquad a'=\frac{e^2v^2r}{n}(1-f^2), \qquad n>0.$$

Both right-hand sides are positive near the core for a regular positive-flux solution. For $n<0$, use the sign-reversed Bogomolny equations; writing them in terms of $|n|$ gives the corresponding monotone radial profiles.

Continue

Bogomolny Bounds and First-Order Equations derives the critical-coupling equations and bound. Moduli-Space Dynamics and Collective Quantization states when vortex positions can be treated as low-energy coordinates.

References

• Jacobs, Laurence, and Claudio Rebbi. “Interaction Energy of Superconducting Vortices.” Physical Review B 19 (1979): 4486–4494. DOI.
• Manton, Nicholas, and Paul Sutcliffe. Topological Solitons. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 2004, ch. 7, pp. 158–240. DOI.
• Nielsen, Holger Bech, and Poul Olesen. “Vortex-Line Models for Dual Strings.” Nuclear Physics B 61 (1973): 45–61. DOI.
• Tong, David. “TASI Lectures on Solitons: Instantons, Monopoles, Vortices and Kinks.” 2005, lecture 3. arXiv:hep-th/0509216.
• Weinberg, Erick J. “Multivortex Solutions of the Ginzburg–Landau Equations.” Physical Review D 19 (1979): 3008–3012. DOI.