Detection of complex nitrogen-bearing molecule ethyl cyanide towards the hot molecular core G10.47+0.03

Abstract: The studies of the complex organic molecular lines towards the hot molecular cores at millimeter and submillimeter wavelengths provide instructive knowledge about the chemical complexity in the interstellar medium (ISM). We present the detection of the rotational emission lines of the complex nitrogen-bearing molecule ethyl cyanide (C$_{2}$H$_{5}$CN) towards the chemically rich hot molecular core G10.47+0.03 using the Atacama Large Millimeter/Submillimeter Array (ALMA) band 4 observations. The estimated column density of C$_{2}$H$_{5}$CN towards the G10.47+0.03 is (7.7$\pm$0.5)$\times$10$^{16}$ cm$^{-2}$ with the high rotational temperature of 352.9$\pm$66.8 K. The estimated fractional abundance of C$_{2}$H$_{5}$CN with respect to H$_{2}$ towards the G10.47+0.03 is 5.70$\times$10$^{-9}$. We observe that the estimated fractional abundance of C$_{2}$H$_{5}$CN is similar to the existing three-phase warm-up chemical modelling abundance of C$_{2}$H$_{5}$CN. We also discuss the possible formation mechanism of C$_{2}$H$_{5}$CN towards the hot molecular cores, and we claim the barrierless and exothermic radical-radical reaction between CH$_{2}$ and CH$_{2}$CN is responsible for the production of low abundant of C$_{2}$H$_{5}$CN ($\sim$10$^{-9}$) in the grain surface of G10.47+0.03.

