{"id":"291cbac2-c4d6-4715-849a-658d71cea8e4","arxiv_id":"2505.01503","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Solar oscillations constrain a compact dark matter core in the Sun to below one hundred-thousandth of a solar mass, and a one-thousandth-solar-mass core improves helioseismic agreement by mimicking a heavy metal core.","lead":"The authors model the Sun with a compact dark matter core and compare the model Suns to neutrino and oscillation data. They find the Sun's vibrations constrain such cores far better than neutrinos, and a one-thousandth-solar-mass core surprisingly improves the fit, which they interpret as mimicking a metal-rich core.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dark-core mass limits are degenerate with the assumed inner-boundary radius: the paper itself states that varying this radius mimics changing the core mass, so the constraints are not uniquely on M0.","rationale":"The paper does what it sets out to do: it computes calibrated MESA solar models with a central point mass, publishes the code, and compares neutrino fluxes and p-mode frequencies to observations. The neutrino result is robust in the sense that only the 10^-2 model produces measurable flux changes. The p-mode sensitivity claim, however, rests on the fixed-Bondi-radius boundary prescription. The authors explicitly report that the Bondi radius changes by an order of magnitude over the evolution and that changing the inner boundary radius has a similar structural effect to changing the dark core mass. That statement, taken at face value, means the observable frequency shifts and chi-squared values are not a function of M0 alone; they are a function of (M0, r_B). Because the class of candidates considered is defined only by being smaller than the Bondi radius, the derived mass thresholds and the claimed 10^-3 improvement lack a unique interpretation until the degeneracy is quantified. This is a correctness risk, not merely a modeling detail: the same observations could be reproduced by a different mass if the core radius or truncation is adjusted. The proposed test—varying r_B around its fiducial value for a fixed M0—would directly show whether the conclusions shift. If they do not shift, the concern is resolved; if they do, the paper should report constraints in the (M0, r_B) plane or restrict to a specific physical model of the compact object. The mixed-mode tension noted in Section 4.1 is a further reason the 10^-3 improvement is not a detection, but it does not replace the boundary degeneracy as the primary concern. This supports the reader's conditional verdict rather than changing it.","tokens_in":20561,"tokens_out":13266,"duration_ms":144556,"concrete_test":"Recompute the calibrated models and Gyre frequencies for at least one fixed dark core mass, say M0 = 10^-3 Msun, with the inner boundary placed at r_B/3, r_B, and 3 r_B, and, if feasible, with a time-dependent r_B using the evolving central sound speed. Compare the resulting radial p-mode frequency differences and chi-squared values for the same observed mode set and surface correction. If the spread across boundary choices is comparable to the difference between adjacent grid masses in Figure 5, or if the 10^-5 threshold shifts by more than a factor of a few, the constraints should be restated as a joint constraint on (M0, r_B) rather than on M0 alone.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claims—the 10^-5 Msun p-mode threshold, the 10^-3 Msun improvement, and the 10^-7 Msun g-mode projection—are all computed with the MESA inner boundary fixed at the Bondi radius r_B = 2GM0/c_s^2 (Section 2). The authors report that r_B can vary by about an order of magnitude during the evolution as the central sound speed changes, and that varying the inner radius boundary condition has 'a very similar effect on the resulting solar structure as when changing the dark core mass.' If the boundary location is as influential as the mass itself, then the observables do not constrain M0 alone; they constrain a degenerate combination of M0 and the assumed core radius/truncation radius. Since the paper's candidate class is defined only by spatial extent less than the Bondi radius, a candidate of mass M0 with a different physical radius—or a time-dependent Bondi radius—would masquerade as a different core mass in the frequency shifts, chi-squared comparisons, and g-mode spacings. The reported mass limits and the claimed 10^-3 improvement are therefore not robust statements about dark core mass until this degeneracy is quantified. The non-accreting assumption is a separate modeling restriction and is not the issue here; the boundary-radius ambiguity affects even the restricted non-accreting case.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses MESA solar evolution models with a point-mass central object (a 'dark core') of mass 10^-8 to 10^-2 M_sun, calibrated to solar luminosity, radius, and surface composition. The inner boundary is fixed at the Bondi radius. The authors compute neutrino fluxes and GYRE pulsation frequencies, then compare with observed solar neutrino fluxes and helioseismic data. They conclude that neutrino fluxes only rule out M0 >= 10^-2 M_sun, p-mode frequencies rule out M0 > 10^-5 M_sun, a 10^-3 M_sun dark core improves agreement with helioseismic data (chi2_r = 308 vs 2633), and future g-mode spacings could probe down to 10^-7 M_sun. They attribute the apparent improvement to a metal-rich core that may emulate star-formation effects rather than dark matter.","tokens_in":20842,"tokens_out":6419,"duration_ms":68189,"significance":"If the central claims hold, the paper provides a new, model-dependent route to constrain macroscopic dark matter using the Sun as a laboratory. It is clearly written, uses publicly available codes (MESA, GYRE, and a public GitHub repository), and makes falsifiable predictions (mixed modes, g-mode period spacings) that distinguish the dark-core scenario from standard solar physics. The qualitative result that a compact central mass can alter the temperature and composition structure so as to mimic a metal-rich core is interesting and connects to discussions of the solar abundance problem. However, the quantitative limits are compromised by an unquantified degeneracy between M0 and the assumed inner boundary radius, the neutrino analysis is based on only two flux components, and the claimed 10^-3 improvement is based on a selected subset of oscillation data while the paper itself notes that mixed modes may rule out that mass. These issues are central to the paper's main claims.","major_comments":[{"comment":"The paper states that varying the inner radius boundary condition 'had a very similar effect on the resulting solar structure as when changing the dark core mass.' Since the inner boundary is set to the Bondi radius r_B = 2GM0/c_s^2 and then held fixed, the constraints reported in Sections 4 and 5 are not on M0 alone but on the joint configuration (M0, r_B). The paper does not quantify how the p-mode limit at ~10^-5 M_sun, the g-mode projection at ~10^-7 M_sun, or the 10^-3 M_sun improvement shift when r_B is varied independently or when the Bondi radius evolves by the order-of-magnitude variation the authors mention. This degeneracy is load-bearing because the claimed mass limits and the apparent improvement could be reinterpreted as constraints on the boundary location rather than on the dark-core mass. A quantitative exploration of the M0-r_B degeneracy, or a physical argument fixing r_B uniquely for each candidate class, is required before the stated limits can be accepted.","section":"Section 2, fixed inner boundary"},{"comment":"The neutrino constraints are derived using only the pp and CNO fluxes. The 8B flux, and to a lesser extent the 7Be and pep fluxes, depend much more steeply on the central temperature (8B approximately as T_c^20) and could be far more sensitive to the presence of a compact core than the pp flux, which is nearly fixed by the solar luminosity. The paper does not report these fluxes or compare them with the observational constraints that are already available. Therefore the conclusion that neutrino measurements only rule out dark-core masses above ~10^-2 M_sun is not established by the presented analysis. The complete set of solar neutrino fluxes should be computed and compared with measurements before claiming that neutrinos provide the weakest constraints.","section":"Section 3, Figure 3"},{"comment":"The chi-squared improvement reported for the 10^-3 M_sun model is computed from radial p-modes only, after a surface correction. The paper itself notes that in this model some mixed modes take the place of ordinary non-radial p-modes and 'may rule out such a massive dark core.' A comparison based on a selected subset of the helioseismic data, while the same data contain modes that are in high tension with the model, does not by itself establish that the 10^-3 model is a better representation of the Sun. A quantitative assessment including non-radial modes, or a clear justification for excluding them, is necessary before the apparent improvement can be presented as a meaningful finding. This does not negate the interest of the result, but it does affect the strength of the central claim.","section":"Section 4 and Section 4.1"}],"minor_comments":[{"comment":"The phrase 'we calibrate standard solar evolution models' may be misleading, because models with a dark core are not standard; consider saying 'we calibrate solar evolution models' or 'solar models.'","section":"Abstract"},{"comment":"The label 'Dark - Standard' in the left panel is not defined in the caption; please state explicitly that this is the theoretical frequency difference between a dark-core model and the standard model.","section":"Figure 4 caption"},{"comment":"The instrument name 'Mesa' appears in the text while the standard capitalization is 'MESA'; likewise 'Gyre' should be 'GYRE'. Please unify the capitalization for consistency.","section":"Throughout"},{"comment":"The non-monotonic behavior of the g-mode period spacing for the 10^-2 M_sun model is explained by a smaller propagation cavity, but it would help to state explicitly that this is a competing effect with the increased buoyancy frequency in the core.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to appeal to both the stellar astrophysics and dark matter communities, and the public code release is a strength. However, the unquantified M0-r_B degeneracy and the incomplete neutrino flux analysis are central to the main quantitative claims. The authors should be encouraged to address these points with additional computations rather than simple discussion. The improvement at 10^-3 M_sun is intriguing but, as the authors themselves acknowledge, is not robust evidence for a dark core; the presentation should more clearly separate this speculative result from the actual constraints."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the punchline: this paper does something genuinely new and does it with transparent, reproducible modeling. It treats the Sun with a non-accreting compact dark core that interacts only gravitationally, calibrates MESA solar models with core masses from 1e-8 to 1e-2 Msun, and derives the first helioseismic and neutrino flux limits for that general class of macroscopic dark matter. The code and data are public, the microphysics is standard, and the comparisons to observations are clear. The neutrino bounds are weak (only about 1e-2 Msun is ruled out), the p-mode acoustic data do the heavy lifting (around 1e-5 Msun), and the g-mode projection (down to 1e-7 Msun) is a reasonable forward look. That core argument holds up.\n\nCredit where it's due: this is a solid new application of existing modeling tools rather than a new framework. The paper also ties the dark core prescription to the solar modeling problem, showing that a 1e-3 Msun core induces a metal-rich center that improves agreement with helioseismic data, while explicitly warning against reading that as dark matter evidence. The literature coverage is fair, and the self-citations to the authors' own PBH work are appropriate.\n\nThe main soft spot is exactly the one stressed: the fixed Bondi radius degeneracy. The paper itself says (Section 2) that varying the inner radius has a very similar effect on the solar structure as changing the dark core mass, and that the Bondi radius can change by an order of magnitude over the evolution. So the quoted mass limits are not uniquely on M0; they constrain a combination of M0 and the assumed core/truncation radius. This does not sink the qualitative ordering—neutrinos weak, p-modes stronger, g-modes most sensitive—but it does mean the headline numbers should be treated as indicative, not precise. The 1e-3 improvement claim is the most fragile. It comes from an a posteriori grid selection, the chi-squared comparison does not include model systematics, and the same model has mixed modes in tension with observations. The authors are appropriately cautious about overinterpreting it, but a referee should ask for variable Bondi radius tests and an explicit discussion of the degeneracy.\n\nMinor: the non-accreting assumption is a real restriction, but it is clearly stated and defensible as a first step. The g-mode projection is speculative but labeled as such.\n\nBottom line: this deserves a serious referee. I would send it to review, with the expectation that the degeneracy and the improvement claim be addressed in revision. It is neither a desk reject nor a clean accept as-is.","headline":"Solid, transparent solar-model study with real first constraints; the mass limits share a degeneracy with the assumed inner boundary radius, so the headline numbers are indicative rather than precise.","tokens_in":21377,"tokens_out":3367,"would_cite":true,"duration_ms":31058,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that current p-mode helioseismology constrains any compact dark core in the Sun to masses below about $10^{-5}$ solar masses—a thousand times tighter than neutrino limits—and that a $10^{-3}$ solar-mass core would…","keywords":["dark matter","compact dark objects","solar core","helioseismology","solar neutrinos","p modes","g modes","solar models"],"falsifier":"Measure the Sun's dipole g-mode period spacings with the precision now achieved for p modes: the standard model predicts approximately uniform spacings of about 25 minutes, a $10^{-3}~\\mathrm{M}_\\odot$ dark core predicts spacings near 1 minute, and a $10^{-7}~\\mathrm{M}_\\odot$ core perturbs the spacings by about 2%. Observing a standard 25-minute pattern across many periods would exclude dark cores above about $10^{-7}~\\mathrm{M}_\\odot$, while a compressed spacing would confirm a massive core; alternatively, a targeted search for the predicted high-inertia mixed modes in the solar oscillation spectrum would, if they are absent at observable amplitudes, rule out dark cores above about $10^{-5}~\\mathrm{M}_\\odot$.","tokens_in":20339,"feed_emoji":"☀️","tokens_out":12088,"duration_ms":95510,"temperature":0.7,"pith_summary":"The paper asks whether the Sun could harbor a compact dark object at its center—something like strange quark matter or a dark-sector macro—that interacts with ordinary matter only through gravity. It calibrates solar evolution models with a central point mass ranging from $10^{-8}$ to $10^{-2}$ solar masses and compares them with neutrino fluxes and with the acoustic (p-mode) and future gravity (g-mode) oscillation spectra. The result is a hierarchy of constraints: neutrinos rule out only the most massive cores, around 1% of the Sun's mass; current p-mode helioseismology already pushes the limit to about $10^{-5}$ solar masses; and future g-mode period spacings could reach $10^{-7}$ solar masses. A model with a $10^{-3}$ solar-mass dark core substantially improves agreement with the observed oscillation frequencies ($\\chi^2_r = 308$ versus 2633 for the standard solar model), but the authors interpret this as the dark core emulating a metal-rich core caused by enhanced gravitational settling, not as evidence for dark matter.","feed_headline":"Sun's p-mode hum caps any dark core at 10^-5 solar masses","feed_subtitle":"Neutrinos only veto cores near 1% of the Sun's mass; the acoustic spectrum is a thousand times sharper.","key_machinery":"The central object is a point-like dark core of mass $M_0$ at the center of an otherwise standard solar evolution model, added to the hydrostatic equilibrium equation as an extra gravitational term $-GM_0\\rho(r)/r^2$. The integration is truncated at the Bondi radius $r_B = 2GM_0/c_s^2$, whose value in solar radii is numerically close to the dark core mass in solar masses; the largest model considered, $10^{-2}~\\mathrm{M}_\\odot$, therefore removes the innermost 1% of the solar radius from the computational domain. This inner boundary produces a density cusp, drives rapid gravitational settling of heavy elements toward the center, raises the mean molecular weight, and thereby changes the sound speed, the buoyancy frequency, and the p- and g-mode oscillation frequencies that are compared with observations.","core_discovery":"On the paper's own terms, the central result is that a non-accreting, non-luminous compact dark core in the present Sun would be nearly invisible to neutrino detectors but visible to seismology: the absence of measured p-mode frequency shifts excludes dark cores above about $10^{-5}~\\mathrm{M}_\\odot$, neutrino fluxes only exclude cores at or above about $10^{-2}~\\mathrm{M}_\\odot$, and a model with a $10^{-3}~\\mathrm{M}_\\odot$ core fits the helioseismic sound-speed profile and radial-mode frequencies better than the standard solar model. The paper argues that this improvement is driven not by dark matter per se but by the heavy metal core that the strong central gravity creates through enhanced settling of magnesium, oxygen, and neon, and it therefore suggests the dark-core construction may be emulating star-formation effects that are usually neglected in solar models. The models also predict high-inertia mixed modes and a greatly compressed g-mode period spacing, which would observationally distinguish a genuine dark core from a metal-rich core.","pith_inferences":["If the $10^{-3}$ model's improvement is truly a proxy for a metal-rich solar core, the same dark-core construction could be used as a numerical shortcut to generate testable asteroseismic predictions for other main-sequence stars, regardless of whether dark matter is involved.","Extending the calculation to a slowly accreting, time-variable core would map how the mass limits shift and would connect these solar constraints to the primordial-black-hole limits from earlier work, since the present non-accreting assumption deliberately sets that physics aside.","A coarse future g-mode detection may separate the dark-core scenario from the metal-rich-core scenario even before high precision is reached, because the predicted period spacing differs by roughly a factor of 25, far larger than most modeling uncertainties."],"forward_implications":["Current p-mode frequencies and the inferred sound-speed profile already push any compact dark core in the Sun below about $10^{-5}$ solar masses, a factor of a thousand tighter than the neutrino limit.","A model with a $10^{-3}$ solar-mass dark core matches the helioseismic data substantially better than the standard solar model ($\\chi^2_r = 308$ versus 2633) and places the convection-zone base within $1\\sigma$ of its observed value, but the paper attributes this to the metal-rich core, not to dark matter.","The predicted high-inertia mixed modes would be observable signatures of dark cores above about $10^{-5}$ solar masses; their absence in the Sun and in solar-like oscillators generally suggests such cores are absent or rare.","Future solar g-mode measurements, if they reach the precision of current p-mode data, could constrain dark cores down to about $10^{-7}$ solar masses, with a $10^{-3}$ core shortening the period spacing from roughly 25 minutes to about 1 minute."],"supporting_citations":[{"why":"Provides the Bondi radius formula used to set the inner boundary of the solar models.","marker":"Bondi 1952"},{"why":"Supplies the stellar evolution code used to construct and calibrate the dark-core solar models.","marker":"Paxton et al. 2011"},{"why":"Supplies the oscillation code used to compute the p- and g-mode frequencies.","marker":"Townsend & Teitler 2013"},{"why":"Provides the helioseismically inferred solar sound-speed profile against which the models are compared.","marker":"Basu et al. 2009"},{"why":"Provides the measured p-mode frequencies and their uncertainties (BiSON) used to compute $\\chi^2_r$.","marker":"Davies et al. 2014"},{"why":"Provides the observed pp-neutrino flux used to set the neutrino constraint.","marker":"Bergström et al. 2016"},{"why":"Provides the measured CNO-neutrino flux used to set the neutrino constraint.","marker":"Basilico et al. 2023"},{"why":"Proposed the metal-rich solar core hypothesis that the dark-core model may be emulating.","marker":"Zhang et al. 2019"},{"why":"Proposed a formation scenario for a metal-rich solar core that the dark-core model may mimic.","marker":"Kunitomo & Guillot 2021"}],"fun_headline_variants":["Solar hum vetoes dark cores heavier than 10^-5 suns","Neutrinos let dark cores slide; Sun's song pins them down","Sun's hum squeezes dark core mass to 10^-5 M_sun","Dark core emulation: metal settling explains solar hum","Seismology beats neutrinos: dark core limit drops 1000x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the dark core is non-accreting and non-luminous and interacts with the surrounding solar plasma only gravitationally, with the evolution truncated at a fixed Bondi radius; if a real compact dark object accretes or radiates, or if the time-varying Bondi radius (which the authors note can change by about an order of magnitude) shifts the effective core size, the derived mass limits and the apparent improvement at $10^{-3}~\\mathrm{M}_\\odot$ would not apply as stated.","fun_headline_variants_meta":{"raw":{"variants":["Solar hum vetoes dark cores heavier than 10^-5 suns","Neutrinos let dark cores slide; Sun's song pins them down","Sun's hum squeezes dark core mass to 10^-5 M_sun","Dark core emulation: metal settling explains solar hum","Seismology beats neutrinos: dark core limit drops 1000x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000831,"raw_usage":{"total_tokens":3674,"prompt_tokens":1037,"completion_tokens":2637,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":2542}},"tokens_in":653,"tokens_out":2637,"duration_ms":19271,"temperature":1.0,"reasoning_tokens":2542,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:17:42.763149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Sun's dipole g-mode period spacings with the precision now achieved for p modes: the standard model predicts approximately uniform spacings of about 25 minutes, a $10^{-3}~\\mathrm{M}_\\odot$ dark core predicts spacings near 1 minute, and a $10^{-7}~\\mathrm{M}_\\odot$ core perturbs the spacings by about 2%. Observing a standard 25-minute pattern across many periods would exclude dark cores above about $10^{-7}~\\mathrm{M}_\\odot$, while a compressed spacing would confirm a massive core; alternatively, a targeted search for the predicted high-inertia mixed modes in the solar oscillation spectrum would, if they are absent at observable amplitudes, rule out dark cores above about $10^{-5}~\\mathrm{M}_\\odot$.","supporting_citations":[],"review_version":1}