Hot corino chemistry and warm carbon chain chemistry (WCCC) are driven by gas-grain interactions in star-forming cores: radical-radical recombination reactions to form complex organic molecules (COMs) in the ice mantle, sublimation of CH$_4$ and COMs, and their subsequent gas-phase reactions. These chemical features are expected to depend on the composition of ice mantle which is set in the prestellar phase. We calculated the gas-grain chemical reaction network considering a layered ice-mantle structure in star-forming cores, to investigate how the hot corino chemistry and WCCC depend on the physical condition of the static phase before the onset of gravitational collapse. We found that WCCC becomes more active, if the temperature is lower, or the visual extinction is lower in the static phase, or the static phase is longer. Dependence of hot corino chemistry on the static-phase condition is more complex. While CH$_3$OH is less abundant in the models with warmer static phase, some COMs are formed efficiently in those warm models, since there are various formation paths of COMs. If the visual extinction is lower, photolysis makes COMs less abundant in the static phase. Once the collapse starts and visual extinction increases, however, COMs can be formed efficiently. Duration of the static phase does not largely affect COM abundances. Chemical diversity between prototypical hot corinos and hybrid sources, in which both COMs and carbon chains are reasonably abundant, can be explained by the variation of prestellar conditions. Deficiency of gaseous COMs in prototypical WCCC sources is, however, hard to reproduce within our models.

The deuterium fractionation of water can serve as a tracer for the chemical and physical evolution of water during star formation and can constrain the origin of water in Solar System bodies. We determine the HDO/H$_2$O ratio in the inner warm gas toward three low-mass Class 0 protostars selected to be in isolated cores, i.e., not associated with any cloud complexes. Previous sources for which the HDO/H$_2$O ratio have been established were all part of larger star-forming complexes. Targeting these isolated protostars allows comparison of the water chemistry in isolated and clustered regions to determine the influence of local cloud environment. We present ALMA observations of the HDO $3_{1,2}$-$2_{2,1}$ and $2_{1,1}$-$2_{1,2}$ transitions at 225.897 GHz and 241.562 GHz along with the H$_2^{18}$O $3_{1,3}$-$2_{2,0}$ transition at 203.407 GHz. The high angular resolution (0\farcs3-1\farcs3) allow the study of the inner warm envelope gas. Model-independent estimates for the HDO/H$_2$O ratios are obtained and compared with previous determinations in the warm gas toward low-mass protostars. We detect the targeted water transitions toward the three sources with S/N > 5. We determine the HDO/H$_2$O ratio toward L483, B335 and BHR71-IRS1 to be ($2.2\pm0.4$)$\times 10^{-3}$, ($1.7\pm0.3$)$\times 10^{-3}$, and ($1.8\pm0.4$)$\times 10^{-3}$, respectively, assuming $T_\mathrm{ex} = 124$ K. The degree of water deuteration of these isolated protostars are a factor of 2-4 higher relative to Class 0 protostars that are members of known nearby clustered star-forming regions. The results indicate that the water deuterium fractionation is influenced by the local cloud environment. This effect can be explained by variations in either collapse timescales or temperatures, which depends on local cloud dynamics and could provide a new method to decipher the history of young stars.

We investigate the evolution of the ortho-to-para ratio of overall (gas + ice) H$_2$ via the nuclear spin conversion on grain surfaces coated with water ice under physical conditions that are relevant to star- and planet-forming regions. We utilize the rate equation model that considers adsorption of gaseous H$_2$ on grain surfaces which have a variety of binding sites with a different potential energy depth, thermal hopping, desorption, and the nuclear spin conversion of adsorbed H$_2$. It is found that the spin conversion efficiency depends on the H$_2$ gas density and the surface temperature. As a general trend, enhanced H$_2$ gas density reduces the efficiency, while the temperature dependence is not monotonic; there is a critical surface temperature at which the efficiency is the maximum. At low temperatures, the exchange of gaseous and icy H$_2$ is inefficient (i.e., adsorbed H$_2$ does not desorb and hinders another gaseous H$_2$ to be adsorbed), while at warm temperatures, the residence time of H$_2$ on surfaces is too short for the spin conversion. Additionally, the spin conversion becomes more efficient with lowering the activation barriers for thermal hopping. We discuss whether the spin conversion on surfaces can dominate over that in the gas-phase in star- and planet-forming regions. Finally, we establish a simple but accurate way to implement the H$_2$ spin conversion on grain surfaces in existing astrochemical models.

Sulfur is one of the most abundant elements in the Universe, with important roles in astro-, geo-, and biochemistry. Its main reservoirs in planet-forming disks have previously eluded detection: gaseous molecules only account for $<1\,$\% of total elemental sulfur, with the rest likely in either ices or refractory minerals. Mechanisms such as giant planets can filter out dust from gas accreting onto disk-hosting stars. For stars above 1.4 solar masses, this leaves a chemical signature on the stellar photosphere that can be used to determine the fraction of each element that is locked in dust. Here, we present an application of this method to sulfur, zinc, and sodium. We analyse the accretion-contaminated photospheres of a sample of young stars and find $(89\pm8)\,$\% of elemental sulfur is in refractory form in their disks. The main carrier is much more refractory than water ice, consistent with sulfide minerals such as FeS.

