---
title: 'GOTHAM Survey: TMC-1 Aromatic Chemistry'
url: https://www.emergentmind.com/topics/gotham-large-program-survey
type: topic
---

# GOTHAM Survey: TMC-1 Aromatic Chemistry

Searching arXiv for the cited GOTHAM survey papers and related overview work.
The **GOTHAM Large Program Survey**—**GBT Observations of TMC-1: Hunting for Aromatic Molecules**—is a deep, broadband centimeter-wavelength spectral line survey of the cyanopolyyne peak in the Taurus Molecular Cloud 1 (TMC-1) conducted with the 100 m Green Bank Telescope. It was launched to investigate the depth and breadth of aromatic chemistry in the interstellar medium at the earliest stages of star formation, following the detection of benzonitrile in TMC-1, and has since developed into a broader program for constraining molecular inventories, line identifications, formation pathways, and astrochemical model deficiencies in a cold, chemically rich dark cloud [2008.12349]. Across the cited works, GOTHAM is characterized by large spectral bandwidth, high sensitivity, high spectral resolution, laboratory-supported spectroscopy, Bayesian line analysis, velocity stacking, and matched filtering; these together have enabled both secure detections of new interstellar molecules and stringent upper limits on chemically informative non-detections [2008.12349].

## 1. Survey design, target selection, and observational regime

GOTHAM targets the **cyanopolyyne peak** of **TMC-1**, a cold, quiescent, chemically rich dark cloud that is nearby, well studied, and frequently treated as a benchmark source for astrochemical models [2008.12349]. The program rationale is tied to four scientific questions stated in the overview paper: whether aromatic species beyond benzonitrile exist in TMC-1, what precursors may build aromatic species through bottom-up chemistry, whether this chemistry is unique to TMC-1 or widespread in dark clouds, and whether such species survive subsequent star and planet formation [2008.12349].

The observational architecture is based on the Green Bank Telescope with X-, K-, and Ka-band coverage. In the early overview, the survey spans **8.0–11.6 GHz**, **18.0–27.5 GHz**, and **26.0–39.5 GHz**, using the VEGAS backend and ON-OFF position-switched observations with **1.4 kHz per channel** spectral resolution [2008.12349]. Later GOTHAM analyses describe nearly continuous coverage across **7.9–36.4 GHz** or comparable ranges, with uniform frequency resolution of **1.4 kHz**, corresponding to roughly **0.05–0.01 km s\(^{-1}\)** in velocity, and high sensitivity reaching RMS noise of a few mK in parts of the band [2303.12221; 2409.16435]. A later inventory paper extends the survey description to **3.9–36.4 GHz**, with **\(\sim 29\) GHz** of spectral bandwidth and **1.43 kHz** resolution [2509.06256].

This instrumental regime is central to the program’s scientific character. TMC-1 exhibits narrow lines and low excitation temperatures, so high frequency precision and fine channelization are necessary not merely for sensitivity but for reliable rest-frequency matching, multi-component decomposition, and statistical recovery of molecules whose individual transitions remain below the single-line detection threshold [2008.12349; 2409.16435]. A plausible implication is that GOTHAM is not only a discovery survey but also a precision spectroscopy platform for centimeter-wave astrochemistry.

## 2. Data reduction, calibration, and statistical analysis framework

The overview paper emphasizes detailed observing and reduction procedures, including archival data reuse, gbtidl-based processing, Doppler corrections, conversion to atmosphere-corrected \(T_A^*\), polarization averaging where applicable, manual RFI cleaning, baseline fitting with polynomial orders between 1 and 20, and noise-weighted averaging [2008.12349]. Later work on the TMC-1 inventory reports a **Python-based, fully automated data reduction pipeline**—the **GOTHAM Spectral Pipeline**—for single-beam, position-switched observations, with automated RFI flagging, per-channel Doppler tracking correction, dynamic zenith opacity correction using real-time weather data, and rejection of integrations with abnormal system temperatures [2509.06256]. This indicates an evolution from early heterogeneous reduction toward a scalable, reproducible calibration framework suited to a survey of very large bandwidth and integration time.

A consistent analysis layer across the program is **Bayesian forward modeling with Markov Chain Monte Carlo (MCMC)**. In different papers, the fit parameters include column density, excitation temperature, linewidth, source size, and one to four velocity components depending on the species and model choice [2008.12349; 2303.12221; 2409.16435; 2509.06256]. The later inventory work gives the likelihood explicitly as

