We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($\alpha_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($\tau\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11{\deg}. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

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

Context. The relationship between outflow launching and formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) has carried out long-baseline observations towards a handful of sources, revealing connections between outflows and the inner regions of disks. Aims. Here we aim to determine the small-scale kinematic and morphological properties of the outflow from the isolated protostar B335 for which no Keplerian disk has, so far, been observed on scales down to 10 au. Methods. We use ALMA in its longest-baseline configuration to observe emission from CO isotopologs, SiO, SO$_2$ and CH$_3$OH. The proximity of B335 provides a resolution of ~3 au (0.03''). We also combine our long-baseline data with archival data to produce a high-fidelity image covering scales up to 700 au (7''). Results. $^{12}$CO has a X-shaped morphology with arms ~50 au in width that we associate with the walls of an outflow cavity, similar to what is observed on larger scales. Long-baseline continuum emission is confined to <7 au of the protostar, while short-baseline continuum emission follows the $^{12}$CO outflow and cavity walls. Methanol is detected within ~30 au of the protostar. SiO is also detected in the vicinity of the protostar, but extended along the outflow. Conclusions. The $^{12}$CO outflow shows no clear signs of rotation at distances $\gtrsim$30 au from the protostar. SiO traces the protostellar jet on small scales, but without obvious rotation. CH$_3$OH and SO$_2$ trace a region <16 au in diameter, centred on the continuum peak, which is clearly rotating. Using episodic, high-velocity, $^{12}$CO features, we estimate the launching radius of the outflow to be <0.1 au and dynamical timescales on the order of a few years.

The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question.

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

(Abridged) The purpose of this paper is to explore and compare the physical and chemical structure of Class I low-mass protostellar sources on protoplanetary disc scales. We present a study of the dust and gas emission towards a representative sample of 12 Class I protostars from the Ophiuchus molecular cloud with the Atacama Large Millimeter/submillimeter Array (ALMA). The continuum at 0.87 mm and molecular transitions from C17O, C34S, H13CO+, CH3OH, SO2 , and C2H were observed at high angular resolution (0.4", ~60 au diameter) towards each source. Disc and stellar masses are estimated from the continuum flux and position-velocity diagrams, and six of the sources show disc-like structures. Towards the more luminous sources, compact emission and large line widths are seen for transitions of SO2 that probe warm gas (Eu ~200 K). In contrast, C17O emission is detected towards the least evolved and less luminous systems. No emission of CH3OH is detected towards any of the continuum peaks, indicating an absence of warm CH3OH gas towards these sources. A power-law relation is seen between the stellar mass and the bolometric luminosity, corresponding to a mass accretion rate of (2.4 +/- 0.6) x 10^-7 Msun/year for the Class I sources. This mass accretion rate is lower than the expected value if the accretion is constant in time and rather points to a scenario of accretion occurring in bursts. The differentiation between C17O and SO2 suggests that they trace different physical components: C17O traces the densest and colder regions of the disc-envelope system, while SO2 may be associated with regions of higher temperature, such as accretion shocks. The lack of warm CH3OH emission suggests that there is no hot-core-like region around any of the sources and that the CH3OH column density averaged over the disc is low.

We report new Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 3 observations at 2.75 mm of the TD around SR 24S with an angular resolution of $\sim$0.11''$\times$ 0.09'' and a peak signal-to-noise ratio of $\sim24$. We detect an inner disk and a mostly symmetric ring-like structure that peaks at $\sim$0.32'', that is $\sim$37 au at a distance of $\sim$114.4 pc. The full width at half maximum of this ring is $\sim$28 au. We analyze the observed structures by fitting the dust continuum visibilities using different models for the intensity profile, and compare with previous ALMA observations of the same disk at 0.45 mm and 1.30 mm. We qualitatively compare the results of these fits with theoretical predictions of different scenarios for the formation of a cavity or large gap. The comparison of the dust continuum structure between different ALMA bands indicates that photoevaporation and dead zone can be excluded as leading mechanisms for the cavity formation in SR 24S disk, leaving the planet scenario (single or multiple planets) as the most plausible mechanism. We compared the 2.75 mm emission with published (sub-)centimeter data and find that the inner disk is likely tracing dust thermal emission. This implies that any companion in the system should allow dust to move inwards throughout the gap and replenish the inner disk. In the case of one single planet, this puts strong constraints on the mass of the potential planet inside the cavity and the disk viscosity of about $\lesssim$5 $M_{\rm{Jup}}$ and $\alpha\sim10^{-4}-10^{-3}$, respectively.

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters.

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet.

Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

We report the detection of methanol in the disk around the young outbursting star V883 Ori with the Atacama Large Millimeter/submillimeter Array (ALMA). Four transitions are observed with upper level energies ranging between 115 and 459 K. The emission is spatially resolved with the 0.14" beam and follows the Keplerian rotation previously observed for C$^{18}$O. Using a rotational diagram analysis, we find a disk-averaged column density of $\sim10^{17}$ cm$^{-2}$ and a rotational temperature of $\sim90-100$ K, suggesting that the methanol has thermally desorbed from the dust grains. We derive outer radii between 120 and 140 AU for the different transitions, compared to the 360 AU outer radius for C$^{18}$O. Depending on the exact physical structure of the disk, the methanol emission could originate in the surface layers beyond the water snowline. Alternatively, the bulk of the methanol emission originates inside the water snowline, which can then be as far out as ~100 AU, instead of 42 AU as was previously inferred from the continuum opacity. In addition, these results show that outbursting young stars like V883 Ori are good sources to study the ice composition of planet forming material through thermally desorbed complex molecules, which have proven to be hard to observe in more evolved protoplanetary disks.

VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $\mu$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $\mu$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $\lambda~<$ 450 $\mu$m. These results, however, do not explain the lack of circumstellar gas emission and puzzling nature of VLA 1623$-$2417 B.

The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars.

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.

Much attention has been placed on the dust distribution in protostellar envelopes, but there are still many unanswered questions regarding the structure of the gas. We aim to start identifying the factors that determine the chemical structure of protostellar regions, by studying and comparing low-mass embedded systems in key molecular tracers. The cold and warm chemical structures of two embedded Class 0 systems, IRAS16293 and VLA1623 are characterized through interferometric observations. DCO+, N2H+ and N2D+ are used to trace the spatial distribution and physics of the cold regions of the envelope, while c-C3H2 and C2H from models of the chemistry are expected to trace the warm (UV-irradiated) regions. Both sources show a number of striking similarities and differences. DCO+ consistently traces the cold material at the disk-envelope interface, where gas and dust temperatures are lowered due to disk shadowing. N2H+ and N2D+, also tracing cold gas, show low abundances towards VLA1623, but for IRAS16293, the distribution of N2D+ is consistent with the same chemical models that reproduce DCO+. c-C3H2 and C2H show different spatial distributions for the two systems. For IRAS16293, c-C3H2 traces the outflow cavity wall, while C2H is found in the envelope material but not the outflow cavity wall. In contrast, toward VLA1623 both molecules trace the outflow cavity wall. Finally, hot core molecules are abundantly observed toward IRAS16293 but not toward VLA1623. We identify temperature as one of the key factors in determining the chemical structure of protostars as seen in gaseous molecules. More luminous protostars, such as IRAS16293, will have chemical complexity out to larger distances than colder protostars, such as VLA1623. Additionally, disks in the embedded phase have a crucial role in controlling both the gas and dust temperature of the envelope, and consequently the chemical structure.

A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $\alpha_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

[Abridged] Protoplanetary disks have been studied extensively, both physically and chemically, to understand the environment in which planets form. However, the first steps of planet formation are likely to occur already when the protostar and disk are still embedded in their natal envelope. The initial conditions for planet formation may thus be provided by these young embedded disks, of which the physical and chemical structure is poorly characterized. We aim to constrain the midplane temperature structure, one of the critical unknowns, of the embedded disk around L1527. In particular, we set out to determine whether there is an extended cold outer region where CO is frozen out, as is the case for Class II disks. We use archival ALMA data to directly observe the midplane of the near edge-on L1527 disk. Optically thick $^{13}$CO ($J=2-1$) and C$^{18}$O ($J=2-1$) emission is observed throughout the disk and inner envelope, while N$_2$D$^+ (J=3-2$), which can only be abundant when CO is frozen out, is not detected. Both CO isotopologues have brightness temperatures $\gtrsim$ 25 K along the midplane. Disk and envelope emission can be disentangled kinematically, because the largest velocities are reached in the disk. A power law radial temperature profile constructed using the highest midplane temperature at these velocities suggest that the temperature is above 20 K out to at least 75 AU, and possibly throughout the entire 125 AU disk. Radiative transfer models show that a model without CO freeze-out in the disk matches the C$^{18}$O observations better than a model with the CO snowline at $\sim$70 AU. In addition, there is no evidence for a large (order of magnitude) depletion of CO. The disk around L1527 is likely to be warm enough to have CO present in the gas phase throughout the disk, suggesting that young embedded disks can indeed be warmer than the more evolved Class II disks.

(Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations.

Snowlines are key ingredients for planet formation. Providing observational constraints on the locations of the major snowlines is therefore crucial for fully connecting planet compositions to their formation mechanism. Unfortunately, the most important snowline, that of water, is very difficult to observe directly in protoplanetary disks due to its close proximity to the central star. Based on chemical considerations, HCO$^+$ is predicted to be a good chemical tracer of the water snowline, because it is particularly abundant in dense clouds when water is frozen out. This work maps the optically thin isotopologue H$^{13}$CO$^+$ ($J=3-2$) toward the envelope of the low-mass protostar NGC1333-IRAS2A (observed with NOEMA at ~0.9" resolution), where the snowline is at larger distance from the star than in disks. The H$^{13}$CO$^+$ emission peaks ~2" northeast of the continuum peak, whereas the previously observed H$_2^{18}$O shows compact emission on source. Quantitative modeling shows that a decrease in H$^{13}$CO$^+$ abundance by at least a factor of six is needed in the inner ~360 AU to reproduce the observed emission profile. Chemical modeling predicts indeed a steep increase in HCO$^+$ just outside the water snowline; the 50% decrease in gaseous H$_2$O at the snowline is not enough to allow HCO$^+$ to be abundant. This places the water snowline at 225 AU, further away from the star than expected based on the 1D envelope temperature structure for NGC1333-IRAS2A. In contrast, DCO$^+$ observations show that the CO snowline is at the expected location, making an outburst scenario unlikely. The spatial anticorrelation of the H$^{13}$CO$^+$ and H$_2^{18}$O emission provide a proof of concept that H$^{13}$CO$^+$ can be used as a tracer of the water snowline.

Planets form in disks around young stars. The planet formation process may start when the protostar and disk are still deeply embedded within their infalling envelope. However, unlike more evolved protoplanetary disks, the physical and chemical structure of these young embedded disks are still poorly constrained. We have analyzed ALMA data for $^{13}$CO, C$^{18}$O and N$_2$D$^+$ to constrain the temperature structure, one of the critical unknowns, in the disk around L1527. The spatial distribution of $^{13}$CO and C$^{18}$O, together with the kinetic temperature derived from the optically thick $^{13}$CO emission and the non-detection of N$_2$D$^+$, suggest that this disk is warm enough ($\gtrsim$ 20 K) to prevent CO freeze-out.

(Abridged) Through spectrally unresolved observations of high-J CO transitions, Herschel-PACS has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components using spectrally resolved Herschel-HIFI data. Observations are presented of the highly excited CO line J=16-15 with Herschel-HIFI toward 24 low-mass protostellar objects. The spectrally resolved profiles show two distinct velocity components: a broad component with an average FWHM of 20 km/s, and a narrower component with a FWHM of 5 km/s that is often offset from the source velocity. The average rotational temperature over the entire profile, as measured from comparison between CO J=16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2x10^-6. Two distinct velocity profiles observed in the HIFI line profiles suggest that the CO ladder observed with PACS consists of two excitation components. The warm component (300 K) is associated with the broad HIFI component, and the hot component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The ubiquity of the warm and hot CO components suggests that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects.

[Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources.

Young stars are associated with prominent outflows of molecular gas. The ejection of gas via these outflows is believed to remove angular momentum from the protostellar system, thus permitting young stars to grow by accretion of material from the protostellar disk. The underlying mechanism for outflow ejection is not yet understood, but is believed to be closely linked to the protostellar disk. Assorted scenarios have been proposed to explain protostellar outflows; the main difference between these models is the region where acceleration of material takes place: close to the protostar itself ('X-wind', or stellar wind), in a larger region throughout the protostellar disk (disk wind), or at the interface between. Because of the limits of observational studies, outflow launching regions have so far only been probed by indirect extrapolation. Here we report observations of carbon monoxide toward the outflow associated with the TMC1A protostellar system. These data show that gas is ejected from a region extending up to a radial distance of 25 astronomical units from the central protostar, and that angular momentum is removed from an extended region of the disk. This demonstrates that the outflowing gas is launched by an extended disk wind from a Keplerian disk. Hence, we rule out X-wind and stellar wind launching scenarios as the source of the emission on the scales we observe.

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.

We present results from a deep sub-millimeter survey in the Serpens Main and Serpens/G3-G6 clusters, conducted with the Submillimetre Common-User Bolometer Array (SCUBA-2) at the James Clerk Maxwell Telescope. We have combined Herschel PACS far-infrared photometry, sub-millimeter continuum and molecular gas line observations, with the aim to conduct a detailed multi-wavelength characterization of `proto-brown dwarf' candidates in Serpens. We have performed continuum and line radiative transfer modeling, and have considered various classification schemes to understand the structure and the evolutionary stage of the system. We have identified four proto-brown dwarf candidates, of which the lowest luminosity source has an Lbol ~0.05 Lsun. Two of these candidates show characteristics consistent with Stage 0/I systems, while the other two are Stage I-T/Class Flat systems with tenuous envelopes. Our work has also revealed a ~20% fraction of mis-identified Class 0/I/Flat sources that show characteristics consistent with Class II edge-on disk systems. We have set constraints on the mass of the central object using the measured bolometric luminosities and numerical simulations of stellar evolution. Considering the available gas+dust mass reservoir and the current mass of the central source, three of these candidates are likely to evolve into brown dwarfs.

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

Solving the mystery of the origin of chondrules is one of the most elusive goals in the field of meteoritics. Recently the idea of planet(esimal) collisions releasing splashes of lava droplets, long considered out of favor, has been reconsidered as a possible origin of chondrules by several papers. One of the main problems with this idea is the lack of quantitative and simple models that can be used to test this scenario by directly comparing to the many known observables of chondrules. In Paper I of this series we presented a simple thermal evolution model of a spherically symmetric expanding cloud of molten lava droplets that is assumed to emerge from a collision between two planetesimals. In the present paper, number II of this series, we use this model to calculate whether or not volatile elements such as Na and K will remain abundant in these droplets or whether they will get depleted due to evaporation. The high density of the droplet cloud (e.g. small distance between adjacent droplets) causes the vapor to quickly reach saturation pressure and thus shutting down further evaporation. We show to which extent, and under which conditions, this keeps the abundances of these elements high, as is seen in chondrules. We find that for most parameters of our model (cloud mass, expansion velocity, initial temperature) the volatile elements Mg, Si and Fe remain entirely in the chondrules. The Na and K abundances inside the droplets will initially stay mostly at their initial values due to the saturation of the vapor pressure, but at some point start to drop due to the cloud expansion. However, as soon as the temperature starts to decrease, most or all of the vapor recondenses again. At the end the Na and K elements retain most of their initial abundances, albeit occasionally somewhat reduced, depending on the parameters of the expanding cloud model.

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) The physical structure of deeply-embedded low-mass protostars (Class 0) on scales of less than 300 AU is still poorly constrained. Determining this is crucial for understanding the physical and chemical evolution from cores to disks. In this study two models of the emission, a Gaussian disk intensity distribution and a parametrized power-law disk model, are fitted to sub-arcsecond resolution interferometric continuum observations of five Class 0 sources, including one source with a confirmed Keplerian disk. For reference, a spherically symmetric single power-law envelope is fitted to the larger scale ($\sim$1000 AU) emission and investigated further for one of the sources on smaller scales. A thin disk model can approximate the emission and physical structure in the inner few 100 AU scales of the studied deeply-embedded low-mass protostars and paves the way for analysis of a larger sample with ALMA. While the disk radii agree with previous estimates the masses are different for some of the sources studied. Assuming a typical temperature distribution, the fractional amount of mass in the disk above 100 K varies in between 7% to 30%. Kinematic data are needed to determine the presence of any Keplerian disk. Using previous observations of p-H$_2^{18}$O, we estimate the relative gas phase water abundances roughly an order of magnitude higher than previously inferred when both warm and cold H$_2$ was used as reference. A spherically symmetric single power-law envelope model fails to simultaneously reproduce both the small and large scale emission.

(Abridged*) Models of the young solar nebula assume a hot initial disk with most volatiles are in the gas phase. The question remains whether an actively accreting disk is warm enough to have gas-phase water up to 50 AU radius. No detailed studies have yet been performed on the extent of snowlines in an embedded accreting disk (Stage 0). Quantify the location of gas-phase volatiles in embedded actively accreting disk system. Two-dimensional physical and radiative transfer models have been used to calculate the temperature structure of embedded protostellar systems. Gas and ice abundances of H$_2$O, CO$_2$, and CO are calculated using the density-dependent thermal desorption formulation. The midplane water snowline increases from 3 to 55 AU for accretion rates through the disk onto the star between $10^{-9}$-$10^{-4} \ M_{\odot} \ {\rm yr^{-1}}$. CO$_2$ can remain in the solid phase within the disk for $\dot{M} \leq 10^{-5} \ M_{\odot} \ {\rm yr^{-1}}$ down to $\sim 20$ AU. Most of the CO is in the gas phase within an actively accreting disk independent of disk properties and accretion rate. The predicted optically thin water isotopolog emission is consistent with the detected H$_2^{18}$O emission toward the Stage 0 embedded young stellar objects, originating from both the disk and the warm inner envelope (hot core). An accreting embedded disk can only account for water emission arising from $R < 50$ AU, however, and the extent rapidly decreases for low accretion rates. Thus, the radial extent of the emission can be measured with ALMA observations and compared to this limit. Volatiles sublimate out to 50 AU in young disks and can reset the chemical content inherited from the envelope in periods of high accretion rates. A hot young solar nebula out to 30 AU can only have occurred during the deeply embedded Stage 0, not during the T-Tauri phase of our early solar system.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

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.

Abridged: Recent simulations have explored different ways to form accretion disks around low-mass stars. We aim to present observables to differentiate a rotationally supported disk from an infalling rotating envelope toward deeply embedded young stellar objects and infer their masses and sizes. Two 3D magnetohydrodynamics (MHD) formation simulations and 2D semi-analytical model are studied. The dust temperature structure is determined through continuum radiative transfer RADMC3D modelling. A simple temperature dependent CO abundance structure is adopted and synthetic spectrally resolved submm rotational molecular lines up to $J_{\rm u} = 10$ are simulated. All models predict similar compact components in continuum if observed at the spatial resolutions of 0.5-1$"$ (70-140 AU) typical of the observations to date. A spatial resolution of $\sim$14 AU and high dynamic range ($> 1000$) are required to differentiate between RSD and pseudo-disk in the continuum. The peak-position velocity diagrams indicate that the pseudo-disk shows a flatter velocity profile with radius than an RSD. On larger-scales, the CO isotopolog single-dish line profiles are similar and are narrower than the observed line widths of low-$J$ lines, indicating significant turbulence in the large-scale envelopes. However a forming RSD can provide the observed line widths of high-$J$ lines. Thus, either RSDs are common or a higher level of turbulence ($b \sim 0.8 \ {\rm km \ s^{-1}}$ ) is required in the inner envelope compared with the outer part. Multiple spatially and spectrally resolved molecular line observations are needed. The continuum data give a better estimate on disk masses whereas the disk sizes can be estimated from the spatially resolved molecular lines observations. The general observable trends are similar between the 2D semi-analytical models and 3D MHD RSD simulations.

(Abridged) Organic molecules are important constituents of protoplanetary disks. Their ro-vibrational lines observed in the near- and mid-infrared are commonly detected toward T Tauri disks. These lines are the only way to probe the chemistry in the inner few au where terrestrial planets form. To understand this chemistry, accurate molecular abundances have to be determined. This is complicated by excitation effects. Most analyses so far have made the assumption of local thermal equilibrium (LTE). Starting from estimates for the collisional rate coefficients of HCN, non-LTE slab models of the HCN emission were calculated to study the importance of different excitation mechanisms. Using a new radiative transfer model, the HCN emission from a full two-dimensional disk was then modeled to study the effect of the non-LTE excitation, together with the line formation. We ran models tailored to the T Tauri disk AS 205 (N) where HCN lines in both the 3 {\mu}m and 14 {\mu}m bands have been observed by VLT-CRIRES and the Spitzer Space Telescope. Reproducing the observed 3 {\mu}m / 14 {\mu}m flux ratios requires very high densities and kinetic temperatures ($n > 10^{14}$ cm$^{-3}$ and $T > 750$ K), if only collisional excitation is accounted for. Radiative pumping can, however, excite the lines easily out to considerable radii $\sim$ 10 au. Consequently, abundances derived from LTE and non-LTE models do not differ by more than a factor of about 3. Models with both a strongly enhanced abundance within $\sim$ 1 au (jump abundance) and constant abundance can reproduce the current observations, but future observations with the MIRI instrument on JWST and METIS on the E-ELT can easily distinguish between the scenarios and test chemical models. Depending on the scenario, ALMA can detect rotational lines within vibrationally excited levels.

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.

Context: The R CrA cloud hosts a handful of Class 0/I low-mass young stellar objects. The chemistry and physics at scales $>500$ AU in this cloud are dominated by the irradiation from the nearby Herbig Be star R CrA. The luminous large-scale emission makes it necessary to use high-resolution spectral imaging to study the chemistry and dynamics of the inner envelopes and discs of the protostars. Aims: We aim to better understand the structure of the inner regions of these protostars and, in particular, the interplay between the chemistry and the presence of discs. Methods: Using Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution spectral imaging interferometry observations, we study the molecular line and dust continuum emission at submillimetre wavelengths. Results: We detect dust continuum emission from four circumstellar discs around Class 0/I objects within the R CrA cloud. Towards IRS7B we detect C$^{17}$O emission showing a rotation curve consistent with a Keplerian disc with a well-defined edge that gives a good estimate for the disc radius at 50 AU. We derive the central object mass to $2.3M_{\odot}$ and the disc mass to $0.024M_{\odot}$. The observations are also consistent with a model of material infalling under conservation of angular momentum; however, this model provides a worse fit to the data. We also report a likely detection of faint CH$_3$OH emission towards this point source, as well as more luminous CH$_3$OH emission in an outflow orthogonal to the major axis of the C$^{17}$O emission. Conclusions: The faint CH$_3$OH emission seen towards IRS7B can be explained by a flat density profile of the inner envelope caused by the disc with a radius $\lesssim50$ AU. We propose that the regions of the envelopes where complex organic molecules are present in Class 0/I young stellar objects can become quenched as the disc grows.

(Abridged) The water deuterium fractionation (HDO/H$_2$O abundance ratio) has traditionally been used to infer the amount of water brought to Earth by comets. Measuring this ratio in deeply-embedded low-mass protostars makes it possible to probe the critical stage when water is transported from clouds to disks in which icy bodies are formed. We present sub-arcsecond resolution observations of HDO in combination with H$_2^{18}$O from the PdBI toward the three low-mass protostars NGC 1333-IRAS 2A, IRAS 4A-NW, and IRAS 4B. The resulting HDO/H$_2$O ratio is $7.4\pm2.1\times10^{-4}$ for IRAS 2A, $19.1\pm5.4\times10^{-4}$ for IRAS 4A-NW, and $5.9\pm1.7\times10^{-4}$ for IRAS 4B. Derived ratios agree with radiative transfer models within a factor of 2-4 depending on the source. Our HDO/H$_2$O ratios for the inner regions (where $T>100$ K) of four young protostars are only a factor of 2 higher than those found for pristine, solar system comets. These small differences suggest that little processing of water occurs between the deeply embedded stage and the formation of planetesimals and comets.

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

Low-mass protostars have been suggested to show highly variable accretion rates through-out their evolution. Such changes in accretion, and related heating of their ambient envelopes, may trigger significant chemical variations on different spatial scales and from source-to-source. We present images of emission from C17O, H13CO+, CH3OH, C34S and C2H toward the low-mass protostar IRAS 15398-3359 on 0.5" (75 AU diameter) scales with the Atacama Large Millimeter/submillimeter Array (ALMA) at 340 GHz. The resolved images show that the emission from H13CO+ is only present in a ring-like structure with a radius of about 1-1.5" (150-200 AU) whereas the CO and other high dipole moment molecules are centrally condensed toward the location of the central protostar. We propose that HCO+ is destroyed by water vapor present on small scales. The origin of this water vapor is likely an accretion burst during the last 100-1000 years increasing the luminosity of IRAS 15398-3359 by a factor of 100 above its current luminosity. Such a burst in luminosity can also explain the centrally condensed CH3OH and extended warm carbon-chain chemistry observed in this source and furthermore be reflected in the relative faintness of its compact continuum emission compared to other protostars.

Context: Rotationally supported disks are critical in the star formation process. The questions of when do they form and what factors influence or hinder their formation have been studied but are largely unanswered. Observations of early stage YSOs are needed to probe disk formation. Aims: VLA1623 is a triple non-coeval protostellar system, with a weak magnetic field perpendicular to the outflow, whose Class 0 component, VLA1623A, shows a disk-like structure in continuum with signatures of rotation in line emission. We aim to determine whether this structure is in part or in whole a rotationally supported disk, i.e. a Keplerian disk, and what are its characteristics. Methods: ALMA Cycle 0 Early Science 1.3 mm continuum and C$^{18}$O (2-1) observations in the extended configuration are presented here and used to perform an analysis of the disk-like structure using PV diagrams and thin disk modelling with the addition of foreground absorption. Results: The PV diagrams of the C$^{18}$O line emission suggest the presence of a rotationally supported component with a radius of at least 50 AU. Kinematical modelling of the line emission shows that the disk out to 180 AU is actually rotationally supported, with the rotation being well described by Keplerian rotation out to at least 150 AU, and the central source mass to be $\sim$0.2 M$_{sun}$ for an inclination of 55$^{\circ}$. Pure infall and conserved angular momentum rotation models are excluded. Conclusions: VLA1623A, a very young Class 0 source, presents a disk with an outer radius $R_{\rm out}$ = 180 AU with a Keplerian velocity structure out to at least 150 AU. The weak magnetic fields and recent fragmentation in this region of rho Ophiuchus may have played a lead role in the formation of the disk.

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

(Abridged) Star and planet formation theories predict an evolution in the density, temperature, and velocity structure as the envelope collapses and forms an accretion disk. The aim of this work is to model the evolution of the molecular excitation, line profiles, and related observables during low-mass star formation. Specifically, the signatures of disks during the deeply embedded stage are investigated. Semi-analytic 2D axisymmetric models have been used to describe the evolution of the density, stellar mass, and luminosity from the pre-stellar to the T-Tauri phase. A full radiative transfer calculation is carried out to accurately determine the time-dependent dust temperatures and CO abundance structure. We present non-LTE near-IR, FIR, and submm lines of CO have been simulated at a number of time steps. In contrast to the dust temperature, the CO excitation temperature derived from submm/FIR lines does not vary during the protostellar evolution, consistent with C18O observations obtained with Herschel and from ground-based telescopes. The near-IR spectra provide complementary information to the submm lines by probing not only the cold outer envelope but also the warm inner region. The near-IR high-J (>8) absorption lines are particularly sensitive to the physical structure of the inner few AU, which does show evolution. High signal-to-noise ratio subarcsec resolution data with ALMA are needed to detect the presence of small rotationally supported disks during the Stage 0 phase and various diagnostics are discussed.

(Abridged) Water is a key tracer of dynamics and chemistry in low-mass protostars, but spectrally resolved observations have so far been limited in sensitivity and angular resolution. In this first systematic survey of spectrally resolved water emission in low-mass protostellar objects, H2O was observed in the ground-state transition at 557 GHz with HIFI on Herschel in 29 embedded Class 0 and I protostars. Complementary far-IR and sub-mm continuum data (including PACS data from our program) are used to constrain the spectral energy distribution of each source. H2O intensities are compared to inferred envelope and outflow properties and CO 3-2 emission. H2O emission is detected in all objects except one. The line profiles are complex and consist of several kinematic components. The profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not the quiescent envelope. Several sources show multiple shock components in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in 7 sources (6 Class 0, 1 Class I) indicative of infalling envelopes, and regular P-Cygni profiles in 4 sources (3 Class I, 1 Class 0) indicative of expanding envelopes. "Bullets" moving at >50 km/s are seen in 4 Class 0 sources; 3 of these are new detections. In the outflow, the H2O/CO abundance ratio as a function of velocity is nearly the same for all sources, increasing from 10^-3 at <5 km/s to >10^-1 at >10 km/s. The H2O abundance in the outer envelope is low, ~10^-10. The different H2O profile components show a clear evolutionary trend: in the Class 0 sources, emission is dominated by outflow components originating inside an infalling envelope. When the infall diminishes during the Class I phase, the outflow weakens and H2O emission disappears.

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.

Gravitational instability plays an important role in driving gas accretion in massive protostellar discs. Particularly strong is the global gravitational instability, which arises when the disc mass is of order 0.1 of the mass of the central star and has a characteristic spatial scale much greater than the disc's vertical scale-height. In this paper we use three-dimensional numerical hydrodynamics to study the development of gravitational instabilities in a disc which is embedded in a dense, gaseous envelope. We find that global gravitational instabilities are the dominant mode of angular momentum transport in the disc with infall, in contrast to otherwise identical isolated discs. The accretion torques created by low-order, global modes of the gravitational instability in a disc subject to infall are larger by a factor of several than an isolated disc of the same mass. We show that this global gravitational instability is driven by the strong vertical shear at the interface between the disc and the envelope, and suggest that this process may be an important means of driving accretion on to young stars.

We summarize the results from analyzing six clusters of galaxies at 0.14 < z < 0.40 observed with the Hubble Space Telescope Advanced Camera for Surveys. We derive deep composite luminosity functions in B,g,V,r,i and z down to an absolute magnitude of -14 +5 log h mag. The luminosity functions are fitted by a single Schechter function with M^*=-19.8,-20.9,-21.9, -22.0,-21.7 and -22.3 i B,g,V.r,i,z respectively and alpha=-1.3 for all bands. The data suggest that the red sequence dominates the luminosity function down to more than 6 mag. below L*, the dwarf spheroidal regime. Hence, at least at z=0.3 the red sequence is well established and galaxies down to dwarf spheroidals are fully assembled within these clusters. We do not detect the faint-end upturn (M > -16) that is observed in lower redshift clusters. If this is real, the faint-end population has originated since z = 0.3.

We present deep composite luminosity functions in $B$, $g$, $V$, $r$, $i$ and $z$ for six clusters at $0.14 < z < 0.40$ observed with the Hubble Space Telescope Advanced Camera for Surveys. The luminosity functions reach to absolute magnitude of $\sim -14 + 5\log h$ mag. and are well fitted by a single Schechter function with $M^*_{BgVriz}=-19.8, -20.9 -21.9, -22.0, -21.7, -22.3$ mag. and $\alpha \sim -1.3$ (in all bands). The observations suggest that the galaxy luminosity function is dominated by objects on the red sequence to at least 6 mags. below the $L^*$ point. Comparison with local data shows that the red sequence is well established at least at $z \sim 0.3$ down to $\sim 1/600^{th}$ of the luminosity of the Milky Way and that galaxies down to the regime of dwarf spheroidals have been completely assembled in clusters at this redshift. We do not detect a steepening of the luminosity function at $M > -16$ as is observed locally. If the faint end upturn is real, the steepening of the luminosity function must be due to a newly infalling population of faint dwarf galaxies.

We derive deep luminosity functions (to $M_z=-15$) for galaxies in Abell 1835 ($z=0.25$) and AC 114 ($z=0.31$) and compare these with the local $z'$ luminosity function for 69 clusters. The data show that the faint-end upturn, the excess of galaxies above a single Schechter function at $M_z < -17$, does not exist in the higher redshift clusters. This suggests that the faint-end upturn galaxies have been created recently, by infall into clusters of star-forming field populations or via tidal disruption of brighter objects.^M