Low-Luminosity AGN: Accretion & Feedback

• Low-luminosity AGN are supermassive black holes with bolometric luminosities ≤10^42–10^43 erg s⁻¹ and low Eddington ratios, powered by radiatively inefficient accretion flows such as ADAF.
• They exhibit unique multiwavelength signatures including hard X-ray spectra, strong compact MIR/NIR emissions, and elevated CO(2–1)/CO(1–0) ratios, setting them apart from classical AGN.
• Their feedback is dominated by kinetic outflows and shocks that modify the circumnuclear ISM, influencing nuclear star formation and host galaxy evolution.

Low-luminosity active galactic nuclei (LLAGN) are accreting supermassive black holes (SMBHs) characterized by bolometric luminosities orders of magnitude below those found in luminous Seyferts, quasars, or classical AGN, typically with $L_\mathrm{bol} \lesssim 10^{42{-}43}$ erg s⁻¹ and Eddington ratios $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$–$10^{-2}$. LLAGN constitute the dominant population of AGN in the local universe and are found in diverse host environments, from late-type disks to massive ellipticals, spanning SMBH masses from $10^4$ to $10^{9}\,M_{\odot}$. Their low luminosities arise from low accretion rates and/or radiatively inefficient accretion modes, with spectral energy distributions, emission line diagnostics, and feedback mechanisms distinct from those in high-luminosity AGN.

1. Defining Characteristics and Accretion Physics

LLAGN are most robustly defined in terms of their Eddington-scaled accretion power and (rest-frame) X-ray and bolometric luminosities. Typical criteria include:

• Bolometric luminosity $L_\mathrm{bol} \lesssim 10^{42{-}43}$ erg s⁻¹ and $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$ erg s⁻¹ (Dumont et al., 2019, Annuar et al., 2020, Asmus et al., 2011).
• Eddington ratio $\lambda_\mathrm{Edd} \sim 10^{-6}$–$10^{-2}$ (Lusso et al., 24 Feb 2025, Zhang et al., 2022).

The physics of accretion in LLAGN is dominated by radiatively inefficient accretion flows (RIAFs), such as advection-dominated accretion flows (ADAF) (Nemmen et al., 2011). In these systems, low-density, hot, optically thin gas flows onto the SMBH, often exhibiting mass-loss via winds, and a significant portion of gravitational energy is advected or converted into kinetic power rather than radiated. The accretion structure is typically a truncated thin disk at radii $R_\mathrm{tr} \sim 30{-}300\,R_S$, with an inner ADAF and an outflowing relativistic jet (Nemmen et al., 2011, Yan et al., 12 Mar 2025).

Radiative efficiency drops to $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$0–$$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$1, and kinetic efficiency in jets $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$2–$$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$3, with jet kinetic power often exceeding bolometric luminosity ($$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$4) (Nemmen et al., 2011). The corona, responsible for X-ray emission in luminous AGN, is generally weak or absent at the lowest $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$5 (Yuk et al., 29 Jan 2025).

2. Multiwavelength Signatures: Spectral Energy Distribution and Emission Components

The SEDs of LLAGN deviate systematically from those of classical AGN (Dumont et al., 2019, Mason et al., 2012, Gilbert et al., 3 Aug 2025, Nemmen et al., 2011):

• The X-ray power-law continuum is often hard (photon index $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$6–2.3), lacking the “big blue bump” of the standard thin disk (Lusso et al., 24 Feb 2025).
• The radio-to-mm regime is dominated by jet synchrotron emission, with parabolic-to-conical jet morphologies seen in VLBI imaging (Yan et al., 12 Mar 2025).
• Strong, compact nuclear mid-infrared (MIR, $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$7m) and near-infrared (NIR, $$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$8m) emission is detected in most LLAGN, correlating tightly with X-ray luminosity ($$\lambda_\mathrm{Edd} \equiv L_\mathrm{bol}/L_\mathrm{Edd} \sim 10^{-6}$$9 with $10^{-2}$0) (Asmus et al., 2011, Mason et al., 2012, Dumont et al., 2019). The MIR–X-ray relation appears independent of $10^{-2}$1 down to $10^{-2}$2 (Asmus et al., 2011).
• At the lowest Eddington ratios ($10^{-2}$3), NIR emission can exceed X-ray output, with $10^{-2}$4 rising toward low $10^{-2}$5 (Dumont et al., 2019).
• Strong silicate emission is present even when the canonical torus is weak or absent, indicative of optically thin, possibly clumpy dust geometries or low dust-to-gas ratios (Mason et al., 2012).
• Star-formation contributions to nuclear NIR/MIR flux are generally minor within tens of parsecs of the SMBH (Asmus et al., 2011, Mason et al., 2012), though some high-resolution SEDs show NIR excesses from young nuclear star clusters.