We performed sensitive observations of the N15ND+(1-0) and 15NND+(1-0) lines toward the prestellar core L1544 using the IRAM 30m telescope. The lines are not detected down to 3 sigma levels in 0.2 km/s channels of around 6 mK. The non-detection provides the lower limit of the 14N/15N ratio for N2D+ of ~700-800, which is much higher than the elemental abundance ratio in the local ISM of ~200-300. The result indicates that N2 is depleted in 15N in the central part of L1544, because N2D+ preferentially traces the cold dense gas, and because it is a daughter molecule of N2. In-situ chemistry is unlikely responsible for the 15N depletion in N2; neither low-temperature gas phase chemistry nor isotope selective photodissociation of N2 explains the 15N depletion; the former prefers transferring 15N to N2, while the latter requires the penetration of interstellar FUV photons into the core center. The most likely explanation is that 15N is preferentially partitioned into ices compared to 14N via the combination of isotope selective photodissociation of N2 and grain surface chemistry in the parent cloud of L1544 or in the outer regions of L1544 which are not fully shielded from the interstellar FUV radiation. The mechanism is the most efficient at the chemical transition from atomic to molecular nitrogen. In other words, our result suggests that the gas in the central part of L1544 already went trough the transition from atomic to molecular nitrogen in the earlier evolutionary stage, and that N2 is currently the primary form of gas-phase nitrogen.

Recent measurements carried out at comet 67P/C-G with the ${\it Rosetta}$ probe revealed that molecular oxygen, O$_2$, is the fourth most abundant molecule in comets. Models show that O$_2$ is likely of primordial nature, coming from the interstellar cloud from which our Solar System was formed. However, gaseous O$_2$ is an elusive molecule in the interstellar medium with only one detection towards quiescent molecular clouds, in the $\rho$ Oph A core. We perform a deep search for molecular oxygen, through the $2_1 - 0_1$ rotational transition at 234 GHz of its $^{16}$O$^{18}$O isotopologue, towards the warm compact gas surrounding the nearby Class 0 protostar IRAS 16293--2422 B with the ALMA interferometer. The targeted $^{16}$O$^{18}$O transition is surrounded by two brighter transitions at $\pm 1$ km s$^{-1}$ relative to the expected $^{16}$O$^{18}$O transition frequency. After subtraction of these two transitions, residual emission at a 3$\sigma$ level remains, but with a velocity offset of $0.3 - 0.5$ km s$^{-1}$ relative to the source velocity, rendering the detection "tentative". We derive the O$_2$ column density for two excitation temperatures $T_{\rm ex}$ of 125 and 300 K, as indicated by other molecules, in order to compare the O$_2$ abundance between IRAS16293 and comet 67P/C-G. Assuming that $^{16}$O$^{18}$O is not detected and using methanol CH$_3$OH as a reference species, we obtain a [O$_2$]/[CH$_3$OH] abundance ratio lower than $2-5$, depending on the assumed $T_{\rm ex}$, a three to four times lower abundance than the [O$_2$]/[CH$_3$OH] ratio of $5-15$ found in comet 67P/C-G. Such a low O$_2$ abundance could be explained by the lower temperature of the dense cloud precursor of IRAS16293 with respect to the one at the origin of our Solar System that prevented an efficient formation of O$_2$ in interstellar ices.

Cometary studies suggest that the organic composition of the early Solar Nebula was rich in complex nitrile species such a CH$_3$CN. Recent ALMA detections in protoplanetary disks suggest that these species may be common during planet and comet formation, but connecting gas phase measurements to cometary abundances first requires constraints on formation chemistry and distributions of these species. We present here the detection of seven spatially resolved transitions of CH$_3$CN in the protoplanetary disk around the T-Tauri star TW Hya. Using a rotational diagram analysis we find a disk-averaged column density of N$_T$=1.45$^{+0.19}_{-0.15}\times10^{12}$ cm$^{-2}$ and a rotational temperature of T$_{rot}$=32.7$^{+3.9}_{-3.4}$ K. A radially resolved rotational diagram shows the rotational temperature to be constant across the disk, suggesting that the CH$_3$CN emission originates from a layer at z/r$\sim$0.3. Through comparison of the observations with predictions from a disk chemistry model, we find that grain-surface reactions likely dominate CH$_3$CN formation and that in situ disk chemistry is sufficient to explain the observed CH$_3$CN column density profile without invoking inheritance from the protostellar phase. However, the same model fails to reproduce a Solar System cometary abundance of CH$_3$CN relative to H$_2$O in the midplane, suggesting that either vigorous vertical mixing or some degree of inheritance from interstellar ices occurred in the Solar Nebula.

