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.

We have observed the Class I protostar L1489 IRS with the Atacama Millimeter/submillimeter Array (ALMA) in Band 6. The C$^{18}$O $J=$2-1 line emission shows flattened and non-axisymmetric structures in the same direction as its velocity gradient due to rotation. We discovered that the C$^{18}$O emission shows dips at a radius of ~200-300 au while the 1.3 mm continuum emission extends smoothly up to r~400 au. At the radius of the C$^{18}$O dips, the rotational axis of the outer portion appears to be tilted by ~15 degrees from that of the inner component. Both the inner and outer components with respect to the C$^{18}$O dips exhibit the $r^{-0.5}$ Keplerian rotation profiles until r~600 au. These results not only indicate that a Keplerian disk extends up to ~600 au but also that the disk is warped. We constructed a three dimensional warped disk model rotating at the Keplerian velocity, and demonstrated that the warped disk model reproduces main observed features in the velocity channel maps and the PV diagrams. Such a warped disk system can form by mass accretion from a misaligned envelope. We also discuss a possible disk evolution scenario based on comparisons of disk radii and masses between Class I and Class II sources.

We present the results of 0.1-pc-scale observations in 250 GHz and 350GHz towards a newly-discovered hot molecular core in a nearby low-metallicity galaxy, the Large Magellanic Cloud (LMC), with the Atacama Large Millimeter/submillimeter Array. A variety of C/N/O/Si/S-bearing molecules are detected towards the high-mass young stellar object, ST16. A rotating protostellar envelope is for the first time detected outside our Galaxy by SO2 and 34SO lines. An outflow cavity is traced by CCH and CN. The isotope abundance of sulfur in the source is estimated to be 32S/34S = 17 and 32S/33S = 53 based on SO, SO2, and CS isotopologues, suggesting that both 34S and 33S are overabundant in the LMC. Rotation diagram analyses show that the source is associated with hot gas (>100K) traced by high-excitation lines of CH3OH and SO2, as well as warm gas (~50K) traced by CH3OH, SO2, 34SO, OCS, CH3CN lines. A comparison of molecular abundances between LMC and Galactic hot cores suggests that organic molecules (e.g., CH3OH, a classical hot core tracer) show a large abundance variation in low metallicity, where the present source is classified into an organic-poor hot core. Our astrochemical simulations suggest that different grain temperature during the initial ice-forming stage would contribute to the chemical differentiation. In contrast, SO2 shows similar abundances within all the known LMC hot cores and the typical abundance roughly scales with the LMC's metallicity. Nitrogen-bearing molecules are generally less abundant in LMC hot cores, except for NO. The present results suggest that chemical compositions of hot cores do not always simply scale with the metallicity.

EC53 is an embedded protostar with quasi-periodic emission in the near-IR and sub-mm. We use ALMA high-resolution observations of continuum and molecular line emission to describe the circumstellar environment of EC 53. The continuum image reveals a disk with a flux that suggests a mass of 0.075 Msun, much less than the estimated mass in the envelope, and an in-band spectral index that indicates grain growth to centimeter sizes. Molecular lines trace the outflow cavity walls, infalling and rotating envelope, and/or the Keplerian disk. The rotation profile of the C17O 3-2 line emission cannot isolate the Keplerian motion clearly although the lower limit of the protostellar mass can be calculated as 0.3 +- 0.1 Msun if the Keplerian motion is adopted. The weak CH3OH emission, which is anti-correlated with the HCO+ 4-3 line emission, indicates that the water snow line is more extended than what expected from the current luminosity, attesting to bygone outburst events. The extended snow line may persist for longer at the disk surface because the lower density increases the freeze-out timescale of methanol and water.

We have observed the submillimeter continuum condensations SMM2, SMM4, SMM9, and SMM11 in the star forming cluster Serpens Main using the Atacama Large Millimeter/submillimeter Array during Cycle 3 in the 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at an angular resolution of ~0.55 (240 au). Sixteen sources have been detected in the 1.3 mm continuum, which can be classified into three groups. Group 1 consists of six sources showing extended continuum emission and bipolar/monopolar 12CO outflows. Although all the Group 1 members are classified as Class 0 protostars, our observations suggest evolutionary trends among them in terms of 12CO outflow dynamical time, SO emission distribution, C18O fractional abundance, and continuum morphology. Group 2 consists of four sources associated with a continuum filamentary structure and no 12CO outflows. Central densities estimated from the 1.3 mm continuum intensity suggest that they are prestellar sources in a marginally Jeans unstable state. Group 3 consists of six Spitzer sources showing point-like 1.3 mm continuum emission and clumpy 12CO outflows. These features of Group 3 suggest envelope dissipation, preventing disk growth from the present size, r <~ 60 au. The Group 3 members are protostars that may be precursors to the T Tauri stars associated with small disks at tens-au radii identified in recent surveys.

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.

We have observed the Class I protostellar source Elias 29 with Atacama Large Millimeter/submillimeter Array (ALMA). We have detected CS, SO, $^{34}$SO, SO$_2$, and SiO line emissions in a compact component concentrated near the protostar and a ridge component separated from the protostar by 4\arcsec\ ($\sim 500$ au). The former component is found to be abundant in SO and SO$_2$ but deficient in CS. The abundance ratio SO/CS is as high as $3^{+13}_{-2} \times 10^2$ at the protostar, which is even higher than that in the outflow-shocked region of L1157 B1. However, organic molecules (HCOOCH$_3$, CH$_3$OCH$_3$, CCH, and c-C$_3$H$_2$) are deficient in Elias 29. We attribute the deficiency in organic molecules and richness in SO and SO$_2$ to the evolved nature of the source or the relatively high dust temperature (\protect\raisebox{-0.7ex}{$\:\stackrel{\textstyle >}{\sim}\:$} 20 K) in the parent cloud of Elias 29. The SO and SO$_2$ emissions trace rotation around the protostar. Assuming a highly inclined configuration ($i \geq 65$\degr; 0\degr\ for a face-on configuration) and Keplerian motion for simplicity, the protostellar mass is estimated to be (0.8 -- 1.0) \Msun. The $^{34}$SO and SO$_2$ emissions are asymmetric in their spectra; the blue-shifted components are weaker than the red-shifted ones. Although this may be attributed to the asymmetric molecular distribution, other possibilities are also discussed.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

Complex organic molecules (COMs), which are the seeds of prebiotic material and precursors of amino acids and sugars, form in the icy mantles of circumstellar dust grains but cannot be detected remotely unless they are heated and released to the gas phase. Around solar-mass stars, water and COMs only sublimate in the inner few au of the disk, making them extremely difficult to spatially resolve and study. Sudden increases in the luminosity of the central star will quickly expand the sublimation front (so-called snow line) to larger radii, as seen previously in the FU Ori outburst of the young star V883 Ori. In this paper, we take advantage of the rapid increase in disk temperature of V883 Ori to detect and analyze five different COMs, methanol, acetone, acetonitrile, acetaldehyde, and methyl formate, in spatially-resolved submillimeter observations. The COMs abundances in V883 Ori is in reasonable agreement with cometary values. This result suggests that outbursting young stars can provide a special opportunity to study the ice composition of material directly related to planet formation.

We investigate the chemical evolution of a collapsing core that starts from a hydrostatic core and finally form a low-mass protostar. New multiphase gas-grain models that include bulk diffusion and photon penetration are simulated by the macroscopic Monte Carlo method in order to derive the chemical evolution. There are two types of species in the ice bulk in the new multiphase models. Interstitial species can diffuse and sublime at their own sublimation temperatures while normal species are locked in the ice bulk. Photodissociation rates of icy species are reduced by the exponential decay of UV flux within the ice mantle. Two-phase models and basic multiphase models without bulk diffusion and photon penetration are also simulated for comparison. Our physical model for the collapsing core is base on a one-dimensional radiation hydrodynamics model. Abundant icy radicals are produced at around 10 K in the new multiphase models. Interstitial radicals can diffuse inside ice mantles to form complex organic molecules (COMS) upon warming-up. Thus, COMs produced by radical recombination at higher temperatures in the new multiphase models are more than one order of magnitude higher than those in the two-phase and basic multiphase models. Moreover, COMs produced at around 10 K in the new multiphase models are about one order of magnitude higher than those in the two-phase model. Our model shows a reasonable agreement with observations toward low-mass protostars. Moreover, molecular oxygen abundances predicted by our new multiphase models agree reasonably well with that found in cometary materials.

