{"id":"36c9175c-14bc-49bc-874b-7bbfcea7b134","arxiv_id":"2501.16492","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"Bulk flow speeds and directions are fitted to supernova data, and the neutrino-coupled f(R) model gives the largest fitted velocities, but the claim is a comparison of fit parameters rather than a prediction.","lead":"This paper fits bulk flow velocities and directions to Pantheon supernova data in three gravity models: f(R), perturbed f(R), and perturbed f(R) coupled to neutrinos. It reports that adding neutrinos raises the fitted bulk flow speeds and aligns the flow with cosmic superclusters and the dark energy dipole.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed neutrino-induced bulk-flow enhancement is within 1σ of the paper's own quoted errors, and V_bulk is a fitted parameter, so the central result is not statistically supported.","rationale":"The reader's weakest assumption, the ad hoc neutrino-coupling parameterization, is a real and serious concern, and I agree that it could invalidate the model if wrong. However, I would place the most load-bearing problem one step earlier: the paper's own quoted error bars show that the reported velocity differences between the f(R), perturbed f(R), and neutrino-coupled models are all within about 1σ. For the key bin 0.8 < z < 1.4, the difference is 306 km/s with σ_Δ ≈ 401 km/s, so the headline 'substantial increase' is not established by the data presented. This is independent of whether the coupling form is physically correct. The reader's rationale did mention overlapping error bars and the fact that V_bulk is a free parameter in Eq. (56), but the formal weakest_assumption focused on the coupling form rather than on this more elementary statistical failure. Therefore I partially agree with the reader's framing while strengthening the argument: even a physically perfect model cannot support the central claim from these fits without a significance test. The absence of code, data, and significance tests strengthens the case for rejection, though the statistical overlap alone is sufficient. My concrete test—a covariance-aware bootstrap refit—would settle the issue computationally; if the 68% interval of the velocity difference contains zero, the claim is unsupported.","tokens_in":23614,"tokens_out":6114,"duration_ms":62495,"concrete_test":"Bootstrap the Pantheon distance moduli using their full covariance, with the same redshift cuts and the same fitting procedure as Eq. (56), and refit V_bulk and direction independently for the f(R) and neutrino-coupled models. Then compute the distribution of ΔV = V_neu − V_fR in each redshift bin; in particular, check whether the 68% interval of ΔV contains zero for 0.8 < z < 1.4, as the quoted 306 ± 401 km/s suggests. If zero is inside the interval in every bin, the claimed neutrino-induced enhancement is not statistically significant and the central claim collapses.","verdict_should_be":"REJECT","load_bearing_attack":"The load-bearing problem is that the paper's headline enhancement is not significant even in its own numbers. For 0.8 < z < 1.4, the neutrino-coupled model gives V_bulk = 3086 ± 286 km/s versus 2780 ± 282 km/s for f(R); the difference is 306 km/s with σ_Δ ≈ 401 km/s, i.e. about 0.8σ. Other bins are similar: 0.4 < z < 0.6 gives Δ ≈ 40 ± 327 km/s; 0.1 < z < 0.2 gives Δ ≈ 123 ± 202 km/s; local bins are below 1.4σ. No significance test or covariance propagation is presented, and the χ² differences between models are small (e.g. Δχ² ≈ 1.9 for the 0.8 < z < 1.4 bin). Moreover, V_bulk is a free parameter fitted per redshift bin in Eq. (56); the dipole term is d_L^dipole = (1+z)/H (n·v_Bulk), so the model dependence enters only through H(z) and d0_L(z). Thus the reported 'substantial increase' is not a dynamical prediction of neutrino coupling but a fitted amplitude whose difference is consistent with noise. The ad hoc coupling form Γ = u^μ∇_μ f_R in Eqs. (19) and (23) is a further concern, but it is not the most decisive issue: the central claim fails on statistical grounds before the model's physical validity is even assessed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes bulk flow in f(R) gravity, perturbed f(R) gravity, and perturbed f(R) gravity with a neutrino coupling, using Pantheon Type Ia supernova data and the Bonvin dipole formalism. For each redshift bin, the bulk flow velocity and direction are obtained by minimizing a chi-square over the Pantheon magnitude residuals. The authors report that the neutrino-coupled model yields the largest bulk flow velocities, exceeding 3000 km/s at 0.8<z<1.4, and that its direction aligns with the Sloan Great Wall, the King Ghidorah Supercluster, and the dark energy dipole. The central claim is that neutrino interactions with modified gravity substantially enhance bulk flow and align it with large-scale structure.","tokens_in":23971,"tokens_out":4671,"duration_ms":46154,"significance":"If the central claim were a genuine model prediction, it would be of interest because it would connect neutrino-modified gravity to an observable large-scale velocity field. However, the bulk flow velocity and direction are free parameters fitted separately in each redshift bin, so the reported 'increase' is a comparison of best-fit amplitudes, not a prediction of the neutrino coupling. Moreover, the differences between models are within the quoted 1-sigma uncertainties in every bin, and no significance test is presented. The direction results also rest on an apparent error in the galactic coordinate transformation. The paper does provide a useful compilation of redshift-binned bulk-flow results and a clear description of the fitting procedure, but those strengths do not overcome the statistical and methodological problems with the headline claim.","major_comments":[{"comment":"The bulk flow velocity and direction are free parameters minimized independently in each redshift bin via Eq. (56). The model dependence enters only through the Hubble parameter and the background luminosity distance in the dipole term of Eq. (51). Therefore, the statement that incorporating neutrinos 'results in a substantial increase in bulk flow velocities' is a statement about fitted values, not a dynamical prediction. To support the claim, the paper would need to show that the perturbation equations of Sec. 4 predict a larger bulk flow amplitude in the neutrino-coupled model without using the Pantheon dipole as input.","section":"Sec. 6, Eq. (56), Tables 5-7"},{"comment":"The reported enhancements are not statistically significant even within the paper's own errors. For 0.8<z<1.4, the neutrino model gives 3086±286 km/s versus 2780±282 km/s for f(R); the difference is 306 km/s with combined error about 401 km/s, i.e., roughly 0.8 sigma. For 0.4<z<0.6 the difference is 40±326 km/s, and for 0.1<z<0.2 it is 123±202 km/s. No bin exceeds 1 sigma, and the chi-square differences between models are small (e.g., Delta chi^2 about 1.9 for the highest bin). The conclusion that neutrinos substantially increase bulk flow is therefore not supported by the analysis as presented.","section":"Tables 5-7 and Sec. 7.1"},{"comment":"The Cartesian unit vector for galactic coordinates appears to be incorrect. For a source at galactic longitude l and latitude b, the standard expression is (cos b cos l, cos b sin l, sin b), but Eq. (53) uses (cos l sin b, sin l sin b, cos b), and Eq. (55) repeats the same convention. As written, sources with b=0 are all assigned to the pole, and the fitted bulk flow directions, including the claimed alignments with the Sloan Great Wall, King Ghidorah Supercluster, and dark energy dipole, are therefore suspect.","section":"Sec. 6, Eqs. (53)-(55)"},{"comment":"The neutrino coupling is introduced as Q_nu = -Gamma rho_nu with Gamma = u^mu grad_mu f_R, but no derivation from the interaction Lagrangian in Eq. (17) is given, and the sign and functional form are simply assumed. The manuscript then uses Gamma = 0.6 ± 0.25, quoted as a best fit to Pantheon+ data in Sec. 5, without explaining how this value is obtained or whether it is re-fit in Eq. (56). If Gamma is a global best fit to the same Pantheon sample used for the bulk-flow dipole, the subsequent enhancement is at least partly circular.","section":"Secs. 3 and 4, Eqs. (19), (23)-(25)"},{"comment":"The constraint on the sum of neutrino masses, Sigma m_nu < 0.142 eV, and the quoted Gamma = 0.6 ± 0.25 are not supported by any presented likelihood, priors, dataset, or fitting procedure. Equations (42)-(45) define a dimensionless density variable but do not by themselves yield a mass constraint. Since Gamma is the only new parameter in the neutrino-coupled model, this missing support is load-bearing for the model comparison.","section":"Sec. 5, Eqs. (42)-(45)"}],"minor_comments":[{"comment":"Several equations contain apparent typographical or OCR errors, such as the term '-g_mu_nu 2 f_R(R)' in Eq. (2) and similar expressions in Eq. (9), which likely should involve the d'Alembertian operator. The notation for derivatives is also inconsistent (fR, f'_R, f''_R) and should be defined uniformly.","section":"Secs. 2-3, Eqs. (2)-(9)"},{"comment":"The text says 'Inserting this in Eq. (21)' when referring to the dipole formula; the intended equation number appears to be Eq. (51) or (52).","section":"Sec. 6, text following Eq. (50)"},{"comment":"The direction uncertainties in Tables 5-7 are large, e.g., b = 69 ± 18 degrees for the 0.4<z<0.6 bin, so statements such as 'near-perfect congruence' with the dark energy dipole or superclusters overstate the precision of the measurement.","section":"Sec. 7 and Appendix"},{"comment":"The CMB power-spectrum discussion is not connected to the bulk-flow fitting procedure, and no Boltzmann solver or parameter choices are described. As presented, the figure and appendix do not add quantitative support to the main claims.","section":"Appendix, Fig. 11"},{"comment":"The reference list contains many duplicated and incomplete entries, including multiple Kashlinsky entries, inconsistent spellings such as 'Watkin' versus 'Watkins', and garbled citations. The paper would need a careful editorial pass before publication.","section":"References"}],"recommendation":"reject","confidential_remarks":"The central result is not supported by the paper's own statistical errors, and the main observable is a fitted parameter rather than a model prediction. The coordinate-transformation error further undermines the direction claims. These are load-bearing issues that cannot be fixed by minor revisions; the manuscript would need to reframe its claims as a null or marginal result and add a proper significance analysis, which is beyond a standard revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take: the headline claim doesn't survive contact with the paper's own numbers. The bulk flow velocities are derived by minimizing χ² in Eq. (56) over V_bulk and dipole direction for each redshift bin, so the comparison across models is a comparison of best-fit values, not a prediction of the model. And the differences are not significant: in the 0.8<z<1.4 bin, the neutrino-coupled model gives 3086±286 km/s versus 2780±282 km/s for f(R); the difference is about 0.8σ. Other bins are similar. The alignment with the dark energy dipole and superclusters is asserted without any significance test.\n\nWhat's new: applying the Bonvin dipole method to a neutrino-coupled perturbed f(R) model, with redshift tomography on the Pantheon dataset, is a legitimate extension of the authors' earlier work. The research question—whether modified gravity plus neutrino coupling leaves a detectable imprint on large-scale velocities—is reasonable, and they include error bars and compare with previous bulk flow studies.\n\nSoft spots, in order of weight. First, the statistical issue above is load-bearing; the abstract's \"substantial increase\" is not supported by the quoted uncertainties. Second, the coupling Γ = u^μ∇_μ f_R is introduced as a parameterization with no derivation; the reader is asked to take it on faith, and it carries the entire neutrino effect. Third, the perturbation equations as typeset contain several malformed terms, so I could not verify the derivation without going back to the source files. Fourth, Γ was fitted to Pantheon+ in the authors' earlier paper and then used here on the Pantheon sample; that is a carry-over rather than a circular fit, but it deserves more scrutiny than it gets. No code or data is released, which makes the results hard to check.\n\nBottom line: the paper is for someone interested in bulk flow in modified gravity, but the central claim fails on the paper's own statistics. I would not send this to review in its current form; the appropriate outcome is a desk reject or a major rewrite that reframes the analysis as a fitting exercise and tests whether the model-to-model differences are significant.","headline":"The claimed neutrino-driven bulk flow enhancement is a fitted-parameter comparison that is not significant in the paper's own error bars.","tokens_in":24511,"tokens_out":3818,"would_cite":false,"duration_ms":34872,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Neutrino coupling in perturbed f(R) gravity amplifies cosmic bulk flow to over 3000 km/s at high redshift.","keywords":["bulk flow","f(R) gravity","neutrino coupling","Hu-Sawicki model","Pantheon supernovae","dark energy dipole","redshift tomography","large-scale structure"],"falsifier":"A direct test is to refit the same Pantheon distance moduli with the full covariance matrix and check whether the $0.8<z<1.4$ dipole still has amplitude near $3086\\,\\mathrm{km/s}$ pointing at $(l,b)=(330^\\circ,-16^\\circ)$; a null or misdirected dipole would falsify the central claim.","tokens_in":23359,"feed_emoji":"🌌","tokens_out":15812,"duration_ms":127360,"temperature":0.7,"pith_summary":"This paper argues that adding a neutrino coupling to perturbed $f(R)$ gravity changes the predicted large-scale velocity field of the universe. Fitting Type Ia supernova distances from the Pantheon catalog in redshift bins, the authors report that the neutrino-coupled model raises the bulk flow velocity in every bin, with values exceeding $3000\\,\\mathrm{km/s}$ in the $0.8<z<1.4$ range. In the same model the flow direction at $z>0.4$ closely follows the dark energy dipole, while at $0.1<z<0.2$ it points toward the Sloan Great Wall and at $0.4<z<0.6$ toward the King Ghidorah Supercluster. If correct, this gives a testable signature of neutrino interactions with modified gravity that standard cosmology does not produce, and it ties neutrino physics to the largest observed cosmic velocities. The paper also reports a neutrino mass sum bound $\\sum m_\\nu < 0.142\\,\\mathrm{eV}$ and a coupling strength $\\Gamma = 0.6 \\pm 0.25$ consistent with earlier work.","feed_headline":"Neutrinos in modified gravity push cosmic flows past 3000 km/s","feed_subtitle":"Redshift-binned supernova data link the high-redshift flow to the dark energy dipole","key_machinery":"The engine of the calculation is a modified neutrino continuity equation in an $f(R)$ background, $\\rho'_\\nu + 3H(\\rho_\\nu + P_\\nu) = -\\Gamma\\rho_\\nu$, where $\\Gamma = u^\\mu\\nabla_\\mu f_R$ measures how the additional scalar degree of freedom of $f(R)$ gravity interacts with the neutrino fluid. The paper evaluates this coupling in the Hu-Sawicki form of $f(R)$, rewrites the linear perturbation equations as a first-order autonomous system in variables $\\xi_1,\\dots,\\xi_8$, and feeds the resulting distance-redshift relation into a dipole formula for the luminosity distance, $d_L^{\\rm(dipole)}(z) = \\frac{1+z}{H}\\,(\\mathbf{n}\\cdot\\mathbf{v}_{\\rm Bulk})$. A $\\chi^2$ fit to Pantheon supernova distance moduli, binned by redshift, then returns the amplitude and direction of the bulk flow for each model.","core_discovery":"The paper's central claim is that the bulk flow, the coherent motion of matter averaged over large volumes, is a discriminating probe of modified gravity, and that neutrinos make it visibly stronger. In the Hu-Sawicki $f(R)$ model, a modified-gravity model whose $f(R)$ is a rational function of the Ricci scalar, adding scalar perturbations and then a neutrino coupling of the form $Q_\\nu = -\\Gamma\\rho_\\nu$ with $\\Gamma = u^\\mu \\nabla_\\mu f_R$ produces a systematic increase in the fitted bulk flow velocity in every redshift bin studied. The highest bin, $0.8<z<1.4$, gives $V_{\\rm bulk} = 3086 \\pm 286\\,\\mathrm{km/s}$ directed at $(l,b) = (330^\\circ \\pm 15^\\circ, -16^\\circ \\pm 17^\\circ)$, which the authors identify with the dark energy dipole; at $0.1<z<0.2$ the flow points toward the Sloan Great Wall, and at $0.4<z<0.6$ toward the King Ghidorah Supercluster, a massive supercluster at $z\\sim0.5$. The authors interpret these alignments and the velocity boost as evidence that neutrinos interacting with the modified gravity sector shape cosmic flows and influence cosmic acceleration.","pith_inferences":["If the proposed coupling is real, the same dipole should survive a full-covariance re-analysis of the Pantheon+ sample; that re-analysis is a direct and inexpensive check.","The claimed high-redshift alignment with the dark energy dipole, if confirmed by independent data, would connect bulk-flow measurements to the long-standing dark-flow and CMB-frame anomalies and give modified gravity a single observable that addresses both.","The linear coupling $Q_\\nu=-\\Gamma\\rho_\\nu$ is a parameterization rather than a derived Lagrangian; testing other functional forms would show whether the $>3000\\,\\mathrm{km/s}$ prediction is generic to neutrino-modified gravity or specific to this choice."],"forward_implications":["In the neutrino-coupled perturbed $f(R)$ model, high-redshift bulk flow exceeds $3000\\,\\mathrm{km/s}$ in $0.8<z<1.4$, a signature that future peculiar-velocity surveys could test.","At $z>0.4$ the bulk flow direction tracks the dark energy dipole, meaning the same sector that drives cosmic acceleration is claimed to steer large-scale velocities.","At lower redshifts the flow points toward the Sloan Great Wall ($0.1<z<0.2$) and the King Ghidorah Supercluster ($0.4<z<0.6$), so the model connects the velocity field to specific observed superclusters.","The neutrino coupling raises bulk flow even in the local universe (from $147\\,\\mathrm{km/s}$ to $173\\,\\mathrm{km/s}$ in $0.001<z<0.016$), leaving a local kinematic imprint of the coupling.","The fit yields $\\sum m_\\nu < 0.142\\,\\mathrm{eV}$ at 95% confidence and $\\Gamma = 0.6\\pm0.25$, giving concrete parameters for future modified-gravity analyses."],"supporting_citations":[{"why":"Supplies the luminosity-distance dipole formula that converts supernova distance moduli into bulk flow velocity and direction.","marker":"Bonvin,Durrer & Gasparini (2006)"},{"why":"Supplies the Pantheon catalog of 1701 Type Ia supernovae used for the redshift-tomography fits.","marker":"Scolnic D. M., 2018"},{"why":"Establishes the neutrino coupling parameterization and the best-fit $\\Gamma = 0.6 \\pm 0.25$ adopted in the perturbed model.","marker":"Yarahmadi et al. (2025)"},{"why":"Defines the King Ghidorah Supercluster's mass, distance, and direction used for the high-redshift alignment comparison.","marker":"Shimakawa et al. (2023)"},{"why":"Documents large observed bulk flows that motivate testing whether modified gravity plus neutrinos can account for them.","marker":"Kashlinsky et al. (2008,2009,2010,2011,2012)"},{"why":"Provides the $f(R)$ perturbation equations and theoretical background on which the perturbed models are built.","marker":"De Felice & Tsujikawa (2010)"},{"why":"Earlier constraints on neutrino mass and bulk flow in $f(R)$ gravity that the paper extends and compares against.","marker":"Yarahmadi & Salehi (2023a,b)"}],"fun_headline_variants":["Neutrinos amplify cosmic bulk flow in f(R) gravity past 3000 km/s","Bulk flow reaches 3086 km/s in neutrino-coupled f(R) gravity","Modified gravity with neutrinos sets cosmic flow toward dark energy dipole","Neutrino-modified gravity boosts bulk flow and aligns with superclusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result rests on the assumption that neutrinos interact with the modified-gravity scalar exactly through the proposed coupling $Q_\\nu = -\\Gamma\\rho_\\nu$ with $\\Gamma\\approx0.6$; if that interaction has a different form or strength, the reported bulk-flow enhancement collapses.","fun_headline_variants_meta":{"raw":{"variants":["Neutrinos amplify cosmic bulk flow in f(R) gravity past 3000 km/s","Bulk flow reaches 3086 km/s in neutrino-coupled f(R) gravity","Modified gravity with neutrinos sets cosmic flow toward dark energy dipole","Neutrino-modified gravity boosts bulk flow and aligns with superclusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000421,"raw_usage":{"total_tokens":2262,"prompt_tokens":1144,"completion_tokens":1118,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":760,"completion_tokens_details":{"reasoning_tokens":1035}},"tokens_in":760,"tokens_out":1118,"duration_ms":10261,"temperature":1.0,"reasoning_tokens":1035,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T12:53:13.172644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test is to refit the same Pantheon distance moduli with the full covariance matrix and check whether the $0.8<z<1.4$ dipole still has amplitude near $3086\\,\\mathrm{km/s}$ pointing at $(l,b)=(330^\\circ,-16^\\circ)$; a null or misdirected dipole would falsify the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Pantheon catalog of 1701 Type Ia supernovae used for the redshift-tomography fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the King Ghidorah Supercluster's mass, distance, and direction used for the high-redshift alignment comparison."}],"review_version":1}