We have observed complex organic molecules (COMs) and deuterated species toward a hot core/corino (HC) associated with the infrared dark cloud (IRDC) clump G34.43+00.24 MM3 with the Atacama Large Millimeter/submillimeter Array (ALMA). We have detected six normal-COMs (CH3OH, CH3CHO, CH3CH2CN, CH3OCH3, HCOOCH3, and NH2CHO), one deuterated-COM (CH2DCN), and two deuterated fundamental molecules (D2CO and DNC) toward G34.43+00.24 MM3 HC. None of these lines, except for CH3OH, are detected toward the shocked regions in our data, which suggests that COMs do not originate in shocks. The abundance of the COMs relative to CH3OH in G34.43+00.24 MM3 HC is found to be similar to those in high-mass hot cores, rather than those in hot corinos in low-mass star-forming regions. This result suggests that the physical conditions of the warm-up phase of G34.43+00.24 MM3 HC are similar to those of high-mass sources. On the other hand, the D2CO abundance relative to CH$_3$OH in G34.43+00.24 MM3 HC is higher than that of other hot cores, and seems to be comparable to that of hot corinos. The relatively high D2CO/CH3OH ratio of G34.43+00.24 MM3 HC implies a long cold starless phase of G34.43+00.24 MM3 HC.

We investigate nitrogen isotope fractionation in forming and evolving molecular clouds using gas-ice astrochemical simulations. We find that the bulk gas can become depleted in heavy nitrogen (15N) due to the formation of 15N-enriched ices. Around the chemical transition from atomic nitrogen to N2, N15N is selectively photodissociated, which results in the enrichment of 15N in atomic nitrogen. As 15N-enriched atomic nitrogen is converted to ammonia ice via grain surface reactions, the bulk gas is depleted in 15N. The level of 15N depletion in the bulk gas can be up to a factor of two compared to the elemental nitrogen isotope ratio, depending on the photodesorption yield of ammonia ice. Once the nitrogen isotopes are differentially partitioned between gas and solids in a molecular cloud, it should remain in the later stages of star formation (e.g., prestellar core) as long as the sublimation of ammonia ice is inefficient. Our model suggests that all the N-bearing molecules in the cold gas of star-forming regions can be depleted in 15N, which is at least qualitatively consistent with the observations toward prestellar core L1544. In our models, icy species show both 15N and deuterium fractionation. The fractionation pattern within ice mantles is different between 15N and deuterium, reflecting their fractionation mechanisms; while the concentration of deuterium almost monotonically increases from the lower layers of the ice mantles to the upper layers, the concentration of 15N reaches the maximum at a certain depth and declines towards the surface.

We investigate deuterium chemistry coupled with the nuclear spin-state chemistry of H$_2$ and H$_3^+$ in protoplanetary disks. Multiple paths of deuterium fractionation are found; exchange reactions with D atoms, such as HCO$^+$ + D, are effective in addition to those with HD. In a disk model with grain sizes appropriate for dark clouds, the freeze-out of molecules is severe in the outer midplane, while the disk surface is shielded from UV radiation. Gaseous molecules, including DCO$^+$, thus become abundant at the disk surface, which tends to make their column density distribution relatively flat. If the dust grains have grown to millimeter size, the freeze-out rate of neutral species is reduced, and the abundances of gaseous molecules, including DCO$^+$ and N$_2$D$^+$, are enhanced in the cold midplane. Turbulent diffusion transports D atoms and radicals at the disk surface to the midplane, and stable ice species in the midplane to the disk surface. The effects of turbulence on chemistry are thus multifold; while DCO$^+$ and N$_2$D$^+$ abundances increase or decrease depending on the regions, HCN and DCN in the gas and ice are much reduced at the innermost radii, compared with the model without turbulence. When cosmic rays penetrate the disk, the ortho-to-para ratio (OPR) of H$_2$ is found to be thermal in the disk, except in the cold ($\lesssim 10$ K) midplane. We also analyze the OPR of H$_3^+$ and H$_2$D$^+$, as well as the main reactions of H$_2$D$^+$, DCO$^+$, and N$_2$D$^+$ to analytically derive their abundances in the cold midplane.

Partitioning of elemental nitrogen in star-forming regions is not well constrained. Most nitrogen is expected to be partitioned among atomic nitrogen, molecular nitrogen (N2), and icy N-bearing molecules, such as ammonia (NH3) and N2. Atomic nitrogen is not directly observable in the cold gas. In this paper, we propose an indirect way to constrain the amount of atomic nitrogen in the cold gas of star-forming clouds, via deuteration in ammonia ice, the [ND2H/NH2D]/[NH2D/NH3] ratio. Using gas-ice astrochemical simulations, we show that if atomic nitrogen remains as the primary reservoir of nitrogen during cold ice formation stages, the [ND2H/NH2D]/[NH2D/NH3] ratio is close to the statistical value of 1/3 and lower than unity, whereas if atomic nitrogen is largely converted into N-bearing molecules, the ratio should be larger than unity. Observability of ammonia isotopologues in the inner hot regions around low-mass protostars, where ammonia ice has sublimated, is also discussed. We conclude that the [ND2H/NH2D]/[NH2D/NH3] ratio can be quantified using a combination of VLA and ALMA observations with reasonable integration times, at least toward IRAS 16293-2422 where high molecular column densities are expected.

