We have analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) cycle 2 data of band 6 toward the G345.4938+01.4677 massive young protostellar object (G345.5+1.47 MYSO) in the IRAS 16562--3959 high-mass star-forming region with an angular resolution of $\sim 0.3"$, corresponding to $\sim 760$ au. We spatially resolve the central region which consists of three prominent molecular emission cores. A hypercompact (HC) H$_{\rm {II}}$ region (Core A) and two molecule-rich cores (Core B and Core C) are identified using the moment zero images of the H30$\alpha$ line and a CH$_{3}$OH line, respectively. Various oxygen-bearing complex organic molecules (COMs), such as (CH$_{3}$)$_{2}$CO and CH$_{3}$OCHO, have been detected toward the positions of Core B and Core C, while nitrogen-bearing species, CH$_{3}$CN, HC$_{3}$N and its $^{13}$C isotopologues, have been detected toward all of the cores. We discuss the formation mechanisms of H$_{2}$CO by comparing the spatial distribution of C$^{18}$O with that of H$_{2}$CO. The $^{33}$SO emission, on the other hand, shows a ring-like structure surrounding Core A, and it peaks on the outer edge of the H30$\alpha$ emission region. These results imply that SO is enhanced in a shock produced by the expanding motion of the ionized region.

The cometary materials are thought to be the reservoir of primitive materials in the solar system. The detection of glycine and CH$_3$NH$_2$ by the ROSINA mass spectrometer in the coma of 67P/Churyumov-Gerasimenko suggests that amino acids and their precursors may have been formed in an early evolutionary phase of the Solar System. We investigate if we can explain "Glycine/CH$_3$NH$_2$" ratio by assuming that glycine and CH$_3$NH$_2$ are inherited from the parentual molecular cloud, and the importance of the concentration process of non-volatile species. We use our astrochemical model representative of physio-chemical evolution during the low-mass (or Solar-type) star-formation to investigate if we can explain the observed "Glycine/CH$_3$NH$_2$" ratio. In addition, we develop a Simplified Cometary Concentration Model (SCCM) to discuss this concentration process with the hypothesis that glycine concentrates on the cometary surface during the numerous passages of perihelion after the formation of the Solar System. The predicted molecular ratio of "Glycine/CH$_3$NH$_2$" by our astrochemical model is less than 10$^{-2}$, while the ROSINA measurement has shown this value to be unity. On the other hand, SCCM results suggest that non-volatile species highly likely concentrate the cometary surface due to its high binding energy. If a small amount of non-volatile species is ejected via comet activity, "Glycine/CH$_3$NH$_2$" ratio would be highly enhanced.

In this work, we present the results of our investigation into the chemistry of Z- and E-cyanomethanimine (HNCHCN), both of which are possible precursors to the nucleobase adenine. Ab initio quantum chemical calculations for a number of reactions with atomic hydrogen were carried out. We find that the reaction H + Z/E-HNCHCN leading both to H-addition as well as H$_2$-abstraction proceed via similar short-range barriers with bimolecular rate coefficients on the order of $\sim10^{-17}$ cm$^{3}$ s$^{-1}$. These results were then incorporated into astrochemical models and used in simulations of the giant molecular cloud G+0.693. The calculated abundances obtained from these models were compared with previous observational data and found to be in good agreement, with a predicted [Z/E] ratio of $\sim3$ - somewhat smaller than the previously derived value of $6.1\pm2.4$. We find that the [Z/E] ratio in our simulations is due mostly to ion-molecule destruction rates driven by the different permanent dipoles of the two conformers. Based on these results, we propose a general rule-of-thumb for estimating the abundances of isomers in interstellar environments.

