Context. High-resolution imaging reveals a large morphological variety of protoplanetary disks. To date, no constraints on their global evolution have been found from this census. An evolutionary ...classification of disks was proposed based on their IR spectral energy distribution, with the Group I sources showing a prominent cold component ascribed to an earlier stage of evolution than Group II. Aims. Disk evolution can be constrained from the comparison of disks with different properties. A first attempt at disk taxonomy is now possible thanks to the increasing number of high-resolution images of Herbig Ae/Be stars becoming available. Methods. Near-IR images of six Group II disks in scattered light were obtained with VLT/NACO in polarimetric differential imaging, which is the most efficient technique for imaging the light scattered by the disk material close to the stars. We compare the stellar/disk properties of this sample with those of well-studied Group I sources available from the literature. Results. Three Group II disks are detected. The brightness distribution in the disk of HD 163296 indicates the presence of a persistent ring-like structure with a possible connection with the CO snowline. A rather compact (<100 AU) disk is detected around HD 142666 and AK Sco. A taxonomic analysis of 17 Herbig Ae/Be sources reveals that the difference between Group I and Group II is due to the presence or absence of a large disk cavity (≳5 AU). There is no evidence supporting the evolution from Group I to Group II. Conclusions. Group II disks are not evolved versions of the Group I disks. Within the Group II disks, very different geometries exist (both self-shadowed and compact). HD 163296 could be the primordial version of a typical Group I disk. Other Group II disks, like AK Sco and HD 142666, could be smaller counterparts of Group I unable to open cavities as large as those of Group I.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims.
We present MCFOST-art, a new non-local thermodynamic equilibrium radiative transfer solver for multilevel atomic systems. The code is embedded in the 3D radiative transfer code MCFOST and is ...compatible with most of the MCFOST modules. The code is versatile and designed to model the close environment of stars in 3D.
Methods.
The code solves for the statistical equilibrium and radiative transfer equations using the Multilevel Accelerated Lambda Iteration method. We tested MCFOST-art on spherically symmetric models of stellar photospheres as well as on a standard model of the solar atmosphere. We computed atomic level populations and outgoing fluxes and compared these values with the results of the TURBOspectrum and RH codes. Calculations including expansion and rotation of the atmosphere were also performed. We tested both the pure local thermodynamic equilibrium and the out-of-equilibrium problems.
Results.
In all cases, the results from all codes agree within a few percent at all wavelengths and reach the sub-percent level between RH and MCFOST-art. We still note a few marginal discrepancies between MCFOST-art and TURBOspectrum as a result of different treatments of background opacities at some critical wavelength ranges.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims: We aim to define a small and large chemical network which can be used for the quantitative simultaneous analysis of molecular emission from the near-IR to the submm. We also aim to revise ...reactions of excited molecular hydrogen, which are not included in UMIST, to provide a homogeneous database for future applications. Methods: We have used the thermo-chemical disk modeling code ProDiMo and a standard T Tauri disk model to evaluate the impact of various chemical networks, reaction rate databases and sets of adsorption energies on a large sample of chemical species and emerging line fluxes from the near-IR to the submm wavelength range. Results: We find large differences in the masses and radial distribution of ice reservoirs when considering freeze-out on bare or polar ice coated grains. Most strongly the ammonia ice mass and the location of the snow line (water) change. As a consequence molecules associated to the ice lines such as N2H+ change their emitting region; none of the line fluxes in the sample considered here changes by more than 25% except CO isotopologues, CN and N2H+ lines. The three-body reaction N+H2+M plays a key role in the formation of water in the outer disk. Besides that, differences between the UMIST 2006 and 2012 database change line fluxes in the sample considered here by less than a factor of two (a subset of low excitation CO and fine structure lines stays even within 25%); exceptions are OH, CN, HCN, HCO+ and N2H+ lines. However, different networks such as OSU and KIDA 2011 lead to pronounced differences in the chemistry inside 100 au and thus affect emission lines from high excitation CO, OH and CN lines. H2 is easily excited at the disk surface and state-to-state reactions enhance the abundance of CH+ and to a lesser extent HCO+. For sub-mm lines of HCN, N2H+ and HCO+, a more complex larger network is recommended. Conclusions: More work is required to consolidate data on key reactions leading to the formation of water, molecular ions such as HCO+ and N2H+ as well as the nitrogen chemistry. This affects many of the key lines used in the interpretation of disk observations. Differential analysis of various disk models using the same chemical input data will be more robust than the interpretation of absolute fluxes.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims.
We aim to study vertical settling and radial drift of dust in protoplanetary disks from a different perspective: an edge-on view. An estimation of the amplitude of settling and drift is highly ...relevant to understanding planet formation.
Methods.
We analyze a sample of 12 HST-selected edge-on protoplanetary disks (i.e., seen with high inclinations) for which the vertical extent of the emission layers can be constrained directly. We present ALMA high angular resolution continuum images (~0.1′′) of these disks at two wavelengths, 0.89 and 2.06 mm (respectively ALMA bands 7 and 4), supplemented with archival band 6 data (1.33 mm) where available.
Results.
Several sources show constant brightness profiles along their major axis with a steep drop at the outer edges. Two disks have central holes with additional compact continuum emission at the location of the central star. For most sources, the millimeter continuum emission is more compact than the scattered light, both in the vertical and radial directions. Six sources are resolved along their minor axis in at least one millimetric band, providing direct information on the vertical distribution of the millimeter grains. For the second largest disk of the sample, Tau 042021, the significant difference in vertical extent between band 7 and band 4 suggests efficient size-selective vertical settling of large grains. Furthermore, the only Class I object in our sample shows evidence of flaring in the millimeter. Along the major axis, all disks are well resolved. Four of them are larger in band 7 than in band 4 in the radial direction, and three have a similar radial extent in all bands. These three disks are also the ones with the sharpest apparent edges (between 80% and 20% of the peak flux, Δ
r
∕
r
~ 0.3), and two of them are binaries. For all disks, we also derive the millimeter brightness temperature and spectral index maps. We find that all edge-on disks in our sample are likely optically thick and that the dust emission reveals low brightness temperatures in most cases (brightness temperatures ≤10 K). The integrated spectral indices are similar to those of disks at lower inclination.
Conclusions.
The comparison of a generic radiative transfer disk model with our data shows that at least three disks are consistent with a small millimeter dust scale height, of a few au (measured at
r
= 100 au). This is in contrast with the more classical value of
h
g
~ 10 au derived from scattered light images and from gas line measurements. These results confirm, by direct observations, that large (millimeter) grains are subject to significant vertical settling in protoplanetary disks.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Aims.We present a new continuum 3D radiative transfer code, MCFOST, based on a Monte-Carlo method. MCFOST can be used to calculate (i) monochromatic images in scattered light and/or thermal emission; ...(ii) polarisation maps; (iii) interferometric visibilities; (iv) spectral energy distributions; and (v) dust temperature distributions of protoplanetary disks. Methods.Several improvements to the standard Monte Carlo method are implemented in MCFOST to increase efficiency and reduce convergence time, including wavelength distribution adjustments, mean intensity calculations, and an adaptive sampling of the radiation field. The reliability and efficiency of the code are tested against a previously-defined benchmark, using a 2D disk configuration. No significant difference (no more than 10% and usually much less) is found between the temperatures and SEDs calculated by MCFOST and by other codes included in the benchmark. Results. A study of the lowest disk mass detectable by Spitzer, around young stars, is presented and the colours of “representative” parametric disks compared to recent IRAC and MIPS Spitzer colours of solar-like young stars located in nearby star-forming regions.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. The observation of planets in their formation stage is a crucial but very challenging step in understanding when, how, and where planets form. PDS 70 is a young pre-main sequence star ...surrounded by a transition disk, in the gap of which a planetary-mass companion has recently been discovered. This discovery represents the first robust direct detection of such a young planet, possibly still at the stage of formation. Aims. We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods. We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 yr, which allowed us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96–3.8 μm). We use different atmospheric models covering a large parameter space in temperature, log g, chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results. PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range 1000–1600 K and log g no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 RJ with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions. This study provides a comprehensive data set on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical of young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planetflux.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of ...continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of 21 ± 2 K. In the outer disc (>300 au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to ≈15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. Magnetospheric accretion has been thoroughly studied in young stellar systems with full non-evolved accretion disks, but it is poorly documented for transition disk objects with large inner ...cavities. Aims. We aim at characterizing the star-disk interaction and the accretion process onto the central star of LkCa 15, a prototypical transition disk system with an inner dust cavity that is 50 au wide. Methods. We obtained quasi-simultaneous photometric and spectropolarimetric observations of the system over several rotational periods. We analyzed the system light curve and associated color variations, as well as changes in spectral continuum and line profile to derive the properties of the accretion flow from the edge of the inner disk to the central star. We also derived magnetic field measurements at the stellar surface. Results. We find that the system exhibits magnetic, photometric, and spectroscopic variability with a period of about 5.70 days. The light curve reveals a periodic dip, which suggests the presence of an inner disk warp that is located at the corotation radius at about 0.06 au from the star. Line profile variations and veiling variability are consistent with a magnetospheric accretion model where the funnel flows reach the star at high latitudes. This leads to the development of an accretion shock close to the magnetic poles. All diagnostics point to a highly inclined inner disk that interacts with the stellar magnetosphere. Conclusions. The spectroscopic and photometric variability on a timescale of days to weeks of LkCa 15 is remarkably similar to that of AA Tau, the prototype of periodic dippers. We therefore suggest that the origin of the variability is a rotating disk warp that is located at the inner edge of a highly inclined disk close to the star. This contrasts with the moderate inclination of the outer transition disk seen on the large scale and thus provides evidence for a significant misalignment between the inner and outer disks of this planet-forming transition disk system.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We report 885 m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging ...from 1.0 to 55.7 mJy. The two nondetections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from ∼0.3 to 20 M⊕. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in Mdust to Mstar, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early-M stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas-to-dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M dwarfs.
Context.
In recent years, protoplanetary disks with spiral structures have been detected in scattered light, millimeter continuum, and CO gas emission. The mechanisms causing these structures are ...still under debate. A popular scenario to drive the spiral arms is the one of a planet perturbing the material in the disk. However, if the disk is massive, gravitational instability isusually the favored explanation. Multiwavelength studies could be helpful to distinguish between the two scenarios. So far, only a handful of disks with spiral arms have been observed in both scattered light and millimeter continuum.
Aims.
We aim to perform an in-depth characterization of the protoplanetary disk morphology around WaOph 6 analyzing data obtained at different wavelengths, as well as to investigate the origin of the spiral features in the disk.
Methods.
We present the first near-infrared polarimetric observations of WaOph 6 obtained with SPHERE at the VLT and compare them to archival millimeter continuum ALMA observations. We traced the spiral features in both data sets and estimated the respective pitch angles. We discuss the different scenarios that can give rise to the spiral arms in WaOph 6. We tested the planetary perturber hypothesis by performing hydrodynamical and radiative transfer simulations to compare them with scattered light and millimeter continuum observations.
Results.
We confirm that the spiral structure is present in our polarized scattered light
H
-band observationsof WaOph 6, making it the youngest disk with spiral arms detected at these wavelengths. From the comparison to the millimeter ALMA-DSHARP observations, we confirm that the disk is flared. We explore the possibility of a massive planetary perturber driving the spiral arms by running hydrodynamical and radiative transfer simulations, and we find that a planet of minimum 10
M
Jup
outside of the observed spiral structure is able to drive spiral arms that resemble the ones in the observations. We derive detection limits from our SPHERE observations and get estimates of the planet’s contrast from different evolutionary models.
Conclusions.
Up to now, no spiral arms had been observed in scattered light in disks around K and/or M stars with ages <1 Myr. Future observations of WaOph 6 could allow us to test theoretical predictions for planet evolutionary models, as well as give us more insightinto the mechanisms driving the spiral arms.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK