{"id":"f8ba325c-675d-4e82-ba51-bd3f82aabf1a","arxiv_id":"2607.29314","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A two-component virial framework with a fitted per-system interaction term reproduces the radial acceleration relation from galaxies to clusters and weak lensing, but the key coupling is fit, not predicted.","lead":"This paper adds one extra acceleration term to ordinary gravity and uses it to match the well-known radial acceleration relation across galaxies, clusters, and weak-lensing data. The catch: the extra term's strength is fitted to each dataset rather than predicted from dark-matter structure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"VIM's local interaction term (Eqs. A33–A34) is not an acceleration: units are m^-1 s^-2, so the weak-lensing VIM fits cannot support the central claim.","rationale":"I read the paper as claiming that the 2VT/CIA/VIM framework provides a physically motivated ΛCDM explanation of the RAR across galaxies, clusters, and weak-lensing data. The reader's weakest assumption about the CIA's constant-DM approximation is legitimate, but the more decisive problem is in the VIM, which the reader did not identify. The VIM is the component explicitly designed for the low-acceleration weak-lensing regime and is included in the paper's central claim. Its defining local acceleration, Eq. A34, is dimensionally inconsistent: it has units m^-1 s^-2 rather than m s^-2. This is an internal correctness defect, not an interpretive or consensus disagreement. Because the VIM fits in Fig. 4 and Table 1 are based on this formula, those results cannot currently be interpreted as evidence for the framework. The CIA comparisons with public data are transparent, and the paper honestly notes that G is fitted and degenerate, so those parts remain conditionally useful. But the central claim as stated spans both models, and the VIM would need to be re-derived and re-fit, not merely parameter-tweaked. I therefore move the verdict from CONDITIONAL to REJECT for the current manuscript.","tokens_in":20644,"tokens_out":12837,"duration_ms":139448,"concrete_test":"Evaluate Eq. A34 in SI units for the NFW best-fit parameters of Table 1 (ρ0=5.00e2 M_sun kpc^-3, rs=0.656 kpc) at r=50 kpc with M_B=1e10.96 M_sun, and compare to the true Newtonian DM acceleration G M_D(<r)/r^2. If the output is not m/s^2, or if it differs from G M_D/r^2 by orders of magnitude with a different radial scaling, the VIM cannot be a physical explanation of the weak-lensing RAR.","verdict_should_be":"REJECT","load_bearing_attack":"Eq. A33 defines g_BD(r) = [1/(4πr^2 M_B(<r))] |dW_BD(<r)/dr|. W_BD is an energy, so dW_BD/dr has dimensions M L T^-2 (force units). Dividing by 4πr^2 gives M L^-1 T^-2 (pressure units), and dividing by M_B leaves L^-1 T^-2. Eq. A34, g_BD = 4παG ρ_D M_D/(M_B r), has exactly these units: in SI, m^-1 s^-2, not m s^-2. The VIM's 'local interaction acceleration' is therefore not an acceleration. The VIM curves in Fig. 4 and Table 1 are fits to a quantity with the wrong physical dimensions. Because the VIM is introduced specifically to handle the weak-lensing regime and is explicitly included in the abstract/conclusions as part of the 'physically motivated framework,' the VIM component of the central claim fails. This is an internal inconsistency, not a disagreement with ΛCDM. The CIA may remain conditionally plausible, but the strongest claim spans both models.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The authors compile radial-acceleration-relation (RAR) measurements from early-type galaxies, brightest cluster galaxies, galaxy clusters, and isolated-galaxy weak lensing, and interpret them with two models: the Constant-Interaction Approximation (CIA), g_CIA = g_B + g_int with g_int = G a_M, and the Virial-Motivated Interaction Model (VIM), g_VIM = g_B + g_BD(r). The two-component virial theorem (2VT) is invoked to set the scale of g_int and to motivate the local interaction term. The paper reports that the CIA reproduces the RAR for galaxies, BCGs, and clusters with a scale-dependent G, and that the VIM reproduces the low-acceleration weak-lensing RAR for an NFW-like profile. It concludes that the RAR is an emergent, environment-dependent outcome within a baryon-DM coupling framework. The central interpretive claim is that the RAR's small scatter, characteristic acceleration scale, and baryonic predictive power follow from the 2VT without modified gravity or exotic microphysics.","tokens_in":20992,"tokens_out":9424,"duration_ms":103342,"significance":"The dataset compilation and transparent chi-square fits (e.g., NFW VIM: chi^2/dof = 1.68, Table 1) are useful, and the paper is honest about the CIA breakdown at very low accelerations. If the 2VT link were made quantitative, a scale-dependent baryon-DM coupling parameter would be an interesting target for future lensing and kinematic surveys. As it stands, however, the VIM contains a dimensional inconsistency and the CIA's characteristic scale is inserted as a fitted constant, so the central interpretive claims are not established. The paper is better viewed as an empirical additive-acceleration description than as a physical derivation from the 2VT.","major_comments":[{"comment":"Dimensional check: W_BD has dimensions of energy, so dW_BD/dr has dimensions of force (M L T^-2). Dividing by 4 pi r^2 and by M_B(<r) yields L^-1 T^-2, i.e., m^-1 s^-2, not m s^-2. The explicit formula (A34) has the same units. Thus g_BD is not an acceleration and cannot be added to g_B in Eq. (A35). The VIM fits to the Mistele et al. RAR are therefore dimensionally inconsistent. Correcting the definition (e.g., removing the 1/(4 pi r^2) factor) would change the fits and requires a full refit.","section":"App. A.5, Eqs. (A33)-(A34); Table 1; Fig. 4"},{"comment":"The CIA interaction term is specified as g_int = G a_M with G fitted. The 2VT identification g_int = R g_D(r_B) is not used quantitatively: neither R nor g_D(r_B) is computed from the baryonic or DM mass distributions. The paper itself notes the degeneracy in separating G. Hence g_B + g_int is a generic additive fit, not a test of the 2VT. Moreover, the same constant-density approximation gives g_D(r) proportional to r (Eq. A19), not a constant; the constant-interaction step therefore needs independent support.","section":"Sec. 3.4 and App. A.3, Eqs. (9), (A20), (A21)"},{"comment":"The horizontal tail and the transition scale are built in: g_CIA = g_B + g_int asymptotes to g_int when g_B goes to 0. Calling a_M 'emergent' from Eq. (9) is circular, since g_int is normalized to a_M by choosing G. The small-scatter claim is likewise not demonstrated: one G value per population is assigned post hoc, and no predicted scatter from variations in baryonic structure or halo properties is calculated.","section":"Sec. 5, first bullet; App. A.4; Eq. (A20)"},{"comment":"The CIA approximation assumes g_int is constant within individual systems over the sampled radial range. For intracluster data, however, the fits assign a different G to each radial aperture (R = 100, 200, 400, 600 kpc), effectively averaging over systems at fixed radius. The resulting G(R) trend is therefore a restatement of the data, not a test of constancy within each system. To support the cluster interpretation, the paper should check the constancy of g_D(r) within each cluster's baryonic region using the fitted halo parameters.","section":"Sec. 4.1, Fig. 2 right panel; Fig. 3; App. A.3"}],"minor_comments":[{"comment":"G denotes both the gravitational constant (Eqs. A5, A12) and the CIA parameter (Eq. 9); R denotes both the 2VT geometrical parameter and a cluster-centric radius. These should be distinguished throughout.","section":"Notation"},{"comment":"The axis labels and fitted equations contain rendered placeholder symbols for exponents; the 'virtual connection' line should be removed or clearly labeled as purely illustrative, with no quantitative weight.","section":"Fig. 3"},{"comment":"No confidence bands or parameter covariances are reported for the VIM fits. With only 15 data points and 3 free parameters, the quoted differences in alpha among profiles would be more convincing with a likelihood contour or covariance matrix.","section":"Sec. 3.5"},{"comment":"The Burkert profile returns alpha = 1.0000 exactly, which suggests a boundary hit; state whether bounds were imposed and whether the boundary is saturated.","section":"Table 1"},{"comment":"The 'small intrinsic scatter' is claimed as explained, but no scatter measure or comparison is provided. Quantify the scatter predicted by the model, or soften the claim.","section":"Sec. 5, first bullet"}],"recommendation":"major_revision","confidential_remarks":"The paper has a useful empirical component, but the central claim is currently not supported: the VIM's interaction term has wrong physical dimensions, and the CIA's '2VT link' reduces to a fitted constant normalized to a_M. I recommend major revision, asking the authors to refit a corrected VIM and to either compute R and g_D from explicit baryonic/DM profiles or explicitly reposition the paper as empirical RAR fitting. If the dimensional error cannot be resolved, the VIM portion should be removed from the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nThe stress-test is correct: Eq. A34 isn't an acceleration. I worked through the dimensions. W_BD is an energy, so its r-derivative is a force; dividing by 4πr^2 M_B leaves L^-1 T^-2, not L T^-2. Eq. A34 gives the same. So the VIM's g_BD has units of m^-1 s^-2. That means the weak-lensing fits in Fig. 4 are comparing a non-acceleration quantity to acceleration data. It's not a modeling choice; it's an internal inconsistency. Since the VIM is explicitly introduced to cover the weak-lensing regime, the paper's main claim is broken.\n\nWhat does survive? The CIA, which adds a constant gint to the baryonic acceleration, is a simple phenomenological formula. It genuinely reproduces the RAR for early-type galaxies, BCGs, and clusters, with different fitted G values. The authors show G varies with scale—low for dwarf spheroidals, higher for clusters. That empirical trend may be a useful benchmark, but it's a fit, not a prediction.\n\nThe problems are threefold. First, the CIA's gint is not derived from the 2VT; the paper never computes R or gD(rB) from actual baryonic and halo structure. The product is just fitted, so the virial theorem connection is decorative. Second, the characteristic scale aM is inserted via gint = G aM, then called 'emergent'; that's circular. Third, no code or data is provided, and there's no comparison with ΛCDM simulations, which would be the natural test. The paper itself admits several of these limitations, which is honest, but they don't go away.\n\nSo my verdict is less optimistic than the reader's. The VIM is not salvageable as written; the CIA is a useful fitting formula with an overreaching interpretation. I would not send this to peer review as is. A revised version that either fixes the VIM's units (which may not be possible without changing its physical meaning) or drops the VIM and recasts the CIA as an empirical parametrization could be a minor contribution. But the current manuscript does not support the claim of a physically motivated framework.\n\nBest,\n[You]","headline":"The VIM's 'local interaction acceleration' has the wrong units (m^-1 s^-2), so the paper's central claim fails; the CIA is a useful fitting formula but not a derived physical model.","tokens_in":21516,"tokens_out":4622,"would_cite":false,"duration_ms":48845,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that the radial acceleration relation in galaxies, clusters, and weak-lensing data emerges from standard dark matter via a two-component virial theorem with a scale-dependent baryon–dark-matter coupling.","keywords":["radial acceleration relation","two-component virial theorem","dark matter halos","baryon-dark matter coupling","galaxy clusters","weak gravitational lensing","acceleration scale","ΛCDM"],"falsifier":"A direct measurement of the dark-matter contribution across the baryonic region of an individual galaxy — for example, from resolved rotation-curve decomposition — showing that the dark-matter acceleration changes by more than a factor of a few where the baryonic acceleration changes by orders of magnitude would break the constant-interaction approximation. Likewise, weak-lensing RAR data that continue to flatten rather than turn downward below 10^-14 m/s^2 would falsify the VIM's cuspy-halo prediction.","tokens_in":20527,"feed_emoji":"🌌","tokens_out":4866,"duration_ms":49673,"temperature":0.7,"pith_summary":"The paper is trying to establish that the radial acceleration relation (RAR) — the tight correlation between observed and baryon-only acceleration — can be reproduced across galaxies, brightest cluster galaxies, galaxy clusters, and weak-lensing measurements without modified gravity or exotic baryon–dark-matter physics. The proposed mechanism is the two-component virial theorem (2VT), which adds an interaction term to the usual virial relation for a baryonic system embedded in a dark-matter halo. At galaxy and cluster scales, the dark-matter contribution is approximated as a constant additive acceleration (the CIA model), which yields the RAR's small scatter, its characteristic acceleration scale, and the predictive role of baryons. At the very low accelerations probed by weak-lensing data, that constant approximation fails and a radially resolved version (the VIM) takes over, fitting the data with a standard dark-matter density profile. If right, the RAR is an emergent, environment-dependent outcome of ΛCDM rather than evidence for new gravity.","feed_headline":"Dark matter term reproduces the radial acceleration relation","feed_subtitle":"One scale-dependent coupling matches galaxies, clusters, and weak-lensing data without modified gravity","key_machinery":"The two-component virial theorem (2VT), extended from the scalar virial theorem, defines a global acceleration scale g2VT(rB) = gvir,B(rB) + R gD(rB) at the baryonic gravitational radius, where R is a dimensionless geometrical parameter weighting the dark-matter contribution by the baryonic mass distribution. The CIA promotes the R gD term to a per-system constant gint = G aM, giving a one-parameter additive model; the VIM derives a local interaction acceleration from the radial derivative of the baryon–dark-matter interaction energy, producing a radius-dependent term that depends on the chosen dark-matter density profile. The 2VT sets the global scale, the CIA supplies the leading-order con","core_discovery":"The central claim is that the observed acceleration in a baryonic system embedded in a dark-matter halo is the Newtonian baryonic acceleration plus a geometric, baryon-weighted dark-matter contribution, and that this sum, evaluated either as a global constant or as a local radial term, reproduces the RAR across many orders of magnitude in acceleration. The constant-interaction model writes gobs = gB + gint, with gint = G aM, and fits early-type galaxies, late-type galaxies, dwarf spheroidals, brightest cluster galaxies, and intracluster regions; the virial-motivated interaction model replaces gint with a local term that depends on the enclosed dark-matter mass and local dark-matter density,","pith_inferences":["If G tracks the halo's central density and concentration, then at fixed baryonic mass, galaxies with more concentrated halos should sit systematically above the mean RAR; this correlation could be searched for in large samples with independent concentration estimates.","The framework implies that the RAR is not strictly universal: systems in dense environments should show a different normalization than field galaxies, a prediction testable with cluster weak-lensing and satellite-galaxy kinematics.","The VIM's low-acceleration downturn, if confirmed, would discriminate against scale-free modified-gravity formulations that predict continued flattening, since the downturn here stems from the local structure of the dark-matter halo."],"forward_implications":["The RAR's small scatter follows because, within each system, the dark-matter contribution is nearly constant, so the baryonic acceleration alone predicts the total acceleration up to a fixed offset.","The characteristic acceleration scale is not a fundamental constant in this framework; it emerges as the normalization of the baryon–dark-matter coupling and can differ between system types.","The cluster RAR's slope-1/2 power law is an emergent projection of a family of CIA curves with scale-dependent G, not a change in the gravitational law.","The VIM predicts a mild downturn in the RAR below roughly 10^-14 m/s^2 for a cuspy halo profile, offering a nontrivial test that future weak-lensing surveys can confirm or refute.","G varies with environment: it is higher in brightest cluster galaxies and clusters than in isolated galaxies, and decreases with radius inside clusters, so the RAR carries information about the dark-matter halo rather than a universal law."],"fun_headline_variants":["One dark-matter term unifies galaxy and cluster acceleration data","Dark matter coupling reproduces radial acceleration relation across scales","Single dark matter term explains acceleration in galaxies and clusters","Virial framework reproduces radial acceleration without modified gravity"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central assumption is that the dark-matter contribution to the observed acceleration is nearly constant across the radii sampled within each galaxy, backed by an approximation of constant dark-matter density inside the baryonic region; if the real halo contribution varies strongly there, the CIA's single-constant fits become fitting artifacts rather than a test of the two-component virial theorem.","fun_headline_variants_meta":{"raw":{"variants":["One dark-matter term unifies galaxy and cluster acceleration data","Dark matter coupling reproduces radial acceleration relation across scales","Single dark matter term explains acceleration in galaxies and clusters","Virial framework reproduces radial acceleration without modified gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000793,"raw_usage":{"total_tokens":3382,"prompt_tokens":851,"completion_tokens":2531,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":2475}},"tokens_in":595,"tokens_out":2531,"duration_ms":16444,"temperature":1.0,"reasoning_tokens":2475,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:22:24.820678+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the dark-matter contribution across the baryonic region of an individual galaxy — for example, from resolved rotation-curve decomposition — showing that the dark-matter acceleration changes by more than a factor of a few where the baryonic acceleration changes by orders of magnitude would break the constant-interaction approximation. Likewise, weak-lensing RAR data that continue to flatten rather than turn downward below 10^-14 m/s^2 would falsify the VIM's cuspy-halo prediction.","supporting_citations":[],"review_version":1}