Pith. sign in

REVIEW 2 major objections 5 minor 38 references

PQ-ball and its Real Scalar Analogue in an Expanding Universe

T0 review · 2 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper shows that PQ-balls — quasi-stable scalar lumps in a spontaneously broken U(1) theory — form from a coherently rotating field in an expanding universe and decay on timescales set by the symmetry-breaking scale, with longer…

desk verdict Useful 3+1D demonstration of PQ-ball formation/decay, but the v^4 lifetime scaling claim isn't supported by the published numbers. read the letter →

arxiv 2507.01082 v1 pith:5M3VTMKD submitted 2025-07-01 hep-ph astro-ph.CO

Kazunori Nakayama , Masaki Yamada This is my paper
classification hep-phastro-ph.CO
keywords PQ-ballQ-balloscillonspontaneoussymmetrybreakinglatticesimulationaxionkineticmisalignmentAffleck-Dinemechanismexpandinguniverse
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper claims that in a complex scalar field theory with a spontaneously broken $U(1)$ symmetry, a coherently rotating field with small perturbations develops localized, quasi-stable lumps called PQ-balls in an expanding universe, and that these lumps later decay because the symmetry-breaking potential lets charge leak away. The claim matters for axion kinetic misalignment models, where the lifetime of these lumps determines when the Peccei-Quinn field thermalizes and whether the scenario remains viable. The evidence comes from $3+1$-dimensional classical lattice simulations starting from a rotating field of amplitude $M$ with a logarithmically running mass, together with a real scalar analogue in which the same formation-and-decay story holds for oscillons.

What carries the argument

The load-bearing object is the logarithmic running-mass potential $V(\Phi)=m^2[1+k_1\log(|\Phi|^2/M^2)-k_2\log((|\Phi|^2+v^2)/M^2)]|\Phi|^2$ with $k_2>k_1>0$, which is shallower than quadratic at large field values and therefore supports Q-ball-like solutions while the $v$ term softly breaks the $U(1)$ symmetry. The analytic core is the surface charge flux estimate $dQ/dt \sim -4\pi v_q \omega R_Q^2 \Phi_\infty^2$, giving $(1/Q)\,dQ/dt \sim m\sqrt{k_2-k_1}(M/\Phi_0)^2(v/M)^{2k_2/k_1}\exp(\cdots)$; this is the formula the simulations test. The numerical machinery is a sixth-order symplectic integrator on a $256^3$ lattice in a comoving box of side $L=20$, with matter-dominated expansion $a\propto t^{2/3}$ and quantum fluctuations seeded with a momentum cutoff.

What would settle it

Repeat the lattice runs at higher resolution and larger box (for example $N=512$, $L=40$) or at additional SSB scales such as $v=0.15,0.25,0.35$ and compare disappearance times: if lifetimes shift by more than order one under resolution changes, or if the measured scaling deviates from $(v/M)^{2k_2/k_1}=v^4$ beyond the scatter in the three current points, the formation-and-decay picture would be in doubt. A direct check would measure the outgoing charge flux at the surface of one isolated PQ-ball and compare it with Eq. (2.11).

Watch

Extended reading notes

Core claim

The central finding is that PQ-ball formation and decay occur generically in an expanding, spontaneously broken $U(1)$ theory without assuming spherical symmetry: starting from a rotating complex scalar of amplitude $M$ with quantum seed perturbations, charge condenses into localized lumps by $mt \sim 100$, those lumps remain stable when the symmetry-breaking scale $v$ is set to zero, and for $v>0$ they eventually dissolve. The full decay time decreases monotonically with $v$, with lifetimes of order $10^4$ for $v=0.2$, $4000$ for $v=0.3$, and $1500$ for $v=0.4$, in marginal agreement with the analytic scaling $(1/Q)\,dQ/dt \sim (v/M)^{2k_2/k_1}$ inherited from the spherical model. The same qualitative behavior holds for a real scalar field, where the analogous objects are oscillons: they form, live long for small $v$, and decay faster as $v$ grows, even though the real field traverses the symmetry-breaking barrier every oscillation.

