Gravitating Abelian-Higgs Cosmic Strings

• Gravitating Abelian-Higgs cosmic strings are finite-width field configurations emerging from a broken U(1) symmetry coupled with gravity, exhibiting a localized vortex core and conical asymptotics.
• They display diverse behaviors across gravity models—from standard Einstein to modified theories like Starobinsky and Rastall—affecting the vortex structure, angular deficit, and dynamic evolution.
• Recent studies highlight their roles in gravitational collapse, wave emission, lensing phenomena, and network cosmology, expanding their relevance in both astrophysics and theoretical models.

Searching arXiv for recent and foundational papers on gravitating Abelian-Higgs cosmic strings. Gravitating Abelian–Higgs cosmic strings are line-like, finite-width field configurations of a spontaneously broken local $U(1)$ gauge theory coupled to gravity, whose matter sector is the Nielsen–Olesen vortex or one of its close generalizations. In the minimal setting, the fields are a complex Higgs scalar and an Abelian gauge field with quantized winding and magnetic flux, while gravity is incorporated either through the full Einstein equations, through modified-gravity extensions such as $f(R)$ or Rastall gravity, or, in some dynamical contexts, through tensor perturbations sourced by the string stress tensor. Across these settings, the central questions are the structure of the vortex core, the asymptotic geometry—often conical in the four-dimensional Einstein case—the existence of self-dual/BPS sectors, the role of higher-dimensional or sigma-model generalizations, and the gravitational signatures of strings ranging from deficit angles and geodesic effects to gravitational-wave emission (Graça, 2015, Slagter et al., 2012, Helfer et al., 2018, Silva et al., 3 Jun 2026).

1. Field-theoretic definition and vortex structure

The standard gravitating Abelian–Higgs string is built from a complex scalar Higgs field and a $U(1)$ gauge field with spontaneous symmetry breaking. In the Starobinsky study, the full action is written as

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$

with

$f(R)=R+\eta R^2,$

$$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$

(Graça, 2015). In the full general-relativistic loop-collapse study, the Einstein–Abelian-Higgs action is

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$

with

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$

(Helfer et al., 2018).

For straight strings, the Nielsen–Olesen ansatz appears in several formulations. In the four-dimensional Starobinsky model one uses

$$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$

(Graça, 2015). In the warped five-dimensional braneworld analysis, the brane vortex ansatz is

$f(R)$0

(Slagter et al., 2012). In the optical study of finite-width gravitating strings, the static string along the $f(R)$1-axis is described by

$f(R)$2

$f(R)$3

with the paper restricting to $f(R)$4 (Silva et al., 3 Jun 2026).

These ansätze encode the standard vortex mechanism. The scalar amplitude vanishes at the axis, the phase winds by $f(R)$5, and the gauge field screens the angular gradient so that finite energy per unit length is obtained. The finite-energy conditions are correspondingly

$f(R)$6

for the unit-winding case (Graça, 2015). In the preheating study, this topological structure is stated explicitly in the form

$f(R)$7

with flux quantization

$f(R)$8

(Dufaux et al., 2010).

A common control parameter is the ratio of gauge to scalar mass scales. In the preheating analysis,

$f(R)$9

(Dufaux et al., 2010). In the Starobinsky model, the dimensionless parameter is

$U(1)$0

(Graça, 2015). In the optical analysis,

$U(1)$1

(Silva et al., 3 Jun 2026). This suggests that the detailed width and local gravitational response of the vortex are governed by the same scalar–vector mass hierarchy that controls type-I/type-II behavior in the flat-space theory.

2. Einstein gravity, conical asymptotics, and self-duality

In four-dimensional Einstein gravity, the canonical picture is that a straight Abelian–Higgs string has a localized core and an asymptotically conical exterior. The Starobinsky paper makes this reference point explicit by noting that asymptotically one has

$U(1)$2

with angular deficit

$U(1)$3

(Graça, 2015). The optical paper uses the same asymptotic structure,

$U(1)$4

and identifies the deficit angle as

$U(1)$5

(Silva et al., 3 Jun 2026). In the analysis of supermassive strings, the deficit angle is written as

$U(1)$6

and, more generally,

$U(1)$7

(Blanco-Pillado et al., 2013).

At critical coupling, the gravitating system admits a Bogomolny reduction. In the supermassive-string analysis, the BPS limit is

$U(1)$8

and the first-order equations are

$U(1)$9

with

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$0

taken as $$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$1 (Blanco-Pillado et al., 2013). In the compact-surface treatment of gravitating vortices, the Einstein–Bogomol'nyi equations are

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$2

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$3

for $$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$4 and $$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$5, and these are explicitly identified with Nielsen–Olesen cosmic strings in the Bogomol'nyi phase (Álvarez-Cónsul et al., 2015).

