🔭 astrophysics

BL Lac host galaxies: how to systematically characterise them in optical-NIR spectroscopy

This paper proposes a systematic method using the QSFit software to distinguish between elliptical and spiral host galaxies of BL Lacertae objects by analyzing optical-NIR spectra, demonstrating that the technique can effectively classify sources with jets fainter than $L_\gamma\sim10^{46}$ erg/s while acknowledging potential biases for certain galaxy types.

Original authors: Gaia Delucchi, Tullia Sbarrato, Giorgio Calderone, Chiara Righi, Silvano Tosi, Boris Sbarufatti

Published 2026-03-18

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Gaia Delucchi, Tullia Sbarrato, Giorgio Calderone, Chiara Righi, Silvano Tosi, Boris Sbarufatti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to identify the type of car parked in a very dark garage. But there's a catch: sitting right on top of the car is a blindingly bright spotlight that shines so intensely it washes out everything else. You can't see the car's shape, color, or features; you only see the glare of the light.

This is exactly the problem astronomers face when studying BL Lac objects.

The Cosmic Problem: The Spotlight vs. The Car

BL Lac objects are a special kind of super-massive black hole at the center of a galaxy. They shoot a jet of energy (like a cosmic laser) straight at Earth. This jet is so bright that it drowns out the light from the galaxy itself.

For decades, astronomers have assumed these galaxies are all Elliptical (big, round, "old" galaxies with no new stars). But because the jet is so loud, they couldn't hear the galaxy's voice to confirm this. They were guessing based on the silence.

The big question is: Are these galaxies actually Elliptical, or could some of them be Spiral (like our Milky Way, with swirling arms and new stars)?

The Solution: A "Spectral" Detective Kit

The authors of this paper (Gaia Delucchi and her team) wanted to build a simple, reliable way to figure this out without needing a super-powerful telescope for every single object. They created a digital "detective kit" using a free software tool called QSFit.

Here is how their method works, using a simple analogy:

1. Building a "Fake Universe" (The Training Ground)

Before testing their method on real stars, they had to teach their software what to look for. They created 3,500 fake galaxies on a computer.

They mixed different types of "cars" (Elliptical and Spiral galaxies).
They added different sizes of "spotlights" (BL Lac jets) of varying brightness.
They even added "static noise" to the signal, just like a bad radio connection, to make it realistic.

2. The Two-Step Test (The "Taste Test")

The software looks at the light from a galaxy and tries to fit it into two different molds:

Mold A: "This is an Elliptical galaxy with a jet."
Mold B: "This is a Spiral galaxy with a jet."

The software calculates a "score" (a statistical number called $\chi^2$ ) for how well the real data fits each mold.

If the data fits Mold A much better, the galaxy is likely Elliptical.
If the data fits Mold B much better, the galaxy is likely Spiral.

3. The Magic Ratio (The "Decision Line")

The authors found a simple trick. They take the score from the Elliptical mold and divide it by the score from the Spiral mold. Let's call this the R-Ratio.

