## Dust-Gas Coupling in Turbulence- and MHD Wind-Driven Protoplanetary Disks: Implications for Rocky Planet Formation

Teng Ee Yap 1 and Konstantin Batygin 1

1 Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125, USA. Corresponding author: tyap@caltech.edu

## ABSTRACT

The degree of coupling between dust particles and their surrounding gas in protoplanetary disks is quantified by the dimensionless Stokes number. The Stokes number ( St ) governs particle size and spatial distributions, in turn establishing the dominant mode of planetary accretion in different disk regions. In this paper, we model the characteristic St of particles across time in disks evolving under both turbulent viscosity and magnetohydrodynamic (MHD) disk winds. In both turbulenceand wind-dominated disks, we find that collisional fragmentation is the limiting mechanism of particle growth. The water-ice sublimation line constitutes a critical transition point between dust settling, drift, and size regimes. For a fiducial value disk evolution parameter ˜ α ≃ 10 -3 , silicate particles interior to the ice-line in our model are characterized by low St ( ≲ 10 -2 ) and sizes in the sub-mm- to 1cm-scale. Icy particle/boulders beyond the ice-line are characterized by high St ( ≳ 10 -2 ) and sizes in the cm to dm size range. Hence, icy particles settle into a thin layer at the outer disk midplane and drift inward at velocities exceeding the gaseous accretionary flow due to substantial headwind drag. Silicate particles in the inner disk remain relatively well dispersed and are to a large extent advected inward with their surrounding gas.

The St dichotomy across the ice-line translates to distinct planet formation pathways between the inner and outer disk. While pebble accretion proceeds slowly for rocky embryos within the ice-line (across most of parameter space), it does so rapidly for volatile-rich embryos beyond it, allowing for the growth of giant planet cores before disk dissipation. Through simulations of rocky planet growth, we evaluate the competition between pebble accretion and classical pairwise collisions between planetesimals. We conclude that the dominance of pebble accretion can only be realized in disks that are driven by MHD winds, slow-evolving (˜ α ≲ 10 -3 5 . ), and devoid of pressure maxima that may concentrate solids and give rise of planetesimal rings in which classical growth is enhanced. Such disks are extremely quiescent, with Shakura-Sunyaev turbulence parameters α ν ≲ 10 -4 . We conclude that for most of parameter space corresponding to values of α ν reflected in observations of protoplanetary disks ( ≳ 10 -4 ), pairwise planetesimal collisions constitute the dominant pathway of rocky planet accretion. Our results are discussed in the context of super-Earth origins, and lend support to the emerging view that they formed in planetesimal rings. Moreover, these results argue against a significant contribution ( ≳ 10%) of outer disk, carbonaceous, material to the proto-Earth in the form of pebbles, in agreement with chemical and isotopic investigations of Earth's accretion history.

Keywords : Protoplanetary Disks, MHD, Pebble Accretion, Planets and Satellites: Formation.

## 1. INTRODUCTION

The evolution of the size and spatial distributions of dust particles in protoplanetary disks is intimately linked with the coupling of these particles to the H 2 -He gas in which they are embedded, and ultimately determines the large-scale architecture of planetary systems that emerge upon disk dissipation. Improvements in our understanding of the interplay between dust and gas in such disks thus translates to a more refined picture of the planet formation process, yielding insights into how our Solar System (SS) evolved to be a comparatively uncommon outcome of planet formation within the Galaxy. Indeed, the SS is devoid of the archetypal close-in rocky planets characteristic of many systems in the galactic exoplanetary census ( e.g., Petigura et al., 2013; Mulders et al., 2018; Batygin &amp; Morbidelli, 2023), yet hosts gas/ice giants on long-period orbits, found only around

∼ 15% of Sun-like stars ( e.g., Rosenthal et al., 2022).

The degree of dust-gas coupling is quantified by the dimensionless 'Stokes number,' defined as the ratio of the particle stopping time to the turnover time of the largest turbulent eddy in the disk, the latter customarily taken as the inverse of the Keplerian angular velocity (Cuzzi et al., 2001; Ormel &amp; Cuzzi, 2007; Sengupta et al., 2024). As such, the Stokes number ( St ) represents the fraction of an orbit (given by the ratio of the stopping time to the orbital period) in radians required for a dust particle to equilibrate with its surrounding gas flow. The St governs both the efficiency of planetesimal formation ( i.e., the onset and growth rate of resonant drag instabilities; Youdin &amp; Goodman, 2005; Yang et al., 2017; Squire &amp; Hopkins, 2018) and the mode of planetary accretion favored in different regions of disks ( e.g., Ormel, 2017).

Recently, Batygin &amp; Morbidelli (2022) (BM22) presented a self-consistent model for St within a steadystate, 'actively' heated viscous disk. In the context of this model, the change in dust composition from purely silicate to icy across the water-ice sublimation line (hereafter denoted simply as the 'ice-line') manifests as a large ( &gt; an order of magnitude) difference in St between the inner and outer disk. This dichotomy in St reflects the difference in fragmentation threshold velocities between silicate and ice grains (Blum &amp; M¨nch, 1993; u Ormel &amp; Cuzzi, 2007; G¨ttler et al., 2010), and underu lies their key conclusion that pebble accretion (Ormel &amp; Klahr, 2010; Lambrechts &amp; Johansen, 2012) can only operate efficiently beyond the ice-line, in the outer disk. This suggests planets in the inner disk primarily grew via stochastic, pairwise collisions between planetesimals ( e.g., Chambers &amp; Wetherill, 1998; Safronov, 1972). While this conclusion is likely robust to leading order, the standard viscous disk theory on which the model is built invites layers of complexity that can yield additional insights into where and how planets form.

Despite the historic appeal of turbulent viscosity as the driver of disk evolution, the extent to which disks are truly turbulent remains unknown. It is also unclear if mechanisms thought to source turbulence ( e.g., the magnetorotational instability (MRI), Balbus &amp; Hawley, 1991; the vertical shear instability (VSI), Nelson et al. 2013) can do so with sufficient vigor to explain observed stellar mass accretion rates ( e.g., Trapman et al., 2022). Magnetohydrodynamic (MHD) disk winds are garnering increasing attention as an alternative avenue for facilitating disk evolution. In essence, given that disks are sufficiently ionized at the surface, magnetic field lines threading the disk allow for the extraction of mass and angular momentum along disk-wind streamlines, driving accretion (Blandford &amp; Payne, 1982; Armitage, 2020).

Tabone et al. (2022) introduced a phenomenological treatment of MHD disk winds, inspired by the functional α -prescription for turbulent viscosity from Shakura &amp; Sunyaev (1973) and the self-similar solutions for disk surface density evolution from Lynden-Bell &amp; Pringle (1974). Using their model, we revisit the exercise of modeling the characteristic St of dust particles across protoplanetary disks, considering other processes beyond fragmentation that may set the characteristic St in different regions of the disk, such as radial drift (Weidenschilling, 1977) and bouncing of colliding particles (Zsom et al., 2010; Windmark et al., 2012). Accordingly, our model illustrates how St evolves both in time and in transition between turbulence-dominated and winddominated disks.

The paper is structured as follows. In Section 2, we provide an overview of the disk model adopted and its parameters. We also define model constants that yield disk surface density and temperature profiles most consistent with heating from both viscous shear and stellar irradiation, as well as the initial disk mass prescribed. In Section 3, we derive characteristic St profiles for two fiducial disks: one turbulence-dominated and the other wind-dominated. We show how these profiles translate to dust properties ( i.e., vertical distribution, radial drift rate, physical size) and planetary accretion regimes across these disks in Section 4. We also assess the favored pathway of planetary accretion within and without the ice-line across parameter space, with a focus on the former, through consideration of mass-doubling timescales and simulations of rocky planet growth. Our results are discussed in the context of super-Earth origins and the Earth's accretion. We consider the potential implications our assumptions have on St across disks in Section 5 before giving our concluding remarks in Section 6.

## 2. DISK MODEL

## 2.1. Overview

We begin by providing an overview of the formulation devised by Tabone et al. (2022) to describe the evolution of the disk surface density Σ( r, t ) under the influence of both turbulence and MHD disk winds (Fig. 1). We defer the reader to their work for an in-depth derivation of the equations to be introduced. Aside from the stellar mass accretion rate ˙ M , a proxy for time ( e.g., Hartmann et al., 1998), and the α turbulence parameter (henceforth denoted α ν ), the inclusion of disk winds necessitates the addition of two parameters to specify the evolution of Σ( r, t ). The first, denoted α DW , is analogous to α ν and quantifies the torque exerted on the disk by winds. The

## MHD Disk Winds

Figure 1. Schematic of a protoplanetary disk evolving under both turbulence and MHD disk winds. Turbulence manifests in an effective viscosity, and facilitates angular momentum transport. MHD disk winds arise from magnetic field lines that thread the disk, and facilitate angular momentum and mass loss. No solids exist interior to the silicate sublimation line, and the ice-line demarcates the boundary between the inner and outer disk, rocky and icy dust particles. The midplane temperature of the flared disk is set by a balance between 'active' viscous heating and 'passive' stellar irradiation, and radiative loss. See Sections 2.1 &amp; 2.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-099feae4-760c-4308-8ef6-87211ccc6418.jpg)

total torque exerted on the disk is encapsulated in ˜, α given by

$$\tilde { \alpha } = \alpha _ { \nu } + \alpha _ { D W }. \text{ \quad \ \ } ( 1 )$$

k b being the Boltzmann constant and µ the mean molecular mass of the disk gas ( ∼ 2.4 proton masses). Note that the Keplerian orbital velocity v K = Ω r K , where r is the radial distance from the star.

The relative contribution of turbulence and disk winds to ˜ is quantified by the ratio of α α DW to α ν , denoted ψ :

$$\psi = \alpha _ { D W } / \alpha _ { \nu }. \text{ \quad \ \ } \text{(2)}$$

Accordingly, disks characterized by high and low ψ are wind- and turbulence-dominated, respectively. The second additional parameter is the magnetic lever arm parameter λ (Blandford &amp; Payne, 1982), the angular momentum per unit mass extracted by disk-wind streamlines. The greater λ is, the smaller the extracted mass required per unit time to explain the contribution to ˙ M from disk winds ( i.e., to yield a given value of α DW ). Thus, λ ultimately sets the rate at which mass is extracted from the disk, given a specification for α DW (or ˜ and α ψ ). The mass loss rate due to disk winds Σ ˙ DW is given by

$$\dot { \Sigma } _ { D W } = \frac { 3 \alpha _ { D W } c _ { s } ^ { 2 } \Sigma } { 4 ( \lambda - 1 ) \Omega _ { K } r ^ { 2 } }, \quad \quad ( 3 ) \quad \ \ 1$$

where r is the radial/cylindrical distance from the star of mass M , Ω K = √ GM/r 3 is the Keplerian angular velocity, and c s = √ k T/µ b the isothermal sound speed, with

With α DW and λ defined, the evolution of Σ( r, t ), given by mass and angular momentum continuity, is expressed as

$$\frac { \partial \Sigma } { \partial t } & = \frac { 3 } { r } \frac { \partial } { \partial r } \left [ \frac { 1 } { r \Omega } \frac { \partial } { \partial r } ( r ^ { 2 } \alpha _ { \nu } \Sigma c _ { s } ^ { 2 } ) \right ] \\ & + \frac { 3 } { 2 r } \frac { \partial } { \partial r } \left [ \frac { \alpha _ { D W } \Sigma c _ { s } ^ { 2 } } { \Omega } \right ] - \dot { \Sigma } _ { D W }$$

Above, the first and second terms on the RHS of the equation represent viscous and wind-driven accretion, respectively. Assuming α ν and α DW are constant in time and space, as well as a temperature profile scaling with r as T ∼ r γ -3 / 2 (see Appendix C in Tabone et al., 2022), the ansatz to the above equation takes the form

$$\Sigma ( r, t ) = \Sigma _ { c } ( t ) \left ( \frac { r } { r _ { c } ( t ) } \right ) ^ { \xi - \gamma } e x p \left [ - ( \frac { r } { r _ { c } ( t ) } ) ^ { 2 - \gamma } \right ], \ \ ( 5 )$$

underpinned by the analytical solution from LyndenBell &amp; Pringle (1974) for the evolution of Σ( r, t ) in a viscous disk with constant α ν . The variable ξ is the socalled 'mass ejection index' describing the flattening of the steady-state Σ profile ( i.e. , that without the exponential decay term, such that Σ ∼ r ξ -γ ; ξ &gt; 0) due to

)