It is important to understand the origin of molecular line intensities and chemical composition in the molecular-cloud scale in the Galactic sources because it serves as a benchmark to compare with the chemical compositions of extragalactic sources. Recent observations of the 3-mm spectra averaged over the 10-pc scale show similar spectral pattern among sources for molecular lines HCN, HCO$^+$, CCH, HNC, HNCO, c-C$_3$H$_2$, CS, SO, N$_2$H$^+$, and CN. To constrain the average physical property emitting such spectral pattern, we model molecular spectra using a time-dependent gas-grain chemical model followed by a radiative transfer calculation. We use a grid of physical parameters such as the density $n=3 \times 10^2 - 3\times 10^4$ cm$^{-3}$, the temperature, $T=10-30$ K, the visual extinction $A_{\rm V} = 2,4,10$ mag, the cosmic-ray ionization rate $\zeta = 10^{-17} - 10^{-16}$ s$^{-1}$, and the sulfur elemental abundance $S/H = 8\times 10^{-8} - 8\times 10^{-7}$. Comparison with the observed spectra indicates that spectra are well reproduced with the relatively low density of $n=(1-3) \times 10^3\,$cm$^{-3}$, $T=10\,$K, $\zeta = 10^{-17}$ s$^{-1}$, and the short chemistry timescale of $10^5$ yrs. This short chemistry timescale may indicate that molecular clouds are constantly affected by the turbulence, and exposed to low-density, low $A_{\rm V}$ regions that "refreshes" the chemical clock by UV radiation. The relatively low density obtained is orders of magnitude lower than the commonly-quoted critical density in the optically thin case. Meanwhile, this range of density is consistent with results from recent observational analysis of molecular-cloud-scale mapping.

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.

We report the results of subparsec-scale submillimeter observations towards an embedded high-mass young stellar object in the Small Magellanic Cloud (SMC) with ALMA. Complementary infrared data obtained with the AKARI satellite and the Gemini South telescope are also presented. The target infrared point source is spatially resolved into two dense molecular cloud cores; one is associated with a high-mass young stellar object (YSO core), while another is not associated with an infrared source (East core). The two cores are dynamically associated but show different chemical characteristics. Emission lines of CS, C33S, H2CS, SO, SO2, CH3OH, H13CO+, H13CN, SiO, and dust continuum are detected from the observed region. Tentative detection of HDS is also reported. The first detection of CH3OH in the SMC has a strong impact on our understanding of the formation of complex organic molecules in metal-poor environments. The gas temperature is estimated to be ~10 K based on the rotation analysis of CH3OH lines. The fractional abundance of CH3OH gas in the East core is estimated to be (0.5-1.5) x 10^(-8), which is comparable with or marginally higher than those of similar cold sources in our Galaxy despite a factor of five lower metallicity in the SMC. This work provides observational evidence that an organic molecule like CH3OH, which is largely formed on grain surfaces, can be produced even in a significantly lower metallicity environment compared to the solar neighborhood. A possible origin of cold CH3OH gas in the observed dense core is discussed.

We have observed the submillimeter continuum condensation SMM4 in Serpens Main using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 3 in 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at angular resolutions of ~0.55" (240 au). The 1.3 mm continuum emission shows that SMM4 is spatially resolved into two protostars embedded in the same core: SMM4A showing a high brightness temperature, 18 K, with little extended structure and SMM4B showing a low brightness temperature, 2 K, with compact and extended structures. Their separation is ~2100 au. Analysis of the continuum visibilities reveals a disk-like structure with a sharp edge at r~240 au in SMM4A, and a compact component with a radius of 56 au in SMM4B. The 12CO emission traces fan-shaped and collimated outflows associated with SMM4A and SMM4B, respectively. The blue and red lobes of the SMM4B outflow have different position angles by ~30 deg. Their inclination and bending angles in the 3D space are estimated at i_b~ 36 deg, i_r~70 deg, and alpha~40 deg, respectively. The SO emission traces shocked regions, such as cavity walls of outflows and the vicinity of SMM4B. The C18O emission mainly traces an infalling and rotating envelope around SMM4B. The C18O fractional abundance in SMM4B is ~50 times smaller than that of the interstellar medium. These results suggest that SMM4A is more evolved than SMM4B. Our studies in Serpens Main demonstrate that continuum and line observations at millimeter wavelengths allow us to differentiate evolutionary phases of protostars within the Class 0 phase.

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.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.

We report the discovery of an extremely young protostar, SMM11, located in the associated submillimeter condensation SMM11 in the Serpens Main cluster using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 3 at 1.3 mm and an angular resolution of ~0.5"~210 AU. SMM11 is a Class 0 protostar without any counterpart at 70 um or shorter wavelengths. The ALMA observations show 1.3 mm continuum emission associated with a collimated 12CO bipolar outflow. Spitzer and Herschel data show that SMM11 is extremely cold (T_bol=26 K) and faint (L_bol<~0.9 Lsun). We estimate the inclination angle of the outflow to be ~80 deg, almost parallel to the plane of the sky, from simple fitting using wind-driven-shell model. The continuum visibilities consist of Gaussian and power-law components, suggesting a spherical envelope with a radius of ~600 AU around the protostar. The estimated low C18O abundance, X(C18O)=1.5-3x10^-10, is also consistent with its youth. The high outflow velocity, a few 10 kms^-1 at a few 1000 AU, is much higher than theoretical simulations of first hydrostatic cores and we suggest that SMM11 is a transitional object right after the second collapse of the first core.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

In order to study a molecular-cloud-scale chemical composition, we have conducted a mapping spectral line survey toward the Galactic molecular cloud W3(OH), which is one of the most active star forming regions in the Perseus arm, with the NRO 45 m telescope. We have observed the area of 16' $\times$ 16', which corresponds to 9.0 pc $\times$ 9.0 pc. The observed frequency ranges are 87--91, 96--103, and 108--112 GHz. We have prepared the spectrum averaged over the observed area, in which 8 molecular species CCH, HCN, HCO$^+$, HNC, CS, SO, C$^{18}$O, and $^{13}$CO are identified. On the other hand, the spectrum of the W3(OH) hot core observed at a 0.17 pc resolution shows the lines of various molecules such as OCS, H$_2$CS CH$_3$CCH, and CH$_3$CN, in addition to the above species. In the spatially averaged spectrum, emission of the species concentrated just around the star-forming core such as CH$_3$OH and HC$_3$N is fainter than in the hot core spectrum, whereas emission of the species widely extended over the cloud such as CCH is relatively brighter. We have classified the observed area into 5 subregions according to the integrated intensity of $^{13}$CO, and have evaluated the contribution to the averaged spectrum from each subregion. The CCH, HCN, HCO$^+$, and CS lines can be seen even in the spectrum of the subregion with the lowest $^{13}$CO integrated intensity range ($< 10$ K km s$^{-1}$). Thus, the contributions of the spatially extended emission is confirmed to be dominant in the spatially averaged spectrum.