\[
l = \sum_{i} \left[ -\ln \left(\sigma_i \sqrt{2\pi}\right) - \frac{1}{2}\frac{(x_i-m_i)^2}{\sigma_i^2} \right],
\]

with \(x_i\) the observed intensity, \(m_i\) the model intensity, and \(\sigma_i\) including noise and systematic uncertainties [2509.06256].

Radiative-transfer-based column density modeling is likewise explicit in the program literature. For example, the HCCCHO study gives the optically thin line-analysis expression

\[
N = \frac{3k}{8\pi^3 \nu S \mu^2} \cdot \frac{Q(T_{ex})}{g_u} \exp\left(\frac{E_u}{kT_{ex}}\right) \cdot \int T_{mb}\,dv
\]

for emission-line column density estimation [2409.16435]. The inventory study further provides expressions for optical depth and observed brightness temperature, together with beam dilution and partition function corrections [2509.06256].

For weak molecules, GOTHAM systematically combines **velocity stacking** and **matched filtering**. Spectra centered on predicted transitions are aligned in velocity, stacked with weights based on predicted intensity and RMS noise, and then cross-correlated against a simulated stacked spectrum to yield an impulse response used as a detection statistic [2008.12349; 2101.05858; 2104.15117; 2209.06851; 2301.07760; 2303.12221]. The repeated use of this framework across chemically diverse targets suggests that statistical signal extraction is a defining methodological contribution of the survey rather than an auxiliary technique.

## 3. Molecular detections and the expansion of the TMC-1 inventory

GOTHAM has been associated with a sequence of first detections and targeted identifications spanning unsaturated chains, substituted nitriles, pure hydrocarbons, sulfur-bearing species, and molecular anions.

The survey overview reported the **first interstellar detection of propargyl cyanide** \(\mathrm{HCCCH_2CN}\), supported by new laboratory spectroscopy with over **110 hyperfine-split features** measured between **5–40 GHz** at \(\sim 2\) kHz uncertainty, and an **18\(\sigma\)** matched-filter detection in TMC-1 with a total column density of \(9.2^{+1.3}_{-0.7}\times10^{11}\ \mathrm{cm}^{-2}\) [2008.12349].

Early GOTHAM analyses also reported the discovery of **trans-(E)-cyanovinylacetylene** and **vinylcyanoacetylene**, with matched-filter significances of **8.0\(\sigma\)** and **5.5\(\sigma\)**, respectively. Their MCMC-derived column densities are \((2.90^{+0.41}_{-0.40}) \times 10^{11}\ \mathrm{cm}^{-2}\) for trans-(E)-cyanovinylacetylene and \((1.87^{+0.37}_{-0.38}) \times 10^{11}\ \mathrm{cm}^{-2}\) for vinylcyanoacetylene; the trans-(Z) isomer remained undetected with an upper limit of \(<2\times10^{11}\ \mathrm{cm}^{-2}\) [2101.05858].

The survey subsequently yielded the **first interstellar detection of indene** \((c\text{-}\mathrm{C_9H_8})\), described as the **first interstellar detection of a pure hydrocarbon PAH**, with a **5.7\(\sigma\)** matched-filter detection and total column density \(9.6^{+4.3}_{-1.6}\times10^{12}\ \mathrm{cm}^{-2}\) [2104.15117]. A later study detected **2-cyanoindene** at **6.3\(\sigma\)**, with \(N_T = 2.10^{+0.60}_{-0.46}\times10^{11}\ \mathrm{cm}^{-2}\), while refining the indene column density to \(9.04^{+0.96}_{-0.96}\times10^{12}\ \mathrm{cm}^{-2}\) [2209.06851]. These measurements enabled the first direct hydrocarbon:CN-substituted comparison for a PAH-related system in the same source.

GOTHAM also reported the **first astronomical detection of the interstellar anion \(\mathrm{C_{10}H^-}\)**. Using velocity-stacked data and matched filtering, the anion was detected at **\(>9\sigma\)** confidence with column density \(4.04^{+10.67}_{-2.23}\times10^{11}\ \mathrm{cm}^{-2}\), while the neutral radical \(\mathrm{C_{10}H}\) was only tentatively detected at **\(\sim 3.2\sigma\)** with \(2.02^{+0.68}_{-0.82}\times10^{11}\ \mathrm{cm}^{-2}\) [2301.07760].

