REVIEW 1 major objections 3 minor 6 cited by
The scaling density of axion strings
T0 review · 1 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper argues that the axion string density approaches a constant, $\zeta_0 = 1.19 \pm 0.20$, so recent reports of logarithmic growth are artifacts of initial conditions, and the implied axion dark matter mass rises by about 50%.
desk verdict New simulations and a promising diagnostic, but the claim that logarithmic growth is an artifact of initial conditions rests on a degeneracy the paper does not close. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the string length density parameter $\zeta$ (equivalently the mean string separation $\xi = \sqrt{V/\ell_s} = t/\sqrt{\zeta}$), and the load-bearing identity is the standard scaling fit $\xi = 2\beta(t - t_0)$, with $\zeta_0 = 1/4\beta^2$. The time offset $t_0$ is the piece that carries the argument: it is the t-axis intercept of the linear fit and is interpreted as the entire memory of the initial conditions. The paper shows that for nonzero $t_0$ the raw estimator $\zeta = t^2/\xi^2$ approaches $\zeta_0$ slowly, behaving like $\zeta_0(1 + 2t_0/t)$, which a logarithmic fit misreads as growth with coefficient $\alpha_*(t_f) \simeq -2\zeta_0 (t_0/t_f)(1 - t_0/t_f)^{-1}$. In this way, the apparent logarithmic and inverse-logarithmic corrections are tied to the ratio $t_0/t_f$, making them properties of the initial setup rather than of the scaling network.
What would settle it
A longer simulation with $t_f/t_0 \gtrsim 20$ using the same shifted estimator would settle the claim: if the residual $\hat\zeta - \zeta_0$ stays flat within $|\alpha| < 0.94 \times 10^{-2}$, standard scaling stands; if it drifts upward at the rate $\alpha_* \simeq 0.2$ reported by the logarithmic-growth groups, the artifact explanation fails.
Extended reading notes
Core claim
The paper's central claim is that the dimensionless string length density parameter $\zeta$, defined by writing the string length per unit volume as $\zeta/t^2$, reaches a constant $\zeta_0 = 1.19 \pm 0.20$ in the scaling regime, rather than growing logarithmically as claimed in Refs. [24–27]. The decisive evidence is that the mean string separation $\xi = t/\sqrt{\zeta}$ follows a straight line in time, $\xi = 2\beta(t - t_0)$, where the intercept $t_0$ encodes the initial string configuration. Fits of the logarithmic and inverse-logarithmic models give coefficients consistent with zero after the ratio $t_0/t_f$ is extrapolated to zero, and a direct fit to the residuals of the linear model gives $\alpha_r = (-0.04 \pm 0.30) \times 10^{-2}$. The paper concludes that no logarithmic correction larger than about 0.5 can occur between the Peccei-Quinn scale and the QCD scale, that constant-$\zeta$ extrapolation from simulations to cosmology is justified, and that the axion dark matter mass estimate should be revised upward by roughly 50%.
Load-bearing premise
The conclusion rests on the assumption that one fitted time offset in the linear string separation fully captures the effect of the initial conditions, so that any residual logarithmic behaviour is just the disguised approach to scaling; if a genuine logarithmic drift could hide inside that offset over the fitted time window, the claim that $\zeta$ is constant would collapse.
Editorial extensions
If this is right
- The string network can be extrapolated at constant $\zeta$ from the Peccei-Quinn transition to the QCD transition, so no large enhancement of string density near the QCD transition is expected.
- The $3\sigma$ bound $|\alpha| < 0.94 \times 10^{-2}$ on the residual logarithmic coefficient limits any drift in $\zeta$ to less than about 0.5 over the full cosmological range.
- The axion dark matter mass in the post-inflationary singlet scenario should be increased by about 50% relative to estimates based on the larger string densities.
- For network models, the energy loss rate per unit length of string must rise at the same rate as the effective string tension, otherwise constant $\zeta$ would be unsustainable.
- Frustrated-string simulations with $\kappa \simeq 50$ produce $\zeta \simeq 20$, far above the $O(1)$ value allowed here, so their large-separation extrapolations of axion number density should not be used.
Reading between the lines
- Editorial extension: the $t_0/t_f$ diagnostic provides a cheap retest of previously published data; if those runs are re-fitted with a shifted estimator and their reported $\alpha_* \simeq 0.2$ survives the extrapolation to $t_0/t_f \to 0$, the artifact explanation would be wrong.
- Editorial extension: the same residual analysis could be applied to local (gauge) cosmic string networks, whose tension does not grow logarithmically, to test whether initial-condition offsets are the generic source of apparent scaling violations.
- Editorial extension: replacing the paper's figure-based reading of $n_{\rm ax}/n_{\rm mis} \simeq 0.5$ at $\kappa = 1$ with a direct simulation measurement would turn the 50% axion mass revision into a sharper prediction.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents new lattice simulations of global axion strings and compares three phenomenological models for the string-density parameter ζ: standard scaling with a constant ζ0, a logarithmic correction ζ=ζ0*+α*log(ηt), and an inverse-logarithmic correction. Fitting the mean string separation ξ to ξ=2β(t−t0), the authors find ζ0=1.19±0.20 and report a residual logarithmic coefficient αr=(−0.04±0.30)×10^-2, concluding that the apparent logarithmic growth seen in earlier work is an artifact of initial conditions and that a constant-ζ extrapolation from the Peccei-Quinn scale to the QCD scale is justified. The paper also re-analyzes previous axion-number-density estimates and suggests that the axion dark-matter mass should be revised upward by about 50%.
Significance. If the central conclusion holds, the paper resolves an active controversy: recent simulations reporting logarithmic growth of ζ would imply substantially larger axion string densities and a different axion dark-matter mass, whereas a constant ζ0≈1.2 supports the traditional scaling extrapolation. The paper's strengths include a careful decomposition of statistical and systematic uncertainties, multiple box sizes and initial correlation lengths, an infinite-volume extrapolation, a physical-width versus constant-comoving-width cross-check, and explicit comparison of the three competing models. The result is also falsifiable: the quoted 3σ bound |α|<0.94×10^-2 gives a concrete target for independent simulations.
major comments (1)
- [§Measurements and results, Eqs. (5), (9)–(10), Table III] The conclusion that the residual logarithmic coefficient is tightly bounded rests on the assumption that the free intercept t0 in Eq. (5) absorbs only initial-condition effects. This assumption is not tested, and the concern lands: if the true density parameter were ζ(t)=ζ0*+α log(ηt), then the estimator in Eq. (9) gives ζhat(t)=(1−t0/t)^2[ζ0*+α log(ηt)] ≈ ζ0*+α log(ηt)−2ζ0*t0/t for small |t0/t|. Over a finite fitting window, a fitted value t0≈−α t_m/(2ζ0*), where t_m is the midpoint of the window, cancels the logarithmic trend to second order. Thus the small αr in Table III is exactly what a genuine log correction with α≈0.2 would produce under this pipeline. Equation (8) demonstrates only the direction from t0 to an apparent logarithmic coefficient; the reverse direction is not excluded. Since the central claims of standard scaling (Eq. (11)) and of a constant-ζ extrapolation to the QCD scale both depend on this non-degeneracy, I ask for an explicit test: fit a nested model that allows both t0 and a logarithmic correction, or perform an injection-recovery test in which a known α log(ηt) term is added to the scaling model and the residual procedure is required to recover it.
minor comments (3)
- [Table III and accompanying text] The statement that the direct log-model fits are 'apparently inconsistent' with α*≃0.2 is stronger than the numbers in Table III support: for the 4k box α*=0.11±0.11, so 0.2 lies within 1σ. The tight exclusion of α*≈0.2 comes mainly from the residual bound in Eq. (12), which is the degenerate test flagged in the major comment. Please present the direct-fit and residual evidence separately when characterizing the constraint.
- [Eq. (5) and footnote 1] The notation is confusing because Eq. (5) is written in terms of physical time t and ξ, while footnote 1 defines β as the slope of comoving separation ξ/a against conformal time τ=2t/a. Please state explicitly which time variable enters the plotted and fitted quantities, and define all symbols consistently.
- [Conclusions, paragraph beginning 'An implication...'] The suggested upward revision of the axion dark-matter mass by about 50% is based on an estimate read from a figure in Ref. [30], with no propagated uncertainty. Please label this as indicative or provide an error estimate, since it appears in the abstract as a quantitative conclusion.
Circularity Check
No significant circularity: the scaling claim is an empirical model comparison with independently fitted parameters, and the cited self-work is not load-bearing.
full rationale
The central claim (standard scaling with ζ0 = 1.19 ± 0.20) rests on direct linear fits to the measured mean string separation ξ(t) with two free parameters (β, t0), Eq. (5). The alternative logarithmic and inverse-logarithmic models in Eq. (6) are fitted separately with their own parameters, and their coefficients are found to be consistent with zero. The apparent-α conversion in Eq. (8) is an algebraic identity explaining how a non-zero t0 can masquerade as a logarithmic term; the paper tests this against simulations with different initial conditions, showing that α* tracks t0/tf and vanishes as t0/tf → 0. The residual test in Eq. (10) is a standard goodness-of-fit check: it subtracts the fitted ζ0 and tests for a remaining trend, but αr is not forced to zero by construction because a genuine logarithmic trend would survive in the residuals unless absorbed by t0, and the finite-window degeneracy is addressed empirically by varying initial conditions and fit ranges. Self-citations [21] and [36] are consistency/method references, not load-bearing assumptions. No equation reduces to its input by definition; the conclusion is an empirical measurement and model comparison, not a derivation from the premise. The skeptical degeneracy concern is a statistical identifiability caveat, not circularity.
Assumptions & free parameters
free parameters (6)
- t0 (string evolution intercept) =
Table I: about -8.9 to -34.5 for 4k runs depending on fit range
- β (slope of ξ(t) linear fit) =
Table II: 0.499±0.042 (1k), 0.486±0.036 (2k), 0.467±0.037 (4k)
- α* (log-growth coefficient in alternative model) =
4k: 0.11±0.11
- α' (inverse-log coefficient in alternative model) =
4k: -3.6±3.8
- α_r (residual log coefficient) =
(-0.04±0.30)×10^-2
- ζ*0 and ζ'0 (intercepts of alternative models) =
4k: 0.42±0.59 and 1.66±0.68
assumptions (6)
- domain assumption The axion field obeys the U(1) symmetric action Eq. (1) and equations of motion Eq. (2) in a radiation-dominated FLRW background with a ∝ τ.
- ad hoc to paper Gaussian random initial fields with power spectrum PΦ(k)=A[1+(k lφ)^2]^-1 and a diffusion preparation phase drive the network into the scaling regime within the simulation time.
- domain assumption The physical string length is correctly measured by counting phase-winding plaquettes with a 2/3 Manhattan correction.
- domain assumption The single time offset t0 in Eq. (5) fully accounts for initial-condition effects in the standard scaling fit.
- domain assumption The QCD axion dark matter scenario with post-inflationary Peccei-Quinn breaking, strings, QCD domain walls, and axion number density proportional to ζ applies to the real universe.
- domain assumption Linear extrapolation in 1/L to infinite volume is valid.
Cite this review
Pith. "Pith review of The scaling density of axion strings." pith.science (2026). https://pith.science/paper/7PY2LQY7
@misc{pith2026190803522,
author = {Pith},
title = {Pith review of: The scaling density of axion strings},
year = {2026},
howpublished = {\url{https://pith.science/paper/7PY2LQY7}},
note = {Machine review of arXiv:1908.03522}
}
abstract
In the QCD axion dark matter scenario with post-inflationary Peccei-Quinn symmetry breaking, the number density of axions, and hence the dark matter density, depends on the length of string per unit volume at cosmic time $t$, by convention written $\zeta/t^2$. The expectation has been that the dimensionless parameter $\zeta$ tends to a constant $\zeta_0$, a feature of a string network known as scaling. It has recently been claimed that in larger numerical simulations $\zeta$ shows a logarithmic increase with time, while theoretical modelling suggests an inverse logarithmic correction. Either case would result in a large enhancement of the string density at the QCD transition, and a substantial revision to the axion mass required for the axion to constitute all of the dark matter. With a set of new simulations of global strings we compare the standard scaling (constant-$\zeta$) model to the logarithmic growth and inverse-logarithmic correction models. In the standard scaling model, by fitting to linear growth in the mean string separation $\xi = t/\sqrt{\zeta}$, we find $\zeta_0 = 1.19 \pm 0.20$. We conclude that the apparent corrections to $\zeta$ are artefacts of the initial conditions, rather than a property of the scaling network. The residuals from the constant-$\zeta$ (linear $\xi$) fit also show no evidence for logarithmic growth, restoring confidence that numerical simulations can be simply extrapolated from the Peccei-Quinn symmetry-breaking scale to the QCD scale. Re-analysis of previous work on the axion number density suggests that recent estimates of the axion dark matter mass in the post-inflationary symmetry-breaking scenario we study should be increased by about 50%.
Figures
Forward citations
Cited by 6 Pith papers
-
CosmoLattice 2.0
CosmoLattice v2.0 extends lattice cosmology simulations with non-minimal scalars, ALP–gauge couplings, defect networks, low-storage RK integrators, optimized GWs, and O(10) GPU speedups.
-
Nambu-Goldstone emissions from the cosmological evolution of global monopoles
First quantitative lattice measurement of NG boson emission from global monopoles: the spectrum peaks at the Hubble scale, the number density grows linearly with H, and the resulting pseudo-NG bosons can be dark matter.
-
Formation and scaling of $\mathbb{Z}_N$ strings for global $\mathrm{SU}(N)/\mathbb{Z}_N$ symmetry
In a global SU(N)/Z_N scalar model, Z_N-string networks with baryon-vertex-like junctions reach a scaling regime for N=2,3,4,5,8, with string density proportional to N^2-1.
-
MOSAIC: Magnonic Observations of Spin-dependent Axion-like InteraCtions
MOSAIC is a proposed scalable array of YIG magnon spheres coupled to electron-on-neon qubits that could search for electron-coupled axion dark matter with projected sensitivity beyond current ferromagnetic haloscopes.
-
QCD Axion on Hilltop by a Phase Shift of $\pi$
A heavy axion inflaton can shift the QCD axion potential by π, placing the QCD axion at the hilltop and allowing f_a ≳ 3×10^9 GeV to explain all dark matter.
-
The Spectrum of Global Axion Strings
Lattice simulations of global axion string networks, corrected for discretization and oscillation effects, are extrapolated to predict the axion dark matter mass in the 95 to 450 micro-eV band.
Reference graph
Works this paper leans on
-
[1]
R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977); Phys. Rev. D16, 1791 (1977)
work page 1977
-
[2]
S. Weinberg, Phys. Rev. Lett.40, 223 (1978); F. Wilczek, Phys. Rev. Lett. 40, 279 (1978)
work page 1978
-
[3]
J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shif- man, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B166, 493 (1980)
work page 1979
-
[4]
A. R. Zhitnitsky, Sov. J. Nucl. Phys.31, 260 (1980), [Yad. Fiz.31,497(1980)]; M. Dine, W. Fischler, and M. Sred- nicki, Phys. Lett. 104B, 199 (1981)
work page 1980
-
[5]
J. Preskill, M. B. Wise, and F. Wilczek, Phys. Lett. B120, 127 (1983); L. F. Abbott and P. Sikivie, Phys. Lett. B120, 133 (1983); M. Dine and W. Fischler, Phys. Lett. B120, 137 (1983)
work page 1983
-
[6]
R. L. Davis, Phys. Lett. B180, 225 (1986)
1986
-
[7]
M. Hindmarsh and T. Kibble, Rept.Prog.Phys. 58, 477 (1995), arXiv:hep-ph/9411342 [hep-ph]
arXiv 1995
-
[8]
Vilenkin and E
A. Vilenkin and E. P. S. Shellard, Cosmic Strings and Other Topological Defects (Cambridge University Press, 2000)
2000
Show all 45 references
-
[9]
Sikivie, Phys
P. Sikivie, Phys. Rev. Lett. 48, 1156 (1982)
1982
-
[10]
Georgi and M
H. Georgi and M. B. Wise, Phys. Lett. 116B, 123 (1982)
1982
-
[11]
E. P. S. Shellard and R. A. Battye, Sources and de- tection of dark matter in the universe. Proceedings, 3rd International Symposium, and Workshop on Primordial Black Holes and Hawking Radiation, Marina del Rey, USA, February 17-20, 1998, Phys. Rept.307, 227 (1998), arXiv:astr...
1998 arXiv
-
[12]
Sikivie, Axions: Theory, cosmology, and experi- mental searches
P. Sikivie, Axions: Theory, cosmology, and experi- mental searches. Proceedings, 1st Joint ILIAS-CERN- CAST axion training, Geneva, Switzerland, November 30-December 2, 2005, Lect. Notes Phys. 741, 19 (2008), [,19(2006)], arXiv:astro-ph/0610440 [astro-ph]
2008 arXiv
-
[13]
J. E. Kim and G. Carosi, Rev. Mod. Phys.82, 557 (2010), arXiv:0807.3125 [hep-ph]
2010 arXiv
-
[14]
Yamaguchi, M
M. Yamaguchi, M. Kawasaki, and J. Yokoyama, Phys. Rev. Lett. 82, 4578 (1999), arXiv:hep-ph/9811311 [hep- ph]
1999 arXiv
-
[15]
Yamaguchi, M
M. Yamaguchi, M. Kawasaki, and J. Yokoyama, in Proceedings, 7th International Symposium on Particles, Strings and Cosmology (PASCOS 99): Lake Tahoe, Cal- ifornia, December 10-16, 1999 (1999) pp. 291–294
1999
-
[16]
Yamaguchi, M
M. Yamaguchi, M. Kawasaki, and J. Yokoyama, in Pro- ceedings, 3rd International Workshop on The identifica- tion of dark matter (IDM 2000): York, UK, September 18-22, 2000 (2000) pp. 297–304
2000
-
[17]
Yamaguchi and J
M. Yamaguchi and J. Yokoyama, Phys. Rev. D67, 103514 (2003), arXiv:hep-ph/0210343 [hep-ph]
2003 arXiv
-
[18]
Hiramatsu, M
T. Hiramatsu, M. Kawasaki, T. Sekiguchi, M. Yam- aguchi, and J. Yokoyama, Phys. Rev. D83, 123531 (2011), arXiv:1012.5502 [hep-ph]
2011 arXiv
-
[19]
Hiramatsu, M
T. Hiramatsu, M. Kawasaki, K. Saikawa, and T. Sekiguchi, Phys. Rev. D85, 105020 (2012), [Erratum: Phys. Rev.D86,089902(2012)], arXiv:1202.5851 [hep-ph]
2012 arXiv
-
[20]
Kawasaki, K
M. Kawasaki, K. Saikawa, and T. Sekiguchi, Phys. Rev. D91, 065014 (2015), arXiv:1412.0789 [hep-ph]
2015 arXiv
-
[21]
Lopez-Eiguren, J
A. Lopez-Eiguren, J. Lizarraga, M. Hindmarsh, and 6 J. Urrestilla, JCAP 1707, 026 (2017), arXiv:1705.04154 [astro-ph.CO]
2017 arXiv
-
[22]
C. J. A. P. Martins and E. P. S. Shellard, Phys. Rev. D54, 2535 (1996), arXiv:hep-ph/9602271 [hep-ph]
1996 arXiv
-
[23]
Martins and E
C. Martins and E. Shellard, Phys.Rev. D65, 043514 (2002), arXiv:hep-ph/0003298 [hep-ph]
2002 arXiv
-
[24]
Gorghetto, E
M. Gorghetto, E. Hardy, and G. Villadoro, JHEP 07, 151 (2018), arXiv:1806.04677 [hep-ph]
2018 arXiv
-
[25]
Kawasaki, T
M. Kawasaki, T. Sekiguchi, M. Yamaguchi, and J. Yokoyama, PTEP 2018, 091E01 (2018), arXiv:1806.05566 [hep-ph]
2018 arXiv
-
[26]
Vaquero, J
A. Vaquero, J. Redondo, and J. Stadler, JCAP 1904, 012 (2019), arXiv:1809.09241 [astro-ph.CO]
2019 arXiv
-
[27]
Buschmann, J
M. Buschmann, J. W. Foster, and B. R. Safdi, (2019), arXiv:1906.00967 [astro-ph.CO]
2019 arXiv
- [28]
-
[29]
V. B. Klaer and G. D. Moore, JCAP 1710, 043 (2017), arXiv:1707.05566 [hep-ph]
2017 arXiv
-
[30]
V. B. Klaer and G. D. Moore, JCAP 1711, 049 (2017), arXiv:1708.07521 [hep-ph]
2017 arXiv
-
[31]
C. T. Hill, A. L. Kagan, and L. M. Widrow, Phys. Rev. D38, 1100 (1988)
1988
-
[32]
C. J. A. P. Martins, Phys. Lett. B788, 147 (2019), arXiv:1811.12678 [astro-ph.CO]
2019 arXiv
-
[33]
Daverio, M
D. Daverio, M. Hindmarsh, and N. Bevis, ArXiv e-prints (2015), arXiv:1508.05610 [physics.comp-ph]
2015 arXiv
-
[34]
W. H. Press, B. S. Ryden, and D. N. Spergel, Astrophys. J. 347, 590 (1989)
1989
-
[35]
Vachaspati and A
T. Vachaspati and A. Vilenkin, Phys. Rev. D30, 2036 (1984)
1984
-
[36]
Bevis, M
N. Bevis, M. Hindmarsh, M. Kunz, and J. Urrestilla, Phys.Rev. D82, 065004 (2010), arXiv:1005.2663 [astro- ph.CO]
2010 arXiv
-
[37]
K. J. Bae, J.-H. Huh, and J. E. Kim, JCAP 0809, 005 (2008), arXiv:0806.0497 [hep-ph]
2008 arXiv
-
[38]
Wantz and E
O. Wantz and E. P. S. Shellard, Phys. Rev. D82, 123508 (2010), arXiv:0910.1066 [astro-ph.CO]
2010 arXiv
-
[39]
Borsanyi et al
S. Borsanyi et al. , Nature 539, 69 (2016), arXiv:1606.07494 [hep-lat]. 7 SUPPLEMENTAL MATERIAL Infinite volume extrapolation of fit parameters In Fig. S1 we show the central values of the fit parame- tersζ0 andαr (see Eq. 10) as well as the 1-σ uncertainties for 1k, 2k and 4k si...
2016 arXiv
-
[40]
Simulations which are scaling in the standard sense (ξ ∝ t) have a slow evolution in the ζ - log(ηt) plane
-
[41]
Simulations converge toζ≃ 1 from both above and below
-
[42]
Convergence from below can look like a logarithmic increase in ζ, as observed in Refs. [24–27]
-
[43]
The coefficients of the logarithm and inverse loga- rithm have a strong dependence on the ratio t0/tf, consistent with their being a feature of the initial conditions
-
[44]
The coefficients of the logarithm and inverse loga- rithm are consistent with zero at t0/tf = 0
-
[45]
The constant terms in logarithm and inverse loga- rithm models are consistent with standard scaling at t0/tf = 0. 8 0 50 100 150 t 0 50 100 150ξ 0 50 100 150 200 250 300 t 0 100 200 300ξ 0 100 200 300 400 500 t 0 100 200 300 400 500ξ 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 log(ηt) 0.0...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.