{"id":"f631c43c-1d90-46fe-a25e-6f4e4c75b35c","arxiv_id":"2505.20545","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Integrating two-fluid TOV equations with MIT bag quark matter and BEC dark matter shows that admixed configurations are generally more compact and less massive than pure quark stars, with the effect depending on the sign of pressure anisotropy.","lead":"The paper numerically solves Einstein's equations for hypothetical stars containing two fluids, quark matter and condensed dark matter, and shows that adding dark matter changes how massive and compact the stars can be. It is a parameter study that could inform interpretation of neutron star mass-radius observations if such mixed stars exist.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Internal contradiction: abstract claims DM always yields more compact, less massive stars, but §5.2 reports the opposite for positive anisotropy.","rationale":"I agree with the reader that the paper has a significant correctness risk and merits a CONDITIONAL verdict. However, I identify a different load-bearing concern than the reader's weakest_assumption. The common-surface boundary condition (p(R)=0 for both fluids at one radius) is a modeling assumption that could shift the M-R curves, but the authors could defend it as a deliberate choice (e.g., if both fluids are confined by the same total pressure). In contrast, the internal contradiction between the abstract/summary and the Case II results is a logical flaw that directly undermines the central claim as stated. The paper cannot simultaneously assert that DM always makes quark stars more compact and less massive while reporting the opposite for positive anisotropy. This must be resolved before any comparison with observations is meaningful. The proposed test—recomputing the Case II M-R curves—would settle the contradiction by determining which statements match the actual numerical solutions. Since the underlying numerical framework appears standard and the issue is correctable by revising the conclusions (or the code, if the text is right), CONDITIONAL rather than REJECT is appropriate. The reader's rationale already noted this contradiction, so my agreement is partial: we agree on the existence of the contradiction but disagree on which concern is most load-bearing.","tokens_in":19420,"tokens_out":4623,"duration_ms":45723,"concrete_test":"Recompute the M-R sequences for Case II (κ1=+0.1, f=0.58, 0.60, 0.62) and for pure quark matter using the same integration scheme and EoS parameters (B=57.64 MeV/fm^3, K=0.01/B, κ2=-0.2). Extract M_max, the radius at M_max, and the compactness M/R at a common central density (e.g., ρc,QM/ρs ≈ 2.75). If the DM-admixed curve has lower M/R and higher M_max than the pure quark curve, then §5.2 is correct and the abstract's universal claim is false; if it has higher M/R and lower M_max, then §5.2 is wrong. Either outcome determines which statement must be corrected.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim, stated in the abstract and repeated in the conclusions ('regardless of whether the anisotropic quark matter satisfies ∆QM > 0 or < 0, the presence of dark matter... results in a more compact and lighter object'), is directly contradicted by the results presented for Case II. In §5.2 the authors write that for positive quark anisotropy (κ1=+0.1) the hybrid configurations are 'less compact yet more massive structures compared to pure quark stars', and Figure 2 is cited as showing this. The summary itself opens with the opposite statement ('Conversely, for positive quark anisotropy... less compact yet more massive'). Thus the universal claim in the abstract is not supported by the paper's own data; the sign of the DM effect on compactness and maximum mass depends on the sign of κ1. This is not a wording issue: the central comparison with pure quark stars, and hence the alignment with gravitational-wave observations, is ambiguous. The reader's verdict identified this contradiction in the rationale, though the weakest_assumption focused on the common-surface boundary condition. The latter is a modeling choice, whereas the contradiction is an internal logical inconsistency that must be resolved before the M-R sequences can be interpreted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript numerically integrates the two-fluid Tolman-Oppenheimer-Volkoff equations for static, spherically symmetric compact stars containing quark matter described by the MIT Bag model and dark matter described by a Bose-Einstein condensate polytropic EoS, with an anisotropy ansatz Δ = κ p (2m/r) for each fluid. The authors fix B = 57.64 MeV/fm^3, K = 0.01/B, κ2 = -0.2, κ1 = ±0.1, and vary the central density ratio α = f/(1-f) = ρ_c,DM/ρ_c,QM. They present mass-radius curves, compactness factors, quark mass fractions, and an adiabatic-index stability check for negative and positive quark anisotropy, overlaying pulsar mass bands and NICER regions. The paper's headline conclusion is that dark-matter-admixed quark stars are more compact and less massive than pure quark stars, but Section 5.2 reports the opposite behavior for positive quark anisotropy.","tokens_in":19681,"tokens_out":7819,"duration_ms":84031,"significance":"If the technical issues are resolved, this parameter study would provide a useful benchmark family of two-fluid equilibrium sequences with transparent EoS choices and explicit observational overlays. The framework is standard, the input parameters are stated concretely, and the comparison with PSR J1614-2230, PSR J0348+0432, PSR J0740+6620, and NICER constraints is a useful way to frame the results. The main limitation is that the advertised universal conclusion is not supported by the paper's own Case II results, and the common-surface and stability assumptions require scrutiny before the quoted masses and radii can be interpreted as physical predictions.","major_comments":[{"comment":"The abstract and the final paragraph of §6 state as a universal result that dark matter-admixed quark stars are 'more compact and lighter' than pure quark stars 'regardless of whether the anisotropic quark matter satisfies ΔQM>0 or ΔQM<0'. This is directly contradicted by §5.2, where for positive quark anisotropy (κ1=+0.1) the hybrid configurations are described as 'less compact yet more massive structures compared to pure quark stars', with the maximum mass slightly increasing as α increases. Since this universal statement is the paper's headline claim and the basis for the claimed alignment with gravitational-wave observations, it must be corrected to a sign-dependent statement or restricted to the negative-anisotropy case.","section":"Abstract; §5.2; §6"},{"comment":"The boundary condition p(R)=0 in Eq. (14), with p=pQM+pDM from Eq. (9), assumes that both fluids terminate at one common radius and that the interior can be matched to a single Schwarzschild exterior. Because the quark and dark matter fluids have independent equations of state and interact only gravitationally, there is no a priori reason for pQM and pDM to vanish at the same radius; one fluid could extend beyond the other, and continuing the MIT-bag EoS beyond pQM=0 would produce negative quark pressure. The manuscript does not verify simultaneous vanishing of the two pressures. Since the reported M, R, and compactness all depend on this surface choice, the authors should either demonstrate numerically that the common-surface condition is satisfied for the parameter sets used or reformulate the matching with separate surfaces.","section":"§3, Eqs. (9), (14), (16)"},{"comment":"Stability is asserted using the approximate adiabatic-index criterion Γ≥Γcr=4/3+(19/21)(M/R), and §5.2 argues that stability of the most massive pure quark configuration 'guarantees that all other configurations must be stable as well' because their masses and compactnesses are lower. This inference is not valid for two-fluid stars, as the paper itself acknowledges via Refs. [135,136]: the maximum-mass configuration need not coincide with the last stable configuration, and extended stable branches can exist. The stability claims should be downgraded to consistency with the approximate Γ criterion, or supported by a proper radial-oscillation analysis of the two-fluid configurations.","section":"§5, Eq. (31), Fig. 3"}],"minor_comments":[{"comment":"Several typographical errors should be corrected: 'posotive' after Eq. (1), 'Alternativelly' in Section 1, 'thorugh' in the Acknowledgments, and 'eigenvalue value problem' in Section 5.2.","section":"Throughout"},{"comment":"The caption lists 'the light HESS compact object (purple region)' twice, attributing it to Refs. [146] and [147]; Ref. [147] is a NICER paper on PSR J0437-4715, so the region labels and references should be reconciled.","section":"Fig. 1 caption"},{"comment":"The quantity plotted as 'Mass fraction' is not defined in the text; please state explicitly whether it is M_QM/M_total and how it is obtained from the integrated density profiles.","section":"§5.1, Fig. 1 bottom-right panel"},{"comment":"The units of K are not specified; since B is given in MeV/fm^3 and the TOV equations are integrated in geometric units, please state the conversion used for K and the resulting units of the polytropic EoS.","section":"§4.1, Eqs. (20), (27), (28)"},{"comment":"The sentence 'As the anisotropic factor of dark matter condensate is always negative [107]' should be qualified, since with the ansatz of Eq. (21) the sign is set by the input parameter κ2, which is chosen as -0.2 here.","section":"§5, paragraph before Eq. (25)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a straightforward parameter study with a standard numerical setup; the referee report above contains the substantive issues. In particular, the internal contradiction between the abstract and Section 5.2 needs to be resolved in revision, and the surface-matching and stability assumptions need to be checked or carefully qualified before the quoted masses, radii, and compactnesses can be used for astrophysical comparisons."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read of arXiv:2505.20545. The paper is a two-fluid TOV parameter scan: MIT Bag quark matter plus a BEC polytrope dark matter, with an anisotropy ansatz proportional to κ p (2m/r). The new content is the systematic map over the DM fraction f and the sign of the quark anisotropy κ1, along with the quark mass-fraction curves. The equations and integration are standard, the text is transparent about the model, and the numerics look reproducible. Credit where due: the authors explicitly note that their adopted singlet-scalar DM model has no DM-nucleon coupling, so the usual capture/thermalization discussion doesn't apply. That is honest framing.\n\nThe problem is that the paper's own headline does not survive contact with its results. The abstract and the final summary state that, regardless of the sign of ΔQM, adding DM produces a more compact and lighter star. But Section 5.2 says the opposite for positive anisotropy: the hybrid configurations are 'less compact yet more massive' than pure quark stars, and Figure 2 shows that. The concluding section contains both statements within a few paragraphs. This is not a wording slip; it makes the central claim ambiguous and the GW/NICER comparison misleading.\n\nThe common-surface boundary condition flagged by the reader is also fair. The two fluids interact only gravitationally, so there is no a priori reason their pressures vanish at the same radius. The paper imposes p(R)=0 and matches a single Schwarzschild exterior without checking which component has the larger radius. If DM forms an outer layer, the reported masses and radii shift. That needs a sentence of justification or a relaxed-surface calculation.\n\nMinor issues: the Γ-based stability argument is explicitly acknowledged as approximate; the paper cites radial-oscillation studies for multi-fluid stars and notes the maximum-mass configuration need not be the last stable one. The alignment with 'gravitational wave observations' is asserted rather than demonstrated quantitatively.\n\nIn short: the numerics are probably fine, and the parameter scan is a modest but usable extension of known results. But the internal contradiction in the headline must be fixed before anyone cites this for the sign of the DM effect. A referee should ask for that fix, plus the surface check, plus a clearer stability statement.\n\nI'd send it to peer review—it's fixable and the underlying computation seems solid—but I wouldn't cite it in its current form.","headline":"A standard two-fluid TOV parameter scan whose abstract contradicts its own results for positive anisotropy; the numerics look solid but the headline needs fixing before citation.","tokens_in":20192,"tokens_out":2633,"would_cite":false,"duration_ms":29257,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C55","85A15"],"pacs":["04.40.-b","95.35.+d","97.60.Jd"],"model":"deepseek-v4-flash","headline":"Dark matter-admixed quark stars are more compact yet less massive than pure quark stars, with the effect controlled by quark anisotropy.","keywords":["anisotropic compact stars","quark stars","dark matter admixed stars","two-fluid Tolman-Oppenheimer-Volkoff equations","MIT bag model","Bose-Einstein condensate dark matter","mass-radius relation","compactness"],"falsifier":"Integrate the same two-fluid TOV system but stop each fluid where its own pressure alone vanishes and compare the two radii; if $R_{DM} \\neq R_{QM}$ for the quoted $f$ and $\\kappa$ values, the single-surface matching assumed here is invalid and the reported masses, radii, and compactness values are not physical. On the observational side, a confirmed pulsar with mass at or above $2\\,M_\\odot$ and a radius in the range predicted by the negative-anisotropy models would rule out those parameter combinations.","tokens_in":19217,"feed_emoji":"🌌","tokens_out":11081,"duration_ms":91644,"temperature":0.7,"pith_summary":"This paper numerically computes the structure of compact stars made of two non-interacting fluids, quark matter and condensed dark matter, by solving the two-fluid Tolman-Oppenheimer-Volkoff equations. Its central claim is that dark matter-admixed quark stars are more compact yet less massive than pure quark stars, and that this effect is controlled by two dials: the central dark matter fraction and the sign of the star's pressure anisotropy. With negative quark anisotropy the computed maximum masses often fall below the two-solar-mass band set by known pulsars, while with positive anisotropy the trend reverses. The paper presents these mass-radius curves as testable predictions that can be compared with X-ray, radio timing, and gravitational-wave constraints on neutron star structure.","feed_headline":"Dark matter makes quark stars more compact, less massive","feed_subtitle":"A two-fluid model says dark matter fraction and quark anisotropy set the curves that X-ray and pulsar timing can test.","key_machinery":"The central object is the two-fluid TOV system: $m' = 4\\pi r^2 \\rho$, $p'_{QM} = -(\\rho_{QM}+p_{QM})\\,g + 2\\Delta_{QM}/r$, $p'_{DM} = -(\\rho_{DM}+p_{DM})\\,g + 2\\Delta_{DM}/r$, with $g(r) = (m + 4\\pi r^3 p)/(r^2(1-2m/r))$, closed by the quark matter relation $p_{QM} = (\\rho_{QM} - 4B)/3$ and the dark matter condensate relation $p_{DM} = K\\rho_{DM}^2$. The anisotropy is encoded by $\\Delta_{QM} = \\kappa_1 p_{QM}(2m/r)$ and $\\Delta_{DM} = \\kappa_2 p_{DM}(2m/r)$. The integration is normalized by the central ratio $\\alpha = f/(1-f)$ with $f = \\rho_{c,DM}/(\\rho_{c,DM}+\\rho_{c,QM})$, and the star is defined by the single-surface matching conditions $p(R)=0$, $m(R)=M$, joined to a Schwarzschild exterior.","core_discovery":"Working with the MIT bag equation of state $p_{QM} = (\\rho_{QM} - 4B)/3$, a polytropic condensate equation of state $p_{DM} = K\\rho_{DM}^2$, and anisotropy terms $\\Delta_{QM} = \\kappa_1 p_{QM}(2m/r)$, $\\Delta_{DM} = \\kappa_2 p_{DM}(2m/r)$, the paper integrates the two-fluid TOV system outward from central densities related by $\\alpha = f/(1-f)$, with $f = 0.58, 0.60, 0.62$. For negative quark anisotropy ($\\kappa_1 = -0.1$), increasing the dark matter fraction produces stars that are more compact yet less massive than pure quark stars, with maximum masses that typically fall below the $\\sim 2\\,M_\\odot$ constraint from PSR J1614-2230 and PSR J0348+0432. For positive quark anisotropy ($\\kappa_1 = +0.1$), the same fractions produce less compact yet more massive configurations, occasionally exceeding the pure-quark maximum mass. The quark mass fraction peaks near $\\rho_{c,QM} \\approx 2.75\\,\\rho_s$ and is shifted downward as $f$ grows, and the maximally massive configurations satisfy the adiabatic index stability criterion $\\Gamma \\ge \\Gamma_{cr} = 4/3 + (19/21)(M/R)$. The authors take these results to show that dark matter content and anisotropy jointly shape the observational signature of quark-matter cores and to align with recent theoretical predictions and gravitational wave observations.","pith_inferences":["If the single-surface assumption is relaxed, dark matter would naturally form an outer halo around the quark core; the exterior would no longer be a single Schwarzschild spacetime, and the inferred mass-radius relation could shift, potentially bringing the negative-anisotropy models back above the two-solar-mass threshold.","The reversal of the dark matter effect with the sign of $\\kappa_1$ suggests that anisotropy, not just dark matter content, is the controlling lever; a single accurate mass-radius measurement of a compact object could therefore distinguish between the two anisotropy mechanisms and, by extension, between the underlying microphysical sources.","A natural follow-up, which the paper states is in progress, is to compute tidal deformabilities for these models; enforcing the GW170817 bound $\\Lambda \\le 800$ would further shrink the allowed $(f,\\kappa_1)$ region and could turn the predicted curves into a sharper test.","If such stars exist, their compactness encodes both the boson self-interaction strength of the dark condensate and the bag constant of quark matter, so precise radius measurements would constrain fundamental constants of both sectors at once."],"forward_implications":["If the negative-anisotropy models are correct, the presence of dark matter lowers the maximum mass of quark stars, so a confirmed two-solar-mass pulsar would rule out the largest dark matter fractions studied here.","For positive quark anisotropy the trend reverses, so measuring a single precise mass and radius for a candidate quark star could indicate which sign of anisotropy is realized in nature.","The quark mass fraction peaks near $\\rho_{c,QM} \\approx 2.75\\,\\rho_s$ and shifts downward as $f$ increases, giving a direct relation between the inferred central density and the dark matter content.","Because the paper's particle model has no dark matter-nucleon interaction, the large dark mass fractions shown are not in conflict with capture-based limits, so the predicted curves are legitimate targets for observation."],"supporting_citations":[{"why":"Supplies the two-fluid Tolman-Oppenheimer-Volkoff equations used to integrate the coupled quark and dark matter fluids.","marker":"[87, 88]"},{"why":"Provides the MIT bag model equations of state for quark matter, fixing $p_{QM} = (\\rho_{QM} - 4B)/3$.","marker":"[89, 90]"},{"why":"Shows that the self-interacting scalar field equation of state reduces to a pressure proportional to density squared at low density, justifying the dark matter polytrope.","marker":"[105]"},{"why":"Establishes that condensate dark matter anisotropy is negative, supporting the $\\kappa_2 = -0.2$ choice.","marker":"[107]"},{"why":"Provides the adiabatic index stability criterion $\\Gamma_{cr} = 4/3 + (19/21)(M/R)$ used to declare the configurations stable.","marker":"[115]"},{"why":"Gives the dark matter induced anisotropy compact star framework whose compactness trends the present results are said to align with.","marker":"[84]"},{"why":"Supplies the PSR J1614-2230 mass measurement used as a $\\sim 2\\,M_\\odot$ constraint band.","marker":"[119]"},{"why":"Supplies the PSR J0348+0432 mass measurement used as a second $\\sim 2\\,M_\\odot$ constraint band.","marker":"[120]"},{"why":"Provides the X-ray radius measurement of PSR J0740+6620 used to restrict the allowed mass-radius regions.","marker":"[122]"}],"fun_headline_variants":["Dark matter shrinks quark stars but trims their mass","Two-fluid stars: dark matter compacts quark cores","Quark stars get denser, lighter with dark matter","Dark matter-admixed quark stars: compact yet light","Anisotropy and dark matter reshape quark star curves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes both fluids terminate at exactly the same radius, with the total pressure $p_{QM}+p_{DM}$ vanishing there, and never checks whether the two pressures separately reach zero at that point; if the dark matter extended further out, the star's mass, radius, and compactness would all change.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter shrinks quark stars but trims their mass","Two-fluid stars: dark matter compacts quark cores","Quark stars get denser, lighter with dark matter","Dark matter-admixed quark stars: compact yet light","Anisotropy and dark matter reshape quark star curves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000154,"raw_usage":{"total_tokens":1253,"prompt_tokens":1028,"completion_tokens":225,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":644,"completion_tokens_details":{"reasoning_tokens":147}},"tokens_in":644,"tokens_out":225,"duration_ms":2773,"temperature":1.0,"reasoning_tokens":147,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:52:28.388440+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the same two-fluid TOV system but stop each fluid where its own pressure alone vanishes and compare the two radii; if $R_{DM} \\neq R_{QM}$ for the quoted $f$ and $\\kappa$ values, the single-surface matching assumed here is invalid and the reported masses, radii, and compactness values are not physical. On the observational side, a confirmed pulsar with mass at or above $2\\,M_\\odot$ and a radius in the range predicted by the negative-anisotropy models would rule out those parameter combinations.","supporting_citations":[{"cited_title":"Colpi, S","cited_arxiv_id":null,"evidence_quote":"Shows that the self-interacting scalar field equation of state reduces to a pressure proportional to density squared at low density, justifying the dark matter polytrope."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that condensate dark matter anisotropy is negative, supporting the $\\kappa_2 = -0.2$ choice."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the adiabatic index stability criterion $\\Gamma_{cr} = 4/3 + (19/21)(M/R)$ used to declare the configurations stable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the dark matter induced anisotropy compact star framework whose compactness trends the present results are said to align with."},{"cited_title":"Shapiro Delay Measure- ment of A Two Solar Mass Neutron Star","cited_arxiv_id":null,"evidence_quote":"Supplies the PSR J1614-2230 mass measurement used as a $\\sim 2\\,M_\\odot$ constraint band."},{"cited_title":"A Massive Pulsar in a Compact Relativistic Binary","cited_arxiv_id":null,"evidence_quote":"Supplies the PSR J0348+0432 mass measurement used as a second $\\sim 2\\,M_\\odot$ constraint band."},{"cited_title":"Riley et al","cited_arxiv_id":null,"evidence_quote":"Provides the X-ray radius measurement of PSR J0740+6620 used to restrict the allowed mass-radius regions."}],"review_version":1}