The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

AGB stars are, together with supernovae, the main contributors of stellar dust to the interstellar medium (ISM). Dust grains formed by AGB stars are thought to be large. However, as dust nucleation and growth within their outflows are still not understood, the dust-grain size distribution (GSD) is unknown. This is an important uncertainty regarding our knowledge of the chemical and physical history of interstellar dust, as AGB dust forms $\sim$ 70% of the starting point of its evolution. We expand on our chemical kinetics model, which uniquely includes a comprehensive dust-gas chemistry. The GSD is now allowed to deviate from the commonly assumed canonical Mathis et al. (1977) distribution. We find that the specific GSD can significantly influence the dust-gas chemistry within the outflow. Our results show that the level of depletion of gas-phase species depends on the average grain surface area of the GSD. Gas-phase abundance profiles and their possible depletions can be retrieved from observations of molecular emission lines when using a range of transitions. Due to degeneracies within the prescription of GSD, specific parameters cannot be retrieved, only (a lower limit to) the average grain surface area. Nonetheless, this can discriminate between dust composed of predominantly large or small grains. We show that when combined with other observables such as the spectral energy distribution and polarised light, depletion levels from molecular gas-phase abundance profiles can constrain the elusive GSD of the dust delivered to the ISM by AGB outflows.

Chemical modelling of AGB outflows is typically focused on either non-thermodynamic equilibrium chemistry in the inner region or photon-driven chemistry in the outer region. We include, for the first time, a comprehensive dust-gas chemistry in our AGB outflow chemical kinetics model, including both dust-gas interactions and grain-surface chemistry. The dust is assumed to have formed in the inner region, and follows an interstellar-like dust-size distribution. Using radiative transfer modelling, we obtain dust temperature profiles for different dust types in an O-rich and a C-rich outflow. We calculate a grid of models, sampling different outflow densities, drift velocities between the dust and gas, and dust types. Dust-gas chemistry can significantly affect the gas-phase composition, depleting parent and daughter species and increasing the abundance of certain daughter species via grain-surface formation followed by desorption/sputtering. Its influence depends on four factors: outflow density, dust temperature, initial composition, and drift velocity. The largest effects are for higher density outflows with cold dust and O-rich parent species, as these species generally have a larger binding energy. At drift velocities larger than $\sim 10$ km s$^{-1}$, ice mantles undergo sputtering; however, they are not fully destroyed. Models with dust-gas chemistry can better reproduce the observed depletion of species in O-rich outflows. When including colder dust in the C-rich outflows and adjusting the binding energy of CS, the depletion in C-rich outflows is also better reproduced. To best interpret high-resolution molecular line observations from AGB outflows, dust-gas interactions are needed in chemical kinetics models.

We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$\sigma$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$\mu$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

Episodic accretion has been used to explain the wide range of protostellar luminosities, but its origin and influence on the star forming process are not yet fully understood. We present an ALMA survey of N$_2$H$^+$ ($1-0$) and HCO$^+$ ($3-2$) toward 39 Class 0 and Class I sources in the Perseus molecular cloud. N$_2$H$^+$ and HCO$^+$ are destroyed via gas-phase reactions with CO and H$_2$O, respectively, thus tracing the CO and H$_2$O snowline locations. A snowline location at a much larger radius than that expected from the current luminosity suggests that an accretion burst has occurred in the past which has shifted the snowline outward. We identified 18/18 Class 0 and 9/10 Class I post-burst sources from N$_2$H$^+$, and 7/17 Class 0 and 1/8 Class I post-burst sources from HCO$^+$.The accretion luminosities during the past bursts are found to be $\sim10-100~L_\odot$. This result can be interpreted as either evolution of burst frequency or disk evolution. In the former case, assuming that refreeze-out timescales are 1000 yr for \ce{H2O} and 10,000 yr for CO, we found that the intervals between bursts increases from 2400 yr in the Class 0 to 8000 yr in the Class I stage. This decrease in the burst frequency may reflect that fragmentation is more likely to occur at an earlier evolutionary stage when the young stellar object is more prone to instability.

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

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

[Abridged] With a growing number of molecules observed in many comets, and an improved understanding of chemical evolution in protoplanetary disk midplanes, comparisons can be made between models and observations that could potentially constrain the formation histories of comets. A $\chi^{2}$-method was used to determine maximum likelihood surfaces for 14 different comets that formed at a given time (up to 8 Myr) and place (out to beyond the CO iceline) in the pre-solar nebula midplane. This was done using observed volatile abundances for the 14 comets and the evolution of volatile abundances from chemical modelling of disk midplanes. Considering all parent species (ten molecules) in a scenario that assumed reset initial chemistry, the $\chi^{2}$ likelihood surfaces show a characteristic trail in the parameter space with high likelihood of formation around 30 AU at early times and 12 AU at later times for ten comets. This trail roughly traces the vicinity of the CO iceline in time. The formation histories for all comets were thereby constrained to the vicinity of the CO iceline, assuming that the chemistry was partially reset early in the pre-solar nebula. This is found, both when considering carbon-, oxygen-, and sulphur-bearing molecules (ten in total), and when only considering carbon- and oxygen-bearing molecules (seven in total). Since these 14 comets did not previously fall into the same taxonomical categories together, this chemical constraint may be proposed as an alternative taxonomy for comets. Based on the most likely time for each of these comets to have formed during the disk chemical evolution, a formation time classification for the 14 comets is suggested.

Observations of different molecular lines in protoplanetary disks provide valuable information on the gas kinematics, as well as constraints on the radial density and temperature structure of the gas. With ALMA we have detected H13CO+ (J=4-3) and HC15N (J=4-3) in the HD97048 protoplanetary disk for the first time. We compare these new detections to the ringed continuum mm-dust emission and the spatially resolved CO (J=3-2) and HCO+ (J=4-3) emission. The radial distributions of the H13CO+ and HC15N emission show hints of ringed sub-structure whereas, the optically thick tracers, CO and HCO+, do not. We calculate the HCO+/H13CO+ intensity ratio across the disk and find that it is radially constant (within our uncertainties). We use a physio-chemical parametric disk structure of the HD97048 disk with an analytical prescription for the HCO+ abundance distribution to generate synthetic observations of the HCO+ and H13CO+ disk emission assuming LTE. The best by-eye fit models require radial variations in the HCO+/H13CO+ abundance ratio and an overall enhancement in H13CO+ relative to HCO+. This highlights the need to consider isotope selective chemistry and in particular low temperature carbon isotope exchange reactions. This also points to the presence of a reservoir of cold molecular gas in the outer disk (T < 10K, R > 200au). Chemical models are required to confirm that isotope-selective chemistry alone can explain the observations presented here. With these data, we cannot rule out that the known dust substructure in the HD97048 disk is responsible for the observed trends in molecular line emission, and higher spatial resolution observations are required to fully explore the potential of optically thin tracers to probe planet-carved dust gaps. We also report non-detections of H13CO+ and HC15N in the HD100546 protoplanetary disk.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