Further detections include **\(E\)-1-cyano-1,3-butadiene**, identified at **5.1\(\sigma\)** with total column density \(3.8^{+1.0}_{-0.9}\times10^{10}\ \mathrm{cm}^{-2}\) [2303.12221], and **cyclopropenethione** \((c\text{-}\mathrm{C_3H_2S})\), detected with \(N_T = 5.72^{+2.65}_{-1.61}\times10^{10}\ \mathrm{cm}^{-2}\) and \(T_{ex} = 4.7^{+1.3}_{-1.1}\ \mathrm{K}\), completing the three low-energy isomers of the \(\mathrm{C_3H_2S}\) family in TMC-1 [2501.06343].

The later inventory analysis synthesizes these and other results into a survey-scale census, reporting constrained column densities for **102 molecular species** in TMC-1, comprising **75 main isotopic species**, **20 \(^{13}\)C-substituted species**, and **7 deuterium-substituted species**, with the detected gas-phase inventory populated primarily by **unsaturated hydrocarbons** [2509.06256].

## 4. Aromatics, cyano proxies, and bottom-up chemical inference

A major theme of GOTHAM is the use of nitrile-functionalized molecules as observational proxies for otherwise radio-inactive hydrocarbons and PAHs. This strategy follows from the difficulty that fully symmetric pure hydrocarbons may lack a permanent dipole moment and therefore a detectable rotational spectrum, whereas CN substitution imparts strong dipole moments and rotational observability [2209.06851].

The indene and 2-cyanoindene pair provides a central case. The measured abundance ratio

\[
\frac{N_T(\mathrm{C_9H_8})}{N_T(\mathrm{2\text{-}C_9H_7CN})} \approx 43
\]

with uncertainty range **27–61** is described as the first direct observation of the ratio of a cyano-substituted PAH to its pure hydrocarbon counterpart in the same source [2209.06851]. That study concludes that while astrochemical models are not yet sufficiently accurate to reproduce the absolute abundances of these species, they do a good job at predicting hydrocarbon:CN-substituted ratios, thereby supporting **-CN tagged species as excellent proxies** for fully symmetric counterparts [2209.06851].

The survey has also used nitriles as proxies for non-polar chain precursors. The detection of \(E\)-1-cyano-1,3-butadiene constrains the abundance of **1,3-butadiene**, a non-polar hydrocarbon that cannot be directly detected by radio observations. Using a modeled CN-functionalization ratio, the work infers

\[
N_T(\ce{CH_2CHCHCH_2}) \approx 3.94\times10^{11}\ \mathrm{cm}^{-2},
\]

and discusses 1,3-butadiene as a key precursor in laboratory, theoretical, and astrophysical models for bottom-up aromatic-ring formation [2303.12221]. The same paper gives the reaction

\[
\ce{CH_2CHCHCH_2 + C_2H \rightarrow C_6H_6 + H}
\]

as a route to benzene and related products [2303.12221].

The indene paper frames pure hydrocarbon PAHs as chemically decisive because they provide a direct handle on bottom-up and top-down aromatic formation scenarios in cold clouds [2104.15117]. It states that indene is the most abundant organic ring detected in TMC-1 to date and that its abundance exceeds NAUTILUS predictions by several orders of magnitude [2104.15117]. This suggests that existing aromatic formation networks remain incomplete even when precursor chains and cyano derivatives are better constrained.

A later statistical inventory strengthens the aromatic context by reporting **ten individual aromatic molecules** in the GOTHAM observations and quantifying their contribution to the gas-phase carbon budget: **0.011%** when including CO and **6%** when excluding CO [2509.06256]. Nine of the ten are CN-functionalized aromatics, which is presented as evidence of the practical role of nitriles as radio proxies [2509.06256].

## 5. Spectroscopic precision, laboratory support, and catalog revision

A defining feature of GOTHAM is its dependence on laboratory spectroscopy and on the revision of public spectral catalogs when survey data reveal inconsistencies. The HCCCHO study is exemplary. It combined astronomical frequencies measured at sufficiently high spectral resolution with new laboratory measurements at **\(\sim 2\) kHz** resolution to produce a new global fit to the rotational spectrum of propynal [2409.16435]. The fit comprised **2674 unique transitions**—or **3230 with blended lines**—covering **\(J=1\) to \(J=100\)** and **\(K_a = 0\)–15**, and achieved astronomical frequency matches **within \(<1\) kHz** for the key \(J=2\text{--}1\) and \(3\text{--}2\) transitions [2409.16435].