We have newly observed the Class 0/I protostar L1527 IRS using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 1 in 220 GHz dust continuum and C18O (J=2-1) line emissions with a ~2 times higher angular resolution (~0.5") and ~4 times better sensitivity than our ALMA Cycle 0 observations. Continuum emission shows elongation perpendicular to the associated outflow, with a deconvolved size of 0.53"x0.15". C18O emission shows similar elongation, indicating that both emissions trace the disk and the flattened envelope surrounding the protostar. The velocity gradient of the C18O emission along the elongation due to rotation of the disk/envelope system is re-analyzed, identifying Keplerian rotation proportional to r^-0.5 more clearly than the Cycle 0 observations. The Keplerian-disk radius and the dynamical stellar mass are kinematically estimated to be ~74 AU and ~0.45 Mo, respectively. The continuum visibility is fitted by models without any annulus averaging, revealing that the disk is in hydrostatic equilibrium. The best-fit model also suggests a density jump by a factor of ~5 between the disk and the envelope, suggesting that disks around protostars can be geometrically distinguishable from the envelope from a viewpoint of density contrast. Importantly, the disk radius geometrically identified with the density jump is consistent with the radius kinematically estimated. Possible origin of the density jump due to the mass accretion from the envelope to the disk is discussed. C18O observations can be reproduced by the same geometrical structures derived from the dust observations, with possible C18O freeze-out and localized C18O desorption.

We used the Atacama Large (sub-)Millimeter Array (ALMA) and the IRAM Plateau de Bure Interferometer (PdBI) to image, with an angular resolution of 0.5$''$ (120 au) and 1$''$ (235 au), respectively, the emission from 11 different organic molecules in the protostellar binary NGC1333 IRAS 4A. We clearly disentangled A1 and A2, the two protostellar cores present. For the first time, we were able to derive the column densities and fractional abundances simultaneously for the two objects, allowing us to analyse the chemical differences between them. Molecular emission from organic molecules is concentrated exclusively in A2 even though A1 is the strongest continuum emitter. The protostellar core A2 displays typical hot corino abundances and its deconvolved size is 70 au. In contrast, the upper limits we placed on molecular abundances for A1 are extremely low, lying about one order of magnitude below prestellar values. The difference in the amount of organic molecules present in A1 and A2 ranges between one and two orders of magnitude. Our results suggest that the optical depth of dust emission at these wavelengths is unlikely to be sufficiently high to completely hide a hot corino in A1 similar in size to that in A2. Thus, the significant contrast in molecular richness found between the two sources is most probably real. We estimate that the size of a hypothetical hot corino in A1 should be less than 12 au. Our results favour a scenario in which the protostar in A2 is either more massive and/or subject to a higher accretion rate than A1, as a result of inhomogeneous fragmentation of the parental molecular clump. This naturally explains the smaller current envelope mass in A2 with respect to A1 along with its molecular richness.

We have conducted a mapping spectral line survey toward the Galactic giant molecular cloud W51 in the 3 mm band with the Mopra 22 m telescope in order to study an averaged chemical composition of the gas extended over a molecular cloud scale in our Galaxy. We have observed the area of $25' \times 30'$, which corresponds to 39 pc $\times$ 47 pc. The frequency ranges of the observation are 85.1 - 101.1 GHz and 107.0 - 114.9 GHz. In the spectrum spatially averaged over the observed area, spectral lines of 12 molecular species and 4 additional isotopologues are identified. An intensity pattern of the spatially-averaged spectrum is found to be similar to that of the spiral arm in the external galaxy M51, indicating that these two sources have similar chemical compositions. The observed area has been classified into 5 sub-regions according to the integrated intensity of $^{13}$CO($J=1-0$) ($I_{\rm ^{13}CO}$), and contributions of the fluxes of 11 molecular lines from each sub-region to the averaged spectrum have been evaluated. For most of molecular species, 50 % or more of the flux come from the sub-regions with $I_{\rm ^{13}CO}$ from 25 K km s$^{-1}$ to 100 K km s$^{-1}$, which does not involve active star forming regions. Therefore, the molecular-cloud-scale spectrum observed in the 3 mm band hardly represents the chemical composition of star forming cores, but mainly represents the chemical composition of an extended quiescent molecular gas. The present result constitutes a sound base for interpreting the spectra of external galaxies at a resolution of a molecular cloud scale ($\sim10$ pc) or larger.

Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations.

We have resolved for the first time the radial and vertical structure of the almost edge-on envelope/disk system of the low-mass Class 0 protostar L1527. For that, we have used ALMA observations with a spatial resolution of 0.25$^{\prime\prime}$$\times$0.13$^{\prime\prime}$ and 0.37$^{\prime\prime}$$\times$0.23$^{\prime\prime}$ at 0.8 mm and 1.2 mm, respectively. The L1527 dust continuum emission has a deconvolved size of 78 au $\times$ 21 au, and shows a flared disk-like structure. A thin infalling-rotating envelope is seen in the CCH emission outward of about 150 au, and its thickness is increased by a factor of 2 inward of it. This radius lies between the centrifugal radius (200 au) and the centrifugal barrier of the infalling-rotating envelope (100 au). The gas stagnates in front of the centrifugal barrier and moves toward vertical directions. SO emission is concentrated around and inside the centrifugal barrier. The rotation speed of the SO emitting gas is found to be decelerated around the centrifugal barrier. A part of the angular momentum could be extracted by the gas which moves away from the mid-plane around the centrifugal barrier. If this is the case, the centrifugal barrier would be related to the launching mechanism of low velocity outflows, such as disk winds.

The Class 0 protostar, L483, has been observed in various molecular lines in the 1.2 mm band at a sub-arcsecond resolution with ALMA. An infalling-rotating envelope is traced by the CS line, while a very compact component with a broad velocity width is observed for the CS, SO, HNCO, NH$_2$CHO, and HCOOCH$_3$ lines. Although this source is regarded as the warm carbon-chain chemistry (WCCC) candidate source at a 1000 au scale, complex organic molecules characteristic of hot corinos such as NH$_2$CHO and HCOOCH$_3$ are detected in the vicinity of the protostar. Thus, both hot corino chemistry and WCCC are seen in L483. Although such a mixed chemical character source has been recognized as an intermediate source in previous single-dish observations, we here report the first spatially-resolved detection. A kinematic structure of the infalling-rotating envelope is roughly explained by a simple ballistic model with the protostellar mass of 0.1--0.2 $M_\odot$ and the radius of the centrifugal barrier (a half of the centrifugal radius) of 30--200 au, assuming the inclination angle of 80\degr\ (0\degr\ for a face-on). The broad line emission observed in the above molecules most likely comes from the disk component inside the centrifugal barrier. Thus, a drastic chemical change is seen around the centrifugal barrier.

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.

Transitional disks show a lack of excess emission at infrared wavelengths due to a large dust cavity, that is often corroborated by spatially resolved observations at ~ mm wavelengths. We present the first spatially resolved ~ mm-wavelength images of the disk around the Herbig Ae/Be star, HD 97048. Scattered light images show that the disk extends to ~640 au. The ALMA data reveal a circular-symmetric dusty disk extending to ~350 au, and a molecular disk traced in CO J=3-2 emission, extending to ~750 au. The CO emission arises from a flared layer with an opening angle ~ 30 deg - 40 deg. HD 97048 is another source for which the large (~ mm-sized) dust grains are more centrally concentrated than the small (~ {\mu}m-sized) grains and molecular gas, likely due to radial drift. The images and visibility data modelling suggests a decrement in continuum emission within ~50 au, consistent with the cavity size determined from mid-infrared imaging (34 +/- 4 au). The extracted continuum intensity profiles show ring-like structures with peaks at ~50, 150, and 300 au, with associated gaps at ~100 and 250 au. This structure should be confirmed in higher-resolution images (FWHM ~ 10 - 20 au). These data confirm the classification of HD 97048 as a transitional disk that also possesses multiple ring-like structures in the dust continuum emission. Additional data are required at multiple and well-separated frequencies to fully characterise the disk structure, and thereby constrain the mechanism(s) responsible for sculpting the HD 97048 disk.

We report the first evidence of a hot corino in a Bok globule. This is based on the ALMA observations in the 1.2 mm band toward the low-mass Class 0 protostar IRAS 19347+0727 in B335. Saturated complex organic molecules (COMs), CH$_3$CHO, HCOOCH$_3$, and NH$_2$CHO, are detected in a compact region within a few 10 au around the protostar. Additionally, CH$_3$OCH$_3$, C$_2$H$_5$OH, C$_2$H$_5$CN, and CH$_3$COCH$_3$ are tentatively detected. Carbon-chain related molecules, CCH and c-C$_3$H$_2$, are also found in this source, whose distributions are extended over a few 100 au scale. On the other hand, sulfur-bearing molecules CS, SO, and SO$_2$, have both compact and extended components. Fractional abundances of the COMs relative to H$_2$ are found to be comparable to those in known hot-corino sources. Though the COMs lines are as broad as 5-8 km s$^{-1}$, they do not show obvious rotation motion in the present observation. Thus, the COMs mainly exist in a structure whose distribution is much smaller than the synthesized beam (0."58 x 0."52).

We have conducted a spectral line survey observation in the 3 mm band toward the low-metallicity dwarf galaxy IC10 with the 45 m radio telescope of Nobeyama Radio Observatory to explore its chemical composition at a molecular-cloud scale (~80 pc). The CS, SO, CCH, HCN, HCO+, and HNC lines are detected for the first time in this galaxy in addition to the CO and 13CO lines, while c-C3H2, CH3OH, CN, C18O, and N2H+ lines are not detected. The spectral intensity pattern is found to be similar to those observed toward molecular clouds in the Large Magellanic Cloud, whose metallicity is as low as IC10. Nitrogen-bearing species are deficient in comparison with the Galactic molecular clouds due to a lower elemental abundance of nitrogen. CCH is abundant in comparison with Galactic translucent clouds, whereas CH3OH may be deficient. These characteristic trends for CCH and CH3OH are also seen in the LMC, and seem to originate from photodissociation regions more extended in peripheries of molecular clouds due to the lower metallicity condition.

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 \sigma$), reaching a peak signal-to-noise of $5.5\sigma$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation.

