# Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the Submillimeter Array

A. Zernickel, P. Schilke, A. Schmiedeke, D. C. Lis, C. L. Brogan, C. Ceccarelli, C. Comito, M. Emprechtinger, T. R. Hunter, T. Möller

► To cite this version:

A. Zernickel, P. Schilke, A. Schmiedeke, D. C. Lis, C. L. Brogan, et al., Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the Submillimeter Array. Astronomy and Astrophysics - A&A, 2012, 546, 10.1051/0004-6361/201219803. insu-03612405

HAL Id: insu-03612405

https://insu.hal.science/insu-03612405

Submitted on 18 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the Submillimeter Array★★★

A. Zernickel<fnref n="1" />, P. Schilke<fnref n="1" />, A. Schmiedeke<fnref n="1" />, D. C. Lis<fnref n="3" />, C. L. Brogan<fnref n="4" />, C. Ceccarelli<fnref n="5" />,
C. Comito<fnref n="2" />, M. Emprechtinger<fnref n="3" />, T. R. Hunter<fnref n="4" />, and T. Möller<fnref n="1" />

Received 12 June 2012 / Accepted 25 August 2012

## ABSTRACT

**Aims.** We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions.

**Methods.** In the framework of the *Herschel* guaranteed time key program CHESS (Chemical *Herschel* Surveys of Star forming regions), NGC 6334I is investigated by using the Heterodyne Instrument for the Far-Infrared (HIFI) aboard the *Herschel* Space Observatory. A spectral line survey is carried out in the frequency range 480–1907 GHz, and further auxiliary interferometric data from the Submillimeter Array (SMA) in the 230 GHz band provide spatial information for disentangling the different physical components contributing to the HIFI spectrum. The spectral lines in the processed *Herschel* data are identified with the aid of former surveys and spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum, assuming local thermal equilibrium, and best fit parameters are derived using a model optimization package.

**Results.** A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High-energy levels ($E_u$ > 1000 K) of the dominant emitter methanol and vibrationally excited HCN ($v_2$ = 1) are detected. The number of unidentified lines remains low with 75, or <2% of the lines detected. The modeling suggests that several spectral features need two or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible in the SiO and H$_2$O emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules identified in SMA1 and O-bearing molecules in SMA2.

**Conclusions.** Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances provide further constraints for astrochemical models.

**Key words.** submillimeter: ISM – line: identification – astrochemistry – ISM: molecules – ISM: abundances – stars: formation

## 1. Introduction

The astrochemical composition of the interstellar gas depends on the physical environment and the evolutionary state. In the cold dense regions molecules are unsaturated, while in contrast saturated molecules are much more abundant in hot cores (Herbst & van Dishoeck 2009). Unsaturated molecules are carbon rich with long carbon chains, like radicals as $C_nH$ ($n$ = 2–8) or cyanopolyynes $HC_nN$ ($n$ = 3, 5, 7, 9, 11), whereas saturated organic molecules are hydrogen rich with single bonds of carbon like $CH_3OCH_3$. Molecular ions are predominantly unsaturated in the cold phase, because most hydrogenation reactions are endothermic or hindered by potential barriers (Herbst 1989). Ices on dust grains begin to form in the cold phase and are observed by broad absorption bands at mid-infrared wavelengths. Sequential hydrogenation of oxygen or nitrogen forms

water or ammonia, respectively. Embedded in the water-ices are relevant amounts of CO, and through further hydrogenation $CO \rightarrow HCO \rightarrow H_2CO \rightarrow H_2COH \rightarrow CH_3OH$. Time-dependent astrochemical models suggest that abundances of certain molecules, such as $H_2S$ and $SO_2$, might be used as an indicator of the evolutionary phase (Nomura & Millar 2004; Herpin et al. 2009), and greater chemical complexity can be expected for longer evolution timescales (Garrod et al. 2008). Approximately 165 molecules have been detected in the ISM or circumstellar shells<fnref n="1" />.

Van Dishoeck (2001) emphasized the importance of spectral line surveys for the research of star-forming regions. Line surveys give the possibility to obtain a census of all atoms and molecules and give insights into their thermal excitation conditions and dynamics by studying line intensities and profiles, which allows one to separate different physical components. The discovery of new species or new excited lines from known molecules is possible. Furthermore, the wide range of energy

<footnote n="1">¹ I. Physikalisches Institut der Universität zu Köln, Zülpicher Straße 77, 50937 Köln, Germany
e-mail: zernickel@phl.uni-koeln.de</footnote>

<footnote n="2">² Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany</footnote>

<footnote n="3">³ California Institute of Technology, MC 301-17, Pasadena, CA 91125 USA</footnote>

<footnote n="4">⁴ NRAO, 520 Edgemont Road, Charlottesville, VA 22903, USA</footnote>

<footnote n="5">⁵ IPAG, Université Joseph Fourier, CNRS, BP 53, 38041 Grenoble Cedex 9, France</footnote>

<footnote>* Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
** Appendices are available in electronic form at http://www.aanda.org</footnote>

<footnote n="1">¹ A current list can be found at http://www.astro.uni-koeln.de/cdms/molecules</footnote>
<table><caption>**Table 1.** Mass and temperature estimates for NGC 6334I.</caption><thead><tr><th>Source</th><th><i>M</i>(<i>M</i><sub>☉</sub>)</th><th><i>T</i><sub>d</sub> (K)</th><th>β</th><th><i>n</i> (cm<sup>-3</sup>)</th><th><i>L</i><sub>bol</sub> (<i>L</i><sub>☉</sub>)</th></tr></thead><tbody><tr><td>Total<sup>a</sup></td><td>200 ± 100</td><td>100</td><td>1.7 ± 0.3</td><td>1.2 × 10<sup>7</sup></td><td>2.6 × 10<sup>5</sup></td></tr><tr><td>SMA1<sup>b</sup></td><td>17 ± 50%</td><td>100</td><td></td><td></td><td></td></tr><tr><td>SMA2<sup>b</sup></td><td>11 ± 50%</td><td>100</td><td></td><td></td><td></td></tr><tr><td>SMA3<sup>b</sup></td><td>33 ± 50%</td><td>60</td><td></td><td></td><td></td></tr><tr><td>SMA4<sup>b</sup></td><td>36 ± 50%</td><td>33</td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Values taken from (a) Sandell (2000) and (b) Hunter et al. (2006). *M* is the mass, *T*<sub>d</sub> the average dust temperature, *β* the dust emissivity index, *n* the hydrogen density and *L*<sub>bol</sub> the bolometric luminosity.

levels and eventually opacities of one species allows one a better determination of the temperature and abundance. The uniform recording with the same spectrometer in one observation reduces the calibration uncertainties.

The giant molecular cloud NGC 6334, also known as the “Cat’s Paw Nebula”, lies in the Carina-Sagittarius Arm of the Milky Way at a distance of (1.7 ± 0.3) kpc (Russeil et al. 2012). Studies at infrared wavelengths revealed sites of massive star formation (Persi et al. 2005). NGC 6334I was first explored by Emerson et al. (1973) as the brightest infrared source in the northern part of the cloud, containing an OH maser and an ultracompact H II region. According to Emerson et al. (1973), the pumping of the OH maser requires an energy source of a massive collapsing protostar (*M* > 30 *M*<sub>☉</sub>). H II regions are associated with high-mass stars (spectral classes O and B) where the UV radiation is strong enough to ionize hydrogen. Russeil et al. (2010) studied the densest cores in NGC 6334 in molecular line and dust continuum emission. For NGC 6334I, they derive a total mass of *M* = 1206 *M*<sub>☉</sub> from the continuum flux density at 1.2 mm, an average density of *n*<sub>H₂</sub> = 1.2 × 10<sup>6</sup> cm<sup>-3</sup> and an upper limit for the virial mass of *M*<sub>vir</sub> = 525 *M*<sub>☉</sub>. The total mass estimate is relatively high, because a low value of the dust temperature of only 20 K was assumed, as compared to the value of *T* = 100 K derived by (Sandell 2000), which leads to a factor of six difference of the mass estimate. While an average dust temperature of *T* = 25 K is appropriate for most cores in the NGC 6334 complex (Matthews et al. 2008), the dust temperature is known to be much higher in NGC 6334I.

Interferometric studies revealed that NGC 6334I consists of several cores, labeled SMA1 to SMA4 (Hunter et al. 2006). The separation distance between SMA1 and SMA2 is about 3.5" or 6000 AU. SMA1 and SMA2 show rich spectra of molecular transitions and SMA3 coincides with the H II region which exhibits free-free and dust emission and is excited by the near-infrared source IRS 1E (De Buizer et al. 2002). SMA4 only shows dust emission but no line emission, which is still unexplained. Table 1 summarizes characteristic values for NGC 6334I, derived from the dust continuum emission at different wavelengths, which classify the source as a low-luminosity high-mass hot core. The overall size of NGC 6334I is 10" × 8", where 1" equals to 1700 AU. SMA1 and SMA2 are so-called hot molecular cores, where certain attributes apply according to Cesaroni (2005): a high temperature (*T* ≥ 100 K), sizes smaller than 0.1 pc, large masses (10–1000 *M*<sub>☉</sub>) and luminosities exceeding 10<sup>4</sup> *L*<sub>☉</sub>. Furthermore, they are often associated with water masers and ultracompact H II regions. Hot cores are heated by high-mass protostars which are embedded in a dusty gas envelope. Their rich molecular spectra originate from the evaporation of ice mantles covering the dust grains, which release at *T* ~ 100 K molecules to the gas phase. Molecules released

from grain mantles then begin to drive a fast, high-temperature gas-phase chemistry, forming complex organic species like acetonitrile, beginning from precursors like methanol. Their rotational, vibrational and radiative excitation leads to a characteristic emission spectrum. This explains why the spectrum is poorer in earlier stages like the colder protostellar core NGC 6334I(N) (Brogan et al. 2009; Walsh et al. 2010), 2' or 1 pc north of NGC 6334I.

Multiple spectral line surveys have been carried out for NGC 6334I: McCutcheon et al. (2000) 334–348 GHz, Schilke et al. (2006) 459–461 GHz and 817–819 GHz, Thorwirth et al. (2003) 88–115 GHz and 218–267 GHz, Walsh et al. (2010) 84–116 GHz and the latest from Kalinina et al. (2010) in the range 81–242 GHz. Besides, several maser lines have been observed from NH<sub>3</sub> (Beuther et al. 2007), CH<sub>3</sub>OH (Walsh et al. 1998) and H<sub>2</sub>O (Migenes et al. 1999) as well as bipolar outflows in CO (Leurini et al. 2006; Qiu et al. 2011) and SiO (Gibb et al. 2007). This source is part of the CHESS program (Chemical HERSchel Surveys of Star forming regions), a key program of the Herschel Space Observatory (Ceccarelli et al. 2010). The aim is to study eight different sites and phases of star formation varying in mass, luminosity, evolutionary state, astrochemical composition etc. by conducting and comparing their spectra. In this article, the entire spectrum (500–1900 GHz) of NGC 6334I is analyzed which was observed by Herschel in 2010. Previously, some small sections of the spectrum were studied to find transitions of new molecules not observed before in the interstellar medium (ISM): Oxidaniumyl (H<sub>2</sub>O<sup>+</sup>) by Ossenkopf et al. (2010) and chloronium (H<sub>2</sub>Cl<sup>+</sup>) by Lis et al. (2010). A detailed analysis of water has been made by Emprechtinger et al. (2010) and of methylidyne (CH) by van der Wiel et al. (2010).

One difficulty is that if a large beam encompasses a complex, unresolved source structure, the spectrum can become complex due to different overlapping components. The spatial information from the SMA is complementary. The Herschel beam is larger than NGC 6334I so that the source is spatially unresolved, while the SMA beam can resolve the individual cores. The SMA maps thus reveal the source morphology and give insights into the distribution of various molecules. This helps to decompose the contributions from different components present in the HIFI spectrum and to differentiate between the two hot cores SMA1 and SMA2, the extended emission and outflows, see Fig. 1.

## 2. Observations

### 2.1. Observations with Herschel

Observations were conducted with the high spectral resolution receiver system HIFI (Heterodyne Instrument for the Far-Infrared; de Graauw et al. 2010) aboard the Herschel Space Observatory (Pilbratt et al. 2010). HIFI is divided into seven bands with different receivers which overlap at the band edges. Each band is separated again into two parts a and b using different polarizations and local oscillator chains. The full frequency range is from 480–1907 GHz, with a gap between 1280 and 1410 GHz. The first five bands employ SIS (superconductor-insulator-superconductor) mixers with an intermediate frequency (IF) bandwidth of 4 GHz. Bands 6 and 7 are based on Hot-Electron Bolometer (HEB) mixers and have an IF bandwidth of 2.4 GHz. The HPBW of the 3.5 m Herschel telescope varies with wavelength from 41" to 12". Table 2 gives a log of the HIFI observations of NGC 6334I and Fig. 2 shows the whole calibrated spectrum. The continuum level increases
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Integrated SMA line intensity map over all species. Overlaid in black contours is the continuum emission. The contour levels are from 90% to 10% of the maximum in steps of 10%. The large and small white circles represent the FWHM *Herschel* beam size for band 1 and band 7, respectively.</figcaption></figure>

