# Investigating the hot molecular core, G10.47+0.03: A pit of nitrogen-bearing complex organic molecules

Suman Kumar Mondal, Wasim Iqbal, Prasanta Gorai, Bratati Bhat, Valentine Wakelam, Ankan Das

► To cite this version:

Suman Kumar Mondal, Wasim Iqbal, Prasanta Gorai, Bratati Bhat, Valentine Wakelam, et al.. Investigating the hot molecular core, G10.47+0.03: A pit of nitrogen-bearing complex organic molecules. Astronomy & Astrophysics - A&A, 2023, 669, pp.A71. 10.1051/0004-6361/202243802. hal-04292079

HAL Id: hal-04292079

https://hal.science/hal-04292079v1

Submitted on 18 Nov 2023

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 CC BY 4.0 - Attribution - International License
# Investigating the hot molecular core, G10.47+0.03: A pit of nitrogen-bearing complex organic molecules

Suman Kumar Mondal<fnref n="2" />, Wasim Iqbal<fnref n="1" />, Prasanta Gorai<fnref n="3" />,<fnref n="2" />, Bratati Bhat<fnref n="2" />, Valentine Wakelam<fnref n="4" />, and Ankan Das<fnref n="5" />

Received 17 April 2022 / Accepted 11 August 2022

## ABSTRACT

**Context.** Recent observations have shown that Nitrogen-bearing complex organic species are present in large quantities in star-forming regions. Thus, investigating the N-bearing species in a hot molecular core, such as G10.47+0.03, is crucial to understanding the molecular complexity in star-forming regions. They also allow us to investigate the chemical and physical processes that determine the many phases during the structural and chemical evolution of the source in star-forming regions.

**Aims.** The aim of this study is to investigate the spatial distribution and the chemical evolution states of N-bearing complex organic molecules in the hot core G10.47+0.03.

**Methods.** We used the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the hot molecular core G10.47+0.03. The extracted spectra were analyzed assuming local thermodynamic equilibrium (LTE). LTE methods are used to estimate the column density of observed species. Furthermore, robust methods such as Markov chain Monte Carlo (MCMC) and rotational diagram methods are implemented for molecules for which multiple transitions were identified to constrain the temperature and column density. Finally, we used the Nautilus gas-grain code to simulate the nitrogen chemistry in the hot molecular core. We carried out both 0D and 1D simulations of the source. We compared the simulated abundances with observational results.

**Results.** We report various transitions of nitrogen-bearing species (NH$_2$CN, HC$_3$N, HC$_5$N, C$_2$H$_3$CN, C$_2$H$_5$CN, and H$_2$NCH$_2$CN) together with some of their isotopologues and isomers. Besides this, we also report the identification of CH$_3$CCH and one of its isotopologues. We present detailed chemical simulation results to investigate the possible N-bearing chemistry in the source.

**Conclusions.** In this study, various transitions of nitrogen-bearing molecules are identified and discussed. The emissions originating from vinyl cyanide, ethyl cyanide, cyanoacetylene, and cyanamide are compact, which could be explained by our astrochemical modeling. Our 0D model shows that the chemistry of certain N-bearing molecules can be very sensitive to initial local conditions such as density or dust temperature. In our 1D model, simulated higher abundances of species such as HCN, HC$_3$N, and HC$_5$N toward the inner shells of the source confirm the observational findings.

**Key words.** astrochemistry – submillimeter: ISM – line: profiles – methods: data analysis – methods: observational – ISM: abundances

## 1. Introduction

Numerous molecules (about 250) are identified in the interstellar medium (ISM) or circumstellar shells (CDMS<fnref n="1" />, McGuire 2018). Among them, 96 molecules contain at least one nitrogen atom. Of these, 33 have cyanide (–CN) and 12 isocyanide (–NC). CN was the first nitrogen-bearing molecule that was observed in 1940 (McKellar 1940). It is the key precursor of cyanopolynes, some hydrocarbons, and bio-molecules. Isomerization can play a key role in the conversion of cyanide to isocyanide and vice versa. The study of nitrogen-bearing complex molecules with CN bonds is of prime importance in prebiotic chemistry. Protein synthesis occurs through a peptide bond formation; for example, CN bonds formed between –COOH and –NH$_2$ groups from different amino acids (Goldman et al. 2010). Nitrogen-bearing molecules play a crucial role in the synthesis of biomolecules (Balucani 2009). Complex organic nitrogen-bearing molecules were detected in various regions of our galaxy, such as star-forming regions, comets, and

atmospheres of planets and the Moon (Johnson et al. 1977; Irvine & Schloerb 1984; Kawaguchi et al. 1992; Jiménez-Serra et al. 2016; Ziurys et al. 1999). Some nitrogen-containing ions such as HC$_3$NH$^+$, NH$_4^+$, H$_2$NCO$^+$, NCCNH$^+$, CN$^-$, C$_3$N$^-$, C$_5$N$^-$, and HCNH$^+$ have been identified in the star forming regions (Kawaguchi et al. 1994; Cernicharo et al. 2013, 2008; Marcelino et al. 2018; Agúndez et al. 2015, 2010; Thaddeus et al. 2008; Ziurys & Turner 1986).

The study of complex organics in star-forming region can shed some light on the region’s evolutionary history. Dense, hot, and chemically enriched inner regions of high-mass star-forming regions are hot molecular cores (HMCs). Many simple and complex species were identified in these regions (van Dishoeck & Blake 1998; Belloche et al. 2019; Jørgensen et al. 2020). The detection of complex organic molecules is an important tool for explaining the physical and chemical structure of high-mass star-forming regions. G10.47+0.03 (hereafter G10) is a massive hot molecular core located at a distance of 8.6 kpc (Sanna et al. 2014). The luminosity of this source is ~7 × 10$^5$ L$_\odot$ (Cesaroni et al. 2010).

<footnote n="1">¹ South-Western Institute for Astronomy Research (SWIFAR), Yunnan University (YNU), Kunming 650500, PR China
e-mail: wasim@ynu.edu.cn</footnote>

<footnote n="2">² Indian Centre for Space Physics, 43 Chalantika, Garia Station Road, Kolkata 700084, India</footnote>

<footnote n="3">³ Department of Space, Earth & Environment, Chalmers University of Technology, 412 96 Gothenburg, Sweden</footnote>

<footnote n="4">⁴ Laboratoire d'astrophysique de Bordeaux, CNRS, Univ. Bordeaux, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="5">⁵ Institute of Astronomy Space and Earth Science, AJ 316, Salt lake Sector II, Kolkata 700091, West Bengal, India
e-mail: ankan.das@gmail.com</footnote>

<footnote n="1">¹ https://cdms.astro.uni-koeln.de/classic/molecules</footnote>
Interstellar dust is used as a catalyst for the synthesis of com-
plex organic molecules. For example, vinyl cyanide (C₂H₃CN)
was primarily detected toward Sgr B2 in 1975 (Gardner &
Winnewisser 1975). Later, this molecule was identified in var-
ious regions of the ISM, such as cold cores, Orion KL, the cir-
cumstellar envelope of the carbon-rich AGB star CW Leo, which
is also known as IRC+10216, and so on (Agúndez et al. 2008;
López et al. 2014). The observed abundance of vinyl cyanide
in IRC+10216 is two times lower than in Orion KL. Agúndez
et al. (2008) found a good agreement between their modeling
and observational results for IRC+10216. Ethyl cyanide was first
detected in the hot core OMC-1 and Sgr B2 (Johnson et al. 1977).
Later it was also observed in different hot cores (G327.3-0.6,
Orion IRC2, Orion KL, Sgr B2(N), Sgr B2 (M)) and in hot cori-
nos (IRAS 16293-2422). Along with the main isotopic species,
singly and doubly ¹³C-substituted ethyl cyanide was identified in
the ISM (Margulès et al. 2016). Cyanoallene was observed with
the 100 m Green Bank Telescope in the cold core of TMC by
Lovas et al. (2006). Cyanamide is a potential prebiotic molecule
(Ibanez et al. 1971) that was first detected several decades ago
in the hot core, Sgr B2 (Turner et al. 1975). To date, the most
complex cyanide molecule observed in the ISM is the propyl
cyanide (n-C₃H₇CN) and its isomer, first branch chain molecule,
isopropyl cyanide (i-C₃H₇CN; Belloche et al. 2009, 2014).

Rolffs et al. (2011b) observed G10 in three different fre-
quency ranges (200 GHz, 345 GHz, and 690 GHz) with the sub-
millimeter array (SMA). They identified many oxygen-bearing,
nitrogen-bearing, and sulfur-bearing species in this source. Fur-
thermore, many complex organic molecules and some highly
excited lines of HC₃N, HCN, and HC¹⁵N observed in this
source reveal that it is an active site for chemical enrich-
ment (Rolffs et al. 2011a,b). Here, we analyze the Atacama
Large Millimeter/submillimeter Array (ALMA) archival data
of G10 in the ALMA band 4 frequency region. This paper
focuses on nitrogen-bearing molecules, their distributions, com-
plex chemistry, and fractional abundances. Various nitrogen-
related molecules such as HCN, CH₃CN, CH₃NH₂, CH₂NH,
C₂H₃CN, C₂H₅CN, HC₃N, NH₂CN, CH₃NCO, NH₂CHO,
and HNCO have been observed (Rolffs et al. 2011b; Fontani et al.
2007; Ohishi et al. 2019; Suzuki et al. 2016; Gorai et al. 2020;
Wyrowski et al. 1999) in G10. In this work, we detect some new
nitrogen-bearing molecules such as HC₅N, H₂NCH₂CN, some
isotopologues of vinyl cyanide, ethyl cyanide, and HC₃N.

This paper is organized as follows. First, in Sect. 2, observa-
tional details are given. Astrochemical modeling and its results
are discussed in Sects. 3 and 4, respectively. Finally, in Sect. 5,
we conclude.

## 2. Observational data and analysis

Here, ALMA Cycle 4 data of G10 is used for the analysis. The phase center of the observation is α (J2000) = 18h08m38s232s, δ (J2000) = -19°51′50.4″. The systematic velocity of the source is 68 km s⁻¹. The flux calibrator was J1733-1304, the phase calibrator was J1832-2039, and the bandpass calibrator was J1924-2914. This observation was performed with the four spectral bands, covering the sky frequencies of 129.50–131.44 GHz, 147.50–149.43 GHz, 153.00–154.93 GHz, and 158.49–160.43 GHz and corresponding angular resolution at four different frequencies are 1.5" (12 900 AU), 1.72" (14 792 AU), 1.66" (14 276 AU), and 1.76" (15 136 AU), considering the distance to the source is 8.6 kpc. Summaries of the observational details (frequency range, channel width, baseline, etc.) are already presented in Gorai et al. (2020). For data analysis and imaging we used CASA

(version 4.7.2; McMullin et al. 2007), and for the line identifica-
tion and analysis we used the CASSIS software<fnref n="2" /> (developed by
IRAP-UPS/CNRS).

2.1. Images

