{"id":"8faf92e6-d1a1-44f0-82a0-1012a15187bb","arxiv_id":"2505.11591","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Dust can form stable horizonless compact objects in the MEMe modified gravity model, offering a dark matter candidate and a mechanism that suppresses black hole formation below about 10^-11 solar masses.","lead":"This paper shows that a modified gravity theory called MEMe can support clouds of ordinary dust as stable, horizonless 'dust stars' instead of collapsing into black holes. If this is right, these objects could form part of the dark matter, and the theory would also suppress the formation of very small black holes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability analysis for the only viable dust-star class is explicitly heuristic; a missed negative mode below the scanned range would overturn the dark-matter claim.","rationale":"The paper's central astrophysical claim is not merely the formal existence of some TOV solution; it is that positive-q dust stars form a stable population of horizonless MACHOs that can constitute a significant dark-matter fraction. This requires the first-region solutions (rho_hat_0 < 2/q) to be radially stable. The authors are honest that the radial-stability framework is not established at r=0 and explicitly defer a rigorous treatment; the numerical eigenvalues in Table I are sparse, and the shooting search starts at a finite negative value. Those two facts mean the statement that all first-region solutions are stable is currently a conjecture supported by one representative eigenvalue calculation rather than a demonstrated result. I do not see a comparably load-bearing flaw in the density-bound and M_min derivation: Eqs. (21)-(23) are internally consistent, Eq. (26) follows from them, and the Buchdahl argument is explicitly labeled as rough. The notational slip in Eq. (24) does not change the numerical factor actually used. The physical-frame ambiguity and the q<0 viability question are worth following up but are secondary to the stability of the proposed dark-matter candidate. The correct disposition is therefore to keep the paper conditional, pending a regularized spectral calculation or an application of the Luz-Carloni formalism to the first-region sequence.","tokens_in":15336,"tokens_out":26919,"duration_ms":299102,"concrete_test":"Reproduce the radial-perturbation spectrum for rho_hat_0 = 1/2q, 1/q, 3/2q, 9/5q, and 199/100q with a spectral method in the regularized coordinate x=r^2 (or the Luz-Carloni formalism of refs. [40-42]), using a search window extending to omega^2 = -10^8 s^-2 and checking convergence under grid refinement. If any omega_n^2 < 0 appears for rho_hat_0 < 2/q, or if the spectrum is not discrete, the stability-based division into a stable first class and unstable higher classes fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing point is the stability of the first-region dust stars (rho_hat_0 < 2/q), because only this class can serve as the claimed stable dark-matter MACHOs. Section VII itself flags that Eq. (30) is singular at r=0 and that the discreteness of the Sturm-Liouville spectrum is assumed, calling the analysis heuristic and in need of a more rigorous approach [40-42]. The numerical part does not close that gap: eigenvalues are shown for only rho_hat_0 = 1/q, 5/q, and 10/q, and the shooting search starts at omega^2 = -10^6 s^-2; a negative mode below that starting value would not be seen. If the singular operator has a continuous unstable sector, or if the effective Einstein-frame EoS does not correctly describe adiabatic perturbations of the Jordan-frame dust, then Table I's positive fundamental mode for one solution does not establish stability of the entire first region. Since removing first-region stability eliminates the only physically viable candidate, this unresolved assumption is the most load-bearing weakness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies static, spherically symmetric compact objects in the Minimal Exponential Measure (MEMe) model, a Type I minimally modified gravity theory with a single free parameter q. For a Jordan-frame perfect fluid with a linear equation of state, the authors rewrite the field equations in an Einstein frame and derive the effective Einstein-frame equation of state. For positive q and Jordan-frame dust, they obtain an implicit Einstein-frame EoS with three branches, prove an upper bound rho_max = 7/(9q) on the Einstein-frame density, and numerically integrate the TOV equations to find horizonless compact objects. They classify solutions by central density, show that only the first branch (rho_hat_0 < 2/q) avoids a negative fundamental radial mode, and derive a minimum black hole mass from a Buchdahl-type argument. They then use constraints on primordial black holes and microlensing to identify a window of q for which these dust stars could be a significant part of dark matter. The paper explicitly labels the radial stability analysis as heuristic because of a singular point at r=0 in the perturbation equations.","tokens_in":15486,"tokens_out":7947,"duration_ms":78101,"significance":"If the results hold, the paper reports a genuinely new mechanism: a minimally modified gravity theory in which pressureless dust can support horizonless compact objects, with masses below roughly 10^-11 solar masses and a parameter range that evades existing MACHO constraints. The analytic derivation of the Einstein-frame EoS (Eq. 20), the density bound (Eqs. 21-23), and the M_min formula (Eq. 26) are clean and transparent. The single parameter q is left free and later constrained by external physical inputs, so the argument is not circular. The paper is also honest about the main weakness, flagging the singular Sturm-Liouville issue in Section VII. However, the dark-matter claim depends on the stability of the first-branch solutions, and that stability is not yet established with the rigor the central claim requires.","major_comments":[{"comment":"The claim that all first-region dust stars (rho_hat_0 < 2/q) are stable is load-bearing for the dark-matter interpretation, but the supporting analysis is explicitly heuristic. The singular point at r=0 in Eq. (30) is assumed not to spoil the discreteness of the Sturm-Liouville spectrum, with no argument or reference supplying that proof. The numerical search starts at omega^2 = -10^6 s^-2, so any unstable mode with a more negative omega^2 would not be seen. Only three representative solutions are shown in Table I, and the statement that the check was performed for all first-region solutions is not accompanied by a scan, a stability criterion, or a plot. Because removing first-region stability eliminates the only viable MACHO candidate, this gap must be closed or the abstract and Section VIII claims must be substantially weakened.","section":"Sec. VII, Eqs. (30)-(31), Table I"},{"comment":"The paper claims that a positive q suppresses the formation of primordial black holes below M_min, but this inference applies a static, equilibrium Buchdahl bound to a dynamical collapse process. The density maximum rho_max = 7/(9q) was derived for static TOV solutions with a monotonically increasing Einstein-frame EoS; it is not shown that the same bound controls the densities reached during gravitational collapse, nor that the Buchdahl divergence of central pressure implies collapse in the presence of the modified matter coupling. A dynamical collapse calculation, or at least a clearly stated additional assumption, is needed to support the abstract statement that the model provides a mechanism for suppressing PBH formation.","section":"Sec. V, Eqs. (24)-(26)"},{"comment":"The adiabatic index used in the perturbation equations is computed from the effective Einstein-frame EoS, but the paper does not derive the perturbation mapping between the Jordan-frame dust, which has zero pressure, and the Einstein-frame perfect fluid. The statement that the unmodified radial oscillation equations apply in the Einstein frame requires that the effective fluid's perturbations faithfully represent perturbations of the underlying physical matter in the MEMe model. Without such a derivation, the stability conclusions for the first branch remain an assumption about the effective description rather than a demonstrated property of the Jordan-frame dust stars.","section":"Sec. VII, Eq. (32)"}],"minor_comments":[{"comment":"The displayed definition rho_c := M/(4 pi R_Buch^3) is inconsistent with the numerical value 16/(243 pi M^2); the uniform-density relation is rho_c = 3M/(4 pi R_Buch^3), which yields the stated value.","section":"Sec. V, Eq. (24)"},{"comment":"The AM-GM step would be easier to verify if the identity 4 - q(3 p - rho) = 3(1 - q p) + (1 + q rho) were stated before applying the inequality.","section":"Sec. IV, Eqs. (21)-(23)"},{"comment":"There is a typo, 'an relatively simple', which should read 'a relatively simple'.","section":"Sec. VIII"},{"comment":"The phrase 'smaller from the minimum black hole mass' should read 'smaller than the minimum black hole mass'.","section":"Fig. 7 and Sec. VIII"},{"comment":"The bound 1/|q| lesssim 10^33 GeV/fm^3 is explicitly conditional on attributing the absence of PBH evaporation signals to the lower mass bound; the paper should state more prominently that this is an assumption, not a direct observational upper limit on q.","section":"Sec. V, Eq. (27)"}],"recommendation":"major_revision","confidential_remarks":"The authors are to be commended for explicitly flagging the heuristic nature of the stability analysis and for providing a clearly written derivation of the equation of state and density bound. The main issue for publication is the gap between the strong claims in the abstract and Section VIII and the unresolved stability assumption in Section VII. I would ask either for a rigorous radial-stability treatment along the lines of the cited Luz-Carloni work, or for a substantial softening of the MACHO and PBH-suppression claims so that they are presented as conditional on the stability assumption. The paper is otherwise technically sound and of interest to the gr-qc community."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, genuinely new result in a niche theory, with one load-bearing caveat that the authors themselves flag. The dust star solutions exist; whether they are stable enough to be dark matter is not settled.\n\nWhat's new: first detailed TOV study of spherically symmetric solutions in the MEMe model, the positive-q branch, the transformed Einstein-frame EoS for dust, the density bound 7/(9q) from an AM-GM argument, and the minimum black hole mass. The analytic derivation in Sec. IV is clean, and the bound is rigorous for any monotonic physical EoS. The three-region classification of the EoS and the M-R curves are useful. Credit where due: the paper does not oversell its stability analysis; it explicitly calls it heuristic and cites the Luz-Carloni work.\n\nSoft spots, in proportion: (1) The radial stability analysis is the biggest gap. Eq. (30) is singular at r=0, and the authors assume the Sturm-Liouville spectrum stays discrete. The numerics show eigenvalues for only three central densities, and the assertion that all first-region solutions are stable is unaccompanied by data. A missed negative mode in the first region would remove the viable dark matter candidate. This is a real load-bearing weakness, but the paper already labels the analysis as heuristic, so readers are warned. (2) The PBH suppression mechanism uses a uniform-density Buchdahl estimate, not a collapse simulation; that is fine for a rough bound, but the abstract's wording is a bit stronger than the argument supports. (3) Eq. (24) has a notational slip: it writes M/(4πR^3) where the subsequent formula uses 3M/(4πR^3). Minor. (4) No code or data are shipped, which makes the 'all first region' stability claim hard to audit.\n\nThe dark matter window is conditional on the single free parameter q, and the paper says so. The existence claim does not depend on the stability analysis, so even if the stability question remains open, the core result stands.\n\nWho this is for: people working on minimally modified gravity theories and alternative compact objects. It deserves a serious referee. The heuristic stability section should be flagged for revision, but this is not a desk-reject paper. I'd send it to review, ask for the stability analysis to be tightened or its range of validity honestly narrowed, and request the code or at least fuller eigenvalue tables.","headline":"Genuinely new dust-supported horizonless compact objects in the MEMe model, with a clean analytic core and an explicitly heuristic stability analysis that is the main thing to fix before the dark matter claim hardens.","tokens_in":16068,"tokens_out":3480,"would_cite":false,"duration_ms":32692,"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":"The paper claims that in the MEMe model of gravity, dust can form horizonless compact objects with masses below about $10^{-11}$ solar masses, making them viable MACHO dark-matter candidates, and that this requires the model's single free…","keywords":["minimal exponential measure","modified gravity","dust stars","compact objects","MACHOs","dark matter","primordial black holes","equation of state"],"falsifier":"A rigorous radial-perturbation calculation (treating the $r=0$ singular term in Eq. (30) with a proper boundary condition, or a nonlinear simulation) that finds a growing mode for any first-branch solution with $\\hat\\rho_0<2/q$ would falsify the stability claim; alternatively, constructing a static horizonless dust solution with $q<0$ would falsify the claimed necessary condition $q>0$.","tokens_in":15058,"feed_emoji":"⭐","tokens_out":15982,"duration_ms":142561,"temperature":0.7,"pith_summary":"The paper tries to establish that the Minimal Exponential Measure (MEMe) model, a modification of gravity that changes how matter couples to geometry without adding new dynamical fields, admits horizonless compact objects supported by dust. In the model's Einstein-frame description, pressureless dust acquires a nonlinear effective pressure, so the stellar-structure equilibrium equations have finite-mass solutions instead of the infinite-mass solutions that dust-like equations of state produce in general relativity. For positive values of the model's single free parameter $q$, these dust-star solutions exist with masses below about $10^{-11}\\,M_\\odot$, placing them in the MACHO window for dark matter. The same positive-$q$ density bound $\\rho_{\\max}=7/(9q)$ sets a minimum mass for black holes formed by collapse, so the model can suppress primordial black hole formation below that scale. Only the solutions with central physical density below $2/q$ are stable under radial perturbations, and the paper treats that stability result as heuristic.","feed_headline":"Dust clouds can form horizonless dark-matter stars","feed_subtitle":"In a minimal modified-gravity model, dust gains an effective pressure and forms stable MACHOs below 10^-11 solar masses.","key_machinery":"The load-bearing object is the exact algebraic map between the Jordan-frame fluid variables $(\\hat\\rho,\\hat p)$ and the Einstein-frame variables $(\\rho,p)$, given by Eqs. (9)-(10) through the determinant $|A|$ of the auxiliary tensor. For a perfect fluid this produces the implicit Einstein-frame equation of state $F(\\rho,p)=0$ (Eq. (20)). The argument then runs through the inequality $\\rho\\le 7/(9q)$, obtained by applying the arithmetic-geometric mean inequality to Eq. (21), and the Buchdahl-limit argument that converts the density cap into $M_{\\min}$. The stellar-structure system (28) is integrated in the Einstein frame, and stability is classified with the standard radial-perturbation equations (30)-(31).","core_discovery":"On the paper's own terms, the central discovery is that the MEMe field equations for a perfect fluid can be rewritten exactly as Einstein equations sourced by a transformed fluid, and that for dust the transformed equation of state is the implicit relation $F(\\rho,p)=0$ of Eq. (20). The transformed density is bounded above by $7/(9q)$, an inequality saturated only by dust, and this cap turns a would-be collapse into a horizonless compact object rather than a black hole. Numerically solving the stellar-structure equations in the Einstein frame gives three classes of solutions, of which only the first class (initial physical density $\\hat\\rho_0<2/q$) is stable; these stable dust stars have masses below about $10^{-11}\\,M_\\odot$ and radii controlled by the same energy scale $1/q$. The paper also derives the minimum black hole mass $M_{\\min}=(c^4/G^{3/2})(4/9)\\sqrt{3q/(7\\pi)}$ from the Buchdahl limit, and identifies $q>0$ as a necessary condition for all of this.","pith_inferences":["The paper leaves implicit that stable dust stars, being horizonless, could be distinguished from black holes by direct observation: gravitational-wave ringdowns, accretion shadows, or tidal signatures would differ even at the same mass, making the MACHO window a live target for next-generation surveys.","The $q>0$ requirement for dust stars is in tension with the bounce-cosmology scenario in earlier MEMe work, which needs $q<0$; the multi-MEMe extension proposed in the appendix could reconcile them, but within a single sector both mechanisms cannot operate. This is a consistency question a future calculation could settle.","A sharp prediction of the suppression mechanism is a mass cutoff in the black hole mass function; if future surveys find collapse-formed black holes below $M_{\\min}$ for an independently measured $q$, the mechanism is falsified rather than merely constrained."],"forward_implications":["The MEMe model predicts a new population of horizonless, dust-supported compact objects with masses up to about $10^{-11}\\,M_\\odot$ that can act as MACHO dark matter and evade current microlensing limits when $1/|q|\\gtrsim 10^{18}$ GeV/fm$^3$.","Black holes formed from the gravitational collapse of matter with a monotonic equation of state must be heavier than $M_{\\min}$, so primordial black hole formation below that mass is suppressed; this offers an alternative to other explanations for the lack of evaporation signals from low-mass PBHs.","The stable dust-star branch is the one with central physical density $\\hat\\rho_0<2/q$; the other branches, though ultra-compact and possibly smaller than their photon spheres, are unstable and should not be observable as enduring objects.","The upper density bound $\\rho_{\\max}=7/(9q)$ applies to every monotonically increasing physical equation of state, so the minimum-mass floor is a generic feature of the $q>0$ model rather than an artifact of the dust choice."],"supporting_citations":[{"why":"It defines the MEMe action and its perfect-fluid Einstein-frame reduction.","marker":"[14]"},{"why":"It shows the general-relativity obstruction that motivates the search: isothermal equations of state give infinite total mass.","marker":"[22]"},{"why":"It supplies primordial-black-hole evaporation constraints used to bound $1/q$.","marker":"[27]"},{"why":"It provides the compiled PBH/MACHO constraint curves used to identify the allowed mass ranges.","marker":"[33]"},{"why":"It is the numerical solver used to integrate the stellar-structure equations and produce the dust-star solutions.","marker":"[34]"},{"why":"It is the source of the radial-oscillation equations used for the stability analysis.","marker":"[37]"},{"why":"It completes the radial-perturbation framework used for the eigenvalue problem.","marker":"[39]"},{"why":"It is the rigorous radial-stability treatment the paper says would be needed to make the analysis conclusive.","marker":"[40]"},{"why":"It provides the microlensing observations that set the MACHO mass constraints the dust stars must evade.","marker":"[53]"}],"fun_headline_variants":["Dust forms horizonless MACHOs in MEMe gravity","MEMe dust stars: horizonless dark matter candidates","Modified gravity lets dust avoid black hole collapse","Horizonless dust stars as dark matter MACHOs","MEMe model: dust makes stable horizonless stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The radial-stability analysis assumes that the singular point at $r=0$ in the oscillation equation (30) does not break the discreteness of the Sturm-Liouville spectrum, so the standard perturbation framework applies; the paper labels this heuristic, and if the assumption fails the conclusion that only the $\\hat\\rho_0<2/q$ dust stars are stable collapses.","fun_headline_variants_meta":{"raw":{"variants":["Dust forms horizonless MACHOs in MEMe gravity","MEMe dust stars: horizonless dark matter candidates","Modified gravity lets dust avoid black hole collapse","Horizonless dust stars as dark matter MACHOs","MEMe model: dust makes stable horizonless stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001861,"raw_usage":{"total_tokens":7302,"prompt_tokens":934,"completion_tokens":6368,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":6291}},"tokens_in":550,"tokens_out":6368,"duration_ms":40358,"temperature":1.0,"reasoning_tokens":6291,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:51:17.703192+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A rigorous radial-perturbation calculation (treating the $r=0$ singular term in Eq. (30) with a proper boundary condition, or a nonlinear simulation) that finds a growing mode for any first-branch solution with $\\hat\\rho_0<2/q$ would falsify the stability claim; alternatively, constructing a static horizonless dust solution with $q<0$ would falsify the claimed necessary condition $q>0$.","supporting_citations":[{"cited_title":"Junction conditions and sharp gradients in generalized coupling theories","cited_arxiv_id":"2203.00011","evidence_quote":"It shows the general-relativity obstruction that motivates the search: isothermal equations of state give infinite total mass."},{"cited_title":"Chandrasekhar,An introduction to the study of stellar structure(Chicago Univ","cited_arxiv_id":null,"evidence_quote":"It supplies primordial-black-hole evaporation constraints used to bound $1/q$."},{"cited_title":"Chen and R","cited_arxiv_id":null,"evidence_quote":"It is the numerical solver used to integrate the stellar-structure equations and produce the dust-star solutions."},{"cited_title":"Chandrasekhar, Astrophys","cited_arxiv_id":null,"evidence_quote":"It completes the radial-perturbation framework used for the eigenvalue problem."},{"cited_title":"Chandrasekhar, Phys","cited_arxiv_id":null,"evidence_quote":"It is the rigorous radial-stability treatment the paper says would be needed to make the analysis conclusive."},{"cited_title":"Aharonov, A","cited_arxiv_id":null,"evidence_quote":"It provides the microlensing observations that set the MACHO mass constraints the dust stars must evade."}],"review_version":1}