We present the observational result of a glycine precursor, methylamine (CH$_3$NH$_2$), together with methanol (CH$_3$OH) and methanimine (CH$_2$NH) towards high-mass star-forming regions, NGC6334I, G10.47+0.03, G31.41+0.3, and W51~e1/e2 using ALMA. The molecular abundances toward these sources were derived using the rotational diagram method and compared with our state-of-the-art chemical model. We found that the observed ratio of "CH$_3$NH$_2$/CH$_3$OH" is in between 0.11 and 2.2. We also found that the observed "CH$_3$NH$_2$/CH$_3$OH" ratio agrees well with our chemical model by considering the formation of CH$_3$NH$_2$ on the grain surface via hydrogenation process to HCN. This result clearly shows the importance of hydrogenation processes to form CH$_3$NH$_2$. NGC63343I MM3, where CH$_3$NH$_2$ was not detected in this study and showed "CH$_3$NH$_2$/CH$_3$OH" ratio of less than 0.02, is clearly distinguished from the other cores.

Despite low elemental abundance of atomic deuterium in interstellar medium (ISM), observational evidences suggest that several species in gas-phase and in ices could be heavily fractionated. We explore various aspects of deuterium enrichment by constructing a chemical evolution model in gas and grain phases. Depending on various physical parameters, gas and grains are allowed to interact with each other through exchange of their chemical species. It is known that HCO+ and N2H+ are two abundant gas phase ions in ISM and their deuterium fractionation are generally used to predict degree of ionization in various regions of a molecular cloud. To have a more realistic estimation, we consider a density profile of a collapsing cloud. We present radial distributions of important interstellar molecules along with their deuterated isotopomers. We carry out quantum chemical simulation to study effects of isotopic substitution on spectral properties of these important interstellar species. We calculate vibrational (harmonic) frequency of the most important deuterated species (neutral & ions). Rotational and distortional constants of these molecules are also computed to predict rotational transitions of these species. We compare vibrational (harmonic) and rotational transitions as computed by us with existing observational, experimental and theoretical results. We hope that our results would assist observers in their quest of several hitherto unobserved deuterated species.

Protoplanetary disks are challenging objects for astrochemical models due to strong density and temperature gradients and due to the UV photons 2D propagation. In this paper, we have studied the importance of several model parameters on the predicted column densities of observed species. We considered: 1) 2-phase (gas and homogeneous grains) or 3-phase (gas, surface, and bulk of grains) models, 2) several initial compositions, 3) grain growth and dust settling, and 4) several cosmic-ray ionization rates. Our main result is that dust settling is the most crucial parameter. Including this effect renders the computed column densities sensitive to all the other model parameters, except cosmic-ray ionization rate. In fact, we found almost no effect of this parameter for radii larger than 10 au (the minimum radius studied here) except for N2H+. We also compared all our models with all the column densities observed in the protoplanetary disk around DM Tau and were not able to reproduce all the observations despite the studied parameters. N2H+ seems to be the most sensitive species. Its observation in protoplanetary disks at large radius could indicate enough N2 in the gas-phase (inhibited by the 3-phase model, but boosted by the settling) and a low electron abundance (favored by low C and S elemental abundances).

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

Methyl cyanide (CH3CN) and propyne (CH3CCH) are two molecules commonly used as gas thermometers for interstellar gas. They are detected in several astrophysical environments and in particular towards protostars. Using data of the low-mass protostar IRAS 16293-2422 obtained with the IRAM 30m single-dish telescope, we constrained the origin of these two molecules in the envelope of the source. The line shape comparison and the results of a radiative transfer analysis both indicate that the emission of CH3CN arises from a warmer and inner region of the envelope than the CH3CCH emission. We compare the observational results with the predictions of a gas-grain chemical model. Our model predicts a peak abundance of CH3CCH in the gas-phase in the outer part of the envelope, at around 2000 au from the central star, which is relatively close to the emission size derived from the observations. The predicted CH3CN abundance only rises at the radius where the grain mantle ices evaporate, with an abundance similar to the one derived from the observations.