Cesaroni et al. (2010) resolved the hot core G10 and found three distinct HII regions inside the core, which are A (Ultra Compact HII region), B1, and B2 (Hyper compact HII regions). Rolffs et al. (2011b) also observed this source with the submillimeter Array (SMA). They observed the continuum at three different frequencies (201/211 GHz, 345/355 GHz, 681/691 GHz). The beam sizes were insufficient to resolve the continuum emission at 201/211 and 681/691 GHz. However, at 345/355 GHz, the beam size was sufficient (~0.3") to resolve this source. Here, continuum maps of G10 were observed at four different wavelengths, 1.88 mm (159.45 GHz), 1.94 mm (153.96 GHz), 2.02 mm (148.51 GHz), and 2.3 mm (130.50 GHz), respectively. Summaries of continuum images (synthesized beam, peak flux, integrated flux, etc.) are discussed in detail in Gorai et al. (2020). Gorai et al. (2020) calculated the H₂ column density for each continuum map and obtained an average value of ~1.35 × 10²⁵ cm⁻<fnref n="2" />. Gorai et al. (2020) also estimated the dust opacity of each continuum map and found all are optically thin (τ_v ~ 0.135).

2.2. Analysis

All lines are extracted from the 2.0" (17 200 AU) diameter
region centered at G10 (α(J2000)=18h08m38.232s, δ(J2000)=
-19°51′50.4"). Line identification has been carried out with CAS-
SIS together with the spectroscopic database, ‘Cologne Database
for Molecular Spectroscopy’ (CDMS, Müller et al. 2001, 2005)<fnref n="3" />,
and the Jet Propulsion Laboratory (JPL, Pickett et al. 1998)<fnref n="4" />
database. The observed transitions of various nitrogen-bearing
molecules together with their quantum numbers, upper state
energies (E_u), V_LSR, line parameters such as peak brightness
temperature (I_max), line width (full width at half maximum;
FWHM), and the integrated intensity (∫T_mbdV) are noted in
Table A.1. The line parameters (I_max, FWHM, ∫T_mbdV) are
obtained using a single Gaussian fit to the observed spectral pro-
file of each unblended transition. To see the spatial distribution
of the observed species, moment zero maps of each transition
are shown in Figs. B.1–B.4. The emitting regions of all the tran-
sitions are calculated using a two-dimensional Gaussian fitting
of their moment maps.

2.3. Rotational diagram analysis

Among all the observed nitrogen-bearing species toward G10,
we only find multiple unblended transitions (more than 2)
for four molecules (see Table A.1), which are vinyl cyanide
(C₂H₃CN), ¹³C isotopologue of vinyl cyanide (¹³CH₂CHCN),
ethyl cyanide (C₂H₅CN), and ¹³C isotopologue of ethyl cyanide
(¹³CH₃CH₂CN). In addition to these four nitrogen-bearing
molecules, multiple transitions are observed for methyl acety-
lene (CH₃CCH). Thus, a rotational diagram analysis was carried
out for these five species to measure their excitation temperatures
and the total column densities.

<footnote n="2">² http://cassis.irap.omp.eu (Vastel et al. 2015)</footnote>

<footnote n="3">³ https://cdms.astro.uni-koeln.de/</footnote>

<footnote n="4">⁴ http://spec.jpl.nasa.gov/</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. Rotational diagram of (a) vinyl cyanide (C₂H₃CN) and (b) its one isotopologue (¹³CH₂CHCN). The best fitted rotational temperature and column density are given in each panel. Filled blue squares are the data points and the blue line represents the error bar. The error bar of each transition is calculated from rms noise (varies 117 mK to 306 mK) and calibration error (considering 5%).</figcaption></figure>

For the optically thin transitions, the column density ($N_u^{thin}$) at local upper state can be expressed as

$$ \frac{N_u^{thin}}{g_u} = \frac{3k_B \int T_{mb} dV}{8\pi^3 v S \mu^2}, \qquad (1) $$

where $g_u$ is the degeneracy of the upper state, $k_B$ is the Boltzmann constant, $\int T_{mb} dV$ is the integrated intensity, $v$ is the rest frequency, $\mu$ is the electric dipole moment, and $S$ is the transition line strength (Goldsmith & Langer 1999). Under the local thermodynamic equilibrium (LTE) condition and for optically thin transitions, the total column density ($N_{total}$) of a molecular species can be prescribed as

$$ \frac{N_u^{thin}}{g_u} = \frac{N_{total}}{Q(T_{rot})} \exp(-E_u/k_B T_{rot}), \qquad (2) $$

where $T_{rot}$ is the rotational temperature, $E_u$ is the upper state energy, $Q(T_{rot})$ is the partition function at rotational temperature. This can be rearranged as

$$ \ln\left(\frac{N_u^{thin}}{g_u}\right) = -\left(\frac{1}{T_{rot}}\right)\left(\frac{E_u}{k}\right) + \ln\left(\frac{N_{total}}{Q(T_{rot})}\right). \qquad (3) $$

There is a linear relationship between the upper state energy and column density at the upper level when all lines are optically thin and all lines are probing the same physical conditions. The column density and rotational temperature are extracted from the rotational diagram. Rotational diagrams of all the five species are shown in Figs. 1, 3, and 6.

## 2.4. Markov chain Monte Carlo model

Though the rotational diagram is used for the species for which multiple transitions are observed, another method, Markov chain Monte Carlo (MCMC), is used to constrain the excitation temperature and column density for the comparison. The MCMC algorithm was initially developed for dealing with the probabilistic behavior of the collection of atomic particles. However, it was challenging to do this analytically. Thus, an iterative process was implemented to solve it through simulation. This stochastic model describes a sequence of possible events in which the

probability of each event depends only on the state attained in the previous one. The MCMC method is an interactive process that goes through all line parameters (e.g., molecular column density, excitation temperature, source size, line width) with a random walk and heads into the solution's space. $\chi^2$ minimization gives the final solution.

For MCMC fitting, the FWHM of all the transitions of a species are kept constant at their average measured value, $v_{LSR}$ is kept constant at 68 km s⁻¹, and source size is kept constant at 2”. Then, the obtained results of these species are compared with those obtained with the rotation diagram method. Obtained results with these two methods are presented in Table 1 and discussed in Sect. 3. The obtained line parameters of the MCMC fitting of five species are noted in Table A.2. MCMC fitted spectra are shown in Figs. C.1–C.3.

## 3. Column densities and spatial distributions of all molecules

The column densities of all species are estimated by generating the synthetic spectra under an LTE approximation. The blended transitions are also identified by carefully examining the synthetic spectra of two nearby transitions. As the input parameters of the LTE model, the source size of ~2” (half-maximum diameter of our beam size), an excitation temperature of 150 K, and molecular hydrogen column density $N(H_2) \sim 1.35 \times 10^{25} \text{ cm}^{-2}$ (Gorai et al. 2020) are considered. The FWHM of a species is chosen by taking an average FWHM of the detected, unblended, optically thin transitions. The column density of a species is varied to obtain a satisfactory fit with the observed data. Figures 2, 4, 5, and 7 show the LTE-fit spectra (see captions of these figures for more details). Table 1 provides the best-fit column densities of all species. We also generated another set of LTE spectra (called LTE2) by keeping all other parameters the same but considering a different source size. For this fitting, the measured emitting regions of an unblended transition and the average value of emitting regions for multiple transitions are used as an input source size. These LTE spectra (LTE2) are shown in Appendices D.1–D.4, and the obtained column densities are noted in Table 1. We find that the estimated column density slightly increased by using measured emitting regions as
<table><caption>**Table 1.** Estimated temperature and column density of the observed species.</caption><thead><tr><th rowspan="2">Species</th><th colspan="4">Temperature (K)</th><th colspan="3">Column density (cm<sup>-2</sup>)</th></tr><tr><th>Other observation</th><th>LTE/LTE2</th><th>RD</th><th>MCMC</th><th>Other observation</th><th>LTE/LTE2<sup>(*)</sup></th><th>RD</th><th>MCMC</th></tr></thead><tbody><tr><td colspan="9" style="text-align:center; font-weight:bold;">Vinyl cyanide and its isotopologues</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>CN</td><td>200<sup>(a)</sup><br>176 ± 35<sup>(b)</sup></td><td>150</td><td>155 ± 53</td><td>209.82 ± 18.02</td><td>7.0 × 10<sup>17</sup><sup>(a)</sup><br>2.0 × 10<sup>14</sup>/6.1 × 10<sup>16</sup></td><td>2.9 × 10<sup>17</sup>/6.0 × 10<sup>17</sup>(1.07)</td><td>(1.7 ± 0.4) × 10<sup>17</sup></td><td>(3.4 ± 0.44) × 10<sup>17</sup></td></tr><tr><td><sup>13</sup>CH<sub>2</sub>CHCN</td><td>—</td><td>150</td><td>119 ± 18</td><td>125.86 ± 13.67</td><td>—</td><td>6.5 × 10<sup>15</sup>/9.60 × 10<sup>15</sup>(1.07)</td><td>(3.5 ± 0.4) × 10<sup>15</sup></td><td>(5.3 ± 0.46) × 10<sup>15</sup></td></tr><tr><td>CH<sub>3</sub><sup>13</sup>CHCN</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>2.8 × 10<sup>15</sup>/5.30 × 10<sup>15</sup>(0.78)</td><td>—</td><td>—</td></tr><tr><td>CH<sub>2</sub>CH<sup>13</sup>CN</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>3.0 × 10<sup>15</sup>/5.60 × 10<sup>15</sup>(0.70)</td><td>—</td><td>—</td></tr><tr><td>CH<sub>2</sub>CHC<sup>15</sup>N</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>9.0 × 10<sup>14</sup></td><td>—</td><td>—</td></tr><tr><td>CH<sub>2</sub>CDCN</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>3.3 × 10<sup>15</sup>/1.00 × 10<sup>16</sup>(0.67)</td><td>—</td><td>—</td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Ethyl cyanide with its isotopologues and isomer</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>CN</td><td>200<sup>(a)</sup><br>103 ± 12<sup>(b)</sup></td><td>150</td><td>92 ± 18</td><td>104.11 ± 9.98</td><td>9 × 10<sup>17</sup><sup>a</sup><br>3.6 × 10<sup>14</sup>/1.1 × 10<sup>17</sup></td><td>2.3 × 10<sup>17</sup>/3.80 × 10<sup>17</sup>(1.35)</td><td>(1.2 ± 0.3) × 10<sup>17</sup></td><td>(1.7 ± 0.10) × 10<sup>17</sup></td></tr><tr><td><sup>13</sup>CH<sub>3</sub>CH<sub>2</sub>CN</td><td>—</td><td>150</td><td>220 ± 112</td><td>130.79 ± 23.17</td><td>—</td><td>8.5 × 10<sup>15</sup>/1.5 × 10<sup>16</sup>(1.16)</td><td>(9.7 ± 2.3) × 10<sup>15</sup></td><td>(7.8 ± 1.12) × 10<sup>15</sup></td></tr><tr><td>CH<sub>3</sub><sup>13</sup>CH<sub>2</sub>CN</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>8.0 × 10<sup>15</sup></td><td>—</td><td>—</td></tr><tr><td>CH<sub>3</sub>CH<sub>2</sub><sup>13</sup>CN</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>3.2 × 10<sup>15</sup>/4.70 × 10<sup>15</sup>(1.22)</td><td>—</td><td>—</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>C<sup>15</sup>N</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>3.1 × 10<sup>15</sup></td><td>—</td><td>—</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>NC</td><td>—</td><td>150</td><td>—</td><td>—</td><td>—</td><td>4.0 × 10<sup>15</sup>/6.40 × 10<sup>15</sup>(1.05)</td><td>—</td><td>—</td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Cyanoacetylene and its isotopologues</td></tr><tr><td>HC<sub>3</sub>N</td><td>200<sup>(a)</sup>/500<sup>(a)</sup></td><td>150</td><td></td><td></td><td>5.0 × 10<sup>16</sup><sup>a</sup>/10<sup>18</sup><sup>a</sup></td><td>4.1 × 10<sup>15</sup>/4.0 × 10<sup>15</sup>(2.18)</td><td></td><td></td></tr><tr><td>HC<sup>13</sup>CCN</td><td></td><td>150</td><td></td><td></td><td></td><td>1.4 × 10<sup>15</sup>/2.30 × 10<sup>15</sup>(1.17)</td><td></td><td></td></tr><tr><td>HCC<sup>13</sup>CN</td><td></td><td>150</td><td></td><td></td><td></td><td>1.7 × 10<sup>15</sup>/2.80 × 10<sup>15</sup>(1.16)</td><td></td><td></td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Cyanodiacetylene and its isotopologues</td></tr><tr><td>HC<sub>5</sub>N</td><td></td><td>150.0</td><td></td><td></td><td></td><td>7.7 × 10<sup>14</sup>/1.70 × 10<sup>15</sup>(1.20)</td><td></td><td></td></tr><tr><td>HCC<sup>13</sup>CCCN</td><td></td><td>150.0</td><td></td><td></td><td></td><td>1.0 × 10<sup>15</sup>/2.35 × 10<sup>15</sup>(0.95)</td><td></td><td></td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Cyanamide</td></tr><tr><td>NH<sub>2</sub>CN</td><td>200<sup>(a)</sup></td><td>150.0</td><td></td><td></td><td>2.0 × 10<sup>16</sup><sup>(a)</sup></td><td>4.8 × 10<sup>15</sup>/1.30 × 10<sup>16</sup>(0.86)</td><td></td><td></td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Aminoacetonitrile</td></tr><tr><td>H<sub>2</sub>NCH<sub>2</sub>CN</td><td></td><td>150.0</td><td></td><td></td><td></td><td>5.5 × 10<sup>15</sup>/9.30 × 10<sup>15</sup>(1.07)</td><td></td><td></td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Methylacetylene and its isotopologues</td></tr><tr><td>CH<sub>3</sub>CCH</td><td></td><td>150.0</td><td>188.9 ± 16</td><td>219.12 ± 14.80</td><td></td><td>7.0 × 10<sup>16</sup>/8.30 × 10<sup>16</sup>(1.67)</td><td>(6.0 ± 0.5) × 10<sup>16</sup></td><td>(8.7 ± 0.75) × 10<sup>16</sup></td></tr><tr><td>CH<sub>3</sub><sup>13</sup>CCH</td><td></td><td>150</td><td></td><td></td><td></td><td>8.0 × 10<sup>15</sup>/9.80 × 10<sup>15</sup>(1.17)</td><td></td><td></td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">Other nitrogen-bearing molecules</td></tr><tr><td>HCN</td><td>200<sup>(a)</sup>/500<sup>(a)</sup></td><td></td><td></td><td></td><td>1.0 × 10<sup>18</sup>/3.0 × 10<sup>18</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>NH<sub>2</sub></td><td>64<sup>(c)</sup></td><td></td><td></td><td></td><td>6.3 × 10<sup>15</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>2</sub>NH</td><td>84/84<sup>(d)</sup></td><td></td><td></td><td></td><td>2.1 × 10<sup>17</sup>/4.7 × 10<sup>15</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>CN</td><td>408/82<sup>(e)</sup><br>240/180<sup>(f)</sup></td><td></td><td></td><td></td><td>5.1 × 10<sup>17</sup>/4.1 × 10<sup>14</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>NCO</td><td>248 ± 19<sup>(g)</sup></td><td></td><td></td><td></td><td>6.0 × 10<sup>16</sup>/3.6 × 10<sup>15</sup></td><td></td><td></td><td></td></tr><tr><td>HNCO</td><td>317 ± 25<sup>(g)</sup></td><td></td><td></td><td></td><td>1.2 × 10<sup>17</sup></td><td></td><td></td><td></td></tr><tr><td>NH<sub>2</sub>CHO</td><td>439 ± 100<sup>(g)</sup></td><td></td><td></td><td></td><td>3.8 × 10<sup>16</sup></td><td></td><td></td><td></td></tr><tr><td></td><td>200<sup>(a)</sup></td><td></td><td></td><td></td><td>2.0 × 10<sup>17</sup></td><td></td><td></td><td></td></tr></tbody></table>

Notes. (a) Rolffs et al. (2011b), consider source size 1.5", H₂ column density 5.0 × 10²⁴ cm⁻², and an additional component of 0.5" and 500 K, which are used for HC₃N and HCN; (b) Fontani et al. (2007), low value column density is beam averaged (23") and high value column density is source averaged (1.3"), consider H₂ column density 1.2 × 10²⁵ cm⁻²; (c) Ohishi et al. (2019), considering source size 10" and H₂ column density 2.6 × 10²³ cm⁻²; (d) Suzuki et al. (2016), for high column density, considering source size 1.4" and H₂ column density 6.69 × 10²⁴ cm⁻² and low column density obtained, considering source size 10" and H₂ column density 1.3 × 10²³ cm⁻²; (e) Hernández-Hernández et al. (2014), a high rotational temperature and column density are for compact component ~0.6" and low for extended components ~4.0", consider H₂ column density 6.69 × 10²⁴ cm⁻²; (f) Olmi et al. (1996), consider two component structure:a dense compact core ~1"–2" and the halo extended about 10 times larger than core; (g) Gorai et al. (2020), using same observation data; (h) The average source size (in arcsec) that value is used for LTE modeling (LTE2) is noted in parentheses.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Black lines represent the observed spectra of vinyl cyanide and its isotopologues ($^{13}$CH$_2$CHCN, CH$_2$$^{13}$CCN, CH$_2$CH$^{13}$CN, CH$_2$CNC$^{15}$N, and CH$_2$CDCN). Unblended transitions are noted in red, whereas blended transitions are indicated in black. The solid red lines represent the Gaussian fit of the unblended transitions. LTE spectra are shown in blue by considering N(H$_2$) = 1.35 × 10$^{25}$ cm$^{-2}$, excitation temperature 150 K, source size of 2'', and the average FWHM (obtained from the unblended optically thin transitions) of the species. The column density of the species is varied until eye-estimated good fits are obtained. The red dashed lines show the systematic velocity (V$_{LSR}$) of the source at ~68 km s$^{-1}$. In addition, the name of the species and their respective transitions (in GHz) are provided in each panel.</figcaption></figure>

source size instead of constant value (~2''). The estimated values of column densities and temperatures with the rotational diagram, MCMC, and LTE methods are summarized in Table 1. All the estimated parameters agree well among themselves. However, the rotational diagram and MCMC method provide more realistic values because the MCMC method is an alternative method to the rotational diagram. Therefore, for comparison and to better constrain the estimated parameters, we employed the MCMC approach. In the LTE model, we have only one free parameter: column density. However, if all other fixed parameters (e.g., FWHM, source size, excitation temperature) are known with certainty, the simple LTE model also provides an outstanding result. Here, the rotational diagram method uses excitation temperature measurements for five molecules. Therefore, we know the gas's excitation temperature. Furthermore, the FWHM and source size of each molecule are also known from our measurement. To test the reliability of the simple LTE method, we apply the LTE model to five species for which we made a rotational diagram. Interestingly, the obtained column densities with the LTE method are comparable to that obtained with the rotational diagram (see Table 1). Therefore, unless otherwise stated, for the rest of the paper in discussions or comparison with modeling results, we use measured values of column densities obtained using a rotational diagram. For all other molecules (if

the result with the rotational diagram method is unavailable), column densities obtained with the LTE model are used. We discuss individual species in the following sections.

## 3.1. Vinyl cyanide and its isotopologues

### 3.1.1. C$_2$H$_3$CN

Altogether, twelve transitions of vinyl cyanide are identified toward G10, which are noted in Table A.1 and shown in Fig. 2. Among them, nine are unblended transitions, and three are blended (147.86523 GHz, 148.028678 GHz, and 159.753858 GHz). A rotational diagram (see Fig. 1a) is performed by considering only the unblended, optically thin transitions. It yields a rotational temperature of 155 ± 53 K and column density of (1.7 ± 0.4) × 10$^{17}$ cm$^{-2}$, respectively (see Table 1). Fontani et al. (2007) found a rotational temperature of vinyl cyanide of 176 ± 35 K, which is in good agreement with our measured value. Comparing our results with Rolffs et al. (2011b), we obtain a similar temperature, while their column density is about two to four times larger than ours. The optically thin unblended transitions of vinyl cyanide have an emitting region of about 0.72''–1.11'' (see Fig. B.1 and Table 2). The 131.168 GHz and 158.773 GHz transitions are found to be
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Rotational diagram of (a) ethyl cyanide ($C_2H_5CN$) and (b) its one isotopologue ($^{13}CH_3CH_2CN$). The best-fit rotational temperature and column density are given in each panel. Filled blue squares are the data points and blue line represents the error bar. Error bar of each transition is calculated from rms noise (varying from 70 mK to 155 mK) and calibration error (considering 5%).</figcaption></figure>

very compact (0.72–0.77"), whereas the other transitions are comparatively little extended (1.04–1.11"). Though these transitions are marginally resolved, our analysis indicates that vinyl cyanide is possibly emitting from the inner regions of the source.

### 3.1.2. $^{13}CH_2CHCN$

Nine transitions of the $^{13}CH_2CHCN$ (isotopologue of vinyl cyanide) are observed toward G10 (see Table A.1). Among them, five are blended. The four unblended transitions are optically thin, and a column density of (3.5 ± 0.4) × 10$^{15}$ cm$^{-2}$ and rotational temperature of 119 ± 18 K are obtained from rotational diagram. The 159.954 GHz transition of $^{13}CH_2CHCN$ is found (see panels 10–13 of Fig. B.1, and Table 2) to be very compact (~0.73"), whereas the other transitions are slightly extended (1.16–1.20").

### 3.1.3. $CH_2^{13}CHCN$ and $CH_2CH^{13}CN$

Another two $^{13}C$ isotopologues of vinyl cyanide ($CH_2^{13}CHCN$ and $CH_2CH^{13}CN$) are tentatively observed toward G10. Only one transition of each of these isotopologues is clearly identified (see Fig. 2). Column densities of 2.8 × 10$^{15}$ cm$^{-2}$ and 3.0 × 10$^{15}$ cm$^{-2}$ for $CH_2^{13}CHCN$ and $CH_2CH^{13}CN$, respectively, depict a best LTE fit to the observed data (see solid blue lines of panels 22–29 of Fig. 2). Both the transitions are found to be very compact (0.70–0.78"; see panels 14 and 15 of Fig. B.1, and Table 2).

### 3.1.4. $CH_2CHC^{15}N$

Our simple LTE model found two possible transitions of $^{15}N$ isotopologue of vinyl cyanide ($CH_2CHC^{15}N$) (see panels 30–31 of Fig. 2). However, the transition at 148.898 GHz is partially contaminated with ethylene glycol, and the other transition at 159.436 GHz is blended with other species. So, no clear transitions are detected toward G10. We derive an upper limits column density of ~9.0 × 10$^{14}$ cm$^{-2}$.

### 3.1.5. $CH_2CDCN$

Altogether, seven transitions of a deuterated form of vinyl cyanide ($CH_2CDCN$) are observed toward G10 (Table A.1). Among them, five are blended. The 130.645 GHz transition is possibly blended with propynal and $C_2H_5C^{15}N$. 130.661 GHz is probably blended with its other transitions and $CH_3NCO$. 130.725 GHz is probably blended with $CH_3OCH_3$. 130.637 GHz is probably blended with its other transitions, and 154.309 GHz is probably blended with $CH_3OCH_3$. Since only two transitions are well detected, we consider that deuterated vinyl cyanide is tentatively detected from our observation. The best LTE fit (see panels 32–38 of Fig. 2) yields a column density of 3.3 × 10$^{15}$ cm$^{-2}$. The emitting region of these two transitions (see Table 2, and panels 11–12 of Fig. B.1) depict that they are comparatively more compact (0.64–0.69") than vinyl cyanide and other isotopologues of vinyl cyanide.

### 3.1.6. Fractionation of vinyl cyanide

The results obtained with the rotational diagram method, MCMC method, and simple LTE estimation reflect the $^{12}C/^{13}C$ isotopic ratio of vinyl cyanide, which varies from 45–64 (see Table 1). Despite a high elemental atomic H/D ratio in the ISM (Linsky 1998, 6.667 × 10$^4$), we see a substantial deuterium enrichment (~88 with simple LTE estimation) for vinyl cyanide. However, the $^{12}C/^{13}C$ isotopic ratio (Clayton & Nittler 2004, ~89) is roughly conserved in vinyl cyanide (45–64, with all the methods implemented).

## 3.2. Ethyl cyanide with its isotopologues and isomer

### 3.2.1. $C_2H_5CN$

In total, nineteen transitions of ethyl cyanide are identified toward G10 (see Table A.1). Among them, four are blended. Within the fifteen unblended transitions, six are optically thick, and nine are optically thin. From rotational diagram analysis (see Fig. 3a), we obtained a column density of 1.2 × 10$^{17}$ cm$^{-2}$ and rotational temperature of 92 K. Our derived rotational temperature is similar to the rotational temperature (103 ± 12 K) of ethyl cyanide derived by Fontani et al. (2007).
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Gaussian-fit and modeled LTE spectra (same as Fig. 2) of the unblended transitions of ethyl cyanide, its isotopologues (¹³CH₃CH₂CN, CH₃¹³CH₂CN, CH₃CH₂¹³CN, C₂H₅C¹⁵N), and one of its isomers (C₂H₅NC).</figcaption></figure>

Rolffs et al. (2011b) estimated a higher column density of ethyl cyanide (9 × 10¹⁷ cm⁻²) than vinyl cyanide (7.0 × 10¹⁷ cm⁻²). However, with all our methods (simple LTE, rotational diagram, and MCMC), a higher column density of vinyl cyanide is obtained compared to ethyl cyanide (see Table 1). Fontani et al. (2007) found vinyl cyanide is ~10 times more abundant than ethyl cyanide. However, when they considered all the transitions (strong, weak and vibrationally excited lines), they obtained the best fit parameters: column density ~1.7 × 10¹⁷ cm⁻², rotational temperature ~140 K, source size ~1”. The emitting regions of the optically thin transitions vary from 1.08”–1.48” (see panels 1–9 of Fig. B.2, and Table 2).

### 3.2.2. ¹³CH₃CH₂CN

Thirteen transitions of the ¹³C isotopologue of ethyl cyanide are observed toward this source (Table A.1). Among them, five transitions are unblended. The rotational diagram yielded a rotational temperature and a column density of 220 K and 9.7 × 10¹⁵ cm⁻², respectively (see Fig. 3b). The emitting region of the unblended and optically thin transitions of this isotopologue is found (see Fig. B.2) to be more compact (0.93–1.29”) than the transitions of ethyl cyanide.

### 3.2.3. CH₃¹³CH₂CN

Two transitions of this ¹³C isotopologues are detected toward G10 (Table A.1) but both transitions are found to be blended.

A column density upper limit of 8.0 × 10¹⁵ cm⁻² is estimated by the LTE approximation.

### 3.2.4. CH₃CH₂¹³CN

Five transitions of this species are observed toward G10 (Table A.1). Our LTE analysis indicates that three of them are blended. Since only two transitions are clearly detected, we consider this tentatively detected molecule with a column density of 2.9 × 10¹⁵ cm⁻². The emitting region’s unblended transitions of CH₃CH₂¹³CN vary from 1.35–1.22” (panels 15 and 16 of Fig. B.2).

### 3.2.5. C₂H₅C¹⁵N

Any transitions of a ¹⁵N isotopologue of ethyl cyanide (C₂H₅C¹⁵N) are not identified toward G10 with full confidence. Our LTE model found only two possible transitions (see panels 40–41 of Fig. 4), with an upper-limit column density of ~3.1 × 10¹⁵ cm⁻².

### 3.2.6. CH₃CH₂NC

CH₃CH₂NC is tentatively detected toward G10. Only two unblended transitions of the isomer of ethyl cyanide are observed. Our best-fit LTE spectra of these transitions yield a column density of 4.0 × 10¹⁵ cm⁻² (blue lines in panels 38–39 of Fig. 4). The emitting region of the transition at 159.061 GHz is
found to be very compact (0.75"), whereas the other transition at 149.12 GHz is found to be comparatively not very extended (1.34").

3.2.7. Fractionation of ethyl cyanide

Multiple transitions of ethyl cyanide and one of its ¹³C iso-
topologues (¹³CH₃CH₂CN) are identified. Based on the vari-
ous methods (simple LTE, rotational diagram, and MCMC), a
¹²C/¹³C ratio of 12–27 is estimated. Two transitions of another
¹³C isotopologue (CH₃CH₂¹³CN) of ethyl cyanide are identi-
fied. With the simple LTE method, a carbon fractionation of
72 is estimated. From the simple LTE fitting, an isomeric ratio
between the cyanide (C₂H₅CN) and isocyanide (C₂H₅NC) of
57.5 is estimated.

3.3. Cyanopolyynes

Cyanopolyynes are a group of species with the chemical for-
mula HC_nN (n = 3, 5, 7...). In addition, it contains the cyano
group (CN). Previously, Rolffs et al. (2011b) identified HCN,
HNC, H¹³CN, and HC₃N in G10. Here, the identifications of
cyanoactylene (HC₃N) and its isotopologues and cyanodiacety-
lene (HC₅N) and its isotopologue are reported.

3.3.1. Cyanoacetylene and its isotopologue

HC₃N is observed toward G10 with one transition at 154.657 GHz. Additionally, transitions of its two ¹³C isotopologues, HC¹³CCN and HCC¹³CN are detected at 154.001 GHz and 154.016 GHz, respectively. We obtained best-fit synthetic spectra with column densities of 4.1 × 10¹⁵ cm⁻², 1.4 × 10¹⁵ cm⁻², and 1.7 × 10¹⁵ cm⁻², for HC₃N, HC¹³CCN, and HCC¹³CN, respectively (see blue lines in panels 1–3 of Fig. 7). The transition of HC₃N is comparatively more extended (~2.18") than the other species reported in this article (see Fig. B.3). The isotopologues of HC₃N are found to be relatively compact (1.16"–1.17"). Rolffs et al. (2011b) obtained a column density of 5.0 × 10¹⁶ cm⁻² when they used an excitation temperature of 200 K. Additionally, with the 500 K and 0.5" they estimated a very high column density (~10¹⁸ cm⁻²) of HC₃N. Here, with our simple LTE analysis, we estimate a column density of HC₃N is 4.1 × 10¹⁵ cm⁻². This mismatch is possibly due to our low spatial resolution observation compared to that of Rolffs et al. (2011b). With the simple LTE calculation, a ¹²C/¹³C fractionation is found to be 2.41–2.93.

3.3.2. Cyanodyacetylene

We clearly detected one transition of HC₅N and one transi-
tion of its isotopologue (HCC¹³CCCN). Both the transitions
yield best LTE fits with the column densities of 7.7 × 10¹⁴ cm⁻²
and 1.0 × 10¹⁵ cm⁻², respectively. The spatial distribution (see
Fig. B.3 and Table 2) of HC₅N is found to be more compact
(1.2") than HC₃N (~2.18"). Furthermore, the HC₅N line we
observed has an upper state energy of ~156 K. In contrast, for
HC₃N, it is ~67 K. The observed spatial distribution of HC₃N
suggests that it traces the envelopes regions, whereas HC₅N may
originate from the warm inner part of the hot molecular core. The
emitting areas (see Table 2) of the isotopologue (¹³C) of HC₅N
depicts that it is more compact (0.95") than their parent isotopo-
logues. The same compactness is also observed for the observed
isotopologues (¹³C) of HC₃N. With the simple LTE calculation,
a ¹²C/¹³C fractionation of 0.77 is obtained.

<figure><img src="image_6.png" /><figcaption>Fig. 5. Gaussian-fit spectra and LTE spectra (same as Fig. 2) of CH₃CCH and its isotopologues.</figcaption></figure>

3.4. Cyanamide

Six transitions of cyanamide (NH₂CN) are observed toward G10 (Table A.1). However, two of them are clearly detected. Our simple LTE model yields a good fit (Fig. 7) to the observed spectra with a column density of 4.8 × 10¹⁵ cm⁻². The emitting region (see Fig. B.3 and Table 2) of the two unblended transitions are compact (~0.91").

3.5. Aminoacetonitrile

Aminoacetonitrile (H₂NCH₂CN) is a well-known precursor of glycine. It hydrolyses to form glycine. Danger et al. (2011) studied the NH₂CH₂CN formation under the interstellar condition by thermal processing of ices containing methanimine, ammonia, and hydrogen cyanide. Five transitions of aminoacetonitrile are observed toward G10 (Table A.1), but four of them are blended with other species. Only the transition at 154.51748 GHz with upper state energy of 86 K is successfully identified. Therefore, we consider that this species was tentatively detected in this observation. A column density of 5.5 × 10¹⁵ cm⁻² yields a good fit to the observed spectra (see panel 16 of Fig. 7). The emitting diameter of aminoacetonitrile is 1.07". The moment zero map of this transition is shown in panel 8 of Fig. B.3.

3.6. Methyl acetylene

Snyder & Buhl (1973) first identified methyl acetylene (CH₃CCH) in Sgr B2. Though methyl acetylene is not a nitrogen-bearing species, being a symmetric top molecule, its identification is reported here because of its utility as a tracer of the physical condition of the star-forming region. Three transitions of methyl acetylene are clearly identified. All the transitions correspond to *J* = 9–8. The three transitions have *k*<sub>a</sub> = 6, 3, and 2, respectively. This is the first identification of CH₃CCH in G10. All three unblended transitions are found to be optically thin. Rotational diagram analysis of the CH₃CCH (Fig. 6) depicts a column density of 6.0 × 10¹⁶ cm⁻² and a rotational temperature of 189 K. Panels 1–3 of Fig. B.4 show moment map of CH₃CCH for different *k*<sub>a</sub> components of *J* = 9–8 transitions. Interestingly, it shows that the emitting regions of CH₃CCH transitions are
<figure><img src="image_7.png" /><figcaption>Fig. 6. Rotational diagram of methyl acetylene (CH₃CCH; same as Figs. 1 and 3). Error bar of each transition is calculated from rms noise (varies 70 mK to 136 mK) and calibration error (considering 5%).</figcaption></figure>

related to their corresponding $k_a$ components. For example, the emitting region of $J = 9-8$ transitions with $k_a = 2$ is 1.8'' ($E_u = 65.8$ K), and for $k_a = 3$ it is 1.9'' ($E_u = 101.7$ K). Since with our spatial resolution these two transitions are marginally resolved, we did not find significant differences for close-lying energy states. However, for $k_a = 6$ ($E_u = 296.9$ K), a compact emission (1.28'') is obtained. With the increase in the $k_a$ component, the upper state energy increases, which yields a more compact emitting region. All three detected transitions have $\Delta k_a = 0$. The ratio of lines between the same $J$-state but other $k_a$-states can be used as temperature probe (Mangum & Wootten 1993). To estimate the kinetic temperature ($T_K$), we chose some selected line ratios ($\frac{J_1k_a - J_2k_a}{J_3k'_a - J_4k'_a}$). These ratios are selected based on the following criteria (Mangum & Wootten 1993; Das et al. 2019): 1) $\Delta J = 1$, $J_1=J_3$, and $J_2=J_4$, 2) $\Delta k_a = 0$ for both the transitions, 3) $k'_a > k_a$, 4) both the transitions are optically thin, 5) The ratio between the difference in excitation temperature belongs to different $k_a$ values.

With the LTE approximation, the line ratio ($R$) between the two transitions satisfying the above criteria can be expressed as (Mangum & Wootten 1993)

$$R = S_R \exp(D/T_K), \qquad (4)$$

where $D = E(J3, k'_a) - E(J1, k_a)$ and $S_R = \frac{S_{J1k_a}}{S_{J3k'_a}}$. Mangum & Wootten (1993) derived this relation for H₂CO, which is also applicable to any slightly asymmetric or symmetric rotor species.

According to the above selection rules, in between the three observed unblended transitions, $\frac{9_2-8_2}{9_3-8_3}$, $\frac{9_2-8_2}{9_6-8_6}$, and $\frac{9_3-8_3}{9_6-8_6}$ are considered for estimating the excitation temperature. The black solid curve in Fig. 8a-c depicts the expected line ratio with the excitation temperature. The observed line ratio of these three transitions is shown with the solid red line. Interestingly, the observed line ratio intersects the curve between 133 and 193 K, closer to the estimated hot core temperature.

Additionally, two transitions ($J = 9-8$ with $k_a = 0, 3$) of ¹³C isotopologue of methyl acetylene (CH₃¹CCH) are observed toward G10. Our LTE synthetic spectra fit the observed spectra well, with a column density of $8.0 \times 10^{15}$ cm⁻². A similar trend of emitting regions is also obtained here. For $k_a = 0$, the emitting region is 1.31'', whereas it is 1.01'' for $k_a = 3$. With

the simple LTE calculation, a carbon fractionation ($\frac{12C}{13C}$ ratio) of 8.8 is obtained between CH₃CCH and CH₃¹CCH. By considering the column density obtained with the rotational diagram and MCMC methods (for CH₃CCH), a $\frac{12C}{13C}$ ratio of 7.5–10.9 is obtained.

A similar process for estimating the kinetic temperature from the transitions obtained with the CH₃¹CCH is carried out. Following the selection rule, $\frac{9_0-8_0}{9_3-8_3}$ is used to estimate the kinetic temperature. Figure 8d depicts that the observed line ratio intersects the expected line ratio at 102 K.

We identify several transitions of nitrogen-bearing species along with their various ¹³C isotopologues. Only the unblended transitions are used for further analysis. In the case of the ¹³C isotopologue of vinyl cyanide, we identify four transitions of ¹³CH₂CHCN, one of CH₂¹CHCN, and one of CH₂CH¹³CN. We obtain two transitions from one of the ¹³C isotopologues of ethyl cyanide (CH₃CH₂¹CCH). Two transitions of the ¹³C isotopologue of cyanoacetylene (HC¹³CCN and HCC₁₃CN) and one transition of the ¹³C isotopologue of cyanodiacetylene (HCC¹³CCCN) are identified. Furthermore, we identify two transitions of a ¹³C isotopologue of methyl acetylene (CH₃¹CCH). We use the rotational diagram analysis and MCMC method to estimate the column density of these species if more than two transitions were identified. Otherwise, we consider the simple LTE approximation. Despite all these issues, we found that the obtained $\frac{12C}{13C}$ ratio of all these species varies from 0.8 to 72. In the case of methyl acetylene, vinyl, and ethyl cyanide, it is 7.5–72, which is generally consistent with the elemental $\frac{12C}{13C}$ ratio (~89). For the cyanopolyynes (cyanoacetylene and cyanodiacetylene), we obtain a small ¹³C enrichment (0.77–2.93). Since not many transitions were obtained for the cyanopolyynes, we restrict our comments based on the upper limit of the column densities obtained.

## 4. Modeling

Modeling the time evolution of chemical composition of high mass sources such as G10 requires knowledge of the history of their physical structure's evolution over the whole hot core formation stage, including the early diffuse stage, collapse stage, and the warm-up stage. However, the physical evolution during the initial stages of massive hot cores and its timescales are still poorly known, and the difficult task of reconstructing the evolutionary history of the physical structure of massive hot sources such as G10 is beyond the scope of this work. Nevertheless, to better understand the factors that can have a major impact on the chemical history of the hot cores such as G10, we did extensive chemical modeling.

In the present work, following Garrod (2013), we adopted a three-stage physical model to simulate the time-dependent chemical evolution of N-bearing species detected in G10. The three stages describing the evolution of physical structure of G10-like sources are 1) an isothermal cloud collapse stage, 2) a relatively short warm-up stage, post cloud collapse, and 3) a hot and dense post-warm-up stage. Here, the three-stage model is only explained briefly as this was already explained in detail in Sil et al. (2018); Das et al. (2019); Gorai et al. (2020); Bhat et al. (2021). In the first stage, the cloud collapses iso-thermally from a diffuse state (initial H density, $n_H = 10^3$ cm⁻³) to a dense state ($n_H = 10^6$ cm⁻³ or $10^7$ cm⁻³ depending on the model). During the cloud collapse stage, visual extinction ($A_V$) also varies from its minimum to maximum value. The first stage is followed by a warm-up stage during which the gas temperature ($T_g$) and the
<figure><img src="image_8.png" /><figcaption>Fig. 7. Gaussian-fit spectra and LTE spectra (same as Fig. 2) of cyanoacetylene, cyanodiacetylene, cyanamide, and aminoacetonitrile.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 8. Black line represents the expected line intensity ratio variation with the kinetic temperature, whereas the red line represents the observed intensity ratio. The ratios of CH<sub>3</sub>CCH transitions are plotted in panels a, b, and c, and that of CH<sub>3</sub><sup>13</sup>CCH is given in panel d.</figcaption></figure>