3. Emission Line Diagnostics and Molecular Gas Content

LLAGN are challenging to identify via standard optical emission-line diagnostics, especially in galaxies with strong star formation. Conventional BPT diagrams ([O III]/Hβ vs. [N II]/Hα) tend to miss LLAGN embedded in rapidly star-forming disks, as H II regions dominate the spectrum (Bär et al., 2016). More sensitive indices such as He II $10^{-2}$6/Hβ are required to uncover AGN photoionization in high-SFR hosts, yielding a $10^{-2}$7 increase in faint AGN identifications (Bär et al., 2016). X-ray confirmation (power-law spectra with $10^{-2}$8–2.1, $10^{-2}$9 erg s⁻¹) substantiates the AGN nature in weak-lined cases, ruling out pure starburst origins.

In molecular gas studies, LLAGN hosts show:

• CO(1–0) and CO(2–1) emission detected in nuclei, though integrated intensities and total gas masses ($10^4$0) are modest (Boeker et al., 2011).
• Elevated CO(2–1)/CO(1–0) line ratios ($10^4$1) in face-on LLAGN, compared to $10^4$2 in quiescent spirals—a signature of AGN-induced molecular gas excitation, not explained by beam filling (Boeker et al., 2011).
• A scarcity of dense molecular gas (n(H$10^4$3) $10^4$4 cm⁻³), as evidenced by upper limits on HCN and HCO$10^4$5 emission, distinguishing LLAGN from powerful AGN with dense star-forming cores (Boeker et al., 2011).

4. Variability, Duty Cycle, and Identification Techniques

LLAGN identification via amplitude statistics of nuclear variability is highly effective. Optical and X-ray variability selection has uncovered populations missed by X-ray, MIR, and line-ratio surveys (Yuk et al., 29 Jan 2025, Ding et al., 2018, Young et al., 2012, Yuk et al., 2022). In high-cadence time-domain optical surveys, variability-selected LLAGN compose $10^4$62% of local bright galaxies for $10^4$7 (Yuk et al., 29 Jan 2025, Yuk et al., 2022), with a spectroscopic AGN confirmation rate of $10^4$8%. These sources often exhibit soft X-ray spectra (mean $10^4$9) and are rarely MIR- or X-ray-selected (Yuk et al., 29 Jan 2025).

X-ray variability metrics (normalized excess variance, power spectral density break scaling as $10^{9}\,M_{\odot}$0) show that LLAGN variability amplitudes plateau at low luminosity, reproducing observed anti-correlation breaks in $10^{9}\,M_{\odot}$1–$10^{9}\,M_{\odot}$2 relations (Ding et al., 2018, Young et al., 2012).

“Changing-look” LLAGN, such as NGC 4614, demonstrate transitions in optical type and X-ray flux coincident with large drops in $10^{9}\,M_{\odot}$3, providing evidence for state transitions (thin disk to RIAF) and the disappearance of broad line regions at critical low accretion rates (Lusso et al., 24 Feb 2025). Faded changing-look quasars at $10^{9}\,M_{\odot}$4 show SED evolution consistent with the low/hard states of X-ray binaries, indicating physical continuity in accretion flow transformation across mass scales (Gilbert et al., 3 Aug 2025).

5. Circumnuclear Environment and Feedback

LLAGN exhibit distinct feedback and ISM-processing signatures compared to their higher-luminosity cousins. Multiwavelength observations resolve:

