Pith. sign in

REVIEW 3 major objections 3 minor 2 cited by

Bubble-wall collisions produce heavy particles at a rate set by a single universal coefficient, V′(2vφ), with a χ⁻⁴ spectrum.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-01 21:03 UTC pith:6IU5SRFH

load-bearing objection Universal χ^{-4} tail with coefficient V'(2vφ) is a sharp, likely-correct correction to earlier bubble-collision particle-production spectra, but the universal coefficient rests on an unvalidated self-similar ansatz; referee it and ask for a full-field check. the 3 major comments →

arxiv 2607.16398 v1 pith:6IU5SRFH submitted 2026-07-17 hep-ph astro-ph.CO

Particle productions during collisions of highly boosted bubble walls

Haipeng An , Hongyi Jiang This is my paper
classification hep-ph astro-ph.CO
keywords bubble collisionsfirst-order phase transitionparticle productionultraviolet spectrumsuperheavy dark mattertrapping equationspectral leakagebaryogenesis
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Collisions of highly boosted bubble walls during a first-order phase transition can produce particles far heavier than the transition's own energy scale, a proposed source of superheavy dark matter and of baryogenesis. This paper claims that the production rate is governed by a universal high-frequency tail: the spectral density falls as χ⁻⁴, with a coefficient fixed solely by V′(2vφ), the derivative of the scalar potential at twice the true-vacuum expectation value. The result holds for both elastic and inelastic collisions, and the paper verifies it with high-precision numerical solutions of the trapping equation while carefully suppressing spectral leakage from finite integration. From this tail it derives analytic rates for heavy scalars and fermions, cosmological yields, and an order-one suppression when the finite bubble radius in 3+1 dimensions is included. If correct, the prediction sharpens—and in places revises—estimates for superheavy dark matter and baryogenesis from bubble collisions.

Core claim

The paper establishes that in the deep-ultraviolet regime χ = ω²−k² ≫ M_h², the Fourier-transformed scalar field profile after two highly boosted bubble walls collide is φ~(χ) = −2V′(2vφ) χ⁻² + O(χ⁻³), for both elastic and inelastic collisions. Consequently the spectral density used in particle-production formulas, F(χ), scales as [V′(2vφ)]² χ⁻⁴, with a logarithmic factor from the kinematic cutoff. Physically, superheavy particles are produced predominantly at the instant the walls first touch—the moment the field reaches φ ≈ 2vφ—so the rest of the potential's shape is irrelevant at leading order. A rigorous derivation via repeated integration by parts in the Bessel-operator representation,

What carries the argument

The trapping equation for the post-collision profile h(s) inside the future light cone, (∂_s² + s⁻¹∂_s)h = −V′(vφ+vφ h)/vφ, together with the eigenfunction identity of the modified Bessel function K0(−is√χ) under the Bessel operator L_B = ∂_s² + s⁻¹∂_s. This identity lets the Fourier transform be expressed as a boundary term plus a series in inverse powers of χ; the boundary term at s = 0 produces the coefficient V′(2vφ). The decay h̄(s) ~ s^{−1/2} cos(M_h s + φ) guarantees the remainder is suppressed.

Load-bearing premise

The analysis assumes the post-collision field inside the forward light cone is exactly the self-similar trapping-equation solution with initial value φ(0)=2vφ at the moment of contact; if a real collision sets a different field value at contact or violates this profile, the universal coefficient V′(2vφ) is not guaranteed.

What would settle it

