{"id":"97bd4747-ae86-4d84-bdef-8f8c732f10c1","arxiv_id":"2608.13557","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For a large class of holographic confining gauge theories far from conformality, the maximum supercooling equals half the squared speed of sound at the critical temperature, up to 1/D corrections.","lead":"In holographic models of strongly coupled gauge theories, the authors derive a general relation between the maximum possible supercooling of the confinement transition and the speed of sound in the deconfined phase. If correct, the result means supercooling is generically small away from conformality, which matters for predicting gravitational wave signals from cosmological phase transitions.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central formula depends on an unproven deep-IR localization condition in Appendix B that is not required by the 1/D expansion and is not tested by the exponential examples; if violated, Eq. (3.19) and ε_sc = c_s^2(T_c)/2 no longer follow.","rationale":"The paper's near-horizon 1/D expansion and the algebra leading from Eq. (3.13) to Eqs. (3.22)–(3.24) are internally consistent. Given the localized free-energy integral, the derivation of ε_sc = c_s^2(T_c)/2 follows: the conditions at φ_min fix T_min, the Taylor expansion around φ_min gives T_c − T_min ∝ δ/D^2, and the same expansion in d log T/d log s gives c_s^2(T_c) = 2δ/D^2. The numerical check in Fig. 4 supports the relation for the exponential superpotential family, although no code or error bars are provided. The single load-bearing weak point is the localization of the free-energy integral, exactly as the reader identified. Eq. (3.19) and the resulting φ_c − φ_min = −√(2D)/D are only valid if the deep-IR contribution I2 in Appendix B is exponentially suppressed. The bound in Eq. (B.9) makes the required conditions explicit: |dA/dφ| must have a positive infimum and dT/dφ must not grow exponentially in the IR. These are asserted, not derived from the slow-roll conditions or from the Gubser good-singularity criterion that selects the allowed potentials in §4. The examples in §4 all have W → e^{γφ} at large φ, so A(φ) is asymptotically linear and the bound is automatically satisfied; they therefore cannot test the assumption. A counterexample with A~−c log φ would shift the dominant contribution away from φ_min and invalidate the universal formula, while a proof that good singularities force the required IR behavior would close the gap. Because the concern is concrete, testable, and explicitly acknowledged in Appendix B, the appropriate disposition remains CONDITIONAL; the present read does not change the reader's verdict.","tokens_in":17439,"tokens_out":22158,"duration_ms":213248,"concrete_test":"Numerically construct the full black-brane solution for a one-parameter family of superpotentials that satisfies the near-T_min conditions (3.13) but interpolates from the exponential IR asymptotics W∼e^{γφ} to a faster-than-exponential asymptotics W∼e^{aφ^2} in the deep IR, for fixed D (say D = 5 and D = 10). For each member, compute T_min, T_c, and c_s^2(T_c), and evaluate I2 in Eq. (B.1) directly from the numerical A(φ) and T(φ). If the branch where I2 becomes comparable to I1 also shows ε_sc deviating from c_s^2(T_c)/2, then the Appendix B localization assumption is load-bearing and must be stated as an explicit condition; if every good-singularity member of the family keeps I2 exponentially suppressed, the gap is closed and the formula survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (3.19), which fixes φ_c via φ_c − φ_min = −√(2D)/D and hence yields ε_sc = δ/D^2 and ε_sc = c_s^2(T_c)/2, is obtained by evaluating the free-energy integral (3.18) as if it were dominated by a window of width 1/D around φ_min. Appendix B splits the integral at φ_⋆ and bounds the deep-IR piece I2 by Eq. (B.9): |I2| ≤ e^{(D−1)A(φ_min)} e^{−(D−1)(φ_⋆−φ_min)/√(2D)} (1/((D−1)c)) |dT/dφ(φ_0)|, with c = inf_{φ>φ_⋆}|dA/dφ|. This bound is exponentially small only if c is not exponentially small and dT/dφ does not grow exponentially. The paper states these as very mild conditions on the IR asymptotics and explicitly assumes away extra extrema of T(φ) (we assume a minimal case...), but does not derive them from the slow-roll conditions (3.1)–(3.2) or from the Gubser good-singularity criterion used for the examples. Superpotentials with faster-than-exponential IR growth (e.g. W∼e^{aφ^2}) give A∼−const×log φ, so c = 0 and I2 is not suppressed; T(φ) can also develop additional extrema in the deep IR. The numerical checks in §4 all use W = 1 + e^{γφ}, which has linear A(φ) at large φ, so they do not probe this assumption. If the IR conditions fail, φ_c is no longer forced to lie within 1/D of φ_min and the universal relation ε_sc = c_s^2/2 is not a consequence of the 1/D expansion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the thermal confinement/deconfinement transition in Einstein-scalar holography by generalising the bulk to D+1 dimensions and working at leading order in the 1/D expansion. It constructs the black brane solution analytically in a boundary-layer approximation, derives the Hawking temperature as a function of the horizon scalar value, and computes the critical temperature by comparing free energies of the black brane and thermal dAdS geometries. The central claim is that the maximum possible supercooling, epsilon_sc = 1 - T_min/T_c, satisfies epsilon_sc = c_s^2(T_c)/2 at leading order, with epsilon_sc generically suppressed as 1/D^2 and independent of the detailed form of the scalar potential. The claim is checked against numerical solutions for the exponential superpotential W = 1 + e^{gamma phi}, against published improved-holography results, and against N=4 super Yang-Mills on a sphere.","tokens_in":17828,"tokens_out":18739,"duration_ms":161692,"significance":"If the central relation holds, it is a striking and potentially universal statement: the maximal supercooling of the deconfined phase is fixed by the thermodynamic speed of sound at T_c rather than by detailed features of the scalar potential. The paper's analytic construction of the localised black brane in the 1/D expansion is a useful technical contribution, and the numerical check in Section 4 for the exponential superpotential is a genuine, if model-specific, independent test. The result is also falsifiable: a counterexample in the stated class with epsilon_sc different from c_s^2(T_c)/2 at leading order would refute it. However, the universality claim is currently conditional on deep-IR assumptions that are stated but not proven, so the advertised breadth of the result goes somewhat beyond what is demonstrated.","major_comments":[{"comment":"The derivation of Eq. (3.19), and hence of Eqs. (3.22) and (3.24), relies on the bound |I2| <= e^{(D-1)A(phi_min)} e^{-(D-1)(phi_star-phi_min)/sqrt(2D)} (1/((D-1)c)) |dT/dphi(phi_0)|. This bound is exponentially small only if c = inf_{phi>phi_star}|dA/dphi| is not exponentially small and |dT/dphi(phi_0)| does not grow exponentially. These conditions are not derived from the slow-roll conditions (3.1)-(3.2) or from the Gubser criterion used in Section 4; the text asserts them as 'very mild conditions on the IR asymptotics' without proof. They can fail for plausible potentials, for example W ~ exp(a phi^2), which gives A ~ -const*log(phi) and hence c=0, and T(phi) can develop additional extrema in the deep IR. Because all numerical checks in Section 4 use W = 1 + e^{gamma phi}, for which A is asymptotically linear, the examples do not probe this assumption. Without a proof, or an explicit restriction of the theorem to potentials satisfying these IR conditions, Eq. (3.19) and the universal relation epsilon_sc = c_s^2(T_c)/2 are not established.","section":"Appendix B, Eqs. (B.8)-(B.9)"},{"comment":"The free-energy comparison assumes a two-branch structure with a unique minimum of T(phi_h) and no additional extrema of T in the deep IR; the text explicitly says 'we assume a minimal case in which such a situation is excluded'. This branch structure is an input to the calculation, not a consequence of the 1/D expansion. If a thermodynamically stable small-black-brane branch exists, or if T(phi_h) has extra extrema, the identification of T_min and the derivation of Eq. (3.21) break down. Please state these assumptions as part of the theorem and verify them explicitly in the examples used for comparison.","section":"Section 3.2 and Appendix B"}],"minor_comments":[{"comment":"The sentence 'the ratio delta/D can be replaced by c_s^2(T_c)/c_s,CFT^2' is off by a factor of 2 relative to the preceding definitions: from Eq. (3.23), c_s^2/c_s,CFT^2 = 2 delta (D-1)/D^2, so delta/D = (1/2) c_s^2/c_s,CFT^2 + O(1/D). The final relation (3.24) is unaffected, but this sentence should be corrected.","section":"Section 3.3, after Eq. (3.23)"},{"comment":"The footnote correctly notes that the sound speed at a first-order transition is formally ill-defined. Please state more explicitly that the relation applies to the smooth continuation of the deconfined-phase equation of state, and define how far below T_c this continuation is being used; as written, a reader may interpret c_s^2(T_c) as the physical sound speed in the mixed-phase region.","section":"Section 3.3, footnote 3"},{"comment":"No numerical details or error estimates are given for the black brane solutions shown in Fig. 4. A brief description of the numerical method, grid resolution, convergence with D, and estimated uncertainties would make the comparison quantitative. The note about ChatGPT in the caption is out of place in a scientific paper.","section":"Section 4, Fig. 4"},{"comment":"The claim of agreement with N=4 SYM should be framed as a consistency check rather than as an independent test of the 1/D expansion, because Eq. (4.4) follows from the same thermodynamic definitions c_s^2 = dlogT/dlogs and T_min/T_c. Relatedly, the abstract's phrase 'independent of the details of the scalar potential' overstates the scope; the result depends on the existence of a minimum of T(phi_h) and on the IR localization assumptions. A phrasing such as 'independent of the form of the potential within the stated class' would be more accurate.","section":"Abstract and Section 4, Eq. (4.4)"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially interesting and the 1/D construction is a genuine technical step forward. My main concern is the unproved deep-IR localization condition in Appendix B; this is the load-bearing point for the claimed universality. The numerical exponential-superpotential check is the strongest independent evidence, but two of the comparison cases in Section 4 (N=4 SYM and the lattice/SU(N) estimates) draw on the authors' own earlier work or on algebraic identities. Please ensure the final version clearly distinguishes between a proof under explicit assumptions and a conjecture for a broader class of potentials."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this paper derives, at leading order in a 1/D expansion, a relation between maximum supercooling and sound speed in Einstein-scalar holographic models away from conformality, and it gives an explicit analytic black brane solution in that regime. The result is genuinely new: nobody had connected epsilon_sc to c_s^2(T_c)/2 before, and the derivation is clean and internally consistent. The paper is worth taking seriously.\n\nWhat it does well: the 1/D machinery is applied carefully. The black brane localization is standard in spirit but the authors go beyond the pure AdS case to a large class of superpotentials, and the leading-order matching is correct. The relation epsilon_sc = c_s^2(T_c)/2 follows from a short chain of steps: minimal temperature from the superpotential conditions (3.13), free energy localized near phi_min (3.19), and the sound speed expansion (3.23). The comparison with the exponential superpotential family in Sec. 4 shows real numerical agreement, and the caveat that the approximation fails near conformality is stated honestly. The paper also correctly flags that the sound speed at a first-order transition is formally ill-defined; it treats c_s as the smooth continuation, which is the right thing to do.\n\nSoft spots, in order of importance:\n\n1. The localization assumption in Appendix B is the load-bearing step and it is not proven. Equation (B.9) is exponentially suppressed only if the IR warp factor decays at least linearly and dT/dphi does not grow exponentially. The authors call these 'very mild conditions' and explicitly assume away extra extrema of T(phi), but they do not derive these conditions from the slow-roll conditions (3.1)-(3.2) or from the Gubser criterion used for the examples. The exponential superpotential W = 1 + e^{gamma phi}, which is the only numerical test in Sec. 4, has linear A(phi) at large phi, so it cannot probe this assumption. A superpotential like W ~ e^{a phi^2} would give A ~ -const * log phi, making c = 0 and the bound useless. This is a real gap, but it is also addressable: the authors could either prove the conditions for a general class of good singularities or find a counterexample. The central relation should be stated as conditional on this IR behavior.\n\n2. The N=4 SYM comparison is a stretch. The paper correctly notes that on a sphere there is no standard sound speed, so it defines c_s via dlog T / dlog s. That is fine as an analogy, but calling it 'striking agreement' with (4.4) is generous—it is a match of a specific functional form, not a universal check. This is a minor overstatement, not a fatal one.\n\n3. The numerics in Sec. 4 lack code and error bars. The reader can reproduce the qualitative agreement from the figure, but a published claim of this strength should have a reproducible script or at least a statement of the numerical method and accuracy. The acknowledgment that one figure was generated by ChatGPT is amusing but not a substitute for reproducibility.\n\n4. The lattice comparison is brief and indirect; the paper uses a lattice sound speed from SU(3) YM to suggest sub-percent supercooling. This is suggestive but not a check of the relation itself.\n\nOn the citation pattern: the self-citations [69,70] are used as confirmatory examples of small supercooling in large-N YM. That is appropriate because those papers are directly on point. No sign of citation padding.\n\nWho this is for: anyone working on holographic models of confinement, cosmological phase transitions, or gravitational wave phenomenology from strongly coupled sectors. It gives a simple, testable relation that can be checked on the lattice and in other holographic models.\n\nReader's verdict: the reader said CONDITIONAL, and I agree. The central derivation holds up under the stated assumptions; the assumptions are not fully justified, but they are plausible and the failure mode is explicit. The paper deserves a serious referee. It should not be desk-rejected.","headline":"A clean 1/D derivation of a new universal supercooling relation in Einstein-scalar gravity, with a load-bearing IR assumption that the paper states but does not prove.","tokens_in":18372,"tokens_out":2005,"would_cite":true,"duration_ms":16878,"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":"For a wide class of holographic gauge theories far from conformality, the maximum supercooling of the confinement transition is fixed by the speed of sound at the critical temperature: $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$.","keywords":["holographic confinement transition","1/D expansion","supercooling","Einstein-scalar gravity","Hawking-Page transition","speed of sound","black brane thermodynamics","large D limit"],"falsifier":"Construct a holographic model satisfying the stated slow-variation conditions but with an IR warp-factor plateau, or with a thermodynamically stable small-black-brane branch, and compute $\\epsilon_{\\rm sc}$ and $c_s^2(T_c)$ numerically: a deviation from $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$ would falsify the claim. Alternatively, on the lattice, measure both $T_{\\min}/T_c$ and the sound speed just above $T_c$ in SU(3) Yang-Mills and check whether $1 - T_{\\min}/T_c$ equals $c_s^2(T_c)/2$.","tokens_in":17244,"feed_emoji":"🕳️","tokens_out":15524,"duration_ms":124003,"temperature":0.7,"pith_summary":"This paper tries to establish a universal relation for the thermal deconfinement transition in a wide class of holographic gauge theories far from conformality: the maximum possible supercooling of the metastable deconfined phase is set, at leading order, by the thermodynamic speed of sound at the critical temperature, $\\epsilon_{\\rm sc} = 1 - T_{\\rm min}/T_c = c_s^2(T_c)/2$. The derivation works in $D+1$-dimensional Einstein-scalar gravity, treating the number of spacetime dimensions $D$ as large. At leading order in $1/D$ the black-brane geometry is analytically tractable because its effects are localized near the horizon, and the result is independent of the detailed shape of the scalar potential. If true, this replaces a model-by-model numerical problem with a single model-independent prediction that can be checked, for example, on the lattice.","feed_headline":"Holographic supercooling equals half the sound speed squared","feed_subtitle":"In the 1/D expansion, supercooling is fixed by the sound speed, not by the scalar potential.","key_machinery":"The central object is the superpotential $W(\\phi)$, the single function whose derivatives encode the scalar potential in Einstein-scalar gravity, together with the $1/D$ boundary-layer structure it controls. In $D+1$ spacetime dimensions the blackening factor $f(u)$ changes from $0$ at the horizon to $1$ over a distance $L_{\\rm eff}/D$, much shorter than the local AdS curvature scale $L_{\\rm eff} = |\\dot A(u_h)|^{-1}$; this separation lets the paper match a near-horizon layer to the far-region superpotential flow equations $\\dot A = -W$ and $\\dot\\phi = 2(D-1)W'$. The thermodynamic engine is the Hawking temperature formula $T = (D/4\\pi) e^{A(u_h)}|\\dot A(u_h)|$, whose minimum defines $T_{\\rm min}$ through the two conditions $(W'/W)^2 = 1/(2D)$ and $W''/W = (1+\\delta)/(2D)$. The free-energy integral localizes near $\\phi_{\\min}$ because of the factor $e^{(D-1)A}$, and a Taylor expansion around $\\phi_{\\min}$ produces the master identity $\\epsilon_{\\rm sc} = \\delta/D^2 = c_s^2(T_c)/2$.","core_discovery":"The central claim is that, for any Einstein-scalar holographic model whose superpotential obeys the slow-variation conditions $(W'/W)^2 \\ll D/[2(D-1)]$ and $W''/W \\ll D/[2(D-1)]$ near the transition, the black brane dual to the deconfined phase can be constructed analytically in a $1/D$ expansion. The blackening function equals $1 - \\exp(D(u-u_h)/L_{\\rm eff})$, so horizon effects are confined to a layer of size $L_{\\rm eff}/D$. The Hawking temperature as a function of the horizon scalar $\\phi_h$ then has a minimum $T_{\\rm min}$, and the free-energy difference between the black-brane and deformed-AdS phases vanishes at $T_c = T_{\\rm min}(1 + \\delta/D^2)$. The same combination $\\delta/D^2$ appears in the deconfined sound speed, $c_s^2(T_c) = 2\\delta/D^2$, giving the identity $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$, with corrections of order $1/D$ and $(c_s/c_{s,\\rm CFT})^2$. The paper verifies this identity numerically for an exponential superpotential and reports agreement with improved holography and with $\\mathcal{N}=4$ super Yang-Mills on a sphere.","pith_inferences":["A testable extension would be to compute both $T_{\\rm min}/T_c$ and $c_s^2(T_c)$ independently in a non-exponential superpotential that satisfies the slow-variation conditions; the paper's logic predicts that the two determinations satisfy $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$ at leading order.","If the identity survives beyond the leading $1/D$ approximation, the ratio $(c_s/c_{s,\\rm CFT})^2$ becomes the natural control parameter for organizing corrections to the transition thermodynamics, much as slow-roll parameters organize inflationary predictions.","For cosmological hidden sectors, the relation implies that strongly coupled theories far from conformality reheat almost immediately to the critical temperature, which would weaken the gravitational-wave signal expected from such confinement transitions compared with estimates that allow large supercooling.","The analogy with slow-roll parameters suggests that the minimal-temperature conditions define a surface in superpotential space; scanning superpotentials on this surface would map out the allowed region where the formula is valid and where it breaks down."],"forward_implications":["A lattice measurement of the sound speed just above $T_c$ in SU(3) Yang-Mills, where $c_s^2(T_c) \\simeq 0.013$, translates through the identity into sub-percent maximum supercooling.","The maximum supercooling is generically suppressed as $1/D^2$, so in $D=4$ models far from conformality the deconfined phase can cool only a little below $T_c$ before the transition completes.","The result is independent of the detailed form of the scalar potential, so the same relation should hold across exponential-superpotential, improved-holography, and similar holographic constructions within the stated validity window.","The formula does not apply to near-conformal theories, where the expansion parameter $(c_s/c_{s,\\rm CFT})^2$ is of order one and supercooling can be large; the paper explicitly excludes that regime.","The $1/D^2$ suppression in the examples, together with lattice large-$N$ results, is presented as evidence that small supercooling may be a general property of strongly coupled gauge theories far from conformality."],"supporting_citations":[{"why":"Supplies the relation between the horizon value of the scalar and the temperature, and the free-energy branch structure used to compute $T_c$.","marker":"[31]"},{"why":"Introduces the superpotential flow equations that determine the vacuum dAdS geometry matched onto the near-horizon layer.","marker":"[50–53]"},{"why":"Provides the large-$D$ expansion of general relativity that makes the black-brane solution analytically tractable.","marker":"[55–58]"},{"why":"Reviews the large-$D$ limit and the boundary-layer methods used to localize the horizon.","marker":"[59]"},{"why":"Establishes the thermal AdS/black-hole phase transition as the holographic model of confinement and deconfinement.","marker":"[17]"},{"why":"Gives the exponential Chamblin-Reall superpotential solution used as an explicit test of the general formula.","marker":"[77]"},{"why":"Supplies the criterion for acceptable IR singularities that bounds the exponential-superpotential parameter range.","marker":"[78]"},{"why":"Provides the lattice sound speed in SU(3) Yang-Mills near $T_c$ used to infer sub-percent supercooling from the identity.","marker":"[81]"},{"why":"Provides an independent large-$N$ lattice prediction of small supercooling with which the $1/D^2$ scaling is compared.","marker":"[70]"}],"fun_headline_variants":["Supercooling pinned to sound speed in holography","1/D expansion reveals universal supercooling rule","Sound speed sets maximum holographic supercooling","Supercooling = half sound speed squared holographically","Holographic transition: supercooling tied to sound speed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the deep-infrared geometry is tame in a specific sense: the warp factor keeps decreasing steadily to minus infinity with a nonzero asymptotic slope, and the temperature as a function of the horizon scalar has no extra extrema, so the far-infrared part of the free-energy integral is exponentially suppressed; if this fails, $T_c$ is not within $1/D$ of $T_{\\min}$ and the relation $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$ is not established.","fun_headline_variants_meta":{"raw":{"variants":["Supercooling pinned to sound speed in holography","1/D expansion reveals universal supercooling rule","Sound speed sets maximum holographic supercooling","Supercooling = half sound speed squared holographically","Holographic transition: supercooling tied to sound speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":3026,"prompt_tokens":1085,"completion_tokens":1941,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":701,"completion_tokens_details":{"reasoning_tokens":1863}},"tokens_in":701,"tokens_out":1941,"duration_ms":14697,"temperature":1.0,"reasoning_tokens":1863,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:11:39.457135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct a holographic model satisfying the stated slow-variation conditions but with an IR warp-factor plateau, or with a thermodynamically stable small-black-brane branch, and compute $\\epsilon_{\\rm sc}$ and $c_s^2(T_c)$ numerically: a deviation from $\\epsilon_{\\rm sc} = c_s^2(T_c)/2$ would falsify the claim. Alternatively, on the lattice, measure both $T_{\\min}/T_c$ and the sound speed just above $T_c$ in SU(3) Yang-Mills and check whether $1 - T_{\\min}/T_c$ equals $c_s^2(T_c)/2$.","supporting_citations":[],"review_version":1}