{"id":"9862e480-6a91-41c2-9c0f-3994041a32c0","arxiv_id":"2501.01493","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"H15A hosts a (2.16+0.23-0.18) x 10^10 solar mass black hole, making it one of the two most massive SMBHs known in the local universe, based on triaxial orbit modeling of Keck KCWI stellar kinematics.","lead":"Using Keck KCWI spectra, the authors measure a supermassive black hole of about 22 billion solar masses in the galaxy Holmberg 15A, along with the galaxy's three-dimensional shape. The result ties the largest black holes known locally and revises an earlier, larger estimate.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Kinematic discrepancy with MUSE is unresolved; the central sigma decline from 340 to 280–300 km/s within 5 arcsec is the main lever on MBH and its origin is not independently established.","rationale":"The reader's weakest_assumption identifies the KCWI/MUSE kinematic discrepancy as the main vulnerability. My independent reading of the paper confirms this: the central σ decline is the kinematic feature that most directly constrains MBH, and while the authors perform extensive internal robustness tests, none of them provides an external anchor for the KCWI kinematics. The discrepancy with MUSE is documented and discussed but not resolved. I agree with the CONDITIONAL verdict: the measurement is plausible and internally consistent, but the unresolved kinematic conflict means the central claim cannot be fully accepted yet. I do not see a more load-bearing concern than this. Other issues (e.g., the effective DOF estimate, MGE modifications) are secondary and would affect the error bars or model comparison rather than the central MBH value as directly as the kinematics.","tokens_in":22747,"tokens_out":1836,"duration_ms":15426,"concrete_test":"Re-extract the stellar kinematics from the publicly available raw KCWI data cubes using an independent pipeline (e.g., a different sky-subtraction and spectral-fitting code) and compare the resulting σ profile to the paper's Figure 7. If the independent extraction reproduces the decline from ~340 to ~280 km/s within R<5 arcsec, the concern is resolved; if it instead produces a flat profile closer to the MUSE values, the central MBH inference is not robust. Alternatively, run the paper's TriOS models with the MUSE kinematics (or with a flat σ profile inserted in the central bins) and check whether the best-fit MBH moves outside the quoted 68% interval.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim, MBH = (2.16+0.23−0.18) × 10^10 M_sun, rests on the KCWI stellar kinematics, especially the central velocity dispersion profile. The paper's own Section 5.1 shows that MUSE-based kinematics (Mehrgan et al. 2019) find a flat σ ~335 km/s over the same radii, whereas KCWI finds σ declining from ~340 to 280–300 km/s by R=5 arcsec. The MBH constraint comes predominantly from the central LOSVD moments, so if the true σ profile were flat as MUSE reports, the inferred MBH would shift substantially. Sections 5.1.1–5.1.4 argue that template choice, GH truncation order, spectral coverage, and parametric versus non-parametric extraction do not explain the difference. However, these tests are all internal to the KCWI data; they do not independently validate the KCWI measurement against an external reference, nor do they reproduce the MUSE measurement on MUSE data. The tests in 5.1.3 are illustrative but cannot settle which instrument/analysis yields the correct kinematics. A second, independent extraction of the KCWI data (e.g., by a different group using a different pipeline and template library, or a reanalysis of the raw MUSE data with the KCWI pipeline choices) is needed. Until that is done, the possibility that the central σ decline is an artifact of KCWI reduction or fitting remains open, and the quoted statistical uncertainties do not cover this systematic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents new Keck KCWI integral-field spectroscopy of the brightest cluster galaxy Holmberg 15A (H15A) covering a roughly 100″ × 100″ contiguous field, binned into 313 spectra from which the first eight Gauss-Hermite moments of the stellar line-of-sight velocity distributions are measured, yielding 2504 kinematic constraints. The authors model these data with the triaxial Schwarzschild code TriOS, running about 40,000 galaxy models over the six-dimensional parameter space of black hole mass MBH, stellar mass-to-light ratio, dark matter mass within 50 kpc, and three intrinsic shape parameters. They infer MBH = (2.16+0.23−0.18) × 10^10 M⊙, luminosity-weighted intrinsic axis ratios p = 0.89 and q = 0.645, a total stellar mass of 2.9 × 10^12 M⊙, and a kinematic axis misaligned by ~62° from the photometric major axis. An axisymmetric reanalysis of the same data gives a worse fit (Δχ2 ≈ 330) and MBH = (2.55 ± 0.20) × 10^10 M⊙. The paper concludes that H15A, along with NGC 4889, hosts one of the two most massive black holes known in the local universe and is a strong outlier relative to the MBH–σ relation, while remaining about a factor of two below the previously reported MUSE-based value of (4.0 ± 0.8) × 10^10 M⊙.","tokens_in":23195,"tokens_out":18243,"duration_ms":164327,"significance":"If the central measurement is correct, this is an important result: H15A would be the first galaxy with a dynamically measured black hole above 2 × 10^10 M⊙ whose mass is determined with a fully triaxial stellar potential, and the mass adds a valuable data point at the extreme high-mass end of the MBH scaling relations. The paper's strengths include the unusually large kinematic dataset (313 bins, eight LOSVD moments), the large orbit-model search (~40,000 models), the split-half consistency test, and the transparency of the comparison with the previous MUSE study. It is a further strength that the orbit-modeling pipeline was validated on simulated galaxies in Pilawa et al. (2024), and that the black-hole mass is derived from a forward model rather than from any scaling relation, so there is no circularity. The paper also honestly quantifies that the assumed symmetry (axisymmetric versus triaxial) changes MBH by only ~18%, isolating stellar kinematics as the dominant source of the disagreement with Mehrgan et al. (2019).","major_comments":[{"comment":"The central kinematic discrepancy between KCWI and MUSE is the main lever on MBH and is not independently resolved. MUSE finds a flat σ ≈ 335 km/s over the central arcseconds while KCWI declines from ≈340 km/s to 280–300 km/s at R = 5″ (Figure 7), and the MBH constraint comes predominantly from the central LOSVD moments. The tests in Sections 5.1.1–5.1.4 show that the KCWI kinematics are stable to template choice, GH truncation order, spectral coverage, and parametric versus non-parametric extraction, but these tests are all internal to the KCWI data; the red-only experiment in Section 5.1.3 does not reproduce the MUSE measurement from MUSE data. The paper itself concludes that \"further tests on MUSE data would be useful.\" As written, the quoted MBH = (2.16+0.23−0.18) × 10^10 M⊙ contains no systematic term for the alternative, flat central σ profile, so the headline number rests on an unresolved instrument-level discrepancy. I recommend reanalyzing the archival MUSE data with the KCWI-style pipeline and template library, or explicitly adding a systematic error to MBH calibrated to the MUSE profile, or demonstrating quantitatively that MBH is insensitive to the central σ shape.","section":"Section 5.1 and Figure 7"},{"comment":"The paper's handling of the low reduced χ2 is qualitative, and the arithmetic does not close the issue. The best-fit χ2 = 1410 with 2504 constraints gives a naive reduced χ2 of 0.57. The paper argues, by analogy with NGC 2693 (Pilawa et al. 2022, 2024), that the effective number of model parameters is about 200, which would raise the reduced χ2 to at most ~0.8 (1410/(2504−200)) — still below unity. A reduced χ2 below 1 indicates that the measurement errors are overestimated, the model is overfitting, or the effective DOF are larger still; in any case the quoted 68% intervals on MBH inherit whatever the resolution is. The split-half test reassures the central value but does not calibrate the uncertainties. The authors should estimate the effective DOF for this specific model and report the resulting reduced χ2, or discuss explicitly how the low χ2 affects the meaning of the quoted intervals.","section":"Section 4.2"},{"comment":"The MGE light model is fit with the imposed constraint σ′ > 0.96″, and the innermost Gaussian component sits exactly at this lower bound (Table 2). Because the deprojected central luminosity density scales as ν0 ∝ Σ0/σ′, this constraint sets a limit on the central stellar density, and the central stellar density is precisely the mass component that competes with MBH in reproducing the central LOSVD. The paper notes that unconstrained fits produce \"exceptionally large and unconstrained central densities,\" but it does not test how MBH responds to the adopted 0.96″ floor. I request a sensitivity test of MBH to the inner-width bound (e.g., floors of 0.75″ and 1.2″), or an explicit demonstration that the central MGE component is not the driver of the MBH inference.","section":"Appendix C and Table 2"},{"comment":"The quoted precision on the intrinsic axis ratios is dominated by the boundaries of the allowed parameter region rather than by the data. The best-fit values q = 0.645+0.001−0.002 (99.7% interval) and u > 0.999 sit at the maximal-flattening boundary q = u q′ with u = 1 of the deprojection inequalities. The luminosity-weighted q is therefore essentially pinned by the projected axis ratios, and the sub-percent formal uncertainty is not a meaningful measurement of the intrinsic shape. Since the triaxial shape is a headline result, the text should note explicitly that q and u are prior-limited at the boundary, and that the robust inference is the coarse statement of strong triaxiality (T ≈ 0.35) rather than the precise q value. The authors' own test showing only a 10% change in MBH when the viewing angles are moved >15° away indicates that this does not affect the black hole mass, but the boundary issue is not discussed in Section 4.2.","section":"Section 4.2, Table 1, Figure 5"}],"minor_comments":[{"comment":"There are several typographical errors: “T able 1” in Table 1; “Give that” should be “Given that” in Section 4.2; “It is is also” in Section 4.2; “T riaxial” and “T riOS” in the title and affiliation line; “W alsh” in the author list; and “T ypeset” in the draft header.","section":"Throughout"},{"comment":"The caption's final sentence is grammatically incomplete: “then performing dynamic nested sampling that surrogate function” should read “then performing dynamic nested sampling of that surrogate function.”","section":"Figure 10 caption"},{"comment":"The sentence “For these, we find that even at the most negative location in the LOSVD, the mean ratio of the amplitude and the uncertainties at that velocity is only −0.88” should be rephrased to state clearly that this is the mean over the affected bins of (LOSVD value)/(1σ error) at the velocity of the most negative value; as written it is easy to misread.","section":"Section 5.1.4"},{"comment":"The claim that H15A and NGC 4889 are “the galaxies hosting the most massive SMBHs known in the local universe” should be explicitly conditioned on the KCWI kinematics, given that the MUSE-based value of (4.0 ± 0.8) × 10^10 M⊙ from Mehrgan et al. (2019) remains unreconciled; a phrase such as “based on the kinematics adopted here” would prevent overstatement.","section":"Abstract and Section 4.3"},{"comment":"The galaxy name is written inconsistently as “H15” in several places in Section 5.1 (e.g., “the central 60″ × 60″ region of H15”) and as “H15A” elsewhere; the paper should use “H15A” consistently.","section":"Section 5.1 and elsewhere"},{"comment":"The caption states that the data have been point-symmetrized for comparison with the model, but does not specify the symmetry convention for each moment (e.g., V and h3 are antisymmetric while σ and h4 are symmetric under point reflection); specifying this would aid the reader.","section":"Figure 9 caption"}],"recommendation":"major_revision","confidential_remarks":"I concur with the stress-test assessment that the unresolved MUSE–KCWI discrepancy in the central velocity dispersion is the decisive issue for this paper. The paper is otherwise unusually thorough, the modeling pipeline is state-of-the-art and externally validated on simulations, and there is no circularity in the mass derivation. The requested revision is substantial but well scoped: either an external cross-check of the central kinematics (e.g., re-extraction of the archival MUSE data with the KCWI-style templates and fitting choices) or an explicit systematic term in the MBH error budget. I do not see grounds for rejection, nor for acceptance as is. The manuscript fits the journal's scope, and the comparison with the prior MUSE study is conducted in good faith; there are no citation or novelty concerns that would affect the verdict."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. This is a careful, technically serious measurement of the SMBH mass in H15A, and it is probably closer to the truth than the prior MUSE-based value. It is also a paper whose headline mass could shift by more than the quoted error bars until the kinematic disagreement with MUSE is resolved by someone outside this group.\n\nWhat is actually new: new deep KCWI spectroscopy over a ~100x100 arcsec field, the first triaxial Schwarzschild modeling for this galaxy, and a mass estimate of MBH = (2.16+0.23-0.18)e10 Msun — roughly half the 4.0e10 Msun from Mehrgan et al. (2019). The galaxy comes out strongly triaxial with luminosity-weighted axis ratios p=0.89 and q=0.65, and axisymmetric fits are worse and give 2.55e10. The TriOS machinery was already tested on simulated galaxies in Pilawa et al. (2024), which is real evidence that the recovery works. Section 5 is unusually candid: the authors lay out the MUSE discrepancy, run honest checks on template choice, GH order, spectral coverage, and parametric vs non-parametric extraction, and state that further tests on MUSE data would be useful. That is how this kind of paper should be written.\n\nSoft spots, in proportion. The main one is exactly what the stress test says. MUSE sees a flat sigma around 335 km/s over the central few arcseconds; KCWI sees it drop from 340 to 280-300 km/s by 5 arcsec. The black hole mass is constrained mainly by those central moments, and all the robustness tests are internal to KCWI — they do not re-reduce MUSE, nor do they validate the KCWI kinematics against a second external instrument. So the discrepancy is unresolved. The quoted uncertainties are statistical only and do not include this systematic. I do not think this sinks the paper; the internal tests are plausible and the triaxial shape is consistent with the observed kinematic misalignment. But 2.16e10 should be treated as provisional until an independent extraction of either dataset appears.\n\nMinor: reduced chi2 of 0.57 with hand-estimated effective degrees of freedom is a real weakness, and the splitting test does not replace a proper covariance treatment. The MGE modifications in Appendix C are post hoc, though described transparently.\n\nBottom line: this paper is for people working on the most massive black holes, IFU kinematics, and triaxial orbit modeling. It deserves a serious referee. I would accept it for review and ask the referee to concentrate on the kinematic systematic rather than the orbital machinery. I would cite it as the new dynamical measurement for H15A while flagging the MUSE discrepancy, and I would bring it to reading group as a clean example of how IFU comparison studies should be presented.","headline":"Careful triaxial KCWI measurement gives H15A a 2.16e10 Msun black hole, probably closer to right than the prior MUSE value, but the unresolved kinematic discrepancy means the mass is provisional.","tokens_in":23668,"tokens_out":3225,"would_cite":true,"duration_ms":32287,"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":"Holmberg 15A hosts one of the two most massive black holes known in the local universe, weighing 22 billion suns.","keywords":["supermassive black holes","galaxy dynamics","triaxial Schwarzschild modeling","Holmberg 15A","stellar kinematics","integral field spectroscopy","scaling relations"],"falsifier":"A re-observation of the inner few arcseconds of H15A with an independent high-resolution spectrograph of comparable or better signal-to-noise, analyzed with both parametric and non-parametric LOSVD fitting, would settle whether the central dispersion drop is real; a flat profile near 335 km/s would falsify the reported black hole mass.","tokens_in":22496,"feed_emoji":"🕳️","tokens_out":6791,"duration_ms":58738,"temperature":0.7,"pith_summary":"Holmberg 15A, the brightest galaxy of the Abell 85 cluster, has an exceptionally dim central core that has made its black hole hard to weigh. This paper reports the most detailed dynamical measurement to date, using Keck KCWI integral-field spectra of 313 spatial bins and triaxial Schwarzschild orbit modeling across roughly 40,000 galaxy models. The authors find a black hole mass of (2.16+0.23-0.18) x $10^{10}$ solar masses, tying NGC 4889 for the most massive black hole known in the local universe. The result matters because it tests how black hole mass scales with galaxy velocity dispersion and stellar mass at the extreme high end, where the scatter is large. It also shows that assuming the galaxy is axisymmetric, rather than triaxial, inflates the inferred black hole mass by about 18 percent.","feed_headline":"22-billion-solar-mass black hole weighed in Holmberg 15A","feed_subtitle":"Triaxial orbit modeling cuts the galaxy's black hole mass nearly in half relative to earlier estimates.","key_machinery":"The central tool is TriOS, a triaxial Schwarzschild orbit-superposition code. Schwarzschild modeling builds a galaxy model from a library of stellar orbits in a given gravitational potential, then assigns non-negative weights to the orbits so their superposition reproduces the observed light distribution and the first eight Gauss-Hermite moments of the stellar velocity distribution at every spatial bin. The paper also uses the observed misalignment between the photometric and kinematic axes as the kinematic signature that forces a triaxial rather than axisymmetric potential, since symmetry in an axisymmetric galaxy requires the two axes to align. The code simultaneously constrains the black hole mass, the stellar mass-to-light ratio, the dark matter halo mass enclosed within 50 kpc, and the three intrinsic shape parameters.","core_discovery":"Using the first eight Gauss-Hermite moments of the stellar line-of-sight velocity distributions from KCWI as constraints, the paper determines the mass and intrinsic shape of H15A simultaneously. The best-fit triaxial model gives MBH = (2.16+0.23-0.18) x $10^{10}$ M_sun, a luminosity-weighted middle-to-long axis ratio p=0.89 and short-to-long ratio q=0.65, a triaxiality parameter T=0.35, and a stellar mass-to-light ratio M*/L_r' = 4.80 solar units. The galaxy's kinematic axis is misaligned from its photometric major axis by about 62 degrees at large radii, which axisymmetric models cannot reproduce. Re-running the analysis with an axisymmetrized version of the orbit code yields a worse fit and raises MBH to (2.55 ± 0.20) x $10^{10}$ M_sun, still far below the (4.0 ± 0.8) x $10^{10}$ M_sun reported from previous MUSE data with axisymmetric modeling. The paper argues through a series of tests that the KCWI kinematics, particularly the central decline in velocity dispersion from ~340 to ~280-300 km/s at 5 arcseconds, are robust to template choice, spectral fitting parameters, spectral coverage, and parametric versus non-parametric LOSVD extraction.","pith_inferences":["If the KCWI kinematics are correct, the earlier MUSE-based mass of 4.0e10 solar masses is likely biased high by a combination of kinematic differences and the axisymmetry assumption; the paper's own axisymmetric fit to KCWI data (2.55e10) isolates the kinematic contribution as roughly the larger part of the gap.","Other ultramassive black hole masses measured with single-aperture or axisymmetric models of cored ellipticals may deserve similar triaxial re-analysis, especially where kinematic misalignment is present or suspected.","A triaxial stellar potential implies the gravitational potential is non-axisymmetric, so the dark matter halo may also be triaxial; modeling the halo as triaxial could shift MBH slightly, though the paper's tests suggest the mass parameters change by only about 10 percent when viewing angles move away from best fit.","The very low central surface brightness of H15A means the black hole's dynamical signature sits on a faint stellar background; if even deeper spectra become available, the central sigma drop could be measured at higher angular resolution and sharpen the MBH constraint further."],"forward_implications":["H15A joins NGC 4889 as the most massive black holes known in the local universe, and both sit well above the mean MBH–sigma relation: H15A is a factor of 10.1 above the McConnell & Ma (2013) relation, a 3-sigma outlier.","The black hole mass is consistent, within 1 sigma, with the MBH–core radius and MBH–bulge mass relations, suggesting that for cored massive ellipticals the stellar core size is a more reliable predictor of MBH than velocity dispersion.","Axisymmetric orbit modeling overestimates MBH by about 18 percent for this galaxy when applied to the same KCWI data, and the axisymmetric models fit significantly worse, with chi-square higher by 330.","The triaxial intrinsic shape with q=0.65 and T=0.35 rules out both oblate and prolate axisymmetry at high confidence, consistent with the observed kinematic misalignment."],"supporting_citations":[{"why":"Supplies the prior MUSE kinematics and the (4.0±0.8)x10^10 M_sun axisymmetric black hole mass that this paper directly compares against and challenges.","marker":"Mehrgan et al. 2019"},{"why":"Provides the TriOS triaxial Schwarzschild code that is the core modeling machinery of the paper.","marker":"Quenneville et al. 2022"},{"why":"Describes the KCWI data reduction and LOSVD measurement procedures that the paper adapts for H15A.","marker":"Liepold et al. 2023"},{"why":"Defines the sphere-of-influence radius and the core-size–black-hole-mass scaling relations used to interpret the result.","marker":"Thomas et al. 2016"},{"why":"Provides the MBH–sigma relation against which H15A is found to be a 3-sigma outlier.","marker":"McConnell & Ma 2013"},{"why":"Supplies the pPXF spectral fitting routine used to extract the first eight Gauss-Hermite moments of the stellar LOSVDs.","marker":"Cappellari 2017"},{"why":"Supplies the Voronoi binning scheme used to co-add spectra into the 313 spatial bins.","marker":"Cappellari & Copin 2003"},{"why":"Supplies the Indo-US stellar template library used as the spectral fitting templates.","marker":"Valdes et al. 2004"}],"fun_headline_variants":["22-billion-solar-mass black hole weighed in triaxial Holmberg 15A","Triaxial orbits pin down 22-billion-solar-mass black hole in Holmberg 15A","Holmberg 15A's black hole weighed at 22 billion suns via triaxial modeling","Triaxial analysis corrects Holmberg 15A black hole to 22 billion solar masses","Black hole in Holmberg 15A weighs 22 billion suns, says triaxial model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central result depends on the KCWI measurement of the central velocity dispersion decline, from about 340 to 280-300 km/s within 5 arcseconds, being unbiased; if the MUSE measurements showing a flat ~335 km/s profile are the true ones, the inferred black hole mass would shift substantially.","fun_headline_variants_meta":{"raw":{"variants":["22-billion-solar-mass black hole weighed in triaxial Holmberg 15A","Triaxial orbits pin down 22-billion-solar-mass black hole in Holmberg 15A","Holmberg 15A's black hole weighed at 22 billion suns via triaxial modeling","Triaxial analysis corrects Holmberg 15A black hole to 22 billion solar masses","Black hole in Holmberg 15A weighs 22 billion suns, says triaxial model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001855,"raw_usage":{"total_tokens":7445,"prompt_tokens":1267,"completion_tokens":6178,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":883,"completion_tokens_details":{"reasoning_tokens":6059}},"tokens_in":883,"tokens_out":6178,"duration_ms":38693,"temperature":1.0,"reasoning_tokens":6059,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:27:10.926632+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A re-observation of the inner few arcseconds of H15A with an independent high-resolution spectrograph of comparable or better signal-to-noise, analyzed with both parametric and non-parametric LOSVD fitting, would settle whether the central dispersion drop is real; a flat profile near 335 km/s would falsify the reported black hole mass.","supporting_citations":[],"review_version":1}