Using this revised catalog, MCMC analysis yielded a best-fit total HCCCHO column density of

\[
7.28^{+4.08}_{-1.94} \times 10^{12}\ \mathrm{cm}^{-2}
\]

with an excitation temperature of approximately **3 K**, around **a factor of 5 times larger** than reported in previous studies [2409.16435]. The paper explicitly states that care is needed when using publicly available spectral catalogs to characterize astronomical spectra and emphasizes the growing importance of such catalogs as next-generation facilities come online [2409.16435].

Comparable laboratory support appears throughout the survey. The propargyl cyanide detection required improved hyperfine-resolved rotational measurements to predict rest frequencies within **0.03 km s\(^{-1}\)** [2008.12349]. The indene work relied on **189 transitions** measured with **2 kHz accuracy** over roughly **2–40 GHz** [2104.15117]. The cyanoindene study used cavity Fourier-transform microwave spectroscopy between **6–40 GHz** for five isomers [2209.06851]. The c-\(\mathrm{C_3H_2S}\) detection was enabled by laboratory characterization of a molecule for which no previous experimental spectra existed [2501.06343].

This recurring pattern shows that GOTHAM is structurally coupled to laboratory molecular spectroscopy. A plausible implication is that the survey’s discovery space is determined as much by catalog fidelity and laboratory rest-frequency accuracy as by telescope sensitivity alone.

## 6. Non-detections, upper limits, and constraints on chemical selectivity

GOTHAM’s scientific value is not confined to positive detections. The survey’s sensitivity and Bayesian multi-line analysis have also produced stringent upper limits that constrain chemical formation and destruction pathways.

The most systematic non-detection study in the cited set is the search for **nitrogen-, oxygen-, and sulfur-bearing heterocycles** in TMC-1 [2204.04744]. Using deep GOTHAM observations over **7.906–35.434 GHz** and Bayesian forward modeling in **molsim**, the search targeted **19 heterocycles** and found **no evidence for the presence of any heterocyclic species** [2204.04744]. The MCMC approach used all lines covered by the survey for each species and reported the **97.5th percentile of the posterior distribution** as a **2\(\sigma\)** upper limit [2204.04744]. The paper states that the MCMC multi-line analysis systematically yielded more stringent upper limits than traditional single-line frequentist methods [2204.04744].

The same work interprets the absence of detectable heterocycles despite the confirmed presence of cyclic and polycyclic hydrocarbons as evidence that heterocycles are **significantly less abundant** than pure carbocyclic counterparts in TMC-1 [2204.04744]. Specific examples given are that pyridine is at most half as abundant as benzonitrile, while pyrrole, thiophene, and furan are depleted by one to two orders of magnitude relative to cyclopentadiene [2204.04744]. The paper further hypothesizes inefficient low-temperature formation pathways, lack of suitable precursors, and possible destruction effects, though it notes that dense clouds such as TMC-1 are UV-shielded [2204.04744].

Upper limits are also chemically informative in chain and nitrile chemistry. The methyl-chain study found no significant emission from **methyltetraacetylene \(\mathrm{CH_3C_8H}\)** and set a **95% confidence** upper limit of \(N_T < 9.8\times10^{10}\ \mathrm{cm}^{-2}\), interpreting this as a sharp abundance drop-off after \(\mathrm{CH_3C_6H}\) [2110.13152]. The trans-(Z)-cyanovinylacetylene isomer likewise remained undetected with \(N_{\rm col} < 2\times10^{11}\ \mathrm{cm}^{-2}\) [2101.05858]. Ethyl cyanide, considered as a hydrogenation product of vinyl cyanide, showed only **4.17\(\sigma\)** significance and was treated as an upper limit of \(<4\times10^{11}\ \mathrm{cm}^{-2}\), an order of magnitude lower than previous limits toward TMC-1 [2101.05858].

These non-detections collectively indicate that GOTHAM is sensitive not just to molecular complexity per se, but to **chemical selectivity**: why some unsaturated carbonaceous species are abundant in TMC-1 while others, including heterocyclic analogues or more saturated products, are not.