Load-bearing premise

The central claim rests on a single benchmark simulation ($k_1=0.1$, $k_2=0.2$, $m/M=0.1$, $L=20$, $N=256$, time step $0.01$) being representative: no convergence checks on lattice spacing, box size, or time step are reported, and the claimed $v^4$ lifetime scaling is checked at only three values of $v$ and found only marginally consistent.

Editorial extensions

If this is right

  • In axion kinetic misalignment models, the PQ field should not be treated as a homogeneous rotating condensate: it fragments into PQ-balls whose eventual decay sets the thermalization epoch of the PQ sector, and lowering the SSB scale postpones that epoch.
  • The decay rate scales roughly as $(v/M)^{2k_2/k_1}$, so for small $v$ a PQ-ball can outlive the moment its surrounding background falls below threshold; its lifetime is then set by intrinsic surface flux, not by Hubble expansion.
  • Because small PQ-balls decay earlier and release charge that can be reabsorbed by larger survivors, the total charge inside identified PQ-balls shows temporary rises during the decay phase rather than a monotonic decline.
  • In the real scalar theory, oscillons form under the same large-amplitude, logarithmic-potential conditions and their lifetime is shortened by the SSB term, so the qualitative phenomenon does not depend on the field being complex.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Inference: the reported simulations establish a proof of principle, not a universal scaling law; a convergence study would be the natural next step, and if resolution artifacts affect the measured lifetimes then Eq. (2.13)'s $v^{2k_2/k_1}$ dependence could be an artifact of the single benchmark.
  • Inference: the same lattice setup could be extended to the QCD axion kinetic misalignment model with the periodic QCD potential included, in which case the decay products of PQ-balls would be axions, possibly observable as dark radiation or a stochastic axion background.
  • Inference: the oscillon result suggests a new decay channel for oscillons in any real scalar theory with radiatively induced symmetry breaking, which could shorten oscillon lifetimes in preheating scenarios where they are otherwise assumed extremely long-lived.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper studies a complex scalar field with the spontaneously-broken-U(1) potential V(Φ) = m^2[1 + k1 log(|Φ|^2/M^2) − k2 log((|Φ|^2+v^2)/M^2)]|Φ|^2 (k2 > k1 > 0) and uses 3+1-dimensional classical lattice simulations in a matter-dominated expanding universe. Starting from a large-amplitude coherently rotating field with small quantum fluctuations, the authors observe the formation of localized, U(1)-charged configurations ('PQ-balls') that remain stable in the v → 0 limit and decay for non-zero SSB scale v after the background amplitude falls below a threshold Φ∞ (Eq. 2.7). They report disappearance times t_d ≈ 10^4, 4×10^3, 1.5×10^3 for v = 0.2, 0.3, 0.4 and state that these are 'marginally consistent' with the analytic lifetime scaling τ ∝ v^{2k2/k1} = v^4 from Eq. (2.13). The paper also presents an analogous real-scalar simulation in which oscillons form and acquire a v-dependent decay channel.

Significance. If the quantitative claims hold, the paper provides an independent 3+1D confirmation that Q-ball-like configurations can form in a spontaneously broken U(1) theory with a running mass, which is directly relevant to the axion kinetic misalignment mechanism and PQ cosmology. The use of identical initial fluctuations across v values and the v = 0 control simulation are good methodological features, and the real-scalar analogue extends the phenomenology to oscillons. The numerical setup is not merely a fit to the analytic formulas of Ref. [12], so the qualitative formation/decay picture is a useful check. However, the central quantitative claim—the v^4 lifetime scaling—is not established by the presented diagnostics, because the reported disappearance times mix the background threshold time with the intrinsic decay time, condition (2.14) is never checked, and no convergence or uncertainty estimates are given. The qualitative statement that larger v leads to shorter lifetimes is credible, but the specific power-law index is not.

