Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks.

The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanets will sallow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine using N-body simulations how the properties of planetary systems are determined during the final stages of assembly. While accretion is a chaotic process, trends in the ensemble properties of planetary systems provide a memory of the initial distribution of solid mass around a star prior to accretion. We also use EPOS, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, either via planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes has a profound impact on the observed planetary architectures, and thus on the formation and delivery of volatiles to possible habitable worlds.

Planet yield calculations may be used to inform the target selection strategy and science operations of space observatories. Forthcoming and proposed NASA missions, such as the Wide-Field Infrared Survey Telescope (WFIRST), the Habitable Exoplanet Imaging Mission (HabEx), and the Large UV/Optical/IR Surveyor (LUVOIR), are expected to be equipped with sensitive coronagraphs and/or starshades. We are developing a suite of numerical simulations to quantify the extent to which ground-based radial velocity (RV) surveys could boost the detection efficiency of direct imaging missions. In this paper, we discuss the first step in the process of estimating planet yields: generating synthetic planetary systems consistent with observed occurrence rates from multiple detection methods. In an attempt to self-consistently populate stars with orbiting planets, it is found that naive extrapolation of occurrence rates (mass, semi-major axis) results in an unrealistically large number-density of Neptune-mass planets beyond the ice-line ($a \gtrsim 5$au), causing dynamic interactions that would destabilize orbits. We impose a stability criterion for multi-planet systems based on mutual Hill radii separation. Considering the influence of compact configurations containing Jovian-mass and Neptune-mass planets results in a marked suppression in the number of terrestrial planets that can exist at large radii. This result has a pronounced impact on planet yield calculations particularly in regions accessible to high-contrast imaging and microlensing. The dynamically compact configurations and occurrence rates that we develop may be incorporated as input into joint RV and direct imaging yield calculations to place meaningful limits on the number of detectable planets with future missions.

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos

The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.

The frequency of Earth-size planets in the habitable zone of Sun-like stars, hereafter $\eta_\oplus$, is a key parameter to evaluate the yield of nearby Earth analogues that can be detected and characterized by future missions. Yet, this value is poorly constrained as there are no reliable exoplanet candidates in the habitable zone of Sun-like stars in the Kepler field. Here, we show that extrapolations relying on the population of small ($< 1.8\,R_\oplus$) short-period ($< 25\,$days) planets bias $\eta_\oplus$ to large values. As the radius distribution at short orbital periods is strongly affected by atmospheric loss, we re-evaluate $\eta_\oplus$ using exoplanets at larger separations. We find that $\eta_\oplus$ drops considerably, to values of only $\sim 5-10$%. Observations of young ($< 100$ Myr) clusters can probe short-period sub-Neptunes that still retain most of their envelope mass. As such, they can be used to quantify the contamination of sub-Neptunes to the population of Kepler short-period small planets and aid in more reliable estimates of $\eta_\oplus$.

The detection of transiting exoplanets in time-series photometry requires the removal or modeling of instrumental and stellar noise. While instrumental systematics can be reduced using methods such as pixel level decorrelation, removing stellar trends while preserving transit signals proves challenging. Due to vast archives of light curves from recent transit surveys, there is a strong need for accurate automatic detrending, without human intervention. A large variety of detrending algorithms are in active use, but their comparative performance for transit discovery is unexplored. We benchmark all commonly used detrending methods against hundreds of Kepler, K2, and TESS planets, selected to represent the most difficult cases for systems with small planet-to-star radius ratios. The full parameter range is explored for each method to determine the best choices for planet discovery. We conclude that the ideal method is a time-windowed slider with an iterative robust location estimator based on Tukey's biweight. This method recovers 99% and 94% of the shallowest Kepler and K2 planets, respectively. We include an additional analysis for young stars with extreme variability and conclude they are best treated using a spline-based method with a robust Huber estimator. All stellar detrending methods explored are available for public use in wotan, an open-source Python package on GitHub (see https://github.com/hippke/wotan).

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

Two studies utilizing sparse aperture masking (SAM) interferometry and $H_{\rm \alpha}$ differential imaging have reported multiple jovian companions around the young solar-mass star, LkCa 15 (LkCa 15 bcd): the first claimed direct detection of infant, newly-formed planets ("protoplanets"). We present new near-infrared direct imaging/spectroscopy from the SCExAO system coupled with the CHARIS integral field spectrograph and multi-epoch thermal infrared imaging from Keck/NIRC2 of LkCa 15 at high Strehl ratios. These data provide the first direct imaging look at the same wavelengths and in the same locations where previous studies identified the LkCa 15 protoplanets and thus offer the first decisive test of their existence. The data do not reveal these planets. Instead, we resolve extended emission tracing a dust disk with a brightness and location comparable to that claimed for LkCa 15 bcd. Forward-models attributing this signal to orbiting planets are inconsistent with the combined SCExAO/CHARIS and Keck/NIRC2 data. An inner disk provides a more compelling explanation for the SAM detections and perhaps also the claimed $H_{\alpha}$ detection of LkCa 15 b. We conclude that there is currently no clear, direct evidence for multiple protoplanets orbiting LkCa 15, although the system likely contains at least one unseen jovian companion. To identify jovian companions around LkCa 15 from future observations, the inner disk should be detected and its effect modeled, removed, and shown to be distinguishable from planets. Protoplanet candidates identified from similar systems should likewise be clearly distinguished from disk emission through modeling.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

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

The orbital distribution of giant planets is crucial for understanding how terrestrial planets form and predicting yields of exoplanet surveys. Here, we derive giant planets occurrence rates as a function of orbital period by taking into account the detection efficiency of the Kepler and radial velocity (RV) surveys. The giant planet occurrence rates for Kepler and RV show the same rising trend with increasing distance from the star. We identify a break in the RV giant planet distribution between ~2-3 au -- close to the location of the snow line in the Solar System -- after which the occurrence rate decreases with distance from the star. Extrapolating a broken power-law distribution to larger semi-major axes, we find good agreement with the ~ 1% planet occurrence rates from direct imaging surveys. Assuming a symmetric power law, we also estimate that the occurrence of giant planets between 0.1-100 au is 26.6 +7.5 -5.4% for planets with masses 0.1-20MJ and decreases to 6.2 +1.5 -1.2% for planets more massive than Jupiter. This implies that only a fraction of the structures detected in disks around young stars can be attributed to giant planets. Various planet population synthesis models show good agreement with the observed distribution, and we show how a quantitative comparison between model and data can be used to constrain planet formation and migration mechanisms.

