Pith. sign in

REVIEW 4 major objections 4 minor 12 references

The paper argues that ratios between near-infrared hydrogen molecule emission lines can identify actively accreting intermediate-mass black holes in dwarf galaxies, even when optical strong lines and coronal lines are undetectable, for hydr

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 00:32 UTC pith:TED54FOG

load-bearing objection Worth reading: a clean model-based case that H2 line ratios can expose 10^3 Msun IMBHs in low-metallicity dwarfs, but the shock degeneracy means the diagnostic is not yet a classifier. the 4 major comments →

arxiv 2607.14971 v1 pith:TED54FOG submitted 2026-07-16 astro-ph.GA

The Goldilocks Molecule: H₂ Emission Lines Can Identify Elusive Dwarf AGN

Morgan Micharski , Chris Richardson This is my paper
classification astro-ph.GA
keywords dwarf galaxiesintermediate-mass black holesactive galactic nucleiH2 emission linesphotoionization modelsnear-infrared spectroscopyexcitation diagnosticslower mass gap
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors are trying to establish that H2 line ratios can serve as a diagnostic for dwarf active galactic nuclei, specifically active intermediate-mass black holes. They run photoionization models anchored to a low-metallicity dwarf with an observed H2 spectrum and find that for hydrogen densities log nH ≳ 4.0, the ratio between H2 emission lines separates AGN from star-forming regions. This separation holds even for a black hole mass of 10^3 solar masses, where optical strong lines and near-infrared coronal lines would go undetected. If true, this gives observers a new way to find the most elusive intermediate-mass black holes in the lower mass gap.

Core claim

The central claim is that the commonly used H2 2.12 micron/Br gamma ratio misclassifies low-metallicity dwarf AGN as star-forming, but a diagram of ratios between different H2 emission lines places these objects in the AGN zone when the hydrogen density is high enough, log nH ≳ 4.0. In the models, this happens even for a 10^3 solar-mass black hole, where the H2 spectrum would be detectable with a near-infrared spectrograph at the distance of the template dwarf, while [O III] 5007 and coronal lines are too weak for standard optical and near-infrared instruments. The authors conclude that the optical spectrum of such a dwarf could be dominated by starlight while an active intermediate-mass bla

What carries the argument

The central object is the set of H2 emission-line ratios, particularly the excitation diagram that plots one H2 line ratio against another and uses an isothermal gas distribution line to separate AGN from star-forming regions. These ratios carry the argument because the H2 lines respond to the AGN radiation field even in low-metallicity gas where dust-catalyzed H2 formation is inefficient, and unlike the H2/Br gamma ratio they do not get diluted by star formation. The photoionization models include the full set of H2 energy levels, and the geometry fixes cloud mass and column density from the observed radius of the template dwarf J1201, so the predicted line intensities are tied to concrete

Load-bearing premise

The diagnostic's AGN region is also populated by low-density shock models, and the paper does not include shocks in its own grid; the method's specificity rests on the assumption that shocks are absent or separable from AGN excitation.

What would settle it

Take a sample of dwarf galaxies that are optically classified as star-forming but show infrared AGN indicators, measure their H2 line ratios with a near-infrared spectrograph, and check against shock-sensitive diagnostics; if a substantial fraction of objects with clear shock signatures also land in the H2-ratio AGN region, the diagnostic loses its specificity. Alternatively, a single dwarf whose H2 ratios mimic the AGN zone but whose spectra show strong shock emission would directly challenge the clean AGN interpretation.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Near-infrared spectroscopy of nearby dwarf galaxies can search for the H2-ratio AGN signature without needing optical AGN lines, expanding the parameter space for finding intermediate-mass black holes.
  • Many dwarf galaxies previously classified as star-forming on the basis of optical strong lines may harbor hidden active intermediate-mass black holes if their H2 ratios fall in the AGN zone.
  • The diagnostic is density-sensitive: only models with hydrogen density log nH ≳ 4.0 separate cleanly, so observers should target regions of high gas density.
  • The method fails for low-mass clouds (around 10^4 solar masses) because the H2 2.12 micron line is predicted to be too faint to detect, limiting the search to more massive clouds.
  • The technique can reach black hole masses as low as 10^3 solar masses, a regime where optical and coronal diagnostics are blind, offering a direct probe of the lower mass gap.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the shock degeneracy can be resolved—for example, by combining H2 ratios with shock-sensitive line ratios or velocity-resolved observations—the H2-ratio diagram could become a standard tool in dwarf surveys with near-infrared instruments.
  • The same ratio approach might extend to higher-redshift low-metallicity galaxies where optical diagnostics are impractical, potentially revealing a population of intermediate-mass black holes across cosmic time.
  • A stacked analysis of H2 spectra from many optically star-forming dwarfs with infrared AGN indicators could statistically reveal the H2-ratio AGN signature even when individual objects are too faint.
  • The model's implication that the AGN can dominate H2 excitation while starlight dominates optical lines suggests that purely optical searches systematically miss the lowest-mass active black holes, so near-infrared surveys are a necessary complement.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper proposes that near-infrared H2 emission-line ratios can diagnose accreting intermediate-mass black holes (IMBHs) in low-metallicity dwarf galaxies. Using Cloudy 17.02 photoionization models with a spherical-shell geometry anchored to the JWST NIRSpec observations of the dwarf AGN J120122.30+021108.3, the authors compute a grid that varies black hole mass (10^3–10^5 Msun), total cloud mass (10^4–10^5 Msun), hydrogen density (log n_H = 3–5 cm^-3), and metallicity (0.05 and 0.1 Zsun). They find that the standard H2 2.12um/Brγ ratio often misclassifies dwarf AGN as star-forming, whereas a diagram of H2 emission-line ratios separates models with log n_H >= 4 into an AGN region, including models with M_BH = 10^3 Msun for which the optical strong lines and coronal lines would be undetectable. The paper concludes that H2 emission-line ratios are a promising route to identifying otherwise elusive IMBHs.

