Violation of Constant Scaling in Networks of Global Monopoles ✛ ✎ ✦

Original: "Constant Scaling Fails for Global Monopole Networks"
· Wakutaka Nakano, Wen Yin
arXiv:2607.05517 · 2026-07-06 · arXiv non-exclusive · independent analysis · 5 min
Large-scale numerical simulations have shown that the density of global monopoles in the early Universe does not plateau but slowly grows logarithmically, changing predictions for dark matter and gravitational waves.

Context

The study of topological defects such as global monopoles is an important bridge between particle physics and cosmology. These objects arise from the spontaneous breaking of the global symmetry O(3) → O(2) in the hot early Universe; the theoretical foundations were laid by Gerard 't Hooft and A. M. Polyakov. The evolution of defects is closely tied to the expansion rate, described by Hubble's law and the Friedmann equations, and the initial field fluctuations could have been set during the epoch of inflation, whose theoretical framework was developed by Alan Guth. If the monopole density per Hubble volume had remained constant, their contribution to cosmological observables would be rigidly fixed. However, recent hints of logarithmic corrections in cosmic strings have forced a reassessment of this postulate. Understanding the true scaling is critically important for dark matter models, in which monopoles or related particles — axions and dark photons — make up the hidden sector.

Methods

To test the hypothesis of constant scaling, the authors conducted a series of numerical experiments based on a modified version of the 3DCosmoLattice code. They simulated a scalar field with O(3) symmetry in an expanding Universe with radiation- and matter-dominated backgrounds. To maintain the resolution of the monopole core on a fixed lattice over long times, a 'fat' monopole technique was used: the coupling constant λ was effectively changed as λ/a², where a is the scale factor. This allowed a dynamic range in time of up to 10⁶. The defect abundance was estimated by counting topological charges on cubes after a smoothing flow procedure, followed by clustering and assignment to Hubble cells. The energy contributions from gradient and total density were also analyzed.

Results

The main result: the parameter ξ, which characterizes the number of monopoles plus antimonopoles per Hubble volume, does not settle to a constant value. Instead, in all investigated initial fluctuation spectra (from red to blue) a positive logarithmic drift is observed: ξ = A + B log η, where η is conformal time. For example, in the model with spectral index p=0 in the matter era, B = 0.087 ± 0.029, and in the radiation era, B = 0.584 ± 0.148. To quantify the sensitivity to the physical hierarchy mᵣ/H (mᵣ is the radial mode mass), a fractional response γ = d log ξ / d log(mᵣ/H) was introduced. It turned out that for blue spectra (p ≥ -1) γ weakly depends on the background equation of state and lies in the range 0.4–0.5. Thus, for p=0, γ = 0.393 ± 0.194 in matter and γ = 0.462 ± 0.161 in radiation. This means that ξ grows as a power law of the hierarchy, with an exponent around 0.5. Such behavior is not described by the standard instantaneous annihilation model and points to a logarithmic correction in the lifetime of the monopole-antimonopole pair.

Implications

The discovered violation of constant scaling requires a revision of analytical and semi-analytical models of defect evolution. The positive logarithmic correction means that the abundance of monopoles at later stages turns out to be higher than previously thought. This directly enhances the predicted signals: the production of dark matter (both in the form of the monopoles themselves and as axions or dark photons), the formation rate of primordial black holes, and the amplitude of the stochastic gravitational wave background. In particular, the gravitational wave spectrum generated in the radiation era may acquire a blue tilt, which is important for interpreting pulsar timing data. In addition, constraints from the anisotropy of the cosmic microwave background should be recalculated taking into account the slowly growing number of defects.

Future development

Future research should focus on the microscopic justification of the logarithmic correction. Why does the pair annihilation time depend on the logarithm of the distance to the core size? Possible reasons include angular momentum loss or the radiation of Goldstone waves. It is also necessary to study the transition to realistic gauge theories, where the long-range field is screened, and to understand how scaling changes as the monopole era approaches its end. In parallel, simulations should be extended to larger lattices and include the effects of inflationary initial conditions for more accurate predictions of cosmological signatures.

Impact

The results will impact the physics of the early Universe, the theory of dark matter, the astrophysics of black holes, gravitational wave astronomy, and the interpretation of cosmic microwave background observations.

Next steps

Conduct simulations with a fixed coupling constant (without the 'fat' monopole approximation) to verify the obtained γ values, and also calculate specific predictions for dark matter particle spectra and gravitational waves in realistic scenarios.

Key open problems

The work is directly connected to fundamental questions: what is the nature of dark matter? Can topological defects explain the recently detected nanohertz stochastic background of gravitational waves? Moreover, logarithmic corrections to scaling are a common feature of systems with long-range interactions, which echoes the problem of infrared divergences in quantum field theory and could shed light on universal aspects of non-equilibrium dynamics.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors performed numerical simulations of global monopoles in an expanding universe using a modified 3DCosmoLattice code with a 'fat' monopole technique to preserve resolution. They found that the parameter ξ does not settle to a constant but exhibits a logarithmic drift over time, with a positive coefficient B that depends on the spectrum of initial fluctuations and the background equation of state. They also introduced a fractional response γ, showing a power-law dependence of ξ on the hierarchy mᵣ/H with an exponent of about 0.5.

Where its strength is

The work systematically investigates the deviation from standard scaling across a wide range of initial conditions and backgrounds, introducing a quantitative measure γ that allows comparison of different scenarios. This provides a basis for constructing refined models of monopole evolution.

Where it could go next

Worth trying

  • Compare the drift of ξ in your simulations with predictions from models with logarithmic corrections for cosmic strings to understand the universality of the effect.
  • Estimate how the obtained dependence of ξ on mᵣ/H changes the expected gravitational wave background from monopoles, using methods from the work on strings.
  • Check whether the logarithmic growth of ξ affects constraints on dark matter models with monopoles, for example, on annihilation rates.

Why this may matter

Refining the scaling of global monopoles is important for cosmological constraints on grand unified theories and for predictions of gravitational wave signals from topological defects.

If the monopole network strictly kept a constant number per Hubble volume, only a few such objects would remain in the visible Universe today. However, the logarithmic growth means there could be hundreds or thousands of times more — enough to leave a trace in gravitational lensing or polarization of the cosmic microwave background, making them potentially detectable.

E(R) \simeq 4\pi v^2 R
the monopole energy grows linearly with the distance R to the antimonopole due to gradients of the Goldstone boson field
\tau_{\rm ann} = D \left( c_0 + c_1 \ln\frac{D}{r_{\rm core}} \right)
The pair annihilation time depends on the ratio of the distance D to the core size r_core, which gives rise to a logarithmic correction
\gamma = \frac{d\log \xi}{d\log(m_r/H)}
Characterizes the relative growth of monopole abundance with changes in the hierarchy between the radial mode mass and the Hubble parameter

Key numbers

  • dynamic range in time: up to 10⁶
  • γ for p=0 in the matter era: 0.393 ± 0.194
  • ξ_count at η=91 for p=0 (matter): 0.31 ± 0.06
  • monopole mass at v=10¹¹ GeV: 3.4×10¹⁰ g
Concepts
galaxy growth mechanisms eternal inflation structure growth magnetic monopole Dollo's law scalable Cope's rule large-scale structure of the universe magnetic monopole dark matter Gause's competitive exclusion principle quantum cosmological inhomogeneities galaxy distribution cosmic expansion cosmological constant problem topological pattern formation inflationary cosmology multi field inflation scale factor galactic evolution processes hubble constant