We present high-resolution millimeter continuum imaging of the disk surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disk system. The system has ...extended mm emission on which are superposed three prominent annular gaps at radii ∼13, 39, and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disk; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exoplanets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius. Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage.
Abstract
It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly ...form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragment's interaction with the disc, but also with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment–fragment interactions while the disc is still present. We find that fragment–fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations – this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment–fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free-floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.
ABSTRACT
Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs, and companion stars. Competition between three processes ...– migration, gas accretion, and tidal disruption – establishes what grows from a given seed. Numerical simulations and population synthesis calculations published to date, however, do not always agree on the outcome. Here, we investigate if the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG, and FARGO give the same answer when faced with the same migrating clump setup. Four test runs with varying assumptions about the initial clump mass and gas accretion on to it are performed. We find that the codes disagree in the clump migration rate by between 10 per cent to ∼50 per cent, depending on the test, but always arrive in the same qualitative picture. Specifically, with gas accretion turned off, planets migrate through the whole effective computational domain. In contrast, for the run with the most massive seed and gas accretion on, the planet opens a deep gap and stalls at separation of order 80 AU. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. We also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations.
Context. Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet ...formation. Dissimilarities in the spatial distribution of small (μm−size) and large (mm−size) dust grains have recently been pointed out. Aims. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks (with particular emphasis on the inner 10−50 AU), thus increasing the scarce sample of high-resolution images of these objects. Methods. We obtained VLT/NACO near-IR high-resolution polarimetric differential imaging observations of SAO 206462 (HD 135344B). This technique allows one to image the polarized scattered light from the disk without any occulting mask and to reach an inner working angle of ~0.1″. Results. A face-on disk is detected in H and Ks bands between 0.1″ and 0.9″. No significant differences are seen between the H and Ks images. In addition to the spiral arms, these new data allow us to resolve for the first time an inner disk cavity for small dust grains. The cavity size (≃28 AU) is much smaller than what is inferred for large dust grains from (sub-)mm observations (39 to 50 AU). This discrepancy cannot be ascribed to any resolution effect. Conclusions. The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas (and, thus, the small grains) at 28 AU, and generating a pressure bump at larger radii (39 AU), which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet (with a mass between 5 and 15 MJ) at 17 to 20 AU from the star is consistent with the observed cavity sizes.
Context. In 2D-simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the ...gravitational potential or force. The softening prescription used in 2D needs to be adjusted properly to correctly resemble the realistic case of vertically extended 3D disks. Aims. We analyze the embedded planet and the self-gravity case and provide a method to evaluate the required smoothing in 2D simulations of thin disks. Methods. Starting from the averaged hydrodynamic equations and using a vertically isothermal disk model, we calculated the force to be used in 2D simulations. We compared our results to the often used Plummer form of the potential, which runs as ∝ 1/(r2 + ϵ2)1/2. For that purpose we computed the required smoothing length ϵ as a function of distance r to the planet or to a disk element within a self-gravitating disk. Results. We find that for longer distances ϵ is determined solely by the vertical disk thickness H. For the planet case we find that outside r ≈ H a value of ϵ = 0.7H describes the averaged force very well, while in the self-gravitating disk the value needs to be higher, ϵ = 1.2H. For shorter distances the smoothing needs to be reduced significantly. Comparing torque densities of 3D and 2D simulations we show that the modification to the vertical density stratification as induced by an embedded planet needs to be taken into account to obtain agreeing results. Conclusions. It is very important to use the correct value of ϵ in 2D simulations to obtain a realistic outcome. In disk fragmentation simulations the choice of ϵ can determine whether a disk will fragment or not. Because a wrong smoothing length can change even the direction of migration, it is very important to include the effect of the planet on the local scale height in 2D planet-disk simulations. We provide an approximate and fast method for this purpose that agrees very well with full 3D simulations.
Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet formation. ...Dissimilarities in the spatial distribution of small ( mu m-size) and large (mm-size) dust grains have recently been pointed out. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks, thus increasing the scarce sample of high-resolution images of these objects. The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas at 28 AU, and generating a pressure bump at larger radii, which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet at 17 to 20 AU from the star is consistent with the observed cavity sizes.
We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD 219666 (M⋆ = 0.92 ± 0.03 M⊙, R⋆ = 1.03 ± 0.03 R⊙, τ⋆ = 10 ± 2 Gyr). With a ...mass of Mb = 16.6 ± 1.3 M⊕, a radius of Rb = 4.71 ± 0.17 R⊕, and an orbital period of Porb ≃ 6 days, HD 219666 b is a new member of a rare class of exoplanets: the hot-Neptunes. The Transiting Exoplanet Survey Satellite (TESS) observed HD 219666 (also known as TOI-118) in its Sector 1 and the light curve shows four transit-like events, equally spaced in time. We confirmed the planetary nature of the candidate by gathering precise radial-velocity measurements with the High Accuracy Radial velocity Planet Searcher (HARPS) at ESO 3.6 m. We used the co-added HARPS spectrum to derive the host star fundamental parameters (Teff = 5527 ± 65 K, log g⋆ = 4.40 ± 0.11 (cgs), Fe/H= 0.04 ± 0.04 dex, log R′HK $\log R^{\prime}_{\textrm{HK}}$ log R HK ′ = −5.07 ± 0.03), as well as the abundances of many volatile and refractory elements. The host star brightness (V = 9.9) makes it suitable for further characterisation by means of in-transit spectroscopy. The determination of the planet orbital obliquity, along with the atmosphericmetal-to-hydrogen content and thermal structure could provide us with important clues on the formation mechanisms of this class of objects.
Context. Within the collision growth scenario for planetesimal formation, the growth step from centimetre-sized pre-planetesimals to kilometre-sized planetesimals remains unclear. The formation of ...larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. Aims. We wish to develop a new classification scheme broad enough to encompass all events with sticking, bouncing, and fragmentation contributions, accurate enough to capture the important collision outcome nuances, and at the same time simple enough to be implementable in global dust coagulation simulations. We furthermore wish to demonstrate the reliability of our numerical smoothed particle hydrodynamics (SPH) model and the applicability of our new collision outcome classification to previous results as well as our simulation results. Methods. We propose and apply a scheme based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. On the basis of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous SiO2 dust aggregates. Results. We demonstrate that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. Intermediate categories are found in our simulations that are difficult to map to existing qualitative categorisations. We show that the four-population model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. In addition, the numerical porosity model successfully passes another benchmark test: the correct simulation of the entire list of collision outcome types yielded by laboratory experiments. As a demonstration of the applicability and the power of the four-population model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.
ABSTRACT We perform two-dimensional hydrodynamical simulations to quantitatively explore the torque balance criterion for gap-opening (as formulated by Crida et al.) in a variety of disks when ...considering a migrating planet. We find that even when the criterion is satisfied, there are instances when planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as an additional condition on the gap-opening timescale, , versus the crossing time, , i.e., the time it takes the planet to cross the region which it is carving out. While this point has been briefly made in the previous literature, our results quantify it for a range of protoplanetary disk properties and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has important implications for the survival of planets formed by core accretion in low mass disks as well as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particularly important for planets with intermediate masses susceptible to Type III-like migration. For some observed transition disks or disks with gaps, we expect that estimates on the potential planet masses based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of this additional timescale criterion theoretical studies may find a reduced planet survivability or that planets may migrate further inwards before opening a gap.
We present long-baseline Atacama Large Millimeter/submillimeter Array observations of the 870 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 Jy beam−1, and reveal an asymmetric ring between ∼20 and 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 (<2 au 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.44 M⊕ and 0.44 au in the optically thin and thick regimes, respectively, for the case of the previously directly imaged protoplanet candidate at ∼55 au (HD 100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD 100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 MJ 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 could 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.