{"id":"746db087-9f13-4b60-b4b0-a2bf607365af","arxiv_id":"1909.02585","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Gallo and Sesana estimate the local black hole mass function between 10^5 and 10^6 solar masses by folding GAMA galaxy counts together with X-ray occupation fraction constraints, finding a normalization uncertainty of about one dex.","lead":"This paper combines the GAMA galaxy census with X-ray measurements of the fraction of nearby galaxies that host a black hole to estimate how many black holes with masses between 100,000 and one million suns exist locally. The estimate matters for LISA, whose gravitational-wave sensitivity peaks in exactly this black hole mass range.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The <1 dex claim is conditioned on an occupation fraction measured only in early-type galaxies, and the paper's own Section 4 calls this a serious drawback; late types dominate the relevant stellar mass range, so the uncertainty band excludes the largest known systematic.","rationale":"The central numerical claim is not that a measured BHMF exists, but that the modeled range has <1 dex normalization uncertainty in [10^5,10^6] M⊙. For that claim to be true, the uncertainty budget in Eq. 3 must include all systematics that can shift the normalization by an order of magnitude. The largest such systematic explicitly acknowledged in the paper is the morphology dependence of the occupation fraction: Section 4 states that using early-type galaxies 'represents a serious drawback for our estimates, especially considering the increased prominence of late types at low stellar masses.' The M⋆ range that maps to MBH≈10^5–10^6 M⊙ under the adopted combined relation is roughly 10^9.2–10^9.8 M⊙, where late types dominate the GAMA mass function. At those masses the early-type λ_occ is near unity (~0.9–0.99), so even a modest factor-of-3 reduction in late-type occupation moves the median Φ by ~0.5 dex; a factor of 10 reduction would exceed the quoted 1 dex band. The paper's Figure 3 right panel shows how a future 5,000-galaxy late-type sample would improve λ_occ, but that is a projection, not a correction included in the current uncertainty. The scaling-relation extrapolation is a second important caveat, but the paper at least brackets it with three variants of the Reines & Volonteri relation; the morphology-dependent occupation is not bracketed at all. Thus the most load-bearing concern is exactly the reader's weakest assumption, and the recommended response is to keep the CONDITIONAL verdict until a morphology-resolved occupation measurement or a bracketing calculation tests the shift.","tokens_in":8419,"tokens_out":10043,"duration_ms":114916,"concrete_test":"Recompute Eq. 3 using a morphology-resolved version of the GAMA stellar mass function: split Φ(M⋆) into early- and late-type components (Wright et al. 2017), keep the early-type λ_occ from Section 2, and assign late types an independently measured occupation function (e.g., from X-ray selected AGN/dwarf nuclei at the same M⋆); if no such sample exists, adopt bracketing extremes λ_occ,late=1.0 and 0.1. Compare the median and 68% band at MBH=3×10^5 M⊙ with Figure 2. If the median shifts by more than ~0.5 dex, or the band widens beyond 1 dex, the sub-dex claim fails and the conclusion must be reworded to state the conditioning on early-type occupation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline normalization claim (abstract; Section 4) is that the modeled BHMF has ≲1 dex 68% uncertainty in [10^5,10^6] M⊙. This number comes from Eq. 3, which convolves the GAMA stellar mass function with λ_occ(M⋆) from Eq. 1. The λ_occ posterior is measured from 326 early-type galaxies only (Section 2). For the adopted combined MBH–M⋆ relation (α=8.13, β=1.72), the M⋆ range that maps to MBH∈[10^5,10^6] M⊙ is roughly 10^9.2–10^9.8 M⊙. In that range late-type galaxies dominate the GAMA mass function, so the central claim silently assumes the early-type occupation fraction applies to them. The paper explicitly labels this 'a serious drawback for our estimates' in Section 4 but does not propagate any morphology-dependent systematic into the yellow 68% band of Figure 2; Figure 3's right panel only simulates future late-type constraints. Because λ_occ(early) is already ~0.9–0.99 at these M⋆, a factor-of-3 lower late-type occupation shifts the median Φ by ~0.5 dex, and a factor-of-10 lower occupation moves it outside the quoted band. The sub-dex uncertainty is therefore an uncertainty over fitted parameters conditioned on the early-type sample, not an uncertainty on the true local BHMF.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter estimates the local black hole mass function (BHMF) below 10^6 solar masses by convolving the GAMA galaxy stellar mass function with a black hole occupation fraction derived from Chandra X-ray observations of 326 nearby early-type galaxies and with a log MBH–log Mstar scaling relation fit to the combined AGN plus inactive galaxy sample of Reines & Volonteri (2015). The central quantitative claim is that the 68% normalization uncertainty of the resulting BHMF is less than 1 dex over the range [10^5, 10^6] M_sun, the range most relevant for LISA. The paper also compares with AGN-only and inactive-only scaling relations, provides a Schechter-like fitting formula for the median BHMF, and illustrates how future X-ray surveys with about 5,000 galaxies could shrink the occupation-fraction uncertainty.","tokens_in":8713,"tokens_out":6888,"duration_ms":72966,"significance":"If the result holds, this is one of the first observationally grounded estimates of the BHMF below 10^6 M_sun, with direct implications for LISA event-rate predictions and for seeding models of supermassive black holes. The convolution framework in Equation (3) is transparent, and the use of an empirically constrained occupation fraction rather than an assumed 100% occupation is a genuine step forward. The paper also deserves credit for propagating the posterior uncertainties of the fitted parameters and for providing a simple analytic approximation to the median BHMF. However, the headline uncertainty claim is conditioned on an assumption that the authors themselves describe as a serious drawback, namely that the early-type occupation fraction applies to the full galaxy population used in the convolution; this systematic is not included in the reported uncertainty band.","major_comments":[{"comment":"In Section 4 the authors write that using the early-type occupation fraction 'represents a serious drawback for our estimates,' and yet the 68% band shown in Figure 2 and the less-than-1-dex claim in the abstract are computed from Equation (3) using only the lambda_occ posterior derived from 326 early-type galaxies. In the stellar mass range that maps to MBH in [10^5, 10^6] M_sun for the adopted scaling relation (roughly Mstar from 10^9.2 to 10^9.8 M_sun), late-type galaxies dominate the GAMA mass function, so a systematically lower late-type occupation fraction would shift the median Phi by about 0.5 dex for a factor-of-3 difference and by more than the quoted band for a factor-of-10 difference. Because no morphology-dependent systematic is included in the yellow band, the headline normalization uncertainty is conditional on the early-type occupation fraction rather than being an uncertainty on the true local BHMF.","section":"Section 4, Figure 2"},{"comment":"Equation (3) adopts a single linear log MBH–log Mstar relation fitted to the Reines & Volonteri sample and then applies it down to the GAMA mass limit, Mstar approximately 10^7.5 M_sun, with no allowance for a change of slope, intercept, or scatter in the dwarf-galaxy regime. The posterior range of alpha, beta, and sigma is propagated, but the unvalidated extrapolation itself is not included in the systematic budget of Figure 2. Since the central result is the normalization in the 10^5–10^6 M_sun range, where this extrapolation is directly relevant, the paper should either justify the extrapolation with low-mass black hole mass estimates and upper limits or add an explicit systematic term.","section":"Section 3, Equation (3)"}],"minor_comments":[{"comment":"Equation (3) writes MBH in the exponent, but the scaling relation is defined for log MBH; please make the logarithm explicit (for example, log MBH - (alpha + beta log Mstar)) to avoid ambiguity.","section":"Section 3, Equation (3)"},{"comment":"The typeset argument of the tanh in Equation (1) is missing a division: it should be 2.5 |8.9 - log Mstar,0| log(Mstar/Mstar,0), not 2.5 |8.9 - log Mstar,0| log Mstar Mstar,0.","section":"Section 2, Equation (1)"},{"comment":"The caption of Figure 2 contains a typo ('Gair et at. 2010' should be 'Gair et al. 2010') and also uses 'BHFM' instead of 'BHMF'; both should be corrected.","section":"Figure 2 caption"},{"comment":"In the discussion of the right panel of Figure 3, the phrase 'the corresponding BHMF uncertainty would shrink to less than about 2' lacks units; please specify whether this is 2 dex or another quantity.","section":"Section 4, right panel of Figure 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript makes a useful contribution, but the abstract and conclusion overstate the robustness of the sub-dex uncertainty. The authors themselves identify the early-type-only occupation fraction as a serious drawback in Section 4, which supports a major revision. I would be satisfied if the headline claim is explicitly conditioned on the early-type occupation fraction and if the morphological systematic is either propagated into the uncertainty band or clearly separated from the statistical uncertainty in the presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one new thing here is that they actually fold an observationally measured occupation fraction into a local BHMF calculation and push it down to 1e5–1e6 Msun. That is a first, and it matters for LISA event-rate forecasts and seeding models. The arithmetic is transparent: Eq. 3 is a standard convolution, the uncertainties are propagated from the posterior draws, and the comparison with earlier BHMFs (Shankar, Gair, Barausse) is useful. The paper also gives a simple Schechter-like fit to their median, which is convenient for people who want to plug it into something. Credit where due: the authors do not oversell the method, and they are explicit about the main caveats.\n\nThe soft spots are real but not hidden. The occupation fraction comes from 326 early-type galaxies only, while late types dominate the stellar mass range (Mstar ~ 1e9–1e10) that maps to the black hole masses of interest. The stress-test note is correct: a factor-of-3 lower late-type occupation shifts the median by ~0.5 dex, and a factor-of-10 lower occupation moves it outside the quoted 68% band. So the \"<1 dex normalization uncertainty\" in the abstract is not an uncertainty on the true local BHMF; it is a conditional statement about the fitted parameters given the early-type sample. The paper even calls this \"a serious drawback\" in Section 4, but does not propagate any morphology-dependent systematic into the yellow band. That is the one place where the claim in the abstract goes a bit further than what the error bars actually represent. The second soft spot is the MBH–Mstar relation: they fit at higher masses and extrapolate down, with scatter 0.61 dex, and they use a combined AGN+inactive sample that Reines & Volonteri themselves argued should not be fit with a single relation. The authors know this and justify it as bracketing uncertainty, which is a defensible choice, but it means the result is a plausible first constraint rather than a closed measurement.\n\nBottom line: this is a useful, clearly stated Letter for anyone working on low-mass black hole demographics or LISA extreme-mass-ratio inspiral rates. It deserves a serious referee—the central idea is sound, the weaknesses are acknowledged, and the result is likely to be cited as the first occupation-corrected BHMF below 1e6. I would accept it to refereeing, with the main request being that the morphology caveat be explicitly reflected in the abstract's uncertainty claim. For my own work, I would cite it with the caveat noted.","headline":"A transparent, honest first estimate of the low-mass black hole mass function with a real weakness—the occupation fraction is early-type only—but the authors say so plainly and the sub-dex claim should be read as conditional on that morphology assumption.","tokens_in":9301,"tokens_out":1184,"would_cite":true,"duration_ms":14493,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the local number density of black holes between 100,000 and one million solar masses is pinned to within a factor of ten once galaxy counts and X-ray occupation constraints are combined.","keywords":["black hole mass function","occupation fraction","dwarf galaxies","intermediate-mass black holes","galaxy stellar mass function","X-ray surveys","LISA gravitational wave detector"],"falsifier":"Measure the nuclear black hole occupation fraction in a large, uniformly observed sample of late-type and dwarf galaxies with $M_\\star$ between roughly $10^9$ and $10^{10}\\,M_\\odot$ (for example, a deep X-ray survey of about 5,000 such galaxies with sub-arcsecond resolution). If the occupation fraction at those masses differs substantially from the early-type-derived tanh ramp, the low-mass normalization of the BHMF shifts by the difference and the claimed sub-dex precision in the $10^5$–$10^6$ $M_\\odot$ band would be contradicted.","tokens_in":8097,"feed_emoji":"🕳️","tokens_out":14584,"duration_ms":128949,"temperature":0.7,"pith_summary":"The paper tries to establish that the local number density of black holes with masses between $10^5$ and $10^6$ solar masses can be pinned down to within about a factor of ten, despite the absence of a clean way to weigh such black holes. The method combines the GAMA galaxy stellar mass function with an X-ray-based estimate of the fraction of galaxies that actually host a nuclear black hole, then convolves the two through a broad $M_{\\rm BH}$–$M_\\star$ scaling relation. If the estimate is right, LISA's most sensitive mass band has a known target density, which sharpens predictions for black hole merger rates and gives a local anchor for models of how the first supermassive black holes were seeded.","feed_headline":"Black hole census below a million solar masses within one dex","feed_subtitle":"A new combination of galaxy counts and X-ray occupation data gives LISA its low-mass target.","key_machinery":"The load-bearing element is the convolution in Equation 3, which builds the BHMF as $\\Phi(M_{\\rm BH}) = \\int \\Phi(M_\\star)\\,\\lambda_{\\rm occ}(M_\\star)\\,\\frac{1}{\\sqrt{2\\pi\\sigma^2}}\\exp[-(M_{\\rm BH}-(\\alpha+\\beta M_\\star))^2/(2\\sigma^2)]\\,dM_\\star$. Here $\\Phi(M_\\star)$ is the GAMA galaxy stellar mass function, $\\lambda_{\\rm occ}$ is the occupation fraction of Equation 1 — a smooth tanh ramp that is nearly zero below $10^7\\,M_\\odot$ and nearly one above $10^{10}\\,M_\\odot$, with a transition mass $M_{\\star,0}$ constrained by X-ray observations of 326 nearby early-type galaxies — and the Gaussian carries the scatter in the $M_{\\rm BH}$–$M_\\star$ relation. This object does the work: it turns a galaxy census into a black hole census while suppressing low-mass galaxies that probably lack black holes, which is why the low-mass end of the BHMF is no longer set by the scaling relation alone.","core_discovery":"The paper's central claim is that the low-mass end of the local black hole mass function is empirically accessible despite the lack of a reliable black hole mass proxy. By folding an X-ray-derived occupation fraction into the convolution between galaxy stellar mass and black hole mass, the authors obtain a BHMF whose 68% normalization uncertainty is less than one dex across $10^5$–$10^6$ $M_\\odot$, the band where LISA is most sensitive. The median best-guess mass function is approximated by $\\log\\Phi = -2.13 - 0.098\\log(M_{\\rm BH}/M_\\odot) - 0.00011\\,(M_{\\rm BH}/M_\\odot)$. The paper further reports that different choices of the $M_{\\rm BH}$–$M_\\star$ scaling relation change the inferred BHMF by only a factor of about three below $10^7$ $M_\\odot$, so the occupation fraction, rather than the scaling relation, is the main lever on the low-mass normalization.","pith_inferences":["If late-type and dwarf galaxies turn out to have a lower nuclear occupation fraction at $M_\\star \\sim 10^9$–$10^{10}\\,M_\\odot$ than the early-type-based ramp, the claimed sub-dex normalization in the $10^5$–$10^6\\,M_\\odot$ band is likely optimistic, since a factor-of-two change in $\\lambda_{\\rm occ}$ shifts the BHMF normalization by roughly the same factor.","The same convolution could be run with the redshift-evolving galaxy stellar mass function to produce a $z>0$ BHMF, giving LISA a direct prediction for the redshift distribution of massive black hole mergers; the paper itself only builds the local function.","The paper's median BHMF, inserted into binary population models, yields concrete predictions for local merger and extreme-mass-ratio inspiral rates that LISA can test once it flies.","A cleaner test of the low-mass census would come from dynamical black hole mass measurements in dwarf galaxies, which would replace the broad $M_{\\rm BH}$–$M_\\star$ scatter with actual masses and directly validate the Gaussian convolution."],"forward_implications":["The local volume density of black holes near a few times $10^5\\,M_\\odot$ is bracketed to within about a factor of ten, so rate predictions for low-mass black hole mergers in LISA no longer float freely over this band.","Within $10^5$–$10^6\\,M_\\odot$ the 68% normalization uncertainty stays below one dex for each of the scaling relations considered, so the remaining disagreement among black hole mass proxies matters less than feared for the census.","Below $10^7\\,M_\\odot$ the occupation-corrected BHMFs built from different scaling relations agree to within a factor of about three, making the low-mass end comparatively stable against the proxy choice.","The median BHMF is approximated by $\\log \\Phi = -2.13 - 0.098\\log(M_{\\rm BH}/M_\\odot) - 0.00011\\,(M_{\\rm BH}/M_\\odot)$, a compact form that population-synthesis calculations can adopt directly.","A future X-ray survey of roughly 5,000 nearby galaxies with uniform, sub-arcsecond coverage could reduce occupation-fraction uncertainties to a few percent and tighten the BHMF further across the LISA band."],"supporting_citations":[{"why":"supplies the GAMA galaxy stellar mass function down to about 10^7.5 solar masses that enters Equation 3 as the host galaxy census.","marker":"Wright et al. (2017)"},{"why":"provides the Bayesian formalism and the base sample of 194 early-type galaxies used to derive occupation fraction constraints.","marker":"Miller et al. (2015)"},{"why":"adds 132 early-type galaxies and supplies the updated Mstar,0 posterior plus the simulated target counts for future occupation measurements.","marker":"Gallo et al. (2019)"},{"why":"assembles the AGN and inactive galaxy samples and the black-hole-to-stellar-mass relations whose combined fit seeds the mass function convolution.","marker":"Reines & Volonteri (2015)"},{"why":"provides the LINMIX_ERR regression machinery used to fit the intercept, slope, and scatter of the adopted black-hole-to-stellar-mass relation.","marker":"Kelly (2007)"},{"why":"provides a high-mass local BHMF estimate used as a comparison anchor above 10^6 solar masses.","marker":"Shankar et al. (2009)"},{"why":"provides updated high-mass BHMF constraints with revised scaling relations and selection corrections used for comparison.","marker":"Shankar et al. (2016)"},{"why":"provides the heavy-plus-light seed semi-analytic BHMF that brackets the high-normalization end of LISA rate predictions.","marker":"Barausse (2012)"},{"why":"provides the low-normalization BHMF prescription used as the pessimistic bracket for LISA rate predictions.","marker":"Gair et al. (2010)"}],"fun_headline_variants":["Black hole census below a million suns within one dex","LISA's target black holes mapped to one-dex precision","Local black hole masses pinned down below a million suns","One-dex black hole count for the LISA era","Low-mass black hole census sharpened for LISA"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the occupation fraction measured from 326 early-type galaxies applies to all galaxies at the same stellar mass, even though late-type and dwarf galaxies outnumber early types below $M_\\star \\sim 10^{10}\\,M_\\odot$; the authors themselves flag this as a serious drawback in Section 4.","fun_headline_variants_meta":{"raw":{"variants":["Black hole census below a million suns within one dex","LISA's target black holes mapped to one-dex precision","Local black hole masses pinned down below a million suns","One-dex black hole count for the LISA era","Low-mass black hole census sharpened for LISA"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000641,"raw_usage":{"total_tokens":2919,"prompt_tokens":885,"completion_tokens":2034,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":1953}},"tokens_in":501,"tokens_out":2034,"duration_ms":16107,"temperature":1.0,"reasoning_tokens":1953,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:46:32.209025+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the nuclear black hole occupation fraction in a large, uniformly observed sample of late-type and dwarf galaxies with $M_\\star$ between roughly $10^9$ and $10^{10}\\,M_\\odot$ (for example, a deep X-ray survey of about 5,000 such galaxies with sub-arcsecond resolution). If the occupation fraction at those masses differs substantially from the early-type-derived tanh ramp, the low-mass normalization of the BHMF shifts by the difference and the claimed sub-dex precision in the $10^5$–$10^6$ $M_\\odot$ band would be contradicted.","supporting_citations":[{"cited_title":"P., Gallo, E., Greene, J","cited_arxiv_id":null,"evidence_quote":"provides the Bayesian formalism and the base sample of 194 early-type galaxies used to derive occupation fraction constraints."},{"cited_title":"Towards a high accuracy measurement of the local black hole occupation fraction in low mass galaxies","cited_arxiv_id":"1903.06629","evidence_quote":"adds 132 early-type galaxies and supplies the updated Mstar,0 posterior plus the simulated target counts for future occupation measurements."},{"cited_title":"E., & V olonteri, M","cited_arxiv_id":null,"evidence_quote":"assembles the AGN and inactive galaxy samples and the black-hole-to-stellar-mass relations whose combined fit seeds the mass function convolution."},{"cited_title":"H., & Miralda-Escud´e, J","cited_arxiv_id":null,"evidence_quote":"provides a high-mass local BHMF estimate used as a comparison anchor above 10^6 solar masses."},{"cited_title":"R., Tang, C., & V olonteri, M","cited_arxiv_id":null,"evidence_quote":"provides the low-normalization BHMF prescription used as the pessimistic bracket for LISA rate predictions."}],"review_version":1}