The dominant reservoirs of elemental nitrogen in protoplanetary disks have not yet been observationally identified. Likely candidates are HCN, NH$_3$ and N$_2$. The relative abundances of these carriers determine the composition of planetesimals as a function of disk radius due to strong differences in their volatility. A significant sequestration of nitrogen in carriers less volatile than N$_2$ is likely required to deliver even small amounts of nitrogen to the Earth and potentially habitable exo-planets. While HCN has been detected in small amounts in inner disks ($<10$ au), so far only relatively insensitive upper limits on inner disk NH$_3$ have been obtained. We present new Gemini-TEXES high resolution spectroscopy of the 10.75 $\mu$m band of warm NH$_3$, and use 2-dimensional radiative transfer modeling to improve previous upper limits by an order of magnitude to $\rm [NH_3/H_{nuc}]<10^{-7}$ at 1 au. These NH$_3$ abundances are significantly lower than those typical for ices in circumstellar envelopes ($[{\rm NH_3/H_{nuc}}]\sim 3\times 10^{-6}$). We also consistently retrieve the inner disk HCN gas abundances using archival Spitzer spectra, and derive upper limits on the HCN ice abundance in protostellar envelopes using archival ground-based 4.7 $\mu$m spectroscopy ([HCN$_{\rm ice}$]/[H$_2$O$_{\rm ice}$]$<1.5-9$\%). We identify the NH$_3$/HCN ratio as an indicator of chemical evolution in the disk, and use this ratio to suggest that inner disk nitrogen is efficiently converted from NH$_3$ to N$_2$, significantly increasing the volatility of nitrogen in planet-forming regions.

Methanol (CH$_3$OH) is at the root of organic ice chemistry in protoplanetary disks. However, its weak emission has made detections difficult. To date, gas-phase CH$_3$OH has been detected in only one Class II disk, TW Hya. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to search for a total of four CH$_3$OH emission lines in bands 6 and 7 toward the disk around the young Herbig Ae star HD 163296. The disk-averaged column density of methanol and its related species formaldehyde (H$_2$CO) are estimated assuming optically thin emission in local thermodynamic equilibrium. We compare these results to the gas-phase column densities of the TW Hya disk. No targeted methanol lines were detected individually nor after line stacking. The 3$\sigma$ disk-integrated intensity upper limits are $< 51$ mJy km s$^{-1}$ for the band 6 lines and $< 26$ mJy km s$^{-1}$ for the band 7 lines. The band 7 lines provide the strictest 3$\sigma$ upper limit on disk-averaged column density with $N_{\mathrm{avg}} < 5.0 \times 10^{11}$ cm$^{-2}$. The methanol-to-formaldehyde ratio is CH$_3$OH/H$_2$CO $< 0.24$ in the HD 163296 disk compared to a ratio of $1.27$ in the TW Hya disk. Differences in the stellar irradiation of Herbig disks compared to T Tauri disks likely influence the gaseous methanol and formaldehyde content. Possible reasons for the lower HD 163296 methanol-to-formaldehyde ratio include: a higher than expected gas-phase formation of H$_2$CO in the HD 163296 disk, uncertainties in the grain surface formation efficiency of CH$_3$OH and H$_2$CO, and differences in the disk structure and/or CH$_3$OH and H$_2$CO desorption processes that release the molecules from ice mantles back into the gas phase. These results provide observational evidence that the gas-phase chemical complexity found in disks may be strongly influenced by the spectral type of the host star.

The bulk composition of the Earth is dramatically carbon poor compared to that of the interstellar medium, and this phenomenon extends to the asteroid belt. To interpret this carbon deficit problem, the carbonaceous component in grains must have been converted into the gas-phase in the inner regions of protoplanetary disks prior to planetary formation. We examine the effect of carbon grain destruction on the chemical structure of disks by calculating the molecular abundances and distributions using a comprehensive chemical reaction network. When carbon grains are destroyed and the elemental abundance of the gas becomes carbon-rich, the abundances of carbon-bearing molecules, such as HCN and carbon-chain molecules, increase dramatically near the midplane, while oxygen-bearing molecules, such as H$_2$O and CO$_2$, are depleted. We compare the results of these model calculations with the solid carbon fraction in the solar system. Although we find a carbon depletion gradient, there are some quantitative discrepancies: the model shows a higher value at the position of asteroid belt and a lower value at the location of the Earth. In addition, using the obtained molecular abundances distributions, coupled with line radiative transfer calculations, we make predictions for ALMA to potentially observe the effect of carbon grain destruction in nearby protoplanetary disks. The results indicate that HCN, H$^{13}$CN and c-C$_3$H$_2$ may be good tracers.

Simulations of massive star formation predict the formation of discs with significant substructure, such as spiral arms and clumps due to fragmentation. Here we present a semi-analytic framework for producing synthetic observations of discs with substructure, in order to determine their observability in interferometric observations. Unlike post-processing of hydrodynamical models, the speed inherent to our approach permits us to explore a large parameter space of star and disc parameters, and thus constrain properties for real observations. We compute synthetic dust continuum and molecular line observations probing different disc masses, distances, inclinations, thermal structures, dust distributions, and number and orientation of spirals and fragments. With appropriate spatial and kinematic filtering applied, our models predict that ALMA observations of massive YSOs at <5 kpc distances should detect spirals in both gas and dust in strongly self-gravitating discs (i.e. discs with up to two spiral arms and strong kinematic perturbations). Detecting spirals will be possible in discs of arbitrary inclination, either by directly spatially resolving them for more face-on discs (inclinations up to ~50 degrees), or through a kinematic signature otherwise. Clumps resulting from disc fragmentation should be detectable in the continuum, if the clump is sufficiently hotter than the surrounding disc material.

