{"id":"39560aed-510a-4d52-a7a2-574b605daa26","arxiv_id":"2607.16398","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Bubble-wall collisions produce ultra-heavy particles with a universal spectrum ∝ [V'(2vφ)]²/χ⁴, localized at the collision instant.","lead":"This paper derives a universal high-energy law for particles produced when fast-moving bubble walls collide during a cosmological phase transition: the production spectrum falls as the fourth power of the invariant mass, controlled only by the slope of the scalar potential at twice the vacuum value. The result sharpens predictions for superheavy dark matter and baryogenesis from bubble collisions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central UV coefficient rests on unvalidated self-similar collision ansatz Eq. (2.4); numerical check solves the same ODE and cannot test it.","rationale":"The reader's weakest_assumption identifies the self-similar ansatz Eq. (2.4) as the load-bearing premise, and I agree. The analytical derivation of Eq. (2.37) is internally coherent: the integration-by-parts machinery and Appendix B give a systematic \\chi^{-2} tail for the trapping-equation solution, and the numerical section carefully handles spectral leakage and confirms the ODE-level prediction. However, the jump from the reduced ODE to the actual two-bubble collision is not tested in this paper. The numerical simulations solve the same ODE (2.22) that defines the ansatz, so they cannot validate Eq. (2.4) or the initial value \\phi(0)=2v_\\phi. A full 1+1D simulation would directly test whether the post-collision field is actually a function of s only and whether the UV coefficient matches -2V'(2v_\\phi). If it does not, the central claim's universality and potential-only dependence are unsupported. This does not invalidate the mathematical core, but it justifies the conditional verdict. I also note the Figure 2/Table 1 inelastic/elastic label swap as a separate but non-central issue. Overall, the reader's conditional verdict is appropriate and no adjustment is needed.","tokens_in":29614,"tokens_out":15943,"duration_ms":146304,"concrete_test":"Run a 1+1D lattice simulation of Eq. (2.1) with two boosted kink initial profiles (e.g., tanh walls with Lorentz factor \\gamma=10 and 100) and potential Eq. (2.2), for both inelastic (a=1,3,5,7) and elastic (a=13,17,21,25) parameters. At several fixed t>0, extract the field along x, invert s=\\sqrt{t^2-x^2} to obtain h(s)=[\\phi(t,x)/v_\\phi - 1], and compare with the solution of the trapping ODE (2.3) with \\phi(0)=2v_\\phi. Independently compute the Fourier transform \\tilde\\phi(\\omega,k) on the same lattice, bin by \\chi=\\omega^2-k^2, and fit \\chi^2\\tilde\\phi(\\chi) in the UV; check that the fitted constant equals -2V'(2v_\\phi) for both regimes. If the extracted h(s) or the fitted coefficient deviates significantly from the ansatz prediction, the universality claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result Eq. (2.37) — \\tilde\\phi(\\chi) = -2V'(2v_\\phi)\\chi^{-2} + O(\\chi^{-3}) — is derived from the self-similar ansatz Eq. (2.4): inside the forward light cone the field is taken to be v_\\phi[1+h(s)], with h(s) solving the reduced trapping ODE (2.3) subject to h(0)=1, h'(0)=0. The coefficient -2V'(2v_\\phi) is exactly the s=0 boundary term of this ansatz. Section 3 verifies the \\chi^{-2} tail by numerically solving the same ODE (2.22) and evaluating the same integral (2.30); it never simulates the full 1+1D field equation (2.1) for two colliding walls. Thus the numerical agreement is evidence only for the ODE asymptotics, not for the physical premise that a real collision sets \\phi(0)=2v_\\phi and evolves self-similarly. At finite boost, finite wall width and nonlinear interactions near the collision point could shift the effective field value and source strength, changing the coefficient of \\chi^{-4} without altering the scaling. The paper imports the ansatz from Ref. [24] but does not demonstrate its validity for the potentials and boosts used here. If Eq. (2.4) is not an accurate description of actual bubble collisions, the universal prediction F(\\chi) \\propto [V'(2v_\\phi)]^2\\chi^{-4} lacks direct support.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":29963,"tokens_out":13331,"duration_ms":105803,"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":[{"comment":"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.","section":"§2.1, Eq. (2.4), and §3"},{"comment":"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.","section":"Eq. (2.36)"},{"comment":"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.","section":"Appendix B, Corollary 2 (B.50)"}],"minor_comments":[{"comment":"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.","section":"Figure 2"},{"comment":"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.","section":"§5.2, Eq. (5.26)"},{"comment":"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.","section":"§4.3"}],"recommendation":"major_revision","confidential_remarks":"The paper's main result is appealing and the ODE-level derivation is coherent. My recommendation is driven by the gap between the reduced model and the full collision problem; the numerical verification does not close that gap. If the ansatz can be validated (or its regime of validity sharply characterized), I would support acceptance. Also please fix the Eq. (2.36) factor and the Appendix B factor/sign before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things: the paper's central asymptotic claim is probably right within its model, and the model-to-physics bridge is not tested. The χ^{-4} tail with coefficient V'(2vφ) is a real correction to earlier treatments (χ^{-2} for elastic collisions in Falkowski–No), and the spectral-leakage explanation for the discrepancy with Mansour–Shakya is plausible and well documented.\n\nWhat's genuinely good: the derivation of Eq. (2.37) is coherent; Appendix B is a real asymptotic proof, not a gloss; the numerics are careful about endpoint leakage with the iε window, and the grid-refinement overlap is a legitimate consistency check. The threshold-dominated production integrals in Sec. 4 are worked out cleanly, and the fermion rate matching Ref. [23] up to O(1) is a nice cross-check. The 3+1 geometric suppression is heuristic but a reasonable order-of-magnitude estimate.\n\nThe soft spot is load-bearing. The universal coefficient V'(2vφ) follows from the ansatz φ(s)=vφ[1+h(s)] inside the forward light cone with h(0)=1, h'(0)=0. That ansatz is imported from Ref. [24] and is not independently validated here: the numerical section solves the same reduced trapping ODE, so the agreement confirms only the ODE asymptotics, not that a real collision produces φ(0)=2vφ and evolves self-similarly. If the actual field value at the collision point differs at O(1) for finite wall width or finite boost, the coefficient changes even though the χ^{-4} scaling may survive. So the strong 'model-independent coefficient' statement is conditional. A full 1+1D simulation of Eq. (2.1), even at moderate boost, would largely settle this, and the paper should at least show such a test or argue why the UV tail can't see the bulk dynamics.\n\nMinor: Figure 2's inelastic/elastic labels are swapped relative to Table 1 (a=13–25 appear as 'inelastic' in the figure but 'elastic' in the table). Cosmetic, but it should be fixed before publication.\n\nWho this is for: people computing superheavy DM or leptogenesis from bubble collisions. It deserves a serious referee; the central claim is important and checkable, and the caveats are fixable with a modest amount of additional work. I'd send it to review and ask for an independent full-field simulation as a condition of publication.","headline":"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.","tokens_in":30458,"tokens_out":2711,"would_cite":true,"duration_ms":24307,"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":"Bubble-wall collisions produce heavy particles at a rate set by a single universal coefficient, V′(2vφ), with a χ⁻⁴ spectrum.","keywords":["bubble collisions","first-order phase transition","particle production","ultraviolet spectrum","superheavy dark matter","trapping equation","spectral leakage","baryogenesis"],"falsifier":"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.","tokens_in":29504,"feed_emoji":"⚛️","tokens_out":5070,"duration_ms":40411,"temperature":0.7,"pith_summary":"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.","feed_headline":"One coefficient sets all heavy-particle yields from bubble collisions","feed_subtitle":"Particle yields from supercooled phase transitions depend only on a single derivative of the scalar potential.","key_machinery":"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.","core_discovery":"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,","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["All heavy particle yields from bubble collisions trace to one constant","Single potential derivative sets every superheavy particle yield","Bubble-wall collisions: one number predicts all heavy particle output","Universal χ⁻⁴ scaling ties heavy particles to a single derivative","One coefficient governs heavy particle production in bubble collisions"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["All heavy particle yields from bubble collisions trace to one constant","Single potential derivative sets every superheavy particle yield","Bubble-wall collisions: one number predicts all heavy particle output","Universal χ⁻⁴ scaling ties heavy particles to a single derivative","One coefficient governs heavy particle production in bubble collisions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1327,"prompt_tokens":813,"completion_tokens":514,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":434}},"tokens_in":557,"tokens_out":514,"duration_ms":5685,"temperature":1.0,"reasoning_tokens":434,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T21:03:21.632139+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}