{"id":"2abb4080-73c3-4f97-b267-731b5961277f","arxiv_id":"2507.04392","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"LISA could detect roughly 10^4 intermediate-mass black hole mergers in four years and AMIGO 10^2 to 10^3 in three years, with linear f(R,T) gravity altering these numbers by less than five percent at allowed couplings.","lead":"Space-based gravitational-wave detectors could see thousands of intermediate-mass black hole mergers: LISA up to about 10,000 over four years and the proposed AMIGO mission hundreds to thousands over three years. A simple modified-gravity extension, f(R,T) = R + 2 lambda T, changes these forecasts by only a few percent, so the paper frames the missions as a weak consistency check on modified gravity and a planning guide for the detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline event counts do not follow from the population model: Section 3.6 substitutes a constant R0=100 and hand-calibrated Gaussian detection functions for the Section 2 rate, so the 1e4/1e2-1e3 forecasts and the f(R,T) percentages rest on an external normalization.","rationale":"The reader's weakest_assumption identifies the same load-bearing weakness: the claim that LISA would see ~10^4 IMBH mergers and AMIGO ~10^2-10^3 is produced by multiplying a constant, externally chosen R0 by ad hoc Gaussian detection probabilities rather than by evaluating the population model developed in Section 2. I agree with that reading, and the manuscript's own text supports it: Section 3.6 explicitly calls the Gaussians parametric approximations, and no equation in the paper relates R0=100 to the derived rate in Eq. (2.6) or to the scenario parameters in Section 2.4. The linear dependence of all quoted counts on R0 means the quantitative conclusions are unsecured unless the integration is actually performed. The f(R,T) portion is also affected: Eqs. (4.10)-(4.11) stretch the same mock Gaussians by α(λ), so the claimed sub-percent redshift shift and sub-five-percent event-count change are relative statements about those mock functions, not about a full SNR recomputation. The qualitative statement that a tiny coupling produces a tiny effect is plausible, but the paper's central quantitative claims are not supported by the model as presented. This does not change the reader's rejection; it reinforces it.","tokens_in":14105,"tokens_out":7056,"duration_ms":79384,"concrete_test":"Recompute the LISA and AMIGO event counts from Eq. (3.1), using the intrinsic rate R(M1,M2,z) of Eq. (2.6) with the fiducial parameters of Section 2.4 and Pdet = 1 when the SNR from Eq. (3.3) exceeds 8, evaluated with the published LISA and AMIGO noise curves of Sections 3.1-3.3. Integrate over M1 ∈ [10^2,10^5] M_sun, q ∈ [0.1,1], z ∈ [0,10], and compare with the counts obtained from Eq. (3.4) with R0=100. If the two totals differ by more than a factor of two for either detector, the headline event-rate forecasts do not follow from the model. Repeating the same integral with R0=10 and R0=1000 quantifies the linear sensitivity of the headline numbers.","verdict_should_be":"REJECT","load_bearing_attack":"Section 3.6, Eq. (3.4), computes the detectable redshift distribution from a constant comoving merger rate density R0 = 100 Gpc^-3 yr^-1 and hand-calibrated Gaussian detection probabilities (3.6)-(3.7), explicitly described as 'mock functions' rather than derived detector responses. This is the step that produces the abstract's ~10^4 LISA and 10^2-10^3 AMIGO counts; Section 3.5's reference to the Section 2 mass function is not backed by a shown integration. The Section 2 model does contain an intrinsic rate built from Φ(M1,z)Φ(M2,z) fbin ∫ p(a0)/τ da0 (Eq. 2.6), with scenario parameters spanning fbin = 0.001-0.05 and ηmerge = 0.01-0.5, but no derivation links R0=100 to those parameters. Since Ndet is linear in R0, any factor-10 error in the assumed rate changes every headline number by a factor of 10. The f(R,T) claims inherit the same weakness: Eqs. (4.10)-(4.11) rescale the Gaussian centers and widths by α(λ) instead of recomputing SNR, so a mis-modeled Gaussian would mis-model the modified-gravity percentages even if the α rescaling itself is correct. Thus the load-bearing quantitative content of the paper rests on an external rate normalization, not on the population model advertised in Section 2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a semi-analytic forecast for the detection of intermediate-mass black hole binaries with the space-based detectors LISA and AMIGO, based on an extrapolated IMBH mass function, and extends the analysis to the minimal matter-coupled modified gravity model f(R,T)=R+2λT. The abstract reports event rates of about 10^4 mergers for LISA in four years and 10^2–10^3 for AMIGO in three years, and claims that for |λ|≤2×10^-3 the peak detection redshift shifts by less than one percent and the total event number changes by less than five percent. The paper identifies optimal detection bands in mass and redshift and discusses complementarity between the two detectors.","tokens_in":14575,"tokens_out":10309,"duration_ms":96919,"significance":"If the quantitative forecasts were reliable, the paper would provide useful mission-planning information for LISA and AMIGO and a simple population-level test of f(R,T) gravity. The authors are candid that the redshift-dependent detection probabilities are 'mock functions' intended as illustrative tools rather than derived detector responses. However, the central event-rate numbers and the modified-gravity percentages are not computed from the mass-function model advertised in Section 2, and the f(R,T) amplitude rescaling contains algebraic and dimensional errors. The qualitative conclusion that LISA and AMIGO are complementary is plausible, but the specific quantitative claims require a substantial reworking.","major_comments":[{"comment":"The differential event rate dNbin/dz is computed from an assumed constant comoving rate density R0 = 100 Gpc^-3 yr^-1 and Gaussian detection probabilities Pdet(z) that the text explicitly calls 'mock functions' and 'illustrative tools rather than physically rigorous detection probabilities.' No derivation connects R0 to the mass-function model of Section 2; the scenario parameters fbin, eta_merge, and <a0> defined in Section 2.4 do not appear in Eq. (3.4), and the paper states that the Gaussians were 'chosen to reflect' the full SNR maps rather than derived from them. Since Ndet is linear in R0, the headline event counts in the abstract and Section 3.5 inherit this external normalization, and the f(R,T) percentages computed in Section 4.3 rescale these same Gaussians. The authors should either compute R0 from the binary mass function or reframe the forecasts as explicitly conditional on an assumed rate density.","section":"Section 3.6, Eq. (3.4)"},{"comment":"The binary mass function and coalescence rate are defined in three mutually inconsistent forms. Eq. (2.3) gives dNbin/(dM1 dM2 dz dV) = Φ(M1) fbin p(q) with no 1/M1 Jacobian, Eq. (2.6) instead uses a product Φ(M1)Φ(M2) and an integral over p(a0)/τmerge, and Eq. (3.8) reintroduces the 1/M1 factor. The merger efficiency ηmerge introduced in Section 2.3 and used in the scenario definitions of Section 2.4 does not appear in Eq. (2.6). The paper does not state which expression is inserted into the numerical integration that produces the event counts, so the quantitative results are not reproducible.","section":"Section 2.2-2.3 and 3.7, Eqs. (2.3), (2.6), (3.8)"},{"comment":"The working field equation (4.4), Gµν = 8πG Tµν + 2λ Tµν + λ T gµν, does not follow from Eq. (4.3) when Θµν = -2Tµν for dust; the substitution gives Gµν = (8πG + 2λ)Tµν with no λTgµν term. In addition, the Newtonian limit Eq. (4.5), Geff = G(1 - λ/8π), is dimensionally inconsistent unless G is set to unity and has a sign opposite to what follows from Eq. (4.4) with a (-,+,+,+) signature. Because the f(R,T) strain rescaling α(λ) = (1 - λ/8π)^(5/6) used throughout Section 4 is derived from this Geff, the modified-gravity conclusions are founded on an incorrect derivation.","section":"Section 4.1, Eqs. (4.3)-(4.5)"},{"comment":"The transformation (z0, σ) → (α z0, α σ) is not the correct transformation of a detection probability that was defined by an SNR threshold. The SNR is not proportional to 1/DL(z) ∝ 1/z at all redshifts, since the redshifted chirp mass and the detector noise enter the SNR integral nontrivially, and rescaling a Gaussian's center and width does not reproduce the effect of multiplying the SNR by α. The paper asserts that these closed forms 'reproduce the full numerical SNR rescaling to better than 2%' but gives no numerical comparison or derivation; this claim needs to be substantiated or removed.","section":"Section 4.3, Eqs. (4.9)-(4.11)"}],"minor_comments":[{"comment":"The event-rate prediction is given only as 'preliminary results suggest', with no indication of which mass-function version from Section 2 was integrated or over what domain; please specify the exact computation.","section":"Section 3.5"},{"comment":"The text writes ln(10^3) ≈ 6.9 but the displayed expression uses log(10^5/10^2); clarify the base or use a single notation.","section":"Section 3.7, Eq. (3.10)"},{"comment":"The statement that GW observatories could verify 'a difference in the propagation speed of waves with different frequencies' is not supported by the model, which has no massive mode at linear order; either connect this to a concrete observable or remove it.","section":"Section 4.4"},{"comment":"These references are cited to support the use of Gaussian mock detectability functions, but those references concern population synthesis and hierarchical mergers, not this specific approximation; please provide a more direct citation or note that the Gaussian is a simplifying assumption.","section":"References [41] and [42]"},{"comment":"The paper claims to be the 'first side-by-side forecast' of IMBH mergers for LISA and AMIGO, but does not compare with existing IMBH rate estimates (e.g., Sesana 2007; Fragione et al. 2018); a brief quantitative comparison would help the reader calibrate the results.","section":"Section 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely problem and the qualitative claim of LISA–AMIGO complementarity is plausible. However, the quantitative forecasts and the f(R,T) perturbation analysis need to be reworked from a consistent model. The authors should also verify the field equations in Section 4.1 and provide a valid derivation of the detection-probability rescaling. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a readable side-by-side forecast of IMBH merger detection for LISA and AMIGO, with a detour into f(R,T) gravity. What's new is the package: the two detectors in one population framework, with event counts and redshift distributions. The qualitative conclusion that Solar-System-compatible λ in linear f(R,T) barely changes detection statistics is plausible and likely robust.\n\nThe problem is that the headline numbers don't come from the population model. Section 2 builds a redshift-dependent binary mass function with pairing fractions and merger timescales, but Section 3.6 switches to a constant comoving rate R0=100 Gpc^-3 yr^-1 and hand-calibrated Gaussian detection probabilities, explicitly called 'mock functions'. The abstract's 10^4 LISA events and 10^2-10^3 AMIGO events are linear in R0. No derivation links R0 to the Section 2 model; the paper doesn't show the numerical integration it refers to in Section 3.5. So the central quantitative claims rest on an external assumption, not on the advertised model.\n\nThere are also internal inconsistencies: the binary mass function appears in three forms (Eqs. 2.3, 2.6, 3.8) that differ by a Jacobian and by whether they include both Phi(M1) and Phi(M2). The f(R,T) part rescales Gaussian centers and widths by alpha(lambda) and claims better-than-2% agreement with a full SNR calculation, but that agreement is asserted, not shown.\n\nWhat the paper does well: it is honest—the mock-function caveat is stated plainly—and the general picture of LISA/AMIGO complementarity is correct. The f(R,T) strain rescaling factor (1 - lambda/8pi)^(5/6) is a standard result and the hierarchy-of-corrections discussion is sensible.\n\nWho is this for? Mission planners and people interested in IMBH population forecasts. It's not a methods paper. As it stands, I wouldn't trust the specific numbers for planning. A revision that either derives R0 from the Section 2 model or clearly frames the paper as a parameter study around an assumed R0 would be much more useful. If the authors did that, this could be a reasonable JCAP contribution.\n\nFor peer review: I'd send it to referees rather than desk-reject, because the topic is timely and the flaws are fixable. But I'd expect a major revision before it's publishable.\n\nBest","headline":"A readable side-by-side LISA/AMIGO IMBH forecast whose headline event counts rest on an assumed rate and mock detection functions, not on the Section 2 population model; the f(R,T) conclusion is plausible but the quantitative claims are not yet supported.","tokens_in":15069,"tokens_out":2937,"would_cite":false,"duration_ms":31972,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83D05","83C57"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper forecasts that LISA will detect about $10^4$ intermediate-mass black hole binary mergers in four years and AMIGO hundreds to thousands in three, with minimal f(R,T) gravity changing those totals by under five percent.","keywords":["intermediate-mass black holes","gravitational wave astronomy","LISA","AMIGO","detection band forecasts","f(R,T) gravity","modified gravity","event rates"],"falsifier":"A redshift-resolved catalog from LISA's first four years would settle it: the model predicts a peak in $dN/dz$ near $z \\approx 4$--$5$ with about $10^4$ total events, so an observed peak shifted by more than roughly $\\Delta z \\approx 1$, a total count differing from $10^4$ by more than about a factor of two or three, or a LISA/AMIGO count ratio that deviates by more than $5\\%$ from the predicted ratio would contradict the fiducial model and the f(R,T) claim.","tokens_in":13890,"feed_emoji":"🌌","tokens_out":15405,"duration_ms":137845,"temperature":0.7,"pith_summary":"The paper tries to establish that intermediate-mass black hole binaries — pairs of black holes between roughly $10^2$ and $10^5$ solar masses — will be a high-yield target for the space-based gravitational wave observatories LISA and AMIGO. It builds a redshift-dependent population model from extrapolated black hole mass functions, flat mass-ratio pairing, and a distribution of orbital separations, then folds in the standard gravitational-wave merger timescale and each detector's sensitivity curve to map the detectable mass–redshift bands. The central numbers are about $10^4$ LISA mergers in four observing years and between $10^2$ and $10^3$ AMIGO mergers in three, with LISA seeing heavier systems at higher redshift and AMIGO lighter systems nearby. The paper then asks whether the minimal matter-coupled modified gravity model $f(R,T) = R + 2\\lambda T$ changes this picture; it finds that for couplings compatible with Solar System tests, $|\\lambda| \\le 2\\times10^{-3}$, the peak detection redshift shifts by less than one percent and the total event count by less than five percent.","feed_headline":"LISA could net 10,000 black-hole mergers in four years","feed_subtitle":"The forecast pairs LISA's millihertz band with AMIGO's decihertz window on intermediate-mass black holes.","key_machinery":"The carrying machinery is the semi-analytic redshift-dependent binary mass function of Section 2 together with the analytic rescaling identity of Section 4. The binary mass function $\\Phi(M_1,z)\\,f_{\\rm bin}(M_1,M_2,z)\\,p(q)$ starts from an extrapolated IMBH mass function, a flat mass-ratio distribution $p(q)\\propto q^0$, and a logarithmically flat distribution of initial orbital separations integrated against the standard gravitational-wave inspiral timescale; this produces the intrinsic coalescence rate. For forecasting, the detector response is compressed into Gaussian redshift profiles $P_{\\rm det}(z)$, and the modified-gravity part is carried entirely by the factor $\\alpha(\\lambda)=(1-\\lambda/8\\pi)^{5/6}$, which rescales the Gaussian centers and widths. That identity converts the f(R,T) strain change into a transparent shift of the event-rate forecast without recomputing signal-to-noise ratios, and is stated to reproduce the full numerical SNR rescaling to better than $2\\%$.","core_discovery":"The central claim is that intermediate-mass black hole binaries are a rich multi-detector population. Starting from a constant comoving merger rate density $R_0 = 100\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$ and Gaussian detection probabilities $P^{\\rm LISA}_{\\rm det}(z) = \\exp[-(z-5)^2/(2.5)^2]$ and $P^{\\rm AMIGO}_{\\rm det}(z) = \\exp[-(z-1.5)^2/(0.75)^2]$, multiplied by the $\\Lambda$CDM comoving volume and a $(1+z)^{-1}$ time-dilation factor, the model yields about $10^4$ events for LISA in four years and $10^2$–$10^3$ for AMIGO in three years. In the linear model $f(R,T) = R + 2\\lambda T$, the gravitational-wave strain is rescaled by $\\alpha(\\lambda) = (1 - \\lambda/8\\pi)^{5/6}$ while the phase is unchanged at leading order, so the detection-probability Gaussians are shifted from $(z_0,\\sigma)$ to $(\\alpha z_0,\\alpha\\sigma)$. For $|\\lambda| \\le 2\\times10^{-3}$, this rescaling moves the peak detection redshift by less than $1\\%$ and the integrated event count by less than $5\\%$; the paper concludes that, at population level, this minimal matter-coupled modified gravity is observationally indistinguishable from general relativity and can serve mainly as a consistency check.","pith_inferences":["A natural extension is to replace the constant $R_0 = 100\\,\\mathrm{Gpc}^{-3}\\,\\mathrm{yr}^{-1}$ with the redshift- and mass-dependent merger rate derived from the paper's own mass-function model, which would give the event counts a fully internal normalization rather than an external input.","Because the f(R,T) forecast is made by rescaling the Gaussian detection functions rather than recomputing SNR maps, a full numerical rerun with realistic noise curves could find larger differences near the detection threshold, particularly for AMIGO's narrower band.","If both detectors fly, the ratio of LISA to AMIGO event counts becomes a useful test that partially cancels population uncertainties, since many astrophysical factors affect both instruments the same way.","The same rescaling technique could be extended to nonlinear f(R,T) variants that are equipped with screening, where population-level shifts of order ten percent are plausible even while Solar System bounds are satisfied."],"forward_implications":["If LISA runs for four years as assumed, it should see on the order of $10^4$ IMBH binary mergers, making intermediate-mass black holes a statistically rich population rather than a handful of special events.","AMIGO would open a complementary decihertz window on lighter binaries at $z \\sim 1.5$--$2$, and together with LISA the two observatories would cover total masses from about $10^2$ to $10^6\\,M_\\odot$ continuously.","Within current Solar System bounds, $f(R,T)=R+2\\lambda T$ cannot be distinguished from general relativity by IMBH event counts; the GR forecasts therefore remain the reference numbers for both detectors.","The amplitude rescaling $\\alpha(\\lambda)$ provides a consistency test: after marginalizing astrophysical uncertainties, a systematic few-percent offset between the predicted and observed merger rates would point to a varying effective gravitational constant or a similar uniform rescaling.","The estimate of about $6\\times10^6$ IMBH binaries in the observable universe sets an upper bound on the total merger yield available to any future detector."],"supporting_citations":[{"why":"defines the f(R,T) action and field equations from which the linear $R+2\\lambda T$ model is taken.","marker":"[19]"},{"why":"provides the AMIGO mission concept and decihertz sensitivity curve used in the AMIGO forecasts.","marker":"[12]"},{"why":"provides the LISA mission sensitivity curve used in the detectability maps and event-rate integrals.","marker":"[13]"},{"why":"supplies the analytic LISA noise power spectral density used for SNR evaluation.","marker":"[40]"},{"why":"gives the gravitational-wave merger timescale for circular binaries that converts the binary mass function into a coalescence rate.","marker":"[36]"},{"why":"derives the weak-field limit of f(R,T) gravity, giving the effective gravitational constant and the Solar System bound on the coupling.","marker":"[43]"},{"why":"supports the leading-order rescaling of gravitational-wave strain in f(R,T) gravity used in the forecast.","marker":"[47]"},{"why":"provides the Solar System bound on $|\\Delta G/G|$ that sets the allowed range of the coupling $\\lambda$.","marker":"[51]"}],"fun_headline_variants":["LISA forecast: ~10,000 IMBH mergers in four years","AMIGO adds decihertz window to IMBH detection","Modified gravity barely changes IMBH event rates","LISA + AMIGO: complementary bands for IMBH binaries","10,000 IMBH mergers for LISA, hundreds for AMIGO"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline numbers assume a constant comoving merger rate of 100 intermediate-mass black hole binaries per cubic gigaparsec per year at every redshift, and model each detector's sensitivity as a fixed bell-shaped curve in redshift; if the true rate is ten times smaller or larger, or the sensitivity shape differs, every headline count changes by that factor.","fun_headline_variants_meta":{"raw":{"variants":["LISA forecast: ~10,000 IMBH mergers in four years","AMIGO adds decihertz window to IMBH detection","Modified gravity barely changes IMBH event rates","LISA + AMIGO: complementary bands for IMBH binaries","10,000 IMBH mergers for LISA, hundreds for AMIGO"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000364,"raw_usage":{"total_tokens":2047,"prompt_tokens":1115,"completion_tokens":932,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":731,"completion_tokens_details":{"reasoning_tokens":849}},"tokens_in":731,"tokens_out":932,"duration_ms":8630,"temperature":1.0,"reasoning_tokens":849,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:51:38.291954+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A redshift-resolved catalog from LISA's first four years would settle it: the model predicts a peak in $dN/dz$ near $z \\approx 4$--$5$ with about $10^4$ total events, so an observed peak shifted by more than roughly $\\Delta z \\approx 1$, a total count differing from $10^4$ by more than about a factor of two or three, or a LISA/AMIGO count ratio that deviates by more than $5\\%$ from the predicted ratio would contradict the fiducial model and the f(R,T) claim.","supporting_citations":[{"cited_title":"Harko, F.S.N","cited_arxiv_id":null,"evidence_quote":"defines the f(R,T) action and field equations from which the linear $R+2\\lambda T$ model is taken."},{"cited_title":"Ni,Gravitational wave detection in space, Int","cited_arxiv_id":null,"evidence_quote":"provides the AMIGO mission concept and decihertz sensitivity curve used in the AMIGO forecasts."},{"cited_title":"Peters,Gravitational radiation and the motion of two point masses, Phys","cited_arxiv_id":null,"evidence_quote":"gives the gravitational-wave merger timescale for circular binaries that converts the binary mass function into a coalescence rate."},{"cited_title":"Alvarenga, A","cited_arxiv_id":null,"evidence_quote":"derives the weak-field limit of f(R,T) gravity, giving the effective gravitational constant and the Solar System bound on the coupling."},{"cited_title":"Sharma, A.K","cited_arxiv_id":null,"evidence_quote":"supports the leading-order rescaling of gravitational-wave strain in f(R,T) gravity used in the forecast."},{"cited_title":"Bertotti, L","cited_arxiv_id":null,"evidence_quote":"provides the Solar System bound on $|\\Delta G/G|$ that sets the allowed range of the coupling $\\lambda$."}],"review_version":1}