(

10

-3

(

Figure 2. Rasters of the deviation in the retrieved vs. assumed temperature power law index δ for (a) ˜ = 10 α -4 and ψ = 0 1 . , (b) ˜ = 10 α -3 and ψ = 1 0 . , and (c) ˜ = 10 α -2 and ψ = 10 ( i.e., diagonal across ˜ α -ψ parameter space). It is clear that most of the variation in δ is attributed to γ . The green point along the center of all three plots, at λ = 3 5 adopted as a constant in our . disk model, represents the optimal values of γ , which changes with ˜ and α ψ . These optimal values yield δ ≲ 0 1. . See Section 2.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-f995fade-22fd-4afb-9cd3-7ba0ff1757bd.jpg)

the extraction of mass by disk winds, and is given by

$$\xi = \frac { 1 } { 4 } ( \psi + 1 ) \left ( \sqrt { 1 + \frac { 4 \psi } { ( \lambda - 1 ) ( \psi + 1 ) ^ { 2 } } } - 1 \right ). \quad ( 6 ) \quad \begin{array} { c } \mathfrak { a } \mathfrak { g } \\ \mathfrak { d i s } \\ \mathfrak { e q } \end{array}$$

timescale t acc, 0 is defined as the time it takes to accrete a gas parcel across a distance of r c (0) / 2, keeping its radial velocity constant at its initial value v r,gas,c (0). It is expressed as

Accordingly, γ is the power law index of Σ( r, t ) in the limit of ψ ≪ 1 ( i.e., ξ ≃ 0; no winds) and small r (where the exponential term is insignificant). For γ = 1 in the absence of disk winds, Σ at small r scales roughly as r -1 and T scales as r -1 / 2 as in the passively irradiated blackbody disk of Chiang &amp; Goldreich (1997).

Additionally, Σ ( ) represents the characteristic surc t face density, and r c ( ) t the characteristic outer edge of the disk at which Σ( r , t c ) = Σ ( ) c t /e . The former is given by

$$\Sigma _ { c } ( t ) = \Sigma _ { c } ( 0 ) \left ( 1 + \frac { t } { ( 1 + \psi ) t _ { a c c, 0 } } \right ) ^ { ( 5 + 2 \xi + \psi ) / [ 2 ( 2 - \gamma ) ] }, \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{
  \text{ \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \end{cases}$$

and the latter by

$$r _ { c } ( t ) = r _ { c } ( 0 ) \left ( 1 + \frac { t } { ( 1 + \psi ) t _ { a c c, 0 } } \right ) ^ { 1 / ( 2 - \gamma ) }. \quad ( 8 ) \quad \begin{array} { c } \ddot { M } \\ \end{array} \vdots \end{array}$$

We identify t = 0 with the end of the disk's infall stage, at which point the characteristic surface density Σ (0) c can be expressed with the initial disk mass M disk (0) and radial extent of the disk r c (0) as

$$\Sigma _ { c } ( 0 ) = \beta \frac { M _ { d i s k } ( 0 ) } { 2 \pi r _ { c } ^ { 2 } ( 0 ) }, \quad \quad ( 9 )$$

where we have introduced a factor of order unity β used to ensure the integration of mass across the disk ( i.e., ∫ ∞ 0 1 . AU 2 πr Σ( r, 0) dr , where 0 1 . AU is taken as the disk inner edge) retrieves M disk (0). The initial accretion

$$t _ { a c c, 0 } = \frac { r _ { c } ( 0 ) } { 2 v _ { r, g a s, c } ( 0 ) } = \frac { r _ { c } ( 0 ) v _ { K, c } } { 3 ( 2 - \gamma ) ^ { 2 } c _ { s, c } ^ { 2 } \tilde { \alpha } }, \quad ( 1 0 )$$

where we have substituted v r,gas,c (0)= (3 / 2)(2 -γ ) 2 ˜ αc 2 s,c v -1 K,c in the second equality. The variables c s,c and v K,c represent the sound speed and Keplerian orbital velocity at r = r c (0).

To summarize, defining a Σ( r, t ) profile of a disk evolving under both turbulence and disk winds requires specification of eight parameters: ˜, α ψ λ γ , , , M disk (0), r c (0), β , and t . In our model, we take ˜, α ψ , and t as free parameters, letting ˜ range between 10 α -4 and 10 -2 ( e.g., Andrews et al., 2009; Rosotti, 2023), and ψ between 0.1 (turbulence-dominated) and 10 (wind-dominated). Note that α ν ∼ 10 -2 corresponds to the maximum attainable value by MRI-generated turbulence ( i.e., in the ideal MHD limit; Armitage, 2020), while α ν ∼ 10 -4 approximates the minimum attainable value by VSI-generated turbulence (Nelson et al., 2013). Below, we discuss how we derive the values for λ and γ which yield Σ( r, t ) and T r, t ( ) profiles most consistent with both viscous and irradiation heating for the range of ˜ and α ψ explored. We set M disk (0) to ≃ 0 05 . M ⊙ and r c (0) to ≃ 12 AU . These values roughly yield Σ( r = 1 AU ) close to that of the 'minimum mass solar nebula' model of Hayashi (1981) ( i.e., ≃ 2 × 10 4 kg/m ) across an accretion timescale 2 in both turbulenceand wind-dominated disks for an intermediate, fiducial value of ˜ α ≃ 10 -3 (Fig. 4). For reference, disks around Class-II young stellar objects like T-Tauri stars are inferred to have masses in the range

Log

=3.5,

∼ 3 × 10 -5 to 0 3 . M ⊙ from measurements of dust emission (Boss &amp; Ciesla, 2014 and references therein; Ansdell et al., 2017; Miotello et al., 2022 ). We also found that a β value of ≃ 1 45 yields Σ profiles that reproduce . M disk (0) to within 15% for the optimal λ and γ found across ˜ α -ψ parameter space.

## 2.2. Optimal Values for λ and γ

To constrain the power law index of the disk's temperature profile ( i.e., γ -3 / 2), we impose that the disk is optically thick and heated by both viscous shear between disk annuli and stellar irradiation. The mid-plane temperature of such a disk T r, t ( ) is established via equilibrium between heat generation and radiative loss, expressed as Armitage, 2020)

$$\sigma T ^ { 4 } \simeq \frac { 3 } { 4 } \tau F _ { v i s c } + \sigma T ^ { 4 } _ { i r r a d }, \quad \quad ( 1 1 )$$

where F visc is the heating rate per unit area from viscous shear, T irrad the interior temperature derived from models of passive radiative equilibrium disks (as opposed to blackbody disks), σ the Stefan-Boltzmann constant, and τ the optical depth. The former is given by (Nakamoto &amp; Nagakawa, 1994)

$$F _ { v i s c } = \frac { 1 } { 2 } \Sigma \nu ( r \frac { \partial \Omega _ { K } } { \partial r } ) ^ { 2 } = \frac { 9 } { 8 } \Sigma \alpha _ { \nu } c _ { s } ^ { 2 } \Omega _ { K }, \quad ( 1 2 ) \quad \text{d} \text{ } \\ \text{ as }$$

where the turbulent viscosity ν = α c h ν s (Shakura &amp; Sunyaev, 1973), and the hydrostatic scale height of the disk (assumed to be vertically isothermal at all r ) h = c / s Ω K . The functional form of T irrad depends on the behavior of the geometrical aspect ratio h/r with r , namely whether it remains virtually constant as in 'flat' disks, or increases as in 'flared' disks. The latter allows for greater interception of stellar irradiation in the outer disk, owing to its curvature away from the midplane. Analyses of spectral energy distributions from TTauri stars indicate that many disks are flared (Kenyon &amp; Hartmann, 1987), prompting us to adopt the expression for T irrad provided by Chiang &amp; Goldreich (1997) for flared radiative equilibrium disks, valid in the domain 0 4 . AU ≲ r ≲ 50 AU :

$$T _ { i r r a d } = 1 5 0 ( \frac { r } { 1 A U } ) ^ { - 3 / 7 } K. \quad \quad ( 1 3 ) \quad \ \gamma _ { l }$$

This power-law fit is obtained for a T-Tauri star with effective temperature T ∗ ≃ 4000 K , mass M ∗ ≃ 0 5 . M ⊙ , and radius R ≃ 2 5 . R ⊙ . For reference, the sun has an effective T ⊙ ≃ 5800 K .

At temperatures below ∼ 1500K ( i.e., the silicate sublimation line), dust constitutes the dominant source of opacity, and the optical depth τ takes the form (Bitsch et al., 2014)

$$\tau = \frac { 1 } { 2 } f _ { d } \Sigma \kappa _ { d },$$

where f d ≃ 0 01 is . the dust-to-gas ratio and κ d ≃ 30 m /kg the dust opacity. Bearing in mind the uncertain2 ties and complications associated with the dependence of κ d on temperature, dust size distribution, and dust composition ( e.g., Chambers, 2009), we refrain from adopting a specific model for variable κ d . We note, nonetheless, that an analytic piecewise model for κ d has been implemented by Stepinski (1998), based on fits to opacities appropriate for dust size distributions centered below the tens of microns and of solar composition (Pollack et al., 1985). The applicability of these fits, given that grain growth of mm- to cm-sized grains via hit-and-stick collisions occurs rapidly, is likely unsubstantiated.

Keeping in mind the relation c s = √ k T/µ b , Eqs. 1114 yield a quartic equation in the mid-plane temperature with coefficients given by

$$p _ { 4 } & = 1 ; \ p _ { 3 } = 0 ; \ p _ { 2 } = 0 \\ p _ { 1 } & = - \frac { 2 7 f _ { d } \kappa _ { d } \alpha _ { \nu } \Omega _ { K } k _ { b } \Sigma ^ { 2 } } { 6 4 \sigma \mu } ; \quad \quad ( 1 5 ) \\ p _ { 0 } & = - ( 7. 3 \times 1 0 ^ { 2 7 } ) \ r ^ { - 1 2 / 7 },$$

where p i corresponds to the i th-order term in T . At small r ( i.e., ≲ r c ), the temperature profile T r, t ( ) is dominated by viscous heating, and can be approximated as

$$T = \left ( \frac { 2 7 f _ { d } \kappa _ { d } \alpha _ { \nu } k _ { b } \Omega _ { K } \Sigma ^ { 2 } } { 6 4 \sigma \mu } \right ) ^ { 1 / 3 }. \quad \ \ ( 1 6 )$$

Indeed, for intermediate values of ˜ α ≃ 10 -3 and ψ ≃ 1 ( α ν = α DW ), T r ( = 10 AU ) as calculated from Eq. 16 deviates from that obtained by solving the quartic equation by only ∼ 12%. This relative deviation drops to ∼ 1% at 5 AU . At large r ( i.e., ≳ r c ), stellar irradiation constitutes the dominant heat generation mechanism, in which case T r, t ( ) can be approximated by Eq. 13 alone. For the same disk, T r ( = 20 AU ) as calculated with Eq. 13 deviates from the quartic root by ∼ 1%.

Having specified ˜ and α ψ , we seek a pair of λ and γ that would yield a Σ profile which, upon substitution into the quartic equation (Eq. 15), yields a T profile with a power law index concordant with the value of γ used ( i.e., an index of γ -3 / 2; recall T ∼ r γ -3 / 2 ). As γ enters directly into the power law index of Σ (Eq. 5), while λ does so through ξ (and weakly; Eq. 6), the former is expected to be the more crucial parameter in this exercise. For a given ˜ and α ψ , we create a raster of Σ, and thus T , profiles using values of γ ranging from 0.4 to 0.9 (corresponding to T power law indices of -0 6 . to -1 1), . and λ ranging from 2 to 5 (Tabone et al. 2022). The deviation of the power law index retrieved via a fit to one of such T profiles from its respective, assumed value of γ -3 / 2 constitute a measure

)

)

(

Log

10(

Log in Myr)

p

) (

,

p PW

/

,

p PA

10(

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-7e3f0c93-641d-4998-a8c7-64e9e863ac3e.jpg)

Log

Figure 3. Rasters of (a) γ opt and (b) the deviation (in %) in the retrieved vs. assumed temperature power law index δ across ˜ α -ψ parameter space, assuming a constant λ = 3 5 . . In (a), values of γ opt in the top-left region of parameter space ( i.e., high ˜ and α ψ ≲ 1) bounded by the dashed black lines are extrapolated (see text for discussion). Across all parameter space, γ opt ranges from ∼ 0 50 . -0 83, and . δ values are ≲ 3%, indicating a high degree of self-consistency. See Section 2.2.

of self-consistency, which we denote δ λ, γ ( ) for the sake of clarity. The values of δ λ, γ ( ) are compared, with the minimum designating the optimal pair of λ and γ for the choice of ˜ and α ψ . The results of this procedure for three cases spanning the range explored for ˜ and α ψ are shown in Fig. 2.

As expected, most of the variation in δ (given ˜ and α ψ ) is attributed to γ , motivating us to adopt an intermediate value of λ = 3 5 regardless of . ˜ and α ψ for simplicity. While it appears there is minimal variation in the optimal γ (henceforth denoted γ opt ) between Figs. 2a, 2b, &amp; 2c ( ≲ 0 05), . γ opt can vary appreciably between turbulence-dominated disks (low ψ ) of high ˜ α and wind-dominated disks (high ψ ) of low ˜. α This can be seen in Fig. 3a, which showcases γ opt across parameter space, with each pair of ˜ and α ψ (each pixel, if you will) a product of the computational exercise described above-simplified, of course, by having λ constant at 3 5. . Values of γ opt range from ≃ 0 50 . to 0 83, . corresponding to T power law indices of ≃ -0 67 to . -1 00. . In the top-left region of Fig. 3a, bounded by the black dashed lines, γ opt values are extrapolated because our algorithm yielded erratic and spurious results therein. This likely stems from inconsistencies in solving the quartic equation for T . In Fig. 3b, we depict δ values (expressed in %) obtained using respective γ opt values from Fig. 3a, all of which lie ≲ 3%.

Before moving on, we note that excluding the bottom right corner of Fig. 3a where α ν = ˜ ( α/ ψ +1) ∼ 10 -5 , γ opt ranges from ≃ 0 50 . to 0 61, . corresponding to T

indices of ≃ -0 89 to . -1 00. . This is similar to the T index of -0 9 . for the steady-state, 'actively' heated viscous disk of BM22, for which the power law index of Σ is ≃ -0 6. . Recalling that the steady-state (not exponential) term in Eq. 5 scales as r ξ -γ , it is clear that for ψ ≪ 1, the present model naturally reduces to the viscous disk model of BM22 for a wide range of ˜, as the α power law index of Σ then becomes ( ξ -γ ) ≃ -γ ∼ -0 6. . This smooth transition in the turbulent regime despite the inclusion of passive irradiation, along with the retrieval of γ opt ∼ 0 6 . for high ψ and ˜ α hints at (i) the significance of the disk at small r ( i.e., where the steady-state term in Σ and viscous heating dominate) in establishing the power law indices of Σ and T , and (ii) reaffirms the notion that the form of Σ( r, t ) is primarily controlled by γ (the 'mass ejection index' ξ is small for a wide range of ψ so long as λ ≳ 2; see Fig. 3 in Tabone et al., 2022). It appears only for α ν ≲ 10 -4 (see Eq. 16) does irradiation begin to influence the form of T r, t ( ) (higher γ corresponds to flatter T profile; Fig. 3a).

## 2.3. Fiducial Σ and T Profiles

With all parameters defined, we depict the Σ( r, t ) and T r, t ( ) profiles for two fiducial cases corresponding to (i) a turbulence-dominated disk (˜ α ≃ 10 -3 , ψ ≃ 0.1) (Figs. 4a, c) and (ii) a MHD wind-dominated disk (˜ α ≃ 10 -3 , ψ ≃ 10) (Figs. 4b, d). The obtained γ opt for these disks are ≃ 0 55 and . ≃ 0 69, . respectively, corresponding to T power law indices of ≃ 0 95 and . ≃ 0 81. . The gray dashed lines in Fig. 4a and 4b represent the Σ

Figure 4. Profiles of the surface density Σ( ) r and mid-plane temperature T r ( ) for the turbulence-dominated (a,c) and wind-dominated (b,d) fiducial disks as a function of time. Time steps (moving from blue to red) are in units of the initial accretion timescale as calculated with Eq. 10. The turbulence-dominated disk is characterized by viscous spreading and the faster inward drift of the ice-line relative to the wind-dominated disk, in which viscous spreading is virtually absent. The flared disk profile facilitates interception of stellar irradiation which prevents T r, t ( ) from plummeting with Σ( r, t ). See Section 2.3.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-ea26a9ff-63e6-4ad9-a9af-bf9bbb606582.jpg)

profiles from the model of BM22 for different values of ˙ M . Turbulence facilitates angular momentum transport between disk materials, while disk winds facilitate mass and angular momentum removal from the disk. Consequently, viscous spreading accompanies stellar accretion in the presence of turbulence, and is virtually absent is a wind-dominated disk. The extent of viscous spreading is controlled by Eq. 8, where ψ governs the evolution rate of r c ( ) with every step of t t acc, 0 . The Σ profiles in the turbulence-dominated disk match those of the steadystate viscous disk for ˙ M between 10 -8 to 10 -7 M ⊙ /yr, deviating from them past r ≃ 10 AU due to the physical cutoff of the disk implemented via the exponential decay term in Eq. 5. In the wind-dominated disk, flattening of the Σ profile at small r due to the increase in the mass ejection index ξ is resisted by higher γ opt . The increase in Σ for a given ˙ M (from BM22; driven only by α ν ), as well as the overall decrease in temperatures in the viscously-heated, smallr portion of the disk (by a factor of ∼ 2 for t = 0) reflects the decrease in α ν by close to an order of magnitude. This in turn yields ice-lines at shorter distances from the star and an increase in t Acc, 0 . The slower evolution as evidenced by the closer proximity of Σ and T profiles across time is attributed to high ψ , and can be understood by inspection of Eq. 7. This, and to a lesser extent the higher t Acc, 0 , culminates in

an inward drift rate of the ice-line that is much slower ( ∼ 1 4 AU/Myr) than that in the turbulence-dominated . disk ( ∼ 7 7 AU/Myr). .

The inclusion of stellar irradiation ensures T r, t ( ) beyond r c ( ) does not plummet with Σ( t r, t ), as would be the case if only viscous heating was present (Eq. 16). As higher temperatures correspond to greater turbulence (Eq. 17) and particle drift velocities (Eq. 18) in the Shakura-Sunyaev prescription, irradiation indirectly stalls particle growth in the outer disk (see Section 3.5).

## 3. THE STOKES NUMBER

The Stokes number is intrinsically related to particle size, and its characteristic value at any given region in the disk is defined through the competition between several barriers to particle growth therein (Birnstiel et al., 2012). Each of these barriers predicts a maximal St , and thus size, attainable by particles undergoing hit-andstick collisions, and the characteristic St is given by the minimum of these maximal values. In other words, particle growth starting at the micron-scale is halted by the first barrier it encounters, and the mechanism enforcing that barrier is location-dependent. Here, we provide an overview of the said barriers, and construct the characteristic St r, t ( ) profiles for our fiducial disks. We emphasize that the characteristic particles at a given location in the disk are those that hold most of the solid mass, not necessarily those that are most numerous. The persistence of optically-thick disks over millions of years suggests the prolonged preeminence of micron-sized dust (Williams &amp; Cieza, 2011), in turn hinting at the prevalence of particle collisions/fragmentation in disks.

## 3.1. Turbulent Fragmentation

Underlying our work here is the assumption that hitand-stick collisional growth operates at a sufficiently high efficiency such that most particles have St close to the maximal attainable value. Numerical simulations have shown that this is indeed the case when St is limited by the most generic of barriers: turbulent fragmentation (Birnstiel et al., 2011). The maximal St in this regime, denoted St tf , is expressed though the equivalence between the average relative velocity of similar-sized dust particles in turbulent motion ∆ v t = √ 3 α Stc ν s (valid for St ≲ 1) and the the fragmentation threshold velocity of those particles v f (Ormel &amp; Cuzzi, 2007):

$$S t _ { t f } = \frac { 1 } { 3 } \frac { v _ { f } ^ { 2 } } { \alpha _ { \nu } c _ { s } ^ { 2 } }$$

Laboratory experiments suggest that the onset of silicate dust fragmentation occurs at v f,rock ≃ 1 m/s (Blum &amp; M¨ unch, 1993; G¨ttler et al., 2010), and we adopt this u value for particles inward of the ice-line. Icy particles beyond the ice-line, constituting ice-rock aggregates, are characterized by surface energies greater than their silicate counterparts, and thus expected to be 'sticker,' fragmenting at higher collisional velocities (Gundlach et al., 2011; Gundlach &amp; Blum, 2015). We adopt v f,ice ≃ 5 m/s for these particles, a factor of two lower than that used by BM22. This is because in our model, v f,ice ≃ 10 m/s yields icy boulders up to the dm- to m-scale in the turbulence- and wind-dominated fiducial disks, respectively, at odds with astronomical observations suggesting dust layers composed of abundant cm-sized particles ( e.g., Testi et al., 2003; Ricci et al., 2010). The viscous disk model of BM22 retrieved cm-sized icy particles because it employed a steady-state Σ profile, that is, one without an exponential decay term that would otherwise decrease Σ and T (solely from viscous heating; Eq. 16), decrease ∆ v t , and thus increase St tf and particle sizes beyond the cm-scale. We note that v f,ice ≃ 5 m/s is substantially lower than the ≃ 15 to 20 m/s obtained assuming cm-sized particles composed of 'weak ice' from the hydrodynamic collision simulations of Leinhardt &amp; Stewart (2009). This discrepancy is insignificant to the key conclusions of our work, which pertain to the inner disk (see Section 4.3). Moreover, the assumption that ice fragments at higher velocities than silicates has been called into question by recent laboratory experiments at temperatures appropriate for most of the outer protoplanetary disk ( ≲ 150 K ), which suggests icy particles therein are comparable to their silicate counterparts in tensile strength (Gundlach et al., 2018; Musiolik &amp; Wurm, 2019; Pinilla et al., 2021). Our choice of a lower v f,ice thus reflects a step in the right direction. It also bears mention that v f,rock is expected to increase under high ( ≳ 1000 K ) temperatures (Pillich et al., 2011), and v f in general depends on porosity, composition (itself related to temperature), and size ( e.g., Benz &amp; Asphaug, 1999; Leinhardt &amp; Stewart, 2009; Beitz et al., 2011). We omit these complications here for simplicity, but briefly discuss them in the context of particle sizes obtained in Section 4.1.2. We consider the potential implications of the temperature dependence more deeply in Section 5.

## 3.2. Radial Drift

Drift towards the Sun, or host star, poses another barrier to particle growth. The H -He gas composing 2 the bulk of protoplanetary disks is partially supported against gravity by a radial pressure gradient, and thus orbits the central star at sub-Keplerian azimuthal velocities. The larger a particle grows (for St ≲ 1), the more susceptible it is to drag from headwind induced by the sub-Keplerian gas flow around it (Weidenschilling,

1977). The increase in drag ( i.e., angular momentum loss per unit time) accelerates the particle's drift inward. Thus, a characteristic particle size (and St ) for a local dust population may be defined through a balance between hit-and-stick collisional growth and radial drift. Consideration of the radial and azimuthal equations of motion for an orbiting particle undergoing drag forces leads to the well-known expression for its radial drift velocity, given by (Nagakawa et al., 1986)

$$v _ { r } = \frac { S t ^ { - 1 } v _ { r, g a s } + | \epsilon | ( c _ { s } ^ { 2 } / v _ { K } ) } { S t + S t ^ { - 1 } }, \quad \ \ ( 1 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where | ϵ | is the power law index of the disk's midplane pressure profile P = ρ c g 2 s ∼ r ϵ , ρ g = Σ ( / √ 2 πh ) being the gas density. The two terms on the RHS represent the contribution to v r from advection with the accretionary gas flow, and the sub-Keplerian headwind, respectively. Recall that v r,gas = (3 2)(2 / -γ ) 2 ˜ αc v 2 s -1 K , and note that for (aerodynamically) small particles that are tightly coupled to the surrounding gas ( St &lt;&lt; 1), v r ≃ v r,gas meaning these particles are virtually advected towards the central star with the gas. At very large sizes ( St &gt;&gt; 1), keeping in mind acceleration from drag scales with the inverse of particle size/radius s ( ∼ 1 /s ), particles/boulders decouple from the gas with v r ≃ 0. In between these two extremes, v r rises to a maximal value of ≃ (1 / 2) | ϵ c | 2 s v -1 K , which occurs at St ≃ 1. In the turbulence- and wind-dominated disks at t = 0 and r = 1 AU, this evaluates to ≃ 140 m/s and ≃ 60 m/s, respectively. Physically, the peak in v r describes particles that are sufficiently coupled to the gas to travel at a sub-Keplerian velocity, and yet have enough mass such that the centrifugal force acting on them fails to counteract gravity. This peak constitutes the so-called 'meter-size barrier' (Weidenschilling &amp; Cuzzi, 1993; Blum &amp; Wurm, 2008), which depletes the local disk region of St ≃ 1 particles rapidly relative to the growth timescale ( τ grow ∼ s/s ˙) of those particles. This is made clear with an example. Assuming particles are in the Epstein drag regime, applicable for sizes s ≲ 9 / 4 times the mean free path of gas molecules λ mfp (Armitage, 2020; see Section 4.1.2 below), the growth timescale takes the simple form τ grow = 1 Ω / K d f , independent of s (Birnstiel et al., 2012). The 'meter-size' drift timescale τ drift , assuming a 'locality' in the disk has a radial (bin) width of ≃ 1 AU, is simply the ratio of 1 AU to the peak v r . Evaluation of τ drift /τ grow at r ≃ 20 AU yields ∼ 0 1 . for both the turbulence- and wind-dominated disks, meaning that St ≃ 1 particles therein are doubling their mass about ten times slower than the time it takes for them to drift across a hundred million kilometers. For simplicity, we set the maximal St in the drift-limited regime, denoted St rd , to unity.

## 3.3. Fragmentation from Relative Drift Velocities

In disks with low turbulence, fragmentation of colliding particles may be facilitated by their relative drift velocities (Birnstiel et al., 2012). The maximal St in this alternate fragmentation-limited regime, denoted St df , is expressed through an equivalence between some characteristic relative drift velocity between similar-sized particles ∆ v r and the fragmentation threshold velocity v f . The radial drift velocity of a particle in the disk is given by Eq. 18. Assuming the smaller of the colliding particles has St 2 = NSt 1 , where N &lt; 1, their relative drift (or collision) velocity is

$$\Delta v _ { r } & = \frac { c _ { s } ^ { 2 } } { v _ { K } } \left ( \frac { 3 \tilde { \alpha } ( 2 - \gamma ) ^ { 2 } / 2 + | \epsilon | S t _ { 1 } } { S t _ { 1 } ^ { 2 } + 1 } \\ - & \frac { 3 \tilde { \alpha } ( 2 - \gamma ) ^ { 2 } / 2 + | \epsilon | N S t _ { 1 } } { N ^ { 2 } S t _ { 1 } ^ { 2 } + 1 } \right ). \\ \bar { \ } _ { \ } _ { \ } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array} \quad ( 1 9 )$$

Setting ∆ v r = v f and rearranging for St 1 ( i.e., St df ), Eq. 26 yields a quartic equation with coefficients given by

$$\nu _ { j } \\ & p _ { 4 } = \frac { v _ { f } v _ { K } N ^ { 2 } } { c _ { s } ^ { 2 } } ; \ p _ { 3 } = | \epsilon | ( N - N ^ { 2 } ) \\ p _ { 2 } = \frac { v _ { f } v _ { K } ( N ^ { 2 } + 1 ) } { c _ { s } ^ { 2 } } + \frac { 3 \tilde { \alpha } ( 2 - \gamma ) ^ { 2 } ( 1 - N ^ { 2 } ) } { 2 } \quad ( 2 0 ) \\ & p _ { 1 } = | \epsilon | ( N - 1 ) ; \ p _ { 0 } = \frac { v _ { f } v _ { K } } { c _ { s } ^ { 2 } } \\ \text{where $n$ corresponds to the third order term in $St$} \,. \ \, \text{We}$$

where p i corresponds to the i th-order term in St df . We set N = 0.5 (Birnstiel et al., 2012), and obtain the solution to St df at each location in the disk numerically as done for T r, t ( ) (Eq. 5). As in the case of turbulent fragmentation, we adopt v f,rock ≃ 1 m/s and v f,ice ≃ 5 m/s for particles inward and outward of the ice-line, respectively. To ensure the quartic equation at each location in our fiducial disks is solved, we substituted the St df ( r, t ) obtained back into Eq. 19 and checked that ∆ v r ( r, t ) ≃ v f ( r, t ).

## 3.4. Bouncing Collisions

A host of laboratory experiments have demonstrated that a possible outcome of particle collisions, in addition to sticking and fragmentation, is bouncing (Blum &amp; M¨ unch, 1993; G¨ uttler et al., 2010). Dust collision models suggests that bouncing can halt growth before particle velocities are sufficiently high for fragmentation to be relevant (Zsom et al., 2010). Despite experimental support for such a 'bouncing barrier,' its significance to dust evolution in protoplanetary disks remains nebulous. Numerical simulations indicate that bouncing may only occur during collisions involving compact aggregates with relatively high coordination numbers ( i.e.,

Figure 5. Profiles of the maximal Stokes number predicted by competing barriers to particle growth at t = 0 (a &amp; b) and the characteristic Stokes number over time (c &amp; d) for the two fiducial disks. In (a) and (b), the grey band defines the minimum of the three maximal St profiles ( i.e., the characteristic St profile). As in Fig. 4, time steps in units of t Acc, 0 are represented in blue to red. Turbulent fragmentation is the dominant growth-limiting mechanism in the turbulence-dominated disk and inner wind-dominated disk. In the latter, fragmentation from relative drift velocities limits growth beyond ∼ 3 1 AU. The ice-line constitutes a major transition point in . St , owing to the higher fragmentation threshold velocity of ice ( ≃ 5 m/s) assumed. In (c), the flared disk profile keeps temperatures high in the outer disk, prohibiting a rapid rise in St tf (Eq. 17). In (d), low temperatures, and hence turbulent velocities, result in higher St values in the inner disk. Points on St profiles represent the transition from St tf to St df at each time step. In both disks, highSt particles drift inward with the cooling of the disk. See Section 3.5.1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-a0a7c89a-b99a-4ce3-b7b3-3311258a00f0.jpg)

the number of neighboring contacts per particles in the aggregate), unrepresentative of dust particles characterized by low bulk densities ( ≲ 100 kg/m ) as generally 3 expected in disks (Wada et al., 2011). It is important to note, however, that the bulk densities of such particles remains poorly constrained. It is possible, for instance, that the majority of particles were sub-mm- to mm-sized chondrules, or at least chondrule-like in their densities ( ρ ◦ ≳ 3000 kg/m ; 2 e.g., Friedrich et al., 2015 and references therein) such that bouncing was significant. Such a proposition is supported by the abundance of chon- drules in primitive meteorites from both the inner and outer SS. The caveat here is that chondrules may not be as ubiquitous in the SS as the meteorite record suggests. In particular, it has been suggested that most of the material in the SS is CI-(Carbonaceous Ivuna)like, with densities and tensile strengths too low to survive passage through Earth's atmosphere (Sears, 1998; Yap &amp; Tissot, 2023). This view is supported by the similarity between CI chondrites and asteroid Ryugu samples returned by Hayabusa2 (Yokoyama et al., 2022), as well as in-situ thermal inertia measurements of asteroid Bennu's sur-

)