CO has long been thought to be the best tracer to measure gas masses as it is readily detected at (sub)mm wavelengths in many disks. Inferred gas masses from CO in recent ALMA observations of large samples of disks seem inconsistent with their inferred dust masses. The derived gas-to-dust mass ratios from CO are 1-2 orders of magnitude lower than $\sim$100 even if photodissociation and freeze-out are included. Herschel measurements of HD line emission imply gas masses in line with gas-to-dust mass ratios of 100. This suggests that at least one additional mechanism is removing gaseous CO. Here we test the suggestion that the bulk of the CO is chemically processed and that the carbon is sequestered into less volatile species such as CO2, CH3OH and CH4 in the dense, shielded midplane regions of the disk. Using our gas-grain chemical code we perform a parameter exploration and follow the CO abundance evolution over a range of conditions representative of shielded disk midplanes. We find that no chemical processing of CO takes place on 1-3 Myr timescales for low cosmic-ray ionisation rates, $< 5\times 10^{-18}$ s$^{-1}$. Assuming an ionisation rate of $10^{-17}$ s$^{-1}$, more than 90% of the CO is destroyed, but only in the parts of the disk below 30 K. This order of magnitude destruction of CO is robust against the choice of grain-surface reaction rate parameters except for temperatures below 20 K, where there is a strong dependence on the H tunnelling barrier. The low temperatures needed for CO chemical processing indicate that the exact disk temperature structure is important, with warm disks around luminous Herbig stars expected to have little to no CO conversion. For cold disks around T Tauri stars, a large fraction of the emitting CO layer is affected unless the disks are young (< 1 Myr). This can lead to inferred gas masses that are up to two orders of magnitude too low.

[Abridged] The detection of abundant O$_{2}$ ice at 1-10% with respect to H$_{2}$O ice in the comae of comets 1P/Halley and 67P/Churyumov-Gerasimenko motivated attempts to explain the origin of the high O$_{2}$ ice abundance. Recent chemical modelling of the outer, colder regions of a protoplanetary disk midplane has shown production of O$_{2}$ ice at the same abundance as that measured in the comet. An updated chemical kinetics code is utilised here to evolve chemistry under pre-Solar nebula midplane conditions. Four different chemical starting conditions, and the effects of various chemical parameters are tested. The parameter space investigation revealed a sweet spot for production of O$_{2}$ ice at an abundance matching those in 67P and 1P, and O$_{3}$ and H$_{2}$O$_{2}$ ices abundances matching those in 67P. This means that there is a radial region in the pre-Solar nebula from 120-150 AU, within which O$_{2}$ could have been produced in-situ via ice chemistry on grain surfaces. However, it is apparent that there is a high degree of sensitivity of the chemistry to the assumed chemical parameters (e.g. binding energy, activation barrier width, and quantum tunnelling barrier). Hence, because the more likely scenario starting with a percentage of elemental oxygen locked in O$_{2}$ also reproduces the O$_{2}$ ice abundance in 67P at early stages, this supports previous suggestions that the cometary O$_{2}$ ice could have a primordial origin.

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

In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$.

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.

In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dust temperature, and gas density. We performed gas-grain astrochemical modeling experiments with exponentially increasing (coupled) gas and dust temperature rising from 10 to 500 K at constant H$_2$ densities of 10$^7$, 10$^8$, and 10$^9$ cm$^{-3}$. We tested the effect of varying the initial ice composition, cosmic-ray ionization rate, warm-up time (over 50, 200, and 1000 kyr), and initial (10, 15, and 25 K) and final temperatures (300 and 500 K). Varying the initial ice compositions within the observed and expected ranges does not noticeably affect the modeled abundances indicating that the chemical make-up of hot cores is determined in the warm-up stage. Complex cyanides vinyl and ethyl cyanide (CH$_2$CHCN and C$_2$H$_5$CN, respectively) cannot be produced in abundances (versus H$_2$) greater than 5x10$^{-10}$ for CH$_2$CHCN and 2x10$^{-10}$ for C$_2$H$_5$CN with a fast warm-up time (52 kyr), while the lower limit for the observed abundance of C$_2$H$_5$CN toward source B3 is 3.4x10$^{-10}$. Complex cyanide abundances are reduced at higher initial temperatures and increased at higher cosmic-ray ionization rates. Reproducing the observed abundances toward G35.20-0.74 Core B3 requires a fast warm-up at a high cosmic-ray ionization rate (1x10$^{-16}$ s$^{-1}$) at a high gas density (>10$^9$ cm$^{-3}$). G35.20-0.74 source B3 only needs to be about 2000 years older than B1/B2 for the observed chemical difference to be present. (This abstract has been shortened)

Modern radio interferometers enable observations of spectral lines with unprecedented spatial resolution and sensitivity. In spite of these technical advances, many lines of interest are still at best weakly detected and therefore necessitate detection and analysis techniques specialized for the low signal-to-noise ratio (SNR) regime. Matched filters can leverage knowledge of the source structure and kinematics to increase sensitivity of spectral line observations. Application of the filter in the native Fourier domain improves SNR while simultaneously avoiding the computational cost and ambiguities associated with imaging, making matched filtering a fast and robust method for weak spectral line detection. We demonstrate how an approximate matched filter can be constructed from a previously observed line or from a model of the source, and we show how this filter can be used to robustly infer a detection significance for weak spectral lines. When applied to ALMA Cycle 2 observations of CH3OH in the protoplanetary disk around TW Hya, the technique yields a ~53% SNR boost over aperture-based spectral extraction methods, and we show that an even higher boost will be achieved for observations at higher spatial resolution. A Python-based open-source implementation of this technique is available under the MIT license at https://github.com/AstroChem/VISIBLE.

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii <30 AU and a broad, ring-like structure from ~50-230 AU with a peak at 100 AU just beyond the millimeter grain edge. We modeled DCO$^+$ emission in HD 169142 with a physical disk structure adapted from the literature, and employed a simple deuterium chemical network to investigate the formation of DCO$^+$ through the cold deuterium fractionation pathway via H$_2$D$^+$. Contributions from the warm deuterium fractionation pathway via CH$_2$D$^+$ are approximated using a constant abundance in the intermediate disk layers. Parameterized models show that alterations to the midplane gas temperature and CO abundance of the literature model are both needed to recover the observed DCO$^+$ radial intensity profile. The best-fit model contains a shadowed, cold midplane in the region z/r < 0.1 with an 8 K decrease in gas temperature and a factor of five CO depletion just beyond the millimeter grain edge, and a 2 K decrease in gas temperature for r > 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines.