\section{Introduction}
\label{sec:intro}
In the ISM, or circumstellar shells, more than 270 interstellar complex and prebiotic molecules were identified at millimeter/submillimeter wavelengths\footnote{\href{https://cdms.astro.uni-koeln.de/classic/molecules}{https://cdms.astro.uni-koeln.de/classic/molecules}}. The identification of the molecular lines in the ISM is a key step in understanding the chemical evolution from simple molecular species to molecules of biological relevance \citep{her09}. The hot molecular cores are known as the early processes of the high-mass star-formation region, which play a key role in increasing the chemical complexity in the ISM \citep{shi21}. At this stage, the complex organic molecules are ejected from the icy surfaces of dust grains \citep{her09}. The hot molecular cores are characterised by their warm temperature ($\gtrsim$100 K), small source size ($\leq$0.1 pc), and high gas density ($\gtrsim$10$^{6}$ cm$^{-3}$) \citep{kur00, van98}. The time-scale of the hot molecular core lies between $\sim$10$^{5}$ years to $\sim$10$^{6}$ years \citep{van98, gar06}. The hot molecular cores contain the high-velocity water (H$_{2}$O) maser emission, which is located near the ultra-compact (UC) H II regions \citep{meh04, man22c}. The hot core phase is identified by chemically rich molecular emission spectrums with many complex organic molecules, including methyl formate (CH$_{3}$CN) and methanol (CH$_{3}$OH) \citep{al17}. These molecules may be created in the hot molecular cores at a high temperature ($\gtrsim$150 K) via endothermic chemical reactions \citep{al17}. Using high spectral and spatial resolution data like ALMA, VLA, etc., we can identify the different complex organic molecules and the spatial distribution of those complex molecules in the hot molecular cores. The identification of disk candidates in the hot molecular cores is very rare, which suggests the correlation between disks and hot molecular core chemistry \citep{al17}. The study of the molecular lines in the disk candidate hot molecular cores can help us to understand the high-mass star-formation process and chemical evaluation on the small physical scales ($\leq$0.05--0.1 pc) \citep{al17}.
\begin{table*}{}

	\centering

	\caption{Observation summary of G10.47+0.03}
	\begin{adjustbox}{width=1.01\textwidth}
		\begin{tabular}{cccccccccccc}
			\hline
			Observation date &	On-source time&Frequency range&Spectral resolution &Sensitivity (10 Km s$^{-1}$)& Angular resolution\\
			(yyyy-dd-mm)     &	(hh:mm)       &(GHz)            & (kHz)                 &(mJy beam$^{-1}$)&($^{\prime\prime}$)\\
			\hline
			2017-03-05       &	01:53.40      &129.50--130.44&1128.91               &	0.61&1.67$^{\prime\prime}$ (14362 AU) 	\\
			--              &   --            &130.09--130.56 &~488.28                            &0.61&\\
			--               &   --           &130.54--131.01 &~488.28     &0.62&\\
			--              &   --            &130.50--131.44&1128.91   &0.62&\\

			\hline
			2017-01-28       &	00:33.06     &147.50--148.43&1128.91               &	0.50&1.52$^{\prime\prime}$ (13072 AU) 	\\
			--              &   --           &147.55--148.01 &~488.28                  &0.50&\\
			--              &   --           &148.00--148.47&~488.28   &0.53&\\
			--               &   --          &148.50--149.43 &1128.91     &0.53&\\

			\hline
			2017-03-06       &	01:03.50      &153.00--153.93&1128.91               &	0.66&1.66$^{\prime\prime}$ (14276 AU)	\\
			--              &   --            &153.36--153.83 &~488.28                            &0.66&\\
			--              &   --            &153.82--154.29 &~488.28   &0.67&\\
			--               &   --           &154.00--154.93 &1128.91      &0.67&\\

			\hline
			2017-03-07       &	00:28.72      &158.49--159.43 &1128.91               &	0.78&1.76$^{\prime\prime}$ (15136 AU) 	\\
			--              &   --           &159.18--159.65 &~488.28                            &0.80&\\
			--              &   --           &159.49--160.43 &1128.91   &0.79&\\
			--               &   --           &159.64--160.11 &~488.28     &0.79&\\

			\hline
		\end{tabular}
	\end{adjustbox}

	\label{tab:data}

\end{table*}

The complex organic molecules that contain --C$\equiv$N functional group play a key role in prebiotic chemistry because this functional group is responsible for the formation of peptides, nucleic acids, amino acids, and nucleobases of DNA and RNA in the Universe \citep{bel09}. The complex nitrogen-bearing molecule ethyl cyanide (C$_{2}$H$_{5}$CN) is also known as propanenitrile, and this is a well-known interstellar molecule that is found in the hot molecular cores in the massive star-formation regions \citep{meh04}. The dipole moments of a-conformer and b-conformer of the asymmetric top molecule C$_{2}$H$_{5}$CN was $\mu_{a}$ = 3.85 D and $\mu_{b}$ = 1.23 D respectively \citep{he74}. The high dipole moment of \ce{C2H5CN} indicated that it exhibits a high intensity and dense rotational spectrum. The \ce{C2H5CN} molecule forms on the dust grains of high density, warm inner region hot molecular cores \citep{meh04}. The emission lines of C$_{2}$H$_{5}$CN were first detected towards the Orion molecular cloud (OMC-1) and Sgr B2 with estimated column densities of 1.8$\times$10$^{14}$ cm$^{-2}$ and 1.6$\times$10$^{14}$ cm$^{-2}$, respectively \citep{jon77}. Later, hundreds of the rotational emission lines of C$_{2}$H$_{5}$CN were detected towards the high-mass star-formation regions Sgr B2, Orion, and W51 \citep{mio97, li01}. The column density of C$_{2}$H$_{5}$CN towards Sgr B2(N) large molecule heimat source or Sgr B2(N-LMH) hot core reaches 1$\times$10$^{17}$ cm$^{-2}$ \citep{mio97}. The higher excited transitional lines of C$_{2}$H$_{5}$CN were also detected in the hot molecular core Orion KL and Sgr B2(N) \citep{dal13, bel13}. The vibrationally excited transition lines of C$_{2}$H$_{5}$CN were also detected from Sgr B2 and in W51 e2 \citep{meh04, dem08}. Earlier, three $^{13}$C isotopologues of ethyl cyanide like $^{13}$CH$_{3}$CH$_{2}$CN, CH$_{3}$$^{13}$CH$_{2}$CN, and CH$_{3}$CH$_{2}$$^{13}$CN were identified from the Orion hot molecular cloud in the frequency range 80--40 GHz and 160--360 GHz \citep{dem07}. The evidence of the C$_{2}$H$_{5}$CN were also found in the low-mass protostars NGC 1333-IRAS 4A, NGC 1333-IRAS 2A, and IRAS 16293--2422 \citep{taq15, caz03}. Recently, \citet{cor14} found evidence of the C$_{2}$H$_{5}$CN in the atmosphere of Titan using the ALMA.

The hot molecular core candidate G10.47+0.03 is known as a UC H II region, which is located at a distance of 8.6 kpc with a luminosity of 5$\times$10$^{5}$ L$_{\odot}$ \citep{san14, man22, man22a, gor20, ce10}. Previously, the submillimeter rotational emission lines of methyl isocyanate (CH$_{3}$NCO), formamide (NH$_{2}$CHO), glycolaldehyde (HOCH$_{2}$CHO), methyl isocyanate (HNCO), ethylene glycol ((CH$_{2}$OH)$_{2}$), methanol (CH$_{3}$OH), propanal (CH$_{3}$CH$_{2}$CHO), methyl formate (CH$_{3}$OCHO), dimethyl ether (CH$_{3}$COCH$_{3}$), acetone (CH$_{3}$OCH$_{3}$), and acetaldehyde (CH$_{3}$CHO) were identified towards G10.47+0.03 \citep{ohi19,suz16, gor20, mondal21}. The rotational emission lines of possible glycine precursor molecules amino acetonitrile ({NH$_{2}$CH$_{2}$CN}), methylamine ({CH$_{3}$NH$_{2}$}), and methenamine ({CH$_{2}$NH}) were also detected from the G10.47+0.03 using the ALMA band 4 and Nobeyama 45 m radio telescope \citep{man22, ohi19, suz16}. Recently, the rotational emission lines of the amide-type molecule cyanamide ({NH$_{2}$CN}) were detected from G10.47+0.03 using the ALMA with an estimated column density of (6.60$\pm$0.1)$\times$10$^{15}$ cm$^{-2}$ \citep{man22a}. The amide-related molecule {NH$_{2}$CN} was previously identified in the ISM as a possible precursor of urea {(NH$_{2}$CONH$_{2}$)} \citep{man22a}.

\begin{figure*}
	\centering
	\scriptsize
	\includegraphics[width=1.0\textwidth]{G10CONTI1.pdf}
	\includegraphics[width=1.0\textwidth]{G10CONTI2.pdf}
	\caption{Submillimeter continuum images of hot molecular core G10.47+0.03 obtained with ALMA band 4 in the frequency range 130.2348 GHz--160.1556 GHz. The contour levels started at 3$\sigma$, where $\sigma$ is the RMS of each continuum image, and the contour levels increased by a factor of $\surd$2. The red circles indicate the synthesised beam of the continuum images. The corresponding synthesised beam size and RMS of all continuum images were presented in Table.~\ref{tab:cont}.}
	\label{fig:continuum}
\end{figure*}

\begin{table*}{}
	\centering
	\caption{Summary of the continuum image properties of G10.47+0.03.
	}
	\begin{adjustbox}{width=1.0\textwidth}
		\begin{tabular}{cccccccccccc}
			\hline
			Observation date&		Frequency&Integrated flux&Peak flux &Beam size &Position angle& RMS\\
			(yyyy-dd-mm) &			(GHz)   &(Jy)    &(Jy beam$^{-1}$)&($^{\prime\prime}$$\times$$^{\prime\prime}$)&($^{\circ}$)&(mJy beam$^{-1}$)\\
			\hline
			2017-03-05&			130.2348&1.368$\pm$0.018&1.089$\pm$0.089&2.386$\times$1.544&71.357&8.669\\
			&130.3207&1.319$\pm$0.016&1.026$\pm$0.079&2.363$\times$1.542&72.412&8.502\\
			&130.7735&1.304$\pm$0.017&1.026$\pm$0.084&2.370$\times$1.525&71.670&5.245 \\
			&130.9490&1.513$\pm$0.018&1.165$\pm$0.089&2.359$\times$1.523&72.154&2.722 \\
			\hline
			2017-01-28 &147.6924&1.962$\pm$0.035&1.378$\pm$0.016&2.027$\times$1.614&70.064&12.816 \\
			&147.9929&2.069$\pm$0.036&1.447$\pm$0.013&2.015$\times$1.601&69.842&14.209 \\
			&148.2379&2.224$\pm$0.033&1.568$\pm$0.015&2.014$\times$1.604&69.471&~~9.749 \\
			&149.0449&2.325$\pm$0.037&1.596$\pm$0.016&2.073$\times$1.604&69.733&16.231 \\
			\hline
			2017-03-06 &			153.4867&2.320$\pm$0.029&1.761$\pm$0.014&2.101$\times$1.383&77.921&4.614 \\
			&153.5975&2.112$\pm$0.027&1.620$\pm$0.018&2.114$\times$1.391&77.891&9.087 \\
			&154.0386&2.085$\pm$0.030&1.585$\pm$0.014&2.173$\times$1.393&80.869&9.924 \\
			&154.4826&2.204$\pm$0.025&1.724$\pm$0.012&2.086$\times$1.392&78.130&2.861 \\
			\hline
			2017-03-07&			158.9425&2.439$\pm$0.039&1.912$\pm$0.019&2.420$\times$1.446&77.480&15.593 \\
			&159.4149&2.396$\pm$0.042&1.873$\pm$0.020&2.465$\times$1.423&75.547&~~7.673 \\
			&159.9283&2.393$\pm$0.040&1.879$\pm$0.018&2.415$\times$1.452&78.002&~~5.158 \\
			&160.1556&2.716$\pm$0.046&2.099$\pm$0.028&2.417$\times$1.446&77.748&15.593 \\
			\hline
		\end{tabular}
	\end{adjustbox}

	\label{tab:cont}

\end{table*}

In this article, we presented the identification of the rotational emission lines of the complex nitrile-bearing molecule ethyl cyanide (C$_{2}$H$_{5}$CN) towards the hot molecular core G10.47+0.03 using the ALMA. The observations and data reduction of the raw data of G10.47+0.03 were presented in Section~\ref{obs}. The result of the identification and spatial distribution of C$_{2}$H$_{5}$CN was shown in Section~\ref{res}. The discussion and conclusion of the detection of the C$_{2}$H$_{5}$CN towards G10.47+0.03 were shown in Section~\ref{dis} and Section~\ref{conclu}.

\begin{figure*}
	\centering
	\includegraphics[width=0.96\textwidth]{ETH1.pdf}
\end{figure*}
\begin{figure*}
	\centering
	\includegraphics[width=0.96\textwidth]{ETH2.pdf}
	\caption{Rotational emission lines of C$_{2}$H$_{5}$CN with different molecular transitions towards the G10.47+0.03 in the frequency ranges of 129.50 GHz--160.43 GHz. The blue spectra indicated the observed submillimeter spectra of G10.47+0.03, while the red synthetic spectra indicated the best fitting Gaussian model over the observed spectra of C$_{2}$H$_{5}$CN. The black lines indicated the peak position of the C$_{2}$H$_{5}$CN transitions.}
	\label{fig:emission}
\end{figure*}

\begin{table*}

	\centering

	\caption{Summary of the Gaussian fitting line properties of the C$_{2}$H$_{5}$CN  towards G10.47+0.03}
	\begin{adjustbox}{width=0.95\textwidth}
		\begin{tabular}{ccccccccccccccccc}
			\hline
			Observed frequency &Transition & $E_{u}$ & $A_{ij}$ &Peak intensity&S$\mu^{2}$&FWHM & V$_{LSR}$ & $\rm{\int T_{mb}dV}$ &Remark\\
			(GHz) &(${\rm J^{'}_{K_a^{'}K_c^{'}}}$--${\rm J^{''}_{K_a^{''}K_c^{''}}}$) &(K)&(s$^{-1}$) &(K)&(Debye$^{2}$)$^{\dagger}$& (km s$^{-1}$) &(km s$^{-1}$)&(K km s$^{-1}$) & \\
			\hline
			129.7681&14(1,14)--13(0,13)&~~44.86&1.33$\times$10$^{-5}$ &~~7.686&~~15.113&6.896$\pm$0.26&68.596$\pm$0.69&~~65.473$\pm$12.39&Blended with \ce{CH3CDO}\\
			129.7957&15(1,15)--14(1,14)&~~51.09&1.82$\times$10$^{-4}$ &23.278&221.206&6.628$\pm$1.29&68.505$\pm$0.29&207.844$\pm$30.39&Partially blended\\
			130.4905&~18(7,12)--19(6,13)$^{*}$&127.98&2.17$\times$10$^{-6}$ &~~4.213&~~~~3.108&6.698$\pm$0.98&68.539$\pm$0.35&~~~~36.852$\pm$9.69&Blended with \ce{HOCO}\\
			130.6939&21(3,19)--21(2,20)&109.42&9.84$\times$10$^{-6}$ &~~4.118&~~16.282&6.414$\pm$0.13&68.409$\pm$0.16&~~~~27.539$\pm$2.96&Blended with \ce{CH3COOH}\\
			130.9039&15(0,15)--14(0,14)&~~50.79&1.87$\times$10$^{-4}$ &21.130&221.586&6.918$\pm$0.29&68.509$\pm$0.22&209.401$\pm$22.69&Non blended\\
			131.0136&30(4,26)--30(3,27)&218.54&1.30$\times$10$^{-5}$ &21.806&~~30.209&6.798$\pm$0.69&68.523$\pm$0.25&234.189$\pm$36.29&Blended with \ce{SO2}\\
			147.7567&17(0,17)--16(0,16)&~~64.57&2.69$\times$10$^{-4}$ &26.124&251.138&6.976$\pm$0.73&68.509$\pm$0.98&121.432$\pm$22.69&Non blended\\
			148.2939&24(4,20)--24(3,21)&147.04&1.63$\times$10$^{-5}$ &~~5.855&~~21.063&6.129$\pm$0.35&68.498$\pm$0.23&~~~~31.233$\pm$2.41&Non blended\\
			148.3628&~~~~~~10(2,9)--9(1,8)&~~28.05&9.75$\times$10$^{-6}$ &~~4.053&~~~~5.384&6.612$\pm$0.57&68.513$\pm$0.16&~~~~23.991$\pm$3.41&Non blended\\
			148.3749&33(2,31)--33(1,32)&249.49&1.24$\times$10$^{-5}$ &10.968&~~21.915&7.012$\pm$1.69&68.596$\pm$0.87&125.118$\pm$18.51&Blended with \ce{CH3NCO}\\
			148.6776&25(2,24)--25(1,25)&143.02&7.74$\times$10$^{-6}$ &~~8.447&~~10.312&6.987$\pm$0.94&68.523$\pm$0.24&122.341$\pm$12.41&Blended with \ce{H2CO}\\
			149.3949&~~~~~~8(2,6)--7(1,7)&~~19.98&7.34$\times$10$^{-6}$ &~~4.115&~~~~3.216&6.703$\pm$0.12&68.498$\pm$0.10&~~~~36.843$\pm$1.36&Non blended\\
			153.2722&17(3,14)--16(3,13)&~~75.97&2.92$\times$10$^{-4}$ &30.143&244.110&6.911$\pm$0.16&68.503$\pm$0.96&237.711$\pm$36.36&Non blended\\
			153.5895&28(3,16)--28(2,27)&184.63&1.44$\times$10$^{-5}$ &~~6.146&~~19.430&6.99$\pm$0.70&68.516$\pm$0.16&~~~~36.137$\pm$2.98&Non blended\\
			154.2422&17(1,26)--16(1,15)&~~68.17&3.06$\times$10$^{-4}$ &32.813&250.591&6.423$\pm$0.12&68.463$\pm$0.69&354.881$\pm$50.39&Non blended\\
			154.4760&21(4,17)--21(3,18)&117.24&1.76$\times$10$^{-5}$ &~~5.352&~~17.647&6.154$\pm$0.12&68.246$\pm$0.94&~~~~36.491$\pm$3.78&Non blended\\
			154.5574&26(2,25)--26(1,26)&154.02&8.40$\times$10$^{-6}$ &~~2.568&~~10.360&6.605$\pm$0.51&68.587$\pm$0.26&~~~~11.946$\pm$0.69&Blended with \ce{c-H$^{13}$CO}\\
			158.7424&20(8,13)--21(7,14)$^{*}$&161.40&3.73$\times$10$^{-6}$ &~~2.199&~~~~3.284&6.605$\pm$0.39&68.506$\pm$0.96&~~~~14.509$\pm$0.14&Non blended\\
			159.3308&17(4,13)--17(3,14)&~~83.61&1.84$\times$10$^{-5}$ &~~6.726&~~13.660&7.225$\pm$0.46&68.501$\pm$0.99&~~~~39.241$\pm$4.69&Non blended\\
			159.7855&18(1,18)--17(0,17)&~~72.24&2.71$\times$10$^{-5}$ &10.364&~~21.145&6.304$\pm$1.29&68.489$\pm$0.63&~~~~83.257$\pm$6.97&Non blended\\
			159.8888&18(2,17)--17(2,16)&~~77.61&3.39$\times$10$^{-4}$ &22.090&263.285&7.887$\pm$0.21&68.498$\pm$0.12&186.139$\pm$35.96&Non blended\\
			160.0501&16(4,12)--16(3,13)&~~76.29&1.84$\times$10$^{-5}$ &~~6.645&~~12.721&7.543$\pm$0.29&68.536$\pm$0.65&~~~~49.218$\pm$6.87&Non blended\\

			\hline
		\end{tabular}
	\end{adjustbox}
	\\
	${*}$ -- Transitions contain double with a frequency difference of less than 100 kHz. The second molecular transition of \ce{C2H5CN} is not shown.\\
	$^{\dagger}$--Values taken from \href{https://splatalogue.online//advanced1.php}{Splatalogue}.

	\label{tab:MOLECULAR DATA}

\end{table*}

\section{Observations and data reductions}
\label{obs}

We have used the cycle 4 archival data of G10.47+0.03, which was observed using the high-resolution Atacama Large Millimeter/submillimeter Array (ALMA) (\#2016.1.00929.S., PI: Ohishi, Masatoshi). The G10.47+0.03 was located at the observed phase centre of ($\alpha,\delta$)$_{\rm J2000}$ = (18:08:38.232, --19:51:50.400).  The observation of G10.47+0.03 was carried out on 28 January 2017, 05 march 2017, 06 march 2017, and 07 march 2017 using the thirty-nine, forty, forty-one, and thirty-nine antennas respectively. During the observations, the minimum baseline of the antennas was 15 m and the maximum baseline of antennas was 331 m. The observations were made with ALMA band 4 receivers with spectral ranges 129.50 GHz--160.43 GHz with the corresponding angular resolution of 1.67$^{\prime\prime}$--1.76$^{\prime\prime}$. The corresponding spectral resolutions of the observed data were 1128.91 kHz and 488.28 kHz. The flux calibrator was taken as J1733--1304, the bandpass calibrator was taken as J1924--2914, and the phase calibrator was taken as J1832--2039. The observation details are shown in Table~\ref{tab:data}.

We used the Common Astronomy Software Application ({\tt CASA 5.4.1}) for the reduction of the ALMA data \citep{m07}. We used the Perley-Butler 2017 flux calibrator model to scale the continuum flux density of the flux calibrator using the CASA task {\tt SETJY} with 5\% accuracy \citep{pal17}. We used the task {\tt MSTRANSFORM} to split the target data with all rest frequencies after the initial data reduction. After the splitting of the target data, we created the continuum images of G10.47+0.03 in all rest frequencies using CASA task {\tt TCLEAN} with {\tt hogbom} deconvolar. Now we perform the continuum subtraction procedure from the UV plane using the task {\tt UVCONTSUB}. After the continuum subtraction, we used the {\tt TCLEAN} task with a Briggs weighting robust parameter of 0.5 to make the spectral data cubes of G10.47+0.03. Finally, we applied the {\tt IMPBCOR} task in CASA for the primary beam pattern correction in the synthesised image. Our data analysis procedure is similar to \citet{man22, man22a}.

\section{Result}
\label{res}
\subsection{Submillimeter continuum emission towards G10.47+0.03}
We presented the submillimeter continuum emission maps towards the hot molecular core G10.47+0.03 at frequencies 130.2348 GHz, 130.3207 GHz, 130.7735 GHz, 130.9490 GHz, 147.6924 GHz, 147.9929 GHz, 148.2379 GHz, 149.0449 GHz, 153.4867 GHz, 153.5975 GHz, 154.0386 GHz, 154.4826 GHz, 158.9425 GHz, 159.4149 GHz, 159.9283 GHz, and 160.1556 GHz in Figure~\ref{fig:continuum}, where the surface brightness colour scale has the unit of Jy beam$^{-1}$. We estimated the integrated flux density, peak flux density, synthesised beam size, position angle, and RMS of G10.47+0.03 to fit a 2D Gaussian over the continuum images using CASA task {\tt IMFIT}. The estimated continuum image properties of G10.47+0.03 is shown in Table.~\ref{tab:cont}. After fitting the 2D Gaussian, we observed that the synthesised beam size of the continuum image of G10.47+0.03 is not sufficient to resolve the continua. Previously, \citet{rol11} detected the submillimeter continuum emission from the G10.47+0.03 between the frequency range of 201--691 GHz with the variation of flux density as 6--95 Jy corresponding to the spectral index 2.8.

\subsection{Identification of C$_{2}$H$_{5}$CN towards G10.47+003}
The observed submillimeter chemically rich spectra of G10.47+0.03 are extracted from the continuum subtracted spectral datacubes of G10.47+0.03 to create a circular region with a diameter of 2.5$^{\prime\prime}$ centred at RA (J2000) = 18$^{h}$08$^{m}$38$^{s}$.232, Dec (J2000) = --19$^\circ$51$^{\prime}$50$^{\prime\prime}$.440. The systematic velocity ($V_{LSR}$) of the G10.47+0.03 was $\sim$67.50--68.50 km s$^{-1}$ \citep{rol11}. We used the {\tt Line Analysis} module in CASSIS for the analysis of the emission spectra and identification of the molecular lines from the submillimeter emission spectra of G10.47+0.03 \citep{vas15}. After the spectral analysis using the {\tt Line Analysis} module, we detected the rotational emission lines of the complex nitrile-bearing molecule C$_{2}$H$_{5}$CN towards G10.47+0.03. We used the Cologne Database for Molecular Spectroscopy (CDMS) \citep{mu05} database to identify the emission lines of C$_{2}$H$_{5}$CN towards G10.47+0.03. We identified a total of twenty-two rotational emission lines of C$_{2}$H$_{5}$CN towards G10.47+0.03 within the observable frequency ranges. Among the twenty-two rotational emission lines, we identified fourteen emission lines of \ce{C2H5CN} in the un-blended conditions and eight emission lines of \ce{C2H5CN} are blended with other nearby molecular transitions. After fitting the Gaussian model over the detected spectra of \ce{C2H5CN} using the Levenberg-Marquardt algorithm in CASSIS, we estimated the full-width half-maximum (FWHM), molecular transitions ({${\rm J^{'}_{K_a^{'}K_c^{'}}}$--${\rm J^{''}_{K_a^{''}K_c^{''}}}$}), upper state energy ($E_u$), Einstein coefficients ($A_{ij}$), peak intensity and integrated intensity ($\rm{\int T_{mb}dV}$). There are no missing transitions of C$_{2}$H$_{5}$CN in the observable frequency ranges. The summary of the line properties of C$_{2}$H$_{5}$CN is shown in Table~\ref{tab:MOLECULAR DATA}. The detected spectra of C$_{2}$H$_{5}$CN with the best fitting Gaussian model are shown in Figure~\ref{fig:emission}.

\begin{table}{}

	\centering
	\caption{Estimated emitting regions of \ce{C2H5CN} towards G10.47+0.03
	}
	\begin{adjustbox}{width=0.48\textwidth}
		\begin{tabular}{cccccccccccc}
			\hline
			Frequency&Transition&$E_{up}$&Emitting region\\

			(GHz)&(${\rm J^{'}_{K_a^{'}K_c^{'}}}$--${\rm J^{''}_{K_a^{''}K_c^{''}}}$)&(K)&[$^{\prime\prime}$]\\
			\hline
			130.9039&15(0,15)--14(0,14)&~~50.79&1.523 \\
			147.7567&17(0,17)--16(0,16)&~~64.57&1.496 \\
			148.2939&24(4,20)--24(3,21)&147.04&1.510\\
			148.3628&~~~~~~10(2,9)--9(1,8)&~~28.05&1.546 \\
			149.3949&~~~~~~~~8(2,6)--7(1,7)&~~19.98&1.398\\
			153.2722&17(3,14)--16(3,13)&~~75.97&1.364\\
			153.5895&28(3,16)--28(2,27)&184.63&1.482\\
			154.2422&17(1,26)--16(1,15)&~~68.17&1.527\\
			154.4760&21(4,17)--21(3,18)&117.24&1.407\\
			158.7424&20(8,13)--21(7,14)&161.40&1.329\\
			159.3308&17(4,13)--17(3,14)&~~83.61&1.494\\
			159.7855&18(1,18)--17(0,17)&~~72.24&1.524\\
			159.8888&18(2,17)--17(2,16)&~~77.61&1.472\\
			160.0501&16(4,12)--16(3,13)&~~76.29&1.394\\
			\hline
		\end{tabular}
	\end{adjustbox}

	\label{tab:spatial}

\end{table}

\begin{figure*}
	\centering
	\includegraphics[width=1.0\textwidth]{map.pdf}
	\caption{ Integrated emission maps of unblended transitions of C$_{2}$H$_{5}$CN towards G10.47+0.03. The integrated emission map of C$_{2}$H$_{5}$CN was overlaid with the 130.773 GHz (2.29 mm) continuum emission map (red contour). The contour levels are at 20\%, 40\%, 60\%, and 80\% of the continuum peak flux. The blue circle indicated the synthesised beam of the C$_{2}$H$_{5}$CN integrated emission map.}
	\label{fig:map}
\end{figure*}

\subsection{Spatial distribution of \ce{C2H5CN}}
We created the integrated emission maps of unblended transitions of C$_{2}$H$_{5}$CN towards G10.47+0.03. To create the integrated emission maps, we used the task {\tt IMMOMENTS} in CASA. During the run of the task {\tt IMMOMENTS}, we defined the channel range of the data cubes where the emission lines of C$_{2}$H$_{5}$CN were detected. The integrated emission maps of C$_{2}$H$_{5}$CN is shown in Figure~\ref{fig:map}, which was overlaid on the 2.29 mm continuum emission map of G10.47+0.03. After the overlaid continuum emission map over the integrated emission maps of C$_{2}$H$_{5}$CN, we observed that the emission maps of C$_{2}$H$_{5}$CN have a peak at the position of the continuum. The resultant integrated emission maps indicated that the different transitions of C$_{2}$H$_{5}$CN molecule arises from the highly dense warm inner hot core region of G10.47+0.03 due to the different upper state energy ($E_{up}$) of all unblended transitions. To estimate the emitting regions of C$_{2}$H$_{5}$CN, we applied the task {\tt IMFIT} for fitting the 2D Gaussian over the integrated emission maps of C$_{2}$H$_{5}$CN. The deconvolved synthesised beam size of the C$_{2}$H$_{5}$CN emitting regions was estimated by the following equation:\\

\begin{equation}
	\theta_{S}=\sqrt{\theta^2_{50}-\theta^2_{\text{beam}}}
\end{equation}
where, $\theta_{\text{beam}}$ was the half-power width of the synthesised beam of the C$_{2}$H$_{5}$CN integrated emission maps and $\theta_{50} = 2\sqrt{A/\pi}$ was the diameter of the circle whose area ($A$) was surrounding $50\%$ line peak of C$_{2}$H$_{5}$CN \citep{man22c, man22, man22a}. The estimated emitting regions of C$_{2}$H$_{5}$CN is shown in Table~\ref{tab:spatial}. The emitting region of C$_{2}$H$_{5}$CN was observed in the range of 1.329$^{\prime\prime}$--1.546$^{\prime\prime}$. We noticed that the estimated emitting regions of C$_{2}$H$_{5}$CN are comparable to or slightly greater than the synthesised beam size of the integrated emission maps, which indicates the observed transitions of C$_{2}$H$_{5}$CN were not well spatially resolved or, at best, marginally resolved. As a result, no conclusions can be drawn about the morphology of the spatial distribution of \ce{C2H5CN} towards the G10.47+0.03.

\begin{figure*}
	\includegraphics[width=0.48\textwidth]{ROTD.pdf}
	\caption{Rotational diagram of \ce{C2H5CN} towards G10.47+0.03. The blue-filled squares indicated the optically thin approximation data points and blue vertical lines represented the error bars. The best-fit column density and rotational temperature were mentioned inside the image.}
	\label{fig:rotd}
\end{figure*}

\subsection{Rotational diagram analysis of \ce{C2H5CN}}
\label{rotd}
We used the rotational diagram model for the estimation of the column density and rotational temperature of the detected unblended rotational emission lines of \ce{C2H5CN}. Initially, we considered that the detected molecular spectrums of \ce{C2H5CN} were optically thin and obeyed the local thermodynamic equilibrium (LTE) conditions. The LTE approximation is appropriate in the G10.47+0.03 environment because the dust optical depth ($\tau$) of G10.47+0.03 was 0.135, indicating that G10.47+0.03 was optically thin in the submillimeter wavelengths \citep{gor20}. The column density of the optically thin molecular emission lines can be written as \citep{gold99},

\begin{equation}
	{N_u^{thin}}=\frac{3{g_u}k_B\int{T_{mb}dV}}{8\pi^{3}\nu S\mu^{2}}
\end{equation}
where, $g_u$ indicated the degeneracy of the upper energy ($E_u$), $k_B$ represented the Boltzmann constant, $\rm{\int T_{mb}dV}$ indicated the integrated intensity, $\mu$ was the electric dipole moment, $S$ denoted the strength of the detected molecular emission lines, and $\nu$ was the rest frequency of identified emission lines of \ce{C2H5CN}. The equation for total column density under LTE conditions is as follows:

\begin{equation}
	\frac{N_u^{thin}}{g_u} = \frac{N_{\text{total}}}{Q(T_{\text{rot}})}\exp(-E_u/k_BT_{\text{rot}})
\end{equation}
where ${Q(T_{\text{rot}})}$ denotes the partition function as a function of rotational temperature ($T_{\text{rot}}$) and $E_{u}$ denotes the identified molecule's upper energy. The partition function of \ce{C2H5CN} at 75 K, 150 K, and 300 K is 4667.9361, 13209.5867, and 37424.5763, respectively. Equation 4 can be rewritten as,

\begin{equation}

	ln\left(\frac{N_u^{thin}}{g_u}\right) = ln\left(\frac{N}{Q}\right)-\left(\frac{E_u}{k_BT_{rot}}\right)
\end{equation}

Equation 4 indicated a linear connection between the values of $\ln(N_{u}/g_{u}$) and E$_{u}$ of the detected emission lines of \ce{C2H5CN}. The values $\ln(N_{u}/g_{u}$) were calculated from the equation 2. We used the {\tt Rotational Diagram} module in CASSIS to draw the rotational diagram from the parameters of the spectral line, which were estimated by fitting the Gaussian model over the observed transitions of \ce{C2H5CN} using the Levenberg-Marquardt algorithm in CASSIS. The spectral line parameters of \ce{C2H5CN} are shown in Table.~\ref{tab:MOLECULAR DATA}. The column density and rotational temperature of \ce{C2H5CN} can be estimated by fitting a straight line over the values of $ln$(N$_{u}$/g$_{u}$), which is plotted as a function of E$_{u}$. After fitting the straight line, the rotational temperature was estimated from the inverse of the slope, and the column density was estimated from the intercept of the slope. During rotational diagram analysis, we used a total of fourteen non-blended transition lines of \ce{C2H5CN} because the rotational diagram method was applicable only for non-blended rotational molecular emission lines. The generated rotational diagram of \ce{C2H5CN} using the {\tt Rotational Diagram} module in CASSIS was shown in Figure~\ref{fig:rotd}. The blue error bars in Figure~\ref{fig:rotd} represent the error bar of $\ln(N_{u}/g_{u}$), which is calculated from the estimated error of $\rm{\int T_{mb}dV}$ by fitting a Gaussian model over the observed emission lines of \ce{C2H5CN}. The obtained column density of \ce{C2H5CN} using the rotational diagram model was (9.5$\pm$0.1)$\times$10$^{16}$ cm$^{-2}$ with a rotational temperature of 223.8$\pm$4.3 K. Our estimated rotational temperature indicated that the emission lines of \ce{C2H5CN} originated from the hot core region of G10.47+0.03 because \citet{rol11} claimed the temperature of the hot core of G10.47+0.03 was above 100 K. The fractional abundance of H$_{2}$CNCN with respect to {H$_{2}$} was calculated to be 1.90$\times$10$^{-8}$, and the hydrogen column density towards G10.47+0.03 was calculated to be 5.0$\times$10$^{24}$ cm$^{-2}$ \citep{rol11}.

\section{Discussion}
\label{dis}
\subsection{Comparison with observation and simulated abundances of C$_{2}$H$_{5}$CN}
We compared the estimated abundance of \ce{C2H5CN} with the three-phase warm-up chemical modelling results of \citet{gar13}, which was applied in the environment of the hot molecular cores. During the chemical simulation, \citet{gar13} assumed there would be an isothermal collapse phase after a static warm-up phase. In the first phase, the gas density increased from 3$\times$10$^{3}$ to 10$^{7}$ cm$^{-3}$ under the free-fall collapse, and the dust temperature decreased from 16 K to 8 K. In the second phase, the gas density remains constant at $\sim$10$^{7}$ cm$^{-3}$ but the dust temperature fluctuates from 8 K to 400 K. The gas density ($n_{H}$) of G10.47+0.03 was 7$\times$10$^{7}$ cm$^{-3}$ \citep{rol11} and temperature of the warm region was $\sim$150 K \citep{rol11, man22, man22a}. So, the three-phase warm-up chemical simulation of \cite{gar13}, which is applied towards hot molecular cores, is suitable for understanding the chemical evolution of \ce{C2H5CN} towards G10.47+0.03. In the three-phase warm-up simulation, \citet{gar13} used the fast (5$\times$10$^{4}$--7.12$\times$10$^{4}$ years), medium (2$\times$10$^{5}$--2.85$\times$10$^{5}$ years), and slow (1$\times$10$^{6}$--1.43$\times$10$^{6}$ years) warm-up models based on the time scales. \cite{gar13} estimated the abundance of \ce{C2H5CN} based on a three-phase warm-up model. After the chemical simulation, \cite{gar13} estimated that the abundance of \ce{C2H5CN} was 5.7$\times$10$^{-9}$ for first warm-up model, 7.7$\times$10$^{-8}$ for medium warm-up model, and 3.0$\times$10$^{-8}$ for slow warm-up model.  We estimated that the abundance of \ce{C2H5CN} towards G10.47+0.03 was 1.90$\times$10$^{-8}$, which is nearly similar to the simulated abundance of \ce{C2H5CN} in the slow warm-up model. During the comparison of the observed abundances and simulated abundances of \ce{C2H5CN}, we noticed the tolerance factor $\sim$1.5 which is a reasonably good agreement between the observational and theoretical results. Earlier, \citet{man22, man22c} and \citet{ohi19} also claimed that the three-phase warm-up chemical simulation of \citet{gar13} is sufficient to understand chemical environment of G10.47+0.03.

\subsection{Possible formation mechanism of C$_{2}$H$_{5}$CN in the hot molecular cores}
The complex nitrogen-bearing molecule \ce{C2H5CN} was produced on the grain surface of the hot molecular cores \citep{meh04, gar17, gar13, gar22}. In the first stage (free-fall collapse phase), the peak abundance of \ce{C2H5CN} was $\sim$10$^{-9}$ \citep{gar13}. In the free-fall collapse stage, the addition of radical \ce{CH2} with the radical \ce{CH2CN} produced the radical \ce{CH2CH2CN} and again hydrogenetion of the radical \ce{CH2CH2CN} produced the low abundance of \ce{C2H5CN} ($\sim$10$^{-9}$) in the grain surface of hot molecular cores \citep{gar13}. The chemical reaction is as follows:\\\\
\ce{CH2}+\ce{CH2CN}$\rightarrow$\ce{CH2CH2CN}~~~~~~~~~~~~(1) \\\\
\ce{CH2CH2CN}+H$\rightarrow$\ce{C2H5CN}~~~~~~~~~~~~~~(2) \\\\
Reactions 1 indicated that the addition of radical \ce{CH2} and radical \ce{CH2CN} (alternatively, radical-radical reactions) is barrierless and exothermic \citep{sin21}. In the second stage (warm-up phase), the abundance of \ce{C2H5CN} are enhanced at $\sim$10$^{-8}$ at high temperature ($\sim$40--80 K) through the surface chemistry and release of HCN in the gas-phase \citep{gar13}. Our estimated abundance of \ce{C2H5CN} towards G10.47+0.03 was 1.90$\times$10$^{-8}$, which indicated that reaction 1 and 2 is responsible for the production of the \ce{C2H5CN} in the grain surface of G10.47+0.03 during the warm-up phases. \cite{gar17} and \cite{gar22} also used that reaction to simulate the abundance of \ce{C2H5CN} in the environment of hot molecular cores in the free-fall collapse and warm-up phases.

\section{Conclusion}
\label{conclu}
We analysed the ALMA band 4 data of the hot molecular core G10.47+0.03 and we extracted the sub-millimeter rotational molecular lines. The main conclusions of this work are as follows: \\\\
1. We detected a total of twenty-two rotational emission lines of \ce{C2H5CN} towards G10.47+0.03 using the ALMA band 4 observation. \\\\
2. The estimated column density of \ce{C2H5CN} towards G10.47+0.03 was 9.50$\pm$0.1$\times$10$^{16}$ cm$^{-2}$ with a rotational temperature of 223.8$\pm$4.3 K. The derived abundance of \ce{C2H5CN} with respect to \ce{H2} towards G10.47+0.03 was 1.90$\times$10$^{-8}$. \\\\
3. We also compared our estimated abundance of \ce{C2H5CN} with the three-phase warm-up models of \cite{gar13} which were applied towards particularly hot molecular cores. After the comparison, we found that the derived abundance of \ce{C2H5CN} is nearly similar to the simulated abundance of \ce{C2H5CN} under the slow warm-up conditions. \\\\
4. We also discuss the possible formation mechanism of \ce{C2H5CN} towards G10.47+0.03 and we claimed the barrierless and exothermic radical-radical reaction between \ce{CH2} and \ce{CH2CN} is responsible for the production of the highly abundant \ce{C2H5CN} towards G10.47+0.03 during the warm-up phases.

\section*{ACKNOWLEDGEMENTS}{ The plots within this paper and other findings of this study are available from the corresponding author upon reasonable request 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 co-operation with the Republic of Chile. Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ.}

\bibliographystyle{aasjournal}

Astrophysics of Galaxies (astro-ph.GA); Chemical Physics (physics.chem-ph)