major comments (2)
  1. [§3; Eqs. (2.7), (2.13), (2.14)] The reported 'lifetimes' t_d = 10^4, 4×10^3, 1.5×10^3 are total disappearance times, not intrinsic PQ-ball decay times. Decay begins only after the rolling background amplitude satisfies ⟨Φ(t)⟩ ≃ Φ∞(v), i.e. after t_th ∼ t0/Φ∞ ∝ v^{-k2/k1} = v^{-2} for the benchmark k2/k1 = 2 in the matter-dominated background. For v = 0.2, 0.3, 0.4, the expected threshold-time ratios are 2.25 and 1.78, while the measured t_d ratios are 2.5 and 2.7, bracketed between the v^{-2} expectations (2.25, 1.78) and the v^{-4} expectations (5.06, 3.16). Condition (2.14), which would justify neglecting the threshold contribution, is never evaluated for the benchmark parameters. Thus the data do not discriminate the claimed v^{2k2/k1} intrinsic scaling from a threshold-dominated scaling. The authors should either measure the decay time after t_th, explicitly verify and state (2.14), or restrict the claim to qualitative v-dependence.
  2. [§3, benchmark parameters and Fig. 3] The quantitative lifetime claim rests on a single benchmark simulation (k1 = 0.1, k2 = 0.2, m/M = 0.1, L = 20, N = 256, Δt = 0.01) without any convergence test in lattice spacing, box size, or time step, and without uncertainties on the disappearance times. The non-monotonic features in Fig. 3, which the authors attribute to decay products being reabsorbed by larger PQ-balls in a finite box, show that finite-volume effects can affect the charge-in-ball curves. Since the claimed v^4 scaling is inferred from only three v values and the paper itself describes the agreement as 'marginally consistent,' resolution or finite-volume effects could change the inferred power-law index. An estimate of systematic errors, or at least one convergence check, is needed to support the quantitative claim.
minor comments (5)
  1. [§1] The sentence 'we discuss the cosmological implications of our finidings' contains a typo: 'finidings' should be 'findings'.
  2. [§3, first sentence] The phrase 'FLR W Universe' appears to be a typo for 'FLRW Universe'.
  3. [Eq. (2.13)] The displayed prefactor m/√(k2−k1) does not follow from substituting R_Q = √2/(m√(k2−k1)) into Eq. (2.12), which would give m√(k2−k1)/√2. The difference is presumably absorbed in the stated 'O(1) numerical factors,' but for the benchmark k2−k1 = 0.1 the discrepancy is a factor of 10 and deserves a clarifying comment.
  4. [Figs. 3 and 6] The vertical axes of Figs. 3 and 6 are not labeled or described in the captions; please state the plotted quantities (Q_in and ρ_in) and their units in the code's normalization.
  5. [§3, initial fluctuations] The statement that 'the precise form of the initial fluctuations does not qualitatively affect our results' is asserted without a supporting test; given that the formation timescale depends logarithmically on the fluctuation amplitude, a brief sensitivity study would strengthen this claim.

Circularity Check

0 steps flagged · score 2.0 of 10

No material circularity: the v^4 lifetime scaling from the authors' prior work is tested against an independent 3+1D lattice simulation, not fitted or defined into existence.

full rationale

The paper's central quantitative claim is that PQ-ball disappearance times scale as v^{2k2/k1}=v^4, following Eq. (2.13), which was derived in the authors' previous paper Ref. [12]. This is a self-citation, but it is not circular: the present paper performs fresh, self-contained 3+1D classical lattice simulations (Section 3) with a fixed benchmark setup (k1=0.1, k2=0.2, m/M=0.1, L=20, N=256) and varying v. The disappearance times are read off from the simulations (t_d ~ 10^4, 4e3, 1.5e3 for v=0.2, 0.3, 0.4), not obtained by substituting simulation outputs into Eq. (2.13). No parameter of Eq. (2.13) is fitted to the simulation data, so the comparison is an external test rather than a tautology. The identification criterion for PQ-balls (charge density exceeding ten times the average) is an operational definition and does not by itself enforce the v^4 scaling. The unresolved issues noted by the skeptic—single benchmark, no convergence checks, the fact that t_d includes the v-dependent threshold time before decay begins, and condition (2.14) never being verified—are correctness or robustness concerns, not circularity. They reduce confidence in the quantitative index but do not make the derivation equivalent to its inputs. Accordingly, the only reason for a nonzero score is the moderate reliance on the same authors' earlier analytic formulas; the central numerical demonstration remains independent.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper's central claim rests on a model and analytic formulas from the authors' prior work (Refs [12,20-22]), plus standard Q-ball formation results from the literature. The parameters k1, k2, m/M are hand-chosen benchmarks. No new particles or entities are introduced.