We have carried out observations in the 42$-$46 and 82$-$103 GHz bands with the Nobeyama 45-m radio telescope, and in the 338.2$-$339.2 and 348.45$-$349.45 GHz bands with the ASTE 10-m telescope toward three high-mass star-forming regions containing massive young stellar objects (MYSOs), G12.89+0.49, G16.86$-$2.16, and G28.28$-$0.36. We have detected HC$_{3}$N including its $^{13}$C and D isotopologues, CH$_{3}$OH, CH$_{3}$CCH, and several complex organic molecules (COMs). Combining our previous results of HC$_{5}$N in these sources, we compare the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratios in the three observed sources. The ratio in G28.28$-$0.36 is derived to be $0.091^{+0.109}_{-0.039}$, which is higher than that in G12.89+0.49 by one order of magnitude, and that in G16.86$-$2.16 by a factor of $\sim 5$. We investigate the relationship between the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratio and the $N$(CH$_{3}$CCH)/$N$(CH$_{3}$OH) ratio. The relationships of the two column density ratios in G28.28$-$0.36 and G16.86$-$2.16 are similar to each other, while HC$_{5}$N is less abundant when compared to CH$_{3}$CCH in G12.89+0.49. These results imply a chemical diversity in the lukewarm ($T \sim 20-30$ K) envelope around MYSOs. Besides, several spectral lines from complex organic molecules, including very-high-excitation energy lines, have been detected toward G12.89+0.49, while the line density is significantly low in G28.28$-$0.36. These results suggest that organic-poor MYSOs are surrounded by a carbon-chain-rich lukewarm envelope (G28.28$-$0.36), while organic-rich MYSOs, namely hot cores, are surrounded by a CH$_{3}$OH-rich lukewarm envelope (G12.89+0.49 and G16.86$-$2.16).

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems.

The study of the chemical evolution of glycine in the interstellar medium is one of challenging topics in astrochemistry. Here, we present the chemical modeling of glycine in hot cores using the state-of-the-art three-phase chemical model NAUTILUS, which is focused on the latest glycine chemistry. For the formation process of glycine on the grain surface, we obtained consistent results with previous studies that glycine would be formed via the reactions of COOH with CH$_2$NH$_2$. However, we will report three important findings regarding the chemical evolution and the detectability of interstellar glycine. First, with the experimentally obtained binding energy from the temperature programmed thermal desorption (TPD) experiment, a large proportion of glycine was destroyed through the grain surface reactions with NH or CH$_3$O radicals before it fully evaporates. As a result, the formation process in the gas phase is more important than thermal evaporation from grains. If this is the case, NH$_2$OH and CH$_3$COOH rather than CH$_3$NH$_2$ and CH$_2$NH would be the essential precursors to the gas phase glycine. Secondly, since the gas phase glycine will be quickly destroyed by positive ions or radicals, early evolutionary phase of the hot cores would be the preferable target for the future glycine surveys. Thirdly, we suggest the possibility that the suprathermal hydrogen atoms can strongly accelerate the formation of COOH radicals from CO$_2$, resulting in the dramatic increase of formation rate of glycine on grains. The efficiency of this process should be investigated in detail by theoretical and experimental studies in the future.

The different spatial distributions of N-bearing and O-bearing species, as is well known towards Orion~KL, is one of the long-lasting mysteries. We conducted a survey observation and chemical modeling study to investigate if the different distributions of O- and N-bearing species are widely recognized in general star-forming regions. First, we report our observational results of complex organic molecules (COMs) with the 45~m radio telescope at the Nobeyama Radio Observatory towards eight star-forming regions. Through our spectral survey ranging from 80 to 108~GHz, we detected CH$_3$OH, HCOOCH$_3$, CH$_3$OCH$_3$, (CH$_3$)$_2$CO, CH$_3$CHO, CH$_3$CH$_2$CN, CH$_2$CHCN, and NH$_2$CHO. Their molecular abundances were derived via the rotation diagram and the least squares methods. We found that N-bearing molecules, tend to show stronger correlations with other N-bearing molecules rather than O-bearing molecules. While G10.47+0.03 showed high fractional abundances of N-bearing species, those in NGC6334F were not so rich, being less than 0.01 compared to CH$_3$OH. Then, the molecular abundances towards these sources were evaluated by chemical modeling with NAUTILUS three-phase gas-grain chemical code. Through the simulations of time evolutions for the abundances of COMs, we suggest that observed correlations of fractional abundances between COMs can be explained by the combination of the different temperature structures inside the hot cores and the different evolutionary phase. Since our modeling could not fully explain the observed excitation temperatures, it is important to investigate the efficiency of grain surface reactions and their activation barriers, and the binding energy of COMs to further promote our understanding.

