Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions.

We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale.

The accuracy of masses of pre-main sequence (PMS) stars derived from their locations on the Hertzsprung-Russell Diagram (HRD) can be tested by comparison with accurate and precise masses determined independently. We present 29 single stars in the Taurus star-forming region (SFR) and 3 in the Ophiuchus SFR with masses measured dynamically to a precision of at least $10 \%$. Our results include 9 updated mass determinations and 3 that have not had their dynamical masses published before. This list of stars with fundamental, dynamical masses, M$_{dyn}$, is drawn from a larger list of 39 targets in the Taurus SFR and 6 in the Ophiuchus SFR. Placing the stars with accurate and precise dynamical masses on HRDs that do not include internal magnetic fields underestimates the mass compared to M$_{dyn}$ by about $30 \%$. Placing them on an HRD that does include magnetic fields yields mass estimates in much better agreement with M$_{dyn}$, with an average difference between M$_{dyn}$ and the estimated track mass of $0.01\pm0.02$~\msun. The ages of the stars, 3--10 MY on tracks that include magnetic fields, is older than the 1--3 MY indicated by the non-magnetic models. The older ages of T Tauri stars predicted by the magnetic models increase the time available for evolution of their disks and formation of the giant gas exoplanets. The agreement between our M$_{dyn}$ values and the masses on the magnetic field tracks provides indirect support for these older ages.

Context. With its high complexity, large size, and close distance, the ringworld around GG Tau A is an appealing case to study the formation and evolution of protoplanetary disks around multiple star systems. However, investigations with radiative transfer models are usually neglecting the influence of the circumstellar dust around the individual stars. Aims. We investigate how circumstellar disks around the stars of GG Tau A are influencing the emission that is scattered at the circumbinary disk and if constraints on these circumstellar disks can be derived. Methods. We perform radiative transfer simulations with the code POLARIS to obtain spectral energy distributions and emission maps in the H-Band (near-infrared). Subsequently, we compare them with observations to achieve our aims. Results. We studied the ratio of polarized intensity at different locations in the circumbinary disk and conclude that the observed scattered-light near-infrared emission is best reproduced, if the circumbinary disk lies in the shadow of at least two co-planar circumstellar disks surrounding the central stars. This implies that the inner wall of the circumbinary disk is strongly obscured around the midplane, while the observed emission is actually dominated by the most upper disk layers. In addition, the inclined dark lane ("gap") on the western side of the circumbinary disk, which is a stable (non rotating) feature since ~20yr, can only be explained by the self-shadowing of a misaligned circumstellar disk surrounding one of the two components of the secondary close-binary star GG Tau Ab.

Protoplanetary disks are challenging objects for astrochemical models due to strong density and temperature gradients and due to the UV photons 2D propagation. In this paper, we have studied the importance of several model parameters on the predicted column densities of observed species. We considered: 1) 2-phase (gas and homogeneous grains) or 3-phase (gas, surface, and bulk of grains) models, 2) several initial compositions, 3) grain growth and dust settling, and 4) several cosmic-ray ionization rates. Our main result is that dust settling is the most crucial parameter. Including this effect renders the computed column densities sensitive to all the other model parameters, except cosmic-ray ionization rate. In fact, we found almost no effect of this parameter for radii larger than 10 au (the minimum radius studied here) except for N2H+. We also compared all our models with all the column densities observed in the protoplanetary disk around DM Tau and were not able to reproduce all the observations despite the studied parameters. N2H+ seems to be the most sensitive species. Its observation in protoplanetary disks at large radius could indicate enough N2 in the gas-phase (inhibited by the 3-phase model, but boosted by the settling) and a low electron abundance (favored by low C and S elemental abundances).

Observations of multiple rotational transitions from a single molecule allow for unparalleled constraints on the physical conditions of the emitting region. We present an analysis of CS in TW~Hya using the $J=7-6$, $5-4$ and $3-2$ transitions imaged at $\sim 0.5^{\prime\prime}$ spatial resolution, resulting in a temperature and column density profile of the CS emission region extending out to 230~au, far beyond previous measurements. In addition, the 15~kHz resolution of the observations and the ability to directly estimate the temperature of the CS emitting gas, allow for one of the most sensitive searches for turbulent broadening in a disk to date. Limits of $v_{\rm turb} \lesssim 0.1 c_s$ can be placed across the entire radius of the disk. We are able to place strict limits of the local H$_2$ density due to the collisional excitations of the observed transitions. From these we find that a minimum disk mass of $3 \times 10^{-4}~M_{\rm sun}$ is required to be consistent with the CS excitation conditions and can uniquely constrain the gas surface density profile in the outer disk.

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems.

