Source author record

Sijme-Jan Paardekooper

Sijme-Jan Paardekooper appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

17 published item(s)

preprint2022arXiv

Dynamics of Dusty Vortices II: Stability of 2D dust laden vortices

Vortices have long been speculated to play a role in planet formation, via the collection of dust in the pressure maxima that arise at the cores of vortices in protoplanetary discs. The question remains however: as dust collects in the core of a vortex, when does that vortex remain stable and able to collect further dust, and when and why does it break up? We study this question by running high resolution 2D simulations of dust laden vortices. By using the terminal velocity approximation in a local shearing box it was possible to efficiently run simulations of back-reacting dust in a gas at high resolution. Our results show how the stability of 2D dust laden vortices in protoplanetary discs depends on their size relative to the disc scale height, as well as the dust coupling. We find small vortices with semiminor axis much smaller than the scale height to be stable for the duration of the simulations ($t>2000$ orbits). Larger vortices, with semiminor axis smaller than but of order of the scale height, exhibit a drag instability after undergoing a long period of contraction where the core becomes progressively more dust rich. The lifetime of these vortices depends on the dust size, with larger dust grains causing the instability to occur sooner. For the size ranges tested in this paper, $μ$m to mm sized grains, vortices survived for several hundreds of orbits. The result implies that the stability of vortices formed by vertical shear instability and zombie vortex instability, or the breakup of larger vortices through hydrodynamic instabilities, is affected by the presence of dust in the disc. The lifetimes observed in this paper, while shortened by the presence of dust for larger vortices, were still long enough to lead to considerable dust enrichment in the vortex cores.

preprint2021arXiv

Polydisperse Streaming Instability II. Methods for solving the linear stability problem

Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The Polydispserse Streaming Instability is a generalisation of the Streaming Instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretisation of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying Python package psitools.

preprint2021arXiv

Polydisperse Streaming Instability III. Dust evolution encourages fast instability

Planet formation via core accretion requires the production of km-sized planetesimals from cosmic dust. This process must overcome barriers to simple collisional growth, for which the Streaming Instability (SI) is often invoked. Dust evolution is still required to create particles large enough to undergo vigorous instability. The SI has been studied primarily with single size dust, and the role of the full evolved dust distribution is largely unexplored. We survey the Polydispserse Streaming Instability (PSI) with physical parameters corresponding to plausible conditions in protoplanetary discs. We consider a full range of particle stopping times, generalized dust size distributions, and the effect of turbulence. We find that, while the PSI grows in many cases more slowly with a interstellar power-law dust distribution than with a single size, reasonable collisional dust evolution, producing an enhancement of the largest dust sizes, produces instability behaviour similar to the monodisperse case. Considering turbulent diffusion the trend is similar. We conclude that if fast linear growth of PSI is required for planet formation, then dust evolution producing a distribution with peak stopping times on the order of 0.1 orbits and an enhancement of the largest dust significantly above the single power-law distribution produced by a fragmentation cascade is sufficient, along with local enhancement of the dust to gas volume mass density ratio to order unity.

preprint2020arXiv

Low-mass planet migration in three dimensional wind-driven inviscid discs: a negative corotation torque

We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.

preprint2019arXiv

Constraining the masses of planets in protoplanetary discs from the presence or absence of vortices -- Comparison with ALMA observations

A massive planet in a protoplanetary disc will open a gap in the disc material. A steep gap edge can be hydrodynamically unstable, which results in the formation of vortices that can act as tracers for the presence of planets in observational results. However, in a viscous disc, the potential formation of these vortices is dependent on the timescale over which the massive planet accretes mass and with a sufficiently long timescale it is possible for no vortices to form. Hence, there is a connection between the presence of vortices and the growth timescale of the planet and it may therefore be possible to exclude a planetary interpretation of observed structure from the absence of vortices. We have investigated the effect of the planet growth timescale on vortex formation for a range of planet masses and viscosities and have found an approximate relation between the planet mass, viscosity and planet growth timescale for which vortices are not formed within the disc. We then interpret these results in the light of recent observations. We have also found that planets do not need to be close to a Jupiter mass to form vortices in the disc if these discs have low viscosity, as these can be caused by planets as small as a few Neptune masses.

preprint2016arXiv

