BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Anders Johansen

Anders Johansen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-phastro-ph.IMastro-ph.SRastro-ph.COastro-ph.GAastro-ph.HEcond-matphysics.comp-ph

Catalog footprint

What is connected

49works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

A low accretion efficiency of planetesimals formed at planetary gap edges

Observations and models of giant planets indicate that such objects are enriched in heavy elements compared to solar abundances. The prevailing view is that giant planets accreted multiple Earth masses of heavy elements after the end of core formation. Such late solid enrichment is commonly explained by the accretion of planetesimals. Planetesimals are expected to form at the edges of planetary gaps, and here we address the question of whether these planetesimals can be accreted in large enough amounts to explain the inferred high heavy element contents of giant planets. We perform a series of N-body simulations of the dynamics of planetesimals and planets during the planetary growth phase, taking into account gas drag as well as the enhanced collision cross-section caused by the extended envelopes. We consider the growth of Jupiter and Saturn via gas accretion after reaching the pebble isolation mass and we include their migration in an evolving disk. We find that the accretion efficiency of planetesimals formed at planetary gap edges is very low: less than 10% of the formed planetesimals are accreted even in the most favorable cases, which in our model corresponds to a few Earth-masses. When planetesimals are assumed to form beyond the feeding zone of the planets, extending to a few Hill radii from a planet, accretion becomes negligible. Furthermore, we find that the accretion efficiency increases when the planetary migration distance is increased and that the efficiency does not increase when the planetesimal radii are decreased. Based on these results we conclude that it is difficult to explain the large heavy element content of giant planets with planetesimal accretion during the gas accretion phase. Alternative processes most likely are required, e.g. accretion of vapor deposited by drifting pebbles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Growth after the streaming instability: The radial distance dependence of the planetary growth

Streaming instability is hypothesized to be triggered at particular protoplanetary disk locations where the volume density of the solid particles is enriched comparable to that of the gas. A ring of planetesimals thus forms when this condition is fulfilled locally. These planetesimals collide with each other and accrete inward drifting pebbles from the outer disk to further increase masses. We investigate the growth of the planetesimals that form in a ring-belt at various disk radii. Their initial mass distributions are calculated based on the formula summarized from the streaming instability simulations. We simulate the subsequent dynamical evolution of the planetesimals with a protoplanetary disk model based either on the minimum mass solar nebula (MMSN) or on the Toomre stability criterion. For the MMSN model, both pebble accretion and planetesimal accretion are efficient at a close-in orbit of $0.3$ AU, resulting in the emergence of several super-Earth mass planets after $1$ Myr. For comparison, only the most massive planetesimals undergo substantial mass growth when they are born at $r{=}3$ AU, while the planetesimals at $r{=}30$ AU experience little or no growth. On the other hand, in the denser Toomre disk, the most massive forming planets can reach Earth mass at $t{=}1$ Myr and reach a mass between that of Neptune and that of Saturn within $3$ Myr at $30$ AU and $100$ AU. Both the pebble and planetesimal accretion rate decrease with disk radial distance. Nevertheless, planetesimal accretion is less pronounced than pebble accretion at more distant disk regions. Taken together, the planets acquire higher masses when the disk has a higher gas density, a higher pebble flux, and/or a lower Stokes number of pebbles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mean motion resonance capture in the context of type-I migration

Capture into mean motion resonance (MMR) is an important dynamical mechanism as it shapes the final architecture of a planetary system. We simulate systems of two or three planets undergoing migration with varied initial parameters such as planetary mass and disk surface density and analyse the resulting resonant chains. In contrast to previous studies, our results show that the disk properties have the dominant impact on capture into mean motion resonance, while the total planetary mass barely affects the final system configuration as long as the planet does not open a gap in the disk. We confirm that the adiabatic resonant capture is the correct framework to understand the conditions leading to MMR formation, since its predictions are qualitatively similar to the numerical results. However, we find that the eccentricity damping can facilitate the capture in a given resonance. We find that under typical disk conditions, planets tend to be captured into 2:1 or 3:2 MMRs, which agrees well with the observed exoplanet MMRs. Our results predict two categories of systems: those that have uniform chains of wide resonances (2:1 or 3:2 MMRs) and those that have a more compact inner pair than the outer pair such as 4:3:2 chains. Both categories of resonant chains are present in observed exoplanet systems. On the other hand, chains with a wider inner pair than the outer one are very rare and emerge from stochastic capture. Our work here can be used to link current configuration of exoplanetary systems to the formation conditions within protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Natural separation of two primordial planetary reservoirs in an expanding solar protoplanetary disk

Meteorites display an isotopic composition dichotomy between non-carbonaceous (NC) and carbonaceous (CC) groups, indicating that planetesimal formation in the solar protoplanetary disk occurred in two distinct reservoirs. The prevailing view is that a rapidly formed Jupiter acted as a barrier between these reservoirs. We show a fundamental inconsistency in this model: if Jupiter is an efficient blocker of drifting pebbles, then the interior NC reservoir is depleted by radial drift within a few hundred thousand years. If Jupiter lets material pass it, then the two reservoirs will be mixed. Instead, we demonstrate that the arrival of the CC pebbles in the inner disk is delayed for several million years by the viscous expansion of the protoplanetary disk. Our results support that Jupiter formed in the outer disk (>10 AU) and allowed a considerable amount of CC material to pass it and become accreted by the terrestrial planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nucleation and growth of iron pebbles explains the formation of iron-rich planets akin to Mercury