We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,\mu$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$\mu$Jy/beam, and reveal an asymmetric ring between $\sim$20-40 au with largely optically thin dust continuum emission. This ring is well fit by two concentric and overlapping Gaussian rings of different widths and a Vortex. In addition, an unresolved component is detected at a position consistent with the central star, which may trace the central inner disk ($<$2au in radius). We report a lack of compact continuum emission at the positions of both claimed protoplanets. We use this result to constrain the circumplanetary disk (CPD) mass and size of 1.44M$_{\rm Earth}$ and 0.44au in the optically thin and thick regime, respectively, for the case of the previously directly imaged protoplanet candidate at $\sim$55 au (HD100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 M$_J$ by using the relation between planet, circumstellar, and circumplanetary masses derived from numerical simulations. Finally, the CO integrated intensity map shows a possible spiral arm feature that would match the spiral features identified in Near-Infrared scattered light polarized emission, which suggests a real spiral feature in the disk surface that needs to be confirmed with further observations.

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

We refine the gap size measurements of the disk surrounding the Herbig Ae star HD 100546 in the N band. Our new mid-infrared interferometric (MIDI) data have been taken with the UT baselines and span the full range of orientations. The correlated fluxes show a wavy pattern in which the minima separation links to a geometrical structure in the disk. We fit each correlated flux measurement with a spline function, deriving the corresponding spatial scale, while assuming that the pattern arises interferometrically due to the bright emission from the inner disk and the opposing sides of the wall of the outer disk. We then fit an ellipse to the derived separations at their corresponding position angles, thereby using the observations to constrain the disk inclination to i =47 +/- 1 degree and the disk position angle to PA =135.0 +/- 2.5 degree East of North, both of which are consistent with the estimated values in previous studies. We also derive the radius of the ellipse to 15.7 +/- 0.8 au. To confirm that the minima separations translate to a geometrical structure in the disk, we model the disk of HD 100546 using a semi-analytical approach taking into account the temperature and optical depth gradients. Using this model, we simultaneously reproduce the level and the minima of the correlated fluxes and constrain the gap size of the disk for each observation. The values obtained for the projected gap size in different orientations are consistent with the separation found by the geometrical model.

ALMA surveys of nearby star-forming regions have shown that the dust mass in the disk is correlated with the stellar mass, but with a large scatter. This scatter could indicate either different evolutionary paths of disks or different initial conditions within a single cluster. We present ALMA Cycle 3 follow-up observations for 14 Class II disks that were low S/N detections or non-detections in our Cycle 2 survey of the $\sim 2$ Myr-old Chamaeleon I star-forming region. With 5 times better sensitivity, we detect millimeter dust continuum emission from six more sources and increase the detection rate to 94\% (51/54) for Chamaeleon I disks around stars earlier than M3. The stellar-disk mass scaling relation reported in \citet{pascucci2016} is confirmed with these updated measurements. Faint outliers in the $F_{mm}$--$M_*$ plane include three non-detections (CHXR71, CHXR30A, and T54) with dust mass upper limits of 0.2 M$_\oplus$ and three very faint disks (CHXR20, ISO91, and T51) with dust masses $\sim 0.5$ M$_\oplus$. By investigating the SED morphology, accretion property and stellar multiplicity, we suggest for the three millimeter non-detections that tidal interaction by a close companion ($<$100 AU) and internal photoevaporation may play a role in hastening the overall disk evolution. The presence of a disk around only the secondary star in a binary system may explain the observed stellar SEDs and low disk masses for some systems.

The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the latest Kepler data release, DR25. We estimate that at least 42% of sun-like stars have nearly coplanar planetary systems with 7 or more exoplanets. The fraction of stars with at least one planet within 1 au could be as high as 100% depending on assumptions about the distribution of single transiting planets. We estimate an occurrence rate of planets in the habitable zone around sun-like stars of eta_earth=36+-14%. The innermost planets in multi-planet systems are clustered around an orbital period of 10 days (0.1 au), reminiscent of the protoplanetary disk inner edge or could be explained by a planet trap at that location. Only a small fraction of planetary systems have the innermost planet at long orbital periods, with fewer than ~8% and ~3% having no planet interior to the orbit of Mercury and Venus, respectively. These results reinforce the view that the solar system is not a typical planetary system, but an outlier among the distribution of known exoplanetary systems. We predict that at least half of the habitable zone exoplanets are accompanied by (non-transiting) planets at shorter orbital periods, hence knowledge of a close-in exoplanet could be used as a way to optimize the search for Earth-size planets in the Habitable Zone with future direct imaging missions.

Exoplanets around different types of stars provide a window into the diverse environments in which planets form. This chapter describes the observed relations between exoplanet populations and stellar properties and how they connect to planet formation in protoplanetary disks. Giant planets occur more frequently around more metal-rich and more massive stars. These findings support the core accretion theory of planet formation, in which the cores of giant planets form more rapidly in more metal-rich and more massive protoplanetary disks. Smaller planets, those with sizes roughly between Earth and Neptune, exhibit different scaling relations with stellar properties. These planets are found around stars with a wide range of metallicities and occur more frequently around lower mass stars. This indicates that planet formation takes place in a wide range of environments, yet it is not clear why planets form more efficiently around low mass stars. Going forward, exoplanet surveys targeting M dwarfs will characterize the exoplanet population around the lowest mass stars. In combination with ongoing stellar characterization, this will help us understand the formation of planets in a large range of environments.

Transitional disks around the Herbig Ae/Be stars are fascinating targets in the contexts of disk evolution and also planet formation. Oph IRS 48 is one of such Herbig Ae stars, which shows an inner dust cavity and azimuthally lopsided large dust distribution. We present new images of Oph IRS 48 at eight mid-infrared (MIR) wavelengths from 8.59 to 24.6\,$\mu {\rm m}$ taken with the COMICS mounted on the 8.2\,m Subaru Telescope. The N-band (7 to 13\,$\mu {\rm m}$) images show that the flux distribution is centrally peaked with a slight spatial extent, while the Q-band (17 to 25\,$\mu {\rm m}$) images show asymmetric double peaks (east and west). Using 18.8 and 24.6\,$\mu$m images, we derived the dust temperature at both east and west peaks to be 135$\pm$22 K. Thus, the asymmetry may not be attributed to a difference in the temperature. % thus other reason is necessary to explain the asymmetry of east and west peaks in Q-band. Comparing our results with previous modeling works, we conclude that the inner disk is aligned to the outer disk. A shadow casted by the optically thick inner disk has a great influence on the morphology of MIR thermal emission from the outer disk.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

We follow the microlensing approach and quantify the occurrence of Kepler exoplanets as a function of planet-to-star mass ratio, q, rather than planet radius or mass. For planets with radii ~1-6 R_earth and periods <100 days, we find that, except for a normalization factor, the occurrence rate vs q can be described by the same broken power law with a break at ~3 x10^{-5} independent of host type for hosts below 1 M_sun. These findings indicate that the planet-to-star mass ratio is a more fundamental quantity in planet formation than planet mass. We then compare our results to those from microlensing for which the overwhelming majority satisfies the M_host<1 M_sun criterion. The break in q for the microlensing planet population, which mostly probes the region outside the snowline, is 3-10 times higher than that inferred from Kepler. Thus, the most common planet inside the snowline is ~3-10 times less massive than the one outside. With rocky planets interior to gaseous planets, the Solar System broadly follows the combined mass-ratio function inferred from Kepler and microlensing. However, the exoplanet population has a less extreme radial distribution of planetary masses than the Solar System. Establishing whether the mass-ratio function beyond the snowline is also host type independent will be crucial to build a comprehensive theory of planet formation.

Future NASA concept missions that are currently under study, like Habitable Exoplanet Imaging Mission (HabEx) & Large Ultra-Violet Optical Infra Red (LUVOIR) Surveyor, would discover a large diversity of exoplanets. We propose here a classification scheme that distinguishes exoplanets into different categories based on their size and incident stellar flux, for the purpose of providing the expected number of exoplanets observed (yield) with direct imaging missions. The boundaries of this classification can be computed using the known chemical behavior of gases and condensates at different pressures and temperatures in a planetary atmosphere. In this study, we initially focus on condensation curves for sphalerite ZnS, H2O, CO2 and CH4. The order in which these species condense in a planetary atmosphere define the boundaries between different classes of planets. Broadly, the planets are divided into rocky (0.5 - 1.0RE), super-Earths (1.0- 1.75RE), sub-Neptunes (1.75-3.5RE), sub-Jovians (3.5 - 6.0RE) and Jovians (6-14.3RE) based on their planet sizes, and 'hot', 'warm' and 'cold' based on the incident stellar flux. We then calculate planet occurrence rates within these boundaries for different kinds of exoplanets, \eta_{planet}, using the community co-ordinated results of NASA's Exoplanet Program Analysis Group's Science Analysis Group-13 (SAG-13). These occurrence rate estimates are in turn used to estimate the expected exoplanet yields for direct imaging missions of different telescope diameter.

We present high resolution ($0.11"\times0.06"$) 3mm ALMA observations of the highly inclined transition disk around the star T Cha. Our continuum image reveals multiple dust structures: an inner disk, a spatially resolved dust gap, and an outer ring. When fitting sky-brightness models to the real component of the 3mm visibilities, we infer that the inner emission is compact ($\le1$au in radius), the gap width is between 18-28 au, and the emission from the outer ring peaks at $\sim36$ au. We compare our ALMA image with previously published 1.6$\mu$m VLT/SPHERE imagery. This comparison reveals that the location of the outer ring is wavelength dependent. More specifically, the peak emission of the 3mm ring is at a larger radial distance than that of the 1.6$\mu$m ring, suggesting that millimeter-sized grains in the outer disk are located further away from the central star than micron-sized grains. We discuss different scenarios to explain our findings, including dead zones, star-driven photoevaporation, and planet-disk interactions. We find that the most likely origin of the dust gap is from an embedded planet, and estimate --- for a single planet scenario --- that T Cha's gap is carved by a $1.2M_{jup}$ planet.

We investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. Dust mass and mass accretion rate in Chamaeleon I are correlated with a slope close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit Mdust-Mstar and Macc-Mstar relations is uncorrelated. Disks with a constant alpha viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass, but over-predict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. 1) The observed scatter in Mdust and Macc is not primoridal, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large time scales affects the mass accretion rates. 2) The observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I due to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, either through direct observations of the gas or spatially resolved multi-wavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

