[Paper Review] An Investigation of Spectral Line Stacking Techniques and Application to the Detection of HC$_{11}$N

Ryan A. Loomis, Andrew M. Burkhardt+11 more|arXiv (Cornell University)|Sep 24, 2020

Astrophysics and Star Formation Studies4 citations

TL;DR

This paper presents a novel hybrid method combining Markov Chain Monte Carlo (MCMC) inference with spectral line stacking and matched filtering to detect weak, rotationally diluted lines of large interstellar molecules in sparse single-dish spectra. Applied to GOTHAM GBT data, the technique enabled the robust detection of HC₁₁N—the largest cyanopolyyne ever found in space—confirming its column density within previous upper limits.

ABSTRACT

As the inventory of interstellar molecules continues to grow, the gulf between small species, whose individual rotational lines can be observed with radio telescopes, and large ones, such as polycyclic aromatic hydrocarbons (PAHs) best studied in bulk via infrared and optical observations, is slowly being bridged. Understanding the connection between these two molecular reservoirs is critical to understanding the interstellar carbon cycle, but will require pushing the boundaries of how far we can probe molecular complexity while still retaining observational specificity. Toward this end, we present a method for detecting and characterizing new molecular species in single-dish observations toward sources with sparse line spectra. We have applied this method to data from the ongoing GOTHAM (GBT Observations of TMC-1: Hunting Aromatic Molecules) Green Bank Telescope (GBT) large program, discovering six new interstellar species. In this paper we highlight the detection of HC$_{11}$N, the largest cyanopolyyne in the interstellar medium.

Motivation & Objective

• To address the challenge of detecting large, rotationally complex interstellar molecules whose individual spectral lines are too weak to be observed due to rotational energy level dilution.
• To improve detection efficiency and parameter characterization of weak molecular emissions in line-sparse, single-dish spectral surveys.
• To develop a robust, statistically grounded method for identifying new interstellar species when no single transition is detectable above noise.
• To apply this method to the GOTHAM GBT survey to identify and characterize complex molecules in TMC-1, a prototypical source for molecular complexity.

Proposed method

• The method integrates MCMC inference to model spectral line parameters (velocity, FWHM, optical depth, column density) while accounting for uncertainties and covariances.
• It applies spectral line stacking by weighting data using a best-fit model derived from MCMC, enhancing signal-to-noise for weak, multiple-line species.
• A matched filter is constructed from the stacked model and applied to the stacked data to assess statistical significance of the detection.
• The technique uses a priori spectral catalogs and telescope response functions, with source parameters (e.g., size, velocity) varied during MCMC fitting.
• Beam dilution and source-size fitting are performed to account for spatial inhomogeneity and beam-averaged emission effects.
• The approach is validated using simulated data and applied to real GBT observations of TMC-1 at X-band (7.8–29.9 GHz).

Experimental results

Research questions

• RQ1 Can MCMC-informed spectral stacking detect interstellar molecules with no individually detectable transitions due to rotational line dilution?
• RQ2 How can the signal-to-noise ratio be maximized for weak, broad, and low-intensity lines of large molecules in sparse spectral surveys?
• RQ3 What is the impact of source size and beam dilution on the inferred column density and emission parameters of large molecules?
• RQ4 Can this method reliably recover known molecular abundances, such as for HC₁₁N, from real observational data?

Key findings

• The method successfully detected HC₁₁N in TMC-1, the largest cyanopolyyne identified in the interstellar medium to date.
• The derived column density of HC₁₁N is consistent with the upper limit previously reported in Loomis et al. (2016), confirming its presence without prior detection.
• The analysis revealed that larger species like HC₁₁N may have broader spatial distributions than smaller molecules, with source components showing non-coincident spatial profiles.
• For optically thin species such as HC₉N, separate component fitting indicated varying source sizes, suggesting non-uniform spatial distribution or beam filling factor effects.
• The technique enabled the detection of six new interstellar species in the GOTHAM survey, demonstrating its effectiveness in probing molecular complexity beyond small molecules.
• The method's robustness was validated through simulations and application to real data, with open-source code publicly available for replication and extension.