## 7. Astrochemical interpretation, model tensions, and benchmark status

A recurring result across GOTHAM papers is that observations often outpace existing astrochemical models in absolute abundance, while some relative trends are reproduced more successfully.

For **propargyl cyanide**, NAUTILUS-v1.1 with the KIDA 2014 network and a dominant formation route

\[
\ce{CN + CH_3CCH \rightarrow HCCCH_2CN + H}
\]

reproduces the observed abundance within a factor of a few, suggesting that the main formation picture for this unsaturated precursor is approximately correct [2008.12349].

By contrast, **indene** presents a much stronger challenge. The indene paper states that NAUTILUS models underpredict the observed abundance by **several orders of magnitude**, with modeled column densities of order \(10^8\)–\(10^9\ \mathrm{cm}^{-2}\) against an observed value near \(10^{13}\ \mathrm{cm}^{-2}\) [2104.15117]. Likewise, the \(E\)-1-cyano-1,3-butadiene study reports that even after expanding the reaction network and increasing modeled 1,3-butadiene abundances by about **two orders of magnitude** relative to previous GOTHAM models, modeled abundances of aromatic species remain underestimated by **1–4 orders of magnitude** [2303.12221].

The methyl-chain study reports that **methylcyanopolyynes** decrease **log-linearly** with chain length and are broadly matched by a NAUTILUS v1.1 three-phase gas-grain model, whereas the sharp observed drop in **methylpolyynes** at \(\mathrm{CH_3C_8H}\) is not reproduced [2110.13152]. That same paper interprets differing \(A/E\) nuclear-spin isomer ratios in methylcyanopolyynes and methylpolyynes as evidence for **separate interstellar formation pathways** [2110.13152].

The \(\mathrm{C_{10}H^-}\) analysis exposes a different class of model discrepancy. The measured ratio

\[
\frac{N(\mathrm{C_{10}H^-})}{N(\mathrm{C_{10}H})} \approx 2.0^{+5.9}_{-1.6}
\]

is described as the highest measured between an anion and neutral species to date and as being at odds with current theories of interstellar anion chemistry [2301.07760]. That work further compares the observed abundances to a gas/grain chemical model and to a machine learning analysis, finding that the machine learning approach matches the detected anion abundances much better than the gas/grain model, which is presented as evidence that current understanding of molecular anion formation chemistry is still highly uncertain [2301.07760].

The c-\(\mathrm{C_3H_2S}\) detection adds an isomeric-chemistry perspective. Its abundance relative to the previously detected CH\(_2\)CCS and HCCCHS follows the **relative dipole principle (RDP)**, a kinetic rule-of-thumb according to which, all other chemistry among a family of isomers being the same, the member with the smallest dipole should be the most abundant [2501.06343]. The paper states that the RDP now holds for both S-bearing and O-bearing counterparts observed in TMC-1, although CH\(_2\)CCO continues to elude detection [2501.06343]. This suggests that kinetic control, including ion-molecule destruction scaling with dipole moment, can dominate over thermodynamic ordering in at least some isomer families.

The 2025 inventory paper formalizes GOTHAM’s benchmark role by presenting a reference set of gas-phase molecular abundances for interstellar clouds [2509.06256]. It states that the inventory accounts for approximately **0.12%** of gas-phase carbon, **0.0015%** of oxygen, **5%** of sulfur, and **~80%** of detected gas-phase nitrogen, while confirming up to **four velocity components** and identifying some molecules—such as CCS, \(\mathrm{H_2CO}\), and \(c\)-\(\mathrm{C_3H_2}\)—as requiring non-LTE or more complex models [2509.06256]. In that sense, GOTHAM functions not only as a discovery program but as an empirical reference against which chemical networks, excitation treatments, isotopic fractionation schemes, and catalog accuracy can be tested.

Taken together, the survey establishes TMC-1 as a chemically rich environment dominated in the detected gas phase by **unsaturated hydrocarbons**, nitriles, and related chain species, in marked contrast to the oxygen-rich organics characteristic of sublimated ices around protostars [2509.06256]. The broader implication is that cold dark-cloud chemistry supports an unexpectedly extensive reservoir of aromatic and pre-aromatic molecules, yet that the mechanistic basis of this reservoir remains only partially captured by current gas-grain models.

Source: https://www.emergentmind.com/topics/gotham-large-program-survey