Significance. If the diagnostic proves robust, it would open a new near-infrared window for finding IMBHs in dwarf galaxies, a population central to black-hole–galaxy co-evolution. The study has clear strengths: it uses a standard, publicly available code (Cloudy); it explores a physically relevant parameter space; and it presents a concrete, falsifiable prediction that the H2-ratio AGN region is occupied only for log n_H >= 4. The use of a real observed template (J1201) is a positive feature. However, the central claim is currently stronger than the evidence: the AGN region is acknowledged to overlap with low-density shock models, the diagram defining the main diagnostic is under-specified, the density grid is described inconsistently, and no uncertainties are propagated. With these gaps, the paper is a useful exploratory study rather than a validated identification method.

major comments (4)
  1. [Sec. 4, Fig. 1d] The paper's own bullet in Sec. 4 states that 'the AGN region is also occupied by low-density shock models (Costa-Souza et al. 2026).' No shock models are included in the present grid, and no quantitative estimate of the contamination is given. Because dwarf galaxies commonly host supernova-driven shocks, a source falling in the AGN region is not uniquely identified as an accreting IMBH. The abstract's claim that H2 ratios 'can reveal active IMBHs' therefore overstates the present evidence. The authors should either add shock models or a literature-based transformation, demonstrate that the high-density H2-ratio regime is not populated by shocks, or explicitly scope the diagnostic to sources for which shocks are independently excluded.
  2. [Sec. 2, Eqs. (1)–(2); Sec. 3] The density-grid description is internally inconsistent. Sec. 2 states: 'For m_tot = 10^4.75 and 10^5.0 Msun, we vary the hydrogen density ... log n_H = 3.0,4.0,5.0 cm^-3, while in least massive m_tot we explore a lower density range log n_H = 3.0,4.0,5.0 cm^-3.' These two sentences are identical. Sec. 3 then says 'the only simulations with m_tot = 10^4.75,10^5.0 Msun not characterized as an AGN have lowest hydrogen densities (n_H = 10^4 cm^-3)', which is not the lowest value in the stated range. Since the central threshold is log n_H >= 4, the exact grid must be specified unambiguously. Please correct the typographical inconsistency and report the full parameter grid.
  3. [Fig. 1d] The text calls Fig. 1d 'particularly promising' without ever defining its axes. It refers to 'H2 emission line ratios' and 'the isothermal gas distribution line' but does not specify which H2 lines are ratioed or how the dividing line is computed. Without these definitions, the diagnostic cannot be reproduced or assessed. Please provide a complete caption and include the relevant formulas in the text.
  4. [Sec. 2, Sec. 3] No uncertainties are propagated, and the sensitivity of the diagnostic to fixed parameters is not explored. The grid fixes the accretion rate at 0.1 mdot_edd and the radiative efficiency at 0.1; H2 excitation can depend on the shape of the ionizing continuum, so varying mdot may change the ratios. Furthermore, the observed J1201 line ratios are not quoted with uncertainties, making the agreement in Fig. 1d only qualitative. Please discuss the dependence on mdot and epsilon_rad and add at least sensitivity limits to the observed point.
