{"id":"877bc4c3-daa6-4611-a115-506f3c042143","arxiv_id":"2412.17946","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Mirror dark matter inside proto-neutron stars reduces maximum mass, radius, and tidal deformability while heating the remnant and raising the speed of sound.","lead":"This paper models a newborn neutron star with a small amount of mirror dark matter in its core and finds that even 1% to 5% dark matter makes the remnant more compact, hotter, and harder to deform. The results suggest that dark matter could leave a measurable fingerprint in gravitational wave signals and neutron star cooling.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing weak point is the admitted assumption that dark and ordinary matter share the same entropy per baryon and lepton fraction at every evolutionary stage; relaxing it could change the magnitude or even the sign of the reported DM-induced compaction, heating, and tidal effects.","rationale":"Both the reader and I identify the same load-bearing assumption: the unverified equality of entropy per baryon and lepton fraction between the dark and visible sectors. This assumption is singled out by the authors themselves as having no established physical basis and being unlikely for gravity-only interactions. It directly controls the dark EoS, and because the two-fluid TOV equations couple the two fluids only through gravity, all stellar quantities are sensitive to it. For example, if the dark sector were colder and more degenerate, the dark core would be more compact, strengthening the gravitational pull and possibly increasing the compactification and heating effects; if it were hotter, the effects would weaken. Thus the quantitative claims in the abstract and Table III cannot be taken as predictions without this input. The sound-speed computation also uses an admitted cubic-spline interpolation rather than a rigorous two-fluid sound speed, but that is a secondary, self-acknowledged approximation affecting only one of the four headline effects. The mass-radius and tidal results are standard two-fluid TOV outcomes and are credible given the inputs. The paper is transparent and uses well-established machinery, and the qualitative message is plausible, so a conditional acceptance is appropriate; no change to the reader's verdict is needed. A concrete sensitivity run with alternative dark-sector thermal states would settle whether the concern is fatal or merely a quantitative caveat.","tokens_in":24502,"tokens_out":9553,"duration_ms":90529,"concrete_test":"Recompute the two-fluid TOV solutions for every stage with at least two alternative dark-sector thermal states while keeping the same visible-sector EoS, DM mass fractions (1%, 5%), and the same total gravitational mass: (i) cold degenerate dark matter with T_D = 0 (s_D = 0), and (ii) a hot dark sector with s_D = 1 (half the adopted entropy) or a different dark lepton fraction, e.g., Y_L,e' = 0. Compare Mmax, R1.4, Λ1.4, and the central temperature increase in Fig. 4. If any of the qualitative trends (DM reduces Mmax, R, and Λ; DM heats matter; DM increases sound speed) reverses or becomes non-monotonic, the equal-entropy assumption is load-bearing and the claims must be re-presented as conditional on an unverified coupling assumption. If the trends persist with similar magnitude, the central claim is robust to this assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative results depend on setting the dark sector's thermodynamic state equal to the visible sector's at each stage. In Sec. I the authors state there is no established reference on how DM thermalizes with OM, that thermal equilibrium is 'unlikely' for gravity-only coupling, yet they adopt the same entropy per baryon as 'a reasonable approximation' and also impose the same lepton fraction. Because the two fluids interact only gravitationally, the dark entropy and lepton fraction are free inputs that set the dark EoS and hence the compactness of the dark core. Changing s_D or the dark lepton fraction alters the DM pressure-density relation, which shifts the two-fluid TOV solutions: the maximum mass, radius, tidal deformability, and the central OM pressure (and therefore the extracted temperature and sound-speed profiles) are all affected. The authors are transparent about this, and the qualitative direction of the effects may survive, but the specific numbers in Table III and Figs. 4-6 are contingent on this single unverified assumption. A non-monotonic entry in Table III (s_B = 1, Y_L,e = 0.4, NH: R_2.1 = 14.04, 13.15, 13.88 km for 0%, 1%, 5% DM) further suggests the numerical two-fluid matching may not be fully robust, reinforcing the need to test the sensitivity of the conclusions to the dark-sector input.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the effects of a mirror dark matter component on proto-neutron stars (supernova remnants) using the DDME2 relativistic mean-field model for ordinary matter (with and without hyperons) and a mirror dark fermion model, within a two-fluid Tolman-Oppenheimer-Volkoff framework. The dark and visible sectors are assumed to have equal entropy per baryon and equal lepton fraction at each evolutionary stage, and 1% and 5% dark matter mass fractions are added to the stellar core. The authors compute mass-radius relations, tidal deformabilities, temperature profiles, particle fractions, and sound speeds for 2.1 solar mass stars. They conclude that dark matter compacts the remnant, reduces its maximum mass, radius, and tidal deformability, heats the matter, decreases isospin asymmetry, and increases sound speed.","tokens_in":24763,"tokens_out":7198,"duration_ms":64742,"significance":"If correct, this is one of the first studies of dark matter effects in hot proto-neutron stars with hyperons, extending the two-fluid DM-admixed neutron star literature to early evolutionary stages. The qualitative trends are plausible and consistent with prior cold-star DM studies, and the paper is transparent about its main assumption—equal entropy and lepton fraction between sectors—while also comparing with pulsar and GW170817 constraints. However, because the quantitative results are conditional on this unverified ansatz and on the fixed DM fractions, the work is exploratory rather than a definitive prediction.","major_comments":[{"comment":"The equal-entropy/lepton-fraction assumption is load-bearing. The authors state that thermal equilibrium is 'unlikely' for gravity-only coupling, yet they set s_D = s_B and Y_L,D = Y_L,e at every stage. Because the dark EoS—and thus the two-fluid TOV solutions for M_max, R, Λ, temperature, and sound speed—depends sensitively on these inputs, the quantitative results in Table III and Figs. 4–6 are contingent on this single unverified choice. Please provide a sensitivity study over a plausible range of s_D and Y_L,D (e.g., s_D = 0.5, 1, 2, 4; Y_L,D = 0.0, 0.2, 0.4) and state how the main conclusions change.","section":"Sec. I; Sec. II B; Eqs. (23)-(24)"},{"comment":"The reported R_2.1 values are 14.04 km (0% DM), 13.15 km (1% DM), and 13.88 km (5% DM). This non-monotonic behavior contradicts the paper's claim that increasing DM mass fraction monotonically reduces the radius. If this is not a typographical error, please explain the physical mechanism or re-examine the numerical matching; otherwise the conclusion 'DM reduces radius' is not robust.","section":"Table III, s_B=1, Y_L,e=0.4, NH row"},{"comment":"The sound speed is not computed from the two-fluid EoS. The authors evaluate c_s^2 from the OM EoS at fixed s_B and then use cubic spline interpolation to map the total pressure P(r) to c_s^2. This neglects the compressibility of the DM fluid. Since the paper claims DM increases the sound speed, please either compute the effective two-fluid sound speed from P = P_OM + P_D and ε = ε_OM + ε_D along the radial profile, or demonstrate that the interpolation is quantitatively accurate for 5% DM.","section":"Sec. IV, Fig. 5"},{"comment":"The text states that the capture rate (about 10^25 GeV/s) over a neutron star lifetime (10^17 s) implies the accumulated DM mass is 'likely insufficient to form a significant fraction of their total mass,' yet the study uses 1% and 5% mass fractions. This is inconsistent. Please clarify whether the 1% and 5% values are intended as upper bounds, parametric choices, or motivated by alternative accumulation scenarios (e.g., asymmetric DM, primordial accumulation, or self-capture), and discuss their astrophysical plausibility.","section":"Sec. V, DM mass fraction motivation"}],"minor_comments":[{"comment":"'Tolman-Oppenheiman-Volkoff' is misspelled; it should be 'Tolman-Oppenheimer-Volkoff'.","section":"Sec. I"},{"comment":"'µ_B′ the total baryon density of the dark sector' should read 'chemical potential' rather than 'density'.","section":"Sec. II B, text after Eq. (19)"},{"comment":"The caption states 'bottom-left: sB = 2 and YL,e = 0', which is inconsistent with the panel label 'sB=2; Yνe=0.0'; please align the caption with the panel.","section":"Fig. 3 caption"},{"comment":"The notation c_s^2 = ∂P/∂ε conflicts with the earlier definition c_s^2 = (dP_OM/dε_OM)|_{s_B}; please define both symbols explicitly to avoid ambiguity.","section":"Sec. IV, Fig. 5"},{"comment":"In reference [41], 'Neuton stars' should be 'Neutron stars'.","section":"Reference list"},{"comment":"The phrase 'For the first time' is a strong claim; consider softening it to 'We use...' or 'We present a...' unless a careful literature check fully supports the novelty.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable incremental contribution to the DM-admixed neutron star literature, extending it to hot proto-neutron stars with hyperons. The authors are transparent about their central assumption, which is commendable, but the quantitative predictions are conditional on the unverified equal-entropy/lepton-fraction ansatz and on the fixed DM fractions. The non-monotonic entry in Table III should be corrected or explained. The work fits the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a genuine first combination of finite-temperature proto-neutron-star (PNS) evolution with hyperons and a two-fluid mirror-DM treatment. The qualitative result — a few percent DM lowers maximum mass, radius, and tidal deformability, and heats the remnant — is plausible and consistent with older cold-star DM literature. The paper is worth reading, but the quantitative numbers rest on an admitted ansatz, so treat them as indicative, not predictive.\n\nWhat is new: the authors apply the two-fluid mirror-DM formalism to hot PNS stages with sB = 1 and 2, trapped and untrapped neutrinos, and T = 0, for both nucleonic and hyperonic DDME2 equations of state. The resulting mass-radius curves, tidal deformabilities, and temperature profiles are new tables and figures. They cite the relevant two-fluid DM and hot-star EoS literature, and the self-citations point to real earlier EoS papers; that is not a problem. They are also transparent about the main assumption — same entropy per baryon and lepton fraction in both sectors — and about the lack of a thermalization model.\n\nSoft spots: the equal-entropy/lepton-fraction assumption is load-bearing. Since the two fluids interact only gravitationally, the dark sector's thermodynamic state is a free input. If the dark entropy differs, the quantitative results shift. The authors acknowledge this but do not test sensitivity. A numerical red flag: Table III has a non-monotonic entry (sB = 1, YL,e = 0.4, NH: R2.1 = 14.04, 13.15, 13.88 km for 0%, 1%, 5% DM), suggesting the two-fluid matching or tabulation needs checking. The sound-speed extraction uses cubic spline interpolation between OM and total pressure; that is a reasonable first approximation, but it is not a fully self-consistent two-fluid sound speed.\n\nNone of this invalidates the qualitative conclusion; the trend that DM compacts and heats the star is likely robust to plausible changes in the dark-sector input. But the specific numbers should not yet be used as observational constraints.\n\nWho this is for: people working on dark matter in compact objects, PNS evolution, and equation-of-state modeling. It deserves a serious referee: the combination is new, the methods are mostly standard, and the limitations are stated. I would send it to review and ask for a sensitivity scan over the dark-sector thermodynamic state and a fix to the table inconsistency.\n\nRecommendation: engage with it as an exploratory model study. A revised version with sensitivity tests around the equal-entropy assumption would be genuinely useful.","headline":"New combination of hot PNS evolution and two-fluid mirror dark matter; the qualitative compaction/heating trends are plausible, but the numbers rest on an admitted equal-entropy ansatz and need sensitivity testing.","tokens_in":25364,"tokens_out":2077,"would_cite":true,"duration_ms":20984,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","97.60.Jd","26.60.-c"],"model":"deepseek-v4-flash","headline":"This paper argues that a few percent of mirror dark matter trapped in a newborn neutron star's core would measurably shrink its radius, lower its maximum mass, reduce its tidal deformability, heat the core, and speed up sound through the…","keywords":["mirror dark matter","proto-neutron stars","supernova remnants","two-fluid TOV","relativistic mean-field","hyperons","tidal deformability","equation of state"],"falsifier":"Re-solve the two-fluid Tolman-Oppenheimer-Volkoff equations with the dark fluid at a different entropy per baryon than the visible fluid (for instance, cold dark matter inside a hot ordinary core); if the predicted extra heating, reduced isospin asymmetry, and smaller tidal deformability vanish or reverse, the central result rests entirely on the equal-entropy assumption. Observationally, a precision mass-radius or tidal-deformability measurement of a compact object that matches the no-dark-matter curve would do the same.","tokens_in":24272,"feed_emoji":"🌌","tokens_out":9990,"duration_ms":88404,"temperature":0.7,"pith_summary":"The paper tries to establish that mirror dark matter, even at the few-percent level, reshapes the observable properties of supernova remnants from their birth as neutrino-rich proto-neutron stars through their cooling into ordinary neutron stars. Treating ordinary and dark matter as two fluids that interact only through gravity, and assuming that both sectors share the same entropy per baryon and lepton fraction at every evolutionary stage, the authors compute separate equations of state and then solve the two-fluid Tolman-Oppenheimer-Volkoff structure equations. They find that 1% and 5% dark-matter mass fractions in the remnant core make the star more compact: maximum mass, radius, and tidal deformability all decrease. They also find that dark matter heats the stellar matter, lowers neutron-proton asymmetry, shifts particle populations toward hyperons, and raises the squared speed of sound. If these results hold, they offer an indirect, gravity-only way to detect dark matter in compact stars and a possible explanation for scatter in neutron-star mass and radius measurements.","feed_headline":"Mirror dark matter shrinks and heats neutron-star remnants","feed_subtitle":"Predicted shifts in mass, radius, and tidal pull are large enough for pulsar and gravitational-wave tests.","key_machinery":"The machinery is a relativistic mean-field description with density-dependent couplings (the DDME2 parameterization) applied separately to ordinary matter and to a mirrored dark sector whose Lagrangian mimics the visible one with dark scalar and vector mesons. The two fluids are coupled thermodynamically by fixing the same entropy per baryon ($s_B$) and lepton fraction ($Y_{L,e}$) at each stage of remnant evolution, which determines temperature and composition as functions of density. Macroscopic quantities come from the two-fluid Tolman-Oppenheimer-Volkoff equations together with the differential equation for the Love number $k_2$, which yields the dimensionless tidal deformability $\\Lambda = (2/3)k_2 C^{-5}$ with compactness $C=M/R$. The same framework produces the pressure and mass profiles used to extract particle fractions, temperature profiles, and the speed of sound.","core_discovery":"The central claim is that dark matter need not annihilate or couple to ordinary matter to leave an imprint on a supernova remnant; gravitational interaction alone is enough. In a two-fluid picture with a mirrored dark sector (dark protons, dark neutrons, dark electrons, and dark neutrinos) coexisting with nucleonic and hyperonic matter, the paper shows that a fixed dark-matter mass fraction of 1% or 5% in the core systematically lowers the maximum gravitational mass (for cold nucleonic stars, from 2.48 to 2.33 $M_\\odot$ at 5% DM), shrinks the corresponding radius, and reduces the dimensionless tidal deformability of a canonical $1.4\\,M_\\odot$ star (for the same cold nucleonic case, from 704 to 356). The same dark matter heats the core by compressing the star through additional gravitational pull, decreases the isospin asymmetry $\\delta=(n_n-n_p)/(n_n+n_p)$ by favoring protons and hyperons over neutrons and leptons, and increases the squared speed of sound $c_s^2$ inside the star. The authors present these as model-dependent predictions that connect microscopic particle fractions to macroscopic observables.","pith_inferences":["Implicit in the results but not pursued in the paper: if the dark sector is colder than the visible sector rather than sharing its entropy, the predicted heating and sound-speed rise would probably weaken, so the quantitative numbers are best read as an upper bound on gravity-only dark-matter effects under the most favorable thermalization assumption.","One extension would be to apply the same two-fluid machinery to merger remnants, which are hotter and more massive than isolated proto-neutron stars; dark-matter-induced compactification could then show up in the post-merger gravitational-wave spectrum of next-generation detectors.","Because mirrored dark matter can form its own compact stars, the two-fluid framework also implies mixed visible-dark binaries; a dark companion would look like an unusually compact object with very low tidal deformability in a gravitational-wave catalog.","The predicted drop in isospin asymmetry could be cross-checked against nuclear symmetry-energy constraints: if terrestrial experiments fix the symmetry-energy slope, the dark-matter fraction needed to produce a given shift in $\\delta$ becomes a testable parameter."],"forward_implications":["A 5% mirror-dark-matter mass fraction lowers the maximum mass of cold nucleonic stars from 2.48 to 2.33 $M_\\odot$ and of cold hyperonic stars from 2.26 to 2.11 $M_\\odot$, shifting where remnants sit in the mass-radius plane.","The canonical $1.4\\,M_\\odot$ tidal deformability drops sharply with dark matter (from 704 to 356 in the cold nucleonic case at 5% DM), moving stars toward the observed binary-merger tidal constraint and changing the expected inspiral signal.","Dark matter heats the core by gravitational compression, most visibly during deleptonization, which would prolong the cooling phase and affect estimates of a young remnant's age.","Dark matter favors hyperon production and lowers isospin asymmetry, making the onset radii of $\\Lambda$, $\\Xi$, and $\\Sigma$ species inside the star sensitive to dark-matter content.","Because the same entropy and lepton-fraction assumptions are applied at every stage, the predictions form a unified evolutionary sequence from the neutrino-rich birth state to the cold catalyzed neutron star."],"supporting_citations":[{"why":"It supplies the DDME2 density-dependent relativistic mean-field parameterization used for both ordinary and dark matter equations of state.","marker":"[65]"},{"why":"It supplies the SU(3)/SU(6) baryon coupling scheme that fixes hyperon couplings and keeps maximum masses above the 2 $M_\\odot$ threshold.","marker":"[66]"},{"why":"It provides the mirror/fermionic dark matter model and its zero-temperature equation of state that this work extends to finite temperature.","marker":"[33]"},{"why":"It supplies the two-fluid fermionic dark matter neutron star framework and the treatment of dark matter parameters used in the present model.","marker":"[22]"},{"why":"It provides the two-fluid Tolman-Oppenheimer-Volkoff equations and boundary conditions for dark-matter-admixed neutron stars.","marker":"[87]"},{"why":"It establishes the proto-neutron-star evolution scheme with fixed entropy per baryon and lepton fraction over the Kelvin-Helmholtz timescale.","marker":"[58]"},{"why":"It provides the hot neutron-star equation of state with exotic baryons from which the temperature and composition profiles are taken.","marker":"[76]"},{"why":"It supplies the GW170817 tidal deformability constraint used to evaluate the dark-matter-admixed results.","marker":"[106]"}],"fun_headline_variants":["Dark mirror sector slims supernova remnants","Mirror dark matter cuts remnant mass, radius, and tidal pull","Gravity-only dark matter leaves its mark on remnants","Mirror dark matter alters remnant cores and sound speed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that dark matter and ordinary matter share the same entropy per baryon and lepton fraction at every evolutionary stage, even though the paper acknowledges that thermal equilibrium between the two sectors is unlikely when they interact only gravitationally.","fun_headline_variants_meta":{"raw":{"variants":["Dark mirror sector slims supernova remnants","Mirror dark matter cuts remnant mass, radius, and tidal pull","Gravity-only dark matter leaves its mark on remnants","Mirror dark matter alters remnant cores and sound speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000781,"raw_usage":{"total_tokens":3455,"prompt_tokens":953,"completion_tokens":2502,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":2438}},"tokens_in":569,"tokens_out":2502,"duration_ms":21442,"temperature":1.0,"reasoning_tokens":2438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:07:54.114746+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-solve the two-fluid Tolman-Oppenheimer-Volkoff equations with the dark fluid at a different entropy per baryon than the visible fluid (for instance, cold dark matter inside a hot ordinary core); if the predicted extra heating, reduced isospin asymmetry, and smaller tidal deformability vanish or reverse, the central result rests entirely on the equal-entropy assumption. Observationally, a precision mass-radius or tidal-deformability measurement of a compact object that matches the no-dark-matter curve would do the same.","supporting_citations":[],"review_version":1}