A terrestrial planet candidate in a temperate orbit around Proxima Centauri

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($α$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

preprint2015arXiv

Spiral arms in scattered light images of protoplanetary discs: Are they the signposts of planets?

One of the striking discoveries of protoplanetary disc research in recent years are the spiral arms seen in several transitional discs in polarised scattered light. An interesting interpretation of the observed spiral features is that they are density waves launched by one or more embedded (proto-)planets in the disc. In this paper we investigate whether planets can be held responsible for the excitation mechanism of the observed spirals. We use locally isothermal hydrodynamic simulations as well as analytic formulae to model the spiral waves launched by planets. Then H-band scattered light images are calculated using a 3D continuum radiative transfer code to study the effect of surface density and pressure scale height perturbation on the detectability of the spirals. We find that a relative change of about 3.5 in the surface density is required for the spirals to be detected with current telescopes in the near-infrared for sources at the distance of typical star-forming regions (140pc). This value is a factor of eight higher than what is seen in hydrodynamic simulations. We also find that a relative change of only 0.2 in pressure scale height is sufficient to create detectable signatures under the same conditions. Therefore, we suggest that the spiral arms observed to date in protoplanetary discs are the results of changes in the vertical structure of the disc (e.g. pressure scale height perturbation) instead of surface density perturbations.

preprint2014arXiv

Dynamical corotation torques on low-mass planets

We study torques on migrating low-mass planets in locally isothermal discs. Previous work on low-mass planets generally kept the planet on a fixed orbit, after which the torque on the planet was measured. In addition to these static torques, when the planet is allowed to migrate it experiences dynamical torques, which are proportional to the migration rate and whose sign depends on the background vortensity gradient. We show that in discs a few times more massive than the Minimum Mass Solar Nebula, these dynamical torques can have a profound impact on planet migration. Inward migration can be slowed down significantly, and if static torques lead to outward migration, dynamical torques can take over, taking the planet beyond zero-torque lines set by saturation of the corotation torque in a runaway fashion. This means the region in non-isothermal discs where outward migration is possible can be larger than what would be concluded from static torques alone.

preprint2014arXiv

Forming Circumbinary Planets: N-body Simulations of Kepler-34

Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors towards high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in-situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34 it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.

preprint2013arXiv

Planet or Brown Dwarf? Inferring the Companion Mass in HD 100546 from the Wall Shape using Mid-Infrared Interferometry

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least α >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and α = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least α = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap.

preprint2013arXiv

Planet-Disc Interactions and Early Evolution of Planetary Systems

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

preprint2013arXiv

The formation of systems with closely spaced low-mass planets and the application to Kepler-36

The Kepler-36 system consists of two planets that are spaced unusually close together, near the 7:6 mean motion resonance. While it is known that mean motion resonances can easily form by convergent migration, Kepler-36 is an extreme case due to the close spacing and the relatively high planet masses of 4 and 8 times that of the Earth. In this paper, we investigate whether such a system can be obtained by interactions with the protoplanetary disc. These discs are thought to be turbulent and exhibit density fluctuations which might originate from the magneto-rotational instability. We adopt a realistic description for stochastic forces due to these density fluctuations and perform both long term hydrodynamical and N-body simulations. Our results show that planets in the Kepler-36 mass range can be naturally assembled into a closely spaced planetary system for a wide range of migration parameters in a turbulent disc similar to the minimum mass solar nebula. The final orbits of our formation scenarios tend to be Lagrange stable, even though large parts of the parameter space are chaotic and unstable.

preprint2012arXiv

How not to build Tatooine: the difficulty of in situ formation of circumbinary planets Kepler 16b, Kepler 34b and Kepler 35b

We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepler planets hostile to planetesimal accretion. We then present results from simulations including planetesimal formation and dust accretion, and show that even in the most favourable case of 100% efficient dust accretion, in situ growth starting from planetesimals smaller than ~10 km is difficult for Kepler 16b, Kepler 34b and Kepler 35b. These planets were likely assembled further out in the disk, and migrated inward to their current location.

preprint2012arXiv

Numerical convergence in self-gravitating shearing sheet simulations and the stochastic nature of disc fragmentation

We study numerical convergence in local two-dimensional hydrodynamical simulations of self-gravitating accretion discs with a simple cooling law. It is well-known that there exists a steady gravito-turbulent state, in which cooling is balanced by dissipation of weak shocks, with a net outward transport of angular momentum. Previous results indicated that if cooling is too fast (typical time scale 3/Omega, where Omega is the local angular velocity), this steady state can not be maintained and the disc will fragment into gravitationally bound clumps. We show that, in the two-dimensional local approximation, this result is in fact not converged with respect to numerical resolution and longer time integration. Irrespective of the cooling time scale, gravito-turbulence consists of density waves as well as transient clumps. These clumps will contract because of the imposed cooling, and collapse into bound objects if they can survive for long enough. Since heating by shocks is very local, the destruction of clumps is a stochastic process. High numerical resolution and long integration times are needed to capture this behaviour. We have observed fragmentation for cooling times up to 20/Omega, almost a factor 7 higher than in previous simulations. Fully three-dimensional simulations with a more realistic cooling prescription are necessary to determine the effects of the use of the two-dimensional approximation and a simple cooling law.

preprint2011arXiv

Numerical convergence in self-gravitating disc simulations: initial conditions and edge effects

We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $Ω^{-1}$, where $Ω$ is the local angular velocity.

preprint2011arXiv

Rapid inward migration of planets formed by gravitational instability

The observation of massive exoplanets at large separation from their host star, like in the HR 8799 system, challenges theories of planet formation. A possible formation mechanism involves the fragmentation of massive self-gravitating discs into clumps. While the conditions for fragmentation have been extensively studied, little is known of the subsequent evolution of these giant planet embryos, in particular their expected orbital migration. Assuming a single planet has formed by fragmentation, we investigate its interaction with the gravitoturbulent disc it is embedded in. Two-dimensional hydrodynamical simulations are used with a simple prescription for the disc cooling. A steady gravitoturbulent disc is first set up, after which simulations are restarted including a planet with a range of masses approximately equal to the clump's initial mass expected in fragmenting discs. Planets rapidly migrate inwards, despite the stochastic kicks due to the turbulent density fluctuations. We show that the migration timescale is essentially that of type I migration, with the planets having no time to open a gap. In discs with aspect ratio ~ 0.1 at their forming location, planets with a mass comparable to, or larger than Jupiter's can migrate in as short as 10000 years, that is, about 10 orbits at 100 AU. Massive planets formed at large separation from their star by gravitational instability are thus unlikely to stay in place, and should rapidly migrate towards the inner parts of protoplanetary discs, regardless of the planet mass.