minor comments (4)
  1. [Sec. 3] Physically, 'associate detachment' should be 'associative detachment' (the H– + H reaction).
  2. [Sec. 3] The sentence about Figures 1e–g is ambiguous: 'Our models are classified as star-forming except in the [O I]/Halpha plot. The simulations with m_tot = 10^4.75,10^5.0 Msun and M_BH = 10^3 Msun do not appear on these diagrams because their weak [O III] emission...' If they do not appear, how are they classified? Please clarify which models appear in which panel.
  3. [Abstract and throughout] The notation 'logn H' should be 'log n_H' with proper formatting.
  4. [Figure 1] The figure is not shown in the manuscript text; ensure that the file is included and that each panel has a clear caption, especially panel (d) where the axes are missing.

Circularity Check

0 steps flagged

No circularity: H2-line-ratio predictions are forward-model outputs from Cloudy, not fits to J1201; admitted shock degeneracy is a validity caveat, not a circular reduction.

full rationale

The central claim is generated by running Cloudy photoionization models with M_BH, m_tot, n_H, and Z as explicit inputs, then comparing predicted H2 line ratios against empirical AGN regions defined externally (Riffel et al. 2013; Costa-Souza et al. 2026). J1201 is used to set fiducial geometric constraints (Eqs. 1-2) and as a template, not to fit the line-ratio diagnostic; no fitted parameter is renamed as a prediction. The AGN classifications in Fig. 1d are not self-definitional because the 'isothermal gas distribution line' criterion comes from independent observational literature, not from the models themselves. The self-citation to Richardson et al. (2025) concerns photoionization methodology only and is modified here, so it is not load-bearing for the claimed result; citing Polimera et al. (2022) for standard optical diagnostics is likewise not circular. Section 4 explicitly admits 'the AGN region is also occupied by low-density shock models (J. H. Costa-Souza et al. 2026)'; this is a real specificity/false-positive limitation for the proposed diagnostic, but it does not reduce the model output to the model input. Therefore no circular step can be exhibited and the paper earns a 0.

Axiom & Free-Parameter Ledger

7 free parameters · 5 axioms · 0 invented entities

The central prediction is a threshold effect (nH ≥ 10^4 cm^-3) within a multi-dimensional model grid. All varied quantities are inputs, none fitted to the H2 data, so the ledger is transparent but the prediction inherits uncertainty from the assumptions above.

free parameters (7)
  • cloud total mass m_tot = 10^4, 10^4.75, 10^5 Msun (grid)
    Sets whether H2 emission is detectable; lower limit justified from BH-accreted mass (Beckmann et al. 2023), not fitted to H2 data.
  • hydrogen density n_H at illuminated face = 10^3–10^5 cm^-3 (grid)
    Central claim requires log nH ≥ 4.0; varied by hand.
  • black hole mass M_BH = 10^3, 10^4, 10^5 Msun (grid)
    Sets L_bol with fixed mdot; 10^3 case anchors headline prediction.
  • metallicity Z/Zsun = 0.05, 0.1
    Chosen to represent low-metallicity dwarfs; affects gas cooling and H2 formation.
  • accretion rate mdot = 0.1 mdot_Edd
    Fixed assumption; with epsilon=0.1 sets bolometric luminosity.
  • radiative efficiency epsilon_rad = 0.1
    Fixed assumption converting accretion rate to L_bol.
  • helium mass fraction Y = 0.25
    Fixed standard value in Cloudy.
axioms (5)
  • domain assumption Spherical shell geometry in Eqs. (1)–(2) relates observed angular size to physical radius and column density.
    Cloud geometry is assumed, not observed; total radius R≈6.6 pc is derived from small-angle approximation.
  • domain assumption Small-angle approximation using J1201 redshift and NIRSpec 0.1 arcsec aperture gives the physical cloud radius.
    Standard geometric approximation applied to a specific observational setup; enters the geometry and column density.
  • domain assumption The emitting cloud mass m_tot is approximately the mass accreted by the BH over cosmic time, setting the 10^4 Msun lower limit.
    Imported from Beckmann et al. (2023); determines which cloud masses are considered.
  • domain assumption Cloudy 17.02 with full H2 and FeII level sets adequately treats H2 excitation and formation at low metallicity.
    The central predictions are generated by this code; no independent verification is provided.
  • domain assumption The H2 ratio AGN/star-forming classification boundaries from Riffel et al. (2013) and Costa-Souza et al. (2026) transfer to low-metallicity dwarf AGN.
    The diagnostic diagram is calibrated on massive AGN and shock models, not on dwarf IMBHs; the paper applies it directly.

pith-pipeline@v1.3.0-alltime-deepseek · 3282 in / 10872 out tokens · 110151 ms · 2026-08-02T00:32:21.271640+00:00 · methodology

0 comments
read the original abstract

