{"id":"528976b3-22c9-4704-8a1b-aa739bd80831","arxiv_id":"2607.04744","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"R-squared gravity enlarges the equilibrium domain of rotating fermion-boson stars and raises static and Keplerian maximum masses relative to GR while remaining compatible with current compact-object constraints.","lead":"The authors build the first static and rotating mixed neutron-star–boson-star models in R-squared gravity and show that the extra scalar force raises maximum masses and widens the allowed mass–radius band relative to general relativity. The models stay compatible with NICER and gravitational-wave bounds and can ease the tension with the heavy secondary in GW190814.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"RPV truncation of the bosonic tail is the load-bearing assumption for the reported mass-radius bands and maximum-mass gains.","rationale":"The paper’s strongest claim is a clean numerical result: the scalaron systematically raises static and Keplerian maximum masses and enlarges the admissible (M,R) domain while remaining compatible with NICER/GW bounds. That claim is supported by the sequences that are actually computed. The only place where the claim can fail without any coding error is if the RPV truncation systematically misrepresents the bosonic contribution for the very configurations that fill the high-mass, high-¯w parts of Fig. 4. The reader already identified this as the weakest assumption; I agree it is also the single most load-bearing one. No other internal inconsistency (coordinate choice, conformal-frame mapping, or single-EOS choice) threatens the claim as directly. Because the paper is transparent about the approximation and the equilibrium catalog remains useful even if the absolute numbers shift modestly, the CONDITIONAL verdict with medium correctness risk is appropriate and needs no adjustment. The concrete test above would settle the issue with one targeted full-KG recomputation.","tokens_in":34312,"tokens_out":694,"duration_ms":6417,"concrete_test":"Recompute a representative high-¯w halo model (e.g., ¯w=0.09, M_T=2 M⊙, J_T=1 M⊙^{2}, a=10^4) by solving the full complex Klein-Gordon equation (Eq. 23) with the same metric potentials, retaining the exponential tail out to several times R_e, and re-evaluate the Komar integrals (Eqs. 42–43). If |ΔM_b|/M_b or |ΔJ_b|/J_b exceeds ~5–10 %, the mass-radius bands and maximum-mass gains of Fig. 4 must be re-assessed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (enlarged equilibrium domain and higher max masses relative to GR, still observationally viable) rests on the RPV high-coupling approximation of Sec. III A (Eqs. 29–33). The bosonic field is set algebraically by the max{0, ...} expression, the exponentially decaying tail is discarded, and T^b_μν is replaced by a compact-support perfect fluid. No independent eigenvalue problem for w is solved. The paper itself notes (Sec. VI) that the neglected tail can become non-negligible for halo-dominated configurations (¯w ≳ 0.08–0.09), which populate a substantial fraction of the colored bands in Fig. 4 and the high-¯w columns of Figs. 1–3. If the truncated mass/angular-momentum integrals (Eqs. 42–43) under- or over-estimate the true ADM quantities once the tail is restored, the quoted static/Keplerian maxima (2.44→2.62 M⊙ and 2.89→3.29 M⊙ for a=10^4) and the claimed compatibility with GW190814/NICER are not reliable. The reader correctly flags this as the weakest assumption; it is also the single most load-bearing one for the headline claim.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper constructs static and uniformly rotating equilibrium configurations of mixed fermion-boson stars in R-squared gravity, f(R)=R+aR², formulated in the Einstein frame as a scalar-tensor theory. The fermionic sector uses the tabulated AkmalPR EOS; the bosonic sector is a self-interacting complex field treated in the RPV high-coupling approximation, so that its energy-momentum tensor is replaced by an effective perfect fluid of compact support. Equilibrium models are obtained with a modified RNS/KEH code for GR and two representative values a=10 and a=10^4. The main claims are that the scalaron modifies the spatial distribution of both components, enlarges the domain of admissible equilibria, and raises static and Keplerian maximum masses relative to GR (e.g. static 2.44→2.62 M⊙ and Keplerian 2.89→3.29 M⊙ for a=10^4), while the sequences remain compatible with NICER, GW170817, and GW190814 constraints for bosonic mass fractions ≲10%.","tokens_in":34637,"tokens_out":1676,"duration_ms":21588,"significance":"This is, to my knowledge, the first self-consistent construction that combines a mixed fermion-boson matter sector, uniform rotation, and a viable f(R) model with a massive scalaron. Prior work treated rotating fermion-boson stars in GR or static mixed stars in f(R), but not the three ingredients together. The explicit Einstein-frame field equations (Appendices A–B), the transparent RPV implementation, and the mass–radius maps with external multimessenger bounds make the contribution useful for the strong-field modified-gravity and dark-matter-admixed compact-star communities. If the reported mass gains and enlarged solution space survive beyond the RPV truncation, the models offer a concrete, observationally testable channel for objects near the low-mass black-hole gap without requiring an extremely stiff nuclear EOS.","major_comments":[{"comment":"Section III A, Eqs. (29)–(33) and the mass/angular-momentum integrals (42)–(43): the headline mass–radius bands and maximum-mass gains in Fig. 4 and Sec. V C rest on the RPV truncation of the bosonic tail and the algebraic replacement of T^b_μν by a compact-support perfect fluid. The paper itself notes (Sec. VI) that the neglected tail can become non-negligible for halo-dominated configurations (¯w ≳ 0.08–0.09), which populate a substantial fraction of the colored regions in Fig. 4 and the high-¯w columns of Figs. 1–3. Without a quantitative estimate of the truncated ADM mass/angular momentum for the adopted (μ_b, λ, ¯w, η) values—or an explicit restriction of the observational claims to the compact/core-like regime where RPV is controlled—the quoted static/Keplerian maxima and the GW190814 compatibility statement are not yet fully reliable. A short validation (e.g. order-of-magnitude ta","section":"Section III A, Eqs. (29)–(33); Fig. 4; Sec. V C"},{"comment":"Section V C and the GW190814 discussion: the claim that a=10^4 allows a near-static mixed star with M_b/M_T ≲ 10% to sit inside the GW A band relies on the RPV sequences and on the single nuclear EOS AkmalPR. The paper correctly notes that stiff EOSs can reach ~2.6 M⊙ in GR, but the quantitative advantage attributed to R-squared gravity is not separated from (i) the RPV truncation and (ii) the fixed 10% bosonic-fraction ceiling. At minimum, the text should state how sensitive the static maximum (2.62 M⊙) is to modest changes in the bosonic fraction ceiling and to restoring a non-zero tail, and should avoid presenting the GW190814 resolution as robust until that sensitivity is shown.","section":"Section V C; Fig. 4"},{"comment":"Units and physical scale of a: results are reported for a=10 and a=10^4 with no explicit statement of the unit system for a (length² in geometrized units). Prior R-squared NS literature typically quotes a in km² or relative to a gravitational radius. Without this, the claimed “near-saturation / Brans–Dicke-like” regime for a=10^4 cannot be mapped to a scalaron mass or to Solar-System/strong-field bounds. Please specify the units of a consistently with the dimensionless code variables and with the literature values used for comparison.","section":"Section V; comparison with Refs. [11, 16]"}],"minor_comments":[{"comment":"Abstract and Introduction: “enlarges the domain of admissible equilibrium solutions” is used without a precise definition (parameter volume, mass–radius area, or existence of new topologies). A one-sentence operational definition would help.","section":"Abstract; Sec. I"},{"comment":"Eq. (36) and surrounding text: the compactified coordinate uses s, and μ=cos θ is introduced with a footnote warning against confusion with μ_b; still, several source-term symbols (S^x_T) reuse subscripts that collide with total quantities. A short notation table would reduce ambiguity.","section":"Section III B"},{"comment":"Fig. 1–3: relative-percentage panels are useful, but the vertical scales differ across columns; stating the peak percentage in the caption (as done for δφ²_max in Fig. 2) for the fermionic panels would aid comparison.","section":"Figs. 1–3"},{"comment":"Appendix A: the static metric is Schwarzschild-like while the rotating metric is quasi-isotropic; the text already warns that the static equations are not a direct coordinate limit of Eq. (24). Cross-referencing that warning again when static profiles are extracted as J_T→0 of the rotating code would avoid reader confusion.","section":"Appendix A; Sec. V A"},{"comment":"Typographical: “dimentionless” (Sec. V B); “taht” in the Introduction’s appendix description; “inequivalent definitions” list is fine but “axion couplings [25]” could be checked for consistency with the reference list numbering after any revision.","section":"Sec. V B; Introduction"},{"comment":"Bosonic microphysical parameters (μ_b ≃ 10^{-16} MeV, λ=100) are stated once in Sec. III B; a brief remark on why this ultralight scale is appropriate for the dark-matter interpretation (or that it is chosen purely to sit in the RPV regime) would help non-specialist readers.","section":"Section III B"}],"recommendation":"major_revision","confidential_remarks":"The work is a natural and well-executed extension of the authors’ prior GR rotating fermion-boson-star paper and of existing R-squared NS literature. Novelty is real but incremental; the main risk for the journal is over-claiming multimessenger compatibility while the RPV truncation remains unquantified for the halo part of the solution space. If the authors supply a tail-mass estimate or restrict the observational discussion, the paper should be publishable. No concerns about citation pattern or ethics."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is the first self-consistent construction of uniformly rotating fermion-boson stars in R-squared gravity. That is the real news. Static mixed stars in f(R) and rotating mixed stars in GR already exist; putting rotation, a complex bosonic sector, and the scalaron together in one KEH/RNS pipeline is new and useful for the subfield.\n\nThey do the technical work carefully. Einstein-frame equations, conformal factors, Komar integrals, and the source terms for the elliptic system are written out. Fixed-(M_T, J_T, w̄) profiles and free-parameter mass-radius maps for GR, a=10, and a=10^4 are clear, with percentage deviations and the usual NICER/GW overlays. The headline numbers (static max mass 2.44\to2.62 M⊙, Keplerian 2.89\to3.29 M⊙ for a=10^4) come from actual sequences, not from a normalization that already encodes the answer. Citation pattern is honest about prior GR rotating FBS and static f(R) mixed work.\n\nThe soft spot is the RPV high-coupling truncation: the bosonic tail is discarded, T^b is replaced by a compact-support perfect fluid, and w is not an eigenvalue. The paper itself notes that the tail can matter for halo-like models (higher w̄), which occupy part of the colored bands. That is a genuine limitation on the absolute mass-radius numbers and on the GW190814 claim for those configurations; it is not a hidden flaw. No stability analysis, one nuclear EOS, and no public code are the usual next-step gaps, not load-bearing cracks.\n\nWho it is for: people who build or use mixed compact-object models in modified gravity and need a concrete rotating catalog. A serious referee should see it. I would engage with the equilibria and cite the sequences when I need a rotating mixed baseline in R^{2}; I would not treat the high-w̄ halo masses as final until the tail is restored.","headline":"First rotating fermion-boson stars in R^{2} gravity; solid numerical catalog with a real but well-flagged RPV limitation.","tokens_in":35293,"tokens_out":499,"would_cite":true,"duration_ms":5793,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"R-squared gravity plus a dark bosonic component raises the maximum mass of rotating mixed stars and widens the allowed mass-radius band relative to general relativity, while remaining compatible with current NICER and gravitational-wave bou","keywords":["fermion-boson stars","R-squared gravity","rotating compact stars","scalar degree of freedom","dark matter admixed neutron stars","mass-radius relation","GW190814"],"falsifier":"A numerical-relativity evolution of a high-mass, low-frequency model that retains the full bosonic tail would show whether the configuration remains long-lived or disperses; if it disperses while the truncated-tail equilibrium predicts stability, the maximum-mass claims fail.","tokens_in":35173,"feed_emoji":"⭐","tokens_out":653,"duration_ms":5921,"temperature":0.7,"pith_summary":"The paper constructs static and uniformly rotating equilibrium stars that contain both ordinary nuclear matter and a self-interacting complex bosonic field, all inside the simplest quadratic modification of Einstein gravity, f(R)=R+aR^{2}. The extra scalar degree of freedom that this theory carries redistributes both the bosonic amplitude and the fermionic pressure, systematically enlarges the space of allowed solutions, and lifts the maximum supported mass (static and mass-shedding) above the general-relativity values. For the largest a considered, the static maximum reaches 2.62 solar masses and the Keplerian maximum 3.29 solar masses. Because these sequences still intersect the NICER mass-radius ellipses and the secondary-mass window of GW190814, the models offer a concrete way to accommodate compact objects near the low-mass black-hole gap without forcing an unrealistically stiff nuclear equation of state.","feed_headline":"R-squared gravity lifts mixed-star mass limits past 3 solar masses","feed_subtitle":"Scalar corrections plus a dark bosonic component widen the mass-radius band while still matching NICER and GW bounds","key_machinery":"The RPV high-coupling approximation that replaces the complex bosonic field by an effective perfect-fluid source of compact support, solved together with the Einstein-frame scalar-tensor equations of R-squared gravity inside a modified RNS self-consistent-field code.","core_discovery":"In R-squared gravity the scalar degree of freedom modifies the spatial profiles of both the bosonic field and the fermionic energy density, expands the domain of admissible mixed-star equilibria, and raises both the static and Keplerian maximum masses relative to pure general relativity, all while the resulting sequences remain compatible with current astrophysical and gravitational-wave constraints.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["R² gravity lifts fermion-boson star max masses above GR limits","Scalar mode expands mixed-star equilibria and raises mass caps","Rotating fermion-boson stars support higher masses in R² gravity","R-squared gravity boosts static and Keplerian mixed-star masses","Dark boson plus R² scalar widens admissible star mass sequences"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The calculation discards the exponentially decaying outer tail of the bosonic field and treats the bosonic sector as a perfect fluid of compact support; if that tail or the strong-self-interaction assumption is not valid, the reported mass-radius bands are unreliable.","fun_headline_variants_meta":{"raw":{"variants":["R² gravity lifts fermion-boson star max masses above GR limits","Scalar mode expands mixed-star equilibria and raises mass caps","Rotating fermion-boson stars support higher masses in R² gravity","R-squared gravity boosts static and Keplerian mixed-star masses","Dark boson plus R² scalar widens admissible star mass sequences"]},"model":"grok-4.5","effort":"low","cost_usd":0.003568,"raw_usage":{"total_tokens":1126,"prompt_tokens":797,"num_sources_used":0,"completion_tokens":93,"cost_in_usd_ticks":35680000,"prompt_tokens_details":{"text_tokens":797,"audio_tokens":0,"image_tokens":0,"cached_tokens":0},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":236,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":797,"tokens_out":93,"duration_ms":2621,"temperature":1.0,"reasoning_tokens":236,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T14:12:21.788798+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A numerical-relativity evolution of a high-mass, low-frequency model that retains the full bosonic tail would show whether the configuration remains long-lived or disperses; if it disperses while the truncated-tail equilibrium predicts stability, the maximum-mass claims fail.","supporting_citations":[],"review_version":1}