K2-19 hosts a planetary system composed of two outer planets, b and c, with size of $7.0\pm 0.2~R_\oplus$ and $4.1\pm0.2~R_\oplus$ , and an inner planet, d, with a radius of $1.11\pm 0.05 R_\oplus$. A recent analysis of Transit-Timing Variations (TTVs) suggested b and c are close to but not in 3:2 mean motion resonance (MMR) because the classical resonant angles circulate. Such an architecture challenges our understanding of planet formation. Indeed, planet migration through the protoplanetary disk should lead to a capture into the MMR. Here, we show that the planets are in fact, locked into the 3:2 resonance despite circulation of the conventional resonant angles and aligned periapses. However, we show that such an orbital configuration cannot be maintained for more than a few hundred million years due to the tidal dissipation experienced by planet d. The tidal dissipation remains efficient because of a secular forcing of the innermost planet eccentricity by planets b and c. While the observations strongly rule out an orbital solution where the three planets are on close to circular orbits, it remains possible that a fourth planet is affecting the TTVs such that the four planet system is consistent with the tidal constraints.

The dynamical stability of tightly packed exoplanetary systems remains poorly understood. While for a two-planet system a sharp stability boundary exists, numerical simulations of three and more planet systems show that they can experience instability on timescales up to billions of years. Moreover, an exponential trend between the planet orbital separation measured in units of Hill radii and the survival time has been reported. While these findings have been observed in numerous numerical simulations, little is known of the actual mechanism leading to instability. Contrary to a constant diffusion process, planetary systems seem to remain dynamically quiescent for most of their lifetime before a very short unstable phase. In this work, we show how the slow chaotic diffusion due to the overlap of three-body resonances dominates the timescale leading to the instability for initially coplanar and circular orbits. While the last instability phase is related to scattering due to two-planet mean motion resonances (MMR), for circular orbits the two-planets MMR are too far separated to destabilize systems initially away from them. We develop an analytical model to generalize the empirical trend obtained for equal mass and equally-spaced planets to general systems. We obtain an analytical estimate of the survival time consistent with simulations over four orders of magnitude for the planet to star mass ratio $\epsilon$, and 6 to 8 orders of magnitude for the instability time. We also confirm that measuring the orbital spacing in terms of Hill radii is not adapted and that the right spacing unit scales as $\epsilon^{1/4}$. We predict that beyond a certain spacing, the three-planet resonances are not overlapped, which results in an increase of the survival time. We finally discuss the extension of our result to more general systems, containing more planets on initially non circular orbits.

Single-binary scattering may lead to an exchange where the single object captures a component of the binary, forming a new binary. This has been well studied in encounters between a star--planet pair and a single star. Here we explore the application of the exchange mechanism to a planet--satellite pair and another planet in the gravitational potential of a central star. As a case study, we focus on encounters between a satellite-bearing object and Neptune. We investigate whether Neptune can capture satellites from that object and if the captured satellites have orbits analogous to the Neptunian moons Triton and Nereid. Using $N$-body simulations, we study the capture probability at different encounter distances. Post-capture, we use a simple analytical argument to estimate how the captured orbits evolve under collisional and tidal effects. We find that the average capture probability reaches $\sim$$10\%$ if Neptune penetrates the donor planet's satellite system. Most moons grabbed by Neptune acquire highly eccentric orbits. Post-capture, around half of those captured, especially those on tight orbits, can be circularised, either by tides only or by collisions+tides, turning into Triton-like objects. Captures further out, on the other hand, stay on wide and eccentric orbits like that of Nereid. Both moon types can be captured in the same encounter and they have wide distributions in orbital inclination. Therefore, Triton naturally has a $\sim$50\% chance of being retrograde. A similar process potentially applies to an exoplanetary system, and our model predicts that exomoons can jump from one planet to another during planetary scattering. Specifically, there should be two distinct populations of captured moons: one on close-in circular orbits and the other on far-out eccentric orbits. The two populations may have highly inclined prograde or retrograde orbits.

The first stages of planet formation take place in protoplanetary disks that are largely made up of gas. Understanding how the gas affects planetesimals in the protoplanetary disk is therefore essential. In this paper, we discuss whether or not gas flow can erode planetesimals. We estimate how much shear stress is exerted onto the planetesimal surface by the gas as a function of disk and planetesimal properties. To determine whether erosion can take place, we compare this with previous measurements of the critical stress that a pebble-pile planetesimal can withstand before erosion begins. If erosion takes place, we estimate the erosion time of the affected planetesimals. We also illustrate our estimates with two-dimensional numerical simulations of flows around planetesimals using the lattice Boltzmann method. We find that the wall shear stress can overcome the critical stress of planetesimals in an eccentric orbit within the innermost regions of the disk. The high eccentricities needed to reach erosive stresses could be the result of shepherding by migrating planets. We also find that if a planetesimal erodes, it does so on short timescales. For planetesimals residing outside of $1 \ \rm{au}$, we find that they are mainly safe from erosion, even in the case of highly eccentric orbits.

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$).

Recent surveys have revealed that protoplanetary disks typically have dust masses that appear to be insufficient to account for the high occurrence rate of exoplanet systems. We demonstrate that this observed dust depletion is consistent with the radial drift of pebbles. Using a Monte Carlo method we simulate the evolution of a cluster of protoplanetary disks, using a 1D numerical method to viscously evolve each gas disk together with the radial drift of dust particles that have grown to 100 $\mu$m in size. For a 2 Myr old cluster of stars, we find a slightly sub-linear scaling between the gas disk mass and the gas accretion rate ($M_\mathrm{g}\propto\dot{M}^{0.9}$). However, for the dust mass we find that evolved dust disks have a much weaker scaling with the gas accretion rate, with the precise scaling depending on the age at which the cluster is sampled and the intrinsic age spread of the disks in the cluster. Ultimately, we find that the dust mass present in protoplanetary disk is on the order of 10-100 Earth masses in 1-3 Myr old star-forming regions, a factor of 10 to 100 depleted from the original dust budget. As the dust drains from the outer disk, pebbles pile up in the inner disk and locally increase the dust-to-gas ratio by a factor of up to 4 above the initial value. In these high dust-to-gas ratio regions we find conditions that are favourable for planetesimal formation via the streaming instability and subsequent growth by pebble accretion. We also find the following scaling relations with stellar mass within a 1-2 Myr old cluster: a slightly super-linear scaling between the gas accretion rate and stellar mass ($\dot{M}\propto M_\star^{1.4}$), a slightly super-linear scaling between the gas disk mass and the stellar mass ($M_\mathrm{g}\propto M_\star^{1.4}$) and a super-linear relation between the dust disk mass and stellar mass ($M_\mathrm{d}\propto M_\star^{1.4-4.1}$).