dust temperature ($T_d$) vary. The third stage is called the post-warm-up stage. We use the above-described three-stage physical model in our 0D and 1D simulations. First, we use 0D simulation for its simplicity to test if our adopted physical model can explain the observed abundances of N-bearing species in G10. 0D simulations do not give us any information about the possible spatial distribution of species abundances. So, we try to find the possible spatial distribution of simulated species in our 1D simulations. We explain our 0D and 1D simulations in Sects. 4.1 and 4.2, respectively. Results for both types of simulations are presented in Sects. 5.1 and 5.2, respectively.

### 4.1. 0D physical model

We started with the simple 0D simulation. We simulate the source as a single point and do not consider any spatial distribution of physical parameters ($n_H$, $A_V$, $T_g$, $T_d$). Following Gorai et al. (2020), we carried out our first simulation with generally accepted parameters for a source such as G10 (initially $n_H = 10^3$ cm$^{-3}$, which goes up to $10^7$ cm$^{-3}$ during cloud collapse; initial $A_V = 2$, which goes up to 250; initial $T_g = 40$ K, which goes up to 200 during the warm-up stage; and initial $T_d = 20$ K, which goes up to 200 K). Since G10 is massive hot
<table><caption>**Table 2.** Emitting regions of observed molecules.</caption><thead><tr><th>Molecules</th><th>Frequency (GHz)</th><th>E<sub>u</sub> (K)</th><th>Diameter (arcsec)</th></tr></thead><tbody><tr><td rowspan="7">C<sub>2</sub>H<sub>3</sub>CN</td><td>130.763576</td><td>18.2</td><td>1.04 ± 0.033</td></tr><tr><td>131.168737</td><td>209.5</td><td>0.77 ± 0.030</td></tr><tr><td>131.267473<sup>(*)</sup></td><td>47.5</td><td>1.07 ± 0.022</td></tr><tr><td>147.561701<sup>(*)</sup></td><td>62.5</td><td>1.28 ± 0.059</td></tr><tr><td>153.42175</td><td>44.13</td><td>1.11 ± 0.015</td></tr><tr><td>153.518944<sup>(*)</sup></td><td>70.9</td><td>1.48 ± 0.016</td></tr><tr><td>154.72454<sup>(*)</sup></td><td>65.4</td><td>1.10 ± 0.016</td></tr><tr><td rowspan="3"><sup>13</sup>CH<sub>2</sub>CHCN</td><td>147.927607</td><td>164.8612</td><td>1.20 ± 0.054</td></tr><tr><td>147.986686</td><td>79.5</td><td>1.17 ± 0.045</td></tr><tr><td>149.423715</td><td>69.2</td><td>1.16 ± 0.049</td></tr><tr><td>159.954635</td><td>71.3</td><td>0.73 ± 0.048</td></tr><tr><td>CH<sub>2</sub><sup>13</sup>CHCN</td><td>153.925606</td><td>65.0</td><td>0.78 ± 0.02</td></tr><tr><td>CH<sub>2</sub>CH<sup>13</sup>C-N</td><td>154.034285</td><td>65.1</td><td>0.70 ± 0.029</td></tr><tr><td rowspan="2">CH<sub>2</sub>CDCN</td><td>130.720066</td><td>74.1</td><td>0.69 ± 0.059</td></tr><tr><td>159.38385</td><td>83.9</td><td>0.64 ± 0.071</td></tr><tr><td rowspan="12">C<sub>2</sub>H<sub>5</sub>CN</td><td>129.76814</td><td>44.9</td><td>1.23 ± 0.023</td></tr><tr><td>129.795661<sup>(*)</sup></td><td>51.1</td><td>1.47 ± 0.021</td></tr><tr><td>130.693882</td><td>109.4</td><td>1.13 ± 0.025</td></tr><tr><td>130.903902<sup>(*)</sup></td><td>50.8</td><td>1.63 ± 0.052</td></tr><tr><td>148.293988</td><td>147.03</td><td>1.32 ± 0.024</td></tr><tr><td>147.756711<sup>(*)</sup></td><td>64.5</td><td>1.63 ± 0.052</td></tr><tr><td>148.36276</td><td>28.1</td><td>1.48 ± 0.019</td></tr><tr><td>153.272198<sup>(*)</sup></td><td>75.9</td><td>1.59 ± 0.018</td></tr><tr><td>153.589524</td><td>184.6</td><td>1.10 ± 0.017</td></tr><tr><td>154.244297<sup>(*)</sup></td><td>68.2</td><td>1.67 ± 0.016</td></tr><tr><td>154.47598</td><td>117.2</td><td>1.22 ± 0.017</td></tr><tr><td>154.55736</td><td>154.0</td><td>1.14 ± 0.015</td></tr><tr><td>158.742366</td><td>161.4</td><td>1.08 ± 0.049</td></tr><tr><td>159.3308</td><td>83.6</td><td>1.22 ± 0.026</td></tr><tr><td>159.888831<sup>(*)</sup></td><td>77.6</td><td>1.46 ± 0.027</td></tr><tr><td rowspan="5"><sup>13</sup>C-H<sub>3</sub>CH<sub>2</sub>CN</td><td>131.155447</td><td>77.8</td><td>1.21 ± 0.023</td></tr><tr><td>131.161428</td><td>183.4</td><td>0.93 ± 0.040</td></tr><tr><td>131.233679</td><td>67.9</td><td>1.17 ± 0.025</td></tr><tr><td>131.249984</td><td>67.9</td><td>1.29 ± 0.026</td></tr><tr><td>159.10002</td><td>74.2</td><td>1.18 ± 0.031</td></tr><tr><td>CH<sub>3</sub>CH<sub>2</sub><sup>13</sup>C-N</td><td>154.679342</td><td>71.8</td><td>1.22 ± 0.022</td></tr><tr><td rowspan="2">C<sub>2</sub>H<sub>5</sub>NC</td><td>149.120363</td><td>62.1</td><td>1.34 ± 0.040</td></tr><tr><td>159.061678</td><td>69.60</td><td>0.75 ± 0.060</td></tr><tr><td>HC<sub>3</sub>N</td><td>154.657284</td><td>66.8</td><td>2.18 ± 0.031</td></tr><tr><td>HC<sup>13</sup>CCN</td><td>154.001217</td><td>66.5</td><td>1.17 ± 0.019</td></tr><tr><td>HCC<sup>13</sup>CN</td><td>154.016078</td><td>66.5</td><td>1.16 ± 0.017</td></tr><tr><td>HC<sub>5</sub>N</td><td>130456.437</td><td>156.5</td><td>1.20 ± 0.027</td></tr><tr><td>HCC<sup>13</sup>C-CCN</td><td>130321.351</td><td>156.3</td><td>0.95 ± 0.030</td></tr><tr><td rowspan="2">NH2CN</td><td>158.891146</td><td>48.8</td><td>0.80 ± 0.109</td></tr><tr><td>160.104805</td><td>164.9</td><td>0.91 ± 0.026</td></tr><tr><td>NH<sub>2</sub>CH<sub>2</sub>CN</td><td>154.51748</td><td>86.4</td><td>1.07 ± 0.018</td></tr><tr><td rowspan="3">CH<sub>3</sub>CCH</td><td>153.71152</td><td>296.1</td><td>1.28 ± 0.021</td></tr><tr><td>153.790769</td><td>101.7</td><td>1.92 ± 0.037</td></tr><tr><td>153.805457</td><td>65.71</td><td>1.80 ± 0.06</td></tr><tr><td rowspan="2">CH<sub>3</sub><sup>13</sup>CCH</td><td>153.72701</td><td>101.7</td><td>1.01 ± 0.10</td></tr><tr><td>153.75337</td><td>36.9</td><td>1.31 ± 0.034</td></tr></tbody></table>