We propose a new simple computational model to estimate adsorption energies of atoms and molecules to low-temperature amorphous water ice, and we present the adsorption energies of carbon (3P), nitrogen (4S), and oxygen (3P) atoms based on quantum chemistry calculations. The adsorption energies were estimated to be 14100 +- 420 K for carbon, 400 +- 30 K for nitrogen, and 1440 +-160 K for oxygen. The adsorption energy of oxygen is well consistent with experimentally reported value. We found that the binding of a nitrogen atom is purely physisorption, while that of a carbon atom is chemisorption in which a chemical bond to an O atom of a water molecule is formed. That of an oxygen atom has a dual character both physisorption and chemisorption. The chemisorption of atomic carbon also implies a possibility of further chemical reactions to produce molecules bearing a C-O bond, while it may hinder the formation of methane on water ice via sequential hydrogenation of carbon atoms. These would be of a large impact to the chemical evolution of carbon species in interstellar environments. We also investigated effects of the newly calculated adsorption energies onto chemical compositions of cold dense molecular clouds with the aid of gas-ice astrochemical simulations. We found that abundances of major nitrogen-bearing molecules, such as N2 and NH3, are significantly altered by applying the calculated adsorption energy, because nitrogen atoms can thermally diffuse on surfaces even at 10 K.

The level of isotopic fractionation in molecules provides insights into their formation environments and how they formed. In this article, we review hydrogen and nitrogen isotopic fractionation in low-mass star-forming regions. Interstellar molecules are significantly enriched in deuterium. The importance of the nuclear spin states of light species on deuterium fractionation and the usefulness of singly and doubly deuterated molecules as chemical tracers are discussed. Observations have revealed that molecules in prestellar cores are enriched in or depleted in 15N depending on molecules. Compared with deuterium fractionation chemistry, our understanding of 15N fractionation chemistry is not well established. We briefly discuss potential 15N fractionation routes, i.e., isotopic-exchange reactions and isotopic selective photodissociation of N2. In addition, the selective freeze-out of 15N atoms onto dust grains around the transition between N atoms and N2 is discussed as a potential mechanism that causes the depletion of 15N in the gas phase.

Molecular oxygen, O$_2$, was recently detected in comet 67P by the ROSINA instrument on board the Rosetta spacecraft with a surprisingly high abundance of 4 % relative to H$_2$O, making O$_2$ the fourth most abundant in comet 67P. Other volatile species with similar volatility, such as molecular nitrogen N$_2$, were also detected by Rosetta, but with much lower abundances and much weaker correlations with water. Here, we investigate the chemical and physical origin of O$_2$ and other volatile species using the new constraints provided by Rosetta. We follow the chemical evolution during star formation with state-of-the-art astrochemical models applied to dynamical physical models by considering three origins: i) in dark clouds, ii) during forming protostellar disks, and iii) during luminosity outbursts in disks. The models presented here favour a dark cloud (or "primordial") grain surface chemistry origin for volatile species in comets, albeit for dark clouds which are slightly warmer and denser than those usually considered as solar system progenitors.

The deuterium enrichment of molecules is sensitive to their formation environment. Constraining patterns of deuterium chemistry in protoplanetary disks is therefore useful for probing how material is inherited or reprocessed throughout the stages of star and planet formation. We present ALMA observations at $\sim0.6"$ resolution of DCO$^+$, H$^{13}$CO$^+$, DCN, and H$^{13}$CN in the full disks around T Tauri stars AS 209 and IM Lup, the transition disks around T Tauri stars V4046 Sgr and LkCa 15, and the full disks around Herbig Ae stars MWC 480 and HD 163296. We also present ALMA observations of HCN in the IM Lup disk. DCN, DCO$^+$, and H$^{13}$CO$^+$ are detected in all disks, and H$^{13}$CN in all but the IM Lup disk. We find efficient deuterium fractionation for the sample, with estimates of disk-averaged DCO$^+$/HCO$^+$ and DCN/HCN abundance ratios ranging from $\sim0.02-0.06$ and $\sim0.005-0.08$, respectively, which is comparable to values reported for other ISM environments. The relative distributions of DCN and DCO$^+$ vary between disks, suggesting that multiple formation pathways may be needed to explain the diverse emission morphologies. In addition, gaps and rings observed in both H$^{13}$CO$^+$ and DCO$^+$ emission provide new evidence that DCO$^+$ bears a complex relationship with the location of the midplane CO snowline.