, 0

Acc

, 0)

Acc

t

t

Time (units of

Time (units of

)

, 0

Acc

t

Time (units of

1.0 0

1.

0.8 8

0.

0.6 6

0.

0.4 4

0.

0.2 2

0.

0.0 0

1.0

0.

0.8

0.6

0.4

0.2

0.0

face, indicating abundant high-porosity boulders (Rozitis et al., 2020). Even with bouncing in effect, it is unclear how easily it may be overcome. Windmark et al. (2012) postulated that cm-sized particles may catalyze growth beyond the mm/sub-mm-scale to which particles are confined by bouncing collisions, and recently, Steinpilz et al. (2020) argued that the bouncing barrier is insignificant in consideration of triboelectric charging of colliding particles.

Given the uncertainty surrounding the bouncing barrier as well as its postulated precedence over fragmentation and radial drift, we do not include it in our primary analysis, but assess its implications for our results separately and only in the turbulence-dominated disk. Here, we derive the maximal St attainable in the bouncinglimited regime (denoted St b ) by equating the turbulent velocity dispersion of dust particles ∆ v t = √ 3 α Stc ν s to the mass-dependent bouncing threshold velocity ∆ v b deduced by Weidling et al. (2012) from collision experiments. With the target mass ≳ the mass of the projectile m proj , ∆ v b (in m/s) is given by

$$\Delta v _ { b } = 0. 0 1 \left ( \frac { m _ { p r o j } } { m _ { b } } \right ) ^ { - 5 / 1 8 }, \quad \quad ( 2 1 ) \quad \ \ 3$$

where the normalizing constant is m b is calibrated at ≃ 3 3 . × 10 -6 kg. We assume that particles limited by bouncing are confined to the Epstein drag regime and confirm this assumption a posteriori . In this regime, the relationship between St and s is given by

$$S t ( s ) _ { E p s } = \sqrt { \frac { \pi } { 8 } } \frac { s \rho _ { \circ } \Omega _ { K } } { \rho _ { g } c _ { s } }. \quad \quad ( 2 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Here, ρ ◦ is the material density of the particle, for which we assumed typical values for compact rock and waterice aggregates ( ρ ◦ ,rock ≃ 3300 kg/m ; 3 ρ ◦ ,ice ≃ 1000 kg/m ) given the aforementioned uncertainties in this 3 parameter. Note that ρ ◦ does not enter into the expressions derived for St tf (Eq. 17) nor St df (Eq. 20). Assuming spherical particles with radii s , St b takes the form

$$S t _ { b } \simeq 0. 0 1 \left ( \frac { m _ { b } \rho _ { o } ^ { 2 } \Omega _ { K } ^ { 3 } } { \rho _ { g } ^ { 3 } \alpha _ { \nu } ^ { 9 / 5 } c _ { s } ^ { 3 3 / 5 } } \right ) ^ { 5 / 2 4 }. \quad \ \ ( 2 3 ) \quad \ n _ { \mathrm d }$$

$$3. 5. \ S t o k e s \ N u m b e r \ P r o f i l e s \\ 3. 5. 1. \ W i t h o u t \ B o u n c i n g$$

The St α, ψ, r (˜ ) profiles constructed for our fiducial cases are shown in Fig. 5. As BM22 have shown, St rises dramatically across the ice-line owing to the higher fragmentation velocities of icy relative to silicate particles. In the turbulence-dominated case, it is clear

1.0

Figure 6. Profiles of the characteristic Stokes number in the bouncing-limited turbulent disk over time. The bouncing barrier limits particles to much lower St than turbulent fragmentation. Here, the transition across the ice-line reflects the change in particle material density, as opposed to fragmentation threshold velocities as in Fig. 5. Viscous spreading increases ρ g in the outer disk over time (Eq. 23), leading to lower St . See Section 3.5.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-6a29966e-74e9-449e-8694-f8441ec8511d.jpg)

that turbulent fragmentation sets the characteristic St throughout the disk (Fig. 5a). While particles interior to the ice-line are confined to St &lt; 10 -3 , those beyond it possess St values ≳ 10 -2 , indicating strong dust-gas coupling in the former, and weak coupling in the latter. Stellar irradiation of the flared outer disk prohibits temperatures from plummeting with Σ and hence St tf from rising dramatically past r c ( ). t The lack of a decrease in St over time for r ≳ r c ( ) t (Fig. 5c), as would be expected from viscous spreading (which increases Σ and T in the outer disk, translating to to faster turbulent velocities in the Shakura-Sunyaev prescription), reaffirms the dominance of irradiation over viscous heating therein.

The maximal St predicted by the two fragmentation barriers increase as turbulence gives way to disk winds (Fig. 5b, d). In the wind-dominated disk, particles possess characteristic St values about an order of magnitude greater than those in the turbulence-dominated disk. Silicate particles in the inner disk are confined to St &lt; 10 -2 , while icy particles in the outer disk possess St ∼ 0 1. . Although turbulent fragmentation remains the St -limiting process for the former, fragmentation due to relative drift velocities dominates beyond ∼ 3 1 . AU and limits the growth of the latter. In both fiducial disks, highSt icy particles drift inwards over time with the ice-line, owing to the cooling of the disk accompa-

)

, 0

Acc

t

Time (units of

1.0

0.8

0.6

0.4

0.2

0.0

nying the decrease in Σ (Eq. 7).

## 3.5.2. With Bouncing

The St α, ψ, r (˜ ) profile for the turbulence-dominated disk with the inclusion of the bouncing barrier is shown in Fig. 6. As expected, bouncing sets the characteristic St across the entire disk, confining particles up to r c ∼ 20 AU below St ≃ 10 -4 and limiting growth far before turbulent fragmentation takes effect. The slight depression in St across the sublimation line reflects the change in particle density, and viscous spreading leads to lower St over time in the outer disk, owing to the increase in Σ and thus ρ g (Eq. 23). The parametrization of the bouncing threshold velocity given in Eq. 21 is based on experiments on millimeter-sized aggregates composed of micron-sized poly-disperse SiO 2 particles. Assuming icy particles are 'stickier' than silicates, it is likely that ∆ v b,ice ≳ ∆ v b,rock . This would result in a dichotomy in St b akin to that for St tf and St df in our two fiducial disks.

## 4. RESULTS &amp; DISCUSSION

## 4.1. Dust Properties

## 4.1.1. Dust Scale Height &amp; Drift Velocities

The St α, ψ, r (˜ ) profiles derived above permit determination of the vertical distribution, radial drift rate, and physical size of the characteristic dust particles across their disks. Here, we compare how these three features vary between the turbulence- and wind-dominated disks absent of the bouncing barrier, before extending our analysis to the bouncing-limited, turbulent disk.

At the macroscopic level, dust in protoplanetary disks can be treated as a trace fluid separate from the bulk H -He gas, modulo the Schmidt number ( 2 S c ≃ 1) and regions where pressure bumps may lead to ρ ◦ ≳ ρ g . By virtue of a balance between turbulent diffusion and gravitational settling, dust particles establish a solid subdisk with a St -dependent scale height h ◦ ( ≲ h ) given by (Dubrulle et al., 1995)

$$h _ { \circ } = \frac { h } { \sqrt { 1 + S t / \alpha _ { \nu } } }. \quad \quad ( 2 4 ) \quad \colon$$

As mentioned above, most of the dust mass is assumed to be held in particles with St close to the maximal value allowed by the growth barriers described. However, the profile of h 0 as calculated using the characteristic St values from above is not that which dominates the observational appearance of disks. The optical thickness of a disk is primarily attributed to micron-sized dust, which, tightly coupled to the gas as they are (with St ≪ 1), have h ◦ ≃ h . Fig. 7 depicts the profiles of the normalized scale height h /h ◦ for our two fiducial cases. Plotted as well are profiles of the normalized radial drift rate | v /v r K | as calculated using Eq. 18.

In both the turbulenceand wind-dominated disks, the water ice sublimation line constitutes a critical transition point in the degree of dust settling and particle drift regimes. Generally, silicate dust particles in the inner disk are well coupled to their surrounding gas, and form a relatively dispersed and well-mixed layer ( h ◦ ≳ 0 1 . h ) about the disk midplane. Their inward drift is predominantly facilitated by advection with the accretionary gas flow, and so proceeds at rather slow velocities ( &lt; 10 -4 v K ). In the turbulence-dominated disk (Fig. 7a), strong dust-gas coupling ( St &lt; 10 -3 ) in the inner disk manifests as a virtual equivalence between h 0 and h , and between v r and v r,gas . The increase in disk winds, and hence St , leads to greater settling ( h 0 ≲ 0 5 . h ) and faster (by a factor of a few) inward drift relative to the gas (Fig. 7b).

Icy particles in the outer disk, being larger and loosely coupled to the gas, settle into a thin layer ( h 0 ≲ 0 1 . h ) at the mid-plane. Their inward drift is driven primarily by headwind drag, and proceeds at velocities that are substantially greater fractions ( &gt; 10 -4 v K ). As with silicate particles, icy particles in the wind-dominated disk settle more readily, and drift faster, than those in the turbulence-dominated disk.

As St df tapers off at ∼ 0 1 in both fiducial disks (Figs. . 5a &amp; 5b), only past the point where there is no solution to St df can St approach unity ( i.e., St rd ) in the outer disk. While this point may be too far in the outskirts of the disk to be meaningful, it is bound to exist since c 2 s /v K falls with r (Eq. 19), unless v f falls equally fast or faster (say due to the aforementioned temperature dependence). Where St rd becomes the dominant growth limiting process, v r will take on its peak value of ∼ (1 / 2) | ϵ c | 2 s /v K .

Figs. 8a &amp; 8b shows the profiles of h /h 0 and | v /v r K | for the bouncing-limited turbulent disk. As expected from the small St s / (Figs. 6, 8c) to which the bouncing barrier relegates particles across the disk, the scale height and drift velocity of dust is essentially identical to those of the gas. Only beyond the outer edge at ∼ r c do particles begin to settle and drift faster than their surrounding gas. Even then, h ◦ ≳ 0 5 . h and v r is only a factor of a few greater than v r,gas .

## 4.1.2. Dust Particle Sizes

The relationship between St and particle size/radius s is dependent on the drag law in effect. In the Epstein drag regime, valid for s ≲ 9 / 4 λ mfp , St s ( ) is given by Eq.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-64744f4c-2ee5-457a-802c-c77ce8b5059c.jpg)

)

, 0

)

, 0

Acc

t

Acc

t

Time (units of

Time (units of

10 0

10 1

Heliocentric distance r (AU) Figure 7. Profiles of the normalized scale height and radial drift velocity of the characteristic dust particles in the (a) turbulence- and (b) wind-dominated disks. Grey bands represent the upper and lower bounds to the drift velocity for all six time steps, the former corresponding to St ≃ 1. The degree of dust settling as well as the headwind contribution to the dust drift velocity change dramatically across the ice-line. In general, silicate particles in the inner disk are vertically dispersed and drift inward with at velocities close to that of the surrounding gas. Icy particles concentrate in a thin sub-disk at the mid-plane, and drift inwards faster than the surrounding gas. See Section 4.1.1.

Heliocentric distance r (AU)

22. In the Stokes drag regime, for which s ≳ 9 / 4 λ mfp , St s ( ) is instead given by (Armitage, 2020)

$$S t ( s ) _ { S t s } = \frac { 8 } { 3 } \frac { s \rho _ { \circ } \Omega _ { K } } { C _ { D } \rho _ { g } v _ { r e l } }, \quad \quad ( 2 5 ) \quad \ \ 1$$

where v rel = v K [1 -(1 - | | ϵ c 2 s /v 2 K ) 1 / 2 ] is the relative velocity between dust and gas ( i.e., the headwind velocity), and C D the drag coefficient dependent on the fluid Reynolds number Re , given by

$$\nu \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } { \dots } \stackrel { \circ } \dots } { \dots } \stackrel { \circ } \dots } R e = \frac { \sqrt { 2 \pi } s v _ { r e l } n \sigma _ { m o l } } { c _ { s } }. \quad \ \ ( 2 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad R e = \frac { \text{ol} } { \text{is} } \quad \text{ol} }$$

In the second equality, we have replaced λ mfp with ( nσ mol ) -1 , n = ρ /µ g being the number density of gas molecules, and σ mol ∼ 2 × 10 -19 m 2 the cross-section for ( H 2 ) molecular collisions (Armitage, 2020). For Re ≲ 1 (Stokes Regime I), C D ≃ 24 Re -1 , while for 1 ≲ Re ≲ 800 (Stokes Regime II), C D ≃ 24 Re -0 6 . . For Re ≳ 800 (Stokes Regime III), C D ≃ 0 44 . (Whipple 1972; Weidenschilling 1977). In constructing s r, t ( ) profiles for our two fiducial disks, depicted in Figs. 9a &amp; 9b, care was taken to ensure self-consistency in the drag law used at each location r . The s r, t ( ) profiles for the bouncing-limited turbulent disk is depicted in Fig. 8c.

In the turbulence-dominated disk, all particles are confined to the Epstein drag regime. While the inner disk is characterized by sub-mm-scale silicate particles, the outer disk hosts icy particles millimeters to centimeters in size. Particles in the wind-dominated disk are universally larger by ∼ two orders of magnitude, with ∼ one-cm particles and cm-to-decimeter-scale particle/boulders in inner and outer disk, respectively. Here, with the exception of the region around r ∼ 2 5 . AU to ∼ 5 AU, wherein particles are in the Stokes regime, all other particles are in the Epstein regime. Particles around r ∼ 3 5 . AU transition from Stokes Regime II to I over an initial accretion timescale t Acc, 0 of ∼ 0 5 . Myr (Fig. 4b). Particles in the bouncinglimited turbulent case are, of course, deep within the Epstein regime, with sizes in the tens of microns. Although not shown, bouncing-limited icy particles in the wind-dominated disk reach sizes of around a hundred microns.

The large, decimeter-sized icy boulders predicted right past the ice-line in the wind-dominated disk may be an artifact of assigning a fixed fragmentation threshold velocity v f ≃ 5 m/s to all icy particles beyond the ice-line, regardless of locally variable disk conditions. Stated differently, they reflect negligence of the effects temperature and particle size have on v f . We note, nonetheless, that experimental and theoretical efforts to understand particle growth in disks have shown that dust grains can feasibly grow up to decimeters in size by hit-and-stick collisions (Blum &amp; Wurm, 2008).

With respect to particle size, for colliding particles ≲ 100 meters in size, the threshold for catastrophic dis-

1.0

1.0

= 0.1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0.0

0.0

)

, 0

Acc

t