A full (1+1)-dimensional numerical solution of the two-bubble field equation (without imposing the self-similar ansatz), computing |φ~(ω,k)|² on the hyperboloid ω²−k² = χ for χ ≫ M_h², would settle the claim: the tail must scale as [V′(2vφ)]² χ⁻⁴ with the logarithmic prefactor, and any exponent differing from −4 or coefficient differing from V′(2vφ) would falsify it.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Heavy-particle production rates (m ≫ M_h) are threshold-dominated and depend on the scalar potential only through V′(2vφ), with a universal log(2Λ) factor.
  • The elastic-collision result revises earlier treatments: the field relaxes from 2vφ over a finite time, so the UV tail is χ⁻⁴, not χ⁻².
  • For fermion pairs, the production per unit area is N_ψ/A ≈ [V′(2vφ)]² y_ψ² log(2Λ_ψ) / (35π³ (4m_ψ²)²).
  • Cosmological yields for superheavy fermions follow Y_ψ ∝ [(M_f²+2M_t²)/(3m_ψ²)]² (v_φ/M_Pl) log(E_max/m_ψ).
  • In 3+1 dimensions, the finite bubble radius contributes only an order-one suppression (~0.65) relative to the parallel-wall approximation, so the planar computation remains the leading estimate.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the universal coefficient holds, predictions for superheavy dark matter depend only on one number per potential, V′(2vφ); model comparisons reduce to computing that derivative, which is directly measurable in any given vacuum.
  • The χ⁻⁴ tail is steeper than some earlier claims (χ⁻² in elastic collisions), so abundance estimates for superheavy particles could shift by orders of magnitude depending on threshold and cutoff—worth re-evaluating in existing dark-matter and leptogenesis models.
  • The localization at the collision instant suggests a possible burst-like signature in multi-particle final states; one testable extension is to compute the two-particle correlation or angular distribution at high masses, which the single-field formula does not fully specify.
  • A direct falsification could come from a full (1+1)-dimensional lattice simulation of the original field equation (2.1) without the self-similar ansatz; if the measured UV tail differs from χ⁻⁴, the trapping-equation reduction is the culprit.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 3 minor

Summary. This paper studies the ultraviolet tail of the scalar field spectrum in collisions of highly boosted bubble walls. Working in 1+1 dimensions, the authors adopt the self-similar ansatz (2.4) in which the post-collision field inside the future light cone is φ(s)=vφ(1+h(s)) with h solving the trapping equation (2.3). They derive the asymptotic expansion (2.35)-(2.37), obtaining φ~(χ) = -2V'(2vφ)χ^{-2}+O(χ^{-3}) for χ≡ω^2-k^2≫M_h^2, and hence a spectral density F(χ) ∝ [V'(2vφ)]^2 χ^{-4} for both elastic and inelastic collisions. They verify this against numerical solutions of the trapping ODE with careful spectral-leakage suppression, derive production rates for heavy scalars and fermions with cosmological yields, and extend the result to 3+1 dimensions, finding an O(1) suppression from finite bubble radius.

Significance. If correct, the result is significant: it replaces earlier scaling assumptions in the literature, makes a parameter-free prediction for the leading high-mass production amplitude depending only on V'(2vφ), and has direct phenomenological consequences for superheavy dark matter and leptogenesis. The paper is transparent about its ODE-level numerical check, includes a formal (if not entirely clean) proof of the asymptotic expansion in Appendix B, and carefully handles spectral leakage. The main risk is not internal consistency but the physical validity of the self-similar ansatz, which is imported from Ref. [24] and not tested against the full 1+1D field equation.

major comments (3)
  1. [§2.1, Eq. (2.4), and §3] The universal coefficient in Eq. (2.37) is exactly the s=0 boundary term of the self-similar ansatz. The numerical 'verification' in Section 3 solves the same trapping ODE (2.22) and evaluates the same integral (2.30); it never simulates the full field equation (2.1) for two colliding walls. Thus it is a consistency check of the ODE, not an independent test of the physical premise that a real high-boost collision sets φ(0)=2vφ and evolves self-similarly. Finite wall width, finite γ_w, and nonlinearities near the collision point could shift the coefficient without changing the χ^{-4} scaling. Please add a full 1+1D simulation of Eq. (2.1) for finite-boost walls, or a controlled derivation of the ansatz, or state the result as explicitly conditional.
  2. [Eq. (2.36)] The term -2∆(\bar h(0)) is missing the factor vφ that appears in the decomposition (2.19), where V'(φ)=M_h^2 vφ \bar h + vφ ∆(\bar h). As printed, with -2∆(\bar h(0)), the simplification to Eq. (2.37) does not follow unless vφ=1. The correct coefficient is -2vφ∆(\bar h(0)), giving Eq. (2.37). This is a typo, but in the central derivation it should be corrected.
  3. [Appendix B, Corollary 2 (B.50)] The stated expansion I(χ)=Σ_{n=0}^{N-1} (-1)^n [(\partial_s^2+s^{-1}\partial_s)^n ∆]_{s=0}/χ^{n+1}+O(χ^{-N-1}) is inconsistent with Eq. (2.35) and with Lemma 4 when α=i√χ: the leading term should be -2∆(\bar h(0))/χ and the general term should be -2(-1)^n [\hat L_B^n ∆]_{s=0}/χ^{n+1}. As written, the appendix's 'rigorous proof' supports a different series than the main text. Please correct the factor/sign and check the subsequent equations.
minor comments (3)
  1. [Figure 2] The sub-panel labels appear swapped: panel (a) labeled 'Inelastic collisions' lists a=13,17,21,25, while Table 1 and the text identify these as elastic parameters; panel (b) has the inelastic set. Please verify the labels.
  2. [§5.2, Eq. (5.26)] The replacement of the oscillatory double integral by the R≃R' contribution is stated without an explicit error estimate. For R0≫√(r1 r2) the stationary-phase approximation may need justification.
  3. [§4.3] The conversion from production per area to number density uses A/V≃3/R* with R* from Ref. [20]; the uncertainty from this geometric approximation is not discussed. A brief comment would help.

Circularity Check

0 steps flagged

No circularity found; the central UV coefficient is a derived consequence of the stated trapping equation and boundary conditions, and the numerical check is an explicit consistency test of the same ODE.

full rationale

The paper derives Eq. (2.37) by taking the stated self-similar ansatz (2.4) and the trapping equation (2.3) with initial conditions φ(0)=2vφ, then performing an asymptotic expansion of the Fourier-space integral. The coefficient −2V′(2vφ) arises from the boundary term at s=0 through the equation of motion, not from fitting any parameter. The numerical section explicitly solves the same trapping equation (2.22) and evaluates the same integral (2.30), so the agreement is a consistency check of the asymptotic expansion rather than an independent test of the physical ansatz; the paper labels this clearly. Ref. [24] is an external citation that supplies the trapping-equation framework, not a self-citation, and Refs. [3,7,8] are self-citations used only for background, not as load-bearing justification. No step reduces a prediction to an input by construction, and no fitted quantity is renamed as a prediction. The main limitation—that the ansatz itself is not validated against a full 1+1D lattice simulation of Eq. (2.1)—is a correctness/validity concern, not a circularity in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The central UV tail has no fitted constants; the coefficient V'(2vφ) is fixed by the potential. The numerical parameters (N=2^24, grid spacings, ε=100/L) are computational choices, not physical inputs; the result is stable across three grids. The main assumptions are the self-similar trapping ansatz, the O(s⁻¹/²) decay lemmas, polynomial expansion of V, the spectral ansatz for Im Γ^(2), and the step-function bubble idealization in 3+1D.

axioms (5)
  • domain assumption The post-collision field in the forward light cone is exactly of the self-similar form Eq. (2.4), with h(s) solving the trapping equation (2.3) and initial condition φ(0)=2vφ, ∂_sφ(0)=0.
    This is the model-to-physics bridge. It is taken from Ref. [24] and is not tested in this paper against a full 1+1D field-theory collision; if real collisions deviate from Eq. (2.4), the derived χ⁻⁴ tail is not universal.
  • standard math Asymptotic decay \bar h(s), \bar h'(s)=O(s⁻¹/²) (Appendix A) and the B₄-algebra estimates (Appendix B) that justify the asymptotic expansion of I₀(χ).
    These are mathematical lemmas used to drop subleading terms in the integration-by-parts expansion; they rely on known results for asymptotically autonomous systems and elliptic regularity, and the proof is only sketched in places.
  • domain assumption V admits a polynomial expansion around the vacua and around 2vφ with V'(0)=V'(vφ)=0 (Eq. 2.2 is used for numerics).
    The universal formula (2.37) needs V' to have a Taylor expansion at φ(0)=2vφ; the paper states this requirement in Sec. 2.4.
  • domain assumption The absorptive part Im Γ^(2)(χ) is parametrized as Γ₀ χ^{αΓ}(1−M²/χ)^{βΓ} Θ(χ−M²) when computing production rates (Eq. 4.3).
    This is a standard spectral-function ansatz, not derived here; it sets the threshold behavior used in Eqs. (4.7), (4.18).
  • domain assumption In 3+1D, bubble walls are idealized as characteristic functions / step functions, with expansion-part singularities suppressed by 1/R (Sec. 5, Appendix C).
    The geometric factor and order-one suppression in Sec. 5 rely on this idealization; deviations for finite wall width are not quantified.

pith-pipeline@v1.3.0-alltime-deepseek · 29267 in / 16182 out tokens · 169308 ms · 2026-08-01T21:03:21.632139+00:00 · methodology

0 comments
read the original abstract

We investigate the production of particles much heavier than the characteristic scale of a cosmological first-order phase transition through collisions of highly boosted bubble walls. Using the scalar order-parameter field, we derive the ultraviolet behavior of its Fourier-space profile for both elastic and inelastic collisions. In the regime $\chi\equiv\omega^2-\mathbf{k}^2\gg M_h^2$, we find the universal result $\tilde{\phi}(\chi) = -2V^\prime(2v_\phi)\chi^{-2}+O(\chi^{-3}),$ implying that the spectral density scales as $F(\chi)\propto [V^\prime(2v_\phi)]^2\chi^{-4}$. Thus, heavy-particle production is localized near the instant of collision and, at leading order, depends on the scalar potential only through $V^\prime(2v_\phi)$. We verify this behavior using high-precision numerical solutions of the trapping equation, carefully suppressing spectral leakage from the finite integration domain, and obtain agreement over a broad ultraviolet range. We then derive analytical production rates for general heavy-particle thresholds and for fermion pairs, together with their cosmological number density and yield. Finally, we extend the analysis to $(3+1)$ dimensions and incorporate the finite bubble radius, finding an order-one suppression relative to the parallel-wall approximation. Our results revise the ultraviolet scaling used in previous treatments and have direct implications for superheavy dark-matter production and baryogenesis from bubble collisions.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Can the universe be matter-dominated after a supercooled first-order phase transition?

    hep-ph 2026-07 conditional novelty 7.0

    After a supercooled first-order phase transition, the scalar field's equation of state is set by the bubble-wall Lorentz factor γ*, and matter domination is delayed until a/a* ≃ γ* in the free-streaming limit.

  2. Particle Production via Rippled Bubble Walls

    hep-ph 2026-07 conditional novelty 6.0

    A rippled bubble wall produces heavy particles resonantly when the momentum transfer matches the ripple frequency, potentially raising dark-matter abundance by orders of magnitude.

Reference graph

Works this paper leans on