The self-dual structure can survive significant generalization of the scalar sector. In the gravitating $$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$6 sigma-model study, the BPS tension bound is

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$7

so that

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$8

with first-order equations

$$S=\int d^4x\,\sqrt{|g|}\left( \frac12 D_\mu\Phi^* D^\mu\Phi -\frac{\lambda}{4}(\Phi^*\Phi-\nu^2)^2 -\frac14 F_{\mu\nu}F^{\mu\nu} +\frac{1}{16\pi G}f(R) \right),$$9

$f(R)=R+\eta R^2,$0

(Alonso-Izquierdo et al., 2020). Likewise, in the generalized Abelian Higgs model with dielectric factor $f(R)=R+\eta R^2,$1, the self-dual equations are

$f(R)=R+\eta R^2,$2

subject to

$f(R)=R+\eta R^2,$3

(Cao et al., 2023). These generalizations preserve the gravitating-vortex mechanism while altering the detailed core structure and target-space geometry.

A common misconception is that “gravitating cosmic string” always means only the conical metric of an infinitesimally thin string. The surveyed papers show a broader usage: the phrase also covers finite-width self-gravitating vortices with regular cores, compact-surface BPS solutions, supermassive compactifying defects, and fully dynamical Einstein–matter evolutions (Álvarez-Cónsul et al., 2015, Blanco-Pillado et al., 2013, Helfer et al., 2018).

3. Beyond minimal Einstein gravity

Several papers examine how the standard gravitating Abelian–Higgs string is modified when the gravitational sector is changed.

In the Starobinsky model,

$f(R)=R+\eta R^2,$4

and the metric field equations become

$f(R)=R+\eta R^2,$5

with trace equation

$f(R)=R+\eta R^2,$6

(Graça, 2015). The paper’s central conclusion is that increasing the Starobinsky coupling lowers the angular deficit: $f(R)=R+\eta R^2,$7 It also raises the critical value $f(R)=R+\eta R^2,$8 at which $f(R)=R+\eta R^2,$9, thereby enlarging the range of symmetry-breaking scales for which regular strings exist (Graça, 2015). The matter core remains close to the Einstein case, while the geometric response is softened.

In Rastall gravity, the field equations are modified to

$$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$0

$$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$1

(Mello et al., 2014). The matter equations are not derived from a total action and are chosen by a “minimal deformation” prescription. The resulting strings still asymptote to a conical spacetime, but lowering the Rastall parameter from its Einstein value increases the deficit angle and can drive the system to a supermassive regime where $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$2 at finite radius (Mello et al., 2014). The paper also observes a “would-be-BPS bound” in which $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$3 can be recovered for tuned parameters, though it is explicitly stated not to be connected to an underlying BPS structure (Mello et al., 2014).

In the warped five-dimensional braneworld model, the Abelian–Higgs matter sector on the brane remains the standard Nielsen–Olesen vortex, but the induced four-dimensional field equations are

$$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$4

(Slagter et al., 2012). The major conclusion is that the usual conical asymptotics are lost on the brane: the induced asymptotic metric is neither conical nor standard Kasner in the ordinary four-dimensional string sense (Slagter et al., 2012). This is a precise example of a gravitating Abelian–Higgs string for which the familiar deficit-angle interpretation fails, not because the vortex disappears, but because bulk/brane corrections change the asymptotic geometry.

A plausible implication is that “gravitating Abelian–Higgs string” is best treated as a family of coupled gauge–scalar–gravity systems rather than a single universal metric type. The matter core can remain recognizably Nielsen–Olesen while the gravitational exterior changes qualitatively across theories (Graça, 2015, Mello et al., 2014, Slagter et al., 2012).

4. Supermassive strings, compactification, and moduli

When the gravitational backreaction becomes strong enough that the deficit angle reaches or exceeds $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$5, the usual open conical interpretation breaks down. In the four-dimensional Einstein–Abelian–Higgs model, the supermassive-string paper defines the criterion

$$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$6

with supermassive strings satisfying $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$7 (Blanco-Pillado et al., 2013). Then $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$8 becomes negative asymptotically, so the circumference decreases and the transverse space closes at finite proper distance $$F_{\mu\nu}=\nabla_\mu A_\nu-\nabla_\nu A_\mu,\qquad D_\mu=\nabla_\mu-i e A_\mu$$9 (Blanco-Pillado et al., 2013). A single supermassive string generically ends at a second conical singularity with

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$0

(Blanco-Pillado et al., 2013).

The same paper shows that globally regular compactifications can be obtained by replacing that would-be singularity with a second vortex. Matching two string solutions across North and South patches leads to the metric and matter conditions

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$1

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$2

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$3

