{"id":"940b5db2-cf55-4988-b71a-5bb983ade111","arxiv_id":"2509.03947","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In the improved holographic EMD model, two hairy black hole phases are separated by a U-shaped boundary whose lower branch is first-order and upper branch third-order, meeting at (mu_B,T)=(765.51,86.54) MeV.","lead":"This paper numerically maps the phase diagram of a holographic model of dense quark matter and finds that the boundary between two types of 'hairy' black holes contains both a first-order and a subtle third-order transition line meeting at a critical point. It matters because it predicts richer structure in a widely used holographic QCD model, which may guide searches for the QCD critical point.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (37) has the wrong sign for dF at fixed mu_B and the free energy is never normalized; the Type-I/II transition-order and critical-point claims are therefore unestablished.","rationale":"I read the paper as a numerical extension of the improved EMD model; the existence of two hairy branches and a first-order line roughly matching [54,55] is plausible and independently supported by prior literature. The genuinely new claims are the third-order line and the critical point at the turning point of the boundary. Those claims are not supported by the thermodynamic argument as written. The reader's weakest assumption, the missing free-energy integration constant, is correct; I add that Eq. (37) also has a sign inconsistency with the paper's own definition F=-P. Since the first-order transition is determined by equality of absolute free energies and the third-order classification by derivatives across the two branches, both depend on a normalization that is never given. This is a load-bearing gap, not a disagreement with the model or with holographic EMD conventions. The requested check is concrete and could validate or refute the critical point. My recommendation remains conditional: the paper should not be accepted until the normalization and sign are addressed and numerical convergence/error estimates are supplied.","tokens_in":16848,"tokens_out":8163,"duration_ms":85785,"concrete_test":"Recompute the Type-I and Type-II free energies from the Euclidean on-shell action of Eq. (1) with holographic renormalization, or at minimum state F(T,mu)=F_0(mu)-integral_{T0}^T s(T',mu)dT' with F_0(mu) fixed by dF/dmu=-rho at fixed T and by the mu=0 lattice-matched limit. Then redraw Figs. 6-9 and re-extract the transition temperatures and orders. If the crossing temperatures or the location of the d^2 s/dT^2 discontinuity shift by more than the reported precision (e.g., 86.54 MeV), the critical-point claim is not robust. Also verify the sign by checking that, at fixed mu, F decreases with T when F=-P.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The thermodynamic classification rests entirely on Eq. (37), dF(T, fixed mu_B) = s dT, used to construct F(T) for each hairy branch and to locate first-order crossings (Fig. 6) and third-order discontinuities (Fig. 9). Two defects undermine it. First, with F = -P (Eq. 36) the Gibbs-Duhem relation gives dP = s dT + rho dmu, hence dF = -s dT at fixed mu; the sign in Eq. (37) is wrong. If the plotted F curves were produced by integrating +s dT, they are not the grand potential, and the swallowtail argument is not connected to the stated thermodynamics. Second, even after the sign is corrected, dF = -s dT fixes F only up to a mu-dependent integration constant. The paper never specifies the holographically renormalized on-shell action, a reference normalization, or any condition fixing F_0(mu); Eq. (36) is only a definition. The first-order transition temperature is determined by equality of the two branches' free energies, so without this constant the claimed transition temperatures, the third-order classification (which requires continuity of F and s and discontinuity of d^2 s/dT^2 across the same branches), and the critical point (765.51, 86.54) MeV are not determined by the text. Numerical derivatives without error bars or convergence data further impede resolving a 'subtle' third-order jump.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the improved holographic Einstein-Maxwell-Dilaton model calibrated to lattice QCD at zero baryon chemical potential (Refs. [54,55]) and classifies static hairy black hole solutions by the sign of the near-horizon scalar derivative psi1_h in Eq. (21). It maps a U-shaped phase boundary in the (mu_B,T) plane between Type-I and Type-II hairy phases. The paper claims that the lower branch of this boundary is a first-order phase transition line consistent with earlier work, the upper branch is a subtle third-order line, and the two branches meet at a critical point (mu_B^crit,T^crit)=(765.51,86.54) MeV that coincides with the turning point of the boundary curve.","tokens_in":17230,"tokens_out":6367,"duration_ms":68906,"significance":"The systematic numerical scan and the gravitational classification of two distinct hairy black hole solutions are useful contributions, and the consistency of the first-order segment with previous EMD results is a valuable cross-check. If the claimed third-order line and critical point were firmly established, they would be a nontrivial addition to holographic black hole thermodynamics and could inform model building. However, the thermodynamic derivation on which these new claims rest is incomplete, so the central novelty is not currently supported.","major_comments":[{"comment":"The free energy whose derivatives define the transition order is never properly constructed. With F = -P in Eq. (36), the Gibbs-Duhem relation dP = s dT + rho dmu gives dF = -s dT at fixed mu, not +s dT as printed. If the F(T) curves in Fig. 6 were generated by integrating Eq. (37) literally, they are not the grand potential. Moreover, the differential fixes F only up to a mu-dependent integration constant; the paper does not specify this constant, a reference normalization, or the holographically renormalized on-shell action. Since the first-order transition temperatures are obtained from equality of F branches and the third-order classification requires continuity of F and s across branches, both are sensitive to this unspecified construction. Thus the first-order line, the third-order line, and the critical point (765.51,86.54) MeV are not derived by the text as it stands.","section":"Sec. IV.B, Eq. (37)"},{"comment":"The third-order claim rests on resolving a 'subtle' discontinuity in the second temperature derivative of the entropy density. No error bars, radial grid convergence tests, or derivative-stencil details are provided. The plotted curves in Fig. 9 appear noisy, e.g. the mu_B=900 panel shows d^2s/dT^2 jumping from about -50 to +50 over roughly 1 MeV, and the mu_B=950 panel shows sharp spikes. Without a convergence study or a quantitative estimate of the jump and its uncertainty, a finite discontinuity cannot be distinguished from numerical noise. Please provide such evidence.","section":"Sec. IV.B, Figs. 8-9"},{"comment":"The phase boundary is defined geometrically by the vanishing of the near-horizon coefficient psi1_h, and the U-shaped curve in Fig. 5 is then identified as the thermodynamic phase boundary. However, in equilibrium the coexistence curve should be determined by equality of the relevant thermodynamic potentials, not by the sign of a near-horizon expansion coefficient. The lower branch is checked against previous first-order results, but the upper branch has no independent thermodynamic determination. This reinforces the need for a properly defined, normalized free energy before the phase diagram can be accepted.","section":"Sec. III.B, Eq. (21) and Sec. IV.A, Fig. 5"}],"minor_comments":[{"comment":"The notation A'_t(r) is used, but no A_t field is defined in the ansatz; the gauge field is denoted phi(r). Please clarify or correct the typo.","section":"Eq. (16)"},{"comment":"The notation beta_+-prime is confusing; define explicitly that the prime denotes differentiation with respect to psi0_h, and write the derivatives as partial derivatives for clarity.","section":"Eqs. (22)-(24)"},{"comment":"The integration is started at r_start=10^-8 and truncated at r_end=10, but no convergence test with respect to these cutoffs is reported. A brief check (e.g. varying r_end by a factor of two) would strengthen the numerical claims.","section":"Sec. III.A"},{"comment":"The horizontal axis label 'Phi_cutoff_1 / Phi_max_1' is ambiguous; please spell out the ratio and explain how the cutoff value was chosen.","section":"Fig. 1, left panel"},{"comment":"In the introductory sentence, 'thorough Gibbs conditions' should be 'through Gibbs conditions.'","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a topic within the journal's scope, and the gravitational-side analysis is potentially interesting. However, the central thermodynamic classification and the critical-point claim require a full derivation of the free energy from the on-shell action, including renormalization and the integration constant, and a convergence study for the reported third-order discontinuity. The sign error in Eq. (37) is concrete and should be corrected before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nHere's my read on Guo et al. The genuinely new thing is the claim that the lattice-calibrated improved EMD model has two distinct hairy black hole branches whose phase boundary contains a first-order segment and a subtle third-order segment meeting at a critical point that sits exactly at the turning point of the boundary curve. The first-order line is consistent with Critelli and Grefa, and the two-branch classification is a reasonable extension of the simplified-model analysis in their ref. [62]. The numerical work looks standard: they solve the ODE system with horizon shooting, extract the holographic coefficients, and map a large region of parameter space. That part is plausible and probably reproducible.\n\nThe soft spot is the thermodynamics. The paper defines the free energy as F = -P (Eq. 36), which is correct. Then Eq. (37) states dF = s dT at fixed mu. That is the wrong sign: from F = -P and Gibbs-Duhem, dP = s dT + rho dmu, you get dF = -s dT. If the plotted free energy curves were actually obtained by integrating +s dT, they are not the grand potential, and the swallowtail argument is not connected to the stated thermodynamics. Even with the sign fixed, dF = -s dT determines F only up to a mu-dependent integration constant. The paper never gives the holographically renormalized on-shell action, a reference normalization, or any condition that fixes F_0(mu). Without that, you cannot compare the free energies of the two branches, which is exactly what determines the first-order transition temperature and whether the solid boundary is a genuine transition at all. The third-order classification, which requires continuity of F and s and a discontinuous second derivative of s across the same branches, is therefore not established by the text. The numerical derivatives in Fig. 9 also have no error bars or convergence checks, and the claimed 'subtle' jump is hard to assess from the plots alone.\n\nI want to stress that the circularity burden is low: the model parameters were fitted to zero-density lattice QCD, not to the claimed phase boundary, so the result is not a fit to the target. And the bulk classification of Type I and Type II hairs is independent of the free energy issue, so part of the paper stands on its own.\n\nBottom line: this is a paper worth reading for the bulk solution taxonomy, but the central thermodynamic claim needs a careful fix. I'd send it to peer review -- a competent referee can sort out the sign error and ask for the missing normalization and data. As it stands, I wouldn't cite the critical point or third-order line.\n\nRecommendation: engage with it, but treat the thermodynamic results as provisional.","headline":"A credible numerical study of two hairy black hole branches, but the claimed third-order transition and critical point rest on a sign error and an unspecified free energy normalization.","tokens_in":17703,"tokens_out":4029,"would_cite":false,"duration_ms":39278,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In the improved holographic Einstein–Maxwell–Dilaton model, the phase boundary between the two hairy black hole types is a first-order line plus a third-order line, meeting at (765.51, 86.54) MeV.","keywords":["holographic QCD phase diagram","Einstein-Maxwell-Dilaton model","hairy black holes","third-order phase transition","first-order phase transition","critical endpoint","baryon chemical potential","gauge/gravity duality"],"falsifier":"Compute the free energy from an explicit holographic on-shell action, or fix the integration constant in dF = s dT by a stated normalization, and re-examine the upper phase boundary; if the discontinuity in the second temperature derivative of the entropy density vanishes or shifts away from (765.51, 86.54) MeV, the claimed third-order line and the critical point are artifacts of the unspecified free-energy construction.","tokens_in":16750,"feed_emoji":"🕳️","tokens_out":7606,"duration_ms":65282,"temperature":0.7,"pith_summary":"This paper tries to establish that the phase boundary between two kinds of hairy black holes in a lattice-calibrated holographic Einstein–Maxwell–Dilaton model has two distinct parts: a first-order transition line and a much weaker third-order transition line. The two lines meet at a critical point at baryon chemical potential 765.51 MeV and temperature 86.54 MeV, which is exactly the turning point of the boundary curve. The claim matters because holographic models of this kind are used to simulate the QCD phase diagram at high baryon density, where direct lattice QCD calculations fail. If the claim is right, the standard first-order line is only half the story: there is a second, mild transition on the same boundary, and the endpoint is fixed by the geometry of the boundary curve.","feed_headline":"At 86.54 MeV, two hairy black hole phases meet at a critical point","feed_subtitle":"The holographic QCD phase boundary is a first-order line plus a subtle third-order line, both ending at one point.","key_machinery":"The sign of the near-horizon expansion coefficient ψ₁ʰ in Eq. (21) decides whether the scalar hair decays monotonically (Type I) or first grows near the horizon (Type II); its vanishing marks the boundary between the two hairy phases. Thermodynamically, the transition order is read from the free energy F(T) at fixed μ_B through dF = s dT and the Ehrenfest classification: a finite jump in entropy gives the first-order line, while a finite discontinuity in ∂²s/∂T² gives the third-order line. The U-shaped boundary curve's turning point is where the two branches meet.","core_discovery":"The paper reports two distinct hairy black hole solutions in the improved holographic EMD model. Type-I hairy black holes are governed by the scalar potential, with a scalar profile that decays monotonically away from the horizon; Type-II solutions are governed by the nonminimal coupling to the U(1) gauge field, with a scalar profile that first grows near the horizon and then decays. In the (μ_B, T) plane the boundary between the two phases forms a U-shaped curve. The lower segment coincides with the first-order transition previously identified in holographic QCD models; the upper segment, analyzed through the Gibbs free energy and its temperature derivatives, shows a finite discontinuity on","pith_inferences":["Because the free-energy normalization is left unspecified, the classification of the upper branch as third-order is only as strong as that choice; an explicit on-shell action could confirm or move the line.","The same black-hole-physics analysis should apply to other lattice-calibrated EMD models; the coincidence of the critical point with the turning point may be a generic feature of the calibrated parameter set, not a structural necessity.","A practical way to search for this weak transition in other holographic models is to look at the nonmonotonicity of the scalar hair near the horizon (the Type II barrier) rather than at boundary free-energy kinks.","Tuning f(ψ) and V(ψ) should move the critical point along the boundary; this suggests a design rule for engineering phase diagrams with a desired first-order line length and a third-order endpoint."],"forward_implications":["For any μ_B above 765.51 MeV, raising the temperature first crosses a sharp first-order transition and then a mild third-order transition on the same phase boundary.","The upper third-order line is subtle: free energy, entropy, and ∂s/∂T are continuous, so this transition would be missed by conventional holographic thermodynamic probes.","The gravitational origin of the two phases is tied to the competition between V′(ψ) and f′(ψ) in the near-horizon expansion, giving a bulk criterion for which boundary phase dominates.","Because the critical point coincides with the turning point of the phase boundary, the endpoint can be read directly from the geometry of the boundary curve rather than from a separate thermodynamic calculation.","Targeted modifications of the scalar potential and gauge coupling should allow deliberate engineering of the phase diagram's first-order line and third-order endpoint."],"supporting_citations":[{"why":"Introduced the five-dimensional EMD model and its holographic critical point, the gravitational setup this work builds on.","marker":"[44]"},{"why":"Extended the EMD construction to critical dynamics, justifying the holographic dictionary used here.","marker":"[45]"},{"why":"Calibrated the improved EMD model to lattice QCD and located the first-order transition line that the paper's lower boundary reproduces.","marker":"[54]"},{"why":"Supplied the parameter values for f(ψ), V(ψ), and κ₅², plus the equation-of-state matching, that define the model solved here.","marker":"[55]"},{"why":"Provided the rescaling of numerical bulk data to physical μ_B and T used to map the phase diagram.","marker":"[56]"},{"why":"The authors' earlier simplified EMD study that introduced the Type I/Type II hairy black hole classification and first identified a third-order transition beyond the critical point.","marker":"[62]"}],"fun_headline_variants":["Critical point where Type-I and Type-II hairy black holes meet","Hairy black hole phases share a critical point on μ_B-T diagram","Improved EMD model: two hairy phases end at one critical point","First-order and third-order lines converge at hairy black hole criticality","At 86.54 MeV, hairy black hole phase boundary hits a critical point"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The classification assumes the free energy density is uniquely determined by its temperature differential dF = s dT at fixed chemical potential, but the paper does not specify how the integration constant is fixed.","fun_headline_variants_meta":{"raw":{"variants":["Critical point where Type-I and Type-II hairy black holes meet","Hairy black hole phases share a critical point on μ_B-T diagram","Improved EMD model: two hairy phases end at one critical point","First-order and third-order lines converge at hairy black hole criticality","At 86.54 MeV, hairy black hole phase boundary hits a critical point"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1360,"prompt_tokens":750,"completion_tokens":610,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":516}},"tokens_in":494,"tokens_out":610,"duration_ms":6499,"temperature":1.0,"reasoning_tokens":516,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:31:25.196485+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the free energy from an explicit holographic on-shell action, or fix the integration constant in dF = s dT by a stated normalization, and re-examine the upper phase boundary; if the discontinuity in the second temperature derivative of the entropy density vanishes or shifts away from (765.51, 86.54) MeV, the claimed third-order line and the critical point are artifacts of the unspecified free-energy construction.","supporting_citations":[{"cited_title":"The underlying black hole phase transitions in an Einstein-Maxwell-dilaton model with a holographic critical point","cited_arxiv_id":"2410.05065","evidence_quote":"The authors' earlier simplified EMD study that introduced the Type I/Type II hairy black hole classification and first identified a third-order transition beyond the critical point."}],"review_version":1}