Time (units of

ruption Q D ∼ v 2 f ∼ St tf /St df = f ( s ) decreases with s within the so-called 'strength dominated' regime (Benz &amp;Asphaug, 1999; Leinhardt &amp; Stewart, 2009; Armitage, 2020). That is, the larger a particle grows, the less resistant it is to fragmentation ( i.e., lower v f ). The most self-consistent way to evaluate s at each location in the disk (when growth is fragmentation-limited) is to seek out a value of v f ∼ √ Q D there which, along with local disk conditions, yields a value of s concordant with the Q D used. The implementation of this procedure, of course, assumes a reliable relationship between Q D and s has been established. Given (i) the uncertainties that still remain concerning the nature and behavior of colliding dust particles in protoplanetary disks and (ii) the observation that Q D seems to change by only a factor of a few with a change in s over two orders of magnitude (Leinhardt &amp; Stewart, 2009), we had omitted the size dependence of v f . More crucial to v f and thus the dubious decimeter boulders we predict for the winddominated disk is the temperature dependence. If v f,ice evolves quickly to ∼ v f,rock past the ice-line, most of the icy particles at r ≲ 20 AU would be mm- and cm-sized in the turbulence- and wind-dominated disks, respectively. While a dichotomy in St would be absent (apart from a 'bump' across the ice-line; see Section 5), a mild dichotomy in s would still persist, reflecting a change in dust material density ρ ◦ . In this case, however, particle sizes fall to below a millimeter rather rapidly, between ∼ 20 to 30 AU given our choice of r c (0).

Turning our attention to silicate particles in the inner disk, we note that the sizes obtained in the turbulence-dominated disk (upper hundreds of microns) closely match those of BM22 for a viscous disk with α ν ≃ 10 -3 and ˙ M ≃ 10 -8 M ⊙ /yr. In their model, silicate particles are just short of a millimeter in size, a factor of order unity higher than our predictions. This slight discrepancy reflects the lower Σ values in their model for the choice of ˙ M (Fig. 5a), which results in lower T and thus higher St tf (Eq. 17).

The growth of silicate particles is limited by turbulent fragmentation, which is expected to operate at higher St ( i.e., larger particle sizes) than the bouncing barrier. However, the collision outcome model of Windmark et al. (2012) predicts that for similar-sized colliding particles, perfect bouncing ( i.e., no sticking whatsoever) is achieved for particles hundreds of microns to a millimeter in size, with fragmentation setting in beyond the mm-scale. This seems to suggests the sizes we have derived are underestimated. The value of v f,rock ≃ 1 m/s is based on collision experiments performed on millimeter-sized grains composed of ∼ 0 5 micron par-. ticles ZrSiO 4 in free-fall and vacuum (Blum &amp; M¨ unch,

1993), and is possibly unrepresentative of particles in actual disks. Due to sintering effects at high ( T ∼ 1000 K ) temperatures in the inner disk (Poppe, 2003), accompanied by dehydration and possible phase changes (Pillich et al., 2021, 2023), silicate particles therein may have v f,rock &gt; 1 m/s. Moreover, v f,rock could increase if progressive compaction of silicate particles ( i.e., evolution in their material density) occurs (Blum &amp; Wurm, 2008), and does so on a meaningful timescale, say, less than the radial drift timescale. We revisit the discussion of the v f temperature dependence in Section 5. Before moving on, note that regardless of whether particles are in the Epstein (Eq. 22) or Stokes (Eqs. 25 &amp; 26) regimes, higher ρ 0 lead to smaller s . Thus, if ρ ◦ ,rock &lt; 3300 kg/m 3 as assumed, silicate particles would be larger than presently predicted.

The particle sizes predicted in our bouncing-limited turbulent disk, in the tens of microns, also appear inconsistent with the expectation that perfect bouncing occurs for particles exceeding a hundred microns. Nonetheless, there is little cause for concern due to two main reasons. First, our results are obtained for a disk with ˜ α ≃ 10 -3 and ψ ≃ 0 1. . As can be deduced from Eq. 23 and the observation that particles sizes increase ψ (due to the decrease in temperatures; Fig. 9), a decrease in ˜ and/or a transition to a α more wind-dominated disk will lead to higher bouncinglimited particle sizes. Second, in the model of Windmark et al. (2012), similar-sized colliding particles between ten to a hundred microns are not strictly within the sticking regime. Experiments suggests there is no sharp boundary between sticking and bouncing regimes, but a so-called 'sticking-to-bouncing' transition in parameter space, where collision outcomes are modeled with a probability distribution (Weidling et al., 2012). This indicates at least that collision outcomes and the parameters that control them remain uncertain. We note that the tens of microns we have obtained are consistent with the fragmentation-limited particle sizes from the fiducial turbulence-dominated disk, in the upper hundreds of microns.

Inspection of Fig. 9 suggests the closest match to inferences in particle sizes from astronomical observations may be found for ψ ∼ 1, in between the two extreme cases we have considered. For ψ ≃ 1 and ˜ α ≃ 10 -3 , silicate and icy particles up to r ≃ 40 AU are confined between a millimeter and decimeter.

## 4.1.3. A Note on Planetesimal Formation

The low characteristic St of silicate particles within the ice-line ( &lt; 10 -2 ) do not appear conducive to efficient planetesimal formation via streaming (or more generally,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-beb4bf3b-962b-42ea-aa48-6d1ed762eaac.jpg)

)

, 0

Acc

t

Time (units of

䡥汩潣敮瑲楣⁤楳瑡湣攠爠⡁唩

Figure 8. Profiles of the (a) normalized dust scale height, (b) normalized dust radial drift velocity, and (c) characteristic particle size in the bouncing-limited turbulent disk. Interior to the outer edge of the disk at r = r c , the dust particles tens of microns in size are completely coupled to the gas. The scale height of the dust sub-disk is equivalent to the hydrostatic scale height of the gas, and dust particles drift inward with the accretionary gas flow. These particles sit deep within the Epstein drag regime. See Sections 4.1.1 &amp; 4.1.2.

resonant drag) instabilities, which develop readily at higher St (Youdin &amp; Goodman, 2005; Squire &amp; Hopkins, 2018). Our model, however, does not consider the potential for radial concentration of dust by pressure bumps ( e.g., associated with the silicate sublimation line; Morbidelli et al., 2022), and other substructures that would hasten the growth rate of instabilities, which depend on both St and the dust-to-gas ratio (metallicity) f d . An increase in f d beckons an increase in St as turbulence becomes less capable, of sufficiently energetic, to combat the settling of dust towards the midplane (see Section 4 in Batygin &amp; Morbidelli, 2020). Settled dust particles grow to larger sizes/ St owing to a smaller turbulent velocity dispersion ∆ v t . In Section 4.3.3 below, we discuss how the buildup of dust at the silicate sublimation line can result in planetesimal 'rings' in which rocky planets may be efficiently accreted ( e.g., Batygin &amp; Morbidelli, 2023).

## 4.2. 2D vs. 3D Pebble Accretion

Pebble accretion (PA) is the process by which the gravity of a planetary body detaches particles with 10 -4 ≲ St ≲ 1 from the background flow of gas, inducing considerable drag which causes them to spiral into the body. In the ensuing analysis and discussion, we use the term 'planetary embryo' to denote both kmto 100-km-sized bodies formed via the gravitational collapse of pebble clouds ( i.e., planetesimals) and larger, 1000-km-sized bodies which have outgrown the planetesimal stage ( i.e., protoplanets). A common misconception about PA, perhaps stemming from its success in resolving the so-called 'timescale conflict' in gas/ice- giant core accretion (Lambrechts &amp; Johansen, 2012), is that it is always efficient. After all, the collisional cross sections characteristic of PA are generally perceived to exceed those in the classical, Safronov-type model for planetary growth involving pairwise planetesimal collisions (Safronov, 1969; Chambers et al., 1998). In reality, the mass accretion rate ˙ m p of an planetary embryo undergoing PA is strongly dependent on the vertical distribution of the surrounding pebbles, namely whether they are concentrated in a thin layer, or diffusely distributed about the disk mid-plane (Ormel, 2017). In the former case, PA proceeds efficiently in the 2D regime, where (in the headwind/Bondi sub-regime; see below)

$$\dot { m } _ { p, 2 D, B o n d i } \simeq f _ { d } \Sigma \sqrt { 8 G m _ { p } S t v _ { r e l } / \Omega _ { K } }, \quad ( 2 7 )$$

while in the latter case, PA proceeds more slowly in the 3D regime, where

$$\dot { \ m } _ { p, 3 D } \simeq 6 \pi f _ { d } \Sigma R _ { H } ^ { 3 } S t \Omega _ { K } / \sqrt { 2 \pi } h _ { \circ }, \quad ( 2 8 )$$

and the planetary Hill radius R H = (10 r Stm / M p 3 ∗ ) 1 / 3 , M ∗ being the mass of the central star (Ormel, 2017; Fig. 10). Recall from above that v rel = v K [1 -(1 -| ϵ c | 2 s /v 2 K ) 1 / 2 ]. The critical planetary-to-stellar mass ratio for the transition between the two regimes is obtained by setting ˙ m p, 2 D,Bondi = ˙ m p, 3 D , which yields

$$\left ( \frac { m _ { p } } { M _ { * } } \right ) _ { 2 D - 3 D } \simeq \frac { 4 \alpha _ { \nu } | \epsilon | h ^ { 4 } } { \pi r ^ { 4 } S t ( S t + \alpha _ { \nu } ) }. \quad \ \ ( 2 9 )$$

Above this critical ratio, PA operates in the 2D regime. 2D pebble accretion is further divided into two subregimes, depending on the relevant length scale over

1.0

ㄮ

0.8

〮

0.6

〮

0.4

〮

0.2

〮

0.0

〮

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-e3181580-5183-48db-977a-b405d9f9e608.jpg)

)

, 0

Acc

t

1.0

0.8

Time (units of

0.6

0.4

0.2

0.0

Figure 9. Profiles of the characteristic particle size for the (a) turbulence- and (b) wind-dominated disks. Grey bands represent the boundary between the Epstein and Stokes drag regimes across the six time steps, defined by s = (9 / 4) λ mfp . In (a), silicate and icy particles across the disk are in the Epstein drag regime. Sizes of the former are confined to the sub-mm scale, while those of the latter range from millimeters to centimeters in size. In (b), only particles between r ∼ 2 5 . AU to ∼ 5 AU enter the Stokes drag regime. A small portion of this region begins in Stokes Regime II, but evolves into Stokes Regime I within the initial accretion timescale. All other particles fall within the Epstein regime. Silicate particles in the inner disk are millimeters to a centimeter in size, while icy particles/boulders in the outer disk are in the cm- to dm-scales. See Section 4.1.2.

which pebble capture ensues. In the headwind or Bondi regime, (Eq. 28) this length scale is the Bondi radius R B = Gm c p s √ StM /m ∗ p | ϵ /v | 3 K , whereas in the (generally) more efficient shear or Hill regime, it is the Hill radius. For the latter, ˙ m p, 2 D is given by (Ormel, 2017)

$$\dot { m } _ { p, 2 D, H i l l } \simeq 2 f _ { d } \Sigma R _ { H } ^ { 2 } \Omega _ { K } S t ^ { 2 / 3 }. \quad \ \ ( 3 0 )$$

accreting protoplanet becomes comparable to the (local) hydrostatic scale height of the disk ( R H ∼ h ) (Lambrechts et al., 2014). At this stage, the protoplanet is sufficiently large to perturb the disk structure, carving out a gap of high pressure that halts the inward drift of pebbles. This 'pebble isolation mass' takes the form:

The critical mass ratio demarcating the boundary between these two sub-regimes is obtained by setting R B = R H , yielding:

$$\left ( \frac { m _ { p } } { M _ { * } } \right ) _ { B - H } \simeq \frac { 1 0 0 | \epsilon | ^ { 3 } c _ { s } ^ { 6 } } { 9 S t v _ { K } ^ { 6 } }. \quad \quad ( 3 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

As the planetary embryo grows, its Bondi radius will eventually surpass its Hill radius, at which point PA transitions from the headwind to shear regime ( e.g., Ida et al., 2016; Ormel, 2017). Moving forward, it is important to keep in mind that the PA rates introduced serve only as approximations, and overestimates for that matter. Indeed, they assume a maximal pebble impact parameter, obtained by equating the pebble settling time ( i.e., the time it takes for a pebble to settle onto the embryo from gravity and drag) and the pebble-embryo encounter time ( i.e., the characteristic time over which the embryo exerts the greatest gravitational tug on the pebble).

Pebble accretion does not proceed indefinitely. It is self-limiting, and terminates when the Hill radius of the

$$\left ( \frac { m _ { p } } { M _ { * } } \right ) _ { I s o } \simeq \frac { 1 } { 2 } \left ( \frac { h } { r } \right ) ^ { 3 } \text{ \quad \quad } ( 3 2 )$$

For regions of a disk with ( m /M p ∗ ) Iso &lt; ( m /M p ∗ ) 2 D -3 D , PA is confined to the 3D regime for all planetary masses.

Profiles of these transitory mass scales for our two fiducial disks are shown in Fig. 11, where ( m /M p ∗ ) 2 D -3 D is represented by solid lines, ( m /M p ∗ ) B -H in faded lines, and ( m /M p ∗ ) Iso in dotted lines. Here and in Section 4.3, a Sun-like star ( i.e., M ∗ ≃ 2 × 10 30 kg) is assumed. The results for the bouncing-limited turbulent case are not shown, as they are readily inferred from the conclusions drawn for silicate particles in the turbulence-dominated disk (Fig. 11a). A salient feature of Fig. 10 is that, regardless of whether disk evolution is driven by turbulence or MHD winds, planetary embryos interior to the ice-line are prohibited from partaking in 2DPA as they would reach the pebble isolation mass ( m /M p ∗ ∼ 10 -5 to 10 -4 ) before the transition point between the 3D an 2D regimes ( m /M p ∗ &gt; 10 -4 ). The situation changes dramatically

)

, 0

Acc

t

Time (units of

1.0

0.8

0.6

0.4

0.2

0.0

Figure 10. Illustration of pebble accretion in the (relatively) inefficient 3D (left) and efficient 2D (right) regimes, relevant for low and high pebble St , respectively. The relevant pebble capture length scale in the former is the Hill radius, while that in the latter is either the Bondi or Hill radius. See Section 4.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-7bdf2449-2338-4b43-9f55-da502a86d745.jpg)

across the ice-line, beyond which sufficiently massive bodies can undergo Bondi and Hill 2DPA.

Across the ice-line in the wind-dominated disk (Fig. 11b), Earth-sized embryos beyond ≃ 3 AU ( i.e., the iceline) are ensured 2DPA from time zero, and for r ≲ 5 AU, transition from the Bondi to Hill regime in ∼ a cycle of t Acc, 0 ( ≃ 0.5 Myr; Fig. 4b). Mars-sized embryos beyond ≃ 3 AU can partake in Bondi 2DPA from time zero, while Moon-sized embryos do so after ∼ half a cycle of t Acc, 0 . Embryos of mass m p ∼ 10 -5 M ∗ are firmly situated in the Hill regime past ≃ 4 AU from time zero. Conditions for 2DPA are far more stringent in the turbulence-dominated disk (Fig. 11a). Only embryos larger than Earth with m p ≳ 10 -5 M ∗ can partake in 2DPA within the first accretion timescale ( ≃ 0.3 Myr; Fig. 4a), and even so, only at r ≳ 10 AU and in the Bondi regime. The result that embryos within the iceline are confined to 3DPA remains true even after many ( ≳ 5) cycles of t Acc, 0 .

r ), in which viscous heating dominates (see Section 2.2), it comes at no surprise that our key results here corroborate theirs. Namely, rocky embryos within the ice-line are confined to 3DPA due to considerable vertical dispersion of silicate particles, the growth of which is limited by turbulent fragmentation. Given the inefficiency of 3DPA relative to 2DPA, BM22 interpreted this as indicating the accretion of such embryos is primarily driven by mutual pairwise collisions. In Section 4.3.2 below, we show that this is true only in disks that are not driven by MHD winds ( ψ ≳ 1) and characterized by low ˜ ( α ≲ 10 -3 5 . ). That is, disks that are not quiescent, with α ν ≲ 10 -4 . The lack of 2DPA alone is an insufficient criterion for the dominance of classical growth.

As the disk cools and the ice-line migrates inwards over time, ( m /M p ∗ ) 2 D -3 D at any location in the disk decreases owing to the increase in St . However, the rate of ice-line migration far exceeds that at which ( m /M p ∗ ) 2 D -3 D falls towards ( m /M p ∗ ) Iso . In the turbulence-dominated disk, Earth- and Mars-sized embryos (at the present-day locations of those planets in the SS) never enter the 2D regime, even after the iceline passes 1 AU at ≳ 10 t Acc, 0 (the efficiency of 3DPA increases across the ice-line nonetheless; see below). In the wind-dominated disk, lower ( m /M p ∗ ) 2 D -3 D values enable those same embryos to partake in 2DPA after ≳ 3 cycles of t Acc, 0 . This, however, does not constitute 2DPA within the ice-line-larger, more settled icy particles have simply moved into the feeding zone of those planets (see Section 4.3.1).

Considering the similarity in Σ( r, t ) between our fiducial disks and the steady-state viscous disk of BM22(Figs. 4a &amp; 4b) for the inner disk ( i.e., small

Bearing in mind the fact that particles across the entire bouncing-limited turbulent disk are in tens of microns (Fig. 8c), it is clear that 2DPA (and virtually PA in general) cannot be realized therein. The existence of gas and ice giants in our SS indicate that, at least exterior to the ice-line, the bouncing barrier was insignificant if not wholly absent. This follows from the aforementioned 'timescale conflict' associated with gas/ice giant planet formation: pairwise collisions, at large distances from the host star, are incapable of facilitating the rapid core growth necessary for the accretion of gaseous envelopes prior to disk dissipation ( e.g., Rafikov, 2011) (see Fig. 13b). Cm-sized (or larger) pebbles must have been present in the outer disk to fuel the growth of, say, proto-Jupiter. Related to this point, our results in Fig. 11 hint at an outer SS disk with α ν = ˜ ( α/ ψ +1) between ∼ 10 -4 and 10 -3 , such that giant planet cores in the outer disk could have partaken in 2DPA fairly early in SS history. This range is in broad agreement with α ν values inferred from telescopic [i.e., Atacama Large Millimeter/sub-millimeter Array (ALMA)] measurements of molecular line broadening and disk h /h 0 ratios (Rosotti, 2023 &amp; references therein). Nonetheless, it should be taken with a grain of salt. Pebble accretion

Figure 11. Profiles of the critical planetary-stellar mass ratio for transitioning between the 2D to 3D pebble accretion regime (solid lines) and between 2D Bondi and Hill pebble accretion (faded lines), as well as the pebble isolation mass (dotted lines), over time for the (a) turbulence- and (b) wind-dominated disks. In both scenarios, planetary embryos interior to the ice-line are prohibited from partaking in 2DPA. Embryos beyond the ice-line are granted 2DPA given masses that surpass a threshold value which decreases in both time, and with increased contribution from disk winds. See Section 4.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-d89b54c8-57fd-445d-b862-b5bbb21de150.jpg)

in the 2D regime is not necessarily required for rapid growth in the outer disk, where high St and low h 0 serve to enhance the 3DPA rate (Eq. 28). Indeed, for a disk with intermediate ˜ α ≃ 10 -3 and ψ ≃ 1, the massdoubling timescale for a Mars-sized planetary embryo just beyond the ice-line ( i.e., close to Jupiter's inferred formation region; Kruijer et al., 2017a) is on the order of ∼ few times 0 01 . Myr despite undergoing 3DPA. This indicates the few Myr lifespan of the SS disk is ample time for giant planet core accretion, especially considering PA rates increase with embryo size (see Section 4.3.2) and the potential for accelerated growth by elevated solid densities induced by principal ( e.g., water-ice, CO) sublimation lines (see Section 4.1.3 &amp; discussion on planetesimal rings in Section 4.3.3).

## 4.3. Rocky Planet Accretion

In this section, we explore the competition between pairwise collisions and PA in rocky planet formation across parameter space. Our results culminate in conditions ( i.e., ˜, α ψ ) under which each growth pathway is dominant, and are discussed in the context of superEarth origins and Earth's accretion history. The motivating questions for each subsection are: Is there a region of parameter space in which rocky embryos can partake in 2DPA (Section 4.3.1)? How do the massdoubling timescales from pairwise collisions and PA

compare for rocky embryos (Section 4.3.2)? How do the contributions from the two competing processes to the final planet mass compare, and how do they change under different envisioned models for rocky planet formation (Section 4.3.3)? Finally, what do our results imply for the origin of super-Earths (Section 4.4), and how do they fare with isotopic constraints on the building blocks of our Earth from meteoritics (Section 4.5)?

## 4.3.1. Possibility of 2DPA in the Inner Disk?

We have demonstrated, albeit under simplifying assumptions in place ( e.g., no temperature dependence on v f , constant α parameters; see Section 5), that 2DPA is not attainable inward of the ice-line for ˜ α ≃ 10 -3 . Nonetheless, there may be region of ˜ α -ψ parameter space in which sufficiently massive rocky embryos can partake in 2DPA. The shrunken gap between ( m /M p ∗ ) 2 D -3 D and ( m /M p ∗ ) Iso in the inner disk accompanying the transition from turbulenceto wind dominance (Fig. 11) suggests this region lies at high ψ . Moreover, the region likely resides towards low ˜ as α this translates to lower relative turbulent velocities ( i.e., ∆ v t ), higher characteristic St , and thus lower values of ( m /M p ∗ ) 2 D -3 D . Here, we systematically seek out pairs of ˜ and α ψ that yield ( m /M p ∗ ) 2 D -3 D &lt; m /M ( p ∗ ) Iso at r = 1 AU for different times following disk infall.

At each time step ( i.e., 0, 1, and 3 Myr), we gener-

)

, 0

Acc

t

