{"id":"dae5f6a9-0399-4a77-b2b1-cf2afe553f9d","arxiv_id":"2507.01082","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"3+1D lattice simulations show PQ-balls form from a rotating complex scalar and decay faster for stronger symmetry breaking, with lifetimes roughly scaling as v^4.","lead":"This paper uses computer simulations to show that rotating scalar fields in the early universe can clump into localized objects called PQ-balls, which later decay. The same behavior appears for real scalar fields as oscillons, which matters for axion and dark matter models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed v^4 lifetime scaling is not established: the measured 'all PQ-balls disappear' times combine the v^{-2} threshold time with the intrinsic v^{-4} decay, and the three data points are bracketed between these scalings.","rationale":"The paper's main novel result is the numerical demonstration of PQ-ball formation in 3+1D and the v-dependence of their decay. The formation and qualitative decay are convincingly shown: v = 0 preserves charge, v > 0 leads to decay, and larger v decays faster. However, the quantitative claim that the lifetime follows Eq. (2.13)'s v^{2k2/k1} scaling is the load-bearing part of the abstract and Section 3. The reported observable—the time when all PQ-balls in the box have disappeared—is not the same as the single-ball intrinsic decay time in Eq. (2.13). It includes the expansion time until the background density drops below Φ∞, which scales as v^{-2} in the benchmark case, and the subsequent decay, which scales as v^{-4}. With only three v values, the data are consistent with a slope between -2 and -4, so the v^4 index is not established. This is a correctness risk for the headline claim, not a mere lack of convergence checks. A targeted analysis separating t_th from the intrinsic decay rate, as described in the concrete test, would settle it. I agree with the reader's CONDITIONAL verdict; no change is needed.","tokens_in":8245,"tokens_out":8175,"duration_ms":200373,"concrete_test":"For each v, separately measure the threshold time t_th (when the volume-averaged |Φ| equals Φ∞ from Eq. 2.7) and the subsequent decay of a tracked individual PQ-ball: fit its charge Q_in(t) to an exponential after t_th and extract the intrinsic rate Γ. Then (i) verify Γ is independent of t_th, (ii) check Γ against Eq. (2.13) and its v^4 scaling using the measured Φ0, and (iii) confirm condition (2.14), Γ << H(t_th), holds; if not, t_d is threshold-dominated. Additionally, add v = 0.15 and v = 0.5 runs and fit t_d vs v: a slope near -2 supports threshold dominance; a slope near -4 supports intrinsic decay.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 reports total disappearance times t_d = 10^4, 4×10^3, 1.5×10^3 for v = 0.2, 0.3, 0.4 and states they are 'marginally consistent' with the v^{2k2/k1} = v^4 scaling of Eq. (2.13). However, t_d is not the intrinsic decay time of a PQ-ball. Decay begins only when the background amplitude falls below the v-dependent threshold Φ∞ (Eq. 2.7). In a matter-dominated universe with Φ ∝ t^{-1}, the threshold time is t_th ~ t0/Φ∞ ∝ v^{-k2/k1} = v^{-2} (for k2/k1 = 2). The observed ratios t_d(0.2)/t_d(0.3) = 2.5 and t_d(0.3)/t_d(0.4) = 2.7 lie between the v^{-2} expectation (2.25, 1.78) and the v^{-4} expectation (5.06, 3.16), so the data do not discriminate the intrinsic decay scaling from a threshold-dominated scaling. Condition (2.14), which would justify neglecting the threshold contribution, is never checked for the benchmark parameters. Thus the central quantitative claim—lifetime scaling as v^{2k2/k1}—is not yet supported by the published diagnostics; the qualitative v-dependence is, but the power-law index is unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8598,"tokens_out":8393,"duration_ms":91094,"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":[{"comment":"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.","section":"§3; Eqs. (2.7), (2.13), (2.14)"},{"comment":"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.","section":"§3, benchmark parameters and Fig. 3"}],"minor_comments":[{"comment":"The sentence 'we discuss the cosmological implications of our finidings' contains a typo: 'finidings' should be 'findings'.","section":"§1"},{"comment":"The phrase 'FLR W Universe' appears to be a typo for 'FLRW Universe'.","section":"§3, first sentence"},{"comment":"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.","section":"Eq. (2.13)"},{"comment":"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.","section":"Figs. 3 and 6"},{"comment":"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.","section":"§3, initial fluctuations"}],"recommendation":"major_revision","confidential_remarks":"The paper's qualitative claims—formation and subsequent decay of PQ-balls in 3+1D, and analogous oscillon behavior in a real scalar theory—are plausible and publishable after revision. The main technical issue is that the quantitative v^4 lifetime scaling is not supported by the data: the reported disappearance times include the threshold time, condition (2.14) is not checked, and there are no convergence or uncertainty estimates. The authors should either provide a cleaner intrinsic-lifetime measurement or explicitly downgrade the power-law claim to qualitative. No concerns about novelty or attribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the genuinely new thing: this is the first 3+1D lattice simulation of PQ-ball formation and decay in an expanding universe, without spherical symmetry. The authors start from a rotating complex scalar with small perturbations, see Q-ball-like objects form, and see them decay when v≠0 while remaining stable at v=0. That is a clean numerical demonstration. The real-scalar oscillon analogue is also new, and the argument that it is nontrivial is correct: the field crosses the SSB barrier every cycle, so the similarity to the complex case is not obvious.\n\nThe main soft spot is quantitative. The paper claims the lifetime scales as v^{2k2/k1}=v^4 from Eq. (2.13), but the data are total-disappearance times, which include the time for the background amplitude to fall to the threshold Φ∞. That threshold time scales as v^{-2}. The reported ratios (2.5 and 2.7) are between the v^{-2} and v^{-4} expectations, so the three points do not actually test the intrinsic decay scaling. Condition (2.14), which would let you ignore the threshold contribution, is never checked for the benchmark parameters. The authors themselves only say 'marginally consistent', which is honest, but it means the central quantitative result is unverified.\n\nThe other gaps are standard but important: no convergence tests on lattice spacing, box size, or time step; lifetimes reported without uncertainties; no code or data release. These are fixable. The circularity concern is minor—the analytic formulas come from the authors' own prior work, but the simulation is an independent test, not a fit.\n\nWhat I like beyond the headline: the non-monotonic charge evolution, with small PQ-balls decaying and their charge being absorbed by larger ones, is a real dynamical effect and worth understanding. The oscillon section is brief but reasonable.\n\nWho should read this: anyone working on axion kinetic misalignment, Q-balls, or oscillon formation. It deserves a serious referee, but the referee should ask for a direct measurement of the intrinsic decay rate (or a check of Eq. 2.14) and at least one resolution/box-size cross-check. I'd cite it for the formation demo, not for the scaling.","headline":"Useful 3+1D demonstration of PQ-ball formation/decay, but the v^4 lifetime scaling claim isn't supported by the published numbers.","tokens_in":9088,"tokens_out":3097,"would_cite":true,"duration_ms":34700,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["PQ-ball","Q-ball","oscillon","spontaneous symmetry breaking","lattice simulation","axion kinetic misalignment","Affleck-Dine mechanism","expanding universe"],"falsifier":"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).","tokens_in":8066,"feed_emoji":"🌀","tokens_out":7557,"duration_ms":76101,"temperature":0.7,"pith_summary":"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.","feed_headline":"PQ-balls form and decay; lifetimes shrink as symmetry scale grows","feed_subtitle":"3D lattice runs trace rotating complex scalars into localized lumps that later dissolve, with lifetimes set by v^4.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Provides the spherically symmetric PQ-ball solution and the analytical decay-rate estimate (Eqs. 2.11-2.13) that the 3D simulations are compared against.","marker":"[12]"},{"why":"Coleman's Q-ball paper defines the stable, charge-conserving localized configurations whose existence and formation the PQ-ball generalizes.","marker":"[1]"},{"why":"Affleck-Dine mechanism supplies the physical setup: a coherently rotating scalar field that generates a global charge.","marker":"[2]"},{"why":"Kasuya-Kawasaki established numerically that rotating fields with perturbations fragment into Q-balls, the mechanism the paper's simulations verify for PQ-balls.","marker":"[9]"},{"why":"Kinetic misalignment mechanism is the cosmological motivation; PQ-ball formation and decay determine the thermalization timescale of the PQ field.","marker":"[16]"},{"why":"Supplies the lattice simulation method for seeding initial quantum fluctuations used in the runs.","marker":"[23]"},{"why":"Gives the oscillon formalism for real scalar fields that the paper's real-scalar analogue builds on.","marker":"[31]"}],"fun_headline_variants":["PQ-balls and oscillons form, then decay faster as symmetry breaking grows","Rotating scalars form PQ-balls and oscillons; bigger symmetry scale shortens lifetimes","PQ-balls and oscillons: formation and decay tied to symmetry-breaking scale","Symmetry breaking size sets PQ-ball and oscillon lifetimes","From rotation to decay: PQ-balls and oscillons in an expanding universe"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["PQ-balls and oscillons form, then decay faster as symmetry breaking grows","Rotating scalars form PQ-balls and oscillons; bigger symmetry scale shortens lifetimes","PQ-balls and oscillons: formation and decay tied to symmetry-breaking scale","Symmetry breaking size sets PQ-ball and oscillon lifetimes","From rotation to decay: PQ-balls and oscillons in an expanding universe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000882,"raw_usage":{"total_tokens":3778,"prompt_tokens":881,"completion_tokens":2897,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2797}},"tokens_in":497,"tokens_out":2897,"duration_ms":20544,"temperature":1.0,"reasoning_tokens":2797,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:00:10.925222+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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).","supporting_citations":[],"review_version":1}