Context. Several sulfur-bearing molecules are observed in the interstellar medium and in comets, in strong contrast to protoplanetary disks where only CS, H$_2$CS and SO have been detected so far. Aims. We combine observations and chemical models to constrain the sulfur abundances and their sensitivity to physical and chemical conditions in the DM Tau protoplanetary disk. Methods. We obtained $0.5^{"}$ ALMA observations of DM Tau in Bands 4 and 6 in lines of CS, SO, SO$_2$, OCS, CCS, H$_2$CS and H$_2$S, achieving a $\sim 5$ mJy sensitivity. Using the non-LTE radiative transfer code RADEX and the forward-modeling tool DiskFit, disk-averaged CS column densities and upper limits for the other species were derived. Results. Only CS was detected with a derived column density of $\sim 2-6 \times 10^{12}$ cm$^{-2}$. We report a first tentative detection of SO$_2$ in DM Tau. The upper limits range between $\sim 10^{11}$ and $10^{14}$ cm$^{-2}$ for the other S-bearing species. The best-fit chemical model matching these values requires a gas-phase C/O ratio of > 1 at $r \sim 50-100$ au. With chemical modeling we demonstrate that sulfur-bearing species could be robust tracers of the gas-phase C/O ratio, surface reaction rates, grain size and UV intensities. Conclusions. The lack of detections of a variety of sulfur-bearing molecules in DM Tau other than CS implies a dearth of reactive sulfur in the gas phase, either through efficient freeze-out or because most of the elemental sulfur is in other large species, as found in comets. The inferred high CS/SO and CS/SO$_2$ ratios require a non-solar C/O gas-phase ratio of > 1, consistent with the recent observations of hydrocarbon rings in DM Tau. The stronger depletion of oxygen-bearing S-species compared to CS is likely linked to the low observed abundances of gaseous water in DM Tau and points to a removal mechanism of oxygen from the gas.

Observations by the Atacama Large Millimeter/sub-millimeter Array of the dust continuum and $^{13}$CO(3-2) millimeter emissions of the triple stellar system GG Tau A are analysed, giving evidence for a rotating gas disc and a concentric and coplanar dust ring. The present work complements an earlier analysis (Tang et al. 2016) by exploring detailed properties of the gas disc. A 95% confidence level upper limit of 0.24 arcsec (34 au) is placed on the disc scale height at a distance of 1 arcsec (140 au) from the central stars. Evidence for Keplerian rotation of the gas disc is presented, the rotation velocity reaching ~3.1 kms$^{-1}$ at 1 arcsec from the central stars, and a 99% confidence level upper limit of 9% is placed on a possible in-fall velocity relative contribution. Variations of the intensity across the disc area are studied in detail and confirm the presence of a hot spot in the south-eastern quadrant. However several other significant intensity variations, in particular a depression in the northern direction, are also revealed. Variations of the intensity are found to be positively correlated to variations of the line width. Possible contributions to the measured line width are reviewed, suggesting an increase of the disc temperature and opacity with decreasing distance from the stars.

We report new dynamical masses for 5 pre-main sequence (PMS) stars in the L1495 region of the Taurus star-forming region (SFR) and 6 in the L1688 region of the Ophiuchus SFR. Since these regions have VLBA parallaxes these are absolute measurements of the stars' masses and are independent of their effective temperatures and luminosities. Seven of the stars have masses $<0.6$ solar masses, thus providing data in a mass range with little data, and of these, 6 are measured to precision $< 5 \%$. We find 8 stars with masses in the range 0.09 to 1.1 solar mass that agree well with the current generation of PMS evolutionary models. The ages of the stars we measured in the Taurus SFR are in the range 1-3 MY, and $<1$ MY for those in L1688. We also measured the dynamical masses of 14 stars in the ALMA archival data for Akeson~\&~Jensen's Cycle 0 project on binaries in the Taurus SFR. We find that the masses of 7 of the targets are so large that they cannot be reconciled with reported values of their luminosity and effective temperature. We suggest that these targets are themselves binaries or triples.

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions.

Accurate temperature calculations for circumstellar disks are particularly important for their chemical evolution. Their temperature distribution is determined by the optical properties of the dust grains, which, among other parameters, depend on their radius. However, in most disk studies, only average optical properties and thus an average temperature is assumed to account for an ensemble of grains with different radii. We investigate the impact of subdividing the grain radius distribution into multiple sub-intervals on the resulting dust temperature distribution and spectral energy distribution (SED). These quantities were computed for two different scenarios: (1) Radius distribution represented by 16 logarithmically distributed radius intervals, and (2) radius distribution represented by a single grain species with averaged optical properties (reference). Within the considered parameter range, i.e., of grain radii between 5 nm and 1 mm and an optically thin and thick disk with a parameterized density distribution, we obtain the following results: In optically thin disk regions, the temperature spread can be as large as ~63% and the relative grain surface below a certain temperature is lower than in the reference disk. With increasing optical depth, the difference in the midplane temperature and the relative grain surface below a certain temperature decreases. Furthermore, below ~20K, this fraction is higher for the reference disk than for the case of multiple grain radii, while it shows the opposite behavior for temperatures above this threshold. The thermal emission in the case of multiple grain radii at short wavelengths is stronger than for the reference disk. The freeze-out radius is a function of grain radius, spanning a radial range between the coldest and warmest grain species of ~30AU.

We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.

We present ALMA Cycle~2 observations at $0.5^{\prime\prime}$ resolution of TW Hya of CS $J=5-4$ emission. The radial profile of the integrated line emission displays oscillatory features outwards of $1.5^{\prime\prime}$ ($\approx 90$ au). A dip-like feature at $1.6^{\prime\prime}$ is coincident in location, depth and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modelling and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modelling performed for scattered light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at $3.1^{\prime\prime}$ is identified as a region of intensified desorption due to enhanced penetration of the interstellar FUV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices.

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection.

The increased sensitivity of millimeter-wave facilities now makes possible the detection of low amounts of gas in debris disks. Some of the gas-rich debris disks harbor peculiar properties, with possible pristine gas and secondary generated dust. The origin of the gas in these hybrid disks is strongly debated and the current sample is too sparse to understand this phenomenon. More detections are necessary to increase the statistics on this population. Lying at the final stages of evolution of proto-planetary disks and at the beginning of the debris disk phase, these objects could provide new insight into the processes involved in the making of planetary systems. We carried out a deep survey of the 12CO(2-1) and 12CO(3-2) lines with the APEX and IRAM radiotelescopes in young debris disks selected according to hybrid disk properties. The survey is complemented with a bibliographic study of the ratio between the emission of the gas and the continuum (S_CO/F_cont) in CTTS, Herbig Ae, WTTS, hybrid, and debris disks. Our sub-mm survey comprises 25 stars, including 17 new targets, and we increase the sensitivity limit by a factor 2 on eight sources compared to similar published studies. We report a 4sigma tentative detection of a double-peaked 12CO(2-1) line around HD23642; an eclipsing binary located in the Pleiades. We also reveal a correlation between the emission of the CO gas and the dust continuum from CTTS, Herbig Ae and few debris disks. The observed trend of the gas to dust flux ratio suggests a concurrent dissipation of the dust and gas components. Hybrid disks systematically lie above this trend, suggesting that these systems may witness a transient phase, when the dust has evolved more rapidly than the gas, with a flux ratio S_CO/F_cont enhanced by a factor of between 10 and 100 compared to standard (proto-)planetary disks.

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN.

We conducted IRAM-30m C18O (2-1) and SMA 1.3mm continuum, 12CO (2-1), and C18O (2-1) observations toward the Class 0/I protostar L1455 IRS1 in Perseus. The IRAM results show L1455 IRS1 located in a dense core of 0.05 pc in size with a mass of 0.54M_sun, connecting to a filamentary structure. Inside the dense core, compact components of 350 AU and 1500 AU are detected in the SMA 1.3 mm continuum and C18O, with a velocity gradient in the latter one perpendicular to a bipolar outflow in 12CO, likely tracing a rotational motion. We measure a rotational velocity profile ~ r^-0.75 that becomes shallower at a turning radius of ~200 AU which is approximately the radius of the 1.3 mm continuum component. These results hint the presence of a Keplerian disk with a radius <200 AU around L1455 IRS1 with a protostellar mass of about 0.28 M_sun. We derive a core rotation that is about one order of magnitude faster than expected. A significant velocity gradient along a filament towards IRS1 indicates that this filament is dynamically important, providing a gas reservoir and possibly responsible for the faster-than-average core rotation. Previous polarimetric observations show a magnetic field aligned with the outflow axis and perpendicular to the associated filament on a 0.1~pc scale, while on the inner 1000 AU scale, the field becomes perpendicular to the outflow axis.This change in magnetic field orientations is consistent with our estimated increase in rotational energy from large to small scales that overcomes the magnetic field energy, wrapping the field lines and aligning them with the disk velocity gradient. These results are discussed in the context of the interplay between filament, magnetic field, and gas kinematics from large to small scale. Possible emerging trends are explored with a sample of 8 Class 0/I protostars.

We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

(abridged) We used the IRAM 30-m to perform a sensitive wideband survey of 30 protoplanetary disks in the Taurus Auriga region. We simultaneously observed HCO$^+$(3-2), HCN(3-2), C$_2$H(3-2), CS(5-4), and two transitions of SO. We combine the results with a previous survey which observed $^{13}$CO (2-1), CN(2-1), two o-H$_2$CO lines and one of SO. We use available interferometric data to derive excitation temperatures of CN and C$_2$H in several sources. We determine characteristic sizes of the gas disks and column densities of all molecules using a parametric power-law disk model. Our study is mostly sensitive to molecules at 200-400 au from the stars. We compare the derived column densities to the predictions of an extensive gas-grain chemical disk model, under conditions representative of T Tauri disks. This survey provides 20 new detections of HCO$^+$ in disks, 18 in HCN, 11 in C$_2$H, 8 in CS and 4 in SO. HCO$^+$ is detected in almost all sources, and its J=3-2 line is essentially optically thick, providing good estimates of the disk radii. The other transitions are (at least partially) optically thin. Variations of the column density ratios do not correlate with any specific property of the star or disk. Disks around Herbig Ae stars appear less rich in molecules than those around T Tauri stars, although the sample remains small. SO is only found in the (presumably younger) embedded objects, perhaps reflecting an evolution of the S chemistry due to increasing depletion with time. Overall, the molecular column densities, and in particular the CN/HCN and CN/C$_2$H ratios, are well reproduced by gas-grain chemistry in cold disks. This study provides a census of simple molecules in disks of radii $> 200-300$ au. Extending that to smaller disks, or searching for less abundant or more complex molecules requires a much more sensitive facility, i.e. NOEMA and ALMA.

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed.

We aim to unveil the observational imprint of physical mechanisms that govern planetary formation in the young, multiple system GG Tau A. We present ALMA observations of $^{12}$CO and $^{13}$CO 3-2 and 0.9 mm continuum emission with 0.35" resolution. The $^{12}$CO 3-2 emission, found within the cavity of the circumternary dust ring (at radius $< 180$ AU) where no $^{13}$CO emission is detected, confirms the presence of CO gas near the circumstellar disk of GG Tau Aa. The outer disk and the recently detected hot spot lying at the outer edge of the dust ring are mapped both in $^{12}$CO and $^{13}$CO. The gas emission in the outer disk can be radially decomposed as a series of slightly overlapping Gaussian rings, suggesting the presence of unresolved gaps or dips. The dip closest to the disk center lies at a radius very close to the hot spot location at $\sim250-260$~AU. The CO excitation conditions indicate that the outer disk remains in the shadow of the ring. The hot spot probably results from local heating processes. The two latter points reinforce the hypothesis that the hot spot is created by an embedded proto-planet shepherding the outer disk.

Dust determines the temperature structure of protoplanetary disks. However, dust temperature determinations almost invariably rely on a complex modeling of the Spectral Energy Distribution. We attempt a direct determination of the temperature of large grains emitting at mm wavelengths.} We observe the edge-on dust disk of the Flying Saucer, which appears in silhouette against the CO J=2-1 emission from a background molecular cloud in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity allows us to directly measure the %absorbing dust temperature. The dust opacity can then be derived from the emitted continuum radiation. The dust disk absorbs the radiation from the CO clouds at several velocities. We derive very low dust temperatures, 5 to 7 K at radii around 100 au, which is much lower than most model predictions. The dust optical depth is $> 0.2$ at 230 GHz, and the scale height at 100 au is at least 8 au (best fit 13 au). However, the dust disk is very flat (flaring index -0.35), which is indicative of dust settling in the outer parts.

Debris disks are usually thought to be gas-poor, the gas being dissipated by accretion or evaporation during the protoplanetary phase. HD141569A is a 5 Myr old star harboring a famous debris disk, with multiple rings and spiral features. We present here the first PdBI maps of the 12CO(2-1), 13CO(2-1) gas and dust emission at 1.3 mm in this disk. The analysis reveals there is still a large amount of (primordial) gas extending out to 250 au, i. e. inside the rings observed in scattered light. HD141569A is thus a hybrid disk with a huge debris component, where dust has evolved and is produced by collisions, with a large remnant reservoir of gas.

The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks and the formation of planetary system. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre-stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model representative of T Tauri stars. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less volatile forms (principally s-CO$_2$, and to a lesser extent s-CH$_4$) instead of being evaporated over a wide range of temperature. The canonical gas-phase abundance of 10$^{-4}$ is only reached above about 30-35 K. The dominant Carbon bearing entity depends on the temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content of disks.

We attempt to determine the molecular composition of disks around young low-mass stars in the $\rho$ Oph region and to compare our results with a similar study performed in the Taurus-Auriga region. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-1 in 29 T Tauri stars located in the $\rho$ Oph and upper Scorpius regions. $^{13}$CO J=2-1 is observed simultaneously to provide an indication of the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO, and the C$^{17}$O J=2-1 line, which provides complementary information on the nature of the emission. Contamination by molecular cloud in $^{13}$CO and even C$^{17}$O is ubiquitous. The CN detection rate appears to be lower than for the Taurus region, with only four sources being detected (three are attributable to disks). H$_2$CO emission is found more frequently, but appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in the $\rho$ Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures in $\rho$ Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile forms. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at high enough resolution with ALMA are required to probe the gas disk population.