We investigate the chemistry in a radiation-hydrodynamics model of star-forming core which evolves from a cold ($\sim 10$ K) prestellar core to the main accretion phase in $\sim 10^5$ yr. A rotationally-supported gravitationally unstable disk is formed around a protostar. We extract the temporal variation of physical parameters in $\sim 1.5 \times 10^3$ SPH particles which end up in the disk, and perform post-processing calculations of the gas-grain chemistry adopting a three-phase model. Inside the disk, the SPH particles migrate both inward and outward. Since a significant fraction of volatiles such as CO can be trapped in the water-dominant ice in the three-phase model, the ice mantle composition depends not only on the current position in the disk but also on whether the dust grain has ever experienced higher temperatures than the water sublimation temperature. Stable molecules such as H$_2$O, CH$_4$, NH$_3$ and CH$_3$OH are already abundant at the onset of gravitational collapse and simply sublimated as the fluid parcels migrate inside the water snow line. On the other hand, various molecules such as carbon chains and complex organic molecules (COMs) are formed in the disk. COMs abundance sensitively depends on the outcomes of photodissociation and diffusion rates of photofragments in bulk ice mantle. As for S-bearing species, H$_2$S ice is abundant in the collapse phase. In the warm regions in the disk, H$_2$S is sublimated to be destroyed, while SO, H$_2$CS, OCS and SO$_2$ become abundant.

We investigate the delivery of regular and deuterated forms of water from prestellar cores to circumstellar disks. We adopt a semi-analytical axisymmetric two-dimensional collapsing core model with post-processing gas-ice astrochemical simulations, in which a layered ice structure is considered. The physical and chemical evolutions are followed until the end of the main accretion phase. When mass averaged over the whole disk, a forming disk has a similar H2O abundance and HDO/H2O abundance ratio as their precollapse values (within a factor of 2), regardless of time in our models. Consistent with previous studies, our models suggest that interstellar water ice is delivered to forming disks without significant alteration. On the other hand, the local vertically averaged H2O ice abundance and HDO/H2O ice ratio can differ more, by up to a factor of several, depending on time and distance from a central star. Key parameters for the local variations are the fluence of stellar UV photons en route into the disk and the ice layered structure, the latter of which is mostly established in the prestellar stages. We also find that even if interstellar water ice is destroyed by stellar UV and (partly) reformed prior to disk entry, the HDO/H2O ratio in reformed water ice is similar to the original value. This finding indicates that some caution is needed in discussions on the prestellar inheritance of H2O based on comparisons between the observationally derived HDO/H2O ratio in clouds/cores and that in disks/comets. Alternatively, we propose that the ratio of D2O/HDO to HDO/H2O better probes the prestellar inheritance of H2O. It is also found that icy organics are more enriched in deuterium than water ice in forming disks. The differential deuterium fractionation in water and organics is inherited from the prestellar stages.

Molecular oxygen has been confirmed as the fourth most abundant molecule in cometary material O$_2$/H$_2$O $\sim 4$ %) and is thought to have a primordial nature, i.e., coming from the interstellar cloud from which our solar system was formed. However, interstellar O$_2$ gas is notoriously difficult to detect and has only been observed in one potential precursor of a solar-like system. Here, the chemical and physical origin of O$_2$ in comets is investigated using sophisticated astrochemical models. Three origins are considered: i) in dark clouds, ii) during forming protostellar disks, and iii) during luminosity outbursts in disks. The dark cloud models show that reproduction of the observed abundance of O$_2$ and related species in comet 67P/C-G requires a low H/O ratio facilitated by a high total density ($\geq 10^5$ cm$^{-3}$), and a moderate cosmic ray ionisation rate ($\leq 10^{-16}$ s$^{-1}$) while a temperature of 20 K, slightly higher than the typical temperatures found in dark clouds, also enhances the production of O$_2$. Disk models show that O$_2$ can only be formed in the gas phase in intermediate disk layers, and cannot explain the strong correlation between O$_2$ and H$_2$O in comet 67P/C-G together with the weak correlation between other volatiles and H$_2$O. However, primordial O$_2$ ice can survive transport into the comet-forming regions of disks. Taken together, these models favour a dark cloud (or "primordial") origin for O$_2$ in comets, albeit for dark clouds which are warmer and denser than those usually considered as solar system progenitors.

Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few % of water ice) at large disc radii (R > 30 AU). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 AU and both are nearly solar in the case of pure infall, but both are sub-solar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically-significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the Solar Nebula.