<table><caption>**Table 2.** Overview over all HIFI bands for the observation of NGC 6334I.</caption><thead><tr><th>Band</th><th>Freq. range (GHz)</th><th>Beam size (")</th><th>Date</th><th>Obs. time (h)</th><th>rms (K)</th></tr></thead><tbody><tr><td>1a</td><td>480–560</td><td rowspan="2">41</td><td>28.02.2010</td><td>3.6</td><td>0.015</td></tr><tr><td>1b</td><td>554–636</td><td>02.03.2010</td><td>2.2</td><td>0.03</td></tr><tr><td>2a</td><td>626–726</td><td rowspan="2">29</td><td>09.10.2010</td><td>3.7</td><td>0.04</td></tr><tr><td>2b</td><td>714–801</td><td>19.03.2010</td><td>3.3</td><td>0.05</td></tr><tr><td>3a</td><td>799–860</td><td rowspan="2">25</td><td>09.10.2010</td><td>1.8</td><td>0.05</td></tr><tr><td>3b</td><td>858–961</td><td>19.03.2010</td><td>3.9</td><td>0.07</td></tr><tr><td>4a</td><td>949–1061</td><td rowspan="2">21</td><td>09.10.2010</td><td>4.8</td><td>0.09</td></tr><tr><td>4b</td><td>1052–1122</td><td>03.03.2010</td><td>2.9</td><td>0.10</td></tr><tr><td>5a</td><td>1108–1244</td><td rowspan="2">18.5</td><td>13.10.2010</td><td>4.2</td><td>0.37</td></tr><tr><td>5b</td><td>1227–1280</td><td>12.03.2011</td><td>1.7</td><td>0.37</td></tr><tr><td>6a</td><td>1410–1575</td><td rowspan="2">14.5</td><td>14.10.2010</td><td>0.9</td><td rowspan="2">1.20</td></tr><tr><td>6b</td><td>1575–1700</td><td>14.10.2010</td><td>3.6</td></tr><tr><td>7a</td><td>1703–1799</td><td rowspan="2">12.5</td><td>14.10.2010</td><td>2.6</td><td>1.27</td></tr><tr><td>7b</td><td>1788–1907</td><td>13.10.2010</td><td>3.2</td><td>1.27</td></tr></tbody></table>

continuously with frequency, except for band 5b, where a small offset is seen. The line emission in band 1a contributes only 3% to the total flux because the HPBW at 520 GHz is four times bigger than the source size so that in addition the continuum from the surrounding area is observed. The emission lines with the highest peak intensity, $T > 15$ K, are those of CO and $^{13}$CO. The telescope was pointed toward the coordinates $\alpha(\text{J2000}): 17^{\prime\prime}20^{\prime\prime}53.32^{\prime\prime}$, $\delta(\text{J2000}): -35^{\circ}46'58.5''$, between SMA1 and SMA2.

The observations were carried out mainly in two sessions in March and October 2010, in dual beam switch mode with an internal chopper mirror switching between the source and two reference positions 3' away from the target to provide flat baselines and calibrate the continuum. The total observing time was ~42 h. The wide band spectrometer provided a spectral resolution of 1.1 MHz, resulting in a velocity resolution of 0.7–0.2 km s$^{-1}$. The double sideband spectra were processed with HIPE (Herschel Interactive Processing Environment; Ott 2010) version 5, deconvolved to single-sideband spectra following the algorithms described in (Comito & Schilke 2002), and

exported to CLASS (see Sect. 3.1) for subsequent analysis. All spectra are corrected for the beam efficiency and the intensity scale is the main-beam brightness temperature in K. The spectra presented here are equally weighted averages of the horizontal and vertical polarizations, in order to reduce the noise. Band 4a has some sinusoidal baseline features from standing waves. For band 6a only a small portion around a $H\alpha$ transition was observed, but the data are corrupted and not discussed here. The band coverage is unbiased, meaning that all frequency ranges are included with the band edges where the noise increases. Different sideband gains can lead to discrepancies between line intensities in the two bands in the overlap region, see Fig. B.1g for an example in the $J = 6-5$ transition of $N_2H^+$ at 558.966 GHz. The discrepancies amount to 10–20% in the bands 1–3 and become negligible in the higher bands. In the further analysis, this uncertainty can be avoided and has marginal effects when several lines of a given molecule are present.

## 2.2. Auxiliary observations from the Submillimeter Array

Supporting SMA data at 230 GHz are kindly provided by Brogan and Hunter. The continuum map is identical to that published in Hunter et al. (2006), whereas the spectral image cubes are unpublished. Detailed information about the observation and calibration can be found therein. The data were calibrated again in 2007 in Miriad (described in Brogan et al. 2009). The separation of the line and continuum emission was performed in the uv plane. For the purpose of modeling, the interferometric line and continuum emission were recombined in CLASS for the portion of the spectrum containing the foreground CN absorption feature. The observed frequency range covers two 2 GHz wide bands, 216.6–218.6 GHz in the lower sideband and 226.6–228.6 GHz in the upper sideband, with a frequency resolution of ~0.8 MHz. The map covers a 64" × 64" field and the synthesized beam is 2.2"×1.3" in the USB. The line-to-continuum ratio is 12% for SMA1 and 21% for SMA2 in the given bands.

# 3. Analysis: line identification and modeling

## 3.1. Simulating spectra with XCLASS

A spectral line is identified by comparing the center frequency with a measured frequency of a known chemical species (atom, molecule, ion) in the laboratory. Two large databases relevant for molecular astrophysics, from the radio to the far-infrared region, have been built up in the last decades: The Cologne Database for Molecular Spectroscopy (CDMS; Müller et al. 2005) and the Jet Propulsion Laboratory (JPL) catalog (Pickett et al. 1998). If the exact source velocity is unknown or several candidates are present for the identification of a spectral line, secondary criteria such as the upper energy level or the Einstein A coefficient must be taken into account. Although many lines in the spectra of hot cores are blended, the identification in spectral surveys is more confident because more lines of a species are present, which can then be checked for consistency, i.e. all lines of a molecule which are expected to be there should be found with the expected line strength, line shape and velocity.

The main program used in this work is XCLASS, which was described in Comito et al. (2005) and applied to a ground-based Orion KL line survey. It is an extension of CLASS (Continuum and Line Analysis Single-dish Software), which is part of the GILDAS<fnref n="2" /> package. XCLASS includes entries from the CDMS,

<footnote n="2">² http://www.iram.fr/IRAMFR/GILDAS</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 2. Complete HIFI spectrum of NGC 6334I. B-bands are colored in blue and bands 5–7 are smoothed by averaging over 5–10 channels.</figcaption></figure>

JPL and a private catalog, which contain all necessary information: rest frequency $v_l$, integrated intensity, lower state energy $E_l$, upper state degeneracy $g_{up}$, quantum numbers and the partition function $Q$. With this database, it is possible to simulate a spectrum of a molecule, by providing parameters such as the excitation temperature and abundance, and compare the simulated spectrum with the observed one. Similar software tools developed in the last years are *Weeds* (Maret et al. 2011), a package for *CLASS*, and *CASSIS* (Caux et al. 2011). The line identification for NGC 6334I, which transfers to a molecule identification in XCLASS, is made easier because 20 detected molecules are reported in the line survey from Kalinina et al. (2010) at lower frequencies and in other surveys mentioned in Sect. 1.

The function which models a spectrum in XCLASS is a solution of the radiative transfer equation for an isothermal object in one dimension, called detection equation (Stahler & Palla 2005). Furthermore, the finite source size and dust attenuation are considered. The model function is:

$$T(v) = \sum_{m} \sum_{c} \eta(\theta_{m,c}) [J(T_{\text{ex}}^{m,c}) - J(T_{\text{bg}})] (1 - e^{-\tau(v)^{m,c}}) \cdot e^{-\tau_d}, \quad (1)$$

with the optical depth

$$\tau(v)^{m,c} = \sum_{l} \tau(v_l)^{m,c} \quad (2)$$

$$\tau(v_l)^{m,c} = \frac{c^3}{8\pi v^3} A_{ul} N_{\text{tot}}^{m,c} \frac{g_l e^{-E_l/kT_{\text{ex}}^{m,c}}}{Q(m, T_{\text{ex}}^{m,c})} (1 - e^{-hv/kT_{\text{ex}}^{m,c}}) \phi(v)^{m,c}, \quad (3)$$

the dust optical depth

$$\tau_d = N_H \cdot \kappa_{1.3\,\text{mm}} \cdot \left( \frac{v}{230\,\text{GHz}} \right)^\beta \cdot 3.3 \times 10^{-26}, \quad (4)$$

the Gaussian line profile

$$\phi(v)^{m,c} = \frac{2 \sqrt{\ln 2}}{\sqrt{\pi} \Delta v^{m,c}} \exp\left(-\frac{(v - (v_l + v_{\text{LSR}}^{m,c}))^2}{2\sigma^2}\right) \quad (5)$$

and the beam filling factor

$$\eta(\theta_{m,c}) = \frac{\theta_{m,c}^2}{\theta_{m,c}^2 + \theta_t^2}. \quad (6)$$

The radiation temperature $J(T)$ is defined as $J(T) = \frac{hv}{k} \frac{1}{e^{hv/kT} - 1}$. The sum goes over the indices $m$ for molecule, $c$ for component and $l$ for spectral line transition. $T_{bg}$ is the background continuum temperature and $\theta_t$ the telescope beam size. The line profile function is normalized so that $\int_0^\infty dv \phi(v) = 1$. The factor $e^{-\tau_d}$ takes the dust attenuation into account where the dust mass opacity $\kappa$ is taken from Ossenkopf & Henning (1994). For NGC 6334I, we adopt a value of $\beta = 2$ and $\kappa_{1.3\,\text{mm}} = 0.42\,\text{cm}^2\,\text{g}^{-1}$, which is between the values for non-coagulated dust grains with thin ice mantles ($\kappa_{1.3\,\text{mm}} = 0.51\,\text{cm}^2\,\text{g}^{-1}$) and without ones ($\kappa_{1.3\,\text{mm}} = 0.31\,\text{cm}^2\,\text{g}^{-1}$). The total hydrogen column density $N_H = N(H)+2N(H_2)$ is $3 \times 10^{24}\,\text{cm}^{-2}$, which is derived from the average hydrogen density given in Table 1 and a source size of $10''$ and is in agreement with values from previous articles (van der Wiel et al. 2010; Lis et al. 2010) considering the uncertainties.

Overall, there are five input parameters for one component of a molecule:

$\theta_{m,c}$ source size in arcsec

$T_{\text{ex}}$ excitation temperature in K

$N_{\text{tot}}$ total column density in cm$^{-2}$

$\Delta v$ velocity width (FWHM) in km s$^{-1}$

$v_{\text{LSR}}$ source velocity in km s$^{-1}$.

Components can be identified as spatially distinct sources within the beam, such as clumps, hot dense cores, colder envelopes or outflows, and can usually be distinguished by different radial velocities. They do not interact with each other radiatively, but are simply superposed in the model. However, for absorption lines the emission from other components is considered first by calculating the background emission spectrum, which is then used as a new “continuum” for absorption lines. By fitting all species
and their components at once, line blending and optical depth effects are taken into account. Because of the high densities in the hot cores ($n_{H_2} > 10^7$ cm$^{-3}$), it is assumed that most molecules are thermalized at the gas temperature, so that the assumption of LTE is reasonable. With the LTE condition, every population of a level i is known and given by the Boltzmann distribution $n_i = \frac{N}{Q}g_i e^{-E_i/(kT)}$. The modeling of a spectrum speeds up line identification enormously, since several hundred lines of a molecule like methanol are predicted at once in a HIFI band.

The modeling can be done together with isotopologues (and higher vibrational states) of a molecule assuming a conversion factor stored in a separate database. All parameters are expected to be the same except the abundance. For carbon, Milam et al. (2005) derive an average gradient of $^{12}$C/$^{13}$C = 6.21(1.00) $\cdot D_{GC} + 18.71(7.37)$, where $D_{GC}$ is the distance to the Galactic center in kpc. The carbon isotopic ratio for NGC 6334I would be 57(10). Eight isotopic species are found in this survey, besides the main isotopes. The ratios used for NGC 6334I are: $[^{35}Cl]/[^{37}Cl] = 3$, $[^{32}S]/[^{34}S] = 23$, $[^{32}S]/[^{33}S] = 127$, $[^{12}Cl]/[^{13}Cl] = 60$, $[^{14}N]/[^{15}N] = 300$, $[^{16}O]/[^{18}O] = 500$ and $[^{16}O]/[^{17}O] = 1500$. Deuterium is found also, but here no common ratio is taken because the chemical fractionation is very high. The ortho to para ratio for some molecules is assumed to be the statistical equilibrium value (high temperature limit).

## 3.2. Model fitting with MAGIX

Once an initial guess is set up, the parameters can be fitted using the program MAGIX (Modeling and Analysis Generic Interface for eXternal numerical codes; Bernst et al. 2011; Möller et al. 2012). It delivers a package of different optimization algorithms and an interface to XCLASS. For every of the five fit parameters a range can be set over which to search for the best fit, or the parameter can be kept constant. Used here is the Levenberg-Marquardt (LM) algorithm with the method of least squares, which tries to minimize the squares of the difference between the fit function $y(x_i|a_0...a_{M-1})$ with the M fit parameters a and the data. The sum over all n data points ($x_i, y_i$) is defined as:

$$ \chi^2 = \sum_{i=0}^{n-1} \frac{(y_i - y(x_i|a_0...a_{M-1}))^2}{\sigma_i^2} \qquad (7) $$

The uncertainty for $y_i$ is $\sigma_i$. The iteration stops when a minimum is found and the algorithm has converged, or when the change in $\chi^2$ from one iteration step to another is very small, e.g. $\Delta\chi^2 < 10^{-6}$. Since this function is a priori unknown, it becomes clear that the algorithm can converge into a local and not a global minimum. A bad fit for all parameter values can be a hint for non-LTE populations, like HCN lines with high critical densities, or maser activity. The automated fitting has to be corrected manually sometimes, since different problems can occur, which are described in the following.

**Starting values:** the choice of the initial values for the parameters, is fundamentally important in the LM algorithm. Bad or random starting values lead to an unnecessary amount of computation time. The starting values are set manually by the user and should be physically reasonable. Because of their fit behavior, the five fit parameters in XCLASS are divided into two groups: $\theta$, $T_{ex}$, N and $\Delta\nu$, $v_{LSR}$. The last two are generally well fixed and do not change anymore, so they are held constant. For high opacities, the line width is broadened which is known as optical depth broadening. To check for a local or global minimum, MAGIX was iterated again with another set of starting values and then

checked if it still converged to the same point. The disadvantage is that this takes computation time with several repetitions and the variation possibilities increase dramatically with more velocity components for a given molecule. This check has been done regularly with molecules with one component, but infrequently for more complicated cases. The best fit parameters for the species with multiple overlapping components may thus not be uniquely determined.

**Calculation time:** the time increases with the number of parameters and data points. To reduce the calculation time, especially for “weeds”, not all transitions from every band are included but a variety of optical thin and thick lines, lowest and highest frequencies and accordingly a big range of different energy states. Suitable for this are Q-branches of methanol or dimethyl ether. Weeds is the colloquial term for interstellar molecules, which show a plethora of spectral lines. They are all organic asymmetric rotors with a methyl group (-CH$_3$), except of SO$_2$. The classification is not strict, but these five are considered to be the most important ones: methanol, dimethyl ether, ethanol, methyl formate and sulfur dioxide. In methanol, CH$_3$OH, the methyl group can tunnel the potential barrier at every 120° relatively to the hydroxyl group (-OH), leading to an internal rotation or torsion about the C-O bond. The torsional levels are split into sub-levels (A and E symmetry), and the interaction of torsional and rotational motions yields multiple energy states and a complex spectrum. Furthermore, heavier and larger molecules have a larger moment of inertia, which is inverse to the rotational constants, so that they produce denser spectra. The determination of the excitation temperature with XCLASS is much faster than the so called Boltzmann plot or rotational diagram method, where the upper state abundance $N_u$ of a line transition can be derived from the integrated area in the optical thin limit. This method has the disadvantages that it is unpractical for many lines and misleading for line blends and high optical depths.

**Uniqueness:** the model function Eq. (1) has two limits where degeneracy occurs: in the optically thick limit $\tau \gg 1$, the last term vanishes and $J(T_{ex})$ and $\eta(\theta)$ enter linearly into the equation; in the optically thin limit $\tau \ll 1$, $(1 - e^{-\tau}) \approx \tau$ and $\eta(\theta)$ and N enter linearly. This results in the fact that these two parameter pairs cannot be determined independently. In the fitting procedure this can be noticed in the behavior that $\chi^2$ changes only slightly while one of the parameter increases and the other decreases. It is not always obvious if a spectral line is optically thick, especially when several components add together. A solution to the problem of degeneracy are isotopologues, if present. With the ratios mentioned above, a consistent model can be fitted with constraints on N and also for $\Delta\nu$. In some cases the rare isotopologues are fitted well, but the main isotopologue is optically thick and too weak to explain the high observed intensities; this is the case for CH$_3$OH, SO, H$_2$CO and HCO$^+$.
This can be a hint that a second component is present when no signs for that could be deduced from the line shape. However, every additional emission component should be justified, since this makes the model more complex and increases the possibility of redundancy. In the case of NGC 6334I, two near overlapping components are justified by the SMA observations which show that both cores SMA1 and SMA2 show line emission at different velocities.

**Weighting:** $\chi^2$ is the sum over all modeled spectral lines. One weighting could be done by the uncertainty of $y_i$, $\sigma_i$, which is currently not implemented in MAGIX. But this is not a severe limitation, because the noise or rms is nearly constant across a given HIFI band. It is relevant when the line intensity of the
<figure><img src="image_4.png" /><figcaption>Fig. 3. a) Fit example around 497 GHz showing the data (black) and the fit (red). The decomposition into the contributions of individual species is shown above in b).</figcaption></figure>

main isotopologue is much higher than for the less abundant isotopologue. Stronger lines are then more weighted compared to weak lines because the change in $\chi^2$ is larger. This can lead to a model that is biased toward fitting the stronger lines. Another problem is the fitting in higher-frequency bands with larger noise (band 5–7). Here the algorithm tends to increase the line width because it smooths the line wings better. Line widths from transitions in lower bands reveal then that this is an overestimation. To fix this, the line width must be held constant and the best fit value is taken from transitions in lower bands.

Error estimate and confidence limits: it is assumed that in general the standard deviation of $\theta$, T and N is at least 10%, because a small higher or lower value does not change the fit significantly. Additionally, there are other errors not accounted for in the fit like calibration uncertainties, uncertain isotopic ratios, N(H₂), etc., so that in the worst case only the magnitude of N is roughly known.

Because lines of one species are blended with lines of other molecules, the fitting procedure is done in three major steps:

1. Fitting intense lines and weeds. By modeling weeds one filters them out and takes into account the line blendings. For intense lines of diatomic molecules like SO, the effect of line blending is small. Blended lines are first excluded from the fitting range.

2. Fitting weaker and absorption lines. For weak lines, the noise becomes important so that the fitting procedure is less reliable. Too weak lines, below $2\sigma$, of a molecule are therefore excluded from the fitting range, except for cases where they can serve to give an upper limit for T or N. Absorption lines are not expected to be associated with the hot cores in NGC 6334I because they come mainly from cold foreground gas. Apart from that are ground state transitions with low energy levels from H₂O or NH₃.

3. Repeating all in a second iteration. The difference now is that the simulated spectrum of all molecules are used as a new background. The parameters are saved globally in a file and read in by MAGIX, but only the molecule of interest is

fitted again. This time, blended lines of a species are taken into account and more components added, if necessary. A simultaneous fit of all parameters of all components of all molecules is in principle possible, but far too time consuming and inefficient. In the end, the result is an improved and consistent fit, see the extract in Fig. 3. This could be done iteratively several times, but it was noticed that the changes are marginal and larger only for molecules with very few transitions present.

### 3.3. Analysis of interferometric data

The presence of several components with small radial velocity differences which overlap makes the assignment to one certain physical source difficult. The advantage of interferometer data is having the spectral and spatial information together. On the other hand, the SMA frequency range of 4 GHz is small compared to HIFI and extended emission can be missed, because it is filtered out by the interferometer. Two spectra were extracted, one from the peak emission toward SMA1 and the other toward the peak emission of SMA2. It must be noticed that for some molecules like methanol, the emission overlaps from both cores. A direct comparison between them reveals that firstly SMA1 is constantly red-shifted compared to SMA2, secondly the line width in SMA1 is broader and thirdly some lines are only present in SMA1 or SMA2. The intensities in Jy/beam are converted to Kelvin by using the formula