We propose using H$_{2}$ emission lines as a novel diagnostic to identify dwarf AGN by running photoionization models incorporating active intermediate-mass black holes and using an observed WISE dwarf AGN as a template. Though many dwarf AGN would be incorrectly classified as star-forming if the H$_{2}$~2.12~$\mu$m/Br$\gamma$ ratio is used, ratios between H$_{2}$ emission lines can reveal active IMBHs in cases where log $n_{\text{H}} \gtrsim 4.0$. This includes the case where $M_{\text{BH}} = 10^3~M_{\odot}$, even though the optical strong lines and NIR coronal lines would likely go undetected. We conclude that H$_{2}$ emission lines show promise in detecting the most elusive IMBHs.

Figures

Figures reproduced from arXiv: 2607.14971 by Chris Richardson, Morgan Micharski.

Figure 1
Figure 1. Figure 1: An illustration of the geometric setup of our simulations (a), followed by the detectability of the H2 2.12 µm line with JWST NIRSpec at z = 0.0034 (b), and several excitation diagnostics in the NIR and optical (c-g). The diagram in (d) is particularly promising for finding elusive IMBHs. calculate the interior radii (rin) and column densities (NH) for each input value, with NH as the boundary condition fo… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

12 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    Nature Astronomy , keywords =

    Variability as a new discovery channel for intermediate-mass black holes in the time-domain era. Nature Astronomy , keywords =. doi:10.1038/s41550-025-02759-5 , archivePrefix =. 2601.13311 , primaryClass =

    arXiv

  2. [2]

    Near-infrared diagnostic diagrams of the gas ionization sources in nearby galaxies: a JWST NIRSpec view

    Near-infrared diagnostic diagrams of the gas ionization sources in nearby galaxies: a JWST NIRSpec view. arXiv e-prints , keywords =. doi:10.48550/arXiv.2605.04925 , archivePrefix =. 2605.04925 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2605.04925

  3. [3]

    , keywords =

    Population statistics of intermediate-mass black holes in dwarf galaxies using the NEWHORIZON simulation. , keywords =. doi:10.1093/mnras/stad1544 , archivePrefix =. 2211.13301 , primaryClass =

    Pith/arXiv arXiv

  4. [4]

    , keywords =

    Intermediate-Mass Black Holes. , keywords =. doi:10.1146/annurev-astro-032620-021835 , archivePrefix =. 1911.09678 , primaryClass =

    Pith/arXiv arXiv 1911

  5. [5]

    , keywords =

    Local Analogs of Primordial Galaxies: In Search of Intermediate-mass Black Holes with JWST NIRSpec. , keywords =. doi:10.3847/1538-4357/adcd59 , archivePrefix =. 2408.04774 , primaryClass =

    Pith/arXiv arXiv

  6. [6]

    Low-ionization nuclear emission-line region and star-forming galaxies

    Molecular hydrogen and [Fe II] in active galactic nuclei - III. Low-ionization nuclear emission-line region and star-forming galaxies. , keywords =. doi:10.1093/mnras/stt026 , archivePrefix =. 1301.0491 , primaryClass =

    Pith/arXiv arXiv

  7. [7]

    , keywords =

    Molecular hydrogen and [Fe II] in Active Galactic Nuclei. , keywords =. doi:10.1051/0004-6361:20034285 , archivePrefix =. astro-ph/0406402 , primaryClass =

    Pith/arXiv arXiv

  8. [8]

    , keywords =

    An Upper Limit to the Degree of Evolution between Supermassive Black Holes and Their Host Galaxies. , keywords =. doi:10.1086/508055 , archivePrefix =. astro-ph/0608091 , primaryClass =

    Pith/arXiv arXiv

  9. [9]

    , keywords =

    The 2017 Release Cloudy. , keywords =. doi:10.48550/arXiv.1705.10877 , archivePrefix =. 1705.10877 , primaryClass =

    Pith/arXiv arXiv doi:10.48550/arxiv.1705.10877 2017

  10. [10]

    BASS. L. Near-infrared Data Release 3: A Spectral Atlas and Characterization of Active Galactic Nuclei. , keywords =. doi:10.3847/1538-4365/ae2d58 , archivePrefix =. 2512.20702 , primaryClass =

    arXiv

  11. [11]

    , keywords =

    RESOLVE and ECO: Finding Low-metallicity z 0 Dwarf AGN Candidates Using Optimized Emission-line Diagnostics. , keywords =. doi:10.3847/1538-4357/ac6595 , archivePrefix =. 2204.03633 , primaryClass =

    Pith/arXiv arXiv

  12. [12]

    , keywords =

    Emission-line Diagnostics for IMBHs in Dwarf Galaxies: Accounting for Black Hole Seeding and Ultraluminous X-Ray Source Excitation. , keywords =. doi:10.3847/1538-4357/adf215 , archivePrefix =. 2505.07749 , primaryClass =

    Pith/arXiv arXiv