We report the first detection of a hot molecular core outside our Galaxy based on radio observations with ALMA toward a high-mass young stellar object (YSO) in a nearby low metallicity galaxy, the Large Magellanic Cloud (LMC). Molecular emission lines of CO, C17O, HCO+, H13CO+, H2CO, NO, SiO, H2CS, 33SO, 32SO2, 34SO2, and 33SO2 are detected from a compact region (0.1 pc) associated with a high-mass YSO, ST11. The temperature of molecular gas is estimated to be higher than 100 K based on rotation diagram analysis of SO2 and 34SO2 lines. The compact source size, warm gas temperature, high density, and rich molecular lines around a high-mass protostar suggest that ST11 is associated with a hot molecular core. We find that the molecular abundances of the LMC hot core are significantly different from those of Galactic hot cores. The abundances of CH3OH, H2CO, and HNCO are remarkably lower compared with Galactic hot cores by at least 1-3 orders of magnitude. We suggest that these abundances are characterized by the deficiency of molecules whose formation requires the hydrogenation of CO on grain surfaces. In contrast, NO shows a high abundance in ST11 despite the notably low abundance of nitrogen in the LMC. A multitude of SO2 and its isotopologue line detections in ST11 imply that SO2 can be a key molecular tracer of hot core chemistry in metal-poor environments. Furthermore, we find molecular outflows around the hot core, which is the second detection of an extragalactic protostellar outflow. In this paper, we discuss physical and chemical characteristics of a hot molecular core in the low metallicity environment.

Sub-arcsecond images of the rotational line emission of CS and SO have been obtained toward the Class I protostar IRAS 04365$+$2535 in TMC-1A with ALMA. A compact component around the protostar is clearly detected in the CS and SO emission. The velocity structure of the compact component of CS reveals infalling-rotating motion conserving the angular momentum. It is well explained by a ballistic model of an infalling-rotating envelope with the radius of the centrifugal barrier (a half of the centrifugal radius) of 50 AU, although the distribution of the infalling gas is asymmetric around the protostar. The distribution of SO is mostly concentrated around the radius of the centrifugal barrier of the simple model. Thus a drastic change in chemical composition of the gas infalling onto the protostar is found to occur at a 50 AU scale probably due to accretion shocks, demonstrating that the infalling material is significantly processed before being delivered into the disk.

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.

Spectral line survey observations of 7 molecular clouds in the Large Magellanic Cloud (LMC) have been conducted in the 3 mm band with the Mopra 22 m telescope to reveal chemical compositions in low metallicity conditions. Spectral lines of fundamental species such as CS, SO, CCH, HCN, HCO+, and HNC are detected in addition to those of CO and 13CO, while CH3OH is not detected in any source and N2H+ is marginally detected in two sources. The molecular-cloud scale (10 pc scale) chemical composition is found to be similar among the 7 sources regardless of different star formation activities, and hence, it represents the chemical composition characteristic to the LMC without influences of star formation activities. In comparison with chemical compositions of Galactic sources, the characteristic features are (1) deficient N-bearing molecules, (2) abundant CCH, and (3) deficient CH3OH. The feature (1) is due to a lower elemental abundance of nitrogen in the LMC, whereas the features (2) and (3) seem to originate from extended photodissociation regions and warmer temperature in cloud peripheries due to a lower abundance of dust grains in the low metallicity condition. In spite of general resemblance of chemical abundances among the seven sources, the CS/HCO+ and SO/HCO+ ratios are found to be slightly higher in a quiescent molecular cloud. An origin of this trend is discussed in relation to possible depletion of sulfur along molecular cloud formation.

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.

We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form in the building blocks of comets, icy bodies, and giant planets. Methane clathrates may play an important role in the evolution of planetary atmospheres. On Earth, the production of methane in clathrates is essentially biological, and these compounds are mostly found in permafrost regions or in the sediments of continental shelves. On Mars, methane would more likely derive from hydrothermal reactions with olivine-rich material. If they do exist, martian methane clathrates would be stable only at depth in the cryosphere and sporadically release some methane into the atmosphere via mechanisms that remain to be determined.

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 have observed the Class I protostar TMC-1A with Atacama Millimeter/submillimeter Array (ALMA) in 12CO and C18O (J=2-1), and 1.3-mm dust continuum emission. Continuum emission with a deconvolved size of 0.50"x0.37", perpendicular to the 12CO outflow, is detected. It most likely traces a circumstellar disk around TMC-1A, as previously reported. In contrast, a more extended structure is detected in C18O although it is still elongated with a deconvolved size of 3.3"x2.2", indicating that C18O traces mainly a flattened envelope surrounding the disk and the central protostar. C18O shows a clear velocity gradient perpendicular to the outflow at higher velocities, indicative of rotation, while an additional velocity gradient along the outflow is found at lower velocities. The radial profile of the rotational velocity is analyzed in detail, finding that it is given as a power-law \propto r^{-a} with an index of ~0.5 at higher velocities. This suggests that the rotation at higher velocities can be explained as Keplerian rotation orbiting a protostar with a dynamical mass of 0.68 Mo (inclination-corrected). The additional velocity gradient of C18O along the outflow is considered to be mainly infall motions in the envelope. Position-Velocity diagrams made from models consisting of an infalling envelope and a Keplerian disk are compared with the observations, revealing that the observed infall velocity is ~0.3 times smaller than the free fall velocity yielded by the dynamical mass of the protostar. Magnetic fields could be responsible for the slow infall velocity. A possible scenario of Keplerian disk formation is discussed.

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.

In this Letter we model the chemistry of DCO$^{+}$ in protoplanetary disks. We find that the overall distribution of the DCO$^{+}$ abundance is qualitatively similar to that of CO but is dominated by thin layer located at the inner disk surface. To understand its distribution, we investigate the different key gas-phase deuteration pathways that can lead to the formation of DCO$^{+}$. Our analysis shows that the recent update in the exothermicity of the reaction involving CH$_2$D$^{+}$ as a parent molecule of DCO$^{+}$ favors deuterium fractionation in warmer conditions. As a result the formation of DCO$^{+}$ is enhanced in the inner warm surface layers of the disk where X-ray ionization occurs. Our analysis points out that DCO$^{+}$ is not a reliable tracer of the CO snow line as previously suggested. We thus predict that DCO$^{+}$ is a tracer of active deuterium and in particular X-ray ionization of the inner disk.

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.

