{"id":"25810200-09f8-4e83-a3ea-bfe9af0d0c63","arxiv_id":"1908.08247","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The authors use SPARC galaxy rotation curves to claim an upper bound m ≤ 10^{-31} eV on the graviton mass in quasi-dilaton massive gravity.","lead":"This paper derives a modified galaxy rotation curve from quasi-dilaton massive gravity and claims an upper bound on the graviton mass of 10^{-31} eV. The result is based on fits to a small selection of SPARC galaxies, with several approximations and no full statistical treatment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The rotation-curve bound relies on an unjustified 'quintic-dominant' truncation; at galactic scales the discarded terms are comparable, so Eq. (26) and m ≤ 10^-31 eV are not established.","rationale":"I read the paper as aiming to derive a theoretically-motivated rotation-curve template for QDMG and to use it to set an upper bound on the graviton mass. For the central claim to hold, the template must be a faithful consequence of QDMG in the regime of galaxies. The weakest point is the transition from the full quintic (39) to the approximate equation (22). The authors simply state that the quintic term is dominant, but do not check that the lower-order terms are small at the scales of interest. My scaling estimate indicates they are not: with typical galaxy parameters and m ≈ 10^-31 eV, the dimensionless source a = A/Λ3^3 is of order unity, and the polynomial coefficients are also order unity, so the solution sits in a regime where all terms contribute. If the approximation fails, the exponent 3/5 and the parameter γ in eqs. (26)-(27) are not the correct prediction, and the bound m ≤ 10^-31 eV could be an artifact. This concern is more fundamental than the reader's weakest assumption about externally-fixed NFW parameters: even if the halo parameters were re-derived jointly with γ, the template itself would be wrong. The paper's own caveat section does not acknowledge this truncation as an approximation, which makes it more likely to be an unexamined step. The numerical check I propose is straightforward and would settle the issue. I do not think the paper should be rejected outright, because the full equations are present in the appendix and the approximation may be salvageable for some parameter choices, but the present derivation does not support the headline bound. Hence the verdict remains CONDITIONAL, with the condition being to solve the full quintic or justify the truncation.","tokens_in":10920,"tokens_out":22107,"duration_ms":183173,"concrete_test":"Take the best-fit values of ω, α3, α4 from ref. [25] (the paper does not quote them) and a representative galaxy from Figure 1, e.g. NGC7814 with its SPARC data. Solve the full quintic (39) numerically for x(R) at each radius, compute the rotation curve v^2(R) = r^2 z(x(R)) from eq. (18), and compare it to eq. (26) with the same halo parameters. If the two curves differ by more than the data error bars (or the effective slope in R differs from the R·f^{3/5} form), then eq. (22) is invalid and the bound must be recomputed. A simpler diagnostic: evaluate the magnitudes of the neglected terms c1 y, c2 y^2, c3 y^3, c4 y^4 at the solution y of the truncated quintic; if any is ≥ 10% of c5 y^5, the truncation is uncontrolled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section B, after deriving the full quintic for x = φ'/r (Appendix eq. 39), the authors replace it by eq. (22), keeping only the x^5 term: -3/Λ3^12 [2(α3+4α4)^2 + (1/ω)(α3-4α4)^2] x^5 = A. No justification is given for neglecting the linear through quartic terms. A scaling estimate shows the issue: for a typical SPARC galaxy (M ~ 10^11 M_sun, r ~ 10 kpc) and m ~ 10^-31 eV, Λ3^3 = m^2 M_Pl ~ 1.2e-61 GeV^3 and A = M/(8π M_Pl r^3) ~ 1e-61 GeV^3, so a = A/Λ3^3 ~ O(1). In dimensionless form the full quintic is c1 y + c2(a) y^2 + c3 y^3 + c4 y^4 + c5 y^5 = a, with coefficients built from O(1) parameters; the root y then has all terms of comparable magnitude, not a single dominant power. This means eq. (22) is not a controlled approximation, and the resulting velocity profile (26), with its characteristic R [ln(1+R)-R/(1+R)]^{3/5} correction and the parameter γ (27), is not the generic prediction of QDMG in this regime. Since the bound m ≤ 10^-31 eV is extracted from fitting exactly this truncated template, the central claim is not established unless the full quintic is solved or the truncation is shown to be valid for the galaxies used. The paper explicitly lists other approximations but does not flag this one.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies quasi-dilaton massive gravity (QDMG), a scalar-extended dRGT theory, and derives constraints on the graviton mass. The first part estimates the decay of gravitational-wave tensor modes into scalar modes and finds the decay width negligible. The second part works in the decoupling limit, derives equations for the Bardeen potentials and the scalar perturbations, truncates the full quintic equation for x to its x^5 term, and fits the resulting modified Navarro-Frenk-White rotation-curve template to SPARC data. The central claim is an upper bound m ≤ 10^-31 eV on the graviton mass in QDMG, with the conclusion that the dark-energy-motivated value m ~ 10^-33 eV remains viable.","tokens_in":11356,"tokens_out":8764,"duration_ms":88032,"significance":"If the bound were established, it would be a useful astrophysical constraint on QDMG and would complement existing LIGO and Solar System bounds. The paper is clearly written and makes a genuine attempt to connect the decoupling-limit equations to galaxy data; it also correctly notes several of its own approximations and the need for a future MCMC treatment. However, the main result depends on an uncontrolled truncation of the equation of motion, and the statistical treatment is too weak to give a robust quantitative bound. The gravitational-wave depletion argument is a reasonable order-of-magnitude null result, but it is not the main claim.","major_comments":[{"comment":"The replacement of the full quintic (39) by the x^5-only equation (22) is not controlled. The full equation contains terms proportional to x, x^2, x^3, x^4 and x^5, and no argument is given that the discarded terms are small at the relevant galactic radii. For a typical SPARC galaxy with M ~ 10^11 M_sun and r ~ 10 kpc, the dimensionless combination A/Lambda_3^3 is of order one at m ~ 10^-31 eV, so the coefficients of the lower powers are not parametrically suppressed relative to the quintic term. The root of the full quintic need not be close to the root of Eq. (22). Since the velocity profile (26), the fitted parameter gamma (27), and the bound (29) all follow from Eq. (22), the central claim is not established unless the full quintic is solved or a controlled dominance argument is supplied for the galaxies actually used.","section":"III B, Eq. (22), Appendix Eq. (39)"},{"comment":"The bound relies on adopting Upsilon_disk, Upsilon_bulge, r_s and rho_s from the GR-based NFW fits of de Almeida et al. [39]. In QDMG the rotation curve is modified, so the inferred halo and stellar parameters should be re-derived jointly with the QDMG parameter gamma. Using GR-fit halo parameters can shift gamma and therefore the inferred graviton mass. The paper acknowledges this and defers a full MCMC, but as it stands the quoted inequality m <= 10^-31 eV is conditional on the GR halo model. In addition, no uncertainty is quoted for gamma or m, so the bound lacks a statistical definition.","section":"III B, after Eq. (28)"},{"comment":"The applicability of the decoupling limit to galactic rotation curves is asserted rather than demonstrated. The text states that the decoupling limit is valid for 'typical scales bigger than 1/m'; for the masses considered here, 1/m ~ 10^24 m (tens of megaparsecs), so galactic scales are not in that regime. If the intended criterion is instead r << 1/m, that should be stated explicitly and the ordering of scales checked. As written, the domain of validity of Eq. (3) does not clearly cover the SPARC data used in Section III B, which directly affects the derivation of the central bound.","section":"II and III B"}],"minor_comments":[{"comment":"The quantity A in Eq. (17) is defined with M(r), but in Eq. (22) it appears as a source term; please clarify whether M(r) is the total enclosed mass or only the dark matter mass entering the NFW profile.","section":"III B, Eq. (17)"},{"comment":"Only two representative galaxies are shown, and no error bars are visible on the data points. The text states that 'the galaxies that we take show consistency' in gamma, but a table or plot reporting gamma for every galaxy used, with uncertainties, is needed for reproducibility.","section":"III B, Figure 1"},{"comment":"The paper oscillates between 'MPl' and 'M_Pl' for the Planck mass; please unify the notation.","section":"Notation"}],"recommendation":"major_revision","confidential_remarks":"The main obstacle is the uncontrolled truncation of the quintic in Eq. (22), which is load-bearing for the claimed bound. The authors should either solve the full quintic or provide a controlled expansion for the actual galaxy sample. A joint fit of gamma with the halo and stellar parameters would also be needed to make the bound quantitative. The gravitational-wave decay section is a minor part of the paper and appears sound as an order-of-magnitude estimate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's novel bit is the QDMG-specific rotation-curve template: a modified NFW velocity profile with a characteristic R [ln(1+R)-R/(1+R)]^{3/5} correction and a parameter γ tied to the graviton mass. Applying the known beyond-Horndeski fitting method to QDMG is a legitimate extension, and the authors are upfront that a full MCMC treatment is missing. The GW decay estimate is a null result that agrees with earlier Horndeski work; that part is fine and not overclaimed.\n\nThe soft spot is load-bearing. In the appendix they give the full quintic for x (eq. 39), then keep only the x^5 term in eq. (22). The truncation is not justified. For a typical SPARC galaxy and m around the claimed bound, Λ3^3 and A are the same order, so the dimensionless source a is order one. In the dimensionless quintic all powers contribute comparably; there is no regime where the x^5 term dominates. I checked the scaling and it holds. That means eq. (26) and the subsequent γ fit are not the generic QDMG prediction in the galaxy regime. The bound m ≤ 10^-31 eV is derived by fitting a template that may not follow from the theory. The authors list several approximations but do not flag this one.\n\nA second, smaller issue: they take halo parameters and stellar mass-to-light ratios from a GR-based fit in ref [39] and do not re-fit them jointly with γ. They acknowledge this, so it is not a hidden flaw, but it makes the bound conditional on the GR halo inference. It is a minor-to-moderate concern relative to the truncation problem.\n\nWhat is good: the paper is honestly written, the derivation up to the full quintic is transparent, and the null GW decay result is consistent with prior literature. The citation pattern is reasonable; using best-fit α and ω coefficients from [25] is appropriate, though it inherits their uncertainties.\n\nWho should read this: anyone working on massive gravity phenomenology. It is a plausible first step, not a definitive constraint. I would send it to a referee, but with the expectation that the truncation issue gets fixed or at least convincingly bounded. As it stands, I would not cite the bound in my own work.","headline":"A sensible application of rotation-curve methods to QDMG, but the claimed bound rests on an unjustified quintic-dominant truncation, so m ≤ 10^-31 eV is not yet established.","tokens_in":11854,"tokens_out":5979,"would_cite":false,"duration_ms":54352,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","98.62.Dm","04.30.-w"],"model":"deepseek-v4-flash","headline":"Galaxy rotation curves cap the graviton mass in quasi-dilaton massive gravity at $10^{-31}$ eV.","keywords":["quasi-dilaton massive gravity","graviton mass bound","galaxy rotation curves","SPARC catalogue","Navarro-Frenk-White profile","Bardeen potentials","decoupling limit","dark energy"],"falsifier":"Re-fit the SPARC galaxies used here (or a larger sample) with a Markov Chain Monte Carlo that varies $\\gamma$, $r_s$, $\\rho_s$, and the stellar mass-to-light ratios simultaneously; if the best-fit $\\gamma$ corresponds to $m > 10^{-31}\\,{\\rm eV}$ while still matching the rotation curves, the paper's bound is false.","tokens_in":2151,"feed_emoji":"🌌","tokens_out":2954,"duration_ms":102063,"temperature":0.7,"pith_summary":"Quasi-dilaton massive gravity (QDMG), a scalar extension of ghost-free dRGT massive gravity, adds a quasi-dilaton field whose extra scalar degrees of freedom alter the Bardeen potentials of galaxies. The paper derives the resulting modification to dark-matter rotation curves and fits them to high-quality SPARC galaxies, obtaining an upper bound on the graviton mass of $m \\leq 10^{-31}\\,{\\rm eV}$. If this bound holds, QDMG remains viable as a dark-energy candidate at $m \\approx 10^{-33}\\,{\\rm eV}$, while ruling out heavier QDMG gravitons up to the previous $10^{-22}\\,{\\rm eV}$ limit. It also shows that tensor-to-scalar gravitational-wave decay is far too weak to deplete observable signals, so gravitational-wave observations cannot easily falsify the theory. The main claim is that QDMG's astrophysical predictions are consistent with existing rotation-curve data only if the graviton is extremely light.","feed_headline":"Quasi-dilaton massive gravity needs gravitons below 10^-31 eV","feed_subtitle":"Rotation-curve data keep the theory alive for dark energy but put space-based detection out of reach.","key_machinery":"The central object is the parameter $\\gamma$ defined in Eq. (27), which packages the graviton mass and the QDMG coupling coefficients into the single quantity controlling how much the dark-matter rotation curve deviates from the pure Navarro-Frenk-White/general-relativity prediction. It is derived from the decoupling-limit equations for the two Bardeen potentials and two galileon-type scalars; after assuming spherical symmetry the system reduces to a quintic for $x = \\phi'/r$, and keeping the dominant $x^5$ term produces the modified velocity profile. A second piece of machinery is the $1\\to 3$ decay vertex $h(\\partial^2\\pi)^3/\\Lambda_3^6$ used to estimate tensor-to-scalar gravitational-wave depletion, with the phase-space integral evaluated following the standard three-body decay treatment.","core_discovery":"The paper's central claim is that quasi-dilaton massive gravity, taken in its decoupling limit up to cubic galileon order and with the quintic term dominant, predicts a dark-matter velocity profile of the form $v^2_{dm} = 4\\pi G r_s^2 \\rho_s \\left[ \\frac{1}{R}\\left(\\ln(1+R) - \\left(1+\\frac{1}{R}\\right)^{-1}\\right) - \\gamma R^{1/5} \\left(\\ln(1+R) - \\left(1+\\frac{1}{R}\\right)^{-1}\\right)^{3/5} \\right]$, where $\\gamma$ encodes the graviton mass $m$ through $\\gamma = (4\\pi G\\rho_s)^{-2/5} m^{4/5} (\\alpha_3+4\\alpha_4) \\left[6(\\alpha_3+4\\alpha_4)^2 + \\frac{3}{\\omega}(\\alpha_3-4\\alpha_4)^2\\right]^{-3/5}$. Fitting this profile to high-quality SPARC galaxies, using Navarro-Frenk-White halo parameters and stellar mass-to-light ratios taken from GR-based fits, yields consistency only for $m \\leq 10^{-31}\\,{\\rm eV}$. This bound is stronger than the previous $10^{-22}\\,{\\rm eV}$ level from gravitational waves and Solar System tests, but it still permits the dark-energy-motivated value $m \\approx 10^{-33}\\,{\\rm eV}$. The paper also computes that the $1\\to 3$ decay of tensor to scalar modes has a width $\\Gamma \\approx 10^{-225}\\,{\\rm eV}$, far too small to observably deplete gravitational-wave signals.","pith_inferences":["A joint fit that varies $\\gamma$, the halo parameters, and the stellar mass-to-light ratios simultaneously could shift the inferred bound; if QDMG changes the effective halo, the bound might move by orders of magnitude.","Extending the calculation beyond the cubic galileon or including nonzero background values of the scalar fields could add terms that compete with the quintic, possibly altering the predicted velocity profile and the mass bound.","The same rotation-curve test could be applied to other scalar-extended massive-gravity theories with stable cosmological solutions, not just QDMG.","The bound is derived from high-mass, high-luminosity galaxies with the best-quality rotation curves; testing lower-mass or dwarf galaxies, where the NFW profile fits more poorly, could expose deviations or strengthen the constraint."],"forward_implications":["If the bound is right, QDMG gravitons heavier than $10^{-31}\\,{\\rm eV}$ are ruled out, shrinking the allowed mass window to $m \\lesssim 10^{-31}\\,{\\rm eV}$.","The dark-energy target $m \\approx 10^{-33}\\,{\\rm eV}$ survives, so QDMG remains a viable explanation of cosmic acceleration.","Space-based gravitational-wave observatories with sensitivity up to $10^{-25}\\,{\\rm eV}$ will not reach this bound, so direct detection of the QDMG graviton mass is out of reach.","Depletion of gravitational-wave signals by tensor-to-scalar decay is negligible, so this decay channel cannot falsify QDMG observationally.","Since the Bardeen potentials satisfy $\\Psi = \\Phi$ as in general relativity, QDMG is indistinguishable from GR in gravitational lensing tests, making rotation curves the discriminating probe."],"supporting_citations":[{"why":"Supplies the NFW halo parameters and stellar mass-to-light ratios used in the QDMG rotation-curve fits.","marker":"[39]"},{"why":"Provides the SPARC galaxy rotation-curve data and quality flags used for the comparison.","marker":"[36]"},{"why":"Supplies the best-fit QDMG coefficients $\\omega$ and $\\alpha_n$ used to evaluate the parameter $\\gamma$ and hence the mass bound.","marker":"[25]"},{"why":"Defines quasi-dilaton massive gravity as the dRGT extension with the quasi-dilaton field.","marker":"[12]"},{"why":"Establishes dRGT ghost-free massive gravity, the base theory whose scalar extension QDMG is.","marker":"[11]"},{"why":"Defines the Navarro-Frenk-White halo profile assumed for the dark matter distribution.","marker":"[38]"},{"why":"Provides the three-body phase-space treatment used for the tensor-to-scalar decay width estimate.","marker":"[35]"},{"why":"Supports the conclusion that higher-derivative corrections leave gravitational-wave signals undepleted, matching the decay-width result.","marker":"[32]"},{"why":"Records previous graviton-mass bounds at the $10^{-22}\\,{\\rm eV}$ level that the new rotation-curve bound improves for QDMG.","marker":"[31]"},{"why":"Lists existing Solar System and gravitational-wave graviton-mass sensitivities used to compare the new bound.","marker":"[41]"}],"fun_headline_variants":["Rotation curves constrain quasi-dilaton graviton mass to 10^-31 eV","Graviton mass bound sharpened to 10^-31 eV by galaxy data","Quasi-dilaton gravity passes galaxy rotation tests with tiny graviton","SPARC galaxies set tightest graviton mass limit in QDMG","Quasi-dilaton gravity survives, but graviton too light for LISA"],"cache_read_input_tokens":13824,"weakest_assumption_plain":"The argument assumes that the NFW halo parameters ($r_s$, $\\rho_s$) and stellar mass-to-light ratios obtained from GR-based fits in ref. [39] can be used unchanged when fitting the QDMG parameter $\\gamma$; if QDMG alters rotation curves, those parameters should be re-derived jointly, and the inferred graviton mass would change.","fun_headline_variants_meta":{"raw":{"variants":["Rotation curves constrain quasi-dilaton graviton mass to 10^-31 eV","Graviton mass bound sharpened to 10^-31 eV by galaxy data","Quasi-dilaton gravity passes galaxy rotation tests with tiny graviton","SPARC galaxies set tightest graviton mass limit in QDMG","Quasi-dilaton gravity survives, but graviton too light for LISA"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000761,"raw_usage":{"total_tokens":3441,"prompt_tokens":1073,"completion_tokens":2368,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":689,"completion_tokens_details":{"reasoning_tokens":2266}},"tokens_in":689,"tokens_out":2368,"duration_ms":16224,"temperature":1.0,"reasoning_tokens":2266,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:45:15.546695+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the SPARC galaxies used here (or a larger sample) with a Markov Chain Monte Carlo that varies $\\gamma$, $r_s$, $\\rho_s$, and the stellar mass-to-light ratios simultaneously; if the best-fit $\\gamma$ corresponds to $m > 10^{-31}\\,{\\rm eV}$ while still matching the rotation curves, the paper's bound is false.","supporting_citations":[{"cited_title":"Navarro, Carlos S","cited_arxiv_id":null,"evidence_quote":"Supplies the NFW halo parameters and stellar mass-to-light ratios used in the QDMG rotation-curve fits."},{"cited_title":"Review of particle physics","cited_arxiv_id":null,"evidence_quote":"Provides the SPARC galaxy rotation-curve data and quality flags used for the comparison."},{"cited_title":"Tate Deskins, John T","cited_arxiv_id":null,"evidence_quote":"Supplies the best-fit QDMG coefficients $\\omega$ and $\\alpha_n$ used to evaluate the parameter $\\gamma$ and hence the mass bound."},{"cited_title":"Implications of the neutron star merger gw170817 for cosmological scalar- tensor theories","cited_arxiv_id":null,"evidence_quote":"Establishes dRGT ghost-free massive gravity, the base theory whose scalar extension QDMG is."},{"cited_title":"Generalized framework for testing gravity with gravitational-wave propagation","cited_arxiv_id":null,"evidence_quote":"Defines the Navarro-Frenk-White halo profile assumed for the dark matter distribution."},{"cited_title":"Gravitational Rainbows: LIGO and Dark Energy at its Cutoff","cited_arxiv_id":null,"evidence_quote":"Provides the three-body phase-space treatment used for the tensor-to-scalar decay width estimate."},{"cited_title":"Tate Deskins, Andrew J","cited_arxiv_id":null,"evidence_quote":"Supports the conclusion that higher-derivative corrections leave gravitational-wave signals undepleted, matching the decay-width result."},{"cited_title":"McGaugh, Arianna Di Cintio, Chris B","cited_arxiv_id":null,"evidence_quote":"Lists existing Solar System and gravitational-wave graviton-mass sensitivities used to compare the new bound."}],"review_version":1}