Recent interferometer observations have found that the D2O/HDO abundance ratio is higher than that of HDO/H2O by about one order of magnitude in the vicinity of low-mass protostar NGC 1333-IRAS 2A, where water ice has sublimated. Previous laboratory and theoretical studies show that the D2O/HDO ice ratio should be lower than the HDO/H2O ice ratio, if HDO and D2O ices are formed simultaneously with H2O ice. In this work, we propose that the observed feature, D2O/HDO > HDO/H2O, is a natural consequence of chemical evolution in the early cold stages of low-mass star formation: 1) majority of oxygen is locked up in water ice and other molecules in molecular clouds, where water deuteration is not efficient, and 2) water ice formation continues with much reduced efficiency in cold prestellar/protostellar cores, where deuteration processes are highly enhanced due to the drop of the ortho-para ratio of H2, the weaker UV radiation field, etc. Using a simple analytical model and gas-ice astrochemical simulations tracing the evolution from the formation of molecular clouds to protostellar cores, we show that the proposed scenario can quantitatively explain the observed HDO/H2O and D2O/HDO ratios. We also find that the majority of HDO and D2O ices are likely formed in cold prestellar/protostellar cores rather than in molecular clouds, where the majority of H2O ice is formed. This work demonstrates the power of the combination of the HDO/H2O and D2O/HDO ratios as a tool to reveal the past history of water ice formation in the early cold stages of star formation and when the enrichment of deuterium in the bulk of water occurred. Further observations are needed to explore if the relation, D2O/HDO > HDO/H2O, is common in low-mass protostellar sources.

We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.

In a protoplanetary disk, a combination of thermal and non-thermal desorption processes regulate where volatiles are liberated from icy grain mantles into the gas phase. Non-thermal desorption should result in volatile-enriched gas in disk-regions where complete freeze-out is otherwise expected. We present ALMA observations of the disk around the young star IM Lup in 1.4 mm continuum, C18O 2-1, H13CO+ 3-2 and DCO+ 3-2 emission at ~0".5 resolution. The images of these dust and gas tracers are clearly resolved. The DCO+ line exhibits a striking pair of concentric rings of emission that peak at radii of ~0".6 and 2" (~90 and 300 AU, respectively). Based on disk chemistry model comparison, the inner DCO+ ring is associated with the balance of CO freeze-out and thermal desorption due to a radial decrease in disk temperature. The outer DCO+ ring is explained by non-thermal desorption of CO ice in the low-column-density outer disk, repopulating the disk midplane with cold CO gas. The CO gas then reacts with abundant H2D+ to form the observed DCO+ outer ring. These observations demonstrate that spatially resolved DCO+ emission can be used to trace otherwise hidden cold gas reservoirs in the outmost disk regions, opening a new window onto their chemistry and kinematics.

We investigate the chemistry of ion molecules in protoplanetary disks, motivated by the detection of N$_2$H$^+$ ring around TW Hya. While the ring inner radius coincides with the CO snow line, it is not apparent why N$_2$H$^+$ is abundant outside the CO snow line in spite of the similar sublimation temperatures of CO and N$_2$. Using the full gas-grain network model, we reproduced the N$_2$H$^+$ ring in a disk model with millimeter grains. The chemical conversion of CO and N$_2$ to less volatile species (sink effect hereinafter) is found to affect the N$_2$H$^+$ distribution. Since the efficiency of the sink depends on various parameters such as activation barriers of grain surface reactions, which are not well constrained, we also constructed the no-sink model; the total (gas and ice) CO and N$_2$ abundances are set constant, and their gaseous abundances are given by the balance between adsorption and desorption. Abundances of molecular ions in the no-sink model are calculated by analytical formulas, which are derived by analyzing the full-network model. The N$_2$H$^+$ ring is reproduced by the no-sink model, as well. The 2D (R-Z) distribution of N$_2$H$^+$, however, is different among the full-network model and no-sink model. The column density of N$_2$H$^+$ in the no-sink model depends sensitively on the desorption rate of CO and N$_2$, and the flux of cosmic ray. We also found that N$_2$H$^+$ abundance can peak at the temperature slightly below the CO sublimation, even if the desorption energies of CO and N$_2$ are the same.

Observations of comets and asteroids show that the Solar Nebula that spawned our planetary system was rich in water and organic molecules. Bombardment brought these organics to the young Earth's surface, seeding its early chemistry. Unlike asteroids, comets preserve a nearly pristine record of the Solar Nebula composition. The presence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect to water, is of special interest because of the importance of C-N bonds for abiotic amino acid synthesis. Comet-like compositions of simple and complex volatiles are found in protostars, and can be readily explained by a combination of gas-phase chemistry to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances complexity[3]. Simple volatiles, including water and HCN, have been detected previously in Solar Nebula analogues - protoplanetary disks around young stars - indicating that they survive disk formation or are reformed in situ. It has been hitherto unclear whether the same holds for more complex organic molecules outside of the Solar Nebula, since recent observations show a dramatic change in the chemistry at the boundary between nascent envelopes and young disks due to accretion shocks[8]. Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary disk around the young star MWC 480. We find abundance ratios of these N-bearing organics in the gas-phase similar to comets, which suggests an even higher relative abundance of complex cyanides in the disk ice. This implies that complex organics accompany simpler volatiles in protoplanetary disks, and that the rich organic chemistry of the Solar Nebula was not unique.