**Notes.** Spatial resolution at four different frequencies are 1.5"(129.5–131.5 GHz), 1.72"(147.5–149.4 GHz), 1.66"(153–154.9 GHz), and 1.76"(158.5–160.43 GHz), (*) optically thick line.

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Time evolution of $n_H$, $A_v$, and $T_d = T_g$ considered from van der Tak et al. (2013) for ten different grids. Red and blue vertical dashed lines show the end of collapse and warm-up phases, respectively.</figcaption></figure>

core, we selected a collapse time of $10^5$ yr and a warm-up time of $5 \times 10^4$ yr.

Knowing that chemistry is very sensitive to certain physical parameters such as temperature and density, we are required to test a number of 0D models by varying the values of most relevant physical parameters. Hence, first we repeat our simulation with a different peak core density of $10^6$ cm$^{-3}$. We also simulated a number of different models where we used different values of peak gas temperature (400, 300, and 200 K) during the warm-up phase and initial dust temperature (15, 20, and 25 K) during the collapse phase to cover the possible physical evolution history of hot cores such as G10. The warm-up phase is a very important phase in chemical evolution history as the dust grain chemistry is very sensitive to dust temperature. This makes the warm-up duration a very critical parameter. Therefore, we used three different values for the total duration ($T_w = 10^4$, $5 \times 10^4$, and $10^5$ yr) of the warm-up phase in order to properly investigate the dependence of simulated abundances of N-bearing species on the warm-up time of the collapsed cloud.

Depending on the values of physical parameters used during the time evolution of these three-stages, a number of 0D models are implemented to study the chemical evolution in G10-like dense cores. In Table 3, we list values of different physical parameters used in different models. More details about these models can be found in Sect. 5.1, where we discuss results for the 0D model.

## 4.2. 1D physical model

To study the spatial distribution of the chemical abundances of G10, we considered the density and temperature distribution obtained by SED fitting by van der Tak et al. (2013). Very recently, Bhat et al. (2021) constructed a similar model for
<table><caption>**Table 3.** Simulated models and physical parameters that describe the 0D model.</caption><thead><tr><th>Simulated models</th><th>t<sub>c</sub> (yr)</th><th>t<sub>w</sub> (yr)</th><th>t<sub>pw</sub> (yr)</th><th>n<sub>H</sub> (cm<sup>-3</sup>)</th><th>A<sub>V</sub></th><th>T<sub>g</sub> (K)</th><th>T<sub>d</sub> (K)</th></tr></thead><tbody><tr><td>N7-Tg200-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–200</td><td>20–200</td></tr><tr><td>N7-Tg300-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–300</td><td>20–300</td></tr><tr><td>N7-Tg400-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–400</td><td>20–400</td></tr><tr><td rowspan="4">Set 1</td><td>N7-Tg200-Td15</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–400</td><td>15–200</td></tr><tr><td>N7-Tg200-Td25</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–400</td><td>25–200</td></tr><tr><td>N7-Tg200-Td20-slow-wup</td><td>10<sup>5</sup></td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–400</td><td>20–200</td></tr><tr><td>N7-Tg200-Td20-fast-wup</td><td>10<sup>5</sup></td><td>10<sup>4</sup></td><td>1.4 × 10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>7</sup></td><td>2–250</td><td>40–400</td><td>20–200</td></tr><tr><td rowspan="6">Set 2</td><td>N6-Tg200-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–200</td><td>20–200</td></tr><tr><td>N6-Tg300-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–300</td><td>20–300</td></tr><tr><td>N6-Tg400-Td20</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–400</td><td>20–400</td></tr><tr><td>N6-Tg200-Td15</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–400</td><td>15–200</td></tr><tr><td>N6-Tg200-Td25</td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>5</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–400</td><td>25–200</td></tr><tr><td>N6-Tg200-Td20-slow-wup</td><td>10<sup>5</sup></td><td>10<sup>5</sup></td><td>5 × 10<sup>4</sup></td><td>10<sup>3</sup>–10<sup>6</sup></td><td>2–250</td><td>40–400</td><td>20–200</td></tr></tbody></table>

**Notes.** During cloud collapse period $t_c$: $T_g$ and $T_d$ values remain constant at their initial values and $n_H$ and $A_V$ vary from initial to final values as given above. During the warm-up period $t_w$: $T_d = T_g$, which varies from the initial to final value as given above, while $n_H$ and $A_V$ values remain unchanged at their final values. During post-warm-up period $t_{pw}$: all $n_H$, $A_V$, $T_g$, and $T_d$ remain constant at their final values. Total time of simulation in each model is $t_c + t_w + t_{pw} = 2.5 \times 10^5$ yr.

G31.041+0.31. They divided the cloud envelope into 23 logarithmic spacing shells. Here, for G10, we divided the cloud envelope into ten logarithmic spacing shells. The time evolutions of the physical parameters ($n_H$, $A_V$, $T_d = T_g$) for these ten shells are shown in Fig. 9.

For each shell, we consider that, at $t = 0$ yr, all the shells have $n_H = 10^3$ cm$^{-3}$. In the first stage, every shell is allowed to evolve up to its spatial density distribution obtained by van der Tak et al. (2013). For example, for the fifth shell's (at ~1717 AU) density may evolve up to $2.08 \times 10^7$ cm$^{-3}$ at the end of $10^5$ yr. Beyond that (warm-up and post-warm-up stages), the density of each shell is kept constant at its maximum value. Since dust temperature is a critical parameter, we used two profiles (15 and 20 K) for the initial dust temperature. For simplicity, we kept the gas temperature $T_g$ equal to the dust temperature $T_d$. In the first stage, $T_d = T_g$ was kept constant at initial value (15 or 20 K depending on model). In the warm-up stage, $T_d = T_g$ evolves from its initial value to final value respective of the spatial temperature distribution as reported by van der Tak et al. (2013; for example, 207 K for the fifth shell). In the post-warm-up stage, gas and dust temperatures are kept constant at their highest values respective to each shell. To calculate $A_V$, we used the relation $A_V = A_{V,0} \times (n_H/n_{H,0})^{(2/3)}$, as given in Sect. 3.3 of Garrod & Pauly (2011), where $A_{V,0}$ is the initial value of visual extinction and $n_{H,0}$ is the initial total number density. We used $A_{V,0} = 2$ and $n_{H,0} = 10^3$ cm$^{-3}$.

### 4.3. Chemical model

To explain the observed abundances of chemical species in G10, a three-phase (gas phase, grain surface, and grain mantle) version of the Nautilus gas-grain chemical code (Ruaud et al. 2016) is used. This code is based on the rate equation approximation (Hasegawa et al. 1992; Hasegawa & Herbst 1993). In all our simulations, the surface of the dust grain is defined as the top two monolayers of species. All ice species under the surface are treated as mantles. When a species desorbs from the surface, another species from the mantle comes to the top to form the

new surface layer, and, similarly, with the accretion of a new species on the surface, another surface species moves to the mantle to form a new mantle layer. In our model, standard physical processes such as accretion, diffusion, recombination, and thermal desorption of species are considered. The cosmic-ray-induced desorption, UV (direct and indirect) photo-desorption, and chemical desorption are also considered. These physical processes are explained in detail in Ruaud et al. (2016); Iqbal et al. (2018); Das et al. (2019); Wakelam et al. (2020). We used a common cosmic-ray ionization rate of $1.3 \times 10^{-17}$ s$^{-1}$. In Table 4, we list the initial abundances of all the gas and the grain species we used in our all chemical models.

### 4.4. Chemical network

Our gas-phase chemistry is mostly based on the kida.uva.2014 public network (Wakelam et al. 2015). We updated our gas phase chemistry by adding all relevant reactions from KIDA and UMIST databases. Apart from this we also followed Tables 8 and 9 of Gorai et al. (2020), and Table 5 of Sil et al. (2018) to complete our gas-phase nitrogen chemistry. Our surface chemistry is mostly based on Iqbal et al. (2018). To complete the surface chemistry, we added all relevant formation and destruction reactions in our chemical network from Table 14 of Belloche et al. (2009) and Table 8 of Gorai et al. (2020). We also added all surface reactions, related to the formation and destruction of $C_2H_3CN$ and $C_2H_5CN$, given in Tables B.3 and B.4 of Garrod et al. (2017).

## 5. Chemical simulation results

We present our results in two sections. First, we present our simulation results for 0D models shown in Table 3. Results from these simulations were considered to select the most suitable parameters for our 1D simulations. Results for 1D simulations are presented in Sect. 5.2. To find the best-fit-simulated model, we used a simple method in which we calculate the mean
<table><caption>**Table 4.** Elemental abundances and initial abundances.</caption><thead><tr><th>Element</th><th>Abundance relative to H</th><th>References</th></tr></thead><tbody><tr><td>H<sub>2</sub></td><td>0.5</td><td></td></tr><tr><td>He</td><td>0.09</td><td>(a)</td></tr><tr><td>N</td><td>6.2 × 10<sup>-5</sup></td><td>(b)</td></tr><tr><td>O</td><td>2.4 × 10<sup>-4</sup></td><td>(c)</td></tr><tr><td>C<sup>+</sup></td><td>1.7 × 10<sup>-4</sup></td><td>(b)</td></tr><tr><td>S<sup>+</sup></td><td>8.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>Si<sup>+</sup></td><td>8.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>Fe<sup>+</sup></td><td>3.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>Na<sup>+</sup></td><td>2.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>Mg<sup>+</sup></td><td>7.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>P<sup>+</sup></td><td>2.0 × 10<sup>-10</sup></td><td>(d)</td></tr><tr><td>Cl<sup>+</sup></td><td>1.0 × 10<sup>-9</sup></td><td>(d)</td></tr><tr><td>Ice</td><td>0</td><td></td></tr><tr><td>All other gas species</td><td>0</td><td></td></tr></tbody></table>

**References.** (a)Wakelam & Herbst (2008), (b)Jenkins (2009), (c)Hincelin et al. (2011), (d)Graedel et al. (1982).

distance of disagreement between the simulated and observed values. This method is explained in detail in Iqbal et al. (2018).

## 5.1. Results obtained with OD model

Density is one of the fundamental properties of a core that determines the future state of the core. It is challenging to consider the final density of a collapsed core just after its isothermal collapse without having comprehensive knowledge of its collapse. Hence, in our simulations, we adopt two final density profiles for G10 to study their impact on chemical abundances at the end of the simulation.