We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is $30\rm~ms^{-1}$. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula (MMSN) regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset ($<10^5$ yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. (Abridged)

Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting.

The New Horizons flyby of the cold classical Kuiper Belt object MU69 showed it to be a contact binary. The existence of other contact binaries in the 1-10km range raises the question of how common these bodies are and how they evolved into contact. Here we consider that the lobes of MU69 formed as a binary in the Solar nebula, and calculate its orbital evolution in the presence of gas. We find that the sub-Keplerian wind of the disk brings the drag timescales for 10km bodies under 1 Myr for quadratic-velocity drag. In the Kuiper belt, however, the drag is linear with velocity and the effect of the wind cancels out as the angular momentum gained in half an orbit is lost in the other half: the drag timescales for 10 km bodies remain over 10 Myr. In this situation we find that a combination of drag and Kozai-Lidov oscillations is a promising channel for collapse. We analytically solve the hierarchical three-body problem with drag and implement it into a Kozai plus tidal friction model. The permanent quadrupoles of the lobes make the Kozai oscillations stochastic, and as drag shrinks the semimajor axis it more easily allows the fluctuations to bring the system into contact. Evolution to contact happens very rapidly (within $10^4$ yr) in the pure, double-average quadrupole, Kozai region between $\approx 85-95^\circ$, and within 3 Myr in the drag-assisted region beyond it. The synergy between $J_2$ and gas drag widens the window of contact to 80-100$^\circ$ initial inclination, over a larger range of semimajor axes. As such, the model predicts an initial contact binary fraction of about 10\% for the cold classicals in the Kuiper belt. The speed at contact deviates from the escape velocity only because of the oblateness. For MU69, the contact velocity should be 3.3-4.2 m/s, in line with the observational evidence from the lack of deformation features and estimate of the tensile strength.

The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${\sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1\%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${\sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.

Nearly-axisymmetric gaps and rings are commonly observed in protoplanetary discs. The leading theory regarding the origin of these patterns is that they are due to dust trapping at the edges of gas gaps induced by the gravitational torques from embedded planets. If the concentration of solids at the gap edges becomes high enough, it could potentially result in planetesimal formation by the streaming instability. We test this hypothesis by performing global 1-D simulations of dust evolution and planetesimal formation in a protoplanetary disc that is perturbed by multiple planets. We explore different combinations of particle sizes, disc parameters, and planetary masses, and find that planetesimals form in all these cases. We also compare the spatial distribution of pebbles from our simulations with protoplanetary disc observations. Planets larger than one pebble isolation mass catch drifting pebbles efficiently at the edge of their gas gaps, and depending on the efficiency of planetesimal formation at the gap edges, the protoplanetary disc transforms within a few 100,000 years to either a transition disc with a large inner hole devoid of dust or to a disc with narrow bright rings. For simulations with planetary masses lower than the pebble isolation mass, the outcome is a disc with a series of weak ring patterns but no strong depletion between the rings. Lowering the pebble size artificially to 100 micrometer-sized "silt", we find that regions between planets get depleted of their pebble mass on a longer time-scale of up to 0.5 million years. These simulations also produce fewer planetesimals than in the nominal model with millimeter-sized particles and always have at least two rings of pebbles still visible after 1 Myr.

The four major satellites of Jupiter, known as the Galilean moons, and Saturn's most massive satellite, Titan, are believed to have formed in a predominantly gaseous circum-planetary disk, during the last stages of formation of their parent planet. Pebbles from the protoplanetary disk are blocked from flowing into the circumplanetary disk by the positive pressure gradient at the outer edge of the planetary gap, so the gas drag assisted capture of planetesimals should be the main contributor to the delivery of solids onto circum-planetary disks. However, a consistent framework for the subsequent accretion of the moons remains to be built. Here we use numerical integrations to show that most planetesimals being captured within a circum-planetary disk are strongly ablated due to the frictional heating they experience, thus supplying the disk with small dust grains, whereas only a small fraction 'survives' their capture. We then construct a simple model of a circum-planetary disk supplied by ablation, where the flux of solids through the disk is at equilibrium with the ablation supply rate, and investigate the formation of moons in such disks. We show that the growth of satellites is driven mainly by accretion of the pebbles that coagulate from the ablated material. The pebble-accreting protosatellites rapidly migrate inward and pile up in resonant chains at the inner edge of the circum-planetary disk. We propose that dynamical instabilities in these resonant chains are at the origin of the different architectures of Jupiter's and Saturn's moon systems. The assembly of moon systems through pebble accretion can therefore be seen as a down-scaled manifestation of the same process that forms systems of super-Earths and terrestrial-mass planets around solar-type stars and M-dwarfs.

The formation of cold gas giants similar to Jupiter and Saturn in orbit and mass is a great challenge for planetesimal-driven core accretion models because the core growth rates far from the star are low. Here we model the growth and migration of single protoplanets that accrete planetesimals and gas. We integrated the core growth rate using fits in the literature to $N$-body simulations, which provide the efficiency of accreting the planetesimals that a protoplanet migrates through. We take into account three constraints from the solar system and from protoplanetary discs: (1) the masses of the terrestrial planets and the comet reservoirs in Neptune's scattered disc and the Oort cloud are consistent with a primordial planetesimal population of a few Earth masses per AU, (2) evidence from the asteroid belt and the Kuiper belt indicates that the characteristic planetesimal diameter is 100 km, and (3) observations of protoplanetary discs indicate that the dust is stirred by weak turbulence; this gas turbulence also excites the inclinations of planetesimals. Our nominal model built on these constraints results in maximum protoplanet masses of $0.1$ Earth masses. Ignoring constraint (1) above, we show that even a planetesimal population of 1,000 Earth masses, corresponding to 50 Earth masses per AU, fails to produce cold gas giants (although it successfully forms hot and warm gas giants). We conclude that a massive planetesimal reservoir is in itself insufficient to produce cold gas giants. The formation of cold gas giants by planetesimal accretion additionally requires that planetesimals are small and that the turbulent stirring is very weak, thereby violating all three above constraints.

Jupiter Trojans are are characterized by dark photometric colors, high inclinations and an asymmetry in number of bodies between the two swarms. Different models have been proposed to explain the high inclination of the Trojans and to interpret their relation with the TNOs, but none of them can also satisfactorily explain the asymmetry ratio. Recently it has been found that the asymmetry can arise if Jupiter has migrated inwards by at least a few au during its growth. The asymmetry and the dark colors of the Trojans are natural outcomes of this model, but simulations with massless unperturbed disc particles led to a flat distribution of the Trojan inclinations and a final total mass that was 3-4 orders of magnitude larger than the current one. In our work, we investigate the possible origin of the peculiar inclination distribution of the Trojans in the scenario where Jupiter migrates inwards. We analyze different possibilities: (a) the secular evolution of an initially flat Trojan population, (b) the presence of planetary embryos among the Trojans and (c) capture of the Trojans from a pre-stirred planetesimal population. We find that the secular evolution of the Trojans and Saturn do not affect the Trojan inclinations appreciably, nor is there any significant mass depletion. Embryos embedded in the swarms, in contrast, can stir the Trojan inclinations and can also deplete the swarms efficiently, but it turns out that it is very difficult to get rid of all of the massive bodies. We propose that the disc where Jupiter's core was forming was already excited by the presence of other embryos competing in Jupiter's core's feeding zone. We show that the trapped Trojans preserve their high inclination through the gas phase of the disc and that Saturn's perturbations are more effective on highly inclined Trojans, leading to a lower capture efficiency and to a substantial depletion of the swarms.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a number of these parameters affect that measured accretion rate of a Saturn-mass planet. We present a global 3D radiative hydrodynamics framework using the FARGOCA-code. The planet is represented by a gravitational potential with a smoothing length at the location of the planet. No mass or energy sink is used, instead luminosity and gas accretion rates are self-consistently computed. We find that the gravitational smoothing length must be resolved by at least 10 grid cells to obtain converged measurements of the gas accretion rates. Secondly, we find gas accretion rates into planetary envelopes compatible with previous studies, and continue to explain those via the structure of our planetary envelopes and their luminosities. Our measured gas accretion rates are formally in the stage of Kelvin-Helmholtz contraction due to the modest entropy loss that can be obtained over the simulation time-scale, but our accretion rates are compatible with those expected during late run-away accretion. Our detailed simulations of the gas flow into the envelope of a Saturn-mass planet provides a framework for understanding the general problem of gas accretion during planet formation and highlights circulation features that develop inside the planetary envelopes. Those circulation features feedback into the envelope energetics and can have further implications for transporting dust into the inner regions of the envelope.

Solid particles in protoplanetary discs can grow by direct vapour deposition outside of ice lines. The presence of microscopic silicate particles may nevertheless hinder growth into large pebbles, since the available vapour is deposited predominantly on the small grains that dominate the total surface area. Experiments on heterogeneous ice nucleation, performed to understand ice clouds in the Martian atmosphere, show that the formation of a new ice layer on a silicate surface requires a substantially higher water vapour pressure than the deposition of water vapour on an existing ice surface. In this paper, we investigate how the difference in partial vapour pressure needed for deposition of vapour on water ice versus heterogeneous ice nucleation on silicate grains influences particle growth close to the water ice line. We developed and tested a dynamical 1D deposition and sublimation model, where we include radial drift, sedimentation, and diffusion in a turbulent protoplanetary disc. We find that vapour is deposited predominantly on already ice-covered particles, since the vapour pressure exterior of the ice line is too low for heterogeneous nucleation on bare silicate grains. Icy particles can thus grow to centimetre-sized pebbles in a narrow region around the ice line, whereas silicate particles stay dust-sized and diffuse out over the disc. The inhibition of heterogeneous ice nucleation results in a preferential region for growth into planetesimals close to the ice line where we find large icy pebbles. The suppression of heterogeneous ice nucleation on silicate grains may also be the mechanism behind some of the observed dark rings around ice lines in protoplanetary discs, as the presence of large ice pebbles outside ice lines leads to a decrease in the opacity there.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

We develop a pebble-driven model to study the formation and evolution of planets around stars in the mass range of 0.08 and 1 solar mass. The growth and migration of a large number of individual protoplanetary embryos are simulated in a population synthesis manner. We test two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplanetary disks. Two types of disks with different turbulent viscous parameters alpha of 1e-3 and 1e-4 are investigated, to shed light on the role of outward migration of protoplanets. The forming planets are compared with the observed exoplanets in terms of masses, semimajor axes, metallicities, and water contents. We find that gas giant planets are likely to form when the characteristic disk sizes are larger, the disk accretion rates are higher, the disks are more metal-rich and/or their stellar hosts are more massive. Our model shows that 1) the characteristic mass of super-Earth is set by the pebble isolation mass. Super-Earth masses increase linearly with the mass of its stellar host, corresponding to one Earth mass around a late M-dwarf star and 20 Earth masses around a solar-mass star. 2) The low-mass planets up to 20 Earth masses can form around stars with a wide range of metallicities, while massive gas giant planets are preferred to grow around metal-rich stars. 3) Super-Earth planets that are mainly composed of silicates, with relatively low water fractions can form from protoplanetary embryos at the water ice line in weakly turbulent disks where outward migration is suppressed. However, if the embryos are formed over a wide range of radial distances, the super-Earths would end up having a distinctive, bimodal composition in water mass. Our model succeeds in quantitatively reproducing several important observed properties of exoplanets and correlations with their stellar hosts.

Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r<1AU) and form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explore slower migration rates which allow the formation of gas giants from embryos originating from the 5-10AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. Giant planet growth requires a sufficiently-high pebble flux to enable growth to out-compete migration. Even higher pebble fluxes produce systems with multiple gas giants. We show that planetary embryos starting interior to 5AU do not grow into gas giants, even if migration is slow and the pebble flux is large. Instead they grow to the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is independent of the pebble flux. Additionally we show that the long term evolution of our formed planetary systems can produce systems with hot super-Earths and outer gas giants as well as systems of giants on eccentric orbits (abridged).

Streaming instability is a key mechanism in planet formation, clustering pebbles into planetesimals. It is triggered at a particular disk location where the local volume density of solids exceeds that of the gas. After their formation, planetesimals can grow by feeding from other planetesimals in the birth ring as well as by accreting inwardly drifting pebbles from the outer disk. To investigate the growth of planetesimals at a single location by the streaming instability, we test the conditions under which super-Earths are able to form within the lifetime of the gaseous disk. We modify the \texttt{Mercury} N-body code to trace the growth and dynamical evolution of a swarm of planetesimals at the ice line for a solar-mass star. Three distributions of planetesimal sizes are investigated: (i) a mono-dispersed population of 400 km radius planetesimals, (ii) a poly-dispersed populations of planetesimals from 200 km up to 1000 km, (iii) a bimodal distribution with a single runaway body and a swarm of smaller, 100 km size planetesimals. The mono-disperse population of 400 km size planetesimals cannot form $\gtrsim$ Earth mass protoplanets. Their velocity dispersions are quickly excited, which suppresses both planetesimal and pebble accretion. Planets can form from the poly-dispersed and bimodal distributions. In these circumstances, the two-component nature damps the random velocity of the large embryo by small planetesimals' dynamical friction, allowing the embryo to accrete pebbles efficiently when it approaches $10^{-2}$ Earth mass. We find that super-Earth planets are preferred to form when the pebble mass flux is higher, the disk turbulence is lower, or the Stokes number of the pebbles is higher.

At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes ($\lesssim15M_{\oplus}$) and super-massive planetary cores potentially able to become gas giant planets ($\gtrsim15M_{\oplus}$). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk's inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to $\sim$95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction ($\gtrsim95\%$) of unstable systems with a sprinkling ($\lesssim5\%$) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (${\rm N\geq2}$) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)

Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.

Pebble accretion is an efficient mechanism able to build up the core of the giant planets within the lifetime of the protoplanetary disc gas-phase. The core grows via this process until the protoplanet reaches its pebble isolation mass and starts to accrete gas. During the growth, the protoplanet undergoes a rapid, large-scale, inward migration due to the interactions with the gaseous protoplanetary disc. In our work, we investigate how this early migration would have affected the minor body populations in our solar system. In particular, we focus on the Jupiter Trojans and the Hildas asteroids. We found that a massive and eccentric Hilda group is captured during the migration from a region between 5 and 8 au and subsequently depleted during the late instability of the giant planets. Our simulations also show that inward migration of the giant planets always produces a Jupiter Trojans' leading swarm more populated than the trailing one, with a ratio comparable to the current observed Trojan asymmetry ratio. The in situ formation of Jupiter, on the other hand, produces symmetric leading/trailing swarms. The reason for the asymmetry is the relative drift between the migrating planet and the particles in the coorbital resonance. The capture happens during the growth of Jupiter's core and Trojan asteroids are afterwards carried along during the giant planet's migration to their final orbits. The asymmetry and eccentricity of the captured Trojans correspond well to observations, but their inclinations are near zero and their total mass is 3-4 orders of magnitude higher than the current population. Future modelling will be needed to understand whether the dynamical evolution of the Trojans over billions of years will raise the inclinations and deplete the masses to observed values.

Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from the star) that become increasingly horizontal (migration-dominated) with rising mass of the protoplanet. It has been shown recently that the migration torque on a protoplanet is reduced proportional to the relative height of the gas gap carved by the growing planet. Here we show from 1-D simulations of planet-disc interaction that the mass at which a planet carves a 50% gap is approximately 2.3 times the pebble isolation mass. Our measurements of the pebble isolation mass from 1-D simulations match published 3-D results relatively well, except at very low viscosities where the 3-D pebble isolation mass is significantly higher, possibly due to gap edge instabilities not captured in 1-D. The pebble isolation mass demarks the transition from pebble accretion to gas accretion. Gas accretion to form gas-giant planets therefore takes place over a few astronomical units of migration after reaching first the pebble isolation mass and, shortly after, the 50% gap mass. Our results demonstrate how planetary growth can outperform migration, both during core accretion and during gas accretion, even when the Stokes number of the pebbles is small, St~0.01, and the pebble-to-gas flux ratio in the protoplanetary disc is in the nominal range of 0.01-0.02. We find that planetary growth is very rapid in the first million years of the protoplanetary disc and that the probability for forming gas-giant planets increases with the initial size of the protoplanetary disc and with decreasing turbulent diffusion.

Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability (SI) is one of the most promising. This work tests the physical assumptions of models that lead to SI in protoplanetary disks (PPDs). We conduct laboratory experiments in which we track the three-dimensional motion of spherical solid particles fluidized in a low-pressure, laminar, incompressible, gas stream. The particle sizes span the Stokes-Epstein drag regime transition and the overall dust-to-gas mass density ratio is close to unity. Lambrechts et al. (2016) established the similarity of the laboratory flow to a simplified PPD model flow. We obtain experimental results suggesting an instability due to particle-gas interaction: i) there exist variations in particle concentration in the direction of the mean drag forces; ii) the particles have a tendency to 'catch up' to one another when they are in proximity; iii) particle clumping occurs on very small scales, which implies local enhancements above the background dust-to-gas mass density ratio by factors of several tens; v) the presence of these density enhancements occurs for a mean dust-to-gas mass density ratio approaching or greater than 1; v) we find evidence for collective particle drag reduction when the local particle number density becomes high and when the background gas pressure is high so that the drag is in the continuum regime. The experiments presented here are precedent-setting for observing SI under controlled conditions and may lead to a deeper understanding of how it operates in nature.