The protonated form of CO2, HOCO+, is assumed to be an indirect tracer of CO2 in the millimeter/submillimeter regime since CO2 lacks a permanent dipole moment. Here, we report the detection of two rotational emission lines (4 0,4-3 0,3) and (5 0,5-4 0,4) of HOCO+ in IRAS 16293-2422. For our observations, we have used EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. The observed abundance of HOCO+ is compared with the simulations using the 3-phase NAUTILUS chemical model. Implications of the measured abundances of HOCO+ to study the chemistry of CO2 ices using JWST-MIRI and NIRSpec are discussed as well.

Methyl isocyanate (CH3NCO) is one of the important complex organic molecules detected on the comet 67P/Churyumov-Gerasimenko by Rosetta's Philae lander. It was also detected in hot cores around high-mass protostars along with a recent detection in the solar-type protostar IRAS 16293-2422. We propose here a gas-grain chemical model to form CH3NCO after reviewing various formation pathways with quantum chemical computations. We have used NAUTILUS 3-phase gas-grain chemical model to compare observed abundances in the IRAS 16293-2422. Our chemical model clearly indicates the ice phase origin of CH3NCO.

Cyclic-C3HD (c-C3HD) is a singly deuterated isotopologue of c-C3H2, which is one of the most abundant and widespread molecules in our Galaxy. We observed IRAS 16293-2422 in the 3 mm band with a single frequency setup using the EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. We observed seven lines of c-C3HD and three lines of c-C3H2. Observed abundances are compared with astrochemical simulations using the NAUTILUS gas-grain chemical model. Our results clearly show that c-C3HD can be used as an important supplement for studying chemistry and physical conditions for cold environments. Assuming that the size of the protostellar envelope is 3000 AU and same excitation temperatures for both c-C3H2 and c-C3HD, we obtain a deuterium fraction of $14_{-3}^{+4}\%$.

Propargyl alcohol (HC2CH2OH, PA) has yet to be observed in the interstellar medium (ISM) although one of its stable isomers, propenal (CH2CHCHO), has already been detected in Sagittarius B2(N) with the 100-meter Green Bank Telescope in the frequency range $18-26$ GHz. In this paper, we investigate the formation of propargyl alcohol along with one of its deuterated isotopomers, HC2CH2OD(OD-PA), in a dense molecular cloud. Various pathways for the formation of PA in the gas and on ice mantles surrounding dust particles are discussed. We use a large gas-grain chemical network to study the chemical evolution of PA and its deuterated isotopomer. Our results suggest that gaseous HC2CH2OH can most likely be detected in hot cores or in collections of hot cores such as the star-forming region Sgr B2(N). A simple LTE (Local thermodynamic equilibrium) radiative transfer model is employed to check the possibility of detecting PA and OD-PA in the millimeter-wave regime. In addition, we have carried out quantum chemical calculations to compute the vibrational transition frequencies and intensities of these species in the infrared for perhaps future use in studies with the James Webb Space Telescope (JWST).

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection.

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi & Kaifu and the values reported in the review from Agundez & Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN.

We conducted survey observations of a glycine precursor, methanimine or methylenimine (CH$_2$NH), with the NRO 45 m telescope and the SMT telescope towards 12 high-mass and two low-mass star-forming regions in order to increase number of CH$_2$NH sources and to better understand the characteristics of CH2NH sources. As a result of our survey, CH$_2$NH was detected in eight sources, including four new sources. The estimated fractional abundances were ~10$^8$ in Orion KL and G10.47+0.03, while they were ~10$^9$ towards the other sources. Our hydrogen recombination line and past studies suggest that CH$_2$NH-rich sources have less evolved HII regions. The less destruction rates by UV flux from the central star would be contributed to the high CH$_2$NH abundances towards CH$_2$NH-rich sources. Our gas-grain chemical simulations suggest that CH$_2$NH is mostly formed in the gas-phase by neutral-neutral reactions rather than the product of thermal evaporation from the dust surfaces.

