Modern experiments, such as pulsar timing arrays (PTAs) and ground-based interferometers (like the LIGO-Virgo-KAGRA network built with contributions from Rainer Weiss), aim to detect the stochastic gravitational-wave background. One promising cosmological source is domain walls — two-dimensional defects that arise from spontaneous breaking of a discrete symmetry in the early Universe shortly after the Big Bang. If the symmetry is exact, the wall network quickly enters a scaling regime and begins to dominate, which is incompatible with observed cosmology described by Hubble's law and the Einstein field equations. Fortunately, even a tiny symmetry breaking — either explicit in the potential or built into the initial conditions — causes the walls to annihilate. So far, attention has focused on the first mechanism, but this work is the first to explore the second in detail: a population bias, where one vacuum is initially slightly more populous.
The field dynamics were tracked using classical 3+1 lattice simulations in an expanding Universe governed by the Friedmann equations, employing a modified version of the CosmoLattice code. A key innovation was the 'fattening' technique, which artificially keeps the physical width of the walls constant, overcoming resolution loss at late times. The initial bias was set by a homogeneous shift in the mean field, with its magnitude encoded in the bias parameter $B_s$. The simulations spanned epochs from wall formation at temperatures around their mass scale to late-time annihilation, with a specialized friction phase at the start accelerating the approach to the scaling regime and suppressing spurious scalar waves.
The evolution of the 'false' vacuum fraction (the initially less populous one) is beautifully described by the function F ~ (1/2) exp[–((η–$η_s$)/Δη_ann)^{p + α (η–$η_s$)/Δη_ann}]. The simulations yield the parameters: exponent p ≈ 1.37, correction α ≈ 0.30, and the characteristic annihilation time Δη_ann ∝ $B_s^{-0.8}$. Annihilation is delayed compared to the potential-bias case: for the same collapse time, the population bias decays more gradually. The gravitational-wave spectrum, averaged over realizations, is well approximated by a single broken power law: a steep infrared rise (∝ $f^3$), a gentle ultraviolet fall-off (∝ $f^{-1}$), and a peak at a frequency roughly 2 times the Hubble rate at the time of generation. The integrated radiation efficiency is ϵ_{gw} ≈ 0.06, only half that of the potential bias. Importantly, the spectral shape is nearly independent of the initial bias magnitude, as confirmed by analyzing the scalar field power spectrum at equal false-vacuum fractions.
The derived gravitational-wave spectra from population-biased networks fill an important gap: now any theoretical scenario with initial vacuum asymmetry can be directly compared with data from radio astronomy PTAs and LVK. Moreover, the work demonstrates that population bias is a viable mechanism on its own, requiring no explicit symmetry breaking in the potential. This is especially valuable for models where the symmetry is gauged (prohibiting a potential bias) or, for instance, post-inflationary QCD axions with $N_{dw}$>1, where the potential bias is constrained by the strong CP problem — a problem that Roberto Peccei and Helen Quinn proposed their mechanism to solve. Thus, our results extend the phenomenological viability of spontaneously broken discrete symmetries.
The immediate observational prospects involve searching PTA data (NANOGrav, EPTA, PPTA) and the fourth observing run of LIGO-Virgo-KAGRA for the described signal. On the theoretical side, it is important to push into the regime of even smaller initial biases to more reliably confirm the universality of the spectrum and clarify the behavior in the infrared tail. Of separate interest is the inclusion of gauge models and the connection to primordial black hole formation from collapsing closed walls — preliminary estimates suggest potential observability.
Observational cosmology and astrophysics: refinement of the parameters of the stochastic gravitational-wave background. Fundamental particle physics: constraints on models with discrete symmetry breaking.
Running longer simulations with improved resolution of the infrared part of the spectrum. Developing optimal strategies for extracting the signal from population-biased walls in real data from gravitational-wave observatories.
The problem of cosmic domain walls (why we don't see them today) and the strong CP problem (axions) — both may be connected through mechanisms for wall annihilation without exact symmetry.
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 paper investigates the evolution of a network of domain walls with an initial population bias in the early Universe. Using lattice simulations with a wall 'inflation' technique, the authors traced the annihilation dynamics and calculated the gravitational wave spectrum. They found that the annihilation time follows a simple power law, and the spectral shape is almost independent of the magnitude of the bias.
Where its strength is
For the first time, a scenario with an initial population bias, rather than explicit symmetry breaking in the potential, has been studied in detail. A wall 'inflation' technique has been proposed and applied, allowing correct modeling of the late-time evolution without loss of resolution.
Where it could go next
Worth trying
Why this may matter
Understanding how initial conditions affect the gravitational wave signal from topological defects is important for interpreting future observations of the background. If a population bias leaves a unique signature, it would open a new way to probe the earliest moments of the Universe, independent of the details of the potential.
Nearby works in our archive
The closest work in the archive sits at 0.69 by meaning — the usual level for our corpus.
If domain walls stretched across the visible Universe, crossing one would create a density jump capable of tearing galaxies apart. Fortunately, a tiny initial imbalance — just 0.01% — makes them vanish long before the first stars form.
The notion of separated vacua echoes the worlds of Roger Zelazny's 'Chronicles of Amber,' where Shadows exist between realities and crossing them takes effort. Here, however, the transitions between vacua turn out to be unstable and collapse on their own.