The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MHD) effects meets this requirement. However, simulations have shown that the midplane of a dead zone is not completely quiescent. In order to examine the concentration of solids in such an environment, we use the local-shearing-box approximation to simulate a particle-gas system with an Ohmic dead zone including mutual drag force between the gas and the solids. We systematically compare the evolution of the system with ideal or non-ideal MHD, with or without back-reaction drag force from particles on gas, and with varying solid abundances. Similar to previous investigations of dead zone dynamics, we find that particles of dimensionless stopping time $\tau_s=0.1$ do not sediment appreciably more than those in ideal magneto-rotational turbulence, resulting in a vertical scale height an order of magnitude larger than in a laminar disk. Contrary to the expectation that this should curb the formation of planetesimals, we nevertheless find that strong clumping of solids still occurs in the dead zone when solid abundances are similar to the critical value for a laminar environment. This can be explained by the weak radial diffusion of particles near the mid-plane. The results imply that the sedimentation of particles to the mid-plane is not a necessary criterion for the formation of planetesimals by the streaming instability.

The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-jupiter-mass planets are piled up at > 1 au, while jupiter/sub-jupiter-mass planets are broadly distributed from ~0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al. (2018), by comparing the migration, growth timescales of gas giants, and disk lifetime and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. (2018) proposed that type II migration speed is no other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent MHD simulations. Population synthesis simulations show that super-jupiter-mass planets remain at > 1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at > 1 au would be theoretically solved, if the new formula is confirmed and wind-driven disk accretion dominates.

The radial drift and diffusion of dust particles in protoplanetary disks affect both the opacity and temperature of such disks as well as the location and timing of planetesimal formation. In this paper, we present results of numerical simulations of particle-gas dynamics in protoplanetary disks that include dust grains with various size distributions. We consider three scenarios in terms of particle size ranges, one where the Stokes number $\tau_{\rm{s}} = 10^{-1} - 10^0$, one where $\tau_{\rm{s}} = 10^{-4} - 10^{-1}$ and finally one where $\tau_{\rm{s}} = 10^{-3} - 10^{0}$. Moreover, we consider both discrete and continuous distributions in particle size. In accordance with previous works we find in our multi-species simulations that different particle sizes interact via the gas and as a result their dynamics changes compared to the single-species case. The larger species trigger the streaming instability and create turbulence that drives the diffusion of the solid materials. We measure the radial equilibrium velocity of the system and find that the radial drift velocity of the large particles is reduced in the multi-species simulations and that the small particle species move on average outwards. We also vary the steepness of the size distribution, such that the exponent of the solid number density distribution, $\rm{d}\it{N}/\rm{d}\it{a} \propto \it{a^{-q}}$, is either $q = 3$ or $q = 4$. We measure the scale height of the particles and observe that small grains are stirred up well above the sedimented midplane layer where the large particles reside. Our measured diffusion and drift parameters can be used in coagulation models for planet formation as well as to understand relative mixing of the components of primitive meteorites (matrix, chondrules and CAIs) prior to inclusion in their parent bodies.

Binary stars form from the same parent molecular cloud and thus have the same chemical composition. Forming planets take building material (solids) away from the surrounding protoplanetary disc. Assuming that the disc's accretion onto the star is the main process that clears the disc, the atmosphere of the star will show abundance reductions caused by the material accreted by the forming planet(s). If planets are only forming around one star of a binary system, the planet formation process can result in abundance differences in wide binary stars, if their natal protoplanetary discs do not interact during planet formation. Abundance differences in the atmospheres of wide binaries hosting giant planets have already been observed and linked to the formation location of giant planets. Here, we model how much building material is taken away for super-Earth planets that form inside/outside of the water ice line as well as ice giants forming inside/outside of the CO ice line. Our model predicts a significant abundance difference $\Delta$[X/H] in the stellar atmospheres of the planet-hosting binary component. Our model predicts that super-Earths that form inside the water ice line ($r<r_{\rm H_2O}$) will result in an $\Delta$[Fe/H]/$\Delta$[O/H] abundance difference in the their host star that is a factor of 2 larger than for super-Earths formed outside the water ice line ($r>r_{\rm H_2O}$) in the water rich parts of the disc. Additionally, our model shows that the $\Delta$[Fe/H]/$\Delta$[C/H] abundance difference in the host star is at least a factor of 3 larger for ice giants formed at $r<r_{\rm CO}$ compared to ice giants formed far out in the protoplanetary disc ($r>r_{\rm CO}$). Future observations of wide binary star systems hosting super-Earths and ice giants could therefore help to constrain the migration pathway of these planets and thus constrain planet formation theories.

Planet-planet collisions are a common outcome of instability in systems of transiting planets close to the star, as well as occurring during in-situ formation of such planets from embryos. Previous N-body studies of instability amongst transiting planets have assumed that collisions result in perfect merging. Here, we explore the effects of implementing a more realistic collision prescription on the outcomes of instability and in-situ formation at orbital radii of a few tenths of an au. There is a strong effect on the outcome of the growth of planetary embryos, so long as the debris thrown off in collisions is rapidly removed from the system (which happens by collisional processing to dust, and then removal by radiation forces) and embryos are small ($< 0.1\mathrm{\,M}_\oplus$). If this is the case, then systems form fewer detectable ($\gtrsim1\mathrm{\,M}_\oplus$) planets than systems evolved under the assumption of perfect merging in collisions. This provides some contribution to the "Kepler Dichotomy": the observed over-abundance of single-planet systems. The effects of changing the collision prescription on unstable mature systems of super-Earths are less pronounced. Perfect mergers only account for a minority of collision outcomes in such systems, but most collisions resulting in mass loss are grazing impacts in which only a few per cent. of mass is lost. As a result, there is little impact on the final masses and multiplicities of the systems after instability when compared to systems evolved under the assumption that collisions always result in perfect merging.

Jupiter played an important role in determining the structure and configuration of the Solar System. Whereas hot-Jupiter type exoplanets preferentially form around metal-rich stars, the conditions required for the formation of planets with masses, orbits and eccentricities comparable to Jupiter (Jupiter analogues) are unknown. Using spectroscopic metallicities, we show that stars hosting Jupiter analogues have an average metallicity close to solar, in contrast to their hot-Jupiter and eccentric cool Jupiter counterparts, which orbit stars with super-solar metallicities. Furthermore, the eccentricities of Jupiter analogues increase with host star metallicity, suggesting that planet-planet scatterings producing highly eccentric cool Jupiters could be more common in metal-rich environments. To investigate a possible explanation for these metallicity trends, we compare the observations to numerical simulations, which indicate that metal-rich stars typically form multiple Jupiters, leading to planet-planet interactions and, hence, a prevalence of either eccentric cool Jupiters or hot-Jupiters with circularized orbits. Although the samples are small and exhibit variations in their metallicities, suggesting that numerous processes other than metallicity affect the formation of planetary systems, the data in hand suggests that Jupiter analogues and terrestrial-sized planets form around stars with average metallicities close to solar, whereas high metallicity systems preferentially host eccentric cool Jupiter or hot-Jupiters, indicating higher metallicity systems may not be favorable for the formation of planetary systems akin to the Solar System.

The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $\alpha$. We find that small pebbles, coupled to the gas, with Stokes number $\tau_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $\alpha$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System.

Unexpected clustering in the orbital elements of minor bodies beyond the Kuiper belt has led to speculations that our solar system actually hosts nine planets, the eight established plus a hypothetical "Planet Nine". Several recent studies have shown that a planet with a mass of about 10 Earth masses on a distant eccentric orbit with perihelion far beyond the Kuiper belt could create and maintain this clustering. The evolutionary path resulting in an orbit such as the one suggested for Planet Nine is nevertheless not easily explained. Here we investigate whether a planet scattered away from the giant-planet region could be lifted to an orbit similar to the one suggested for Planet Nine through dynamical friction with a cold, distant planetesimal belt. Recent simulations of planetesimal formation via the streaming instability suggest that planetesimals can readily form beyond 100au. We explore this circularisation by dynamical friction with a set of numerical simulations. We find that a planet that is scattered from the region close to Neptune onto an eccentric orbit has a 20-30% chance of obtaining an orbit similar to that of Planet Nine after 4.6Gyr. Our simulations also result in strong or partial clustering of the planetesimals; however, whether or not this clustering is observable depends on the location of the inner edge of the planetesimal belt. If the inner edge is located at 200au the degree of clustering amongst observable objects is significant.

Two basic routes for planetesimal formation have been proposed over the last few decades. One is a classical "slow-growth" scenario. Another one is particle concentration models, in which small pebbles are concentrated locally and then collapse gravitationally to form planetesimals. Both types of models make certain predictions for the size spectrum and internal structure of newly-born planetesimals. We use these predictions as input to simulate collisional evolution of debris discs left after the gas dispersal. The debris disc emission as a function of a system's age computed in these simulations is compared with several Spitzer and Herschel debris disc surveys around A-type stars. We confirm that the observed brightness evolution for the majority of discs can be reproduced by classical models. Further, we find that it is equally consistent with the size distribution of planetesimals predicted by particle concentration models - provided the objects are loosely bound "pebble piles" as these models also predict. Regardless of the assumed planetesimal formation mechanism, explaining the brightest debris discs in the samples uncovers a "disc mass problem." To reproduce such discs by collisional simulations, a total mass of planetesimals of up to ~1000 Earth masses is required, which exceeds the total mass of solids available in the protoplanetary progenitors of debris discs. This may indicate that stirring was delayed in some of the bright discs, that giant impacts occurred recently in some of them, that some systems may be younger than previously thought, or that non-collisional processes contribute significantly to the dust production.

The processes that led to the formation of the planetary bodies in the Solar System are still not fully understood. Using the results obtained with the comprehensive suite of instruments on-board ESA's Rosetta mission, we present evidence that comet 67P/Churyumov-Gerasimenko likely formed through the gentle gravitational collapse of a bound clump of mm-sized dust aggregates ("pebbles"), intermixed with microscopic ice particles. This formation scenario leads to a cometary make-up that is simultaneously compatible with the global porosity, homogeneity, tensile strength, thermal inertia, vertical temperature profiles, sizes and porosities of emitted dust, and the steep increase in water-vapour production rate with decreasing heliocentric distance, measured by the instruments on-board the Rosetta spacecraft and the Philae lander. Our findings suggest that the pebbles observed to be abundant in protoplanetary discs around young stars provide the building material for comets and other minor bodies.

Chondrules in primitive meteorites likely formed by recrystallisation of dust aggregates that were flash-heated to nearly complete melting. Chondrules may represent the building blocks of rocky planetesimals and protoplanets in the inner regions of protoplanetary discs, but the source of ubiquitous thermal processing of their dust aggregate precursors remains elusive. Here we demonstrate that escape of positrons released in the decay of the short-lived radionuclide $^{26}$Al leads to a large-scale charging of dense pebble structures, resulting in neutralisation by lightning discharge and flash-heating of dust and pebbles. This charging mechanism is similar to a nuclear battery where a radioactive source charges a capacitor. We show that the nuclear battery effect operates in circumplanetesimal pebble discs. The extremely high pebble densities in such discs are consistent with conditions during chondrule heating inferred from the high abundance of sodium within chondrules. The sedimented mid-plane layer of the protoplanetary disc may also be prone to charging by the emission of positrons, if the mass density of small dust there is at least an order of magnitude above the gas density. Our results imply that the decay energy of $^{26}$Al can be harvested to drive intense lightning activity in protoplanetary discs. The total energy stored in positron emission is comparable to the energy needed to melt all solids in the protoplanetary disc. The efficiency of transferring the positron energy to the electric field nevertheless depends on the relatively unknown distribution and scale-dependence of pebble density gradients in circumplanetesimal pebble discs and in the protoplanetary disc mid-plane layer.

Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time $\tau_\mathrm{s} = 10^{-2}$ and $10^{-3}$ -- representing mm- and cm-sized particles interior of the water ice line -- concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of $\tau_\mathrm{s} \ll 1$. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability.

We study the collapse of pebble clouds with a statistical model to find the internal structure of comet-sized planetesimals. Pebble-pebble collisions occur during the collapse and the outcome of these collisions affect the resulting structure of the planetesimal. We expand our previous models by allowing the individual pebble sub-clouds to contract at different rates and by including the effect of gas drag on the contraction speed and in energy dissipation. Our results yield comets that are porous pebble-piles with particle sizes varying with depth. In the surface layers there is a mixture of primordial pebbles and pebble fragments. The interior, on the other hand, consists only of primordial pebbles with a narrower size distribution, yielding higher porosity there. Our results imply that the gas in the protoplanetary disc plays an important role in determining the radial distribution of pebble sizes and porosity inside planetesimals.

