Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.

We present ALMA observations of organic molecules towards five low-mass Class 0/I protostellar disk candidates in the Serpens cluster. Three sources (Ser-emb 1, Ser-emb 8, and Ser-emb 17) present emission of CH3OH as well as CH3OCH3, CH3OCHO, and CH2CO, while NH2CHO is detected in just Ser-emb 8 and Ser-emb 17. Detecting hot corino-type chemistry in three of five sources represents a high occurrence rate given the relative sparsity of these sources in the literature, and this suggests a possible link between protostellar disk formation and hot corino formation. For sources with CH3OH detections, we derive column densities of 10^{17}-10^{18} cm^{-2} and rotational temperatures of ~200-250 K. The CH3OH-normalized column density ratios of large, oxygen-bearing COMs in the Serpens sources and other hot corinos span two orders of magnitude, demonstrating a high degree of chemical diversity at the hot corino stage. Resolved observations of a larger sample of objects are needed to understand the origins of chemical diversity in hot corinos, and the relationship between different protostellar structural elements on disk-forming scales.

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

Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morphology that cannot be accounted for in traditional, spherical 1D models of the envelope. Aims. The purpose of this paper is to study the environment of the two components of the binary through 3D radiative transfer modeling and to compare with data from the Atacama Large Millimeter/submillimeter Array. Such comparisons can be used to constrain the protoplanetary disk structures, the luminosities of the two components of the binary and the chemistry of simple species. Methods. We present 13CO, C17O and C18O J=3-2 observations from the ALMA Protostellar Interferometric Line Survey (PILS), together with a qualitative study of the dust and gas density distribution of IRAS 16293-2422. A 3D dust and gas model including disks and a dust filament between the two protostars is constructed which qualitatively reproduces the dust continuum and gas line emission. Results and conclusions. Radiative transfer modeling of source A and B, with the density solution of an infalling, rotating collapse or a protoplanetary disk model, can match the constraints for the disk-like emission around source A and B from the observed dust continuum and CO isotopologue gas emission. If a protoplanetary disk model is used around source B, it has to have an unusually high scale-height in order to reach the dust continuum peak emission value, while fulfilling the other observational constraints. Our 3D model requires source A to be much more luminous than source B; LA ~ 18 $L_\odot$ and LB ~ 3 $L_\odot$.

FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars.

A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.

H2CO is one of the most readily detected organic molecules in protoplanetary disks. Yet its distribution and dominant formation pathway(s) remain largely unconstrained. To address these issues, we present ALMA observations of two H2CO lines (3_{12}-2_{11} and 5_{15}-4_{14}) at 0".5 (~30 au) spatial resolution toward the disk around the nearby T Tauri star TW Hya. Emission from both lines is spatially resolved, showing a central depression, a peak at 0".4 radius, and a radial decline at larger radii with a bump at ~1", near the millimeter continuum edge. We adopt a physical model for the disk and use toy models to explore the radial and vertical H2CO abundance structure. We find that the observed emission implies the presence of at least two distinct H2CO gas reservoirs: (1) a warm and unresolved inner component (<10 au), and (2) an outer component that extends from ~15 au to beyond the millimeter continuum edge. The outer component is further constrained by the line ratio to arise in a more elevated disk layer at larger radii. The inferred H2CO abundance structure agrees well with disk chemistry models, which predict efficient H2CO gas-phase formation close to the star, and cold H2CO grain surface formation, through H additions to condensed CO, followed by non-thermal desorption in the outer disk. The implied presence of active grain surface chemistry in the TW Hya disk is consistent with the recent detection of CH3OH emission, and suggests that more complex organic molecules are formed in disks, as well.

