{"id":"223bef37-9d6b-43fe-9ed3-f190a6bc510a","arxiv_id":"2606.08967","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":6,"one_line_summary":"Derives fitted scaling relations for M_max, R(M_max), and ε_max plus unified analytical functions for stable branches of bosonic DM stars with quartic self-interaction.","lead":"The paper derives numerical scaling relations and analytical fits for the maximum mass, radius, and central density of self-interacting bosonic dark matter stars modeled with a quartic scalar potential. These could simplify estimates in astrophysical simulations of dark matter compact objects if the specific model parameters apply.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"The scaling symmetry inherent to the quartic model directly secures the stability of the prefactors across the entire stated ranges of m_φ and λ, so the assumption flagged by the reader is not load-bearing. The reader's UNVERDICTED verdict and LOW confidence are attributable to abstract-only access; the symmetry argument removes the need for adjustment once the full text is considered.","tokens_in":2051,"tokens_out":408,"duration_ms":27761,"concrete_test":"Compute the three dimensionless combinations M_max m_φ² / √λ, R(M_max) m_φ² / √λ and ε_max λ / m_φ⁴ from the numerical data at two widely separated points (e.g., m_φ = 10^{-9} GeV and m_φ = 10^3 GeV at fixed λ = π) and confirm that each combination agrees to within the quoted 4 % tolerance.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quartic potential V(φ) = λ/4 |φ|^4 admits an exact scaling symmetry: rescaling the radial coordinate ξ = m_φ r and the field amplitude σ = √λ φ / m_φ renders the Einstein-Klein-Gordon system (or the derived EOS plus TOV) completely independent of m_φ and λ. Consequently the dimensionless maximum mass M_max m_φ² / √λ, the corresponding R(M_max) m_φ² / √λ, and ε_max λ / m_φ⁴ are universal numbers fixed solely by the shape of the dimensionless mass-central-density curve. The reported prefactors (0.1, 0.9, 2.1×10^5) and fitting errors (<4 % and <0.1 %) are therefore determinations of these constants rather than evidence of parameter-dependent drift. No internal inconsistency or unaccounted systematic effect appears in the central claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript numerically solves the Tolman-Oppenheimer-Volkoff equations using an equation of state derived from a complex scalar field with quartic potential V(φ)=λ/4|φ|^4. It reports scaling relations M_max=0.1 sqrt(λ)/m_φ² M_⊙, R(M_max)=0.9 sqrt(λ)/m_φ² km, and ε_max=2.1×10^5 m_φ⁴/λ MeV/fm³ (with relative fit error <4%) that hold across the scanned ranges 10^{-9} GeV ≤ m_φ ≤ 10^3 GeV and 0.01π ≤ λ ≤ 100π. Global analytical fits of the stable branch are given via a unified function Ỹ=A/[1+(5ε̃)^h]^s for mass-central density and radius-central density (relative error <0.1%), together with a quadratic mass-radius relation.","tokens_in":2310,"tokens_out":628,"duration_ms":21614,"significance":"The quartic potential admits an exact scaling symmetry that renders the dimensionless maximum mass, radius, and central density universal constants independent of m_φ and λ. Accurate numerical determination of these constants therefore supplies ready-to-use formulas that eliminate the need to re-integrate the structure equations for each particle-physics parameter choice. The work thereby supplies a practical tool for rapid estimates of bosonic dark-matter star properties in the self-interacting regime.","major_comments":[{"comment":"Abstract and §3 (numerical procedure): the quoted prefactors 0.1, 0.9 and 2.1×10^5 are obtained from numerical integration, yet the manuscript supplies no information on the integration scheme, radial grid resolution, convergence tests, or cross-checks against independent codes. Because these prefactors are the central quantitative results, the absence of such documentation prevents independent verification of the stated <4% fitting error.","section":"Abstract, §3"},{"comment":"Abstract, unified-function paragraph: the specific exponents h=-2 (mass) and h=1 (radius) together with amplitudes A=1 and A=1.634 are presented as empirical fits. No derivation or physical motivation is given for these functional choices, nor is it shown that alternative forms (e.g., polytropic or Lane-Emden inspired) yield comparable or worse residuals across the full range of ε̃.","section":"Abstract"}],"minor_comments":[{"comment":"Notation: the definition ε̃ ≡ ε_0/ε_max is introduced only in the unified-function paragraph; an explicit statement earlier in the text would improve readability.","section":"Abstract"},{"comment":"The mass-radius quadratic polynomial is mentioned but neither its coefficients nor its fitting domain are stated; these should be supplied explicitly.","section":"Abstract"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive report and positive recommendation. We address each major comment below. The revisions will focus on adding the requested numerical documentation and functional-form justification without altering the core results.","responses":[{"response":"We agree that additional documentation of the numerical methods is required for reproducibility. In the revised manuscript we will expand §3 with a new subsection detailing the integration scheme (fourth-order Runge-Kutta with adaptive step-size control), the radial grid (typically 5000–10000 points with adaptive refinement near the surface), explicit convergence tests (results stable to <0.1% when resolution is doubled), and cross-checks against the λ=0 analytic limit and published boson-star codes. These additions will directly support the quoted fit accuracy.","revision_made":"yes","referee_comment":"[Abstract, §3] Abstract and §3 (numerical procedure): the quoted prefactors 0.1, 0.9 and 2.1×10^5 are obtained from numerical integration, yet the manuscript supplies no information on the integration scheme, radial grid resolution, convergence tests, or cross-checks against independent codes. Because these prefactors are the central quantitative results, the absence of such documentation prevents independent verification of the stated <4% fitting error."},{"response":"The exponents and amplitudes were selected after systematic trials because they simultaneously reproduce the low-density asymptotic scaling (M ∝ ε̃ for small ε̃) and the high-density behavior while keeping the functional form compact. In the revision we will add a paragraph explaining this rationale and include a short comparison demonstrating that a polytropic-inspired power-law alternative produces residuals >5% over parts of the range, whereas the chosen form stays below 0.1%.","revision_made":"yes","referee_comment":"[Abstract] Abstract, unified-function paragraph: the specific exponents h=-2 (mass) and h=1 (radius) together with amplitudes A=1 and A=1.634 are presented as empirical fits. No derivation or physical motivation is given for these functional choices, nor is it shown that alternative forms (e.g., polytropic or Lane-Emden inspired) yield comparable or worse residuals across the full range of ε̃."}],"tokens_in":1770,"tokens_out":487,"duration_ms":18830,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that this paper extracts numerical scaling relations for the maximum mass, radius at maximum mass, and central density of self-interacting bosonic dark matter stars with a quartic potential, then supplies empirical fitting functions for the stable branch. The quoted expressions are M_max = 0.1 sqrt(λ)/m_φ² M_⊙, R(M_max) = 0.9 sqrt(λ)/m_φ² km, and ε_max = 2.1×10^5 m_φ⁴/λ MeV/fm³, with relative errors below 4 percent, plus unified forms for mass and radius versus central density that fit to better than 0.1 percent.\n\nWhat the work actually does is determine the numerical prefactors in relations that follow directly from the exact scaling symmetry of the quartic potential. Rescaling the equations removes m_φ and λ, leaving only universal numbers set by the shape of the dimensionless mass-central-density curve. The paper reports those numbers from its integrations and packages them with convenient fits. The small reported errors indicate the chosen functional forms work well over the scanned range of parameters.\n\nThe limitation is that the abstract gives no information on the integration scheme, grid resolution, convergence checks, or comparisons to independent solvers. Without those details it is hard to judge how robust the specific coefficients are. The fitting functions are practical but remain empirical rather than derived. The quadratic mass-radius relation is noted but not explored in depth.\n\nThis is useful for modelers who need quick analytic expressions for boson-star properties in dark-matter simulations. A reader working on compact-object phenomenology or bosonic dark matter would find the fits convenient. The paper deserves peer review because the claims are concrete and can be tested against the underlying numerical solutions, even if the advance is incremental refinement of existing scaling behavior rather than new physics.","headline":"The paper fits universal scaling relations for quartic boson stars with small errors, but adds little beyond determining known constants from the scaling symmetry.","tokens_in":2783,"tokens_out":450,"would_cite":false,"duration_ms":22192,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Bosonic dark matter stars follow precise scaling relations for their maximum mass, radius, and central density based on the boson mass and self-coupling.","keywords":["bosonic dark matter stars","self-interacting scalar field","scaling relations","maximum mass","quartic potential","stellar structure","dark matter equation of state"],"falsifier":"A set of numerical stellar models at a boson mass and coupling value within the studied ranges whose maximum mass deviates from the predicted scaling by more than 4 percent would falsify the claimed precision of the relations.","tokens_in":2952,"feed_emoji":"🌌","tokens_out":753,"duration_ms":24887,"temperature":0.7,"pith_summary":"The paper systematically studies the properties of stars made from bosonic dark matter with self-interactions described by a quartic potential. Numerical solutions of the structure equations across a broad range of particle masses and couplings yield scaling relations for the maximum mass, critical radius, and critical central density. These relations are presented with explicit prefactors and hold with relative errors below 4 percent. The authors also supply unified analytical fits for how mass and radius vary with central density on the stable branch, accurate to 0.1 percent. This approach replaces repeated numerical integrations with direct formulas that depend only on the dark matter particle parameters.","feed_headline":"Scaling laws fix max mass of bosonic dark matter stars","feed_subtitle":"Formulas give mass, radius and density from boson mass and coupling with errors below 4 percent","key_machinery":"Scaling relations derived from numerical integration of the stellar structure equations using the equation of state for a complex scalar field with quartic self-interaction potential.","core_discovery":"Numerical solutions of the stellar structure equations with the equation of state from a complex scalar field with quartic potential produce the scaling relations M_max = 0.1 sqrt(λ)/m_φ² solar masses, R(M_max) = 0.9 sqrt(λ)/m_φ² km, and ε_max = 2.1×10^5 m_φ⁴/λ MeV/fm³ where m_φ is in GeV, with fitting relative error less than 4 percent. The mass-central density and radius-central density relations on the stable branch are described by a single functional form with parameters chosen separately for mass and radius, achieving fitting relative error less than 0.1 percent. A simple quadratic polynomial mass-radius relation is also identified.","pith_inferences":["The scaling relations could constrain the allowed range of dark matter particle parameters if any bosonic stars are observed.","The unified fitting function might apply to other interaction potentials or be derived analytically in limiting cases.","Comparison with the mass-radius relations of neutron stars or other exotic compact objects could distinguish bosonic dark matter stars observationally."],"forward_implications":["Maximum mass and size of bosonic dark matter stars can be calculated directly from the boson mass and coupling without solving the differential equations each time.","The stable branch configurations collapse onto universal curves when expressed in terms of the critical values.","The mass-radius relation takes a simple quadratic form that follows from the underlying equation of state.","These properties hold uniformly for boson masses spanning twelve orders of magnitude and couplings from 0.01π to 100π."],"fun_headline_variants":["Scaling relations fix bosonic dark matter star max mass","Precise scaling links boson mass to star radius and density","Analytical fits describe stable bosonic dark matter branches","Unified functions fit bosonic DM star mass and radius curves"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The numerical results produce stable prefactors in the scaling relations that can be captured by the chosen functional forms without additional systematic effects across the parameter ranges examined.","fun_headline_variants_meta":{"raw":{"variants":["Scaling relations fix bosonic dark matter star max mass","Precise scaling links boson mass to star radius and density","Analytical fits describe stable bosonic dark matter branches","Unified functions fit bosonic DM star mass and radius curves"]},"model":"grok-4.3","cost_usd":0.005277,"raw_usage":{"total_tokens":2630,"prompt_tokens":985,"num_sources_used":0,"completion_tokens":62,"cost_in_usd_ticks":52765500,"prompt_tokens_details":{"text_tokens":985,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1583,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":985,"tokens_out":62,"duration_ms":10649,"temperature":1.0,"reasoning_tokens":1583,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T15:51:40.826850+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A set of numerical stellar models at a boson mass and coupling value within the studied ranges whose maximum mass deviates from the predicted scaling by more than 4 percent would falsify the claimed precision of the relations.","supporting_citations":[],"review_version":1}