Super-Earths in the TW Hya disc

We test the hypothesis that the sub-millimetre thermal emission and scattered light gaps seen in recent observations of TW Hya are caused by planet-disc interactions. We perform global three-dimensional dusty smoothed particle hydrodynamics simulations, comparing synthetic observations of our models with dust thermal emission, CO emission and scattered light observations. We find that the dust gaps observed at 24 au and 41 au can be explained by two super-Earths ($\sim 4 \mathrm{M}_{\oplus}$). A planet of approximately Saturn-mass can explain the CO emission and the depth and width of the gap seen in scattered light at 94 au. Our model produces a prominent spiral arm while there are only hints of this in the data. To avoid runaway growth and migration of the planets we require a disc mass of $\lesssim 10^{-2}\,\mathrm{M}_{\odot}$ in agreement with CO observations but 10$-$100 times lower than the estimate from HD line emission.Comment: 6 pages, 5 figures, accepted for publication in MNRA

The Foggy Disks Surrounding Herbig Ae Stars: a Theoretical Study of the H2O Line Spectra

Water is a key species in many astrophysical environments, but it is particularly important in proto-planetary disks. So far,observations of water in these objects have been scarce, but the situation should soon change thanks to the Herschel satellite. We report here a theoretical study of the water line spectrum of a proto-planetary disk surrounding Ae stars. We show that several lines will be observable with the HIFI instrument onboard the Herschel Space Observatory. We predict that some maser lines could also be observable with ground telescopes and we discuss how the predictions depend not only on the adopted physical and chemical model but also on the set of collisional coefficients used and on the H2 ortho to para ratio through its effect on collisional excitation. This makes the water lines observations a powerful, but dangerous -if misused- diagnostic tool.Comment: Accepted for publication in ApJ Letter

On the stratified dust distribution of the GG Tau circumbinary ring

Our objective is to study the vertical dust distribution in the circumbinary ring of the binary system GG Tau and to search for evidence of stratification, one of the first steps expected to occur during planet formation. We present a simultaneous analysis of four scattered light images spanning a range of wavelength from 800 nm to 3800 nm and compare them with (i) a parametric prescription for the vertical dust stratification, and (ii) with the results of SPH bi-fluid hydrodynamic calculations. The parametric prescription and hydrodynamical calculations of stratification both reproduce the observed brightness profiles well. These models also provide a correct match for the observed star/ring integrated flux ratio. Another solution with a well-mixed, but ``exotic'', dust size distribution also matches the brightness profile ratios but fails to match the star/ring flux ratio. These results give support to the presence of vertical stratification of the dust in the ring of GG Tau and further predict the presence of a radial stratification also.Comment: 9 pages, 11 figures. Accepted for publication in A&

Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.Comment: 15 pages, 10 figures, accepted for publication in A&

The accumulation and trapping of grains at planet gaps: effects of grain growth and fragmentation

We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap's detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.Comment: Accepted for publication in Planetary and Space Science. 13 pages, 4 figure

The Taurus Boundary of Stellar/Substellar (TBOSS) Survey I: far-IR disk emission measured with Herschel

With Herschel/PACS 134 low mass members of the Taurus star-forming region spanning the M4-L0 spectral type range and covering the transition from low mass stars to brown dwarfs were observed. Combining the new Herschel results with other programs, a total of 150 of the 154 M4-L0 Taurus members members have observations with Herschel. Among the 150 targets, 70um flux densities were measured for 7 of the 7 ClassI objects, 48 of the 67 ClassII members, and 3 of the 76 ClassIII targets. For the detected ClassII objects, the median 70um flux density level declines with spectral type, however, the distribution of excess relative to central object flux density does not change across the stellar/substellar boundary in the M4-L0 range. Connecting the 70um TBOSS values with the results from K0-M3 ClassII members results in the first comprehensive census of far-IR emission across the full mass spectrum of the stellar and substellar population of a star-forming region, and the median flux density declines with spectral type in a trend analogous to the flux density decline expected for the central objects. SEDs were constructed for all TBOSS targets covering the optical to far-IR range and extending to the submm/mm for a subset of sources. Based on an initial exploration of the impact of different physical parameters; inclination, scale height and flaring have the largest influence on the PACS flux densities. From the 24um to 70um spectral index of the SEDs, 5 new candidate transition disks were identified. The steep 24um to 70um slope for a subset of 8 TBOSS targets may be an indication of truncated disks in these systems.Two examples of mixed pair systems that include secondaries with disks were measured. Finally, comparing the TBOSS results with a Herschel study of Ophiuchus brown dwarfs reveals a lower fraction of disks around the Taurus substellar population.Comment: 64 pages, 33 figures, 12 tables, accepted for publication in A&

Nonlinear system-identification of the filling phase of a wet-clutch system

The work presented illustrates how the choice of input perturbation signal and experimental design improves the derived model of a nonlinear system, in particular the dynamics of a wet-clutch system. The relationship between the applied input current signal and resulting output pressure in the filling phase of the clutch is established based on bandlimited periodic signals applied at different current operating points and signals approximating the desired filling current signal. A polynomial nonlinear state space model is estimated and validated over a range of measurements and yields better fits over a linear model, while the performance of either model depends on the perturbation signal used for model estimation

