Probing the turbulent mixing strength in protoplanetary disks across the stellar mass range: no significant variations

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.Comment: 13 pages, 16 figures, accepted by A&

Where are the Water Worlds? Identifying the Exo-water-worlds Using Models of Planet Formation and Atmospheric Evolution

Planet formation models suggest that the small exoplanets that migrate from beyond the snowline of the protoplanetary disk likely contain water-ice-rich cores ($\sim 50\%$ by mass), also known as the water worlds. While the observed radius valley of the Kepler planets is well explained with the atmospheric dichotomy of the rocky planets, precise measurements of mass and radius of the transiting planets hint at the existence of these water worlds. However, observations cannot confirm the core compositions of those planets owing to the degeneracy between the density of a bare water-ice-rich planet and the bulk density of a rocky planet with a thin atmosphere. We combine different formation models from the Genesis library with atmospheric escape models, such as photo-evaporation and impact stripping, to simulate planetary systems consistent with the observed radius valley. We then explore the possibility of water worlds being present in the currently observed sample by comparing them with the simulated planets in the mass-radius-orbital period space. We find that the migration models suggest $\gtrsim 10\%$ and $\gtrsim 20\%$ of the bare planets, i.e. planets without primordial H/He atmospheres, to be water-ice-rich around G- and M-type host stars respectively, consistent with the mass-radius distributions of the observed planets. However, most of the water worlds are predicted to be outside a period of 10 days. A unique identification of water worlds through radial velocity and transmission spectroscopy is likely to be more successful when targeting such planets with longer orbital periods.Comment: Accepted for publication in ApJ, a csv file containing analyzed observational data is attache

The Exoplanet Population Observation Simulator. I - The Inner Edges of Planetary 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.Comment: Accepted in AAS journals, code available on githu

Earths in Other Solar Systems N-body simulations: the Role of Orbital Damping in Reproducing the Kepler Planetary Systems

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.Comment: Resubmitted to ApJ. Planet formation models available online at http://eos-nexus.org/genesis-database

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

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/eposComment: Accepted in Ap

Resolved images of the protoplanetary disk around HD 100546 with ALMA

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.Comment: 6 pages, 4 figures. Accepted in ApJ

Constraints from Dust Mass and Mass Accretion Rate Measurements on Angular Momentum Transport in Protoplanetary Disks

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.Comment: See also the paper by Lodato et a

Hints for a Turnover at the Snow Line in the Giant Planet Occurrence Rate

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.Comment: 16 pages, 10 figure

An ALMA Survey of faint disks in the Chamaeleon I star-forming region: Why are some Class II disks so faint?

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.Comment: ApJ accepte