### 5.1.1. Impact of peak H density

In Fig. 10, we show simulated abundances (gas and total ice) as a function of time for nine N-bearing species and CH<sub>3</sub>CCH (observed in G10) during the collapse stage, where total ice abundance includes the contribution of both surface and mantle. In all plots, solid lines are for the model N7-Tg200-Td20 (final $n_H = 10^7$ cm<sup>-3</sup>) and dotted lines are for the model N6-Tg200-Td20 (final $n_H = 10^6$ cm<sup>-3</sup>). In Fig. 10, we notice a general trend that is followed by all species. In denser clouds (plotted with solid lines), as the cloud collapses and the gas density increases, it takes less than $2 \times 10^4$ yr for all species to attain their maximum abundance values in the gas phase. This is followed by a rapid depletion of these species before a relatively stable state is reached; this state is maintained until the start of warm-up phase. At this stage, the abundances of different species continue to evolve, but at a very low rate.

During the isothermal collapse, the dust temperature is only 20K, which is too low for thermal desorption of these complex organic species from the dust surface. We further checked to verify that nonthermal desorption does not contribute to this rapid growth in gas phase species. In fact, a rapid increase in the gas phase abundances is almost entirely driven by the gas-phase formation of these species. During this phase, HC<sub>5</sub>N is mostly produced in the gas phase through H<sub>2</sub>C<sub>5</sub>N<sup>+</sup> + e<sup>-</sup> → H + HC<sub>5</sub>N and N + C<sub>6</sub>H → C + HC<sub>5</sub>N. Similarly, two major reactions

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Impact of peak H density (solid lines ($10^7$ cm<sup>-3</sup>) versus dotted lines ($10^6$ cm<sup>-3</sup>)) on simulated abundances of selected species during the cloud-collapse phase. Plots on the left and the right are for the ice and gas-phase abundances, respectively. Legends in the top left plot are for both plots at the top, and legends in the bottom left plot are for both plots at the bottom. All solid lines are for model N7-Tg200-Td20, and all dotted lines are for model N6-Tg200-Td20 (see Table 3 for details of each model).</figcaption></figure>

contributing to the rapid growth of HC<sub>3</sub>N in the gas phase are N + C<sub>4</sub>H → C + HC<sub>3</sub>N and C + H<sub>2</sub>CCN → H + HC<sub>3</sub>N.

As expected, in model N6-Tg200-Td20, with a lower peak H density (plotted with dotted lines), gas phase formation and destruction rates are relatively small compared to model N7-Tg200-Td20, and it takes longer (up to $6 \times 10^4$ yr) to reach a similar state of chemical evolution to that seen in model N7-Tg200-Td20.

During the warm-up stage (see Fig. 11), the gas-phase abundances rise rapidly due to desorption from the dust surface as surface temperature rises. HNCO is the first to desorb (just above 60 K) from the surface. All other species start rapid desorption when the surface temperature rises above 100 K. The best fit time for model with lower peak H density also lies during the early stages of warm-up phase when abundances of most species lies in the range of $10^{-10}$ to $10^{-9}$.

We note that the model with lower peak H density ($10^6$ cm<sup>-3</sup>) results in higher abundances of the considered species at the end of the cloud-collapse stage. However, this difference rapidly disappears during the warm-up stage (except for C<sub>2</sub>H<sub>3</sub>CN and C<sub>2</sub>H<sub>5</sub>CN). We also note that the gas-phase abundance is affected by considering a lower value of peak H density, while the ice abundances remain very much unaffected.

### 5.1.2. Impact of initial dust temperature

Surface chemistry is very sensitive to small changes in surface temperature. Hence, we used our models (N7-Tg200-Td20 and N6-Tg200-Td20) presented in Sect. 5.1.1 and ran them again with initial dust temperatures of 15 and 25 K. Thus, in total we ran four additional models: N7-Tg200-Td15 and N7-Tg200-Td25
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Same as Fig. 10, but for warm-up and post-warm-up stages. Plots on the left show the ice abundance (warm-up phase only). Plots in the middle and right show the gas-phase abundances during the warm-up and the post-warm-up phases, respectively.</figcaption></figure>

with a peak H density of $10^7$ cm$^{-3}$, and N6-Tg200-Td15 and N6-Tg200-Td25 with a peak H density of $10^6$ cm$^{-3}$.

In Fig. 12, we show the chemical evolution of selected species as a function of time in the gas phase during all three stages of physical structure evolution of a simulated cloud. In this figure, the pink horizontal bar represents the observed abundance. We note that, in general, the species have the highest gas-phase abundance during the cloud-collapse stage and the warm-up stage in the model with the highest dust temperature (25 K), and the lowest abundance when the dust temperature is lowest (15 K), with a few exceptions such as CH$_3$NCO and CH$_3$CCH. In fact, CH$_3$NCO and CH$_3$CCH are produced efficiently to cross the limit of observed abundance in the model with initial dust temperature of only 15 K. In the other two models, abundances of these species are always lower than observed abundances by a small factor. We see that for most of the species, the best-fit time falls during the early to mid-warm-up phase of the cloud.

We can conclude that initial dust temperature has a major impact on the time evolution of the selected species, with a higher dust temperature resulting in higher abundance values in general. A model with an initial dust temperature of 15 K produces a difference of few orders of magnitude in the abundances of a number of species when compared to a model with an initial dust temperature of 20 K. This difference is higher during the cloud-collapse stage and lowest during the post-warm-up stage. When we increase the dust temperature to 25 K, the difference remains small compared to the 20 K model. Hence, it is critical to select the initial dust temperature during the cloud-collapse stage, and it may result in a difference of few orders of magnitude if a temperature under 20 K is used.

### 5.1.3. Sensitivity to initial gas temperature and the peak gas and dust temperature

In all models discussed so far, we used an initial gas temperature of 40 K and the peak gas and dust temperature of 200 K. We

found that chemistry is not very sensitive to the initial gas temperature if the temperature is below 100 K. So, we do not show results for different initial gas temperatures here.

However, we found that variations in peak gas temperature do have noticeable impact on chemistry, especially during the warm-up phase. We know that in a collapsed cloud, gas and dust temperatures are coupled due to extreme density, such that $T_d = T_g$. Basically, the peak gas and dust temperatures are the same, and both change with the models presented in this section. Again, we used our reference models (N7-Tg200-Td20 and N6-Tg200-Td20) presented in Sect. 5.1.1 and ran four additional models to show the impact of variations in the peak gas and dust temperatures on chemical abundances of N-bearing species. We ran new models with the peak gas and dust temperatures of 300 K and 400 K; these models are N7-Tg300-Td20 and N7-Tg400-Td20, respectively, with a peak H density of $10^7$ cm$^{-3}$, and N6-Tg300-Td20 and N6-Tg400-Td20, respectively, with a peak H density of $10^6$ cm$^{-3}$.

In Fig. 13, we show simulated abundances as a function of time for selected species in the gas phase during all the three stages of physical structure evolution of simulated cloud. We see that the variation in peak gas and dust temperature only shifts the steady state abundance time for different species during warm-up stage. This shift in time is on the order of $10^4$ yr. We also note that, during the warm-up phase, the steady-state condition is reached earlier, and this happens at a higher gas and dust temperature if we increase the peak value of the gas and dust temperatures in different models. Once the gas and dust temperature goes above 200 K, there is little difference observed in abundances between different models. Our best-fit time falls again during the warm-up stage and is only shifted to the left by a factor of $10^4$ yr depending on the peak gas temperature.

### 5.1.4. Impact of warm-up duration

We implemented three different values for the total duration of the warm-up stage: $t_w = 10^4$, $5 \times 10^4$, and $10^5$ yr, respectively,
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Impact of initial dust temperature (black (15 K), red (20 K), and blue (25 K) lines) and final H density (dashed lines ($10^7$ cm$^{-3}$), and dotted lines ($10^6$ cm$^{-3}$)) on simulated abundances of selected species. For all plots, black, red, and blue dashed lines are for models N7-Tg200-Td15, N7-Tg200-Td20, and N7-T2400-Td25, respectively, and black, red, and blue dotted lines are for models N6-Tg200-Td15, N6-Tg200-Td20, and N6-Tg200-Td25, respectively. Pink horizontal bar represents the observed value in G10 (see Table 3 for details of each model).</figcaption></figure>

for fast, mid, and slow cases. In Fig. 14, we show the impact of different warm-up durations on the time-dependent chemical evolution of selected species. In Fig. 14, red, blue, and green vertical dotted lines mark the end of the warm-up period and the start of the post-warm-up stage for fast, mid, and slow models,

respectively. During the warm-up stage, the initial dust temperature is 20 K and the gas temperature is 40 K. During warm-up, the dust temperature is coupled with the gas temperature and climbs to its peak value (200 K). Due to a gradual increase in dust temperature, thermal desorption becomes very efficient
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Impact of peak gas and dust temperature (black, red, and blue lines) and peak H density (dashed lines ($10^7$ cm$^{-3}$) versus dotted lines ($10^6$ cm$^{-3}$)) on simulated abundances of selected species. For all plots, black, red, and blue dashed lines are for models N7-Tg200-Td20, N7-Tg300-Td20, and N7-Tg400-Td20, respectively, and black, red, and blue dotted lines are for models N6-Tg200-Td20, N6-Tg300-Td20, and N6-Tg400-Td20. Pink horizontal bar represents the observed value in G10 (see Table 3 for details of each model).</figcaption></figure>

even for heavier species with a binding energy above 4000 K. Due to efficient thermal desorption, almost all the ice species are evaporated within $10^4$ yr once the dust temperature goes above 100 K and continues to rise. Soon after that, abundances of gas species reach a steady-state value. In Fig. 14, we can see that

abundances of all species rise fastest in a fast warm-up model, as the rise in temperature is steeper, followed by mid and slow warm-up models. However, there is no significant difference in the steady state abundance values in different models after $2 \times 10^5$ yr, except for $C_2H_3CN$ and $C_2H_5CN$, as our simulation
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Impact of warm-up time (fast, mid, and slow cases are plotted in red, blue, and green, respectively) and final H density (dashed line (10⁷ cm⁻³) versus dotted line (10⁶ cm⁻³)) on simulated abundances of selected species. Red, blue, and green dashed lines are for models N7-Tg200-Td20-fast-wup, N7-Tg200-Td20, and N7-Tg200-Td20-slow-wup, respectively, and red, blue, and green dotted lines are for models N6-Tg200-Td20-fast-wup, N6-Tg200-Td20, and N6-Tg200-Td20-slow-wup, respectively. All plots are for both warm-up and post-warm-up stages, with both stages separated by red, blue, and green vertical dotted lines for fast, mid, and slow cases, respectively. Pink horizontal bar represents the observed value in G10 (see Table 3 for details of each model).</figcaption></figure>

shows that these two species are sensitive to the warm-up duration. In all three models, simulated abundances of all selected species fit well with observed values during the warm-up stage.

CH₃CCH is slightly under-produced in all our models. CH₃CCH is efficiently produced only when the initial dust temperature is 15 K (see Fig. 12).
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Time evolution of abundances of selected species during the warm-up and post-warm-up phases for 1D physical model with initial dust and gas temperature of 15 K. Red dashed vertical line at $1.5 \times 10^5$ yr shows the end of the warm-up phase. Pink horizontal bar represents observed abundance range.</figcaption></figure>

## 5.2. Results obtained with 1D model

Here, we show the results for 1D simulations, which are essential to studying the possible spatial distribution of N-bearing species in the source G10. We ran two simulations (based on initial dust temperatures of 15 and 20 K) of the 1D model. In each simulation, we have ten grids representing ten spatial points in the source. Details on these spatial points are given in Sect. 4.2.

The results showing the time evolution of the abundances of selected species for all ten grids are shown in Fig. 15. In this simulation, the initial gas and dust temperature was set at 15 K. The same plot for an initial dust and gas temperature of 20 K is given in Fig. E.1. The observed abundances are marked with horizontal bars. Grid 1 is the innermost grid, whereas Grid 10 is the outermost. Our model shows that the abundances of all species vary by a few orders of magnitude between the inner and outer grids. Some species such as HC₃N and HC₅N are most abundant at inner grids, while others such as NH₂CN and CH₃CCH are most abundant in the outer grids.

For a better comparison with values presented in Rolffs et al. (2011b), we show the spatial distribution of peak abundances of selected species in Fig. 16. Again, this plot is for models with initial dust temperatures of 15 K. A similar plot with an initial dust

temperature of 20 K is given in Fig. E.2. The dotted and dashed green vertical lines represent the observed spatial resolutions of Rolffs et al. (2011b), whereas the dashed pink vertical line represents the spatial resolution of our observation. Peak abundances (black curves) are extracted by noting the abundances achieved during warm-up and post-warm-up periods. To avoid the time uncertainty of the peak values, we also show the abundances obtained just after the warm-up phase ($1.5 \times 10^5$ yr) and final time ($2.5 \times 10^5$ yr).

