{"id":"8f400a70-0c97-4d26-84ed-889d62169684","arxiv_id":"2509.03364","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulated ELT/HARMONI observations of VCC 1861 show that a 2.5e5 solar mass IMBH (5% of the nuclear cluster mass) can be recovered from stellar kinematics at Virgo distance.","lead":"This paper simulates mock ELT/HARMONI observations of a faint dwarf galaxy in the Virgo Cluster to test whether an intermediate-mass black hole weighing 5% of the nuclear star cluster can be seen. The simulations recover the injected black hole mass, arguing the ELT will be able to probe such black holes at 16.5 Mpc.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Extrapolated central surface brightness profile is the key unvalidated input; alternative core-Sersic slopes must be tested before the 5% IMBH detectability claim is accepted.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the HST-to-ELT extrapolation of the central surface brightness profile with gamma=0.1 is unvalidated and directly controls the BH detectability. The paper's central claim, that ELT/HARMONI can measure M_BH=2.5e5 Msun (5% of the NSC mass) in Virgo Cluster dwarfs, rests on mock observations built from this extrapolated model. The simulation is internally consistent: the MGE is fitted, the JAMcyl model generates the noiseless cube, HSIM adds realistic noise and PSF, pPXF extracts kinematics, and adamet recovers the injected BH mass with small uncertainties. This is a valid demonstration of pipeline precision, but it does not by itself establish robustness to model misspecification. The most direct missing test is to vary the central slope gamma and see whether the recovery remains unbiased and the signal remains detectable. This is not a fatal flaw in the pipeline; it is a missing robustness check on the central astrophysical input. Since the reader already returned CONDITIONAL and the concern matches, the verdict should remain UNCHANGED. I would caution that if a future test with gamma=0.5 or 1.0 shows large biases, the paper's conclusion would need to be downgraded.","tokens_in":14734,"tokens_out":5468,"duration_ms":59857,"concrete_test":"Rerun the simulation with alternative core-Sersic slopes gamma = 0.0, 0.5, and 1.0. For each gamma, re-fit the HST surface brightness profile, construct a new MGE, regenerate the M_BH=2.5e5 and no-BH mock cubes at the same HARMONI H-high settings, and repeat the pPXF + adamet recovery. If any alternative gamma yields a median log M_BH differing by more than ~0.3 dex from the input, or a 3-sigma interval that includes log M_BH = 0 (i.e., no detection), then the claimed detectability at 5% of NSC mass is not robust to plausible central density profiles. Report the recovered log M_BH and its 3-sigma uncertainty for each gamma.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the stellar density profile at radii below HST resolution. In Section 3.1, the HST/ACS image (FWHM 0.11 arcsec, pixel 0.05 arcsec) is fitted with a core-Sersic profile with the power-law index fixed at gamma=0.1, 'as assumed by [41]', then extrapolated 'deeply toward the center' to the ELT/HARMONI pixel scale. The resulting MGE (Table 1) has a central Gaussian with sigma=0.008 arcsec, about 10 times smaller than the HST PSF and comparable to the planned 10 mas spaxels. This extrapolation sets the central stellar density, which directly controls the contrast between the M_BH=2.5e5 Msun and no-BH kinematics shown in Section 4.1. If the true NSC cusp has a different gamma (e.g., 0.5 or 1.0), the stellar mass distribution inside the sphere of influence (R_soi ~17 mas) changes, altering the recovered BH mass and the detectability threshold. The paper explicitly acknowledges the need to extrapolate but does not test alternative slopes or quantify how the central MGE component is constrained by the HST data. Because the mock data are generated and fitted with the same JAMcyl model family, the exercise is a closed-loop consistency check that inherits whatever bias the assumed gamma introduces. Thus the feasibility forecast is conditional on an unmeasured structural parameter.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a simulation study assessing whether ELT/HARMONI high-resolution integral-field stellar kinematics can measure intermediate-mass black holes (IMBHs) in Virgo Cluster dwarf galaxies. The authors build a stellar mass model of VCC 1861 from HST/ACS imaging, extrapolated to sub-HST scales with a core-Sersic profile (γ=0.1), construct noiseless mock cubes using JAMcyl, pass them through the HSIM instrument simulator for the H-high and K-long high-resolution gratings with LTAO, extract kinematics with Voronoi binning and pPXF, and recover the BH mass via JAMcyl+MCMC. For an injected IMBH of log M_BH=5.4, they recover 5.38 (H-high) and 5.34 (K-long); for the no-BH case they place upper limits. They conclude that ELT/HARMONI can detect IMBHs with masses about 5% of the host nuclear star cluster mass at 16.5 Mpc.","tokens_in":15096,"tokens_out":7967,"duration_ms":76977,"significance":"If the central claim holds, the paper provides a valuable feasibility forecast: extending IMBH dynamical mass measurements from the Local Volume (≲3.5 Mpc) to the Virgo Cluster distance would dramatically increase the sample of accessible IMBH candidates. The study uses a realistic instrument simulation pipeline (HSIM, LTAO, Voronoi binning, pPXF with a 13-template library) and the recovery of the injected BH mass is clean. However, the test is a closed-loop exercise—the same JAMcyl model family, the same MGE, and the same template library are used for both mock generation and fitting—and the central stellar density profile is extrapolated below HST resolution using a fixed, untested core-Sersic slope. These two issues make the headline detectability claim conditional on unvalidated structural assumptions. The paper also offers a concrete, falsifiable prediction about HARMONI's capabilities, which is useful for planning ELT observations.","major_comments":[{"comment":"The central MGE component has σ=0.008″ (Table 1), about 10 times smaller than the HST/ACS PSF FWHM (0.11″). It is derived from a core-Sersic fit with γ=0.1 fixed 'as assumed by [41]' and extrapolated below the HST resolution to the 10 mas spaxel scale. This component sets the stellar density inside the BH sphere of influence (R_soi≈17 mas) and directly controls the contrast between the BH and no-BH Vrms profiles that underlies the detectability claim in Sections 4.1 and 4.2. The authors acknowledge the need to 'extrapolate deeply toward the center' but do not test alternative γ values. If the actual NSC cusp is steeper (e.g., γ=0.5 or 1.0), the central stellar mass would be larger and could mimic or mask a 2.5×10^5 M⊙ BH, changing the 5% detectability threshold. Please add sensitivity tests with at least two alternative inner slopes (or a non-parametric MGE) and quantify the effect on th","section":"Section 3.1"},{"comment":"The mock data are generated with JAMcyl using the same MGE, β_z=0.02, M/L=1.4, and the same MARCS library that are used in the pPXF extraction and the JAMcyl fit. The recovery of log M_BH=5.38/5.34 is therefore primarily a check of internal consistency of the MCMC and the instrument pipeline. The claim of 'the ELT's capability to detect IMBHs' overstates the evidence: real galaxies have unknown anisotropy, possible M/L gradients, and unmeasured central density slopes. Please either (i) run a robustness test with mock data generated at different β_z and with a different MGE inner slope, verifying that the recovery remains unbiased, or (ii) soften the abstract and conclusion to state that the detectability is conditional on the assumed stellar mass profile. As written, this is a load-bearing caveat for the central claim.","section":"Sections 3.2 and 4.2"}],"minor_comments":[{"comment":"Affiliation contains the typo 'Physisc' for 'Physics'.","section":"Section 1"},{"comment":"The PSF is described as 'HST/WFPC2 F850LP'; the observations are ACS/WFC, so this should be corrected.","section":"Section 3.1"},{"comment":"Text says 'pixel size of 5 × 5 mas2. This pixel size is smaller than the 10 × 10 arcsec² of ELT/HARMONI'—the units should be '10 × 10 mas²'.","section":"Section 3.2"},{"comment":"The text refers to inclination spanning '300 − 900' and a 'best fit at 300'; these should be '30°–90°' and '30°' (or, given the input of 44° stated in Section 3.2, the statement 'aligning with our input parameters' is incorrect).","section":"Section 4.2"},{"comment":"The chi-squared equation is mangled: '∑ (V_rms,i − V̅_rms,i)^n' should be 'Σ (V_rms,i − V̅_rms,i)^2 / σ^2', and the meaning of σ in that expression should be clarified.","section":"Section 4.2"},{"comment":"The caption says 'medium resolution gratings,' but the paper uses the H-high and K-long high-resolution gratings; this appears to be a copy-paste error.","section":"Figure 5 caption"},{"comment":"The 'Input Parameters' column lists prior ranges, not the true input values. Please add a separate column (or row) with the injected values (i=44°, M/L=1.4, log M_BH=5.4) so the recovery can be directly compared.","section":"Table 2"},{"comment":"The statement 'β_z is also well-constrained with the difference of Δβ_z < 0.1' is not supported by Table 2: the no-BH K-long row reports β_z = 0.007 ± 0.240. This claim should be qualified or removed.","section":"Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The reader's stress-test concern about the extrapolated central surface brightness profile is valid and is the main driver of the major-revision recommendation. The paper is otherwise a reasonable feasibility study. In addition to the requested sensitivity tests, the authors should carefully proofread the manuscript: there are many typos, inconsistent inclination values (44° vs. 30°), and at least one figure caption that appears copied from another paper. The journal should also ensure that the table and figure content agree with the text. If the sensitivity tests confirm the detectability claim, the paper would be a useful contribution to the ELT/HARMONI planning literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper gets the mechanics right and shows that ELT/HARMONI can recover an injected IMBH in a Virgo dwarf under the assumptions they make. The new piece is honest: VCC 1861 at 16.5 Mpc, BH 2.5e5 Msun (5% of the NSC), both H-high and K-long. They build mock cubes with HSIM/LTAO, bin to S/N 20, extract kinematics with pPXF, and recover log MBH = 5.38/5.34 vs input 5.4. The sigma maps show a clear central peak for the BH case and a central drop for the no-BH case, consistent with expected Jeans behavior.\n\nWhere I would push back mirrors the stress-test concern. The central stellar density profile is inferred from a core-Sersic fit to HST data and extrapolated to 10 mas scales with gamma fixed at 0.1, taken 'as assumed by [41]'. That is a load-bearing assumption: it sets the stellar contribution inside the sphere of influence and therefore the contrast that makes the BH detectable. Alternative gamma values are not explored, and the paper doesn't demonstrate that HST alone pins gamma down. The authors explicitly note they need to extrapolate deeply toward the center, but they don't quantify the uncertainty this introduces.\n\nThe second issue is the closed-loop character. JAMcyl generates the cubes and JAMcyl fits them. The recovery is then a consistency check, not an independent validation of the dynamical model. I don't think this is a fatal flaw—the pipeline is realistic, and JAMcyl is a standard tool—but it limits the strength of the forecast. A non-JAM forward model (even axisymmetric Schwarzschild) or a statement about why that is unnecessary would make the claim robust.\n\nThird, no detection significance or completeness is given. The paper shows a kinematic signature but does not say how many spaxels/bin outliers could mimic it, nor what fraction of similar targets would show a measurable peak.\n\nMinor but easy: text says average inclination 44°, while best fits in Table 2 are 31-37° and the text says they align with input; that's not consistent. Also typos and figure captions that label high-resolution as medium.\n\nWho it's for: ELT/HARMONI planning and IMBH search strategy. It deserves serious refereeing; I'd ask for the gamma sensitivity test, a significance calculation, and code release before accepting.","headline":"A believable closed-loop recovery, but the detectability claim rests on an untested central cusp slope; deserves refereeing with conditions.","tokens_in":15648,"tokens_out":3179,"would_cite":true,"duration_ms":33413,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations indicate ELT/HARMONI can detect 2.5e5 solar-mass black holes in Virgo dwarfs.","keywords":["intermediate-mass black holes","nuclear star clusters","ELT/HARMONI","stellar kinematics","Jeans Anisotropic Modeling","Virgo Cluster","mock integral-field spectroscopy","black hole mass measurement"],"falsifier":"Re-running the same mock pipeline on VCC 1861 with a different inner power-law index (for example gamma=0.5 instead of 0.1) would test whether the central velocity-dispersion peak of the 2.5e5 solar-mass IMBH remains distinguishable from the no-black-hole case; if the profiles become indistinguishable, the central claim fails. Equivalently, a real 4.5-hour ELT/HARMONI observation of VCC 1861 showing no central sigma peak would contradict the predicted detectability.","tokens_in":14609,"feed_emoji":"🔭","tokens_out":9664,"duration_ms":90368,"temperature":0.7,"pith_summary":"Intermediate-mass black holes are the likely seeds of supermassive black holes, but almost none are confirmed beyond the Local Group because their gravitational spheres of influence are too small for current telescopes. This paper tries to establish that the forthcoming ELT/HARMONI can change that at Virgo Cluster distances of 16.5 Mpc. Using VCC 1861, a faint dwarf elliptical with a nuclear star cluster, the authors build a realistic mass model from HST imaging, generate mock HARMONI observations with and without a 2.5e5 solar-mass central black hole, and ask whether the black-hole mass can be recovered from the simulated stellar kinematics. They find that the high-resolution H-high and K-long gratings at 10 mas spaxels produce a clear central peak in velocity dispersion when the IMBH is present and a central drop when it is not, and that Bayesian recovery returns the input mass to within about 0.05-0.09 dex. If this holds, ELT/HARMONI will be able to measure IMBHs in Virgo dwarfs and begin mapping their occupation fraction, directly testing seeding scenarios for supermassive black holes.","feed_headline":"Mock ELT observations weigh a 250,000-solar-mass black hole","feed_subtitle":"HARMONI's high-resolution gratings resolve the central velocity peak of an IMBH in Virgo dwarf VCC 1861.","key_machinery":"The load-bearing observable is the inner velocity-dispersion profile: within the central roughly 0.1 arcsec, the black hole's gravity creates a sigma peak, whereas a cluster without a black hole shows a sigma drop. Whether that contrast survives depends on the central stellar density, which the authors extrapolate from 0.11 arcsec HST photometry to the 0.01 arcsec HARMONI pixel scale using a core-Sersic profile with a fixed inner power-law index gamma=0.1. The simulation chain then links a Multi-Gaussian Expansion mass model, JAMcyl anisotropic Jeans models, HSIM mock-cube generation, pPXF spectral fitting, and adaptive-Metropolis MCMC to recover masses. The main degeneracy is an anti-correl","core_discovery":"The paper reports that ELT/HARMONI's high-resolution H-high and K-long gratings at 10 mas spaxels can detect and measure a 2.5e5 solar-mass IMBH in the nuclear star cluster of VCC 1861 at 16.5 Mpc. The simulated black hole creates a central velocity-dispersion peak of about 24 km/s inside 0.1 arcsec, while the no-black-hole model produces a central drop, so the two cases are clearly distinguishable. Bayesian Jeans Anisotropic Modeling recovers log(M_BH/Msun)=5.382±0.049 in H-high and 5.337±0.090 in K-long against an input of 5.4, with the mass-to-light ratio near the input value of 1.4. The no-black-hole cubes yield upper limits of 3.18 and 2.54 in log(M_BH/Msun), so the same pipeline can bo","pith_inferences":["The same pipeline could be pointed at other Virgo dwarf galaxies with nuclear star clusters; a survey of detections and non-detections would measure the IMBH occupation fraction in dwarfs, which the paper does not attempt.","The fixed gamma=0.1 inner slope is the least tested input; if ELT resolves the true nuclear cusp, comparing the assumed and observed profiles would directly test the detectability calculation.","Because the field of view covers only the nuclear star cluster, inclination is poorly constrained; adding the galaxy's outer disk in future simulations could reduce the mass-to-light versus black-hole-mass degeneracy and tighten the recovery.","A central sigma-peak-versus-drop dichotomy, if generic, would make cheap IMBH screening possible from short observations before expensive full MCMC modeling."],"forward_implications":["ELT/HARMONI high-resolution gratings at 10 mas spaxels can distinguish a 2.5e5 solar-mass IMBH from a no-black-hole nuclear star cluster at 16.5 Mpc.","In the no-black-hole case the method returns tight upper limits, so it can also rule out IMBHs at the 5% mass contrast.","The H-high and K-long gratings yield consistent kinematics with differences below 5%, giving an internal cross-check on the recovered masses.","The recovered black-hole mass is accurate to about 0.05-0.09 dex, and the mass-to-light ratio is recovered near its input value.","This extends IMBH searches from the Local Group out to Virgo Cluster distance, substantially enlarging the volume over which IMBH demographics can be probed."],"supporting_citations":[{"why":"The earlier ELT/HARMONI IMBH simulation for Local Group galaxies that this work extends to Virgo distance.","marker":"[41]"},{"why":"Provides the stellar mass and age of VCC 1861's nuclear star cluster, setting the input cluster mass of 5.2e6 solar masses.","marker":"[46]"},{"why":"Source of the HST ACS/WFC Virgo survey image and nuclei photometry used for the surface brightness and mass model.","marker":"[51]"},{"why":"The HSIM pipeline used to generate the mock HARMONI integral-field data cubes.","marker":"[52]"},{"why":"The core-Sersic profile used to fit and extrapolate the central surface brightness to the ELT pixel scale.","marker":"[57]"},{"why":"The axisymmetric Jeans Anisotropic Modeling formalism used to construct the kinematic maps.","marker":"[65]"},{"why":"The JAMcyl solver used to compute the modeled root-mean-square velocity maps.","marker":"[66]"},{"why":"The pPXF full-spectrum fitting method used to extract stellar kinematics from the mock cubes.","marker":"[74]"},{"why":"The adaptive-Metropolis MCMC algorithm used to recover black-hole mass and uncertainties from the kinematic data.","marker":"[76]"}],"fun_headline_variants":["ELT mock data weigh 250k-solar-mass IMBH in VCC 1861","Simulated ELT data weigh a 250k-solar-mass IMBH in Virgo","ELT's HARMONI to unmask a 250,000-solar-mass black hole","HARMONI mock cubes reveal Virgo dwarf's hidden 250k-solar black hole","ELT mock observations spot 250k-solar-mass black hole in Virgo dwarf"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The central stellar density profile is extrapolated from 0.11 arcsec Hubble resolution down to 0.01 arcsec ELT pixels using a core-Sersic model with a fixed inner slope gamma=0.1, and if the real nuclear cusp is steeper or shallower the claimed 5% detectability could disappear.","fun_headline_variants_meta":{"raw":{"variants":["ELT mock data weigh 250k-solar-mass IMBH in VCC 1861","Simulated ELT data weigh a 250k-solar-mass IMBH in Virgo","ELT's HARMONI to unmask a 250,000-solar-mass black hole","HARMONI mock cubes reveal Virgo dwarf's hidden 250k-solar black hole","ELT mock observations spot 250k-solar-mass black hole in Virgo dwarf"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001456,"raw_usage":{"total_tokens":5750,"prompt_tokens":850,"completion_tokens":4900,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":4779}},"tokens_in":594,"tokens_out":4900,"duration_ms":31030,"temperature":1.0,"reasoning_tokens":4779,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:56:23.136959+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-running the same mock pipeline on VCC 1861 with a different inner power-law index (for example gamma=0.5 instead of 0.1) would test whether the central velocity-dispersion peak of the 2.5e5 solar-mass IMBH remains distinguishable from the no-black-hole case; if the profiles become indistinguishable, the central claim fails. Equivalently, a real 4.5-hour ELT/HARMONI observation of VCC 1861 showing no central sigma peak would contradict the predicted detectability.","supporting_citations":[{"cited_title":"Galaxy Populations in the Fornax and Virgo Clusters,","cited_arxiv_id":null,"evidence_quote":"Source of the HST ACS/WFC Virgo survey image and nuclei photometry used for the surface brightness and mass model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The axisymmetric Jeans Anisotropic Modeling formalism used to construct the kinematic maps."},{"cited_title":"The E-ELT first light spectrograph HARMONI: capabilities and modes,","cited_arxiv_id":null,"evidence_quote":"The pPXF full-spectrum fitting method used to extract stellar kinematics from the mock cubes."},{"cited_title":"Simulated stellar kinematics studies of high-redshift galaxies with the HARMONI Integral Field Spectrograph,","cited_arxiv_id":null,"evidence_quote":"The adaptive-Metropolis MCMC algorithm used to recover black-hole mass and uncertainties from the kinematic data."}],"review_version":1}