$$T_{\text{mb}}(\text{K}) = \frac{1.22 \times 10^6 \cdot F(\text{Jy})}{[\nu(\text{GHz}) \cdot \theta(\text{''})]^2}, \quad (8)$$

where $\theta$ is the beam size. The conversion factor is then about 12.5 K/Jy for the LSB and 8.5 K/Jy for the USB. The rms is about 2.5 K in the LSB and 2 K in the USB. Maps of the integrated intensity over one spectral line of all detected molecules were made with the toolkit Karma (Gooch 1996) and the source sizes were derived from Gaussian fits. The whole procedure explained in Sect. 3.2 was performed again with the SMA spectra.
The model derived from the HIFI spectrum was extrapolated to the lower frequencies of the SMA bands and compared with the observed one. Deviations can be explained by the fact that the HIFI model includes the outer parts like the envelope and is an average over a range of temperatures, whereas the SMA spectra correspond to the inner hot core with higher temperatures.

The new derived source sizes $\theta$ were adapted for the HIFI spectrum. For some molecules, the source sizes derived from the HIFI spectrum were already in good agreement with the one from the SMA data. Small deviations (≤2′′) are permitted when they lead to a lower $\chi^2$ value. The source size is likely to change when excitation conditions differ, so that the high-J transitions in the HIFI spectra may come from a more compact region than that observed in the SMA maps. For the cases where the lines are optically thin or thick, $\theta$ is held fixed to cancel out the ambiguity.

# 4. Results

## 4.1. General results

A total of 41 molecules with their isotopologues are identified in the HIFI spectrum. The best fit results are given in Table 3. The whole survey would need several hundred pages of figures and tables with line lists. Therefore, to give an impression, only selected fits for each molecule are attached in Appendix B. The modeled spectrum for each molecule and the full observed and modeled spectrum for NGC 6334I will be made available online, after reprocessing the HIFI spectra with HIPE 8. See Appendix B for further details. The average source velocity (considering all components in the range −5 to −10 km s⁻¹ associated with NGC 6334I) is −7.3 km s⁻¹ and average line width 4.2 km s⁻¹ with a standard deviation of about ±1 km s⁻¹ for both spectral features. The velocity width is dominated by turbulence, since the thermal Doppler width contributes less than 1 km s⁻¹ for typical values of T = 100 K and m = 20 amu. Some molecules show strong broad wings which are produced by outflows: CO, CS, SiO, H₂O and HCN. H₂CS has an asymmetric line profile with a stronger redshifted wing. Especially the component at −17 to −20 km s⁻¹, observed in C, NH₃ and H₂S with T = 10–30 K is presumably the blue outflow lobe. Two cold (~5 K), extended foreground clouds are detected consistently with components around −2 km s⁻¹ and +7 km s⁻¹ in the following species: C and C⁺, CH and CH⁺, H₂O and H₂O⁺, HF, H₂Cl⁺ and OH⁺. This is in agreement with an observed OH absorption component at +6 km s⁻¹ toward several NGC 6334 cores reported by Brooks & Whiteoak (2001).

The best fit results for the SMA spectra are presented in Table A.1, the list of detected lines in Table A.2 and the spectra for the LSB and USB in Fig. A.1 in Appendix A. Overall, 20 molecules are identified, of which 15 are also present in the HIFI observations. The additional five molecules are all organic: acetone (CH₃CH₃CO), ethylene glycol ((CH₂OH)₂), formic acid (HCOOH), acrylonitrile (C₂H₃CN) and propionitrile (C₂H₅CN). The latter two only have low intensity features (<3σ) in the HIFI band 1a spectrum and are therefore considered to be not detected. Also, 3 isotopologues of OCS are found. In general, the fits give higher temperatures than those derived from the HIFI data, up to 200–300 K, because the emission originates from the center of the hot cores, observed in a smaller beam. Although not true for all molecules, SMA2 has on average higher excitation temperatures, lower column densities, abundances and smaller line widths than SMA1. On average, the source velocity for SMA1 is −5.6 ± 0.1 km s⁻¹ and

for SMA2 −7.6 ± 0.1 km s⁻¹, so that the relative radial velocity difference is 2 km s⁻¹. When applied e.g. to ethanol (C₂H₅OH), which is only seen toward SMA2, the source velocity derived from the SMA spectra is consistent with that from HIFI. In the case of C₂H₄O or NH₂CHO, a small discrepancy of 0.8 km s⁻¹ is seen, which can be explained by the fact that in the SMA bands only a few transitions are detected, some of which are blended, so that the precise velocity determination is hindered.

The chemical differences between both cores are determined from the SMA maps and the velocity components in the HIFI spectra. Those molecules that show dominant emission or only emission in SMA2 are HCOOH, (CH₂OH)₂, NH₂CHO, C₂H₄O and C₂H₅OH. For SMA1, the dominant species are HC₃N, C₂H₃CN, C₂H₅CN, H₂S, CN, HCN and HNC. It is noticeable that more oxygen-bearing species, including water (Emprechtinger et al. 2012), arise in SMA2 and more nitrogen-bearing species in SMA1. This variation in the peak abundance is also observed in other star-forming regions, such as Orion KL between the hot core and the compact ridge (Comito et al. 2005).

## 4.2. Statistics

A spectral line is considered to be detected when the peak intensity exceeds the 3σ level of the random noise. Furthermore, the line width should be comparable to that of all other detected lines. In this way one excludes single channel spikes and sets a lower limit for the integrated line intensity. A lower limit of 2σ would be acceptable when this limit is not taken for a single line of a molecule but rather for many lines since the probability that this is a coincidence decreases when having many lines observed and a good model fit. Nevertheless, some weak features are included in certain cases in the fitting range because they can give useful upper limits for the model.

Overall, 49 chemical species have been found in NGC 6334I so far, listed in Table 4, ranging from fine structure transitions of C, through over mostly di- and triatomic molecules, to organic molecules with 10 atoms. About half of them were found in previous line surveys. Methyl acetylene (CH₃CCH) has been identified by Kalinina et al. (2010), but is too weak in HIFI (<1σ). We find 31 isotopologues in our data. The number of lines in all HIFI bands derived from the XCLASS fits is 4217, see Table 5. Hyperfine structure splitting is not considered and counted as one line. For the number of unidentified lines see next section. An automated line counting on the observational spectra gives the same number in orders of magnitude, but is too inaccurate at the band edges and for blended lines. The integrated line intensities are derived from the simulated spectrum, and for bad fits like CO, CS, HCN etc. from Gaussian fits. Because the fits are not perfect, the presented values should be viewed with caution. The line density is 5.5/GHz in the first 5 HIFI bands, ranging from 17/GHz in 1a to <1/GHz in 5a. Bands 6 and 7 have a high rms and less than 10 lines each. It is evident that the five weeds mentioned in Sect. 3.2 account for over 75% of all lines, and that methanol is the dominant emitter besides CO. Methyl formate has a very dense spectrum, but is less intense at HIFI frequencies than in both SMA bands where it has the highest number of lines.

## 4.3. Unidentified features

Table 6 lists the frequencies and strengths of the 75 unidentified lines in the HIFI spectra and the 27 unidentified lines in the SMA spectra. The strength of the lines are derived from
<table><caption>**Table 3.** Fit results for the HIFI spectrum of NGC 6334I.</caption><thead><tr><th>Species</th><th>Θ<br>(")</th><th>T<sub>ex</sub><br>(K)</th><th>N<br>(cm<sup>-2</sup>)</th><th>Δv<br>(km s<sup>-1</sup>)</th><th>v<sub>LSR</sub><br>(km s<sup>-1</sup>)</th><th>Notes</th><th>Species</th><th>Θ<br>(")</th><th>T<sub>ex</sub><br>(K)</th><th>N<br>(cm<sup>-2</sup>)</th><th>Δv<br>(km s<sup>-1</sup>)</th><th>v<sub>LSR</sub><br>(km s<sup>-1</sup>)</th><th>Notes</th></tr></thead><tbody><tr><td rowspan="5">C</td><td>40</td><td>56</td><td>2.4 × 10<sup>18</sup></td><td>4.3</td><td>-7.3</td><td rowspan="5">b</td><td>HCS<sup>+</sup></td><td>4.0</td><td>68</td><td>7.6 × 10<sup>14</sup></td><td>4.3</td><td>-7.1</td><td rowspan="5">a</td></tr><tr><td>29</td><td>68</td><td>2.1 × 10<sup>17</sup></td><td>4</td><td>0.2</td><td>NO</td><td>3.5</td><td>50</td><td>2.7 × 10<sup>17</sup></td><td>3.9</td><td>-7.5</td></tr><tr><td>ext</td><td>6</td><td>2.3 × 10<sup>17</sup></td><td>1</td><td>-2</td><td>SiO</td><td>3.4</td><td>30</td><td>3.0 × 10<sup>15</sup></td><td>5</td><td>-7.4</td></tr><tr><td>ext</td><td>2.7</td><td>4.7 × 10<sup>15</sup></td><td>1.3</td><td>6.5</td><td>HF</td><td>ext</td><td>2.7</td><td>1.2 × 10<sup>14</sup></td><td>8.6</td><td>-5</td></tr><tr><td>48</td><td>26</td><td>4.0 × 10<sup>16</sup></td><td>3</td><td>-20.6</td><td>ext</td><td>2.7</td><td>1.2 × 10<sup>13</sup></td><td>1.9</td><td>0.3</td><td></td></tr><tr><td rowspan="3">C<sup>+</sup></td><td>ext</td><td>112</td><td>1.8 × 10<sup>18</sup></td><td>6</td><td>-9</td><td rowspan="3">1</td><td>HCl</td><td>22</td><td>39</td><td>1.9 × 10<sup>14</sup></td><td>4.3</td><td>-6</td><td rowspan="3">6</td></tr><tr><td>ext</td><td>4</td><td>6.0 × 10<sup>17</sup></td><td>3</td><td>-0.5</td><td>H<sub>2</sub>Cl<sup>+</sup></td><td>ext</td><td>11</td><td>1.5 × 10<sup>14</sup></td><td>8.6</td><td>-10</td></tr><tr><td>ext</td><td>5</td><td>6.0 × 10<sup>17</sup></td><td>4.6</td><td>6.6</td><td>SH<sup>+</sup></td><td>ext</td><td>2.7</td><td>2.5 × 10<sup>12</sup></td><td>4</td><td>-5.7</td></tr><tr><td rowspan="2">CH<sup>+</sup></td><td>ext</td><td>13</td><td>1.5 × 10<sup>14</sup></td><td>5</td><td>-9.8</td><td rowspan="2">1</td><td>H<sub>2</sub>S</td><td>3</td><td>95</td><td>2.5 × 10<sup>17</sup></td><td>4</td><td>-6.1</td><td rowspan="2">8</td></tr><tr><td>ext</td><td>8</td><td>1.6 × 10<sup>14</sup></td><td>6.6</td><td>-2</td><td>ext</td><td>20</td><td>1.0 × 10<sup>15</sup></td><td>5</td><td>-7.8</td></tr><tr><td rowspan="2">CH</td><td>13.4</td><td>95</td><td>1.9 × 10<sup>14</sup></td><td>2.9</td><td>-8.2</td><td rowspan="2">2</td><td></td><td>18.5</td><td>11</td><td>8.6 × 10<sup>14</sup></td><td>17</td><td>-17.5</td><td rowspan="2">b</td></tr><tr><td>64</td><td>4</td><td>6.7 × 10<sup>13</sup></td><td>1.9</td><td>6.7</td><td>CS</td><td>2.4</td><td>87</td><td>3.6 × 10<sup>17</sup></td><td>4.5</td><td>-6.8</td><td rowspan="2">9</td></tr><tr><td rowspan="2">NH</td><td>39</td><td>2.7</td><td>6.0 × 10<sup>13</sup></td><td>2</td><td>-6</td><td rowspan="2">2</td><td></td><td>2.2</td><td>73</td><td>1.4 × 10<sup>17</sup></td><td>3.8</td><td>-7.8</td><td rowspan="2">a</td></tr><tr><td>ext</td><td>2.7</td><td>9.0 × 10<sup>13</sup></td><td>10</td><td>-4.6</td><td>ext</td><td>103</td><td>6.5 × 10<sup>13</sup></td><td>6</td><td>-6.8</td><td></td></tr><tr><td rowspan="2">NH<sub>2</sub></td><td>3.4</td><td>108</td><td>3.6 × 10<sup>15</sup></td><td>4</td><td>-8.5</td><td rowspan="2">2</td><td>NS</td><td>3</td><td>89</td><td>3.4 × 10<sup>15</sup></td><td>5</td><td>-7</td><td rowspan="2">a</td></tr><tr><td>ext</td><td>12</td><td>2.9 × 10<sup>13</sup></td><td>2</td><td>-6.1</td><td>OCS</td><td>2.5</td><td>102</td><td>1.2 × 10<sup>18</sup></td><td>4</td><td>-7.5</td><td></td></tr><tr><td rowspan="2">NH<sub>3</sub></td><td>8</td><td>28</td><td>2.3 × 10<sup>16</sup></td><td>4.2</td><td>-7.7</td><td rowspan="2">b</td><td>SO</td><td>2.7</td><td>70</td><td>2.1 × 10<sup>17</sup></td><td>4.1</td><td>-7.4</td><td rowspan="2">10</td></tr><tr><td>ext</td><td>9</td><td>6.0 × 10<sup>13</sup></td><td>2.9</td><td>-6.9</td><td>SO<sub>2</sub></td><td>1.8</td><td>100</td><td>9.5 × 10<sup>16</sup></td><td>4.7</td><td>-7.8</td><td></td></tr><tr><td rowspan="2">NH<sub>2</sub>D</td><td>8</td><td>28</td><td>2.3 × 10<sup>14</sup></td><td>3.8</td><td>-8.6</td><td rowspan="2">b</td><td>H<sub>2</sub>CS</td><td>3.4</td><td>106</td><td>2.7 × 10<sup>16</sup></td><td>5</td><td>-7.9</td><td rowspan="2">11</td></tr><tr><td>ext</td><td>9</td><td>2.2 × 10<sup>14</sup></td><td>16</td><td>-17</td><td>CH<sub>3</sub>OH</td><td>2.7</td><td>99</td><td>1.4 × 10<sup>19</sup></td><td>2.9</td><td>-7.9</td><td rowspan="2">12</td></tr><tr><td rowspan="2">N<sub>2</sub>H<sup>+</sup></td><td>31</td><td>29</td><td>4.3 × 10<sup>13</sup></td><td>2.5</td><td>-6.8</td><td rowspan="2">b</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>2.1</td><td>87</td><td>3.0 × 10<sup>18</sup></td><td>3.2</td><td>-7.4</td><td rowspan="2">13</td></tr><tr><td>ext</td><td>20</td><td>35</td><td>3.5 × 10<sup>13</sup></td><td>2.5</td><td>-9.5</td><td>CH<sub>3</sub>OCHO</td><td>5.0</td><td>69</td><td>7.0 × 10<sup>16</sup></td><td>2.7</td><td>-8.1</td><td></td></tr><tr><td rowspan="2">OH<sup>+</sup></td><td>ext</td><td>2.7</td><td>2.0 × 10<sup>13</sup></td><td>4</td><td>3.3</td><td rowspan="2">3</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>3</td><td>115</td><td>7.1 × 10<sup>17</sup></td><td>3.5</td><td>-7.65</td><td rowspan="2">8</td></tr><tr><td>ext</td><td>2.7</td><td>2.8 × 10<sup>13</sup></td><td>6.5</td><td>-2</td><td>H<sub>2</sub>CO</td><td>2</td><td>117</td><td>1.4 × 10<sup>17</sup></td><td>3.9</td><td>-8</td><td></td></tr><tr><td rowspan="2">H<sub>2</sub>O</td><td>4.2</td><td>59</td><td>2.1 × 10<sup>18</sup></td><td>4.3</td><td>-8.3</td><td rowspan="2">a</td><td></td><td>2.7</td><td>78</td><td>1.5 × 10<sup>17</sup></td><td>4.5</td><td>-7.4</td><td rowspan="2">14</td></tr><tr><td>ext</td><td>4.4</td><td>440</td><td>1.2 × 10<sup>15</sup></td><td>5.3</td><td>-7.7</td><td>CH<sub>3</sub>CN</td><td>1.9</td><td>154</td><td>2.8 × 10<sup>16</sup></td><td>4.7</td><td>-7.3</td><td></td></tr><tr><td rowspan="2">H<sub>2</sub>O<sup>+</sup></td><td>ext</td><td>24</td><td>3.7 × 10<sup>13</sup></td><td>8.3</td><td>-0.7</td><td rowspan="2">4</td><td>C<sub>2</sub>H<sub>4</sub>O</td><td>2.1</td><td>80</td><td>5.4 × 10<sup>16</sup></td><td>3</td><td>-7.6</td><td rowspan="2">a</td></tr><tr><td>ext</td><td>8.7</td><td>32</td><td>4.3 × 10<sup>15</sup></td><td>5</td><td>-6.9</td><td>H<sub>2</sub>CCO</td><td>3</td><td>117</td><td>6.0 × 10<sup>15</sup></td><td>3</td><td>-7.7</td><td></td></tr><tr><td rowspan="2">C<sub>2</sub>H</td><td>14.8</td><td>22</td><td>1.1 × 10<sup>15</sup></td><td>5</td><td>-6.4</td><td rowspan="2">5</td><td>HCO<sup>+</sup></td><td>5.1</td><td>73</td><td>1.3 × 10<sup>15</sup></td><td>4.7</td><td>-7.5</td><td rowspan="2">15</td></tr><tr><td>ext</td><td>32</td><td>60</td><td>3.7 × 10<sup>19</sup></td><td>4.6</td><td>-6.8</td><td>CH<sub>2</sub>NH</td><td>2.4</td><td>75</td><td>6.0 × 10<sup>15</sup></td><td>5</td><td>-7.0</td><td rowspan="2">a</td></tr><tr><td rowspan="2">CO</td><td>3.5</td><td>32</td><td>2.6 × 10<sup>18</sup></td><td>6</td><td>-5.6</td><td rowspan="2">5</td><td>HNCO</td><td>2.0</td><td>284</td><td>1.5 × 10<sup>16</sup></td><td>5</td><td>-7.7</td><td rowspan="2">a</td></tr><tr><td>ext</td><td>5.0</td><td>48</td><td>3.0 × 10<sup>16</sup></td><td>5</td><td>-7.5</td><td>NH<sub>2</sub>CHO</td><td>2</td><td>117</td><td>4.7 × 10<sup>15</sup></td><td>4.2</td><td>-7.7</td><td></td></tr><tr><td rowspan="2">DCN</td><td>3.5</td><td>54</td><td>3.7 × 10<sup>14</sup></td><td>6</td><td>-5.6</td><td rowspan="2">5</td><td>HC<sub>3</sub>N</td><td>2.8</td><td>134</td><td>4.2 × 10<sup>15</sup></td><td>9</td><td>-5.6</td><td rowspan="2">a</td></tr><tr><td>ext</td><td>3.5</td><td>50</td><td>1.1 × 10<sup>15</sup></td><td>6</td><td>-5.6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">HNC</td><td>5.0</td><td>40</td><td>3.3 × 10<sup>14</sup></td><td>4.4</td><td>-8.2</td><td rowspan="2">5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ext</td><td>5.0</td><td>40</td><td>3.3 × 10<sup>14</sup></td><td>4.4</td><td>-8.2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Several rows for one species denote multiple components. Extended source sizes are abbreviated with ext. (1) Including $^{13}$CH$^{+}$; (2) including $^{15}$NH$_{3}$; (3) only isotopologues $^{18}$O, $^{17}$O; (4) only isotopologues $^{13}$CO, $^{18}$O, $^{13}$C$^{18}$O, $^{17}$O; (5) only isotopologues $^{13}$CN, $^{15}$N, $^{13}$C$^{15}$N; (6) including $^{37}$Cl; (7) including $^{37}$Cl$^{+}$; (8) including $^{34}$S; (9) only isotopologues $^{34}$S, $^{33}$S, $^{13}$CS, $^{13}$C$^{34}$S; (10) including $^{34}$SO, $^{33}$SO; (11) including $^{2}$C$^{34}$S; (12) including $^{13}$CH$_{3}$OH; (13) including $^{13}$CH$_{3}$OCH$_{3}$; (14) including $^{13}$CO; (15) including $^{13}$CO$^{+}$. (a) τ small. (b) outflow component.