We have observed the clump G34.43+00.24 MM3 associated with an infrared dark cloud in DNC $J$=3--2, HN$^{13}$C $J$=3--2, and N$_2$H$^+$ $J$=3--2 with the Atacama Large Millimeter/submillimeter Array (ALMA). The N$_2$H$^+$ emission is found to be relatively weak near the hot core and the outflows, and its distribution is clearly anti-correlated with the CS emission. This result indicates that a young outflow is interacting with cold ambient gas. The HN$^{13}$C emission is compact and mostly emanates from the hot core, whereas the DNC emission is extended around the hot core. Thus, the DNC and HN$^{13}$C emission traces warm regions near the protostar differently. The DNC emission is stronger than the HN$^{13}$C emission toward most parts of this clump. The DNC/HNC abundance ratio averaged within a $15^{\prime\prime} \times 15^{\prime\prime}$ area around the phase center is higher than 0.06. This ratio is much higher than the value obtained by the previous single-dish observations of DNC and HN$^{13}$C $J$=1--0 ($\sim$0.003). It seems likely that the DNC and HNC emission observed with the single-dish telescope traces lower density envelopes, while that observed with ALMA traces higher density and highly deuterated regions. We have compared the observational results with chemical-model results in order to investigate the behavior of DNC and HNC in the dense cores. Taking these results into account, we suggest that the low DNC/HNC ratio in the high-mass sources obtained by the single-dish observations are at least partly due to the low filling factor of the high density regions.

We have observed a molecular clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the CH3OH J_K=9_{-1}--8_0 E, 5_0-4_0 E, and 5_{-1}-4_{-1} E lines at sub-arcsecond resolution by using the Atacama Large Millimeter/submillimeter Array. By comparing the CH3OH J_K=9_{-1}-8_0 E emission with the CH3OH 5_0-4_0 E and 5_{-1}-4_{-1} E emission, we have found that the CH3OH J_K=9_{-1}-8_0 E emission is masing. We have clearly shown that the CH3OH J_K=9_{-1}-8_0 masers arise from the post shocked gas in the interacting regions between the outflows and ambient dense gas. Toward the strongest peak of the CH3OH maser emission, SiO J=6-5 emission is very weak. This indicates that the CH3OH maser emission traces relatively old shocks or weak shocks.

We study the influence of the turbulent transport on ice chemistry in protoplanetary disks, focusing on carbon and nitrogen bearing molecules. Chemical rate equations are solved with the diffusion term, mimicking the turbulent mixing in the vertical direction. Turbulence can bring ice-coated dust grains from the midplane to the warm irradiated disk surface, and the ice mantles are reprocessed by photoreactions, thermal desorption, and surface reactions. The upward transport decreases the abundance of methanol and ammonia ices at r < 30 AU, because warm dust temperature prohibits their reformation on grain surfaces. This reprocessing could explain the smaller abundances of carbon and nitrogen bearing molecules in cometary coma than those in low-mass protostellar envelopes. We also show the effect of mixing on the synthesis of complex organic molecules (COMs) are two ways: (1) transport of ices from the midplane to the disk surface and (2) transport of atomic hydrogen from the surface to the midplane. The former enhances the COMs formation in the disk surface, while the latter suppresses it in the midplane. Then, when mixing is strong, COMs are predominantly formed in the disk surface, while their parent molecules are (re)formed in the midplane. This cycle expands the COMs distribution both vertically and radially outward compared with that in the non-turbulent model. We derive the timescale of the sink mechanism by which CO and N2 are converted to less volatile molecules to be depleted from the gas phase, and find that the vertical mixing suppresses this mechanism in the inner disks.

Chemical models predict that the deuterated fraction (the column density ratio between a molecule containing D and its counterpart containing H) of N2H+, Dfrac(N2H+), is high in massive pre-protostellar cores and rapidly drops of an order of magnitude after the protostar birth, while that of HNC, Dfrac(HNC), remains constant for much longer. We tested these predictions by deriving Dfrac(HNC) in 22 high-mass star forming cores divided in three different evolutionary stages, from high-mass starless core candidates (HMSCs, 8) to high-mass protostellar objects (HMPOs, 7) to Ultracompact HII regions (UCHIIs, 7). For all of them, Dfrac (N2H+) was already determined through IRAM-30m Telescope observations, which confirmed the theoretical rapid decrease of Dfrac(N2H+) after protostar birth (Fontani et al. 2011). Therefore our comparative study is not affected by biases introduced by the source selection. We have found average Dfrac(HNC) of 0.012, 0.009 and 0.008 in HMSCs, HMPOs and UCHIIs, respectively, with no statistically significant differences among the three evolutionary groups. These findings confirm the predictions of the chemical models, and indicate that large values of Dfrac(N2H+) are more suitable than large values of Dfrac(HNC) to identify cores on the verge of forming high-mass stars, likewise what found in the low-mass regime.