We investigate the origin of the period distribution of giant planets. We fit the bias-corrected distribution of gas-giant planets inside 300 days found by Santerne et al. (2016) using a planet formation model based on pebble accretion. We investigate two possible initial conditions: a linear distribution of planetary seeds, and seeds injected exclusively on the water and CO icelines. Our simulations exclude the linear initial distribution of seeds with a high degree of confidence. Our bimodal model based on snowlines give a more reasonable fit to the data, with the discrepancies reducing significantly if we assume the water snowline to be a factor 3-10 less efficient at producing planets. This model moreover performs better on both the warm/hot Jupiters ratio and a Gaussian mixture model as comparison criteria. Our results hint that the gas-giant exoplanets population inside 300 days is more compatible with planets forming preferentially at special locations.

Atmospheric chemical abundances of giant planets lead to important constraints on planetary formation and migration. Studies have shown that giant planets that migrate through the protoplanetary disk can accrete substantial amounts of oxygen-rich planetesimals, leading to super-solar metallicities in the envelope and solar or sub-solar C/O ratios. Pebble accretion has been demonstrated to play an important role in core accretion and to have growth rates that are consistent with planetary migration. The high pebble accretion rates allow planetary cores to start their growth beyond 10 AU and subsequently migrate to cold (>~ 1 AU), warm (~0.1 AU- 1AU) or hot (<~ 0.1 AU) orbits. In this work we investigate how the formation of giant planets via pebble accretion influences their atmospheric chemical compositions. We find that under the standard pebble accretion scenario, where the core is isolated from the envelope, the resulting metallicities (O/H and C/H ratios) are sub-solar, while the C/O ratios are super-solar. Planets that migrate through the disk to become hot Jupiters accrete substantial amounts of water vapour, but still acquire slightly sub-solar O/H and super-solar C/O of 0.7-0.8. The metallicity can be substantially sub-solar (~0.2-0.5x solar) and the C/O can even approach 1.0 if the planet accretes its envelope mostly beyond the CO2 ice line, i.e. cold Jupiters or hot Jupiters that form far out and migrate in by scattering. Allowing for core erosion yields significantly super-solar metallicities and solar or sub-solar C/O, which can also be achieved by other means, e.g. photo-evaporation and late-stage planetesimal accretion.

The star EPIC 210894022 has been identified from a light curve acquired through the K2 space mission as possibly orbited by a transiting planet. Our aim is to confirm the planetary nature of the object and derive its fundamental parameters. We combine the K2 photometry with reconnaissance spectroscopy and radial velocity (RV) measurements obtained using three separate telescope and spectrograph combinations. The spectroscopic synthesis package SME has been used to derive the stellar photospheric parameters that were used as input to various stellar evolutionary tracks in order to derive the parameters of the system. The planetary transit was also validated to occur on the assumed host star through adaptive imaging and statistical analysis. The star is found to be located in the background of the Hyades cluster at a distance at least 4 times further away from Earth than the cluster itself. The spectrum and the space velocities of EPIC 210894022 strongly suggest it to be a member of the thick disk population. We find that the star is a metal poor ([Fe/H]=-0.53+/-0.05 dex) and alpha-rich somewhat evolved solar-like object of spectral type G3 with Teff=5730+/-50 K, logg=4.15+/-0.1 (cgs), radius of 1.3+/-0.1 R_Sun, and mass of 0.88+/-0.02 M_Sun. The RV detection together with the imaging confirms with a high level of significance that the transit signature is caused by a super-Earth orbiting the star EPIC 210894022. We measure a mass of 8.6+/-3.9 M_Earth and a radius of 1.9+/-0.2 R_Earth. A second more massive object with a period longer than about 120 days is indicated by a long term linear acceleration. With an age of > 10 Gyrs this system is one of the oldest where planets is hitherto detected. Further studies of this planetary system is important since it contains information about the planetary formation process during a very early epoch of the history of our Galaxy.

The distribution of exoplanet masses is not primordial. After the initial stage of planet formation is complete, the gravitational interactions between planets can lead to the physical collision of two planets, or the ejection of one or more planets from the system. When this occurs, the remaining planets are typically left in more eccentric orbits. Here we use present-day eccentricities of the observed exoplanet population to reconstruct the initial mass function of exoplanets before the onset of dynamical instability. We developed a Bayesian framework that combines data from N-body simulations with present-day observations to compute a probability distribution for the planets that were ejected or collided in the past. Integrating across the exoplanet population, we obtained an estimate of the initial mass function of exoplanets. We find that the ejected planets are primarily sub-Saturn type planets. While the present-day distribution appears to be bimodal, with peaks around $\sim 1 M_{\rm J}$ and $\sim 20 M_\oplus$, this bimodality does not seem to be primordial. Instead, planets around $\sim 60 M_\oplus$ appear to be preferentially removed by dynamical instabilities. Attempts to reproduce exoplanet populations using population synthesis codes should be mindful of the fact that the present population has been depleted of intermediate-mass planets. Future work should explore how the system architecture and multiplicity might alter our results.

Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_\oplus$) planetesimal belts beyond 100 au and up to $\sim 20 M_\oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_\oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.

We study how close-in systems such as those detected by Kepler are affected by the dynamics of bodies in the outer system. We consider two scenarios: outer systems of giant planets potentially unstable to planet--planet scattering, and wide binaries that may be capable of driving Kozai or other secular variations of outer planets' eccentricities. Dynamical excitation of planets in the outer system reduces the multiplicity of Kepler-detectable planets in the inner system in $\sim20-25\%$ of our systems. Accounting for the occurrence rates of wide-orbit planets and binary stars, $\approx18\%$ of close-in systems could be destabilised by their outer companions in this way. This provides some contribution to the apparent excess of systems with a single transiting planet compared to multiple, however, it only contributes at most $25\%$ of the excess. The effects of the outer dynamics can generate systems similar to Kepler-56 (two coplanar planets significantly misaligned with the host star) and Kepler-108 (two significantly non-coplanar planets in a binary). We also identify three pathways to the formation of eccentric warm Jupiters resulting from the interaction between outer and inner systems: direct inelastic collision between an eccentric outer and an inner planet, secular eccentricity oscillations that may "freeze out" when scattering resolves in the outer system; and scattering in the inner system followed by "uplift", where inner planets are removed by interaction with the outer planets. In these scenarios, the formation of eccentric warm Jupiters is a signature of a past history of violent dynamics among massive planets beyond $\sim1$ au.

The streaming instability is a mechanism to concentrate solid particles into overdense filaments that undergo gravitational collapse and form planetesimals. However, it remains unclear how the initial mass function of these planetesimals depends on the box dimensions of numerical simulations. To resolve this, we perform simulations of planetesimal formation with the largest box dimensions to date, allowing planetesimals to form simultaneously in multiple filaments that can only emerge within such large simulation boxes. In our simulations, planetesimals with sizes between 80 km and several hundred kilometers form. We find that a power law with a rather shallow exponential cutoff at the high-mass end represents the cumulative birth mass function better than an integrated power law. The steepness of the exponential cutoff is largely independent of box dimensions and resolution, while the exponent of the power law is not constrained at the resolutions we employ. Moreover, we find that the characteristic mass scale of the exponential cutoff correlates with the mass budget in each filament. Together with previous studies of high-resolution simulations with small box domains, our results therefore imply that the cumulative birth mass function of planetesimals is consistent with an exponentially tapered power law with a power-law exponent of approximately -1.6 and a steepness of the exponential cutoff in the range of 0.3-0.4.

The study of cosmology, galaxy formation and exoplanets has now advanced to a stage where a cosmic inventory of terrestrial planets may be attempted. By coupling semi-analytic models of galaxy formation to a recipe that relates the occurrence of planets to the mass and metallicity of their host stars, we trace the population of terrestrial planets around both solar-mass (FGK type) and lower-mass (M dwarf) stars throughout all of cosmic history. We find that the mean age of terrestrial planets in the local Universe is $7\pm{}1$ Gyr for FGK hosts and $8\pm{}1$ Gyr for M dwarfs. We estimate that hot Jupiters have depleted the population of terrestrial planets around FGK stars by no more than $\approx 10\%$, and that only $\approx 10\%$ of the terrestrial planets at the current epoch are orbiting stars in a metallicity range for which such planets have yet to be confirmed. The typical terrestrial planet in the local Universe is located in a spheroid-dominated galaxy with a total stellar mass comparable to that of the Milky Way. When looking at the inventory of planets throughout the whole observable Universe, we argue for a total of $\approx 1\times 10^{19}$ and $\approx 5\times 10^{20}$ terrestrial planets around FGK and M stars, respectively. Due to light travel time effects, the terrestrial planets on our past light cone exhibit a mean age of just $1.7\pm 0.2$ Gyr. These results are discussed in the context of cosmic habitability, the Copernican principle and searches for extraterrestrial intelligence at cosmological distances.

Previous work on protoplanetary dust growth shows halt at centimeter sizes owing to the occurrence of bouncing at velocities of $\geq$ 0.1 $ms^{-1}$ and fragmentation at velocities $\geq$ 1 $ms^{-1}$. To overcome these barriers, spatial concentration of cm-sized dust pebbles and subsequent gravitational collapse have been proposed. However, numerical investigations have shown that dust aggregates may undergo fragmentation during the gravitational collapse phase. This fragmentation in turn changes the size distribution of the solids and thus must be taken into account in order to understand the properties of the planetesimals that form. To explore the fate of dust pebbles undergoing fragmenting collisions, we conducted laboratory experiments on dust-aggregate collisions with a focus on establishing a collision model for this stage of planetesimal formation. In our experiments, we analysed collisions of dust aggregates with masses between 1.4 g and 180 g, mass ratios between target and projectile from 125 to 1 at a fixed porosity of 65%, within the velocity range of 1.5 $-$ 8.7 $ms^{-1}$ , at low atmospheric pressure of $\sim 10^{-3}$ mbar and in free-fall conditions. We derived the mass of the largest fragment, the fragment size/mass distribution, and the efficiency of mass transfer as a function of collision velocity and projectile/target aggregate size. Moreover, we give recipes for an easy-to-use fragmentation and mass-transfer model for further use in modeling work. In a companion paper, we utilize the experimental findings and the derived dust-aggregate collision model to investigate the fate of dust pebbles during gravitational collapse.

Some scenarios for planetesimal formation go through a phase of collapse of gravitationally bound clouds of mm-cm-sized pebbles. Such clouds can form for example through the streaming instability in protoplanetary disks. We model the collapse process with a statistical model to obtain the internal structure of planetesimals with solid radii between 10 and 1,000 km. In the collapse, pebbles collide and, depending on relative speed, collisions have different outcomes. A mixture of particle sizes inside a planetesimal leads to better packing capabilities and higher densities. In this paper we apply results from new laboratory experiments of dust aggregate collisions (presented in a companion paper) to model collision outcomes. We find that the internal structure of a planetesimal is strongly dependent on both its mass and the applied fragmentation model. Low-mass planetesimals have no/few fragmenting pebble collisions in the collapse phase and end up as porous pebble-piles. The amount of fragmenting collisions increases with increasing cloud mass, resulting in wider particle size distributions and higher density. The collapse is nevertheless "cold" in the sense that collision speeds are damped by the high collision frequency. This ensures that a significant fraction of large pebbles survive the collapse in all but the most massive clouds. Our results are in broad agreement with the observed increase in density of Kuiper belt objects with increasing size as exemplified by the recent characterization of the highly porous comet 67P/Churyumov-Gerasimenko.

Many observed giant planets lie on eccentric orbits. Such orbits could be the result of strong scatterings with other giant planets. The same dynamical instability that produces these scatterings may also cause habitable planets in interior orbits to become ejected, destroyed, or be transported out of the habitable zone. We say that a habitable planet has resilient habitability if it is able to avoid ejections and collisions and its orbit remains inside the habitable zone. Here we model the orbital evolution of rocky planets in planetary systems where giant planets become dynamically unstable. We measure the resilience of habitable planets as a function of the observed, present-day masses and orbits of the giant planets. We find that the survival rate of habitable planets depends strongly on the giant planet architecture. Equal-mass planetary systems are far more destructive than systems with giant planets of unequal masses. We also establish a link with observation; we find that giant planets with present-day eccentricities higher than 0.4 almost never have a habitable interior planet. For a giant planet with an present-day eccentricity of 0.2 and semimajor axis of 5 AU orbiting a Sun-like star, 50% of the orbits in the habitable zone are resilient to the instability. As semimajor axis increases and eccentricity decreases, a higher fraction of habitable planets survive and remain habitable. However, if the habitable planet has rocky siblings, there is a significant risk of rocky planet collisions that would sterilize the planet.