Time (units of

1.0

0.8

0.6

0.4

0.2

0.0

ate a raster of ( m /M p ∗ ) 2 D -3 D at 1 AU for all values of ˜ α and ψ explored. This is compared to a similar raster of ( m /M p ∗ ) Iso , and the pairs (pixels) that yield ( m /M p ∗ ) 2 D -3 D &lt; m /M ( p ∗ ) Iso are determined. It is important to differentiate the potential to partake in 2DPA within the ice-line from that arising from ice-line migration past 1 AU. This is especially relevant for disks that evolve relatively rapidly due to large ˜ ( α i.e., ≳ 10 -2 5 . ) and thus low t Acc, 0 (scaling as 1 /αc ˜ 2 s,c ; Eq. 10). As a related point, it is meaningful to define a proxy for disk dissipation arising from the E/FUV photoevaporation of disk gas ( e.g., St¨rzer &amp; Hollenbach, 1999; Ercolano o &amp; Pascucci, 2017; Concha-Ram´ ırez et al., 2019). Here, we take the time at which the stellar mass accretion rate ˙ M = 2 πr Σ v r,gas ≃ 10 -10 M ⊙ /yr at 1 AU as heralding the end of the disk's life.

Results from this exercise are shown in Fig. 12. At t = 0 (Fig. 12a) and for high ˜ and low α ψ (top left corner), no solids can condense at and within 1 AU, which lies interior to the silicate sublimation line at T = 1500 K . For most of parameter space, rocky embryos are confined to 3DPA by the pebble isolation mass. Only in the bottom right ( i.e., low ˜ and high α ψ as predicted), beyond the solid black curve, can 2DPA be granted for a sufficiently high mass [( m /M p ∗ ) 2 D -3 D ≳ 10 -6 -10 -5 ]. The region in which 2DPA is accessible increases in size with increasing m /M p ∗ (the dashed red curve for Earth sits above that for Mars), reflecting more ˜ α -ψ pairs for which ( m /M p ∗ ) 2 D -3 D &lt; m /M p ∗ . After 1 Myr (Fig. 12b), the silicate sublimation line has migrated within 1 AU for all parameter space, and a region wherein ( m /M p ∗ ) 2 D -3 D &lt; m /M ( p ∗ ) Iso has developed in the top right ( i.e., high ˜ and high α ψ ). This results from iceline migration past 1 AU, and as such, does not represent the possibility of 2DPA within the ice-line. In fact, for all disks with ˜ α ≳ 10 -2 5 . , the ice-line lies interior to 1 AU by this time. The reason most embryos there are still confined to 3DPA can be understood by comparing Figs. 11a &amp; 11b. Greater turbulence ( i.e., lower ψ ; higher α ν for a given ˜) α results in higher temperatures across the disk due to viscous heating (see Section 2.2), in turn leading to lower St and higher ( m /M p ∗ ) 2 D -3 D . For high ˜ and moderately α low ψ ( ≲ 5), the entire ( m /M p ∗ ) 2 D -3 D profile lies above ( m /M p ∗ ) Iso . The rapidity of ice-line migration for these disks arises not only from their large ˜, but stellar irraα diation of the flared disk, which keeps T (or c s ) higher than it would be at r c, 0 (where t Acc, 0 is defined; see Eq. 10) in its absence (Fig. 4c, d). After 3 Myr (Fig. 12c), disks at the upper right corner of parameter space have dissipated, and the region in which embryos at 1 AU are prohibited from 2DPA has shrunken considerably, owing mostly to ice-line migration. In quiescent disks ( α ν ≲ few times 10 -4 ) at the bottom right, the iceline still sits beyond 1 AU. By now, this region where ( m /M p ∗ ) 2 D -3 D &lt; ( m /M p ∗ ) Iso interior to the ice-line has expanded only slightly due to large t Acc, 0 .

We briefly note that t = 0 in our model is not anchored to SS time zero, conventionally defined by the condensation of Calcium-Aluminum-rich Inclusions (CAIs). Moreover, the definition of t Acc, 0 is rather arbitrary, such that its values should not be given much weight [especially for high ˜, where α t Acc, 0 is extremely low (on the order of ∼ 0 01 Myr)]. . It mainly serves to show how the morphology of parameter space evolves with time. Clearly, so long as embryos in even moderately turbulent disks sit within the ice-line ( i.e., accrete 'dry'), they are confined to 3DPA.

## 4.3.2. Embryo Mass-Doubling Timescales

Our investigation of parameter space indicates that for slow-evolving ( i.e., low ˜) disks driven by MHD winds α ( i.e., high ψ ), rocky embryos can accrete pebbles in the efficient 2D regime. While such embryos are relegated to 3DPA for most ˜ α -ψ pairs, it would be unseemly to assume classical growth by pairwise collisions thus constitutes their primary mode of accretion. Instead, a quantitative statement about the favored pathway for rocky planet formation can be made via consideration of the planetary mass-doubling timescale τ p ≃ m /m p ˙ p . The accretion rates ˙ m p for the different regimes of PA are provided above by Eqs. 27, 28, &amp; 30, while that for pairwise collisional growth, in account of gravitational (Safronov) focusing, is given by (Lissauer, 1993; Kokubo &amp; Ida, 2000)

$$\dot { m } _ { p w } \simeq f _ { d } \Sigma \pi R _ { p } ^ { 2 } \Omega _ { K } ( 1 + \Theta ) / \sqrt { 2 \pi }. \quad ( 3 3 )$$

Here, we assume a Mars-sized embryo with m p ≃ 6 4 . × 10 23 kg, and R p ≃ 3 4 . × 10 6 m. The Safronov number Θ = ( v esc /σ v ) 2 quantifies the degree of gravitational focusing and is generally taken to be of order a few. The embryo escape velocity v esc = √ 2 Gm /R p p and σ v is the local velocity dispersion of embryos participating in mutual collisions.

A raster of τ p can be generated for PA assuming the appropriate regime (dependent on r and t ) for each pair of ˜ and α ψ . We denote values in this raster τ p,PA to differentiate them from those in a similar raster generated assuming pairwise collisional growth, denoted τ p,PW . In our calculation of τ p,PA and τ p,PW , we assume the total solid mass f d Σ in the feeding zone ( r ≃ 1 AU ) is distributed evenly between pebbles and embryos. That is, the fraction of solid mass held in pebbles f Peb = 1 -f Emb ≃ 0 5 (see below). . We use the ratio

Figure 12. Rasters of the critical planetary-to-stellar mass ratio for transitioning from 3DPA to 2DPA at r = 1 AU and after (a) 0 Myr, (b) 1 Myr, and (c) 3 Myr of disk evolution. For most of parameter space, planetary embryos in the inner disk are confined to 3DPA by the pebble isolation mass. Only in wind-dominated disks with low ˜ α can sufficiently massive embryos interior to the ice-line partake in 2DPA. The red dashed and dotted lines denote the Earth-to-Sun and Mars-to-Sun mass ratios, respectively. See Section 4.3.1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-96209f39-8261-491d-8d31-76024ff1926b.jpg)

τ p,PA /τ p,PW as a metric for evaluating the competition between PA and classical growth, neglecting for now the dependence of τ p on time (see Section 4.3.3). Assuming Θ ≃ 5, our results for t = 0 and two locations r = 1 AU and r = 20 AU, corresponding to within and without the ice-line, are shown in Fig. 13.

Interior to the ice-line (Fig. 13a), rocky Mars-sized embryos in disks characterized by relatively low ˜ ( α ≲ 10 -3 5 . ; below the black dashed line) have τ p,PA /τ p,PW &lt; 1, indicating PA constitutes their primary mode of accretion. At higher ˜, PA also dominates growth in disks α driven by MHD winds ( ψ ≳ a few). A key observation is that for most of the pebble-dominated parameter space, PA operates in the 3D regime, suggesting the inaccessibility of 2PDA does not necessarily entail the dominance of classical growth. At ψ ≃ 1 and ˜ α ≃ 10 -3 5 . (keeping in mind our choice of f Peb ≃ 0 5 and . Θ ≃ 5) 3DPA is ∼ 10 times more efficient than growth by pairwise collisions. In accordance with Fig. 12, only in wind-dominated ( ψ ≃ 10) disks of low ˜ ( α ≃ 10 -4 ) are the rocky Mars-sized embryos accreting pebbles in the 2D regime. For a turbulence-dominated disk with ˜ α ≃ 10 -3 , such as that of BM22, pairwise collisions remain ∼ 5 times more efficient than PA.

Beyond the ice-line (Fig. 13b), τ p,PA /τ p,PW &lt; 1 across all of parameter space, indicating embryos therein grow primarily from pebbles. Notably, PA is ∼ 10,000 times more efficient than pairwise collisions. The 2DPA region has expanded owing to greater settling of icy particles, though embryos in turbulence-driven disks at relatively high ˜ are still α confined to 3DPA (see Section 4.3.1). For all disks wherein Mars-sized embryos at r = 20AU (and 1 AU) partake in 2DPA, they do so in the Bondi regime. As can be deduced from inspection of Fig. 11, only embryos more massive than Earth in relatively wind-dominated disks can enter the Hill regime.

In both Figs. 13a &amp; 13b, the lack of a gradient in the 2DPA field reflects the effect of the square root in Eq. 27, which dampens the variations in the characteristic St with ˜ and α ψ . Consideration of τ p,PA in going from the bottom right to top left of parameter space ( i.e., deeper into the 3DPA regime) makes the inefficiency of 3D- relative to 2DPA apparent. Within the ice-line (Fig. 13a), for instance, τ p,PA is roughly consistent at ≃ 0.003 Myr within the 2D regime, and remains so across the 2D3D boundary. Once in the 3D regime, however, τ p,PA increases rapidly, reaching ≳ 5 Myr at intermediate values ˜ α ∼ 10 -3 and ψ ∼ 1, and &gt; 10 Myr farther up.

At this stage, it is worthwhile to consider how our results for the inner disk in Fig. 13a change with f Peb and Θ. The scenario explored so far, with f Peb ≃ 0 5, . is likely 'pebble optimistic.' While the reservoir of peb-

bles in the planetary feeding zone may be constantly replenished by their inward drift from farther disk regions, the emergence of embryos from the collapse of concentrated pebble clouds followed by their subsequent growth necessitates a decline in f Peb (and increase in f d ). This process is inferred to be rapid and extensive in the SS. Indeed, Hf-W chronology of iron meteorites coupled with thermal modeling of planetesimals internally heated by 26 Al decay indicate that planetesimals ≳ 100 km in size accreted within ≃ 0.5 to 1 Myr following CAI condensation (preceding time zero in our model), both within and without the ice-line (Kruijer et al., 2014, 2017a; Hunt et al., 2018; Spitzer et al., 2021). Rapid ( ≲ 1 Myr) planetesimal growth in the inner disk is further supported by dust coagulation models ( e.g., Izidoro et al., 2021a; Morbidelli et al., 2022). By the time Mars-sized embryos are present and ubiquitous in the disk, f Peb has likely fallen more closely to ≃ 0.1 or 0.2, if not lower. As for the Safronov number Θ, we note that even if v esc exceeds σ v by a factor of a few due to say, gas drag, collisional damping, and dynamical friction, Θ increases to ∼ 10.

As f Peb decreases and Θ increases, the pebbledominated region of parameter shrinks towards the bottom right (high ψ and low ˜). α For f Peb ≃ 0 5 . and Θ ≃ 5 (Fig. 13a), ∼ 42% of parameter space (excluding the region where T ≳ 1500 K ) is dominated by classical growth. This fraction increases to ∼ 66% for f Peb ∼ 0.1, and to ∼ 72% with an additional increase in Θ to ∼ 10. Accounting for the fact that the equations for ˙ m p, 2 D and ˙ m p, 3 D are overestimates (see Section 4.2), this suggests the growth of Mars-sized embryos to ∼ twice their mass proceeds primarily via pairwise collisions across parameter space. Time, however, has not been considered.

It is important to realize that as accretion progresses, while the disk is still present (see below), PA grows in its dominance. This can be understood as a result of both the increase of the embryo mass m p and more crucially the St of particles with time. Regarding to the former, ˙ m p, 3 D (Eq. 28) scales as ∼ R 3 Hill ∼ m p while ˙ m pw (Eq. 33) scales as ∼ R 2 p ∼ m 2 / 3 p . As m p grows, then, the ratio ˙ m p, 3 D /m ˙ pw does so too as ∼ m 1 / 3 p . Note that in the less applicable 2D regime, ˙ m p, 2 D,Bondi /m ˙ pw ∼ m -1 / 2 p (Eq. 27) and ˙ m p, 2 D,Hill /m ˙ pw ∼ m 0 p (Eq. 30). As time progresses, pebble St increases owing to the decrease in Σ and thus T . This leads to an increase in ˙ m p, 3 D due to both (i) the increase in the PA impact parameter (or equivalently, the decrease in the pebble settling time; see Section 4.2) (Ormel, 2017) and (ii) the decrease in the scale height of the solid sub-disk h 0 ∼ St -1 / 2 (note that h 0 also decreases with h ∼ √ T ). For rocky embryos at r ≃ 1 AU (Fig. 13a), the pebble-dominated region ( i.e., where τ p,PA /τ p,PW &lt; 1) grows rapidly with time, and slightly with mass. The rapid rise in St is largely due to the short initial accretion timescale t Acc, 0 across parameter space, owing to stellar irradiation at r c, 0 , the parametrization for which is independent of ˜ α and ψ (Eq. 13). Short t Acc, 0 also results in a rapid decline in Σ, which slows accretion. Thus, for embryos much larger than Mars, while τ p,PA /τ p,PW (evaluated at t = 0) may be &gt; 1 for most ˜ α -ψ pairs, this dominance of pairwise collisions across parameter space is not necessarily reflected in their contribution to the final planetary mass. This motivates our work in Section 4.3.3, wherein the contributions from pebbles and embryos throughout the full accretion history of a planet are compared. We emphasize that time zero in our model is physically undefined and may well succeed the onset of planetary accretion, at which point higher Σ and T yielded lower St , translating to a greater dominance of classical growth. Moreover, t Acc, 0 values are not necessarily meaningful (see end of Section 4.3.1), prompting us to consider accretion in steady-state disks in Section 4.3.3 below.

To summarize, PA can constitute the dominant pathway of rocky planetary accretion regardless of whether it operates in the 2D or 3D regime . Pairwise collisions between embryos dominate the initial growth of Mars-sized embryos for most of parameter space, given low f Peb and accounting for the fact that ˙ m p, 2 D and ˙ m p, 3 D are upper limits. This dominance may not apply across the full planetary accretion process given the rapid decline in Σ (and thus T ) with time due to short disk evolution timescales governed by t Acc, 0 . The rapid decline in Σ increases the time required to accrete a given mass, and as more time passes, the rise in pebble St and decline in the disk hydrostatic scale height h increase the efficiency of PA relative to pairwise collisions. The physical cutoff to this process is disk dissipation, beyond which pairwise collisions are the sole contributors to planetary growth . The post-CAI timing of SS disk dissipation is constrained to ∼ 4 Myr from paleomagnetic analyses of angrites and CO chondrites (Wang et al., 2017; Borlina et al., 2022) as well as extinct radionuclide systematics of CAIs and chondrules (Dauphas &amp; Chaussidon, 2011). This is concordant with disk dissipation timescales inferred around T-Tauri stars (few to ∼ 10 Myr; e.g., Podosek &amp; Cassen, 1994; Brice˜o et al., 2001). n

As a prelude to Section 4.3.3, we note that the dominance of PA in the growth of a ∼ 100 km-sized embryo to an Earth- or super-Earth-sized planet necessitates (i) most of the growth to be completed within ∼ few Myr, and (ii) the dominance of PA across virtually the entire growth period.

-2

Ratio between

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-dc9c8c1c-2d5b-4912-a7f1-3d173d4c3e6b.jpg)

parameters

(

DW

/

)

Figure 13. Rasters of Log 10 ( τ p,PA /τ p,PW ) at t = 0 for a Mars-sized embryo at (a) r = 1 AU and (b) r = 20 AU. In (a), values corresponding to disks with low ˜ and to a lesser extent high α ψ are &lt; 1, indicating the dominance of pebble accretion in the mass-doubling growth of the embryo. In (b), values are &lt; 1 across all parameter space. See Section 4.3.2.

4.3.3. Building a Super-Earth: Pebble vs. Planetesimal Contributions

Time is clearly a crucial parameter in evaluating the competition between classical growth and PA. In our model, it determines the density of solids (and hence the rapidity of planetary accretion) as well as the pebble St (and hence the efficiency of PA) at any location in the disk. Furthermore, it enables the account of disk dissipation, and thus the cessation of PA, in the planetary growth process. To include the impact of time in our analysis, we simulate the growth of a ∼ 500 kmsized embryo towards ∼ 8 M ⊕ as inferred for most 'superEarths' (Wu, 2019) (see Section 4.4), truncating growth at the pebble isolation mass which decreases over time ( i.e., through the hydrostatic scale height h ; see Eq. 32). In doing so, we keep track of the contribution of pebbles and other embryos to the final planetary mass, denoting them by M PA and M PW , respectively. We use time steps of ∆ t = 0 05 Myr, adding the contributions . ˙ m p,PA ( )∆ t t and ˙ m p,pw ( )∆ t t to the embryo mass at each step. The characteristic size of planetesimals formed via streaming instabilities is estimated to be on the order of 100 km from both theory and numerical simulations (Youdin &amp; Goodman, 2005; Sch¨ afer et al., 2017; Klahr &amp; Schreiber, 2020; Gerbig &amp; Li, 2023). Thus, a starting size of 500 km is within reason. We emphasize that our ensuing analyses, while capable of providing insight into the key factors influencing planetary growth, are fundamentally approximate and do not constitute a replacement for N-body simulations.

Given that our embryo is a 'so-called' first-generation planetesimal, t = 0 here is presumably close to (say, ≲ 1 Myr) the timing of CAI formation. While it is instructive to simulate planetary growth in an evolving disk, the uncertainty in t Acc, 0 values (setting the clock for disk evolution; see Eqs. 7 &amp; 8), beckon a consideration of accretion assuming steady state disks as well. Similar to the exercise in Section 4.3.2, we generate a raster of M PA /M PW for all pairs of ˜ and α ψ at r = 1 AU. We explore how M PA /M PW vary across parameter space, assuming both time evolution and steady state. Note that while the pebble St remains constant in a steady state disk, accretion rates from both pairwise collisions and PA still account for the increasing planetary mass m p . We first consider, as done in Section 4.3.2, disks with smooth solid densities. That is, disks without substructures of concentrated solids ( e.g., dust traps) induced by, say, pressure maxima. In these 'smooth' disks, we set f Peb = 0 1, as expected to be more representative of . the full accretion history compared to the conservative value of 0.5 used in Section 4.3.2. The Safronov number Θ is kept at 5. The dust/solid-to-gas ratio f d 'seen' by the growing embryo ( i.e., excluding its own mass) at 1 AU is assumed to remain ∼ constant throughout accretion, being replenished from the inward drift of solids. Note that the solid density f d Σ applies only before disk dissipation. Following the photoevaporation of disk gas, while f d Σ falls to zero, embryos can persist and proceed to partake in pairwise collisions for tens of Myr. Moreover, f Emb d f Σ underestimates the solid mass held in embryos over time. As embryos proceed to form by the collapse of dust clouds and mutual collisions, f Emb will

approach unity even for a constant pebble abundance, and f d will naturally grow even in the absence of disk substructures.

The raster described for disks evolving in time and in steady state are shown in Figs. 14a &amp; 14b, respectively. Recall from Section 4.3.2 that disk dissipation timescales are on the order of a few Myr following planetesimal formation. That said, ˜ α -ψ pairs which yield planetary accretion times exceeding 10 Myr (to be conservative) are discarded. In the time-evolving case (Fig. 14a), these disks severely overestimate the pebble contribution to the final mass, owing to the increase in St with time (see Section 4.3.2). Within the discarded region, planetary masses at 10 Myr range from ∼ few times 10 -4 to 2 M ⊕ from high ˜ (and to a lesser extent low α ψ ) to low ˜ with most below 1 α M ⊕ . This slow rate of accretion reflects the rapid fall in Σ over time (due to short t Acc, 0 ), leading to small ˙ m pw (Eq. 33) and ˙ m p, 3 D (Eq. 28). Corresponding ratios of M PA /M PW evaluated at t = 10 Myr range from ∼ 0 1 to 120. . The pebble isolation mass across the entire discarded region is greater than 2 M ⊕ . Thus, for most disks wherein the accretion time exceeds disk dissipation, growth is far from completion, and with high ˜, dominated by classical growth α thus far. This suggests that if rocky Mars-sized and larger planets reside in such disks, they grew primarily by pairwise collisions persisting after disk dissipation. For all ˜ α -ψ pairs for which accretion takes place in &lt; 10 Myr, PA dominates planetary growth. These disks are, as expected, those characterized by relatively low ˜ α and driven by MHD winds. The final planetary masses reached in these disks are their respective pebble isolation masses at 1 AU, ranging from ∼ 2 to 7 M ⊕ with increasing ψ . While the initial temperature at 1 AU is higher with greater turbulence, translating to a higher pebble isolation mass, t Acc, 0 and the (3D) PA accretion rate decrease with turbulence. This means that the pebble isolation mass falls faster towards a slower growing planet in turbulence-driven disks, explaining the smaller final planetary masses in them relative to wind-driven disks. As our work in Section 4.3.1 illustrates (Fig. 12), the region of parameter space within which 2DPA is accessible expands towards low ψ with time, owing mostly to the drift of the ice-line past 1 AU. This phenomenon is exhibited in Fig. 14a by the higher of the two regions wherein 2DPA took place.