The nature of the main sulphur reservoir in star forming regions is a long standing mystery. The observed abundance of sulphur-bearing species in dense clouds is only about 0.1 per cent of the same quantity in diffuse clouds. Therefore, the main sulphur species in star forming regions of the interstellar medium are still unknown. IRAS 16293-2422 is one of the regions where production of S-bearing species is favourable due to its conditions which allows the evaporation of ice mantles. We carried out observations in the 3 mm band towards the solar type protostar IRAS 16293-2422 with the IRAM 30m telescope. We observed a single frequency setup with the EMIR heterodyne 3 mm receiver with an Lower Inner (LI) tuning frequency of 89.98 GHz. Several lines of the complex sulphur species CH3SH were detected. Observed abundances are compared with simulations using the NAUTILUS gas-grain chemical model. Modelling results suggest that CH3SH has the constant abundance of 4e-9 (compared to H2) for radii lower than 200 AU and is mostly formed on the surfaces. Detection of CH3SH indicates that there may be several new families of S-bearing molecules (which could form starting from CH3SH) which have not been detected or looked for yet.

Work on the chemical evolution of pre-biotic molecules remains incomplete since the major obstacle is the lack of adequate knowledge of rate coefficients of various reactions which take place in interstellar conditions. In this work, we study the possibility of forming three pyrimidine bases, namely, cytosine, uracil and thymine in interstellar regions. Our study reveals that the synthesis of uracil from cytosine and water is quite impossible under interstellar circumstances. For the synthesis of thymine, reaction between uracil and :CH2 is investigated. Since no other relevant pathways for the formation of uracil and thymine were available in the literature, we consider a large gas-grain chemical network to study the chemical evolution of cytosine in gas and ice phases. Our modeling result shows that cytosine would be produced in cold, dense interstellar conditions. However, presence of cytosine is yet to be established. We propose that a new molecule, namely, C4N3OH5 could be observable in the interstellar region. C4N3OH5 is a precursor (Z isomer of cytosine) of cytosine and far more abundant than cytosine. We hope that observation of this precursor molecule would enable us to estimate the abundance of cytosine in interstellar regions. We also carry out quantum chemical calculations to find out the vibrational as well as rotational transitions of this precursor molecule along with three pyrimidine bases.

We carry out Monte-Carlo simulation to study deuterium enrichment of interstellar grain mantles under various physical conditions. Based on the physical properties, various types of clouds are considered. We find that in diffuse cloud regions, very strong radiation fields persists and hardly a few layers of surface species are formed. In translucent cloud regions with a moderate radiation field, significant number of layers would be produced and surface coverage is mainly dominated by photo-dissociation products such as, C,CH_3,CH_2D,OH and OD. In the intermediate dense cloud regions (having number density of total hydrogen nuclei in all forms ~ 2 x 10^4 cm^-3), water and methanol along with their deuterated derivatives are efficiently formed. For much higher density regions (~ 10^6 cm^-3), water and methanol productions are suppressed but surface coverage of CO,CO_2,O_2,O_3 are dramatically increased. We find a very high degree of fractionation of water and methanol. Observational results support a high fractionation of methanol but surprisingly water fractionation is found to be low. This is in contradiction with our model results indicating alternative routes for de-fractionation of water. Effects of various types of energy barriers are also studied. Moreover, we allow grain mantles to interact with various charged particles (such as H^+, Fe^+,S^+ and C^+) to study the stopping power and projected range of these charged particles on various target ices.