R < 1: The Elliptical mold fits better. (It's an Elliptical galaxy).
R > 1: The Spiral mold fits better. (It's a Spiral galaxy).
R is right in the middle: The signal is too weak or the galaxy is too similar to the other type to tell for sure. This is the "Confusion Zone."

The Catch: How Bright is the Spotlight?

The method works beautifully, but only under one condition: The spotlight (the jet) cannot be too blinding.

If the jet is too powerful (specifically, if it's brighter than a certain threshold), it completely drowns out the galaxy, and the software can't tell the difference. The authors found that for about 80% of known BL Lac objects, the jet isn't quite bright enough to hide the galaxy completely. For these 80%, their method works like a charm.

Why This Matters

Previously, astronomers had to take very expensive, high-resolution photos to see the galaxy shape. This new method only requires a standard spectrum (a rainbow of light) and two quick computer runs.

The Good News: It confirms that most BL Lacs are indeed in Elliptical galaxies, but it opens the door to finding the rare "Spiral" ones that were previously hidden.
The Future: With new telescopes like JWST and Euclid coming online, this method will allow astronomers to quickly sort through thousands of galaxies, finally hearing the "voice" of the galaxy behind the blinding jet.

In a Nutshell

The authors built a digital filter that separates the "noise" of a super-bright jet from the "signal" of the host galaxy. By running the data through two different filters (one for Elliptical, one for Spiral) and comparing the results, they can now systematically tell us what kind of galaxy is hiding behind the cosmic spotlight. It's like finally being able to see the car in the dark garage, even with the spotlight on top of it.

1. Problem Statement

Active Galactic Nuclei (AGN), particularly BL Lacertae (BL Lac) objects, are characterized by relativistic jets aligned with the observer's line of sight. This alignment causes the non-thermal jet emission to dominate the electromagnetic spectrum (from radio to $\gamma$ -rays), often completely obscuring the thermal emission from the host galaxy.

The Challenge: Determining the morphological type of the host galaxy (elliptical vs. spiral) is crucial for understanding the co-evolution of Supermassive Black Holes (SMBHs) and their environments. However, current knowledge is sparse and biased. Historically, BL Lacs have been assumed to reside exclusively in giant elliptical galaxies, supported by limited observations.
The Gap: There is no systematic method to distinguish between elliptical and spiral hosts in BL Lac spectra, especially when the jet flux overwhelms the host. The only known spiral candidate (PKS 1413+135) is debated, and previous studies often assume ellipticity to derive photometric redshifts, potentially introducing biases.

2. Methodology

The authors propose a systematic, automated spectroscopic analysis method using the open-source software QSFit. The methodology involves three main stages:

A. Synthetic Data Generation

To test and calibrate the method, the authors created a large sample of 3,500 synthetic spectra covering the UV-to-NIR range (1000–10,000 Å).

Components: Each spectrum is a combination of:
- Jet Emission: Phenomenological Spectral Energy Distributions (SEDs) from Ghisellini et al. (2017) representing 5 different $\gamma$ -ray luminosity bins.
- Host Galaxies: 7 templates from the Swire collection (M. Polletta et al. 2007), including 3 elliptical galaxies (ages: 2, 5, 13 Gyr) and 4 spiral types (Sa, Sb, Sc, Sd).
Parameters: The sample spans redshifts $z = 0.1$ to $2.0$ and various host-to-jet luminosity ratios ( $L_{host}^R / L_{BL Lac}^R$ ).
Noise: Gaussian noise (standard deviation $\sigma = 10\%$ of flux density) was added to simulate observational conditions. Tests were also run with 1%, 5%, and 15% noise levels.

B. Modified QSFit Analysis

The standard QSFit software (designed for quasars with emission lines) was modified for BL Lacs:

Removal of Emission Lines: All broad/narrow emission lines and pseudo-continua (Fe, Balmer) were removed, as BL Lacs typically lack them.
Continuum Model: The jet emission was modeled as a single power-law.
Dual-Recipe Approach: Each synthetic spectrum was analyzed twice using two distinct "recipes":
1. Recipe 1: Power-law + Ell5 template (5 Gyr old elliptical galaxy).
2. Recipe 2: Power-law + Sb template (Sb spiral galaxy).
Output: For each run, the software returns the reduced chi-square statistic ( $\chi^2_{red}$ ) and normalization parameters.

C. Classification Metric

The core of the method is the ratio of the fit statistics from the two recipes:
$R = \frac{\chi^2_{Ell5}}{\chi^2_{Sb}}$

$R < 1$ : The elliptical template provides a better fit.
$R > 1$ : The spiral template provides a better fit.
$R \approx 1$ : The fit is ambiguous (host is diluted by the jet or templates are indistinguishable).

3. Key Contributions

Systematic Algorithm: Development of a simple, reproducible algorithm requiring only two QSFit runs per spectrum to classify host morphology.
Quantification of Uncertainty: Definition of a "Confusion Region" bounded by limits $\alpha$ and $\beta$ , where classification is statistically unreliable.
Luminosity Thresholds: Establishment of specific $\gamma$ -ray and R-band luminosity limits below which host galaxy characterization is feasible.
Validation on Real Data: Application of the method to a sample of 240 BL Lacs from the Fermi 4LAC-DR3 and SDSS-DR17 catalogs.

4. Key Results

A. The Classification Metric ( $R$ )

The ratio $R$ effectively discriminates host types based on the host-to-jet luminosity ratio ( $L_{host}^R / L_{BL Lac}^R$ ):

Dominant Jet ( $L_{host} \ll L_{jet}$ ): $R \approx 1$ . The host is undetectable; classification is impossible.
Comparable/Dominant Host: As the host becomes visible, $R$ $R$ deviates from 1.
- Elliptical Hosts: Tend to yield $R < 1$ .
- Spiral Hosts: Tend to yield $R > 1$ .

B. The Confusion Region

The authors defined a region of uncertainty where misclassification occurs:

Limits: $\beta = 0.97$ (lower) and $\alpha = 1.09$ (upper) for 10% noise.
Asymmetry: The region is asymmetric. Elliptical galaxies (specifically the 13 Gyr template, Ell13) are more likely to be misclassified as spirals (falling in the $R > 1$ region) than vice versa. This is attributed to the excess blue-UV flux in older elliptical templates mimicking spiral signatures.
Coverage: Approximately 70% of the synthetic sample falls within this confusion region, primarily due to low host-to-jet ratios.

C. Luminosity Limits

The method is only effective for BL Lacs with jets fainter than a specific threshold:

$\gamma$ -ray Limit: $L_\gamma \lesssim 10^{46}$ erg/s. Above this, the jet completely swamps the host.
R-band Limit: For a typical host ( $M_R = -22.9$ ), the host can be identified if $L_{BL Lac}^R \lesssim 4.4 \times 10^{45}$ erg/s.
Feasibility: About 79-80% of BL Lacs in the 3LAC sample fall below these limits, making them candidates for this analysis.

D. Real Data Validation

Applied to 240 SDSS BL Lacs ( $L_\gamma < 10^{46}$ erg/s):

14 sources fell outside the confusion region.
11 sources were classified as Elliptical; 3 as Spiral.
Cross-Check: Sources near the spiral boundary were found to have irregular/interacting morphologies or bluer stellar populations, consistent with the method's logic that bluer hosts shift $R$ toward the spiral region.

5. Significance and Future Outlook

Paradigm Shift: This work challenges the assumption that all BL Lacs reside in elliptical galaxies, providing a tool to systematically test for spiral hosts.
Readiness for New Facilities: The method is designed for upcoming data from JWST, Euclid, and LSST.
- Euclid NISP: The authors note that Euclid's Near-Infrared Spectro-Photometer (NISP) could utilize 2D spectroscopy to spatially separate the jet (center) from the host (edges), potentially overcoming the spectral confusion entirely.
Conservative Approach: The authors emphasize that the "Confusion Region" must be respected. Misclassifications are asymmetric, with faint spirals being harder to identify than faint ellipticals.
Impact: By enabling the distinction of host types, this method opens the door to studying the evolutionary tracks of AGN, the role of accretion efficiency, and the environmental factors governing jet formation.

In summary, the paper provides a robust, statistically grounded framework to extract host galaxy information from the overwhelming glare of BL Lac jets, defining clear boundaries for when such analysis is possible and reliable.