QUIJOTE Q-Band Line Survey in TMC-1

• QUIJOTE line survey is a comprehensive Q-band molecular survey that maps weak rotational transitions in the cold dark cloud TMC-1 using ultralow noise observations.
• It employs advanced techniques including frequency-switching and sophisticated data reduction methods to detect rare molecules and characterize their spectral signatures.
• The survey benchmarks astrochemical models by quantifying new hydrocarbons, PAHs, ions, and isotopologues, providing crucial insights into interstellar chemistry.

The QUIJOTE line survey is a highly sensitive, broadband molecular-line survey conducted in the Q-band (31.0–50.3 GHz) using the Yebes 40 m radio telescope, primarily targeting the cold dark cloud TMC-1. It has been a key driver in identifying an unprecedented inventory of interstellar molecules—spanning hydrocarbon chains, rings, radicals, isotopologue systems, ions, and early polycyclic aromatic hydrocarbons (PAHs)—and in benchmarking the physical and chemical network of cold molecular clouds. The survey achieves deep integration, ultralow rms noise (down to 0.06–0.2 mK), and high spectral resolution (38.15 kHz), enabling detection of species with low dipole moments and weak line intensities. The following sections detail its technical parameters, observing modes, data analysis methodologies, key discoveries, and astrochemical impact.

1. Technical Setup and Survey Parameters

QUIJOTE is based on the Yebes 40 m radio telescope, equipped with a dual-polarization, cryogenically cooled HEMT receiver covering the Q-band (31.0–50.3 GHz) in a single tuning. The signal is split into eight 2.5 GHz sub-bands per polarization, each analyzed by FFT spectrometers, providing 38.15 kHz channel spacing (corresponding to ~0.25 km s⁻¹ at 45 GHz). Main-beam efficiency η_MB varies from ~0.60 at 32 GHz to ~0.43 at 50 GHz, with a half-power beamwidth shrinking from 54.4″ at 32.4 GHz to 36.4″ at 48.4 GHz. System temperatures improved during the survey, reaching 16 K at 32 GHz and 30 K at 50 GHz in 2021. Achievable rms noise is ~0.08 mK at 32 GHz and 0.25 mK at 49.5 GHz after typical on-source integrations of 500–1500 h.

Observations are targeted at the cyanopolyyne peak of TMC-1 at (α, δ)_J2000 = (04ʰ 41ᵐ 41.9ˢ, +25° 41′ 27.0″), with a uniform source size of ∼80″ adopted in analyses for beam-dilution corrections. The survey is about 50–100 times more sensitive than previous Q-band surveys of TMC-1, uncovering sub-milliKelvin, weak rotational signatures of highly abundant, rare, and previously undetectable species (Fuentetaja et al., 2023, Fuentetaja et al., 19 Jan 2026, Fuentetaja et al., 2024, Cernicharo et al., 2024, Fuentetaja et al., 17 Sep 2025).

Parameter	Value/Characteristic	Reference
Telescope	Yebes 40 m, IGN (Spain)	(Fuentetaja et al., 2023, Fuentetaja et al., 19 Jan 2026)
Frequency Coverage	31.0–50.3 GHz (Q-band)	(Fuentetaja et al., 2023, Cernicharo et al., 2024)
Spectral Resolution	38.15 kHz (~0.25 km s⁻¹)	(Fuentetaja et al., 2023)
rms Noise	0.06–0.25 mK	(Fuentetaja et al., 19 Jan 2026, Fuentetaja et al., 2024)
Beam Size (HPBW)	54.4″ (32.4 GHz)→36.4″ (48.4 GHz)	(Fuentetaja et al., 2023)
On-source Integration	500–1500 h	(Fuentetaja et al., 19 Jan 2026)

2. Observing Strategy and Data Reduction

Data are collected in frequency-switching mode with throws of 8 or 10 MHz and multiple LO tunings to suppress and blank negative artifacts and avoid spurious signals. Calibration is performed with hot/cold load techniques and the ATM atmospheric model, maintaining an overall absolute intensity uncertainty of ≤10%. The GILDAS/CLASS and MADEX software frameworks serve for reduction, baseline subtraction, folding, blanking of artifacts, and line profile fitting. Integrated spectra can achieve rms noise well below 0.2 mK, approaching the confusion limit for TMC-1 (Fuentetaja et al., 19 Jan 2026, Cernicharo et al., 2023).