Recent ALMA (Atacama Large Millimeter/submillimeter Array) observations of young protostellar objects detected warm SO emission, which could be associated with a forming protostellar disk. In order to investigate if such warm gas can be produced by accretion shock onto the forming disk, we calculate the sputtering and thermal desorption of various grain surface species in one dimensional shock waves. We find that thermal desorption is much more efficient than the sputtering in the post-shock region. While H$_{2}$O can be thermally desorbed, if the accretion velocity is larger than 8 km s$^{-1}$ with the pre-shock gas number density of 10$^{9}$ cm$^{-3}$, SO is desorbed, if the accretion velocity $\gtrsim$ 2 km s$^{-1}$ and $\gtrsim$ 4km s$^{-1}$, with the pre-shock density of 10$^{9}$ cm$^{-3}$ and 10$^{8}$ cm$^{-3}$, respectively. We also find that the column density of hydrogen nuclei in warm post-shock gas is $N_{{\rm warm}} \sim 10^{21}$ cm$^{-2}$.

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.

IRAS 16293-2422 is a well studied low-mass protostar characterized by a strong level of deuterium fractionation. In the line of sight of the protostellar envelope, an additional absorption layer, rich in singly and doubly deuterated water has been discovered by a detailed multiline analysis of HDO. To model the chemistry in this source, the gas-grain chemical code Nautilus has been used with an extended deuterium network. For the protostellar envelope, we solve the chemical reaction network in infalling fluid parcels in a protostellar core model. For the foreground cloud, we explored several physical conditions (density, cosmic ionization rate, C/O ratio). The main results of the paper are that gas-phase abundances of H2O, HDO and D2O observed in the inner regions of IRAS16293-2422 are lower than those predicted by a 1D dynamical/chemical (hot corino) model in which the ices are fully evaporated. The abundance in the outer part of the envelope present chaotic profiles due to adsorption/evaporation competition, very different from the constant abundance assumed for the analysis of the observations. We also found that the large abundances of gas-phase H2O, HDO and D2O observed in the absorption layer are more likely explained by exothermic surface reactions rather than photodesorption processes.

We report Atacama Large Millimeter/submillimeter Array (ALMA) cycle 0 observations of C$^{18}$O ($J=2-1$), SO ($J_N= 6_5-5_4$) and 1.3mm dust continuum toward L1527 IRS, a class 0 solar-type protostar surrounded by an infalling and rotating envelope. C$^{18}$O emission shows strong redshifted absorption against the bright continuum emission associated with L1527 IRS, strongly suggesting infall motions in the C$^{18}$O envelope. The C$^{18}$O envelope also rotates with a velocity mostly proportional to $r^{-1}$, where $r$ is the radius, while the rotation profile at the innermost radius (54 AU) may be shallower than $r^{-1}$, suggestive of formation of a Keplerian disk around the central protostar of 0.3 Mo in dynamical mass. SO emission arising from the inner part of the C$^{18}$O envelope also shows rotation in the same direction as the C$^{18}$O envelope. The rotation is, however, rigid-body like which is very different from the differential rotation shown by C$^{18}$O. In order to explain the line profiles and the position-velocity (PV) diagrams of C$^{18}$O and SO observed, simple models composed of an infalling envelope surrounding a Keplerian disk of 54 AU in radius orbiting a star of 0.3 Mo are examined. It is found that in order to reproduce characteristic features of the observed line profiles and PV diagrams, the infall velocity in the model has to be smaller than the free-fall velocity yielded by a star of 0.3 Mo. Possible reasons for the reduced infall velocities are discussed.

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

We have conducted ALMA observations in the 1.3 mm continuum and 12CO (2-1), C18O (2-1) and SO (5_6-4_5) lines toward L1489 IRS, a Class I protostar surrounded by a Keplerian disk and an infalling envelope. The Keplerian disk is clearly identified in the 12CO and C18O emission, and its outer radius (~700 AU) and mass (~0.005 Msun) are comparable to those of disks around T Tauri stars. The protostellar mass is estimated to be 1.6 Msun with the inclination angle of 66 deg. In addition to the Keplerian disk, there are blueshifted and redshifted off-axis protrusions seen in the C18O emission pointing toward the north and the south, respectively, adjunct to the middle part of the Keplerian disk. The shape and kinematics of these protrusions can be interpreted as streams of infalling flows with a conserved angular momentum following parabolic trajectories toward the Keplerian disk, and the mass infalling rate is estimated to be ~5E-7 Msun/yr. The specific angular momentum of the infalling flows (~2.5E-3 km/s*pc) is comparable to that at the outer radius of the Keplerian disk (~4.8E-3 km/s*pc). The SO emission is elongated along the disk major axis and exhibits a linear velocity gradient along the axis, which is interpreted as that the SO emission primarily traces a ring region in the flared Keplerian disk at radii of ~250-390 AU. The local enhancement of the SO abundance in the ring region can be due to the accretion shocks at the centrifugal radius where the infalling flows fall onto the disk. Our ALMA observations unveiled both the Keplerian disk and the infalling gas onto the disk, and the disk can further grow by accreting material and angular momenta from the infalling gas.

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.

While recent studies of the solar-mass protostar IRAS16293-2422 have focused on its inner arcsecond, the wealth of Herschel/HIFI data has shown that the structure of the outer envelope and of the transition region to the more diffuse ISM is not clearly constrained. We use rotational ground-state transitions of CH (methylidyne), as a tracer of the lower-density envelope. Assuming LTE, we perform a $\chi^2$ minimization of the high spectral resolution HIFI observations of the CH transitions at ~532 and ~536 GHz in order to derive column densities in the envelope and in the foreground cloud. We obtain column densities of (7.7$\pm$0.2)$\times10^{13}$ cm$^{-2}$ and (1.5$\pm$0.3)$\times10^{13}$ cm$^{-2}$, respectively. The chemical modeling predicts column densities of (0.5-2)$\times10^{13}$ cm$^{-2}$ in the envelope (depending on the cosmic-ray ionization rate), and 5$\times10^{11}$ to 2.5$\times10^{14}$ cm$^{-2}$ in the foreground cloud (depending on time). Both observed abundances are reproduced by the model at a satisfactory level. The constraints set by these observations on the physical conditions in the foreground cloud are however weak. Furthermore, the CH abundance in the envelope is strongly affected by the rate coefficient of the reaction H+CH$\rightarrow$C+H$_2$ ; further investigation of its value at low temperature would be necessary to facilitate the comparison between the model and the observations.

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-4 that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, 1e-12 - 1e-7. Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards several T Tauri star-disk systems. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

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 performed one-dimensional hydrodynamic simulations with detailed cooling, heating and chemical processes to examine the thermal stability of shocked gas in cold neutral medium (CNM) and molecular clouds. We find that both CNM and molecular clouds can be thermally unstable in the cooling layer behind the shock wave. The characteristic wavelength of the thermal instability ranges from 10^-5 pc to 0.1 pc in the CNM, and from 10^-7 pc to 0.1 pc in the molecular clouds. This coincides with the size of observed tiny scale structures in the CNM and molecular clouds, indicating that the thermal instability in the shocked gas could be a formation mechanism of these tiny structures in the interstellar medium. We have also calculated the e-folding number of the thermal instability to estimate the amplification of the density fluctuation in the shocked gas. Density perturbations in the CNM grow by a factor of exp(5)~150, whereas the perturbations in the molecular clouds grow only by a factor of a few behind a high Mach number shock. The amplification factor is larger at lower densities and higher velocities. Formation of very small-scale structures by thermal instability in shocked gas is more effective in lower densities.