The leading order observations from Fig. 14a are reprised in the steady state case (Fig. 14b). Here, the inward drift of the ice-line is prohibited such that only embryos in wind-dominated disks at low ˜ can partake α in 2DPA (Fig. 12). Without the aid of time evolution (or more specifically, the accompanying increase in St ), more disks are characterized by accretion times exceeding 10 Myr. Across the PA dominated region at low ˜ α and high ψ , M PA /M PW values are comparable to those in the time-evolving case, owing to relatively large t Acc, 0 . In other words, these disks evolve slowly so the removal of time evolution hardly impacts embryo growth within them. The final planetary masses in these disks range from ∼ 6 to 8 M ⊕ . In the discarded region, planetary masses at 10 Myr increase from high ˜ and low α ψ towards the boundary with the non-discarded region, going from ∼ 10 -3 to 8 M ⊕ . Corresponding ratios of M PA /M PW at 10 Myr range from ∼ 10 -3 to 10. As in the time-evolving case, it is safe to infer that rocky Mars-sized and larger embryos in such disks grow primarily through pairwise collisions.

The paradigm of planet formation in disks with 'smooth' solid densities have recently been challenged by models suggesting rapid embryo growth in planetesimal rings arising at pressure maxima in protoplanetary disks, notably associated with principal ( i.e., silicate and water) sublimation lines (Cai et al., 2022; Izidoro et al., 2022; Morbidelli et al., 2022; Batygin &amp; Morbidelli, 2023). These models are buttressed by observations of ringed structures in T-Tauri disks by ALMA ( e.g., Andrews et al., 2016; Dullemond et al., 2018). As an alternative to accretion in 'smooth' disks, we simulate the growth of the 500 km-sized embryo to super-Earth mass in this emerging framework, once again with and without time evolution. In our model, the planetesimal ring that develops around 1 AU constitutes an annulus with a radial width of ∼ 0 3 . AU hosting ∼ 2 5 . M ⊕ . This translates to a solid density (held in embryos) of ∼ 500 kg/m 3 as assumed by Batygin &amp; Morbidelli (2023). Within such a ring, inelastic collisions coupled with aerodynamic drag dampens the embryo velocity dispersion, enhancing gravitational focusing. This motivates us to adopt a slightly higher Safronov number of Θ = 10. The pebble flux into the ring is taken to be ˙ M peb ∼ 100 M ⊕ /Myr (Bitsch et al., 2019; Lenz et al., 2019), such that the solid (pebble) density Σ peb ≃ ˙ M peb / 2 πrv r , with r = 1 AU and v r evaluated there accordingly. Results from our simulations are shown in Figs. 14c &amp; 14d.

In the time-evolving case (Figs. 14c), growth is completed within 10 Myr across parameter space. The dominance of PA ( i.e., M PA /M PW &gt; 1) in all but turbulence-driven disks at high and moderately low ˜ α ( i.e., &gt; 10 -2 5 . and &lt; 10 -3 ) is attributed to the efficient 2D regime, the widespread accessibility to which results from rapid ice-line migration past 1 AU for high ˜ ( α ≳ 10 -3 ) (Fig. 12). Recall from Section 4.3.1 that for high ˜, the ice-line sits within 1 AU after only 1 Myr has α

Figure 14. Rasters of Log 10 ( M PA /M PW ) for the growth of a 500 km-sized rocky embryo towards 8 M ⊕ at r = 1 AU, assuming accretion in (a) 'smooth', time-evolving disks, (b) 'smooth', steady state disks, (c) a planetesimal ring within time-evolving disks, and (d) a planetesimal ring within steady state disks. In disks with 'smooth' solid densities (a &amp; b), PA dominates growth for high ψ and low ˜ ( α i.e., bottom right of parameter space). Pairwise planetesimal collisions dominate growth at lower ψ and higher ˜, as PA ceases with disk dissipation. α In planetesimal rings (c &amp; d), it appears a substantial region of parameter space is still PA-dominated. However, M PA /M PW in this region is at most ∼ 3. Given the overestimation of PA rates, the omission of 'pebble filtering' effects, etc. (see text for full discussion), pairwise collisions likely dominate growth across virtually all parameter space. See Section 4.3.3.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-1f651cf4-1b63-441b-ac0c-34c6123fd57a.jpg)

passed. Nonetheless, for low ψ ( i.e., substantial turbulence), the transitory stellar-planet mass ratio between the 3D- and 2DPA ( i.e., ( m /M p ∗ ) 2 D -3 D ) sits above the pebble isolation mass ratio ( m /M p ∗ ) Iso . It is due to this confinement in 3DPA that enables the domination of pairwise collisions. Growth to the terminal mass of 8 M ⊕ can still occur in &lt; 10 Myr in these disks owing to both enhanced collisional growth in the planetesimal ring and the sharp increase in St accompanying the introduction of icy pebbles past the ice-line. Keeping in mind inferences on the volatile-poor nature of super-

Earths ( e.g., Owen &amp; Wu, 2017; Gupta &amp; Schlichting, 2019; see Section 4.4), it is unclear if the substantial accretion of icy particles in high ˜ disks as suggested here α (and in Fig. 4a) is meaningful, motivating consideration of the steady state case in which ice-line migration is prohibited. Before moving on, we note that while PA prevails over classical growth across a broad swath of parameter space, it does so by only a factor of a few. Considering the overestimated PA accretion rates and probable reductions in the pebble flux over time due to say, 'pebble filtering' by planets forming at farther

radii ( e.g., Izidoro et al., 2021a), the significance of PA is dubious. It is likely that in most disks (˜ α -ψ pairs), rocky embryos in the planetesimal ring grew primarily via pairwise collisions. Only in wind-dominated disks could the contribution of pebbles and other embryos to growth be comparable.

Interestingly, the final planetary masses ( i.e., pebble isolation masses) in the planetesimal ring model decrease with increasing ψ , a trend opposite of that in our model of a 'smooth' disk solid density (see discussion above for Fig. 14a). This reflects the faster growth rates in planetesimal rings, which enable embryos in turbulence-driven disks to exceed pebble isolation masses in wind-driven disks before pebble isolation masses in their own disks truncate their growth.

In the steady state case (Fig. 14d), the prohibition of ice-line migration changes the morphology of parameter space dramatically. Here, only rocky embryos in wind-dominated disks at low ˜ can access 2DPA, and α grow primarily via PA. As in the time-evolving case, the cumulative contribution from pebbles surpasses that from collisions with other embryos by only a small factor ( &lt; 3). Given the aforementioned advantages conferred on PA in our model, it is safe to conclude that pairwise collisions contribute a comparable if not greater amount of mass to planetary growth, constituting the primary mode of accretion even in such quiescent disks. Across the rest of parameter space, rocky embryos in planetesimal rings owe virtually all their growth to pairwise collisions. It is worth noting that in for high ˜ and α low ψ ( i.e., the discarded region), the planetary mass at 10 Myr far exceeds 1 M ⊕ , ranging from ∼ 6 7 . M ⊕ to just short of 8 M ⊕ . Thus, rapid classical growth is facilitated across all parameter space by planetesimal rings (see Section 4.4).

In reality, the evolution of most disks likely falls somewhere between the steady state and (rapidly) timeevolving cases explored. The 10 Myr cutoff to PA is likely an overestimate, given that disks are inferred to dissipate only a few Myr following planetesimal formation. The application of our model to the SS, for instance, requires implementation of a dissipation time &lt; 4 Myr for consistency with paleomagnetic meteorite analyses (see Section 4.3.2). Dissipation timescales &lt; 10 Myr would lend support to the dominance of classical growth, as less of the final planetary mass can be attributed to PA. Additional support for classical growth, to summarize key points from above, include (i) the overestimation of PA rates, (ii) the negligence of 'pebble filtering' effects from planet formation beyond 1 AU, which would decrease the pebble flux at 1 AU over time, (iii) the fact that the density of solids stored in embryos ought to increase with time (not decrease as f Emb d f Σ does in our model, since f Emb and f d are kept constant), and finally, (iv) the prevalence of disk substructures such as concentric rings and gaps in a host of spatially resolved observations of protoplanetary disks, suggestive of dust trapping and planet formation in planetesimal rings. Even in the absence of rings, (iii) would accelerate accretion by pairwise collisions, which is clearly underestimated in high ˜ disks of Figs. 14a &amp; 14b. α All things considered, it appears the dominance of PA in rocky planet formation can only be realized in disks (i) driven by MHD winds ( i.e., ψ ≳ 1), (ii) undergoing low torques ( i.e., ˜ α ≲ 10 -3 5 . ; slow-evolving), and (iii) absent of pressure maxima that evolve rapidly into planetesimal rings within which classical growth is strongly enhanced. Criteria (i) and (ii) imply disks with α ν ≲ 10 -4 , significantly lower than the range inferred for protoplanetary disks by astronomical observations (Rosotti, 2023), which probe α ν far beyond the ice-line ( α ν is likely higher closer in to the star; see Section 5). This further supports the notion that rocky planets grow primarily by pairwise planetesimal collisions.

## 4.4. Origin of Rocky Super-Earths in Planetesimal Rings

Exoplanet surveys from the Kepler space telescope have revealed the ubiquity of low mass (several M ⊕ ) planets orbiting close (¡¡ 1 AU) to their host stars (Petigura et al., 2013; Mulders et al., 2018), and that such planets exhibit a so-called 'radius valley' in their population between ∼ 1 5 to 2 . R ⊕ (Fulton et al., 2017). It is conventionally accepted that the smaller 'super-Earths' ( R ≲ 1 5 . R ⊕ ) originated as the larger 'sub-Neptunes' ( R ≳ 2 R ⊕ ), having lost their primordial H 2 atmospheres (Bean et al., 2021). The provenance of these planets in their protoplanetary disks remains enigmatic, largely reflecting the uncertainty in their bulk compositions. The two prevailing models for how sub-Neptunes evolve into super-Earths-photoevaporation ( e.g., Owen et al., 2013; Owen &amp; Wu, 2017) and core-powered mass loss (Gupta &amp; Schlichting, 2019)-constrain sub-Neptune core compositions to be volatile-poor. Along with predictions that super-Earths originating in the outer disk will be enriched in ices (Izidoro et al., 2021b), these constraints point to their accretion within the ice-line. This picture of 'dry' accretion has been challenged by increasingly detailed analyses of Kepler planet size distributions, suggesting a subset of super-Earths may in fact be water-rich ( e.g., Teske et al., 2018). Adding to the complexity, initially water-rich super-Earths ( ≲ 3 M ⊕ ) can lose virtually all their water ice to photoevaporation close to their host stars ( ≲ 0 03 AU) (Kurosaki et .

al., 2014). Moreover, accretion of rocky pebbles during the inward migration of icy cores can result in final, largely volatile-poor compositions (Bitsch et al., 2019).

Although it is unclear whether super-Earths ( i.e., sub-Neptune cores) formed interior to the ice-line ( e.g., Chatterjee &amp; Tan, 2013; Ida &amp; Lin, 2010) or across a broad swath of the disk ( e.g., Coleman &amp; Nelson, 2016), there is general consensus that they did not form in-situ ( i.e., at their present-day locations) (Inamdar &amp; Schlichting, 2015; Ogihara et al., 2015). This sets a timescale constraint for rocky super-Earth accretion in the inner disk, as growth must be sufficiently rapid as to allow for disk-driven (Type-I) inward migration prior to disk dissipation. While such rapidity can be achieved by volatile-rich cores undergoing PA in the outer disk, it cannot, for most of parameter space, be realized by rocky cores in the 3DPA-confined inner disk assuming a 'smooth' solid density (see Section 4.3.3). Indeed, only in relatively slow-evolving and wind-dominated 'smooth' disks (Figs. 14a &amp; 14b) may embryos grow sufficiently fast to drift inwards before disk dissipation (recall that disk dissipation timescales are a few Myr, not the conservative 10 Myr assumed). In both the time-evolving and steady state cases (Figs. 14a, b), only 500 km-sized rocky embryos in disks with ˜ α ≲ 10 -3 and ψ ≳ 1 reach the size of Mars in &lt; 4 Myr. While PA in such wind-dominated disks can facilitate rapid growth, we stress that this assumes pebbles are in constant supply. Furthermore, as discussed in Section 4.5, the terrestrial isotopic composition does not support a considerable contribution of pebbles sourced from the outer SS to the Earth's accretion.

Planetesimal rings offer a more compelling scenario for rapid planetary growth. In the time-evolving (Fig. 14c) and steady state (Fig. 14d) cases, rocky embryos grow to Mars size in &lt; . 1 5 Myr and &lt; 2 5 Myr across parameter . space, respectively. In the model for super-Earth formation recently posited by Batygin &amp; Morbidelli (2023), accretion by mutual planetesimal collisions unfolds in a planetesimal ring centered on the silicate sublimation line. The ring arises from rocky grains that condense from silicate vapor diffusing outward across the said line early in the disk's life, which participate in gravitohydrodynamic instabilities to form planetesimals. While the ring may instead owe its origin to the inward drift and pile-up of rocky pebbles at pressure bumps ( e.g., Flock et al., 2019), the exact mechanism by which it forms is hardly as important as the benefits it reaps. Not only is the planetesimal ring advantageous from a timescale perspective, it lends itself nicely to explaining the galactic abundance and intra-system uniformity of super-Earths. Regarding the latter point, growth in the ring proceeds until planetary masses are sufficiently large for removal by Type-I migration. In this way, the terminal masses of super-Earths 'manufactured' in a planetesimal ring, within a given disk, are regulated. In the framework of Batygin &amp; Morbidelli (2023), the dearth of super-Earths in our SS is simply a consequence of excessive turbulence, which prevented the construction of a ring sufficiently massive to induce growth in time for disk-driven inward migration.

To conclude, the emerging paradigm of planetesimal rings is supported by both theastronomically observed ubiquity of disk cavities and dust traps, and its ability to explain the statistics of a substantial exoplanet population. If super-Earths formed 'dry', planetesimal rings offer an avenue for accelerated growth in the inner disk, allowing for disk-driven migration close to their host stars before disk dissipation.

## 4.5. Earth's Accretion from Planetesimal Collisions

In Section 4.3.3, we concluded that PA can dominate rocky planet formation only in quiescent disks characterized by 'smooth' solid densities. The prevalence of classical growth suggested by this result, from a statistical and observational standpoint, has two key implications for Earth's accretion. First, it argues against the rapid accretion of the proto-Earth ( i.e., pre-Moon-forming Giant Impact; MGI) by pebbles, as suggested previously by several authors ( e.g., Johansen et al., 2015; Schiller et al., 2020). It does not, however, preclude the possibility that the proto-Earth accreted rapidly by pairwise planetesimal collisions before the dissipation of the SS disk in a few Myr. As illustrated in Figs. 14c &amp; 14d, and discussed above, classical growth can be greatly accelerated in planetesimal rings, and dominates rocky embryo growth across a broad span of parameter space. If the proto-Earth accreted within such a ring, growth must still have been slow enough such that the SS disk dissipated before substantial inward migration occurred. Interestingly, the relatively low abundance of 182 W in the present-day mantle indicates a series of large impacts ( i.e., planetesimal collisions) necessarily accompanies rapid early accretion. This reflects the fact that an early ( i.e., pre-dissipation) proto-Earth mantle would be enriched in 182 Wtoday from the decay of 182 Hf, lest large impacts facilitate extensive core-mantle equilibration (Olson &amp; Sharp, 2023).

A significant contribution to Earth's accretion from pairwise collisions is indirectly supported by the MGI, which may be considered the final major step in forming Earth. The timing inferred for the MGI ( ≳ 50 Myr after CAIs) from Hf-W isotope systematics of lunar materials (Barboni et al., 2017; Kruijer et al., 2017b; Thiemens et

al., 2021) postdates disk dissipation by tens of Myr, and is consistent with a prolonged accretion of the protoEarth.

Another implication of our work is that the protoEarth did not incorporate a substantial fraction of outer SS ( i.e., carbonaceous; CC) material in the form of pebbles. This is in agreement with analyses of nucleosynthetic isotope anomalies in multi-element space ( e.g., Burkhardt et al., 2021) that suggests the Earth accreted ≲ 10% of CC material by mass, assuming enstatite chondrites (which bear the closest isotopic resemblance to the Earth; e.g., Dauphas, 2017), constitute the Earth's non-carbonaceous (NC) component. Chemical considerations also support an inner disk, predominantly NC origin of the Earth ( e.g., Sossi et al., 2022; Tissot et al., 2022; Liu et al., 2023). Notably, our results contradict the model for Earth's accretion recently proposed by Onyett et al. (2023) on the basis of Si isotope anomalies. This model considers only Si anomalies, and assumes the primary building blocks of Earth are derived from a mixture of NC achondrite material and CI (carbonaceous Ivuna) material. As a consequence, it suggests that ≳ 25% of the proto-Earth mass is attributed to CIlike pebbles from the outer SS. The assumption of NC achondrite-CI mixing is founded upon an inferred temporal evolution in the isotopic composition of inner disk materials from the introduction of CI-like pebbles, itself based on the achondrite-chondrite dichotomy observed in Si anomalies. This dichotomy contrasts with the conventional NC-CC dichotomy observed for all other elements ( e.g., Burkhardt et al., 2019; Kleine et al., 2020) and is possibly due to the inadequate correction for natural mass dependent fractionation with the exponential law, especially considering the relatively small dispersion of Si anomalies across SS materials (Dauphas et al., 2023). We emphasize that pebbles, given a prolonged accretion of Earth, are not expected to contribute significantly to rocky planet formation (Fig. 14).

The heliocentric, s-process trend observed for NC meteorites in multi-element space (Burkhardt et al., 2021) point to a missing, isotopically Earth-like endmember in the innermost disk that has yet to be isolated. Recently, Yap &amp; Tissot (2023) suggested that this missing component may represent the bulk composition of the SS parent molecular cloud. Regardless of its true identity, the inclusion of this component to mixing models for the Earth serves to diminish the fraction of CC material required to explain the terrestrial composition.

## 5. IMPLICATIONS OF OMITTED DETAILS

Having dedicated much of the last section to rocky planet formation, we now return to the disk-wide pic-

Figure 15. Profile of the fragmentation threshold velocity as a function of temperature, having implemented the increase in surface energy at high temperatures due to extensive drying ( i.e., the evaporation of the monolayer of water on grain surfaces/interstices. At T ≃ 1300 K v , f,rock rises to ≃ 9 m/s, close to v f,ice ≃ 10 m/s. See Section 5.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00159v1-figure-fb9463f7-722c-40e1-8ce3-5a38301a12f0.jpg)

ture of the characteristic St and its implications, as addressed in Sections 3, 4.1, &amp; 4.2. Specifically, we discuss here the potential implications of several significant assumptions made in our disk and dust-gas coupling model. As we have shown, fragmentation is the dominant process limiting growth in the absence of the bouncing barrier, resulting in the dramatic rise in St across the ice-line. This dichotomy in St translates to that in characteristic particle sizes, radial drift velocities, and degrees of vertical settling, the latter controlling the dominant mode of planetary accretion. This result reflects the fragmentation threshold velocities v f we have prescribed to silicate and icy particles within and without the ice-line, respectively. That said, it is worthwhile if not proper to consider the temperature dependence of v f , neglected above for simplicity. In particular, v f,rock is expected to increase towards the silicate sublimation line at T ≃ 1500 K , while v f,ice is expected to decrease at temperatures lower than T ∼ 200 K . Indeed, Pillich et al. (2021) found that at temperatures exceeding ∼ 1200 K , the evaporation of water at the interstices and surfaces of silicate dust grains (coupled with phase changes such as the reduction of Fe in silicates to metallic Fe; Pillich et al., 2023) leads to an increase in surface energies, and thus adhesion, by orders of magnitude. Sintering may also contribute to higher v f,rock (Poppe, 2003). Musiolik &amp; Wurm (2019) found that the surface energy of water ice falls by a factor of ≃ 60 simply in going from T &gt; 200 K to T ≃ 175 K , suggesting that icy particles a few AU or less beyond the ice-line will not have an advantage over silicate particles in sticking.

