Context. The nearby and young β Pictoris system hosts a well resolved disk, a directly imaged massive giant planet orbiting at '9 au, as well as an inner planet orbiting at '2.7 au, which was... Show moreContext. The nearby and young β Pictoris system hosts a well resolved disk, a directly imaged massive giant planet orbiting at '9 au, as well as an inner planet orbiting at '2.7 au, which was recently detected through radial velocity (RV). As such, it offers several unique opportunities for detailed studies of planetary system formation and early evolution. Aims. We aim to further constrain the orbital and physical properties of β Pictoris b and c using a combination of high contrast imaging, long base-line interferometry, and RV data. We also predict the closest approaches or the transit times of both planets, and we constrain the presence of additional planets in the system. Methods. We obtained six additional epochs of SPHERE data, six additional epochs of GRAVITY data, and five additional epochs of RV data. We combined these various types of data in a single Markov-chain Monte Carlo analysis to constrain the orbital parameters and masses of the two planets simultaneously. The analysis takes into account the gravitational influence of both planets on the star and hence their relative astrometry. Secondly, we used the RV and high contrast imaging data to derive the probabilities of presence of additional planets throughout the disk, and we tested the impact of absolute astrometry. Results. The orbital properties of both planets are constrained with a semi-major axis of 9.8 ± 0.4 au and 2.7 ± 0.02 au for b and c, respectively, and eccentricities of 0.09 ± 0.1 and 0.27 ± 0.07, assuming the HIPPARCOS distance. We note that despite these low fitting error bars, the eccentricity of β Pictoris c might still be over-estimated. If no prior is provided on the mass of β Pictoris b, we obtain a very low value that is inconsistent with what is derived from brightness-mass models. When we set an evolutionary model motivated prior to the mass of β Pictoris b, we find a solution in the 10–11 MJup range. Conversely, β Pictoris c’s mass is well constrained, at 7.8 ± 0.4 MJup, assuming both planets are on coplanar orbits. These values depend on the assumptions on the distance of the β Pictoris system. The absolute astrometry HIPPARCOS-Gaia data are consistent with the solutions presented here at the 2σ level, but these solutions are fully driven by the relative astrometry plus RV data. Finally, we derive unprecedented limits on the presence of additional planets in the disk. We can now exclude the presence of planets that are more massive than about 2.5 MJup closer than 3 au, and more massive than 3.5 MJup between 3 and 7.5 au. Beyond 7.5 au, we exclude the presence of planets that are more massive than 1–2 MJup. Conclusions. Combining relative astrometry and RVs allows one to precisely constrain the orbital parameters of both planets and to give lower limits to potential additional planets throughout the disk. The mass of β Pictoris c is also well constrained, while additional RV data with appropriate observing strategies are required to properly constrain the mass of β Pictoris b. Show less
Muro-Arena, G.A.; Benisty, M.; Ginski, C.; Dominik, C.; Facchini, S.; Villenave, M.; ... ; Wildi, F. 2020
Context. Shadows in scattered light images of protoplanetary disks are a common feature and support the presence of warps or misalignments between disk regions. These warps are possibly caused by... Show moreContext. Shadows in scattered light images of protoplanetary disks are a common feature and support the presence of warps or misalignments between disk regions. These warps are possibly caused by an inclined (sub-)stellar companion embedded in the disk.Aims. We aim to study the morphology of the protoplanetary disk around the Herbig Ae star HD 139614 based on the first scattered light observations of this disk, which we model with the radiative transfer code MCMax3D.Methods. We obtained J- and H-band observations that show strong azimuthal asymmetries in polarized scattered light with VLT/SPHERE. In the outer disk, beyond similar to 30 au, a broad shadow spans a range of similar to 240 deg in position angle, in the east. A bright ring at similar to 16 au also shows an azimuthally asymmetric brightness, with the faintest side roughly coincidental with the brightest region of the outer disk. Additionally, two arcs are detected at similar to 34 and similar to 50 au. We created a simple four-zone approximation to a warped disk model of HD 139614 in order to qualitatively reproduce these features. The location and misalignment of the disk components were constrained from the shape and location of the shadows they cast.Results. We find that the shadow on the outer disk covers a range of position angles too wide to be explained by a single inner misaligned component. Our model requires a minimum of two separate misaligned zones - or a continuously warped region - to cast this broad shadow on the outer disk. A small misalignment of similar to 4 degrees between adjacent components can reproduce most of the observed shadow features.Conclusions. Multiple misaligned disk zones, potentially mimicking a warp, can explain the observed broad shadows in the HD 139614 disk. A planetary mass companion in the disk, located on an inclined orbit, could be responsible for such a feature and for the dust-depleted gap responsible for a dip in the SED. Show less
Mesa, D.; Keppler, M.; Cantalloube, F.; Rodet, L.; Charnay, B.; Gratton, R.; ... ; Wildi, F. 2019
