{"id":"ae39876c-dcf2-4776-9bd0-c794228be300","arxiv_id":"2412.09588","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Holographic simulations show supercooled bubble walls in a QCD-like fluid reach at most about 0.13 times light speed, while superheated walls reach only about 0.03.","lead":"Physicists simulated the formation and growth of bubbles in a dense, strongly interacting fluid whose behavior is designed to mimic the conjectured phase changes of quark matter. They found the bubble walls move slowly, especially for superheated bubbles, which affects predictions for gravitational waves from neutron star mergers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No convergence tests are reported for the new Maxwell-extended Jecco code, so the central numerical values (v_w^max ≈ 0.126 and 0.032) and the superheated/supercooled asymmetry are not yet secured against discretization error.","rationale":"The reader's verdict is CONDITIONAL, and the main numerical-convergence concern supports that conditionality. The reader's stated weakest assumption is the qualitative nature of the holographic model; that is a real limitation, but it is explicitly disclosed in Section 2 and does not invalidate the paper's stated goal of studying a QCD-like phase diagram. The more actionable weakness is the absence of any convergence evidence for the new Maxwell-extended code, which directly bears on the reliability of the quantitative centerpiece: the reported maximum velocities and the supercooled/superheated asymmetry. Since the manuscript already lacks public code and data, a targeted resolution study is the natural way to settle this concern. A possible secondary issue is an apparent internal inconsistency in Section 6, where the text says one would expect v_w to decrease with δ after stating the large-jump approximation gives v_w ∝ δ for superheated bubbles; this affects the comparison with Ref. [27] but not the primary bubble-velocity result. If the proposed convergence tests pass, the qualitative findings (slow walls, superheated walls slower than supercooled ones) are likely robust; if they fail, the numerical values would need revision. Either way the verdict remains CONDITIONAL, so no change to the reader's recommendation is needed.","tokens_in":21157,"tokens_out":5695,"duration_ms":60094,"concrete_test":"Re-run two representative steady-state simulations (for example supercooled c3 and superheated h3 of Fig. 11) at twice the radial and transverse resolution and at half the time step, using identical initial data; extract v_w from the inflection-point trajectory over the same self-similar ξ interval and report the relative change in v_w. Independently, repeat one of these runs with the initial Gaussian perturbation amplitude changed by a factor of two to verify the claimed insensitivity to initial data. If either test shifts v_w by more than approximately 10%, the quantitative values in Eqs. (6.1)-(6.2) should be revised or presented with explicit uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result is quantitative: the paper reports three-digit maximum wall velocities in Eqs. (6.1)-(6.2) and a monotonic increase with supercooling/superheating in Figs. 17-19. However, no resolution study, no error bars, and no convergence checks are presented for the Einstein-scalar-Maxwell system introduced in Section 2 and implemented in Appendix A. The new ingredient relative to the previously validated Jecco code of Ref. [41] is the Maxwell field; Appendix A shows that the equations for Φ_t and F become coupled, as do those for dot ϕ, dot A_x and dot B, so the established discretization is not automatically reliable for this system. Because superheated walls are very slow (v_w ≲ 0.032), numerical dissipation can in principle bias their apparent speed, possibly exaggerating or even creating the asymmetry with supercooled bubbles. The self-similar collapse of snapshots in Figs. 15-16 is encouraging but does not quantify the error in v_w. Without a grid-refinement check, the statement 'within our set of simulations, the maximum velocity...' characterizes a particular numerical setup rather than a robust physical result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses a five-dimensional Einstein-scalar-Maxwell holographic model to study bubble dynamics in a theory whose equilibrium phase diagram contains a first-order transition line ending at a critical point, qualitatively resembling the conjectured QCD phase diagram at finite baryon density. The authors first construct the thermodynamic phase diagram, identify metastable regions, and then evolve planar bubbles for both supercooled and superheated cases, extracting terminal wall velocities from the steady-state profiles. They report maximum velocities v_w^max ≈ 0.126 for supercooled and ≈0.032 for superheated bubbles, claim a monotonic increase of the velocity with the amount of supercooling/superheating, and compare their results with a hydrodynamic large-jump approximation and with a linear relation between velocity and pressure difference. Implications for gravitational wave emission in neutron star mergers are discussed.","tokens_in":21363,"tokens_out":6493,"duration_ms":63463,"significance":"If the numerical results are robust, this is the first microscopic, strongly coupled computation of bubble wall velocities in a holographic model with a QCD-like phase diagram, thereby extending Ref. [30] to nonzero baryon chemical potential and to superheated bubbles. The qualitative messages — that the walls remain slow compared to the speed of light, that the velocity increases away from the transition line, and that superheated bubbles are markedly slower than supercooled ones — are interesting and potentially relevant for gravitational wave phenomenology. Strengths of the paper include the explicit derivation of the new Maxwell-sector equations in Appendix A, the use of equilibrium data in the comparisons, the self-similar-profile checks, and the systematic comparison with independent approximations. The principal weakness is that the central quantitative outputs are not yet supported by convergence tests, error estimates, or an equally thorough initial-condition independence check as in Ref. [30].","major_comments":[{"comment":"The central quantitative results — v_w^max ≈ 0.126 and 0.032, and the velocities shown in Figs. 17–20 — are reported without any resolution study, convergence test, or error estimate for the new Einstein-scalar-Maxwell system. The Maxwell extension changes the nested structure of the equations of motion (A.1): the equations for Ψ_t and F become coupled, and the dot variables φ̇, Ȧ_x and Ḃ are coupled, so the validation of the neutral Jecco code in Ref. [41] does not automatically transfer to this case. This is especially concerning for the superheated velocities, which are of order 0.03 and could be significantly biased by numerical dissipation or by finite-domain effects. I ask for at least a grid-refinement test and a domain-size test for representative supercooled and superheated simulations, with the resulting error bars propagated to all plotted velocities.","section":"Sec. 6, Eqs. (6.1)–(6.2); App. A"},{"comment":"The statement that the terminal velocity depends only on the metastable state outside the bubble is used to organize all the results in Section 6, but in this paper it is not tested with the same thoroughness as in Ref. [30]. The text explicitly says that the expected independence has been checked only in 'initial investigations' and that an analysis as exhaustive as Ref. [30] was not performed. Since the model now includes a charged sector and the perturbations are introduced in both a4 and At,2, the initial-data independence should be demonstrated for at least a few representative states by varying the amplitude and width of the Gaussian perturbation; otherwise the extracted v_w values could carry an undetermined dependence on the initial condition.","section":"Sec. 5, initial conditions"},{"comment":"The wording 'the wall velocity increases monotonically with the amount of supercooling or superheating' is stronger than the evidence presented. The support consists of a small number of points without error bars, ordered by a Euclidean distance in (E/Λ^4, ρ/Λ^3) that the authors themselves describe as having no intrinsic physical meaning, and Fig. 20 shows considerable scatter around the linear fits. The monotonicity claim should either be presented with quantitative tolerances, including the scatter of all simulations and the sensitivity to the chosen ordering coordinate, or be softened to a monotonic trend.","section":"Sec. 6, Figs. 18–19"}],"minor_comments":[{"comment":"The sentence 'this multivalued region coincides exactly coincides with the region where metastable states exist' contains a duplicated word; it should read 'coincides exactly with'.","section":"Sec. 3"},{"comment":"The caption contains 'as as a function of the self-similar parameter'; remove the duplicated 'as'.","section":"Fig. 16 caption"},{"comment":"Refs. [30] and [38] list the same paper and should be merged into a single reference.","section":"References"},{"comment":"The phrase 'highly sensitively to the velocity' should be 'highly sensitive to the velocity'.","section":"Abstract"},{"comment":"The horizontal axis of both panels in Fig. 20 is not labeled; please add the definition of the plotted variable to the figure.","section":"Sec. 6, Fig. 20"},{"comment":"The suggestion that the large-jump approximation becomes accurate for QCD because the error decreases with δ is an extrapolation: the data only reach δ ≈ 0.3, while the QCD estimate quoted is δ ≈ 0.1. Please phrase this as an extrapolation or add supporting data at smaller δ.","section":"Sec. 6, Fig. 21"}],"recommendation":"major_revision","confidential_remarks":"This is a competent and interesting paper, and the main request for major revision is numerical robustness rather than conceptual soundness. I would not reject on the current evidence; the authors should be able to supply convergence tests and initial-data checks within a normal revision cycle. The fit to the journal is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the finite-density, superheated extension of the holographic bubble wall program, and it is a real step forward. For the first time you get microscopic wall velocities in a QCD-like phase diagram with nonzero baryon number: about 0.126 max for supercooled bubbles, 0.032 for superheated, with monotonic dependence on metastability depth. The authors also show that the large-jump hydrodynamic approximation misses the exact values by 35–70% for supercooled bubbles and fails for superheated ones. That second result is probably the most durable part of the paper.\n\nWhat is genuinely new: Ref. [30] did zero chemical potential and only supercooled bubbles. Here the bulk is Einstein-scalar-Maxwell, chosen to have a critical point and a line of first-order transitions, and the Jecco code is extended to evolve the Maxwell field. The physics output is model-level, not quantitative QCD: Section 2 says this explicitly, and the gravitational-wave numbers are order-of-magnitude. The paper deserves credit for not overstating that.\n\nSoft spots, proportional. The stress-test concern is correct. The two headline velocities have no error bars and no grid-refinement study. The Maxwell extension is precisely the new part of the code, and the coupled equations in Appendix A mean the old neutral-code validation does not automatically transfer. \"Within our set of simulations\" is an honest qualifier, but it also means the maximum is a grid-search result rather than a controlled extremum. There is also a softer statement on initial-condition independence than in Ref. [30]. None of that undermines the monotonic-scaling story, which is supported by many data points and by the self-similar profiles. It does mean the precise magnitudes should not be quoted until a convergence test appears. The self-citations are mostly to their own code and prior results, which is appropriate; the missing validation is not hidden by them.\n\nWho benefits: gravitational-wave modelers needing a concrete strong-coupling benchmark, and holographers working on phase transitions. I would send this to a serious referee. The fix is straightforward: one convergence/monitoring study, error estimates on the velocities, and ideally release of the evolution data.","headline":"A genuinely new finite-density holographic bubble dynamics computation with a clean superheated/supercooled asymmetry; referee it, but require convergence checks before the velocity magnitudes are quoted.","tokens_in":21932,"tokens_out":3085,"would_cite":true,"duration_ms":30633,"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 reports the first microscopic, holographic computation of bubble wall velocities in a phase diagram qualitatively mirroring QCD, finding slow walls with maximum speeds of about 0.126 c (supercooled) and 0.032 c (superheated)…","keywords":["holography","bubble wall velocity","first-order phase transition","QCD phase diagram","neutron star mergers","gravitational waves","supercooled bubbles","superheated bubbles"],"falsifier":"A first-principles simulation of the same type in a holographic model quantitatively matched to QCD thermodynamics—or a future lattice determination of the QCD transition at high baryon density—that found superheated walls moving faster than supercooled walls at comparable metastability, or wall speeds above about half the speed of light in the most metastable states, would contradict the paper's central results.","tokens_in":20919,"feed_emoji":"🫧","tokens_out":8870,"duration_ms":74352,"temperature":0.7,"pith_summary":"This paper reports the first microscopic computation of bubble-wall velocities in a holographic theory whose phase diagram qualitatively mirrors the conjectured QCD phase diagram at nonzero baryon density. It simulates bubbles of the stable phase growing inside metastable states in both directions—supercooled and superheated—and finds that in steady state the wall velocity is always small: at most about 0.126 times the speed of light for supercooled bubbles and 0.032 for superheated ones. The velocity rises monotonically with the amount of supercooling or superheating and with the distance from the critical point, and superheated walls are systematically slower than supercooled ones for comparable metastability. These results matter because the gravitational-wave signal expected from neutron star mergers is highly sensitive to the wall velocity, so the calculation sharpens predictions for what future detectors might see.","feed_headline":"QCD-like bubble walls crawl under 13% of light speed","feed_subtitle":"First holographic clocking: supercooled walls reach 0.126 c, superheated only 0.032 c.","key_machinery":"The central object is the five-dimensional Einstein-scalar-Maxwell holographic model, in which a scalar field dual to a dimension-three operator and a Maxwell field dual to baryon number are chosen so the boundary theory has the desired phase diagram: a crossover at low chemical potential that turns into a line of first-order phase transitions at $\\mu \\gtrsim 1.09\\Lambda$, ending at a critical point at $(T,\\mu)=(0.34,1.09)\\Lambda$, with well-defined metastable regions found from the Hessian of the free energy. The dynamics are produced by seeding a homogeneous metastable state with a localized Gaussian perturbation in the energy and charge densities, whose sign selects a superheated or supercooled bubble, and then evolving the full bulk equations in ingoing Eddington-Finkelstein coordinates with a nested time-integration scheme. The steady-state wall velocity is read off from the motion of the inflection point of the energy-density profile once the flow becomes self-similar. This setup is what makes the first microscopic, parameter-free extraction of $v_w$ possible in a QCD-like phase diagram.","core_discovery":"The paper's central result is that, in a five-dimensional Einstein-scalar-Maxwell model chosen so the boundary theory has a line of first-order phase transitions ending in a critical point—like the conjectured QCD line at finite baryon chemical potential—expanding bubbles of the stable phase settle into slow, steady-state motion. For every metastable state studied, the wall reaches a terminal velocity, extracted from the inflection point of the energy-density profile, and this velocity grows monotonically with the degree of supercooling or superheating, vanishing at the first-order line and growing with distance from the critical point. The largest velocities found are $v_w^{\\max} \\simeq 0.126$ for supercooled bubbles and $v_w^{\\max} \\simeq 0.032$ for superheated bubbles, so superheated walls are markedly slower. Both types of bubbles are deflagrations: supercooled bubbles push fluid outward and leave an overdense shell ahead of the wall, while superheated bubbles absorb energy from the outside and leave an underdense shell. The paper also compares the simulated velocities with existing estimates, finding that the large-jump-in-degrees-of-freedom approximation underpredicts supercooled velocities by 35–70% and becomes unreliable for superheated bubbles in the regime probed.","pith_inferences":["If real QCD walls are as slow as these, the gravitational-wave strain from neutron star mergers could be smaller than estimates that assume relativistic walls, and the difference between the two directions implies the signal may distinguish heating from cooling phases of the merger.","The approximate linear relation between $v_w$ and $\\Delta P/E$ (or its baryon-density generalization) holds only roughly in the supercooled data; the scatter suggests that a quantitative relation for superheated walls may require genuinely microscopic input rather than equilibrium quantities alone.","Because the model was selected for numerical tractability rather than quantitative QCD matching, the qualitative ordering—supercooled faster than superheated, monotonic growth with metastability—is the transferable content; the specific speeds would likely shift in a model calibrated to QCD thermodynamics.","A natural testable extension is to repeat the same holographic evolution for spherical bubbles and for bubble collisions, since the paper's planar isolated-bubble setup isolates $v_w$ but leaves the sound-wave and collision contributions to the gravitational-wave spectrum uncomputed."],"forward_implications":["Maximum steady-state bubble wall speeds are $v_w^{\\max}\\simeq 0.126$ for supercooled and $v_w^{\\max}\\simeq 0.032$ for superheated bubbles, so walls in this strongly coupled theory are slow compared with the speed of light.","Wall velocity increases monotonically with the amount of supercooling or superheating and with distance from the critical point, and it goes to zero at zero metastability.","Superheated bubbles are slower than supercooled bubbles for comparable metastability, and the two types produce opposite fluid-flow patterns: overdense shells for supercooled, underdense shells for superheated.","The large-jump-in-degrees-of-freedom approximation underestimates supercooled wall velocities by 35–70% in this model, with the error decreasing as the jump parameter $\\delta$ decreases; the same approximation is not applicable to the simulated superheated regime.","Applied to neutron star mergers, the results imply gravitational-wave emission dominated by slow walls, with the peak frequency in the MHz range and an amplitude sensitive to the computed velocities."],"supporting_citations":[{"why":"The original holographic computation of the bubble wall velocity that this paper extends to superheated bubbles and nonzero baryon chemical potential.","marker":"[30]"},{"why":"Supplies the Einstein-scalar-Maxwell construction and the holographic critical point whose phase diagram is mirrored here.","marker":"[31]"},{"why":"Provides the superpotential family and the scalar potential used to define the model.","marker":"[32]"},{"why":"Identifies neutron star mergers as sources of megahertz gravitational waves from QCD-like phase transitions, the observational target of the calculation.","marker":"[5]"},{"why":"Gives the hydrodynamic description of relativistic superheated bubbles that the paper compares with its microscopic results.","marker":"[16]"},{"why":"Provides the large-jump-in-degrees-of-freedom prediction for the wall velocity that is tested against the simulations.","marker":"[27]"},{"why":"Derives the approximately linear relation between wall velocity and normalized pressure difference that the paper examines.","marker":"[26]"},{"why":"Earlier holographic study of spinodal gravitational waves that found the linear velocity-versus-pressure-difference behavior used as a comparison.","marker":"[40]"}],"fun_headline_variants":["Holographic sim: QCD bubble walls top out at 0.126c","Supercooled bubble walls reach 0.126c, superheated only 0.032c","First holographic bubble dynamics: max wall speed 0.126c","QCD bubble expansion computed: wall velocity max 0.126c","Holography measures QCD bubble wall speeds: 0.126c max"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything said about QCD rests on the assumption that this particular holographic model—chosen for numerical tractability, with its specific superpotential and gauge coupling—mirrors the real QCD phase diagram closely enough that its bubble dynamics are representative, despite not being quantitatively matched to QCD.","fun_headline_variants_meta":{"raw":{"variants":["Holographic sim: QCD bubble walls top out at 0.126c","Supercooled bubble walls reach 0.126c, superheated only 0.032c","First holographic bubble dynamics: max wall speed 0.126c","QCD bubble expansion computed: wall velocity max 0.126c","Holography measures QCD bubble wall speeds: 0.126c max"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001098,"raw_usage":{"total_tokens":4573,"prompt_tokens":925,"completion_tokens":3648,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":541,"completion_tokens_details":{"reasoning_tokens":3540}},"tokens_in":541,"tokens_out":3648,"duration_ms":25338,"temperature":1.0,"reasoning_tokens":3540,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:54:54.340756+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A first-principles simulation of the same type in a holographic model quantitatively matched to QCD thermodynamics—or a future lattice determination of the QCD transition at high baryon density—that found superheated walls moving faster than supercooled walls at comparable metastability, or wall speeds above about half the speed of light in the most metastable states, would contradict the paper's central results.","supporting_citations":[],"review_version":1}