The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies 1 . How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem 2 : boulders stick together poorly 3 , and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas 4 . Gravitational collapse of the solid component has been suggested to overcome this barrier 1,5,6 . Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer 2, 7 , but turbulence must be present to explain observed gas accretion in protostellar discs 8 . Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas 9 , and these concentra-1 arXiv:0708.3890v1 [astro-ph] 29 Aug 2007 tions are augmented a further order of magnitude by a streaming instability [10][11][12] driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.Planet formation models typically treat turbulence as a diffusive process that opposes the gravitational sedimentation of solids to a high density midplane layer in circumstellar discs 7,13 . Recent models of solids moving in turbulent gas reveal that the turbulent motions not only mix them, but also concentrate metre-sized boulders in the transient gas overdensities 9 formed in magnetorotational turbulence 14 , in giant gaseous vortices 15,16 , and in spiral arms of self-gravitating discs 17 .Short-lived eddies at the dissipation scale of forced turbulence concentrate smaller millimetre-sized solids 18 . Some simulations mentioned above 9,11,12 were performed with the Pencil Code, which solves the magnetohydrodynamic (MHD) equations on a three-dimensional grid for a gas that interacts through drag forces with boulders. Boulders are represented as superparticles with independent positions and velocities, each having the mass of a huge number of boulders but the aerodynamic behaviour of a single boulder. We have now further developed the Pencil Code to include a fully parallel solver for the gravitational potential of the particles (see Supplementary Information). The particle density is mapped on the grid using the Triangular Shaped Cloud assignment scheme 19 and the gravitational potential of the solids is found using a Fast Fourier Transform method 20 .
2This allows us, for the first time, to simulate the dynamics of self-gravitating solid particles in magnetised, three-dimensional turbulence.We model a corotating, local box with linearised Keplerian shear that straddles the protoplanetary disc midplane and orbits the young star ...