We present near infrared spectroscopic observations of 19 molecular clouds made using the AKARI satellite, and the data reduction pipeline written to analyse those observations. The 2.5 --~5 $\mu$m spectra of 30 objects -- 22 field stars behind quiescent molecular clouds and eight low mass YSOs in cores -- were successfully extracted using the pipeline. Those spectra are further analysed to calculate the column densities of key solid phase molecular species, including H$_2$O, CO$_2$, CO, and OCN$^-$. The profile of the H$_2$O ice band is seen to vary across the objects observed and we suggest that the extended red wing may be an evolutionary indicator of both dust and ice mantle properties. The observation of 22 spectra with fluxes as low as $<$~5 mJy towards background stars, including 15 where the column densities of H$_2$O, CO and CO$_2$ were calculated, provides valuable data that could help to benchmark the initial conditions in star-forming regions prior to the onset of star formation.

We report the first interstellar detection of c-C3D2. The doubly deuterated cyclopropenylidene, a carbene, has been detected toward the starless cores TMC- 1C and L1544 using the IRAM 30m telescope. The J(Ka,Kc) = 3(0,3)-2(1,2), 3(1,3)-2(0,2), and 2(2,1)-1(1,0) transitions of this species have been observed at 3 mm in both sources. The expected 1:2 intensity ratio has been found in the 3(0,3)-2(1,2) and 3(1,3)-2(0,2) lines, belonging to the para and ortho species respectively. We also observed lines of the main species, c-C3H2, the singly deuterated c-C3HD, and the species with one 13C off of the principal axis of the molecule, c-H13CC2H. The lines of c-C3D2 have been observed with high signal to noise ratio, better than 7.5 sigma in TMC-1C and 9 sigma in L1544. The abundance of doubly deuterated cyclopropenylidene with respect to the normal species is found to be (0.4 - 0.8)% in TMC-1C and (1.2 - 2.1)% in L1544. The deuteration of this small hydrocarbon ring is analysed with a comprehensive gas-grain model, the first including doubly deuterated species. The observed abundances of c-C3D2 can be explained solely by gas-phase processes, supporting the idea that c-C3H2 is a good indicator of gas-phase deuteration.

We study the evolution of phosphorous-bearing species in one-dimensional C-shock models. We find that the abundances of P-bearing species depend sensitively on the elemental abundance of P in the gas phase and on the abundance of N atoms in the pre-shock gas. The observed abundance of PN and the non-detection of PO towards L1157 B1 are reproduced in C-shock models with shock velocity v=20km s-1 and pre-shock density n(H2) =10^4 - 10^5, if the elemental abundance of P in the gas phase is 10^-9 and the N-atom abundance is n(N)/nH - 10^-5 in the pre-shock gas. We also find that P-chemistry is sensitive to O- and N-chemistry, because N atoms are destroyed mainly by OH and NO. We identify the reactions of O-bearing and N-bearing species that significantly affect P chemistry.

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

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.

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

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 investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a layer, Z/R <~ 0.1 - 0.2, in the disk midplane, within R \approx 200 AU.

To investigate the composition and evolution of circumstellar ice around low-mass YSOs, we observed ice absorption bands in the near infrared (NIR) towards eight YSOs ranging from class 0 to class II, among which seven are associated with edge-on disks. We performed slit-less spectroscopic observations using the grism mode of the Infrared Camera (IRC) on board AKARI, which enables us to obtain full NIR spectra from 2.5 $\mu$m to 5 $\mu$m. The spectra were fitted with polynomial baselines to derive the absorption spectra. The molecular absorption bands were then fitted with the laboratory database of ice absorption bands, considering the instrumental line profile and the spectral resolution of the grism dispersion element. Towards the class 0-I sources (L1527, IRC-L1041-2, and IRAS04302), absorption bands of H$_2$O, CO$_2$, CO, and XCN are clearly detected. Column density ratios of CO$_2$ ice and CO ice relative to H$_2$O ice are 21-28% and 13-46%, respectively. If XCN is OCN$^-$, its column density is as high as 2-6% relative to H$_2$O ice. The HDO ice feature at 4.1 $\mu$m is tentatively detected towards the class 0-I sources and HV Tau. Non-detections of the CH-stretching mode features around 3.5 $\mu$m provide upper limits to the CH$_3$OH abundance of 26% (L1527) and 42% (IRAS04302) relative to H$_2$O. We tentatively detect OCS ice absorption towards IRC-L1041-2. Towards class 0-I sources, the detected features should mostly originate in the cold envelope, while CO gas and OCN$^-$ could originate in the region close to the protostar, where there are warm temperatures and UV radiation. We detect H$_2$O ice band towards ASR41 and 2MASSJ1628137-243139, which are edge-on class II disks. We also detect H$_2$O ice and CO$_2$ ice towards HV Tau, HK Tau, and UY Aur, and tentatively detect CO gas features towards HK Tau and UY Aur.

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.

We have observed dense gas around the Very Low-Luminosity Ob jects (VeLLOs) L1521F-IRS and IRAM 04191+1522 in carbon-chain and organic molecular lines with the Nobeyama 45 m telescope. Towards L1521F-IRS, carbon-chain lines of CH3CCH (50-40), C4H (17/2-15/2), and C3H2 (212-101) are 1.5 - 3.5 times stronger than those towards IRAM 04191+1522, and the abundances of the carbon-chain molecules towards L1521F-IRS are 2 to 5 times higher than those towards IRAM 04191+1522. Mapping observations of these carbon-chain molecular lines show that in L1521F the peak positions of these carbon-chain molecular lines are different from each other and there is no emission peak towards the VeLLO position, while in IRAM 04191+1522 these carbon-chain lines are as weak as the detection limits except for the C3H2 line. The observed chemical differentiation between L1521F and IRAM 04191+1522 suggests that the evolutionary stage of L1521F-IRS is younger than that of IRAM 04191+1522, consistent with the extent of the associated outflows seen in the 13CO (1-0) line. The non-detection of the organic molecular lines of CH3OH (6-2-7-1 E) and CH3CN (60-50) implies that the warm (~ 100 K) molecular-desorbing region heated by the central protostar is smaller than ~ 100 AU towards L1521F-IRS and IRAM 04191+1522, suggesting the young age of these VeLLOs. We propose that the chemical status of surrounding dense gas can be used to trace the evolutionary stages of VeLLOs.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

We have studied the dependence of the profiles of molecular abundances and line emission on the accretion flow in the hot ($\ga 100$K) inner region of protoplanetary disks. The gas-phase reactions initiated by evaporation of the ice mantle on dust grains are calculated along the accretion flow. We focus on methanol, a molecule that is formed predominantly through the evaporation of warm ice mantles, to show how the abundance profile and line emission depend on the accretion flow. Our results show that some evaporated molecules keep high abundances only when the accretion velocity is large enough, and that methanol could be useful as a diagnostic of the accretion flow by means of ALMA observations at the disk radius of $\la 10$AU.

We present Submillimeter Array (SMA) observations of several deuterated species in the disk around the classical T Tauri star TW Hydrae at arcsecond scales, including detections of the DCN J=3-2 and DCO+ J=3-2 lines, and upper limits to the HDO 3(1,2)-2(2,1), ortho-H2D+ 1(1,0)-1(1,1) and para-D2H+ 1(1,0)-1(0,1) transitions. We also present observations of the HCN J=3-2, HCO+ J=3-2 and H13CO+ J=4-3 lines for comparison with their deuterated isotopologues. We constrain the radial and vertical distributions of various species in the disk by fitting the data using a model where the molecular emission from an irradiated accretion disk is sampled with a 2D Monte Carlo radiative transfer code. We find that the distribution of DCO+ differs markedly from that of HCO+. The D/H ratios inferred change by at least one order of magnitude (0.01 to 0.1) for radii <30 AU to >70 AU and there is a rapid falloff of the abundance of DCO+ at radii larger than 90 AU. Using a simple analytical chemical model, we constrain the degree of ionization, x(e-)=n(e-)/n(H2), to be ~10^-7 in the disk layer(s) where these molecules are present. Provided the distribution of DCN follows that of HCN, the ratio of DCN to HCN is determined to be 1.7\pm0.5 \times 10^-2; however, this ratio is very sensitive to the poorly constrained vertical distribution of HCN. The resolved radial distribution of DCO+ indicates that {\it in situ} deuterium fractionation remains active within the TW Hydrae disk and must be considered in the molecular evolution of circumstellar accretion disks.