We investigate water and deuterated water chemistry in turbulent protoplanetary disks. Chemical rate equations are solved with the diffusion term, mimicking turbulent mixing in vertical direction. Water near the midplane is transported to the disk atmosphere by turbulence and destroyed by photoreactions to produce atomic oxygen, while the atomic oxygen is transported to the midplane and reforms water and/or other molecules. We find that this cycle significantly decreases column densities of water ice at r < 30 AU, where dust temperatures are too high to reform water ice effectively. The radial extent of such region depends on the desorption energy of atomic hydrogen. Our model indicates that water ice could be deficient even outside the sublimation radius. Outside this radius, the cycle decreases the D/H ratio of water ice from 2x10^-2, which is set by the collapsing core model, to 10^-4-10^-2 in 10^6 yr, without significantly decreasing the water ice column density. The resultant D/H ratios depend on the strength of mixing and the radial distance from the central star. Our finding suggests that the D/H ratio of cometary water (10^-3-10^-4) could be established (i.e. cometary water could be formed) in the solar nebula, even if the D/H ratio of water ice delivered to the disk was very high (10^-2).

We have observed a cluster forming clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the 1.3 mm continuum and the CH3OH, CS, 13CS, SiO, CH3CH2CN, and HCOOCH3 lines with the Atacama Large Millimeter/submillimeter Array and in K-band with the Keck telescope. We have found a young outflow toward the center of this clump in the SiO, CS, and CH3OH lines. This outflow is likely driven by a protostar embedded in a hot core, which is traced by the CH3CH2CN, HCOOCH3, 13CS, and high excitation CH3OH lines. The size of the hot core is about 800 x 300 AU in spite of its low mass (<1.1 M_sun), suggesting a high accretion rate or the presence of multiple star system harboring a few hot corinos. The outflow is highly collimated, and the dynamical timescale is estimated to be less than 740 yr. In addition, we have also detected extended emission of SiO, CS, and CH3OH, which is not associated with the hot core and the outflow. This emission may be related to past star formation activity in the clump. Although G34.43+00.24 MM3 is surrounded by a dark feature in infrared, it has already experienced active formation of low-mass stars in an early stage of clump evolution.

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

We have observed the HN13C J=1-0 and DNC J=1-0 lines toward 18 massive clumps, including infrared dark clouds (IRDCs) and high-mass protostellar objects (HMPOs), by using the Nobeyama Radio Observatory 45 m telescope. We have found that the HN13C emission is stronger than the DNC emission toward all the observed sources. The averaged DNC/HNC ratio is indeed lower toward the observed high-mass sources (0.009\pm0.005) than toward the low-mass starless and star-forming cores (0.06). The kinetic temperature derived from the NH3 (J, K) = (1, 1) and (2, 2) line intensities is higher toward the observed high-mass sources than toward the low-mass cores. However the DNC/HNC ratio of some IRDCs involving the Spitzer 24 {\mu}m sources is found to be lower than that of HMPOs, although the kinetic temperature of the IRDCs is lower than that of the HMPOs. This implies that the DNC/HNC ratio does not depend only on the current kinetic temperature. With the aid of chemical model simulations, we discuss how the DNC/HNC ratio decreases after the birth of protostars. We suggest that the DNC/HNC ratio in star-forming cores depends on the physical conditions and history in their starless-core phase, such as its duration time and the gas kinetic temperature.

We construct the gas-grain chemical network model which includes carbon isotopes (12C and 13C) with an emphasis on isotopomer-exchange reactions. Temporal variations of molecular abundances, the carbon isotope ratios (12CX/13CX) and the isotopomer ratios (12C13CX/13C12CX) of CCH and CCS in cold dense cloud cores are investigated by numerical calculations. We confirm that the isotope ratios of molecules, both in the gas phase and grain surfaces, are significantly different depending on whether the molecule is formed from the carbon atom (ion) or the CO molecule. Molecules formed from carbon atoms have the CX/13CX ratios greater than the elemental abundance ratio of [12C/13C]. On the other hand, molecules formed from CO molecules have the CX/13CX ratios smaller than the [12C/13C] ratio. We reproduce the observed C13CH/13CCH ratio in TMC-1, if the isotopomer exchange reaction, 13CCH + H <-> C13CH + H + 8.1 K, proceeds with the forward rate coefficient kf > 10^-11 cm3 s-1. However, the C13CS/13CCS ratio is lower than that observed in TMC-1. We then assume the isotopomer exchange reaction catalyzed by the H atom, 13CCS + H <-> C13CS + H + 17.4 K. In the model with this reaction, we reproduce the observed C13CS/13CCS, CCS/C13CS and CCS/13CCS ratio simultaneously.