Solving the mystery of the origin of chondrules is one of the most elusive goals in the field of meteoritics. Recently the idea of planet(esimal) collisions releasing splashes of lava droplets, long considered out of favor, has been reconsidered as a possible origin of chondrules by several papers. One of the main problems with this idea is the lack of quantitative and simple models that can be used to test this scenario by directly comparing to the many known observables of chondrules. In Paper I of this series we presented a simple thermal evolution model of a spherically symmetric expanding cloud of molten lava droplets that is assumed to emerge from a collision between two planetesimals. In the present paper, number II of this series, we use this model to calculate whether or not volatile elements such as Na and K will remain abundant in these droplets or whether they will get depleted due to evaporation. The high density of the droplet cloud (e.g. small distance between adjacent droplets) causes the vapor to quickly reach saturation pressure and thus shutting down further evaporation. We show to which extent, and under which conditions, this keeps the abundances of these elements high, as is seen in chondrules. We find that for most parameters of our model (cloud mass, expansion velocity, initial temperature) the volatile elements Mg, Si and Fe remain entirely in the chondrules. The Na and K abundances inside the droplets will initially stay mostly at their initial values due to the saturation of the vapor pressure, but at some point start to drop due to the cloud expansion. However, as soon as the temperature starts to decrease, most or all of the vapor recondenses again. At the end the Na and K elements retain most of their initial abundances, albeit occasionally somewhat reduced, depending on the parameters of the expanding cloud model.

HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~ 30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by preventing close encounters. We simulate planetary systems similar to the HR 8799 planetary system, exploring the parameter space in separation between the orbits, planetary masses and distance from the Sun to the star. We find systems that look like HR 8799 and remain stable for longer than the estimated age of HR 8799. None of our systems are forced into resonances. We find, with nominal masses and in a narrow range of orbit separations, that 5 of 100 systems match the observations and lifetime. Considering a broad range of orbit separations, we find 12 of 900 similar systems. The systems survive significantly longer because of their slightly increased initial orbit separations compared to assuming circular orbits from the observed positions. A small increase in separation leads to a significant increase in survival time. The low eccentricity the orbits develop from gravitational interaction is enough for the planets to match the observations. With lower masses, but still comfortably within the estimated planet mass uncertainty, we find 18 of 100 matching and long-lived systems in a narrow orbital separation range. In the broad separation range, we find 82 of 900 matching systems. Our results imply that the planets in the HR 8799 system do not have to be in strong mean motion resonances.

The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff, particularly because of small particles and/or strong local solid concentration, that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can for example be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

The size distribution of asteroids in the solar system suggests that they formed top-down, with 100-1000 km bodies forming from the gravitational collapse of dense clumps of small solid particles. We investigate the conditions under which solid particles can form dense clumps in a protoplanetary disc. We used a hydrodynamic code to model the solid-gas interaction in disc. We found that particles down to millimeter size can form dense clumps, but only in regions where solids make $\sim$ 8% of the local surface density. More generally, we mapped the range of particle sizes and concentrations that is consistent with the formation of particle clumps.

The behaviour of sedimenting particles depends on the dust-to-gas ratio of the fluid. Linear stability analysis shows that solids settling in the Epstein drag regime would remain homogeneously distributed in non-rotating incompressible fluids, even when dust-to-gas ratios reach unity. However, the non-linear evolution has not been probed before. Here, we present numerical calculations indicating that in a particle-dense mixture solids spontaneously mix out of the fluid and form swarms overdense in particles by at least a factor 10. The instability is caused by mass-loaded regions locally breaking the equilibrium background stratification. The driving mechanism depends on non-linear perturbations of the background flow and shares some similarity to the streaming instability in accretion discs. The resulting particle-rich swarms may stimulate particle growth by coagulation. In the context of protoplanetary discs, the instability could be relevant for aiding small particles to settle to the midplane in the outer disc. Inside the gas envelopes of protoplanets, enhanced settling may lead to a reduced dust opacity, which facilitates the contraction of the envelope. We show that the relevant physical set up can be recreated in a laboratory setting. This will allow our numerical calculations to be investigated experimentally in the future.

The formation of planetary cores must proceed rapidly in order for the giant planets to accrete their gaseous envelopes before the dissipation of the protoplanetary gas disc (<3 Myr). In orbits beyond 10 AU, direct accumulation of planetesimals by the cores is too slow. Fragments of planetesimals could be accreted faster, but planetesimals are likely too large for fragmentation to be efficient, and resonant trapping poses a further hurdle. Here we instead investigate the accretion of small pebbles (mm-cm sizes) that are the natural outcome of an equilibrium between the growth and radial drift of particles. We construct a simplified analytical model of dust coagulation and pebble drift in the outer disc, between 5 AU and 100 AU, which gives the temporal evolution of the solid surface density and the dominant particle size. These two key quantities determine how core growth proceeds at various orbital distances. We find that pebble surface densities are sufficiently high to achieve the inside-out formation of planetary cores within the disc lifetime. The overall efficiency by which dust gets converted to planets can be high, close to 50 % for planetary architectures similar to the Solar System. Growth by pebble accretion in the outer disc is sufficiently fast to overcome catastrophic Type I migration of the cores. These results require protoplanetary discs with large radial extent (~100 AU) and assume a low number of initial seed embryos. Our findings imply that protoplanetary discs with low disc masses, as expected around low-mass stars (<1 M_sun), or with sub-solar dust-to-gas ratios, do not easily form gas-giant planets (M > 100 M_E), but preferentially form Neptune-mass planets or smaller (< 10 M_E). This is consistent with exoplanet surveys which show that gas giants are relatively uncommon around stars of low mass or low metallicity.

We explore the role of dynamics in shaping planetary system multiplicities, focussing on two particular problems. (1) We propose that the lack of close-in super-Earths in hot Jupiter systems is a signature of the migration history of the hot Jupiters and helps to discriminate between different mechanisms of migration. We present N-body simulations of dynamical migration scenarios where proto-hot Jupiters are excited to high eccentricities prior to tidal circularisation and orbital decay. We show that in this scenario, the eccentric giant planet typically destroys planets in the inner system, in agreement with the observed lack of close super-Earth companions to hot Jupiters. (2) We explore the role of the dynamics of outer systems in affecting the multiplicities of close-in systems such as those discovered by Kepler. We consider specifically the effects of planet--planet scattering and Kozai perturbations on an exterior giant planet on the architecture of the inner system, and evaluate the ability of such scenarios to reduce the inner system's multiplicity and contribute to the observed excess of single Kepler planets.

The temperature and density profiles of protoplanetary discs depend crucially on the mass fraction of micrometre-sized dust grains and on their chemical composition. A larger abundance of micrometre-sized grains leads to an overall heating of the disc, so that the water ice line moves further away from the star. An increase in the water fraction inside the disc, maintaining a fixed dust abundance, increases the temperature in the icy regions of the disc and lowers the temperature in the inner regions. Discs with a larger silicate fraction have the opposite effect. Here we explore the consequence of the dust composition and abundance for the formation and migration of planets. We find that discs with low water content can only sustain outwards migration for planets up to 4 Earth masses, while outwards migration in discs with a larger water content persists up to 8 Earth masses in the late stages of the disc evolution. Icy planetary cores that do not reach run-away gas accretion can thus migrate to orbits close to the host star if the water abundance is low. Our results imply that hot and warm super-Earths found in exoplanet surveys could have formed beyond the ice line and thus contain a significant fraction in water. These water-rich super-Earths should orbit primarily around stars with a low oxygen abundance, where a low oxygen abundance is caused by either a low water-to-silicate ratio or by overall low metallicity.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

The formation of planets depends on the underlying protoplanetary disc structure, which influences both the accretion and migration rates of embedded planets. The disc itself evolves on time-scales of several Myr during which both temperature and density profiles change as matter accretes onto the central star. Here we use a detailed model of an evolving disc to determine the growth of planets by pebble accretion and their migration through the disc. Cores that reach their pebble isolation mass accrete gas to finally form giant planets with extensive gas envelopes, while planets that do not reach pebble isolation mass are stranded as ice giants and ice planets containing only minor amounts of gas in their envelopes. Unlike earlier population synthesis models, our model works without any artificial reductions in migration speed and for protoplanetary discs with gas and dust column densities similar to those inferred from observations. We find that in our nominal disc model the emergence of planetary embryos preferably occurs after approximately 2 Myr in order to not exclusively form gas giants, but also ice giants and smaller planets. The high pebble accretion rates ensure that critical core masses for gas accretion can be reached at all orbital distances. Gas giant planets nevertheless experience significant reduction in semi-major axes by migration. Considering instead planetesimal accretion for planetary growth, we show that formation time-scales are too long to compete with the migration time-scales and the dissipation time of the protoplanetary disc. Altogether, we find that pebble accretion overcomes many of the challenges in the formation of ice and gas giants in evolving protoplanetary discs.

Hot Jupiters are giant planets on orbits a few hundredths of an AU. They do not share their system with low-mass close-in planets, despite these latter being exceedingly common. Two migration channels for hot Jupiters have been proposed: through a protoplanetary gas disc or by tidal circularisation of highly-eccentric planets. We show that highly-eccentric giant planets that will become hot Jupiters clear out any low-mass inner planets in the system, explaining the observed lack of such companions to hot Jupiters. A less common outcome of the interaction is that the giant planet is ejected by the inner planets. Furthermore, the interaction can implant giant planets on moderately-high eccentricities at semimajor axes $<1$ AU, a region otherwise hard to populate. Our work supports the hypothesis that most hot Jupiters reached their current orbits following a phase of high eccentricity, possibly excited by other planetary or stellar companions.

Asteroids and meteorites provide key evidence on the formation of planetesimals in the Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation within a population of initially km-scale planetesimals. However, new models challenge this idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the gravitational collapse of small particles which have organised themselves in dense filaments and clusters in the turbulent gas. Particles concentrate passively between eddies down to the smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The streaming instability causes particles to take an active role in the concentration, by piling up in dense filaments whose friction on the gas reduces the radial drift compared to that of isolated particles. In this chapter we review new paradigms for asteroid formation and compare critically against the observed properties of asteroids as well as constraints from meteorites. Chondrules of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. We review recent progress on understanding the incorporation of chondrules into the asteroids, including layered accretion models where chondrules are accreted onto asteroids over millions of years. We highlight in the end ten unsolved questions in asteroid formation where we expect that progress will be made over the next decade.

The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we investigate the conditions under which solid particles can form dense clumps in a protoplanetary disk. We use a hydrodynamic code to model the interaction between solid particles and the gas inside a shearing box inside the disk, considering particle sizes from sub-millimeter-sized chondrules to meter-sized rocks. We find that particles down to millimeter sizes can form dense particle clouds through the run-away convergence of radial drift known as the streaming instability. We make a map of the range of conditions (strength of turbulence, particle mass-loading, disk mass, and distance to the star) which are prone to producing dense particle clumps. Finally, we estimate the distribution of collision speeds between mm-sized particles. We calculate the rate of sticking collisions and obtain a robust upper limit on the particle growth timescale of $\sim$$10^5$ years. This means that mm-sized chondrule aggregates can grow on a timescale much smaller than the disk accretion timescale ($\sim$$10^6 - 10^7$ years). Our results suggest a pathway from the mm-sized grains found in primitive meteorites to fully formed asteroids. We speculate that asteroids may form from a positive feedback loop in which coagualation leads to particle clumping driven by the streaming instability. This clumping, in turn reduces collision speeds and enhances coagulation.} Future simulations should model coagulation and the streaming instability together to explore this feedback loop further.

