{"id":"505d19a3-3331-440d-8b0e-9daa1f99f0f2","arxiv_id":"2411.15094","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Bubble-wall terminal velocity depends on plasma thermalization; non-equilibrium and free-streaming regimes give slower or different walls, and stationary solutions can be bypassed by runaways.","lead":"This paper studies how fast vacuum bubbles grow during a cosmological first-order phase transition, depending on whether the surrounding plasma stays in thermal equilibrium. It finds that non-equilibrium effects slow bubble walls in some regimes and can make bubbles run away even when a steady state exists.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MPC collision scheme cannot reproduce LTE even at τ→0, so the finite-τ wall-velocity interpolation in Fig. 6 and the existence claims in Fig. 7 are not settled.","rationale":"Good-faith reading: the authors are careful and candid about the limitations of the N-body collision algorithm, and they do not overclaim; the qualitative message that LTE may be insufficient and that stationary states are not guaranteed is supported by independent evidence: (i) the LTE hydrodynamical simulations reproduce the matching-method velocities and show runaways even when stationary states exist, (ii) the fully ballistic analytic pressure formula agrees with the large-τ N-body results, and (iii) the analytic intermediate scenario (LTE away from wall, ballistic inside wall) is consistent with entropy-production intuition. My concern is therefore not with the overall direction of the paper but with the specific intermediate-thermalization interpolation and the ballistic existence claims, which are the least secure and are exactly where the reader's weakest assumption is located. The paper itself flags the relevant failure modes (non-locality, once-per-timestep collisions, no LTE at τ = 0, poor shock resolution), so the concern is internal rather than a matter of disagreeing with consensus. A homogeneous calibration and a resolution or variant study of MPC collisions would distinguish physical nonequilibrium effects from algorithmic artifacts. I do not see the concern as fatal: the main qualitative conclusions survive through the analytic and LTE simulation channels, but the quantitative finite-τ curves and existence statements should be treated as conditional pending such a test.","tokens_in":20172,"tokens_out":4424,"duration_ms":41682,"concrete_test":"Rerun the small-τ Fig. 6 points (e.g., Δm/T = 3.4, α ≈ 0.5) with a modified MPC collision step that allows up to 8 particles per cell and/or halves Rcell, keeping N, Δt, and τ fixed. If the terminal velocity shifts by more than ~10% toward the analytic ΔP = ΔV0 curve, or if a homogeneous τ = 0 run fails to relax to the Maxwell–Boltzmann distribution, then the finite-τ interpolation is a collision-scheme artifact rather than evidence about plasma thermalization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix B states that the MPC algorithm 'violates causality due to the finite size of the collision size', that results for vw ≳ 0.7 'become less reliable', and that 'even if τ = 0, the N-body simulation algorithm does not result in LTE in front of the bubble, as collisions in the algorithm happen once per timestep.' These are exactly the conditions under which the paper's central interpolation claim is made: Sec. 4.2 and Fig. 6 use small-τ N-body runs to show that the system 'asymptotically approaches the LTE limit' and that thermalization time affects the terminal velocity, and the τ→0 N-body points lie below the analytic ΔP = ΔV0 limit. The disagreement is attributed to physical nonequilibrium effects, but a collision scheme whose τ→0 limit is not LTE cannot separate physical non-thermal behavior from algorithmic artifacts. The existence statements in Sec. 4.3 / Fig. 7, where σ is found to affect whether a stationary solution is reached in ballistic simulations, rest on the same MPC implementation. Thus the quantitative interpolation and the 'stationary state may not be reached' claim in the ballistic regime are load-bearing on an unvalidated collision model; the supporting LTE hydrodynamics and fully ballistic analytic limits do not cover this intermediate regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the terminal velocity of vacuum bubble walls in a first-order cosmological phase transition, comparing three regimes: local thermal equilibrium (LTE) everywhere, LTE away from the wall with ballistic motion inside the wall, and a fully ballistic fluid. The authors derive terminal velocities from entropy conservation in the LTE case and from a pressure-balance condition in the ballistic cases, and they support the analysis with hydrodynamic lattice simulations and N-body simulations using multi-particle collision dynamics (MPC). The main claims are that the wall velocity depends on the degree of thermalization, that the fully ballistic scenario gives significantly different velocities, and that in both hydrodynamic and ballistic simulations a stationary solution, even when it exists, may not actually be reached.","tokens_in":20339,"tokens_out":8714,"duration_ms":88955,"significance":"If correct, the paper's central message is important: LTE-based predictions of bubble-wall velocities can be unreliable in non-equilibrium regimes, and time-dependent simulations, rather than stationary-state analyses alone, may be needed to determine the final fate of a bubble. The strengths of the paper are the internally consistent analytic derivations (entropy conservation in Sec. 3.1 and pressure balance in Secs. 3.2-3.3), the good quantitative agreement between hydrodynamic simulations and the stationary-state method in the LTE regime (Fig. 2 and Fig. 4, top panel), and the reproduction of the fully ballistic analytic limit in the large-mean-free-path N-body simulations (Fig. 6, tau to infinity limit). The significance is tempered, however, by the reliance of the intermediate-thermalization results on an MPC collision scheme whose limitations are acknowledged in Appendix B.","major_comments":[{"comment":"The small-tau N-body points in Fig. 6 are used to argue that the system asymptotically approaches the LTE limit and that the thermalization time changes the stationary wall velocity. Appendix B states, however, that 'even if tau = 0, the N-body simulation algorithm does not result in LTE in front of the bubble, as collisions in the algorithm happen once per timestep.' Because the tau to 0 limit of the algorithm is not LTE, the observed offset between the N-body points and the analytic Delta-P = Delta-V0 limit cannot be unambiguously attributed to physical nonequilibrium effects, as is done in Sec. 4.2. The authors should either demonstrate convergence of the MPC scheme to a known LTE result in a controlled benchmark (for example, a planar wall with decreasing Delta-t and R_cell) or soften the quantitative interpolation claim in Fig. 6 and the associated conclusion that thermalization time changes the stationary wall velocity.","section":"Sec. 4.2 / Fig. 6 and Appendix B"},{"comment":"The conclusion that, in the ballistic regime, the existence of a stationary solution depends on sigma and on whether thermalization is fast enough rests on the same MPC implementation. Appendix B states that the algorithm 'violates causality due to the finite size of the collision size' and that results for v_w >~ 0.7 'become less reliable', and the right panel of Fig. 7 shows wall velocities extending close to v_w = 1. The authors should quantify how much of the parameter region in Fig. 7 lies in the unreliable regime and check whether the 'no stationary solution' region persists when the causality violation is removed or when a different collision scheme is used. As it stands, the abstract's statement that both simulations show a stationary solution may not be reached is stronger than the ballistic evidence supports.","section":"Sec. 4.3 / Fig. 7 and Appendix B"},{"comment":"The paper does not calibrate the MPC collision probability to a physical transport coefficient, so the mapping between the parameter tau and the physical mean free path or thermalization time is not established. The collision probability contains a factor exp(Delta-t/tau), a density-dependent factor, and a cap of two particles per cell, but no comparison with known transport coefficients (for example, viscosity or mean free path in a Maxwell-Boltzmann gas) is provided. Without this calibration, the horizontal axis of Fig. 6 is only a simulation parameter, and the quantitative statement that thermalization time has an important effect on the stationary wall velocity is not yet tied to a physically meaningful scale.","section":"Appendix B and Fig. 6"}],"minor_comments":[{"comment":"The sentence 'will not make this distinction in terminology' is missing a subject; it should read 'we will not make this distinction'.","section":"Footnote 1"},{"comment":"The phrase 'violates causality due to the finite size of the collision size' should presumably read 'collision cell size'; please correct this and clarify which characteristic scale is meant.","section":"Appendix B"},{"comment":"The text says that the effect of the statistics is 'quantified explicitly in Sec. 3.1', but the quantitative comparison of Maxwell-Boltzmann, Bose-Einstein, and Fermi-Dirac distributions appears in Sec. 4.2 (Fig. 5).","section":"Sec. 2.1"},{"comment":"The notation T_j and R_j in Eq. (3.5) is not defined; the text defines a reflection coefficient R and a transmission coefficient T = 1 - R. Please clarify whether the subscript j labels the side of the wall and how T_j and R_j relate to R.","section":"Sec. 3.2, Eq. (3.5)"},{"comment":"The horizontal axis label 'tau / Delta t infinity' is ambiguous; please indicate explicitly that the right boundary corresponds to the free-streaming limit tau to infinity.","section":"Fig. 6"},{"comment":"The definition 'Delta V0 = V(0,T) - V(phi0,T)' appears inconsistent with the earlier definition of Delta V0 as the bare potential difference in Eq. (2.14) and with the later statement in Sec. 3.2 that Delta V0 omits thermal corrections. Please clarify whether V here denotes the bare potential V0 or the thermal potential V.","section":"Sec. 2.3, Eq. (2.17)"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about the limitations of the MPC collision scheme, and my recommendation is driven by the gap between that honesty and the strength of the abstract and conclusions. The analytic parts and the hydrodynamic simulations are solid and valuable; the request for validation or softening of the intermediate-thermalization claims is, in my view, a normal part of revision rather than a reason to reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper deserves a serious read. It compares bubble-wall dynamics across three thermalization regimes—LTE everywhere, LTE outside the wall with ballistic crossing, and fully ballistic—and shows that the thermalization time changes the terminal wall velocity and can determine whether a stationary state is reached at all. That last point is the substantive one: finding a stationary solution is not the same as the bubble reaching it.\n\nWhat's genuinely new is the unified three-regime map and the finite-thermalization-time scan in Fig. 6, plus the surface-tension effect on reaching a stationary state in the ballistic runs (Fig. 7). The authors are transparent that parts of the non-reachability story already appeared in their earlier work; the novel contribution is the systematic comparison and the explicit time-dependence. The entropy-conservation derivation is clean, and the LTE hydro simulations match the stationary-state matching wherever a steady state is reached. The fully ballistic limit reproduces the analytic formula, which checks the pipeline.\n\nThe main soft spot is the intermediate regime. The MPC collision scheme used in the N-body runs cannot reproduce LTE even as τ→0, violates causality for vw ≳ 0.7, and has poor shock resolution. The authors say this in Appendix B, which I respect, but it means the quantitative interpolation in Fig. 6 and the existence statements based on self-interacting runs are the least secure parts. The τ→0 points lying below the analytic LTE limit could be physical nonequilibrium physics, but with a collision scheme whose τ→0 limit isn't LTE, the attribution is not established. I'd phrase those results as indicative, not quantitative. Missing code and error bars make this harder to check. That said, the qualitative picture is robust: the three regimes give clearly different terminal velocities, and the LTE and fully ballistic limits are well anchored.\n\nI would send this to a serious referee. The topic matters for gravitational-wave predictions and electroweak baryogenesis, and the paper is a clear step toward understanding when equilibrium approximations break down. I'd recommend conditional acceptance: improve or validate the collision model for the intermediate regime, release the code and data, and mark Fig. 6 as proof-of-principle until then. I would cite it for the three-regime comparison even with those caveats.","headline":"A useful three-regime comparison of bubble-wall thermalization; the qualitative message holds, but the intermediate-regime numbers rest on an admittedly imperfect collision scheme.","tokens_in":20960,"tokens_out":4227,"would_cite":true,"duration_ms":40443,"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 terminal speed of an expanding vacuum bubble is set by the degree of thermalization of the surrounding plasma, and stationary solutions that exist in steady-state theory may never be reached in practice.","keywords":["bubble wall velocity","cosmological first-order phase transitions","local thermal equilibrium","ballistic pressure","N-body simulations","entropy conservation","gravitational waves","electroweak baryogenesis"],"falsifier":"Run the same ballistic setup with a smaller collision-cell size and with collision cells allowed to hold more than two particles; if the inferred terminal velocity changes by more than the quoted error, or if the $\\tau \\to 0$ limit does not approach the entropy-conservation terminal velocity as the cell shrinks, then the intermediate-thermalization results are a numerical artefact of the collision scheme.","tokens_in":1607,"feed_emoji":"🫧","tokens_out":1841,"duration_ms":90405,"temperature":0.7,"pith_summary":"This paper asks what limits the growth speed of vacuum bubbles in cosmological first-order phase transitions when the surrounding plasma does not fully thermalize. It argues that the bubble wall's terminal velocity depends on where the particle mean free path lies relative to the wall thickness: entropy conservation fixes the speed when local thermal equilibrium holds everywhere, while matching ballistic pressure against the vacuum-energy difference fixes a different speed when particles cross the wall without thermalizing. The paper finds that the fully ballistic limit gives significantly different wall velocities from the equilibrium predictions, and that intermediate thermalization times can either speed up or slow down the wall. It also finds, in both hydrodynamic and particle-based simulations, that a stationary terminal velocity that exists in steady-state theory may never be reached because the accelerating wall outruns the formation of its heated plasma shell. If this is right, predictions for gravitational-wave signals and electroweak baryogenesis that assume a single LTE terminal velocity carry a systematic uncertainty.","feed_headline":"Thermalization controls vacuum-bubble wall velocity","feed_subtitle":"How much the plasma thermalizes sets the wall speed; assuming equilibrium can miss runaway bubbles.","key_machinery":"The organizing quantity is the particle mean free path compared with the bubble-wall thickness, which selects one of three regimes: LTE everywhere, LTE outside the wall with ballistic transport inside it, or a fully ballistic fluid. The terminal velocity is then fixed by one of two closure conditions: entropy-flux conservation, $s_-\\gamma_- v_- = s_+\\gamma_+ v_+$, in the LTE case, or pressure balance $\\Delta P = \\Delta V_0$, with $\\Delta P$ computed from particle momentum transfer through the wall, in the ballistic cases. For a fully free-streaming Maxwell–Boltzmann fluid, $\\Delta P$ takes the closed form of Eq. (3.8), whose relativistic limit reproduces the standard result. Lattice hydrodynamic simulations implement the LTE closure, while N-body simulations with collision dynamics implement the ballistic and intermediate cases and test whether the stationary state is actually approached.","core_discovery":"On the paper's own terms, the bubble-wall terminal velocity is determined by how well the plasma thermalizes around the wall. In local thermal equilibrium everywhere, the additional condition fixing the wall velocity is entropy conservation across the wall. When equilibrium holds only away from the wall and particles move ballistically through it, the terminal velocity is instead fixed by matching the ballistic pressure with the potential-energy difference driving the bubble; this gives slightly slower walls, because entropy production inside the wall acts as extra friction. When the fluid is fully ballistic, the terminal velocity is set by a closed-form pressure balance and can differ substantially from the other two cases. Beyond these stationary predictions, the paper shows that even when such a solution exists, the actual dynamics matters: bubbles accelerating strongly can run away before the heated plasma shell characteristic of a steady deflagration or hybrid forms, so the terminal state is not automatically realized.","pith_inferences":["If these results carry over to realistic models, gravitational-wave spectra computed from a single LTE wall velocity should be treated as a band rather than a line, and electroweak-baryogenesis calculations that assume steady subsonic walls may miss a population of runaways.","The finding that existing stationary states can be unreachable suggests the useful observable is the distribution of wall velocities at bubble collision, not a single terminal velocity from a steady-state analysis.","Extending the three-limit comparison to models with nontrivial quantum reflection coefficients and soft-particle emission would test whether the paper's step-function reflection and neglected emission bracket the error or underestimate it.","Because the paper's collision scheme violates causality for wall speeds above about $0.7c$ by its own account, the fast-wall branch of the thermalization scan is the least trustworthy; a causal kinetic or particle-in-cell treatment would provide a sharper check of whether the short-mean-free-path limit converges to the LTE result."],"forward_implications":["In the LTE regime, the hydrodynamic simulations and the entropy-conservation matching agree on the wall velocity whenever a steady state is actually reached, so the mapped velocity contours can be used for that regime.","The intermediate scenario, with ballistic transport inside the wall and LTE outside it, predicts slightly lower terminal velocities than full LTE, meaning non-equilibrium effects inside the wall provide extra friction.","The fully ballistic limit predicts substantially different wall velocities, and for strong transitions the wall is typically faster when particles barely interact.","Whether a stationary state is reached can depend on wall tension and on thermalization time, even though the terminal velocity itself does not depend on tension in the large-radius limit; quickly accelerating walls can run away before the heated plasma shell forms.","Because thermalization time affects both the value of the terminal velocity and the existence of a terminal state, time-dependent non-equilibrium simulation is needed to predict the final fate of the bubble."],"supporting_citations":[{"why":"Supplies the N-body simulation method and the ballistic pressure calculation used for bubble-wall motion in a fluid.","marker":"[74]"},{"why":"Provides the ballistic-pressure formula and dynamics of false vacuum bubbles used for the fully non-interacting limit.","marker":"[75]"},{"why":"Gives the relativistic-wall pressure limit that the fully ballistic pressure formula is required to reproduce.","marker":"[87]"},{"why":"Establishes the entropy-conservation condition for bubble wall velocity in local equilibrium, used in the LTE scenario.","marker":"[86]"},{"why":"Reviews phase-transition dynamics and supports the entropy-current conservation used in the LTE matching condition.","marker":"[85]"},{"why":"Provides the model-independent energy budget and matching equations used to construct stationary fluid profiles.","marker":"[42]"},{"why":"Gives the hydrodynamical constraints and numerical setup for the lattice simulations.","marker":"[88]"},{"why":"Reports the earlier observation that not all stationary states are reached, which this paper confirms and extends.","marker":"[89]"}],"fun_headline_variants":["Bubble wall speed hinges on plasma thermalization","Thermalization decides bubble-wall terminal velocity","Non-equilibrium slows vacuum bubble walls","Runaway vacuum bubbles defy thermal equilibrium forecasts"],"cache_read_input_tokens":23040,"weakest_assumption_plain":"The load-bearing premise is that the particle-collision algorithm used in the N-body simulations, with at most two particles per collision cell, density- and time-dependent collision probabilities, and finite-size cells, represents physical thermalization accurately enough to determine both the terminal velocity and whether a stationary state is reached; the paper itself reports that this scheme violates causality for wall speeds above about $0.7c$, making those fast-wall results less reliable.","fun_headline_variants_meta":{"raw":{"variants":["Bubble wall speed hinges on plasma thermalization","Thermalization decides bubble-wall terminal velocity","Non-equilibrium slows vacuum bubble walls","Runaway vacuum bubbles defy thermal equilibrium forecasts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000395,"raw_usage":{"total_tokens":2044,"prompt_tokens":889,"completion_tokens":1155,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":1100}},"tokens_in":505,"tokens_out":1155,"duration_ms":9637,"temperature":1.0,"reasoning_tokens":1100,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:33:13.867341+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same ballistic setup with a smaller collision-cell size and with collision cells allowed to hold more than two particles; if the inferred terminal velocity changes by more than the quoted error, or if the $\\tau \\to 0$ limit does not approach the entropy-conservation terminal velocity as the cell shrinks, then the intermediate-thermalization results are a numerical artefact of the collision scheme.","supporting_citations":[],"review_version":1}