{"id":"68fa1020-0e99-4dc6-b9f5-3299e9ae51b4","arxiv_id":"1908.10415","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Cold, beta-equilibrated charm quark matter from perturbative QCD yields charm star configurations that are dynamically unstable under radial oscillations.","lead":"This paper extends a perturbative QCD equation of state to include charm quarks and tests whether 'charm stars' could exist. It finds that such stars are dynamically unstable under radial oscillations, so they would not be realized in nature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The EoS is built by summing two independent one-massive-flavor systems (Eq. 4), which is not a valid NNLO Nf=2+1+1 thermodynamic potential; the claimed instability may be an artifact of this construction.","rationale":"The reader identified the pairwise flavor summation as an unresolved technical point but treated the pQCD convergence and X>=3 restriction as the weakest assumption. I view the EoS construction as more load-bearing: it is the foundation of the entire stability analysis, and there are concrete algebraic reasons to doubt it. The paper's own statement that Eq. (4) is 'a convenient way of writing the degrees of freedom' is not a derivation; the NNLO thermodynamic potential contains gluonic and ring contributions that are functions of the total flavor content and cannot be decomposed into independent pair sums. No validation against the known Nf=4 limit or a direct two-massive-flavor calculation is provided. If the test reveals that Eq. (4) is incorrect, the central claim of charm-star instability could be an artifact of the flawed EoS; if the test shows the pairwise sum is exact (which is unlikely), the claim stands. Because the paper transparently restricts to the X>=3 band satisfying the Bodmer-Witten hypothesis and reproduces known pulsation frequencies, I do not reject it outright; rather, this concern justifies keeping the conditional verdict and making validation of Eq. (4) a necessary condition for acceptance. This is the same verdict as the reader, so no adjustment is needed.","tokens_in":11463,"tokens_out":14565,"duration_ms":134532,"concrete_test":"Set m_s = m_c = 0 in Eq. (4) for the (u+c)+(d+s) split and compare the resulting pressure to the standard NNLO Nf=4 massless pQCD pressure (Ghisoiu et al., Nucl. Phys. B915 (2017) 102, or the O(alpha_s^2) result from Ref. [18]). The pairwise sum predicts an N_f^2 coefficient of 1^2+1^2=2, whereas the correct Nf=4 coefficient is 4; any disagreement at that order settles whether the pairwise construction is valid. If the pairwise sum fails this limit, the EoS behind the instability claim is not the NNLO pQCD EoS for Nf=2+1+1 and the central result is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, that charm quark stars are unstable under radial oscillations, depends entirely on the equation of state constructed in Sec. III. In Eq. (4), the Nf=2+1+1 potential is written as a sum of two independent (N_l=1)+1 systems, (u+c) and (d+s). Each term in the sum contains the full massless one-flavor thermodynamic potential, including gluonic and Debye-screening contributions. These contributions are not additive in flavor number: the ring term is nonlinear in the total quark density, so computing it separately for each pair and adding the results omits cross-flavor screening and double-counts purely gluonic NNLO terms (e.g., terms whose flavor dependence is N_f^2). The paper provides no proof that such an additivity holds, and the massless limit m_s, m_c -> 0 does not reproduce the known Nf=4 NNLO pressure unless the missing cross terms are restored. Since the EoS determines the TOV mass-radius curve and the radial oscillation eigenfrequencies of the Gondek et al. method, an error here can shift the charm threshold, the maximum-mass configuration, and the sign of Im(f0). The instability conclusion is therefore not established by the calculation as presented.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the perturbative QCD equation-of-state formalism for N_f = N_l + 1 flavors to systems with multiple massive quarks by writing the thermodynamic potential as a sum of independent (N_l + 1) blocks. It applies this construction to beta-equilibrated, electrically neutral charm quark matter, builds the EoS, solves the TOV equations, and studies radial stability using the two first-order equations of Gondek et al. The central claim is that charm quark stars are dynamically unstable under radial oscillations across the renormalization-scale band X >= 3.","tokens_in":11732,"tokens_out":10366,"duration_ms":110608,"significance":"If the result holds, it would provide a first-principles pQCD update to the older MIT-bag-model conclusion that charm stars are unstable, with a falsifiable prediction for a possible new branch of ultra-dense hybrid stars. The paper is careful in using external lattice/PDG inputs, scanning the renormalization scale instead of fitting it, and validating the pulsation code against known EoSs from Ref. [54]. However, the central conclusion rests entirely on an EoS construction in Eq. (4) that is not the genuine NNLO QCD thermodynamic potential for N_f = 2 + 1 + 1. Because the missing cross-flavor and pure-gluon terms can change the EoS, the stability analysis, and the final instability claim, the result is not established by the calculation as presented.","major_comments":[{"comment":"The central thermodynamic potential is written as a sum of independent (N_l + 1) blocks, e.g. (u+c) + (d+s) for N_f = 2 + 1 + 1. This is not the NNLO QCD thermodynamic potential for a system with four active quark flavors. The ring/Debye-screening contribution is a nonlinear function of the total Debye mass built from all active flavors; summing separate one-massive-flavor results therefore omits cross-flavor screening between the two pairs and double-counts pure-gluon NNLO terms. A concrete diagnostic is the massless limit: as m_s, m_c -> 0, Eq. (4) does not reduce to the known massless N_f = 4 NNLO pressure. Since this EoS is the input to the TOV integration and to the Gondek et al. pulsation equations, the instability conclusion in Sec. IV is not supported by the calculation as presented. The authors need either to use a genuine N_f = 2 + 1 + 1 NNLO potential, or to demonstrate quantitatively that the omitted cross terms are negligible in the range mu_s ~ 1.2-1.4 GeV used for the charm threshold.","section":"Sec. II.B, Eq. (4); Sec. III.C"},{"comment":"The pressure is obtained by integrating number densities from a lower limit mu0(X), described in the text as the point of zero pressure. The value of mu0(X) is never given, nor is the criterion for choosing it. The resulting EoS P(epsilon), and hence the mass-radius curves and pulsation eigenfrequencies, depend on this integration constant. Please specify mu0(X) and demonstrate that the instability result is not an artifact of this choice.","section":"Sec. III.C, Eqs. (19)-(22)"},{"comment":"Dynamical instability under radial oscillations is displayed only for X = 3. The statement that the same behavior was obtained for larger values of X is not accompanied by any quantitative evidence. Since the conclusion is stated for the whole band X >= 3, please provide the Im(f0) versus epsilon_c curves for X = 4 and X = 5, or a table of fundamental-mode eigenvalues at the charm-threshold configurations.","section":"Sec. IV, Fig. 6"}],"minor_comments":[{"comment":"The general definition of the renormalization scale, Lambda_bar = X sum_i mu_i / N_f, is not used consistently in Eq. (17b), where the denominator is 3 rather than 4 for the N_f = 4 case. Please clarify whether this is intentional and explain the averaging prescription above the charm threshold.","section":"Sec. II.A and Sec. III.B, Eq. (17b)"},{"comment":"The restriction to X >= 3 is said to satisfy the Bodmer-Witten hypothesis as shown in Ref. [11], but the connection is not explained; a one-sentence statement of the criterion used would help the reader.","section":"Sec. IV"},{"comment":"The notation Omega[N_l^{(i)}] and Omega[1^{(i)}] is hard to follow because these are not defined in terms of the explicit functions of Ref. [11]. Writing the formula for N_f = 2 + 1 + 1 explicitly, or adding a short glossary, would improve readability.","section":"Sec. II.B, Eq. (4)"},{"comment":"The figure would benefit from a panel showing Re(f0) as well, since the text discusses Re(omega0) = 0 as the onset of instability and the reader cannot verify this condition from the current plot.","section":"Sec. IV, Fig. 6"},{"comment":"Reference [3] is cited as an arXiv preprint; if a peer-reviewed version exists, it should be cited instead.","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper addresses an interesting astrophysical question and the numerical implementation appears careful, but the EoS construction in Eq. (4) is not a controlled pQCD result for N_f = 2 + 1 + 1. The central instability claim depends on this construction, so the result is not established. A resubmission with a correct multi-flavor thermodynamic potential, or with the construction clearly labeled as a phenomenological model and its uncertainty quantified, could be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper asks a good question—are charm stars stable if you use a pQCD equation of state?—and it reaches an answer that agrees with old bag-model work, but the calculation as presented does not actually support that answer. The stress-test note is right, and it lands on the load-bearing part of the paper.\n\nWhat is new and good: the application of a pQCD EoS to charm quark matter under beta equilibrium and charge neutrality is new, and the authors are transparent about their assumptions. They restrict the renormalization-scale parameter to X >= 3 to satisfy the Bodmer-Witten hypothesis, show particle fractions and mass-radius curves, and use the Gondek method for radial oscillations. The qualitative conclusion, that charm stars are unstable, is consistent with earlier MIT bag-model analyses, which is reassuring. The citation pattern to the pQCD EoS literature is standard; I see no self-citation problem.\n\nThe problem is in Sec. II B, Eq. (4). The thermodynamic potential for Nf = 2+1+1 is written as a sum of two independent (Nl+1)-flavor blocks, (u+c) and (d+s). That is not the NNLO thermodynamic potential of QCD with four flavors. The ring term, which handles Debye screening, is nonlinear in the total quark content—it behaves like (sum over flavors)^{3/2}, not like a sum of independent blocks. Cross-flavor screening is missed, and purely gluonic NNLO terms are mis-weighted. In the massless limit, Eq. (4) does not reduce to the known Nf = 4 pressure. Since the EoS drives the TOV curves and the oscillation eigenfrequencies, the instability claim rests on an unvalidated construction. This is not a minor technical point; it affects the charm threshold and can change the sign of Im(f0).\n\nSecondary issues are smaller: the zero-pressure integration constant mu0 is not specified, no code or data are provided, and uncertainties from lattice inputs and alpha_s are not propagated. These are not fatal if the main construction were solid, but they add to the need for revision.\n\nBottom line: the paper deserves engagement because charm stars are a real question and the qualitative answer is likely correct. But as written, the central result is not established. If I were the editor, I would send it to peer review rather than desk-reject it, but I would not accept the instability claim without a corrected or explicitly justified EoS.","headline":"The paper's instability conclusion is plausible but not established: the equation of state is built from a pairwise flavor sum that is not the NNLO Nf=2+1+1 thermodynamic potential.","tokens_in":12259,"tokens_out":4775,"would_cite":false,"duration_ms":52342,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t","12.38.Bx","21.65.Qr","12.38.Mh","04.40.Dg"],"model":"deepseek-v4-flash","headline":"Charm stars, hypothetical ultradense objects with charm-quark cores, are dynamically unstable under a first-principles QCD equation of state.","keywords":["charm quark matter","perturbative QCD equation of state","charm stars","quark stars","radial stability","beta equilibrium","compact stars","renormalization scale"],"falsifier":"A stable compact object whose measured mass and radius fall on the charm-star branch, for example near $M \\approx 1.0\\,M_\\odot$ and $R \\approx 7.3$ km for $X = 3$ at the charm threshold, would contradict the instability claim. Alternatively, a recalculation of the same equation of state with next-to-next-to-next-to-leading-order corrections that yields a real fundamental pulsation frequency throughout the charm branch would also falsify it.","tokens_in":11261,"feed_emoji":"⭐","tokens_out":11450,"duration_ms":105071,"temperature":0.7,"pith_summary":"The paper asks whether quark matter that contains charm quarks can form stable compact stars. It extends a perturbative QCD equation of state, previously built for massless quarks plus one massive flavor, to include two massive flavors (strange and charm), and imposes beta equilibrium and electric charge neutrality. Using this equation of state and standard first-order radial pulsation equations, it finds that the stellar configurations with a charm core are dynamically unstable. A sympathetic reader would take this as a first-principles confirmation, within perturbative QCD, of the older bag-model conclusion that charm stars cannot exist.","feed_headline":"Charm stars are unstable, first-principles QCD shows","feed_subtitle":"First-principles QCD rules out a new ultradense branch of compact stars","key_machinery":"The central machinery is the extension of the perturbative QCD thermodynamic potential from $N_f = N_l + 1$ flavors to $N_f = N_l + N_m$ flavors by summing independent massless-plus-massive contributions, so that charm matter is treated as $(u+c)+(d+s)$. The pressure is reconstructed by integrating the number densities from the strangeness threshold, with $\\beta$-equilibrium and charge-neutrality conditions reducing all chemical potentials to functions of the strange chemical potential $\\mu_s$. The renormalization-scale parameter $X$, defined by $\\bar{\\Lambda}=X\\sum_i \\mu_i/N_f$, enters through the running coupling and running masses, and the requirement of a continuous passage through the charm threshold reduces the allowed band to $X > 4/3$; the stability analysis is restricted to $X \\geq 3$. Dynamical stability is tested with the first-order coupled equations from Ref. [25] for the relative radial displacement and the Lagrangian pressure perturbation, where the imaginary part of the fundamental eigenfrequency signals the onset of instability.","core_discovery":"The central claim is that cold, dense, electrically neutral quark matter in $\\beta$ equilibrium, described by a perturbative QCD equation of state that includes charm quarks up to next-to-next-to-leading order, does not admit stable charm stars. In the parameter band $X \\geq 3$, which satisfies the stability hypothesis for strange quark matter, the mass--central-energy-density curve has a second branch where charm quarks appear and the thermodynamic condition $\\partial M/\\partial \\epsilon_c \\geq 0$ holds, but the fundamental radial oscillation mode acquires a positive imaginary frequency there. Because the radial modes obey Sturm-Liouville ordering, once the fundamental mode is unstable all higher modes are unstable as well. The paper therefore concludes that bare charm stars are excluded as a new family of ultradense hybrid compact stars, while still allowing small charm fractions in the cores of heavy hybrid stars.","pith_inferences":["If the instability persists under higher-order corrections, the exclusion of the charm branch is generic to perturbative QCD equations of state; a natural next test is to repeat the radial analysis with published next-to-next-to-next-to-leading-order corrections.","The same heavy-flavor extension applied to bottom quarks would likely also produce unstable bare 'bottom stars', but the more tangible consequence is transient charm production in merger remnants where beta equilibrium is not yet established.","The pressure kink at the charm threshold could serve as a subtle signature: if the instability sets a maximum central density for quark cores, gravitational-wave constraints on the post-merger equation of state might indirectly reveal the charm threshold."],"forward_implications":["If the stability conclusion is right, no pure quark star with a charm core can be dynamically stable within this equation-of-state family, so searches for such objects should place the charm branch in excluded parameter regions.","The charm threshold softens the pressure and reduces the renormalization-scale uncertainty band to $X > 4/3$, making charm-aware quark matter equations of state more tightly constrained at high densities than their three-flavor counterparts.","Bare charm stars are ruled out, but the same framework leaves open a matching between nuclear matter and a quark phase with a small charm contamination in the cores of the heaviest neutron stars via a first-order transition.","A non-negligible charm fraction could contribute to the equation of state during the early stages of neutron-star mergers, when densities exceed the charm threshold."],"supporting_citations":[{"why":"Supplies the three-loop perturbative QCD thermodynamic potential for massless quarks plus one massive flavor that this paper extends to two massive flavors.","marker":"[11]"},{"why":"Provides the first-order radial pulsation equations whose solutions determine whether a stellar configuration is dynamically stable.","marker":"[25]"},{"why":"Gives the TOV structure equations, the thermodynamic stability condition, and the compact-star framework used here.","marker":"[21]"},{"why":"Establishes the earlier bag-model analysis of charm stars whose instability conclusion this paper revisits with perturbative QCD.","marker":"[20]"},{"why":"Supports the statement that higher-order perturbative and nonperturbative contributions modify the equation of state only mildly.","marker":"[19]"},{"why":"Provides the threshold matching conditions that motivate the continuity requirement on the renormalization scale.","marker":"[10]"},{"why":"Supplies the two-loop running quark mass formulas used for the strange and charm masses.","marker":"[35]"},{"why":"Provides the general-relativistic and causality exclusions used to delimit the allowed mass-radius region.","marker":"[45]"}],"fun_headline_variants":["Charm stars ruled out by perturbative QCD","Heavy quarks destabilize charm stars","First-principles QCD: no stable charm stars","Charm stars unstable, perturbative QCD says","No charm stars: QCD equation of state forbids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion depends on perturbative QCD being a reliable description at charm-threshold densities (strange chemical potential near 1.2 to 1.4 GeV, where the running coupling is still about 0.3) and on restricting the renormalization-scale parameter to the band $X \\geq 3$ that satisfies the stability hypothesis for quark matter.","fun_headline_variants_meta":{"raw":{"variants":["Charm stars ruled out by perturbative QCD","Heavy quarks destabilize charm stars","First-principles QCD: no stable charm stars","Charm stars unstable, perturbative QCD says","No charm stars: QCD equation of state forbids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000162,"raw_usage":{"total_tokens":1168,"prompt_tokens":801,"completion_tokens":367,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":417,"completion_tokens_details":{"reasoning_tokens":292}},"tokens_in":417,"tokens_out":367,"duration_ms":3789,"temperature":1.0,"reasoning_tokens":292,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:44:39.625746+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A stable compact object whose measured mass and radius fall on the charm-star branch, for example near $M \\approx 1.0\\,M_\\odot$ and $R \\approx 7.3$ km for $X = 3$ at the charm threshold, would contradict the instability claim. Alternatively, a recalculation of the same equation of state with next-to-next-to-next-to-leading-order corrections that yields a real fundamental pulsation frequency throughout the charm branch would also falsify it.","supporting_citations":[{"cited_title":"Gondek, P","cited_arxiv_id":null,"evidence_quote":"Provides the first-order radial pulsation equations whose solutions determine whether a stellar configuration is dynamically stable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the TOV structure equations, the thermodynamic stability condition, and the compact-star framework used here."},{"cited_title":"Kettner, F","cited_arxiv_id":null,"evidence_quote":"Establishes the earlier bag-model analysis of charm stars whose instability conclusion this paper revisits with perturbative QCD."},{"cited_title":"Gorda, A","cited_arxiv_id":null,"evidence_quote":"Supports the statement that higher-order perturbative and nonperturbative contributions modify the equation of state only mildly."},{"cited_title":"Rodrigo and A","cited_arxiv_id":null,"evidence_quote":"Provides the threshold matching conditions that motivate the continuity requirement on the renormalization scale."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-loop running quark mass formulas used for the strange and charm masses."},{"cited_title":"Haque, A","cited_arxiv_id":null,"evidence_quote":"Provides the general-relativistic and causality exclusions used to delimit the allowed mass-radius region."}],"review_version":1}