Let us first consider the expected increase in v f,rock at high temperatures, using dimensional analysis as a guide to a quantitative determination of how v f,rock may vary towards the disk inner edge. From Pillich et al. (2021), the effective surface energy γ e as a function of temperature is modeled as such:

$$\gamma _ { e } = \gamma _ { e _ { w } } ( f _ { w } e ^ { \frac { 3 0 0 - T } { T _ { w } } } + f _ { s d } e ^ { \frac { T - 3 0 0 } { T _ { s d } } } ), \quad ( 3 4 )$$

where 300 K ≲ T ≲ 1300 K , and the normalizing constants γ e w ≃ 0 07 . J/m , 2 f w ≃ 8, f sd ≃ 2 × 10 -4 , T w ≃ 950 06 . K , and T sd ≃ 67 16 . K . The functional form of γ e sets that of v f ( T &gt; 300 K ), and we relate the two quantities by setting the specific energy for catastrophic disruption (Benz &amp; Asphaug, 1999; Leinhardt &amp; Stewart, 2009) Q = m v / m 1 2 f 2 2 equal to Cγ e where C is a constant of normalization with units of m /kg (such 2 that Cγ e has units of velocity squared):

$$v _ { f } ( T > 3 0 0 K ) = \sqrt { 2 C \gamma _ { e } } = C ^ { \prime } \sqrt { \gamma _ { e } }, \quad ( 3 5 ) \quad \stackrel { \natural } { w }$$

where we have assumed that the colliding particles are of equal mass and grouped the constant √ 2 C into another constant C ′ , which will be used to to normalize v f at T ≃ 300 K to v f,rock ( T = 300 K ) = 1 m/s. The profile of v f ( T ) which results from implementing Eq. 35 is shown in Fig. 15. At T ≃ 1300 K v , f,rock has steeply risen to ≃ 8 5 m/s, higher than the value of . v f,ice assumed. This rise in v f,rock appears to suggest that there are conditions ( i.e., low ˜ and high α ψ ) for which the dry silicate particle in the inner disk can reach relatively high St , settle considerably, and consequently partake in 2DPA. Nonetheless, the image of a thin sub-disk of dust at the disk inner edge is at odds with intuition, namely the notion the inner region ought to be extremely turbulent due to high ionization fractions for which the MRI operates close to the ideal MHD limit (Armitage, 2020). Here, we call attention to another key assumption we have made in our analysis: that ˜ (including both α α ν and α DW themselves) are constant in space and time.

To first-order, ˜ likely peaks close to the inner edge α of the disk, and falls off with distance from the star. This is true for disks dominated by MRI-driven turbulence ( e.g., Penna et al., 2013) and disk winds so long as the strength of the magnetic field threading the disk decreases outward as well. Moreover, measurements of IR molecular line broadening suggests the inner edge of disks are characterized by α ν ≳ 10 -2 ( e.g., Carr et al., 2004), far higher than values estimated by other techniques ( e.g., sub-mm line broadening, dust scale height) for tens or hundreds of AU ( i.e., α ν ∼ 10 -4 to 10 -3 ) (Rosotti, 2023). As ˜ α rises for a given ψ , the characteristic St as set by turbulent fragmentation and fragmentation from relative drift velocities will decrease through both ˜ α explicitly and the isothermal sound speed c s ∼ T 1 / 2 (Eqs. 17 &amp; 19). High ˜( α ≳ 10 -2 ) likely outweighs the increase in v f,rock in the innermost regions of the disk such that 2DPA is only attainable after several units of t Acc, 0 have passed (˜ may also decay with α time; see below), at which point the disk has probably dissipated.

The counteraction of effects from ˜ and α v f may also apply to the outer disk. With ˜ kept constant, the deα crease in v f,ice with T ≲ 200 K implies that the characteristic St of icy particles (or where T ≲ 150 K ; ≳ 10 AU for the turbulence-dominated disk, and even closer for the wind-dominated disk) will plunge towards values close to those of their silicate counterparts in the inner disk, dampening 2DPA therein. In the event that ˜ decays with distance, it would serve to increase α St . If v f,ice ≃ v f,rock for most of the disk at low temperatures ( ≲ 150K), the dramatic rise in St across the ice-line would be limited in extent (resembling a bump if you will). With the continual decrease of ˜ beyond the iceα line, the outer disk will still host particles larger and more readily settled than the inner disk.

Having considered the spatial variation in ˜, α let us consider briefly the implications of its temporal variation. The decay of α ν (or ˜ α for low ψ ) early in the disk's life is expected as infall onto the disk and gravitational instabilities dwindle ( e.g., Morbidelli et al., 2022). Following the short-lived infall stage, the evolution of α ν will largely be dictated by temperature, which controls the extent of ionization (of elements such as K) and thus the vigor of the MRI (Armitage, 2020). The evolution of α DW , similarly, is governed by the strength of magnetic field lines threading the disk, giving rise to MHD winds. We proceed by assuming ˜ globally decays as α time passes, noting that the ensuing points can simply be reversed for a growing ˜. The decrease in ˜ results in α α an increase in characteristic St across the disk. Graphically, this means that the St profiles in Figs. 5c &amp; 5d will not only evolve to the left towards the star, but also up towards higher values. Over time, then, particles across the disk are more prone to gravitational settling and higher drift rates due to increased headwind drag. Moreover, the threshold to transition from 3DPA to 2DPA decreases more rapidly at any given location in the disk, and PA rates increase. The decay rate of ˜ may scale with ˜ itself, such the effects discussed are α α most applicable to disks with high ( ≳ 10 -3 ) ˜. α

Several, second-order, avenues for improving the model at hand include the implementation of a size dependence for v f (see Section 4.1) and ice-lines for other prominent condensates such as CO and CH . We do not 4

foresee these details changing the key takeaways from this work to any significant degree.

## 6. CONCLUDING POINTS

In this work, we model the degree of dust-gas coupling across protoplanetary disks evolving under both turbulence and MHD disk winds. The key conclusions of our work are the following:

- · Absent of the bouncing barrier, fragmentation is the dominant mechanism limiting particle growth, trumping radial drift in setting the characteristic Stokes number of particles across the disk.
- · Assuming icy particles fragment at higher collisional velocities than silicate particles, the waterice sublimation line constitutes a crucial transition point in particle Stokes number, and consequently, dust settling, drift velocity, and planetary accretion regimes. Namely, higher St translates to a thinner solid sub-disk, faster radial drift velocities, and a lower threshold to enter the 2D pebble accretion regime.
- · For ˜ = 10 α -3 , silicate particles interior to the iceline in turbulence-dominated ( ψ ≃ 0 1) and wind-. dominated ( ψ ≃ 10) disks are characterized by St ≲ 10 -3 and St ≲ 10 -2 , respectively, corresponding to sizes in the mm- to 1-cm-scale. Icy particles beyond the ice-line possess Stokes numbers ∼ two orders of magnitude larger than their silicate counterparts, at St ≳ 10 -2 and St ≳ 0 1, . corresponding to the cm- to dm-scale. Bouncinglimited particles are smaller than fragmentationlimited particles, and assume sizes from the tens to hundreds of microns, depending on disk parameters.
- · For most of parameter space, with the exception of a rather small corner at low ˜ ( α i.e., ≲ 10 -3 5 . ) and moderately high ψ ( i.e., ≳ 1), planetary embryos interior to the ice-line ( i.e., at r = 1 AU) are confined to the 3D regime of pebble accretion. Inaccessibility to 2D pebble accretion in the inner disk, however, does not necessarily imply the dominance of classical growth by pairwise collisions between embryos.
- · For pebble accretion to dominate rocky planet formation, the main stage of growth must be (i) completed within a few million years, prior to disk dissipation, and be (ii) attributed mostly to pebbles. This can only be realized in disks whose evolution is (i) driven by MHD winds ( i.e., high ψ )
- and (ii) relatively slow owing to low torques ( i.e., low ˜), α and (iii) in the absence of pressure maxima that can facilitate the concentration of solids in the planetary feeding zone, resulting in planetesimal rings that significantly accelerate classical growth. Such disks are extremely quiescent, with α ν ≲ 10 -4 . Thus, for almost all of parameter space corresponding to α ν values observed for protoplanetary disks ( ≳ few × 10 -4 ), pairwise collisional growth is predicted to dominate pebble accretion interior to the water ice-line, serving as the primary mode of rocky planet accretion.
- · For planetary embryos exterior to the ice-line ( i.e., at r = 20 AU), pebble accretion dominates classical collisional growth for all parameter space, regardless of whether the planetary embryo lies in the 2D or 3D pebble accretion regime. The rapid accretion of giant planet cores suggests the outer SS disk was not turbulence-dominated, such that 2D pebble accretion was accessible therein. Rapid growth can be facilitated even in the 3D regime with the aid of sublimation fronts that elevate local solid densities.
- · Super-Earth formation in planetesimal rings is supported by the rapidity of rocky planet accretion facilitated by collisional growth within such rings, allowing for disk-driven migration to take place before disk dissipation. The planetesimal ring framework is also favored from an observational standpoint, and for its ability to explain the galactic abundance and intra-system uniformity of super-Earths.
- · The dominance of classical growth argues against a significant ( ≳ 10%) contribution of outer disk, carbonaceous material to the Earth in the form of pebbles, in accordance with chemical and isotopic investigations of Earth's accretion history.

We conclude by emphasizing that growth rates associated with different accretionary processes depend on both the overall (solid and gas) disk structure as well as the Stokes number, and quantitative measures must be undertaken to delineate the details of their competition. Furthermore, insofar as our results on rocky planet accretion are concerned, MHD winds can be replaced by any mechanism aside from turbulent viscosity that exerts a torque on the disk, driving stellar accretion. Seeing as fragmentation is key in setting the characteristic Stokes number throughout the disk, future models

should feature a more refined prescription of fragmentation threshold velocities as constrained by laboratory experiments. A prescription of spatially and/or temporally variable α parameters also represents a worthwhile pursuit.

## DECLARATION OF COMPETING INTEREST

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

## ACKNOWLEDGEMENTS

We thank three anonymous reviewers for their constructive reviews on the manuscript, and editor Sean Raymond for prompt and careful editorial handling. We also thank Fran¸ cois L. H. Tissot for advice on improving clarity in communicating our results. This work was supported by a Packard fellowship and National Science Foundation grant (AST 2109276) to KB, as well as the Caltech Center for Comparative Planetary Evolution (3CPE).

## REFERENCES

Ansdell, M., Williams, J. P., Manara, C. F., Miotello, A., Facchini, S., van der Marel, N., ... &amp; van Dishoeck, E. F. (2017). An ALMA survey of protoplanetary disks in the σ Orionis cluster. The Astronomical Journal , 153 (5), 240.

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., &amp; Dullemond, C. P. (2009). Protoplanetary disk structures in Ophiuchus. The Astrophysical Journal, 700(2), 1502.

Andrews, S. M., Wilner, D. J., Zhu, Z., Birnstiel, T., Carpenter, J. M., P´rez, L. M., ... e &amp; Ricci, L. (2016). Ringed substructure and a gap at 1 au in the nearest protoplanetary disk. The Astrophysical Journal Letters , 820 (2), L40.

Armitage, P. J. (2020). Astrophysics of planet formation. Cambridge University Press.

Balbus, S. A., &amp; Hawley, J. F. (1998). Instability, turbulence, and enhanced transport in accretion disks. Reviews of modern physics, 70(1), 1.

Barboni, M., Boehnke, P., Keller, B., Kohl, I. E., Schoene, B., Young, E. D., &amp; McKeegan, K. D. (2017). Early formation of the Moon 4.51 billion years ago. Science advances , 3 (1), e1602365.

Batygin, K., &amp; Morbidelli, A. (2020). Formation of giant planet satellites. The Astrophysical Journal , 894 (2), 143.

Batygin, K., &amp; Morbidelli, A. (2022). Self-consistent model for dust-gas coupling in protoplanetary disks. Astronomy &amp; Astrophysics, 666, A19.

Batygin, K., &amp; Morbidelli, A. (2023). Formation of rocky super-earths from a narrow ring of planetesimals. Nature Astronomy, 7(3), 330-338.

Bean, J. L., Raymond, S. N., &amp; Owen, J. E. (2021). The Nature and Origins of Sub-Neptune Size Planets. Journal of Geophysical Research: Planets, 126(1), e2020JE006639.

Beitz, E., G¨ttler, C., Blum, J., Meisner, T., Teiser, u J., &amp; Wurm, G. (2011). Low-velocity collisions of centimeter-sized dust aggregates. The Astrophysical Journal, 736(1), 34.

Benz, W., &amp; Asphaug, E. (1999). Catastrophic disruptions revisited. Icarus, 142(1), 5-20.

Birnstiel, T., Ormel, C. W., &amp; Dullemond, C. P. (2011). Dust size distributions in coagulation/fragmentation equilibrium: numerical solutions and analytical fits. Astronomy &amp; Astrophysics, 525, A11.

Birnstiel, T., Klahr, H., &amp; Ercolano, B. (2012). A simple model for the evolution of the dust population in protoplanetary disks. Astronomy &amp; Astrophysics, 539, A148.

Bitsch, B., Morbidelli, A., Lega, E., &amp; Crida, A. (2014). Stellar irradiated discs and implications on migration of embedded planets-II. Accreting-discs. Astronomy &amp; Astrophysics, 564, A135.

Bitsch, B., Izidoro, A., Johansen, A., Raymond, S. N., Morbidelli, A., Lambrechts, M., &amp; Jacobson, S. A. (2019). Formation of planetary systems by pebble accretion and migration: growth of gas giants. Astronomy &amp; Astrophysics , 623 , A88.

Bitsch, B., Raymond, S. N., &amp; Izidoro, A. (2019). Rocky super-Earths or waterworlds: The interplay of planet migration, pebble accretion, and disc evolution. Astronomy &amp; Astrophysics, 624, A109.

Borlina, C. S., Weiss, B. P., Bryson, J. F., &amp; Armitage, P. J. (2022). Lifetime of the outer solar system nebula from carbonaceous chondrites. Journal of Geophysical Research: Planets , 127 (7), e2021JE007139.

Boss, A. P., &amp; Ciesla, F. J. (2014). The solar nebula. Planets, Asteroids, Comets and The Solar System, 2, 37-53.

Blandford, R. D., &amp; Payne, D. G. (1982). Hydromagnetic flows from accretion discs and the production of radio jets. Monthly Notices of the Royal Astronomical Society, 199(4), 883-903.

Blum, J., &amp; M¨nch, M. (1993). u Experimental investigations on aggregate-aggregate collisions in the early solar nebula. Icarus, 106(1), 151-167.

Blum, J., &amp; Wurm, G. (2008). The growth mechanisms of macroscopic bodies in protoplanetary disks. Annu. Rev. Astron. Astrophys. , 46 , 21-56.

Briceno, C., Vivas, A. K., Calvet, N., Hartmann, L., Pacheco, R., Herrera, D., ... &amp; Andrews, P. (2001). The CIDA-QUEST large-scale survey of Orion OB1: Evidence for rapid disk dissipation in a dispersed stellar population. Science, 291(5501), 93-96.

Burkhardt, C., Spitzer, F., Morbidelli, A., Budde, G., Render, J. H., Kruijer, T. S., &amp; Kleine, T. (2021). Terrestrial planet formation from lost inner solar system material. Science advances, 7(52), eabj7601.

Cai, M. X., Tan, J. C., &amp; Portegies Zwart, S. (2022). Inside-out planet formation: VI. oligarchic coagulation of planetesimals from a pebble ring?. Monthly Notices of the Royal Astronomical Society , 510 (4), 5486-5499.

Carr, J. S., Tokunaga, A. T., &amp; Najita, J. (2004). Hot H2O emission and evidence for turbulence in the disk of a young star. The Astrophysical Journal , 603 (1), 213.

Chambers, J. E., &amp; Wetherill, G. W. (1998). Making the terrestrial planets: N-body integrations of planetary embryos in three dimensions. Icarus, 136(2), 304-327.

Chambers, J. E. (2009). An analytic model for the evolution of a viscous, irradiated disk. The Astrophysical Journal , 705 (2), 1206.

Chatterjee, S., &amp; Tan, J. C. (2013). Inside-out planet formation. The Astrophysical Journal, 780(1), 53.

Chiang, E. I., &amp; Goldreich, P. (1997). Spectral energy distributions of T Tauri stars with passive circumstellar disks. The Astrophysical Journal, 490(1), 368.

Concha-Ram´ ırez, F., Wilhelm, M. J., Portegies Zwart, S., &amp; Haworth, T. J. (2019). External photoevaporation of circumstellar discs constrains the time-scale for planet formation. Monthly Notices of the Royal Astronomical Society , 490 (4), 5678-5690.

Cuzzi, J. N., Hogan, R. C., Paque, J. M., &amp; Dobrovolskis, A. R. (2001). Size-selective concentration of chondrules and other small particles in protoplanetary nebula turbulence. The Astrophysical Journal, 546(1), 496.

Dauphas, N., &amp; Chaussidon, M. (2011). A perspective from extinct radionuclides on a young stellar object: the Sun and its accretion disk. Annual Review of Earth and Planetary Sciences , 39 , 351-386.

Dauphas, N. (2017). The isotopic nature of the Earth's accreting material through time. Nature, 541(7638), 521-524.

Dauphas, N., Hopp, T., &amp; Nesvorn´ y, D. (2024). Bayesian inference on the isotopic building blocks of Mars and Earth. Icarus , 408 , 115805.

Dubrulle, B., Morfill, G., &amp; Sterzik, M. (1995). The dust subdisk in the protoplanetary nebula. icarus, 114(2), 237-246.

Dullemond, C. P., Birnstiel, T., Huang, J., Kurtovic, N. T., Andrews, S. M., Guzm´ an, V. V., ... &amp; Ricci, L. (2018). The disk substructures at high angular resolution project (DSHARP). VI. Dust trapping in thin-ringed protoplanetary disks. The Astrophysical Journal Letters , 869 (2), L46.

Ercolano, B., &amp; Pascucci, I. (2017). The dispersal of planet-forming discs: theory confronts observations. Royal Society Open Science , 4 (4), 170114.

Flock, M., Turner, N. J., Mulders, G. D., Hasegawa, Y., Nelson, R. P., &amp; Bitsch, B. (2019). Planet formation and migration near the silicate sublimation front in protoplanetary disks. Astronomy &amp; Astrophysics, 630,

A147.

Friedrich, J. M., Weisberg, M. K., Ebel, D. S., Biltz, A. E., Corbett, B. M., Iotzov, I. V., ... &amp;Wolman, M. D. (2015). Chondrule size and related physical properties: A compilation and evaluation of current data across all meteorite groups. Geochemistry , 75 (4), 419-443.

Fulton, B. J., Petigura, E. A., Howard, A. W., Isaacson, H., Marcy, G. W., Cargile, P. A., ... &amp; Hirsch, L. A. (2017). The California-Kepler survey. III. A gap in the radius distribution of small planets. The Astronomical Journal, 154(3), 109.

Gerbig, K., &amp; Li, R. (2023). Planetesimal Initial Mass Functions Following Diffusion-regulated Gravitational Collapse. The Astrophysical Journal , 949 (2), 81.

Gundlach, B., Kilias, S., Beitz, E., &amp; Blum, J. (2011). Micrometer-sized ice particles for planetary-science experiments-I. Preparation, critical rolling friction force, and specific surface energy. Icarus, 214(2), 717-723.

Gundlach, B., &amp; Blum, J. (2014). The stickiness of micrometer-sized water-ice particles. The Astrophysical Journal , 798 (1), 34.