We present resolved Plateau de Bure Array observations of DM Tau in lines of HCO+ (3-2), (1-0) and DCO+ (3-2). A power-law fitting approach allowed a derivation of column densities of these two molecules. A chemical inner hole of ~50 AU was found in both HCO+ and DCO+ with DCO+ emission extending to only 450 AU. An isotopic ratio of R_D = N(DCO+) / N(HCO+) was found to range from 0.1 at 50 AU and 0.2 at 450 AU. Chemical modeling allowed an exploration of the sensitivity of these molecular abundances to physical parameters out with temperature, finding that X-rays were the domination ionization source in the HCO+ molecular region and that R_D also is sensitive to the CO depletion. The ionization fraction, assuming a steady state system, was found to be x(e-) ~ 10$^{-7}$. Modeling suggests that HCO+ is the dominant charged molecule in the disk but its contribution to ionization fraction is dwarfed by atmoic ions such as C+, S+ and H+.

Recent exo-planetary surveys reveal that planets can orbit and survive around binary stars. This suggests that some fraction of young binary systems which possess massive circumbinary disks (CB) may be in the midst of planet formation. However, there are very few CB disks detected. We revisit one of the known CB disks, the UY Aurigae system, and probe 13CO 2-1, C18O 2-1, SO 5(6)-4(5) and 12CO 3-2 line emission and the thermal dust continuum. Our new results confirm the existence of the CB disk. In addition, the circumstellar (CS) disks are clearly resolved in dust continuum at 1.4 mm. The spectral indices between the wavelengths of 0.85 mm and 6 cm are found to be surprisingly low, being 1.6 for both CS disks. The deprojected separation of the binary is 1.26" based on our 1.4 mm continuum data. This is 0.07" (10 AU) larger than in earlier studies. Combining the fact of the variation of UY Aur B in $R$ band, we propose that the CS disk of an undetected companion UY Aur Bb obscures UY Aur Ba. A very complex kinematical pattern inside the CB disk is observed due to a mixing of Keplerian rotation of the CB disk, the infall and outflow gas. The streaming gas accreting from the CB ring toward the CS disks and possible outflows are also identified and resolved. The SO emission is found to be at the bases of the streaming shocks. Our results suggest that the UY Aur system is undergoing an active accretion phase from the CB disk to the CS disks. The UY Aur B might also be a binary system, making the UY Aur a triple system.

Aims: We attempt to determine the masses of single or multiple young T Tauri and HAeBe stars from the rotation of their Keplerian disks. Methods:We used the IRAM PdBI interferometer to perform arcsecond resolution images of the CN N=2-1 transition with good spectral resolution. Integrated spectra from the 30-m radiotelescope show that CN is relatively unaffected by contamination from the molecular clouds. Our sample includes 12 sources, among which isolated stars like DM Tau and MWC 480 are used to demonstrate the method and its accuracy. We derive the dynamical mass by fitting a disk model to the emission, a process giving M/D the mass to distance ratio. We also compare the CN results with higher resolution CO data, that are however affected by contamination. Results: All disks are found in nearly perfect Keplerian rotation. We determine accurate masses for 11 stars, in the mass range 0.5 to 1.9 solar masses. The remaining one, DG Tau B, is a deeply embedded object for which CN emission partially arises from the outflow. With previous determination, this leads to 14 (single) stars with dynamical masses. Comparison with evolutionary tracks, in a distance independent modified HR diagram, show good overall agreement (with one exception, CW Tau), and indicate that measurement of effective temperatures are the limiting factor. The lack of low mass stars in the sample does not allow to distinguish between alternate tracks.

We aim at unveiling the observational imprint of physical mechanisms that govern planetary formation in young, multiple systems. In particular, we investigate the impact of tidal truncation on the inner circumstellar disks. We observed the emblematic system GG Tau at high-angular resolution: a hierarchical quadruple system composed of low-mass T Tauri binary stars surrounded by a well-studied, massive circumbinary disk in Keplerian rotation. We used the near-IR 4-telescope combiner PIONIER on the VLTI and sparse-aperture-masking techniques on VLT/NaCo to probe this proto-planetary system at sub-AU scales. We report the discovery of a significant closure-phase signal in H and Ks bands that can be reproduced with an additional low-mass companion orbiting GG Tau Ab, at a (projected) separation rho = 31.7 +/- 0.2mas (4.4 au) and PA = 219.6 +/- 0.3deg. This finding offers a simple explanation for several key questions in this system, including the missing-stellar-mass problem and the asymmetry of continuum emission from the inner dust disks observed at millimeter wavelengths. Composed of now five co-eval stars with 0.02 <= Mstar <= 0.7 Msun, the quintuple system GG Tau has become an ideal test case to constrain stellar evolution models at young ages (few 10^6yr).

Grain-surface reactions play an essential role in interstellar chemistry, since dust grain catalyses reactions at its surface allowing for the formation of molecules. We used a chemical model in which both gas-phase and grain-surface reactions occur and studied particularly the diffusion mechanisms on the surface of the grains. Surface reactions can occur via thermal hopping when species cross over a potential barrier or via quantum tunneling when species cross through this barrier. We show that the thermal diffusion (hopping) can be much more efficient after a cosmic ray particle collides with a dust grain, heating it to a peak temperature of 70K. We present here the results of numerical simulations after including the quantum tunneling mechanism for species H, H2 and O and considering the effect of cosmic ray particle collision on the surface reactions. As a consequence, the gas-phase and grain-surface abundances are affected and we show that more complex molecules can be formed in molecular clouds.