29 extracted references · 19 linked inside Pith · cited by 2 Pith papers

  1. [1]

    Witten,Cosmic Separation of Phases,Phys

    E. Witten,Cosmic Separation of Phases,Phys. Rev. D30(1984) 272–285

    1984

  2. [2]

    A. D. Sakharov,Violation of CP Invariance, C Asymmetry, and Baryon Asymmetry of the Universe,JETP Lett.5(1967) 24–27. [Pisma Zh. Eksp. Teor. Fiz. 5 (1967) 32]

    1967

  3. [3]

    H. An, X. Tong, and S. Zhou,Superheavy dark matter production from a symmetry-restoring first-order phase transition during inflation,Phys. Rev. D107(2023), no. 2 023522, [arXiv:2208.14857]

    Pith/arXiv arXiv 2023

  4. [4]

    Freese and M

    K. Freese and M. W. Winkler,Dark matter and gravitational waves from a dark big bang,Phys. Rev. D107(2023), no. 8 083522, [arXiv:2302.11579]

    Pith/arXiv arXiv 2023

  5. [5]

    Cheng and F

    Z. Cheng and F. P. Huang,Dark Matter Production from Bubble Collisions during a First-Order Phase Transition at the End of Inflation,arXiv:2605.03758

    Pith/arXiv arXiv

  6. [6]

    Casey and C

    R. Casey and C. Ilie,Dark sector tunneling field potentials for a dark big bang,Phys. Rev. D 110(2024), no. 10 103522, [arXiv:2407.05752]

    Pith/arXiv arXiv 2024

  7. [7]

    An and T

    H. An and T. Li,A New Origin of the Big Bang from Dark-Sector-Induced Vacuum Decay and Its Gravitational-Wave Signal,arXiv:2606.06587. – 41 –

    Pith/arXiv arXiv

  8. [8]

    H. An, T. Li, and C. Yang,Gravitational Waves and Primordial Black Holes produced by Dark Meta Stable Vacuum Decay,arXiv:2601.14366

    arXiv

  9. [9]

    Falkowski and J

    A. Falkowski and J. M. No,Non-thermal Dark Matter Production from the Electroweak Phase Transition: Multi-TeV WIMPs and ’Baby-Zillas’,JHEP02(2013) 034, [arXiv:1211.5615]

    Pith/arXiv arXiv 2013

  10. [10]

    S. Iso, P. D. Serpico, and K. Shimada,QCD-Electroweak First-Order Phase Transition in a Supercooled Universe,Phys. Rev. Lett.119(2017), no. 14 141301, [arXiv:1704.04955]

    Pith/arXiv arXiv 2017

  11. [11]

    Hambye, A

    T. Hambye, A. Strumia, and D. Teresi,Super-cool Dark Matter,JHEP08(2018) 188, [arXiv:1805.01473]

    Pith/arXiv arXiv 2018

  12. [12]

    Baratella, A

    P. Baratella, A. Pomarol, and F. Rompineve,The Supercooled Universe,JHEP03(2019) 100, [arXiv:1812.06996]

    Pith/arXiv arXiv 2019

  13. [13]

    Baldes, Y

    I. Baldes, Y. Gouttenoire, F. Sala, and G. Servant,Supercool composite Dark Matter beyond 100 TeV,JHEP07(2022) 084, [arXiv:2110.13926]

    Pith/arXiv arXiv 2022

  14. [14]

    G. F. Giudice, H. M. Lee, A. Pomarol, and B. Shakya,Nonthermal heavy dark matter from a first-order phase transition,JHEP12(2024) 190, [arXiv:2403.03252]

    Pith/arXiv arXiv 2024

  15. [15]

    W.-Y. Ai, M. Fairbairn, K. Mimasu, and T. You,Non-thermal production of heavy vector dark matter from relativistic bubble walls,JHEP05(2025) 225, [arXiv:2406.20051]. [Erratum: JHEP 05, 021 (2026)]

    arXiv 2025

  16. [16]

    Shakya,A Cosmic Higgs Collider,arXiv:2512.13815

    B. Shakya,A Cosmic Higgs Collider,arXiv:2512.13815

    arXiv

  17. [17]

    Fairbairn and W

    M. Fairbairn and W. S. A. Shellard,Phenomenology of Vector Dark Matter produced by a First Order Phase Transition,arXiv:2604.24528

    Pith/arXiv arXiv

  18. [18]

    Konstandin and G

    T. Konstandin and G. Servant,Natural Cold Baryogenesis from Strongly Interacting Electroweak Symmetry Breaking,JCAP07(2011) 024, [arXiv:1104.4793]

    Pith/arXiv arXiv 2011

  19. [19]

    Katz and A

    A. Katz and A. Riotto,Baryogenesis and Gravitational Waves from Runaway Bubble Collisions,JCAP11(2016) 011, [arXiv:1608.00583]

    Pith/arXiv arXiv 2016

  20. [20]

    Cataldi and B

    M. Cataldi and B. Shakya,Leptogenesis via bubble collisions,JCAP11(2024) 047, [arXiv:2407.16747]

    Pith/arXiv arXiv 2024

  21. [21]

    Datta and A

    A. Datta and A. Sil,Probing Leptogenesis through Gravitational Waves,arXiv:2410.01900

    Pith/arXiv arXiv

  22. [22]

    Mansour and B

    H. Mansour and B. Shakya,Particle production from phase transition bubbles,Phys. Rev. D 111(2025), no. 2 023520, [arXiv:2308.13070]

    Pith/arXiv arXiv 2025

  23. [23]

    Ghoshal, P

    A. Ghoshal, P. Pal, and A. Strumia,Particle production from bubble collisions, arXiv:2607.15279

    Pith/arXiv arXiv

  24. [24]

    Jinno, T

    R. Jinno, T. Konstandin, and M. Takimoto,Relativistic bubble collisions—a closer look,JCAP 09(2019) 035, [arXiv:1906.02588]

    Pith/arXiv arXiv 2019

  25. [25]

    Watkins and L

    R. Watkins and L. M. Widrow,Aspects of Reheating in First Order Inflation,Nucl. Phys. B 374(1992) 446–468

    1992

  26. [26]

    F. J. Harris,On the use of windows for harmonic analysis with the discrete fourier transform, Proceedings of the IEEE66(1978), no. 1 51–83

    1978

  27. [27]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery,Numerical Recipes: The Art of Scientific Computing. Cambridge University Press, New York, NY, USA, 3rd ed., 2007. – 42 –

    2007

  28. [28]

    Markus,II

    L. Markus,II. ASYMPTOTICALLY AUTONOMOUS DIFFERENTIAL SYSTEMS, pp. 17–30. Princeton University Press, Princeton, 2016

    2016

  29. [29]

    Cabot, H

    A. Cabot, H. Engler, and S. Gadat,On the long time behavior of second order differential equations with asymptotically small dissipation,Transactions of the American Mathematical Society361(2009), no. 11 5983–6017. – 43 –

    2009