{"id":"819cdf75-902b-426e-9623-e015d9b62b0e","arxiv_id":"2412.06739","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Dark matter-admixed quarkyonic neutron stars show f-mode frequencies shifted by up to tens of percent, while universal frequency-density and frequency-compactness relations remain approximately valid.","lead":"This paper computes the vibration frequencies of neutron stars made of quarkyonic matter with dark matter mixed in, using two nuclear models and three free parameters. It finds that adding dark matter shifts the frequencies, yet the stars still follow known universal relations linking frequency to density and compactness.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universal-relations claim rests on a self-fit to the same 26 models; no out-of-sample or residual test rules out a systematic DM-induced deviation.","rationale":"The reader's weakest_assumption focuses on the DM abundance and couplings. Those affect the magnitude of the DM-induced frequency shifts, but even if the abundance is uncertain, the paper's central scientific claim is about universality: the abstract is not 'DM shifts f-modes' but 'universal relations hold despite DM.' The load-bearing step for that claim is the demonstration of the universal relation, not the precise value of kf^DM. The current demonstration is a self-fit with no validation, so it cannot exclude the possibility that DM introduces a systematic deviation. This is a correctness risk in the argument, not merely a missing robustness test: if the residual trend is real, the fitted Eq. (34) would overpredict or underpredict f for DM-admixed stars, and the conclusion 'demonstrating their robustness in the presence of dark matter' would be unsupported. The proposed test is cheap and decisive: the paper already has all the data needed to compute residuals against a non-DM fit. It would also address the reader's concern about uncertainties in the fitted constants. I therefore keep the reader's CONDITIONAL verdict, since the qualitative finding that DM changes f-mode frequencies is plausible and the Cowling approximation is standard, but the universal-relation robustness should be conditional on passing the out-of-sample residual test.","tokens_in":16650,"tokens_out":11154,"duration_ms":120485,"concrete_test":"Refit Eqs. (34) and (35) using only the 14 non-DM EOSs (baryonic and quarkyonic with kf^DM=0). For each of the 12 DM-admixed EOSs, compute residuals in f (fixed Mbar, Rbar) and in omega_M (fixed C) relative to these fits. The claim of robustness is supported only if the DM residual RMS is comparable to the non-DM residual RMS and shows no monotonic trend with kf^DM (0.03-0.04 GeV); if the residual RMS doubles or trends upward with DM fraction, the universal relation is not DM-robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that f-mode universal relations remain valid in the presence of dark matter, as stated in the abstract and quantified by Eqs. (34) and (35). The evidence is a linear fit to all 26 EOSs in Figs. 3 and 4, with correlation coefficients of 95% and 99%. This is not sufficient to establish universality. The same data are used both to define the fit and to test it, so the relation is guaranteed to have some degree of fit; no residual plot, scatter measure, or out-of-sample validation is provided. The concern is load-bearing because the paper emphasizes that DM causes 'notable variations' in f-mode frequencies. If the DM-admixed models systematically deviate from the relation defined by the baryonic and non-DM quarkyonic models, then the conclusion that the relation is robust to DM would be false, even though a combined fit can absorb the offset. The correlation of 95% for Eq. (34) is also a loose standard for a universal relation; asteroseismology applications typically require residuals at the few-percent level, not an R^2 of 0.90. Moreover, the plotted points are not independent: each EOS provides a whole M-R sequence, so the effective number of independent constraints is only 26, not the total number of plotted stars. This undercuts the statistical significance implied by the high correlation. A concrete check is to fit the relation on the non-DM subsample and examine the residuals of the DM-admixed models.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes l=2 f-mode oscillation frequencies of non-rotating neutron stars built from quarkyonic equations of state with admixed dark matter, using the relativistic Cowling approximation. The nuclear input combines an E-RMF baryonic sector, a quarkyonic matter model with two parameters (transition density n_t and QCD confinement scale Lambda_cs), and a Higgs-portal neutralino DM component characterized by the DM Fermi momentum k_f^DM. The authors integrate the TOV equations and solve the linearized perturbation equations for two parameter sets (G3 and IOPB-I), varying the three free parameters over a small grid. They report that adding DM lowers the maximum mass and increases f-mode frequencies for a given compactness, and they claim that universal relations among f-mode frequency, average density, and compactness remain valid, fitting Eqs. (34) and (35) with correlation coefficients of 95% and 99%. The paper also presents correlations among macroscopic properties such as mass, radius, tidal deformability, and f-mode frequency.","tokens_in":16933,"tokens_out":5514,"duration_ms":57147,"significance":"If the universal-relation claim is robust, this work would usefully extend f-mode universality to a class of exotic compact stars, which matters for planned gravitational-wave asteroseismology analyses. The manuscript has clear strengths: it uses standard, well-defined TOV and Cowling equations; it tabulates macroscopic properties for 26 EOSs in Table I; and it provides explicit, falsifiable fitting formulas in Eqs. (34)-(35). The qualitative statement that DM shifts f-mode frequencies is supported by the tabulated values. However, the central claim of EOS independence is not yet established, because the evidence is a self-fit to the same data used to define the relation, with no residual analysis, no out-of-sample test, and no comparison against the established Andersson-Kokkotas relation. The Cowling systematic error is acknowledged but not propagated. These are load-bearing gaps rather than presentation issues.","major_comments":[{"comment":"The claim that f-mode universal relations hold across baryonic, quarkyonic, and DM-admixed EOSs is supported only by linear least-squares fits to the same 26 EOSs used to define the fit. No out-of-sample validation, residual plot, or comparison with the Andersson-Kokkotas relation (cited in the Introduction) is provided. If the DM-admixed models systematically deviate from the relation defined by the non-DM models, the combined fit could absorb a DM-induced offset and masquerade as universality. Please fit Eq. (34) to the non-DM subsample and report residuals of the DM-admixed models, or provide an equivalent out-of-sample test.","section":"Sec. III.B, Eq. (34), Fig. 3"},{"comment":"The reported correlation coefficients are not a sufficient statistical basis for universality. For Eq. (34), r = 95% corresponds to R^2 ~ 0.90, implying roughly 10% unexplained scatter; asteroseismology applications typically require residuals at the few-percent level. Additionally, the plotted points are not independent: each EOS produces a continuous M-R sequence, so the effective number of independent constraints is 26 rather than the total number of plotted stars. Please report the RMS fractional residual around each fit and account for the non-independence of points when assessing the significance of the correlation.","section":"Sec. III.B, Figs. 3-4"},{"comment":"The text states that the Cowling approximation error is less than 20% relative to full general-relativistic results, but this systematic uncertainty is not propagated into the fitted universal relations. A frequency error of this size could be comparable to or larger than the scatter around Eqs. (34) and (35), and if the error is density-dependent it could bias the fitted slope. Please quantify the impact of the Cowling approximation on the universal relations, or add an explicit qualification that the reported correlations are within the Cowling approximation only.","section":"Sec. II.E and Sec. III.B"}],"minor_comments":[{"comment":"Typos: \"Mclerran\" should be \"McLerran\", and \"Jhao and Lattimer\" should be \"Zhao and Lattimer\".","section":"Introduction"},{"comment":"The transition Fermi momentum is introduced as k_Ft but appears as kt in Eq. (5); please define kt explicitly and use consistent notation.","section":"Sec. II.B, Eq. (5)"},{"comment":"The caption says \"the numbers show the corresponding p-values,\" but the displayed numbers are correlation coefficients, not p-values; please correct the caption.","section":"Fig. 7 caption"},{"comment":"There is a typo: \"dimesionless\" should be \"dimensionless\".","section":"Sec. III.B"},{"comment":"Typo: \"steller\" should be \"stellar\".","section":"Fig. 4 caption"},{"comment":"The phrase \"higher lower tidal deformability region\" is confusing; please rephrase to state whether the region has higher or lower tidal deformability.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a competent extension of the authors' earlier work on DM-admixed quarkyonic stars, and the f-mode calculation itself appears sound. The main weakness is that the universal-relation claim rests on an in-sample fit with no residual or out-of-sample checks; this is fixable within the scope of the paper. If the authors provide the requested residual analysis and compare with established relations, the paper could become acceptable. The novelty is modest (the DM model and quarkyonic EOS are taken from previous works, while the f-mode computation is new), but it may still be suitable for a specialized journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it actually does compute f-mode frequencies for dark matter–admixed quarkyonic neutron stars over a grid of transition density, confinement scale, and DM Fermi momentum, and that specific combination is new. Second, the central robustness claim—that universal relations survive the presence of DM—is weaker than the abstract implies, because the relations are fitted to and tested on the same 26 EOSs. No out-of-sample check or residual analysis is reported, so the paper does not actually demonstrate that DM-induced deviations are absent.\n\nWhat is genuinely useful: Table I gives a clean set of M_max, R_1.4, Lambda_1.4, and f-mode frequencies for 26 baryonic/quarkyonic/DM-admixed models with two nuclear parametrizations. The qualitative trend—DM raises f_1.4 while lowering M_max—is clear and physically plausible. The TOV and Cowling equations are standard, and the numerical machinery appears reliable. Self-citations to refs. [13,14] and [65] are appropriate here; the quarkyonic model and the Cowling approach come from those works.\n\nThe soft spots are real but not fatal. The universal-relations section (Eqs. 34-35) reports correlation coefficients of 95% and 99% from linear fits to all points combined. That is not evidence of EOS-independence. Each EOS contributes a whole sequence of stars, so the effective number of independent constraints is only 26, and the correlation is inflated. The right check is to fit on the non-DM subsample and look at residuals of DM-admixed models; the paper does not do this. The Cowling error (which the authors themselves say can reach 20%) is never propagated into the fits. The DM parameters, especially k_f^DM = 0.03 GeV from the M_chi/M_NS = 1/6 assumption, are somewhat ad hoc, though that is a modeling choice rather than an error.\n\nWho is this for? Researchers studying f-modes in exotic NS matter, or looking for DM signatures in future gravitational-wave detections. It is a niche but legitimate contribution. I would send it to peer review—the calculation is sound enough, and the parameter scan has value—but I would ask the authors to either perform an out-of-sample test of the universal relations or soften the claim to a statement that the relations are not strongly violated within their model. As it stands, the paper slightly oversells its main conclusion.\n\nRecommendation: engage with it, but in review push on the self-fit issue. It is a fixable weakness, not a fatal one.","headline":"Useful new f-mode numbers for DM-admixed quarkyonic stars, but the universal-relation claim is a self-fit, not a test.","tokens_in":17497,"tokens_out":2166,"would_cite":false,"duration_ms":25408,"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":"This paper claims that dark-matter-admixed quarkyonic neutron stars obey the same f-mode universal relations as ordinary neutron stars: frequency tracks mean density and compactness.","keywords":["f-mode oscillations","quarkyonic neutron stars","dark-matter admixed neutron stars","universal relations","Cowling approximation","relativistic mean field","neutron star equation of state"],"falsifier":"Recompute the f-mode frequencies for the same 26 equations of state with a fully general-relativistic oscillation code that includes metric perturbations; if the scatter around $f = 0.6282 + 2.0476\\sqrt{\\bar{M}/\\bar{R}^3}$ widens beyond the claimed 95% correlation, the universal relation is an artifact of the Cowling approximation. Observationally, a future post-merger gravitational-wave signal that yields an f-mode frequency for a neutron star with independently measured mass and radius would falsify the universal behavior if it falls off the fitted line by more than the scatter quoted here.","tokens_in":16464,"feed_emoji":"🌌","tokens_out":10543,"duration_ms":98581,"temperature":0.7,"pith_summary":"This paper asks whether the low-frequency f-mode oscillations of neutron stars built from quarkyonic matter—a mixed phase where quarks fill low-momentum states and nucleons occupy a shell near the Fermi surface—still follow the same universal scaling laws when dark matter is added. Using the Cowling approximation of linearized general relativity and two relativistic mean-field parameter sets, the authors vary three model parameters: quark–hadron transition density, QCD confinement scale, and dark-matter Fermi momentum. They find that adding dark matter shifts the f-mode frequency by up to a few tenths of a kHz and lowers the maximum mass. The paper's core claim is that the universal relations survive this shift: across all 26 baryonic, quarkyonic, and dark-matter-admixed equations of state, the frequency tracks the square root of mean density with 95% correlation and the scaled frequency tracks compactness with 99% correlation. If true, gravitational-wave measurements of f-modes from such exotic stars would remain interpretable with the same simple formulas used for ordinary neutron stars.","feed_headline":"Dark matter shifts neutron-star f-modes; universal law holds","feed_subtitle":"Quarkyonic stars with DM still obey one frequency-density law across 26 equations of state.","key_machinery":"The load-bearing machinery is the Cowling approximation for $\\ell=2$ f-modes: the spacetime metric is held fixed and only the fluid displacement is perturbed, which the paper notes reproduces fully relativistic frequencies to within about 20%. On top of this sits an effective relativistic mean-field (E-RMF) description of baryonic matter, a quarkyonic equation of state in which nucleons occupy a Fermi shell while quarks fill the low-momentum states, and a neutralino dark-matter component coupled through Higgs exchange with fixed couplings $y=0.07$, $f=0.35$ and a dark-matter Fermi momentum $k_f^{\\rm DM}$ taken as 0, 0.03, or 0.04 GeV. The three free parameters, transition density $n_t$, QCD confinement scale $\\Lambda_{\\rm cs}$, and $k_f^{\\rm DM}$, generate the 26 equations of state over which the two universal fits are drawn.","core_discovery":"The central claim, stated on the paper's own terms, is that f-mode oscillations of quarkyonic neutron stars, with and without a dark-matter admixture, obey the same approximate linear universal relations as ordinary neutron stars. From a grid of 26 equations of state spanning two nuclear parameter sets, three transition densities, two confinement scales, and three dark-matter Fermi momenta, the paper obtains $f\\,(\\mathrm{kHz}) = 0.6282 + 2.0476\\sqrt{\\bar{M}/\\bar{R}^3}$ with a correlation coefficient of about 95%, where $\\bar{M}$ and $\\bar{R}$ are mass and radius in units of $1.4\\,M_\\odot$ and $10$ km, and $\\omega M\\,(\\mathrm{kHz\\,km}) = -4.665 + 199.95\\,C$ with 99% correlation, where $C = M/R$ is compactness. Dark matter changes where a given star sits on these lines, raising the f-mode frequency at fixed microscopic parameters and lowering the maximum mass, but does not push stars off them. This is what the authors mean by saying the universal relations hold in the presence of dark matter.","pith_inferences":["A natural next step the paper does not take is to rerun the Cowling calculation in full general relativity; if the 95% and 99% correlations survive, the fits could serve as priors for gravitational-wave template banks targeting neutron stars with exotic cores.","Because the dark-matter sector enters only through the single parameter $k_f^{\\rm DM}$, the same calculation could be repeated for other dark-matter candidates, such as asymmetric or bosonic dark matter, without changing the formalism; whether universality holds would depend only on the resulting effective equation of state.","The universal fits are drawn over a deliberately sparse grid of parameter combinations; a Monte Carlo marginalization over $n_t$, $\\Lambda_{\\rm cs}$, and $k_f^{\\rm DM}$ would test whether the correlations reflect genuine equation-of-state independence or a coincidence of the chosen grid."],"forward_implications":["A detected f-mode frequency from a post-merger remnant can be converted into a model-independent estimate of mean density or compactness using either universal fit, regardless of whether the star contains quarkyonic matter or dark matter.","Because raising $k_f^{\\rm DM}$ increases $f_{1.4}$ while decreasing tidal deformability $\\Lambda_{1.4}$, a joint measurement of f-mode frequency and tidal deformability could, in principle, separate dark-matter content from transition-density effects in this model.","The 99%-correlation fit $\\omega M = -4.665 + 199.95 C$ ties scaled frequency to compactness tightly enough that a single f-mode detection would pin down the star's compactness to within the fit's scatter.","The fact that both the G3 and IOPB-I parameter sets fall on the same fits indicates that the universal relations are not an artifact of one specific nuclear parametrization."],"supporting_citations":[{"why":"Foundational quarkyonic matter model whose density dependence the paper adopts.","marker":"[85]"},{"why":"Extended quarkyonic model with beta equilibrium, charge neutrality, and chemical equilibrium that generates the equations of state used here.","marker":"[13]"},{"why":"Empirical f-mode universal relations that the paper tests and extends to dark-matter-admixed quarkyonic stars.","marker":"[62]"},{"why":"Original Cowling approximation, the oscillation scheme used for all f-mode calculations.","marker":"[63]"},{"why":"Neutralino-Higgs-portal dark-matter model with relativistic mean-field treatment used for the dark-matter component.","marker":"[40]"},{"why":"Empirical constraints used to fix the dark-matter couplings y=0.07 and f=0.35.","marker":"[23]"},{"why":"Earlier relativistic mean-field analysis of dark-matter influence on quarkyonic stars that supplies the equation-of-state inputs and parameter choices.","marker":"[14]"}],"fun_headline_variants":["DM shifts f-modes but universal law survives","Quarkyonic star f-modes: DM moves, law holds","Dark matter bends f-modes, not universal curve","F-modes shift with DM, universality intact","Neutron-star f-modes: DM changes, law stands"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dark-matter effect is sized by the assumption that the dark-matter Fermi momentum is roughly 0.03 GeV, obtained by taking the neutralino mass to be one sixth of the neutron-star mass and the nucleon number density to be 1000 times the dark-matter density; a smaller dark-matter fraction or weaker Higgs-portal couplings would shrink the frequency shifts and could remove the apparent universality.","fun_headline_variants_meta":{"raw":{"variants":["DM shifts f-modes but universal law survives","Quarkyonic star f-modes: DM moves, law holds","Dark matter bends f-modes, not universal curve","F-modes shift with DM, universality intact","Neutron-star f-modes: DM changes, law stands"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000157,"raw_usage":{"total_tokens":1204,"prompt_tokens":908,"completion_tokens":296,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":214}},"tokens_in":524,"tokens_out":296,"duration_ms":3590,"temperature":1.0,"reasoning_tokens":214,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:19:23.520765+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the f-mode frequencies for the same 26 equations of state with a fully general-relativistic oscillation code that includes metric perturbations; if the scatter around $f = 0.6282 + 2.0476\\sqrt{\\bar{M}/\\bar{R}^3}$ widens beyond the claimed 95% correlation, the universal relation is an artifact of the Cowling approximation. Observationally, a future post-merger gravitational-wave signal that yields an f-mode frequency for a neutron star with independently measured mass and radius would falsify the universal behavior if it falls off the fitted line by more than the scatter quoted here.","supporting_citations":[{"cited_title":"Andersson and K","cited_arxiv_id":null,"evidence_quote":"Empirical f-mode universal relations that the paper tests and extends to dark-matter-admixed quarkyonic stars."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original Cowling approximation, the oscillation scheme used for all f-mode calculations."}],"review_version":1}