$H_2$ is the most abundant interstellar species. Its deuterated forms ($HD$ and $D_2$) are also significantly abundant. Huge abundances of these molecules could be explained by considering the chemistry occurring on the interstellar dust. Because of its simplicity, Rate equation method is widely used to study the formation of grain-surface species. However, since recombination efficiency of formation of any surface species are heavily dependent on various physical and chemical parameters, Monte Carlo method would be best method suited to take care of randomness of the processes. We perform Monte Carlo simulation to study the formation of $H_2$, $HD$ and $D_2$ on interstellar ices. Adsorption energies of surface species are the key inputs for the formation of any species on interstellar dusts but binding energies of deuterated species are yet to known with certainty. A zero point energy correction exists between hydrogenated and deuterated species which should be considered while modeling the chemistry on the interstellar dusts. Following some earlier studies, we consider various sets of adsorption energies to study the formation of these species in diverse physical circumstances. As expected, noticeable difference in these two approaches (Rate equation method and Monte Carlo method) is observed for production of these simple molecules on interstellar ices. We introduce two factors, namely, $S_f$ and $\beta$ to explain these discrepancies: $S_f$ is a scaling factor, which could be used to correlate discrepancies between Rate equation and Monte Carlo methods. $\beta$ factor indicates the formation efficiency under various circumstances. Higher values of $\beta$ indicates a lower production efficiency. We found that $\beta$ increases with a decrease in rate of accretion from gas phase to grain phase.

It is long debated if pre-biotic molecules are indeed present in the interstellar medium. Despite substantial works pointing to their existence, pre-biotic molecules are yet to be discovered with a complete confidence. In this paper, our main aim is to study the chemical evolution of interstellar adenine under various circumstances. We prepare a large gas-grain chemical network by considering various pathways for the formation of adenine. Majumdar et al. (2012) proposed that in the absence of adenine detection, one could try to trace two precursors of adenine, namely, HCCN and NH_2CN. Recently Merz et al. (2014), proposed another route for the formation of adenine in interstellar condition. They proposed two more precursor molecules. But it was not verified by any accurate gas-grain chemical model. Neither was it known if the production rate would be high or low. Our paper fills this important gap. We include this new pathways to find that the contribution through this pathways for the formation of Adenine is the most dominant one in the context of interstellar medium. We propose that observers may look for the two precursors (C_3NH and HNCNH) in the interstellar media which are equally important for predicting abundances of adenine. We perform quantum chemical calculations to find out spectral properties of adenine and its two new precursor molecules in infrared, ultraviolet and sub-millimeter region. Our present study would be useful for predicting abundance of adenine.

Methyl acetate (CH_3COOCH_3) has been recently observed by IRAM 30 m radio telescope in Orion though the presence of its deuterated isotopomers is yet to be confirmed. We therefore study the properties of various forms of methyl acetate, namely, CH_3COOCH_3, CH_2DCOOCH_3 and CH_3COOCH_2D. Our simulation reveals that these species could be produced efficiently both in gas as well as in ice phases. Production of methyl acetate could follow radical-radical reaction between acetyl (CH_3CO) and methoxy (CH_3O) radicals. To predict abundances of CH_3COOCH_3 along with its two singly deuterated isotopomers and its two isomers (ethyl formate and hydroxyacetone), we prepare a gas-grain chemical network to study chemical evolution of these molecules. Since gas phase rate coefficients for methyl acetate and its related species were unknown, either we consider similar rate coefficients for similar types of reactions (by following existing data bases) or we carry out quantum chemical calculations to estimate the unknown rate coefficients. For the surface reactions, we use adsorption energies of reactants from some earlier studies. Moreover, we perform quantum chemical calculations to obtain spectral properties of methyl acetate in infrared and sub-millimeter regions. We prepare two catalog files for the rotational transitions of CH_2DCOOCH_3 and CH_3COOCH_2D in JPL format, which could be useful for their detection in regions of interstellar media where CH_3COOCH_3 has already been observed.

We carried out an experiment in order to obtain the InfraRed (IR) spectra of methyl propionate (CH3CH2COOCH3) in astrochemical conditions and present the IR spectra for future identification of this molecule in the InterStellar Medium (ISM). The experimental IR spectrum is compared with the theoretical spectrum and an attempt was made to assign the observed peak positions to their corresponding molecular vibrations in condensed phase. Moreover, our calculations suggest that methyl propionate must be synthesized efficiently within the complex chemical network of the ISM and therefore be present in cold dust grains, awaiting identification.