and produces the global relation

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$4

(Blanco-Pillado et al., 2013). In the BPS limit $$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$5, some of these compactifications embed into $$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$6 supergravity and are half-BPS (Blanco-Pillado et al., 2013).

A related but distinct modulus phenomenon appears in gravitating strings with flat directions. In semilocal, tachyonic, and axionic BPS models coupled to gravity, the paper proves that the flat-direction zero mode survives gravitational coupling (Hartmann et al., 2012). For semilocal strings, the family of solutions is parameterized by

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$7

while the total energy remains fixed,

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$8

and the deficit angle stays unchanged,

$$S = S_{EH}-\int d^4x\sqrt{-g}\left[(D_{\mu}{\phi})^{*}(D^{\mu}{\phi})+\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+V(\phi) \right],$$9

(Hartmann et al., 2012). Yet the local energy density and the metric profile $$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$0 do change. The paper’s central lesson is that tension degeneracy does not imply local gravitational degeneracy (Hartmann et al., 2012).

This point is important for interpreting finite-width gravitating strings. A common oversimplification is that once the integrated tension is fixed, the gravitational field is fixed. The flat-direction analysis shows this is not generally true: asymptotic conical data can coincide while the local geometry, energy-density distribution, and even the size of a closed transverse space in supermassive cases differ across a modulus family (Hartmann et al., 2012).

5. Geodesics, lensing, and optical structure

Finite-width gravitating Abelian–Higgs strings are not optically equivalent to the ideal infinitely thin conical string. The 2026 optical study analyzes the straight Einstein–Abelian–Higgs vortex directly and finds several effects absent in the idealized model (Silva et al., 3 Jun 2026).

The static metric is

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$1

with asymptotics

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$2

(Silva et al., 3 Jun 2026). The Ricci scalar is

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$3

and the core radius is defined by

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$4

(Silva et al., 3 Jun 2026). Numerically, the string core is a smooth localized region of negative Ricci curvature.

Photon propagation is treated Hamiltonianly, with lens equation

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$5

and magnification

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$6

(Silva et al., 3 Jun 2026). The paper’s most distinctive optical result is finite-core triple imaging: for suitable source positions, two external rays pass outside the core and one central ray crosses the regular interior. The central image is strongly demagnified because near the axis $$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$7 is large, so

$$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$8

(Silva et al., 3 Jun 2026). In the ideal-string limit, the central image becomes effectively inaccessible.

The same paper also identifies an intrinsic Shapiro time delay due to the nontrivial lapse $$S_{EH}=\int d^4x\,\sqrt{-g}\,\frac{R}{16\pi G}, \qquad V(\phi)=\frac{\lambda}{4}\left(|\phi|^2-\eta^2\right)^2$$9,

$$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$0

(Silva et al., 3 Jun 2026). Its sign is controlled by the mass ratio parameter $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$1:

• for $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$2, $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$3, so the core acts as a temporal shortcut;
• for $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$4, $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$5, so the core acts as a temporal barrier;
• at $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$6, $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$7 and $$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$8 (Silva et al., 3 Jun 2026).

Geodesic structure can become even richer when more than one Abelian–Higgs sector is present. In the study of two magnetically interacting Abelian–Higgs strings, the static metric is

$$ds^2=-N^2(r)\,dt^2+dr^2+L^2(r)\,d\phi^2+N^2(r)\,dz^2,$$9

with deficit angle

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$0

(Hartmann et al., 2012). The effective potential for test particles is

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$1

and the paper shows that magnetic interaction between the two string sectors can create bound timelike orbits even where noninteracting strings would not (Hartmann et al., 2012). By contrast, bound null geodesics are shown to be impossible (Hartmann et al., 2012).

These results correct another common misconception: a cosmic-string spacetime is not exhausted by the ideal double-image conical lens. Finite-width regular cores, nontrivial lapse functions, and field-theoretic interactions can produce central images, demagnification, time-delay sign changes, and core-scale bound orbits (Silva et al., 3 Jun 2026, Hartmann et al., 2012).

6. Time-dependent gravity, collapse, and gravitational radiation

Gravitating Abelian–Higgs strings also appear as fully dynamical sources of gravitational radiation. Two distinct regimes emerge in the surveyed work.

The first is genuine nonlinear Einstein–matter dynamics. In the 3+1 numerical-relativity study of collapsing Abelian–Higgs loops, the theory is evolved with the full Einstein equations and the Abelian–Higgs matter system using GRChombo and the BSSN formulation (Helfer et al., 2018). The loops are circular, planar, and initially at rest. At critical coupling,

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$2

(Helfer et al., 2018). The characteristic core thickness is

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$3