The dependence of the mass accretion rate on the stellar properties is a key constraint for star formation and disk evolution studies. Here we present a study of a sample of stars in the Chamaeleon I star forming region carried out using the VLT/X-Shooter spectrograph. The sample is nearly complete down to M~0.1Msun for the young stars still harboring a disk in this region. We derive the stellar and accretion parameters using a self-consistent method to fit the broad-band flux-calibrated medium resolution spectrum. The correlation between the accretion luminosity to the stellar luminosity, and of the mass accretion rate to the stellar mass in the logarithmic plane yields slopes of 1.9 and 2.3, respectively. These slopes and the accretion rates are consistent with previous results in various star forming regions and with different theoretical frameworks. However, we find that a broken power-law fit, with a steeper slope for stellar luminosity smaller than ~0.45 Lsun and for stellar masses smaller than ~ 0.3 Msun, is slightly preferred according to different statistical tests, but the single power-law model is not excluded. The steeper relation for lower mass stars can be interpreted as a faster evolution in the past for accretion in disks around these objects, or as different accretion regimes in different stellar mass ranges. Finally, we find two regions on the mass accretion versus stellar mass plane empty of objects. One at high mass accretion rates and low stellar masses, which is related to the steeper dependence of the two parameters we derived. The second one is just above the observational limits imposed by chromospheric emission. This empty region is located at M~0.3-0.4Msun, typical masses where photoevaporation is known to be effective, and at mass accretion rates ~10^-10 Msun/yr, a value compatible with the one expected for photoevaporation to rapidly dissipate the inner disk.