(abridged) Most Class II sources (of nearby star forming regions) are surrounded by disks with weak millimeter continuum emission. These "faint" disks may hold clues to the disk dissipation mechanism. We attempt to determine the characteristics of such faint disks around classical T Tauri stars, and to explore the link between disk faintness and the proposed disk dispersal mechanisms (accretion, viscous spreading, photo-evaporation, planetary system formation). We performed high-angular resolution (0.3") imaging of a small sample of disks (9 sources) with low 1.3mm continuum flux (mostly <30 mJy) with the IRAM Plateau de Bure interferometer and simultaneously searched for 13CO (or CO) J=2-1 line emission. Using a simple parametric disk model, we determine characteristic sizes of the disks, in dust and gas, and we constrain surface densities in the central 50 AU. All disks are much smaller than the bright disks imaged so far, both in continuum and 13CO lines (5 detections). In continuum, half of the disks are very small, with characteristic radii less than 10AU, but still have high surface density values. Small sizes appear to be the main cause for the low disk luminosity. Direct evidence for grain growth is found for the three disks that are sufficiently resolved. Low continuum opacity is attested in two systems only, but we cannot firmly distinguish between a low gas surface density and a lower dust emissivity resulting from grain growth. We report a tentative discovery of a 20 AU radius cavity in DS Tau, bringing the proportion of transitional disks to a similar value to that of brighter sources, but cavities cannot explain the low mm flux. This study highlights a category of very compact dust disks, still exhibiting high surface densities, which may represent up to 25 % of the whole disk population, but its origin is unclear with the current data alone.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

An outstanding question of astrobiology is the link between the chemical composition of planets, comets, and other solar system bodies and the molecules formed in the interstellar medium. Understanding the chemical and physical evolution of the matter leading to the formation of protoplanetary disks is an important step for this. We provide some new clues to this long-standing problem using three-dimensional chemical simulations of the early phases of disk formation: we interfaced the full gas-grain chemical model Nautilus with the radiation-magnetohydrodynamic model RAMSES, for different configurations and intensities of magnetic field. Our results show that the chemical content (gas and ices) is globally conserved during the collapsing process, from the parent molecular cloud to the young disk surrounding the first Larson core. A qualitative comparison with cometary composition suggests that comets are constituted of different phases, some molecules being direct tracers of interstellar chemistry, while others, including complex molecules, seem to have been formed in disks, where higher densities and temperatures allow for an active grain surface chemistry. The latter phase, and its connection with the formation of the first Larson core, remains to be modeled.

Context: Circumstellar disks are considered to be the environment for the formation of planets. The growth of dust grains in these disks is the first step in the core accretion-gas capture planet formation scenario. Indicators and evidence of disk evolution can be traced in spatially resolved images and the spectral energy distribution (SED) of these objects. Aims: We develop a model for the dust phase of the edge-on oriented circumstellar disk of the Butterfly Star which allows one to fit observed multi-wavelength images and the SED simultaneously. Methods: Our model is based on spatially resolved high angular resolution observations at 1.3 mm, 894 micron, 2.07 micron, 1.87 micron, 1.60 micron, and 1.13 micron and an extensively covered SED ranging from 12 micron to 2.7 mm, including a detailed spectrum obtained with the Spitzer Space Telescope in the range from 12 micron to 38 micron. A parameter study based on a grid search method involving the detailed analysis of every parameter was performed to constrain the disk parameters and find the best-fit model for the independent observations. The individual observations were modeled simultaneously, using our continuum radiative transfer code. Results: We derived a model that is capable of reproducing all of the observations of the disk at the same time. We find quantitative evidence for grain growth up to ~100 micron-sized particles, vertical settling of larger dust grains toward the disk midplane, and radial segregation of the latter toward the central star. Within our best-fit model the large grains have a distribution with a scale height of 3.7 AU at 100 AU and a radial extent of 175 AU compared to a hydrostatic scale height of 6.9 AU at 100 AU and an outer disk radius of 300 AU. Our results are consistent with current theoretical models for the evolution of circumstellar disks and the early stages of planet formation.

Investigating the dynamical evolution of dust grains in proto-planetary disks is a key issue to understand how planets should form. We identify under which conditions dust settling can be constrained by high angular resolution observations at mm wavelengths, and which observational strategies are suited for such studies. Exploring a large range of models, we generate synthetic images of disks with different degrees of dust settling, and simulate high angular resolution (~ 0.05-0.3") ALMA observations of these synthetic disks. The resulting data sets are then analyzed blindly with homogeneous disk models (where dust and gas are totally mixed) and the derived disk parameters are used as tracers of the settling factor. Our dust disks are partially resolved by ALMA and present some specific behaviors on radial and mainly vertical directions, which can be used to quantify the level of settling. We find out that an angular resolution better than or equal to ~ 0.1" (using 2.3 km baselines at 0.8mm) allows us to constrain the dust scale height and flaring index with sufficient precision to unambiguously distinguish between settled and non-settled disks, provided the inclination is close enough to edge-on (i >= 75{\deg}). Ignoring dust settling and assuming hydrostatic equilibrium when analyzing such disks affects the derived dust temperature and the radial dependency of the dust emissivity index. The surface density distribution can also be severely biased at the highest inclinations. However, the derived dust properties remain largely unaffected if the disk scale height is fitted separately. ALMA has the potential to test some of the dust settling mechanisms, but for real disks, deviations from ideal geometry (warps, spiral waves) may provide an ultimate limit on the dust settling detection.

We present the first high-angular resolution (up to 0.7", ~5000 AU) polarization and thermal dust continuum images toward the massive star-forming region W51 North. The observations were carried out with the Submillimeter Array (SMA) in both the subcompact (SMA-SubC) and extended (SMA-Ext) configurations at a wavelength of 870 micron. W51 North is resolved into four cores (SMA1 to SMA4) in the 870 micron continuum image. The associated dust polarization exhibits more complex structures than seen at lower angular resolution. We analyze the inferred morphologies of the plane-of-sky magnetic field (B_bot) in the SMA1 to SMA4 cores and in the envelope using the SMA-Ext and SMA-SubC data. These results are compared with the B_bot archive images obtained from the CSO and JCMT. A correlation between dust intensity gradient position angles (phi_{nabla I}) and magnetic field position angles (phi_B) is found in the CSO, JCMT and both SMA data sets. This correlation is further analyzed quantitatively. A systematically tighter correlation between phi_{nabla I} and phi_B is found in the cores, whereas the correlation decreases in outside-core regions. Magnetic field-to-gravity force ratio (Sigma_B) maps are derived using the newly developed polarization - intensity gradient method by Koch, Tang & Ho 2012. We find that the force ratios tend to be small (Sigma_B <= 0.5) in the cores in all 4 data sets. In regions outside of the cores, the ratios increase or the field is even dominating gravity (Sigma_B > 1). This possibly provides a physical explanation of the tightening correlation between phi_{nabla I} and phi_B in the cores: the more the B field lines are dragged and aligned by gravity, the tighter the correlation is. Finally, we propose a schematic scenario for the magnetic field in W51 North to interpret the four polarization observations at different physical scales.

Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the transport of angular momentum, and the growth of dust particles. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1 ") and spectral resolution (0.126 km/s) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the disk. A low sensitivity CS J=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic linewidth derived from the CS J=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this layer is located near one scale height, in reasonable agreement with chemical model predictions.

We use the IRAM 30-m telescope to perform a sensitive search for CN N=2-1 in 42 T Tauri or Herbig Ae systems located mostly in the Taurus-Auriga region. $^{13}$CO J=2-1 is observed simultaneously to indicate the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO and the C$^{17}$O J=2-1 line which provide complementary information on the nature of the emission. While $^{13}$CO is in general dominated by residual emission from the cloud, CN exhibits a high disk detection rate $> 50$% in our sample. We even report CN detection in stars for which interferometric searches failed to detect $^{12}$CO, presumably because of obscuration by a foreground, optically thick, cloud. Comparison between CN and o-H$_2$CO or SO line profiles and intensities divide the sample in two main categories. Sources with SO emission are bright and have strong H$_2$CO emission, leading in general to [H$_2$CO/CN]$ > 0.5$. Furthermore, their line profiles, combined with a priori information on the objects, suggest that the emission is coming from outflows or envelopes rather than from a circumstellar disk. On the other hand, most sources have [H$_2$CO/CN]$ < 0.3$, no SO emission, and some of them exhibit clear double-peaked profiles characteristics of rotating disks. In this second category, CN is likely tracing the proto-planetary disks. From the line flux and opacity derived from the hyperfine ratios, we constrain the outer radii of the disks, which range from 300 to 600 AU. The overall gas disk detection rate (including all molecular tracers) is $\sim 68%$, and decreases for fainter continuum sources. This study shows that gas disks, like dust disks, are ubiquitous around young PMS stars in regions of isolated star formation, and that a large fraction of them have $R > 300$ AU.

The circumstellar disk of AB Aurigae has garnered strong attention owing to the apparent existence of spirals at a relatively young stage and also the asymmetric disk traced in thermal dust emission. However, the physical conditions of the spirals are still not well understood. The origin of the asymmetric thermal emission is unclear. We observed the disk at 230 GHz (1.3 mm) in both the continuum and the spectral line ^12CO J=2-1 with IRAM 30-m, the Plateau de Bure interferometer, and the Submillimeter Array to sample all spatial scales from 0.37" to about 50". To combine the data obtained from these telescopes, several methods and calibration issues were checked and discussed. The 1.3 mm continuum (dust) emission is resolved into inner disk and outer ring. Molecular gas at high velocities traced by the CO line is detected next to the stellar location. The inclination angle of the disk is found to decrease toward the center. On a larger scale, based on the intensity weighted dispersion and the integrated intensity map of ^12CO J=2-1, four spirals are identified, where two of them are also detected in the near infrared. The total gas mass of the 4 spirals (M_spiral) is 10^-7 < M_spiral < 10^-5 M_sun, which is 3 orders of magnitude smaller than the mass of the gas ring. Surprisingly, the CO gas inside the spiral is apparently counter-rotating with respect to the CO disk, and it only exhibits small radial motion. The wide gap, the warped disk, and the asymmetric dust ring suggest that there is an undetected companion with a mass of 0.03 M_sun at a radius of 45 AU. Although an hypothetical fly-by cannot be ruled out, the most likely explanation of the AB Aurigae system may be inhomogeneous accretion well above or below the main disk plane from the remnant envelope, which can explain both the rotation and large-scale motions detected with the 30-m image.

The radio-quiet quasar BR1202-0725 (z=4.695) is a remarkable source with a bright Northwest (NW) companion detected at submm and radio wavelengths but invisible in the optical. In the absence of amplification by gravitational lensing, BR1202-0725 would be the most luminous binary CO and FIR source in the Universe. In this paper, we report observations with the IRAM Plateau de Bure interferometer of BR1202-0725 in the redshifted emission of the CO(5-4) and (7-6) lines, the [C I](3P2-3P1) line, a high angular resolution (0.3 x 0.8 arcsec) 1.3 mm map of the rest-frame, far-IR dust continuum, and a search for the CO(11-10) line. We compare these results with recent ALMA data in the [C II] line. Both the quasar host galaxy and its NW companion are spatially resolved in the molecular line emission and the dust continuum. The CO profile of the NW companion is very broad with a full width at half maximum of 1000 +/- 130 km/s, compared to 360 +/- 40 km/s for the quasar host galaxy to the Southeast (SE). The difference in linewidths and center velocities, and the absence of any lens candidate or arc-like structure in the field, at any wavelength, show that the obscured NW galaxy and the SE quasar host galaxy cannot be lensed images of the same object. Instead, we find morphological and kinematic evidence for sub-structures in both the NW and SE sources. We interpret these results as strong indications that the BR1202-0725 complex is a group of young, interacting, and highly active starburst galaxies.

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCa 15 and the Herbig Ae star MWC 480. We report the first detection of HC3N at 5 sigma in the GO Tau and MWC 480 disks with the IRAM 30-m, and in the LkCa 15 disk (5 sigma), using the IRAM array, with derived column densities of the order of 10^{12}cm^{-2}. We also obtain stringent upper limits on CCS (N < 1.5 x 10^{12} cm^{-3}). We discuss the observational results by comparing them to column densities derived from existing chemical disk models (computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule observations. The observed column densities of HC3N are typically two orders of magnitude lower than the existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO).

Circumstellar disks are characteristic for star formation and vanish during the first few Myr of stellar evolution. During this time planets are believed to form in the dense midplane by growth, sedimentation and aggregation of dust. Indicators of disk evolution, such as holes and gaps, can be traced in the spectral energy distribution (SED) and spatially resolved images. We aim to construct a self-consistent model of HH 30 by fitting all available continuum observations simultaneously. New data sets not available in previous studies, such as high-resolution interferometric imaging with the Plateau de Bure Interferometer (PdBI) at lambda = 1.3 mm and SED measured with IRS on the Spitzer Space Telescope in the mid-infrared, put strong constraints on predictions and are likely to provide new insights into the evolutionary state of this object. A parameter study based on simulated annealing was performed to find unbiased best-fit models for independent observations made in the wavelength domain lambda ~ 1 micron ... 4 mm. The method essentially creates a Markov chain through parameter space by comparing predictions generated by our self-consistent continuum radiation transfer code MC3D with observations. We present models of the edge-on circumstellar disk of HH 30 based on observations from the near-infrared to mm-wavelengths that suggest the presence of an inner depletion zone with about 45 AU radius and a steep decline of mm opacity beyond 140 AU. Our modeling indicates that several modes of dust evolution such as growth, settling, and radial migration are taking place in this object. High-resolution observations of HH 30 at different wavelengths with next-generation observatories such as ALMA and JWST will enable the modeling of inhomogeneous dust properties and significantly expand our understanding of circumstellar disk evolution.

We present high spatial resolution maps of ro-vibrational molecular hydrogen emission from the environment of the GG Tau A binary component in the GG Tau quadruple system. The H2 v= 1-0 S(1) emission is spatially resolved and encompasses the inner binary, with emission detected at locations that should be dynamically cleared on several hundred-year timescales. Extensions of H2 gas emission are seen to ~100 AU distances from the central stars. The v = 2-1 S(1) emission at 2.24 microns is also detected at ~30 AU from the central stars, with a line ratio of 0.05 +/- 0.01 with respect to the v = 1-0 S(1) emission. Assuming gas in LTE, this ratio corresponds to an emission environment at ~1700 K. We estimate that this temperature is too high for quiescent gas heated by X-ray or UV emission from the central stars. Surprisingly, we find that the brightest region of H2 emission arises from a spatial location that is exactly coincident with a recently revealed dust "streamer" which seems to be transferring material from the outer circumbinary ring around GG Tau A into the inner region. As a result, we identify a new excitation mechanism for ro-vibrational H2 stimulation in the environment of young stars. The H2 in the GG Tau A system appears to be stimulated by mass accretion infall as material in the circumbinary ring accretes onto the system to replenish the inner circumstellar disks. We postulate that H2 stimulated by accretion infall could be present in other systems, particularly binaries and "transition disk" systems which have dust cleared gaps in their circumstellar environments.

We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the $1_{10}-1_{01}$ line of H$_2$S at 168.8 GHz and $2_{23}-1_{12}$ line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. We fail to detect the SO and H$_2$S lines, although CS is detected in LkCa15, DM\,Tau, and GO\,Tau but not in MWC\,480. However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the sulfur chemistry. Our best modeling of disks is obtained for a C/O ratio of 1.2, starting from initial cloud conditions of H density of $2\times 10^5$ cm$^{-3}$ and age of $10^6$ yr. The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce the H$_2$S upper limits. The predicted H$_2$S column densities are too high by at least one order of magnitude. H$_2$S may remain locked onto grain surfaces and react with other species, thereby preventing the desorption of H$_2$S.

The chemistry of proto-planetary disks is thought to be dominated by two major processes: photodissociation near the disk surface, and depletion on dust grains in the disk mid-plane, resulting in a layered structure with molecules located in a warm layer above the disk mid-plane. We attempt here to confront this warm molecular layer model prediction with the distribution of two key molecules for dissociation processes: CN and HCN. Using the IRAM Plateau de Bure interferometer, we obtained high spatial and spectral resolution images of the CN J=2-1 and HCN J=1-0 lines in the disks surrounding the two T-Tauri DM Tau and LkCa 15 and the Herbig Ae MWC 480. Disk properties are derived assuming power law distributions. The hyperfine structure of the observed transitions allows us to constrain the line opacities and excitation temperatures. We compare the observational results with predictions from existing chemical models, and use a simple PDR model (without freeze-out of molecules on grains and surface chemistry) to illustrate dependencies on UV field strength, grain size and gas-to-dust ratio. We also evaluate the impact of Lyman alpha radiation. The temperature ordering follows the trend found from CO lines, with DM Tau being the coldest object and MWC 480 the warmest. Although CN indicates somewhat higher excitation temperatures than HCN, the derived values in the T-Tauri disks are very low (8-10 K). They agree with results obtained from CCH, and are in contradiction with thermal and chemical model predictions. These very low temperatures, as well as geometrical constraints, suggest that substantial amounts of CN and HCN remain in the gas phase close to the disk mid-plane, and that this mid-plane is quite cold. The observed CN/HCN ratio (5-10) is in better agreement with the existence of large grains, and possibly also a substantial contribution of Lyman alpha radiation.

The structure in density and temperature of protoplanetary disks surrounding low-mass stars is not yet well known. The protoplanetary disks mid-planes are expected to be very cold and thus depleted in molecules in gas phase, especially CO. Recent observations of molecules at very low apparent temperature (~ 6 K) challenge this current picture of the protoplanetary disk structures. We aim at constraining the physical conditions, and in particular the gas-phase CO abundance in the mid-plane of protoplanetary disks. The light molecule H2D+,is a tracer of cold and CO-depleted environment. It is therefore a good candidate to explore the disks mid-planes. We performed a deep search for H2D+ in the two well-known disks surrounding TW Hya and DM Tau using the APEX and JCMT telescopes. The analysis of the observations are done with DISKFIT, a radiative transfer code dedicated to disks. In addition, we used a chemical model describing deuterium chemistry to infer the implications of our observations on the level of CO depletion and on the ionization rate in the disks mid-plane. The ortho-H2D+(1(1,0))-1(1,1)) line at 372 GHz was not detected. Although our limit is three times better than previous observations, comparison with the chemical modeling indicates that it remains insufficient to put valuable constraints on the CO abundance in the disk mid-plane. Even with ALMA, the detection of H2D+ may not be straightforward, and H2D+ may not be a sufficiently sensitive tracer of the protoplanetary disks mid-plane

Previous studies have found that Vega is surrounded by an extended debris disc that is very smooth in the far infrared, but displays possible clumpiness at 850micron and dust emission peaks at 1.3mm. We reobserved Vega at 1.3mm with PdBI to constrain its circumstellar dust distribution. Our observations of a three-field mosaic have a factor of two higher sensitivity than previous observations. We detect Vega photosphere with the expected flux, but none of the previously reported emission peaks that should have been detected at the >6sigma level, with a sensitivity <1mK. This implies that the dust distribution around Vega is principally smooth and circularly symmetric. This also means that no planet is needed to account for dust trapped in mean-motion resonnance.

Context: Dark cloud chemical models usually predict large amounts of O2, often above observational limits. Aims: We investigate the reason for this discrepancy from a theoretical point of view, inspired by the studies of Jenkins and Whittet on oxygen depletion. Methods: We use the gas-grain code Nautilus with an up-to-date gas-phase network to study the sensitivity of the molecular oxygen abundance to the oxygen elemental abundance. We use the rate coefficient for the reaction O + OH at 10 K recommended by the KIDA (KInetic Database for Astrochemistry) experts. Results: The updates of rate coefficients and branching ratios of the reactions of our gas-phase chemical network, especially N + CN and H3+ + O, have changed the model sensitivity to the oxygen elemental abundance. In addition, the gas-phase abundances calculated with our gas-grain model are less sensitive to the elemental C/O ratio than those computed with a pure gas-phase model. The grain surface chemistry plays the role of a buffer absorbing most of the extra carbon. Finally, to reproduce the low abundance of molecular oxygen observed in dark clouds at all times, we need an oxygen elemental abundance smaller than 1.6E-4. Conclusion: The chemistry of molecular oxygen in dense clouds is quite sensitive to model parameters that are not necessarily well constrained. That O2 abundance may be sensitive to nitrogen chemistry is an indication of the complexity of interstellar chemistry.

(Abridged) We attempt to characterize the radial distribution of dust in disks around a sample of young stars from an observational point of view, and, when possible, in a model-independent way, by using parametric laws. We used the IRAM PdBI interferometer to provide very high angular resolution (down to 0.4" in some sources) observations of the continuum at 1.3 mm and 3 mm around a sample of T Tauri stars in the Taurus-Auriga region. The sample includes single and multiple systems, with a total of 23 individual disks. We used track-sharing observing mode to minimize the biases. We fitted these data with two kinds of models: a "truncated power law" model and a model presenting an exponential decay at the disk edge ("viscous" model). ect evidence for tidal truncation is found in the multiple systems. The temperature of the mm-emitting dust is constrained in a few systems. Unambiguous evidence for large grains is obtained by resolving out disks with very low values of the dust emissivity index Beta. In most disks that are sufficiently resolved at two different wavelengths, we find a radial dependence of Beta, which appears to increase from low values (as low as 0) at the center to about 1.7 -- 2 at the disk edge. The same behavior could apply to all studied disks. It introduces further ambiguities in interpreting the brightness profile, because the regions with apparent Beta = 0 can also be interpreted as being optically thick when their brightness temperature is high enough. Despite the added uncertainty on the dust absorption coefficient, the characteristic size of the disk appears to increase with a higher estimated star age. These results provide the first direct evidence of the radial dependence of the grain size in proto-planetary disks.

Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.html

Gas and dust dissipation processes of proto-planetary disks are hardly known. Transition disks between Class II (proto-planetary disks) and Class III (debris disks) remain difficult to detect. We investigate the carbon chemistry of the peculiar CQ Tau gas disk. It is likely a transition disk because it exhibits weak CO emission with a relatively strong millimeter continuum, indicating that the disk might be currently dissipating its gas content. We used APEX to observe the two CI lines at 492GHz and 809 GHz in the disk orbiting CQ Tau. We compare the observations to several chemical model predictions. We focus our study on the influence of the stellar UV radiation shape and gas-to-dust ratio. We did not detect the CI lines. However, our upper limits are deep enough to exclude high-CI models. The only available models compatible with our limits imply very low gas-to-dust ratio, of the order of a few, only. These observations strengthen the hypothesis that CQ Tau is likely a transition disk and suggest that gas disappears before dust.

Recent sub-millimetric observations at the Plateau de Bure interferometer evidenced a cavity at ~ 46 AU in radius into the proto-planetary disk around the T Tauri star LkCa15 (V1079 Tau), located in the Taurus molecular cloud. Additional Spitzer observations have corroborated this result possibly explained by the presence of a massive (>= 5 MJup) planetary mass, a brown dwarf or a low mass star companion at about 30 AU from the star. We used the most recent developments of high angular resolution and high contrast imaging to search directly for the existence of this putative companion, and to bring new constraints on its physical and orbital properties. The NACO adaptive optics instrument at VLT was used to observe LkCa15 using a four quadrant phase mask coronagraph to access small angular separations at relatively high contrast. A reference star at the same parallactic angle was carefully observed to optimize the quasi-static speckles subtraction (limiting our sensitivity at less than 1.0). Although we do not report any positive detection of a faint companion that would be responsible for the observed gap in LkCa15's disk (25-30 AU), our detection limits start constraining its probable mass, semi-major axis and eccentricity. Using evolutionary model predictions, Monte Carlo simulations exclude the presence of low eccentric companions with masses M >= 6 M Jup and orbiting at a >= 100 AU with significant level of confidence. For closer orbits, brown dwarf companions can be rejected with a detection probability of 90% down to 80 AU (at 80% down to 60 AU). Our detection limits do not access the star environment close enough to fully exclude the presence of a brown dwarf or a massive planet within the disk inner activity (i.e at less than 30 AU). Only, further and higher contrast observations should unveil the existence of this putative companion inside the LkCa15 disk.

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.

We studied several representative circumstellar disks surrounding the Herbig Ae star MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)millimeter wavelengths in lines of CCH. Our aim is to characterize photochemistry in the heavily UV-irradiated MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative minimization technique, the CCH column densities and excitation conditions are con- strained. Very low excitation temperatures are derived for the T Tauri stars. These values are compared with the results of advanced chemical modeling, which is based on a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical network with surface reactions. Both model and observations suggest that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial dependency and source to source variations of CCH column densities qualitatively agree with the observed values, but the predicted column densities are too low by a factor of several. The chemical model fails to reproduce high concentrations of CCH in very cold disk midplane as derived from the observed low excitation condition for both the (1-0) and (2-1) transitions.

Aims. We attempt to understand the presence of gas phase CO below its freezing temperature in circumstellar disks. We study two promising mechanisms to explain this phenomenon: turbulent mixing and photodesorption. Methods. We compute the chemical evolution of circumstellar disks including grain surface reactions with and without turbulent mixing and CO photodesorption. Results. We show that photodesorption significantly enhances the gas phase CO abundance, by extracting CO from the grains when the visual extinction remains below about 5 magnitudes. However the resulting dependence of column density on radial distance is not consistent with observations so far. We propose that this inconsistency could be the result of grain growth. On the other hand, the influence of turbulent mixing is not found to be straightforward. The efficiency of turbulent mixing depends upon a variety of parameters, including the disk structure. For the set of parameters we chose, turbulent mixing is not found to have any significant influence on the CO column density.

The disk-outflow connection is thought to play a key role in extracting excess angular momentum from a forming proto-star. Though jet rotation has been observed in a few objects, no rotation of molecular outflows has been unambiguously reported so far. We report new millimeter-interferometric observations of the edge-on T Tauri star - disk system in the isolated Bok globule CB26. The aim of these observations was to study the disk-outflow relation in this 1Myr old low-mass young stellar object. The IRAM PdBI array was used to observe 12CO(2-1) at 1.3mm in two configurations, resulting in spectral line maps with 1.5 arcsec resolution. We use an empirical parameterized steady-state outflow model combined with 2-D line radiative transfer calculations and chi^2-minimization in parameter space to derive a best-fit model and constrain parameters of the outflow. The data reveal a previously undiscovered collimated bipolar molecular outflow of total length ~2000 AU, escaping perpendicular to the plane of the disk. We find peculiar kinematic signatures that suggest the outflow is rotating with the same orientation as the disk. However, we could not ultimately exclude jet precession or two misaligned flows as possible origin of the observed peculiar velocity field. There is indirect indication that the embedded driving source is a binary system, which, together with the youth of the source, could provide the clue to the observed kinematic features of the outflow. CB26 is so far the most promising source to study the rotation of a molecular outflow. Assuming that the outflow is rotating, we compute and compare masses, mass flux, angular momenta, and angular momentum flux of disk and outflow and derive disk dispersal timescales of 0.5...1 Myr, comparable to the age of the system.

We study the molecular content and chemistry of a circumstellar disk surrounding the Herbig Ae star AB Aur at (sub-)millimeter wavelengths. Our aim is to reconstruct the chemical history and composition of the AB Aur disk and to compare it with disks around low-mass, cooler T Tauri stars. We observe the AB Aur disk with the IRAM Plateau de Bure Interferometer in the C- and D- configurations in rotational lines of CS, HCN, C2H, CH3OH, HCO+, and CO isotopes. Using an iterative minimization technique, observed columns densities and abundances are derived. These values are further compared with results of an advanced chemical model that is based on a steady-state flared disk structure with a vertical temperature gradient, and gas-grain chemical network with surface reactions. We firmly detect HCO+ in the 1--0 transition, tentatively detect HCN, and do not detect CS, C2H, and CH3OH. The observed HCO+ and 13CO column densities as well as the upper limits to the column densities of HCN, CS, C2H, and CH3OH are in good agreement with modeling results and those from previous studies. The AB Aur disk possesses more CO, but is less abundant in other molecular species compared to the DM Tau disk. This is primarily caused by intense UV irradiation from the central Herbig A0 star, which results in a hotter disk where CO freeze out does not occur and thus surface formation of complex CO-bearing molecules might be inhibited.

The overall properties of disks surrounding intermediate PMS stars (HAe) are not yet well constrained by current observations. The disk inclination, which significantly affect SED modeling, is often unknown. We attempted to resolve the disks around CQ Tau and MWC 758, to provide accurate constraints on the disk parameters, in particular the temperature and surface density distribution. We report arcsecond resolution observations of dust and CO line emissions with the IRAM array. The disk properties are derived using a standard disk model. We use the Meudon PDR code to study the chemistry. The two disks share some common properties. The mean CO abundance is low despite disk temperatures above the CO condensation temperature. Furthermore, the CO surface density and dust opacity have different radial dependence. The CQ Tau disk appears warmer, and perhaps less dense than that of MWC 758. Modeling the chemistry, we find that photodissociation of CO is a viable mechanism to explain the low abundance. The photospheric flux is not sufficient for this: a strong UV excess is required. In CQ Tau, the high temperature is consistent with expectation for a PDR. The PDR model has difficulty explaining the mild temperatures obtained in MWC 758, for which a low gas-to-dust ratio is preferred. A yet unexplored alternative could be that, despite currently high gas temperatures, CO remains trapped in grains, as the models suggest that large grains can be cold enough to prevent thermal desorption of CO. The low inclination of the CQ Tau disk, ~30^\circ, challenges previous interpretations given for the UX Ori - like luminosity variations of this star. We conclude that CO cannot be used as a simple tracer of gas-to-dust ratio, the CO abundance being affected by photodissociation, and grain growth.

We present high spatial resolution maps, obtained with the Plateau de Bure Interferometer, of the blue lobe of the L1157 outflow. We observed four lines at 3 mm, namely CH3OH (2_K-1_K), HC3N (11-10), HCN (1-0) and OCS (7-6). Moreover, the bright B1 clump has also been observed at better spatial resolution in CS (2-1), CH3OH (2_1-1_1)A-, and 34SO (3_2-2_1). These high spatial resolution observations show a very rich structure in all the tracers, revealing a clumpy structure of the gas superimposed to an extended emission. In fact, the three clumps detected by previous IRAM-30m single dish observations have been resolved into several sub-clumps and new clumps have been detected in the outflow. The clumps are associated with the two cavities created by two shock episodes driven by the precessing jet. In particular, the clumps nearest the protostar are located at the walls of the younger cavity with a clear arch-shape form while the farthest clumps have slightly different observational characteristics indicating that they are associated to the older shock episode. The emission of the observed species peaks in different part of the lobe: the east clumps are brighter in HC3N (11-10), HCN (1-0) and CS (2-1) while the west clumps are brighter in CH3OH(2_K-1_K), OCS (7-6) and 34SO (3_2-2_1). This peak displacement in the line emission suggests a variation of the physical conditions and/or the chemical composition along the lobe of the outflow at small scale, likely related to the shock activity and the precession of the outflow. In particular, we observe the decoupling of the silicon monoxide and methanol emission, common shock tracers, in the B1 clump located at the apex of the bow shock produced by the second shock episode.

We analyze the line radiative transfer in protoplanetary disks using several approximate methods and a well-tested Accelerated Monte Carlo code. A low-mass flaring disk model with uniform as well as stratified molecular abundances is adopted. Radiative transfer in low and high rotational lines of CO, C18O, HCO+, DCO+, HCN, CS, and H2CO is simulated. The corresponding excitation temperatures, synthetic spectra, and channel maps are derived and compared to the results of the Monte Carlo calculations. A simple scheme that describes the conditions of the line excitation for a chosen molecular transition is elaborated. We find that the simple LTE approach can safely be applied for the low molecular transitions only, while it significantly overestimates the intensities of the upper lines. In contrast, the Full Escape Probability (FEP) approximation can safely be used for the upper transitions ($J_{\rm up} \ga 3$) but it is not appropriate for the lowest transitions because of the maser effect. In general, the molecular lines in protoplanetary disks are partly subthermally excited and require more sophisticated approximate line radiative transfer methods. We analyze a number of approximate methods, namely, LVG, VEP (Vertical Escape Probability) and VOR (Vertical One Ray) and discuss their algorithms in detail. In addition, two modifications to the canonical Monte Carlo algorithm that allow a significant speed up of the line radiative transfer modeling in rotating configurations by a factor of 10--50 are described.

Characterizing the Martian atmosphere is an essential objective to understand its meteorology and its climate. The lower atmosphere (< 40 km) and middle atmosphere (40-80 km) of Mars appear dynamically coupled at much higher levels than in the case of the Earth. The vertical extension of the weather phenomena is considerable with for example Hadley's cells reaching the top of the neutral atmosphere (120 km). The circulation in the middle atmosphere modifies the meteorology of the lower atmosphere, affecting the transport and climatic processes Observations of the CO rotational lines at millimeter (mm) wavelengths (Clancy et al 1990) have strongly contributed in the study of the vertical distribution of this compound and the thermal profile in the atmosphere of Mars over 0-70 km. Singledish observations of the CO Doppler lineshift have allowed direct wind measurements in the martian middle atmosphere near 50 km altitude (Lellouch et al 1991), but at a low spatial resolution (12'') enabling only an essentially hemispheric resolution of the martian disk. The use of mm interferometry has allowed us to better spatially resolve the Martian disk, in order to obtain wind maps of the middle atmosphere (Moreno et al 2001).

We report sub-arcsec CO observations of the disks around MWC 480, LkCa 15 and DM Tau, and simultaneous measurements of HCO+ J=1-0. We derive the disk properties by fitting a standard disk model, with all parameters power laws of the distance to the star. Biases are explained and discussed. We find evidence for vertical temperature gradient in the disks of MWC 480 and DM Tau, as in AB Aur, but not in LkCa 15. The disks temperature increase with stellar effective temperature. Most of the CO gas is at temperatures smaller than 17 K, the condensation temperature on grains. The scale height of the CO distribution appears larger (by 50%) than the hydrostatic scale height. The more UV luminous stars have more CO, but there is no simple dependency of CO abundance and isotopologue ratio with stellar type. The 13CO/HCO+ ratio is around 600. The temperature behaviour is consistent with expectations, but published chemical models have difficulty reproducing the observed CO quantities. Vertical mixing and photo-dissociation at the disk outer edge seem important chemical agents. The CO data suggest a more complex surface density distribution than assumed in models.

Aims: To constrain the ionization fraction in protoplanetary disks, we present new high-sensitivity interferometric observations of N$_2$H$^+$ in three disks surrounding DM Tau, LkCa 15, and MWC 480. Methods: We used the IRAM PdBI array to observe the N$_2$H$^+$ J=1-0 line and applied a $\chi^2$-minimization technique to estimate corresponding column densities. These values are compared, together with HCO$^+$ column densities, to results of a steady-state disk model with a vertical temperature gradient coupled to gas-grain chemistry. Results: We report two \ndhp detections for LkCa 15 and DM Tau at the $5 \sigma$ level and an upper limit for MWC 480. The column density derived from the data for LkCa 15 is much lower than previously reported. The [N$_2$H$^+$/HCO$^+$] ratio is on the order of 0.02--0.03. So far, HCO$^+$ remains the most abundant observed molecular ion in disks. Conclusions: All the observed values generally agree with the modelled column densities of disks at an evolutionary stage of a few million years (within the uncertainty limits), but the radial distribution of the molecules is not reproduced well. The low inferred concentration of N$_2$H$^+$ in three disks around low-mass and intermediate-mass young stars implies that this ion is not a sensitive tracer of the overall disk ionization fraction.

We performed sub-arcsecond high-sensitivity nterferometric observations of the thermal dust emission at 1.4 mm and 2.8 mm in the disks surrounding LkCa 15 and MWC 480, with the new 750 m baselines of the IRAM PdBI array. This provides a linear resolution of about 60 AU at the Taurus distance. We report the existence of a cavity of about 50 AU radius in the inner disk of LkCa 15. Whereas LkCa 15 emission is optically thin, the optically thick core of MWC 480 is resolved at 1.4 mm with a radius of about 35 AU, constraining the dust temperature. In MWC 480, the dust emission is coming from a colder layer than the CO emission, most likely the disk mid-plane. These observations provide direct evidence of an inner cavity around LkCa 15. Such a cavity most probably results from the tidal disturbance created by a low mass companion or large planet at about 30 AU from the star. These results suggest that planetary system formation is already at work in LkCa 15. They also indicate that the classical steady-state viscous disk model is a too simplistic description of the inner 50 AU of ''proto-planetary'' disks, and that the disk evolution is coupled to the planet formation process. The MWC 480 results indicate that a proper estimate of the dust temperature and size of the optically thick core are essential to determine the dust emissivity index.

HH30 is a well-known Pre-Main-Sequence star in Taurus. HST observations have revealed a flared, edge-on disk driving a highly-collimated optical jet, making this object a case study for the disk-jet-outflow paradigm. We obtained high angular resolution (about 1") observations of the dust continuum at 2.7 and 1.3 mm, and of the 12CO(2-1), 13CO(2-1) & (1-0), C18O(1-0) emissions around HH30. A standard disk model is used to fit the 13CO(2-1) uv-plane visibilities and derive the disk properties, and the stellar mass. It results that HH30 is a low mass TTauri of spectral type around M1 and age 1 to 4 Myrs, surrounded by a medium size Keplerian disk, of mass around 4e-3 Msun and outer radius 420 AU. The disk rotation vector points toward the North-Eastern jet. Using a distance of 140 pc, we deduce a stellar mass of 0.45 Msun. A highly asymmetric outflow originates from the inner parts of the disk. It presents to first order a conical morphology with a 30 degree half opening angle and a constant (12 km/s) radial velocity field. Outflow rotation was searched for but not found. These observations do not enable to assign the origin of the molecular outflow to entrainment by the optical jet or to a disk wind. In the latter case, the upper limit of the outflow rotation velocity implies an origin in the inner 15 AU of the disk.

We report the detection of the J=2-1 line of DCO+ in the proto-planetary disk of DM Tau and re-analyze the spectrum covering the 465 GHz transition of HDO in this source, recently published by Ceccarelli et al. (2005). A modelling of the DCO+ line profile with the source parameters derived from high resolution HCO+ observations yields a DCO+/HCO+ abundance ratio of about 0.004, an order of magnitude smaller than that derived in the low mass cores. The re-analysis of the 465 GHz spectrum, using the proper continuum flux (0.5 Jy) and source systemic velocity (6.05 km/s), makes it clear that the absorption features attributed to HDO and C6H are almost certainly unrelated to these species. We show that the line-to-continuum ratio of an absorption line in front of a Keplerian disk can hardly exceed the ratio of the turbulent velocity to the projected rotation velocity at the disk edge, unless the line is optically very thick (tau > 10 000). This ratio is typically 0.1-0.3 in proto-planetary disks and is about 0.15 in DM Tau, much smaller than that for the alleged absorption features. We also show that the detection of H2D+ in DM Tau, previously reported by these authors, is only a 2-sigma detection when the proper velocity is adopted. So far, DCO+ is thus the only deuterated molecule clearly detected in proto-planetary disks.

We present sub-arcsecond images of AB Auriga obtained with the IRAM Plateau de Bure interferometer in the isotopologues of CO, and in continuum at 3 and 1.3 mm. Instead of being centrally peaked, the continuum emission is dominated by a bright, asymmetric (spiral-like) feature at about 140 AU from the central star. The large scale molecular structure suggests the AB Aur disk is inclined between 23 and 43 degrees, but the strong asymmetry of the continuum and molecular emission prevents an accurate determination of the inclination of the inner parts. We find significant non-Keplerian motion, with a best fit exponent for the rotation velocity law of 0.41 +/- 0.01, but no evidence for radial motions. The disk has an inner hole about 70 AU in radius. The disk is warm and shows no evidence of depletion of CO. The dust properties suggest the dust is less evolved than in typical T Tauri disks. Both the spiral-like feature and the departure from purely Keplerian motions indicates the AB Aur disk is not in quasi-equilibrium. Disk self-gravity is insufficient to create the perturbation. This behavior may be related either to an early phase of star formation in which the Keplerian regime is not yet fully established and/or to a disturbance of yet unknown origin. An alternate, but unproven, possibility is that of a low mass companion located about 40 AU from AB Aur.

We present high angular resolution interferometric observations of the dust continuum at 2.7 and 1.3 mm, and of the HC3N (12-11) and C18O (2-1) emission around L1157-mm, a Class~0 object that drives a spectacular molecular outflow. The millimeter dust emission is clearly resolved into two components, a flattened compact source of ~450 x 250 AU at 1.3 mm, and mass ~0.1 Msun, plus an extended envelope of ~3000 AU at 1.3 mm, and mass ~1.1 Msun. The millimeter spectral index varies throughout the region, with the lower value found toward the compact protostar, possibly indicating grain growth in the denser regions. A strong interaction between the molecular outflow and the close protostellar environment is taking place and affects the structure of the innermost parts of the envelope. This is shown by the spatial coincidence between the molecular outflow and the dust (1.3 mm continuum) and HC3N emission: both tracers show structures associated to the edges of the outflow lobes. Basically, the global picture sketched for the Class 0 object L1157-mm by Gueth et al. (1997) is supported. We find possible evidence of infall, but we do not detect any velocity gradient indicative of a rotating circumstellar disk.

We report high-resolution millimeter mapping with the IRAM Plateau de Bure interferometer of rest-frame 335 micron continuum and CO(3-2) line emission from the z=2.8 submillimeter galaxy SMMJ02399-0136. The continuum emission comes from a ~3" diameter structure whose elongation is approximately east-west and whose centroid is coincident within the astrometric errors with the brightest X-ray and rest-UV peak (L1). The line data show that this structure is most likely a rapidly rotating disk. Its rotation velocity of >420 km/s implies a total dynamical mass of >3x10^11 solar masses within an intrinsic radius of 8 kpc, most of which is plausibly in the form of stars and gas. SMMJ02399-0136 is thus a very massive system, whose formation at z~3 is not easy to understand in current CDM hierarchical merger cosmogonies.

New, high-sensitivity interferometric CO J=2-1 observations of the HH 7-11 outflow show that despite previous doubts, this system is powered by the Class I source SVS 13. The molecular outflow from SVS 13 is formed by a shell with a large opening angle at the base, which is typical of outflows from Class I sources, but it also contains an extremely-high-velocity jet composed of ``molecular bullets'', which is more typical of Class 0 outflows. This suggests that SVS 13 could be a very young Class I, which still keeps some features of the previous evolutionary stage. We briefly discuss the nature of some sources in the SVS 13 vicinity which are emitters of cm-wave continuum, but have no counterpart at mm wavelengths.

We have used the high sensitivity and resolution of the IRAM interferometer to produce sub-arcsecond 12CO 2-1 images of 9 protoplanetary disks surrounding T Tauri stars in the Taurus-Auriga cloud (7 singles and 2 binaries). The images demonstrate the disks are in Keplerian rotation around their central stars. Using the least square fit method described in Guilloteau and Dutrey (1998), we derive the disk properties, in particular its inclination angle and rotation velocity, hence the dynamical mass. Since the disk mass is usually small, this is a direct measurement of the stellar mass. Typically, we reach an internal precision of 10% in the determinations of stellar mass. The over-all accuracy is limited by the uncertainty in the distance to a specific star. In a distance independent way, we compare the derived masses with theoretical tracks of pre-main-sequence evolution. Combined with the mean distance to the Taurus region (140 pc), for stars with mass close to 1 Msun, our results tend to favor the tracks with cooler photospheres (higher masses for a given spectral type). We find that in UZ Tau E the disk and the spectroscopic binary orbit appear to have different inclinations.

The z=2.39 radio galaxy 53W002 lies in a cluster of Ly-alpha emission line objects and may itself be undergoing a major burst of star formation. CO(3--2) emission, at 102 GHz, was detected from 53W002 by Scoville et al. (1997a), who also reported a possible 30 kpc extension and velocity gradient suggesting a rotating gaseous disk. In this paper we present new interferometric CO(3--2) observations which confirm the previous line detection with improved signal-to-noise ratio, but show no evidence for source extension or velocity gradient. The compact nature of the CO source and the molecular mass found in this object are similar to luminous infrared galaxies and other AGNs previously studied.

We report the detection of sub-mm emission from dust at 850 microns and of the 12CO J=4-3 line in the two distant powerful radio galaxies 4C 60.07 (z=3.79) and 6C 1909+722 (z=3.53). In the case of 4C 60.07 the dust emission is also detected at 1.25 mm. The estimated molecular gas masses are large, of the order of ~(0.5-1)x10^{11} Solar. The large FIR luminosities (L_fir ~ 10^{13} Solar) suggest that we are witnessing two major starburst phenomena, while the observed large velocity widths (FWHM > 500 km/sec) are characteristic of mergers. In the case of 4C 60.07 the CO emission extends over ~30 kpc and spans a velocity range of >1000 km/sec. It consists of two distinct features with FWHM of >= 550 km/sec and ~150 km/sec and line centers separated by >=700 km/sec The least massive of these components is probably very gas-rich with potentially >=60% of its dynamical mass in the form of molecular gas. The extraordinary morphology of the CO emission in this object suggests that it is not just a scaled-up version of a local Ultra Luminous Infrared Galaxy, and it may be a formative stage of the elliptical host of the residing radio-loud AGN. Finally we briefly explore the effects of the wide range of gas excitation conditions expected for starburst environments on the luminosity of high-J CO lines. We conclude that in unlensed objects, CO (J+1-->J), J+1>3 lines can be significantly weak with respect to CO J=1-0 and this can hinder their detection even in the presence of substantial molecular gas masses.

The IRAM interferometer has been used to detect the submm source HDF 850.1 found by Hughes et al. (1998) in the Hubble Deep Field. The flux density measured at 1.3mm is 2.2 mJy, in agreement with the flux density measured at the JCMT. The flux densities and upper limits measured at 3.4, 2.8, 1.3, 0.85, and 0.45 mm show that the emission is from dust. We suggest that the 1.3mm dust source is associated with the optical arc-like feature, 3-593.0, that has a photometric redshift of about 1.7. If HDF 850.1 is at this redshift and unlensed, its spectral energy distribution, combined with that of 3-593.0, matches closely that of the ultraluminous galaxy VII Zw 31. Another possibility is that the dust source may be gravitationally lensed by the elliptical galaxy 3-586.0 at a redshift of 1.

We report the detection with the IRAM interferometer at 3mm of J=5-4 CO line in the radio quiet quasar BRI1335-0415 at redshift of z=4.41. After BR1202-0725 at z=4.69 (Ohta etal 1996, Omont etal 1996a), this is the second detection of CO at z>3. The integrated line intensity is 2.8+-0.3 Jy km s^{-1} with a linewidth of 420+-60 km s^{-1}. The dust continuum emission has also been mapped at 1.35mm. The 1.35mm flux is found to be 5.6+-1.1mJy. The ratio of the CO to 1.35mm continuum flux is slightly larger than for BR1202-0725. Contrary to the case of BR1202-0725, there is only marginal evidence of extension of the 1.35mm continuum and 3mm CO emission. In the absence of gravitational lensing, for which there is no a priori evidence, and within the uncertainties of the CO to M(H_2) conversion factors, the mass of molecular gas M(H_2) could be 10^{11} solar masses.

We have detected CO(3-2) line emission from the gravitationally lensed quasar MG 0414+0534 at redshift 2.64, using the IRAM Plateau de Bure Interferometer. The line is broad, with Delta v_FWHM = 580 km/s. The velocity-integrated CO flux is comparable to, but somewhat smaller than, that of IRAS F10214+4724 and the Cloverleaf quasar (H1413+117), both of which are at similar redshifts. The lensed components A1+A2 and B were resolved, and separate spectra are presented for each. We also observed the unlensed radio quiet quasar PG 1634+706 at z=1.33, finding no significant CO emission.

We present a new investigation of the Cloverleaf (z=2.558) based on the combination of archival HST/WFPC2 data, recent IRAM CO(7-6) maps and wide field CFHT/FOCAM images. The deepest WFPC2 observation shows a significant overdensity of I ~ 23--25 galaxies around the Cloverleaf that we interpret as the presence of a cluster of galaxies along the line of sight. The typical magnitude, red color (R-I ~ 0.9) and small angular size of these galaxies suggest that the cluster is very distant and could be associated with the absorption systems observed in the spectra of the quasar spots. The Cloverleaf is probably the result of the lensing effects of a system which includes a single-galaxy and one of the most distant cluster of galaxies. We have modelled the lens using altogether the WFPC2 data and the IRAM/CO(7-6) map, considereing two cases: one in which the mass model is a galaxy and a dark halo at z=1.7, and a second one in which the mass model is the combination of a cluster and an individual galaxy located amid the Cloverleaf, both at z=1.7. The CO source is found to form a ring-like structure orbiting the central engine at ~ 100km/s at a radial distance of ~ 100pc, leading to a black-hole mass of ~ 10^9 M_solar. The cluster component increases significantly the convergence of the lens and this pulls down the requirement on the mass of the lensing galaxy by a factor 2. This may explain the mystery of the lensing galaxy not detected yet.

We report the detection of emission from dust and carbon monoxide in the radio-quiet quasar BR1202--0725 at redshift z=4.7. Maps of these emissions reveal two objects separated by a few arc seconds, which could indicate either the presence of a companion to the quasar or gravitational lensing of the quasar itself. Regardless of the precise interpretation of the maps, the detection of carbon monoxide confirms the presence of a large mass of molecular gas in one of the most distant galaxies known, and shows that conditions conducive to huge bursts of star formation existed in the early Universe.