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.
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.
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.
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 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.
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.
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.
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.
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
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.
Nearby works in our archive
The closest work in the archive sits at 0.67 by meaning — the usual level for our corpus.
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.