The mass of a protoplanetary disk limits the formation and future growth of any planet. Masses of protoplanetary disks are usually calculated from measurements of the dust continuum emission by assuming an interstellar gas-to-dust ratio. To investigate the utility of CO as an alternate probe of disk mass, we use ALMA to survey $^{13}$CO and C$^{18}$O J = $3-2$ line emission from a sample of 93 protoplanetary disks around stars and brown dwarfs with masses from 0.03 -- 2 M$_{\odot}$ in the nearby Chamaeleon I star-forming region. We detect $^{13}$CO emission from 17 sources and C$^{18}$O from only one source. Gas masses for disks are then estimated by comparing the CO line luminosities to results from published disk models that include CO freeze-out and isotope-selective photodissociation. Under the assumption of a typical ISM CO-to-H$_2$ ratios of $10^{-4}$, the resulting gas masses are implausibly low, with an average gas mass of $\sim$ 0.05 M$_{Jup}$ as inferred from the average flux of stacked $^{13}$CO lines. The low gas masses and gas-to-dust ratios for Cha I disks are both consistent with similar results from disks in the Lupus star-forming region. The faint CO line emission may instead be explained if disks have much higher gas masses, but freeze-out of CO or complex C-bearing molecules is underestimated in disk models. The conversion of CO flux to CO gas mass also suffers from uncertainties in disk structures, which could affect gas temperatures. CO emission lines will only be a good tracer of the disk mass when models for C and CO depletion are confirmed to be accurate.

The properties of disks around brown dwarfs and very-low mass stars (hereafter VLMOs) provide important boundary conditions on the process of planet formation and inform us about the numbers and masses of planets than can form in this regime. We use the Herschel Space Observatory PACS spectrometer to measure the continuum and [OI] 63um line emission towards 11 VLMOs with known disks in the Taurus and Chamaeleon I star-forming regions. We fit radiative transfer models to the spectral energy distributions of these sources. Additionally, we carry out a grid of radiative transfer models run in a regime that connects the luminosity of our sources with brighter T~Tauri stars. We find VLMO disks with sizes [1.3--78] au, smaller than typical T~Tauri disks, fit well the spectral energy distributions assuming disk geometry and dust properties are stellar-mass independent. Reducing the disk size increases the disk temperature and we show that VLMOs do not follow previously derived disk temperature-stellar luminosity relationships if the disk outer radius scales with stellar mass. Only 2 out of 11 sources are detected in [OI] despite a better sensitivity than was achieved for T Tauri stars, suggesting that VLMO disks are underluminous. Using thermochemical models we show that smaller disks can lead to the unexpected [OI] non-detections in our sample. The disk outer radius is an important factor in determining the gas and dust observables. Hence, spatially resolved observations with ALMA -- to establish if and how disk radii scale with stellar mass -- should be pursued further.

High-contrast scattered light observations have revealed the surface morphology of several dozens of protoplanetary disks at optical and near-infrared wavelengths. Inclined disks offer the opportunity to measure part of the phase function of the dust grains that reside in the disk surface which is essential for our understanding of protoplanetary dust properties and the early stages of planet formation. We aim to construct a method which takes into account how the flaring shape of the scattering surface of an (optically thick) protoplanetary disk projects onto the image plane of the observer. This allows us to map physical quantities (scattering radius and scattering angle) onto scattered light images and retrieve stellar irradiation corrected (r^2-scaled) images and dust phase functions. We apply the method on archival polarized intensity images of the protoplanetary disk around HD 100546 that were obtained with VLT/SPHERE in R'-band and VLT/NACO in H- and Ks-band. The brightest side of the r^2-scaled R'-band polarized intensity image of HD 100546 changes from the far to the near side of the disk when a flaring instead of a geometrically flat disk surface is used for the r^2-scaling. The decrease in polarized surface brightness in the scattering angle range of ~40-70 deg is likely a result of the dust phase function and degree of polarization which peak in different scattering angle regimes. The derived phase functions show part of a forward scattering peak which indicates that large, aggregate dust grains dominate the scattering opacity in the disk surface. Projection effects of a protoplanetary disk surface need to be taken into account to correctly interpret scattered light images. Applying the correct scaling for the correction of stellar irradiation is crucial for the interpretation of the images and the derivation of the dust properties in the disk surface layer.

The host star metallicity provide a measure of the conditions in protoplanetary disks at the time of planet formation. Using a sample of over 20,000 Kepler stars with spectroscopic metallicities from the LAMOST survey, we explore how the exoplanet population depends on host star metallicity as a function of orbital period and planet size. We find that exoplanets with orbital periods less than 10 days are preferentially found around metal-rich stars ([Fe/H]~ 0.15 +- 0.05 dex). The occurrence rates of these hot exoplanets increases to ~30% for super-solar metallicity stars from ~10% for stars with a sub-solar metallicity. Cooler exoplanets, that resides at longer orbital periods and constitute the bulk of the exoplanet population with an occurrence rate of >~ 90%, have host-star metallicities consistent with solar. At short orbital periods, P<10 days, the difference in host star metallicity is largest for hot rocky planets (<1.7 R_Earth), where the metallicity difference is [Fe/H] =~ 0.25 +- 0.07 dex. The excess of hot rocky planets around metal-rich stars implies they either share a formation mechanism with hot Jupiters, or trace a planet trap at the protoplanetary disk inner edge which is metallicity-dependent. We do not find statistically significant evidence for a previously identified trend that small planets toward the habitable zone are preferentially found around low-metallicity stars. Refuting or confirming this trend requires a larger sample of spectroscopic metallicities.