<table><caption>**Table 4.** List of all detected chemical species in NGC 6334I.</caption><thead><tr><th>1 atom</th><th>2 atoms</th><th>3 atoms</th><th>4 atoms</th><th>5 atoms</th><th>6 atoms</th><th>7 atoms</th><th>8 atoms</th><th>9 atoms</th><th>10 atoms</th></tr></thead><tbody><tr><td>C</td><td>SO, NO</td><td>SO<sub>2</sub>, HCO<sup>+</sup></td><td>H<sub>2</sub>CO</td><td>CH<sub>2</sub>NH</td><td>CH<sub>3</sub>OH</td><td>C<sub>2</sub>H<sub>4</sub>O</td><td>CH<sub>3</sub>OCHO</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>CH<sub>3</sub>CH<sub>3</sub>CO</td></tr><tr><td>C<sup>+</sup></td><td>CS, SiO</td><td>OCS, H<sub>2</sub>O</td><td>H<sub>2</sub>CS</td><td>HC<sub>3</sub>N</td><td>CH<sub>3</sub>CN</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td></td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>(CH<sub>2</sub>OH)<sub>2</sub></td></tr><tr><td></td><td>NS, OH<sup>+</sup></td><td>HCN, NH<sub>2</sub></td><td>HNCO</td><td>HCOOH</td><td>NH<sub>3</sub></td><td>CH<sub>3</sub>CCH</td><td></td><td>C<sub>2</sub>H<sub>5</sub>CN</td><td></td></tr><tr><td></td><td>CO, CH<sup>+</sup></td><td>HNC, HCS<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>CN, CH</td><td>H<sub>2</sub>S, H<sub>2</sub>O<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>HCl, HF</td><td>C<sub>2</sub>H, H<sub>2</sub>Cl<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>NH, SH<sup>+</sup></td><td>N<sub>2</sub>H<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Gaussian fits. 211 unidentified features were detected initially, going up from band 1a to 4b, and none in the HEB Bands. All of them are emission lines except of one weak absorption line at 507 686 MHz below $5\sigma$. Most lines are relatively weak. The

highest peak-to-noise ratio of 14 has the line at 492 325 MHz. Erroneous features are not listed here, but can be checked by comparing the bands in *H* and *V* polarization or at the edges where the bands overlap. In this way, 5 false lines (called ghosts)
<table><caption>**Table 5.** Integrated spectral line intensities of all detected chemical species of the HIFI spectrum for NGC 6334I.</caption><thead><tr><th>Species</th><th>∫Tdv<br>(K km s<sup>-1</sup>)</th><th>Lines</th><th>Species</th><th>∫Tdv<br>(K km s<sup>-1</sup>)</th><th>Lines</th></tr></thead><tbody><tr><td>CH<sub>3</sub>OH</td><td>6016</td><td>2063</td><td>CH<sub>3</sub>CN</td><td>50</td><td>48</td></tr><tr><td>CO</td><td>4715</td><td>31</td><td>HCl</td><td>40</td><td>3</td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>1486</td><td>637</td><td>N<sub>2</sub>H<sup>+</sup></td><td>23</td><td>4</td></tr><tr><td>H<sub>2</sub>O</td><td>1257</td><td>31</td><td>OCS</td><td>22</td><td>16</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>582</td><td>558</td><td>H<sub>2</sub>O<sup>+</sup></td><td>16</td><td>1</td></tr><tr><td>H<sub>2</sub>CO</td><td>428</td><td>105</td><td>NS</td><td>16</td><td>12</td></tr><tr><td>CS</td><td>280</td><td>45</td><td>CN</td><td>15</td><td>6</td></tr><tr><td>C<sup>+</sup></td><td>262</td><td>1</td><td>HNC</td><td>15</td><td>4</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>227</td><td>134</td><td>CH<sub>2</sub>NH</td><td>14</td><td>17</td></tr><tr><td>NH<sub>3</sub></td><td>226</td><td>13</td><td>CH</td><td>14</td><td>1</td></tr><tr><td>SO<sub>2</sub></td><td>207</td><td>109</td><td>OH<sup>+</sup></td><td>12</td><td>2</td></tr><tr><td>HCN</td><td>205</td><td>27</td><td>NO</td><td>11</td><td>8</td></tr><tr><td>NH<sub>2</sub></td><td>181</td><td>6</td><td>H<sub>2</sub>Cl<sup>+</sup></td><td>11</td><td>2</td></tr><tr><td>H<sub>2</sub>S</td><td>168</td><td>32</td><td>C<sub>2</sub>H</td><td>11</td><td>6</td></tr><tr><td>HCO<sup>+</sup></td><td>152</td><td>10</td><td>HNCO</td><td>10</td><td>25</td></tr><tr><td>HF</td><td>146</td><td>1</td><td>H<sub>2</sub>CCO</td><td>4</td><td>9</td></tr><tr><td>H<sub>2</sub>CS</td><td>141</td><td>78</td><td>HCS<sup>+</sup></td><td>4</td><td>4</td></tr><tr><td>SO</td><td>119</td><td>51</td><td>NH<sub>2</sub>CHO</td><td>2</td><td>8</td></tr><tr><td>C</td><td>115</td><td>2</td><td>H<sub>3</sub>N</td><td>2</td><td>4</td></tr><tr><td>CH<sup>+</sup></td><td>108</td><td>3</td><td>SiO</td><td>2</td><td>1</td></tr><tr><td>NH</td><td>103</td><td>3</td><td>SH<sup>+</sup></td><td>1</td><td>1</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>67</td><td>94</td><td></td><td></td><td></td></tr><tr><td>Sum</td><td>17 485</td><td>4217</td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Included are all transitions from isotopologues (except CH<sub>3</sub>OD and CH<sub>3</sub><sup>18</sup>OH) and torsionally or vibrationally excited states.

at 720 GHz are noticed which are present in 2a, but not 2b. They arose presumably from a feedback in the LO and are no longer present after calibration with a newer version of HIPE. Notable are patterns like branches where several lines (5–10) appear at once in a narrow frequency range (1 GHz), namely at 510, 532, 537, 602, 641, 695, 720, 1009 and 1069 GHz. Other clear patterns like a recurring difference between lines could not be discovered. If there are several candidates for an unidentified line, the simulation with XCLASS helps to constrain them. Any anticoincidence, which means that a simulated line in other quantum numbers in other frequency bands does not match with the observed spectrum, is an indication for a false identification. That is why complete line surveys are better suited for confident line identifications than small windows or single line searches. Another constraint from interferometric data is that all lines of a species originate from the same source. Since weeds have the highest number of lines, it is most probable that the unidentified lines originate from them, or one of their isotopologues not yet included in the database. Two papers were checked with reported line frequencies of deuterated methanol CH<sub>3</sub>OD (Walsh et al. 2000) and the isotopologue CH<sub>3</sub><sup>18</sup>OH (Fisher et al. 2007). 8 lines could be assigned to deuterated methanol and 128, especially most of the former reported branches, to CH<sub>3</sub><sup>18</sup>OH. Subtracting them, the remaining number of unidentified lines is 75 lines (including the 1009 GHz and 1069 GHz branches) or about 2% of the total of 4428 lines. In the SMA bands, 144 transitions are counted of which 27 are unidentified. That makes 19% and is much more than in HIFI frequency range. It is probable that they originate from organic molecules in higher excited states. All unidentified lines in the SMA spectra are detected in both cores, but it is remarkable that, with the exception of two, they are all stronger in SMA2 than SMA1.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Derived excitation temperatures in dependence of the source size for molecules associated with NGC 6334I.</figcaption></figure>

# 5. Discussion

## 5.1. General remarks

Schilke et al. (2006) analyzed spectra of NGC 6334I with APEX at 460 and 810 GHz with a bandwidth of 2.5 GHz. 12 molecules were identified and roughly fitted, of which 6 show an acceptable agreement (same order of magnitude in *N*) with our results: CH<sub>3</sub>CN, CH<sub>3</sub>OCH<sub>3</sub>, CH<sub>3</sub>OCHO, HNCO, N<sub>2</sub>H<sup>+</sup> and SO<sub>2</sub>. The deviation for other molecules can be explained by an over- or underestimation of the source size and by setting the temperature in steps of 50 K. In general, there is a good agreement between the HIFI and SMA column densities. The temperatures given in Table A.1 are not well determined, because for many molecules only one transition is present.

Because no error bars are given, it should be pointed out that the given parameters for some molecules are less reliable than for others. The spectra of strong lines (CO, N<sub>2</sub>H<sup>+</sup>, HCN, C<sup>+</sup>...) could be affected by emission in the off-beam position. In some cases, the assumption of LTE does not hold, e.g. for HCN, or the ambiguity leads to an uncertainty. Reliable values are these which are checked consistently with the SMA and have isotopologues, so that better constraints are available for θ and N.

In Fig. 4, the excitation temperature is plotted as a function of the source size for all molecules associated with NGC 6334I. The extended source sizes (θ ≫ beam size) are not included. The data points can be roughly fitted by a function $f(x) = a \cdot x^b + c$ with an exponent of $b = -1.3 \pm 0.5$ and a constant $c = (35 \pm 11)$ K for $\theta > 6''$ in the outer parts or envelope where the temperature gradient becomes flat. This is consistent with models of hot cores which predict a steep increase inwards (Rolffs et al. 2011). The profile could be affected by the fit function, where in the optical thick limit, T and θ are inversely proportional to each other with $T \propto 1/\theta$, but this only concerns few molecules.

