{"id":"bd85e448-79b5-44ad-bf12-ebe1325db1c9","arxiv_id":"2505.04477","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"Hot quark stars built from a two-flavor holographic QCD equation of state reach masses of 2 to 17 solar masses and could mimic black holes.","lead":"The authors use a holographic model of two-flavor quark matter to build theoretical hot quark stars with masses from 2 to 17 times the Sun's. These stars might look like black holes from afar and could fill the 'mass gap' between neutron stars and black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass range rests on treating five fixed-(mu,T) fits as a barotropic stellar EoS; absent a physical temperature profile or proof the curves bracket the true EoS, the 2-17 Msun claim is unsupported.","rationale":"The reader flagged the barotropic extrapolation and the unproven bracketing by five curves as the weakest assumption; my analysis agrees. The central calculation is a straightforward TOV integration, but its physical interpretation as a 'hot' quark star requires a barotropic EoS. The five fits are arbitrary slices of the (mu,T) plane, and the paper provides no argument that a real stellar trajectory, with varying temperature, lies inside their envelope. The CEP coordinate swap is a separate concrete inconsistency that undermines the 'around the critical point' framing. Because the reader already issued CONDITIONAL, this stress-test does not shift the verdict; it sharpens the condition under which the paper should be accepted: the authors should provide either an actual temperature-profile calculation or a quantitative bracketing argument.","tokens_in":15205,"tokens_out":8047,"duration_ms":78729,"concrete_test":"Use the (mu,T)-dependent p and epsilon data behind Figure 7 to construct physically motivated stellar trajectories (e.g., isentropic s/n_B = const, or the T(r) profile from Eq. (4.2) for a chosen Rc), derive the corresponding barotropic EoS, solve the TOV equations for the same central pressures as in Section 3, and compare the M-R curves to Figures 2-6. If the true trajectory leaves the band of the five fits or yields masses outside 2-17 Msun, the bracketing assumption fails. As a minimal check, recompute the 'simple' models with the quark EoS terminated at its lowest fitted pressure p_min and matched to a hadron crust, removing the p-to-zero extrapolation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that Eqs. (2.8)-(2.12), each fitted from a fixed slice of the (mu,T) plane near the claimed CEP, can be integrated in the TOV equations as a single barotropic EoS from the center down to p=0. For a star with central temperatures of order 182 MeV, the pressure depends on both T and mu; no temperature profile is solved in Section 3, and the five fits are simply asserted in Section 5 to 'enclose the real EoS'. Section 4 shows the physical EoS for a star with a thermal profile is T0-dependent, so the five curves are arbitrary candidates, not a proven bracket. The extrapolation below the minimum fitted pressure p_min is also unjustified, since the quark-gluon plasma has nonzero pressure at the phase transition. An internal inconsistency in the CEP coordinates (Introduction assigns (555,105) to the 2-flavor CEP, while Section 2 uses (219,182)) further weakens the claim that these slices are centered on the critical point. Unless one of these issues is resolved, the mass range 2-17 Msun is an artifact of the chosen slices rather than a prediction for hot quark stars.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies a holographic Einstein-Maxwell-dilaton model for two-flavor QCD, with parameters matched to lattice data, to extract equations of state for hot quark-gluon plasma near the claimed critical endpoint. Five EoS curves are fitted for fixed values of mu/T, T, or mu (Eqs. 2.8-2.12), and these are used as barotropic inputs to the TOV equations. For the 'simple' models, the fits are extended to zero pressure; for the 'combined' models, the quark core is matched to a polytropic hadron shell through an energy-density jump parameter m. The resulting static configurations have masses from about 2 to 17 solar masses, radii around 100-275 km, and maximum compactness near 0.22, which the authors present as black-hole-mimicker and mass-gap candidates. The paper also analyzes I-Love-Q-C universal relations, presents two-dimensional (mu,T) heatmaps of energy density and pressure, and discusses a constant-thermal-conductivity temperature profile with an illustrative stellar solution.","tokens_in":15570,"tokens_out":7812,"duration_ms":77166,"significance":"If the EoS treatment were justified, the paper would offer a concrete holographic prediction for hot quark-star masses and a possible interpretation of mass-gap and black-hole-mimicker candidates. The full (mu,T) parameter maps in Section 4 could be a useful resource for future dense-matter studies, and the I-Love-Q-C tests extend universal-relation studies to two-layer models with phase transitions. The authors are honest about several limitations, such as the quark pressure being nonzero at the transition and the simplified nature of the analysis. However, the central astrophysical claim currently rests on unverified assumptions about how five fixed-slice fits bound the true two-dimensional EoS, so the significance is not yet established.","major_comments":[{"comment":"The central mass range of 2-17 solar masses is obtained by treating five fixed-(mu,T) fits as barotropic EoSs, but the actual energy density and pressure depend on both mu and T. The statement in Section 2 that 'the real EoS will be contained within this range' is an assertion, not a derivation: one-dimensional slices through a two-dimensional surface need not bracket the surface, and Section 4 shows explicitly that the stellar EoS depends on the central temperature T0. Please either prove the bracketing property, for example by demonstrating monotonic behavior of epsilon(p) over the enclosed (mu,T) region, or solve the TOV equations with the full two-dimensional EoS and a self-consistent thermal profile.","section":"Section 2-3, Eqs. (2.8)-(2.12)"},{"comment":"The fits are extrapolated to p=0 although the fitting ranges in Eqs. (2.8)-(2.12) begin at p_min around 2-3 x 10^-8 in solar units and the authors acknowledge that the quark-gluon plasma has nonzero pressure at the phase transition. The 'simple' models therefore place the stellar surface in a regime where the fitted EoS has no data and where the quark phase should already have ended. The claimed minimum mass near 2 solar masses and the low-mass parts of the mass-radius curves in Figure 2 depend on this extrapolation; please remove it or show quantitatively that the results are insensitive to the lower pressure cutoff.","section":"Section 3, after Eqs. (2.8)-(2.12)"},{"comment":"The location of the two-flavor critical endpoint is stated inconsistently. Section 2 states that the CEP of two-flavor QCD is at (mu,T)=(219,182) MeV, while the Introduction assigns (555,105) to the two-flavor case and (219,182) to the 2+1-flavor case. Since the five fitted slices are selected 'around the CEP,' this is not a cosmetic issue; please correct the CEP coordinates and state explicitly which reference and which flavor number give the CEP used for the fits in Figures 1-6.","section":"Introduction and Section 2"},{"comment":"The combined-case predictions depend on the ad hoc energy-density jump parameter m and on the choice of the hadronic polytrope epsilon_n = kappa_n p^{0.5}. No physical or observational motivation is given for selecting m in {1, 1.2, 0.8, 0.5}, and the resulting radii differ by more than a factor of two across these choices. Since the abstract quotes a single mass range spanning the combined models, please report the mass and radius ranges separately for each m and discuss how m could be constrained by microphysics or observations.","section":"Section 3, Eq. (3.2), Figures 3-6"},{"comment":"The paper does not provide a quantitative stability analysis for the TOV solutions. The text cites reference [52] for stability judging methods and mentions 'unstable components' for larger phase transitions, but no radial-oscillation modes or turning-point analysis is shown. Without such an analysis, it is not established that the computed sequences correspond to stable stars rather than unstable equilibrium branches, which is especially important given the very high masses and large radii far outside the usual neutron-star range.","section":"Section 3, Figures 2-6"}],"minor_comments":[{"comment":"The blue data points, especially at low pressure, are difficult to discern, and no fit residuals or uncertainties are reported for Eqs. (2.8)-(2.12); adding these would allow the reader to assess the quality of the EoS fits.","section":"Figure 1"},{"comment":"The Laplace equation for temperature assumes a flat, source-free region, while the star is relativistic and the text assumes a central heat source. Please state the approximations under which Eq. (4.2) is taken as the temperature profile.","section":"Section 4, Eq. (4.1)"},{"comment":"The statement that substituting T(r) into p(mu,T) yields mu(p,r) assumes a unique inversion; please specify the numerical method used and how any multi-valued branches are selected.","section":"Section 4"},{"comment":"The quantities I, Love, Q, and C are not defined in the text or figure captions; please define them or cite the conventions of reference [38] so that the plots are self-contained.","section":"Figures 2-6"},{"comment":"The Figure 5 caption contains the duplicated phrase 'the the I-Love-Q-C relations,' and reference [50] is dated 2023 in the arXiv listing but labeled 2025 in the reference list; please correct these typographical errors.","section":"Figure 5 caption and References"}],"recommendation":"major_revision","confidential_remarks":"To the editor: This manuscript is a follow-up of the authors' previous 2+1-flavor work, and the central astrophysical prediction currently rests on an asserted bracketing of the true two-dimensional EoS by five one-dimensional slices. I recommend major revision rather than rejection because the holographic framework and the full (mu,T) heatmaps could support a corrected calculation if the authors either prove the bracket or integrate the full EoS with a self-consistent temperature profile. The internal inconsistency in the CEP coordinates between the Introduction and Section 2 should be fixed before the paper can be considered reliable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: the paper knows its model, but the headline mass range is not a prediction. Equations (2.8)–(2.12) are fits to five fixed slices through the (μ,T) plane near the claimed critical endpoint, and the whole stellar calculation treats each as a single barotropic EoS integrated down to p=0. There is no temperature profile for the 'hot' star—Section 3 uses an EoS whose temperature is fixed along each curve, then Section 4 admits that with a realistic thermal profile each star needs its own EoS. The assertion that the five curves 'enclose the real EoS' is just that: an assertion. So the 2–17 M⊙ envelope is the envelope of the chosen slices, not a bracketing of physically possible stars.\n\nWhat is genuinely useful: the paper produces new two-flavor EMD-derived EoS fits from a model matched to lattice data, gives full ε(μ,T) and p(μ,T) heatmaps, computes M–R and I–Love–Q–C curves for simple and two-layer stars, and adds a constant-conductivity toy model. The TOV and tidal calculations are standard; nothing looks botched. The I–Love–Q universality surviving across their models is a mild, plausible finding.\n\nThe soft spots are real and load-bearing for the astrophysical claims. First, the CEP coordinates: the Introduction assigns (555,105) to 2-flavor and (219,182) to 2+1-flavor, then Section 2 says [27] obtained (219,182) for 2-flavor. That is an internal contradiction in the central anchor of the calculation. Second, the fits are stated without residuals or error bars, and they are extended to p=0 even though the authors note the quark pressure is nonzero at the phase transition; extending to zero is a modeling choice that needs physical justification. Third, the hadron shell is a single polytrope with γ=0.5 and κ chosen to cap NS mass at 2.2 M⊙, and the m parameter in Eq. (3.2) is scanned over 0.5–1.2 without a criterion. None of these are fatal to the paper as an exploratory study, but they are fatal to the abstract's claim that these are BH-mimicker candidates in the mass gap.\n\nWho gets value: someone working on holographic QCD phenomenology or exotic compact-object signatures will want this as a reference for the 2-flavor EoS and heatmaps. I would not cite the mass-range claim in an astrophysics paper until the thermal-structure issue is addressed.\n\nMy recommendation: send it to peer review, but with a major-revision bar. The authors should resolve the CEP inconsistency, provide fit residuals and a clear statement about the p=0 extension, and either solve a temperature profile or label the 2–17 M⊙ range more cautiously as an illustration.","headline":"The 2–17 M⊙ mass range is an artifact of treating five ad hoc (μ,T) slices as a barotropic EoS; the paper is still a competent and useful extension of the authors' earlier holographic quark-star work.","tokens_in":16081,"tokens_out":2721,"would_cite":false,"duration_ms":24662,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.40.Dg","12.38.Mh"],"model":"deepseek-v4-flash","headline":"The holographic two-flavor QCD equation of state yields hot quark stars of 2–17 solar masses that can mimic black holes and fill the mass gap.","keywords":["quark stars","holographic QCD","Einstein-Maxwell-Dilaton model","equation of state","I-Love-Q relations","Tolman-Oppenheimer-Volkoff equations","black hole mimickers","QCD critical endpoint"],"falsifier":"Compute the holographic equation of state along the actual stellar adiabat from center to surface, with temperature declining outward, and integrate the TOV equations; if the resulting maximum mass falls below about $2\\,M_\\odot$ or the mass–radius band no longer reaches $17\\,M_\\odot$, the central claim fails. A second decisive check would be a gravitational-wave detection of a mass-gap compact object whose tidal Love number is measured to be consistent with zero, which would remove the observational need for these stars as black-hole mimickers.","tokens_in":14933,"feed_emoji":"🌟","tokens_out":8653,"duration_ms":75881,"temperature":0.7,"pith_summary":"This paper aims to establish that the holographic Einstein–Maxwell-dilaton model, with parameters fixed by two-flavor lattice QCD, can supply equations of state for quark–gluon plasma near the critical endpoint at $T=182$ MeV that behave like stellar matter. Integrated through the Tolman–Oppenheimer–Volkoff equations, those equations of state produce hot quark stars with masses from about 2 to 17 solar masses and a maximum compactness around 0.22. A sympathetic reader would care because such stars would sit in the black-hole mass gap and mimic stellar-mass black holes while remaining distinguishable by nonzero tidal Love numbers. The paper also claims that the I–Love–Q relations continue to hold when hadron shells are added, as long as the energy-density jump at the quark–hadron interface is not too large.","feed_headline":"Hot quark stars from 2 to 17 solar masses mimic black holes","feed_subtitle":"A two-flavor QCD equation of state yields mass-gap candidates that gravitational waves could distinguish from black holes.","key_machinery":"The central object is the Einstein–Maxwell-dilaton (EMD) action in five dimensions, a holographic model of strongly coupled QCD in which the scalar potential and gauge coupling are fixed by lattice data. Holographic renormalization of this action yields $\\epsilon$ and $p$ as functions of temperature and chemical potential. The load-bearing mechanism for the stellar claim is the set of five fitted barotropic relations $\\epsilon_i(p_i)$ in Eqs. (2.8)–(2.12), each a two-term power law, which are extrapolated to zero pressure and fed into the Tolman–Oppenheimer–Volkoff equations. The I–Love–Q–C relations are computed with two-layer stellar procedures to test whether the resulting stars obey the neutron-star universal relations.","core_discovery":"The authors claim that the two-flavor version of the holographic Einstein–Maxwell-dilaton model, with parameters matched to two-flavor lattice QCD, yields equations of state for quark–gluon plasma near the critical endpoint at $T=182$ MeV that can serve as stellar cores. Fitting five curves around the critical point (Eqs. (2.8)–(2.12)) and integrating the Tolman–Oppenheimer–Volkoff equations, they obtain static hot quark stars with masses from about 2 to 17 $M_\\odot$ and maximum compactness around 0.22, which they call black-hole mimickers and mass-gap candidates. Adding a hadron shell with an energy-density jump factor $m$ preserves the I–Love–Q relations for $m=1$, $1.2$, and $0.8$, while the $m=0.5$ case shows larger deviations in compactness-related relations. The paper also presents the full $\\epsilon(\\mu,T)$ and $p(\\mu,T)$ parameter maps for the quark phase and a constant-thermal-conductivity stellar model with a central heat source, in which the equation of state depends on the central temperature.","pith_inferences":["A natural next step, not carried out in the paper, is to compute the full temperature and chemical-potential profiles inside the star and evaluate the equation of state along that actual adiabat, which would test whether the five fitted curves really bracket the physical stellar equation of state.","The paper assumes the quark–hadron interface strongly reflects radiation to keep the star hot; a quantitative calculation of the interface reflectivity and resulting cooling time would turn the black-hole-mimicker scenario into a testable lifetime prediction.","The constant-thermal-conductivity model with a central heat source is treated in a simplified way; computing thermal conductivity from the same holographic model would replace the assumed Laplace profile and could change the mass–radius predictions.","If the mass-gap claim survives a full equation-of-state calculation, the same EMD framework could be used to predict the merger signatures of two such stars, extending the paper's static analysis to dynamical gravitational-wave signals."],"forward_implications":["Mass-gap events in the range $2.5$–$5\\,M_\\odot$ could be hot quark stars rather than black holes, so gravitational-wave catalogues should treat that mass interval as possibly populated by exotic stars.","Because black holes have zero tidal Love number in general relativity, a measured nonzero tidal deformability for a mass-gap compact object would support the quark-star interpretation.","The I–Love–Q universality extends to quark cores with hadron shells when the interface jump is small, so the relations can be used to extract quark–hadron transition parameters from observations.","The predicted minimum mass near $2\\,M_\\odot$ means hot quark stars occupy a separate mass band above ordinary neutron stars, consistent with quark matter being disfavored inside canonical neutron stars.","The strong flavor dependence of the predicted masses and radii offers a way to infer whether a hot quark star contains only up and down quarks or also strange quarks."],"supporting_citations":[{"why":"Supplies the two-flavor holographic model, its lattice-fixed parameters, and the critical endpoint location used to select the five fitting curves.","marker":"[27]"},{"why":"Provides the holographic renormalization and thermodynamic relations used to extract energy density and pressure from the EMD action.","marker":"[23]"},{"why":"Provides the two-flavor lattice equation-of-state data used to fix the model parameters.","marker":"[28]"},{"why":"Provides lattice results at low chemical potential used to validate the holographic model's baryon-density predictions.","marker":"[29]"},{"why":"Supplies the I-Love-Q calculation procedure and the polytropic neutron-star comparison model used in the universal-relation analysis.","marker":"[38]"},{"why":"Supplies the two-layer stellar procedure used to join quark cores to hadron shells in the 'combined' models.","marker":"[51]"},{"why":"The authors' earlier 2+1-flavor hot quark star work that this paper's mass and radius results are compared with.","marker":"[21]"},{"why":"Provides the gamma_n = 0.5 polytropic neutron-star model whose parameters are used for the hadron shells.","marker":"[50]"}],"fun_headline_variants":["Hot quark stars: black hole mimickers from 2 to 17 solar masses","Mass-gap candidates: hot quark stars from holographic QCD","Holographic hot quark stars mimic black holes across mass gap","Hot quark stars up to 17 solar masses could be black hole mimickers","Black hole mimickers: hot quark stars from 2 to 17 solar masses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that equations of state fitted from high-temperature quark–gluon plasma data around the critical point can be extrapolated down to zero pressure and used in the TOV equations as a barotropic stellar model, without solving the internal temperature structure and without proof that the five chosen curves bracket the true equation of state.","fun_headline_variants_meta":{"raw":{"variants":["Hot quark stars: black hole mimickers from 2 to 17 solar masses","Mass-gap candidates: hot quark stars from holographic QCD","Holographic hot quark stars mimic black holes across mass gap","Hot quark stars up to 17 solar masses could be black hole mimickers","Black hole mimickers: hot quark stars from 2 to 17 solar masses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000945,"raw_usage":{"total_tokens":4037,"prompt_tokens":947,"completion_tokens":3090,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":2992}},"tokens_in":563,"tokens_out":3090,"duration_ms":22281,"temperature":1.0,"reasoning_tokens":2992,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:27:53.638126+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the holographic equation of state along the actual stellar adiabat from center to surface, with temperature declining outward, and integrate the TOV equations; if the resulting maximum mass falls below about $2\\,M_\\odot$ or the mass–radius band no longer reaches $17\\,M_\\odot$, the central claim fails. A second decisive check would be a gravitational-wave detection of a mass-gap compact object whose tidal Love number is measured to be consistent with zero, which would remove the observational need for these stars as black-hole mimickers.","supporting_citations":[{"cited_title":"Phase structure and critical phenomena in two-flavor QCD by holography.Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the two-flavor holographic model, its lattice-fixed parameters, and the critical endpoint location used to select the five fitting curves."},{"cited_title":"Probing QCD critical point and induced gravitational wave by black hole physics.Phys","cited_arxiv_id":null,"evidence_quote":"Provides the holographic renormalization and thermodynamic relations used to extract energy density and pressure from the EMD action."},{"cited_title":"Equation of state of quark-gluon matter from lattice QCD with two flavors of twisted mass Wilson fermions","cited_arxiv_id":null,"evidence_quote":"Provides the two-flavor lattice equation-of-state data used to fix the model parameters."},{"cited_title":"Quark number susceptibilities and equation of state at finite chemical potential in staggered QCD with Nt = 8.Phys","cited_arxiv_id":null,"evidence_quote":"Provides lattice results at low chemical potential used to validate the holographic model's baryon-density predictions."},{"cited_title":"I-love-q relations in neutron stars and their applications to astrophysics, gravitational waves, and fundamental physics.Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the I-Love-Q calculation procedure and the polytropic neutron-star comparison model used in the universal-relation analysis."},{"cited_title":"GW170817 and GW190425 as hybrid stars of dark and nuclear matter.Eur","cited_arxiv_id":null,"evidence_quote":"Supplies the two-layer stellar procedure used to join quark cores to hadron shells in the 'combined' models."},{"cited_title":"Dark I-Love-Q","cited_arxiv_id":"2309.07971","evidence_quote":"Provides the gamma_n = 0.5 polytropic neutron-star model whose parameters are used for the hadron shells."}],"review_version":1}