LkCa 15 hosts a pre-transitional disk as well as at least one accreting protoplanet orbiting in its gap. Previous disk observations have focused mainly on the outer disk, which is cleared inward of ~50 au. The planet candidates, on the other hand, reside at orbital radii around 15 au, where disk observations have been unreliable until recently. Here we present new J-band imaging polarimetry of LkCa 15 with SPHERE IRDIS, yielding the most accurate and detailed scattered-light images of the disk to date down to the planet-hosting inner regions. We find what appear to be persistent asymmetric structures in the scattering material at the location of the planet candidates, which could be responsible at least for parts of the signals measured with sparse-aperture masking. These images further allow us to trace the gap edge in scattered light at all position angles and search the inner and outer disks for morphological substructure. The outer disk appears smooth with slight azimuthal variations in polarized surface brightness, which may be due to shadowing from the inner disk or a two-peaked polarized phase function. We find that the near-side gap edge revealed by polarimetry matches the sharp crescent seen in previous ADI imaging very well. Finally, the ratio of polarized disk to stellar flux is more than six times larger in J-band than in the RI bands.

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets. We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the Pynpoint package. The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate 'c' might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of 'b'. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to 'b'.

Trends in the planet population with host star mass provide an avenue to constrain planet formation theories. We derive the planet radius distribution function for Kepler stars of different spectral types, sampling a range in host star masses. We find that M dwarf stars have 3.5 times more small planets (1.0-2.8 R_Earth) than main-sequence FGK stars, but two times fewer Neptune-sized and larger planets (>2.8 R_Earth). We find no systematic trend in the planet size distribution between spectral types F, G, and K to explain the increasing occurrence rates. Taking into account the mass-radius relationship and heavy-element mass of observed exoplanets, and assuming those are independent of spectral type, we derive the inventory of the heavy-element mass locked up in exoplanets at short orbits. The overall higher planet occurrence rates around M stars are not consistent with the redistribution of the same mass into more, smaller planets. At the orbital periods and planet radii where Kepler observations are complete for all spectral types, the average heavy-element mass locked up in exoplanets increases roughly inversely with stellar mass from 4 M_Earth in F stars to 5 M_Earth in G and K stars to 7 M_Earth in M stars. This trend stands in stark contrast with observed protoplanetary disk masses that decrease towards lower mass stars, and provides a challenge for current planet formation models. Neither models of in situ formation nor migration of fully-formed planets are consistent with these results. Instead, these results are indicative of large-scale inward migration of planetary building blocks --- either through type-I migration or radial drift of dust grains --- that is more efficient for lower mass stars, but does not result in significantly larger or smaller planets.

We present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.

We aim to constrain the structure of the circumstellar material around the post-AGB binary and RV Tauri pulsator AC Her. We want to constrain the spatial distribution of the amorphous as well as of the crystalline dust. We present very high-quality mid-IR interferometric data that were obtained with MIDI/VLTI. We analyse the MIDI data and the full SED, using the MCMax radiative transfer code, to find a good structure model of AC Her's circumbinary disk. We include a grain size distribution and midplane settling of dust self-consistently. The spatial distribution of crystalline forsterite in the disk is investigated with the mid-IR features, the 69~$\mu$m band and the 11.3~$\mu$m signatures in the interferometric data. All the data are well fitted. The inclination and position angle of the disk are well determined at i=50+-8 and PA=305+-10. We firmly establish that the inner disk radius is about an order of magnitude larger than the dust sublimation radius. Significant grain growth has occurred, with mm-sized grains being settled to the midplane of the disk. A large dust mass is needed to fit the sub-mm fluxes. By assuming {\alpha}=0.01, a good fit is obtained with a small grain size power law index of 3.25, combined with a small gas/dust ratio <10. The resulting gas mass is compatible with recent estimates employing direct gas diagnostics. The spatial distribution of the forsterite is different from the amorphous dust, as more warm forsterite is needed in the surface layers of the inner disk. The disk in AC Her is very evolved, with its small gas/dust ratio and large inner hole. Mid-IR interferometry offers unique constraints, complementary to mid-IR features, for studying the mineralogy in disks. A better uv coverage is needed to constrain in detail the distribution of the crystalline forsterite in AC Her, but we find strong similarities with the protoplanetary disk HD100546.

The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.

We imaged circumstellar disks around 22 Herbig Ae/Be stars at 25 \mu m using Subaru/COMICS and Gemini/T-ReCS. Our sample consists of equal numbers of objects belonging to the two categories defined by Meeus et al. (2001); 11 group I (flaring disk) and II (at disk) sources. We find that group I sources tend to show more extended emission than group II sources. Previous studies have shown that the continuous disk is hard to be resolved with 8 meter class telescopes in Q-band due to the strong emission from the unresolved innermost region of the disk. It indicates that the resolved Q-band sources require a hole or gap in the disk material distribution to suppress the contribution from the innermost region of the disk. As many group I sources are resolved at 25 \mu m, we suggest that many, not all, group I Herbig Ae/Be disks have a hole or gap and are (pre-)transitional disks. On the other hand, the unresolved nature of many group II sources at 25 \mu m supports that group II disks have continuous at disk geometry. It has been inferred that group I disks may evolve into group II through settling of dust grains to the mid-plane of the proto-planetary disk. However, considering growing evidence for the presence of a hole or gaps in the disk of group I sources, such an evolutionary scenario is unlikely. The difference between groups I and II may reflect different evolutionary pathways of protoplanetary disks.

Most models of volatile delivery to accreting terrestrial planets assume that the carriers for water are similar in water content to the carbonaceous chondrites in our Solar System. Here we suggest that the water content of primitive bodies in many planetary systems may actually be much higher, as carbonaceous chondrites have lost some of their original water due to heating from short-lived radioisotopes that drove parent body alteration. Using N-body simulations, we explore how planetary accretion would be different if bodies beyond the water line contained a water mass fraction consistent with chemical equilibrium calculations, and more similar to comets, as opposed to the more traditional water-depleted values. We apply this model to consider planet formation around stars of different masses and identify trends in the properties of Habitable Zone planets and planetary system architecture which could be tested by ongoing exoplanet census data collection. Comparison of such data with the model predicted trends will serve to evaluate how well the N-body simulations and the initial conditions used in studies of planetary accretion can be used to understand this stage of planet formation.