Chondrules are millimeter-sized spherules that dominate primitive meteorites (chondrites) originating from the asteroid belt. The incorporation of chondrules into asteroidal bodies must be an important step in planet formation, but the mechanism is not understood. We show that the main growth of asteroids can result from gas-drag-assisted accretion of chondrules. The largest planetesimals of a population with a characteristic radius of 100 km undergo run-away accretion of chondrules within ~3 Myr, forming planetary embryos up to Mars sizes along with smaller asteroids whose size distribution matches that of main belt asteroids. The aerodynamical accretion leads to size-sorting of chondrules consistent with chondrites. Accretion of mm-sized chondrules and ice particles drives the growth of planetesimals beyond the ice line as well, but the growth time increases above the disk life time outside of 25 AU. The contribution of direct planetesimal accretion to the growth of both asteroids and Kuiper belt objects is minor. In contrast, planetesimal accretion and chondrule accretion play more equal roles for the formation of Moon-sized embryos in the terrestrial planet formation region. These embryos are isolated from each other and accrete planetesimals only at a low rate. However, the continued accretion of chondrules destabilizes the oligarchic configuration and leads to the formation of Mars-sized embryos and terrestrial planets by a combination of direct chondrule accretion and giant impacts.

The formation of chondrules is one of the oldest unsolved mysteries in meteoritics and planet formation. Recently an old idea has been revived: the idea that chondrules form as a result of collisions between planetesimals in which the ejected molten material forms small droplets which solidify to become chondrules. Pre-melting of the planetesimals by radioactive decay of 26Al would help producing sprays of melt even at relatively low impact velocity. In this paper we study the radiative cooling of a ballistically expanding spherical cloud of chondrule droplets ejected from the impact site. We present results from a numerical radiative transfer models as well as analytic approximate solutions. We find that the temperature after the start of the expansion of the cloud remains constant for a time t_cool and then drops with time t approximately as T ~ T_0[(3/5)t/t_cool+ 2/5]^(-5/3) for t>t_cool. The time at which this temperature drop starts t_cool depends via an analytical formula on the mass of the cloud, the expansion velocity and the size of the chondrule. During the early isothermal expansion phase the density is still so high that we expect the vapor of volatile elements to saturate so that no large volatile losses are expected.

The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.

The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying structure of the protoplanetary disc, specifically the profiles of temperature, gas scale height and density. The commonly used disc structure of the Minimum Mass Solar Nebular (MMSN) is a simple power law in all these quantities. However, protoplanetary disc models with both viscous and stellar heating show several bumps and dips in temperature, scale height and density caused by transitions in opacity, which are missing in the MMSN model. These play an important role in the formation of planets, as they can act as sweet spots for the formation of planetesimals via the streaming instability and affect the direction and magnitude of type-I-migration. We present 2D simulations of accretion discs that feature radiative cooling, viscous and stellar heating, and are linked to the observed evolutionary stages of protoplanetary discs and their host stars. These models allow us to identify preferred planetesimal and planet formation regions in the protoplanetary disc as a function of the disc's metallicity, accretion rate and lifetime. We derive simple fitting formulae that feature all structural characteristics of protoplanetary discs during the evolution of several Myr. These fits are straightforward to apply for modelling any growth stage of planets where detailed knowledge of the underlying disc structure is required.

Accumulation of dust and ice particles into planetesimals is an important step in the planet formation process. Planetesimals are the seeds of both terrestrial planets and the solid cores of gas and ice giants forming by core accretion. Left-over planetesimals in the form of asteroids, trans-Neptunian objects and comets provide a unique record of the physical conditions in the solar nebula. Debris from planetesimal collisions around other stars signposts that the planetesimal formation process, and hence planet formation, is ubiquitous in the Galaxy. The planetesimal formation stage extends from micrometer-sized dust and ice to bodies which can undergo run-away accretion. The latter ranges in size from 1 km to 1000 km, dependent on the planetesimal eccentricity excited by turbulent gas density fluctuations. Particles face many barriers during this growth, arising mainly from inefficient sticking, fragmentation and radial drift. Two promising growth pathways are mass transfer, where small aggregates transfer up to 50% of their mass in high-speed collisions with much larger targets, and fluffy growth, where aggregate cross sections and sticking probabilities are enhanced by a low internal density. A wide range of particle sizes, from mm to 10 m, concentrate in the turbulent gas flow. Overdense filaments fragment gravitationally into bound particle clumps, with most mass entering planetesimals of contracted radii from 100 to 500 km, depending on local disc properties. We propose a hybrid model for planetesimal formation where particle growth starts unaided by self-gravity but later proceeds inside gravitationally collapsing pebble clumps to form planetesimals with a wide range of sizes.

The first stage of planet formation is the accumulation of dust and ice grains into mm-cm-sized pebbles. These pebbles can clump together through the streaming instability and form gravitationally bound pebble 'clouds'. Pebbles inside such a cloud will undergo mutual collisions, dissipating energy into heat. As the cloud loses energy, it gradually contracts towards solid density. We model this process and investigate two important properties of the collapse: (i) the timescale of the collapse and (ii) the temporal evolution of the pebble size distribution. Our numerical model of the pebble cloud is zero-dimensional and treats collisions with a statistical method. We find that planetesimals with radii larger than 100 km collapse on the free-fall timescale of about 25 years. Lower-mass clouds have longer pebble collision timescales and collapse much more slowly, with collapse times of a few hundred years for 10-km-scale planetesimals and a few thousand years for 1-km-scale planetesimals. The mass of the pebble cloud also determines the interior structure of the resulting planetesimal. The pebble collision speeds in low-mass clouds are below the threshold for fragmentation, forming pebble-pile planetesimals consisting of the primordial pebbles from the protoplanetary disk. Planetesimals above 100 km in radius, on the other hand, consist of mixtures of dust (pebble fragments) and pebbles which have undergone substantial collisions with dust and other pebbles. The Rosetta mission to the comet 67P/Churyumov-Gerasimenko and the New Horizons mission to Pluto will both provide valuable information about the structure of planetesimals in the Solar System. Our model predicts that 67P is a pebble-pile planetesimal consisting of primordial pebbles from the Solar Nebula, while the pebbles in the cloud which contracted to form Pluto must have been ground down substantially during the collapse.

In the Solar System giant planets come in two flavours: 'gas giants' (Jupiter and Saturn) with massive gas envelopes and 'ice giants' (Uranus and Neptune) with much thinner envelopes around their cores. It is poorly understood how these two classes of planets formed. High solid accretion rates, necessary to form the cores of giant planets within the life-time of protoplanetary discs, heat the envelope and prevent rapid gas contraction onto the core, unless accretion is halted. We find that, in fact, accretion of pebbles (~ cm-sized particles) is self-limiting: when a core becomes massive enough it carves a gap in the pebble disc. This halt in pebble accretion subsequently triggers the rapid collapse of the super-critical gas envelope. As opposed to gas giants, ice giants do not reach this threshold mass and can only bind low-mass envelopes that are highly enriched by water vapour from sublimated icy pebbles. This offers an explanation for the compositional difference between gas giants and ice giants in the Solar System. Furthermore, as opposed to planetesimal-driven accretion scenarios, our model allows core formation and envelope attraction within disc life-times, provided that solids in protoplanetary discs are predominantly in pebbles. Our results imply that the outer regions of planetary systems, where the mass required to halt pebble accretion is large, are dominated by ice giants and that gas-giant exoplanets in wide orbits are enriched by more than 50 Earth masses of solids.

The streaming instability is a promising mechanism to overcome the barriers in direct dust growth and lead to the formation of planetesimals. Most previous studies of the streaming instability, however, were focused on a local region of a protoplanetary disk with a limited simulation domain such that only one filamentary concentration of solids has been observed. The characteristic separation between filaments is therefore not known. To address this, we conduct the largest-scale simulations of the streaming instability to date, with computational domains up to 1.6 gas scale heights both horizontally and vertically. The large dynamical range allows the effect of vertical gas stratification to become prominent. We observe more frequent merging and splitting of filaments in simulation boxes of high vertical extent. We find multiple filamentary concentrations of solids with an average separation of about 0.2 local gas scale heights, much higher than the most unstable wavelength from linear stability analysis. This measures the characteristic separation of planetesimal forming events driven by the streaming instability and thus the initial feeding zone of planetesimals.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

Understanding the dominant brown dwarf and giant planet formation processes, and finding out whether these processes rely on completely different mechanisms or share common channels represents one of the major challenges of astronomy and remains the subject of heated debates. It is the aim of this review to summarize the latest developments in this field and to address the issue of origin by confronting different brown dwarf and giant planet formation scenarios to presently available observational constraints. As examined in the review, if objects are classified as "Brown Dwarfs" or "Giant Planets" on the basis of their formation mechanism, it has now become clear that their mass domains overlap and that there is no mass limit between these two distinct populations. Furthermore, while there is increasing observational evidence for the existence of non-deuterium burning brown dwarfs, some giant planets, characterized by a significantly metal enriched composition, might be massive enough to ignite deuterium burning in their core. Deuterium burning (or lack of) thus plays no role in either brown dwarf or giant planet formation. Consequently, we argue that the IAU definition to distinguish these two populations has no physical justification and brings scientific confusion. In contrast, brown dwarfs and giant planets might bear some imprints of their formation mechanism, notably in their mean density and in the physical properties of their atmosphere. Future direct imaging surveys will undoubtedly provide crucial information and perhaps provide some clear observational diagnostics to unambiguously distinguish these different astrophysical objects.

We show that condensation is an efficient particle growth mechanism, leading to growth beyond decimeter-sized pebbles close to an ice line in protoplanetary discs. As coagulation of dust particles is frustrated by bouncing and fragmentation, condensation could be a complementary, or even dominant, growth mode in the early stages of planet formation. Ice particles diffuse across the ice line and sublimate, and vapour diffusing back across the ice line recondenses onto already existing particles, causing them to grow. We develop a numerical model of the dynamical behaviour of ice particles close to the water ice line, approximately 3 AU from the host star. Particles move with the turbulent gas, modelled as a random walk. They also sediment towards the midplane and drift radially towards the central star. Condensation and sublimation are calculated using a Monte Carlo approach. Our results indicate that, with a turbulent alpha-value of 0.01, growth from millimeter to at least decimeter-sized pebbles is possible on a time scale of 1000 years. We find that particle growth is dominated by ice and vapour transport across the radial ice line, with growth due to transport across the atmospheric ice line being negligible. Ice particles mix outwards by turbulent diffusion, leading to net growth across the entire cold region. The resulting particles are large enough to be sensitive to concentration by streaming instabilities, and in pressure bumps and vortices, which can cause further growth into planetesimals. In our model, particles are considered to be homogeneous ice particles. Taking into account the more realistic composition of ice condensed onto rocky ice nuclei might affect the growth time scales, by release of refractory ice nuclei after sublimation. We also ignore sticking and fragmentation in particle collisions. These effects will be the subject of future investigations.

Recent numerical simulations have shown long-lived axisymmetric sub- and super-Keplerian flows in protoplanetary disks. These zonal flows are found in local as well as global simulations of disks unstable to the magnetorotational instability. This paper covers our study of the strength and lifetime of zonal flows and the resulting long-lived gas over- and underdensities as functions of the azimuthal and radial size of the local shearing box. We further investigate dust particle concentrations without feedback on the gas and without self-gravity. Strength and lifetime of zonal flows increase with the radial extent of the simulation box, but decrease with the azimuthal box size. Our simulations support earlier results that zonal flows have a natural radial length scale of 5-7 gas pressure scale heights. This is the first study that combines three-dimensional MHD simulations of zonal flows and dust particles feeling the gas pressure. The pressure bumps trap particles with $\textrm{St} = 1$ very efficiently. We show that $\textrm{St} = 0.1$ particles (of some centimeters in size if at $5\textrm{AU}$ in an MMSN) reach a hundred-fold higher density than initially. This opens the path for particles of $\textrm{St} = 0.1$ and dust-to-gas ratio of 0.01 or for particles of $\textrm{St} \geq 0.5$ and dust-to-gas ratio $10^{-4}$ to still reach densities that potentially trigger the streaming instability and thus gravoturbulent formation of planetesimals.