The millimeter-wave rotational emission lines (4(04)-3(03) and 5(05)-4(04)) of protonated carbon dioxide, HCO2+(HOCO+), has been detected toward the low-mass class 0 protostar IRAS 04368+2557 in L1527 with the IRAM 30 m telescope. This is the first detection of HCO2+ except for the Galactic Center clouds. The column density of HCO2+ averaged over the beam size (29") is determined to be 7.6x10^10 cm^-2, assuming the rotational temperature of 12.3 K. The fractional abundance of gaseous CO2 relative to H2 is estimated from the column density of HCO2+ with an aid of a simplified chemical model. If the HCO2+ emission only comes from the evaporation region of CO2 near the protostar (T>50 K), the fractional abundance of CO2 is estimated to be higher than 6.6x10^-4. This is comparable to the elemental abundance of carbon in interstellar clouds, and hence, the direct evaporation of CO2 from dust grain is unrealistic as a source of gaseous CO2 in L1527. A narrow line width of HCO2+ also supports this. On the other hand, the fractional abundance of CO2 is estimated to be 2.9x10^-7, if the source size is comparable to the beam size. These results indicate that gaseous CO2 is abundant even in the low-mass star-forming region. Possible production mechanisms of gaseous CO2 are discussed.

We investigate molecular evolution in a star-forming core that is initially a hydrostatic starless core and collapses to form a low-mass protostar. The results of a one-dimensional radiation-hydrodynamics calculation are adopted as a physical model of the core. We first derive radii at which CO and large organic species sublimate. CO sublimation in the central region starts shortly before the formation of the first hydrostatic core. When the protostar is born, the CO sublimation radius extends to 100 AU, and the region inside $\lesssim 10$ AU is hotter than 100 K, at which some large organic species evaporate. We calculate the temporal variation of physical parameters in infalling shells, in which the molecular evolution is solved using an updated gas-grain chemical model to derive the spatial distribution of molecules in a protostellar core. The shells pass through the warm region of $10 -100$ K in several $\times$ $10^4$ yr, and fall into the central star $\sim 100$ yr after they enter the region where $T \gtrsim 100$ K. We find that large organic species are formed mainly via grain-surface reactions at temperatures of $20 -40$ K and then desorbed into the gas-phase at their sublimation temperatures. Carbon-chain species can be formed by a combination of gas-phase reactions and grain-surface reactions following the sublimation of CH$_4$. Our model also predicts that CO$_2$ is more abundant in isolated cores, while gas-phase large organic species are more abundant in cores embedded in ambient clouds.

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and ultraviolet (UV) irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains in the disks are numerically computed by solving the coagulation equation for settling dust particles. The level populations and line emission of molecular hydrogen are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk region and in the surface layer, while the far UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-ray induced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the properties of the dust grains in the disks. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient, while the UV radiation fields become stronger due to the decrease of grain opacity. This makes the H2 level populations change from local thermodynamic equilibrium (LTE) to non-LTE distributions, which results in changes to the line ratios of H2 emission. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.

In a disk around DM Tau, previous observation of 13CO (J=2-1 and 1-0 transitions) derived the 13CO gas temperature of \sim 13-20K, which is lower than the sublimation temperature of CO (20 K). We argue that the existence of such cold CO can be explained by a vertical mixing of disk material. As the gas is transported from a warm layer to a cold layer, CO is depleted onto dust grains with a timescale of \sim 10^3 yr. Because of the steep temperature gradient in the vertical direction, an observable amount of CO is still in the gas phase when the fluid parcel reaches the layer of \sim 13 K. Apparent temperature of CO decreases as the maximum grain size increases from micron-size to mm-size.

In this review we re-evaluate our observational and theoretical understanding of the chemical evolution of protoplanetary disks. We discuss how improved observational capabilities have enabled the detection of numerous molecules exposing an active disk chemistry that appears to be in disequilibrium. We outline the primary facets of static and dynamical theoretical chemical models. Such models have demonstrated that the observed disk chemistry arises from warm surface layers that are irradiated by X-ray and FUV emission from the central accreting star. Key emphasis is placed on reviewing areas where disk chemistry and physics are linked: including the deuterium chemistry, gas temperature structure, disk viscous evolution (mixing), ionization fraction, and the beginnings of planet formation.

We calculate the physical structure of protoplanetary disks by evaluating the gas density and temperature self-consistently and solving separately for the dust temperature. The effect of grain growth is taken into account by assuming a power-law size distribution and varying the maximum radius of grains a_max. In our fiducial model with a_max=10um, the gas is warmer than the dust in the surface layer of the disk, while the gas and dust have the same temperature in deeper layers. In the models with larger a_max, the gas temperature in the surface layer is lower than in the fiducial model because of reduced photo-electric heating rates from small grains, while the deeper penetration of stellar radiation warms the gas at intermediate height. A detailed chemical reaction network is solved at outer radii (r \ge 50 AU). Vertical distributions of some molecular species at different radii are similar, when plotted as a function of hydrogen column density Sigma_H from the disk surface. Consequently, molecular column densities do not much depend on disk radius. In the models with larger a_max, the lower temperature in the surface layer makes the geometrical thickness of the disk smaller, and the gaseous molecules are confined to smaller heights. However, if we plot the vertical distributions of molecules as a function of Sigma_H, they do not significantly depend on a_max. The dependence of the molecular column densities on a_max is not significant, either. Notable exceptions are HCO+, H3+ and H2D+, which have smaller column densities in the models with larger a_max.

On the basis of near-infrared imaging observations, we derived visual extinction (Av) distribution toward ten Bok globules through measurements of both the color excess (E_{H-K}) and the stellar density at J, H, and Ks (star count). Radial column density profiles for each globule were analyzed with the Bonnor-Ebert sphere model. Using the data of our ten globules and four globules in the literature, we investigated the stability of globules on the basis of xi_max, which characterizes the Bonnor-Ebert sphere as well as the stability of the equilibrium state against the gravitational collapse. We found that more than half of starless globules are located near the critical state (xi_max = 6.5 +/- 2). Thus, we suggest that a nearly critical Bonnor-Ebert sphere characterizes the typical density structure of starless globules. Remaining starless globules show clearly unstable states (xi_max > 10). Since unstable equilibrium states are not long maintained, we expect that these globules are on the way to gravitational collapse or that they are stabilized by non-thermal support. It was also found that all the star-forming globules show unstable solutions of xi_max >10, which is consistent with the fact that they have started gravitational collapse. We investigated the evolution of a collapsing gas sphere whose initial condition is a nearly critical Bonnor-Ebert sphere. We found that the column density profiles of the collapsing sphere mimic those of the static Bonnor-Ebert spheres in unstable equilibrium. The collapsing gas sphere resembles marginally unstable Bonnor-Ebert spheres for a long time. We found that the frequency distribution of xi_max for the observed starless globules is consistent with that from model calculations of the collapsing sphere.

The gravitational collapse of a spherical cloud core is investigated by numerical calculations. The initial conditions of the core lie close to the critical Bonnor-Ebert sphere with a central density of \sim 10^4 cm^{-3} in one model (alpha=1.1), while gravity overwhelms pressure in the other (alpha=4.0), where alpha is the internal gravity-to-pressure ratio. The alpha=1.1 model shows reasonable agreement with the observed velocity field in prestellar cores. Molecular distributions in cores are calculated by solving a chemical reaction network that includes both gas-phase and grain-surface reactions. When the central density of the core reaches 10^5 cm^{-3}, carbon-bearing species are significantly depleted in the central region of the alpha=1.1 model, while the depletion is only marginal in the other model. The two different approaches encompass the observed variations of molecular distributions in different prestellar cores, suggesting that molecular distributions can be probes of contraction or accumulation time scales of cores. The central enhancement of the NH3/N2H+ ratio, which is observed in some prestellar cores, can be reproduced under certain conditions by adopting recently measured branching fractions for N2H+ recombination. Various molecular species, such as CH3OH and CO2, are produced by grain-surface reactions. The ice composition depends sensitively on the assumed temperature. Multi-deuterated species are included in our most recent gas-grain chemical network. The deuterated isotopomers of H3+ are useful as probes of the central regions of evolved cores, in which gas-phase species with heavy elements are strongly depleted. At 10 K, our model can reproduce the observed abundance ratio of ND3/NH3, but underestimates the isotopic ratios of deuterated to normal methanol.

The molecular evolution that occurs in collapsing prestellar cores is investigated. To model the dynamics, we adopt the Larson-Penston (L-P) solution and analogues with slower rates of collapse. For the chemistry, we utilize the new standard model (NSM) with the addition of deuterium fractionation and grain-surface reactions treated via the modified rate approach. The use of surface reactions distinguishes the present work from our previous model. We find that these reactions efficiently produce H2O, H2CO, CH3OH, N2, and NH3 ices. In addition, the surface chemistry influences the gas-phase abundances in a variety of ways. The current reaction network along with the L-P solution allows us to reproduce satisfactorily most of the molecular column densities and their radial distributions observed in L1544. The agreement tends to worsen with models that include strongly delayed collapse rates. Inferred radial distributions in terms of fractional abundances are somewhat harder to reproduce. In addition to our standard chemical model, we have also run a model with the UMIST gas-phase chemical network. The abundances of gas-phase S-bearing molecules such as CS and CCS are significantly affected by uncertainties in the gas-phase chemical network. In all of our models, the column density of N2H+ monotonically increases as the central density of the core increases during collapse from 3 10^4 cm-3 to 3 10^7 cm-3. Thus, the abundance of this ion can be a probe of evolutionary stage. Molecular D/H ratios in assorted cores are best reproduced in the L-P picture with the conventional rate coefficients for fractionation reactions. If we adopt the newly measured and calculated rate coefficients, the D/H ratios, especially N2D+/N2H+, become significantly lower than the observed values.

Emission from the $2_{12}-1_{11}$ line of H$_2$CO has been detected and marginally resolved toward LkCa15 by the Nobeyama Millimeter Array. The column density of H$_2$CO is higher than that observed in DM Tau and than predicted by theoretical models of disk chemistry; also the line-intensity profile is less centrally peaked than that for CO. A similar behavior is observed in other organic gaseous molecules in the LkCa 15 disk.

A new two-dimensional axi-symmetric ultraviolet radiative transfer code is presented, which is used to calculate photodissociation and ionization rates for use in chemistry models of flaring circumstellar disks. Scattering and absorption of photons from the central star and from the interstellar radiation field are taken into account. The molecules are effectively photodissociated in the surface layer of the disk, but can exist in the intermediate, moderately warm layers. A comparison has been made with an approximate 2D ray-tracing method and it was found that the latter underestimates the ultraviolet field and thus the molecular photodissociation rates below the disk surface. The full 2D results show significantly higher abundances of radicals such as CN and C2H than previous work, partly due to the fact that CO is dissociated to greater depths. Results for different stellar radiation fields are also presented. The CN/HCN ratio shows a strong dependence on the stellar spectrum, whereas other ratios such as HCO+/CO show only little variation.

We have investigated molecular distributions in protoplanetary disks, adopting a disk model with a temperature gradient in the vertical direction. The model produces sufficiently high abundances of gaseous CO and HCO+ to account for line observations of T Tauri stars using a sticking probability of unity and without assuming any non-thermal desorption. In regions of radius R > 10 AU, with which we are concerned, the temperature increases with increasing height from the midplane. In a warm intermediate layer, there are significant amounts of gaseous molecules owing to thermal desorption and efficient shielding of ultraviolet radiation by the flared disk. The column densities of HCN, CN, CS, H2CO, HNC and HCO+ obtained from our model are in good agreement with the observations of DM Tau, but are smaller than those of LkCa15. Molecular line profiles from our disk models are calculated using a 2-dimensional non-local-thermal-equilibrium (NLTE) molecular-line radiative transfer code for a direct comparison with observations. Deuterated species are included in our chemical model. The molecular D/H ratios in the model are in reasonable agreement with those observed in protoplanetary disks.

We have investigated the evolution and distribution of molecules in collapsing prestellar cores via numerical chemical models, adopting the Larson-Penston solution and its delayed analogues to study collapse. Molecular abundances and distributions in a collapsing core are determined by the balance among the dynamical, chemical and adsorption time scales. When the central density n_H of a prestellar core with the Larson-Penston flow rises to 3 10^6 cm^{-3}, the CCS and CO column densities are calculated to show central holes of radius 7000 AU and 4000 AU, respectively, while the column density of N2H+ is centrally peaked. These predictions are consistent with observations of L1544. If the dynamical time scale of the core is larger than that of the Larson-Penston solution owing to magnetic fields, rotation, or turbulence, the column densities of CO and CCS are smaller, and their holes are larger than in the Larson-Penston core with the same central gas density. On the other hand, N2H+ and NH3 are more abundant in the more slowly collapsing core. Therefore, molecular distributions can probe the collapse time scale of prestellar cores. Deuterium fractionation has also been studied via numerical calculations. The deuterium fraction in molecules increases as a core evolves and molecular depletion onto grains proceeds. When the central density of the core is n_H=3 10^6 cm^{-3}, the ratio DCO+/HCO+ at the center is in the range 0.06-0.27, depending on the collapse time scale and adsorption energy; this range is in reasonable agreement with the observed value in L1544.

We have investigated the two-dimensional (R,Z) distribution of deuterated molecular species in circumstellar disks around young stellar objects. The abundance ratios between singly deuterated and normal molecules (``D/H ratios'') in disks evolve in a similar way as in molecular clouds. Fractionation is caused by rapid exchange reactions that are exothermic because of energy differences between deuterated and normal species. In the midplane region, where molecules are heavily depleted onto grain surfaces, the D/H ratios of gaseous molecules are higher than at larger heights. The D/H ratios for the vertical column densities of NH3, H2O, and HCO+ are sensitive to the temperature, and decrease significantly with decreasing radial distance for R < 300 AU. The analogous D/H ratios for CH4 and H2CO, on the other hand, are not very sensitive to the temperature in the range (T=10-50 K) we are concerned with, and do not decrease with decreasing R at R > 50 AU. The D/H column-density ratios also depend on disk mass. In a disk with a larger mass, the ratios of deuterated species to normal species are higher, because of heavier depletion of molecules onto grains. In the second part of the paper, we report molecular column densities for disks embedded in ambient cloud gas. Our results suggest that CN and HCO+ can be tracers of gaseous disks, especially if the central object is a strong X-ray source. Our results also suggest that the radial distributions of CN, C2H, HCN, and H2CO may vary among disks depending on the X-ray luminosity of the central star.

We investigate the evolution of molecular abundance in quiescent protoplanetary disks which are presumed to be around weak-line T Tauri stars. In the region of surface density less than $10^2$ g cm$^{-2}$ (distance from the star $\gtrsim 10$ AU in the minimum- mass solar nebula), cosmic rays are barely attenuated even in the midplane of the disk and produce chemically active ions such as He$^+$ and H$_{3}^+$. Through reactions with these ions CO and N$_2$ are finally transformed into CO$_2$, NH$_3$, and HCN. In the region where the temperature is low enough for these products to freeze onto grains, considerable amount of carbon and nitrogen is locked up in the ice mantle and is depleted from the gas phase in a time scale $\lesssim 3\times 10^6$ yr. Oxidized (CO$_2$) ice and reduced (NH$_3$ and hydrocarbon) ice naturally coexist in this part of the disk. The molecular abundance both in the gas phase and in ice mantle varies significantly with the distance from the central star.