As shown in Fig. 16, the spatial resolution of the observation presented in Rolffs et al. (2011b) is much higher than that of the observations in this paper. We also note that the simulated peak abundances of all species is at its maximum toward the inner grids, though not in the same place. This is important to understand why, for most of the species, Rolffs et al. (2011b) recorded a higher maximum abundance compared to us for the source G10. Our 1D simulations shows a high abundance of HCN and HC₃N deeper inside the cloud (i.e., around the high-temperature region), which could not be resolved with observation presented in this paper. However, HC₃N is more extended compared to HC₅N. This is reflected in our modeling results (Figs. 16 and E.2). They show that HC₅N decays in the outer part of the envelope more rapidly compared to HC₃N. The abundance of vinyl cyanide and ethyl cyanide seems to be
<figure><img src="image_17.png" /><figcaption>**Fig. 16.** Spatial distribution of abundances of selected species at three selected times in simulation. Red, blue, and black curves show the spatial distribution of abundances at ~1.5 × 10⁵ yr, ~2.5 × 10⁵ yr, and the peak abundance, respectively. Peak abundance is the max abundance during the warm-up and post-warm-up phases. Here, we considered the initial dust temperature of 15 K. Green dotted (at 2580 AU), green dashed (at 12 900 AU), and pink dashed (at 17 200 AU) lines represent the highest spatial resolution (0.3") and the moderate resolution (1.5") of Rolffs et al. (2011b), and our resolution (2"), respectively. Pink horizontal bar represents observed abundance range.</figcaption></figure>

heavily dependent on the initial ice temperature. In the model with the initial gas and dust temperature of 20 K, shown in Fig. E.1, abundances of vinyl cyanide and ethyl cyanide drop rapidly compared to the model with the initial gas and dust temperature of 15 K, shown in Fig. 15. More interestingly, inside 3000 AU, it depicts a higher abundance of ethyl cyanide than vinyl cyanide. Rolffs et al. (2011b) also obtained a higher abundance of ethyl cyanide (9 × 10¹⁷ cm⁻²) compared to vinyl cyanide (7 × 10¹⁷ cm⁻²). In contrast, we obtained a higher abundance of vinyl cyanide than ethyl cyanide with our resolution. This is also reflected in our model. So, our model at the initial ice phase temperature of 15 K successfully explains a high abundance of ethyl cyanide over vinyl cyanide deep inside, and of vinyl cyanide over ethyl cyanide at the outer part of the envelope. An abundance of CH₃NH₂ seems roughly invariant across the cloud. The abundances of CH₂NH and CH₃NH₂ are overproduced in our 15 K model, whereas they are underproduced in the 20 K model. Results presented here can explain the observed abundances of the peptide-like bond containing species (HNCO, CH₃NCO, and NH₂CHO). With the initial gas and dust temperature of 15 K, CH₃CCH shows a good match with the observational results (see Fig. 16).

Simulated (obtained in best-fit 1D model) and observed column densities of the identified species are shown in Fig. 17. Various methods are implemented to derive the observed column density of the species. RD (shown in green) represents the column density obtained with the rotational diagram method, and MCMC (shown in blue) illustrates the MCMC Method. LTE1 (shown in red) represents the cases when a fixed source size of ~2" is used, whereas LTE2 represents the circumstances for which the average emitting region of species are used (shown in black). Figure 17 shows that all the methods of deriving the observed column density of a species agree well among themselves. The simulated column densities of species are shown in yellow. The grid number (noted in black) represents the peak column density of that grid, for which we have a good match with the observed results. Additionally, for the cases of HC₃N and HC₅N, we include the results obtained for the same grid at ~2.5 × 10⁵ yr.

Figure 17 summarizes our findings. For example, we obtain a good match between the observed and simulated results for NH₂CH₂CN in the innermost grid and for NH₂CN in the fourth grid. It signifies the presence of these species in the inner part of the hot core. Figure 16 shows that in the case of C₂H₅CN, we
<figure><img src="image_18.png" /></figure>

**Fig. 17.** Observed and simulated column density of the species are shown. LTE1 represents column density obtained with a constant source size of 2", whereas LTE2 represents the cases with the average emitting regions of each species. Simulated column densities (obtained in best-fit 1D model) are shown in yellow. The grid number (noted in black) of the modeled column density represents the peak column density of that grid, which matches well with the observed results. For HC₃N and HC₅N, the modeled column density at ~2.5 × 10⁵ yr is additionally shown for the same grid number noted in red. Error bars obtained with the rotational diagram and MCMC methods are shown in magenta.

obtain a dramatic change in column density between grid numbers 6 and 7. Thus, the column density obtained at the middle of 6 and 7 (i.e., 6.5) is considered from the column density profile of C₂H₅CN. Overall, we note that the C₂H₃CN emission and CH₃CCH emission are more extended (these match well with the ninth grid) than that of C₂H₅CN (~6.5th grid). For HC₃N and HC₅N, we obtain a good match with the observation during the latter stages of evolution (~2.5 × 10⁵) in the ninth grid. However, a higher column density (simulated and observed) of HC₃N compared to HC₅N implies a relatively extended emission of HC₃N.

## 6. Conclusions

This paper reports the identification of some nitrogen-bearing species in the high-mass protostellar envelope of G10 using the archival ALMA data. We also used the Nautilus gas-grain code to investigate the evolution of Nitrogen chemistry in a hot molecular source such as G10. For our chemical modeling, we utilized a physical structure similar to that observed in G10. Our chemical network contains both the gas- and dust-phase chemistry. The main highlights of this paper are presented below.

– We identify numerous transitions of nitrogen-bearing species in G10. Among them, vinyl cyanide and its one ¹³C isotopologue (¹³CH₂CHCN) and ethyl cyanide and its one ¹³C isotopologue (¹³CH₃CH₂CN) are identified. Two other ¹³C isotopologues of vinyl cyanide (CH₂¹³CHCN and CH₂CH¹³CN), one deuterated form of vinyl cyanide (CH₂CDCN), another ¹³C isotopologue of ethyl cyanide (CH₃CH₂¹³CN), an isomer of ethyl cyanide (C₂H₅NC), cyanoacetylene and its two ¹³C isotopologues (HC¹³CCN and HCC¹³CN), cyanodiacetylene and its one ¹³C isotopologue (HCC¹³CCCN), cyanamide, and aminoacetonitrile are tentatively identified. Additionally, one ¹³C isotopologue of (CH₃¹³CCH) is tentatively identified with our investigation.

– Furthermore, we also identify multiple transitions of methyl acetylene in this source, which is a valuable tracer of

physical condition. From our analysis, a kinetic temperature of 133–193 K is estimated.

– Based on our detection of vinyl cyanide and one of its ¹³C isotopologues, a ¹²C/¹³C ratio of 45–64 is obtained. Similarly, with the ethyl cyanide and its one ¹³C isotopologue, a ¹²C/¹³C ratio of 12–27 is obtained. Additionally, a ratio of 57.5 is obtained between ethyl cyanide and ethyl isocyanide.

– Our 0D model with an initial dust temperature of 15 K and a cloud density of 10⁷ cm⁻³ best describes the observed abundances in the source. Our best-fit model is consistent with the source properties of hot and dense sources similar to G10.

– 0D simulations show that the chemistry of certain N-bearing species is very sensitive to initial local conditions such as density or dust temperature or both, and it is not very sensitive to initial gas temperature.

– In our 1D model, simulated higher abundances of species such as HCN, HC₃N, and HC₅N toward the inner shells of the source confirm the observational findings.

– We obtained a lower abundance of vinyl cyanide, ethyl cyanide, cyanoacetylene, and cyanamide compared to Rolffs et al. (2011b). We note that the emission of most of these species is compact. Since Rolffs et al. (2011b) have a higher spatial resolution (highest resolution ~0.3" corresponds to 2580 AU and average resolution ~1.5" corresponds to 12 900 AU) than ours (~2" corresponds to 17 200 AU), they obtained a higher abundance of these species.

– Results obtained with our 1D model, shown in Figs. 16 and E.2, show that the abundance of most of the species increases as we move toward the denser and warmer inner regions of the source. In this work, we have identified various nitrogen-bearing species and investigated the N-bearing chemistry in G10. In addition, we have discussed how the chemistry of N-bearing molecules can vary with initial local conditions and investigated these species' spatial distribution. Since most of the transitions were marginally resolved, higher-angular-resolution observations would be needed to know their spatial distribution more accurately.

*Acknowledgements.* This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.00929.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. S.K.M. acknowledges CSIR fellowship (File no. 09/904(0014)/2018-EMR-I). W.I. acknowledges the fellowship of China Postdoctoral Science Foundation (Grant number: 2020M683386). P.G. acknowledges support from a Chalmers Cosmic Origins postdoctoral fellowship. B.B. gratefully acknowledges the DST, India's Government, for providing financial assistance through the DST-INSPIRE Fellowship [IF170046] scheme.

## References

Agúndez, M., Fonfría, J. P., Cernicharo, J., Pardo, J. R., & Guélin, M. 2008, A&A, 479, 493

Agúndez, M., Cernicharo, J., Guélin, M., et al. 2010, A&A, 517, A2

Agúndez, M., Cernicharo, J., de Vicente, P., et al. 2015, A&A, 579, A10

Balucani, N. 2009, Int. J. Mol. Sci., 10, 2304

Belloche, A., Garrod, R. T., Müller, H. S. P., et al. 2009, A&A, 499, 215

Belloche, A., Garrod, R. T., Müller, H. S. P., & Menten, K. M. 2014, Science, 345, 1584

Belloche, A., Garrod, R. T., Müller, H. S. P., et al. 2019, A&A, 628, A10

Bhat, B., Gorai, P., Mondal, S. K., Chakrabarti, S. K., & Das, A. 2021, Adv. Space Res., 69, 415

Cernicharo, J., Guélin, M., Agúndez, M., McCarthy, M. C., & Thaddeus, P. 2008, ApJ, 688, L83

Cernicharo, J., Tercero, B., Fuente, A., et al. 2013, ApJ, 771, L10

Cesaroni, R., Hofner, P., Araya, E., & Kurtz, S. 2010, A&A, 509, A50

Clayton, D. D., & Nittler, L. R. 2004, ARA&A, 42, 39
Danger, G., Borget, F., Chomat, M., et al. 2011, A&A, 535, A47

Das, A., Gorai, P., & Chakrabarti, S. K. 2019, A&A, 628, A73

Fontani, F., Pascucci, I., Caselli, P., et al. 2007, A&A, 470, 639

Gardner, F. F., & Winnewisser, G. 1975, ApJ, 195, L127

Garrod, R. T. 2013, ApJ, 765, 60

Garrod, R. T., & Pauly, T. 2011, ApJ, 735, 15

Garrod, R. T., Belloche, A., Müller, H. S. P., & Menten, K. M. 2017, A&A, 601, A48

Goldman, N., Reed, E. J., Fried, L. E., William Kuo, I. F., & Maiti, A. 2010, Nat. Chem., 2, 949

Goldsmith, P. F., & Langer, W. D. 1999, ApJ, 517, 209

Gorai, P., Bhat, B., Sil, M., et al. 2020, ApJ, 895, 86

Graedel, T. E., Langer, W. D., & Frerking, M. A. 1982, ApJS, 48, 321

Hasegawa, T. I., & Herbst, E. 1993, MNRAS, 263, 589

Hasegawa, T. I., Herbst, E., & Leung, C. M. 1992, ApJS, 82, 167

Hernández-Hernández, V., Zapata, L., Kurtz, S., & Garay, G. 2014, ApJ, 786, 38

Hincelin, U., Wakelam, V., Hersant, F., et al. 2011, A&A, 530, A61

Ibanez, J. D., Kimball, A. P., & Oro, J. 1971, Science, 173, 444

Iqbal, W., Wakelam, V., & Gratier, P. 2018, A&A, 620, A109

Irvine, W. M., & Schloerb, F. 1984, ApJ, 282, 516

Jenkins, E. B. 2009, ApJ, 700, 1299

Jiménez-Serra, I., Vasyunin, A. I., Caselli, P., et al. 2016, ApJ, 830, L6

Johnson, D. R., Lovas, F. J., Gottlieb, C. A., et al. 1977, ApJ, 218, 370

Jørgensen, J. K., Belloche, A., & Garrod, R. T. 2020, ARA&A, 58, 727

Kawaguchi, K., Ohishi, M., Ishikawa, S.-I., & Kaifu, N. 1992, ApJ, 386, L51

Kawaguchi, K., Kasai, Y., Ishikawa, S.-I., et al. 1994, ApJ, 420, L95

Linsky, J. L. 1998, Space Sci. Rev., 84, 285

López, A., Tercero, B., Kisiel, Z., et al. 2014, A&A, 572, A44

Lovas, F. J., Remijan, A. J., Hollis, J. M., Jewell, P. R., & Snyder, L. E. 2006, ApJ, 637, L37

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

Marcelino, N., Agúndez, M., Cernicharo, J., Roueff, E., & Tafalla, M. 2018, A&A, 612, L10

Margulès, L., Belloche, A., Müller, H. S. P., et al. 2016, A&A, 590, A93

McGuire, B. A. 2018, ApJS, 239, 17

McKellar, A. 1940, PASP, 52, 187

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser., 376, 127

Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, L49

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

Ohishi, M., Suzuki, T., Hirota, T., Saito, M., & Kaifu, N. 2019, PASJ, 71, 86

Olmi, L., Cesaroni, R., Neri, R., & Walmsley, C. M. 1996, A&A, 315, 565

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

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

Rolffs, R., Schilke, P., Zhang, Q., & Zapata, L. 2011b, A&A, 536, A33

Ruaud, M., Wakelam, V., & Hersant, F. 2016, MNRAS, 459, 3756

Sanna, A., Reid, M. J., Menten, K. M., et al. 2014, ApJ, 781, 108

Sil, M., Gorai, P., Das, A., et al. 2018, ApJ, 853, 139

Snyder, L. E., & Buhl, D. 1973, Nat. Phys. Sci., 243, 45

Suzuki, T., Ohishi, M., Hirota, T., et al. 2016, ApJ, 825, 79

Thaddeus, P., Gottlieb, C. A., Gupta, H., et al. 2008, ApJ, 677, 1132

Turner, B. E., Liszt, H. S., Kaifu, N., & Kisliakov, A. G. 1975, ApJ, 201, L149

van der Tak, F. F. S., Chavarría, L., Herpin, F., et al. 2013, A&A, 554, A83

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

Vastel, C., Bottinelli, S., Caux, E., Glorian, J. M., & Boiziot, M. 2015, in SF2A-2015: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics (SF2A), 313

Wakelam, V., & Herbst, E. 2008, ApJ, 680, 371

Wakelam, V., Loison, J.-C., Herbst, E., et al. 2015, ApJS, 217, 20

Wakelam, V., Iqbal, W., Melisse, J. P., et al. 2020, MNRAS, 497, 2309

Wyrowski, F., Schilke, P., & Walmsley, C. M. 1999, A&A, 341, 882

Ziurys, L. M., & Turner, B. E. 1986, ApJ, 302, L31

Ziurys, L. M., Savage, C., Brewster, M. A., et al. 1999, ApJ, 527, L67
## Appendix A: Observed line parameters

<table><caption>Table A.1. Summary of line parameters of observed molecules toward G10.47+0.03.</caption><thead><tr><th>Species</th><th>J<sub>K<sub>1</sub>K<sub>2</sub></sub><sup>+</sup>J<sub>K<sub>1</sub>K<sub>2</sub></sub><sup>-</sup></th><th>Frequency<br>(GHz)</th><th>E<sub>u</sub><br>(K)</th><th>g<sub>up</sub><br>(km s<sup>-1</sup>)</th><th>ΔV<br>(km s<sup>-1</sup>)</th><th>I<sub>max</sub><br>(K)</th><th>V<sub>LSR</sub><br>(km s<sup>-1</sup>)</th><th>Sμ<sup>2</sup><br>(Debye<sup>2</sup>)</th><th>∫ T<sub>mb</sub>dV<br>(K km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="11" style="text-align: center; font-style: italic;">Vinyl cyanide and its isotopologues</td></tr><tr><td rowspan="10">C<sub>2</sub>H<sub>3</sub>CN</td><td>6<sub>2,4</sub> → 6<sub>1,5</sub>,v=0</td><td>130.763576</td><td>18.2</td><td>39</td><td>9.5±0.5</td><td>3.9±0.2</td><td>68.1±0.2</td><td>8.0</td><td>39.2±3.9</td></tr><tr><td>29<sub>2,27</sub> → 29<sub>1,28</sub>,v=0</td><td>131.168737</td><td>209.5</td><td>177</td><td>7.6±0.2</td><td>7.5±0.2</td><td>67.2±0.1</td><td>54.3</td><td>61.2±3.3</td></tr><tr><td>14<sub>0,14</sub> → 13<sub>0,13</sub>,v=0</td><td>131.267473<sup>b</sup></td><td>47.5</td><td>87</td><td>10.3±0.4</td><td>33.5±0.9</td><td>66.7±0.1</td><td>610.5</td><td>370.7±23.2</td></tr><tr><td>16<sub>1,16</sub> → 15<sub>1,15</sub>,v=0</td><td>147.561701<sup>b</sup></td><td>62.5</td><td>99</td><td>9.0±0.1</td><td>21.3±0.1</td><td>66.8±0.1</td><td>695.6</td><td>205.5±1.8</td></tr><tr><td>18<sub>0,18</sub> → 17<sub>1,17</sub>,v=0</td><td>147.865229<sup>a</sup></td><td>77.1</td><td>111</td><td></td><td></td><td></td><td>28.7</td><td></td></tr><tr><td>10<sub>2,9</sub> → 10<sub>1,10</sub>,v=0</td><td>148.028678<sup>a</sup></td><td>33.6</td><td>63</td><td></td><td></td><td></td><td>11.5</td><td></td></tr><tr><td>12<sub>2,11</sub> → 12<sub>1,12</sub>,v=0</td><td>153.421750</td><td>44.1</td><td>75</td><td>8.5±0.7</td><td>4.5±0.3</td><td>67.9±0.3</td><td>13.4</td><td>40.9±5.9</td></tr><tr><td>16<sub>2,14</sub> → 15<sub>2,13</sub>,v=0</td><td>153.518944<sup>b</sup></td><td>70.9</td><td>99</td><td>11.9±0.7</td><td>26.5±1.0</td><td>67.4±0.2</td><td>687.7</td><td>336.1±32.2</td></tr><tr><td>16<sub>1,15</sub> → 15<sub>1,14</sub>,v=0</td><td>154.724540<sup>b</sup></td><td>65.4</td><td>99</td><td>10.2±0.4</td><td>15.1±0.5</td><td>66.8±0.2</td><td>695.5</td><td>167.2±13.1</td></tr><tr><td>17<sub>0,17</sub> → 16<sub>0,16</sub>,v=0</td><td>158.657428<sup>b</sup></td><td>69.1</td><td>105</td><td>10.2±0.2</td><td>40.3±0.5</td><td>67.2±0.1</td><td>741.2</td><td>436.5±13.1</td></tr><tr><td>19<sub>0,19</sub> → 18<sub>1,18</sub>,v=0</td><td>158.773785</td><td>85.5</td><td>117</td><td>7.3±0.6</td><td>8.3±0.5</td><td>66.8±0.2</td><td>31.1</td><td>65.2±9.7</td></tr><tr><td>26<sub>1,25</sub> → 26<sub>0,26</sub>,v=0</td><td>159.753858<sup>a</sup></td><td>164.7</td><td>159</td><td></td><td></td><td></td><td>25.6</td><td></td></tr><tr><td rowspan="10"><sup>13</sup>CH<sub>2</sub>CHCN</td><td>16<sub>5,11</sub> → 15<sub>5,10</sub></td><td>147.907794<sup>a</sup></td><td>113.7</td><td>99</td><td></td><td></td><td></td><td>630.4</td><td></td></tr><tr><td>16<sub>6,11</sub> → 15<sub>6,10</sub></td><td>147.909047<sup>a</sup></td><td>137.2</td><td>99</td><td></td><td></td><td></td><td>600.30</td><td></td></tr><tr><td>16<sub>7,9</sub> → 15<sub>7,8</sub></td><td>147.927607</td><td>164.8</td><td>99</td><td>9.7±0.5</td><td>4.2±0.2</td><td>66.9±0.2</td><td>564.8</td><td>42.9±3.8</td></tr><tr><td>16<sub>8,9</sub> → 15<sub>8,8</sub></td><td>147.957689<sup>a</sup></td><td>196.7</td><td>99</td><td></td><td></td><td></td><td>523.9</td><td></td></tr><tr><td>16<sub>3,14</sub> → 15<sub>3,13</sub></td><td>147.986686</td><td>79.5</td><td>99</td><td>9.1±0.5</td><td>5.9±0.2</td><td>66.0±0.2</td><td>673.9</td><td>57.9±4.7</td></tr><tr><td>16<sub>9,7</sub> → 15<sub>9,6</sub></td><td>147.996640<sup>a</sup></td><td>232.8</td><td>99</td><td></td><td></td><td></td><td>477.5</td><td></td></tr><tr><td>16<sub>2,14</sub> → 15<sub>2,13</sub></td><td>149.423715</td><td>69.2</td><td>99</td><td>9.8±2.5</td><td>5.8±1.1</td><td>67.4±0.9</td><td>687.6</td><td>61.7±26.7</td></tr><tr><td>17<sub>0,17</sub> → 16<sub>0,16</sub></td><td>154.577621<sup>a</sup></td><td>67.2</td><td>105</td><td></td><td></td><td></td><td>741.2</td><td></td></tr><tr><td>17<sub>1,16</sub> → 16<sub>1,15</sub></td><td>159.954635</td><td>71.2</td><td>105</td><td>7.4±0.4</td><td>8.2±0.3</td><td>66.9±0.2</td><td>739.4</td><td>64.7±5.9</td></tr><tr><td rowspan="4">CH<sub>2</sub><sup>13</sup>CHCN</td><td>14<sub>0,14</sub> → 13<sub>0,13</sub>,v=0</td><td>130.481700<sup>a</sup></td><td>47.2</td><td>87</td><td></td><td></td><td></td><td>610.5</td><td></td></tr><tr><td>16<sub>0,16</sub> → 15<sub>0,15</sub>,v=0</td><td>148.643016<sup>a</sup></td><td>61.1</td><td>99</td><td></td><td></td><td></td><td>697.6</td><td></td></tr><tr><td>16<sub>1,15</sub> → 15<sub>1,14</sub>,v=0</td><td>153.925606</td><td>65.1</td><td>99</td><td>6.9±0.1</td><td>3.9±0.1</td><td>67.4±0.2</td><td>695.5</td><td>29.1±1.2</td></tr><tr><td>17<sub>2,16</sub> → 16<sub>2,15</sub>,v=0</td><td>159.872596<sup>a</sup></td><td>77.59</td><td>105</td><td></td><td></td><td></td><td>731.7</td><td></td></tr><tr><td rowspan="4">CH<sub>2</sub>CH<sup>13</sup>CHN</td><td>14<sub>0,14</sub> → 13<sub>0,13</sub>,v=0</td><td>130.700756<sup>a</sup></td><td>47.3</td><td>87</td><td></td><td></td><td></td><td>610.6</td><td></td></tr><tr><td>16<sub>0,16</sub> → 15<sub>0,15</sub>,v=0</td><td>148.916702<sup>a</sup></td><td>61.1</td><td>99</td><td></td><td></td><td></td><td>697.7</td><td></td></tr><tr><td>16<sub>1,15</sub> → 15<sub>1,14</sub>,v=0</td><td>154.034285</td><td>65.1</td><td>99</td><td>7.6±0.1</td><td>3.9±0.1</td><td>66.6±0.1</td><td>695.5</td><td>31.9±0.8</td></tr><tr><td>17<sub>2,16</sub> → 16<sub>2,15</sub>,v=0</td><td>160.051637<sup>a</sup></td><td>77.8</td><td>105</td><td></td><td></td><td></td><td>731.8</td><td></td></tr><tr><td rowspan="2">CH<sub>2</sub>CHC<sup>15</sup>N</td><td>16<sub>2,14</sub> → 15<sub>2,13</sub></td><td>148.8982<sup>a</sup></td><td>61.9</td><td>33</td><td></td><td></td><td></td><td>223.1</td><td></td></tr><tr><td>17<sub>1,16</sub> → 16<sub>1,15</sub></td><td>159.436449<sup>a</sup></td><td>71.2</td><td>35</td><td></td><td></td><td></td><td>240.4</td><td></td></tr><tr><td rowspan="10">CH<sub>2</sub>CDCN</td><td>14<sub>6,9</sub> → 13<sub>6,8</sub></td><td>130.637285<sup>a</sup></td><td>108.3</td><td>29</td><td></td><td></td><td></td><td>166.4</td><td></td></tr><tr><td>14<sub>5,10</sub> → 13<sub>5,9</sub></td><td>130.661515<sup>a</sup></td><td>89.6</td><td>29</td><td></td><td></td><td></td><td>177.8</td><td></td></tr><tr><td>14<sub>9,5</sub> → 13<sub>9,4</sub></td><td>130.645695<sup>a</sup></td><td>184.7</td><td>29</td><td></td><td></td><td></td><td>119.6</td><td></td></tr><tr><td>14<sub>4,11</sub> → 13<sub>4,10</sub></td><td>130.720066</td><td>74.2</td><td>29</td><td>11.2±0.5</td><td>2.0±0.1</td><td>67.6±0.1</td><td>23.6</td><td>187.1±1.6</td></tr><tr><td>14<sub>4,10</sub> → 13<sub>4,9</sub></td><td>130.725273<sup>a</sup></td><td>74.3</td><td>29</td><td></td><td></td><td></td><td>187.1</td><td></td></tr><tr><td>17<sub>0,16</sub> → 16<sub>0,16</sub></td><td>154.309753<sup>a</sup></td><td>67.4</td><td>35</td><td></td><td></td><td></td><td>246.8</td><td></td></tr><tr><td>17<sub>3,14</sub> → 16<sub>3,13</sub></td><td>159.383849</td><td>83.9</td><td>35</td><td>6.7±0.2</td><td>3.2±0.1</td><td>68.2±0.1</td><td>23.2</td><td>239.7±0.9</td></tr><tr><td colspan="11" style="text-align: center; font-style: italic;">Ethyl cyanide and its isotopologues and isomer</td></tr><tr><td rowspan="10">C<sub>2</sub>H<sub>5</sub>CN</td><td>14<sub>1,14</sub> → 13<sub>0,13</sub>,v=0</td><td>129.768140</td><td>44.8</td><td>29</td><td>9.8±0.8</td><td>15.5±0.9</td><td>65.6±0.3</td><td>15.1</td><td>162.8±22.2</td></tr><tr><td>15<sub>1,15</sub> → 14<sub>1,14</sub>,v=0</td><td>129.795661<sup>b</sup></td><td>51.1</td><td>31</td><td>9.4±0.2</td><td>43.1±0.9</td><td>66.9±0.1</td><td>221.2</td><td>431.7±19.4</td></tr><tr><td>21<sub>3,19</sub> → 21<sub>2,20</sub>,v=0</td><td>130.693882</td><td>109.4</td><td>43</td><td>7.9±0.4</td><td>11.7±0.5</td><td>66.5±0.2</td><td>16.3</td><td>99.1±8.9</td></tr><tr><td>15<sub>0,15</sub> → 14<sub>0,14</sub>,v=0</td><td>130.903902<sup>b</sup></td><td>50.8</td><td>31</td><td>10.5±0.3</td><td>44.9±1.1</td><td>67.5±0.1</td><td>221.6</td><td>503.9±28.2</td></tr><tr><td>17<sub>0,17</sub> → 16<sub>0,16</sub>,v=0</td><td>147.756711<sup>b</sup></td><td>64.6</td><td>35</td><td>11.8±1.1</td><td>31.9±2.0</td><td>65.5±0.4</td><td>251.13</td><td>403.5±63.8</td></tr><tr><td>24<sub>4,20</sub> → 24<sub>3,21</sub>,v=0</td><td>148.293988</td><td>147.1</td><td>49</td><td>7.5±0.7</td><td>7.8±0.6</td><td>66.8±0.3</td><td>21.1</td><td>62.1±10.1</td></tr><tr><td>10<sub>2,9</sub> → 9<sub>1,8</sub>,v=0</td><td>148.362760</td><td>28.1</td><td>21</td><td>9.8±1.9</td><td>6.4±0.9</td><td>68.2±0.7</td><td>5.4</td><td>66.5±22.7</td></tr><tr><td>33<sub>2,31</sub> → 33<sub>1,32</sub>,v=0</td><td>148.374933<sup>a</sup></td><td>249.5</td><td>67</td><td></td><td></td><td></td><td>21.9</td><td></td></tr><tr><td>25<sub>2,24</sub> → 25<sub>1,25</sub>,v=0</td><td>148.677605<sup>a</sup></td><td>143.0</td><td>51</td><td></td><td></td><td></td><td>10.3</td><td></td></tr><tr><td>17<sub>3,14</sub> → 16<sub>3,13</sub>,v=0</td><td>153.272198<sup>b</sup></td><td>75.9</td><td>35</td><td>9.3±0.2</td><td>26.6±0.4</td><td>66.8±0.1</td><td>244.1</td><td>266.5±7.9</td></tr><tr><td>28<sub>3,26</sub> → 28<sub>2,27</sub>,v=0</td><td>153.589524</td><td>184.6</td><td>57</td><td>8.7±0.8</td><td>5.4±0.4</td><td>67.6±0.3</td><td>19.4</td><td>50.4±8.5</td></tr><tr><td>17<sub>1,16</sub> → 16<sub>1,15</sub>,v=0</td><td>154.244297<sup>b</sup></td><td>68.2</td><td>35</td><td>10.5±0.1</td><td>28.9±0.2</td><td>67.9±0.1</td><td>250.6</td><td>324.4±5.9</td></tr><tr><td>21<sub>4,17</sub> → 21<sub>3,18</sub>,v=0</td><td>154.475980</td><td>117.2</td><td>43</td><td>7.6±0.2</td><td>5.9±0.1</td><td>66.9±0.1</td><td>17.6</td><td>48.5±1.9</td></tr><tr><td>26<sub>2,25</sub> → 26<sub>1,26</sub>,v=0</td><td>154.557360</td><td>154.0</td><td>53</td><td>7.1±0.2</td><td>3.5±0.1</td><td>66.4±0.1</td><td>10.4</td><td>26.8±1.0</td></tr><tr><td>20<sub>8,13</sub> → 21<sub>7,14</sub>,v=0</td><td>158.742366</td><td>161.4</td><td>41</td><td>6.3±0.2</td><td>4.8±0.2</td><td>66.5±0.1</td><td>3.3</td><td>32.3±2.4</td></tr><tr><td>17<sub>4,13</sub> → 17<sub>3,14</sub>,v=0</td><td>159.330800</td><td>83.6</td><td>35</td><td>7.3±0.2</td><td>13.2±0.2</td><td>66.8±0.1</td><td>13.6</td><td>102.9±4.1</td></tr><tr><td>18<sub>1,18</sub> → 17<sub>0,17</sub>,v=0</td><td>159.78554<sup>a</sup></td><td>72.2</td><td>37</td><td></td><td></td><td></td><td>21.1</td><td></td></tr><tr><td>18<sub>2,17</sub> → 17<sub>2,16</sub>,v=0</td><td>159.888831<sup>b</sup></td><td>77.6</td><td>37</td><td>9.8±0.1</td><td>53.6±0.5</td><td>66.9±0.1</td><td>263.3</td><td>561.5±11.8</td></tr><tr><td>16<sub>4,12</sub> → 16<sub>3,13</sub>,v=0</td><td>160.050058<sup>a</sup></td><td>76.3</td><td>33</td><td></td><td></td><td></td><td>12.7</td><td></td></tr><tr><td rowspan="3"><sup>13</sup>CH<sub>3</sub>CH<sub>2</sub>CN</td><td>15<sub>2,14</sub> → 14<sub>2,13</sub>,v=0</td><td>130.319847<sup>a</sup></td><td>54.5</td><td>31</td><td></td><td></td><td></td><td>214.4</td><td></td></tr><tr><td>15<sub>8,7</sub> → 14<sub>8,6</sub>,v=0</td><td>131.109474<sup>a</sup></td><td>120.8</td><td>31</td><td></td><td></td><td></td><td>156.3</td><td></td></tr><tr><td>15<sub>6,9</sub> → 14<sub>6,8</sub>,v=0</td><td>131.119334<sup>a</sup></td><td>89.9</td><td>31</td><td></td><td></td><td></td><td>183.4</td><td></td></tr></tbody></table>

<footnote><sup>a</sup> Gaussian fits are not performed for the blended transitions. <sup>b</sup> Optically thick line.</footnote>
<table><caption>**Table A.1. Continued.**</caption><thead><tr><th>Species</th><th>J<sub>K<sub>s</sub>K<sub>s</sub></sub><sup>+</sup>J<sub>K<sub>s</sub>K<sub>s</sub></sub><sup>-</sup></th><th>Frequency (GHz)</th><th>E<sub>u</sub> (K)</th><th>g<sub>up</sub></th><th>ΔV (km s<sup>-1</sup>)</th><th>I<sub>max</sub> (K)</th><th>V<sub>LSR</sub> (km s<sup>-1</sup>)</th><th>Sμ<sup>2</sup> (Debye<sup>2</sup>)</th><th>∫ T<sub>mb</sub>dV (K km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="12">CH<sub>3</sub><sup>13</sup>CH<sub>2</sub>CN</td><td>15<sub>5,10</sub> → 14<sub>5,9, v=0</td><td>131.155447</td><td>77.8</td><td>31</td><td>8.2±0.6</td><td>7.8±0.3</td><td>67.9±0.2</td><td>194.2</td><td>68.06±7.3</td></tr><tr><td>15<sub>11,4</sub> → 14<sub>11,3, v=0</td><td>131.161428</td><td>183.4</td><td>31</td><td>8.1±0.1</td><td>3.2±0.1</td><td>66.7±0.1</td><td>100.9</td><td>28.2±0.7</td></tr><tr><td>15<sub>4,12</sub> → 14<sub>4,11, v=0</td><td>131.233679</td><td>67.9</td><td>31</td><td>8.3±0.5</td><td>4.8±0.1</td><td>66.1±0.1</td><td>202.9</td><td>42.9±3.8</td></tr><tr><td>15<sub>4,11</sub> → 14<sub>4,10, v=0</td><td>131.249984</td><td>67.9</td><td>31</td><td>9.9±0.7</td><td>6.1±0.2</td><td>66.4±0.2</td><td>202.8</td><td>63.8±7.4</td></tr><tr><td>17<sub>6,11</sub> → 16<sub>6,10, v=0</td><td>148.625244<sup>a</sup></td><td>103.3</td><td>35</td><td></td><td></td><td></td><td>216.6</td><td></td></tr><tr><td>17<sub>11,6</sub> → 16<sub>11,5, v=0</td><td>148.647252<sup>a</sup></td><td>197.2</td><td>35</td><td></td><td></td><td></td><td>143.8</td><td></td></tr><tr><td>17<sub>5,13</sub> → 16<sub>5,12, v=0</td><td>148.682220<sup>a</sup></td><td>91.7</td><td>35</td><td></td><td></td><td></td><td>226.1</td><td></td></tr><tr><td>17<sub>4,13</sub> → 16<sub>4,12, v=0</td><td>148.833416<sup>a</sup></td><td>81.5</td><td>35</td><td></td><td></td><td></td><td>233.8</td><td></td></tr><tr><td>18<sub>1,17</sub> → 17<sub>1,16, v=0</td><td>159.10002</td><td>74.2</td><td>37</td><td>7.7±0.3</td><td>8.6±0.1</td><td>66.6±0.1</td><td>260.7</td><td>70.2±3.6</td></tr><tr><td>12<sub>4,8</sub> → 12<sub>3,9, v=0</td><td>160.026237<sup>a</sup></td><td>78.2</td><td>25</td><td></td><td></td><td></td><td>275.6</td><td></td></tr><tr><td rowspan="2">15<sub>0,15</sub> → 14<sub>0,14, v=0</sub></td><td>130.052706<sup>a</sup></td><td>50.5</td><td>31</td><td></td><td></td><td></td><td></td><td>217.6</td><td></td></tr><tr><td>159.009721<sup>a</sup></td><td>77.1</td><td>37</td><td></td><td></td><td></td><td></td><td>258.6</td><td></td></tr><tr><td rowspan="5">CH<sub>3</sub>CH<sub>2</sub><sup>13</sup>CN</td><td>15<sub>0,15</sub> → 14<sub>0,14, v=0</td><td>130.280373<sup>a</sup></td><td>50.5</td><td>31</td><td></td><td></td><td></td><td>217.7</td><td></td></tr><tr><td>17<sub>1,16</sub> → 16<sub>1,15, v=0</td><td>153.487081<sup>a</sup></td><td>67.8</td><td>35</td><td></td><td></td><td></td><td>246.1</td><td></td></tr><tr><td>17<sub>2,15</sub> → 16<sub>2,14, v=0</td><td>154.609607<sup>a</sup></td><td>70.7</td><td>35</td><td></td><td></td><td></td><td>244.2</td><td></td></tr><tr><td>18<sub>1,18</sub> → 17<sub>1,17, v=0</td><td>154.679342</td><td>71.8</td><td>37</td><td>7.7±0.7</td><td>3.1±0.1</td><td>66.9±0.2</td><td>261.1</td><td>25.9±3.6</td></tr><tr><td>18<sub>2,17</sub> → 17<sub>2,16, v=0</td><td>159.105544<sup>a</sup></td><td>77.2</td><td>37</td><td></td><td></td><td></td><td>258.6</td><td></td></tr><tr><td rowspan="2">C<sub>2</sub>H<sub>5</sub>C<sup>15</sup>N</td><td>15<sub>10,5</sub> → 14<sub>10,4</sub></td><td>130.567344<sup>a</sup></td><td>161.2</td><td>31</td><td></td><td></td><td></td><td>123.5</td><td></td></tr><tr><td>18<sub>2,16</sub> → 17<sub>2,15</sub></td><td>159.848501<sup>a</sup></td><td>76.8</td><td>37</td><td></td><td></td><td></td><td>236.6</td><td></td></tr><tr><td rowspan="2">C<sub>2</sub>H<sub>5</sub>NC</td><td>16<sub>1,16</sub> → 15<sub>1,15, v=0</sub></td><td>149.120363</td><td>62.2</td><td>33</td><td>11.9±0.1</td><td>4.8±0.1</td><td>69.6±0.1</td><td>228.7</td><td>61.2±1.1</td></tr><tr><td>17<sub>0,17</sub> → 16<sub>0,16, v=0</sub></td><td>159.061678</td><td>69.6</td><td>35</td><td>8.4±0.1</td><td>2.9±0.1</td><td>69.6±0.1</td><td>243.3</td><td>26.7±0.7</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">Cyanocacetylene and its isotopologue</td></tr><tr><td>HC<sub>3</sub>N</td><td>17→16, v=0</td><td>154.657284</td><td>66.8</td><td>35</td><td>13.6±0.2</td><td>36.4±0.4</td><td>66.5±0.2</td><td>236.7</td><td>528.8±23.5</td></tr><tr><td>HC<sup>13</sup>CCN</td><td>17→16, v=0</td><td>154.001217</td><td>66.5</td><td>35</td><td>8.3±0.1</td><td>20.9±0.2</td><td>67.0±0.03</td><td>236.8</td><td>183.4±2.5</td></tr><tr><td>HCC<sup>13</sup>CN</td><td>17→16, v=0</td><td>154.016078</td><td>66.5</td><td>35</td><td>9.3±0.1</td><td>22.38±0.1</td><td>66.9±0.02</td><td>236.8</td><td>220.9±3.1</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">Cyanodiacetylene and its isotopologue</td></tr><tr><td>HC<sub>5</sub>N</td><td>49→48</td><td>130.456437</td><td>156.5</td><td>297</td><td>9.2±0.4</td><td>7.0±0.1</td><td>68.9±0.2</td><td>2755.8</td><td>68.8±5.4</td></tr><tr><td>HCC<sup>13</sup>CCCN</td><td>49→48</td><td>130.321350</td><td>156.4</td><td>99</td><td>8.3±0.2</td><td>8.2±0.1</td><td>69.8±0.1</td><td>918.7</td><td>72.1±3.2</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">Cyanamide</td></tr><tr><td rowspan="6">NH<sub>2</sub>CN</td><td>8<sub>1,8</sub> → 7<sub>1,7, v=0</sub></td><td>158.891146</td><td>48.8</td><td>51</td><td>11.6±1.2</td><td>9.4±0.3</td><td>68.1±0.2</td><td>442.2</td><td>116.2±16.2</td></tr><tr><td>8<sub>2,7</sub> → 7<sub>2,6, v=1</sub></td><td>159.595146<sup>a</sup></td><td>162.3</td><td>51</td><td></td><td></td><td></td><td>401.8</td><td></td></tr><tr><td>8<sub>2,6</sub> → 7<sub>2,5, v=1</sub></td><td>159.607737<sup>a</sup></td><td>162.3</td><td>51</td><td></td><td></td><td></td><td>401.8</td><td></td></tr><tr><td>8<sub>3,6</sub> → 7<sub>3,5, v=1</sub></td><td>159.749861<sup>a</sup></td><td>232.0</td><td>17</td><td></td><td></td><td></td><td>125.2</td><td></td></tr><tr><td>8<sub>0,8</sub> → 7<sub>0,7, v=1</sub></td><td>159.814663<sup>a</sup></td><td>105.8</td><td>51</td><td></td><td></td><td></td><td>438.3</td><td></td></tr><tr><td>8<sub>3,6</sub> → 7<sub>3,5, v=0</sub></td><td>160.104805</td><td>164.9</td><td>51</td><td>10.2±0.5</td><td>13.4±0.4</td><td>68.5±0.1</td><td>391.6</td><td>146.3±11.5</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">Aminoacetonitrile</td></tr><tr><td rowspan="6">H<sub>2</sub>NCH<sub>2</sub>CN</td><td>17<sub>7,11</sub> → 16<sub>7,10, v=0</sub></td><td>154.369232<sup>a</sup></td><td>127.1</td><td>35</td><td></td><td></td><td></td><td>93.7</td><td></td></tr><tr><td>17<sub>8,9</sub> → 16<sub>8,8, v=0</sub></td><td>154.373384<sup>a</sup></td><td>145.5</td><td>35</td><td></td><td></td><td></td><td>87.89</td><td></td></tr><tr><td>17<sub>6,12</sub> → 16<sub>6,11, v=0</sub></td><td>154.382222<sup>a</sup></td><td>111.1</td><td>35</td><td></td><td></td><td></td><td>98.93</td><td></td></tr><tr><td>17<sub>9,8</sub> → 16<sub>9,7, v=0</sub></td><td>154.388904<sup>a</sup></td><td>166.4</td><td>35</td><td></td><td></td><td></td><td>81.25</td><td></td></tr><tr><td>17<sub>4,14</sub> → 16<sub>4,13, v=0</sub></td><td>154.517480</td><td>86.4</td><td>35</td><td>8.7±0.4</td><td>2.3±0.1</td><td>65.8±0.1</td><td>106.6</td><td>21.3±1.5</td></tr><tr><td>9<sub>6</sub> → 8<sub>6</sub></td><td>153.711520</td><td>296.1</td><td>38</td><td>9.2±0.9</td><td>2.5±0.2</td><td>66.6±0.4</td><td>5.6</td><td>24.0±4.3</td></tr><tr><td>9<sub>5</sub> → 8<sub>5</sub></td><td>153.743800<sup>a</sup></td><td>216.9</td><td>19</td><td></td><td></td><td></td><td>3.5</td><td></td></tr><tr><td>9<sub>4</sub> → 8<sub>4</sub></td><td>153.770224<sup>a</sup></td><td>152.1</td><td>19</td><td></td><td></td><td></td><td>4.0</td><td></td></tr><tr><td>9<sub>3</sub> → 8<sub>3</sub></td><td>153.790769</td><td>101.7</td><td>38</td><td>10.8±0.9</td><td>8.7±0.6</td><td>67.2±0.3</td><td>9.0</td><td>100.3±14.9</td></tr><tr><td>9<sub>2</sub> → 8<sub>2</sub></td><td>153.805457</td><td>65.7</td><td>19</td><td>10.0±0.4</td><td>6.8±0.2</td><td>67.1±0.2</td><td>4.8</td><td>72.5±5.4</td></tr><tr><td>9<sub>1</sub> → 8<sub>1</sub></td><td>153.814273<sup>a</sup></td><td>44.9</td><td>19</td><td></td><td></td><td></td><td>4.9</td><td></td></tr><tr><td>9<sub>0</sub> → 8<sub>0</sub></td><td>153.817211<sup>a</sup></td><td>36.1</td><td>19</td><td></td><td></td><td></td><td>5.0</td><td></td></tr><tr><td rowspan="4">CH<sub>3</sub><sup>13</sup>CCH</td><td>9<sub>3</sub> → 8<sub>3</sub></td><td>153.72701</td><td>101.7</td><td>38</td><td>8.6±0.1</td><td>1.9±0.1</td><td>65.9±0.1</td><td>9.0</td><td>17.7±0.4</td></tr><tr><td>9<sub>2</sub> → 8<sub>2</sub></td><td>153.741680<sup>a</sup></td><td>65.7</td><td>19</td><td></td><td></td><td></td><td>4.8</td><td></td></tr><tr><td>9<sub>1</sub> → 8<sub>1</sub></td><td>153.750470<sup>a</sup></td><td>44.1</td><td>19</td><td></td><td></td><td></td><td>5.0</td><td></td></tr><tr><td>9<sub>0</sub> → 8<sub>0</sub></td><td>153.753370</td><td>36.8</td><td>19</td><td>9.6±0.4</td><td>1.8±0.1</td><td>67.7±0.2</td><td>5.1</td><td>18.8±1.3</td></tr></tbody></table>

<footnote><sup>a</sup> Gaussian fits are not performed for the blended transitions. <sup>b</sup> Optically thick line.</footnote>
<table><caption>**Table A.2.** Summary of the best-fit line parameters obtained using MCMC method considering a source size = 2'' and a V<sub>lsr</sub> = 68 km/s.</caption><thead><tr><th>Species</th><th>Quantum numbers</th><th>Frequency (GHz)</th><th>E<sub>u</sub> (K)</th><th>FWHM (Km s<sup>-1</sup>)</th><th>Best fitted column density (cm<sup>-2</sup>)</th><th>Best fitted T<sub>ex</sub> (K)</th><th>Optical depth (τ)</th></tr></thead><tbody><tr><td colspan="8" style="text-align: center; font-style: italic;">Vinyl cyanide and its isotopologues</td></tr><tr><td rowspan="4">C<sub>2</sub>H<sub>3</sub>CN</td><td>6<sub>2,4</sub> → 6<sub>1,5</sub>,v=0</td><td>130.763576</td><td>18.2</td><td rowspan="4">8.2</td><td rowspan="4">(3.4 ± 0.4) × 10<sup>17</sup></td><td rowspan="4">209.8±18.02</td><td>0.02</td></tr><tr><td>29<sub>2,27</sub> → 29<sub>1,28</sub>,v=0</td><td>131.168737</td><td>209.5</td><td>0.05</td></tr><tr><td>12<sub>2,11</sub> → 12<sub>1,12</sub>,v=0</td><td>153.42175</td><td>44.1</td><td>0.03</td></tr><tr><td>19<sub>0,19</sub> → 18<sub>1,18</sub>,v=0</td><td>158.773785</td><td>85.5</td><td>0.06</td></tr><tr><td rowspan="4"><sup>13</sup>CH<sub>2</sub>CHCN</td><td>16<sub>7,9</sub> → 15<sub>7,8</sub></td><td>147.927607</td><td>164.8</td><td rowspan="4">9.0</td><td rowspan="4">(5.3±0.4) × 10<sup>15</sup></td><td rowspan="4">125.8±13.67</td><td>0.03</td></tr><tr><td>16<sub>3,14</sub> → 15<sub>3,13</sub></td><td>147.986686</td><td>79.5</td><td>0.07</td></tr><tr><td>16<sub>2,14</sub> → 15<sub>2,13</sub></td><td>149.423715</td><td>69.2</td><td>0.07</td></tr><tr><td>17<sub>1,16</sub> → 16<sub>1,15</sub></td><td>159.954635</td><td>71.2</td><td>0.08</td></tr><tr><td colspan="8" style="text-align: center; font-style: italic;">Ethyl cyanide with its isotopologues and isomer</td></tr><tr><td rowspan="8">C<sub>2</sub>H<sub>5</sub>CN</td><td>14<sub>1,14</sub> → 13<sub>0,13</sub>,v=0</td><td>129.768140</td><td>44.8</td><td rowspan="8">7.8</td><td rowspan="8">(1.7 ± 0.10) × 10<sup>17</sup></td><td rowspan="8">104.1±9.9</td><td>0.20</td></tr><tr><td>21<sub>3,19</sub> → 21<sub>2,20</sub>,v=0</td><td>130.693882</td><td>109.4</td><td>0.12</td></tr><tr><td>24<sub>4,20</sub> → 24<sub>3,21</sub>,v=0</td><td>148.293988</td><td>147.1</td><td>0.12</td></tr><tr><td>10<sub>2,9</sub> → 9<sub>1,8</sub>,v=0</td><td>148.36276</td><td>28.1</td><td>0.10</td></tr><tr><td>28<sub>3,26</sub> → 28<sub>2,27</sub>,v=0</td><td>153.589524</td><td>184.6</td><td>0.08</td></tr><tr><td>21<sub>4,17</sub> → 21<sub>3,18</sub>,v=0</td><td>154.47598</td><td>117.2</td><td>0.14</td></tr><tr><td>26<sub>2,25</sub> → 26<sub>1,26</sub>,v=0</td><td>154.55736</td><td>154.0</td><td>0.07</td></tr><tr><td>20<sub>8,13</sub> → 21<sub>7,14</sub>,v=0</td><td>158.742366</td><td>161.4</td><td>0.02</td></tr><tr><td>17<sub>4,13</sub> → 17<sub>3,14</sub>,v=0</td><td>159.3308</td><td>83.6</td><td>0.15</td></tr><tr><td rowspan="5"><sup>13</sup>CH<sub>3</sub>CH<sub>2</sub>CN</td><td>15<sub>5,10</sub> → 14<sub>5,9</sub>,v=0</td><td>131.155447</td><td>77.8</td><td rowspan="5">8.4</td><td rowspan="5">(7.8 ± 1.1) × 10<sup>15</sup></td><td rowspan="5">130.8±23.2</td><td>0.05</td></tr><tr><td>15<sub>11,4</sub> → 14<sub>11,3</sub>,v=0</td><td>131.161428</td><td>183.4</td><td>0.01</td></tr><tr><td>15<sub>4,12</sub> → 14<sub>4,11</sub>,v=0</td><td>131.233679</td><td>67.9</td><td>0.06</td></tr><tr><td>15<sub>4,11</sub> → 14<sub>4,10</sub>,v=0</td><td>131.249984</td><td>67.9</td><td>0.06</td></tr><tr><td>18<sub>1,17</sub> → 17<sub>1,16</sub>,v=0</td><td>159.10002</td><td>74.2</td><td>0.08</td></tr><tr><td colspan="8" style="text-align: center; font-style: italic;">Methylacetylene and its isotopologue</td></tr><tr><td rowspan="3">CH<sub>3</sub>CCH</td><td>9<sub>6</sub> → 8<sub>6</sub></td><td>153.71152</td><td>296.1</td><td rowspan="3">10.0</td><td rowspan="3">(8.7 ± 0.7) × 10<sup>16</sup></td><td rowspan="3">219.1±14.8</td><td>0.015</td></tr><tr><td>9<sub>3</sub> → 8<sub>3</sub></td><td>153.790769</td><td>101.7</td><td>0.06</td></tr><tr><td>9<sub>2</sub> → 8<sub>2</sub></td><td>153.805457</td><td>65.7</td><td>0.04</td></tr></tbody></table>
**Appendix B: Line maps**

The integrated intensity distribution of all observed molecules is shown in Figures B.1, B.2, B.3, and B.4. The integrated intensity is obtained in the velocity range roughly between 58 km/s and 78 km/s.

<figure><img src="image_19.png" /><figcaption>**Fig. B.1.** Moment maps of unblended transition of C<sub>2</sub>H<sub>3</sub>CN and its isotopologues and isomer. Contours are drawn at 5%, 15%, 30%, 50%, and 80% of the maximum intensity (in Jy/beam.km/s noted in each panel) obtained for each transition. The black contour represents 50% of the maximum intensity.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Fig. B.2. Moment maps of transitions of C<sub>2</sub>H<sub>5</sub>CN, its isotopologues, and isomer (same as Figure B.1).</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Fig. B.3. Moment maps of unblended transition of HC<sub>3</sub>N, HC<sup>13</sup>CCN, HC<sub>5</sub>N, HCC<sup>13</sup>CCCN, NH<sub>2</sub>CN, and NH<sub>2</sub>CH<sub>2</sub>CN (same as Figure B.1).</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>Fig. B.4. Moment maps of unblended transition of CH<sub>3</sub>CCH and CH<sub>3</sub><sup>13</sup>CCH (same as Figure B.1).</figcaption></figure>
Appendix C: Modeled line parameters

<figure><img src="image_23.png" /><figcaption>Fig. C.1. MCMC fit of observed, unblended, optically thin transitions of vinyl cyanide and its isotopologues (solid green lines). The blue lines represent the modeled spectra, whereas the observed spectra are shown in black. The vertical red dashed line shows the position of $V_{lsr} = 68$ km/s. The green solid lines represent the LTE spectra of some species by considering the rotational temperature and column density obtained from the rotational diagram analysis (see Fig. 1).</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>Fig. C.2. MCMC-fit spectra (same as Figure C.1) of the observed, unblended, optically thin transitions of ethyl cyanide, its one isotopologue.</figcaption></figure>
<figure><img src="image_25.png" /><figcaption>Fig. C.3. MCMC-fit spectra of CH<sub>3</sub>CCH transitions.</figcaption></figure>
Appendix D: LTE spectra (LTE2)

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

**Fig. D.1.** Black lines represent observed spectra of vinyl cyanide and its isotopologues (¹³CH₂CHCN, CH₂¹³CCN, CH₂CH¹³CN, CH₂CNC¹⁵N, and CH₂CDCN). Unblended transitions are noted in red, whereas blended transitions are given in black. The solid red lines represent the Gaussian fit of the unblended transitions. LTE spectra (LTE2) are shown in blue by considering N(H₂) = 1.35 × 10²⁵ cm⁻², an excitation temperature of 150 K, average source size (average of emitting regions of unblended transitions), and the average FWHM (obtained from the unblended, optically thin transitions) of the species. The column density of the species is varied until estimated good fits are obtained by eye. The red dashed lines show the systematic velocity (V<sub>LSR</sub>) of the source at ~ 68 km s⁻¹. In addition, the name of the species and their respective transitions (in GHz) are provided in each panel.
<figure><img src="image_27.png" /><figcaption>**Fig. D.2.** Gaussian-fit and -modeled LTE spectra (LTE2; same as Figure D.1) of the unblended transitions of ethyl cyanide, its isotopologues ($^{13}$CH$_3$CH$_2$CN, CH$_3$$^{13}$CH$_2$CN, CH$_3$CH$_2$$^{13}$CN, C$_2$H$_5$C$^{15}$N), and one of its isomers (C$_2$H$_5$NC).</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>**Fig. D.3.** Gaussian-fit spectra and LTE spectra (LTE2; same as Figure D.1) of cyanoacetylene, cyanodiacetylene, cyanamide, and aminoacetonitrile.</figcaption></figure>
<figure><img src="image_29.png" /><figcaption>**Fig. D.4.** Gaussian-fit spectra and LTE spectra (LTE2; same as Figure D.1) of CH<sub>3</sub>CCH and its isotopologues.</figcaption></figure>
Appendix E: Extra figures for chemical evolution of selected species in 1D models

<figure><img src="image_30.png" /><figcaption>Fig. E.1. Same as Figure 15, but for initial dust and gas temperature of 20 K.</figcaption></figure>
<figure><img src="image_31.png" /><figcaption>Fig. E.2. Same as Figure 16, but for initial dust temperature of 20 K.</figcaption></figure>