In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 114 disks in different evolutionary stages (full, transitional and debris disks) collected from the {\it Spitzer} archive. We focus on the two brighter {H~{\sc i}} lines observed in the {\it Spitzer} spectra, the {H~{\sc i}}(7-6) at 12.37$\mu$m and the {H~{\sc i}}(9-7) at 11.32$\mu$m. We detect the {H~{\sc i}}(7-6) line in 46 objects, and the {H~{\sc i}}(9-7) in 11. We compare these lines with the other most common gas line detected in {\it Spitzer} spectra, the {[Ne~{\sc iii}]} at 12.81$\mu$m. We argue that it is unlikely that the {H~{\sc i}} emission originates from the photoevaporating upper surface layers of the disk, as has been found for the {[Ne~{\sc iii}]} lines toward low-accreting stars. Using the {H~{\sc i}}(9-7)/{H~{\sc i}}(7-6) line ratios we find these gas lines are likely probing gas with hydrogen column densities of 10$^{10}$-10$^{11}$~cm$^{-3}$. The subsample of objects surrounded by full and transitional disks show a positive correlation between the accretion luminosity and the {H~{\sc i}} line luminosity. These two results suggest that the observed mid-IR {H~{\sc i}} lines trace gas accreting onto the star in the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of these lines can be excluded for the vast majority of full and transitional disks.We report for the first time the detection of the {H~{\sc i}}(7-6) line in eight young (< 20~Myr) debris disks. A pure chromospheric origin cannot be ruled out in these objects. If the {H~{\sc i}}(7-6) line traces accretion in these older systems, as in the case of full and transitional disks, the strength of the emission implies accretion rates lower than 10$^{-10}$M$_{\odot}$/yr. We discuss some advantages of extending accretion indicators to longer wavelengths.

The Kepler Spacecraft has discovered a large number of planets up to one-year periods and down to terrestrial sizes. While the majority of the target stars are main-sequence dwarfs of spectral type F, G, and K, Kepler covers stars with effective temperature as low as 2500 K, which corresponds to M stars. These cooler stars allow characterization of small planets near the habitable zone, yet it is not clear if this population is representative of that around FGK stars. In this paper, we calculate the occurrence of planets around stars of different spectral types as a function of planet radius and distance from the star, and show that they are significantly different from each other. We further identify two trends: First, the occurrence of Earth to Neptune-sized planets is successively higher toward later spectral types at all orbital periods probed by Kepler; Planets around M stars occur twice as frequently as around G stars, and thrice as frequently as around F stars. Second, a drop in planet occurrence is evident at all spectral types inward of a 10 day orbital period, with a plateau further out. By assigning to each spectral type a median stellar mass, we show that the distance from the star where this drop occurs is stellar mass dependent, and scales with semi-major axis as the cube root of stellar mass. By comparing different mechanisms of planet formation, trapping and destruction, we find that this scaling best matches the location of the pre-main-sequence co-rotation radius, indicating efficient trapping of migrating planets or planetary building blocks. These results demonstrate the stellar-mass dependence of the planet population, both in terms of occurrence rate and of orbital distribution. The prominent stellar-mass dependence of the inner boundary of the planet population shows that the formation or migration of planets is sensitive to the stellar parameters.

Binaries with circumbinary disks are commonly found among optically bright post-AGB stars. Although clearly linked to binary interaction processes, the formation, evolution and fate of these disks are still badly understood. Due to their compactness, interferometric techniques are required to resolve them. Here, we discuss our high-quality multiwavelength interferometric data of two prototypical yet very different post-AGB binaries, AC and 89 Herculis, as well as the modeling thereof with radiative transfer models. A detailed account of the data and models of both objects is published in three separate papers elsewhere; here we focus on comparing the modeling results for the two objects. In particular we discuss the successes and limitations of the models which were developed for protoplanetary disks around young stars. We conclude that multiwavelength high-angular-resolution observations and radiative transfer disk models are indispensible to understand these complex interacting objects and their place in the grand scheme of the (binary) evolution of low and intermediate mass stars.

The presence of disks and outflows is widespread among post-AGB binaries. In the first paper of this series, a surprisingly large fraction of optical light was found to be resolved in the 89 Her post-AGB system. The data showed this flux to arise from close to the central binary. Scattering off the inner rim of the circumbinary disk, or in a dusty outflow were suggested as two possible origins. With detailed dust radiative transfer models of the disk we aim to discriminate between these two configurations. By including Herschel/SPIRE photometry, we extend the SED such that it now fully covers UV to sub-mm wavelengths. The MCMax radiative transfer code is used to create a large grid of disk models. Our models include a self-consistent treatment of dust settling as well as of scattering. A Si-rich composition with two additional opacity sources, metallic Fe or amorphous C, are tested. The SED is fit together with mid-IR (MIDI) visibilities as well as the optical and near-IR visibilities of Paper I, to constrain the structure of the disk and in particular of its inner rim. The near-IR visibility data require a smooth inner rim, here obtained with a two-power-law parameterization of the radial surface density distribution. A model can be found that fits all the IR photometric and interferometric data well, with either of the two continuum opacity sources. Our best-fit passive models are characterized by a significant amount of mm-sized grains, which are settled to the midplane of the disk. Not a single disk model fits our data at optical wavelengths though, the reason being the opposing constraints imposed by the optical and near-IR interferometric data. A geometry in which a passive, dusty, and puffed-up circumbinary disk is present, can reproduce all the IR but not the optical observations of 89 Her. Another dusty, outflow or halo, component therefore needs to be added to the system.