free parameters (4)
  • k1 = 0.1
    Log slope of the running mass below the scale v; chosen by hand for the benchmark. The lifetime exponent 2k2/k1 depends on it.
  • k2 = 0.2
    Log slope above v; with k1=0.1 gives K=k1-k2=-0.1, the gravity-mediation-like condition for Q-ball formation.
  • m/M = 0.1
    Ratio used to set the amplitude of seeded quantum fluctuations after rescaling; chosen as an example.
  • t0 = 100
    Initial time in units of 1/m at which the rotating field and fluctuations are imposed; hand-chosen.
assumptions (5)
  • domain assumption The potential (2.1) with k2 > k1 > 0 realizes a spontaneously broken vacuum at fa much smaller than v and M via radiative stabilization.
    Taken from Refs [20-22]; the hierarchy is central to the PQ-ball setup.
  • domain assumption For K = k1 - k2 < 0, small perturbations of a rotating scalar field grow and form Q-balls in the v=0 limit.
    Established in Refs [4-11]; motivates the expectation that PQ-balls form.
  • domain assumption The analytic PQ-ball profile, radius, critical amplitude, and decay rate (Eqs. 2.4-2.13) from Ref [12] are valid.
    The numerical results are interpreted through these formulas; Ref [12] is the authors' own prior work.
  • domain assumption Expansion is matter-dominated with a proportional to t^{2/3} throughout the simulation.
    Assumes an inflaton-oscillation-dominated phase, as stated in Section 3.
  • ad hoc to paper The precise form of the initial quantum fluctuations does not qualitatively affect the results.
    Asserted in Section 3 but not demonstrated by varying the fluctuation spectrum.

how reviews work

0 comments
Cite this review

Pith. "Pith review of PQ-ball and its Real Scalar Analogue in an Expanding Universe." pith.science (2026). https://pith.science/paper/5M3VTMKD

@misc{pith2026250701082,
  author       = {Pith},
  title        = {Pith review of: PQ-ball and its Real Scalar Analogue in an Expanding Universe},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5M3VTMKD}},
  note         = {Machine review of arXiv:2507.01082}
}
read the original abstract

We demonstrate the formation of quasi-stable localized scalar configurations in spontaneously symmetry breaking U(1) model by 3+1-dimensional classical lattice simulations. Such configurations are called PQ-balls, as the primary motivation of this kind of configuration is Peccei-Quinn theory under the kinetic misalignment mechanism. Our numerical simulations demonstrate that they can form if the PQ charge is generated through the coherent rotation of a complex scalar field in the complex plane, via dynamics analogous to the Affleck-Dine mechanism. These configurations subsequently decay due to the U(1)-breaking effect induced by spontaneous symmetry breaking. We also demonstrate the formation and decay of oscillons in a similar setup in a real scalar field theory.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

38 extracted references · 10 canonical work pages

  1. [12]

    Kobayashi, K

    J. Kobayashi, K. Nakayama and M. Yamada, Decay rate of PQ-ball , 2504.13510

    arXiv

  2. [1]

    S. R. Coleman, Q-balls, Nucl. Phys. B 262 (1985) 263

    1985

  3. [2]

    Affleck and M

    I. Affleck and M. Dine, A New Mechanism for Baryogenesis , Nucl. Phys. B 249 (1985) 361

    1985

  4. [3]

    M. Dine, L. Randall and S. D. Thomas, Baryogenesis from flat directions of the supersymmetric standard model, Nucl. Phys. B 458 (1996) 291 [ hep-ph/9507453]. – 10 – 0 2000 4000 6000 8000 10 000 0 2000 4000 6000 8000 10 000 Figure 6. Time evolution of total energy density inside oscillons for the case with v = 0, 0.01, 0.1, 0.2, 0.3, from top to bottom

    arXiv 1996

  5. [4]

    Kusenko, Small Q balls , Phys

    A. Kusenko, Small Q balls , Phys. Lett. B 404 (1997) 285 [ hep-th/9704073]

    arXiv 1997

  6. [5]

    Kusenko, Solitons in the supersymmetric extensions of the standard model , Phys

    A. Kusenko, Solitons in the supersymmetric extensions of the standard model , Phys. Lett. B 405 (1997) 108 [ hep-ph/9704273]

    arXiv 1997

  7. [6]

    Kusenko and M

    A. Kusenko and M. E. Shaposhnikov, Supersymmetric Q balls as dark matter , Phys. Lett. B 418 (1998) 46 [ hep-ph/9709492]

    arXiv 1998

  8. [7]

    Enqvist and J

    K. Enqvist and J. McDonald, Q balls and baryogenesis in the MSSM , Phys. Lett. B 425 (1998) 309 [hep-ph/9711514]

    arXiv 1998