The fractional abundances $X = N_{\text{Species}}/N_{\text{H}}$ are distributed widely, ranging over eight orders of magnitude from the nitrogen hydrides ($X = 10^{-12}$) to CO ($X = 10^{-5}$, derived from the isotopologues). The hydrogen column density for NGC 6334I is $N_{\text{H}} = 3 \times 10^{24} \text{ cm}^{-2}$, see Sect. 3.1. The most abundant molecules with the highest column densities after CO are CH<sub>3</sub>OH and CH<sub>3</sub>OCH<sub>3</sub> (both $X = 10^{-6}$), HCN and H<sub>2</sub>O (both $X = 10^{-7}$). The hydrogen column density can be estimated from the derived CO column density of $N = 3.7 \times 10^{19} \text{ cm}^{-2}$ by using a standard fractional abundance of [CO]/[H] = $5 \times 10^{-5}$.
<table><caption>**Table 6.** List of unidentified lines in the SMA and HIFI spectra for NGC 6334I.</caption><thead><tr><th>Frequency (MHz)</th><th>Area (K km s<sup>-1</sup>)</th><th>I (K)</th><th>I/rms</th><th>Width (km s<sup>-1</sup>)</th><th>Frequency (MHz)</th><th>Area (K km s<sup>-1</sup>)</th><th>I (K)</th><th>I/rms</th><th>Width (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="10"><strong>SMA</strong></td></tr><tr><td>216 671</td><td>41.2</td><td>8.1</td><td>3.2</td><td>4.8</td><td>227 128</td><td>68.0</td><td>11.3</td><td>5.7</td><td>5.7</td></tr><tr><td>216 715</td><td>71.7<sup>b</sup></td><td>13.5</td><td>5.4</td><td>5.0</td><td>227 239</td><td>26.0<sup>b</sup></td><td>8.1</td><td>4.1</td><td>3.0</td></tr><tr><td>216 913</td><td>48.0</td><td>12.4</td><td>5.0</td><td>3.6</td><td>227 244</td><td>32.3<sup>b</sup></td><td>7.0</td><td>3.5</td><td>4.3</td></tr><tr><td>216 919</td><td>82.8</td><td>24.1</td><td>9.6</td><td>3.2</td><td>227 534</td><td>61.0</td><td>14.2</td><td>7.1</td><td>4.0</td></tr><tr><td>217 000</td><td>54.9</td><td>11.4</td><td>4.6</td><td>4.5</td><td>227 728</td><td>43.7</td><td>12.8</td><td>6.4</td><td>3.2</td></tr><tr><td>217 056</td><td>90.0</td><td>20.8</td><td>8.3</td><td>4.1</td><td>227 803</td><td>41.6</td><td>11.4</td><td>5.7</td><td>3.4</td></tr><tr><td>217 262</td><td>91.5<sup>b</sup></td><td>22.3</td><td>8.9</td><td>3.9</td><td>227 823</td><td>47.8<sup>b</sup></td><td>7.4</td><td>3.7</td><td>6.1</td></tr><tr><td>217 427</td><td>38.2</td><td>11.6</td><td>4.6</td><td>3.1</td><td>227 841</td><td>99.0<sup>b</sup></td><td>21.1</td><td>10.6</td><td>4.4</td></tr><tr><td>217 496</td><td>70.9</td><td>18.8</td><td>7.5</td><td>3.5</td><td>228 001</td><td>53.0</td><td>7.6</td><td>3.8</td><td>6.6</td></tr><tr><td>217 593</td><td>88.3</td><td>24.5</td><td>9.8</td><td>3.3</td><td>228 141</td><td>112<sup>b</sup></td><td>17.0</td><td>8.5</td><td>6.2</td></tr><tr><td>217 678</td><td>43.5</td><td>9.6</td><td>3.8</td><td>4.3</td><td>228 230</td><td>24.5</td><td>7.5</td><td>3.8</td><td>3.1</td></tr><tr><td>217 689</td><td>119</td><td>27.4</td><td>11.0</td><td>4.1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>217 710</td><td>120</td><td>27.7</td><td>11.1</td><td>4.1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>217 738</td><td>76.1</td><td>18.0</td><td>7.2</td><td>4.0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>218 158</td><td>116<sup>b</sup></td><td>18.8</td><td>7.5</td><td>5.8</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>218 600</td><td>58.8</td><td>11.8</td><td>4.7</td><td>4.7</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="10"><strong>HIFI</strong></td></tr><tr><td>492 313</td><td>0.88</td><td>0.30</td><td>14.3</td><td>2.8</td><td>1 004 810</td><td>2.13</td><td>0.50</td><td>4.2</td><td>4.0</td></tr><tr><td>501 501</td><td>0.40<sup>b</sup></td><td>0.06</td><td>4.0</td><td>6.3</td><td>1 007 185</td><td>2.50</td><td>0.44</td><td>3.1</td><td>5.3</td></tr><tr><td>504 512</td><td>0.14<sup>b</sup></td><td>0.05</td><td>3.3</td><td>3.0</td><td>1 008 602</td><td>2.10</td><td>0.33</td><td>2.5</td><td>6.1</td></tr><tr><td>505 339</td><td>0.62</td><td>0.10</td><td>6.3</td><td>6.2</td><td>1 008 962</td><td>3.16</td><td>0.47</td><td>3.4</td><td>6.3</td></tr><tr><td>507 167</td><td>0.72<sup>b</sup></td><td>0.07</td><td>4.7</td><td>9.8</td><td>1 009 245</td><td>1.61</td><td>0.34</td><td>2.4</td><td>4.4</td></tr><tr><td>507 673</td><td>0.33</td><td>0.07</td><td>4.3</td><td>4.8</td><td>1 009 460</td><td>1.85</td><td>0.54</td><td>4.2</td><td>3.2</td></tr><tr><td>515 619</td><td>0.25</td><td>0.06</td><td>4.3</td><td>3.6</td><td>1 009 618</td><td>2.31</td><td>0.62</td><td>5.6</td><td>3.5</td></tr><tr><td>519 093</td><td>0.70<sup>b</sup></td><td>0.043</td><td>5</td><td>9.1</td><td>1 009 728</td><td>2.78</td><td>0.51</td><td>4.3</td><td>5.1</td></tr><tr><td>528 091</td><td>0.33</td><td>0.06</td><td>4.1</td><td>5.0</td><td>1 009 804</td><td>2.90</td><td>0.53</td><td>4.4</td><td>5.1</td></tr><tr><td>530 769</td><td>0.28<sup>b</sup></td><td>0.09</td><td>6.0</td><td>2.8</td><td>1 009 850</td><td>2.74</td><td>0.58</td><td>4.8</td><td>4.5</td></tr><tr><td>531 303</td><td>0.36<sup>b</sup></td><td>0.062</td><td>4.1</td><td>5.4</td><td>1 009 881</td><td>3.20</td><td>0.50</td><td>4.2</td><td>6.0</td></tr><tr><td>531 386</td><td>0.65</td><td>0.09</td><td>6.0</td><td>6.8</td><td>1 017 993</td><td>4.00</td><td>0.71</td><td>7.6</td><td>5.2</td></tr><tr><td>534 943</td><td>0.63<sup>b</sup></td><td>0.049</td><td>3.3</td><td>12</td><td>1 024 147</td><td>2.78</td><td>0.60</td><td>6.0</td><td>4.3</td></tr><tr><td>536 632</td><td>0.14</td><td>0.06</td><td>3.7</td><td>2.4</td><td>1 040 932</td><td>1.86</td><td>0.52</td><td>3.7</td><td>3.3</td></tr><tr><td>539 134</td><td>0.36</td><td>0.08</td><td>5.3</td><td>4.2</td><td>1 045 510</td><td>2.30</td><td>0.54</td><td>6.0</td><td>4.0</td></tr><tr><td>556 311</td><td>0.44</td><td>0.10</td><td>6.7</td><td>4.1</td><td>1 047 493</td><td>2.00</td><td>0.47</td><td>3.6</td><td>4.1</td></tr><tr><td>602 910</td><td>0.53</td><td>0.10</td><td>3.1</td><td>4.8</td><td>1 047 982</td><td>3.34</td><td>0.62</td><td>4.1</td><td>5.0</td></tr><tr><td>649 022</td><td>0.51</td><td>0.12</td><td>4.0</td><td>4.0</td><td>1 052 768</td><td>4.50<sup>b</sup></td><td>0.75</td><td>3.0</td><td>5.7</td></tr><tr><td>649 093</td><td>0.47</td><td>0.12</td><td>4.0</td><td>3.7</td><td>1 065 139</td><td>3.90</td><td>0.60</td><td>4.6</td><td>6.2</td></tr><tr><td>650 170</td><td>0.74</td><td>0.13</td><td>3.8</td><td>5.4</td><td>1 066 454</td><td>2.00</td><td>0.49</td><td>4.9</td><td>3.8</td></tr><tr><td>692 016</td><td>0.39</td><td>0.10</td><td>2.5</td><td>3.7</td><td>1 069 053</td><td>1.45</td><td>0.30</td><td>3.0</td><td>4.5</td></tr><tr><td>693 184</td><td>0.60</td><td>0.13</td><td>3.3</td><td>4.4</td><td>1 069 207</td><td>1.00</td><td>0.32</td><td>3.2</td><td>2.9</td></tr><tr><td>695 857</td><td>0.41</td><td>0.13</td><td>3.3</td><td>2.9</td><td>1 069 495</td><td>1.35</td><td>0.34</td><td>3.4</td><td>3.7</td></tr><tr><td>697 351</td><td>0.95</td><td>0.12</td><td>3.0</td><td>7.7</td><td>1 069 625</td><td>1.35</td><td>0.32</td><td>3.2</td><td>3.9</td></tr><tr><td>698 919</td><td>1.60<sup>b</sup></td><td>0.27</td><td>6.8</td><td>5.6</td><td>1 069 741</td><td>1.48</td><td>0.34</td><td>3.4</td><td>4.1</td></tr><tr><td>706 552</td><td>0.55</td><td>0.16</td><td>4.0</td><td>3.2</td><td>1 069 851</td><td>2.70<sup>b</sup></td><td>0.43</td><td>3.9</td><td>5.9</td></tr><tr><td>717 678</td><td>0.74</td><td>0.24</td><td>5.3</td><td>2.9</td><td>1 076 527</td><td>3.80</td><td>0.85</td><td>7.7</td><td>4.2</td></tr><tr><td>741 290</td><td>1.24</td><td>0.15</td><td>3.0</td><td>7.9</td><td>1 078 549</td><td>2.63</td><td>0.51</td><td>4.6</td><td>4.8</td></tr><tr><td>743 318</td><td>0.81</td><td>0.15</td><td>3.0</td><td>5.2</td><td>1 088 059</td><td>2.00</td><td>0.55</td><td>4.6</td><td>3.4</td></tr><tr><td>781 479</td><td>1.40</td><td>0.22</td><td>4.4</td><td>6.0</td><td>1 088 302</td><td>1.65</td><td>0.44</td><td>3.7</td><td>3.4</td></tr><tr><td>785 042</td><td>0.85</td><td>0.18</td><td>3.6</td><td>4.4</td><td>1 095 012</td><td>2.86</td><td>0.69</td><td>6.3</td><td>3.9</td></tr><tr><td>790 927</td><td>3.30</td><td>0.52</td><td>5.8</td><td>6.0</td><td>1 095 149</td><td>1.70</td><td>0.45</td><td>4.1</td><td>3.5</td></tr><tr><td>791 560</td><td>1.00</td><td>0.27</td><td>3.4</td><td>3.4</td><td>1 098 546</td><td>2.15</td><td>0.45</td><td>4.5</td><td>4.4</td></tr><tr><td>795 572</td><td>1.60</td><td>0.26</td><td>3.7</td><td>5.7</td><td>1 099 864</td><td>2.90</td><td>0.62</td><td>6.2</td><td>4.4</td></tr><tr><td>840 812</td><td>0.80</td><td>0.21</td><td>4.2</td><td>3.6</td><td>1 099 903</td><td>1.50</td><td>0.32</td><td>3.2</td><td>4.4</td></tr><tr><td>885 659</td><td>1.05</td><td>0.31</td><td>4.4</td><td>3.2</td><td>1 108 349</td><td>1.96</td><td>0.51</td><td>3.6</td><td>3.6</td></tr><tr><td>972 428</td><td>1.70</td><td>0.40</td><td>3.6</td><td>4.0</td><td>1 112 291</td><td>4.00</td><td>0.48</td><td>3.7</td><td>7.7</td></tr><tr><td>1 000 929</td><td>3.90</td><td>0.67</td><td>7.4</td><td>5.4</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Given are the center frequency in the velocity rest frame of SMA2, the line area, the peak intensity, the corresponding signal-to-noise ratio, and the line width (FWHM). (b) blended line.

The result is $N_H = 7.4 \times 10^{23}$ for a source size of $\theta = 32''$ and is lower compared to the value given above where $\theta = 10''$ is assumed. The following differences and similarities are evident when comparing NGC 6334I with an example of a

low-mass protostar IRAS 16293-2422 and the high-mass star-forming region Orion KL, Table 7: The methanol abundance $X = 5 \times 10^{-6}$ in NGC 6334I is higher than in the Orion HC and IRAS 16293, which have the same ratio. As is shown
<table><caption>**Table 7.** Comparison of the fractional abundances in NGC 6334I with IRAS 16293-2422 and the hot core (HC) and compact ridge (CR) of Orion KL.</caption><thead><tr><th>Molecule</th><th>NGC<br>6334I</th><th>IRAS<br>16293-2422<sup>c</sup></th><th>Orion KL<br>HC<sup>e</sup></th><th>Orion KL<br>CR<sup>e</sup></th></tr></thead><tbody><tr><td>CO</td><td>1.2(-5)</td><td>4.0(-5)<sup>d</sup></td><td>9.0(-5)</td><td>1.1(-4)</td></tr><tr><td>CH<sub>3</sub>OH</td><td>4.7(-6)</td><td>3.0(-7)</td><td>1.4(-7)</td><td>4.0(-7)</td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>1.0(-6)</td><td>2.4(-7)</td><td>8.0(-9)</td><td>1.9(-8)</td></tr><tr><td>HCN</td><td>8.7(-7)</td><td>-</td><td>3.0(-7)</td><td>-</td></tr><tr><td>H<sub>2</sub>O</td><td>6.7(-7)</td><td>-</td><td>5.3(-6)</td><td>1.0(-4)</td></tr><tr><td>OCS</td><td>4.0(-7)</td><td>2.5(-7)</td><td>1.1(-8)</td><td>3.0(-8)</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>2.4(-7)</td><td>2.0(-7)</td><td>1.4(-8)</td><td>3.0(-8)</td></tr><tr><td>CS</td><td>1.2(-7)</td><td>3.0(-9)<sup>d</sup></td><td>6.0(-9)</td><td>1.0(-8)</td></tr><tr><td>SO</td><td>8.3(-8)</td><td>2.5(-7)</td><td>1.9(-7)</td><td>3.0(-7)</td></tr><tr><td>H<sub>2</sub>S</td><td>8.3(-8)</td><td>9.0(-8)</td><td>2.7(-6)<sup>f</sup></td><td>-</td></tr><tr><td>H<sub>2</sub>CO</td><td>5.0(-8)</td><td>6.0(-8)</td><td>7.0(-9)</td><td>4.0(-8)</td></tr><tr><td>SO<sub>2</sub></td><td>3.3(-8)</td><td>1.0(-7)</td><td>1.2(-7)</td><td>1.6(-7)</td></tr><tr><td>HCOOH</td><td>1.8(-8)</td><td>6.2(-8)</td><td>8.0(-10)</td><td>1.4(-9)</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>1.4(-8)</td><td>1.2(-8)</td><td>3.0(-9)</td><td>5.0(-9)</td></tr><tr><td>H<sub>2</sub>CS</td><td>1.0(-8)</td><td>5.5(-9)</td><td>8.0(-10)</td><td>1.2(-9)</td></tr><tr><td>CH<sub>3</sub>CN</td><td>9.3(-9)</td><td>1.0(-8)</td><td>4.0(-9)</td><td>5.0(-9)</td></tr><tr><td>HNCO</td><td>5.0(-9)</td><td>9.0(-9)</td><td>5.8(-9)</td><td>-</td></tr><tr><td>C<sub>2</sub>H</td><td>1.4(-9)</td><td>2.2(-10)<sup>d</sup></td><td>5.5(-10)</td><td>5.3(-9)</td></tr><tr><td>HC<sub>3</sub>N</td><td>1.4(-9)</td><td>1.0(-9)</td><td>1.8(-9)</td><td>6.0(-9)</td></tr><tr><td>SiO</td><td>1.0(-9)</td><td>4.5(-9)</td><td>6.0(-9)</td><td>8.0(-9)</td></tr><tr><td>HCO<sup>+</sup></td><td>4.3(-10)</td><td>1.4(-9)<sup>d</sup></td><td>1.1(-9)</td><td>1.0(-9)</td></tr><tr><td>CN</td><td>3.7(-10)</td><td>8.0(-11)<sup>d</sup></td><td>8.0(-10)</td><td>8.0(-10)</td></tr></tbody></table>

**Notes.** Notation a(b) indicates a × 10⁻⁶. (d) Referring to the cooler, less dense envelope.

**References.** (c) Cazaux et al. (2003), Schöier et al. (2002). (e) Sutton et al. (1995), Caselli et al. (1993) and references therein. (f) From Persson et al. (2007).

by time-dependent chemical models of hot cores (Nomura & Millar 2004), the parent molecule methanol (formed on grain surfaces) is stable for timescales of at least 10⁴ yr. The other organic molecules, dimethyl ether and methyl formate, are more abundant in NGC 6334I than in Orion, reflecting perhaps that NGC 6334I is more time evolved. These second generation molecules start to form in the gas phase between evaporated molecules and peak up in abundance at 10⁴⁻⁵ yr. After that, the parent molecules are eventually destroyed. For sulfur-bearing molecules, an anticorrelation between NGC 6334I and Orion is observable: while H₂S is two orders of magnitude less abundant in NGC 6334I, OCS is much more abundant than in Orion. The higher depletion of H₂S could be explained by the fact that this molecule is efficiently destroyed by H₃O⁺ in hot cores and subsequently transformed into SO and SO₂ (Wakelam et al. 2004). Overall, the fractional abundance of the sulfur-bearing molecules and their sequence in NGC 6334I is more comparable to IRAS 16293 than to Orion (except for CS). Some nitrogen-bearing molecules like nitriles have about the same fractional abundances in NGC 6334I as in Orion and IRAS 16293: CH₃CN with X = 1 × 10⁻⁸, CN with X = 4 × 10⁻¹⁰, HC₃N with X = 2 × 10⁻⁹ and C₂H₅CN with X = 1 × 10⁻⁸. When comparing fractional abundances in the literature, it should be pointed out that possible uncertainties are the adopted H₂ column density and the source size. It needs to be further investigated if all these deviations are due to different precursor abundances, varying gas temperatures and densities or are due to age effects solely.

Results for most molecules are discussed briefly below, see Figs. B.1 and B.2 for their spectra. They are grouped by their moment of inertia, ranging from the lightest diatomic to the heaviest complex organic molecules. The electronic state for

molecules is denoted by 2S+1ΛΩ, where S is the total electronic spin, Λ the projection of the total electronic orbital angular momenta along the internuclear axis and Ω = Λ + Σ, where Σ is the projection of S onto the axis.

## 5.2. Atomic and diatomic species

CI and CII, Fig. B.1a. Two fine structure transitions of neutral carbon fall into the submillimeter range, ³P₁⁻³P₀ at 492 161 MHz and ³P₂⁻³P₁ at 809 342 MHz. Up to five velocity components can be observed, with the main emission at -7.3 km s⁻¹. The nondetection of ¹³C is used to set an upper limit for the abundance. The [CII] ²P₃/₂⁻²P₁/₂ fine-structure transition at 1.9 THz is one of the most important cooling lines in the ISM and one of the most intense single lines (∫ Tdv = 262 K km s⁻¹) in the HIFI spectrum besides CO. Three components are found, corresponding to the ones in CI. The large source sizes for all components indicate a broad distribution. The dust attenuation is neglected in the fit of the CII emission line because it is assumed that most of the emission originates from the diffuse outer layer. Therefore, the derived values for N and T should be taken as lower limits, and we notice that some of the spectral features could be affected or caused by contamination in the off-beam.

HF, Fig. B.1s. The ground transition J = 1–0 of hydrogen fluoride at 1 232 476 MHz is seen in absorption (Emprechtinger et al. 2012). Three components at -5, 0.3 and 7 km s⁻¹ are visible. The broad velocity width of the -5 km s⁻¹ component indicates that HF is widely distributed and extended over NGC 6334I. HF is considered to be a good tracer of H₂ (Neufeld et al. 2005), since the formation is directly coupled through the exothermic reaction F + H₂ → HF + H, and the ratio in the diffuse medium is X = 3.6 × 10⁻⁸. This would correspond to N<sub>H</sub> = 3 × 10²¹ cm⁻² in this diffuse layer.

SO, Fig. B.2d. Sulfur monoxide is a radical with two unpaired electrons (ground state ³Σ) and a spin-multiplicity of 3, so that each total angular momentum level J = N + S is split into 3 levels according to the projection of S. The notation is J,N → J - 1, N - 1, where N is the rotational momenta. Beginning from J = 12, the transitions are strong and can be observed up to J = 24. Additionally, the isotopologue ³³SO is present for J = 12 and ³⁴SO for J = 12 up to J = 16. Weaker SO transitions are present 4<sub>3</sub>-1<sub>2</sub> and 5<sub>4</sub>-2<sub>3</sub> in Band 1a and 7<sub>7</sub>-6<sub>7</sub> up to 12<sub>12</sub>-11<sub>12</sub>. The emission comes predominantly from SMA2, as revealed by the velocity of -7.4 km s⁻¹ and the SMA map. The abundance is derived from the isotopologues, but it was necessary to implement a second component which is optically thin.

SiO, Fig. B.1r. Silicon monoxide is viewed as a good tracer for outflows and the J = 5–4 transition was observed by Gibb et al. (2007). NGC 6334I shows a bipolar outflow with a red lobe in the northeast and a blue lobe in the southwest direction. With better resolution and additional transitions in the HIFI survey the abundance and spatial position can be given more securely. The J = 12 and J = 13 transitions in the HIFI data are quite weak. Therefore the parameters were only manually determined. Notable in the SMA data is the narrow -17 km s⁻¹ component of the blue lobe and the wide, more diffuse outflow in the red wing. Note that the massive outflows are better traced by CO maps of other groups. Leurini et al. (2006) analyzed the ratios of several CO transitions and derived 1–3 M<sub>☉</sub> per lobe (dependent on the inclination angle), a dynamical time scale (lobe length divided by v<sub>max</sub>) of 2–4 × 10³ yr and a total mass entrainment rate of M<sub>tot</sub> = 2–15 × 10⁻³ M<sub>☉</sub>/yr.
NH, Fig. B.1d. The ground state transitions $J = 1-0$ of the radical Imidogen with a spin multiplicity of 3 ($^3\Sigma$) and its hyperfine structure are observed. The coupling of the nuclear spin $I = 1$ of $^{14}$N with $J$ leads to the splitting into $F = J + I$ components. To fit it properly, a second broad component had to be included. The $J = 12-0_1$ transition is blended by the $J = 11-10$ transition of HCN.

SH$^+$, Fig. B.1v. Sulfoniumylidene has been detected only recently in Sgr B2(M) (Menten et al. 2011). The $12-0_1$ transition and its hyperfine structure are observed tentatively above $3\sigma$ in absorption and is heavily blended with methanol. The velocity matches that of SMA1. A comparative study of SH$^+$ absorption lines toward other high-mass star-forming regions has been made by Godard et al. (2012).

CH and CH$^+$, Fig. B.1c and B.1b. The hyperfine transitions of methylidine (CH) were analyzed by van der Wiel et al. (2010). We included the blended lines from methanol and find three components: one broad emission from NGC 6334I and two absorption lines from the foreground. The values are in very good agreement with their model I. The two lowest transitions of methylidinium (CH$^+$) are observed including the isotope $^{13}$C (Lis et al., in prep.). The three components are approximately the same as in CH and show the same abundance.

NO, Fig. B.1q. The radical nitrogen monoxide with its doublet structure ($^2\Pi$) is detected in the transitions $^{2}\Pi_{1/2}$, $J = 5.5-4.5$ (e and f) up to $^{2}\Pi_{1/2}$, $J = 8.5-7.5$ (e and f). The hyperfine structure is not resolved.

NS, Fig. B.2b. The electronic ground state of nitrogen sulfide is $^2\Pi$, which splits the rotation levels due to the unpaired electron. All transitions from $J = 11$ to 19 are detected. Noticeable is that the fit for the $\Pi_{1/2}$ states is good, but the $\Pi_{3/2}$ states are underexcited. The optical depth is small, so the source size was fixed at $\theta = 3''$.

CN, Fig. B.1l. Due to the weak emission and small number of transitions, there is a large uncertainty of the abundance depending on the velocity width. A width of 4 km s$^{-1}$ would result in an abundance higher by a factor of 2. Here the best fit is taken. The $J = 2_{0,3}-1_{0,2}$ transition at 226.8 GHz in the SMA data is seen in absorption and resembles well the dust emission. The assignment to SMA1 due to $v_{LSR} = -6.4$ km s$^{-1}$ is confirmed by the map.

CO and CS, Fig. B.1m and B.2a. Carbon monoxide is the most abundant molecule after $H_2$ in the ISM. In the HIFI spectrum, CO has the most intense lines, which are optically thick. The line profile shows strong outflows with widths of over 80 km s$^{-1}$ and self-absorption. The same is true for carbon monosulfide, so no direct abundance determination is possible. Therefore only the four detected isotopologues of CS and CO are considered, with a focus on $C^{18}\text{O}$ and $C^{17}\text{O}$ by Ceccarelli et al. (in prep.). For CS, a multi-component approach was tried with two components from the hot cores and one extended emission. But because the first and third component have the same offset velocity, there might be a redundancy.

### 5.3. Polyatomic linear molecules

$N_2H^+$, Fig. B.1g. Four emission lines from $J = 6$ to 10 of diazenylum are detected which exhibit an asymmetric line profile (Ceccarelli et al., in prep.). The intensities could not be fitted properly with a single component. The two component model shows a much better fit, but the values are very uncertain and a narrow line width had to be taken. Russeil et al. (2010) find two components also by observing the hyperfine structure of

the $10_1-0_{12}$ line, at -7.3 and -10 km s$^{-1}$, but with large differences between them concerning $T$ (127 K to 5 K) and $\Delta v$ (4 to 19 km s$^{-1}$). Another possibility could be self-absorption.

HNC and HCN, Fig. B.1o and B.1n. Hydrogen isocyanide is an isomer of hydrogen cyanide (HCN). Transitions from $J = 6$ to 10 are present for HNC and $J = 6$ to 12 for HCN. As for $N_2H^+$, a two component model was taken. The line profile is changing and affected by self absorption. HCN has a high critical density ($n_{cr} > 10^9$ for $J = 9-8$ at 797 GHz) and shows an outflow in the red wing of the asymmetric line profile with self-absorption, which indicates an infall. Also vibrationally excited lines ($v_2 = 1, 2$) of HCN are observed. Both aspects, infall and hot HCN, are characteristic for evolved cluster-forming regions (Rolffs et al. 2011). A uniform fit with the isotopologues of HCN is included for the completeness of the survey. For HNC, no other isotopologues are found apart from the $HN^{13}C J = 6-5$ line.

HCO$^+$, Fig. B.2p. The rotational transitions of the formyl ion are relatively intense and present up to band 5a (Ceccarelli et al., in prep.). At first, one component including $H^{13}CO^+$ was tried, but the fitting analysis revealed that either the lower-$J$ or the higher-$J$ transitions are consistent but not all together. This is an indication that either LTE is not satisfied, or several components with different excitation conditions overlap. Assuming two components with different optical depths, the fit could be improved. At higher transitions the component with the higher optical depth dominates.

OCS, Fig. B.2c. Carbonyl sulfide is detected in the rotational transitions from $J = 40$ to 53. Because at first the source size was unknown and no isotopologues are detected in the HIFI spectrum, the abundance was underestimated due to a low opacity ($\tau < 0.3$). With the detection of three isotopologues in the SMA bands the column density was corrected upwards. The emission of $O^{13}CS$ comes predominantly from SMA1.

$HC_3N$, Fig. B.2t. Cyanoacetylene lines are weak in the HIFI data, but more intense in the SMA data. The fit for SMA1 gives a large line width of 9 km s$^{-1}$, which is taken for the HIFI model.

### 5.4. Symmetric top molecules

$CH_3CN$, Fig. B.2m. Methyl cyanide is considered to be a good tracer of hot cores (Olmi et al. 1993). Rotational transitions from $J = 27$ to 39 with K-components $K = 0-6$ are detected in the HIFI spectrum. Thorwirth et al. (2003) fitted the $J = 12-11$ transition with an unknown source size resulting in $T = 145$ K and $N = 5.2 \times 10^{15}$ cm$^{-2}$. The temperature is in good agreement with our value of 154 K. Additionally, weak ($<3\sigma$) bending modes $v_8 = 1$ are observed in band 1a. The small source size and high temperature confirms that methyl cyanide probes the central part of hot cores.

$NH_3$, Fig. B.1f. The ground transitions of ammonia are detected mostly in absorption. The overlapping emission and absorption components could not be fitted consistently, and the column density for the emission component associated with NGC 6334I is therefore probably underestimated.

### 5.5. Inorganic asymmetric molecules

$NH_2$, Fig. B.1e. Azanyl is an asymmetric top with a hyperfine structure. Four transitions are detected including the ground state transition. The emission component can be associated with NGC 6334I. As in the case of NH, a broad component ($\Delta v = 30$ km s$^{-1}$) is added to fit the absorption.
$H_2Cl^+$ and HCl, Fig. B.1u and B.1t. Hydrogen chloride and chloronium were analyzed by Lis et al. (2010). The *J* = 2–1 at 1251 GHz in band 5b of HCl is additionally detected, showing a second component in absorption from the colder and more extended envelope at −10 km s⁻¹ with a broad width comparable to $H_2Cl^+$. The emission component can be assigned to the center of NGC 6334I. The abundance for this component differs by 50% compared to (Lis et al. 2010) because the source size is not known and the lines are optically thick and affected by the absorption feature. For $H_2Cl^+$, the blended lines of dimethylether are included in the fit (blue line). Its velocity and line width leads to the conclusion that it originates from a diffuse foreground molecular cloud.

$H_2O^+$, Fig. B.1j. Oxidaniumyl was analyzed by Ossenkopf et al. (2010). Additionally, another weak fine structure transition was observed at 1140 GHz from the *J* = 0.5 level. The derived abundance is above the lower limit given by Ossenkopf et al. (2010). Its radial velocity and line width leads to the conclusion that it is not associated with NGC 6334I.

$H_2O$, Fig. B.1i. Water was analyzed by Emprechtinger et al. (2010, 2012). It has complex line shapes with many components like optically thick emission and broad outflows, therefore no fit was attempted. Instead, the optically thin isotopes were included to get an estimation of the abundance. Most fits including deuterated water can not reproduce the line strength, so that a non-LTE calculation method is necessary to model them properly.

$H_2S$, Fig. B.1w. Hydrogen sulfide shows 3 components. The emission one at −6.1 km s⁻¹ is affected by two absorption lines, of which one is from the outflow. While the emission peak could be constrained through the isotopologues, the absorptions features could not be reproduced well.

## 5.6. Organic asymmetric molecules

$H_2CO$, Fig. B.2l. Formaldehyde is a nearly prolate top and the asymmetry splits every *K* ≠ 0 level into doublets. For astrophysical applications, due to its structure this molecule is a good probe of the kinetic temperature and volume density and traces the inner and outer part of hot cores (Mangum & Wootten 1993). Interstellar formaldehyde is formed by hydrogenation of solid CO: $H + CO \rightarrow HCO + H \rightarrow H_2CO$. Another less efficient pathway in the gas phase is $CH_3+O \rightarrow H_2CO + H$. Besides of the weeds, this organic molecule has the most intense lines throughout the HIFI spectrum from band 1a to 4b. The rotational transitions range from *J* = 7 to 15 with *K* = 0 to *K* = 6, and the isotopologues $H_2^{13}CO$ and $H_2^{18}O$ are observed. The fits with XCLASS show two problems: one component is optically thick, and all the *K* components could not be described by a single rotation temperature. The *K* = 0,1 do not match with *K* = 2–6 components consistently. From the SMA data the $3_0-2_0$ transition shows that the *K* = 0 value is lower in intensity than the *K* = 2 component. This hints to self absorption in SMA1 or optical thickness.

$H_2CCO$ and $H_2CS$, Fig. B.2o and B.2f. Ethenone and thioformaldehyde are nearly prolate tops with a spectrum similar to $H_2CO$. For $H_2CCO$, only the *K* = 1 components from *J* = 24 to 27 are detected, the other *K* = 0 and *K* = 3 components are included to set an upper limit. Since the transitions are weak and the optical depth is small, $\theta$ and $\Delta v$ needed to be fixed. For $H_2CS$ ranging from *J* = 14 to 23 with *K* = 0 to *K* = 5 transitions, a good fit is obtained where all *K* components can be described consistently with one excitation temperature.

$CH_3OH$, Fig. B.2g and B.2h. Methanol is the strongest weed in the HIFI spectrum and accounts for about half of all detected lines. One component is optically thick, which can be derived from the fit of the $^{13}CH_3OH$ isotopologue, and the other component accounts for the broad emission. The model fits the observation quite well, except of the most intense lines. As can be seen in the SMA spectra, the XCLASS LTE model cannot reproduce the methanol lines due to a high optical depth. Furthermore, high energy lines >500 K are underestimated, so that a 3rd component might be necessary.

$NH_2CHO$, Fig. B.2s. Formamide is only detected weakly in band 1a. The optical depth is small, and the source size is determined from the SMA map which shows that the emission comes mainly from SMA2.

$C_2H_4O$, Fig. B.2n. Ethylene oxide is surprisingly only detected in SMA2. As for formamide, many weak lines are observed. The source velocity $v_{LSR}$ in the HIFI data is consistent with the general source velocity of SMA2.

HNCO, Fig. B.2r. Isocyanic acid is almost linear, but asymmetric and shows relatively weak emission. The very high temperature and small source size may trace the hot core SMA2.

$CH_3OCH_3$, Fig. B.2i. Dimethyl ether has two methyl groups, and similarly to methanol, a splitting into 4 sublevels (AA,EE,AE,EA). It is, after methanol, the organic molecule with the most dense and intense lines in the entire HIFI spectrum. The laboratory measured spectrum from the $^{13}CH_3OCH_3$ isotopologue are taken from (Koerber et al. 2010), but are not yet publicly available. The SMA spectra show that the dominant emission arises from SMA1 and that they are weak associations with the outflows.

$CH_3OCHO$, Fig. B.2j. Methyl formate has the densest spectrum of all detected molecules in the lowest HIFI band 1a (Brouillet et al., in prep.) and the highest number of lines in both SMA bands. The SMA maps show an association with the outflows.

## 6. Conclusions

In this article, observational data from the HIFI instrument aboard the Herschel Space Observatory and from the Submillimeter Array are analyzed for the high-mass star-forming region NGC 6334I. We employ a LTE radiative transfer model to get insights into the physical structure of the source. The molecular spectral line survey gives an overview of the chemical inventory of the molecular species and their excitation, whereas the high-resolution spectrometer allows us tracing the kinematics of the different velocity components. The interferometer maps reveal the source morphology and help to distinguish the embedded cores and resolve their sizes. The results fit well with the overall picture of a high-mass young stellar object and hot molecular cores described by van Dishoeck & Blake (1998).

Two main software tools, XCLASS and MAGIX, are used to demonstrate that modeling and handling of the vast HIFI data sets is possible. These tools allow us to decompose the spectrum into contributions of individual molecules and their isotopologues. By assuming LTE condition and employing the LM fitting algorithm, a fast extraction of information is possible. To take line blends into account, a model for every molecule is produced, which is consistent over multiple HIFI bands. It should be clear that for other excitation conditions than LTE and for optically thick lines like CO this method reaches its limitations. Sophisticated 3-d radiative transfer codes such as
RADMC-3D<fnref n="3" /> and LIME (Brinch & Hogerheijde 2010) with the possibility of non-LTE calculations will be used in the future to model NGC 6334I into more detail, taking into account the geometry and non-homogeneous distribution of temperature and abundances. With the improved capabilities of ALMA (Atacama Large Millimeter Array), it will be possible to study the molecular chemistry and kinematics on sub-arcsecond scales (<1000 AU) and to follow the fragmentation of protostellar cores with much improved mass sensitivity.

The molecular survey revealed a variety of chemical species, including saturated organic molecules. The astrochemical question is: how can these abundances be explained and how are they related to each other? The next step is therefore a careful comparison of the column densities presented here with predictions of astrochemical models and detailed comparisons with the molecular abundances and abundance ratios in other CHESS targets.

**Acknowledgements.** This work was carried out within the Collaborative Research Council 956, funded by the Deutsche Forschungsgemeinschaft (DFG). C. Ceccarelli acknowledges the financial support from the French spatial agency CNES and the Agence Nationale pour la Recherche (ANR), France (project FORCOMS, contracts ANR-08-BLAN-022). Support for this work was provided by NASA through an award issued by JPL/Caltech. HIFI has been designed and built by a consortium of institutes and university departments from across Europe, Canada and the United States under the leadership of SRON Netherlands Institute for Space Research, Groningen, The Netherlands and with major contributions from Germany, France and the US. Consortium members are: Canada: CSA, U Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INA F, Osservatorio Astrofisico di Arcetri-INA F; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronómico Nacional (IGN), Centro de Astrobiología (CSIC-INTA). Sweden: Chalmers University of Technology – MC2, RSS & GARD; Onsala Space Observatory; Swedish National Space Board, Stockholm University – Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: Caltech, JPL, NHSC.

## References

Bernst, I., Schilke, P., Moeller, T., et al. 2011, in Astronomical Data Analysis Software and Systems XX, eds. I. N. Evans, A. Accomazzi, D. J. Mink, & A. H. Rots, ASP Conf. Ser., 442, 505

Beuther, H., Walsh, A. J., Thorwirth, S., et al. 2007, A&A, 466, 989

Brinch, C., & Hogerheijde, M. R. 2010, A&A, 523, A25

Brogan, C. L., Hunter, T. R., Cyganowski, C. J., et al. 2009, ApJ, 707, 1

Brooks, K. J., & Whiteoak, J. B. 2001, MNRAS, 320, 465

Caselli, P., Hasegawa, T. I., & Herbst, E. 1993, ApJ, 408, 548

Caux, E., Bottinelli, S., Vastel, C., & Glorian, J. M. 2011, in IAU Symp., 280, 120

Cazaux, S., Tielens, A. G. G. M., Ceccarelli, C., et al. 2003, ApJ, 593, L51

Ceccarelli, C., Bacmann, A., Boogert, A., et al. 2010, A&A, 521, L22

Cesaroni, R. 2005, in Massive Star Birth: A Crossroads of Astrophysics, eds. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley, IAU Symp., 227, 59

Comito, C., & Schilke, P. 2002, A&A, 395, 357

Comito, C., Schilke, P., Phillips, T. G., et al. 2005, ApJS, 156, 127

De Buizer, J. M., Radomski, J. T., Piña, R. K., & Telesco, C. M. 2002, ApJ, 580, 305

de Graauw, T., Helmich, F. P., Phillips, T. G., et al. 2010, A&A, 518, L6

Emerson, J. P., Jennings, R. E., & Moorwood, A. F. M. 1973, ApJ, 184, 401

Emprechtinger, M., Lis, D. C., Bell, T., et al. 2010, A&A, 521, L28

Emprechtinger, M., Lis, D. C., Rolffs, R., et al. 2012, ApJ, submitted

Fisher, J., Paciga, G., Xu, L.-H., et al. 2007, J. Molecul. Spectrosc., 245, 7

Garrod, R. T., Weaver, S. L. W., & Herbst, E. 2008, ApJ, 682, 283

Gibb, A. G., Davis, C. J., & Moore, T. J. T. 2007, MNRAS, 382, 1213

Godard, B., Falgarone, E., Gerin, M., et al. 2012, A&A, 540, A87

Gooch, R. 1996, in Astronomical Data Analysis Software and Systems V, eds. G. H. Jacoby, & J. Barnes, ASP Conf. Ser., 101, 80

Herbst, E. 1989, in The Physics and Chemistry of Interstellar Molecular Clouds – mm and Sub-mm Observations in Astrophysics, eds. G. Winnewisser, & J. T. Armstrong, Lecture Notes in Physics (Berlin: Springer Verlag), 331, 344

Herbst, E., & van Dishoeck, E. F. 2009, ARA&A, 47, 427

Herpin, F., Marseille, M., Wakelam, V., Bontemps, S., & Lis, D. C. 2009, A&A, 504, 853

Hunter, T. R., Brogan, C. L., Megeath, S. T., et al. 2006, ApJ, 649, 888

Kalinina, N. D., Sobolev, A. M., & Kalenskii, S. V. 2010, New Astron., 15, 590

Koerber, M., Endres, C. P., Lewen, F., et al. 2010, in 65th OSU Int. Symp. Mol. Spectrosc., FC02

Leurini, S., Schilke, P., Parise, B., et al. 2006, A&A, 454, L83

Lis, D. C., Pearson, J. C., Neufeld, D. A., et al. 2010, A&A, 521, L9

Mangum, J. G., & Wootten, A. 1993, ApJS, 89, 123

Maret, S., Hily-Blant, P., Pety, J., Bardeau, S., & Reynier, E. 2011, A&A, 526, A47

Matthews, H. E., McCutcheon, W. H., Kirk, H., White, G. J., & Cohen, M. 2008, AJ, 136, 2083

McCutcheon, W. H., Sandell, G., Matthews, H. E., et al. 2000, MNRAS, 316, 152

Menten, K. M., Wyrowski, F., Belloche, A., et al. 2011, A&A, 525, A77

Migenes, V., Horiuchi, S., Slysh, V. I., et al. 1999, ApJS, 123, 487

Milam, S. N., Savage, C., Brewster, M. A., Ziurys, L. M., & Wyckoff, S. 2005, ApJ, 634, 1126

Möller, T., Bernst, I., Panoglou, D., et al. 2012, A&A, submitted

Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215

Neufeld, D. A., Wolfire, M. G., & Schilke, P. 2005, ApJ, 628, 260

Nomura, H., & Millar, T. J. 2004, A&A, 414, 409

Olmi, L., Cesaroni, R., & Walmsley, C. M. 1993, A&A, 276, 489

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943

Ossenkopf, V., Müller, H. S. P., Lis, D. C., et al. 2010, A&A, 518, L111

Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, eds. Y. Mizumoto, K.-I. Morita, & M. Ohishi, ASP Conf. Ser., 434, 139

Persi, P., Tapia, M., Roth, M., Gómez, M., & Marenzi, A. R. 2005, in Massive Star Birth: A Crossroads of Astrophysics, eds. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley, IAU Symp., 227, 291

Persson, C. M., Olofsson, A. O. H., Koning, N., et al. 2007, A&A, 476, 807

Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spec. Radiat. Transf., 60, 883

Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1

Qiu, K., Wyrowski, F., Menten, K. M., et al. 2011, ApJ, 743, L25

Rolffs, R., Schilke, P., Wyrowski, F., et al. 2011, A&A, 527, A68

Russeil, D., Zavagno, A., Motte, F., et al. 2010, A&A, 515, A55

Russeil, D., Zavagno, A., Adami, C., et al. 2012, A&A, 538, A142

Sandell, G. 2000, A&A, 358, 242

Schilke, P., Comito, C., Thorwirth, S., et al. 2006, A&A, 454, L41

Schöier, F. L., Jørgensen, J. K., van Dishoeck, E. F., & Blake, G. A. 2002, A&A, 390, 1001

Stahler, S. W., & Palla, F. 2005, The Formation of Stars (Wiley-VCH)

Sutton, E. C., Peng, R., Danchi, W. C., et al. 1995, ApJS, 97, 455

Thorwirth, S., Winnewisser, G., Megeath, S. T., & Tieftrunk, A. R. 2003, in Galactic Star Formation Across the Stellar Mass Spectrum, eds. J. M. De Buizer, & N. S. van der Bliek, ASP Conf. Ser., 287, 257

van der Wiel, M. H. D., van der Tak, F. F. S., Lis, D. C., et al. 2010, A&A, 521, L43

van Dishoeck, E. F. 2001, in The Promise of the Herschel Space Observatory, eds. G. L. Pilbratt, J. Cernicharo, A. M. Heras, T. Prusti, & R. Harris, ESA SP, 460, 185

van Dishoeck, E. F., & Blake, G. A. 1998, ARA&A, 36, 317

Wakelam, V., Caselli, P., Ceccarelli, C., Herbst, E., & Castets, A. 2004, A&A, 422, 159

Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 301, 640

Walsh, A. J., Thorwirth, S., Beuther, H., & Burton, M. G. 2010, MNRAS, 404, 1396

Walsh, M. S., Xu, L.-H., Lees, R. M., et al. 2000, J. Molecul. Spectrosc., 204, 60

<footnote n="3">³ http://www.it.a.uni-heidelberg.de/~dullemond/software/radmc-3d/</footnote>
Appendix A: SMA spectra and tables

<table><caption>Table A.1. XCLASS fit results from SMA for NGC 6334I.</caption><thead><tr><th>Source</th><th>Species</th><th>Θ<br>(")</th><th>T<sub>ex</sub><br>(K)</th><th>N<br>(cm<sup>-2</sup>)</th><th>Δν<br>(km s<sup>-1</sup>)</th><th>ν<sub>LSR</sub><br>(km s<sup>-1</sup>)</th><th>Notes</th></tr></thead><tbody><tr><td>SMA1</td><td>CN</td><td>2</td><td>3</td><td>0.2 × 10<sup>15</sup></td><td>2.2</td><td>1.3</td><td></td></tr><tr><td>SMA1</td><td>CN</td><td>2</td><td>3</td><td>0.8 × 10<sup>15</sup></td><td>6.2</td><td>-4.6</td><td></td></tr><tr><td>SMA2</td><td>CN</td><td>2</td><td>3</td><td>0.7 × 10<sup>15</sup></td><td>6.9</td><td>-0.6</td><td></td></tr><tr><td>SMA2</td><td>CN</td><td>2</td><td>3</td><td>1.5 × 10<sup>15</sup></td><td>3.2</td><td>-4.5</td><td></td></tr><tr><td>SMA1</td><td>DCN</td><td>3.5</td><td>359</td><td>1.6 × 10<sup>15</sup></td><td>5.8</td><td>-4.9</td><td>2</td></tr><tr><td>SMA2</td><td>DCN</td><td>5.0</td><td>146</td><td>0.6 × 10<sup>15</sup></td><td>6</td><td>-7.8</td><td>2</td></tr><tr><td>SMA1</td><td>SiO</td><td>1.7</td><td>40</td><td>1.0 × 10<sup>15</sup></td><td>8</td><td>-3.2</td><td>2</td></tr><tr><td>SMA2</td><td>SiO</td><td>2.0</td><td>40</td><td>1.0 × 10<sup>15</sup></td><td>6</td><td>-5.2</td><td>2</td></tr><tr><td>SMA1</td><td>H<sub>2</sub>S</td><td>3</td><td>62</td><td>1.4 × 10<sup>17</sup></td><td>5</td><td>-4</td><td>2, 3</td></tr><tr><td>SMA2</td><td>H<sub>2</sub>S</td><td>3</td><td>147</td><td>0.4 × 10<sup>17</sup></td><td>6.5</td><td>-6</td><td>2, 3</td></tr><tr><td>SMA1</td><td>CS</td><td>2.4</td><td>40</td><td>1.6 × 10<sup>18</sup></td><td>5.1</td><td>-3.2</td><td>2, 3, 5</td></tr><tr><td>SMA2</td><td>CS</td><td>2.2</td><td>32</td><td>6.6 × 10<sup>17</sup></td><td>2.7</td><td>-6.4</td><td>2, 3, 5</td></tr><tr><td>SMA1</td><td>OCS</td><td>2.5</td><td>137</td><td>2.3 × 10<sup>18</sup></td><td>4</td><td>-5.6</td><td>1</td></tr><tr><td>SMA2</td><td>OCS</td><td>2.5</td><td>126</td><td>1.1 × 10<sup>18</sup></td><td>4</td><td>-7</td><td>1</td></tr><tr><td>SMA1</td><td>CH<sub>3</sub>OH</td><td>4</td><td>78</td><td>2.4 × 10<sup>19</sup></td><td>3.6</td><td>-5.2</td><td></td></tr><tr><td>SMA2</td><td>CH<sub>3</sub>OH</td><td>4</td><td>104</td><td>2.0 × 10<sup>19</sup></td><td>3.0</td><td>-7.5</td><td></td></tr><tr><td>SMA1</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>3.0</td><td>78</td><td>2.3 × 10<sup>18</sup></td><td>4.9</td><td>-5.6</td><td></td></tr><tr><td>SMA2</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>3.0</td><td>98</td><td>0.7 × 10<sup>18</sup></td><td>5.4</td><td>-7.7</td><td></td></tr><tr><td>SMA1</td><td>CH<sub>3</sub>OCHO</td><td>2.3</td><td>104</td><td>7.7 × 10<sup>17</sup></td><td>4.6</td><td>-6.1</td><td></td></tr><tr><td>SMA2</td><td>CH<sub>3</sub>OCHO</td><td>3</td><td>115</td><td>9.6 × 10<sup>17</sup></td><td>3.8</td><td>-7.6</td><td></td></tr><tr><td>SMA1</td><td>H<sub>2</sub>CO</td><td>3.0</td><td>104</td><td>1.6 × 10<sup>17</sup></td><td>5.7</td><td>-5.5</td><td></td></tr><tr><td>SMA2</td><td>H<sub>2</sub>CO</td><td>3.0</td><td>138</td><td>1.3 × 10<sup>17</sup></td><td>4.6</td><td>-7.5</td><td></td></tr><tr><td>SMA1</td><td>NH<sub>2</sub>CHO</td><td>2</td><td>25</td><td>4.7 × 10<sup>15</sup></td><td>5.2</td><td>-4.7</td><td>2, 3</td></tr><tr><td>SMA2</td><td>NH<sub>2</sub>CHO</td><td>2</td><td>45</td><td>5.0 × 10<sup>15</sup></td><td>4.0</td><td>-6.9</td><td>2, 3</td></tr><tr><td>SMA1</td><td>HC<sub>3</sub>N</td><td>2.8</td><td>208</td><td>4.2 × 10<sup>15</sup></td><td>9.4</td><td>-5.7</td><td></td></tr><tr><td>SMA2</td><td>HC<sub>3</sub>N</td><td>2.4</td><td>223</td><td>1.4 × 10<sup>15</sup></td><td>5.8</td><td>-6.6</td><td></td></tr><tr><td>SMA1</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>2.0</td><td>147</td><td>5.2 × 10<sup>15</sup></td><td>6.0</td><td>-6.3</td><td></td></tr><tr><td>SMA1</td><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>2.2</td><td>280</td><td>4.2 × 10<sup>16</sup></td><td>6.6</td><td>-5.5</td><td></td></tr><tr><td>SMA2</td><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>1.0</td><td>123</td><td>1.6 × 10<sup>16</sup></td><td>4.1</td><td>-6.7</td><td></td></tr><tr><td>SMA1</td><td>CH<sub>3</sub>CH<sub>3</sub>CO</td><td>2.6</td><td>73</td><td>2.1 × 10<sup>16</sup></td><td>5.6</td><td>-5.4</td><td></td></tr><tr><td>SMA2</td><td>CH<sub>3</sub>CH<sub>3</sub>CO</td><td>2.6</td><td>159</td><td>2.9 × 10<sup>16</sup></td><td>4.8</td><td>-7.4</td><td></td></tr><tr><td>SMA2</td><td>HCOOH</td><td>2.0</td><td>258</td><td>5.5 × 10<sup>16</sup></td><td>4.3</td><td>-9</td><td>2</td></tr><tr><td>SMA2</td><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>2.2</td><td>217</td><td>2.6 × 10<sup>16</sup></td><td>3.4</td><td>-8.2</td><td></td></tr><tr><td>SMA2</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>2.5</td><td>77</td><td>1.4 × 10<sup>17</sup></td><td>4.9</td><td>-7.9</td><td></td></tr><tr><td>SMA2</td><td>C<sub>2</sub>H<sub>4</sub>O</td><td>2.1</td><td>335</td><td>6.7 × 10<sup>16</sup></td><td>5</td><td>-6.8</td><td></td></tr><tr><td>SMA2</td><td>SO</td><td>2.5</td><td>70</td><td>1.7 × 10<sup>17</sup></td><td>4.2</td><td>-8.7</td><td>2, 6</td></tr></tbody></table>

Notes. (1) Including only isotopologues <sup>18</sup>OCS, O<sup>13</sup>CS, OC<sup>33</sup>S. (2) 1 or few lines. (3) Blended. (4) Including <sup>13</sup>CH<sub>3</sub>OH. (5) Including only <sup>13</sup>C<sup>34</sup>S. (6) Including only <sup>33</sup>SO.
<figure><img src="image_6.png" /><figcaption>**Fig. A.1.** LSB Spectra (in steps of 500 MHz) of the two hot cores SMA1 (upper) and SMA2 (lower) of NGC 6334L, observed with the Submillimeter Array. Overlaid in red are the XCLASS fits. For more detailed information about each identified spectral line, see Table A.2.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
<table><caption>Table A.2. List of identified lines in the SMA bands of NGC 6334I in order of increasing rest frequency.</caption><thead><tr><th>Molecule</th><th>Transition</th><th>LSB<br>Frequency (MHz)</th><th>Molecule</th><th>Transition</th><th>USB<br>Frequency (MHz)</th></tr></thead><tbody><tr><td>H<sub>2</sub>S</td><td>2 2 0–2 1 1</td><td>216 710.437</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>10 2 9 2–9 1 8 2</td><td>226 661.701</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>16 5 11–16 4 12</td><td>216 741.468</td><td>CN</td><td>2 0 2 2–1 0 1 1</td><td>226 679.311</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>16 6 11–16 5 12</td><td>216 741.712</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 3 17 0–19 4 16 1</td><td>226 692.48</td></tr><tr><td><sup>18</sup>OCS</td><td>19–18</td><td>216 753.488</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 4 17 0–19 3 16 1</td><td>226 692.48</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>32 6 27 1–32 5 28 1</td><td>216 758.969</td><td>CH<sub>3</sub>OCHO</td><td>20 2 19 1–19 2 18 1</td><td>226 713.06</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>32 6 27 0–32 5 28 0</td><td>216 809.014</td><td>CH<sub>3</sub>OCHO</td><td>20 2 19 0–19 2 18 0</td><td>226 718.688</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>18 2 16 2–17 2 15 2</td><td>216 830.197</td><td><sup>13</sup>CH<sub>3</sub>OH</td><td>12 5 8 0–13 4 9 0</td><td>226 716.392</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>18 2 16 0–17 2 15 0</td><td>216 838.889</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 4 17 1–19 3 16 1</td><td>226 784.639</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>19 4 16 0–18 4 15 1</td><td>216 881.337</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 17 1–19 4 16 1</td><td>226 784.639</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>19 3 16 0–18 4 15 1</td><td>216 881.337</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 4 17 1–19 3 16 2</td><td>226 784.67</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>19 4 16 0–18 3 15 1</td><td>216 881.337</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 17 1–19 4 16 2</td><td>226 784.67</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>19 3 16 0–18 3 15 1</td><td>216 881.337</td><td>CH<sub>3</sub>OCHO</td><td>20 1 19 2–19 1 18 2</td><td>226 773.13</td></tr><tr><td>CH<sub>3</sub>OH</td><td>5 1 0 0–4 2 0 0</td><td>216 945.521</td><td>CH<sub>3</sub>OCHO</td><td>20 1 19 0–19 1 18 0</td><td>226 778.786</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>17 3 14 5–16 3 13 5</td><td>216 958.834</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 4 17 0–19 3 16 1</td><td>226 832.055</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>20 0 20 2–19 1 19 1</td><td>216 962.989</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 17 0–19 4 16 1</td><td>226 832.055</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 4 16 1–18 3 15 1</td><td>216 974.426</td><td>CH<sub>3</sub>OCHO</td><td>20 2 19 1–19 1 18 2</td><td>226 856.872</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 3 16 1–18 4 15 1</td><td>216 974.426</td><td>CH<sub>3</sub>OCHO</td><td>20 2 19 0–19 1 18 0</td><td>226 862.257</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 4 16 1–18 3 15 2</td><td>216 974.464</td><td>CN</td><td>2 0 3 3–1 0 2 2</td><td>226 874.191</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 3 16 1–18 4 15 2</td><td>216 974.464</td><td>CN</td><td>2 0 3 4–1 0 2 3</td><td>226 874.781</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 4 16 0–18 3 15 1</td><td>217 022.509</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 4 17 0–19 3 16 0</td><td>226 879.393</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 3 16 0–18 4 15 1</td><td>217 022.509</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 17 0–19 4 16 0</td><td>226 879.393</td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OH</td><td>14 1 13 0–13 2 12 0</td><td>217 044.616</td><td>CN</td><td>2 0 3 2–1 0 2 2</td><td>226 887.42</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 4 16 0–18 3 15 0</td><td>217 070.504</td><td>CN</td><td>2 0 3 3–1 0 2 3</td><td>226 892.128</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>19 3 16 0–18 4 15 0</td><td>217 070.504</td><td>CH<sub>3</sub>OCHO</td><td>6 6 1 4–5 5 1 4</td><td>226 913.736</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>30 4 26 0–30 4 27 0</td><td>217 077.079</td><td>CH<sub>3</sub>OH</td><td>14–2 0 0–13 3 0 0</td><td>226 939.41</td></tr><tr><td>SiO</td><td>5–4</td><td>217 104.98</td><td>CH<sub>3</sub>OCHO</td><td>25 9 17 1–25 8 18 1</td><td>226 974.039</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>35 9 27 1–35 8 28 1</td><td>217 165.366</td><td>CH<sub>3</sub>OCHO</td><td>19 2 17 2–18 2 16 2</td><td>227 019.55</td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>22 4 19 3–22 3 20 3</td><td>217 189.668</td><td>CH<sub>3</sub>OCHO</td><td>19 2 17 0–18 2 16 0</td><td>227 028.121</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>32 9 24 0–32 8 25 0</td><td>217 215.847</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 3 17 1–19 4 16 1</td><td>227 032.073</td></tr><tr><td>DCN</td><td>3 2–2 2</td><td>217 238.631</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 4 17 1–19 3 16 1</td><td>227 032.073</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>30 4 26 0–30 3 27 0</td><td>217 266.549</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 4 17 1–19 3 16 2</td><td>227 043.496</td></tr><tr><td>CH<sub>3</sub>OH</td><td>6 1–1 1–7 2–1 1</td><td>217 299.205</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 3 17 1–19 4 16 2</td><td>227 043.496</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>17 4 13 3–16 4 12 3</td><td>217 312.626</td><td>CH<sub>3</sub>OH</td><td>21 1 0 0–21 0 0 0</td><td>227 094.747</td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OH</td><td>10 2 8 0–9 3 7 0</td><td>217 399.55</td><td>CH<sub>3</sub>OCHO</td><td>25 9 16 2–25 8 18 1</td><td>227 106.67</td></tr><tr><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>24 1 24 0–23 1 23 1</td><td>217 449.995</td><td>CH<sub>3</sub>OH</td><td>12–1 0 0–11 2 0 0</td><td>227 229.511</td></tr><tr><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>24 0 24 0–23 0 23 1</td><td>217 450.27</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>13 3 10 1–12 3 9 1</td><td>227 294.752</td></tr><tr><td>CH<sub>3</sub>OH</td><td>16 1 0 0–15 3 0 0</td><td>217 525.002</td><td><sup>13</sup>C<sup>34</sup>S</td><td>5–4</td><td>227 300.506</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>5 1 4 1–4 0 4 0</td><td>217 548.152</td><td>HC<sub>3</sub>N, v = 0</td><td>25–24</td><td>227 418.905</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>25 3 22 0–24 4 20 1</td><td>217 549.328</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>21 3 19 0–20 3 18 0</td><td>227 560.862</td></tr><tr><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>21 2 19 1–20 2 18 0</td><td>217 587.548</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>21 2 19 0–20 2 18 0</td><td>227 560.862</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>15 4 11–15 3 12</td><td>217 615.673</td><td>CH<sub>3</sub>OCHO</td><td>21 0 21 2–20 0 20 2</td><td>227 561.741</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>15 5 11–15 4 12</td><td>217 615.72</td><td>(CH<sub>2</sub>OH)<sub>2</sub></td><td>22 4 18 0–21 4 17 1</td><td>227 587.074</td></tr><tr><td>CH<sub>3</sub>OH</td><td>15 6 1 1–16 5 1 1</td><td>217 642.678</td><td>CH<sub>3</sub>OCHO</td><td>18 3 15 3–17 3 14 3</td><td>227 599.261</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 2 18 0–19 2 17 0</td><td>217 744.747</td><td>NH<sub>2</sub>CHO</td><td>11 0 11–10 0 10</td><td>227 606.176</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 3 18 0–19 3 17 0</td><td>217 744.747</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>18 5 13 2–18 4 14 2</td><td>227 606.079</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>5 3 3 2–4 2 2 2</td><td>217 803.689</td><td>C<sub>2</sub>H<sub>5</sub>CN, v = 1</td><td>26 1 26 1–25 1 25 1</td><td>227 624.735</td></tr><tr><td><sup>33</sup>SO</td><td>5 6 8–4 5 7</td><td>217 832.642</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>47 11 36 1–46 12 35 1</td><td>227 668.331</td></tr><tr><td>CH<sub>3</sub>OH</td><td>20 1 0 0–20 0 0 0</td><td>217 886.504</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>47 11 37 3–46 12 34 3</td><td>227 675.733</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 2 18 0–19 3 17 1</td><td>218 023.182</td><td>C<sub>2</sub>H<sub>4</sub>CN</td><td>25 3 22–24 3 21</td><td>227 780.972</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>20 3 18 0–19 2 17 1</td><td>218 023.182</td><td>CH<sub>3</sub>OH</td><td>16 1 1 0–15 2 1 0</td><td>227 814.528</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 18 1–19 2 17 2</td><td>218 091.448</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>21 3 19 0–20 3 18 1</td><td>227 836.273</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 2 18 1–19 3 17 2</td><td>218 091.448</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>21 2 19 0–20 2 18 1</td><td>227 836.273</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>17 4 13 5–16 4 12 5</td><td>218 108.438</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 7 18–23 7 17</td><td>227 897.607</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 2 18 0–19 3 17 1</td><td>218 127.207</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 7 17–23 7 16</td><td>227 897.608</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 18 0–19 2 17 1</td><td>218 127.207</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 1–20 3 18 1</td><td>227 903.922</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 3 18 0–19 2 17 0</td><td>218 162.929</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 1–20 3 18 2</td><td>227 903.953</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>20 2 18 0–19 3 17 0</td><td>218 162.929</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 2 19 1–20 2 18 2</td><td>227 903.953</td></tr><tr><td>O<sup>13</sup>CS</td><td>18–17</td><td>218 198.998</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 6 19–23 6 18</td><td>227 906.683</td></tr><tr><td>H<sub>2</sub>CO</td><td>3 0 3–2 0 2</td><td>218 222.192</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 6 18–23 6 17</td><td>227 906.709</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>31 9 23 0–31 8 24 0</td><td>218 259.594</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 2 19 0–20 2 18 1</td><td>227 939.374</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>17 3 14 2–16 3 13 2</td><td>218 280.9</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 2 19 0–20 3 18 1</td><td>227 939.374</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>17 3 14 0–16 3 13 0</td><td>218 297.89</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 0–20 2 18 1</td><td>227 939.374</td></tr><tr><td>HC<sub>3</sub>N, v = 0</td><td>24–23</td><td>218 324.723</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 0–20 3 18 1</td><td>227 939.374</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>14 3 11–14 2 12</td><td>218 335.626</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 0–20 3 18 1</td><td>227 939.374</td></tr></tbody></table>
<table><caption>Table A.2. continued.</caption><thead><tr><th>Molecule</th><th>Transition</th><th>LSB<br>Frequency (MHz)</th><th>Molecule</th><th>Transition</th><th>USB<br>Frequency (MHz)</th></tr></thead><tbody><tr><td>C<sub>2</sub>H<sub>4</sub>O</td><td>14 4 11–14 3 12</td><td>218 335.634</td><td>CH<sub>3</sub>OCHO</td><td>24 9 15 0–24 8 16 0</td><td>227 954.452</td></tr><tr><td>(CH<sub>3</sub>OH)<sub>2</sub></td><td>22 4 19 0–21 4 18 1</td><td>218 371.495</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 5 20–23 5 19</td><td>227 966.032</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>24 3 21–23 3 20</td><td>218 389.97</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 5 19–23 5 18</td><td>227 967.589</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 7 17–22 7 16</td><td>218 398.555</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 3 19 0–20 3 18 0</td><td>227 974.758</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 7 17–22 7 16</td><td>218 398.555</td><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 0</td><td>21 2 19 0–20 2 18 0</td><td>227 974.758</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 7 16–22 7 15</td><td>218 398.555</td><td>HC<sub>3</sub>N,v<sub>7</sub>=1</td><td>25 1–24–1</td><td>227 977.277</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 6 18–22 6 17</td><td>218 402.435</td><td>CH<sub>3</sub>OCHO</td><td>31 4 27 2–31 3 28 2</td><td>227 986.255</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 6 17–22 6 16</td><td>218 402.451</td><td>CH<sub>3</sub>OCHO</td><td>24 9 15 2–24 8 16 2</td><td>227 994.545</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>11 11 1 0–10 10 0 0</td><td>218 402.853</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>13 1 12 1–12 1 11 1</td><td>227 891.911</td></tr><tr><td>CH<sub>3</sub>CH<sub>3</sub>CO, v = 1</td><td>11 11 0 0–10 10 1 0</td><td>218 404.816</td><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>13 3 10 0–12 3 9 0</td><td>228 029.05</td></tr><tr><td>CH<sub>3</sub>OH</td><td>4 2 0 0–3 1 0 0</td><td>218 440.063</td><td>CH<sub>3</sub>OCHO</td><td>31 4 27 0–31 3 28 0</td><td>228 057.899</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 5 19–22 5 18</td><td>218 451.297</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 11 13–23 11 12</td><td>228 087.245</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 5 18–22 5 17</td><td>218 452.357</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 11 14–23 11 13</td><td>228 087.245</td></tr><tr><td>NH<sub>2</sub>CHO</td><td>10 1 9–9 1 8</td><td>218 459.653</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 4 21–23 4 20</td><td>228 104.614</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>5 3 2 2–4 2 3 2</td><td>218 461.226</td><td>OC<sup>33</sup>S</td><td>19–18</td><td>228 151.943</td></tr><tr><td>H<sub>2</sub>CO</td><td>3 2 2–2 2 1</td><td>218 475.632</td><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>24 4 20–23 4 19</td><td>228 160.305</td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>23 3 21 1–23 2 22 1</td><td>218 491.914</td><td>CH<sub>3</sub>OCHO</td><td>24 9 16 1–24 8 17 1</td><td>228 205.832</td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>23 3 21 0–23 2 22 0</td><td>218 494.39</td><td>CH<sub>3</sub>OCHO</td><td>18 3 15 5–17 3 14 5</td><td>228 211.291</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 10 13–22 10 12</td><td>218 519.997</td><td>CH<sub>3</sub>OCHO</td><td>24 9 16 0–24 8 17 0</td><td>228 270.535</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 10 14–22 10 13</td><td>218 519.997</td><td>HC<sub>3</sub>N,v<sub>7</sub>=1</td><td>25–1–24 1</td><td>228 303.174</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>OH</td><td>21 5 16 2–21 4 17 2</td><td>218 554.499</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>26 3 24 0–25 4 21 0</td><td>228 421.737</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 4 20–22 4 19</td><td>218 573.646</td><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>26 3 24 1–25 4 21 1</td><td>228 422.878</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>23 3 21–22 3 20</td><td>218 585.072</td><td>CH<sub>3</sub>OH</td><td>21 0 0 0–20 3 0 0</td><td>228 467.902</td></tr><tr><td>CH<sub>3</sub>OCHO</td><td>27 7 21 0–27 5 22 0</td><td>218 593.341</td><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>25 2 23–24 2 22</td><td>228 483.136</td></tr><tr><td></td><td></td><td></td><td>HCOOH</td><td>10 2 8–9 2 7</td><td>228 544.168</td></tr></tbody></table>
**Appendix B: HIFI fitted curves of all molecules**

Examples of selected, fitted spectra for each molecule. In red is the model, in blue are mostly annotations or blended lines. If only few transitions of a molecule are present, they are all shown in the plot. For weeds or molecules with many lines, no quantum numbers are given and only selected frequency ranges over different bands are presented. Heavy blended lines are excluded for simplicity. All intensities are given in units of Kelvin. The quantum numbers refer to upper state – lower state. For absorption lines, this sequence should be read reversed.

The entire observed and simulated spectrum for each HIFI band will be made available online in the CLASS format on the CHESS website (http://www-laog.obs.ujf-grenoble.fr/heberges/hs3f/), as soon as the HIPE 8 reprocessed data are available. Additionally, the modeled spectrum for each molecule (including isotopologues and vibrationally excited states) can be loaded. It should be noticed that for the simulated spectrum of only one molecule, the blending by other molecules is not taken into account. This results in a deviation for the fits of absorption lines when a nearby emission line contributes to the continuum. In XCLASS, at first all the emission components are calculated and then the absorption components.
<figure><img src="image_8.png" /><figcaption>Fig. B.1. XCLASS fits of chemical species in the HIFI spectrum of NGC 6334I.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>(e) NH₂</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(g) N₂H⁺</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(d) NH</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(f) NH₃</figcaption></figure>

Fig. B.1. continued.
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

Fig. B.1. continued.
<figure><img src="image_17.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Fig. B.1. continued.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Fig. B.2. XCLASS fits of chemical species in the HIFI spectrum of NGC 6334I. continued.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Fig. B.2. continued.</figcaption></figure>
<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /><figcaption>Fig. B.2. continued.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>Fig. B.2. continued.</figcaption></figure>
<figure><img src="image_25.png" /><figcaption>Fig. B.2. continued.</figcaption></figure>