Most major planetary bodies in the Solar System rotate in the same direction as their orbital motion their spin is prograde. Theoretical studies to explain the direction as well as the magnitude of the spin vector have had mixed success. When the accreting building blocks are  kmsize planetesimals  as predicted by the classical model  the accretion process is so symmetric that it cancels out prograde with retrograde spin contributions rendering the net spin minute. For this reason the currentlyfavored model for the origin of planetary rotation is the giant impact model in which a single collision suffices to deliver a spin which magnitude is close to the breakup rotation rate. However the giant impact model does not naturally explain the preference for prograde spin. Similarly an increasing number of spin vector measurement of asteroids also shows that the spin vector of large primordial asteroids is not isotropic. Here we reassess the viability of smaller particles to bestow planetary bodies with a net spin focusing on the pebble accretion model in which gas drag and gravity join forces to accrete small particles at a large cross section. Similar to the classical calculation for planetesimals we integrate the pebble equation of motion and measure the angular momentum transfer at impact. We consider a variety of disk conditions and pebble properties and conduct our calculations in the limits of 2D planar and 3D homogeneous pebble distributions. We find that in certain regions of the parameter space the angular momentum transfer is significant much larger than with planetesimals and on par with or exceeding the current spin of planetary bodies. We link this large net spin delivery to the appearance of asymmetries during the accretion process of pebbles. For example prograde contribution may dominate in certain regions of the parameter space because they originate from trajectories that are preferentially captured. For simplicity our calculations have ignored certain important effects e.g. collisions the backreaction on the gas and formation of atmospheres and do not address how the eventual distribution of spin vectors is obtained for which collisions and postformation processes must have played a role to explain the scatter. Irrespective of these issues pebble accretion is a viable mechanism to not only grow planetary bodies but also to impart them with a significant spin.
['Astrophysics - Earth and Planetary Astrophysics']
Spinning up planetary bodies by pebble accretion
{'Spinning up planetary bodies by pebble accretion': 'R.G. Visser a, ∗ , C.W. Ormel a , C. Dominik a , S. Ida b  \na Anton Pannekoek Institute, University of Amsterdam, Science Park 904, PO box 94249, Amsterdam, The Netherlands  \nb Earth-Life Science Institute (ELSI), Tokyo Institute of Technology, Meguro, Tokyo, 152-8550, Japan', 'Abstract': 'Most major planetary bodies in the solar system rotate in the same direction as their orbital motion: their spin is prograde. Theoretical studies to explain the direction as well as the magnitude of the spin vector have had mixed success. When the accreting building blocks are ∼ km-size planetesimals - as predicted by the classical model - the accretion process is so symmetric that it cancels out prograde with retrograde spin contributions, rendering the net spin minute. For this reason, the currently-favored model for the origin of planetary rotation is the giant impact model, in which a single collision suffices to deliver a spin, which magnitude is close to the breakup rotation rate. However, the giant impact model does not naturally explain the preference for prograde spin. Similarly, an increasing number of spin-vector measurement of asteroids also shows that the spin vector of large (primordial) asteroids is not isotropic. Here, we re-assess the viability of smaller particles to bestow planetary bodies with a net spin, focusing on the pebble accretion model in which gas drag and gravity join forces to accrete small particles at a large cross section. Similar to the classical calculation for planetesimals, we integrate the pebble equation of motion and measure the angular momentum transfer at impact. We consider a variety of disk conditions and pebble properties and conduct our calculations in the limits of 2D (planar) and 3D (homogeneous) pebble distributions. We find that in certain regions of the parameter space the angular momentum transfer is significant, much larger than with planetesimals and on par with or exceeding the current spin of planetary bodies. We link this large net spin delivery to the appearance of asymmetries during the accretion process of pebbles. For example, prograde contribution may dominate (in certain regions of the parameter space) because they originate from trajectories that are preferentially captured. For simplicity, our calculations have ignored certain important effects (e.g., collisions, the backreaction on the gas, and formation of atmospheres) and do not address how the eventual distribution of spin vectors is obtained for which collisions and postformation processes must have played a role to explain the scatter. Irrespective  \nof these issues, pebble accretion is a viable mechanism to not only grow planetary bodies, but also to impart them with a significant spin.', '1. Introduction': "With the exception of Venus, all other planets in the solar system have their spin and orbital angular momentum vectors aligned (prograde rotation). The giant planets were likely spun-up by accretion through a circumplanetary disks (Machida et al., 2008). For planets dominated by solids - the terrestrial planets as well as Uranus and Neptune - the situation is more complex. The classical model of planet formation asserts that planetary cores grew by accretion of km-sized planetesimals (Pollack et al., 1994). However, planetesimals do not provide planetary embryos with a significant amount of angular momentum (Ida & Nakazawa, 1990; Lissauer & Safronov, 1991; Dones & Tremaine, 1993a). On average, the planar components do not bestow angular momentum by symmetry, 〈 l x 〉 = 〈 l y 〉 = 0 where 〈 l i 〉 is the average momentum per unit mass transferred by the impactor in the direction i . To obtain the vertical spin component 〈 l z 〉 numerical integration of accreting planetesimal trajectories are conducted, averaging over impact parameter and phase angles. But these too display nearcancellation behavior; essentially, the point where planetesimals hit the surface of a planet is random, meaning that they are about as likely to impart positive angular momentum as negative angular momentum. In fact, the largest net contribution occurs for a dynamically cold planetesimal disk, for which the spin will be retrograde (Ida & Nakazawa, 1990) (see Figure 1); the maximum prograde contribution (also indicated in Figure 1) occurs when the planetesimals are dynamically moderately excited, such that their scaleheight extends over a scale larger than the Hill radius of the body (Dones & Tremaine, 1993a). This is more realistic, but it fails to reach the observed spins by many magnitudes. The only way for planetesimals to bestow a significant (prograde) spin is to invoke a non-uniform disk (Ohtsuki & Ida, 1998), which may materialize when the planetary embryos grow in isolation, but not in a more general setting where multiple embryos scatter planetesimals around.  \nBecause of the failure to deliver a systematic spin by accreting planetesimals, the current paradigm holds that the primordial spin of planetary bodies was instead stochastically determined by late large impacts (Dones & Tremaine, 1993b). The spin would be stochastic when it is dominated by accretion of a single body that constituted a significant fraction of the target bodies' mass, as for example in the hypothesized formation of the Earth-Moon system. Indeed, the current paradigm for the formation of the terrestrial planets asserts the existence of a giant impact phase, where many Mars-size embryos were perturbed onto crossing orbits after the dispersal of the gaseous disks (Wetherill, 1985). These equal-mass embryos collide at velocities similar to the surface escape velocity, conveying an angular momentum close to the breakup frequency ω crit = √ GM p /R 3 (Kokubo & Ida, 2007; Miguel & Brunini, 2010). Importantly, the giant impact scenario also predicts that the distribution of the spin  \nvector must be isotropic, such that the spin axis is oriented towards low ecliptic latitude, as is the case with Uranus but not with the terrestrial planets. However, this result is not considered a problem for the giant impact hypothesis, as the number of planets is limited and the spin of the terrestrial planets is known to be significantly affected by post-formation processes (Dobrovolskis, 1980). For example, the spin obliquity for Mars is known to be unstable (chaotic) on timescales short as ∼ 10 Myr (Laskar & Robutel, 1993).  \nThe asteroid belt provides us with an opportunity to inform us on the primordial spin state of solar system bodies, as the spin vector of a growing number of bodies is measured (Warner et al., 2009). Small asteroids are believed to be collisional products, which have erased any memory of the primordial spin. In addition, the spin of small asteroids is affected by post-formation processes like YORP (Rubincam, 2000). On the other hand, the spins of asteroids with diameter D > 120 km are thought to have been hardly affected by collisions or post-formation dynamics (Bottke et al., 2005; Steinberg & Sari, 2015). In Figure 1 we plot the spin vector (in terms of its rotation frequency and ecliptic latitude β ) for asteroids of D > 150 km. Although it is obvious that asteroids evolved in a highly collisional environment, Figure 1 clearly displays a preference for a prograde spin compared to what one would expect from an isotropic distribution. A Kolmogorov-Smirnov test indicates that the null-hypothesis of a uniform distribution in sin β can be rejected with a probability of > 99% . This result hints that asteroids assembled through a mechanism that provided them with a systematic (vertical) spin. Similarly to the asteroid belt, recent observations of the mutual orbit spin orientations of trans-neptunian binaries also show a preference for prograde motion (Grundy & et al., 2019) suggesting that prograde spin is the natural starting state of (solar system) planetary bodies.  \nPebble accretion (Ormel & Klahr, 2010; Lambrechts & Johansen, 2012) is an alternative theory to explain the assembly of planetary bodies. In pebble accretion bodies accrete aerodynamically small particles ('pebbles') which drift inwards through the disk by virtue of the sub-Keplerian motion of the gas (Weidenschilling, 1977a). The pebbles are then captured by the combined effects of gravity and gas drag to spiral inwards (settle) towards the protoplanet surface. Pebble accretion has been invoked to help explain the architecture of the solar system: its terrestrial planets (Levison et al., 2015a), its giant planets (Levison et al., 2015b) and its asteroids (Johansen et al., 2015).  \nThe accumulation of planetary spin by pebble accretion has not been deeply studied. Johansen & Lacerda (2010) were the first to investigate the link between accreting small particles and the spin of the accreting body in a gaseous environment. The authors performed local hydrodynamical simulations (a shearing sheet with periodic boundary conditions), involving high pebble densities resulting in instantaneous and very rapid accretion. Because of the high densities particle feedback effects are important, resulting in a prograde spinning circumplanetary disk. While these pioneering simulations are insightful in highlighting the ability of small particles to spin up bodies, they rely on rather specific conditions.  \nIn this paper we explore the implications pebble accretion has for the spin of  \naccreting bodies in a more general fashion. We calculate the angular momentum that is being transfered to the accreting body by integrating the equation of motion of individual pebbles. Summing over the range of impact trajectories we find, in this way, the net AM that is being transfered to the accreting body. Although neglecting feedback effects, our more 'basic' approach has the advantage that we can easily scan a wide parameter space and in this way identify specific trends in the simulated data. Our approach is by-and-large similar to the classical calculations for planetesimals (Ida & Nakazawa, 1990; Lissauer & Safronov, 1991; Dones & Tremaine, 1993a) but then tailored towards pebbles. Specifically, we will show in this paper that spin transfer is connected to the appearance and disappearance of asymmetries in the pebble accretion process. We show that these asymmetries lead to mean spin outcomes which are on a par with what is observed in the asteroid belt and planets.  \nThe paper is structured as follows. In Section 2 we discuss our numerical model, in Section 3 we present the obtained results, the discussion of the results is given in Section 4 and we relate our results to the solar system in Section 5. The key conclusions are presented in Section 6.", '2. Numerical model': 'To determine the amount of rotation pebbles supply to the accreting protoplanet at impact, pebble trajectories are numerically integrated with a RungeKutta-Fehlberg variable step scheme (RKF45) (Fehlberg, 1969), (used by Visser & Ormel (2016)). We use a reference frame co-moving with a protoplanet with mass M p at distance r 0 from the star (Figure 2). The three dimensional equation of motion (EOM) for a test particle with velocity v = ( v x , v y , v z ) and position r = ( x, y, z ) in the local frame is given by (Ida & Nakazawa, 1990; Ormel & Klahr, 2010):  \nd v d t =   2Ω 0 v y +3Ω 2 0 x -2Ω 0 v x 0   -GM p r 3   x y z   + F d , (1)  \nwith G the universal gravity constant, (2Ω 0 v y , -2Ω 0 v x , 0) the Coriolis acceleration, 3Ω 2 0 x the tidal force and Ω 0 the Keplerian frequency around the star. The second and last term contain the protoplanet gravity and gas drag, respectively. The stellar gravity component in the z dimension has been omitted, an approximation that is further motivated in Section 3.3. The gas drag force is given by:  \nF d = -v -v g t s , (2)  \nwith t s the stopping time, the time after which solids are slowed down by the gas. For pebbles of radius s and internal density ρ · s (Whipple, 1972):  \nt s =        ρ · s s ρ g v th Epstein regime: s < 9 4 l mfp 2 ρ · s s 2 9 η d Stokes regime: s ≥ 9 4 l mfp , (3)  \nwhere the Stokes regime is valid for low Reynolds number, Re glyph[lessmuch] 1, with η d the dynamic viscosity, l mfp the molecular mean free path, v th the thermal speed and ρ g the gas density. The gas flow is approximated as the Keplerian shear motion around the planet:  \nv g = ( -v hw -3 2 Ω 0 x )ˆ y , (4)  \nwith v hw the gas headwind which arises because the gas orbits the star with a speed slightly lower than the Keplerian speed v k due to pressure support pointing radially outwards (Weidenschilling, 1977a). A large body in an unperturbed Keplerian orbit, in a co-moving reference frame, hence faces a headwind of magnitude:  \nv hw = ηv k , (5)  \nwhere η measures the deviation from the Kepler velocity (Nakagawa et al., 1986):  \nη = c 2 s 2 v k ∂ log P ∂ log r 0 ∼ ( c s v k ) 2 , (6)  \nand P is the gas pressure in the midplane.  \nThe Hill radius of the accreting body is given by:  \nR H = r 0 ( M p 3 M glyph[star] ) 1 / 3 , (7)  \nwith M glyph[star] the mass of the central star. Normalizing the velocities to the Hill velocity R H Ω 0 , distances to the Hill radius R H and units of time to the inverse orbital frequency of the protoplanet Ω -1 0 in line with Ormel & Klahr (2010), the 3D pebble EOM in Hill units is given by:  \nd d t   v x v y v z   =   2 v y +3 x -3 x/r 3 -2 v x -3 y/r 3 -3 z/r 3   -1 τ s   v x v y + ζ w +3 x/ 2 v z   , (8)  \nvelocities and positions are now represented as dimensionless quantities, with GM p = 3 in Hill units. The dimensionless headwind parameter ζ w is given by (Ormel & Klahr, 2010):  \nζ w ≡ v hw Ω 0 R H , (9)  \nwhich is therefore an indicator of the mass (gravity) of the protoplanet. The physical protoplanet radius is normalized by the Hill radius:  \nα p ≡ R R H . (10)  \nand the Stokes number is given by:  \nτ s = t s Ω 0 . (11)  \nIn these units, the input parameters for our integrations are reduced to ζ w , α p and pebble Stokes number τ s . A direct relation between the main physical and dimensionless quantities is summarized in Table 1.  \nIn general specific values of disk parameters are not known and vary greatly since protoplanetary disks are subject to complex time evolution during their lifetime (e.g. (Bitsch et al., 2015)). Nevertheless, we will occasionally transform our results into physical quantities adopting power-law profiles for both the gas temperature and gas surface density profile following from the the Minimum Mass Solar Nebula (MMSN) (Weidenschilling, 1977b; Hayashi et al., 1985; Nakagawa et al., 1986):  \nT ( r ) = 170 K ( r 0 1 au ) -1 / 2 , (12)  \nΣ( r ) = 1300 g cm -2 ( r 0 1 au ) -3 / 2 . (13)  \nThe disk gas scaleheight is given by:  \nH g = c s Ω 0 . (14)  \nAssuming an isothermal column, the gas distribution in the vertical dimension z is equal to:  \nρ g = Σ H g √ 2 π exp [ -1 2 ( z H g ) 2 ] . (15)  \nThe main input parameters ζ w , α p and pebble Stokes number τ s can now be expressed in terms of the disk parameters as:  \nv hw = 30 ms -1 ( µ 2 . 34 ) -1 ( M glyph[star] M glyph[circledot] ) -1 / 2 , (16)  \nζ w = 12 . 5 ( ρ · 2 . 5 g cm -3 ) -1 / 3 ( v hw 30 m s -1 ) ( R 100 km ) -1 ( r 0 au ) 1 / 2 , (17)  \nα p = 5 . 7 × 10 -3 ( M glyph[star] M glyph[circledot] ) 1 / 3 ( ρ · 2 . 5 g cm -3 ) -1 / 3 ( r 0 au ) -1 , (18)  \nτ s = 0 . 1 ( s 1 cm ) ( ρ · s 2 . 5 g cm -3 )( M glyph[star] M glyph[circledot] ) 1 / 2 ( r 0 au ) -3 / 2 , (19)  \nwith µ the mean molecular weight in atomic mass units and where we applied the ideal gas law.', '2.1. Pebble-to-planet spin supply': 'The spin supply of pebbles to the protoplanet is calculated by tracing individual specific angular momentum (SAM) contributions from pebbles at impact on the protoplanet. An individual pebble heading for the protoplanet bestows a SAM of:  \nl = r × v . (20)  \nThe z -component is the one of interest:  \nl z = ( xv y -yv x ) + Ω 0 ( x 2 + y 2 ) , (21)  \nwhere the latter term is to correct for the rotation of the co-moving coordinate frame (Dones & Tremaine, 1993a). The mean SAM transferred to the protoplanet is found by averaging the SAM contributions over the accretion cross-section. For example, in 2D:  \n〈 l z 〉 = ∫ accret F ( x 0 ) l z ( x 0 )d x 0 ∫ accret F ( x 0 )d x 0 , (22)  \nwith F ( x 0 ) = v y, ∞ Σ the flux of pebbles entering our domain and Σ the pebble surface density. It is useful to normalize the specific mean spin angular momentum 〈 l z 〉 with l z , esc = √ 2 GMR for the analysis of the results. Converting this to a fraction of the break-up (critical) frequency gives:  \nω z ω crit ≈ 3 . 5 〈 l z 〉 l z , esc , (23)  \nwith ω crit = √ GM/R 3 .', '2.1.1. Initial conditions': 'Integration of pebble trajectories start at a position x = ( x 0 , y 0 , z 0 ) (see Figure 2). We take y 0 far from the protoplanet with radius R , to ensure that the planet gravity is initially insignificant. Pebbles are initiated with the unperturbed azimuthal and radial drift velocities (Weidenschilling, 1977a):  \nv x, ∞ = -2 v hw τ s 1 + τ 2 s , (24)  \nv y, ∞ = -v hw 1 + τ 2 s -3 2 Ω 0 x 0 . (25)  \nThe vertical release velocity v z, ∞ = 0. To ensure that the pebbles are initially not influenced by the protoplanet gravity we demand that the gas drag force is larger than the two-body gravitational force of the protoplanet by a safety factor C. Specifically, we fix y 0 = C √ GMt s /v hw , with C = 100 (Ormel & Klahr, 2010). We then scan along the y = y 0 line to determine for which range in x 0 pebbles are accreted. The condition for accretion of a pebble onto the protoplanet is given by:  \nr -R ≤ 0 . (26)  \nIn this way we find x 1 and x 2 , the interior, respectively, exterior edge of the protoplanet impact range.  \nIn the 3D case the collision range is symmetric around z = 0 and we release pebbles over an interval z 0 ∈ [ -b z , b z ] while ensuring that misses are included on both sides of the impact range.', '3. Results': 'The main result of this study are highlighted in Section 3.1, which summarizes the result of our numerical integrations (Figure 3). In Section 3.2 (2D) and Section 3.3 (3D) we provide an in-depth trajectory analysis, which allows us to understand the trends that we observe in the simulated data. We present test cases with earlier work in Section 3.4.', '3.1. Spin contributions': "In Figure 3 the imparted spin as function of pebble stopping time τ s and headwind-to-shear parameter ζ w is plotted for α p = 2 × 10 -3 , which corresponds to a distance of about 2.5 au from a solar-mass star and a planet density of 2.5 g cm -3 . The top panel of Figure 3 provides ω z /ω crit for a planar configuration (pebbles having settled to a midplane), whereas the bottom panel gives ω z /ω crit in the 3D limit, applicable when the pebble scaleheight exceeds the impact radius for pebble accretion. The 3D configuration is therefore more applicable for small pebbles that experience strong turbulence and for small planets. These results show that large prograde spin can be obtained even beyond the levels seen in the solar system's bodies. For example, the spin of Ceres and other asteroids ( ω/ω crit ∼ 0 . 3; thick contour) is well matched by accretion of τ s ≈ 0 . 05 pebbles either in the 2D or 3D mode. The spin of smaller asteroids ( R ∼ 100 km) can also be attributed to pebble accretion, provided that the headwind velocity is less than 30 m s -1 . Retrograde spins are also possible, albeit in small regions of the parameter space. Figure 3b shows 3D accretion of larger pebbles may spin up planets up to their breakup rate. Values ω z > ω crit are however artificial; in reality the pebble will bounce off to be decelerated by gas drag until it is re-accreted at the critical spin frequency, resulting in these bodies to rotate at breakup speeds. In any case, it is more likely that in those cases the interaction geometry follows the 2D limit.  \nPebble accretion is in fact quite efficient to deliver a systematic spin. The degree of asymmetry between prograde and retrograde spin contributions, δ = 〈 l z 〉 / √ 〈 l 2 z 〉 exceeds δ = 0 . 1 over a large part of the parameter space (see Figure 8), whereas for planetesimal accretion δ ∼ α p ∼ 10 -3 (Dones & Tremaine, 1993a), see Section 3.4. A qualitative understanding of this result follows from considering asymmetries in the interaction timescale of pebbles that encounter the planet on an interior vs exterior trajectory, see Figure 4. Pebble accretion first starts to operate on the exterior trajectories, which are flung back in a direction ≈ 45 · against the disk headwind due to the drift angle. Hence these encounters last long, resulting in their capture, while the interior trajectories  \nmomentarily avoid accretion. This explains the prograde band at R = R PA , the radius where pebble accretion commences (Visser & Ormel, 2016). For larger planets, on the other hand, encounters last long for the pebbles entering the Hill sphere on interior orbits due to the importance of the shear motion (Figure 4b). This again results in a prograde net contribution, which is even stronger in the 3D (see Figure 4b and Section 3.3)", '3.2. Trajectory analysis for τ s = 0 . 1 pebbles': 'In the co-moving frame we classify pebbles that cross the x-axis initially from x glyph[lessorsimilar] -R as interior trajectories and pebbles that cross the x-axis initially from x glyph[greaterorsimilar] R as exterior trajectories . There are three directions in which spin can be imparted to the protoplanet by a pebble:  \n- 1. pebbles fall in counter-clockwise to the protoplanet imparting spin aligned with the orbital rotation Ω 0 (prograde)\n- 2. pebbles fall in perpendicular to the protoplanet surface, imparting no net spin\n- 3. pebbles fall in clockwise to the protoplanet imparting spin opposite to the orbital rotation Ω 0 (retrograde)  \nIn the case of symmetry the sum of individual pebble spin contributions over the collision cross-section is zero. However, as it turns out several asymmetries produce non-zero mean spin outcomes.  \nWe use the example model for fixed τ s = 0 . 1 and α p = 0 . 003 (Fig.5, blue curve) to analyze pebble trajectories. We select five protoplanet radii from this curve corresponding to the most extreme transition between prograde and retrograde outcomes throughout the curve. For each protoplanet radius we then select several pebble trajectories ranging from an interior miss, a range of impacts and an exterior miss (Fig. 6a to e, left panel). The right panel of Figure 6a to e shows the specific angular momentum an individual pebble imparts at impact on the protoplanet (every red dot), with respect to the release distance x 0 within the collision cross-section. The impacting trajectories in the left panel correspond to the highlighted dots with equal color in the right panel, going from dark (most interior trajectory) to light (most exterior trajectory). The highlighted dots in the right panel are uniquely labelled for reference in the text.  \nAs protoplanet radii in our simulations range from R = 10 km to R =3000 km there are a range of different accretion regimes which we will use to structure the analysis. These regimes are also shown in Figure 8 where we show the degree of asymmetry between prograde and retrogade contributions for the 2D case (blue curve). Specifically, for our parameter space the following regimes are relevant:  \n- (a) geometrical and ballistic accretion (Safronov, 1972), R glyph[lessorsimilar] 200 km(Fig. 6a)\n- (b) onset pebble accretion (Visser & Ormel, 2016), 200 km glyph[lessorsimilar] R glyph[lessorsimilar] 400 km(Fig. 6b)  \n- (c) pebble accretion (Ormel & Klahr, 2010; Lambrechts & Johansen, 2012), 400 km glyph[lessorsimilar] R glyph[lessorsimilar] 1000 km(Fig. 6c)\n- (d) transition to shear dominated pebble accretion, R ∼ 1000 km(Fig. 6d)\n- (e) shear dominated pebble accretion (Ormel & Klahr, 2010; Lambrechts & Johansen, 2012) R glyph[greaterorsimilar] 1400 km(Fig. 6e)', 'Regime (a): Geometrical accretion and ballistic accretion': 'For both geometrical and ballistic impacts, prograde pebble spin imparted by impacts from the interior band (Fig. 6a, trajectories between a0 and a1) is canceled by the exterior band ((Fig. 6a, trajectories between a1 and a2), which deliver retrograde spin. The sum of all the spin contributions sum to zero and symmetry prevails.', 'Regime (b): Onset pebble accretion': 'Pebble accretion commmences at protoplanet radius (Visser & Ormel, 2016):  \nR PA ≈ 160 km ( v hw 30 m s -1 ) ( ρ · 2 . 5 g cm -3 ) -0 . 36 ( r 0 au ) 0 . 42 ( τ s 0 . 1 ) 0 . 28 , (27)  \nwhich is valid for τ s < 1. From R PA , pebbles at interior trajectories ((Fig. 6b, trajectories between b0 and b1) start to differ from pebbles at exterior trajectories((Fig. 6b, trajectories between b1 and b3). In particular trajectory b0 is still accreting in a ballistic fashion equivalent to trajectory a0. The exterior trajectories on the other hand (b2 and b3) are captured last moment already entering the settling regime, enhancing the exterior accretion band. That (axial) symmetry is broken can also be seen from two trajectories that both reach the same minimum distance r min upon first approach to the protoplanet. The interior blue trajectory escapes, while the red exterior trajectory hits the protoplanet with a prograde contribution (Fig. 4). This asymmetry between interior and exterior encounters can be understood from the angle at which the drifting pebbles approach the protoplanet. For this reason the escaping interior trajectory is deflected more downwards and trajectory b2 is deflected more upwards . Consequently, encounter (crossing) times are longer for trajectory b2, resulting in its capture.  \nNaively it would be expected that the exterior trajectories b2 and b3 accrete clockwise (retrograde) by chasing spirals. Instead they fly across the protoplanet again, reverse orbit and wrap around the planet from below to impart a significant amount of counter-clockwise (prograde) spin. Effectively it can be conluded that the exterior band already experiences a transition from ballistic encounters (too fast to be captured) to settling encounters in which the protoplanet gravity is able to capture the pebble. Because this is the transition to settling, pebbles are captured at the verge of escape, leading to the last moment orbit reversal back to the protoplanet surface, to impart prograde spin. The onset of pebble accretion stated in Equation 27 is therefore essentially only valid for  \nthe exterior band. At somewhat larger protoplanet radius (gravity) R glyph[greaterorsimilar] 300 km, the interior band catches up by also accreting pebbles through the settling mechanism, equivalent to the exterior band between b2 and b3. The interior band imparts spin in a retrograde fashion, restoring the symmetry and bringing mean spin 〈 l z 〉 close to zero again (Fig. 6b-right panel).', 'Regime (c): Pebble accretion': "For R = 400 km pebble trajectories mean spin is dominated by retrogade contributions due to an asymmetry associated with the in-spiralling process. Trajectories c1 and c3 are captured and spiral inwards due to settling. However, before even finishing half a spiral these pebbles already collide with the protoplanet surface. (Fig. 6c-left and right panel). Again the interior trajectory c1 is deflected more downwards than the exterior trajectory c3. Consequently, trajectory c1 is prevented from delivering a large amount of prograde rotation since it is guided toward a perpendicular impact to the protoplanet surface. Trajectory c3 on the other hand is able to wrap around the protoplanet, finally 'grazing' the protoplanet surface at impact to impart significant retrograde spin. With increasing gravity the symmetry is again restored for R ∼ 500 km as pebbles completely spiral in (settle) (like the c0 and c4 trajectories) for both the interior and exterior bands.", 'Regime (d): Transition to shear dominated pebble accretion': 'For protoplanet radius R ∼ 1000 km the background Keplerian shear velocity ( v sh = -1 . 5Ω 0 x ) becomes comparable to the gas headwind velocity v hw since collision cross-sections reach a significant fraction ( ∼ 0 . 5 R H ) of the protoplanet Hill sphere. The x -coordinate for which the gas headwind velocity equals the upwards Keplerian shear layer is given by:  \nx co = 2 3 Ω -1 v hw , (28)  \nwhich is referred to as the co-rotation line. As pebbles drift closer to this line, their y-velocity becomes decreasingly negative to eventually follow horseshoes upwards (gray missing trajectories in Fig. 6d-left panel).  \nInterior pebbles experience longer encounter times as they approach the co-rotation line while drifting inwards. As long as x co < -R H this enhances the interior accretion (the band between trajectory d0 and d2 in Fig. 6d) since gravity has a long time to capture these pebbles. These pebbles accrete prograde because they fall in with counter-clockwise spirals. Exterior pebbles on the other hand have a more reduced encounter time; they approach in shear layers with a higher downwards velocity (the band between trajectory d2 and d4 in Fig. 6d). As a consequence, the interior band (between d0 and d2) is visibly broader than the exterior one (between d2 and d4) as shown in the right panel of Fig. 6d. Since the interior pebbles spiral in counter-clockwise to the protoplanet, the mean spin outcome is prograde. Thus, the transition to the shear-regime implies another asymmetry, which favors prograde contribution.', 'Regime (e): Shear dominated pebble accretion': 'For R ∼ 3000 km the Hill radius extends beyond the co-rotation line (Fig. 6e, left panel); collision cross sections are now entirely determined by the Keplerian shear. Although a small part of the interior band (Fig. 6e, trajectories e0 to e1) still accretes in a prograde fashion a significant part has been lost after crossing the co-rotation line to follow upwards horseshoes (interior gray miss, Fig. 6e, left panel). Consequently, the retrograde contribution from exterior trajectories dominate.  \nA subset of the exterior trajectories, however, now accrete from the back. After having traveled through the Hill sphere once, they cross the co-rotation line due to radial drift. The pebble then travels upwards from -y and enters the Hill sphere a second time (trajectory e4, Fig. 6e, left panel). These backaccreting pebbles will always supply prograde angular momentum, as we can see easily in the following way: The protoplanet orbits with the Keplerian velocity while the pebble is pulled from x ∼ x co to x = R , starting with v y ∼ 0 at the co-rotation line and then, while approaching the planet, moving toward the downward headwind speed and the downwards shear field combined. This forces the pebble to impact the protoplanet in a counter-clock-wise fashion resulting in a maximum prograde contribution to the angular momentum (see Figure 7 for clarification). Nonetheless, the total prograde contributions do not outweigh the retrograde ones and the mean spin outcome is dominantly retrograde.  \nFor even higher protoplanet gravity the Hill sphere extends even further beyond the co-rotation line making the secondary contributions such as trajectory e4 take over the full accretion band. As a consequence mean spin outcomes are effectively prograde again as shown in Figure 5 for R glyph[greaterorsimilar] 4000 km.', '3.3. Analysis of the 3D simulations': "Generally, pebbles are expected to reside in a sub-disk, in which the pebble sedimention is balanced by turbulent diffusion (Dubrulle et al., 1995). The thickness of the particle disk is therefore determined by the aerodynamical properties of the pebbles ( τ s ) as well as the strength of the turbulence (usually characterized in terms of a diffusivity D = α z c s H ). Therefore the pebbles should be modelled according to a Langevin equation (Ormel & Liu, 2018) or with hydrodynamical simulations (Homann et al., 2016). We opt to avoid these formal but complicated models by omitting the vertical stellar gravity component. In the absence of the protoplanet gravity, pebbles would then preserve their height ( z ). In addition, in our 3D calculations we assume an homogeneous distribution of pebbles in z , i.e., we assume that the pebble scale height is significantly larger than the Hill radius of the protoplanet. In the 3D case turbulent trajectories will be ballistic ('straight') over timescales less than the turbulence correlation time. Since the latter is typically ∼ Ω -1 in most turbulence models Cuzzi et al. (2001), longer than the settling timescale, 3D accretion onto small planets/asteroids can be approximated as ballistic. Then, our 2D and 3D simulations constitute the two limits of the general situation.  \nIn Figure 9 3D simulations for α p = 3 × 10 -3 and for τ s = 0 . 1 show that mean spins are dominantly more prograde in the 3D case with respect to the  \n2D case. We provide a similar analysis as applied to the 2D results for the most notable transitions in the 3D spin curves shown in Figure 5. We will show that most of the features in the 2D case can be extended to the 3D case. Heat maps of individual spin contributions l z,i at impact are shown in Figure 10 for initial release distance z 0 on the y-axis and initial release distance x 0 on the x-axis.  \nThe physical processes responsible for the behavior of the 2D spin curve in section 3.2 are equivalent in 3D. The trajectory types in the 2D case described in Figure 6 apply to a mid-plane slice of these heat maps ( z 0 = 0). The trajectory types are however non-changing for non-zero release distance z 0 , while keeping x 0 fixed. This is indicated in Figure 10 with the vertical solid lines with the corresponding unique label shown in the 2D case (Fig. 6).  \nglyph[negationslash]  \nThe main difference in 3D, favouring the prograde spin, is the increase in the encounter time discrepancy between interior and exterior pebbles due to shear. In the shear important regime interior pebbles released at ( x 0 , y 0 , z = 0) drift to the co-rotation line (trajectory e0, for example), gaining small velocities and consequently long encounters. This effect remains unchanged in 3D since the influence of Keplerian shear and drift on a pebble are equivalent for release distance ( x 0 , y 0 , z = 0) (Fig. 10 D , trajectory e0). The long lasting interior encounters for z 0 = 0 lead to timely vertical settling to the mid-plane and a strong l z contribution. In particular the following effects are observed in 3D:  \nglyph[negationslash]  \n- 1. Analogous to the 2D results, the mean spin is effectively zero for the geometrical and ballistic regime(Fig. 10 A ).\n- 2. A larger prograde contribution compared to the 2D case due to the transition from a 1D to a 2D collision cross-section (larger surface ratio, Fig. 10 B ),\n- 3. Shear dominated regime; comparatively longer encounter times for interior pebbles (trajectory d0 to d2 and e0 to e1), which are accreted accross a larger z . Vertical settling is made possible for interior pebbles since they drift closer towards the co-rotation line, promoting long lasting encounters. This acts strongly in favour of prograde spin supply, see Fig 10 C,D .\n- 4. (Conversely) shorter encounter times for exterior pebbles (trajectory d2 to d4 and e1 to e3) ; and a smaller 'accretion height' (lower z ) due to the weakening gravity of the protoplanet with increasing vertical release distance.", '3.4. Accreting larger pebbles and planetesimals': 'Next, we consider varying the size of the accreted particles (in terms of the dimensionless stopping time τ s ) while fixing the size of the body ( ζ w ). We consider a range of τ s ∈ [10 -2 , 10 4 ] in the case of zero inclination and eccentricity, i = e = 0. The high τ s limit effectively implies a vanishing gas-drag force and we expect therefore that our results must converge to the i = 0 , e = 0 limit that have been previously calculated (Ida & Nakazawa, 1990; Dones & Tremaine, 1993a). Our results are presented in Figure 10, where, following  \nDones & Tremaine (1993a), we have normalized with the product of the Hill velocity and the body radius R ; 〈 ˜ l z 〉 = 〈 l z 〉 /R H Ω 0 R to compare results.  \nFor 0 . 01 glyph[lessorsimilar] τ s glyph[lessorsimilar] 0 . 1 pebbles fall in symmetrically both from the interior and exterior bands since drift is negligible. However, the spin outcome is retrograde due to higher exterior flux because of shear. For τ s ∼ 0 . 1 the co-rotation line eats away prograde accretion as discussed in previous sections. Combined with the higher exterior flux the spin outcome is even more retrograde. For 0 . 1 < τ s < 1 back-accretion from the exterior band delivers high amounts of prograde spin due to increasing drift and higher exterior flux. For τ s > 1 drift decreases and chaotic bands are formed, meaning that the collision crosssection can be broken up by bands of missing trajectories (Ida & Nakazawa, 1990). Furthermore, the gas becomes less important and the co-rotation line therefore shifts more to the x = 0 line which supresses prograde contributions as pebbles follow parabolas from the interior band. The exterior pebbles are still moderately focused towards the planet by gas, leading to a strong retrograde infall. As both gas and radial drift effects diminish, the spin converges to the retrograde value found in previous works (Ida & Nakazawa, 1990; Lissauer & Safronov, 1991; Dones & Tremaine, 1993a).', '4.1. Caveats': "The results have been obtained in the approximation of a local shearing sheet patch for which Hill's approximations apply. This approach is perfectly valid for x 0 and y 0 glyph[lessmuch] than the orbital distance to the central body r 0 (Liu & Ormel, 2018). The release distance y 0 increases for larger protoplanet radii due to the growing Hill radius to ensure that pebble trajectories are governed by gas-drag initially. If y 0 reaches a significant fraction of the orbital radius r 0 , the curvature of the protoplanetary disk cannot be neglected anymore.  \nTo ensure the validity of our findings in the shearing sheet approximation, simulations are conducted in the global frame for the model α p = 3 × 10 -3 and τ s = 0 . 1, see Figure 12. We applied the same procedure as Liu & Ormel (2018) who determined pebble accretion efficiency's in the global frame, referring to the global equations of motion of a pebble in the stellar frame and initial release conditions of the pebbles. The results in a global fashion are the same as our local simulations showing that the local frame produces valid outcomes for the intermediate to lower range of Stokes numbers. We do note that results might start to deviate for Stokes numbers higher than τ s = 1 due to the transition to gas-free accretion.  \nCollisions among pebbles within the Hill sphere have been ignored, which may result in the formation of circumplanetary pebble disks (Schlichting & Sari, 2007) . This would in particular be important for small pebbles and large pebble fluxes; if bodies formed from τ s = 0 . 1 particle over a protracted period, however, collisions are unimportant. The timescale on which pebbles with mass  \nm p = (4 / 3) πs 3 ρ · s and radius s collide is given by:  \nt coll = ( nσ ∆ v ) -1 , (29)  \nwith n = ρ p /m p the number density of pebbles, σ = 4 πs 2 the mutual crosssection of pebble colliders and ∆ v the typical velocity dispersion of the pebbles. Comparing this to the settling time of a pebble t sett = b coll /v sett with b coll the impact parameter, gives:  \nN = t sett t coll = 3 b coll ρ p ρ · s s ∆ v v sett , (30)  \nwith N the number of collisions pebbles undergo before settling onto the protoplanet. In terms of the disk parameters this can be written as:  \nN ∼ 0 . 05 ( q p 10 -9 ) 1 / 3 ( τ s 0 . 1 ) 2 / 3 ( h 0 . 05 )( ρ p ρ g )( ∆ v v settl ) , (31)  \nwhere q p = M p /M glyph[star] and h = H/r 0 is the disk aspect ratio. Furthermore we assume that v th ∼ c s = H Ω 0 , pebbles are in the Epstein regime, and that the Hill limit applies: b coll ∼ R H τ 1 / 3 (Ormel & Klahr, 2010). For asteroids N ∼ 0 . 05, unless τ becomes very small. On the other hand, for larger bodies pebbles take increasingly longer to spiral in to the accreting body increasing ρ p /ρ g . This renders it likely that a disk will form around the body analogous to Johansen & Lacerda (2010). Thus, depending on the parameters collisions could become important and the effect on imparted spin should be considered in future work.  \nIn this work pebble feedback on the gas has not been considered. Indeed, if the solid to gas ratio is of the order ρ p /ρ g glyph[greaterorsimilar] 1 this alters dynamics of pebbles and the gas as the latter is forced to accrete with the solids in a prograde disk around the central body (Johansen & Lacerda, 2010). However, in the pebble accretion regime pebbles spiral in rapidly to the protoplanet with the radial velocity v r ∼ 1 /r 2 . This indicates that ρ p should not drastically increase because of the high accretion flux. Furthermore, our study is focused towards a more general description of spin delivery in the pebble accretion framework by considering the net sum of individual spin contributions of pebbles. We therefore have ignored collective effects for simplicity.  \nWe assume no initial rotation of the protoplanet in determining mean spin outcomes. However, bodies with radii R glyph[lessorsimilar] 100 km experience a growth barrier as they are too inefficient to accrete pebbles (Visser & Ormel, 2016). The angular momentum imparted during planetesimal formation could change the total spin outcomes obtained in the pebble accretion framework. Recent numerical results show however that trans-neptunian binaries forming through the streaming instability also preferentially orbit prograde around their mutual center of mass (Nesvorn'y et al., 2019). Furthermore we can assess the importance of the imprint of planetesimal spin by considering the total angular momentum supplied to a body of mass M . This quantity is given by L = ∫ 〈 l z 〉 dm . Rewriting this expression in logarithmic space gives L = ∫ 〈 l z 〉 md log m , showing that the  \nintegration of the mean angular momentum in Figure 3 is weighted by mass. Assuming that most of the mass of asteroids and protoplanets come from subsequent PA, rather than the initial planetesimal formation process, the spin state of these bodies will be given by Figure 3.  \nWhen atmospheres form around the body, the gas flow becomes more complex (Ormel et al., 2015a,b; Cimerman et al., 2017; Lambrechts & Lega, 2017; Kurokawa & Tanigawa, 2018; Popovas et al., 2018). These complicated but more realistic sub-structures of gas streamlines within the protoplanet Hill sphere are not considered in this work. Atmospheric structures do become important if R b /R glyph[greaterorsimilar] 1, where R b = GM/c 2 s is the Bondi radius of the body. For R b /R glyph[greaterorsimilar] 3 we estimate that the atmosphere has become dense enough to quantitatively alter the pebble trajectories and hence the spin contributions. While these effects are important to consider in future work for planets, the spin outcomes we obtain for asteroids are unaffected since R b /R glyph[lessorsimilar] 1.", '5. Solar system bodies': "For asteroids, we can use our general results presented in Figure 3 to constrain the conditions that operated during pebble accretion. First, we require relatively large pebbles τ s ∼ 0 . 1; much smaller pebbles would not impart a significant net angular momentum, since they barely drift and shear across the impact cross section is small. Silicates particles, being less sticky, generally have a much lower τ s . It is generally assumed that icy particles stick more readily and at larger collision velocities than silicate particles, so the need for relatively large pebbles could point us to icy pebbles. However, a set of experiments by Musiolik & Wurm (2019) indicate that the improved sticking only occurs in a relatively narrow temperature range at temperatures above ∼ 170 K , so this enhancement might only be relevant near the snow line if these results turn up to hold in general. Also, a reduction of the disk headwind lowers relative velocities between grains possibly stimulating growth to larger (icy) pebbles. This would be consistent with the finding that (large) asteroids were icy bodies with a significant post-formation processing of the ice, making them dry (Schmidt et al., 2009). ALMA observations also provide indications of the presence of decimeter sized pebbles in relatively young disks (Zhang et al., 2015). Nevertheless, the sticking properties of ice remain controversial and it is important to pursue further experimental and numerical research on the topic in the near future. Second, for pebble accretion to operate on bodies ∼ 100 km we infer that the disk headwind parameter v hw cannot be too high. Our standard value of 30 m s -1 already implies pebble accretion will fail for bodies smaller than 200 km in radius (assuming τ s = 0 . 1). However, the disk headwind is rather uncertain and could have been much smaller, in which case pebble accretion would operate on a much larger range of asteroids. This finding is again consistent with a cold formation environment.  \nFor asteroids, the post-formation sublimation of ice could alter the rotation state of asteroids, although it is unlikely that this affects the rotation rate of large asteroids (Jewitt, 1997). For planets, our results are in good agreement  \nwith the observed spin periods of both the terrestrial planets Earth and Mars, and the ice giants Uranus and Neptune. Not only does this suggest that pebble accretion naturally provides systematic spins to (exo)solar system bodies, it provides an additional validation for pebble accretion as a planet formation mechanism. Undoubtedly, impacts, outgassing, long-term chaotic evolution, and the formation of a planet atmosphere, which leads to different pebble aerodynamics and eventually ablation of pebbles (Brouwers et al., 2018) have shaped the periods and orientation of the spin axes of the solar system's planetary bodies. Nonetheless, the evidence in the solar system renders it likely that bodies started off with a systematic vertical spin provided by pebble accretion.", '6. Conclusions': 'We have calculated the spin angular momentum per unit mass that pebbles supply to a protoplanet during accretion, both in a 2D and 3D fashion. The net angular momentum gained by the protoplanet is obtained by summing individual pebble supply to the body at impact. The main conclusions of the obtained results are as follows:  \n- 1. Pebble accretion can deliver significant systematic rotation to an accreting body for a large part of the parameter space. In our results, spin values reach the value of spin periods typically observed in the Solar System both in magnitude and direction. The results are in good agreement with the observed spin rotation of the larger asteroids such as for example Ceres, Pallas, Vesta, and Hygiea.\n- 2. The cause of the asymmetry in spin accumulation by the accreting body, is a discrepancy in encounter times of pebbles that accrete interior and pebbles that accrete exterior to the collision cross-section. The difference in encounter times follows from the asymmetric infall angle of pebbles due to radial drift around the onset of pebble accretion and the large discrepancy over the collision cross-section due to keplerian shear. The asymmetry favors prograde spin rotation by the accreting body and, in a smaller region of the parameter space, favors retrogade spin rotation.\n- 3. The degree of asymmetry between prograde and retrograde spin contributions, δ = 〈 l z 〉 / √ 〈 l 2 z 〉 reaches values of δ ∼ ± 0 . 3. This is significantly higher than δ = 10 -3 that has been found in gas-free planetesimal accretion.\n- 4. In 3D there is a propensity for prograde rotation, because the pebbles are accreted preferentially from the (prograde) side characterized by long lasting encounters.\n- 5. In generally, small pebbles ( τ s glyph[lessmuch] 0 . 1) bestow much smaller spin than larger pebbles, due to the more symmetric impact geometry. Applying our results to asteroids, we infer that they formed in a cold environment  \nand accreted pebbles with significant ice fractions. Much of the ice was subsequently lost by post-formation processes.', 'References': "Baer, J., Chesley, S., & Britt, D. (2012). Asteroid Masses V3.0. NASA Planetary Data System , 191 , EAR-A-COMPIL-5-ASTMASS-V3.0.  \nBitsch, B., Johansen, A., Lambrechts, M., & Morbidelli, A. (2015). The structure of protoplanetary discs around evolving young stars. Astronomy and Astrophysics , 575 , A28. doi: 10.1051/0004-6361/201424964 . arXiv:1411.3255 .  \nBottke, W. F., Durda, D. D., Nesvorn'y, D., Jedicke, R., Morbidelli, A., Vokrouhlick'y, D., & Levison, H. (2005). The fossilized size distribution of the main asteroid belt. Icarus , 175 , 111-140. doi: 10.1016/j.icarus.2004.10.026 .  \nBrouwers, M. G., Vazan, A., & Ormel, C. W. (2018). How cores grow by pebble accretion. I. Direct core growth. Astronomy and Astrophysics , 611 , A65. doi: 10.1051/0004-6361/201731824 . arXiv:1708.05392 .  \nCimerman, N. P., Kuiper, R., & Ormel, C. W. (2017). Hydrodynamics of embedded planets' first atmospheres - III. The role of radiation transport for super-Earth planets. Monthly Notices of the RAS , 471 , 4662-4676. doi: 10. 1093/mnras/stx1924 . arXiv:1707.08079 .  \nCuzzi, J. N., Hogan, R. C., Paque, J. M., & Dobrovolskis, A. R. (2001). Size-selective Concentration of Chondrules and Other Small Particles in Protoplanetary Nebula Turbulence. Astrophysical Journal , 546 , 496-508. doi: 10.1086/318233 . arXiv:astro-ph/0009210 .  \nDobrovolskis, A. R. (1980). Atmospheric tides and the rotation of Venus. II Spin evolution. Icarus , 41 , 18-35. doi: 10.1016/0019-1035(80)90157-8 .  \nDones, L., & Tremaine, S. (1993a). On the origin of planetary spins. Icarus , 103 , 67-92. doi: 10.1006/icar.1993.1059 .  \nDones, L., & Tremaine, S. (1993b). Why does the earth spin forward? Science , 259 , 350-354. doi: 10.1126/science.259.5093.350 .  \nDubrulle, B., Morfill, G., & Sterzik, M. (1995). The dust subdisk in the protoplanetary nebula. Icarus , 114 , 237-246. doi: 10.1006/icar.1995.1058 .  \nFehlberg, E. (1969). . NASA-TR-R-313. , .  \nGrundy, W. M., & et al. (2019). Mutual orbit orientations of transneptunian binaries. Icarus , . doi: https://doi.org/10.1016/j.icarus.2019.03.035 .  \nHayashi, C., Nakazawa, K., & Nakagawa, Y. (1985). Formation of the solar system. In D. C. Black & M. S. Matthews (Ed.), Protostars and Planets II (pp. 1100-1153). Univ. of Arizona Press, Tuscon.  \nHomann, H., Guillot, T., Bec, J., Ormel, C. W., Ida, S., & Tanga, P. (2016). Effect of turbulence on collisions of dust particles with planetesimals in protoplanetary disks. Astronomy and Astrophysics , 589 , A129. doi: 10.1051/0004-6361/201527344 . arXiv:1602.03037 .  \nIda, S., & Nakazawa, K. (1990). Did rotation of the protoplanets originate from the successive collisions of planetesimals? Icarus , 86 , 561-573. doi: 10.1016/ 0019-1035(90)90233-Y .  \nJewitt, D. (1997). Cometary Rotation: an Overview. Earth Moon and Planets , 79 , 35-53. doi: 10.1023/A:1006272914117 .  \n- Johansen, A., & Lacerda, P. (2010). Prograde rotation of protoplanets by accretion of pebbles in a gaseous environment. Monthly Notices of the RAS , 404 , 475-485. doi: 10.1111/j.1365-2966.2010.16309.x . arXiv:0910.1524 .  \nJohansen, A., Mac Low, M.-M., Lacerda, P., & Bizzarro, M. (2015). Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion. Science Advances , 1 , 15109. doi: 10.1126/sciadv.1500109 . arXiv:1503.07347 .  \nKokubo, E., & Ida, S. (2007). Formation of Terrestrial Planets from Protoplanets. II. Statistics of Planetary Spin. Astrophysical Journal , 671 , 2082-2090. doi: 10.1086/522364 .  \nKurokawa, H., & Tanigawa, T. (2018). Suppression of atmospheric recycling of planets embedded in a protoplanetary disc by buoyancy barrierc. Monthly Notices of the RAS , 479 , 635-648. doi: 10.1093/mnras/sty1498 . arXiv:1806.01695 .  \nLambrechts, M., & Johansen, A. (2012). Rapid growth of gas-giant cores by pebble accretion. Astronomy and Astrophysics , 544 , A32. doi: 10.1051/ 0004-6361/201219127 . arXiv:1205.3030 .  \nLambrechts, M., & Lega, E. (2017). Reduced gas accretion on super-Earths and ice giants. Astronomy and Astrophysics , 606 , A146. doi: 10.1051/0004-6361/ 201731014 . arXiv:1708.00767 .  \nLaskar, J., & Robutel, P. (1993). The chaotic obliquity of the planets. Nature , 361 , 608-612. doi: 10.1038/361608a0 .  \nLevison, H. F., Kretke, K. A., & Duncan, M. J. (2015a). Growing the gasgiant planets by the gradual accumulation of pebbles. Nature , 524 , 322-324. doi: 10.1038/nature14675 .  \nLevison, H. F., Kretke, K. A., Walsh, K. J., & Bottke, W. F. (2015b). Growing the terrestrial planets from the gradual accumulation of sub-meter sized objects. Proceedings of the National Academy of Science , 112 , 14180-14185. doi: 10.1073/pnas.1513364112 . arXiv:1510.02095 .  \nZhang, K., Blake, G. A., & Bergin, E. A. (2015). Evidence of Fast Pebble Growth Near Condensation Fronts in the HL Tau Protoplanetary Disk. Astrophysical Journal Letters , 806 , L7. doi: 10.1088/2041-8205/806/1/L7 . arXiv:1505.00882 .", 'Acknowledgements': 'This work has been financially supported by the Netherlands Organisation for Scientific Research (NWO project no. ALWGO/15-01 and VIDI project 639.042.422)  \nTable 1: Relation between the dimensionfull (Dim.full) and dimensionless (Dim.less) quantities', 'Asteroid Lightcurve Data': 'Diameter, latitude and period are obtained from the Asteroid Lightcurve Data Base (LCDB) V2.0 (Warner et al., 2009). The density is taken from Baer et al. (2012). If unknown ρ = 2 g cm -3 is used.  \n·  \nFigure 1: Spin rate and orientation of planets (Earth, Mars, Uranus, and Neptune indicated with E, M, U, N, respectively) and asteroids of diameter D larger than 150 km (circles) . The spin frequency in terms of the critical breakup frequency ω crit = √ GM p /R 3 is plotted against the orientation of the spin vector to the ecliptic, where β measures ecliptic latitude. For asteroids, the symbol size scales with the diameter. The giant impact scenario predicts a spin magnitude close to breakup and an isotropic orientation of spin vectors (Kokubo & Ida, 2007) (sin β , which is proportional to the x-axis scaling, would then be uniformly distributed). The lower horizontal lines indicate the maximum retrograde and prograde spin rates obtained from systematic planetesimal accretion (Dones & Tremaine, 1993a). Planetary spins are too large for planetesimal accretion and appear too much skewed towards β = 90 · for the giant impact scenario. See Appendix A for the asteroid properties.  \n<!-- image -->  \nFigure 2: A sketch of the shearing sheet domain co-rotating with a protoplanet (center). Pebbles enter from the top at the right side from ( x 0 , y 0 ) due to a negative shear velocity and the downwards pointing headwind. The flux of pebbles increases to the right due to the increasing shear velocity. The highlighted pebble (red dotted) supplies its angular momentum to the protoplanet at impact.  \n<!-- image -->  \nFigure 3: Spin rotation rates obtained from pebble accretion for the planar case (top) and 3D case (bottom) . For each combination of stopping time τ s and headwind-to-shear parameter ζ w the average vertical angular momentum per unit mass 〈 l z 〉 transfered by pebbles to the protoplanet is calculated. This is converted into a normalized spin through the relation ω/ω crit = 〈 l z 〉 /n √ GMR where an inertia factor of n = 0 . 4 has been used. Valid for α p = 0 . 002. The transformation from ζ w to R assumes a disk headwind of v hw = 30 m/s.  \n<!-- image -->  \n10  \n2  \n(b) 3D  \nR  \nPA  \nR [km]  \n10  \n3  \n10  \ns  \n10  \n10  \n0  \n10  \n1  \nw  \n1  \n2  \n2.0  \n1.5  \n1.0  \n0.5  \n0.0  \nt  \ni  \nr  \nc  \n/  \nz  \n0.5  \n1.0  \n1.5  \n2.0  \nFigure 4: Preference of spin asymmetry explained by geometry of pebble accretion . Pebble accreting trajectories (for τ s = 0 . 1 and α p = 3 × 10 -3 are plotted in the frame of the circularly-moving planet (center). In this frame pebbles approach the planet at an angle, due to their radially drift motion. The distance between two adjacent dots on the trajectory indicate the same amount of time. (a) Due to the slanted approach pebbles encountering the planet (black dot) on an exterior orbit (red) will be scattered at a steeper angle with respect to the headwind than those pebbles that approach on the interior side (blue). Consequently, these pebbles face a stronger aerodynamic deceleration and the prolonged encounter duration allows their gravitational capture, whereas the more horizontally-deflected pebbles (blue) escape. In (b) the situation for an accretion range of the order of the planet Hill radius R H (dashed circle) is illustrated, where the red and blue curves indicate the extreme trajectories that lead to capture. As the shear (arrows) is significant for larger planets, pebbles on interior (blue) trajectories approach at a lower velocity, resulting in longer encounters and a higher total capture rate by the planet. In both cases the asymmetry results in a significant net prograde contribution.  \n<!-- image -->  \nFigure 5: Spin results obtained in dimensionfull and dimensionless units . A : The specific period for τ s = 0 . 1 for planetesimal/protoplanet sizes R ∈ [10 km , 7000 km] and a protoplanet density of 1 g cm -3 . The dashed vertical line shows the protoplanet radius where pebble accretion starts (Equation 27). B : the same data in terms of dimensionless quantities ζ w vs the fraction of the break-up frequency.  \n<!-- image -->  \n<!-- image -->  \n]  \nR  \n[  \nY  \n]  \nR  \n[  \nY  \nx  \nco  \nx  \nX  \n[  \nR  \n]  \nco  \nX  \n[  \nR  \n]  \ne1  \nx  \n0 [R]  \nFigure 6: Pebble trajectories (left) and corresponding spin curves (right) . They are selected from the α p = 3 × 10 -3 integrations and τ s = 0 . 1 for indicated protoplanet radii. The spin curves show individual SAM contributions to the protoplanet versus initial release distance x 0 . x 0 ranges from the interior side ( x 0 = x 1 ) to the exterior side ( x 0 = x 2 ) of the collision cross-section. The resulting mean spin 〈 l z 〉 is indicated with the dashed blue line. The extent of the y-axis ranges from ± l z , esc . Shown pebble trajectories consist of two misses (gray trajectories) and a range of impacts. The impacting trajectories correspond from black (interior) to blue (exterior) with the indicated dots in the corresponding spin curve to clarify individual spin supply behavior. Additionally these selected trajectories are uniquely labeled for reference in the main text. For D , and E , the Hill sphere is indicated by the dashed circle and the co-rotation line with the red vertical line.  \n<!-- image -->  \n+  \nl  \nz, esc  \nlz  \n=0  \nl  \nz, esc  \ne0  \ne4  \ne3  \n(e)  R = 3000 km  \n(d)  R = 1000 km  \nlz  \ne2  \nFigure 7: Clarifying illustration of a secondary Hill sphere encounter of a pebble. At first approach the pebble escapes from the Hill sphere and drifts to the co-rotation line at which the shear + headwind (gray arrows) are zero. The headwind velocity is indicated by the black arrow. As the pebble approaches the Hill sphere for a second time, the pebble is pulled inwards by the protoplanet gravity. During the path from -R Hill to the protoplanet surface, the pebble is simultaneously accelerated downwards due to the negative shear velocity and headwind velocity. As a result the pebble is forced to impact counterclockwise below the protoplanet delivering prograde spin rotation.  \n<!-- image -->  \nFigure 8: The degree of asymmetry δ vs the planet radius for both 2D and 3D mean spin outcomes. The relevant accretion regimes are subdivided by the vertical black lines.  \n<!-- image -->  \nFigure 9: Comparison of the 2D results (gray curves) with 3D results (blue curves) for the example model. A : The physical quantities protoplanet radius vs specific spin period. The blue dashed vertical line indicates R PA . B : the same data in terms of dimensionless quantities ζ w vs the fraction of the break-up frequency.  \n<!-- image -->  \n<!-- image -->  \nglyph[negationslash]  \nFigure 10: Heat maps of initial vertical and x-release coordinates for the example model, τ s = 0 . 1. The color bars indicate the sign and magnitude of SAM an individual pebble imparts to the protoplanet. The vertical solid lines with unique labels indicate that the trajectory type is similar to the corresponding 2D trajectory types in Fig. 6 for z 0 = 0 and fixed release distance x 0 . For A , 200 km; the contributions are symmetric and sum up to ∼ 0. For B , 400 km; the retrogade mean spin is slightly dampened with respect to the analogous 2D case. C , 1000 km; the Keplerian shear is equally important in the vertical release distance z 0 , resulting in additional prograde contributions interior to the collision cross-section and dampening in retrogade contributions exterior, explaining the deformation of the elliptic shape to an egg-like shape. D , 1500 km; the co-rotation line interior to the protoplanet Hill sphere causes long lasting encounters and pebbles can settle vertically to the protoplanet in a timely fashion. Conversely, the exterior pebbles fail to be captured for increasing z 0 due to increasingly weaker protoplanet gravity, leading to the rocket-like shape.Figure 11: Mean dimensionless spin normalized according to Dones & Tremaine (1993a) plotted vs Stokes number. For the convergence a fixed α p = 0 . 003 and ζ w = 1 are adopted, corresponding to the zero dispersion ( i = e = 0) and high gravity regime in the gas-free case ( τ s > 100) in Dones & Tremaine (1993a).  \n<!-- image -->  \nFigure 12: Global results vs local results. The shear-to-headwind parameter (x-axis) vs the fraction of the break-up speed (y-axis) for α p = 3 × 10 -3 and τ s = 0 . 1. The grey curve indicate results obtained from the local shearing sheet approximation , while the blue curve is obtained from global simulations performed with the same setup as Liu & Ormel (2018).  \n<!-- image -->'}