Issue 
A&A
Volume 638, June 2020



Article Number  A1  
Number of page(s)  7  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/202037983  
Published online  29 May 2020 
Planet formation by pebble accretion in ringed disks
Département Lagrange, University of Nice – Sophia Antipolis, CNRS,
Observatoire de la Côte d’Azur,
Nice,
France
email: morby@obsnice.fr
Received:
19
March
2020
Accepted:
9
April
2020
Context. Pebble accretion is expected to be the dominant process for the formation of massive solid planets, such as the cores of giant planets and superEarths. So far, this process has been studied under the assumption that dust coagulates and drifts throughout the full protoplanetary disk. However, observations show that many disks are structured in rings that may be due to pressure maxima, preventing the global radial drift of the dust.
Aims. We aim to study how the pebbleaccretion paradigm changes if the dust is confined in a ring.
Methods. Our approach is mostly analytic. We derived a formula that provides an upper bound to the growth of a planet as a function of time. We also numerically implemented the analytic formulæ to compute the growth of a planet located in a typical ring observed in the DSHARP survey, as well as in a putative ring rescaled at 5 AU.
Results. Planet Type I migration is stopped in a ring, but not necessarily at its center. If the entropydriven corotation torque is desaturated, the planet is located in a region with low dust density, which severely limits its accretion rate. If the planet is instead near the ring’s center, its accretion rate can be similar to the one it would have in a classic (ringless) disk of equivalent dust density. However, the growth rate of the planet is limited by the diffusion of dust in the ring, and the final planet mass is bounded by the total ring mass. The DSHARP rings are too far from the star to allow the formation of massive planets within the disk’s lifetime. However, a similar ring rescaled to 5 AU could lead to the formation of a planet incorporating the full ring mass in less than 1/2 My.
Conclusions. The existence of rings may not be an obstacle to planet formation by pebbleaccretion. However, for accretion to be effective, the resting position of the planet has to be relatively near the ring’s center, and the ring needs to be not too far from the central star. The formation of planets in rings can explain the existence of giant planets with core masses smaller than the socalled pebble isolation mass.
Key words: planets and satellites: formation / protoplanetary disks / planetdisk interactions
© A. Morbidelli 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The formation of massive planets (cores of giant planets, superEarths) is not yet fully elucidated. The classic model of oligarchic growth in a disk of planetesimals has difficulties producing such massive bodies within the disk’s lifetime, particularly at distances typical of the giant planets of the solar system or warm jupiters around other stars (Levison et al. 2010; Morbidelli et al. 2015; Johansen & Lambrechts 2017; Johansen & Bitsch 2019). Oligarchic growth in planetesimal disks is also particularly inefficient if the initial planetesimals are mostly ~100 km in size (Fortier et al. 2013), as predicted by the streaming instability model (Johansen et al. 2015; Simon et al. 2017) and suggested by observations of the sizefrequency distribution of the remaining solar system planetesimals (for asteroids see e.g., Morbidelli et al. 2009, Delbo’ et al. 2017; for Kuiperbelt objects see e.g., Morbidelli & Nesvorný 2020).
For these reasons, a new paradigm for planet formation has been developed in the last decade, dubbed pebble accretion (Ormel & Klahr 2010; Lambrechts & Johansen 2012, 2014; Lambrechts et al. 2014; Levison et al. 2015; Ida et al. 2016; Ormel 2017, to quote just a few). In this paradigm, planets grow by accreting small solid material (dust grains, pebblesized objects) as this material drifts radially in the disk due to aerodynamic drag. If the flux of material is sufficiently large, planetary cores of 10–20 Earth masses can form within the lifetime of the disk.
Nevertheless, the models of pebble accretion developed so far are quite simplistic, because they assume a continuous flux of pebbles as if their initial reservoir had an infinite radial extension. However, protoplanetary disks are not infinitely wide. The median observed sizes (in gas) of ClassI and ClassII disks are 100 and 200 AU, respectively, although some exceptional ones can extend to ~1000 AU (Najita & Bergin 2018). Due to observational biases, the real medians are certainly smaller. There is evidence that the protoplanetary disk of the solar system did not exceed 80 AU in radial extension (Kretke et al. 2012). Given the expected rates of growth and radial drift of dust (Brauer et al. 2008; Birnstiel et al. 2010, 2016), in disks of these sizes the flux of pebbles is expected to last much less than the lifetime of the disk. For instance, Sato et al. (2016) estimated that the flux of pebbles should decay sharply after a time $t~2\times {10}^{5}{({r}_{\text{gas}}/100\text{\hspace{0.17em}}\text{AU})}^{3/2}$ yr, where r_{gas} is the radius of the gas disk. As a corollary, one would expect that protoplanetary disks at several AUs from the star become dust poor very rapidly because of dust growth and drift (Birnstiel et al. 2010; Rosotti et al. 2019), but observations show that disks remain dust rich (compared to gas) for at least a few My (Manara et al. 2016), given the good correlation existing between the gas accretion rate onto the central star and the disk mass in dust. There is no evidence that disks are systematically smaller in dust than gas, because the observed ratios (typically r_{gas}∕r_{dust} ~ 2 – Najita & Bergin 2018; Ansdell et al. 2018) could simply be the consequence of optical depth effects (Trapman et al. 2019). For balance, it should be noted that the size of disks in gas is also very uncertain and may depend on the tracer that is used (Anderson et al. 2019).
Another issue with the pebble accretion paradigm is that it is inefficient. To form a core of 10 Earth masses starting from a Moonmass seed, more than 100 Earth masses of pebbles have to pass through the planetary orbit, most of which are not accreted by the growing planet (Morbidelli et al. 2015; Bitsch et al. 2019). The same is true for the formation of a system of superEarths (Lambrechts et al. 2019). These unaccreted pebbles/dust should accumulate at the inner edge of the disk, where the transition from a turbulent to a more quiescent gas is expected to create a local maximum in the pressureradial distribution (Flock et al. 2016, 2017). There, the accumulated material should form massive planets (Flock et al. 2019). The corollary is that every star should have closein massive planets, which is not observed.
These problems could be solved if the dust grows much more slowly than expected in coagulation models, and so drifts toward the star on a longer timescale (Johansen et al. 2019). However, in this case the planets would also grow more slowly (because the pebbles are smaller). Moreover, observations of disks at mmwavelength suggest that mmsize dust does exist in disks. Another possibility is that there are impediments to dustdrift in protoplanetary disks, possibly due to the presence of multiple pressure bumps (Pinilla et al. 2012). Recent highresolution observations with ALMA do indeed show that most protoplanetary disks have ringlike structures (ALMA Partnership 2015; Andrews et al. 2018; Long et al. 2018; Dullemond et al. 2018), indicating that bottlenecks for dustdrift exist at multiple locations.
The origin of these rings is still debated. Magnetohydrodynamical (MHD) simulations show that the surface density distribution of disks can become a radial wavy function (Bethune et al. 2017), leading to the formation of multiple pressure bumps and hence of dust rings (Riols et al. 2019). More subtle dustgas instabilities can also generate dust rings (Tominaga et al. 2019). On the other hand, massive planets can open gaps in the gas distribution and generate pressure bumps (e.g., Zhang et al. 2018; Weber et al. 2017; Yang & Zhu 2020; WafflardFernandez & Baruteau 2020). Thus, whether rings are a prerequisite for planet formation or a consequence of planet formation is not yet known. Similarly, for the solar system, there is evidence that the dust from the outer disk did not penetrate into the inner disk in large quantities, leading to the observed dichotomy in the isotopic properties of inner and outer solar system planetesimals (Kruijer et al. 2017). Whether this separation between the inner and outer disk’s dust was due to the existence of a pressure bump in the disk (Brasser & Mojzsis 2020) or to the formation of Jupiter (Kruijer et al. 2017) is under debate.
This short paper places itself in the first hypothesis (rings exist before planet formation takes place) and explores how the pebble accretion paradigm changes in the presence of rings. If rings are a consequence of planet formation, what follows can be regarded as a description of how a second generation of planets can grow, from the rings induced by the formation of the first generation. Two important assumptions of this work are that (i) the rings are due to pressure maxima, as advocated by Dullemond et al. (2018), and (ii) they are longlived structures which, frankly, is not known. If the rings are ephemeral, or just the result of a modulation of the dust drift speed, they probably don’t have a significant impact on planet accretion.
There are three big changes in the pebbleaccretion narrative if the dust is confined in rings. Firstly, there is no net radial drift of dust. Dust diffuses back and forth within the ring, which gives the planet multiple chances to accrete it, alleviating the inefficiency of pebble accretion upon a single passage of a dust grain through the planet’s orbit. Secondly, the reservoir of solid mass available for planet growth is finite, limited by the mass of the ring. Thirdly, there is no differential azimuthal velocity between the planet and gas/dust near the pressure maximum. Section 2 analyzes all of these aspects, leading to an analytic upper bound of the planet’s growth as a function of ring mass, turbulent diffusion, and the dust’s Stokes number. Section 3 then presents a more detailed computation of the growth of a planet in the rings observed in the DSHARP survey, discussed in Dullemond et al. (2018). The results are discussed in Sect. 4.
2 Dust dynamics in a ring and its accretion by a planet
The formation of a dust ring in the disk is most likely due to the appearance of a pressure bump in the disk, whatever the cause for this pressure bump (e.g., MHD instabilities, dustgas instabilities, the existence of other planets). As in Dullemond et al. (2018), we assume that the radial pressure profile is Gaussian, namely: $$p(r)={p}_{0}\mathrm{exp}\left(\frac{{(r{r}_{0})}^{2}}{2{w}^{2}}\right),$$(1)
where r_{0} is the center of the pressure maximum and w its width.
Because of gas drag, the dust is attracted toward the pressure maximum, with a radial speed that is $${v}_{r}(r)=\frac{{H}^{2}\Omega \tau}{{w}^{2}(1+{\tau}^{2})}(r{r}_{0})$$(2)
(see Formula 44 in Dullemond et al. 2018), where H is the pressure scaleheight of the disk, Ω is the keplerian orbital frequency, and τ is the dust’s Stokes number^{1}.
On the other hand, the dust also undergoes turbulent diffusion. The diffusion coefficient for the dust is (Youdin & Lithwick 2007): $${D}_{\text{d}}=\frac{D}{1+{\tau}^{2}},$$(3)
where D is the turbulent diffusion in the gas, (usually assumed to be equal to its viscosity) which, following the αprescription of Shakura & Sunyaev (1973), gives D = αH^{2}Ω, where α is the turbulent parameter. Observations of the dust distribution in disks suggest that α = 10^{−4}–10^{−3} (Pinte et al. 2016; Dullemond et al. 2018).
The balance between diffusion (which tends to disperse the dust) and the drag force (which tends to bring the dust back towards the pressure maximum) gives the dust the steady state radial distribution: $${\Sigma}_{d}(r)={\Sigma}_{d,0}\mathrm{exp}\left(\frac{{(r{r}_{0})}^{2}}{2{w}_{d}^{2}}\right),$$(4)
with a width w_{d} related to the width w of the pressure bump by the relationship $${w}_{d}=w/\sqrt{1+\tau /\alpha}$$(5)
(Dullemond et al. 2018), which reduces to the oftenquoted relationship ${w}_{d}=w\sqrt{\alpha /\tau}$ for α ≪ τ. Dullemond et al. (2018) showed that there is typically an order of magnitude contrast between the center of the ring and its border (the “gap” between rings), arguing that the description given by Eq. (4) is valid for r − r_{0}≲ 2w_{d}.
The dynamics of the growing planet responds to the surface density distribution of the gas Σ_{g}, which sets its migration speed and direction. The latter is related to the pressure profile by the relationship $$p(r)={c}_{s}^{2}{\Sigma}_{g}/(\sqrt{2\pi}H),$$(6)
where c_{s} = HΩ is the sound speed. Assuming, for simplicity, H ∝ r (flat disk), one gets Σ_{g}(r) ∝ r^{2}p(r), with p(r) given by Eq. (1). It is trivial to see that Σ_{g}(r) is proportional to r^{2} at the location of the pressure maximum r_{0}, whereas it has a maximum at ${r}_{M}={r}_{0}(1+\sqrt{1+8{w}^{2}/{r}_{0}^{2}})/2$. For w∕r_{0} ≪ 1, the location of the maximum of Σ_{g} can be approximated by r_{M} = r_{0} + 2w^{2}∕r_{0}. Approximating Σ_{g} (r) near its maximum by a parabolic profile, we then find that that Σ(r) ∝ 1∕r at r_{1} = r_{0} + 3w^{2}∕r_{0}.
The planet’s radial migration in disks with arbitrary surface density profiles has been studied in detail in Paardekooper et al. (2010, 2011). They demonstrated that the planet feels a torque Γ from the disk that is: $$\begin{array}{lll}\Gamma \hfill & =\hfill & \frac{{\Gamma}_{0}}{\gamma}[(2.51.7b+0.1a)+1.1{f}_{\text{v}}(\nu ,M)\left(\frac{3}{2}a\right)\hfill \\ \hfill & \hfill & +{f}_{\text{E}}(\nu ,\chi ,M)\left(7.9\frac{\xi}{\gamma}\right)],\hfill \end{array}$$(7)
where ${\Gamma}_{0}={[(M/{M}_{*})({r}_{p}/{H}_{p})]}^{2}{\Sigma}_{g}({r}_{p}){r}_{p}^{2}{\Omega}_{p}^{2}$ is the nominal torque (M and M_{*} being the masses of the planet and the star respectively and r_{p} the location of the planet, where the disk has height H_{p} and rotationfrequency Ω_{p}). In Formula (7), the coefficient a is the exponent of the gas surfacedensity profile, approximated as a powerlaw of type ${\Sigma}_{g}(r)={\Sigma}_{g}({r}_{p})/{(r{r}_{p})}^{a}$ around r_{p} ; similarly, b is the exponent of the disk’s temperature profile, written as $T(r)={T}_{0}({r}_{p})/{(r{r}_{p})}^{b}$. Coherently with our assumption that the ratio H∕r in the disk is constant, we assume that b = 1. The coefficient γ in Eq. (7) is the adiabatic index, in general assumed to be equal to 1.4 and ξ = b − (γ − 1)a. In Formula (7), the first term within (.) is due to the socalled Lindblad torque; the second one is due to the vortensitydriven corotation torque, and the last is due to the entropydriven corotation torque.
The functions f_{V}(ν, M) and f_{E} (ν, χ, M) are between 0 and 1 and depend on the disk’s viscosity ν, thermal conductivity χ, and the planet’s mass (Paardekooper et al. 2011). In particular, f_{V}(ν, M) → 0 as ν → 0; in other words, the vortensitydriven corotation torque vanishes in the inviscid limit, which is called torque saturation. In this case, the planet is capable of changing the surface density profile of the disk in its corotation region toward the slope with a = 3∕2, which erases the corotation torque. However, this is only true under the hypothesis that the viscous torque is responsible for the original density profile in the disk. However, in a lowviscosity disk, there must be other processes that are responsible for shaping the pressure maximum, for instance the magnetic stress (Bethune et al. 2017), or the torque from an already existing planet opening a gap between two adjacent rings. It is reasonable to expect that these processes would restore the original disk’s profile in the planet’s corotation region, fighting against the planet’s action. Thus, we may assume that f_{V} (ν, M) ~ 1 in our case, although a specific evaluation would be needed once the actual origin of the pressure maximum is understood.
The function f_{E}(ν, χ, M) → 0 if χ → 0, for example, in the adiabatic limit, where thermal diffusivity is null. Instead, in the isothermal limit (χ →∞), the function f_{E} tends to a limit value ${f}_{\text{E}}^{\infty}(\nu ,M)<1$, and ${f}_{\text{E}}^{\infty}(\nu ,M)\to 0$ if ν → 0. Thus, a variety of values of f_{E} can be envisioned.
There is a planet trap – where the planet stops migrating (Masset et al. 2006a) – where the term within [.] in Formula (7) vanishes. Assuming the limit case f_{V} = f_{E} = 1, this happensat the location where a ~ 1 (precisely: 0.95), meaning ~3w^{2}∕r_{0} beyond the pressure bump. If instead f_{V} = 1, but f_{E} = 0, the solution of the [.] = 0 equation occurs for a = −2.55, meaning ~ 0.5w^{2}∕r_{0} inward of the pressure bump (which corresponds to a = −2). We note that, if the vortensitydriven torque were also strongly saturated (f_{V} ~ 0), the equation [.] = 0 would not have a solution: this means that there would be no planet trap and the planet would migrate out of the ring. Because this study only makes sense if the planet is confined in the ring, we assume that the vortensitydriven corotation torque is strongly desaturated (f_{V} ~ 1, as justified above), so that the planet trap is at a distance δw^{2}∕r_{0} from the pressure bump, with δ between − 0.5 and 3 depending on the actual degree of saturation of the entropydriven corotation torque. Thus, δ is a measure of where the planet trap is relative to the pressure maximum, in units of w^{2} ∕r_{0}. The radial profiles of pressure, gas, and dust surface densities, as well as the range of possible locations of the planet trap and corresponding values of δ are illustrated in Fig. 1. For completeness, had we assumed a flared disk, meaning with H∕r ∝ r^{2∕7} and b = 3∕7, the planet trap would be located at δ = 0.15–2 (again, depending on the saturation of the entropydriven torque), and the maximum of the gas surface density distribution would be at δ = 12∕7.
The planet accretes pebbles/dust at a rate $$\dot{M}=C{\Sigma}_{d}({r}_{p}),$$(8)
where Σ_{d}(r_{p}) is the density of dust at the planet’s location, and the coefficient C depends on the accretion crosssection of the planet and the velocity of the dust relative to the planet, according to the various accretion regimes (Bondi accretion, Hill accretion, 2D, 3D, etc.; see the compendium in Ida et al. 2016).
From Formula (4), the density of the dust at the planet trap is $${\Sigma}_{d}({r}_{p})={\Sigma}_{d,0}\mathrm{exp}\left(\frac{{\delta}^{2}{w}^{4}}{2{r}_{0}^{2}{w}_{d}^{2}}\right)={\Sigma}_{d,0}\mathrm{exp}\left(\frac{{\delta}^{2}{w}^{2}(\alpha +\tau )}{2\alpha {r}_{0}^{2}}\right),$$(9)
where we expressed r_{p} − r_{0} in terms of δ and used Formula (5) to pass from the second to the third expression. We note that, in a lowviscosity disk with $\alpha \ll {\delta}^{2}{w}^{2}\tau /(2{r}_{0}^{2})$, the dust density (Eq. (9)) is exponentially small at the location δ of the planet, so that the planet can accrete only at a slow rate. This is a first big difference with respect to a disk without pressure traps, where the dust density scales smoothly with r (approximately $~1/\sqrt{r}$, assuming τ and the dust massflux to be independent of r), so that, wherever the planet is, it can accrete at a comparable rate.
If, instead, δ is small, Σ_{d}(δ) ~ Σ_{d,0}. But in this case the velocity of the gas relative to the planet becomes small, because it is given by the expression: $$\Delta v=\frac{1}{2}\frac{{H}^{2}\Omega}{r}\frac{\text{d}\mathrm{log}P}{\text{d}\mathrm{log}r}=\frac{1}{2}\frac{\delta {H}^{2}\Omega}{{r}_{0}}.$$(10)
However, the planet’s accretion rate in the 3D Bondi regime, for a given particle’s τ, is independent of Δv, and it is equal to $3/2\pi {R}_{\text{H}}^{3}\Omega {\rho}_{d}$, where R_{H} = r_{p}[(M∕(3M_{*})^{]1∕3} is the Hill radius, and ρ_{d} is the volume density of dust near the midplane (Ida et al. 2016; Ormel 2017). Thus, the accretion rate is the same as in a disk without pressure bump and dust surface density Σ_{d,0}.
A caveat, however, needs to be mentioned. At δ ~ 0, the gas near the planet’s orbit is not pressure supported, and thus it has to follow horseshoe streamlines relative to the planet, even if the planet’s mass is small. The halfwidth of the horseshoe region is ${x}_{s}=1.16{r}_{p}{(M/{M}_{*})}^{1/2}{(H/r)}^{1/2}$ when the planet mass is small and tends to ~2.5R_{H} when M∕M_{*}≳ 10^{−4} (Masset et al. 2006b). If one neglects the perturbations of the planet on the gas, the nominal accretion radius in the Hill regime is ${r}_{\text{acc}}={(\tau /0.1)}^{1/3}{R}_{\text{H}}$. If x_{s} ≪ r_{acc}, the dust reaches the Hill radius of the planet without being substantially affected by the distortions of the gas streamlines, so the nominal accretion radius is the correct radius. But if x_{s} ≳ r_{acc,} the specific flow of the gas becomes important. Only the dust carried along a narrow stream of gas around the separatrix of the coorbital motion is accreted by the planet. The radial width W(M, τ) of this stream is a complicated function of planet mass and the particles’ Stokes number, and is described in Kuwahara & Kurokawa (2020). On the other hand, the shear velocity of this dust relative to the planet is 3∕2x_{s}Ω, which is larger than 3∕2r_{acc}Ω. Using numerical simulations, Kuwahara and Kurokawa showed that, if the dust density is radially uniform, these two effects approximately cancel each other out, and the actual accretion rate onto the planet is almost the same as that computed, a priori incorrectly, using the nominal accretion radius. However, in the case of a ring, the density of dust in Formula (4) is not uniform. Thus, the estimate of the dust accretion rate should be decreased by a factor $\mathrm{exp}({x}_{s}^{2}/2{w}_{d}^{2})$. In most cases, this factor will be ~1, but it can be significantly smaller if the ring is narrow (small w_{d}) and the planet massive (large x_{s}).
In any case, the planet cannot accrete more dust than the amount that the ring can deliver by diffusion and/or advection. Because the dust density profile (Eq. (4)) is set by the equilibrium between diffusion and advection, we can use either of them to evaluate the dust delivery rate to the planet’s orbit. Using the advection term (2) and expression (9), we find that the planet’s accretion rate is capped by (assuming α ≪ τ ≪ 1): $${\dot{M}}_{\text{Max}}=2\times 2\pi {r}_{p}{\Sigma}_{d}({r}_{p}){v}_{r}({r}_{p})=A\tau \delta {\Sigma}_{d,0}\mathrm{exp}\left(\frac{{\delta}^{2}{w}^{2}\tau}{2\alpha {r}_{0}^{2}}\right),$$(11)
with $A=4\pi {r}_{p}{H}_{p}^{2}{\Omega}_{p}/{r}_{0}$. We note the first factor of 2 at the right hand side of the first equals sign, which is due to the fact that dust is delivered from both sides of the planet’s orbit due to diffusion, unlike in a disk dominated by advection, where the material flows only from one side. For the reasons explained above, if δ ~ 0, the dust flux should be computed by setting $\delta =\text{max}({r}_{\text{acc}},{x}_{s}){r}_{0}/{w}_{d}^{2}$. For comparison, the maximal accretion rate in a disk without pressure bump and dust density Σ_{d,0} would be ${\dot{M}}_{\text{Max}}=4\pi {r}_{p}^{2}{\Omega}_{p}\tau \eta {\Sigma}_{d,0}$, which is a factor $~1/\delta \mathrm{exp}[{\delta}^{2}{w}^{2}\tau /(2\alpha {r}_{0}^{2})]$ larger than the one in Eq. (11). Formula (11) has a maximum for $\delta =\sqrt{\alpha /\tau}{r}_{0}/w$, where ${\dot{M}}_{\text{Max}}=A\sqrt{\alpha \tau}{r}_{0}/w{\Sigma}_{d,0}\mathrm{exp}(1/2)$, meaning it is proportional to the square root of the diffusion coefficient in the disk.
Another aspect to consider about pebble accretion in a ring is that the reservoir of solid mass is finite and equal to the dust mass of the ring. As the planet is growing by accreting dust, the total dust mass in the ring decreases as ${M}_{\text{ring}}={M}_{\text{ring}}^{\text{init}}M$, where ${M}_{\text{ring}}^{\text{init}}$ is the initial dust mass in the ring, and M is the mass of the planet. To obtain an upper bound on the planetary mass as a function with time, we assume that the planet always accretes at the maximal possible rate (Eq. (11)). By approximating $${\Sigma}_{d,0}=\frac{{M}_{\text{ring}}}{{(2\pi )}^{3/2}{r}_{0}{w}_{d}}=\frac{{M}_{\text{ring}}}{{(2\pi )}^{3/2}{r}_{0}w}\sqrt{\frac{\tau}{\alpha}}\text{\hspace{0.17em}},$$(12)
we find that the planet’s mass would evolve according to the differential equation: $$\dot{M}=\frac{A}{{(2\pi )}^{3/2}{r}_{0}w}\sqrt{\frac{{\tau}^{3}}{\alpha}}({M}_{\text{ring}}^{\text{init}}M)\delta \mathrm{exp}\left(\frac{{\delta}^{2}{w}^{2}\tau}{2\alpha {r}_{0}^{2}}\right)$$(13)
whose solution is: $M(t)={M}_{\text{ring}}^{\text{init}}[1\mathrm{exp}(Bt)]$, with $$B=\frac{A}{{(2\pi )}^{3/2}{r}_{0}w}\sqrt{\frac{{\tau}^{3}}{\alpha}}\delta \mathrm{exp}\left(\frac{{\delta}^{2}{w}^{2}\tau}{2\alpha {r}_{0}^{2}}\right).$$(14)
Therefore, initially the planet grows at a slow rate if $\alpha {r}_{0}^{2}/\tau {w}^{2}\ll 1$ or τ^{3}∕αδ^{2} ≪ 1 (we remind the reader that it’s an upper bound); then the rate slows down as the mass of the planet approaches the limit ${M}_{\text{ring}}^{\text{init}}$ asymptotically.
Fig. 1 Radial profiles of pressure (black), gas surface density (red), and dust surface density (blue) for the DHSARP ring B77, from Dullemond et al. (2018) – see Sect. 3. The vertical scale of each profile has been renormalized for illustrative purposes. The dashed thin red curves show the profiles r^{2.55} and 1∕r, and their tangential points to the red solid curve shows the locations where the gas surface density profile haslocally these profiles. These two locations delimit the range of possible positions of the planet at the migration trap, depending on the degree of saturation of the vortensitydriven corotation torque. Their distances relative to the location of the pressure maximum are given in terms of the parameter δ, as well as the location of the maximum of the gas surfacedensity distribution. The widths w and w_{d} in Eqs. (1) and (4) are also graphically represented. 
3 An application: growth of a planet in a DSHARP ring
To make a quantitative calculation of planet accretion in a dust ring, we considered the rings studied in Dullemond et al. (2018). In particular, we chose the ring denoted B77 in that paper, as it is the “median ring” in terms of distance from the star and width. The properties of this ring given in Dullemond et al. (2018) are: r_{0} = 75 AU, T = 22 K, w_{d} = 4.5 AU, M_{dust} = 40 M_{⊕}. These properties are deduced directly from the observations. The maximal width of the pressure bump w is estimated from the distance between adjacent rings and is w = 17 AU. The radio w_{d}∕w sets α∕τ= 7.7 × 10^{−2}. Educated guesses on the gas density (Σ_{g} = 10 g cm^{−2}) and dust size then lead to τ = 3 × 10^{−3} and α = 2.3 × 10^{−4}. Even if there are large uncertainties on some of these parameters, we adopted them as fiducial numbers for the sake of example. We also assumed, for simplicity, that the central star has the mass of the Sun, even though the real star is 0.8 M_{⊙}. These parameters were also used to draw Fig. 1, for illustrative purposes.
We considered a planetary seed with M = 0.1 M_{⊕}. We first assumed, for the sake of example, that the planet is fixed at three different locations (δ = 3, −0.5 or 0); then, we followed the expected location of the planet as its mass grows, using the Paardekooper et al. (2011) formulæ for torque saturation.
For the computation of the coefficient C in Eq. (8), we computed the fraction of the local dust density accreted according to the Bondi, Hill, 2D, or 3D accretion regimes, following the formulæ in Ida et al. (2016).
If the planet is at δ = 3, we find that the planet starts (and remains) in the 3D Bondi accretion regime, so that $$C=\frac{\sqrt{8\pi}}{{H}_{d}}{R}_{\text{B}}^{2}\frac{{t}_{s}}{{t}_{B}}\Delta v,$$(15)
where Δv is given in (10), the Bondi radius is ${R}_{\text{B}}=GM/{(\Delta v)}^{2}$, the crossing time is t_{B} = R_{B}∕Δv, and the stopping time is t_{s} = τ∕Ω. The mass evolution of the planet is shown over the first 3 My in Fig. 2 (red curve). We kept the disk properties fixed over this amount of time, for simplicity, but it would not be realistic to continue beyond this timescale, because the disk’s gas and dust density should decay considerably. Also, we have no knowledge of the actual lifetime of dust rings. The planet’s mass seems constant, because it increases by only 12% on the considered timescale, which is invisible at the scale of the plot. This is because, at δ = 3, the dust density is very reduced, as can be seen by extrapolating the blue curve below the lower boundary of Fig. 1.
If the planet is instead at δ = −0.5, the dust density is higher, and the growth is significantly faster. After 1.6 My, the planet has reached 0.27M_{⊕} and the accretion switches to the 3D Hill regime, so that $$C=\frac{3\sqrt{\pi}}{\sqrt{8}{H}_{d}}{R}_{\text{H}}^{3}\frac{\tau}{0.1}\Omega .$$(16)
The planet reaches a final mass of 0.8 M_{⊕} (Fig. 2, blue curve).
If the planet is at the dust trap, it can grow at the maximal speed. In this case, the planet reaches 1.1 M_{⊕} in the end (Fig. 2, cyan curve), with a growth that is basically identical to that of a planet in a disk with dust density Σ_{d,0} without pressure bumps, because ${x}_{s}\ll 2\sqrt{2}{w}_{d}$ throughout the planet’s growth. In all three cases we considered, the growth rate is much smaller than the maximum rate reported in Eq. (13).
To add a bit of realism, we now compute, for each planet mass, the expected location of the planet’s trap accounting for the saturation of the corotation torques, and compute the accretion rate of the planet at that location. To evaluate the saturation of the corotation torques, we used formulæ (52) and (53) from Paardekooper et al. (2011) using the disk parameters reported above, and an opacity κ = 0.03 cm^{2} g^{−1}. The result is illustrated by the black curve, which is very close to the cyan curve for the planet at the pressure maximum. In this case, the planet indeed starts at δ = 0.2, meaning, very close to the pressure maximum, and slowly moves away from it, reaching δ = 0.4 at the end of the simulation, as labeled on the plot. Throughout the simulation, the saturation functions f_{V} and f_{E} are very close to the values corresponding to the linear components of the torques, namely 0.7/1.1 and 1.7/7.9, respectively (Paardekooper et al. 2011). In fact, for the considered planetary masses, the ratios $\nu {r}_{p}/({x}_{s}^{3}\Omega )$ and $\chi {r}_{p}/({x}_{s}^{3}\Omega )$, which enter in the saturation functions, are both much larger than unity, meaning that the system is in the isothermal, highviscosity regime (Paardekooper et al. 2011).
We conclude that planetary seeds are unlikely to become a sizable (observable) planet in this ring (and probably all rings) observed by the DSHARP survey. This is consistent with the fact that these rings are observable. In fact, if the rings had formed a planet containing most of their mass within the age of the disk, they would no longer be visible.
It is therefore interesting to rescale the considered ring into the inner disk. For simplicity, we multiplied r_{0}, w and w_{d} by 0.067, so that r_{0} = 5 AU. We kept the same values of α and τ. We multiplied the mass of the ring by $\sqrt{0.067}$, which is equivalent to assuming a disk’s global surfacedensity profile ∝ 1∕r^{3∕2}. We also reset the temperature to T = 120 K and the aspect ratio to H∕r = 0.05, which are appropriate for a disk at 5 AU.
The mass growths for planets situated at δ = 3, −0.5, and 0 in such a rescaled ring are illustrated in Fig. 3. Both the planets located at δ = 0 and − 0.5 grow at almost indistinguishable rates and show the saturation effect discussed at the end of Sect. 2. In fact, they accrete all the available mass (10.2 M_{⊕}) in less than 1/2 My. The planet at δ = 3 grows muchmore slowly, although faster than in the ring at 75 AU, and in this case reaches 3.5 M_{⊕} at the end of the simulation. The planet that is allowed to migrate as it grows (black curve) shows very interesting behavior. Initially, it is at δ= 0.6, but the planet trap moves farther away from the pressure maximum as the planet grows, reaching δ = 2.7 at M = 3M_{⊕}. This limits the planet’s accretion rate relatively to the cases with δ = 0.5 or 0 shown before, because the density of dust at large δ is lower. Figure 3 labels the value of δ as a function of the planet’s mass. For larger masses, the planet migrates back toward the pressure maximum. This behavior is the effect of the corotation torques; f_{V} and f_{E} are initially close to the values corresponding to the linear components of the torques, then they evolve towards ~ 1 as the planet grows and furthermore they decrease as the ratios $\nu {r}_{p}/({x}_{s}^{3}\Omega )$ and $\chi {r}_{p}/({x}_{s}^{3}\Omega )$ approach unity. As a result of this back and forth migration relative to the center of the ring, the planet reaches the maximal mass of 10.2 M_{⊕} at 2.5 My. However, we warn that the saturation functions f_{V} and f_{E} are very sensitive on the disk’s parameters, and different assumptions may lead to different results and evolutions. Here, we assumed an opacity κ = 3 cm^{2} g^{−1}. A lower value of κΣ_{g} would have kept the planet closer to the pressure maximum, enhancing its accretion rate. Also, we remind the reader that it is unclear how saturation would occur in the presence of additional forces (e.g., magnetic stresses) shaping the disk’s density profile.
Of course, we do not know whether rings of this kind exist in the inner part of protoplanetary disks, but if they do, they may be effective sites to produce giant planet cores, provided that the planet remains in the proximity of the pressure maximum.
Fig. 2 Evolution of the mass of a planet of initially 0.1 M_{⊕} accreting pebbles in the B77 DHARP ring, for three different locations (i.e., values of δ) relative to the position of the pressure maximum as well as a case (black curve) where the planet follows the (massdependent) location of the planet trap as it grows, and therefore δ changes withtime (from ~0.2 to ~ 0.4). 
Fig. 3 Same as Fig. 2, but rescaling the geometry of the B77 DHSARP ring at 5 AU and increasing the dust surface density ∝ 1∕r^{3∕2}. 
4 Conclusions
PierreSimon de Laplace is recognized to have been the first scientist to theorize back in 1796 the formation of a circumstellar disk due to angular momentum conservation during the contraction of gas towards the central star (Emmanuel Kant had also envisioned the formation of a disk, but his reasoning was not scientifically correct). But what is less widely known is that Laplace, to explain the final formation of a discrete set of planets, envisioned that the circumstellar disks would break in rings, each ring giving origin to one planet. Astonishingly, modern resolved observationsof protoplanetary disks (Andrews et al. 2018) indeed show that most of them are structured in rings.
In this paper, we have examined how the now popular paradigm of planet formation by pebble accretion changes if the rings are due to pressure maxima that prevent the global drift of dust in the disk. We discuss in the Introduction that the radial drift of dust through the entire disk would lead to consequences (in terms of dust distribution, dust/gas ratio, formation of planets at the inner disk edge) that are inconsistent with the observations. This suggests that pressure maxima are ubiquitous in protoplanetary disks (Pinilla et al. 2012), probably also in the inner parts not yet observationally resolved by current instruments.
In the analytic part of this paper (Sect. 2), we point out that planet Type I migration is probably stalled in a ring, but the location at which the planet is at rest is not necessarily the pressure maximum. It depends on the degree of saturation of the socalled corotation torques, particularly the one depending on the entropy gradient (Paardekooper et al. 2010, 2011). Thus, the dust accumulation may be offset from the planet’s radial position by a significant amount. In this case, the planet finds itself in a lowdensity region of dust, and its accretion is necessarily slow.
If the planet location is near the pressure maximum instead, the dust density is high, but pebble accretion still proceeds differently than in the classic case where dust has a global radial drift. Firstly, the difference in orbital velocity between the planet and the gas/particles is much smaller (it is null at the exact pressure maximum). Thus, the Bondi radius becomes very large, and the accretion regime already transits to the Hill regime at small planet masses (Ormel 2017). Moreover, in the absence of differential azimuthal velocity imposed by the gas pressure gradient, the motion of dust relative to the planet is governed by threebody coorbital dynamics instead of being a simple flyby+drag process (Kuwahara & Kurokawa 2020). We show, however, that both these aspects have moderate effects on the accretion rate of the planet, which remains similar to the one that the same planet would have in a disk without pressure maximum but with the same dust surface density.
More relevant is the fact that the planet cannot accrete more mass than that advected towards the planet’s orbit in the ring. The density profile in the ring is the result of the balance between the drift of dust towards the pressure maximum and the turbulent diffusion in the disk (Dullemond et al. 2018). Far from the pressure maximum, the drift is fast, but the equilibrium surface density of dust is exponentially low; approaching the pressure maximum, it’s the drift speed that tends to be zero. Thus, the maximal accretion rate is, in general, significantly lower than the maximal accretion rate in a classic, powerlaw disk and it can, in some cases, effectively taper the mass growth of the planet. Moreover, the reservoir of mass available to the planet is finite (i.e., the total mass of the ring) so that the maximal planet’s mass is bounded. In Sect. 2 we obtain an analytic formula of the mass evolution of a planet assuming that it accretes at the maximal rate until exhausting the ring’s dust reservoir. Of course, this formula is derived from very optimistic assumptions and serves just as an upper bound of the real mass evolution.
In Sect. 3, we compute the accretion rate of a planet in a typical ring observed by the DSHARP survey (Dullemond et al. 2018), using a numerical implementation of the analytic formulæ. We show that if the planet feels unsaturated corotation torques, it basically does not grow because it lies too far from the peak of the ring’s dust density. If the planet lies closer to the ring’s center, its accretion is more significant, but nevertheless insufficient for the planet to grow beyond one Earth masses in 3 My. However, it is fair to say that this inefficiency of planet formation is not due to the dust being confined in a ring, but to the long orbital timescales and low dust densities, even at peak values.
This result suggests that it is unlikely that the DSHARP rings themselves are the result of the formation of a first generation of planets. In fact, in a disk with uniform density Σ_{d,0}, the formation of the first generation of planets would still have taken as long as for the planet illustrated by the magenta curve in Fig. 2, reaching at most ~1 M_{⊕}. Also, a planet with this mass is not enough to form a pressure bump in a disk with 10% aspect ratio (Crida et al. 2005). Only planets growing when the disk had a higher dust density could have become more massive (Manara et al. 2018). However, it is important to remember that without preexisting pressure bumps, the dust would have drifted towards the star on a timescale of 5000 orbital periods (assuming a differential velocity of the gas with respect to the Keplerian velocity of 0.1c_{s} and τ = 3 × 10^{−3}). Thus, the growth of the first generation of planets would have rapidly stalled due to the absence of solid material. These conclusions are in apparent conflict with Pinte et al. (2020), who claim the detection of velocity “kinks” in eight of 18 circumstellar disks observed by the DSHARP program, which could be due to the presence of planets. However, we note that, if confirmed, these kinks would require the existence of planets of multiple Jupiter masses; massive planets at such large distances could have been formed by gravitational instability (see e.g., Kratter & Lodato 2016 for a recent review) rather than coreaccretion.
We then “rescaled” the considered ring to 5 AU, reducing its widths and vertical thickness proportionally to its radial location and scaling the dust surface density ∝ 1∕r^{3∕2}. In this case, if the planet is close to the ring’s center, it can accrete the whole mass of the ring within half My and can become a giant planet core. A planet that is free to position itself within the ring according to the massdependent partial saturation of its corotation torques can also accrete the full mass of the ring, although on a longer timescale.
This experiment reveals another difference with respect to the pebbleaccretion paradigm in the absence of pressure bumps. In that case, the growth of a planet continues until it reaches the socalled pebbleisolation mass (Lambrechts et al. 2014; Bitsch et al. 2018), at which point it creates a pressure bump that stops the flux of pebbles. When this happens, there is no additional energy deposited by the accretion of solids, the core can cool off and it can start accreting gas (Lambrechts & Johansen 2014). This predicts that the cores of giant planets should always be of the order of the pebble isolation mass. For a planet in a disk with H∕r = 0.05 (a typical value at 5 AU for a viscous accretion of ~ 10^{−7}M_{⊙} y^{−1}; Bitsch et al. 2014), this should be ~20 M_{⊕}. But the core of Jupiter could be smaller than this estimate by a factor of ~ 3 (Wahl et al. 2017); there could even be no compact core (Debras & Chabrier 2019). A pebbleisolation mass of ~ 5–7 M_{⊕} would require H∕r ≲ 0.03 which, at 5 AU, implies a very late disk with no accretional heating; but in this case, it would be hard to accrete the massive envelope around the core. A suggested solution of this conundrum is that part of the core has been eroded into the envelope (Stevenson 1985; Guillot et al. 2004). However, formation in a ring could be another elegant solution. In fact, accretion stops when the mass of the ring is exhausted, which can happen well before the pebble isolation mass is reached. From that point, accretion of gas could start. This would produce giant planets with small cores.
Acknowledgements
I wish to thank K. Batygin, W. Bethune, B. Bitsch, A. Crida, A. Johanen, M. Lambrechts and E. Lega for their comments that have improved this manuscript. This work is part of the project Gepard, funded by the French ANR and German DFG, focused on migration in lowviscosity disks.
References
 ALMA Partnership (Brogan, C. L., et al.) 2015, ApJ, 808, L3 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, D. E., Blake, G. A., Bergin, E. A., et al. 2019, ApJ, 881, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, ApJ, 869, L41 [NASA ADS] [CrossRef] [Google Scholar]
 Ansdell, M., Williams, J. P., Trapman, L., et al. 2018, ApJ, 859, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Béthune, W., Lesur, G., & Ferreira, J. 2017, A&A, 600, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Birnstiel, T., Dullemond, C. P., & Brauer, F. 2010, A&A, 513, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Birnstiel, T., Fang, M., & Johansen, A. 2016, Space Sci. Rev., 205, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Bitsch, B., Morbidelli, A., Lega, E., et al. 2014, A&A, 564, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bitsch, B., Morbidelli, A., Johansen, A., et al. 2018, A&A, 612, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bitsch, B., Izidoro, A., Johansen, A., et al. 2019, A&A, 623, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brasser, R., & Mojzsis, S. J. 2020, Nat. Astron., 4, 492 [NASA ADS] [CrossRef] [Google Scholar]
 Brauer, F., Dullemond, C. P., & Henning, T. 2008, A&A, 480, 859 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587 [Google Scholar]
 Delbo’, M., Walsh, K., Bolin, B., et al. 2017, Science, 357, 1026 [NASA ADS] [CrossRef] [Google Scholar]
 Debras, F., & Chabrier, G. 2019, ApJ, 872, 100 [NASA ADS] [CrossRef] [Google Scholar]
 Dullemond, C. P., Birnstiel, T., Huang, J., et al. 2018, ApJ, 869, L46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Flock, M., Fromang, S., Turner, N. J., et al. 2016, ApJ, 827, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Flock, M., Fromang, S., Turner, N. J., et al. 2017, ApJ, 835, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Flock, M., Turner, N. J., Mulders, G. D., et al. 2019, A&A, 630, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fortier, A., Alibert, Y., Carron, F., et al. 2013, A&A, 549, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guillot, T., Stevenson, D. J., Hubbard, W. B., et al. 2004, Jupiter. The Planet, Satellites and Magnetosphere (Cambridge: Cambridge University Press), 35 [Google Scholar]
 Ida, S., Guillot, T., & Morbidelli, A. 2016, A&A, 591, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Johansen, A., & Bitsch, B. 2019, A&A, 631, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Johansen, A., & Lambrechts, M. 2017, Ann. Rev. Earth Planet. Sci., 45, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Johansen, A., Mac Low, M.M., Lacerda, P., et al. 2015, Sci. Adv., 1, 1500109 [NASA ADS] [CrossRef] [Google Scholar]
 Johansen, A., Ida, S., & Brasser, R. 2019, A&A, 622, A202 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kratter, K., & Lodato, G. 2016, ARA&A, 54, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Kretke, K. A., Levison, H. F., Buie, M. W., et al. 2012, AJ, 143, 91 [Google Scholar]
 Kruijer, T. S., Burkhardt, C., Budde, G., et al. 2017, Proc. Natl. Acad. Sci., 114, 6712 [NASA ADS] [Google Scholar]
 Kuwahara, A., & Kurokawa, H. 2020, A&A, 633, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lambrechts, M., & Johansen, A. 2012, A&A, 544, A32 [Google Scholar]
 Lambrechts, M., & Johansen, A. 2014, A&A, 572, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lambrechts, M., Johansen, A., & Morbidelli, A. 2014, A&A, 572, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lambrechts, M., Morbidelli, A., Jacobson, S. A., et al. 2019, A&A, 627, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Levison, H. F., Thommes, E., & Duncan, M. J. 2010, AJ, 139, 1297 [NASA ADS] [CrossRef] [Google Scholar]
 Levison, H. F., Kretke, K. A., & Duncan, M. J. 2015, Nature, 524, 322 [CrossRef] [Google Scholar]
 Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Manara, C. F., Rosotti, G., Testi, L., et al. 2016, A&A, 591, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Manara, C. F., Morbidelli, A., & Guillot, T. 2018, A&A, 618, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Masset, F. S., Morbidelli, A., Crida, A., et al. 2006a, ApJ, 642, 478 [NASA ADS] [CrossRef] [Google Scholar]
 Masset, F. S., D’Angelo, G., & Kley, W. 2006b, ApJ, 652, 730 [NASA ADS] [CrossRef] [Google Scholar]
 Morbidelli, A., & Nesvorný, D. 2020, The Transneptunian Solar System (Amsterdam: Elsevier), 25 [Google Scholar]
 Morbidelli, A., Bottke, W. F., Nesvorný, D., et al. 2009, Icarus, 204, 558 [NASA ADS] [CrossRef] [Google Scholar]
 Morbidelli, A., Lambrechts, M., Jacobson, S., et al. 2015, Icarus, 258, 418 [NASA ADS] [CrossRef] [Google Scholar]
 Najita, J. R., & Bergin, E. A. 2018, ApJ, 864, 168 [Google Scholar]
 Ormel, C. W. 2017, Astrophys. Space Sci. Lib., 445, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Ormel, C. W., & Klahr, H. H. 2010, A&A, 520, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Paardekooper,S.J., Baruteau, C., Crida, A., et al. 2010, MNRAS, 401, 1950 [NASA ADS] [CrossRef] [Google Scholar]
 Paardekooper,S.J., Baruteau, C., & Kley, W. 2011, MNRAS, 410, 293 [NASA ADS] [CrossRef] [Google Scholar]
 Pinilla, P., Birnstiel, T., Ricci, L., et al. 2012, A&A, 538, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pinte, C., Dent, W. R. F., Ménard, F., et al. 2016, ApJ, 816, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Pinte, C., Price, D. J., Ménard, F., et al. 2020, ApJ, 890, L9 [NASA ADS] [CrossRef] [Google Scholar]
 Riols, A., & Lesur, G. 2019, A&A, 625, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rosotti, G. P., Tazzari, M., Booth, R. A., et al. 2019, MNRAS, 486, 4829 [NASA ADS] [CrossRef] [Google Scholar]
 Sato, T., Okuzumi, S., & Ida, S. 2016, A&A, 589, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 500, 33 [NASA ADS] [Google Scholar]
 Simon, J. B., Armitage, P. J., Youdin, A. N., et al. 2017, ApJ, 847, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Stevenson, D. J. 1985, Icarus, 62, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Tominaga, R. T., Takahashi, S. Z., & Inutsuka, S..i. 2019, ApJ, 881, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Trapman, L., Facchini, S., Hogerheijde, M. R., et al. 2019, A&A, 629, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 WafflardFernandez, G., & Baruteau, C. 2020, MNRAS, 493, 5892 [NASA ADS] [CrossRef] [Google Scholar]
 Wahl, S. M., Hubbard, W. B., Militzer, B., et al. 2017, Geophys. Res. Lett., 44, 4649 [NASA ADS] [CrossRef] [Google Scholar]
 Weber, P., Pérez, S., BenítezLlambay, P., et al. 2019, ApJ, 884, 178 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, C.C., & Zhu, Z. 2020, MNRAS, 491, 4702 [Google Scholar]
 Youdin, A. N., & Lithwick, Y. 2007, Icarus, 192, 588 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, S., Zhu, Z., Huang, J., et al. 2018, ApJ, 869, L47 [NASA ADS] [CrossRef] [Google Scholar]
The original formula in Dullemond et al. (2018) uses the sound speed c_{s}, which is equal to HΩ.
All Figures
Fig. 1 Radial profiles of pressure (black), gas surface density (red), and dust surface density (blue) for the DHSARP ring B77, from Dullemond et al. (2018) – see Sect. 3. The vertical scale of each profile has been renormalized for illustrative purposes. The dashed thin red curves show the profiles r^{2.55} and 1∕r, and their tangential points to the red solid curve shows the locations where the gas surface density profile haslocally these profiles. These two locations delimit the range of possible positions of the planet at the migration trap, depending on the degree of saturation of the vortensitydriven corotation torque. Their distances relative to the location of the pressure maximum are given in terms of the parameter δ, as well as the location of the maximum of the gas surfacedensity distribution. The widths w and w_{d} in Eqs. (1) and (4) are also graphically represented. 

In the text 
Fig. 2 Evolution of the mass of a planet of initially 0.1 M_{⊕} accreting pebbles in the B77 DHARP ring, for three different locations (i.e., values of δ) relative to the position of the pressure maximum as well as a case (black curve) where the planet follows the (massdependent) location of the planet trap as it grows, and therefore δ changes withtime (from ~0.2 to ~ 0.4). 

In the text 
Fig. 3 Same as Fig. 2, but rescaling the geometry of the B77 DHSARP ring at 5 AU and increasing the dust surface density ∝ 1∕r^{3∕2}. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.