• Powerful parsec-scale mass outflows and jets, producing fast (v$10^{9}\,M_{\odot}$5–$10^{9}\,M_{\odot}$6 km s⁻¹), shock-excited, coronal line—emitting gas. In NGC 1386, the outflow rate is $10^{9}\,M_{\odot}$7 yr⁻¹, with kinetic power $10^{9}\,M_{\odot}$8 erg s⁻¹, often $10^{9}\,M_{\odot}$9\% of $L_\mathrm{bol} \lesssim 10^{42{-}43}$0 (Rodríguez-Ardila et al., 2017).
• Nuclear feedback operates in the “kinetic mode”: LLAGN with $L_\mathrm{bol} \lesssim 10^{42{-}43}$1–$L_\mathrm{bol} \lesssim 10^{42{-}43}$2 exhibit strong shocks, which dominate the excitation of both ionized and molecular phases on $L_\mathrm{bol} \lesssim 10^{42{-}43}$3 pc scales, as established by JWST and Spitzer spectroscopy of diagnostic fine-structure lines ([Ne V], [O IV]) and H$L_\mathrm{bol} \lesssim 10^{42{-}43}$4 bands (Zhang et al., 2022, Zhang et al., 12 Jan 2026).
• The kinetic feedback imprints on the ISM: small PAH molecules are destroyed—evidenced by suppressed PAH $L_\mathrm{bol} \lesssim 10^{42{-}43}$5 and enhanced $L_\mathrm{bol} \lesssim 10^{42{-}43}$6 band ratios, with most surviving PAH population consisting of large, neutral molecules ($L_\mathrm{bol} \lesssim 10^{42{-}43}$7) (Zhang et al., 2022, Zhang et al., 12 Jan 2026).
• H$L_\mathrm{bol} \lesssim 10^{42{-}43}$8 rotational transitions indicate slow, jet-driven molecular shocks at velocities $L_\mathrm{bol} \lesssim 10^{42{-}43}$9 km s⁻¹, incompatible with pure radiative (PDR or XDR) excitation. A few percent of the jet mechanical power suffices to power the observed H$L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$0 emission (Zhang et al., 12 Jan 2026).
• Compared to Seyferts and quasars, which manifest “radiative mode” UV/X-ray feedback and strong suppression of all PAH bands, LLAGN display a kinetic–mechanical regime that directly modifies their immediate molecular environment (Zhang et al., 2022).

6. Demographics, Host Properties, and Environmental Dependence

LLAGN are prevalent in both bulge and disk galaxies across the local universe. Narrow-line, low-mass LLAGN (FWHM Hβ < 2200 km s⁻¹) show sub-Eddington accretion, weak Fe II, strong optical variability, but are physically distinct from classical NLSy1 galaxies, forming a “low-tail” in 4DE1 parameter space and representing a low-accretion analogue within the broader narrow-line AGN family (Sharma et al., 5 Dec 2025).

Variability- and emission-line-selected LLAGN are found in diverse environments, but modern statistical surveys report a $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$1 higher fraction of LLAGN in galaxy clusters relative to the field, contrasting with the tendency for luminous AGN to avoid clusters (Yuk et al., 29 Jan 2025). Proposed triggering mechanisms at low $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$2 include tidal interactions with the cluster potential, cooling flows, and minor mergers rather than major mergers or violent gas inflows.

Multiwavelength studies indicate that nuclear star formation is suppressed in LLAGN hosts, linked to low dense-gas fractions and mechanical heating by weak jets or AGN-driven feedback (Boeker et al., 2011). The overall LLAGN fraction in volume-limited surveys is $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$3 at $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$4, constraining the local AGN duty cycle and black-hole growth scenarios (Yuk et al., 2022).

7. Future Prospects and Observational Diagnostics

The unexpected strength of nuclear K-band and MIR emission in LLAGN, alongside tight correlations with X-ray luminosity, highlight the efficacy of IR selection, particularly with JWST/NIRSpec and MIRI (Dumont et al., 2019). JWST enables spatially resolved identification of hot-dust AGN emission and separation from stellar backgrounds down to $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$5–$L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$6 erg s⁻¹ at distances up to tens of Mpc.