and the loop mass is approximated by

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$4

(Helfer et al., 2018).

The simulations reveal a two-outcome structure: subcritical loops shrink to their core width, unwind, and disperse, whereas supercritical loops form black holes (Helfer et al., 2018). The transition is captured by a hoop-conjecture estimate based on

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$5

which yields

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$6

(Helfer et al., 2018). For black-hole-forming cases, the emitted gravitational wave is dominated by the $$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$7 mode, and for the representative case

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$8

the total radiated efficiency is

$$\Phi(r)=\nu f(r)e^{i\phi}, \qquad A_\mu dx^\mu=\frac{1}{e}[1-P(r)]\,d\phi$$9

(Helfer et al., 2018).

The second regime is perturbative gravitational radiation from dynamically formed Abelian–Higgs strings during preheating. In that work, gravity is not treated as a nonlinear metric sector; instead tensor perturbations satisfy

$f(R)$00

with source

$f(R)$01

(Dufaux et al., 2010). The paper emphasizes that these are not long-lived scaling networks but transient Nielsen–Olesen-like string configurations formed during tachyonic symmetry breaking. Their formation, widening, fragmentation, and decay generate a multipeaked gravitational-wave spectrum with peaks tied to the Higgs and gauge mass scales (Dufaux et al., 2010).

These two regimes delimit a useful distinction. In one, the string is self-gravitating in the strong-field sense and can form horizons (Helfer et al., 2018). In the other, the string acts as a localized anisotropic-stress source of gravitational waves on an essentially fixed background (Dufaux et al., 2010). Both fall within the modern literature on gravitating Abelian–Higgs strings, but they address different gravitational questions.

7. Network phenomenology and broader generalizations

On cosmological scales, Abelian–Higgs strings are often studied not as isolated self-gravitating vortices but as networks with gravitational and particle signatures. The multi-messenger network paper emphasizes that large-scale numerical simulations show classical field radiation is the dominant decay channel for generic Abelian–Higgs loops produced from random initial conditions, while special initial conditions can produce NG-like loops that radiate primarily gravitationally (Hindmarsh et al., 2022). The paper models this by introducing a fraction $f(R)$02 of NG-like loops and finds that explaining the NANOGrav signal while satisfying simulation and CMB constraints requires

$f(R)$03

(Hindmarsh et al., 2022). It also concludes that more than $f(R)$04 of the network energy must end up as dark matter or dark radiation in the viable region (Hindmarsh et al., 2022). This is not a study of individual self-gravitating vortices, but it is directly relevant to the gravitational phenomenology of Abelian–Higgs string networks.

Thermal-inflation strings provide another indirect comparison point. That paper studies a local $f(R)$05 string model with potential

$f(R)$06

and shows that, unlike the usual Abelian–Higgs scaling, the tension is

$f(R)$07

rather than $f(R)$08, while the loop cutoff is governed by cusp annihilation rather than gravitational radiation for the benchmark thermal-inflation parameters (Brandenberger et al., 2023). This paper is not a self-gravitating-vortex study, but it sharpens the distinction between minimal Abelian–Higgs intuition and more general local-string models.

A number of further extensions show how robust the gravitating-vortex framework is under deformation of the matter sector. The compactified six-dimensional gauge-theory model yields an effective Abelian–Higgs string whose deficit angle behaves approximately as

$f(R)$09

while the radius of the extra $f(R)$10 varies with distance from the string (Nakamula et al., 2019). The $f(R)$11 sigma-model study exhibits two species of self-dual gravitating cosmic strings associated with different coordinate charts (Alonso-Izquierdo et al., 2020). The generalized Abelian Higgs model with dielectric gauge factor $f(R)$12 preserves self-duality and yields the gravitating string equation

$f(R)$13

(Cao et al., 2023). The Chern–Simons–Higgs gravity paper, while not Maxwell–Higgs, demonstrates another mathematically sharp codimension-two self-dual system with gravitational existence conditions such as

$f(R)$14

on compact surfaces (Cao et al., 2024).

Taken together, these results support a broad but precise understanding of the subject. Gravitating Abelian–Higgs cosmic strings comprise not only the classic self-gravitating Nielsen–Olesen vortex with conical asymptotics, but also modified-gravity variants, braneworld embeddings, compactifying supermassive defects, semilocal and sigma-model generalizations, fully dynamical collapsing loops, transient string sources during preheating, and cosmological networks whose gravitational signatures depend on both $f(R)$15 and their dominant decay channels (Graça, 2015, Mello et al., 2014, Slagter et al., 2012, Blanco-Pillado et al., 2013, Helfer et al., 2018, Dufaux et al., 2010, Hindmarsh et al., 2022).