Gundlach, B., Schmidt, K. P., Kreuzig, C., Bischoff, D., Rezaei, F., Kothe, S., ... &amp; Stoll, E. (2018). The tensile strength of ice and dust aggregates and its dependence on particle properties. Monthly Notices of the Royal Astronomical Society , 479 (1), 1273-1277.

Gupta, A., &amp; Schlichting, H. E. (2019). Sculpting the valley in the radius distribution of small exoplanets as a by-product of planet formation: the core-powered mass-loss mechanism. Monthly Notices of the Royal Astronomical Society, 487(1), 24-33.

G¨ uttler, C., Blum, J., Zsom, A., Ormel, C. W., &amp; Dullemond, C. P. (2010). The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals?I. Mapping the zoo of laboratory collision experiments. Astronomy &amp; Astrophysics, 513, A56.

Hartmann, L., Calvet, N., Gullbring, E., &amp; D'Alessio, P. (1998). Accretion and the evolution of T Tauri disks. The Astrophysical Journal, 495(1), 385.

Hayashi, C. (1981). Structure of the solar nebula, growth and decay of magnetic fields and effects of magnetic and turbulent viscosities on the nebula. Progress

Hunt, A. C., Cook, D. L., Lichtenberg, T., Reger, P. M., Ek, M., Golabek, G. J., &amp; Sch¨nb¨chler, M. (2018). o a Late metal-silicate separation on the IAB parent asteroid: constraints from combined W and Pt isotopes and thermal modelling. Earth and Planetary Science Letters , 482 , 490-500.

Ida, S., Guillot, T., &amp; Morbidelli, A. (2016). The radial dependence of pebble accretion rates: A source of diversity in planetary systems-I. Analytical formulation. Astronomy &amp; Astrophysics, 591, A72.

Ida, S., &amp; Lin, D. N. C. (2010). Toward a deterministic model of planetary formation. VI. Dynamical interaction and coagulation of multiple rocky embryos and super-Earth systems around solar-type stars. The Astrophysical Journal , 719 (1), 810.

Inamdar, N. K., &amp; Schlichting, H. E. (2015). The formation of super-Earths and mini-Neptunes with giant impacts. Monthly Notices of the Royal Astronomical Society , 448 (2), 1751-1760.

Izidoro, A., Bitsch, B., &amp; Dasgupta, R. (2021a). The effect of a strong pressure bump in the Sun's natal disk: terrestrial planet formation via planetesimal accretion rather than pebble accretion. The Astrophysical Journal , 915 (1), 62.

Izidoro, A., Bitsch, B., Raymond, S. N., Johansen, A., Morbidelli, A., Lambrechts, M., &amp; Jacobson, S. A. (2021b). Formation of planetary systems by pebble accretion and migration-Hot super-Earth systems from breaking compact resonant chains. Astronomy &amp; Astrophysics, 650, A152.

Izidoro, A., Dasgupta, R., Raymond, S. N., Deienno, R., Bitsch, B., &amp; Isella, A. (2022). Planetesimal rings as the cause of the Solar System's planetary architecture. Nature Astronomy , 6 (3), 357-366.

Johansen, A., Low, M. M. M., Lacerda, P., &amp; Bizzarro, M. (2015). Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion. Science Advances, 1(3), e1500109.

Kenyon, S. J., &amp; Hartmann, L. W. (1987). Spectral energy distributions of T Tauri stars-Disk flaring and limits on accretion. Astrophysical Journal, Part 1 (ISSN

0004-637X), vol. 323, Dec. 15, 1987, p. 714-733. , 323 , 714-733.

Klahr, H., &amp; Schreiber, A. (2020). Turbulence sets the length scale for planetesimal formation: Local 2D simulations of streaming instability and planetesimal formation. The Astrophysical Journal , 901 (1), 54.

Kleine, T., Budde, G., Burkhardt, C., Kruijer, T. S., Worsham, E. A., Morbidelli, A., &amp; Nimmo, F. (2020). The non-carbonaceous-carbonaceous meteorite dichotomy. Space Science Reviews, 216, 1-27.

Kokubo, E., &amp; Ida, S. (2000). Formation of protoplanets from planetesimals in the solar nebula. Icarus , 143 (1), 15-27.

Kruijer, T. S., Burkhardt, C., Budde, G., &amp; Kleine, T. (2017a). Age of Jupiter inferred from the distinct genetics and formation times of meteorites. Proceedings of the National Academy of Sciences , 114 (26), 6712-6716.

Kruijer, T. S., &amp; Kleine, T. (2017b). Tungsten isotopes and the origin of the Moon. Earth and Planetary Science Letters , 475 , 15-24.

Kruijer, T. S., Touboul, M., Fischer-G¨dde, o M., Bermingham, K. R., Walker, R. J., &amp; Kleine, T. (2014). Protracted core formation and rapid accretion of protoplanets. Science , 344 (6188), 1150-1154.

Kurosaki, K., Ikoma, M., &amp; Hori, Y. (2014). Impact of photo-evaporative mass loss on masses and radii of water-rich sub/super-Earths. Astronomy &amp; Astrophysics , 562 , A80.

Lambrechts, M., &amp; Johansen, A. (2012). Rapid growth of gas-giant cores by pebble accretion. Astronomy &amp; Astrophysics, 544, A32.

Lambrechts, M., Johansen, A., &amp; Morbidelli, A. (2014). Separating gas-giant and ice-giant planets by halting pebble accretion. Astronomy &amp; Astrophysics, 572, A35.

Leinhardt, Z. M., &amp; Stewart, S. T. (2009). Full numerical simulations of catastrophic small body collisions. Icarus, 199(2), 542-559.

Lenz, C. T., Klahr, H., &amp; Birnstiel, T. (2019). Planetesimal population synthesis: Pebble flux-regulated planetesimal formation. The Astrophysical Journal , 874 (1), 36.

Liu, W., Zhang, Y., Tissot, F. L., Avice, G., Ye, Z., &amp; Yin, Q. Z. (2023). I/Pu reveals Earth mainly accreted from volatile-poor differentiated planetesimals. Science Advances, 9(27), eadg9213.

Lissauer, J. J. (1993). Planet formation. Annual review of astronomy and astrophysics , 31 (1), 129-172.

Lynden-Bell, D., &amp; Pringle, J. E. (1974). The evolution of viscous discs and the origin of the nebular variables. Monthly Notices of the Royal Astronomical Society, 168(3), 603-637.

Miotello, A., Kamp, I., Birnstiel, T., Cleeves, L. I., &amp; Kataoka, A. (2022). Setting the stage for planet formation: measurements and implications of the fundamental disk properties. arXiv preprint arXiv:2203.09818 .

Morbidelli, A., Baillie, K., Batygin, K., Charnoz, S., Guillot, T., Rubie, D. C., &amp; Kleine, T. (2022). Contemporary formation of early Solar System planetesimals at two distinct radial locations. Nature Astronomy, 6(1), 72-79.

Mulders, G. D., Pascucci, I., Apai, D., &amp; Ciesla, F. J. (2018). The exoplanet population observation simulator. I. the inner edges of planetary systems. The Astronomical Journal , 156 (1), 24.

Musiolik, G., &amp; Wurm, G. (2019). Contacts of water ice in protoplanetary disks-laboratory experiments. The Astrophysical Journal, 873(1), 58.

Nakagawa, Y., Sekiya, M., &amp; Hayashi, C. (1986). Settling and growth of dust particles in a laminar phase of a low-mass solar nebula. Icarus, 67(3), 375-390.

Nakamoto, T., &amp; Nakagawa, Y. (1994). Formation, early evolution, and gravitational stability of protoplanetary disks. Astrophysical Journal, Part 1 (ISSN 0004-637X), vol. 421, no. 2, p. 640-650, 421, 640-650.

Nelson, R. P., Gressel, O., &amp; Umurhan, O. M. (2013). Linear and non-linear evolution of the vertical shear instability in accretion discs. Monthly Notices of the Royal Astronomical Society, 435(3), 2610-2632.

Ogihara, M., Morbidelli, A., &amp; Guillot, T. (2015). A reassessment of the in situ formation of close-in super-

Earths. Astronomy &amp; Astrophysics, 578, A36.

Olson, P. L., &amp; Sharp, Z. D. (2023). Hafniumtungsten evolution with pebble accretion during Earth formation. Earth and Planetary Science Letters , 622 , 118418.

Onyett, I. J., Schiller, M., Makhatadze, G. V., Deng, Z., Johansen, A., &amp; Bizzarro, M. (2023). Silicon isotope constraints on terrestrial planet accretion. Nature, 1-6.

Ormel, C. W., &amp; Cuzzi, J. N. (2007). Closed-form expressions for particle relative velocities induced by turbulence. Astronomy &amp; Astrophysics, 466(2), 413420.

Ormel, C. W., &amp; Klahr, H. H. (2010). The effect of gas drag on the growth of protoplanets-analytical expressions for the accretion of small bodies in laminar disks. Astronomy &amp; Astrophysics, 520, A43.

Ormel, C. W. (2017). The emerging paradigm of pebble accretion. Formation, Evolution, and Dynamics of Young Solar Systems, 197-228.

Owen, J. E., &amp; Wu, Y. (2013). Kepler planets: a tale of evaporation. The Astrophysical Journal, 775(2), 105.

Owen, J. E., &amp; Wu, Y. (2017). The evaporation valley in the Kepler planets. The Astrophysical Journal, 847(1), 29.

Penna, R. F., Sadowski, A., Kulkarni, A. K., &amp; Narayan, R. (2013). The Shakura-Sunyaev viscosity prescription with variable α (r). Monthly Notices of the Royal Astronomical Society, 428(3), 2255-2274.

Petigura, E. A., Howard, A. W., &amp; Marcy, G. W. (2013). Prevalence of Earth-size planets orbiting Sunlike stars. Proceedings of the National Academy of Sciences, 110(48), 19273-19278.

Pillich, C., Bogdan, T., Landers, J., Wurm, G., &amp; Wende, H. (2021). Drifting inwards in protoplanetary discs-II. The effect of water on sticking properties at increasing temperatures. Astronomy &amp; Astrophysics, 652, A106.

Pillich, C., Bogdan, T., Tasto, J., Landers, J., Wurm, G., &amp; Wende, H. (2023). Composition and Sticking of Hot Chondritic Dust in a Protoplanetary Hydrogen At- mosphere. The Planetary Science Journal , 4 (10), 195.

Pinilla, P., Lenz, C. T., &amp; Stammler, S. M. (2021). Growing and trapping pebbles with fragile collisions of particles in protoplanetary disks. Astronomy &amp; Astrophysics , 645 , A70.

Podosek, F. A., &amp; Cassen, P. (1994). Theoretical, observational, and isotopic estimates of the lifetime of the solar nebula. Meteoritics, 29(1), 6-25.

Pollack, J. B., McKay, C. P., &amp; Christofferson, B. M. (1985). A calculation of the Rosseland mean opacity of dust grains in primordial solar system nebulae. Icarus , 64 (3), 471-492.

Poppe, T. (2003). Sintering of highly porous silicaparticle samples: analogues of early Solar-System aggregates. Icarus , 164 (1), 139-148.

Rafikov, R. R. (2011). Constraint on the giant planet production by core accretion. The Astrophysical Journal, 727(2), 86.

Ricci, L., Testi, L., Natta, A., Neri, R., Cabrit, S., &amp; Herczeg, G. J. (2010). Dust properties of protoplanetary disks in the Taurus-Auriga star forming region from millimeter wavelengths. Astronomy &amp; Astrophysics , 512 , A15.

Rosenthal, L. J., Knutson, H. A., Chachan, Y., Dai, F., Howard, A. W., Fulton, B. J., ... &amp; Wright, J. T. (2022). The California legacy survey. III. On the shoulders of (some) giants: the relationship between inner small planets and outer massive planets. The Astrophysical Journal Supplement Series, 262(1), 1.

Rosotti, G. P. (2023). Empirical constraints on turbulence in proto-planetary discs. New Astronomy Reviews, 101674.

Rozitis, B., Ryan, A. J., Emery, J. P., Christensen, P. R., Hamilton, V. E., Simon, A. A., ... &amp; Lauretta, D. S. (2020). Asteroid (101955) Bennu's weak boulders and thermally anomalous equator. Science Advances , 6 (41), eabc3699.

Safronov, V. S. (1972). Evolution of the Protoplanetary Cloud and Formation of the Earth and the Planets.

Sch¨ afer, U., Yang, C. C., &amp; Johansen, A. (2017). Initial mass function of planetesimals formed by the

streaming instability. Astronomy &amp; Astrophysics , 597 , A69.

Schiller, M., Bizzarro, M., &amp; Siebert, J. (2020). Iron isotope evidence for very rapid accretion and differentiation of the proto-Earth. Science advances, 6(7), eaay7604.

Sears, D. W. (1998). The case for rarity of chondrules and calcium-aluminum-rich inclusions in the early solar system and some implications for astrophysical models. The Astrophysical Journal , 498 (2), 773.

Shakura, N. I., &amp; Sunyaev, R. A. (1973). Black holes in binary systems. Observational appearance. Astronomy and Astrophysics, Vol. 24, p. 337-355, 24, 337-355.

Sossi, P. A., Stotz, I. L., Jacobson, S. A., Morbidelli, A., &amp; O'Neill, H. S. C. (2022). Stochastic accretion of the Earth. Nature astronomy , 6 (8), 951-960.

Spitzer, F., Burkhardt, C., Nimmo, F., &amp; Kleine, T. (2021). Nucleosynthetic Pt isotope anomalies and the Hf-W chronology of core formation in inner and outer solar system planetesimals. Earth and Planetary Science Letters , 576 , 117211.

Squire, J., &amp; Hopkins, P. F. (2018). Resonant drag instabilities in protoplanetary discs: the streaming instability and new, faster growing instabilities. Monthly Notices of the Royal Astronomical Society , 477 (4), 50115040.

Steinpilz, T., Joeris, K., Jungmann, F., Wolf, D., Brendel, L., Teiser, J., ... &amp; Wurm, G. (2020). Electrical charging overcomes the bouncing barrier in planet formation. Nature Physics, 16(2), 225-229.

Stepinski, T. F. (1998). The solar nebula as a process-An analytic model. Icarus , 132 (1), 100-112.

St¨rzer, o H., &amp; Hollenbach, D. (1999). Photodissociation region models of photoevaporating circumstellar disks and application to the proplyds in Orion. The Astrophysical Journal , 515 (2), 669.

Tabone, B., Rosotti, G. P., Cridland, A. J., Armitage, P. J., &amp; Lodato, G. (2022). Secular evolution of MHD wind-driven discs: analytical solutions in the expanded α -framework. Monthly Notices of the Royal Astronomical Society, 512(2), 2290-2309.

Teske, J. K., Ciardi, D. R., Howell, S. B., Hirsch, L. A., &amp; Johnson, R. A. (2018). The effects of stellar companions on the observed transiting exoplanet radius distribution. The Astronomical Journal , 156 (6), 292.

Testi, L., Natta, A., Shepherd, D. S., &amp; Wilner, D. J. (2003). Large grains in the disk of CQ Tau. Astronomy &amp; Astrophysics , 403 (1), 323-328.

Thiemens, M. M., Tusch, J., OC Fonseca, R., Leitzke, F., Fischer-G¨dde, M., Debaille, V., ... o &amp; M¨ unker, C. (2021). Reply to: No 182W evidence for early Moon formation. Nature Geoscience , 14 (10), 716-718.

Tissot, F. L., Collinet, M., Namur, O., &amp; Grove, T. L. (2022). The case for the angrite parent body as the archetypal first-generation planetesimal: Large, reduced and Mg-enriched. Geochimica et Cosmochimica Acta, 338, 278-301.

Trapman, L., Tabone, B., Rosotti, G., &amp; Zhang, K. (2022). Effect of MHD wind-driven disk evolution on the observed sizes of protoplanetary disks. The Astrophysical Journal, 926(1), 61.

Wada, K., Tanaka, H., Suyama, T., Kimura, H., &amp; Yamamoto, T. (2011). The rebound condition of dust aggregates revealed by numerical simulation of their collisions. The Astrophysical Journal, 737(1), 36.

Wang, H., Weiss, B. P., Bai, X. N., Downey, B. G., Wang, J., Wang, J., ... &amp; Zucolotto, M. E. (2017). Lifetime of the solar nebula constrained by meteorite paleomagnetism. Science , 355 (6325), 623-627.

Weidenschilling, S. J. (1977). Aerodynamics of solid bodies in the solar nebula. Monthly Notices of the Royal Astronomical Society, 180(2), 57-70.

Weidenschilling, S. J., &amp; Cuzzi, J. N. (1993). Formation of planetesimals in the solar nebula. In Protostars and planets III (pp. 1031-1060).

Weidling, R., G¨ uttler, C., &amp; Blum, J. (2012). Free collisions in a microgravity many-particle experiment. I. Dust aggregate sticking at low velocities. Icarus, 218(1), 688-700.

Whipple, F. L. (1972). On certain aerodynamic processes for asteroids and comets. In From plasma to planet (p. 211).

Williams, J. P., &amp; Cieza, L. A. (2011). Protoplanetary disks and their evolution. Annual Review of Astronomy and Astrophysics, 49, 67-117.

Windmark, F., Birnstiel, T., G¨ uttler, C., Blum, J., Dullemond, C. P., &amp; Henning, T. (2012). Planetesimal formation by sweep-up: how the bouncing barrier can be beneficial to growth. Astronomy &amp; Astrophysics, 540, A73.

Wu, Y. (2019). Mass and mass scalings of superEarths. The Astrophysical Journal, 874(1), 91.

Yang, C. C., Johansen, A., &amp; Carrera, D. (2017). Concentrating small particles in protoplanetary disks through the streaming instability. Astronomy &amp; Astrophysics, 606, A80.

Yap, T. E., &amp; Tissot, F. L. (2023). The NC-CC dichotomy explained by significant addition of CAI-like dust to the Bulk Molecular Cloud (BMC) composition. Icarus, 405, 115680.

Yokoyama, T., Nagashima, K., Nakai, I., Young, E. D., Abe, Y., Al´ eon, J., ... &amp; Yurimoto, H. (2022). Samples returned from the asteroid Ryugu are similar to Ivuna-type carbonaceous meteorites. Science , 379 (6634), eabn7850.

Youdin, A. N., &amp; Goodman, J. (2005). Streaming instabilities in protoplanetary disks. The Astrophysical Journal, 620(1), 459.

Zsom, A., Ormel, C. W., G¨ uttler, C., Blum, J., &amp; Dullemond, C. P. (2010). The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals?II. Introducing the bouncing barrier. Astronomy &amp; Astrophysics, 513, A57.
Teng Ee Yap Konstantin Batygin
Dust-Gas Coupling in Turbulence- and MHD Wind-Driven Protoplanetary Disks: Implications for Rocky Planet Formation
The degree of coupling between dust particles and their surrounding gas in
protoplanetary disks is quantified by the dimensionless Stokes number. The
Stokes number (St) governs particle size and spatial distributions, in turn
establishing the dominant mode of planetary accretion in different disk
regions. In this paper, we model the characteristic St of particles across time
in disks evolving under both turbulent viscosity and magnetohydrodynamic (MHD)
disk winds. In both turbulence- and wind-dominated disks, we find that
collisional fragmentation is the limiting mechanism of particle growth, and the
water-ice sublimation line constitutes a critical transition point between dust
settling, drift, and size regimes.
  The St dichotomy across the ice-line translates to distinct planet formation
pathways between the inner and outer disk. While pebble accretion proceeds
slowly for rocky embryos within the ice-line (across most of parameter space),
it does so rapidly for volatile-rich embryos beyond it, allowing for the growth
of giant planet cores before disk dissipation. Through simulations of rocky
planet growth, we evaluate the competition between pebble accretion and
classical pairwise collisions between planetesimals. We conclude that the
dominance of pebble accretion can only be realized in disks that are driven by
MHD winds, slow-evolving, and devoid of pressure maxima that may concentrate
solids and give rise of planetesimal rings. Such disks are extremely quiescent,
with Shakura-Sunyaev turbulence parameters $\alpha_{\nu} \sim 10^{-4}$. We
conclude that for most of parameter space corresponding to values of
$\alpha_{\nu}$ reflected in observations of protoplanetary disks ($\gtrsim
10^{-4}$), pairwise collisions constitute the dominant pathway of rocky planet
accretion. Our results are discussed in the context of super-Earth origins and
Earth's accretion history.