{"id":"c092c90f-160d-4238-81b7-7fe07e7c56c9","arxiv_id":"2509.07076","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"CMB anisotropies from stochastic bubble nucleation constrain late-time phase transitions to release less than ~1% of dark energy when β/H⋆≲25, much tighter than Hubble-budget limits.","lead":"This paper derives new limits on a late-universe first-order phase transition that converts some dark energy into invisible radiation, using the anisotropic imprint it would leave on the cosmic microwave background. The bounds are up to an order of magnitude stronger than from the Hubble expansion rate alone, and rule out certain negative-vacuum-energy phase transitions that would otherwise be allowed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8)'s distance-matching condition is asserted, not derived; if a linear-in-δz_pt response survives, the r≲10^-5(β/H⋆)^2 bound changes.","rationale":"The reader identified the same load-bearing assumption: Eq. (8) is asserted rather than derived, and the central bound depends on the resulting quadratic scaling. I agree that this is the weakest point. The Gaussian truncation in Eq. (A10) is real but less severe, since it only shifts the overall coefficient of the power spectrum, not the functional form of the bound. The concern does not by itself prove the paper wrong: a careful gauge choice or a full perturbative derivation might recover the same result. But the burden is on the authors to show that the first-encounter surface can be treated as a fixed-distance source. Given that the paper otherwise presents a transparent calculation, uses an independently published P_δt [43], and gives order-of-magnitude estimates that can be checked, the appropriate verdict is conditional on this derivation being supplied. Hence I keep the reader's CONDITIONAL verdict, which is represented here as UNCHANGED.","tokens_in":19539,"tokens_out":11251,"duration_ms":116278,"concrete_test":"Compute the CMB temperature perturbation from the same PT-time power spectrum using standard first-order cosmological perturbation theory: introduce the stochastic PT time δt_c(x) as a perturbation to the expansion rate (or as a source in the metric perturbations), integrate the photon geodesic from a fixed comoving last-scattering surface to the observer, and expand to first order in δt_c. If the resulting δT/T contains a term linear in δt_c (equivalently, ∂δz0/∂δz_pt ≠ 0 at δz_pt=0), then Eq. (8) and the quadratic Eq. (9) are not the correct mapping. A simpler analytic check: repeat the Appendix expansion but instead of equating distances from z̄_pt+δz_pt to today, equate distances from z_rec ≈ 1100 to today in the two histories and expand to first order; this directly tests whether the linear term vanishes for the physical source.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central bound (15) follows from δz0 ∝ δz_pt^2 in Eq. (9), which is obtained by imposing Eq. (8): that the comoving distance from a fixed initial redshift z̄_pt+δz_pt to the observer is unchanged by the phase transition. This condition is stated, not derived. CMB photons originate at recombination, whose comoving distance from us is fixed; a spatially varying PT time should source a first-order temperature perturbation (an ISW-type line-of-sight integral) unless a gauge choice removes it. Eq. (8) is not that gauge choice; it equates distances in two different background histories, which is not a consequence of the geodesic equation or of photon number conservation. The Appendix uses Eq. (8) as input, and the fact that the linear term cancels in the expansion is a property of this imposed equality, not of the physical photon redshift. If a linear δz0 ∝ r δz_pt term survives, Pδz0 scales differently—Pδz0 ~ r^2 Pδt rather than r^2(β/H⋆)^{-4} convolutions—and the amplitude and ℓ-shape of D_ℓ^{TT,pt}, hence the r≲10^-5(β/H⋆)^2 bound, change. The Gaussian truncation (A10) is secondary: it affects the coefficient, not the scaling.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers a late-time first-order phase transition in a dark sector during dark-energy domination (z ≲ 0.3), with the latent heat converted into free-streaming dark radiation. It computes CMB temperature anisotropies generated by fluctuations in the redshift at which photons enter the true vacuum, described by a spatially varying δz_pt. Using Planck 2018 data through a simplified three-bin χ² analysis, it derives an approximate upper bound r ≲ 10^{-5}(β/H_⋆)^2 on the released fraction of dark energy, and it discusses transitions to negative vacuum energy (r > 1), concluding that β/H_⋆ ≳ 300 is required and that β/H_⋆ ≲ 500 still leaves at least 14 Gyr before a crunch. The central analytic results are Eq. (9), δz_0 ∝ δz_pt^2, and Eq. (11), with the derivation in Appendix A.","tokens_in":19862,"tokens_out":8814,"duration_ms":108958,"significance":"Late-time FOPTs are a plausible and understudied dark-sector signature, and translating the stochastic completion time into CMB anisotropies is a fresh and potentially useful angle. If the derivation were sound, the bound r ≲ 10^{-5}(β/H_⋆)^2 would be an important, broadly applicable constraint, and the negative-vacuum/crunch discussion is a nice addition. The paper is explicit about several approximations: latent heat as free-streaming radiation, fixed Planck fiducial parameters, a three-bin χ², and the use of a published P_δt(k) from Ref. [43]. It also provides an explicit scaling relation that can be checked numerically. However, the central quantitative claim currently rests on an unproven distance-matching condition and an unjustified Gaussian truncation, so the headline numbers are not yet established.","major_comments":[{"comment":"The equality χ_AO = χ_BO is asserted with the sentence that 'the total comoving distance from the same initial redshift to us must remain the same,' but no derivation is given. In the two histories being compared, H(z) differs over a finite redshift interval; a comoving distance in a modified background is not an invariant quantity, and the photon redshift is determined by the null geodesic equation, not by an independent matching condition. This condition is precisely what removes the linear-in-δz_pt term in the expansion leading to Eq. (9). If a term ∝ r δz_pt survives in δz_0, then P_{δz0}(k) scales as r^2 P_{δt}(k) rather than the convolution in Eq. (12), changing the ℓ-shape and the β/H scaling of the bound in Eq. (15). The appendix uses Eq. (8) as input and therefore does not resolve the issue. The authors should either derive Eq. (8) from a consistent gauge/matching prescription o","section":"Section III, Eq. (8)"},{"comment":"The four-point function is truncated to products of two-point functions with the statement that only cross-correlations between x and y contribute. This is a Gaussian approximation for δz_pt. However, the text in Section III states δz_pt ≥ 0, so δz_pt cannot be an exactly Gaussian field; in a bubble-nucleation process the connected part of the four-point function is generically nonzero and can be comparable to the disconnected part at relevant scales. No estimate or bound on the connected contribution is given. Because the normalization of D_ℓ^{TT,pt}, and hence the 2σ curve in Fig. 3, depends on this truncation, the numerical results are not fully supported.","section":"Appendix A, Eq. (A10)"},{"comment":"The statement that 'for β/H_⋆ ≲ 500, the universe will not crunch for at least 14 Gyr' is obtained by combining the t_end formula with the upper bounds r ≤ 2 and r ≤ 3 from the CMB anisotropy estimate. Those bounds inherit the uncertainties in Eqs. (8) and (A10). The 14-Gyr claim should therefore be presented as conditional on the anisotropy analysis, not as a robust consequence of the phase-transition dynamics alone.","section":"Section IV and Fig. 3"}],"minor_comments":[{"comment":"The text around Eq. (4) contains an inconsistency: the fiducial values are given as ΩΛ = 0.69, Ωm = 0.31 in one place, but the duplicated passage states ΩΛ = 0.31, Ωm = 0.69. The correct values should be used consistently.","section":"Section III, Eq. (4)"},{"comment":"The caption includes stray editorial text ('SK : ¯zpt, also above in Figure with r=0.1') and duplicated figure text. This should be cleaned before submission.","section":"Figure 2 caption"},{"comment":"The three-bin χ² uses Planck error bars that are correlated and derived from a baseline cosmology. The authors acknowledge the limitation, but the resulting upper bound should be described explicitly as an order-of-magnitude estimate rather than a rigorous confidence limit.","section":"Section III, Eq. (14)"},{"comment":"The statement δz_pt ≥ 0 and the subsequent treatment of δz_pt as a Gaussian random field in Appendix A should be reconciled; this is related to the major comment on Eq. (A10).","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The submitted version appears to contain duplicated text blocks and an editorial artifact in the caption of Fig. 2; the editor may want a clean manuscript. The overlap with Ref. [43], which shares an author, is acceptable because that paper is an independent published calculation, but the present manuscript's central new step (the matching condition) is what needs the most scrutiny."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is the first paper I've seen that actually probes first-order phase transitions in the z ≲ 0.3 dark-energy-dominated epoch with CMB anisotropies, and the constraints are not incremental—they beat the Hubble-budget bound by orders of magnitude. The mechanism is genuinely new: stochastic bubble nucleation makes the 'first encounter' of a CMB photon with the true-vacuum region direction-dependent, and the induced redshift shift goes as (δz_pt)². The headline bound, r ≲ 10^-5 (β/H⋆)², would make a completed PT with β/H⋆ ≲ 25 limited to less than 1% of the dark energy—a strong statement and a testable one.\n\nWhat's actually solid: the input P_δt(k) is a published first-principles calculation (ref. [43]). It shares an author with this paper, which you should note, but it's not a fit to the target bound, so the circularity worry is mild. The appendix works the two-point function of the composite operator δz_pt² through explicitly, and the paper is upfront about its approximations: free-streaming dark radiation, fixed Planck parameters, and a three-bin χ² on archive error bars. For an order-of-magnitude exclusion, that's honest and adequate. The r > 1 negative-vacuum section is a nice bonus.\n\nNow the soft spot, and it's the one the referee should attack: Eq. (8) states that the comoving distance from a fixed initial redshift to the observer is unchanged by the PT, and that is what cancels the linear-in-δz_pt response and produces the quadratic scaling. The condition is asserted, not derived from the geodesic equation or from a gauge choice. A spatially varying PT time should normally source a first-order metric perturbation, and an ISW-type integral along the photon path is the natural consequence. The paper's argument that free-streaming radiation homogenizes the density is plausible but doesn't demonstrate that the linear term is absent in a consistent gauge. If a linear term does survive, P_δz0 picks up r² P_δt directly and the bound becomes roughly r ≲ 10^-5 (β/H⋆)—actually stronger, so the qualitative conclusion survives, but the claimed power law and the ℓ-shape don't. The Gaussian truncation in Eq. (A10) is a real but secondary approximation: it affects the coefficient, not the scaling.\n\nWho's this for? Anyone working on dark-sector FOPTs, late-time ISW physics, or DESI-adjacent dark energy models. It deserves a serious referee—send it, and ask specifically for a derivation or explicit justification of Eq. (8). If that condition survives contact with perturbation theory, this is a solid and useful paper.","headline":"Novel constraints on late-time phase transitions via a first-encounter-surface mechanism that are orders of magnitude stronger than the Hubble-budget bound, but the central quadratic scaling rests on an asserted distance-matching condition; send to a referee who will push on it.","tokens_in":20399,"tokens_out":17434,"would_cite":true,"duration_ms":169392,"reading_group":"yes","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 claims current CMB temperature data cap a completed slow dark-sector phase transition at about 1% of the dark energy via anisotropic photon redshifting at a stochastic first-encounter surface.","keywords":["first-order phase transition","dark energy","CMB temperature anisotropy","integrated Sachs-Wolfe effect","bubble nucleation","first-encounter surface","vacuum energy","late-time cosmology"],"falsifier":"One decisive check is to compute or simulate photon geodesics through a stochastic bubble-nucleation history at z≈0.1–0.3 and test the paper's matching condition χ_AO = χ_BO (Eq. 8) at first order in δz_pt: any linear-in-δz_pt contribution to δz0 invalidates the quadratic formula and the derived bound. Observationally, a low-multipole (ℓ<20) CMB temperature measurement with uncertainties smaller than the predicted D_l^TT,pt for a given (r, β/H⋆, z̄_pt) would confirm or exclude the signal; an absence of the predicted r²(β/H⋆)⁻⁴ scaling would likewise settle the claim.","tokens_in":19420,"feed_emoji":"🌌","tokens_out":12787,"duration_ms":134169,"temperature":0.7,"pith_summary":"This paper tries to establish that a first-order phase transition in the dark sector during the last few billion years—when vacuum energy already dominates—is far more constrained by CMB temperature anisotropies than by measurements of the Hubble expansion rate. The mechanism is that bubble nucleation is stochastic, so different lines of sight enter the true vacuum at slightly different redshifts; the resulting corrugated \"first-encounter surface\" makes CMB photons redshift by different amounts across the sky. In the minimal case where the released latent heat becomes free-streaming dark radiation, the authors derive a quadratic relation δz0 ∝ δz_pt² and project the induced redshift power spectrum onto low multipoles. Their central bound is r ≲ 10⁻⁵ (β/H⋆)², where β/H⋆ measures the transition's inverse duration in units of the Hubble rate at completion: a transition with β/H⋆ ≲ 25 can release no more than about 1% of the dark energy. This matters because Hubble-expansion measurements alone miss everything below roughly ten percent of the dark energy, so the CMB turns an otherwise nearly invisible late-time event into a testable signal.","feed_headline":"Slow late-universe phase transitions capped at 1% by the CMB","feed_subtitle":"Dark-energy releases as small as 1% show up as CMB anisotropy through patchy vacuum entry.","key_machinery":"The load-bearing object is the \"first-encounter surface\": the corrugated surface at z̄_pt + δz_pt(θ,φ) where each CMB photon first enters the true vacuum. Two identities carry the argument: the comoving-distance matching condition χ_AO = χ_BO (Eq. 8), which cancels the linear term and leaves δz0 quadratic in δz_pt, and the projection of the resulting redshift power spectrum Pδz0(k)—built from the finite-bubble-statistics spectrum Pδt(k)—onto CMB multipoles via spherical Bessel functions (Eq. 13). The whole effect is a late-time analogue of the integrated Sachs-Wolfe effect: stochastic bubble nucleation sources metric perturbations near z̄_pt, and the second-order relation sets how strongly t","core_discovery":"Late-time first-order phase transitions leave a CMB anisotropy even if the dark radiation they release is homogeneous afterwards. Fluctuations in local transition time enter the observed photon redshift only at second order, δz0 ≈ δz_pt² rΩΛ /[(1+z̄_pt)(ΩΛ+Ωm(1+z̄_pt)³)^{3/2}], because the equality of comoving distances before and after the transition cancels the linear term. Squaring this with the finite-bubble-statistics spectrum Pδt(k) gives an induced CMB power spectrum peaked at ℓ<20 with amplitude ∝ r²(β/H⋆)⁻⁴. Equating that to the observed scalar-perturbation amplitude yields r ≲ 10⁻⁵ (β/H⋆)²: slow completed transitions are capped near 1% of dark energy, negative-vacuum endpoints requ","pith_inferences":["An immediate extension is to compute the same first-encounter-surface signal for transitions that reheat into massive or self-interacting particles; the paper notes such cases keep density perturbations alive, so the bounds presented here should be treated as lower limits on the observable signal.","Because the observable is δz_pt², the induced CMB map is not Gaussian: a skewness or bispectrum search at ℓ<20 could separate a phase-transition origin from ordinary adiabatic integrated Sachs-Wolfe anisotropy.","The same corrugated-entry-surface logic should apply at other redshifts—for example in the 1 eV to 10 keV window where other studies bound dark-sector phase transitions—so the quadratic redshift relation offers a template for a full timeline of first-order phase-transition constraints.","If a late-time transition is incomplete, the anisotropy calculation does not apply as written, but photons would still cross a patchwork of false and true vacuum; quantifying that regime is the natural next step and could close the remaining allowed window."],"forward_implications":["A completed slow phase transition (β/H⋆ ≲ 25) in the dark-energy-dominated era can release at most about 1% of the dark energy; anything larger would show up as an excess in the low-multipole CMB temperature spectrum.","Naive Hubble-expansion constraints would allow r as large as ~0.65 for a transition at z̄_pt = 0.3, while the CMB anisotropy bound is much stronger for β/H⋆ ≲ 200.","Negative-vacuum (r > 1) endpoints are not excluded, but only for fast transitions with β/H⋆ ≳ 300; if β/H⋆ ≲ 500, the universe must survive at least 14 Gyr before crunching.","Because the signal is quadratic in the transition-time fluctuation, the bound scales as r ∝ (β/H⋆)², so faster transitions evade the CMB constraint more easily.","The constraints are minimal: converting latent heat into free-streaming radiation homogenizes density perturbations, and the authors note that slow-moving final products would preserve perturbations and yield even stronger bounds."],"supporting_citations":[{"why":"Supplies the nucleation-rate parametrization Γ = Γ0 e^{-S(t)} with exponent β(t-t_f), the starting point for stochastic transition-time fluctuations.","marker":"[79–81]"},{"why":"Gives the average comoving bubble separation d_b, which sets the scale where Pδt(k) changes branch.","marker":"[80, 81]"},{"why":"Defines the reference time t_f as the moment when an e^{-1} fraction of space remains in the false vacuum, fixing the meaning of β/H⋆.","marker":"[82]"},{"why":"Supplies the dimensionless power spectrum of phase-transition-time fluctuations from finite bubble statistics, the input from which the CMB redshift power spectrum is built.","marker":"[43]"},{"why":"Supplies the Planck 2018 best-fit ΛCDM parameters (H0, ΩΛ, Ωm) used in the Hubble model H_pt(z, z_pt, r).","marker":"[83]"},{"why":"Supplies the 1σ CMB temperature-anisotropy error bars used in the χ² upper bound on r.","marker":"[84]"},{"why":"Provides the integrated Sachs-Wolfe formalism used to project induced redshift perturbations onto CMB temperature multipoles.","marker":"[85]"},{"why":"Provides the baryon-acoustic-oscillation h(z) constraint that defines the naive Hubble-expansion bound the CMB bound improves upon.","marker":"[86]"}],"fun_headline_variants":["CMB caps slow dark-energy phase transitions at 1%","Even homogeneous late phase transitions show in CMB","Patchy vacuum entry leaves CMB signature","Dark-energy transitions must be fast or nearly invisible to CMB","CMB anisotropy from patchy vacuum entry constrains phase transitions"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The limit rests on assuming that a patchy phase transition changes only the expansion rate, not the paths photons take, so every photon that starts from the same redshift covers exactly the same spatial distance; if the transition bends photon trajectories at first order, the quadratic redshift relation and the r≲10⁻⁵(β/H⋆)² bound no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["CMB caps slow dark-energy phase transitions at 1%","Even homogeneous late phase transitions show in CMB","Patchy vacuum entry leaves CMB signature","Dark-energy transitions must be fast or nearly invisible to CMB","CMB anisotropy from patchy vacuum entry constrains phase transitions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000883,"raw_usage":{"total_tokens":3677,"prompt_tokens":799,"completion_tokens":2878,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":2799}},"tokens_in":543,"tokens_out":2878,"duration_ms":24254,"temperature":1.0,"reasoning_tokens":2799,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:51:58.146228+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One decisive check is to compute or simulate photon geodesics through a stochastic bubble-nucleation history at z≈0.1–0.3 and test the paper's matching condition χ_AO = χ_BO (Eq. 8) at first order in δz_pt: any linear-in-δz_pt contribution to δz0 invalidates the quadratic formula and the derived bound. Observationally, a low-multipole (ℓ<20) CMB temperature measurement with uncertainties smaller than the predicted D_l^TT,pt for a given (r, β/H⋆, z̄_pt) would confirm or exclude the signal; an absence of the predicted r²(β/H⋆)⁻⁴ scaling would likewise settle the claim.","supporting_citations":[],"review_version":1}