The disk around the Herbig Ae/Be star HD 100546 has been extensively studied and it is one of the systems for which there are observational indications of ongoing and/or recent planet formation. However, up until now no resolved image of the millimeter dust emission or the gas has been published. We present the first resolved images of the disk around HD 100546 obtained in Band 7 with the ALMA observatory. The CO (3-2) image reveals a gas disk that extends out to 350 au radius at the 3-sigma level. Surprisingly, the 870um dust continuum emission is compact (radius <60 au) and asymmetric. The dust emission is well matched by a truncated disk with outer radius of $\approx$50 au. The lack of millimeter-sized particles outside the 60 au is consistent with radial drift of particles of this size. The protoplanet candidate, identified in previous high-contrast NACO/VLT L' observations, could be related to the sharp outer edge of the millimeter-sized particles. Future higher angular resolution ALMA observations are needed to determine the detailed properties of the millimeter emission and the gas kinematics in the inner region (<2arcsec). Such observations could also reveal the presence of a planet through the detection of circumplanetary disk material.

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15, and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67 [-0.11,+0.18]) and a significantly tapered gap edge ('round wall'), but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27 [-20,+19] mas) is notably lower than that found in our earlier H-band images (Thalmann et al. 2010). Intriguingly, we also find, at high significance, an offset of y = 69 [-25, +49] mas perpendicular to the line of nodes. If confirmed by future observations, this would imply a highly elliptical -- or otherwise asymmetric -- disk gap with an effective eccentricity of e = ~0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk, and appears brighter than the far side because of strong forward scattering.

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

A region of roughly half of the solar system scale around the star HD 100546 is largely cleared of gas and dust, in contrast to the bright outer disc. However, some material is observed in the immediate vicinity of the star. We investigate how the dust is distributed within and outside the gap, and constrain the disc geometry with mid-infrared interferometric observations using VLTI/MIDI. With baseline lengths of 40m, our long baseline observations are sensitive to the inner few AU from the star, and we combined them with observations at shorter, 15m baselines, to probe emission beyond the gap at up to 20AU from the star. We modelled the mid-infrared emission using radial temperature profiles. Our model is composed of infinitesimal concentric annuli emitting as black bodies, and it has distinct inner and outer disc components. We derived an upper limit of 0.7AU for the radial size of the inner disc, from our longest baseline data. This small dusty disc is separated from the edge of the outer disc by a large, roughly 10AU wide gap. Our short baseline data place a bright ring of emission at 11+-1AU, consistent with prior observations of the transition region between the gap and the outer disc, known as the disc wall. The inclination and position angle are constrained by our data to i=53+-8deg and PA=145+-5deg. Compared to the rim and outer disc geometry this suggests co-planarity. Brightness asymmetry is evident in both short and long baseline data, and it is unequivocally discernible from any atmospheric or instrumental effects. The origin of the asymmetry is consistent with the bright disc wall, which we find to be 1-2AU wide. The gap is cleared of micron-sized dust, but we cannot rule out the presence of larger particles and/or perturbing bodies.

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least \alpha{} >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and \alpha{} = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least \alpha{} = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap.

Context: We have analysed Herschel-PACS spectra of 32 circumstellar disks around Herbig Ae/Be and T-Tauri stars obtained within the Herschel key programme DIGIT. In this paper we focus on the 69mu emission band of the crystalline silicate forsterite. Aims: This work provides an overview of the 69mu forsterite bands in the DIGIT sample. We aim to derive the temperature and composition of the forsterite grains. With this information, constraints can be placed on the spatial distribution of the forsterite in the disk and its formation history. Methods: Position and shape of the 69mu band are used to derive the temperature and composition of the dust by comparison to laboratory spectra of that band. We combine our data with existing Spitzer IRS spectra to compare the presence and strength of the 69mu band to the forsterite bands at shorter wavelengths. Results: A total of 32 sources have been observed, 8 of them show a 69mu emission band that can be attributed to forsterite. With the exception of the T-Tauri star AS205, all of the detections are for disks associated with Herbig Ae/Be stars. Most of the forsterite grains that give rise to the 69mu bands are warm (~100-200 K) and iron-poor (less than ~2% iron). Only AB-Aur requires approximately 3-4% of iron. Conclusions: Our findings support the hypothesis that the forsterite grains form through an equilibrium condensation process at high temperatures. The connection between the strength of the 69 and 33mu bands shows that at least part of the emission in these bands originates from the same dust grains. Further, any model that explains the PACS and the Spitzer IRS observations must take the effects of a wavelength dependent optical depth into account. We find indications of a correlation of the detection rate of the 69mu band with the spectral type of the host stars. However, our sample is too small to obtain a definitive result.

Scattered light images of circumstellar disks play an important role in characterizing the planet forming environments around young stars. The characteristic size of the scattering dust grains can be estimated from the observed brightness asymmetry between the front and back side of the disk, for example using standard Mie theory. However such models often overpredict their brightness by one or two orders of magnitude, and have difficulty explaining very red disk colors. We aim to develop a dust model that explains simultaneously the observed disk surface brightness, colors and asymmetry in scattered light, focusing on constraining grain sizes. We use the 2D radiative transfer code MCMax with anisotropic scattering to explore the effects of grain size on synthetic scattered light images of circumstellar disks. We compare the results with surface brightness profiles of the protoplanetary disk HD 100546 in scattered light at wavelengths from 0.4 to 2.2 micron. We find that extreme forward scattering by micron sized particles lowers the /effective/ dust albedo and creates a faint and red disk that /appears/ only slightly forward scattering. For the outer (>~100 AU) disk of HD 100546 we derive a minimum grain size of 2.5 micron, likely present in the form of aggregates. Intermediate sized grains are too bright, whereas smaller grains are faint and scatter more isotropically, but also produce disk colors that are too blue. Observed surface brightness asymmetries alone are not sufficient to constrain the grain size in circumstellar disks. Additional information, such as the brightness and colors of the disk are needed to provide additional constraints.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