We present a study of molecular gas in the inner disk $\left(r < 20 \, \text{AU} \right)$ around RY Lupi, with spectra from HST-COS, HST-STIS, and VLT-CRIRES. We model the radial distribution of flux from hot gas in a surface layer between $r = 0.1-10$ AU, as traced by Ly$\alpha$-pumped H$_2$. The result shows H$_2$ emission originating in a ring centered at $\sim$3 AU that declines within $r < 0.1$ AU, which is consistent with the behavior of disks with dust cavities. An analysis of the H$_2$ line shapes shows that a two-component Gaussian profile $\left(\text{FWHM}_{broad, H_2} = 105 \pm 15 \, \text{km s}^{-1}; \, \text{FWHM}_{narrow, H_2} = 43 \pm 13 \, \text{km s}^{-1} \right)$ is statistically preferred to a single-component Gaussian. We interpret this as tentative evidence for gas emitting from radially separated disk regions $\left( \left \langle r_{broad, H_2} \right \rangle \sim 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, H_2} \right \rangle \sim 3 \pm 2 \, \text{AU} \right)$. The 4.7 $\mu$m $^{12}$CO emission lines are also well fit by two-component profiles $\left( \left \langle r_{broad, CO} \right \rangle = 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, CO} \right \rangle = 15 \pm 2 \, \text{AU} \right)$. We combine these results with 10 $\mu$m observations to form a picture of gapped structure within the mm-imaged dust cavity, providing the first such overview of the inner regions of a young disk. The HST SED of RY Lupi is available online for use in modeling efforts.

Methanol is formed via surface reactions on icy dust grains. Methanol is also detected in the gas-phase at temperatures below its thermal desorption temperature and at levels higher than can be explained by pure gas-phase chemistry. The process that controls the transition from solid state to gas-phase methanol in cold environments is not understood. The goal of this work is to investigate whether thermal CO desorption provides an indirect pathway for methanol to co-desorb at low temperatures. Mixed CH$_{3}$OH:CO/CH$_{4}$ ices were heated under UHV (ultra-high vacuum) conditions and ice contents are traced using RAIRS (reflection absorption IR spectroscopy), while desorbing species were detected mass spectrometrically. An updated gas-grain chemical network was used to test the impact of the results of these experiments. The physical model used is applicable for TW Hya, a protoplanetary disk in which cold gas-phase methanol has recently been detected. Methanol release together with thermal CO desorption is found to be an ineffective process in the experiments, resulting in an upper limit of $\leq$ 7.3 $\times$ 10$^{-7}$ CH3OH molecules per CO molecule over all ice mixtures considered. Chemical modelling based on the upper limits shows that co-desorption rates as low as 10$^{-6}$ CH$_{3}$OH molecules per CO molecule are high enough to release substantial amounts of methanol to the gas-phase at and around the location of the CO thermal desorption front in a protoplanetary disk. The impact of thermal co-desorption of CH3OH with CO as a grain-gas bridge mechanism is compared with that of UV induced photodesorption and chemisorption.