Recent high angular resolution ($\sim$0.2") ALMA observations of the 1.1 mm continuum and of HCO+ J=3-2 and HCN J=3-2 gas towards the binary protostar IRS 43 reveal multiple Keplerian disks which are significantly misaligned ($\gt$ 60$^\circ$), both in inclination and position angle and also with respect to the binary orbital plane. Each stellar component has an associated circumstellar disk while the binary is surrounded by a circumbinary disk. Together with archival VLA measurements of the stellar positions over 25 years, and assuming a circular orbit, we use our continuum measurements to determine the binary separation, a = 74 $\pm$ 4 AU, and its inclination, i $\lt$ 30$^\circ$. The misalignment in this system suggests that turbulence has likely played a major role in the formation of IRS 43.

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

We present SIGAME (SImulator of GAlaxy Millimetre/submillimetre Emission), a new numerical code designed to simulate the 12CO rotational line emission spectrum of galaxies. Using sub-grid physics recipes to post-process the outputs of smoothed particle hydrodynamics (SPH) simulations, a molecular gas phase is condensed out of the hot and partly ionized SPH gas. The gas is subjected to far-UV radiation fields and cosmic ray ionization rates which are set to scale with the local star formation rate volume density. Level populations and radiative transport of the CO lines are solved with the 3-D radiative transfer code LIME. We have applied SIGAME to cosmological SPH simulations of three disc galaxies at z=2 with stellar masses in the range ~(0.5-2)x10^11 Msun and star formation rates ~40-140 Msun/yr. Global CO luminosities and line ratios are in agreement with observations of disc galaxies at z~2 up to and including J=3-2 but falling short of the few existing J=5-4 observations. The central 5 kpc regions of our galaxies have CO 3-2/1-0 and 7-6/1-0 brightness temperature ratios of ~0.55-0.65 and ~0.02-0.08, respectively, while further out in the disc the ratios drop to more quiescent values of ~0.5 and <0.01. Global CO-to-H2 conversion (alpha_CO) factors are ~=1.5 Msun*pc^2/(K km s/1), i.e. ~2-3 times below typically adopted values for disc galaxies, and alpha_CO increases with radius, in agreement with observations of nearby galaxies. Adopting a top-heavy Giant Molecular Cloud (GMC) mass spectrum does not significantly change the results. Steepening the GMC density profile leads to higher global line ratios for J_up>=3 and CO-to-H2 conversion factors [~=3.6 Msun*pc^2/(K km/s)].

Changing physical conditions in the vicinity of protostars allow for a rich and interesting chemistry to occur. Heating and cooling of the gas allows molecules to be released from and frozen out on dust grains. These changes in physics, traced by chemistry, as well as the kinematical information allows us to distinguish between different scenarios describing the infall of matter and the launching of molecular outflows and jets. We aim at determining the spatial distribution of different species, of different chemical origin. This is to examine the physical processes in play in the observed region. From the kinematical information of the emission lines we aim at determining the nature of the infalling and outflowing gas in the system. We also aim at determining the physical properties of the outflow. Maps from the Sub-Millimeter Array reveal the spatial distribution of the gaseous emission toward IRAS15398-3359. The line radiative transfer code LIME is used to construct a full 3D model of the system taking all relevant components and scales into account. CO, HCO+ and N2H+ are detected and are shown to trace the motions of the outflow. For CO, also the circumstellar envelope and the surrounding cloud have a profound impact on the observed line profiles. N2H+ is detected in the outflow, but is suppressed towards the central region, perhaps due to the competing reaction between CO and H3+ in the densest regions as well as destruction of N2H+ by CO. N2D+ is detected in a ridge south-west from the protostellar condensation. The morphology and kinematics of the CO emission suggests that the source is younger than 1000 years. The mass, momentum, momentum rate, mechanical luminosity, kinetic energy and mass-loss rate are also all estimated to be low. A full 3D radiative transfer model of the system can explain all the kinematical and morphological features in the system.

We present S\'IGAME simulations of the [CII]157.7$\mu$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase.

The HH54 shock is a Herbig-Haro object, located in the nearby Chamaeleon II cloud. Observed CO line profiles are due to a complex distribution in density, temperature, velocity, and geometry. Resolving the HH54 shock wave in the far-infrared cooling lines of CO constrain the kinematics, morphology, and physical conditions of the shocked region. We used the PACS and SPIRE instruments on board the Herschel space observatory to map the full FIR spectrum in a region covering the HH54 shock wave. Complementary Herschel-HIFI, APEX, and Spitzer data are used in the analysis as well. The observed features in the line profiles are reproduced using a 3D radiative transfer model of a bow-shock, constructed with the Line Modeling Engine code (LIME). The FIR emission is confined to the HH54 region and a coherent displacement of the location of the emission maximum of CO with increasing J is observed. The peak positions of the high-J CO lines are shifted upstream from the lower J CO lines and coincide with the position of the spectral feature identified previously in CO(10-9) profiles with HIFI. This indicates a hotter molecular component in the upstream gas with distinct dynamics. The coherent displacement with increasing J for CO is consistent with a scenario where IRAS12500-7658 is the exciting source of the flow, and the 180 K bow-shock is accompanied by a hot (800 K) molecular component located upstream from the apex of the shock and blueshifted by -7 km s$^{-1}$. The spatial proximity of this knot to the peaks of the atomic fine-structure emission lines observed with Spitzer and PACS ([OI]63, 145 $\mu$m) suggests that it may be associated with the dissociative shock as the jet impacts slower moving gas in the HH54 bow-shock.

A fundamental question about the early evolution of low-mass protostars is when circumstellar disks may form. High angular resolution observations of molecular transitions in the (sub)millimeter wavelength windows make it possible to investigate the kinematics of the gas around newly-formed stars, for example to identify the presence of rotation and infall. IRAS 16293-2422 was observed with the extended Submillimeter Array (eSMA) resulting in subarcsecond resolution (0.46" x 0.29", i.e. $\sim$ 55 $\times$ 35~AU) images of compact emission from the C$^{17}$O (3-2) and C$^{34}$S (7-6) transitions at 337~GHz (0.89~mm). To recover the more extended emission we have combined the eSMA data with SMA observations of the same molecules. The emission of C$^{17}$O (3-2) and C$^{34}$S (7-6) both show a velocity gradient oriented along a northeast-southwest direction with respect to the continuum marking the location of one of the components of the binary, IRAS16293A. Our combined eSMA and SMA observations show that the velocity field on the 50--400~AU scales is consistent with a rotating structure. It cannot be explained by simple Keplerian rotation around a single point mass but rather needs to take into account the enclosed envelope mass at the radii where the observed lines are excited. We suggest that IRAS 16293-2422 could be among the best candidates to observe a pseudo-disk with future high angular resolution observations.

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.

Interferometers play an increasingly important role for spatially resolved observations. If employed at full potential, interferometry can probe an enormous dynamic range in spatial scale. Interpretation of the observed visibilities requires the numerical compu- tation of Fourier integrals over the synthetic model images. To get the correct values of these integrals, the model images must have the right size and resolution. Insufficient care in these choices can lead to wrong results. We present a new general-purpose scheme for the computation of visibilities of radiative transfer images. Our method requires a model image that is a list of intensities at arbitrarily placed positions on the image-plane. It creates a triangulated grid from these vertices, and assumes that the intensity inside each triangle of the grid is a linear function. The Fourier integral over each triangle is then evaluated with an analytic expression and the complex visibility of the entire image is then the sum of all triangles. The result is a robust Fourier trans- form that does not suffer from aliasing effects due to grid regularities. The method automatically ensures that all structure contained in the model gets reflected in the Fourier transform.

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.

A fundamental part of the study of star formation is to place young stellar objects in an evolutionary sequence. Establishing a robust evolutionary classification scheme allows us not only to understand how the Sun was born but also to predict what kind of main sequence star a given protostar will become. Traditionally, low-mass young stellar objects are classified according to the shape of their SEDs. Such methods are however prone to misclassification due to degeneracy and do not constrain the temporal evolution. More recently, young stellar objects have been classified based on envelope, disk, and stellar masses determined from resolved images of their continuum and line emission at submillimeter wavelengths. Through detailed modeling of two Class I sources, we aim at determining accurate velocity profiles and explore the role of freeze-out chemistry in such objects. We present new Submillimeter Array observations of the continuum and HCO$^+$ line emission at 1.1 mm toward two protostars, IRS 63 and IRS 43 in the Ophiuchus star forming region. The sources are modeled in detail using dust radiation transfer to fit the SED and continuum images and line radiation transfer to produce synthetic position-velocity diagrams. We use a $\chi^2$ search algorithm to find the best model fit to the data and to estimate the errors in the model variables. Our best fit models present disk, envelope and stellar masses as well as the HCO$^+$ abundance and inclination of both sources. We also identify a ring structure with a radius of about 200 AU in IRS 63. We find that freeze-out chemistry is important in IRS 63 but not for IRS 43. We show that the velocity field in IRS 43 is consistent with Keplerian rotation. Due to molecular depletion, it is not possible to draw a similar conclusion for IRS 63. We identify a ring shaped structure in IRS 63 on the same spatial scale as the disk outer radius.

The goal of this paper is to determine the respective importance of solid state vs. gas phase reactions for the formation of dimethyl ether. This is done by a detailed analysis of the excitation properties of the ground state and the torsionally excited states, v11=1 and v15=1, toward the high-mass star-forming region G327.3-0.6. With the Atacama Pathfinder EXperiment 12 m submillimeter telescope, we performed a spectral line survey. The observed spectrum is modeled assuming local thermal equilibrium. CH3OCH3 has been detected in the ground state, and in the torsionally excited states v11=1 and v15=1, for which lines have been detected here for the first time. The emission is modeled with an isothermal source structure as well as with a non-uniform spherical structure. For non-uniform source models one abundance jump for dimethyl ether is sufficient to fit the emission, but two components are needed for the isothermal models. This suggests that dimethyl ether is present in an extended region of the envelope and traces a non-uniform density and temperature structure. Both types of models furthermore suggest that most dimethyl ether is present in gas that is warmer than 100 K, but a smaller fraction of 5%-28% is present at temperatures between 70 and 100 K. The dimethyl ether present in this cooler gas is likely formed in the solid state, while gas phase formation probably is dominant above 100 K. Finally, the v11=1 and v15=1 torsionally excited states are easily excited under the density and temperature conditions in G327.3-0.6 and will thus very likely be detectable in other hot cores as well.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

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.

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

(Abridged) As the stellar X-ray and UV light penetration of a protoplanetary disk depends sensitively on the dust properties, trace molecular species like HCO+, HCN, and CN are expected to show marked differences from photoprocessing effects as the dust content in the disk evolves. We investigate the evolution of the UV irradiation of the molecular gas in a sample of classical T Tauri stars in Taurus that exhibit a wide range in grain growth and dust settling properties. We obtained HCO+ (J=3-2), HCN (J=3-2), and CN (J=2-1) observations of 13 sources with the JCMT. Our sample has 1.3mm fluxes in excess of 75mJy, indicating the presence of significant dust reservoirs; a range of dust settling as traced through their spectral slopes between 6, 13, and 25 microns; and varying degrees of grain growth as extrapolated from the strength of their 10-micron silicate emission features. We compare the emission line strengths with the sources' continuum flux and infrared features, and use detailed modeling based on two different model prescriptions to compare typical disk abundances for HCO+, HCN, and CN with the gas-line observations for our sample. We detected HCO+ (3-2) toward 6 disks, HCN (3-2) from 0 disks, and CN (2-1) toward 4 disks. For the complete sample, there is no correlation between the gas-line strengths or their ratios and either the sources' dust continuum flux or infrared slopes.

We present a new code for solving the molecular and atomic excitation and radiation transfer problem in a molecular gas and predicting emergent spectra. This code works in arbitrary three dimensional geometry using unstructured Delaunay latices for the transport of photons. Various physical models can be used as input, ranging from analytical descriptions over tabulated models to SPH simulations. To generate the Delaunay grid we sample the input model randomly, but weigh the sample probability with the molecular density and other parameters, and thereby we obtain an average grid point separation that scales with the local opacity. Our code does photon very efficiently so that the slow convergence of opaque models becomes traceable. When convergence between the level populations, the radiation field, and the point separation has been obtained, the grid is ray-traced to produced images that can readily be compared to observations. Because of the high dynamic range in scales that can be resolved using this type of grid, our code is particularly well suited for modeling of ALMA data. Our code can furthermore deal with overlapping lines of multiple molecular and atomic species.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

Observations of Serpens have been performed at the JCMT using Harp-B. Maps over a 4.5'x5.4' region were made in a frequency window around 338 GHz, covering the 7-6 transitions of methanol. Emission is extended over each source, following the column density of H2 but showing up also particularly strongly around outflows. The rotational temperature is low, 15-20 K, and does not vary with position within each source. The abundance is typically 10^-9 - 10^-8 with respect to H2 in the outer envelope, whereas "jumps" by factors of up to 10^2 -10^3 inside the region where the dust temperature exceeds 100 K are not excluded. A factor of up to ~ 10^3 enhancement is seen in outflow gas. In one object, SMM4, the ice abundance has been measured to be ~ 3x10^-5 with respect to H2 in the outer envelope, i.e., a factor of 10^3 larger than the gas-phase abundance. Comparison with C18O J=3-2 emission shows that strong CO depletion leads to a high gas-phase abundance of CH3OH not just for the Serpens sources, but for a larger sample of protostars. The observations illustrate the large-scale, low-level desorption of CH3OH from dust grains, extending out to and beyond 7500 AU from each source, a scenario which is consistent with non-thermal (photo-)desorption from the ice. The observations also illustrate the usefulness of CH3OH as a tracer of energetic input in the form of outflows, where methanol is sputtered from the grain surfaces. Finally, the observations provide further evidence of CH3OH formation through CO hydrogenation proceeding on grain surfaces in low-mass envelopes.

Low-mass star formation is described by gravitational collapse of dense cores of gas and dust. At some point during the collapse, a disk is formed around the protostar and the disk will spin up and grow in size as the core contracts because of angular momentum conservation. The question is how early the disk formation process occurs. In this paper we aim to characterize the kinematical state of a deeply embedded, Class 0 young stellar object, NGC1333-IRAS2A, based on high angular resolution (< 1$''$ $\approx$ 200 AU) interferometric observations of HCN and H$^{13}$CN J = 4-3 from the Submillimeter Array, and test whether a circumstellar disk can be detected based on gas kinematic features. We adopt a physical model which has been shown to describe the object well and obtain a fit of a parameterized model of the velocity field, using a two-dimensional axis-symmetric radiation transfer code. The parameterization and fit to the high angular resolution data characterize the central dynamical mass and the ratio of infall velocity to rotation velocity. We find a large amount of infall and very little rotation on all scales. The central object has a relatively low mass of 0.25 M$_\odot$ . As an object with a low stellar mass compared to the envelope mass, we conclude that NGC1333-IRAS2A is consistent with the suggestion that, as a Class 0 object, it represents the earliest stages of star formation. The large amount of infall relative to rotation also suggests that this is a young object. We do however find the need of a central compact component on scales of a few hundred AU based on the continuum data, which suggests that disk formation happens shortly after the initial gravitational collapse. The data do not reveal a distinct velocity field for this 0.1 M$_\odot$ component.

The physical conditions in a collapsing cloud can be traced by observations of molecular lines. To correctly interpret these observations the abundance distributions of the observed species need to be derived. The chemistry in a collapsing molecular cloud is not in a steady state as the density and temperature evolve. We therefore need to follow chemical reactions, both in the gas phase and on dust grains, as well as gas-grain interactions, to predict the abundance distributions. Our aim is to model the abundances of molecules, in the gas phase and on grain mantles in the form of ice, from prestellar core collapse to disk formation. We use a 2-dimensional hydrodynamical simulation as a physical model from which we take the density, temperature, and the flow of the gas. Trace particles, moving along with the gas, are used to follow the chemistry during prestellar core collapse and disk formation. The evolution of the abundances and the composition of ices on grain mantles were compared to observations. We also investigated the initial abundances to be adopted in more detailed modeling of protoplanetary disks by following the chemical evolution of trace particles accreting onto the disk. Fractional abundances of HCO+, N2H+, H2CO, HC3N, and CH3OH from our model with grain surface reactions provide a good match to observations, while abundances of CO, CS, SO, HCN, and HNC show better agreement without grain surface reactions. The observed mantle composition of dust grains is best reproduced when we include surface reactions. The initial chemical abundances to be used for detailed modeling of a protoplanetary disk are found to be different from those in dark interstellar clouds. Ices with a binding energy lower than about 1200 K sublimate before accreting onto the disk, while those with a higher binding energy do not.

Understanding the gas abundance distribution is essential when tracing star formation using molecular line observations. Changing density and temperature conditions cause gas to freeze-out onto dust grains, and this needs to be taken into account when modeling a collapsing molecular cloud. This study aims to provide a realistic estimate of the CO abundance distribution throughout the collapse of a molecular cloud. We provide abundance profiles and synthetic spectral lines which can be compared to observations. We use a 2D hydrodynamical simulation of a collapsing cloud and subsequent formation of a protoplanetary disk as input for the chemical calculations. From the resulting abundances, synthetic spectra are calculated using a molecular excitation and radiation transfer code. We compare three different methods to calculate the abundance of CO. Our models also consider cosmic ray desorption and the effects of an increased CO binding energy. The resulting abundance profiles are compared to observations from the literature and are found to agree well. The resulting abundance profiles agree well with analytic approximations, and the corresponding line fluxes match observational data. Our developed method to calculate abundances in hydrodynamical simulations should greatly aid in comparing these to observations, and can easily be generalized to include gas-phase reaction networks.

When low-mass stars form, the collapsing cloud of gas and dust goes through several stages which are usually characterized by the shape of their spectral energy distributions. Such classification is based on the cloud morphology only and does not address the dynamical state of the object. In this paper we investigate the initial cloud collapse and subsequent disk formation through the dynamical behavior as reflected in the sub-millimeter spectral emission line profiles. If a young stellar object is to be characterized by its dynamical structure it is important to know how accurately information about the velocity field can be extracted and which observables provide the best description of the kinematics. Of particular interest is the transition from infalling envelope to rotating disk, because this provides the initial conditions for the protoplanetary disk, such as mass and size. We use a hydrodynamical model, describing the collapse of a core and formation of a disk, to produce synthetic observables which we compare to calculated line profiles of a simple parameterized model. Because we know the velocity field from the hydrodynamical simulation we can determine in a quantitative way how well our best-fit parameterized velocity field reproduces the original. We use a molecular line excitation and radiation transfer code to produce spectra of both our hydro dynamical simulation as well as our parameterized model. We find that information about the velocity field can reasonably well be derived by fitting a simple model to either single-dish lines or interferometric data, but preferentially by using a combination of the two. Our result shows that it is possible to establish relative ages of a sample of young stellar objects using this method, independently of the details of the hydrodynamical model.

Circumstellar disks are expected to form early in the process that leads to the formation of a young star, during the collapse of the dense molecular cloud core. It is currently not well understood at what stage of the collapse the disk is formed or how it subsequently evolves. We aim to identify whether an embedded Keplerian protoplanetary disk resides in the L1489 IRS system. Given the amount of envelope material still present, such a disk would respresent a very young example of a protoplanetary disk. Using the Submillimeter Array (SMA) we have observed the HCO$^+$ $J=$ 3--2 line with a resolution of about 1$''$. At this resolution a protoplanetary disk with a radius of a few hundred AUs should be detectable, if present. Radiative transfer tools are used to model the emission from both continuum and line data. We find that these data are consistent with theoretical models of a collapsing envelope and Keplerian circumstellar disk. Models reproducing both the SED and the interferometric continuum observations reveal that the disk is inclined by 40$^\circ$ which is significantly different to the surrounding envelope (74$^\circ$). This misalignment of the angular momentum axes may be caused by a gradient within the angular momentum in the parental cloud or if L1489 IRS is a binary system rather than just a single star. In the latter case, future observations looking for variability at sub-arcsecond scales may be able to constrain these dynamical variations directly. However, if stars form from turbulent cores, the accreting material will not have a constant angular momentum axis (although the average is well defined and conserved) in which case it is more likely to have a misalignment of the angular momentum axes of the disk and the envelope.

We investigate the velocity transition in the low-mass protostar L1489 IRS, which is known to be embedded in a flattened, disc-like structure that shows both infall and rotation. We construct a model for L1489 IRS consisting of an flattened envelope and a velocity field that can vary from pure infall to pure rotation. We obtain best-fit parameters by comparison to 24 molecular transitions from the literature, and using a molecular excitation code and a Voronoi optimisation algorithm. We test the model against existing millimeter interferometric observations, near-infrared scattered light imaging, and 12CO ro-vibrational lines.We find that L1489 IRS is well described by a central stellar mass of 1.3M$_\odot$ surrounded by a 0.10M$_\odot$ flattened envelope with approximate scale height h\approx 0.57 R, inclined at 74^\circ. The velocity field is strongly dominated by rotation, with the velocity vector making an angle of 15^\circ with the azimuthal direction. Reproducing low-excitation transitions requires that the emission and absorption by the starless core 1' (8400 AU) east of L1489 IRS is included properly, implying that L1489 IRS is located partially behind this core. We speculate that L1489 IRS was originally formed closer to the center of this core, but has migrated to its current position over the past few times 10^5 yr, consistent with their radial velocity difference of 0.4 kms-1. This suggests that L1489 IRS' unusual appearance may be result of its migration, and that it would appear as a `normal' embedded protostar if it were still surrounded by an extended cloud core. Conversely, we hypothesize that the inner envelopes of embedded protostars resemble the rotating structure seen around L1489 IRS.

AIMS : To characterize the molecular environment of classical T Tauri stars in Lupus observed with the Spitzer Space Telescope and to search for gas-rich disks toward these sources. METHODS : Submillimeter observations of 12CO, 13CO and C18O in the J=3-2 and 2-1 lines were obtained with the James Clerk Maxwell Telescope toward a sample of 21 T Tauri stars with disks in the Lupus molecular clouds. Pointings at the sources and at selected off-positions are presented in order to disentangle material associated with disks from ambient cloud material. RESULTS : One source, IM Lup (Sz 82), has been found with double-peaked 12CO and 13CO profiles characteristic of a large rotating gas disk. The inclination of the disk is found to be ~20 degrees, with an outer radius of 400-700 AU. For most other sources, including GQ Lup with its substellar companion, the single-dish 12CO lines are dominated by extended cloud emission with a complex velocity structure. No evidence for molecular outflows is found. Compact C18O emission due to dense circumstellar material has been detected toward only two sources. Future searches for gas-rich disks in Lupus should either use interferometers or perform very deep single-dish integrations in dense gas tracers to separate the cloud and disk emission.