Main focus of this paper is to explore the possibility of findings two deuterated isotopomers of H2Cl+ (Chloronium) in and around the Interstellar Medium (ISM). Presence of Chloronium ion has recently been confirmed by Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (Neufeld et al., 2012). It observed para-chloronium towards six sources in the Galaxy. Existence of its deuterated isotopomers (HDCl+ & D2Cl+) are till date not discussed in the literature. We find that these deuterated gas phase ions could be destroyed by various ion-molecular reactions, dissociative recombination (DR) and by cosmic rays (CR). We compute all the Ion-molecular (polar) reaction rates by using the parameterized trajectory theory and the Ion-molecular (non-polar) reaction rates by using the Langevin theory. For DR and CR induced reactions, we adopt two well behaved rate formulas. We also include these rate coefficients into our large gas-grain chemical network to study the chemical evolution of these species around the outer edge of the cold dense cloud. In order to study spectral properties of chloronium ion and its two deuterated isotopomers, we have carried out quantum chemical simulations. We calculated ground state properties of these species by employing second order Moller-Plesset perturbation theory (MP2) along with quadruple-zeta correlation consistent (aug-cc-pVQZ) basis set. Infrared and electronic absorption spectra of these species are calculated by using the same level of theory. MP2/aug-cc-pVQZ level of theory is used to report the different spectroscopic constants of these gas phase species. These spectroscopic constants are essential to predict rotational transitions of these species. Our predicted column densities of D2Cl^+, HDCl^+ along with the spectral information may enable their future identification around outer edges of cold dark clouds.

It has long been suggested that CH2CN- might be a carrier of one of the many poorly characterized diffuse interstellar bands. In this paper, our aim is to study various forms of CH2CN in the interstellar medium. Aim of this paper is to predict spectroscopic characteristics of various forms of CH2CN and its deuterated derivatives. Moreover, we would like to model the interstellar chemistry for making predictions for the column densities of such species around dark cloud conditions. A detailed quantum chemical simulations to present the spectral properties of various forms of the CH2CN. MP2 theory along with the aug-CCPVTZ basis set is used to obtain different spectroscopic constants of CH2CN-, CHDCN- and CD2CN- in the gas phase which are essential to predict rotational spectra of these species. We performed quantum chemical calculation to find out energetically the most stable spin states for these species. We have computed IR and electronic absorption spectra for different forms of CH2CN. Moreover, we have also implemented a large gas-grain chemical network to predict the column densities of various forms of the cyanomethyl radical and its related species. In order to mimic physical conditions around a dense cloud region, the variation of the visual extinction parameters are considered with respect to the hydrogen number density of the simulated cloud. Our quantum chemical calculation reveals that the singlet spin state is the most stable form of cyanomethyl anion and its deuterated forms. For the confirmation of the detection of the cyanomethyl anion and its two deuterated forms, namely, CHDCN- and CD2CN-, we present the rotational spectral information of these species in the Appendix. Our chemical model predicts that the deuterated forms of cyanomethyl radicals (specially the anions) are also reasonably abundant around the dense region of the molecular cloud.

Cyanoformaldehyde (HCOCN) molecule has recently been suspected towards the Sagittarius B2(N) by the Green Bank telescope, though a confirmation of this observation has not yet been made. In and around a star forming region, this molecule could be formed by the exothermic reaction between two abundant interstellar species, H2CO and CN. Till date, the reaction rate coefficient for the formation of this molecule is unknown. Educated guesses were used to explain the abundance of this molecule by chemical modeling. In this paper, we carried out quantum chemical calculations to find out empirical rate coefficients for the formation of HCOCN and different chemical properties during the formation of HCOCN molecules. Though HCOCN is stable against unimolecular decomposition, this gas phase molecule could be destroyed by many other means, like: ion-molecular reactions or by the effect of cosmic rays. Ion-molecular reaction rates are computed by using the capture theories. We have also included the obtained rate coefficients into our large gas-grain chemical network to study the chemical evolution of these species in various interstellar conditions. Formation of one of the isotopologue(DCOCN) of HCOCN is also studied. Our study predicts the possibility of finding HCOCN and DCOCN in the ice phase with a reasonably high abundance. In order to detect HCOCN or DCOCN in various interstellar environments, it is necessary to know the spectroscopic properties of these molecules. To this effect, we carried out quantum chemical calculations to find out different spectral parameters of HCOCN for the transition in electronic, infrared and rotational modes. We clearly show how the isotopic substitution (DCOCN) plays a part in the vibrational progressions of HCOCN.