Episodic accretion has been proposed as a solution to the long-standing luminosity problem in star formation; however, the process remains poorly understood. We present observations of line emission from N2H+ and CO isotopologues using the Atacama Large Millimeter/submillimeter Array (ALMA) in the envelopes of eight Very Low Luminosity Objects (VeLLOs). In five of the sources the spatial distribution of emission from N2H+ and CO isotopologues shows a clear anti-correlation. It is proposed that this is tracing the CO snow line in the envelopes: N2H+ emission is depleted toward the center of these sources in contrast to the CO isotopologue emission which exhibits a peak. The positions of the CO snow lines traced by the N2H+ emission are located at much larger radii than those calculated using the current luminosities of the central sources. This implies that these five sources have experienced a recent accretion burst because the CO snow line would have been pushed outwards during the burst due to the increased luminosity of the central star. The N2H+ and CO isotopologue emission from DCE161, one of the other three sources, is most likely tracing a transition disk at a later evolutionary stage. Excluding DCE161, five out of seven sources (i.e., ~70%) show signatures of a recent accretion burst. This fraction is larger than that of the Class 0/I sources studied by J{\o}rgensen et al. (2015) and Frimann et al. (2016) suggesting that the interval between accretion episodes in VeLLOs is shorter than that in Class 0/I sources.

Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5$\sigma$ detection in the stacked line profile. The emission is compact but is spectrally resolved and the line profile has two components. One component peaks at the source velocity and the other is blue-shifted by 5 km s$^{-1}$. The kinematics and spatial distribution of the SO emission are not consistent with that expected from a purely Keplerian disk. We detect additional blue-shifted emission that we attribute to a disk wind. The disk component was simulated using LIME and a physical disk structure. The disk emission is asymmetric and best fit by a wedge of emission in the north east region of the disk coincident with a `hot-spot' observed in the CO $J=3-2$ line. The favoured hypothesis is that a possible inner disk warp (seen in CO emission) directly exposes the north-east side of the disk to heating by the central star, creating locally the conditions to launch a disk wind. Chemical models of a disk wind will help to elucidate why the wind is particularly highlighted in SO emission and whether a refractory source of sulphur is needed. An alternative explanation is that the SO is tracing an accretion shock from a circumplanetary disk associated with the proposed protoplanet embedded in the disk at 50 au. We also report a non-detection of SO in the protoplanetary disk around HD 97048.

We present spatially resolved ALMA images of CO J=3-2 emission from the protoplanetary disk around HD100546. We model the spatially-resolved kinematic structure of the CO emission. Assuming a velocity profile which prescribes a flat or flared emitting surface in Keplerian rotation, we uncover significant residuals with a peak of $\approx7\delta v$, where $\delta v = 0.21$ km s$^{-1}$ is the width of a spectral resolution element. The residuals reveal the possible presence of a severely warped and twisted inner disk extending to at most 100au. Adapting the model to include a misaligned inner gas disk with (i) an inclination almost edge-on to the line of sight, and (ii) a position angle almost orthogonal to that of the outer disk reduces the residuals to $< 3\delta v$. However, these findings are contrasted by recent VLT/SPHERE, MagAO/GPI, and VLTI/PIONIER observations of HD100546 that show no evidence of a severely misaligned inner dust disk down to spatial scales of $\sim 1$au. An alternative explanation for the observed kinematics are fast radial flows mediated by (proto)planets. Inclusion of a radial velocity component at close to free-fall speeds and inwards of $\approx 50$au results in residuals of $\approx 4 \delta v$. Hence, the model including a radial velocity component only does not reproduce the data as well as that including a twisted and misaligned inner gas disk. Molecular emission data at a higher spatial resolution (of order 10au) are required to further constrain the kinematics within $\lesssim 100$au. HD100546 joins several other protoplanetary disks for which high spectral resolution molecular emission shows that the gas velocity structure cannot be described by a purely Keplerian velocity profile with a universal inclination and position angle. Regardless of the process, the most likely cause is the presence of an unseen planetary companion. (Abridged)

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

Gas-phase methanol was recently detected in a protoplanetary disk for the first time with ALMA. The peak abundance and distribution of methanol observed in TW Hya differed from that predicted by chemical models. Here, the chemistry of methanol gas and ice is calculated using a physical model tailored for TW Hya with the aim to contrast the results with the recent detection in this source. New pathways for the formation of larger complex molecules (e.g., ethylene glycol) are included in an updated chemical model, as well as the fragmentation of methanol ice upon photodesorption. It is found that including fragmentation upon photodesorption improves the agreement between the peak abundance reached in the chemical models with that observed in TW Hya ($\sim 10^{-11}$ with respect to \ce{H2}); however, the model predicts that the peak in emission resides a factor of $2-3$ farther out in the disk than the ALMA images. Reasons for the persistent differences in the gas-phase methanol distribution between models and the observations of TW Hya are discussed. These include the location of the ice reservoir which may coincide with the compact mm-dust disk ($\lesssim 60$~au) and sources of gas-phase methanol which have not yet been considered in models. The possibility of detecting larger molecules with ALMA is also explored. Calculations of the rotational spectra of complex molecules other than methanol using a parametric model constrained by the TW Hya observations suggest that the detection of individual emission lines of complex molecules with ALMA remains challenging. However, the signal-to-noise ratio can be enhanced via stacking of multiple transitions which have similar upper energy levels.

(Abridged) Exoplanet atmospheres are thought be built up from accretion of gas as well as pebbles and planetesimals in the midplanes of planet-forming disks. The chemical composition of this material is usually assumed to be unchanged during the disk lifetime. However, chemistry can alter the relative abundances of molecules in this planet-building material. To assess the impact of disk chemistry during the era of planet formation, an extensive kinetic chemistry gas-grain reaction network is utilised to evolve the abundances of chemical species over time. Given a high level of ionisation, chemical evolution in protoplanetary disk midplanes becomes significant after a few times $10^{5}$ yrs, and is still ongoing by 7 Myr between the H$_{2}$O and the O$_{2}$ icelines. Importantly, the changes in the abundances of the major elemental carbon and oxygen-bearing molecules imply that the traditional "stepfunction" for the C/O ratios in gas and ice in the disk midplane (as defined by sharp changes at icelines of H$_{2}$O, CO$_{2}$ and CO) evolves over time, and cannot be assumed fixed. In addition, at lower temperatures (< 29 K), gaseous CO colliding with the grains gets converted into CO$_{2}$ and other more complex ices, lowering the CO gas abundance between the O$_{2}$ and CO thermal icelines. This effect can mimic a CO iceline at a higher temperature than suggested by its binding energy. Chemistry in the disk midplane is ionisation-driven, and evolves over time. In order to reliably predict the atmospheric compositions of forming planets, as well as to relate observed atmospheric C/O ratios of exoplanets to where and how the atmospheres have formed in a disk midplane, chemical evolution needs to be considered and implemented into planet formation models.

Advanced telescopes, such as ALMA and JWST, are likely to show that the chemical universe may be even more complex than currently observed, requiring astrochemical modelers to improve their models to account for the impact of new data. However, essential input information for gas-grain models, such as binding energies of molecules to the surface, have been derived experimentally only for a handful of species, leaving hundreds of species with highly uncertain estimates. We present in this paper a systematic study of the effect of uncertainties in the binding energies on an astrochemical two-phase model of a dark molecular cloud, using the rate equations approach. A list of recommended binding energy values based on a literature search of published data is presented. Thousands of simulations of dark cloud models were run, and in each simulation a value for the binding energy of hundreds of species was randomly chosen from a normal distribution. Our results show that the binding energy of H$_{2}$ is critical for the surface chemistry. For high binding energy, H$_{2}$ freezes out on the grain forming an H$_{2}$ ice. This is not physically realistic and we suggest a change in the rate equations. The abundance ranges found are in reasonable agreement with astronomical ice observations. Pearson correlation coefficients revealed that the binding energy of HCO, HNO, CH$_{2}$, and C correlate most strongly with the abundance of dominant ice species. Finally, the formation route of complex organic molecules was found to be sensitive to the branching ratios of H$_{2}$CO hydrogenation.

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

We report on subarcsecond observations of complex organic molecules (COMs) in the high-mass protostar IRAS20126+4104 with the Plateau de Bure Interferometer in its most extended configurations. In addition to the simple molecules SO, HNCO and H2-13CO, we detect emission from CH3CN, CH3OH, HCOOH, HCOOCH3, CH3OCH3, CH3CH2CN, CH3COCH3, NH2CN, and (CH2OH)2. SO and HNCO present a X-shaped morphology consistent with tracing the outflow cavity walls. Most of the COMs have their peak emission at the putative position of the protostar, but also show an extension towards the south(east), coinciding with an H2 knot from the jet at about 800-1000 au from the protostar. This is especially clear in the case of H2-13CO and CH3OCH3. We fitted the spectra at representative positions for the disc and the outflow, and found that the abundances of most COMs are comparable at both positions, suggesting that COMs are enhanced in shocks as a result of the passage of the outflow. By coupling a parametric shock model to a large gas-grain chemical network including COMs, we find that the observed COMs should survive in the gas phase for about 2000 yr, comparable to the shock lifetime estimated from the water masers at the outflow position. Overall, our data indicate that COMs in IRAS20126+4104 may arise not only from the disc, but also from dense and hot regions associated with the outflow.

Observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

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

[Abridged] Snowlines in protoplanetary disks play an important role in planet formation and composition. Since the CO snowline is difficult to observe directly with CO emission, its location has been inferred in several disks from spatially resolved ALMA observations of DCO+ and N2H+. N2H+ is considered to be a good tracer of the CO snowline based on astrochemical considerations predicting an anti-correlation between N2H+ and gas-phase CO. In this work, the robustness of N2H+ as a tracer of the CO snowline is investigated. A simple chemical network is used in combination with the radiative transfer code LIME to model the N2H+ distribution and corresponding emission in the disk around TW Hya. The assumed CO and N2 abundances, corresponding binding energies, cosmic ray ionization rate, and degree of large-grain settling are varied to determine the effects on the N2H+ emission and its relation to the CO snowline. For the adopted physical structure of the TW Hya disk and molecular binding energies for pure ices, the balance between freeze-out and thermal desorption predicts a CO snowline at 19 AU, corresponding to a CO midplane freeze-out temperature of 20 K. A model with a total, i.e. gas plus ice, CO abundance of 3e-6 with respect to H2 fits the position of the emission peak observed by Qi et al. 2013 for the TW Hya disk. However, the relationship between N2H+ and the CO snowline is more complicated than generally assumed: for the investigated parameters, the N2H+ column density peaks at least 5 AU outside the CO snowline. Moreover, the N2H+ emission can peak much further out, as far as ~50 AU beyond the snowline. Hence, chemical modeling, as done here, is necessary to derive a CO snowline location from N2H+ observations.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also marginally resolved. The most remarkable finding is that the spectral index $\alpha (R)$ between bands 4 and 6 peaks at the 22~au gap. The derived power-law index of the dust opacity $\beta (R)$ is $\sim1.7$ at the 22~au gap and decreases toward the disk center to $\sim0$. The most prominent gap at 22~au could be caused by the gravitational interaction between the disk and an unseen planet with a mass of $\lesssim$1.5 $M_\mathrm{Neptune}$ although other origins may be possible. The planet-induced gap is supported by the fact that $\beta (R)$ is enhanced at the 22~au gap, indicating a deficit of $\sim$mm-sized grains within the gap due to dust filtration by a planet.

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

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

[Abridged] Chemical evolution in the protoplanetary disk midplane can modify the composition of ices and gases. We have investigated if and how chemical evolution affects the abundances and distributions of key volatile species in the midplane of a protoplanetary disk in the 0.2-30 AU range. A full chemical network including gas-phase, gas-grain interactions and grain-surface chemistry is used to evolve chemistry in time, for 1 Myr. Great diversity is observed in the relative abundance ratios of the main considered species: H2O, CO, CO2, CH4, O2, NH3 and N2. The choice of ionisation level, the choice of initial abundances, as well as the extent of chemical reaction types included are all factors that affect the chemical evolution. The only exception is the inheritance scenario with a low ionisation level, which results in negligible changes compared with the initial abundances, regardless of whether grain-surface chemistry is included. The chemical processing changes the C/O ratios for gas and ice significantly with radius and with model assumptions. For high ionisation levels, chemical processing becomes significant after a few times $10^{5}$~yrs. It is concluded that chemical evolution in the disk midplane needs to be considered in the determination of the volatile composition of planetesimals.

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.

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

Comets provide a unique insight into the molecular composition and complexity of the material in the primordial solar nebula. Recent results from the Rosetta mission, currently monitoring comet 67P/Churyumov-Gerasimenko in situ, and ALMA (the Atacama Large Millimeter/submillimeter Array), have demonstrated a tantalising link between the chemical complexity now confirmed in disks (via the detection of gas-phase CH3CN; Oberg et al. 2015) and that confirmed on the surface of 67P (Goesmann et al. 2015), raising questions concerning the chemical origin of such species (cloud or inheritance versus disk synthesis). Results from an astrochemical model of a protoplanetary disk are presented in which complex chemistry is included and in which it is assumed that simple ices only are inherited from the parent molecular cloud. The model results show good agreement with the abundances of several COMs observed on the surface of 67P with Philae/COSAC. Cosmic-ray and X-ray-induced photoprocessing of predominantly simple ices inherited by the protoplanetary disk is sufficient to generate a chemical complexity similar to that observed in comets. This indicates that the icy COMs detected on the surface of 67P may have a disk origin. The results also show that gas-phase CH3CN is abundant in the inner warm disk atmosphere where hot gas-phase chemistry dominates and potentially erases the ice chemical signature. Hence, CH3CN may not be an unambiguous tracer of the complex organic ice reservoir. However, a better understanding of the hot gas-phase chemistry of CH3CN is needed to confirm this preliminary conclusion.

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

(Abridged) Near- to mid-IR observations of protoplanetary disks show that the inner regions (<10AU) are rich in small organic volatiles (e.g., C2H2 and HCN). Trends in the data suggest that disks around cooler stars (~3000K) are potentially more carbon- and molecule-rich than their hotter counterparts. Our aims are to explore the composition of the planet-forming region of disks around stars from M dwarf to Herbig Ae and compare with the observed trends. Models of the disk physical structure are coupled with a gas-grain chemical network to map the abundances in the planet-forming zone. N2 self shielding, X-ray-induced chemistry, and initial abundances, are investigated. The composition in the 'observable' atmosphere is compared with that in the midplane where the planet-building reservoir resides. M dwarf disk atmospheres are relatively more molecule rich than those for T Tauri or Herbig Ae disks. The weak far-UV flux helps retain this complexity which is enhanced by X-ray-induced ion-molecule chemistry. N2 self shielding has only a small effect and does not explain the higher C2H2/HCN ratios observed towards cooler stars. The models underproduce the OH/H2O column density ratios constrained in Herbig Ae disks, despite reproducing the absolute value for OH: H2O self shielding only increases this discrepency. The disk midplane content is sensitive to the initial main elemental reservoirs. The gas in the inner disk is generally more carbon rich than the midplane ices and is most significant for disks around cooler stars. The atmospheric C/O ratio appears larger than it actually is when calculated using observable tracers only because gas-phase O2 is predicted to be a significant oxygen reservoir. The models suggest that the gas in the inner regions of disks around cooler stars is more carbon rich; however, calculations of the molecular emission are necessary to confirm the observed trends.

Context. The disk around the Herbig Ae star, HD 100546, shows structures that suggest the presence of two companions in the disk at $\sim10$ and $\sim70$ AU. The outer companion seems to be in the act of formation. Aims. Our aims are to provide constraints on the age of the planets in HD 100546 and to explore the potential evidence for sequential planet formation in transition disks such as HD 100546. Methods. We compare the recent resolved continuum observations of the disk around HD 100546 with the results of dust evolution simulations using an analytical prescription for the shapes of gaps carved by massive planets. Results. An inner pressure bump must have been present since early in the disk lifetime to have good agreement between the dust evolution models and the continuum observations of HD 100546. This pressure bump may have resulted from the presence of a very massive planet ($\sim20 M_{\rm{Jup}}$), which formed early in the inner disk ($r\sim$10 AU). If only this single planet exists, the disk is likely to be old, comparable to the stellar age ($\sim$5-10 Myr). Another possible explanation is an additional massive planet in the outer disk ($r\sim70$ AU): either a low-mass outer planet ($\lesssim5 M_{\rm{Jup}}$) injected at early times, or a higher mass outer planet ($\gtrsim15 M_{\rm{Jup}}$) formed very recently, traps the right amount of dust in pressure bumps to reproduce the observations. In the latter case, the disk could be much younger ($\sim3.0$ Myr). Conclusions. In the case in which two massive companions are embedded in the disk around HD 100546, as suggested in the literature, the outer companion could be at least $\gtrsim$2.5 Myr younger than the inner companion.

There is much evidence that planet formation is occurring in the disk around the Herbig Be star HD100546. To learn more about the processes occurring in this disk we conducted high resolution imaging at 43/45 GHz with the Australia Telescope Compact Array (ATCA). Multiple array configurations were used, providing a best spatial resolution of $\sim$ 0.15 arcsec, or 15 AU at HD100546's distance of $\sim$ 100 pc. Significant structure is revealed, but its precise form is dependent on the $u-v$ plane sampling used for the image reconstruction. At a resolution of $\leq$ 30 AU we detected an inner gap in the disk with a radius of $\sim$ 25 AU and a position angle approximately along the known disk major axis. With different weighting, and an achieved resolution of $\sim$ 15 AU, emission appears at the centre and the disk takes on the shape of an incomplete ring, much like a horseshoe, again with a gap radius of $\sim$ 25 AU. The position angle of the disk major axis and its inclination from face-on are determined to be $140^{\circ}\pm5^{\circ}$ and $40^{\circ}\pm5^{\circ}$ respectively. The $\sim$ 25 AU gap radius is confirmed by a null in the real part of the binned visibilities at 320$\pm$10 k$\lambda$, whilst the non-axisymmetric nature is also confirmed through significant structure in the imaginary component. The emission mechanism at the central peak is most likely to be free-free emission from a stellar or disk wind. Overall our data support the picture of at least one, but probably several, giant planets orbiting HD100546 within 25 AU.

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent observations suggest that their emission stems, not only from the traditional hot corino, but also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity system is carried out. Wavelength-dependent radiative transfer calculations are performed and a comprehensive gas-grain chemical network is used to simulate the physical and chemical structure. The morphology of the system delineates three distinct regions: the cavity wall layer with time-dependent and species-variant enhancements; a torus rich in complex organic ices, but not reflected in gas-phase abundances; and the remaining outer envelope abundant in simpler solid and gaseous molecules. Strongly irradiated regions, such as the cavity wall layer, are subject to frequent photodissociation in the solid phase. Subsequent recombination of the photoproducts leads to frequent reactive desorption, causing gas-phase enhancements of several orders of magnitude. This mechanism remains to be quantified with laboratory experiments. Direct photodesorption is found to be relatively inefficient. If radicals are not produced directly in the icy mantle, the formation of complex organics is impeded. For efficiency, a sufficient number of FUV photons needs to penetrate the envelope; and elevated cool dust temperatures need to enable grain-surface radical mobility. As a result, a high stellar luminosity and a sufficiently wide cavity favor chemical complexity. Furthermore within this paradigm, complex organics are demonstrated to have unique lifetimes and be grouped into early (formaldehyde, ketene, methanol, formic acid, methyl formate, acetic acid, glycolaldehyde) and late (acetaldehyde, dimethyl ether, ethanol) species.

The HDO/H2O ratio in interstellar gas is often used to draw conclusions on the origin of water in star-forming regions and on Earth. In cold cores and in the outer regions of protoplanetary disks, gas-phase water comes from photodesorption of water ice. We present fitting formulae for implementation in astrochemical models using photodesorption efficiencies for all water ice isotopologues obtained using classical molecular dynamics simulations. We investigate if the gas-phase HDO/H2O ratio reflects that present in the ice or whether fractionation can occur during photodesorption. Probabilities for the top four monolayers are presented for photodesorption of X (X=H,D) atoms, OX radicals, and X2O and HDO molecules following photodissociation of H2O, D2O, and HDO in H2O amorphous ice at temperatures from 10-100 K. Isotope effects are found for all products: (1) H atom photodesorption probabilities from H2O ice are larger than those for D atom photodesorption from D2O ice by a factor of 1.1; the ratio of H and D photodesorbed upon HDO photodissociation is a factor of 2. This process will enrich the ice in deuterium atoms over time; (2) the OD/OH photodesorption ratio upon D2O and H2O photodissociation is on average a factor of 2, but the ratio upon HDO photodissociation is almost constant at unity for all temperatures; (3) D atoms are more effective in kicking out neighbouring water molecules than H atoms. However, the ratio of the photodesorbed HDO and H2O molecules is equal to the HDO/H2O ratio in the ice, therefore, there is no isotope fractionation upon HDO and H2O photodesorption. Nevertheless, the enrichment of the ice in D atoms due to photodesorption can over time lead to an enhanced HDO/H2O ratio in the ice, and, when photodesorbed, also in the gas. The extent to which the ortho/para ratio of H2O can be modified by the photodesorption process is also discussed. (Abridged)

Aims: Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods: A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., FUV field, initial abundances, timescales, and kinematics). Results: The observed water ice abundances with respect to hydrogen nuclei in our sample are 30-80ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T(dust)>15 K. This requires short prestellar core lifetimes of less than about 0.1Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions: The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short pre-collapse time of less than 0.1 Myr at densities of n(H)~1e4 cm-3, or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.

Interstellar methanol is considered to be a parent species of larger, more complex organic molecules. A physicochemical simulation of infalling parcels of matter is performed for a low-mass star-forming system to trace the chemical evolution from cloud to disc. An axisymmetric 2D semi-analytic model generates the time-dependent density and velocity distributions, and full continuum radiative transfer is performed to calculate the dust temperature and the UV radiation field at each position as a function of time. A comprehensive gas-grain chemical network is employed to compute the chemical abundances along infall trajectories. Two physical scenarios are studied, one in which the dominant disc growth mechanism is viscous spreading, and another in which continuous infall of matter prevails. The results show that the infall path influences the abundance of methanol entering each type of disc, ranging from complete loss of methanol to an enhancement by a factor of > 1 relative to the prestellar phase. Critical chemical processes and parameters for the methanol chemistry under different physical conditions are identified. The exact abundance and distribution of methanol is important for the budget of complex organic molecules in discs, which will be incorporated into forming planetary system objects such as protoplanets and comets. These simulations show that the comet-forming zone contains less methanol than in the precollapse phase, which is dominantly of prestellar origin, but also with additional layers built up in the envelope during infall. Such intriguing links will soon be tested by upcoming data from the Rosetta mission.

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.

The envelopes of AGB stars are irradiated externally by ultraviolet photons; hence, the chemistry is sensitive to the photodissociation of N$_2$ and CO, which are major reservoirs of nitrogen and carbon, respectively. The photodissociation of N$_2$ has recently been quantified by laboratory and theoretical studies. Improvements have also been made for CO photodissociation. For the first time, we use accurate N$_2$ and CO photodissociation rates and shielding functions in a model of the circumstellar envelope of the carbon-rich AGB star, IRC +10216. We use a state-of-the-art chemical model of an AGB envelope, the latest CO and N$_2$ photodissociation data, and a new method for implementing molecular shielding functions in full spherical geometry with isotropic incident radiation. We compare computed column densities and radial distributions of molecules with observations. The transition of N$_2$ $\to$ N (also, CO $\to$ C $\to$ C$^+$) is shifted towards the outer envelope relative to previous models. This leads to different column densities and radial distributions of N-bearing species, especially those species whose formation/destruction processes largely depend on the availability of atomic or molecular nitrogen, for example, C$_n$N ($n$=1, 3, 5), C$_n$N$^-$ ($n$=1, 3, 5), HC$_n$N ($n$=1, 3, 5, 7, 9), H$_2$CN and CH$_2$CN. The chemistry of many species is directly or indirectly affected by the photodissociation of N$_2$ and CO, especially in the outer shell of AGB stars where photodissociation is important. Thus, it is important to include N$_2$ and CO shielding in astrochemical models of AGB envelopes and other irradiated environments. In general, while differences remain between our model of IRC +10216 and the observed molecular column densities, better agreement is found between the calculated and observed radii of peak abundance.

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

(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.

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (~400 pc).

We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we also now include, for download, state-specific deuterated rate coefficients, deuterium exchange reactions and a list of surface binding energies for many neutral species. Where possible, we have referenced the original source of all new and existing data. We have tested the main reaction network using a dark cloud model and a carbon-rich circumstellar envelope model. We present and briefly discuss the results of these models.

Prompted by the ongoing interest in Spitzer Infrared Spectrometer spectra of carbon stars in the Large Magellanic Cloud, we have investigated the circumstellar chemistry of carbon stars in low-metallicity environments. Consistent with observations, our models show that acetylene is particularly abundant in the inner regions of low metallicity carbon-rich AGB stars -- more abundant than carbon monoxide. As a consequence, larger hydrocarbons have higher abundances at the metallicities of the Magellanic Clouds than in stars with solar metallicity. We also find the oxygen and nitrogen chemistry is suppressed at lower metallicity, as expected. Finally, we calculate molecular line emission from carbon stars in the Large and Small Magellanic Cloud and find that several molecules should be readily detectable with the Atacama Large Millimeter Array at Full Science operations.

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.

We observed the 9-cm Lambda-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0e-8, to the southeast where it reaches a value of 2.0e-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of "early-type" molecules in this region.

We calculate the chemical evolution of protoplanetary disks considering radial viscous accretion, vertical turbulent mixing and vertical disk winds. We study the effects on the disk chemical structure when different models for the formation of molecular hydrogen on dust grains are adopted. Our gas-phase chemistry is extracted from the UMIST Database for Astrochemistry (Rate06) to which we have added detailed gas-grain interactions. We use our chemical model results to generate synthetic near- and mid-infrared LTE line emission spectra and compare these with recent Spitzer observations. Our results show that if H2 formation on warm grains is taken into consideration, the H2O and OH abundances in the disk surface increase significantly. We find the radial accretion flow strongly influences the molecular abundances, with those in the cold midplane layers particularly affected. On the other hand, we show that diffusive turbulent mixing affects the disk chemistry in the warm molecular layers, influencing the line emission from the disk and subsequently improving agreement with observations. We find that NH3, CH3OH, C2H2 and sulphur-containing species are greatly enhanced by the inclusion of turbulent mixing. We demonstrate that disk winds potentially affect the disk chemistry and the resulting molecular line emission in a similar manner to that found when mixing is included.

We report the detection of microwave emission lines from the hydrocarbon anion C6H- and its parent neutral C6H in the star-forming region L1251A (in Cepheus), and the pre-stellar core L1512 (in Auriga). The carbon-chain-bearing species C4H, HC3N, HC5N, HC7N and C3S are also detected in large abundances. The observations of L1251A constitute the first detections of anions and long-chain polyynes and cyanopolyynes (with more than 5 carbon atoms) in the Cepheus Flare star-forming region, and the first detection of anions in the vicinity of a protostar outside of the Taurus molecular cloud complex, highlighting a wider importance for anions in the chemistry of star formation. Rotational excitation temperatures have been derived from the HC3N hyperfine structure lines, and are found to be 6.2 K for L1251A and 8.7 K for L1512. The anion-to-neutral ratios are 3.6% and 4.1%, respectively, which are within the range of values previously observed in the interstellar medium, and suggest a relative uniformity in the processes governing anion abundances in different dense interstellar clouds. This research contributes towards the growing body of evidence that carbon chain anions are relatively abundant in interstellar clouds throughout the Galaxy, but especially in the regions of relatively high density and high depletion surrounding pre-stellar cores and young, embedded protostars.

We have developed a high resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (sub-milli-arcsecond in the inner disk for objects at the distance of Taurus, ~ 140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas-grain interactions including freeze out and thermal and non-thermal desorption (cosmic-ray induced desorption, photodesorption and X-ray desorption) and a grain-surface network. We find that cosmic-ray induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenise the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3 and CH3OCH3 to potentially observable values (i.e. a fractional abundance of >~ 1.0E-11).

We have investigated the role of molecular anion chemistry in pseudo-time dependent chemical models of dark clouds. With oxygen-rich elemental abundances, the addition of anions results in a slight improvement in the overall agreement between model results and observations of molecular abundances in TMC-1 (CP). More importantly, with the inclusion of anions, we see an enhanced production efficiency of unsaturated carbon-chain neutral molecules, especially in the longer members of the families CnH, CnH2, and HCnN. The use of carbon-rich elemental abundances in models of TMC-1 (CP) with anion chemistry worsens the agreement with observations obtained in the absence of anions.

The recent detection of the hydrocarbon anion C6H- in the interstellar medium has led us to investigate the synthesis of hydrocarbon anions in a variety of interstellar and circumstellar environments. We find that the anion/neutral abundance ratio can be quite large, on the order of at least a few percent, once the neutral has more than five carbon atoms. Detailed modeling shows that the column densities of C6H- observed in IRC+10216 and TMC-1 can be reproduced. Our calculations also predict that other hydrocarbon anions, such as C4H- and C8H-, are viable candidates for detection in IRC+10216, TMC-1 and photon-dominated regions such as the Horsehead Nebula.