{"id":"417e31ba-1571-4d17-aaad-55c102cb6375","arxiv_id":"2502.07879","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In an effective holographic model, expanding and cooling plasma forms 'hot remnants' that stay at the critical temperature while shrinking, then heat up during dissolution.","lead":"This paper tracks what happens when a hot, expanding plasma cools through a first-order phase transition, using a simplified model calibrated to a holographic description of quark confinement. It finds that pockets of hot plasma can survive, stay near the transition temperature, and then heat up while dissolving, instead of simply cooling and turning into bubbles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dissolution-stage scaling and heating above Tc follow from exact entropy conservation in the zero-dissipation model; no test shows these survive in a real plasma with viscosity and order-parameter friction.","rationale":"The paper's central claim has two parts: (i) hot remnants persist at T ≈ Tc during shrinking, and (ii) during dissolution the residual plasma heats up above Tc. Part (i) is supported by fully converged simulations (Figs. 1,3,4,5) and by simple energy conservation; I see no serious threat to it. Part (ii) is derived from two ingredients: exact entropy conservation (Supplemental) and the virial approximation using the potential VTOT extrapolated beyond its fitted range. The reader flagged no-dissipation; I agree that it is the more fundamental of the two because it enters the analytic derivation directly. The potential extrapolation affects the quantitative T(t) curve but not necessarily the existence of heating: any smooth potential with the fitted first T-derivative at Tc will have its interior zero γ*(T) move toward γ=1 as T rises, so the qualitative heating is likely generic. Entropy non-conservation, by contrast, removes the derivation of the 1/√τ scaling entirely; with dissipation the system has an additional entropy-production channel and the late-time behavior is unconstrained by the model. A concrete numerical experiment adding a minimal friction term would settle whether the dissolution heating survives. The paper deserves credit for the fully converged shrinking-stage results and for providing analytic scalings; the conditional verdict is appropriate until such a test is performed or code/data are provided. Thus verdict_should_be is UNCHANGED.","tokens_in":11908,"tokens_out":17757,"duration_ms":165519,"concrete_test":"Run the same boost-invariant dissolution setup (Fig. 6) with a minimal friction term in the γ equation of motion, e.g., replace Eq. (A.8) with c ∂^2γ − V′_TOT = ηγ (u·∂)γ, choosing ηγ so that the terminal domain-wall velocity equals the holographic value quoted in [16], and add a corresponding viscous contribution to T^{μν} with η/s = 1/4π. Record T(τ) and (1−γ)max(τ) in the dissolving regime. If T still rises above Tc and (1−γ)max still decays as τ^{−1/2}, the no-dissipation assumption is not load-bearing. If T drops below Tc or the decay changes qualitatively, the central claim of heating above Tc is an artifact of the conservative model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Supplemental Material derives the exactly conserved entropy current jμ = −∂T VTOT(γ,T) uμ; the letter invokes this conservation law to obtain the dissolution scaling 1−Γ ∼ 1/τ (equivalently 1−γ ∼ 1/√τ) and, via the virial balance VTOT(γ*,T)=0, the counterintuitive temperature rise above Tc shown in Figs. 6 and 9. This conservation law is a direct consequence of the model having no dissipation: the authors state 'we do not have any dissipation terms which are usually added ... as our model by itself reproduces correct bubble wall velocities in holographic theories' (p.2). The shrinking stage at T ≈ Tc is governed by energy conservation and latent heat and is therefore robust. The load-bearing step is the dissolution stage: real plasmas have shear and bulk viscosity and order-parameter friction, so ∂μjμ > 0, and the entropy budget that forces 1−Γ ∼ 1/τ is modified. If a dissipative term of holographically motivated size is added, the deconfined fraction may dilute at a different rate and the virial condition driving T above Tc need not hold; the remnant could instead cool and nucleate bubbles. The authors themselves flag that the onset of dissolution is sensitive to the numerical filter (Supplemental), so the zero-dissipation idealization adds a further unquantified source of uncertainty for the headline claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a first-order confinement/deconfinement phase transition in an expanding medium using an effective boundary model fitted to the holographic Witten model. The model couples boost-invariant (or FRW) hydrodynamics to a scalar order parameter γ that interpolates between deconfined and confined phases. The authors report a generic dynamical outcome: after supercooling and bubble nucleation, trapped regions of deconfined plasma form 'hot remnants' that cool to Tc, shrink while releasing latent heat, and then dissolve, with the remaining deconfined fraction decreasing as 1−γ ~ 1/√τ while the temperature rises above Tc. The shrinking and dissolution stages are reproduced in both flat-space boost-invariant expansion and expanding FRW spacetime. The paper interprets the shrinking stage via energy conservation and the dissolution stage via the exactly conserved entropy current combined with a virial approximation for the scalar profile.","tokens_in":12190,"tokens_out":2488,"duration_ms":32875,"significance":"If the findings are robust, they challenge the standard picture in which a first-order transition in an expanding medium proceeds purely via bubble nucleation and cooling below Tc. The hot-remnant phenomenon and the reheating during dissolution could be relevant for heavy-ion collisions and for early-universe phase transitions. The paper provides transparent physical arguments, two independent numerical setups, and analytic scalings (L ~ 1/τ, 1−γ ~ 1/√τ) that are directly confronted with the numerics. The authors also supply a detailed supplemental description of the numerical scheme, including regularization, grid settings, and convergence checks for several runs. These are concrete strengths. The main limitations are the zero-dissipation idealization, the reliance on a potential fitted only near Tc, and the sensitivity of the dissolution onset to the numerical filter; these affect the headline claim that the phenomenon is 'very generic'.","major_comments":[{"comment":"The dissolution-stage scaling 1−Γ ~ 1/τ (equivalently 1−γ ~ 1/√τ) and the associated temperature rise are derived from the exactly conserved entropy current j^μ = −∂_T V_TOT u^μ (Eq. A.14 in the Supplemental Material). This conservation law requires the absence of dissipation, as the authors explicitly state: 'we do not have any dissipation terms which are usually added'. Real plasmas possess shear and bulk viscosity and order-parameter friction, which would make ∂_μ j^μ > 0. The paper does not quantify how the reported scalings and the reheating would be modified by a dissipative term of physically motivated size. Since the abstract and discussion claim that hot remnants appear 'very generically', this idealization is load-bearing; a test with a small dissipative term, or at least an order-of-magnitude estimate, is needed to support the extrapolation beyond the zero-dissipation model.","section":"Supplemental Material, 'Entropy current' and main text, p. 2"},{"comment":"The 'virial theorem prediction' for the temperature rise during dissolution reads off the zeros of V_TOT(γ*, T)=0, but V_TOT is the same fitted potential whose form is used 'also away from Tc' even though it 'was fitted only in the neighbourhood of Tc = 1 up to its first T-derivative'. The reheating above Tc is therefore not an independent prediction: it follows from an extrapolation of the fitted potential. The paper should either test the sensitivity of the dissolution-stage temperature rise to the form of V_TOT away from Tc (e.g., by varying β(γ) or adding higher-order terms) or soften the claim that the reheating is a generic consequence of the holographic model rather than of the chosen extrapolation.","section":"p. 2, paragraph on the potential; Figs. 6 and 9"},{"comment":"The dissolution-onset timescale is reported to be sensitive to the numerical filter: 'the results in Figs. 6, 7, 8 and 9 suffer from a mild uncertainty in the timescale where the γ field profile starts to dissolve. This timescale turns out to be sensitive to the filter we are using.' Figures 6 and 9 are exactly the figures that establish the counterintuitive reheating stage. Since the filter is a numerical artifact, the uncertainty in the onset time propagates to the duration and quantitative details of the dissolution stage. The paper should provide a quantitative estimate of how much this uncertainty affects the extracted scalings and the temperature-rise curve, rather than only stating that the predictions of the virial theorem are reproduced.","section":"Supplemental Material, 'Some details on the numerical simulations'"}],"minor_comments":[{"comment":"The potential V_TOT(γ, T) is written with a term '+1' at the end; it would help to clarify that this constant shifts the overall zero of the potential and does not affect the equations of motion.","section":"p. 3, Eq. (3)"},{"comment":"The profile ansatz 1−γ ∼ b(τ)/cosh^4(d(τ) q_* x / 4) is introduced with parameters b(τ) and d(τ); the rationale for the specific power 4 and the factor 1/4 in the argument would be clearer if tied directly to the potential approximation in Eq. (14).","section":"p. 4, Eq. (12)"},{"comment":"The text says 'expanding s = −∂_T V_TOT around γ = 1 gives 1−Γ ∼ 1/τ' but does not show the intermediate steps; for a letter this is acceptable, but a one-line derivation in the Supplemental Material would improve reproducibility.","section":"p. 4, after Eq. (14)"},{"comment":"The statement 'we checked that similar results hold also e.g. for a(t) ∼ e^{Ht} with small H' is not accompanied by any quantitative detail; adding a brief comment on the range of H and the observed scalings would strengthen the genericity claim.","section":"p. 5, FRW section"},{"comment":"The reference list would benefit from the page numbers or article numbers for the arXiv:2411.17806 paper (Ref. [23]) and for the published versions of Refs. [27] and [28].","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is interesting and the core 'hot remnant' phenomenon is convincingly demonstrated within the zero-dissipation effective model. My main concern is the gap between the idealized model and the 'very generic' claim in the abstract, especially because the dissolution-stage reheating is tied to the exact entropy conservation that holds only in the absence of dissipation. The numerical-filter sensitivity in the dissolution onset is also understated in the main text. These are fixable with additional tests or a more careful framing, but they are load-bearing for the headline claim, hence major revision rather than minor. I do not see grounds for rejection: the central mechanism is physically transparent and the numerical evidence with the stated caveats is consistent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the identification of a three-stage evolution for deconfined remnants in an expanding medium: reheating to Tc, shrinking at constant temperature, then dissolution with the residual plasma heating above Tc. The shrinking stage is the solid core of the paper. The L ~ 1/tau scaling in boost-invariant flow and L ~ 1/a^3 in FRW follow from simple energy conservation plus latent heat, and they match two independent numerical setups. That part deserves credit and is likely robust.\n\nThe dissolution stage is the soft spot, and the stress-test note lands on it correctly. The 1-gamma ~ 1/sqrt(tau) scaling and the temperature rise above Tc are derived from exact entropy conservation, which is built into the model by having no dissipation terms. Real plasmas have viscosity and order-parameter friction, so the entropy budget that forces these scalings will be modified. The authors are upfront about the absence of dissipation, and they motivate it by the model's success in reproducing holographic bubble wall velocities, but that success does not automatically extend to the late-time dissolution of a small remnant. They also admit the onset of dissolution is sensitive to the numerical filter. So the headline claim that remnants 'heat up' during dissolution is conditional on the zero-dissipation idealization, not an established prediction for real systems.\n\nThe circularity concern is real but mild. The parameters are fit to the AMW domain wall, and the virial-theorem prediction of heating reads off the same fitted potential. To that extent it is not an independent test. But the hot-remnant phenomenon itself is an emergent dynamical outcome, not put in by hand, and the scaling laws in the shrinking stage are derived from conservation principles that don't depend on the fine details of the fit. I would not call the central mechanism circular.\n\nWho is this for? People working on holographic models of phase transitions, early-universe cosmology, and heavy-ion phenomenology. The paper is clearly written and the physics arguments are transparent. It deserves a serious referee: the phenomenon is interesting, the shrinking-stage results are solid, and the dissolution-stage predictions are clearly flagged as model-dependent. The main things I would want from a revision are a quantification of how dissipation would change the dissolution scalings and either code release or more detail on the numerics.\n\nMy recommendation: send it to peer review. It is a legitimate model-based result with a clear new qualitative phenomenon, and the soft spots are addressable rather than fatal.","headline":"A credible model-based discovery of hot remnants in expanding plasma, with robust scaling laws in the shrinking stage and a load-bearing zero-dissipation assumption in the dissolution stage.","tokens_in":12731,"tokens_out":616,"would_cite":true,"duration_ms":8197,"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":"In an expanding medium, a first-order phase transition can leave hot remnants of the high-energy phase that remain at the transition temperature while shrinking, and then heat up further as they dissolve.","keywords":["first-order phase transition","deconfinement","hot remnants","holographic confinement model","boost-invariant expansion","FRW cosmology","entropy conservation","domain wall"],"falsifier":"Rerun the boost-invariant simulation with the dissipative or friction terms that the paper deliberately omits, using the same initial conditions: if a blob of deconfined matter near $T_c$ then cools below $T_c$, nucleates bubbles, or dissolves without reheating above $T_c$, the claim of generic hot remnants is falsified; a weaker test is to measure the entropy production rate, since the scaling laws require exact conservation.","tokens_in":11673,"feed_emoji":"🔥","tokens_out":7120,"duration_ms":58807,"temperature":0.7,"pith_summary":"The paper argues that the standard picture of a first-order phase transition in an expanding medium — supercool, then nucleate bubbles of the low-temperature phase — is incomplete. Using a holography-inspired effective model of the confinement/deconfinement transition, it shows that regions of the high-energy phase can persist as 'hot remnants' whose temperature stays pinned near the critical temperature $T_c$ even as the medium keeps expanding. These remnants shrink rather than nucleate bubbles; once small enough, they dissolve, and the remaining plasma heats up above $T_c$. The same qualitative stages appear in boost-invariant flat-space expansion and in cosmological FRW expansion, suggesting a generic mechanism rather than a peculiarity of a specific geometry. The paper's conclusion matters because hot remnants would alter predictions for heavy-ion collisions and for first-order phase transitions in the early universe.","feed_headline":"Hot plasma remnants survive expansion without cooling","feed_subtitle":"First-order phase transitions may leave long-lived hot remnants that shrink, then heat up.","key_machinery":"The central object is an effective boundary model in which the total energy-momentum tensor is a mixture of a perfect-fluid deconfined phase and a trivial confined phase, controlled by a dynamical order parameter $\\gamma$ (the 'mixing fraction'), plus a domain-wall contribution. The dynamics of $\\gamma$ come from a Lagrangian with double-well potential $V_{TOT}(\\gamma,T)$, whose explicit temperature dependence couples the order parameter to hydrodynamics through the formalism of [13]. This coupling is such that the model needs no dissipation terms and still reproduces holographic domain-wall velocities [16]. The argument is carried by the exactly conserved entropy current $j^\\mu = -\\partial_T V_{TOT}(\\gamma,T)\\,u^\\mu$: entropy conservation plus the virial approximation $c\\,\\gamma'^2=2V_{TOT}(\\gamma,T)$ yields the shrinking and dissolution scalings and the late-time heating.","core_discovery":"The central discovery is that in an expanding medium, first-order confinement/deconfinement transitions do not simply proceed by supercooling and bubble nucleation. Instead, trapped regions of the deconfined high-energy phase remain at a temperature near the critical temperature $T_c$ while the rest of the system cools: latent heat released as the surrounding confined phase expands keeps the remnant hot. These hot remnants shrink, and only when their width reaches the domain-wall scale do they dissolve; during dissolution the residual deconfined plasma is heated further, moving above $T_c$. The paper derives scaling laws for each stage — $L(\\tau)\\sim 1/\\tau$ for boost-invariant expansion and $L(t)\\sim 1/a(t)^3$ in an expanding FRW universe, then $1-\\gamma\\sim1/\\sqrt{\\tau}$ with rising temperature — and shows that all of it follows from exact entropy conservation in a dissipation-free hydrodynamic model with a dynamical mixing fraction.","pith_inferences":["If the mechanism is as generic as claimed, first-order transitions in the early universe could leave long-lived hot regions whose subsequent decay produces gravitational waves and particle abundances different from the standard spherical-bubble scenario; the paper identifies this as the most promising application but does not compute it.","In the holographic dual, the dissolving remnant should correspond to a shrinking black hole that eventually disappears into the confining geometry; the paper notes this is an open problem, and a bulk simulation would provide a sharp test.","Adding a small bulk or shear viscosity to the same model would test whether the heating stage persists; the paper explicitly works in the dissipation-free limit, so this is an extension, not a claim.","The virial-theorem explanation ties the reheating to the shape of the Landau potential, so a different potential shape could reverse the effect; the claimed robustness may therefore be limited to models fitted to the same domain-wall data."],"forward_implications":["In an expanding medium that passes through a first-order transition, bubble nucleation can be suppressed because trapped supercooled plasma is reheated to the transition temperature by latent heat from the shrinking boundary.","The remnant size shrinks as $L(\\tau)\\sim 1/\\tau$ in boost-invariant flow and $L(t)\\sim 1/a(t)^3$ in an FRW universe, giving clean scaling predictions for numerical and possibly experimental searches.","During the final dissolution stage, the residual deconfined plasma heats up above $T_c$ while $1-\\gamma\\sim 1/\\sqrt{\\tau}$, a counterintuitive signature that distinguishes this mechanism from ordinary supercooled droplet collapse.","Because the same stages appear in flat and cosmological geometries, the effect is expected in any slowly expanding medium with a first-order transition and negligible dissipation."],"supporting_citations":[{"why":"Supplies the extended-hydrodynamics model with a dynamical mixing fraction that the paper uses for all simulations.","marker":"[1]"},{"why":"Defines the holographic confinement model whose phase structure the effective model is fitted to.","marker":"[2]"},{"why":"Provides the numerical domain-wall data in the holographic model used to fit the effective potential.","marker":"[3]"},{"why":"Gives the adiabatic-hydrodynamics formalism that generates the temperature coupling and the conserved entropy current.","marker":"[13]"},{"why":"Shows that the dissipation-free model reproduces holographic domain-wall velocities, justifying the omission of friction.","marker":"[16]"},{"why":"Introduces the Bjorken boost-invariant expansion used as the flat-space test setup.","marker":"[17]"},{"why":"Provides the integrated conservation equation used for the remnant-size estimates in the FRW case.","marker":"[26]"}],"fun_headline_variants":["Hot remnants ignore expansion cooling, then heat up","Plasma remnants stay hot as space expands, then surge","Expansion can't cool hot remnants that later reheat","Hot remnants shrink, don't cool, then heat up further","Plasma remnants resist cooling, then reheat on dissolution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the model has no dissipation or friction, so the entropy current is exactly conserved and the remnant's behavior is derived from that conservation; if real plasma viscosity is significant, the remnants could cool, nucleate bubbles, or dissolve away without reheating.","fun_headline_variants_meta":{"raw":{"variants":["Hot remnants ignore expansion cooling, then heat up","Plasma remnants stay hot as space expands, then surge","Expansion can't cool hot remnants that later reheat","Hot remnants shrink, don't cool, then heat up further","Plasma remnants resist cooling, then reheat on dissolution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000575,"raw_usage":{"total_tokens":2644,"prompt_tokens":806,"completion_tokens":1838,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":422,"completion_tokens_details":{"reasoning_tokens":1758}},"tokens_in":422,"tokens_out":1838,"duration_ms":12950,"temperature":1.0,"reasoning_tokens":1758,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:31:55.911125+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the boost-invariant simulation with the dissipative or friction terms that the paper deliberately omits, using the same initial conditions: if a blob of deconfined matter near $T_c$ then cools below $T_c$, nucleates bubbles, or dissolves without reheating above $T_c$, the claim of generic hot remnants is falsified; a weaker test is to measure the entropy production rate, since the scaling laws require exact conservation.","supporting_citations":[{"cited_title":"high energy","cited_arxiv_id":null,"evidence_quote":"Provides the numerical domain-wall data in the holographic model used to fit the effective potential."}],"review_version":1}