ABSTRACT
Many stars are in binaries or higher order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in ...multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper, we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here, we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time-scales.
Abstract
Spiral arms have been observed in nearly a dozen protoplanetary discs in near-infrared scattered light and recently also in the submillimetre continuum. While one of the most compelling ...explanations is that they are driven by planetary or stellar companions, in all but one cases such companions have not yet been detected and there is even ambiguity on whether the planet should be located inside or outside the spirals. Here, we use 3D hydrodynamic simulations to study the morphology of spiral density waves launched by embedded planets taking into account the vertical temperature gradient, a natural consequence of stellar irradiation. Our simulations show that the pitch angle of the spirals in thermally stratified discs is the lowest in the disc mid-plane and increases towards the disc surface. We combine the hydrodynamic simulations with 3D radiative transfer calculations to predict that the pitch angle of planetary spirals observed in the near-infrared is higher than in the submillimetre. We also find that in both cases the spirals converge towards the planet. This provides a new powerful observational method to determine if the perturbing planet is inside or outside the spirals, as well as map the thermal stratification of the disc.
We investigate the minimum planet mass that produces observable signatures in infrared scattered light and submillimetre (submm) continuum images and demonstrate how these images can be used to ...measure planet masses to within a factor of about 2. To this end, we perform multi-fluid gas and dust simulations of discs containing low-mass planets, generating simulated observations at 1.65, 10 and 850 μm. We show that the minimum planet mass that produces a detectable signature is ∼15 M⊕: this value is strongly dependent on disc temperature and changes slightly with wavelength (favouring the submm). We also confirm previous results that there is a minimum planet mass of ∼20 M⊕ that produces a pressure maximum in the disc: only planets above this threshold mass generate a dust trap that can eventually create a hole in the submm dust. Below this mass, planets produce annular enhancements in dust outwards of the planet and a reduction in the vicinity of the planet. These features are in steady state and can be understood in terms of variations in the dust radial velocity, imposed by the perturbed gas pressure radial profile, analogous to a traffic jam. We also show how planet masses can be derived from structure in scattered light and submm images. We emphasize that simulations with dust need to be run over thousands of planetary orbits so as to allow the gas profile to achieve a steady state and caution against the estimation of planet masses using gas-only simulations.
Abstract
Gas giants’ early (≲5 Myr) orbital evolution occurs in a disc losing mass in part to photoevaporation driven by high energy irradiance from the host star. This process may ultimately ...overcome viscous accretion to disperse the disc and halt migrating giants by starving their orbits of gas, imprinting on giant planet separations in evolved systems. Inversion of this distribution could then give insight into whether the stellar FUV, EUV or X-ray flux dominates photoevaporation, constraining planet formation and disc evolution models. We use a 1D hydrodynamic code in population syntheses for gas giants undergoing Type II migration in a viscously evolving disc subject to either a primarily FUV, EUV or X-ray flux from a pre-solar T Tauri star. The photoevaporative mass loss profile’s unique peak location and width in each energetic regime produces characteristic features in the distribution of giant separations: a severe dearth of ≲2 MJ planets interior to 5 au in the FUV scenario, a sharp concentration of ≲3 MJ planets between ≈1.5–2 au in the EUV case and a relative abundance of ≈2–3.5 MJ giants interior to 0.5 au in the X-ray model. These features do not resemble the observational sample of gas giants with mass constraints, although our results do show some weaker qualitative similarities. We thus assess how the differing photoevaporative profiles interact with migrating giants and address the effects of large model uncertainties as a step to better connect disc models with trends in the exoplanet population.
Abstract
We report the robust detection of coherent, localized deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disk around HD 163296. The ...analysis is performed using the
discminer
channel map modeling framework on
12
CO
J
= 2–1 DSHARP data. Not only orbital radius but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at
R
= 94 ± 6 au,
ϕ
= 50° ± 3° (P94), is near the center of one of the gaps in dust continuum emission and is consistent with a planet mass of 1
M
Jup
. The other planet, located at
R
= 261 ± 4 au,
ϕ
= 57° ± 1° (P261), is in the region where a velocity kink was previously observed in
12
CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disk and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles, we detect gas gaps at
R
= 38, 88, and 136 au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at
R
= 88 au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the
discminer
is capable of finding localized, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary disks.
ABSTRACT
Protoplanetary disc surveys conducted with Atacama Large Millimetre Array (ALMA) are measuring disc radii in multiple star-forming regions. The disc radius is a fundamental quantity to ...diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous protoplanetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust mass increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly on to the star, the residual dust retained on Myr time-scales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc flux (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time.
Abstract
We explore how measurements of protoplanetary disc masses and accretion rates provided by surveys of star-forming regions can be analysed via the dimensionless accretion parameter, which we ...define as the product of the accretion rate and stellar age divided by the disc mass. By extending and generalizing the study of Jones et al., we demonstrate that this parameter should be less than or of order unity for a wide range of evolutionary scenarios, rising above unity only during the final stages of outside-in clearing by external photoevaporation. We use this result to assess the reliability of disc mass estimates derived from CO isotopologues and sub-mm continuum emission by examining the distribution of accretion efficiencies in regions that are not subject to external photoevaporation. We find that while dust-based mass estimates produce results compatible with theoretical expectations assuming a canonical dust-to-gas ratio, the systematically lower CO-based estimates yield accretion efficiencies significantly above unity in contrast with the theory. This finding provides additional evidence that CO-based disc masses are an underestimate, in line with arguments that have been made on the basis of chemical modelling of relatively small samples. On the other hand, we demonstrate that dust-based mass estimates are sufficiently accurate to reveal distinctly higher accretion efficiencies in the Trapezium cluster, where this result is expected, given the evident importance of external photoevaporation. We therefore propose the dimensionless accretion parameter as a new diagnostic of external photoevaporation in other star-forming regions.
ABSTRACT
The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of ...viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscous evolution for the MHD wind case. The transport of angular momentum and mass induced by the wind is parametrized by an α-like parameter and by the magnetic lever arm parameter λ. Extensions of the paradigmatic Lynden-Bell and Pringle similarity solutions to the wind case are presented. We show that wind-driven accretion leads to a steeper decrease in the disc mass and accretion rate than in viscous models due to the absence of disc spreading. If the decline of the magnetic field strength is slower than that of the gas surface density, the disc is dispersed after a finite time. The evolution of the disc in the $\dot{M}_*-M_{D}$ plane is sensitive to the wind and turbulence parameters. A disc population evolving under the action of winds can exhibit a correlation between $\dot{M}_*$ and MD depending on the initial conditions. The simplified framework proposed in this work opens to a new avenue to test the effectiveness of wind-driven accretion from the observed disc demographics and constitutes an important step to include wind-driven accretion in planet population synthesis models.
ABSTRACT
When imaged at high resolution, many protoplanetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local ...maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper, we present a method to measure the dust–gas coupling α/St and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission line data cubes. As a proof of concept, we then apply the method to two discs with prominent substructure, HD 163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good (α/St ∼ 0.1). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the α turbulent parameter (α ∼ 10−2). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD 163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1 mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.
Context.
Recent years have seen building evidence that planet formation starts early, in the first ~0.5 Myr. Studying the dust masses available in young disks enables us to understand the origin of ...planetary systems given that mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones.
Aims.
We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets.
Methods.
We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.1–1.3 mm) continuum observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA)
Ka
-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities measured in the image plane.
Results.
We find a strong linear correlation between the ALMA and VLA fluxes, demonstrating that emission at both wavelengths is dominated by dust emission. For a subsample of optically thin sources, we find a median spectral index of 2.5 from which we derive the dust opacity index
β
= 0.5, suggesting significant dust growth. Comparison with ALMA surveys of Orion shows that the Class I dust disk mass distribution between the two regions is similar, but that the Class 0 disks are more massive in Perseus than those in Orion. Using the DIANA opacity model including large grains, with a dust opacity value of
κ
9 mm
= 0.28 cm
2
g
−1
, the median dust masses of the embedded disks in Perseus are 158
M
⊕
for Class 0 and 52
M
⊕
for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47
M
⊕
and 12
M
⊕
for Class 0 and Class I, respectively, obtained using the maximum dust opacity value
κ
1.3 mm
= 2.3 cm
2
g
−1
. The dust masses of young Class 0 and I disks are larger by at least a factor of ten and three, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions.
Conclusions.
The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of ~15%. A higher efficiency of ~30% is necessary if the planet formation is set to start in Class I disks.