The planet candidates discovered by the Kepler mission provide a rich sample to constrain the architectures and relative inclinations of planetary systems within approximately 0.5 AU of their host stars. We use the triple-transit systems from the Kepler 16-months data as templates for physical triple-planet systems and perform synthetic transit observations. We find that all the Kepler triple-transit and double-transit systems can be produced from the triple-planet templates, given a low mutual inclination of around five degrees. Our analysis shows that the Kepler data contains a population of planets larger than four Earth radii in single-transit systems that can not arise from the triple-planet templates. We explore the hypothesis that high-mass counterparts of the triple-transit systems underwent dynamical instability to produce a population of massive double-planet systems of moderately high mutual inclination. We perform N-body simulations of mass-boosted triple-planet systems and observe how the systems heat up and lose planets, most frequently by planet-planet collisions, yielding transits in agreement with the large planets in the Kepler single-transit systems. The resulting population of massive double-planet systems can nevertheless not explain the additional excess of low-mass planets among the observed single-transit systems and the lack of gas-giant planets in double-transit and triple-transit systems. Planetary instability of systems of triple gas-giant planets can be behind part of the dichotomy between systems hosting one or more small planets and those hosting a single giant planet. The main part of the dichotomy, however, is more likely to have arisen already during planet formation when the formation, migration or scattering of a massive planet, triggered above a threshold metallicity, suppressed the formation of other planets in sub-AU orbits.

In this paper we report on the formation of magnetically-levitating accretion disks around supermassive black holes. The structure of these disks is calculated by numerically modelling tidal disruption of magnetized interstellar gas clouds. We find that the resulting disks are entirely supported by the pressure of the magnetic fields against the component of gravitational force directed perpendicular to the disks. The magnetic field shows ordered large-scale geometry that remains stable for the duration of our numerical experiments extending over 10% of the disk lifetime. Strong magnetic pressure allows high accretion rate and inhibits disk fragmentation. This in combination with the repeated feeding of manetized molecular clouds to a supermassive black hole yields a possible solution to the long-standing puzzle of black hole growth in the centres of galaxies.

The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time-scale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals.

Modelling the formation of super-km-sized planetesimals by gravitational collapse of regions overdense in small particles requires numerical algorithms capable of handling simultaneously hydrodynamics, particle dynamics and particle collisions. While the initial phases of radial contraction are dictated by drag forces and gravity, particle collisions become gradually more significant as filaments contract beyond Roche density. Here we present a new numerical algorithm for treating momentum and energy exchange in collisions between numerical superparticles representing a high number of physical particles. We adopt a Monte Carlo approach where superparticle pairs in a grid cell collide statistically on the physical collision time-scale. Collisions occur by enlarging particles until they touch and solving for the collision outcome, accounting for energy dissipation in inelastic collisions. We demonstrate that superparticle collisions can be consistently implemented at a modest computational cost. In protoplanetary disc turbulence driven by the streaming instability, we argue that the relative Keplerian shear velocity should be subtracted during the collision calculation. If it is not subtracted, density inhomogeneities are too rapidly diffused away, as bloated particles exaggerate collision speeds. Local particle densities reach several thousand times the mid-plane gas density. We find efficient formation of gravitationally bound clumps, with a range of masses corresponding to contracted radii from 100 to 400 km when applied to the asteroid belt and 150 to 730 km when applied to the Kuiper belt, extrapolated using a constant self-gravity parameter. The smaller planetesimals are not observed at low resolution, but the masses of the largest planetesimals are relatively independent of resolution and treatment of collisions.

We present high-resolution computer simulations of dust dynamics and planetesimal formation in turbulence generated by the magnetorotational instability. We show that the turbulent viscosity associated with magnetorotational turbulence in a non-stratified shearing box increases when going from 256^3 to 512^3 grid points in the presence of a weak imposed magnetic field, yielding a turbulent viscosity of $\alpha\approx0.003$ at high resolution. Particles representing approximately meter-sized boulders concentrate in large-scale high-pressure regions in the simulation box. The appearance of zonal flows and particle concentration in pressure bumps is relatively similar at moderate (256^3) and high (512^3) resolution. In the moderate-resolution simulation we activate particle self-gravity at a time when there is little particle concentration, in contrast with previous simulations where particle self-gravity was activated during a concentration event. We observe that bound clumps form over the next ten orbits, with initial birth masses of a few times the dwarf planet Ceres. At high resolution we activate self-gravity during a particle concentration event, leading to a burst of planetesimal formation, with clump masses ranging from a significant fraction of to several times the mass of Ceres. We present a new domain decomposition algorithm for particle-mesh schemes. Particles are spread evenly among the processors and the local gas velocity field and assigned drag forces are exchanged between a domain-decomposed mesh and discrete blocks of particles. We obtain good load balancing on up to 4096 cores even in simulations where particles sediment to the mid-plane and concentrate in pressure bumps.

We perform hydrodynamical simulations of the accretion of pebbles and rocks onto protoplanets of a few hundred kilometres in radius, including two-way drag force coupling between particles and the protoplanetary disc gas. Particle streams interacting with the gas far out within the Hill sphere of the protoplanet spiral into a prograde circumplanetary disc. Material is accreted onto the protoplanet due to stirring by the turbulent surroundings. We speculate that the trend for prograde rotation among the largest asteroids is primordial and that protoplanets accreted 10%-50% of their mass from pebbles and rocks during the gaseous solar nebula phase. Our model also offers a possible explanation for the narrow range of spin periods observed among the largest bodies in the asteroid and trans-Neptunian belts, and predicts that 1000 km-scale Kuiper belt objects that have not experienced giant impacts should preferentially spin in the prograde direction.

As accretion in protoplanetary disks is enabled by turbulent viscosity, the border between active and inactive (dead) zones constitutes a location where there is an abrupt change in the accretion flow. The gas accumulation that ensues triggers the Rossby wave instability, that in turn saturates into anticyclonic vortices. It was suggested that the trapping of solids within them leads to a burst of planet formation on very short timescales. We perform two-dimensional global simulations of the dynamics of gas and solids in a non-magnetized thin protoplanetary disk with the Pencil Code. We use multiple particle species of radius 1, 10, 30, and 100 cm, solving for the particles' gravitational interaction by a particle-mesh method. The dead zone is modeled as a region of low viscosity. Adiabatic and locally isothermal equations of state are used. We find that the Rossby wave instability is triggered under a variety of conditions, thus making vortex formation a robust process. Inside the vortices, fast accumulation of solids occurs and the particles collapse into objects of planetary mass in timescales as short as five orbits. Because the drag force is size-dependent, aerodynamical sorting ensues within the vortical motion, and the first bound structures formed are composed primarily of similarly-sized particles. In addition to erosion due to ram pressure, we identify gas tides from the massive vortices as a disrupting agent of formed protoplanetary embryos. We also estimate the collisional velocity history of the particles that compose the most massive embryo by the end of the simulation, finding that the vast majority of them never experienced a collision with another particle at speeds faster than 1 m/s.

We present three-dimensional numerical simulations of particle clumping and planetesimal formation in protoplanetary disks with varying amounts of solid material. As centimeter-size pebbles settle to the mid-plane, turbulence develops through vertical shearing and streaming instabilities. We find that when the pebble-to-gas column density ratio is 0.01, corresponding roughly to solar metallicity, clumping is weak, so the pebble density rarely exceeds the gas density. Doubling the column density ratio leads to a dramatic increase in clumping, with characteristic particle densities more than ten times the gas density and maximum densities reaching several thousand times the gas density. This is consistent with unstratified simulations of the streaming instability that show strong clumping in particle dominated flows. The clumps readily contract gravitationally into interacting planetesimals of order 100 km in radius. Our results suggest that the correlation between host star metallicity and exoplanets may reflect the early stages of planet formation. We further speculate that initially low metallicity disks can be particle enriched during the gas dispersal phase, leading to a late burst of planetesimal formation.

The role of magnetic fields for the formation of planets is reviewed. Protoplanetary disc turbulence driven by the magnetorotational instability has a huge influence on the early stages of planet formation. Small dust grains are transported both vertically and radially in the disc by turbulent diffusion, counteracting sedimentation to the mid-plane and transporting crystalline material from the hot inner disc to the outer parts. The conclusion from recent efforts to measure the turbulent diffusion coefficient of magnetorotational turbulence is that turbulent diffusion of small particles is much stronger than naively thought. Larger particles -- pebbles, rocks and boulders -- get trapped in long-lived high pressure regions that arise spontaneously at large scales in the turbulent flow. These gas high pressures, in geostrophic balance with a sub-Keplerian/super-Keplerian zonal flow envelope, are excited by radial fluctuations in the Maxwell stress. The coherence time of the Maxwell stress is only a few orbits, where as the correlation time of the pressure bumps is comparable to the turbulent mixing time-scale, many tens or orbits on scales much greater than one scale height. The particle overdensities contract under the combined gravity of all the particles and condense into gravitationally bound clusters of rocks and boulders. These planetesimals have masses comparable to the dwarf planet Ceres. I conclude with thoughts on future priorities in the field of planet formation in turbulent discs.

We study the behavior of magnetorotational turbulence in shearing box simulations with a radial and azimuthal extent up to ten scale heights. Maxwell and Reynolds stresses are found to increase by more than a factor two when increasing the box size beyond two scale heights in the radial direction. Further increase of the box size has little or no effect on the statistical properties of the turbulence. An inverse cascade excites magnetic field structures at the largest scales of the box. The corresponding 10% variation in the Maxwell stress launches a zonal flow of alternating sub- and super-Keplerian velocity. This in turn generates a banded density structure in geostrophic balance between pressure and Coriolis forces. We present a simplified model for the appearance of zonal flows, in which stochastic forcing by the magnetic tension on short time-scales creates zonal flow structures with life-times of several tens of orbits. We experiment with various improved shearing box algorithms to reduce the numerical diffusivity introduced by the supersonic shear flow. While a standard finite difference advection scheme shows signs of a suppression of turbulent activity near the edges of the box, this problem is eliminated by a new method where the Keplerian shear advection is advanced in time by interpolation in Fourier space.

Centimeter and meter sized solid particles in protoplanetary disks are trapped within long lived high pressure regions, creating opportunities for collapse into planetesimals and planetary embryos. We study the accumulations in the stable Lagrangian points of a giant planet, as well as in the Rossby vortices launched at the edges of the gap it carves. We employ the Pencil Code, tracing the solids with a large number of interacting Lagrangian particles, usually 100,000. For particles of 1 cm to 10 cm radii, gravitational collapse occurs in the Lagrangian points in less than 200 orbits. For 5 cm particles, a 2 Earth mass planet is formed. For 10 cm, the final maximum collapsed mass is around 3 Earth masses. The collapse of the 1 cm particles is indirect, following the timescale of depletion of gas from the tadpole orbits. In the edges of the gap vortices are excited, trapping preferentially particles of 30 cm radii. The rocky planet that is formed is as massive as 17 Earth masses, constituting a Super-Earth. By using multiple particle species, we find that gas drag modifies the streamlines in the tadpole region around the classical L4 and L5 points. As a result, particles of different radii have their stable points shifted to different locations. Collapse therefore takes longer and produces planets of lower mass. Three super-Earths are formed in the vortices, the most massive having 4.4 Earth masses. We conclude that a Jupiter mass planet can induce the formation of other planetary embryos in the outer edge of its gas gap. Trojan Earth mass planets are readily formed, and although not existing in the solar system, might be common in the exoplanetary zoo.

In the borders of the dead zones of protoplanetary disks, the inflow of gas produces a local density maximum that triggers the Rossby wave instability. The vortices that form are efficient in trapping solids. We aim to assess the possibility of gravitational collapse of the solids within the Rossby vortices. We perform global simulations of the dynamics of gas and solids in a low mass non-magnetized self-gravitating thin protoplanetary disk with the Pencil code. We use multiple particle species of radius 1, 10, 30, and 100 cm. The dead zone is modeled as a region of low viscosity. The Rossby vortices excited in the edges of the dead zone are very efficient particle traps. Within 5 orbits after their appearance, the solids achieve critical density and undergo gravitational collapse into Mars sized objects. The velocity dispersions are of the order of 10 m/s for newly formed embryos, later lowering to less than 1 m/s by drag force cooling. After 200 orbits, 38 gravitationally bound embryos were formed inside the vortices, half of them being more massive than Mars. The embryos are composed primarily of same-sized particles. We conclude that the presence of a dead zone naturally gives rise to a population of protoplanetary cores in the mass range of 0.1-0.6 Earth masses, on very short timescales.