Evidence for a circumplanetary disk around protoplanet PDS 70 b

We present the first observational evidence for a circumplanetary disk around the protoplanet PDS~70~b, based on a new spectrum in the $K$ band acquired with VLT/SINFONI. We tested three hypotheses to explain the spectrum: Atmospheric emission from the planet with either (1) a single value of extinction or (2) variable extinction, and (3) a combined atmospheric and circumplanetary disk model. Goodness-of-fit indicators favour the third option, suggesting circumplanetary material contributing excess thermal emission --- most prominent at $\lambda \gtrsim 2.3 \mu$m. Inferred accretion rates ($\sim 10^{-7.8}$--$10^{-7.3} M_J$ yr$^{-1}$) are compatible with observational constraints based on the H$\alpha$ and Br$\gamma$ lines. For the planet, we derive an effective temperature of 1500--1600 K, surface gravity $\log(g)\sim 4.0$, radius $\sim 1.6 R_J$, mass $\sim 10 M_J$, and possible thick clouds. Models with variable extinction lead to slightly worse fits. However, the amplitude ($\Delta A_V \gtrsim 3$mag) and timescale of variation ($\lesssim$~years) required for the extinction would also suggest circumplanetary material.Comment: 8 pages, 2 figures, 1 table. This is a pre-copyedited, author-produced PDF of an article accepted for publication in ApJL on 2019 May 1

Where is the warm H2 ? A search for H2 emission from disks around Herbig Ae/Be stars

Mid-IR emission lines of H2 are useful probes to determine the mass of warm gas present in the surface layers of disks. Numerous observations of Herbig Ae/Be stars (HAeBes) have been performed, but only 2 detections of mid-IR H2 toward HD97048 and AB Aur have been reported. We aim at tracing the warm gas in the disks of 5 HAeBes with gas-rich environments and physical characteristics close to those of AB Aur and HD97048, to discuss whether the detections toward these 2 objects are suggestive of peculiar conditions for the gas. We search for the H2 S(1) emission line at 17.035 \mu\m with VISIR, and complemented by CH molecule observations with UVES. We gather the H2 measurements from the literature to put the new results in context and search for a correlation with some disk properties. None of the 5 VISIR targets shows evidence for H2 emission. From the 3sigma upper limits on the integrated line fluxes we constrain the amount of optically thin warm gas to be less than 1.4 M_Jup in the disk surface layers. There are now 20 HAeBes observed with VISIR and TEXES instruments to search for warm H2, but only two detections (HD97048 and AB Aur) were made so far. We find that the two stars with detected warm H2 show at the same time high 30/13 \mu\m flux ratios and large PAH line fluxes at 8.6 and 11.3 \mu\m compared to the bulk of observed HAeBes and have emission CO lines detected at 4.7 \mu\m. We detect the CH 4300.3A absorption line toward both HD97048 and AB Aur with UVES. The CH to H2 abundance ratios that this would imply if it were to arise from the same component as well as the radial velocity of the CH lines both suggest that CH arises from a surrounding envelope, while the detected H2 would reside in the disk. The two detections of the S(1) line in the disks of HD97048 and AB Aur suggest either peculiar physical conditions or a particular stage of evolution.Comment: accepted for publication in A&A : 10 pages, 6 figure

Planet gaps in the dust layer of 3D protoplanetary disks. II. Observability with ALMA

[Abridged] Aims: We provide predictions for ALMA observations of planet gaps that account for the specific spatial distribution of dust that results from consistent gas+dust dynamics. Methods: In a previous work, we ran full 3D, two-fluid Smoothed Particle Hydrodynamics (SPH) simulations of a planet embedded in a gas+dust T Tauri disk for different planet masses and grain sizes. In this work, the resulting dust distributions are passed to the Monte Carlo radiative transfer code MCFOST to construct synthetic images in the ALMA wavebands. We then use the ALMA simulator to produce images that include thermal and phase noise for a range of angular resolutions, wavelengths, and integration times, as well as for different inclinations, declinations and distances. We also produce images which assume that gas and dust are well mixed with a gas-to-dust ratio of 100 to compare with previous ALMA predictions, all made under this hypothesis. Results: Our findings clearly demonstrate the importance of correctly incorporating the dust dynamics. We show that the gap carved by a 1 M_J planet orbiting at 40 AU is visible with a much higher contrast than the well-mixed assumption would predict. In the case of a 5 M_J planet, we clearly see a deficit in dust emission in the inner disk, and point out the risk of interpreting the resulting image as that of a transition disk with an inner hole if observed in unfavorable conditions. Planet signatures are fainter in more distant disks but declination or inclination to the line-of-sight have little effect on ALMA's ability to resolve the gaps. Conclusions: ALMA has the potential to see signposts of planets in disks of nearby star-forming regions. We present optimized observing parameters to detect them in the case of 1 and 5 M_J planets on 40 AU orbits.Comment: 15 pages, 21 figures, accepted by Astronomy & Astrophysics, a higher resolution version of the paper is available at http://www-obs.univ-lyon1.fr/labo/perso/jean-francois.gonzalez/Papers/Gaps_ALMA.pd