Identification of molecular lines involves searching for harmonic or near-harmonic series of transitions matching predicted frequencies from laboratory spectroscopy and quantum-chemical calculations. Gaussian fitting of each candidate requires a ≥3σ detection and linewidths consistent with the subsonic turbulence of TMC-1 (typical Δv ≃ 0.6–1.1 km s⁻¹).

3. Line Identification and Quantitative Analysis

Identifications are consolidated using spectroscopic models (Hamiltonians of increasing complexity for asymmetric rotors, including RAM36, Watson A-reduction) and line catalogs (MADEX, CDMS, JPL). Transition frequencies, dipole moments, Einstein A coefficients, and partition functions are taken from laboratory measurements or ab initio calculations. Examples include methyl ketene (CH₃CHCO, μ_a=1.65 D, μ_b=0.33 D), cyclopentindene (c-C₁₁H₈, μ_a=0.60 D, μ_b=0.80 D), and 1,4-pentadiyne (HCCCH₂CCH, μ_a=0.67 D). Detection of low-dipole PAH derivatives is enabled by the survey's sensitivity (Fuentetaja et al., 19 Jan 2026, Cernicharo et al., 2024, Fuentetaja et al., 2024).

The rotational diagram (Boltzmann plot) method is the standard for deriving rotational temperatures and column densities, assuming local thermodynamic equilibrium (LTE) and optically thin emission. The population of the upper state $N_u$ follows:

$$N_u = \frac{8\pi k_B \nu^2}{hc^3 A_{ul}} \int T_A^*\,dv$$

The total column density $N$ and rotational temperature $T_{\rm rot}$ are extracted by fitting:

$$\ln\frac{N_u}{g_u} = \ln\frac{N}{Q_{\rm rot}(T_{\rm rot})} - \frac{E_u}{k_B T_{\rm rot}}$$

Non-LTE effects (e.g., multi-level population structure in radicals like C₆H) are analyzed where collisional-rate coefficients are available (Cernicharo et al., 2023).

4. Chemical Inventory and Key Discoveries

QUIJOTE has enabled the secure identification of a broad spectrum of molecular species in TMC-1. Notable discoveries and quantifications include:

• New hydrocarbons and PAHs: Methyl ketene (CH₃CHCO), 1,4-pentadiyne (HCCCH₂CCH), cyclopentindene (c-C₁₁H₈), fulvenallene (c-C₅H₄CCH₂), benzyne (o-C₆H₄), and cyano derivatives of acenaphthylene (1- and 5-C₁₂H₇CN) (Fuentetaja et al., 2023, Fuentetaja et al., 2024, Fuentetaja et al., 19 Jan 2026, Cernicharo et al., 2021, Cernicharo et al., 2022, Cernicharo et al., 2024).
• Exotic (iso)topologues/isomers: Thiofulminic acid (HCNS), rare S and C isotopologues of the CₙS family, deuterated methyldiacetylene (CH₂DC₄H), and vibrationally excited C₆H (Cernicharo et al., 2024, Fuentetaja et al., 17 Sep 2025, Cabezas et al., 2021, Cernicharo et al., 2023).
• Radicals and ions: HCCCO, C₅O, C₅H⁺, C₃H⁺, NC₃S, HC₃S, and families of O- and N-bearing carbon chains (Cernicharo et al., 2021, Cernicharo et al., 2022, Cernicharo et al., 2024).

Column densities for these species often reach as low as 10¹⁰–10¹² cm⁻². Fractional abundances, derived with $N({\rm H}_2)=10^{22}$ cm⁻², span from $10^{-9}$ (major PAHs/chain radicals) down to $10^{-12}$ (rare cations, isotopologues).

5. Astrochemical Modeling and Isomerism

Astrochemical interpretation relies on pseudo–time–dependent gas-phase networks, typically based on UMIST RATE12 or later, updated with new laboratory or theoretical reaction rates. Key insights from modeling with QUIJOTE constraints include:

• Accurate reproduction of small chain and ring abundances: Most detected cyclic and acyclic species’ abundances are matched within factors of a few by models with updated reaction rates and proper branching ratios (Fuentetaja et al., 2024, Fuentetaja et al., 2023).
• Discrimination among isomers: For C₃H₄O (e.g., methyl ketene, trans- and cis-propenal), C₅H₄, and C₁₂H₈ derivatives, computed column densities and relative energies place constraints on formation mechanisms and kinetic control over thermodynamic equilibria (Fuentetaja et al., 2023, Fuentetaja et al., 2024, Cernicharo et al., 2024).
• Importance of barrierless, neutral–neutral and ion–neutral routes: Many cyclic species (e.g., o-benzyne, cyclopentadiene derivatives) and complex PAHs likely form via rapid radical–radical or ion–neutral reactions at 10 K, evidenced by the necessity of such pathways to reconcile models with observed abundances (Cernicharo et al., 2022).
• Anomalies in isotopic and D/H fractionation: Position-dependent $^{13}$C fractionation in CCS/CCCS is well explained by exothermic H-exchange reactions; deuterium enrichment in methyldiacetylene families is matched only when atomic D–exchange channels are included in network models (Fuentetaja et al., 17 Sep 2025, Cabezas et al., 2021).