Hydromagnetic stresses in accretion discs have been the subject of intense theoretical research over the past one and a half decades. Most of the disc simulations have assumed a small initial magnetic field and studied the turbulence that arises from the magnetorotational instability. However, gaseous discs in galactic nuclei and in some binary systems are likely to have significant initial magnetisation. Motivated by this, we performed ideal magnetohydrodynamic simulations of strongly magnetised, vertically stratified discs in a Keplerian potential. Our initial equilibrium configuration, which has an azimuthal magnetic field in equipartion with thermal pressure, is unstable to the Parker instability. This leads to the expelling of magnetic field arcs, anchored in the midplane of the disc, to around five scale heights from the midplane. Transition to turbulence happens primarily through magnetorotational instability in the resulting vertical fields, although magnetorotational shear instability in the unperturbed azimuthal field plays a significant role as well, especially in the midplane where buoyancy is weak. High magnetic and hydrodynamical stresses arise, yielding an effective $\alpha$-value of around 0.1 in our highest resolution run. Azimuthal magnetic field expelled by magnetic buoyancy from the disc is continuously replenished by the stretching of a radial field created as gas parcels slide in the linear gravity field along inclined magnetic field lines. This dynamo process, where the bending of field lines by the Parker instability leads to re-creation of the azimuthal field, implies that highly magnetised discs are astrophysically viable and that they have high accretion rates.

To treat the problem of growing protoplanetary disc solids across the meter barrier, we consider a very simplified two-component coagulation-fragmentation model that consists of macroscopic boulders and smaller dust grains, the latter being the result of catastrophic collisions between the boulders. Boulders in turn increase their radii by sweeping up the dust fragments. An analytical solution to the dynamical equations predicts that growth by coagulation-fragmentation can be efficient and allow agglomeration of 10-meter-sized objects within the time-scale of the radial drift. These results are supported by computer simulations of the motion of boulders and fragments in 3-D time-dependent magnetorotational turbulence. Allowing however the fragments to diffuse freely out of the sedimentary layer of boulders reduces the density of both boulders and fragments in the mid-plane, and thus also the growth of the boulder radius, drastically. The reason is that the turbulent diffusion time-scale is so much shorter than the collisional time-scale that dust fragments leak out of the mid-plane layer before they can be swept up by the boulders there. Our conclusion that coagulation-fragmentation is not an efficient way to grow across the meter barrier in fully turbulent protoplanetary discs confirms recent results by Brauer, Dullemond, & Henning who solved the coagulation equation in a parameterised turbulence model with collisional fragmentation, cratering, radial drift, and a range of particle sizes. We find that a relatively small population of boulders in a sedimentary mid-plane layer can populate the entire vertical extent of the disc with small grains and that these grains are not first generation dust, but have been through several agglomeration-destruction cycles during the simulations.

We present global 3D MHD simulations of disks of gas and solids, aiming at developing models that can be used to study various scenarios of planet formation and planet-disk interaction in turbulent accretion disks. A second goal is to show that Cartesian codes are comparable to cylindrical and spherical ones in handling the magnetohydrodynamics of the disk simulations, as the disk-in-a-box models presented here develop and sustain MHD turbulence. We investigate the dependence of the magnetorotational instability on disk scale height, finding evidence that the turbulence generated by the magnetorotational instability grows with thermal pressure. The turbulent stresses depend on the thermal pressure obeying a power law of 0.24+/-0.03, compatible with the value of 0.25 found in shearing box calculations. The ratio of stresses decreased with increasing temperature. We also study the dynamics of boulders in the hydromagnetic turbulence. The vertical turbulent diffusion of the embedded boulders is comparable to the turbulent viscosity of the flow. Significant overdensities arise in the solid component as boulders concentrate in high pressure regions.

Proposed mechanisms for the formation of km-sized solid planetesimals face long-standing difficulties. Robust sticking mechanisms that would produce planetesimals by coagulation alone remain elusive. The gravitational collapse of smaller solids into planetesimals is opposed by stirring from turbulent gas. This proceeding describes recent works showing that "particle feedback," the back-reaction of drag forces on the gas in protoplanetary disks, promotes particle clumping as seeds for gravitational collapse. The idealized streaming instability demonstrates the basic ability of feedback to generate particle overdensities. More detailed numerical simulations show that the particle overdensities produced in turbulent flows trigger gravitational collapse to planetesimals. We discuss surprising aspects of this work, including the large (super-Ceres) mass of the collapsing bound cluster, and the finding that MHD turbulence aids gravitational collapse.

This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). 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 (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) 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.

Millimeter interferometry provides evidence for the presence of mm to cm size "pebbles" in the outer parts of disks around pre-main-sequence stars. The observations suggest that large grains are produced relatively early in disk evolution (< 1 Myr) and remain at large radii for longer periods of time (5 to 10 Myr). Simple theoretical estimates of the radial drift time of solid particles, however, imply that they would drift inward over a time scale of less than 0.1 Myr. In this paper, we address this conflict between theory and observation, using more detailed theoretical models, including the effects of sedimentation, collective drag forces and turbulent viscosity. We find that, although these effects slow down the radial drift of the dust particles, this reduction is not sufficient to explain the observationally determined long survival time of mm/cm-sized grains in protoplanetary disks. However, if for some reason the gas to dust ratio in the disk is reduced by at least a factor of 20 from the canonical value of 100 (for instance through photoevaporation of the gas), then the radial drift time scales become sufficiently large to be in agreement with observations.

We present local simulations that verify the linear streaming instability that arises from aerodynamic coupling between solids and gas in protoplanetary disks. This robust instability creates enhancements in the particle density in order to tap the free energy of the relative drift between solids and gas, generated by the radial pressure gradient of the disk. We confirm the analytic growth rates found by Youdin & Goodman (2005) using grid hydrodynamics to simulate the gas and, alternatively, particle and grid representations of the solids. Since the analytic derivation approximates particles as a fluid, this work corroborates the streaming instability when solids are treated as particles. The idealized physical conditions -- axisymmetry, uniform particle size, and the neglect of vertical stratification and collisions -- provide a rigorous, well-defined test of any numerical algorithm for coupled particle-gas dynamics in protoplanetary disks. We describe a numerical particle-mesh implementation of the drag force, which is crucial for resolving the coupled oscillations. Finally we comment on the balance of energy and angular momentum in two-component disks with frictional coupling. A companion paper details the non-linear evolution of the streaming instability into saturated turbulence with dense particle clumps.

We present simulations of the non-linear evolution of streaming instabilities in protoplanetary disks. The two components of the disk, gas treated with grid hydrodynamics and solids treated as superparticles, are mutually coupled by drag forces. We find that the initially laminar equilibrium flow spontaneously develops into turbulence in our unstratified local model. Marginally coupled solids (that couple to the gas on a Keplerian time-scale) trigger an upward cascade to large particle clumps with peak overdensities above 100. The clumps evolve dynamically by losing material downstream to the radial drift flow while receiving recycled material from upstream. Smaller, more tightly coupled solids produce weaker turbulence with more transient overdensities on smaller length scales. The net inward radial drift is decreased for marginally coupled particles, whereas the tightly coupled particles migrate faster in the saturated turbulent state. The turbulent diffusion of solid particles, measured by their random walk, depends strongly on their stopping time and on the solids-to-gas ratio of the background state, but diffusion is generally modest, particularly for tightly coupled solids. Angular momentum transport is too weak and of the wrong sign to influence stellar accretion. Self-gravity and collisions will be needed to determine the relevance of particle overdensities for planetesimal formation.

We study the effect of an imposed vertical magnetic field on the turbulent mass diffusion properties of magnetorotational turbulence in protoplanetary discs. It is well-known that the effective viscosity generated by the turbulence depends strongly on the magnitude of such an external field. In this letter we show that the turbulent diffusion of the flow also grows, but that the diffusion coefficient does not rise with increasing vertical field as fast as the viscosity does. The vertical Schmidt number, i.e. the ratio between viscosity and vertical diffusion, can be close to 20 for high field magnitudes, whereas the radial Schmidt number is increased from below unity to around 3.5. Our results may have consequences for the interpretation of observations of dust in protoplanetary discs and for chemical evolution modelling of these discs.

We perform numerical simulations of the Kelvin-Helmholtz instability in the mid-plane of a protoplanetary disk. A two-dimensional corotating slice in the azimuthal--vertical plane of the disk is considered where we include the Coriolis force and the radial advection of the Keplerian rotation flow. Dust grains, treated as individual particles, move under the influence of friction with the gas, while the gas is treated as a compressible fluid. The friction force from the dust grains on the gas leads to a vertical shear in the gas rotation velocity. As the particles settle around the mid-plane due to gravity, the shear increases, and eventually the flow becomes unstable to the Kelvin-Helmholtz instability. The Kelvin-Helmholtz turbulence saturates when the vertical settling of the dust is balanced by the turbulent diffusion away from the mid-plane. The azimuthally averaged state of the self-sustained Kelvin-Helmholtz turbulence is found to have a constant Richardson number in the region around the mid-plane where the dust-to-gas ratio is significant. Nevertheless the dust density has a strong non-axisymmetric component. We identify a powerful clumping mechanism, caused by the dependence of the rotation velocity of the dust grains on the dust-to-gas ratio, as the source of the non-axisymmetry. Our simulations confirm recent findings that the critical Richardson number for Kelvin-Helmholtz instability is around unity or larger, rather than the classical value of 1/4

We explore the effect of magnetorotational turbulence on the dynamics and concentrations of boulders in local box simulations of a sub-Keplerian protoplanetary disc. The solids are treated as particles each with an independent space coordinate and velocity. We find that the turbulence has two effects on the solids. 1) Meter and decameter bodies are strongly concentrated, locally up to a factor 100 times the average dust density, whereas decimeter bodies only experience a moderate density increase. The concentrations are located in large scale radial gas density enhancements that arise from a combination of turbulence and shear. 2) For meter-sized boulders, the concentrations cause the average radial drift speed to be reduced by 40%. We find that the densest clumps of solids are gravitationally unstable under physically reasonable values for the gas column density and for the dust-to-gas ratio due to sedimentation. We speculate that planetesimals can form in a dust layer that is not in itself dense enough to undergo gravitational fragmentation, and that fragmentation happens in turbulent density fluctuations in this sublayer.

We measure the turbulent diffusion coefficient of dust grains embedded in magnetorotational turbulence in a protoplanetary disc directly from numerical simulations and compare it to the turbulent viscosity of the flow. The simulations are done in a local coordinate frame comoving with the gas in Keplerian rotation. Periodic boundary conditions are used in all directions, and vertical gravity is not applied to the gas. Using a two-fluid approach, small dust grains of various sizes (with friction times up to $\varOmega_0 \tau_{\rm f}=0.02$) are allowed to move under the influence of friction with the turbulent gas. We measure the turbulent diffusion coefficient of the dust grains by applying an external sinusoidal force field acting in the vertical direction on the dust component only. This concentrates the dust around the mid-plane of the disc, and an equilibrium distribution of the dust density is achieved when the vertical settling is counteracted by the turbulent diffusion away from the mid-plane. Comparing with analytical expressions for the equilibrium concentration we deduce the vertical turbulent diffusion coefficient. The vertical diffusion coefficient is found to be lower than the turbulent viscosity and to have an associated vertical diffusion Prandtl number of about 1.5. A similar radial force field also allows us to measure the radial turbulent diffusion coefficient. We find a radial diffusion Prandtl number of about 0.85 and also find that the radial turbulent diffusion coefficient is around 70% higher than the vertical. We also find evidence for trapping of dust grains of intermediate friction time in turbulent eddies.

Local three-dimensional shearing box simulations of the compressible coupled dust-gas equations are used in the fluid approximation to study the evolution of different initial vortex configurations in a protoplanetary disc and their dust-trapping capabilities. The initial conditions for the gas are derived from an analytic solution to the compressible Euler equation and the continuity equation. The solution is valid if there is a vacuum outside the vortex. In the simulations the vortex is either embedded in a hot corona, or it is extended in a cylindrical fashion in the vertical direction. Both configurations are found to survive for at least one orbit and lead to accumulation of dust inside the vortex. This confirms earlier findings that dust accumulates in anticyclonic vortices, indicating that this is a viable mechanism for planetesimal formation.