The pathway to forming the iron-rich planet Mercury remains mysterious. Mercury's core makes up 70% of the planetary mass, which implies a significant enrichment of iron relative to silicates, while its mantle is strongly depleted in oxidized iron. The high core mass fraction is traditionally ascribed to evaporative loss of silicates, e.g. following a giant impact, but the high abundance of moderately volatile elements in the mantle of Mercury is inconsistent with reaching temperatures much above 1,000 K during its formation. Here we explore the nucleation of solid particles from a gas of solar composition that cools down in the hot inner regions of the protoplanetary disc. The high surface tension of iron causes iron particles to nucleate homogeneously (i.e., not on a more refractory substrate) under very high supersaturation. The low nucleation rates lead to depositional growth of large iron pebbles on a sparse population of nucleated iron nano-particles. Silicates in the form of iron-free MgSiO$_3$ nucleate at similar temperatures but obtain smaller sizes due to the much higher number of nucleated particles. This results in a chemical separation of large iron particles from silicate particles with ten times lower Stokes numbers. We propose that such conditions lead to the formation of iron-rich planetesimals by the streaming instability. In this view, Mercury formed by accretion of iron-rich planetesimals with a sub-solar abundance of highly reduced silicate material. Our results imply that the iron-rich planets known to orbit the Sun and other stars are not required to have experienced mantle-stripping impacts. Instead their formation could be a direct consequence of temperature fluctuations in protoplanetary discs and chemical separation of distinct crystal species through the ensuing nucleation process.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A pebble accretion model for the formation of the terrestrial planets in the Solar System

Pebbles of millimeter sizes are abundant in protoplanetary discs around young stars. Chondrules inside primitive meteorites - formed by melting of dust aggregate pebbles or in impacts between planetesimals - have similar sizes. The role of pebble accretion for terrestrial planet formation is nevertheless unclear. Here we present a model where inwards-drifting pebbles feed the growth of terrestrial planets. The masses and orbits of Venus, Earth, Theia (which later collided with the Earth to form the Moon) and Mars are all consistent with pebble accretion onto protoplanets that formed around Mars' orbit and migrated to their final positions while growing. The isotopic compositions of Earth and Mars are matched qualitatively by accretion of two generations of pebbles, carrying distinct isotopic signatures. Finally, we show that the water and carbon budget of Earth can be delivered by pebbles from the early generation before the gas envelope became hot enough to vaporise volatiles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capture of satellites during planetary encounters A case study of the Neptunian moons Triton and Nereid

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Erosion of planetesimals by gas flow

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitoviscous protoplanetary disks with a dust component. III. Evolution of gas, dust, and pebbles

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)

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pebble drift and planetesimal formation in protoplanetary discs with embedded planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pebble-driven Planet Formation around Very Low-mass Stars and Brown Dwarfs

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\%$).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Resonance in the K2-19 system is at odds with its high reported eccentricities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Formation of moon systems around giant planets: Capture and ablation of planetesimals as foundation for a pebble accretion scenario

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the inclinations of the Jupiter Trojans

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Forming the cores of giant planets from the radial pebble flux in protoplanetary discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hot Jupiters and Super-Earths

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How to form asteroids from mm-sized grains

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of the water content in protoplanetary discs on planet migration and formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Initial mass function of planetesimals formed by the streaming instability

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Long-term stability of the HR 8799 planetary system without resonant lock

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spontaneous concentrations of solids through two-way drag forces between gas and sedimenting particles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Survival of habitable planets in unstable planetary systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forming chondrules in impact splashes. I. Radiative cooling model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Growth of asteroids, planetary embryos and Kuiper belt objects by chondrule accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

How to form planetesimals from mm-sized chondrules and chondrule aggregates

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New Paradigms For Asteroid Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The destruction of inner planetary systems during high-eccentricity migration of gas giants

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation of pebble-pile planetesimals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Feeding Zone of Planetesimal Formation by the Streaming Instability

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Separating gas-giant and ice-giant planets by halting pebble accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The multifaceted planetesimal formation process

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The structure of protoplanetary discs around evolving young stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravoturbulent Planetesimal Formation: The Positive Effect of long-lived Zonal Flows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ice condensation as a planet formation mechanism

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Can planetary instability explain the Kepler dichotomy?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetically-levitating disks around supermassive black holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rapid growth of gas-giant cores by pebble accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-resolution simulations of planetesimal formation in turbulent protoplanetary discs

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 $α\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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Prograde rotation of protoplanets by accretion of pebbles in a gaseous environment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Particle Clumping and Planetesimal Formation Depend Strongly on Metallicity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The role of magnetic fields for planetary 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Zonal Flows and Long-Lived Axisymmetric Pressure Bumps in Magnetorotational Turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Protoplanetary Disk Turbulence Driven by the Streaming Instability: Linear Evolution and Numerical Methods

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Rapid planetesimal formation in turbulent circumstellar 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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Large financial crashes

We propose that large stock market crashes are analogous to critical points studied in statistical physics with log-periodic correction to scaling. We extend our previous renormalization group model of stock market prices prior to and after crashes [D. Sornette et al., J.Phys.I France 6, 167, 1996] by including the first non-linear correction. This predicts the existence of a log-frequency shift over time in the log-periodic oscillations prior to a crash. This is tested on the two largest historical crashes of the century, the october 1929 and october 1987 crashes, by fitting the stock market index over an interval of 8 years prior to the crashes. The good quality of the fits, as well as the consistency of the parameter values obtained from the two crashes, promote the theory that crashes have their origin in the collective ``crowd'' behavior of many interacting agents.

0 citations0 reviewsNo code linkedOpen access