PAH band ratios and H$L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$7 excitation provide a robust mid-IR diagnostic to distinguish kinetic-mode feedback in LLAGN from radiative-mode processes in more powerful AGN (Zhang et al., 2022, Zhang et al., 12 Jan 2026). The employment of variability-based, emission-line, and multiwavelength selection, coupled with high-resolution IFU spectroscopy, is essential to correct for misclassification of Compton-thick, obscured, or host-dominated LLAGN and to avoid significant underestimates in low-$L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$8 AGN space density (Lambrides et al., 2020, Ding et al., 2018).

The transition between radiative- and kinetic-dominated feedback as a function of $L_\mathrm{2-10\,keV} \lesssim 10^{40{-}41}$9, the survival and dissipation timescale of the dusty torus, the structure of accretion flows, and the nature of accretion state changes remain central topics for future theoretical and observational work—areas in which deep, high-angular-resolution, multi-epoch IR and X-ray spectroscopy will be decisive.

---

Selected Key Quantitative Results and Correlations in LLAGN:

Observable or Ratio	LLAGN Regime	Typical Value/Trend	Reference
$\lambda_\mathrm{Edd} \sim 10^{-6}$0	LLAGN	$\lambda_\mathrm{Edd} \sim 10^{-6}$1 erg s⁻¹	(Dumont et al., 2019, Annuar et al., 2020)
$\lambda_\mathrm{Edd} \sim 10^{-6}$2	LLAGN	$\lambda_\mathrm{Edd} \sim 10^{-6}$3–$\lambda_\mathrm{Edd} \sim 10^{-6}$4	(Lusso et al., 24 Feb 2025, Zhang et al., 2022)
$\lambda_\mathrm{Edd} \sim 10^{-6}$5	NIR-to-X ratio ($\lambda_\mathrm{Edd} \sim 10^{-6}$6)	$\lambda_\mathrm{Edd} \sim 10^{-6}$7 exceeds $\lambda_\mathrm{Edd} \sim 10^{-6}$8 as $\lambda_\mathrm{Edd} \sim 10^{-6}$9 drops	(Dumont et al., 2019)
MIR–X-ray correlation	log–log slope across $10^{-2}$0 erg/s	$10^{-2}$1; no dependence on $10^{-2}$2	(Asmus et al., 2011)
Silicate ($10^{-2}$3)	High in high $10^{-2}$4 LLAGN	$10^{-2}$5–$10^{-2}$6 cm² (order mag above Seyferts)	(Mason et al., 2012)
$10^{-2}$7	LLAGN (face-on disks)	$10^{-2}$8 ($10^{-2}$9 confidence, elevated)	(Boeker et al., 2011)
Dense H$R_\mathrm{tr} \sim 30{-}300\,R_S$0 (HCN/CO ratio)	LLAGN nuclei	$R_\mathrm{tr} \sim 30{-}300\,R_S$1 (low density)	(Boeker et al., 2011)
Nuclear X-ray photon index	LLAGN (variable)	Mean $R_\mathrm{tr} \sim 30{-}300\,R_S$2 (soft, unobscured)	(Young et al., 2012, Ding et al., 2018)
BLR radius–luminosity scaling	Down to low $R_\mathrm{tr} \sim 30{-}300\,R_S$3	$R_\mathrm{tr} \sim 30{-}300\,R_S$4, $R_\mathrm{tr} \sim 30{-}300\,R_S$5–0.19 dex	(Bentz et al., 2013)
PAH band ratios	LLAGN (shock-dominated)	$R_\mathrm{tr} \sim 30{-}300\,R_S$6, $R_\mathrm{tr} \sim 30{-}300\,R_S$7, $R_\mathrm{tr} \sim 30{-}300\,R_S$8	(Zhang et al., 2022, Zhang et al., 12 Jan 2026)
Optical variability, log σ	Low-mass LLAGN	median log σ ≈ –0.68	(Sharma et al., 5 Dec 2025)

---