{"id":"b599692d-61f0-4949-965b-599480b60155","arxiv_id":"2412.09266","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"By letting the plasma background temperature and velocity vary with position, the authors remove the singularity that plagued previous fluid Ansatz calculations of bubble wall friction.","lead":"This paper fixes a known mathematical blow-up in calculations of how fast vacuum bubbles expand during a cosmological phase transition. The fix matters because bubble wall speed controls the gravitational wave signal and matter-antimatter asymmetry these transitions could produce.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detonation solutions rest on an enforced source cancellation that the paper's own Fig. 8 shows is violated for exactly the detonation branch.","rationale":"The reader's weakest assumption correctly identifies the enforcement of S1 = S2 = 0 as the load-bearing premise. My stress-test sharpens this to the detonation branch: the paper's own Fig. 8 and Sec. 6.1 show that the integrated sources do not vanish for detonations, which are precisely the solutions claimed to exist only in a narrow parameter window. Since the non-singularity of the friction and the algebraic elimination of light fluctuations both rely on Eq. (12) being sourceless, the detonation results are internally unsupported. The concrete test would settle whether the enforcement is a harmless approximation or a source of spurious solutions. I do not propose rejecting the paper outright because the method may still be salvageable by quantifying the residual-source corrections or restricting claims to deflagrations; however, the central claim as stated in the abstract ('all terms that would give rise to a singularity now vanish') is too strong given the demonstrated nonzero residuals for detonations. Hence the reader's CONDITIONAL verdict is appropriate, and my analysis confirms rather than changes it.","tokens_in":21456,"tokens_out":2373,"duration_ms":26766,"concrete_test":"Re-solve the linearized Boltzmann system for Λ in [625, 635] GeV on the detonation branch without enforcing S1 = S2 = 0: keep the actual on-shell residual sources (Eqs. 28-29) on the right-hand side of Eq. (12), and recompute vw and Lw from Eqs. (50)-(51). If no finite solution exists, or if the light-species fluctuation δvlight diverges as vw approaches cs, the claimed non-singular detonation solutions are an artifact of the enforced cancellation.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that introducing a z-dependent background makes the linearized energy-momentum conservation equations (Eq. 12) sourceless, so the singular undamped combination is no longer sourced and the singularity disappears. This requires S1 = S2 = 0, where S1 and S2 are given in Eqs. (28)-(29). However, Sec. 3 states that the sources are enforced to vanish, not derived from first principles: 'we will enforce these sources to vanish and then study their behavior a posteriori'. The a posteriori check in Sec. 6.1 (Fig. 8) shows the cancellation is incomplete in general. For deflagrations the integrated sources approximately vanish, but for the detonation case the paper states explicitly that 'the sources do not integrate to zero in this case'. Since detonations are claimed only in the narrow window Λ ∈ [625, 635] GeV, the inconsistency is not peripheral: the existence and properties of detonation solutions depend on setting S1 = S2 = 0 even though the on-shell background produces nonzero integrated sources there. If the residual sources are included, the undamped equation is sourced again, the algebraic relation used to eliminate light fluctuations (App. C) loses its justification, and the singular behavior may reappear or the terminal velocity and width may shift substantially. The approximations listed in Sec. 6.1 (neglect of mass dependence in Boltzmann coefficients, specific wall Ansatz) are exactly what break the exact cancellation, so this is an internal inconsistency of the method rather than a mere technical detail.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the well-known singularity in linearized fluid-Ansatz treatments of bubble wall friction when the wall velocity approaches the speed of sound. The authors promote the equilibrium background to a space-dependent Tbg(z), vbg(z) determined by nonlinear energy-momentum conservation, and include the resulting derivative terms in the linearized Boltzmann equations. They argue that the two moments representing total energy-momentum conservation then have vanishing sources, so the undamped combination that caused the singularity is no longer excited; they enforce S1=S2=0, eliminate the light fluctuations via an algebraic relation, and solve the Higgs equation of motion for vw and Lw. For the Standard Model with a low cutoff they find deflagrations for essentially all Λ in [600,800] GeV and non-runaway detonations only for Λ∈[625,635] GeV. They also quantify the out-of-equilibrium contribution to the wall pressure and find it non-negligible, giving corrections of 25–45% in vw relative to the equilibrium-only treatment.","tokens_in":21788,"tokens_out":5608,"duration_ms":59441,"significance":"If the method is valid, it provides a semi-analytic, physically interpretable route to bubble wall velocities without the cs singularity, while including heavy gauge bosons and light fluctuations, and it sharpens the current debate on whether equilibrium friction suffices for wall velocity estimates. The paper is unusually explicit about its approximations and provides an a posteriori check in Fig. 8, which is a strength. The main quantitative predictions—deflagrations for almost all Λ and detonations in a narrow Λ window—are concrete and falsifiable. However, the central singularity-removal claim rests on an enforced source cancellation whose on-shell residual is nonzero precisely for the detonation branch, so the current version does not fully establish that headline result.","major_comments":[{"comment":"The singularity removal is implemented by imposing S1 = S2 = 0 rather than by deriving these equations from the solved system. The authors state in Sec. 3: 'we will enforce these sources to vanish and then study their behavior a posteriori.' The a posteriori check in Sec. 6.1 shows that the cancellation is not exact in general, and for the detonation case 'the sources do not integrate to zero in this case.' This is not a peripheral failure: the detonation solutions in Fig. 7 exist only in the narrow window Λ ∈ [625, 635] GeV, and their computation uses the algebraic relation in App. C that is justified only when the energy-momentum moments are truly sourceless. If the residual sources are retained, the undamped combination is sourced again, and the detonation velocity and width can shift substantially. The manuscript should either include the residual sources in a controlled way and show that the detonation solutions survive, or explicitly limit the singularity-resolution claim to the deflagration/hybrid branches and present the detonation results as an illustration of the enforced-S1=S2=0 system rather than as a prediction.","section":"Sec. 3 and Sec. 6.1, Eqs. (28)-(29), Fig. 8"},{"comment":"The on-shell argument in Appendix B demonstrates a cancellation at linear order in the fluctuations, but the exactness of S1 = S2 = 0 is broken by precisely the two approximations listed in Sec. 6.1: neglect of mass dependence in most Boltzmann moments and the tanh wall Ansatz of Eq. (35). The paper acknowledges this, saying 'Both of these effects will lead to the fact that the sources in (26) are not exactly zero.' Given this limitation, the abstract's statement that 'all terms that would give rise to a singularity now vanish' is stronger than what is established. What is established is that, under the enforced source condition and with the stated approximations, the singular source is absent at the order retained. Please rephrase the central claim so that this conditionality is visible, and state explicitly which of the two approximations is responsible for the residual visible in Fig. 8 for strong deflagrations.","section":"Appendix B and Sec. 6.1"},{"comment":"The quantitative conclusion that out-of-equilibrium effects change vw by 25–45% is obtained by comparing the full solutions with 'equilibrium only' solutions obtained by dropping the ffl, flight, gfl, and glight terms in Eqs. (50)-(51). In Fig. 6 the wall width is fixed at LwT+ = 15 for all vw, whereas the actual full solutions in Fig. 7 have Lw varying between roughly 5 and 45. A sensitivity check in which the equilibrium-only case is solved with a self-consistently determined Lw would make the comparison more robust and would clarify how much of the claimed correction is an artifact of the fixed width used in the pressure curves.","section":"Sec. 6, Figs. 6-7"}],"minor_comments":[{"comment":"The quantities δτf, δτb, and δτlight are used in Eq. (47) before being defined; please define δτ in terms of the fluctuation variables introduced in Eq. (4), including the sign convention.","section":"Sec. 5.2, Eq. (47)"},{"comment":"The caption says 'Source momentum' and labels curves as '1st old', '1st new', '1st tot', and '2nd tot', but the legend in the figure also uses line styles; please make explicit which line corresponds to S2 and whether any rescaling has been applied to make the curves comparable.","section":"Fig. 8"},{"comment":"The statement that a term from the Liouville operator cancels against a force term, ensuring that the remaining force term vanishes for zero wall velocity, is given without derivation; a short equation or reference would improve reproducibility.","section":"Sec. 2, footnote 2"},{"comment":"The phrase 'non-singular solutions to the Boltzmann equation' in the title and abstract could be read as applying to the full nonlinear equation, whereas the paper solves the linearized system with a truncated fluid Ansatz; a qualifier such as 'linearized' would avoid overstatement.","section":"Abstract and Sec. 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent about its main weakness, which is commendable, but that weakness is nonetheless central: the detonation branch is a headline quantitative result and it relies on an enforced source cancellation that the paper's own Fig. 8 shows to be violated for detonations. I would ask the authors to either control the residual sources or substantially qualify the detonation claims before acceptance. The relation to WallGo [40] is handled candidly; a numerical comparison of wall velocities with that package would strengthen the revised version if the authors continue to claim differences."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper directly attacks a recognized problem and gets partway there. The deflagration resolution of the speed-of-sound singularity looks credible; the detonation results do not.\n\nWhat's actually new: the z-dependent background in the fluid Ansatz was suggested earlier and implemented by Laurent & Cline with Chebyshev polynomials, but this paper adds the light-species fluctuations and keeps the fluid Ansatz, which preserves the physical interpretation in terms of local temperature, velocity, and chemical potential. That is a genuine step forward. The formalism is clearly presented, the benchmark (SM with a low cutoff) is well chosen, and the quantitative comparison with equilibrium-only friction is useful: out-of-equilibrium effects shift vw by 25–45%, which speaks directly to recent claims that equilibrium suffices.\n\nThe soft spot is the source cancellation. The two linear combinations of the Boltzmann equations that represent energy-momentum conservation, Eqs. (28)–(29), are enforced to vanish rather than derived. Appendix B gives an on-shell argument that they should vanish at linear order, but the numerical check in Fig. 8 shows the cancellation is only approximate: for deflagrations the integrated sources roughly vanish, while for detonations they explicitly do not. Since the detonation solutions, confined to Λ ∈ [625, 635] GeV, are computed under S1 = S2 = 0, and the paper itself says the integrated sources don't vanish there, the algebraic relation that eliminates the light fluctuations (App. C) loses its justification on exactly that branch. The terminal velocities and widths for those detonations could shift substantially if the residual sources are included, and the singularity could reappear. This is not a nitpick: the detonation window is one of the paper's concrete results.\n\nWhat holds up: the deflagration side. The weaker condition that actually removes the singularity is the integrated version (Eq. 13), and for deflagrations that integrated condition is approximately satisfied. So the main resolution of the speed-of-sound problem for subsonic walls looks plausible, and the relative importance of light species is sensibly quantified. The citation pattern is honest—it builds on the authors' earlier singularity paper, engages with Laurent & Cline, and compares with the recent WallGo work. No code or data are shipped, which is a minor annoyance for a computational paper.\n\nWho it's for: anyone computing bubble wall velocities in first-order phase transitions. It deserves a serious referee—the topic is important and the method is a step forward—but the referee should push for the detonation branch to be either fixed or dropped, and for the abstract to match the caveats in Sec. 6.1. I'd send it to review with the expectation of major revision.","headline":"A promising fix for the speed-of-sound singularity in the fluid Ansatz, but the detonation branch rests on an enforced source cancellation that the paper's own check shows to fail.","tokens_in":22272,"tokens_out":5240,"would_cite":true,"duration_ms":50255,"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":"The friction on a cosmological bubble wall no longer diverges at the speed of sound when the plasma background is allowed to vary in space, and in a Standard Model with a low cutoff this yields deflagration solutions for almost all cutoff…","keywords":["cosmological phase transitions","bubble wall velocity","Boltzmann equation","fluid Ansatz","speed of sound singularity","plasma friction","electroweak phase transition","gravitational waves"],"falsifier":"Solve the Boltzmann equation for the same low-cutoff Standard Model on a momentum-space lattice without imposing the fluid Ansatz, and check whether the friction diverges as the wall speed approaches the speed of sound; a divergence would show the singularity is a real feature rather than an artifact of the constant-background linearization.","tokens_in":21281,"feed_emoji":"🌌","tokens_out":12928,"duration_ms":103491,"temperature":0.7,"pith_summary":"The paper claims that the well-known singularity in bubble-wall friction, which appears exactly when the wall velocity equals the speed of sound, is an artifact of assuming a constant plasma background. It shows that allowing the background temperature and fluid velocity to vary with distance to the wall modifies the linearized Boltzmann equation so that all singular source terms vanish. The authors apply this to a Standard Model with a low cutoff, solving for the terminal wall velocity and width. They find non-singular friction, with deflagration solutions for almost every cutoff scale between 600 and 800 GeV and non-runaway detonations only for cutoffs around 625–635 GeV. If correct, this removes a key obstruction to predicting bubble wall velocities and the gravitational waves and baryon asymmetries that depend on them.","feed_headline":"Non-constant plasma background cures bubble-wall friction singularity","feed_subtitle":"Spatially varying background removes speed-of-sound divergence, yielding deflagrations for almost all cutoffs.","key_machinery":"The load-bearing object is the fluid Ansatz with a space-dependent background: out-of-equilibrium distribution functions are expanded around local equilibrium using $T_{\\rm bg}(z)$ and $v_{\\rm bg}(z)$ obtained from the non-linear conservation equations (26). The key identity is the on-shell cancellation of the sources $S_1$ and $S_2$ of the linearized energy-momentum conservation equations: the new terms proportional to $\\partial_z T$ and $\\partial_z v$ cancel the mass-dependent source terms, turning the undamped equations into algebraic relations between light and heavy species fluctuations. This cancellation removes the $1/(1-3v_w^2)$ denominator that caused the singularity in the friction.","core_discovery":"The central discovery is that the speed-of-sound singularity in bubble-wall friction is a consequence of expanding the fluid Ansatz around a constant background. The paper defines a z-dependent background — temperature $T_{\\rm bg}(z)$ and velocity $v_{\\rm bg}(z)$ — by solving the non-linear energy-momentum conservation equations (26) across the wall, including the scalar field contribution. Linearizing the Boltzmann equation around this background introduces new source terms proportional to $\\partial_z T$ and $\\partial_z v$. The authors show that when the Higgs equation of motion is satisfied, these new terms cancel the original mass-dependent sources in the two linear combinations of Boltzmann equations that encode energy-momentum conservation; these undamped equations then become algebraic constraints linking light and heavy fluctuations, and the singular factor $1/(1-3v_w^2)$ no longer appears. Enforcing this cancellation, the paper obtains finite friction across the speed of sound, with a discontinuity only at the Jouguet velocity, and uses the resulting pressures to compute terminal wall velocities in the low-cutoff Standard Model benchmark.","pith_inferences":["The same background-absorption strategy could be applied to other linearized transport problems at bubble walls, such as electroweak baryogenesis, where constant-background expansions produce analogous singular denominators.","The incomplete source cancellation for strong transitions and detonations (shown in Fig. 8) suggests that including next-order corrections could shift wall velocities in those regimes; a direct computation of the residual sources without enforcing $S_1 = S_2 = 0$ would quantify this shift.","A fully variational wall profile, rather than the two-parameter tanh shape, would test whether the source cancellation is robust to the profile choice.","If the narrow detonation window found here is generic across models, gravitational wave predictions based on non-runaway detonations may need to focus on deflagration and hybrid contributions; scanning a broader model space would settle this."],"forward_implications":["The friction on the bubble wall is finite and continuous across the speed of sound; only a discontinuity at the Jouguet velocity remains, where the hydrodynamic boundary conditions change between hybrids and detonations.","Light-species fluctuations now represent only genuine departures from local equilibrium, and their contribution to the wall pressure is small but non-negligible, with opposite sign to the heavy-species contribution.","Out-of-equilibrium corrections shift the terminal wall velocity by roughly 25–45% and the wall width by 10–25% relative to equilibrium-only calculations, so equilibrium backreaction alone does not give accurate predictions.","For the low-cutoff Standard Model benchmark, deflagration solutions exist for almost every cutoff scale in the range 600–800 GeV, while non-runaway detonations occur only for $\\Lambda \\in [625,\\,635]$ GeV and always coexist with a deflagration.","For this benchmark, the deflagration is most likely the stable solution, so detonations would not be realized even where they exist mathematically."],"supporting_citations":[{"why":"Introduced the fluid Ansatz for the out-of-equilibrium distribution, the framework that this paper adapts to a varying background.","marker":"[17]"},{"why":"Established the linearized Boltzmann setup for bubble wall velocities, including the kinetic matrices and source terms used here.","marker":"[18]"},{"why":"Identified the speed-of-sound singularity in the friction and traced it to the constant-background linearization, the problem this paper resolves.","marker":"[20]"},{"why":"First implemented a space-dependent background for wall-velocity calculations, though with Chebyshev expansions rather than the fluid Ansatz; it provides the key comparison point.","marker":"[23]"},{"why":"Introduced the hydrodynamic obstruction (equilibrium backreaction) term that this paper includes alongside non-equilibrium contributions.","marker":"[24]"},{"why":"Provided results on electroweak bubble wall expansion in a Standard Model-like plasma, against which the deflagration and detonation solutions are compared.","marker":"[27]"},{"why":"A recent computation of bubble wall velocities used for comparison; the paper discusses differences in implementation.","marker":"[40]"}],"fun_headline_variants":["Spatially varying fluid kills bubble-wall speed-of-sound singularity","Boltzmann fix: z-dependent background removes friction divergence","Bubble-wall velocity singularity vanishes with non-constant plasma","Deflagrations for nearly all cutoffs after Boltzmann background fix","Fluid Ansatz with spatial variation ends singular bubble friction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's resolution depends on assuming that the two source terms in the linearized energy-momentum conservation equations can be set to zero once the background is chosen; the paper itself shows this cancellation is incomplete for detonations and strong transitions, so if the enforcement fails the singularity could return.","fun_headline_variants_meta":{"raw":{"variants":["Spatially varying fluid kills bubble-wall speed-of-sound singularity","Boltzmann fix: z-dependent background removes friction divergence","Bubble-wall velocity singularity vanishes with non-constant plasma","Deflagrations for nearly all cutoffs after Boltzmann background fix","Fluid Ansatz with spatial variation ends singular bubble friction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000144,"raw_usage":{"total_tokens":1185,"prompt_tokens":963,"completion_tokens":222,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":139}},"tokens_in":579,"tokens_out":222,"duration_ms":3280,"temperature":1.0,"reasoning_tokens":139,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:05:09.124745+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the Boltzmann equation for the same low-cutoff Standard Model on a momentum-space lattice without imposing the fluid Ansatz, and check whether the friction diverges as the wall speed approaches the speed of sound; a divergence would show the singularity is a real feature rather than an artifact of the constant-background linearization.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identified the speed-of-sound singularity in the friction and traced it to the constant-background linearization, the problem this paper resolves."},{"cited_title":"Laurent and J","cited_arxiv_id":null,"evidence_quote":"First implemented a space-dependent background for wall-velocity calculations, though with Chebyshev expansions rather than the fluid Ansatz; it provides the key comparison point."}],"review_version":1}