6. Benchmarks for Sulphur, Oxygen, and Nitrogen Chemistry

The survey has yielded the most complete inventory of S-bearing CₙS, HCS⁺, and related ions (including C¹³C³⁴S, CC³³S, CCC³³S, HC³³S⁺, and HCC³⁴S⁺), with measured isotopic ratios matching the solar system’s S and C isotopic values, except for known anomalies in multipositional carbon chains due to fractionation (Fuentetaja et al., 17 Sep 2025). Similar completeness is achieved in the O-bearing chains (e.g., HCCCO, HC₅O, C₅O), as well as in families of cyano and nitrile derivatives (Cernicharo et al., 2021, Cernicharo et al., 2022). This breadth enables constraints on elemental depletion (S, N, O), benchmarking models of grain-surface vs. gas-phase chemistry, and the identification of missing or underestimated reaction routes.

7. Implications, Limitations, and Future Directions

QUIJOTE’s legacy is multifunctional: (i) providing comprehensive, high-precision molecular inventories for cold dark clouds; (ii) yielding stringent benchmarks for statistical and kinetic chemical models over a wide dynamic range in physical and chemical complexity; and (iii) facilitating the derivation or refinement of rotational constants, aiding future surveys. Non-detections (e.g., certain high-energy isomers, longer S-chains, weak isotopologues) deliver critical upper limits that delimit feasible astrochemical parameter space (Cernicharo et al., 2024).

Challenges remain in accelerating laboratory reaction-rate studies (especially for neutral–neutral and certain DR channels at 10 K), completing collisional rate coefficient databases for excitation modeling, and expanding chemical networks to accommodate the efficient buildup of PAHs and multiply substituted ring systems consistent with QUIJOTE observations (Fuentetaja et al., 19 Jan 2026, Cernicharo et al., 2024). The detection statistics and derived abundance ratios strongly favor bottom-up, kinetic-dominated formation of rings and complex molecules, supporting a scenario in which stepwise addition of small radicals (C₂H, CN) assembles even the earliest aromatic and PAH structures in the cold ISM (Cernicharo et al., 2024, Fuentetaja et al., 19 Jan 2026).

• (Fuentetaja et al., 2023) Discovery of CH3CHCO in TMC-1 with the QUIJOTE line survey
• (Fuentetaja et al., 2024) Discovery of HCCCH2CCH in TMC-1 with the QUIJOTE line survey
• (Fuentetaja et al., 19 Jan 2026) Discovery of 1H-cyclopent[cd]indene (c-C11H8) in TMC-1 with the QUIJOTE line survey
• (Fuentetaja et al., 17 Sep 2025) Analysis of the isotopologues of CS, CCS, CCCS, HCS+, HCCS+, and H2CS in TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2023) Detection of vibrationally excited C6H in the cold prestellar core TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2021) Discovery of benzyne, o-C6H4, in TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2021) Discovery of HCCCO and C5O in TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2022) Discovery of fulvenallene in TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2022) Discovery of five cyano derivatives of propene with the QUIJOTE line survey
• (Cernicharo et al., 2024) Discovery of thiofulminic acid with the QUIJOTE line survey: A study of the isomers of HNCS and HNCO in TMC-1
• (Cabezas et al., 2021) New deuterated species in TMC-1: Detection of CH2DC4H with the QUIJOTE line survey
• (Cernicharo et al., 2022) Discovery of C5H+ and detection of C3H+ in TMC-1 with the QUIJOTE line survey
• (Fuentetaja et al., 2022) Discovery of CH2CCHC4H and a rigorous detection of CH2CCHC3N in TMC-1 with the QUIJOTE line survey
• (Cernicharo et al., 2024) More sulphur in TMC-1: Discovery of the NC3S and HC3S radicals with the QUIJOTE line survey
• (Cernicharo et al., 2024) Discovery of two cyano derivatives of acenaphthylene (C12H8) in TMC-1 with the QUIJOTE line survey