Chemical composition of a molecular cloud is highly sensitive to the physical properties of the cloud. In order to obtain the chemical composition around a star forming region, we carry out a two dimensional hydrodynamical simulation of the collapsing phase of a proto-star. A total variation diminishing scheme (TVD) is used to solve the set of equations governing hydrodynamics. This hydrodynamic code is capable of mimicking evolution of the physical properties during the formation of a proto-star. We couple our reasonably large gas-grain chemical network to study the chemical evolution during the collapsing phase of a proto-star. To have a realistic estimate of the abundances of bio-molecules in the interstellar medium, we include the recently calculated rate coefficients for the formation of several interstellar bio-molecules into our gas phase network. Chemical evolution is studied in detail by keeping grain at the constant temperature throughout the simulation as well as by using the temperature variation obtained from the hydrodynamical model. By considering a large gas-grain network with the sophisticated hydrodynamic model more realistic abundances are predicted. We find that the chemical composition are highly sensitive to the dynamic behavior of the collapsing cloud, specifically on the density and temperature distribution.

One of the stumbling blocks for studying the evolution of interstellar molecules is the lack of adequate knowledge of the rate co-efficients of various reactions which take place in the Interstellar medium and molecular clouds. Some of the theoretical models of rate coefficients do exist in the literature for computing abundances of the complex pre-biotic molecules. So far these have been used to study the abundances of these molecules in space. However, in order to obtain more accurate final compositions in these media, we find out the rate coefficients for the formation of some of the most important interstellar pre-biotic molecules by using quantum chemical theory. We use these rates inside our hydro-chemical model to find out the chemical evolution and the final abundances of the pre-biotic species during the collapsing phase of a proto-star. We find that a significant amount of various pre-biotic molecules could be produced during the collapsing phase of a proto-star. We study extensively the formation these molecules via successive neutral-neutral and radical-radical/radical-molecular reactions. We present the time evolution of the chemical species with an emphasis on how the production of these molecules varies with the depth of a cloud. We compare the formation of adenine in the interstellar space using our rate-coefficients and using those obtained from the existing theoretical models. Formation routes of the pre-biotic molecules are found to be highly dependent on the abundances of the reactive species and the rate coefficients involved in the reactions. Presence of grains strongly affect the abundances of the gas phase species. We also carry out a comparative study between different pathways available for the synthesis of adenine, alanine, glycine and other molecules considered in our network.

We carry out a quantum chemical calculation to obtain the infrared and electronic absorption spectra of several complex molecules of the interstellar medium (ISM). These molecules are the precursors of adenine, glycine & alanine. They could be produced in the gas phase as well as in the ice phase. We carried out a hydro-chemical simulation to predict the abundances of these species in the gas as well as in the ice phase. Gas and grains are assumed to be interacting through the accretion of various species from the gas phase on to the grain surface and desorption (thermal evaporation and photo-evaporation) from the grain surface to the gas phase. Depending on the physical properties of the cloud, the calculated abundances varies. The influence of ice on vibrational frequencies of different pre-biotic molecules was obtained using Polarizable Continuum Model (PCM) model with the integral equation formalism variant (IEFPCM) as default SCRF method with a dielectric constant of 78.5. Time dependent density functional theory (TDDFT) is used to study the electronic absorption spectrum of complex molecules which are biologically important such as, formamide and precursors of adenine, alanine and glycine. We notice a significant difference between the spectra of the gas and ice phase (water ice). The ice could be mixed instead of simple water ice. We have varied the ice composition to find out the effects of solvent on the spectrum. We expect that our study could set the guidelines for observing the precursor of some bio-molecules in the interstellar space.