Dust settling and grain growth are the first steps in the planet-formation process in protoplanetary disks. These disks are observed around stars with different spectral types, and there are indications that the disks around lower mass stars are significantly flatter, which could indicate that they settle and evolve faster, or in a different way. We aim to test this assumption by modeling the median spectral energy distributions (SEDs) of three samples of protoplanetary disks: around Herbig stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing strength to avoid a strong observational bias from disk and stellar properties that depend on stellar mass. We generated SEDs with the radiative transfer code MCMax, using a hydrostatic disk structure and settling the dust in a self-consistent way with the alpha-prescription to probe the turbulent mixing strength. We are able to fit all three samples with a disk with the same input parameters, scaling the inner edge to the dust evaporation radius and disk mass to millimeter photometry. The Herbig stars require a special treatment for the inner rim regions, while the T-Tauri stars require viscous heating, and the brown dwarfs lack a good estimate of the disk mass because only few millimeter detections exist. We find that the turbulent mixing strength does not vary across the stellar mass range for a fixed grain size distribution and gas-to-dust ratio. Regions with the same temperature have a self-similar vertical structure independent of stellar mass, but regions at the same distance from the central star appear more settled in disks around lower mass stars. We find a relatively low turbulent mixing strength of alpha = 10^(-4) for a standard grain size distribution, but our results are also consistent with alpha = 0.01 for a grain size distribution with fewer small grains or a lower gas-to-dust ratio.

Imaging polarimetry is a powerful tool for imaging faint circumstellar material. For a correct analysis of observations we need to fully understand the effects of dust particle parameters, as well as the effects of the telescope, atmospheric seeing, and assumptions about the data reduction and processing of the observed signal. Here we study the major effects of dust particle structure, size-dependent grain settling, and instrumental properties. We performed radiative transfer modeling using different dust particle models and disk structures. To study the influence of seeing and telescope diffraction we ran the models through an instrument simulator for the ExPo dual-beam imaging polarimeter mounted at the 4.2m William Herschel Telescope (WHT). Particle shape and size have a strong influence on the brightness and detectability of the disks. In the simulated observations, the central resolution element also contains contributions from the inner regions of the protoplanetary disk besides the unpolarized central star. This causes the central resolution element to be polarized, making simple corrections for instrumental polarization difficult. This effect strongly depends on the spatial resolution, so adaptive optics systems are needed for proper polarization calibration. We find that the commonly employed homogeneous sphere model gives results that differ significantly from more realistic models. For a proper analysis of the wealth of data available now or in the near future, one must properly take the effects of particle types and disk structure into account. The observed signal depends strongly on the properties of these more realistic models, thus providing a potentially powerful diagnostic. We conclude that it is important to correctly understand telescope depolarization and calibration effects for a correct interpretation of the degree of polarization.

Forsterite is one of the crystalline dust species that is often observed in protoplanetary disks and solar system comets. Being absent in the interstellar medium, it must be produced during the disk lifetime. It can therefore serve as a tracer of dust processing and disk evolution, which can lead to a better understanding of the physical processes occurring in the disk, and possibly planet formation. However, the connection of these processes with the overall disk crystallinity remains unclear. We aim to characterize the forsterite abundance and spatial distribution in the disk of the Herbig Be star HD 100546, to investigate if a connection exists with the large disk gap. We use a 2D radiative transfer code, MCMax, to model the circumstellar dust around HD 100546. We use VISIR Q-band imaging to probe the outer disk geometry and mid-infrared features to model the spatial distribution of forsterite. The temperature-dependent shape of the 69 micron feature observed with Herschel PACS is used as a critical tool to constrain this distribution. We find a crystalline mass fraction of 40 - 60 %, located close to the disk wall between 13 and 20 AU, and possibly farther out at the disk surface. The forsterite is in thermal contact with the other dust species. We put an upper limit on the iron content of forsterite of 0.3 %. Optical depth effects play a key role in explaining the observed forsterite features, hiding warm forsterite from view at short wavelengths. The disk wall acts as a showcase: it displays a localized high abundance of forsterite, which gives rise to a high observed crystallinity, while the overall mass fraction of forsterite is a factor of ten lower.

The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will develop planetary systems. An interesting case is HD142527, for which previous observations revealed a complex circumstellar environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be stars. This suggests that the disk surrounding HD142527 is in an uncommon evolutionary stage. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig Ae/Be star HD142527. We map the composition and spatial distribution of the dust around HD142527. We analyze SEST and ATCA millimeter data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the observables satisfactorily. We find that the disk of HD142527 has three geometrically distinct components separated by a disk gap running from 30 to 130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in small grains of 10^-3 M0 and a vertical height of the inner wall of the outer disk of h = 60 AU. ....

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

With the legacy of Spitzer and current advances in (sub)mm astronomy, a large number of 'transitional' disks has been identified which are believed to contain gaps or have developped large inner holes, some filled with dust. This may indicate that complex geometries may be a key feature in disk evolution that has to be understood and modeled correctly. The disk around LkCa 15 is such a disk, with a gap ranging from ~5 - 50 AU, as identified by Espaillat et al. (2007) using 1+1D radiative transfer modelling. To fit the SED, they propose 2 possible scenarios for the inner (<5 AU) disk - optically thick or optically thin - and one scenario for the outer disk. We use the gapped disk of LkCa 15 as a showcase to illustrate the importance of 2D radiative transfer in transitional disks, by showing how the vertical dust distribution in dust-filled inner holes determines not only the radial optical depth but also the outer disk geometry. We use MCMax, a 2D radiative transfer code with a self-consistent vertical structure, to model the SED. We identify two possible geometries for the inner and outer disk, that are both different from those in Espaillat et al. (2007). An inner disk in hydrostatic equilibrium reprocesses enough starlight to fit the near infrared flux, but also casts a shadow on the inner rim of the outer disk. This requires the outer disk scale height to be large enough to rise out of the shadow. An optically thin inner disk does not cast such a shadow, and the SED can be fitted with a smaller outer disk scale height. For the dust in the inner regions to become optically thin however, the scale height would have to be so much larger than its hydrostatic equilibrium value that it effectively becomes a dust shell. It is currently unclear if a physical mechanism exists which could provide for such a configuration.

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of group II Herbig Ae stars demonstrate a lack of gas. Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately reproduces all the observables. Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a flaring gas surface and is relatively void of 'visible' dust grains. Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk.