{"id":"b781378a-4105-4aa1-87d9-d6a028b09812","arxiv_id":"2505.19584","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using the full one-loop effective potential with the LTE approximation, the authors find that for xSM deflagration, bag-model gravitational wave peak predictions can differ by up to 48% in frequency and 90% in amplitude from full-potential integrated-K results, and that deflagration dominates the…","lead":"The authors compute bubble wall velocities and gravitational wave spectra using the full finite-temperature effective potential under the local thermal equilibrium approximation, instead of the usual bag equation of state. In a scan of the singlet-extended Standard Model, they find deflagration dominates, and peak gravitational wave predictions can differ by up to 48% in frequency and 90% in amplitude from bag-model estimates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fixed-VEV approximation in the wall matching conditions (Sec. 3) is untested and internally inconsistent with the code's temperature-dependent speed of sound, so the reported 48%/90% discrepancies could shift if VEV variation is included.","rationale":"The reader's verdict of CONDITIONAL is appropriate, and the weakest assumption identified is indeed the fixed-VEV approximation. This is load-bearing because the headline 48%/90% discrepancies are dominated by the LTE-derived wall velocity, which enters through the matching conditions that rely on the fixed VEVs. The internal inconsistency with the temperature-dependent speed of sound in the profile integration makes the approximation demonstrably present in the code, not just a hypothetical concern. I considered other potential issues, such as the LTE approximation being an upper bound on the wall velocity or the arbitrary xi_w>0.99 cutoff; these are acknowledged limitations, whereas the fixed-VEV choice is unquantified and testable. The WallGo cross-check at a single benchmark outside the scanned range does not directly test this approximation in the parameter space where the claims are made. The proposed concrete test would determine whether the fixed-VEV choice changes the results by more than a few percent; if it does, the quantitative claims would need revision, but until then the conditional acceptance stands.","tokens_in":20773,"tokens_out":12868,"duration_ms":114066,"concrete_test":"Select 3-5 benchmark points in the scanned parameter space (e.g., ms=200, lambda_hs=1.5; ms=300, lambda_hs=2.0; ms=400, lambda_hs=1.0, with lambda_s=0.2). Recompute both deflagration and detonation solutions using self-consistent VEVs: at each matching iteration, set X_+ and X_- to the phase minima at the current T_+ and T_- (the same valAt(T) routine used for cs_squ) and solve Eqs. (3.2)-(3.4) again. Compare the resulting xi_w, K, and peak Omega_sw h^2 with the fixed-VEV results, and directly report the fractional change of the broken-phase VEV between Tn and the solved T_- for each point. If xi_w or K shifts by more than about 5%, the fixed-VEV approximation is not valid in this region and the reported discrepancies are not stable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claims (up to 48% in frequency, 90% in amplitude) hinge on the LTE-derived wall velocity and integrated kinetic energy fraction, which come from the wall matching conditions (3.1)-(3.4). In those conditions, p_+ and p_- and the energy densities are evaluated at VEVs fixed to their reference-temperature values, justified only by the assertion that 'the VEV typically exhibits weak temperature dependence' (Sec. 3). No test of this approximation is given. Meanwhile, the speed of sound used in the fluid profile integration is computed with valAt(T), i.e., with temperature-dependent VEVs (Appendix 6.1, cs_squ). Thus the wall boundary conditions and the subsequent profile evolution effectively use different equations of state. In the xSM, the false vacuum has a nonzero singlet VEV that can vary substantially with temperature, and in deflagrations T_+ differs from the nucleation temperature due to shock heating. Fixing the VEVs at Tn can therefore change the pressure difference across the wall, directly shifting the solved xi_w and K. The only external cross-check (WallGo) uses a benchmark with ms=120 GeV, outside the scanned range (ms>=150 GeV), so it does not validate this approximation in the region where the headline discrepancies are claimed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a hydrodynamic framework in which the bubble wall velocity is computed from the full one-loop finite-temperature effective potential under the local thermal equilibrium (LTE) approximation, and the kinetic energy fraction K is obtained by direct integration of the fluid profile. The framework is validated against the bag model on an artificial benchmark and against WallGo on one xSM parameter point. The authors then scan the xSM parameter space, classify the allowed hydrodynamic modes, and compare gravitational wave spectra obtained with different equations of state, different methods for computing K, and different choices of wall velocity. Their main quantitative findings are that deflagration is the most prevalent mode in the scanned region and that, in the deflagration regime, spectra from the full effective potential with LTE-derived wall velocity and integrated K differ from bag-model spectra with fitted K by up to 48% in peak frequency and 90% in peak amplitude.","tokens_in":21034,"tokens_out":4537,"duration_ms":44950,"significance":"If the quantitative claims hold, the paper provides a useful step toward reducing the systematic uncertainty in gravitational wave forecasts from cosmological phase transitions: it replaces input wall velocities with an LTE-derived value and computes K from the actual fluid profile rather than from bag-model fitting formulas. The WallGo cross-check at one benchmark agrees to about 1.2%, and the internal bag-model consistency test in Fig. 1 supports the numerical implementation. The main value of the paper is therefore not a new principle but a concrete, model-specific quantification of how much the predicted spectra depend on the treatment of the equation of state, K, and wall velocity. However, the load-bearing fixed-VEV approximation and the arbitrary runaway cutoff need to be tested before the reported 48% and 90% numbers can be regarded as robust.","major_comments":[{"comment":"The fixed-VEV approximation is load-bearing and is not quantified. The matching conditions evaluate p± and ρ± at VEVs frozen at a reference temperature, with the sole justification that 'the VEV typically exhibits weak temperature dependence.' Meanwhile, in Appendix 6.1 the speed of sound is computed using valAt(T), i.e., with temperature-dependent VEVs, so the matching conditions and the subsequent profile integration effectively use two different equations of state. In the xSM the false vacuum has a nonzero singlet VEV that can vary with temperature, and in a deflagration T+ differs from Tn because of shock heating. A VEV shift changes the pressure difference across the wall and therefore directly shifts the solved ξw and K, on which the headline 48%/90% discrepancies rest. Please provide a numerical test, for example by evaluating the matching conditions at Tn, T−, and T+ or by showing that the VEV variation is negligible across the scanned range.","section":"§3, Eqs. (3.1)–(3.4)"},{"comment":"The exclusion of solutions with ξw > 0.99 as 'nonphysical' is arbitrary. The paper discards all runaway detonation candidates with this cutoff, and the scan statistics that lead to the conclusion that deflagration is the most prevalent mode are computed after this cut. No sensitivity to the cutoff value is given, and the only justification is a qualitative reference to Ref. [26]. Please show how the mode classification and the reported mean discrepancies change when the cutoff is varied (e.g., 0.95, 0.99, 0.999) or provide a physical criterion that fixes the cutoff.","section":"§4, Fig. 2 and text after Eq. (3.8)"},{"comment":"The WallGo validation uses a benchmark with ms = 120 GeV, λhs = 0.9, λs = 1.0, which lies outside the scan range ms ≥ 150 GeV stated in Eq. (4.5), and the reconstructed potential uses a different renormalization scheme and omits daisy resummation. The 1.22% agreement is encouraging for the numerical machinery, but it does not validate the fixed-VEV approximation in the region where the central discrepancies are claimed. Please state this limitation explicitly or add a validation point inside the scanned region.","section":"§4, WallGo cross-check paragraph"}],"minor_comments":[{"comment":"The abbreviation 'LET' is used in several places (e.g., 'the LET result' and 'LET approach') and should be 'LTE' for consistency with the rest of the text.","section":"§4 and Fig. 4 caption"},{"comment":"'Lorenz factor' should be 'Lorentz factor'.","section":"§2.1, after Eq. (2.4)"},{"comment":"The table header contains a formatting artifact 'T able 1', and the bold variables that are meant to indicate specified quantities are not visible in the rendered text; please ensure the formatting conveys the intended information.","section":"Table 1"},{"comment":"The scan range in the text states 150 GeV ≤ ms ≤ 500 GeV, but both panels of Fig. 2 show ms only up to about 325 GeV. Please clarify whether the axes are truncated or whether no solutions were found above that mass.","section":"Fig. 2 and Eq. (4.5)"},{"comment":"The parameters Ms, δ, and λ in the artificial bag-model benchmark are dimensionful but no units are given; specifying Ms in GeV and δ in GeV would avoid ambiguity.","section":"§3, Eq. (3.9)"},{"comment":"The abstract contains a stray '/github' token, and Appendix 6.1 shows code snippets but no repository URL. If the framework is intended to be publicly available, please provide a working link.","section":"Abstract and code appendix"}],"recommendation":"major_revision","confidential_remarks":"The paper is not circular and the central numerical comparisons are internally consistent, but the fixed-VEV approximation and the ξw > 0.99 cutoff are exactly the points on which the quantitative claims rest. I would request a direct test of the fixed-VEV approximation and a cutoff sensitivity study before publication. The WallGo benchmark agreement is a genuine strength, but it is outside the scanned parameter region. No citation or novelty concerns beyond the note added regarding Ref. [67]."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nHere's the short version: this is a useful, competent pipeline paper. It computes bubble wall velocities in LTE and kinetic energy fractions K directly from a full one-loop effective potential rather than a bag EOS, applies it systematically to the xSM, and quantifies how much the bag-model shortcut biases gravitational wave forecasts. The headline numbers—up to 48% in peak frequency and 90% in amplitude for deflagration—are the reason to read it.\n\nWhat is actually new is the systematic comparison. The ingredients aren't new: LTE matching from Ai, Garbrecht and Tamarit, full effective potentials, and the Espinosa et al. energy-budget machinery all exist. The contribution is putting them together with direct integration of K and doing a parameter scan rather than cherry-picked benchmarks. The internal checks are genuinely reassuring: an artificial bag model reproduces the expected profiles for all three modes, and the WallGo comparison at one xSM benchmark agrees to about 1.2%. The discussion of which approximation drives the GW error—wall velocity, not the EOS, in most of the scanned space—is clear.\n\nThe soft spots are real but not fatal. Section 3 fixes the VEVs in the two phases to reference-temperature values when evaluating p± and ρ±, with only the assertion that VEVs vary weakly. That approximation is load-bearing because the whole point of the full potential is temperature dependence. Worse, the code in Appendix 6.1 computes the speed of sound with valAt(T), i.e., temperature-dependent VEVs, so the EOS used at the wall is not the same as the EOS used in the profile integration. That inconsistency needs to be quantified. If VEVs shift substantially between Tn and the wall temperature, especially in deflagrations with shock heating, the solved ξw and K will shift and the 48%/90% numbers could move. The WallGo cross-check doesn't cover this either: the benchmark has ms = 120 GeV, while the scan starts at 150 GeV. The ξw > 0.99 cutoff is also arbitrary (cited to Ref. [26]) but not tested for sensitivity.\n\nNet: the central qualitative claim—bag-model forecasts can be off substantially in deflagration, mostly because of the wall velocity—probably survives. The specific percentages should be viewed as conditional until the fixed-VEV approximation is tested and, ideally, the code is released. This deserves a serious referee. I'd send it to review and ask for that test plus code; I'd also cite it in my own work.","headline":"Useful LTE pipeline paper with a credible quantitative claim that bag-model deflagration GW forecasts can be off by up to 48%/90%, but the untested fixed-VEV approximation in the wall matching needs to be quantified before those numbers are quoted.","tokens_in":21558,"tokens_out":3516,"would_cite":true,"duration_ms":68743,"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 full effective potential, not the bag model, should set the bubble wall velocity and kinetic energy fraction: in the xSM this shifts gravitational wave predictions by up to 48% in peak frequency and 90% in amplitude.","keywords":["bubble wall velocity","local thermal equilibrium","first-order phase transition","gravitational waves","effective potential","kinetic energy fraction","bag model","singlet extension of SM"],"falsifier":"Recompute the deflagration benchmark points with temperature-dependent field values inserted directly into the three wall-matching conditions; if the resulting wall velocities and kinetic-energy fractions move enough to erase the claimed 48%/90% differences from the bag model, the central comparison rests on the fixed-VEV approximation rather than on the full potential itself.","tokens_in":20560,"feed_emoji":"🌊","tokens_out":8117,"duration_ms":71146,"temperature":0.7,"pith_summary":"The paper establishes a framework for computing the bubble wall velocity and the kinetic energy fraction $K$ of a cosmological first-order phase transition directly from the full one-loop finite-temperature effective potential, using the local thermal equilibrium (LTE) approximation instead of treating the wall velocity as a free input. Applied to the singlet-extended Standard Model (xSM), the framework finds that deflagration is the most common fluid-motion mode. It then shows that gravitational wave spectra computed this way can differ from the commonly used bag-model-with-fitting-formula approach by up to 48% in peak frequency and 90% in peak amplitude in the deflagration regime, while the detonation regime shows much smaller differences. This matters because current forecasts for space-based gravitational wave detectors rely on exactly those approximate inputs, and the numbers quantify how much the forecasts can be trusted.","feed_headline":"Full potential shifts gravitational wave peaks by up to 90%","feed_subtitle":"Wall velocity from the full potential, not the bag model, changes deflagration predictions by 48-90%.","key_machinery":"The load-bearing mechanism is the LTE closure $\\gamma T = \\text{const}$ (equivalently, entropy conservation across the wall, $s\\gamma v = \\text{const}$), added as a third matching condition to the two standard energy-momentum flux conditions across the bubble wall; this removes one free input variable, so the wall velocity is solved for rather than guessed. The full effective potential supplies the equation of state through the thermodynamic relations for $p$, $\\rho$, and $w$, and a shooting method through the shock front fixes the position of the shock and the temperature profile. The framework is built on the xSM effective potential with one-loop Coleman-Weinberg, finite-temperature, and daisy-resummed corrections, and it discards detonation solutions with $\\xi_w > 0.99$ as unphysical.","core_discovery":"The paper's central claim is that the three wall-matching conditions—energy-flux conservation, momentum-flux conservation, and the LTE condition $\\gamma T = \\text{const}$ from entropy conservation—together with thermodynamic quantities derived from the full effective potential (pressure $p=-V_{\\rm eff}$, energy density $\\rho = V_{\\rm eff}-T\\,\\partial V_{\\rm eff}/\\partial T$, enthalpy $w=-T\\,\\partial V_{\\rm eff}/\\partial T$) determine the bubble wall velocity and the complete fluid profile without any ad hoc velocity input. Within the scanned xSM parameter space the authors find that deflagration is the prevalent steady-state mode, that the bag model approximates the full equation of state well for these benchmark points because the squared speed of sound stays close to $1/3$ in both phases, and that replacing the bag-model efficiency fits with integrated $K$ and LTE-derived wall velocities shifts gravitational wave peak predictions substantially in the deflagration regime (up to 48% in frequency and 90% in amplitude) but only mildly in detonation (6% and 18%). A notable secondary finding is that when the wall velocity is supplied as an external input in deflagration, the bag-model spectra resemble the LTE-based spectra more closely than the full-potential spectra with that same input velocity.","pith_inferences":["If the 48% and 90% discrepancies generalize beyond xSM, gravitational wave forecasts from the bag model carry a systematic error larger than the usual detector-sensitivity uncertainties, so a full-potential LTE computation should become the default for any model claiming a detectable signal.","The framework's fixed-VEV approximation can be tested directly: recomputing benchmark points with temperature-dependent field values in the matching conditions will show whether the quoted discrepancies survive; if they do not, the bag-model comparison needs to be redone.","The same machinery, applied to models with particles whose masses sit near the transition temperature (where the speed of sound deviates from $1/3$), should produce larger bag-model deviations than the xSM scan, giving a targeted prediction for where approximate forecasts most need revision."],"forward_implications":["Within the scanned xSM parameter space, the bag model is a reliable stand-in for the full equation of state for the fluid profiles, since the speed of sound in both phases is close to $1/3$.","The fitted efficiency formulas overestimate the kinetic energy fraction in deflagration (average peak-amplitude deviation about 14.3%) but agree with direct integration in detonation (about 2.9%).","The choice of wall velocity dominates the uncertainty in the predicted gravitational wave spectra: fixing it to a representative value (0.3 for deflagration, 0.9 for detonation) instead of using the LTE result produces deviations that track variations in $K$.","Detonation solutions are unphysical when the LTE wall velocity approaches the speed of light ($\\xi_w > 0.99$); such parameter points are excluded, and LTE is only valid for small transition strengths $\\alpha_N$.","In deflagration scenarios where the wall velocity is treated as an input parameter, mapping the model to the bag model gives spectra closer to the LTE-based calculation than using the full potential with the same input velocity."],"supporting_citations":[{"why":"Supplies the local thermal equilibrium approximation and the entropy-conservation matching condition $\\gamma T = \\mathrm{const}$ that the framework adds to the hydrodynamic equations.","marker":"[26]"},{"why":"Provides the energy-budget formalism, the efficiency-factor fitting formulas, and the wall-matching conditions that the paper compares against.","marker":"[51]"},{"why":"Open-source package used as the benchmark: for the same xSM parameter point it gives wall velocity 0.6203 vs the authors' 0.6127, validating the code.","marker":"[32]"},{"why":"Public phase-transition code used to reconstruct the xSM effective potential for the cross-check with the benchmark package.","marker":"[60]"},{"why":"Supplies the xSM effective potential with one-loop, thermal, and daisy-ring corrections used in the parameter scan and gravitational wave predictions.","marker":"[56]"},{"why":"Provides the sound-wave gravitational wave spectrum formula (amplitude and peak frequency) that the paper evaluates for every parameter point.","marker":"[61]"},{"why":"Supplies the definition of the Jouguet velocity as the minimum detonation velocity, used to classify which fluid motion modes are allowed.","marker":"[49]"},{"why":"Discussed as the reference for effects beyond the bag equation of state, including non-constant temperature profiles and the speed-of-sound corrections to the efficiency factor.","marker":"[52]"}],"fun_headline_variants":["Full potential vs bag model: deflagration GW peaks shift up to 90%","LTE wall velocity plus full potential changes deflagration GW peaks by 90%","Bag model underestimates deflagration GW amplitude by up to 90%","Deflagration prevalent: full potential alters GW predictions up to 90%","Full potential shifts GW peaks 48% in frequency, 90% in amplitude"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The matching conditions assume the vacuum field values in the two phases stay fixed at their reference-temperature values even though the effective potential and speed of sound used elsewhere are temperature-dependent.","fun_headline_variants_meta":{"raw":{"variants":["Full potential vs bag model: deflagration GW peaks shift up to 90%","LTE wall velocity plus full potential changes deflagration GW peaks by 90%","Bag model underestimates deflagration GW amplitude by up to 90%","Deflagration prevalent: full potential alters GW predictions up to 90%","Full potential shifts GW peaks 48% in frequency, 90% in amplitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00039,"raw_usage":{"total_tokens":2101,"prompt_tokens":1040,"completion_tokens":1061,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":954}},"tokens_in":656,"tokens_out":1061,"duration_ms":8839,"temperature":1.0,"reasoning_tokens":954,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:11:58.291206+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the deflagration benchmark points with temperature-dependent field values inserted directly into the three wall-matching conditions; if the resulting wall velocities and kinetic-energy fractions move enough to erase the claimed 48%/90% differences from the bag model, the central comparison rests on the fixed-VEV approximation rather than on the full potential itself.","supporting_citations":[],"review_version":1}