Show all 38 references
  1. [8]

    Enqvist and J

    K. Enqvist and J. McDonald, B - ball baryogenesis and the baryon to dark matter ratio , Nucl. Phys. B 538 (1999) 321 [ hep-ph/9803380]

    1999 arXiv

  2. [9]

    Kasuya and M

    S. Kasuya and M. Kawasaki, Q ball formation through Affleck-Dine mechanism , Phys. Rev. D 61 (2000) 041301 [ hep-ph/9909509]

    2000 arXiv

  3. [10]

    Kasuya and M

    S. Kasuya and M. Kawasaki, Q Ball formation in the gravity mediated SUSY breaking scenario , Phys. Rev. D 62 (2000) 023512 [ hep-ph/0002285]

    2000 arXiv

  4. [11]

    Kasuya and M

    S. Kasuya and M. Kawasaki, A New type of stable Q balls in the gauge mediated SUSY breaking, Phys. Rev. Lett. 85 (2000) 2677 [ hep-ph/0006128]

    2000 arXiv

  5. [13]

    Kawasaki, K

    M. Kawasaki, K. Murai and F. Takahashi, Q-balls Under Spontaneously Broken U(1) Symmetry, 2504.13509

    arXiv

  6. [14]

    Nugaev, A

    E. Nugaev, A. Shkerin and M. Smolyakov, Q-holes, JHEP 12 (2016) 032 [ 1609.05568]

    2016 arXiv

  7. [15]

    E. Y. Nugaev and A. V. Shkerin, Review of Nontopological Solitons in Theories with U (1)-Symmetry, J. Exp. Theor. Phys. 130 (2020) 301 [ 1905.05146]

    2020 arXiv

  8. [16]

    R. T. Co, L. J. Hall and K. Harigaya, Axion Kinetic Misalignment Mechanism , Phys. Rev. Lett. 124 (2020) 251802 [ 1910.14152]

    2020 arXiv

  9. [17]

    Chang and Y

    C.-F. Chang and Y. Cui, New Perspectives on Axion Misalignment Mechanism , Phys. Rev. D 102 (2020) 015003 [ 1911.11885]

    2020 arXiv

  10. [18]

    R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons , Phys. Rev. Lett. 38 (1977) 1440

    1977

  11. [19]

    R. D. Peccei and H. R. Quinn, Constraints Imposed by CP Conservation in the Presence of Instantons, Phys. Rev. D 16 (1977) 1791. – 11 –

    1977

  12. [20]

    N. Abe, T. Moroi and M. Yamaguchi, Anomaly mediated supersymmetry breaking with axion , JHEP 01 (2002) 010 [ hep-ph/0111155]

    2002 arXiv

  13. [21]

    Nakamura, K.-i

    S. Nakamura, K.-i. Okumura and M. Yamaguchi, Axionic Mirage Mediation , Phys. Rev. D 77 (2008) 115027 [ 0803.3725]

    2008 arXiv

  14. [22]

    Moroi, K

    T. Moroi, K. Mukaida, K. Nakayama and M. Takimoto, Axion Models with High Scale Inflation, JHEP 11 (2014) 151 [ 1407.7465]

    2014 arXiv

  15. [23]

    G. N. Felder and I. Tkachev, LATTICEEASY: A Program for lattice simulations of scalar fields in an expanding universe , Comput. Phys. Commun. 178 (2008) 929 [ hep-ph/0011159]

    2008 arXiv

  16. [24]

    Kawasaki and K

    M. Kawasaki and K. Saikawa, Study of gravitational radiation from cosmic domain walls , JCAP 09 (2011) 008 [ 1102.5628]

    2011 arXiv

  17. [25]

    Hiramatsu, M

    T. Hiramatsu, M. Kawasaki and K. Saikawa, On the estimation of gravitational wave spectrum from cosmic domain walls , JCAP 02 (2014) 031 [ 1309.5001]

    2014 arXiv

  18. [26]

    Kamada and M

    A. Kamada and M. Yamada, Gravitational waves as a probe of the SUSY scale , Phys. Rev. D 91 (2015) 063529 [ 1407.2882]

    2015 arXiv

  19. [27]

    Kamada and M

    A. Kamada and M. Yamada, Gravitational wave signals from short-lived topological defects in the MSSM , JCAP 10 (2015) 021 [ 1505.01167]

    2015 arXiv

  20. [28]

    I. L. Bogolyubsky and V. G. Makhankov, On the Pulsed Soliton Lifetime in Two Classical Relativistic Theory Models, JETP Lett. 24 (1976) 12

    1976

  21. [29]

    Segur and M

    H. Segur and M. D. Kruskal, Nonexistence of Small Amplitude Breather Solutions in ϕ4 Theory, Phys. Rev. Lett. 58 (1987) 747

    1987

  22. [30]

    Gleiser, Pseudostable bubbles , Phys

    M. Gleiser, Pseudostable bubbles , Phys. Rev. D 49 (1994) 2978 [ hep-ph/9308279]

    1994 arXiv

  23. [31]

    E. J. Copeland, M. Gleiser and H. R. Muller, Oscillons: Resonant configurations during bubble collapse, Phys. Rev. D 52 (1995) 1920 [ hep-ph/9503217]

    1995 arXiv

  24. [32]

    Kasuya, M

    S. Kasuya, M. Kawasaki and F. Takahashi, I-balls, Phys. Lett. B 559 (2003) 99 [hep-ph/0209358]

    2003 arXiv

  25. [33]

    Gleiser and D

    M. Gleiser and D. Sicilia, Analytical Characterization of Oscillon Energy and Lifetime , Phys. Rev. Lett. 101 (2008) 011602 [ 0804.0791]

    2008 arXiv

  26. [34]

    P. M. Saffin, P. Tognarelli and A. Tranberg, Oscillon Lifetime in the Presence of Quantum Fluctuations, JHEP 08 (2014) 125 [ 1401.6168]

    2014 arXiv

  27. [35]

    Mukaida and M

    K. Mukaida and M. Takimoto, Correspondence of I- and Q-balls as Non-relativistic Condensates, JCAP 08 (2014) 051 [ 1405.3233]

    2014 arXiv

  28. [36]

    Mukaida, M

    K. Mukaida, M. Takimoto and M. Yamada, On Longevity of I-ball/Oscillon , JHEP 03 (2017) 122 [1612.07750]

    2017 arXiv

  29. [37]

    M. Ibe, M. Kawasaki, W. Nakano and E. Sonomoto, Decay of I-ball/Oscillon in Classical Field Theory, JHEP 04 (2019) 030 [ 1901.06130]

    2019 arXiv

  30. [38]

    Zhang, M

    H.-Y. Zhang, M. A. Amin, E. J. Copeland, P. M. Saffin and K. D. Lozanov, Classical Decay Rates of Oscillons , JCAP 07 (2020) 055 [ 2004.01202]. – 12 –

    2020 arXiv

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.