Showing 1–50 of 141 results for author: Kley, W

Viscous circumbinary protoplanetary discs -- I. Structure of the inner cavity

Authors: Anna B. T. Penzlin, Richard A. Booth, Richard P. Nelson, Christoph M. Schäfer, Wilhelm Kley

Abstract: Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficul… ▽ More Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficult problem, for which there is no complete analytical solution, that predicts the shape of the observed disc accurately. We use numeric simulation to evolve circumbinary discs with varying disc viscosities and investigate the size and shape of the inner cavities in such protoplanetary discs. We have simulated over 140 locally isothermal 2D grid-based disc models for > 3e4 binary orbits each and mapped out the parameter space relevant for protoplanetary discs. With this, it becomes possible to create parametrised profiles for individual discs to compare to observation and find limits to their binary eccentricity or internal viscosity from the simulation data. In the long-term simulations larger cavity sizes than previously considered are possible within the parameter space (< 6 binary separations). As an example, we find that the eccentricity of the disc around HD 142527 suggests the impact of the binary dynamics on the disc. However, even considering the larger cavity sizes, the large size of the cavity in HD 142527 remains unexplained by the simulations considering the most recent orbital constraints. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 14 figures (appendix: 2 pages, 3 figures); accepted by mnras

Overstability of the 2:1 mean motion resonance: Exploring disc parameters with hydrodynamic simulations

Authors: Zahra Afkanpour, Sareh Ataiee, Alexandros Ziampras, Anna B. T. Penzlin, Rafael Sfair, Christoph Schäfer, Wilhelm Kley, Hilke Schlichting

Abstract: Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overst… ▽ More Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overstability in a planetary system with two moderate-mass planets in a first-order 2:1 mean motion resonance. We also aim to analyse the impact of key disc parameters, namely the viscosity, surface density, and aspect ratio, on the occurrence of overstability in this planetary system when the masses of the planets are kept constant. We performed 2D locally isothermal hydrodynamical simulations of two planets, with masses of 5 and 10 $M_{\oplus}$, in a 2:1 resonance. Upon identifying the fiducial model in which the system exhibits overstability, we performed simulations with different disc parameters to explore the effects of the disc on the overstability of the system. We observe an overstable planetary system in our hydrodynamic simulations. In the parameter study, we note that overstability occurs in discs characterised by low surface density and low viscosity. Increasing the surface density reduces the probability of overstability within the system. A limit cycle was observed in a specific viscous model with $α_ν = 10^{-3}$. In almost all our models, planets create partial gaps in the disc, which affects both the migration timescale and structure of the planetary system. We demonstrate the existence of overstability using hydrodynamic simulations but find deviations from the analytic approximation and show that the main contribution to this deviation can be attributed to dynamic gap opening. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 16 pages (including appendices), 27 figures, 3 tables, accepted for publication by A&A

Journal ref: A&A 686, A277 (2024)

FargoCPT: A 2D Multi-Physics Code for Simulating the Interaction of Disks with Stars, Planets and Particles

Authors: Thomas Rometsch, Lucas M. Jordan, Tobias W. Moldenhauer, Dennis Wehner, Steven Rendon Restrepo, Tobias W. A. Müller, Giovanni Picogna, Wilhelm Kley, Cornelis P. Dullemond

Abstract: Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, a… ▽ More Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, an updated version of FARGO incorporating other previous enhancements to the code, to provide a simulation environment tailored to study interactions between stars, planets, and disks, ensuring accurate representation of planet systems, hydrodynamics, and dust dynamics with a focus on usability. Methods: The radiation-hydrodynamics part of FargoCPT uses a second-order upwind scheme in 2D polar coordinates supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. Integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapting to drag laws for the Stokes and Epstein regimes. Moreover, FargoCPT provides mechanisms to simulate accretion onto the stars and planets. Results: The code has been tested in practice by its use in various publications. Additionally, it comes with an automated test suite to test the physics modules. Conclusions: FargoCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in the ongoing exploration of planet formation. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 31 pages, 18 figures; accepted for publication in Astronomy and Astrophysics (A&A)

Stability of coorbital planets around binaries

Authors: Stefan Adelbert, Anna B. T. Penzlin, Christoph M. Schäfer, Wilhelm Kley, Billy Quarles, Rafael Sfair

Abstract: In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the para… ▽ More In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the parameters to study general n-body stability of planet pairs in eccentric horseshoe orbits around binaries. The stability is tested by evolving the planet orbits for $10^5$ binary periods with varying initial semi-major axes and planet eccentricities. The unstable region of a single circumbinary planet is used as a comparison to the investigated coorbital configurations in this work. We confirm previous findings on the stability of single planets and find a first order linear relation between orbit eccentricity and pericentre to identify stable orbits for various binary configurations. Such a linear relation is also found for the stability of 1:1 resonant planets around binaries. Stable orbits for eccentric horseshoe configurations exist with a pericentre closer than seven binary separations and, in the case of Kepler-413, the pericentre of the first stable orbit can be approximated by $r_{c,peri} = (2.88 e_p + 2.46) a_{bin}$. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 10 figures, accepted for publication by A&A

Journal ref: A&A 680, A29 (2023)

Hydrodynamic turbulence in disks with embedded planets

Authors: Alexandros Ziampras, Wilhelm Kley, Richard P. Nelson

Abstract: The vertical shear instability (VSI) is a source of hydrodynamic turbulence that can drive vigorous vertical mixing and moderate levels of accretion in protoplanetary disks, and it could be observable in the near future. With high-resolution three-dimensional numerical hydrodynamics simulations, we modeled the behavior of the VSI in protoplanetary disks with and without embedded planets. We then m… ▽ More The vertical shear instability (VSI) is a source of hydrodynamic turbulence that can drive vigorous vertical mixing and moderate levels of accretion in protoplanetary disks, and it could be observable in the near future. With high-resolution three-dimensional numerical hydrodynamics simulations, we modeled the behavior of the VSI in protoplanetary disks with and without embedded planets. We then measured its accretion and mixing capabilities by comparing the full Reynolds stress, which includes the contribution of nonaxisymmetric features, such as spiral arms and vortices, to the Reynolds stress due to the azimuthally averaged velocity field, which can be attributed to good approximation to the VSI. We verified that the VSI can contribute to the accretion stress and showed that, depending on disk conditions, an embedded planet can coexist with or suppress VSI turbulent stress. Specifically, the presence of spiral shocks launched by a planet or planet-generated vortices can interfere with the VSI near the planet's vicinity, with the instability recovering at large enough distances from the planet or vortex. Our results suggest that observations of VSI signatures are unlikely in disks that contain massive, nonaxisymmetric features. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 12 pages, 14 figures, 1 table; accepted in A&A

Journal ref: A&A 670, A135 (2023)

How cooling influences circumbinary discs

Authors: Prakruti Sudarshan, Anna B. T. Penzlin, Alexandros Ziampras, Wilhelm Kley, Richard P. Nelson

Abstract: Circumbinary disc observations and simulations show large, eccentric inner cavities. Recent work has shown that the shape and size of these cavities depend on the aspect ratio and viscosity of the disc, as well as the binary eccentricity and mass ratio. It has been further shown that, for gaps created by planets, the cooling timescale significantly affects the shape and size of the gap. In this st… ▽ More Circumbinary disc observations and simulations show large, eccentric inner cavities. Recent work has shown that the shape and size of these cavities depend on the aspect ratio and viscosity of the disc, as well as the binary eccentricity and mass ratio. It has been further shown that, for gaps created by planets, the cooling timescale significantly affects the shape and size of the gap. In this study, we consider the effect of different cooling models on the cavity shape in a circumbinary disc. We compare locally isothermal and radiatively cooled disc models to ones with a parametrised cooling timescale ($β$-cooling), implemented in 2D numerical simulations for varying binary eccentricities. While the shape of the cavity for radiative and locally isothermal models remains comparable, the inner disc structure changes slightly, leading to a change in the precession rate of the disc. With $β$-cooled models, the shape and size of the cavity changes dramatically towards values of $β$=1. Based on our findings, we introduce a parametrised $β$ model that accounts for the shorter cooling timescale inside the cavity while adequately reproducing the results of the radiative model, and we highlight that accurate treatment of the thermodynamics inside the cavity has a significant impact in modelling circumbinary systems. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: 12 pages, 11 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 664, A157 (2022)

The morphology of CSCha circumbinary disk suggesting the existence of a Saturn-mass planet

Authors: N. T. Kurtovic, P. Pinilla, Anna B. T. Penzlin, M. Benisty, L. Pérez, C. Ginski, A. Isella, W. Kley, F. Menard, S. Pérez, A. Bayo

Abstract: Planets have been detected in circumbinary orbits in several different systems, despite the additional challenges faced during their formation in such an environment. We investigate the possibility of planetary formation in the spectroscopic binary CS Cha by analyzing its circumbinary disk. The system was studied with high angular resolution ALMA observations at 0.87mm. Visibilities modeling and K… ▽ More Planets have been detected in circumbinary orbits in several different systems, despite the additional challenges faced during their formation in such an environment. We investigate the possibility of planetary formation in the spectroscopic binary CS Cha by analyzing its circumbinary disk. The system was studied with high angular resolution ALMA observations at 0.87mm. Visibilities modeling and Keplerian fitting are used to constrain the physical properties of CS Cha, and the observations were compared to hydrodynamic simulations. Our observations are able to resolve the disk cavity in the dust continuum emission and the 12CO J:3-2 transition. We find the dust continuum disk to be azimuthally axisymmetric (less than 9% of intensity variation along the ring) and of low eccentricity (of 0.039 at the peak brightness of the ring). Under certain conditions, low eccentricities can be achieved in simulated disks without the need of a planet, however, the combination of low eccentricity and axisymmetry is consistent with the presence of a Saturn-like planet orbiting near the edge of the cavity. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 18 pages, 13 figures, accepted in A&A. Watch a 3min summary by the first author in https://youtu.be/6j5ROSj2N-U

Journal ref: A&A 664, A151 (2022)

Optical and Near-infrared View of Planet-forming Disks and Protoplanets

Authors: M. Benisty, C. Dominik, K. Follette, A. Garufi, C. Ginski, J. Hashimoto, M. Keppler, W. Kley, J. Monnier

Abstract: In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spi… ▽ More In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spirals, arcs, and shadows, are seen in all imaged infrared-bright disks to date, and are ubiquitous across spectral types. These substructures are believed to be the result of disk evolution processes, and in particular disk-planet interactions. Since small dust grains that scatter light are tightly bound to the disk's gas, these observations closely trace disk structures predicted by hydrodynamical models and serve as observational tests of the predictions of planet formation theories. We argue that the results of current and next-generation high-contrast imaging surveys will, when combined with complementary data from ALMA, lead to a much deeper understanding of the co-evolution of disks and planets, and the mechanisms by which planets form. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Accepted version, before community feedback. 31 pages (37 with references and appendix), 17 figures. This Chapter is dedicated to our co-author Willy Kley. An obituary written by Richard Nelson can be found at the end of the Chapter

Recycling of the first atmospheres of embedded planets: Dependence on core mass and optical depth

Authors: T. W. Moldenhauer, R. Kuiper, W. Kley, C. W. Ormel

Abstract: Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypo… ▽ More Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypothesis and showed that atmosphere-disk recycling is able to fully compensate for radiative cooling and thereby halt Kelvin-Helmholtz contraction to prevent runaway gas accretion. To understand the parameters that determine the efficiency of atmospheric recycling, we extend our earlier studies by exploring the effects of the core mass, the effect of circumstellar gas on sub-Keplerian orbits (headwind), and the optical depth of the surrounding gas on the recycling timescale. Additionally, we analyze their effects on the size and mass of the forming atmosphere. For the explored parameter space, all simulations eventually reach an equilibrium where heating due to hydrodynamic recycling fully compensates radiative cooling. In this equilibrium, the atmosphere-to-core mass ratio stays well below $10 \, \%$, preventing the atmosphere from becoming self-gravitating and entering runaway gas accretion. Higher core masses cause the atmosphere to become turbulent, which further enhances recycling. Even for our highest core mass of $10 \, M_\mathrm{Earth}$, atmosphere-disk recycling is efficient enough to fully compensate for radiative cooling and prevent the atmosphere from becoming self-gravitating. Hence, in-situ formation of hot Jupiters is very unlikely, and migration of gas giants is a key process required to explain their existence. Our findings imply that atmosphere-disk recycling is the most natural explanation for the prevalence of close-in super-Earths and mini-Neptunes. △ Less

Submitted 24 February, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 18 pages, 13 figures, Accepted for publication by A&A on 2/7/2022

Binary orbital evolution driven by a circumbinary disc

Authors: Anna B. T. Penzlin, Wilhelm Kley, Hugo Audiffren, Christoph M. Schäfer

Abstract: The question whether the interaction of a circumbinary disc with the central binary system leads to shrinking or expansion of the binary orbit has attracted considerable interest as it impacts the evolution of binary black holes and stellar binary stars in their formation phase. We performed two-dimensional hydrodynamical simulations of circumbinary discs for a large parameter set of disc viscosit… ▽ More The question whether the interaction of a circumbinary disc with the central binary system leads to shrinking or expansion of the binary orbit has attracted considerable interest as it impacts the evolution of binary black holes and stellar binary stars in their formation phase. We performed two-dimensional hydrodynamical simulations of circumbinary discs for a large parameter set of disc viscosities and thicknesses and two different binary mass ratios for binaries on circular orbits. For those we measured carefully the net angular momentum and mass transfer between disc and binary system, and evaluate the normalised specific angular momentum accretion, $j_\mathrm{s}$ . This is compared to the theoretical, critical specific angular momentum change $j_\mathrm{s,crit}$ that separates contracting from expanding cases which depends on the the binary's mass ratio and the relative accretion onto the two stars. Using finite and infinite disc models we show that the inferred binary evolution is very similar for both setups and confirm that $j_\mathrm{s}$ can be measured accurately with cylindrical simulations that do not include the central binary. However, to obtain the relative accretion onto the stars for non-equal mass binaries, simulations that cover the whole domain including the binary are required. We find that for thick discs with aspect ratio $h = 0.1$ the binaries expand for all viscosities, while discs with $h = 0.05$ lead to an expansion only for larger viscosities with $α$ exceeding $\sim 0.005$. Overall, the regime of binary expansion extends to a much wider parameter space than previously anticipated, but for thin, low viscosity discs the orbits shrink. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 16 pages, 23 figures, accepted by A&A

Journal ref: A&A 660, A101 (2022)

Survival of planet-induced vortices in 2D disks

Authors: Thomas Rometsch, Alexandros Ziampras, Wilhelm Kley, William Béthune

Abstract: Context: Several observations of protoplanetary disks display non-axisymmetric features, often interpreted as vortices. Numerical modeling has repeatedly shown that gap-opening planets are capable of producing large and long-lasting vortices at their outer gap edge, making massive planets popular candidates as the source of such features. Aims: We explore the lifetime of vortices generated by Jupi… ▽ More Context: Several observations of protoplanetary disks display non-axisymmetric features, often interpreted as vortices. Numerical modeling has repeatedly shown that gap-opening planets are capable of producing large and long-lasting vortices at their outer gap edge, making massive planets popular candidates as the source of such features. Aims: We explore the lifetime of vortices generated by Jupiter-sized planets as a function of the thermal relaxation timescale, the level of turbulence, and the effect of disk self-gravity. Methods: We conduct 2D numerical simulations using the hydrodynamics codes PLUTO and FARGO, scanning through several physical and numerical parameters. Vortex properties are automatically extracted from thousands of simulation snapshots. Results: We find that vortices that spawn at the outer gap edge can survive for about 100-3000 planetary orbits, with the shortest lifetimes occurring for moderately efficient dissipation and cooling. However, we also observe a different regime of long-lasting vortices with lifetimes of at least 15 000 orbits for very low viscosity and very short thermal relaxation timescales. Disk self-gravity significantly shortens the lifetime of regular vortices but still allows long-lived ones to survive. Conclusions: Our results suggest that the cooling timescale plays an important role in vortex formation and lifetime and that planet-generated vortices should be observable at large distances from the star for typical thermal relaxation timescales and low turbulence levels. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 17 pages, 13 figures, 1 table; accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 656, A130 (2021)

Disks in close binary stars: γ-Cephei revisited

Authors: Lucas M. Jordan, Wilhelm Kley, Giovanni Picogna, Francesco Marzari

Abstract: Close binaries ($a_{bin} < 20$ au) are known to harbor planets, yet planet formation is unlikely to succeed in such systems. Studying the dynamics of disks in close binaries can help to understand how those planets could have formed. We study the impact that numerical and physical parameters have on the dynamics of disks in close binaries. We use the $γ$-Cephei system as an example and focus on di… ▽ More Close binaries ($a_{bin} < 20$ au) are known to harbor planets, yet planet formation is unlikely to succeed in such systems. Studying the dynamics of disks in close binaries can help to understand how those planets could have formed. We study the impact that numerical and physical parameters have on the dynamics of disks in close binaries. We use the $γ$-Cephei system as an example and focus on disk quantities such as disk eccentricity and the precession rate as indicators for the dynamical state of the disks. We simulate disks in close binaries by performing two-dimensional radiative hydrodynamical simulations using a modified version of the Fargo code. First, we perform a parameter study for different numerical parameters to confirm that our results are robust. In the second part, we study the effects of different masses and different viscosities on the disks' dynamics. Previous studies on radiative disks in close binaries used too low resolutions and too small simulation domains, which impacted the disk's dynamics. We find that radiative disks in close binaries, after an initialization phase, become eccentric with mean eccentricities between 0.06 to 0.27 and display a slow retrograde precession with periods ranging from $4 - 40T_{bin}$ which depends quadratically on the disk's mean aspect ratio. In general, the disks show a coherent, rigid precession which can be broken, however, by changes in the opacity law reducing the overall eccentricity of the disk. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: Accepted in A&A; 14 Pages, 23 figures

Journal ref: A&A 654, A54 (2021)

Pushing planets into an inner cavity by a resonant chain

Authors: S. Ataiee, W. Kley

Abstract: Context. The orbital distribution of exoplanets indicates an accumulation of super-Earth sized planets close to their host stars in compact systems. When an inward disc-driven migration scenario is assumed for their formation, these planets could have been stopped and might have been parked at an inner edge of the disc, or be pushed through the inner disc cavity by a resonant chain. This topic has… ▽ More Context. The orbital distribution of exoplanets indicates an accumulation of super-Earth sized planets close to their host stars in compact systems. When an inward disc-driven migration scenario is assumed for their formation, these planets could have been stopped and might have been parked at an inner edge of the disc, or be pushed through the inner disc cavity by a resonant chain. This topic has not been properly and extensively studied. Using numerical simulations, we investigate the possibility that the inner planets in a resonant chain can be pushed into the disc inner cavity by outer planets. We performed hydrodynamical and N-body simulations of planetary systems embedded in their nascent disc. The inner edge of the disc was represented in two different ways, resembling either a dead zone inner edge (DZ) or a disc inner boundary (IB). The main difference lies in the steepness of the surface density profile. The innermost planet always has a mass of 10 M Earth , with additional outer planets of equal or higher mass. A steeper profile is able to stop a chain of planets more efficiently than a shallower profile. The final configurations in our DZ models are usually tighter than in their IB counterparts, and therefore more prone to instability. We derive analytical expressions for the stopping conditions based on power equilibrium, and show that the final eccentricities result from torque equilibrium. For planets in thinner discs, we found, for the first time, clear signs for over-stable librations in the hydrodynamical simulations, leading to very compact systems. We also found that the popular N-body simulations may overestimate the number of planets in the disc inner cavity. △ Less

Submitted 13 March, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 22 pages (including appendices), Accepted for publication in A & A

Journal ref: A&A 648, A69 (2021)

Steady State by Recycling prevents Premature Collapse of Protoplanetary Atmospheres

Authors: T. W. Moldenhauer, R. Kuiper, W. Kley, C. W. Ormel

Abstract: In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or… ▽ More In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or cooling) timescale, but it has also been suggested that thermodynamical cooling can be prevented by hydrodynamical planet atmosphere-disk recycling. We investigate the efficacy of the recycling hypothesis in preventing the collapse of the atmosphere, check for the existence of a steady state configuration, and determine the final atmospheric mass to core mass ratio. We use three-dimensional radiation-hydrodynamic simulations to model the formation of planetary proto-atmospheres. Equations are solved in a local frame centered on the planet. Ignoring small oscillations that average to zero over time, the simulations converge to a steady state where the velocity field of the gas becomes constant in time. In a steady state, the energy loss by radiative cooling is fully compensated by the recycling of the low entropy gas in the planetary atmosphere with high entropy gas from the circumstellar disk. For close-in planets, recycling naturally halts the cooling of planetary proto-atmospheres, preventing them from contracting toward the runaway regime and collapsing into gas giants. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Comments: 5 pages, 4 figures, accepted to A&A

Journal ref: A&A 646, L11 (2021)

Spiral structures in gravito-turbulent gaseous disks

Authors: W. Béthune, H. Latter, W. Kley

Abstract: Gravitational instabilities can drive small-scale turbulence and large-scale spiral arms in massive gaseous disks under conditions of slow radiative cooling. These motions affect the observed disk morphology, its mass accretion rate and variability, and could control the process of planet formation via dust grain concentration, processing, and collisional fragmentation. We study gravito-turbulence… ▽ More Gravitational instabilities can drive small-scale turbulence and large-scale spiral arms in massive gaseous disks under conditions of slow radiative cooling. These motions affect the observed disk morphology, its mass accretion rate and variability, and could control the process of planet formation via dust grain concentration, processing, and collisional fragmentation. We study gravito-turbulence and its associated spiral structure in thin gaseous disks subject to a prescribed cooling law. We characterize the morphology, coherence, and propagation of the spirals and examine when the flow deviates from viscous disk models. We used the finite-volume code Pluto to integrate the equations of self-gravitating hydrodynamics in three-dimensional spherical geometry. The gas was cooled over longer-than-orbital timescales to trigger the gravitational instability and sustain turbulence. We ran models for various disk masses and cooling rates. In all cases considered, the turbulent gravitational stress transports angular momentum outward at a rate compatible with viscous disk theory. The dissipation of orbital energy happens via shocks in spiral density wakes, heating the disk back to a marginally stable thermal equilibrium. These wakes drive vertical motions and contribute to mix material from the disk with its corona. They are formed and destroyed intermittently, and they nearly corotate with the gas at every radius. As a consequence, large-scale spiral arms exhibit no long-term global coherence, and energy thermalization is an essentially local process. In the absence of radial substructures or tidal forcing, and provided a local cooling law, gravito-turbulence reduces to a local phenomenon in thin gaseous disks. △ Less

Submitted 30 March, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 20 pages, 21 figures; accepted for publication in A&A

Journal ref: A&A 650, A49 (2021)

Giant planet migration during the disc dispersal phase

Authors: Kristina Monsch, Giovanni Picogna, Barbara Ercolano, Wilhelm Kley

Abstract: Transition discs are expected to be a natural outcome of the interplay between photoevaporation (PE) and giant planet formation. Massive planets reduce the inflow of material from the outer to the inner disc, therefore triggering an earlier onset of disc dispersal due to PE through a process known as Planet-Induced PhotoEvaporation (PIPE). In this case, a cavity is formed as material inside the pl… ▽ More Transition discs are expected to be a natural outcome of the interplay between photoevaporation (PE) and giant planet formation. Massive planets reduce the inflow of material from the outer to the inner disc, therefore triggering an earlier onset of disc dispersal due to PE through a process known as Planet-Induced PhotoEvaporation (PIPE). In this case, a cavity is formed as material inside the planetary orbit is removed by PE, leaving only the outer disc to drive the migration of the giant planet. We investigate the impact of PE on giant planet migration and focus specifically on the case of transition discs with an evacuated cavity inside the planet location. This is important for determining under what circumstances PE is efficient at halting the migration of giant planets, thus affecting the final orbital distribution of a population of planets. For this purpose, we use 2D FARGO simulations to model the migration of giant planets in a range of primordial and transition discs subject to PE. The results are then compared to the standard prescriptions used to calculate the migration tracks of planets in 1D planet population synthesis models. The FARGO simulations show that once the disc inside the planet location is depleted of gas, planet migration ceases. This contradicts the results obtained by the impulse approximation, which predicts the accelerated inward migration of planets in discs that have been cleared inside the planetary orbit. These results suggest that the impulse approximation may not be suitable for planets embedded in transition discs. A better approximation that could be used in 1D models would involve halting planet migration once the material inside the planetary orbit is depleted of gas and the surface density at the 3:2 mean motion resonance location in the outer disc reaches a threshold value of $0.01\,\mathrm{g\,cm^{-2}}$. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 16 pages, 11 figures; accepted for publication in A&A

Journal ref: A&A 646, A169 (2021)

Migration of Jupiter mass planets in low viscosity discs

Authors: E. Lega, R. P. Nelson, A. Morbidelli, W. Kley, W. Béthune, A. Crida, D. Kloster, H. Méheut, T. Rometsch, A. Ziampras

Abstract: Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instabili… ▽ More Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instability at the edges of the gap opened by the planet. Migration is then determined by vortex-planet interactions. Our aim is to study migration in low viscosity 3D discs. We performed numerical simulations using 2D (including self-gravity) and 3D codes. After selecting disc masses for which self-gravity is not important, 3D simulations without self-gravity can be safely used. In our nominal simulation only numerical viscosity is present. We then performed simulations with prescribed viscosity to assess the threshold below which the new migration processes appear. We show that for alpha viscosity <= 1.e-5 two migration modes are possible which differ from classical Type-II migration, in the sense that they are not proportional to the disc's viscosity. The first occurs when the gap opened by the planet is not very deep. This occurs in 3D simulations and/or when a big vortex forms at the outer edge of the planetary gap, diffusing material into the gap. We call this type of migration "vortex-driven migration". This migration is very slow and cannot continue indefinitely, because eventually the vortex dissolves. The second migration mode occurs when the gap is deep so that the planet's eccentricity grows to a value ~0.2 due to inefficient eccentricity damping by corotation resonances. This second, faster migration mode appears to be typical of 2D models in discs with slower damping of temperature's perturbations. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 18 pages, 13 figures

Journal ref: A&A 646, A166 (2021)

Modeling the nonaxisymmetric structure in the HD 163296 disk with planet-disk interaction

Authors: Peter J. Rodenkirch, Thomas Rometsch, Cornelis P. Dullemond, Philipp Weber, Wilhelm Kley

Abstract: High resolution ALMA observations revealed a variety of rich substructures in numerous protoplanetary disks. These structures consist of rings, gaps and asymmetric features. It is debated whether planets can be accounted for these substructures in the dust continuum. Characterizing the origin of asymmetries as seen in HD 163296 might lead to a better understanding of planet formation and the under… ▽ More High resolution ALMA observations revealed a variety of rich substructures in numerous protoplanetary disks. These structures consist of rings, gaps and asymmetric features. It is debated whether planets can be accounted for these substructures in the dust continuum. Characterizing the origin of asymmetries as seen in HD 163296 might lead to a better understanding of planet formation and the underlying physical parameters of the system. We test the possibility of the formation of the crescent-shaped asymmetry in the HD 163296 disk through planet-disk interaction. The goal is to obtain constraints on planet masses and eccentricities and disk viscosities. Two dimensional, multi-fluid, hydrodynamical simulations are performed with the FARGO3D code including three embedded planets. Dust is described with the pressureless fluid approach and is distributed over eight size bins. Resulting grids are post-processed with the radiative transfer code RADMC-3D and the CASA software to model synthetic observations. We find that the crescent-shaped asymmetry can be qualitatively modeled with a Jupiter mass planet at a radial distance of 48 au. Dust is trapped preferably in the trailing Lagrange point L5 with a mass of 10 to 15 earth masses. Increased values of eccentricity of the innermost Jupiter mass planet damages the stability of the crescent-shaped feature and does not reproduce the observed radial proximity to the first prominent ring in the system. Generally, a low level of viscosity ($α\leq 2\cdot10^{-3}$) is necessary to allow the existence of such a feature. Including dust feedback the leading point L4 can dominantly capture dust for dust grains with an initial Stokes number $\leq 3.6\cdot 10^{-2}$. The observational results suggest a negligible effect of dust feedback since only one such feature has been detected so far. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A174 (2021)

Parking planets in circumbinary discs

Authors: A. B. T. Penzlin, W. Kley, R. P. Nelson

Abstract: The Kepler space mission discovered about a dozen planets orbiting around binary stars systems. Most of these circumbinary planets lie near their instability boundaries at about 3 to 5 binary separations. Past attempts to match these final locations through an inward migration process were only successful for the Kepler-16 system. Here, we study 10 circumbinary systems and try to match the final p… ▽ More The Kepler space mission discovered about a dozen planets orbiting around binary stars systems. Most of these circumbinary planets lie near their instability boundaries at about 3 to 5 binary separations. Past attempts to match these final locations through an inward migration process were only successful for the Kepler-16 system. Here, we study 10 circumbinary systems and try to match the final parking locations and orbital parameters of the planets with a disc driven migration scenario. We performed 2D locally isothermal hydrodynamical simulations of circumbinary discs with embedded planets and followed their migration evolution using different values for the disc viscosity and aspect ratio. We found that for the six systems with intermediate binary eccentricities ($0.1 \le e_{bin}\le 0.21$) the final planetary orbits matched the observations closely for a single set of disc parameters, specifically a disc viscosity of $α= 10^{-4}$, and an aspect ratio of $H/r \sim 0.04$. For these systems the planet masses were large enough to open at least a partial gap in their discs as they approach the binary, forcing the discs to become circularized and allowing for further migration towards the binary, leading to good agreement with the observed planetary orbital parameters. For systems with very small or large binary eccentricities the match was not as good because the very eccentric discs and large inner cavities in these cases prevented close-in planet migration. In test simulations with higher than observed planet masses better agreement could be found for those systems. The good agreement for 6 out of the 10 modelled systems, where the relative difference between observed and simulated final planet orbit is $\leq 10\%$, strongly supports the idea that planet migration in the disc brought the planets to their present locations. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Journal ref: A&A 645, A68 (2021)

Migration jumps of planets in transition disks

Authors: Thomas Rometsch, Peter J. Rodenkirch, Wilhelm Kley, Cornelis P. Dullemond

Abstract: Transition disks form a special class of protoplanetary disks that are characterized by a deficiency of disk material close to the star. In a subgroup, inner holes in these disks can stretch out to a few tens of au while there is still mass accretion onto the central star observed. We analyse the proposition that this type of wide transition disks is generated by the interaction of the disk with a… ▽ More Transition disks form a special class of protoplanetary disks that are characterized by a deficiency of disk material close to the star. In a subgroup, inner holes in these disks can stretch out to a few tens of au while there is still mass accretion onto the central star observed. We analyse the proposition that this type of wide transition disks is generated by the interaction of the disk with a system of embedded planets. We performed 2D hydrodynamics simulations of a flat disk using either a locally isothermal equation of state or considering also radiative effects. Two 3 to 9 Jupiter mass planets were embedded in the disk and their dynamical evolution due to disk-planet interaction was followed for over 100 000 years. The simulations account for mass accretion onto the star and planets. We included models with parameters geared to the system PDS 70. To assess the observability of features in our models we performed synthetic ALMA observations. For systems with a more massive inner planet there are phases where both planets migrate outward engaged in a 2:1 mean motion resonance via the Masset-Snellgrove mechanism. In sufficiently massive disks the formation of a vortex in the outer disk can trigger rapid outward migration of the outer planet where its distance increases by tens of au within a few thousand years. Later, the outer planet migrates back inwards settling again into resonance with the inner planet. We call this emerging composite phenomenon a 'migration jump'. Outward migration and the migration jumps are accompanied by a high mass accretion rate onto the star. The synthetic images reveal numerous substructures depending on the type of dynamical behaviour. Our results suggest that the outward migration of two embedded planets is a prime candidate for the explanation of the observed high stellar mass accretion rate in wide transition disks. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: 21 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 643, A87 (2020)

High Resolution Parameter Study of the Vertical Shear Instability

Authors: Natascha Manger, Hubert Klahr, Wilhelm Kley, Mario Flock

Abstract: Theoretical models of protoplanetary disks have shown the Vertical Shear Instability (VSI) to be a prime candidate to explain turbulence in the dead zone of the disk. However, simulations of the VSI have yet to show consistent levels of key disk turbulence parameters like the stress-to-pressure ratio $α$. We aim to reconcile these different values by performing a parameter study on the VSI with fo… ▽ More Theoretical models of protoplanetary disks have shown the Vertical Shear Instability (VSI) to be a prime candidate to explain turbulence in the dead zone of the disk. However, simulations of the VSI have yet to show consistent levels of key disk turbulence parameters like the stress-to-pressure ratio $α$. We aim to reconcile these different values by performing a parameter study on the VSI with focus on the disk density gradient $p$ and aspect ratio $h := H/R$. We use full 2$π$ 3D simulations of the disk for chosen set of both parameters. All simulations are evolved for 1000 reference orbits, at a resolution of 18 cells per h. We find that the saturated stress-to-pressure ratio in our simulations is dependent on the disk aspect ratio with a \review{strong} scaling of $α\propto h^{2.6}$, in contrast to the traditional $α$ model, where viscosity scales as $ν\propto αh^2$ with a constant $α$. We also observe consistent formation of large scale vortices across all investigated parameters. The vortices show uniformly aspect ratios of $χ\approx 10$ and radial widths of approximately 1.5 $H$. With our findings we can reconcile the different values reported for the stress-to-pressure ratio from both isothermal and radiation hydrodynamics models, and show long-term evolution effects of the VSI that could aide in the formation of planetesimals. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 15 pages, 18 figures, submitted to MNRAS

Gap, shadows, spirals, streamers: SPHERE observations of binary-disk interactions in GG Tau A

Authors: M. Keppler, A. Penzlin, M. Benisty, R. van Boekel, T. Henning, R. G. van Holstein, W. Kley, A. Garufi, C. Ginski, W. Brandner, G. H. -M. Bertrang, A. Boccaletti, J. de Boer, M. Bonavita, S. Brown Sevilla, G. Chauvin, C. Dominik, M. Janson, M. Langlois, G. Lodato, A. -L. Maire, F. Ménard, E. Pantin, Ch. Pinte, T. Stolker , et al. (9 additional authors not shown)

Abstract: A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the… ▽ More A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated by periodic perturbations by the binary, tearing off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations we find that the main features, in particular the gap size, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring. △ Less

Submitted 4 June, 2020; v1 submitted 18 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 639, A62 (2020)

Solid accretion onto planetary cores in radiative disks

Authors: Apostolos Zormpas, Giovanni Picogna, Barbara Ercolano, Wilhelm Kley

Abstract: The solid accretion rate, necessary to grow gas giant planetary cores within the disk lifetime, has been a major constraint for theories of planet formation. We tested the solid accretion rate efficiency on planetary cores of different masses embedded in their birth disk, by means of 3D radiation-hydrodynamics, where we followed the evolution of a swarm of embedded solids of different sizes. We fo… ▽ More The solid accretion rate, necessary to grow gas giant planetary cores within the disk lifetime, has been a major constraint for theories of planet formation. We tested the solid accretion rate efficiency on planetary cores of different masses embedded in their birth disk, by means of 3D radiation-hydrodynamics, where we followed the evolution of a swarm of embedded solids of different sizes. We found that using a realistic equation of state and radiative cooling, the disk at 5 au is able to cool efficiently and reduce its aspect ratio. As a result, the pebble isolation mass is reached before the core grows to 10 Earth masses, stopping efficiently the pebble flux and creating a transition disk. Moreover, the reduced isolation mass halts the solid accretion before the core reaches the critical mass, leading to a barrier to giant planet formation, and it explains the large abundance of super-Earth planets in the observed population. △ Less

Submitted 22 June, 2020; v1 submitted 3 April, 2020; originally announced April 2020.

Comments: Published in A&A, 8 pages, 9 figures

Journal ref: A&A 638, A97 (2020)

Spirals inside the millimeter cavity of transition disk SR 21

Authors: G. A. Muro-Arena, C. Ginski, C. Dominik, M. Benisty, P. Pinilla, A. J. Bohn, T. Moldenhauer, W. Kley, D. Harsono, T. Henning, R. G. van Holstein, M. Janson, M. Keppler, F. Ménard, L. M. Pérez, T. Stolker, M. Tazzari, M. Villenave, A. Zurlo, C. Petit, F. Rigal, O. Möller-Nilsson, M. Llored, T. Moulin, P. Rabou

Abstract: Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In thi… ▽ More Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11°. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures

Importance of radiative effects in gap opening by planets in protoplanetary disks

Authors: Alexandros Ziampras, Wilhelm Kley, Cornelis P. Dullemond

Abstract: Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100… ▽ More Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100 au). We test this assumption by conducting hydrodynamics simulations of embedded planets in thin locally isothermal and radiative disks that mimic the systems HD 163296 and AS 209 in order to examine the effect of including the energy equation in a seemingly locally isothermal environment as far as planet-disk interaction is concerned. We find that modeling such disks with an ideal equation of state makes a difference in terms of the number of produced rings and the spiral arm contrast in the disk. Locally isothermal disks produce sharper annular or azimuthal features and overestimate a single planet's gap-opening capabilities by producing multiple gaps. In contrast, planets in radiative disks carve a single gap for typical disk parameters. Consequently, for accurate modeling of planets with semimajor axes up to about 100 au, radiative effects should be taken into account even in seemingly locally isothermal disks. In addition, for the case of AS 209, we find that the primary gap is significantly different between locally isothermal and radiative models. Our results suggest that multiple planets are required to explain the ring-rich structures in such systems. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Comments: 10 pages, 12 figures, 2 tables; accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 637, A50 (2020)

The role of disc torques in forming resonant planetary systems

Authors: S. Ataiee, W. Kley

Abstract: The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet's gravity, ignoring its own gravity. For the non-migrat… ▽ More The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet's gravity, ignoring its own gravity. For the non-migrating the neglect of the disc gravity results in a consistent torque calculation while for the migrating case it is inconsistent. We aim to study how much this inconsistent torque calculation can affect the final configuration of a two-planet system. Our focus will be on low-mass planets because most of the multi-planetary systems, discovered by the Kepler survey, have masses around 10 Earth masses. Performing hydrodynamical simulations of planet-disc interaction, we measure the torques on non-migrating and migrating planets for various disc masses as well as density and temperature slopes with and without considering the disc self-gravity. Using this data, we find a relation that quantifies the inconsistency, use it in an N-body code, and perform an extended parameter study modelling the migration of a planetary system with different planet mass ratios and disc surface densities, in order to investigate the impact of the torque inconsistency on the architecture of the planetary system. Not considering disc self-gravity produces an artificially larger torque on the migrating planet that can result in tighter planetary systems. The deviation of this torque from the correct value is larger in discs with steeper surface density profiles. In hydrodynamical modelling of multi-planetary systems, it is crucial to account for the torque correction, otherwise the results favour more packed systems. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 7 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 635, A204 (2020)

Effect of wind-driven accretion on planetary migration

Authors: C. N. Kimmig, C. P. Dullemond, W. Kley

Abstract: Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a r… ▽ More Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a renaissance, as evidence is mounting that PPDs may be less turbulent than previously thought, and 3-D non-ideal MHD modeling of the wind-launching process is maturing. Aim: We wish to investigate how wind-driven accretion may affect migration. We aim for a qualitative exploration of the main effects, rather than a quantitative prediction. Methods: We perform 2-D hydrodynamic simulations with the FARGO3D code in the $(r,φ)$-plane. The vertical coordinate and the launching of the wind are not treated explicitly. Instead, the torque of the wind onto the disk is treated using a simple 2-parameter formula treating the wind mass loss rate and the lever arm. Results: We find that the wind-driven accretion process has a different way of replenishing the co-orbital region than the viscous accretion. The former always injects mass from the outer edge of the co-orbital region and removes mass from the inner edge, while the latter injects or removes mass from the co-orbital region depending on the radial density gradients in the disk. The migration behavior can differ very much and under certain conditions it can drive rapid type-III-like outward migration. We derive an analytic expression for the parameters under which this outward migration occurs. Conclusion: If wind-driven accretion plays a role in the secular evolution of PPDs, migration studies have to include this process as well, because it can strongly affect the resulting migration rate and direction. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A4 (2020)

The impact of planet wakes on the location and shape of the water iceline in a protoplanetary disk

Authors: Alexandros Ziampras, Sareh Ataiee, Wilhelm Kley, Cornelis P. Dullemond, Clément Baruteau

Abstract: Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric natu… ▽ More Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: Accepted by Astronomy and Astrophysics (A&A). 15 pages, 15 figures, 2 tables

Journal ref: A&A 633, A29 (2020)

1:1 orbital resonance of circumbinary planets

Authors: Anna B. T. Penzlin, Sareh Ataiee, Wilhelm Kley

Abstract: The recent detection of the third planet in Kepler-47 has shown that binary stars can host several planets in circumbinary orbits. To understand the evolution of such systems we have performed two-dimensional hydrodynamic simulations of the circumbinary disc with two embedded planets for several Kepler systems. In two cases, Kepler-47 and -413, the planets are captured in a 1:1 mean-motion resonan… ▽ More The recent detection of the third planet in Kepler-47 has shown that binary stars can host several planets in circumbinary orbits. To understand the evolution of such systems we have performed two-dimensional hydrodynamic simulations of the circumbinary disc with two embedded planets for several Kepler systems. In two cases, Kepler-47 and -413, the planets are captured in a 1:1 mean-motion resonance at the planet parking position near the inner edge of the disc. The orbits are fully aligned, have mean eccentricities of about 0.25 to 0.30, and the planets are entangled in a horseshoe type of motion. Subsequent n-body simulations without the disc show that the configurations are stable. Our results point to the existence of a new class of stable resonant orbits around binary stars. It remains to be seen if such orbits exist in reality. △ Less

Submitted 27 August, 2019; originally announced August 2019.

Comments: 5 pages, accepted by A&A Letters

Journal ref: A&A 630, L1 (2019)

Circumbinary discs with radiative cooling and embedded planets

Authors: Wilhelm Kley, Daniel Thun, Anna B. T. Penzlin

Abstract: As of today ten circumbinary planets orbiting solar type main sequence stars have been discovered. Nearly all of them orbit around the central binary very closely to the region of instability where it is difficult to form them in situ. It is assumed that they formed further out and migrated to their observed position. We extend previous studies to a more realistic thermal disc structure and determ… ▽ More As of today ten circumbinary planets orbiting solar type main sequence stars have been discovered. Nearly all of them orbit around the central binary very closely to the region of instability where it is difficult to form them in situ. It is assumed that they formed further out and migrated to their observed position. We extend previous studies to a more realistic thermal disc structure and determine what parameter influence the final parking location of a planet around a binary star. We perform two-dimensional numerical simulations of viscous accretion discs around a central binary that include viscous heating and radiative cooling from the disc surfaces. We vary the binary eccentricity as well as disc viscosity and mass. Concerning the disc evolution we find that it can take well over 100000 binary orbits until an equilibrium state is reached. As seen previously, we find that the central cavity opened by the binary becomes eccentric and precesses slowly in a prograde sense. Embedded planets migrate to the inner edge of the disc. In cases of lower disc viscosity they migrate further in maintaining a circular orbit, while for high viscosity they are parked further out on an eccentric orbit. The final location of an embedded planet is linked to its ability to open a gap in the disc. Gap opening planets separate inner from outer disc, preventing eccentricity excitation in the latter and making it more circular. This allows embedded planets to migrate closer to the binary, in agreement with the observations. The necessary condition for gap opening and the final planet position depend on the planet mass and disc viscosity. △ Less

Submitted 21 May, 2019; originally announced May 2019.

Comments: 15 pages, 16 figures, accepted by Astronomy and Astrophysics

Journal ref: A&A 627, A91 (2019)

Formation of a planetary Laplace resonance through migration in an eccentric disk - The case of GJ876

Authors: Nicolas P. Cimerman, Wilhelm Kley, Rolf Kuiper

Abstract: Orbital mean motion resonances in planetary systems originate from dissipative processes in disk-planet interactions that lead to orbital migration. In multi-planet systems that host giant planets, the perturbation of the protoplanetary disk strongly affects the migration of companion planets. By studying the well-characterized resonant planetary system around GJ 876 we aim to explore which effect… ▽ More Orbital mean motion resonances in planetary systems originate from dissipative processes in disk-planet interactions that lead to orbital migration. In multi-planet systems that host giant planets, the perturbation of the protoplanetary disk strongly affects the migration of companion planets. By studying the well-characterized resonant planetary system around GJ 876 we aim to explore which effects shape disk-driven migration in such a multi-planet system to form resonant chains. We modelled the orbital migration of three planets embedded in a protoplanetary disk using two-dimensional locally isothermal hydrodynamical simulations. We performed a parameter study by varying the disk thickness, $α$ viscosity, mass as well as the initial position of the planets. Moreover, we have analysed and compared simulations with various boundary conditions at the disk's inner rim. We find that due to the high masses of the giant planets, substantial eccentricity can be excited in the disk. This results in large variations of the torque acting on the outer lower-mass planet, which we attribute to a shift of Lindblad and corotation resonances due to disk eccentricity. Depending on disk parameters, the migration of the outer planet can be stopped in a non-resonant state. In other models, the outer planet is able to open a partial gap and to circularize the disk again, later entering a 2:1 resonance with the most massive planet in the system to complete the observed Laplace resonance. Disk-mediated interactions between planets due to spiral waves and excitation of disk eccentricity cause deviations from smooth inward migration of exterior lower mass planets. Self-consistent modelling of the disk-driven migration of multi-planet systems is thus mandatory. Our results are compatible with a late migration of the outermost planet into the resonant chain, when the giant planet pair already is in resonance. △ Less

Submitted 13 August, 2018; originally announced August 2018.

Comments: 20 pages, 14 figures, accepted in A&A

Journal ref: A&A 618, A169 (2018)

Migration of planets in circumbinary discs

Authors: Daniel Thun, Wilhelm Kley

Abstract: The discovery of planets in close orbits around binary stars raises questions about their formation. It is believed that these planets formed in the outer regions of the disc and then migrated through planet-disc interaction to their current location. Considering five different systems (Kepler-16, -34, -35, -38 and -413) we model planet migration through the disc, with special focus on the final o… ▽ More The discovery of planets in close orbits around binary stars raises questions about their formation. It is believed that these planets formed in the outer regions of the disc and then migrated through planet-disc interaction to their current location. Considering five different systems (Kepler-16, -34, -35, -38 and -413) we model planet migration through the disc, with special focus on the final orbital elements of the planets. We investigate how the final orbital parameters are influenced by the disc and planet masses. Using 2D, locally isothermal, viscous hydrodynamical simulations we first model the disc dynamics for all five systems, followed by a study of the migration properties of embedded planets with different masses. For all systems we find that the discs become eccentric and precess slowly. We confirm the bifurcation feature in the precession period -- gap-size diagram for different binary mass ratios. The Kepler-16, -35, -38 and -413 systems lie on the lower branch and Kepler-34 on the upper one. For systems with small binary eccentricity we find a new non-monotonic, loop-like feature. In all systems the planets migrate to the inner edge of the disc cavity. Depending on the planet-disc mass ratio we observe two different regimes. Massive planets can significantly alter the disc structure by compressing and circularising the inner cavity and they remain on nearly circular orbits. Lower mass planets are strongly influenced by the disc, their eccentricity is exited to high values, and their orbits are aligned with the inner disc in a state of apsidal corotation. In our simulations the final locations of the planets are typically too large with respect to the observations because of the large inner gaps of the discs. The migrating planets in the most eccentric discs (around Kepler-34 and -413) show the largest final eccentricity in agreement with the observations. △ Less

Submitted 16 August, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: 19 pages, 19 figures, published in A&A, language edited version

Journal ref: A&A 616, A47 (2018)

Particle accretion onto planets in discs with hydrodynamic turbulence

Authors: Giovanni Picogna, Moritz H. R. Stoll, Wilhelm Kley

Abstract: The growth process of proto-planets can be sped-up by accreting a large number of solid, pebble-sized objects that are still present in the protoplanetary disc. It is still an open question on how efficient this process works in realistic turbulent discs. Here, we investigate the accretion of pebbles in turbulent discs that are driven by the purely hydrodynamical vertical shear instability (VSI).… ▽ More The growth process of proto-planets can be sped-up by accreting a large number of solid, pebble-sized objects that are still present in the protoplanetary disc. It is still an open question on how efficient this process works in realistic turbulent discs. Here, we investigate the accretion of pebbles in turbulent discs that are driven by the purely hydrodynamical vertical shear instability (VSI). For this purpose, we perform global three-dimensional simulations of locally isothermal, VSI turbulent discs with embedded protoplanetary cores from 5 to 100 $M_\oplus$ that are placed at 5.2 au distance from the star. In addition, we follow the evolution of a swarm of embedded pebbles of different size under the action of drag forces between gas and particles in this turbulent flow. Simultaneously, we perform a set of comparison simulations for laminar viscous discs where the particles experience stochastic kicks. For both cases, we measure the accretion rate onto the cores as a function of core mass and Stokes number ($τ_s$) of the particles and compare it to recent MRI turbulence simulations. Overall the dynamic is very similar for the particles in the VSI turbulent disc and the laminar case with stochastic kicks. For the small mass planets (i.e. 5 and 10 $M_\oplus$), well-coupled particles with $τ_s = 1$, which have a size of about one meter at this location, we find an accretion efficiency (rate of particles accreted over drifting inward) of about 1.6-3%. For smaller and larger particles this efficiency is higher. However, the fast inward drift for $τ_s = 1$ particles makes them the most effective for rapid growth, leading to mass doubling times of about 20,000 yr. For masses between 10 and 30 $M_\oplus$ the core reaches the pebble isolation mass and the particles are trapped at the pressure maximum just outside of the planet, shutting off further particle accretion. △ Less

Submitted 16 April, 2018; v1 submitted 23 March, 2018; originally announced March 2018.

Comments: 18 pages, accepted to A&A

Journal ref: A&A 616, A116 (2018)

Growth and evolution of satellites in a Jovian massive disc

Authors: Ricardo A Moraes, Wilhelm Kley, Ernesto Vieira Neto

Abstract: The formation of satellite systems in circum-planetary discs is considered to be similar to the formation of rocky planets in a proto-planetary disc, especially Super-Earths. Thus, it is possible to use systems with large satellites to test formation theories that are also applicable to extrasolar planets. Furthermore, a better understanding of the origin of satellites might yield important inform… ▽ More The formation of satellite systems in circum-planetary discs is considered to be similar to the formation of rocky planets in a proto-planetary disc, especially Super-Earths. Thus, it is possible to use systems with large satellites to test formation theories that are also applicable to extrasolar planets. Furthermore, a better understanding of the origin of satellites might yield important information about the environment near the growing planet during the last stages of planet formation. In this work we investigate the formation and migration of the Jovian satellites through N-body simulations. We simulated a massive, static, low viscosity, circum-planetary disc in agreement with the minimum mass sub-nebula model prescriptions for its total mass. In hydrodynamic simulations we found no signs of gaps, therefore type II migration is not expected. Hence, we used analytic prescriptions for type I migration, eccentricity and inclination damping, and performed N-body simulations with damping forces added. Detailed parameter studies showed that the number of final satellites is strong influenced by the initial distribution of embryos, the disc temperature, and the initial gas density profile. For steeper initial density profiles it is possible to form systems with multiple satellites in resonance while a flatter profile favours the formation of satellites close to the region of the Galilean satellites. We show that the formation of massive satellites such as Ganymede and Callisto can be achieved for hotter discs with an aspect ratio of H/r ~ 0.15 for which the ice line was located around 30 R_J. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 17 pages, 21 figures

Accretion outbursts in massive star formation

Authors: D. M. -A. Meyer, E. I. Vorobyov, R. Kuiper, W. Kley

Abstract: Using the HPC ressources of the state of Baden-Württemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in… ▽ More Using the HPC ressources of the state of Baden-Württemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: Proceedings of the 3rd bwHPC-Symposium, Heidelberg (2016), 2 pages, 2 figures

Journal ref: ISBN 978-3-946531-70-8 - published on 02.10.2017

Forming spectroscopic massive proto-binaries by disk fragmentation

Authors: D. M. -A. Meyer, R. Kuiper, W. Kley, K. G. Johnston, E. Vorobyov

Abstract: The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-… ▽ More The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-hydrodynamic numerical simulations of 100 Mo pre-stellar cores. We find that accretion discs of young massive stars violently fragment without preventing the (highly variable) accretion of gaseous clumps onto the protostars. While acquiring the characteristics of a nascent low-mass companion, some disc fragments migrate onto the central massive protostar with dynamical properties showing that its final Keplerian orbit is close enough to constitute a close massive proto-binary system, having a young high-mass and a low-mass component. We conclude on the viability of the disc fragmentation channel for the formation of such short-period binaries, and that both processes -close massive binary formation and accretion bursts- may happen at the same time. FU-Orionis-type bursts, such as observed in the young high-mass star S255IR-NIRS3, may not only indicate ongoing disc fragmentation, but also be considered as a tracer for the formation of close massive binaries - progenitors of the subsequent massive spectroscopic binaries - once the high-mass component of the system will enter the main-sequence phase of its evolution. Finally, we investigate the ALMA-observability of the disc fragments. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: Accepted for publication by MNRAS, 24 pages, 16 figures

Planet formation and disk-planet interactions

Authors: Wilhelm Kley

Abstract: This review is based on lectures given at the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation' held in March 2015 in Les Diablerets, Switzerland. Starting with an overview of the main characterictics of the Solar System and extrasolar planets, we describe the planet formation process in terms of the sequential accretion scenario. First the growth processes of dust part… ▽ More This review is based on lectures given at the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation' held in March 2015 in Les Diablerets, Switzerland. Starting with an overview of the main characterictics of the Solar System and extrasolar planets, we describe the planet formation process in terms of the sequential accretion scenario. First the growth processes of dust particles to planetesimals and subsequently to terrestrial planets or planetary cores are presented. This is followed by the formation process of the giant planets either by core accretion or gravitational instability. Finally, the dynamical evolution of the orbital elements as driven by disk-planet interaction and the overall evolution of multi-object systems is presented. △ Less

Submitted 22 July, 2017; originally announced July 2017.

Comments: Lecture Notes of the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation', March 2015, 111 Pages

Planet-disc interaction in laminar and turbulent discs

Authors: Moritz Stoll, Giovanni Picogna, Wilhelm Kley

Abstract: In weakly ionized discs turbulence can be generated through the vertical shear instability (VSI). Embedded planets feel a stochastic component in the torques acting on them which can impact their migration. In this work we study the interplay between a growing planet embedded in a protoplanetary disc and the VSI-turbulence. We performed a series of three-dimensional hydrodynamical simulations for… ▽ More In weakly ionized discs turbulence can be generated through the vertical shear instability (VSI). Embedded planets feel a stochastic component in the torques acting on them which can impact their migration. In this work we study the interplay between a growing planet embedded in a protoplanetary disc and the VSI-turbulence. We performed a series of three-dimensional hydrodynamical simulations for locally isothermal discs with embedded planets in the mass range from 5 to 100 Earth masses. We study planets embedded in an inviscid disc that is VSI unstable, becomes turbulent and generates angular momentum transport with an effective $α= 5 \cdot 10^{-4}$. This is compared to the corresponding viscous disc using exactly this $α$-value. In general we find that the planets have only a weak impact on the disc turbulence. Only for the largest planet ($100 M_\oplus$) the turbulent activity becomes enhanced inside of the planet. The depth and width of a gap created by the more massive planets ($30, 100 M_\oplus$) in the turbulent disc equal exactly that of the corresponding viscous case, leading to very similar torque strengths acting on the planet, with small stochastic fluctuations for the VSI disc. At the gap edges vortices are generated that are stronger and longer lived in the VSI disc. Low mass planets (with $M_p \leq 10 M_\oplus$) do not open gaps in the disc in both cases but generate for the turbulent disc an over-density behind the planet that exerts a significant negative torque. This can boost the inward migration in VSI turbulent discs well above the Type I rate. Due to the finite turbulence level in realistic three-dimensional discs the gap depth will always be limited and migration will not stall in inviscid discs. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 10 pages, 7 figures, accepted by A&A

Journal ref: A&A 604, A28 (2017)

The vertical structure of the boundary layer around compact objects

Authors: Marius Hertfelder, Wilhelm Kley

Abstract: Context. Mass transfer due to Roche lobe overflow leads to the formation of an accretion disk around a weakly magnetized white dwarf (WD) in cataclysmic variables. At the inner edge of the disk, the gas comes upon the surface of the WD and has to get rid of its excess kinetic energy in order to settle down on the more slowly rotating outer stellar layers. This region is known as the boundary layer… ▽ More Context. Mass transfer due to Roche lobe overflow leads to the formation of an accretion disk around a weakly magnetized white dwarf (WD) in cataclysmic variables. At the inner edge of the disk, the gas comes upon the surface of the WD and has to get rid of its excess kinetic energy in order to settle down on the more slowly rotating outer stellar layers. This region is known as the boundary layer (BL). Aims. In this work we investigate the vertical structure of the BL, which is still poorly understood. We shall provide details of the basic structure of the two-dimensional (2D) BL and how it depends on parameters such as stellar mass and rotation rate, as well as the mass-accretion rate. We further investigate the destination of the disk material and compare our results with previous one-dimensional (1D) simulations. Methods. We solve the 2D equations of radiation hydrodynamics in a spherical geometry using a parallel grid-based code that employs a Riemann solver. The radiation energy is considered in the two-temperature approach with a radiative flux given by the flux-limited diffusion approximation. Results. The BL around a non-rotating WD is characterized by a steep drop in angular velocity over a width of only 1% of the stellar radius, a heavy depletion of mass, and a high temperature (500 000 K) as a consequence of the strong shear. Variations in Om_star, M_star, and M_dot influence the extent of the changes of the variables in the BL but not the general structure. Depending on Om_star, the disk material travels up to the poles or is halted at a certain latitude. The extent of mixing with the stellar material also depends on Om_star. We find that the 1D approximation matches the 2D data well, apart from an underestimated temperature. △ Less

Submitted 22 May, 2017; originally announced May 2017.

Comments: 15 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 605, A24 (2017)

Numerical Simulations of Regolith Sampling Processes

Authors: Christoph M. Schäfer, Samuel Scherrer, Robert Buchwald, Thomas I. Maindl, Roland Speith, Wilhelm Kley

Abstract: We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation and failure flows for dynamical behaviour of regolith. In the context of projected small body (asteroid or small moons) sample return missions, we investigate… ▽ More We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation and failure flows for dynamical behaviour of regolith. In the context of projected small body (asteroid or small moons) sample return missions, we investigate the efficiency and feasibility of a particular material sampling method: Brushes sweep material from the asteroid's surface into a collecting tray. We analyze the influence of different material parameters of regolith such as cohesion and angle of internal friction on the sampling rate. Furthermore, we study the sampling process in two environments by varying the surface gravity (Earth's and Phobos') and we apply different rotation rates for the brushes. We find good agreement of our sampling simulations on Earth with experiments and provide estimations for the influence of the material properties on the collecting rate. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 15 pages, 11 figures, accepted for publication in Planetary and Space Science

Circumbinary discs: Numerical and physical behaviour

Authors: Daniel Thun, Wilhelm Kley, Giovanni Picogna

Abstract: We study the evolution of circumbinary disks under the gravitational influence of the binary using two-dimensional hydrodynamical simulations to investigate the impact of disk and binary parameters on the dynamical aspects of the disk. To distinguish between physical and numerical effects we apply three hydrodynamical codes. First we analyse in detail numerical issues concerning the conditions at… ▽ More We study the evolution of circumbinary disks under the gravitational influence of the binary using two-dimensional hydrodynamical simulations to investigate the impact of disk and binary parameters on the dynamical aspects of the disk. To distinguish between physical and numerical effects we apply three hydrodynamical codes. First we analyse in detail numerical issues concerning the conditions at the boundaries and grid resolution. We then perform a series of simulations with different binary (eccentricity, mass ratio) and disk parameters (viscosity, aspect ratio) starting from a reference model with Kepler-16 parameters. Concerning the numerical aspects we find that the inner grid radius must be of the order of the binary semi-major axis, with free outflow conditions applied such that mass can flow onto the central binary. A closed inner boundary leads to unstable evolutions. We find that the inner disk turns eccentric and precesses for all investigated physical parameters. The precession rate is slow with periods ($T_\mathrm{prec}$) starting at around 500 binary orbits ($T_\mathrm{bin}$) for high viscosity and large $H/R$ where the inner hole is smaller and more circular. Reducing $α$ and $H/R$ increases the gap size and $T_\mathrm{prec}$ reaches 2500 $T_\mathrm{bin}$. For varying binary mass ratios $q_\mathrm{bin}$ the gap size remains constant whereas $T_\mathrm{prec}$ decreases for increasing $q_\mathrm{bin}$. For varying binary eccentricities $e_\mathrm{bin}$ we find two separate branches in the gap size and eccentricity diagram. The bifurcation occurs at around $e_\mathrm{crit} \approx 0.18$ where the gap is smallest with the shortest $T_\mathrm{prec}$. For $e_\mathrm{bin}$ smaller and larger than $e_\mathrm{crit}$ the gap size and $T_\mathrm{prec}$ increase. Circular binaries create the most eccentric disks. △ Less

Submitted 13 June, 2018; v1 submitted 26 April, 2017; originally announced April 2017.

Comments: 20 pages, 25 figures, published in A&A, language edited version

Journal ref: A&A 604, A102 (2017)

Anisotropic hydrodynamic turbulence in accretion disks

Authors: Moritz H. R. Stoll, Giovanni Picogna, Wilhelm Kley

Abstract: Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether th… ▽ More Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks. △ Less

Submitted 7 February, 2017; v1 submitted 1 February, 2017; originally announced February 2017.

Comments: v2: minor corrections: language and equations

Journal ref: A&A 599, L6 (2017)

The accretion of migrating giant planets

Authors: Christoph Dürmann, Wilhelm Kley

Abstract: Most studies concerning the growth and evolution of massive planets focus either on their accretion or their migration only. In this work we study both processes concurrently to investigate how they might mutually affect each other. We modeled a 2-dimensional disk with a steady accretion flow onto the central star and embed a Jupiter mass planet at 5.2 au. The disk is locally isothermal and viscos… ▽ More Most studies concerning the growth and evolution of massive planets focus either on their accretion or their migration only. In this work we study both processes concurrently to investigate how they might mutually affect each other. We modeled a 2-dimensional disk with a steady accretion flow onto the central star and embed a Jupiter mass planet at 5.2 au. The disk is locally isothermal and viscosity is modeled using a constant $α$. The planet is held on a fixed orbit for a few hundred orbits to allow the disk to adapt and carve a gap. After this period, the planet is released and free to move according to the gravitational interaction with the gas disk. The mass accretion onto the planet is modeled by removing a fraction of gas from the inner Hill sphere, and the removed mass and momentum can be added to the planet. Our results show that a fast migrating planet is able to accrete more gas than a slower migrating planet. Utilizing a tracer fluid we analyzed the origin of the accreted gas which comes predominantly originating from the inner disk for a fast migrating planet. In case of slower migration the fraction of gas from the outer disk increases. We also found that even for very high accretion rates in some cases gas crosses the planetary gap from the inner to the outer disk. Our simulations show that the crossing of gas changes during the migration process as the migration rate slows down. Therefore classical type II migration where the planet migrates with the viscous drift rate and no gas crosses the gap is no general process but may only occur for special parameters and at a certain time during the orbital evolution of the planet. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 9 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 598, A80 (2017)

Bow shock nebulae of hot massive stars in a magnetized medium

Authors: D. M. -A. Meyer, A. Mignone, R. Kuiper, A. Raga, W. Kley

Abstract: A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic field… ▽ More A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic fields, aligned with the motion of the flow, on the formation and emission properties of these circumstellar structures. Our axisymmetric, magneto-hydrodynamical simulations with optically-thin radiative cooling, heating and anisotropic thermal conduction show that the presence of the background ISM magnetic field affects the projected optical emission our bow shocks at Ha and [OIII] lambda 5007 which become fainter by about 1-2 orders of magnitude, respectively. Radiative transfer calculations against dust opacity indicate that the magnetic field slightly diminishes their projected infrared emission and that our bow shocks emit brightly at 60 micron. This may explain why the bow shocks generated by ionizing runaway massive stars are often difficult to identify. Finally, we discuss our results in the context of the bow shock of Zeta Ophiuchi and we support the interpretation of its imperfect morphology as an evidence of the presence of an ISM magnetic field not aligned with the motion of its driving star. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: Accepted by MNRAS, 21 pages, 13 figures

On the existence of accretion-driven bursts in massive star formation

Authors: D. M. -A. Meyer, E. I. Vorobyov, R. Kuiper, W. Kley

Abstract: Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experienc… ▽ More Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experience variable accretion and show that this process is accompanied by luminous outbursts induced by the episodic accretion of gaseous clumps falling from the circumstellar disk onto the protostar. Consequently, the process of accretion-induced luminous flares is also conceivable in the high-mass regime of star formation and we propose to regard this phenomenon as a general mechanism that can affect protostars regardless of their mass and/or the chemical properties of the parent environment in which they form. In addition to the commonness of accretion-driven outbursts in the star formation machinery, we conjecture that luminous flares from regions hosting forming high-mass star may be an observational implication of the fragmentation of their accretion disks. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 5 pages, Accepted for MNRAS Letters

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

Authors: Stefan Kraus, John D. Monnier, Michael J. Ireland, Gaspard Duchene, Catherine Espaillat, Sebastian Hoenig, Attila Juhasz, Chris Mordasini, Johan Olofsson, Claudia Paladini, Keivan Stassun, Neal Turner, Gautam Vasisht, Tim J. Harries, Matthew R. Bate, Jean-Francois Gonzalez, Alexis Matter, Zhaohuan Zhu, Olja Panic, Zsolt Regaly, Alessandro Morbidelli, Farzana Meru, Sebastian Wolf, John Ilee, Jean-Philippe Berger , et al. (53 additional authors not shown)

Abstract: 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 prot… ▽ More 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. △ Less

Submitted 16 August, 2016; v1 submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 5 figures, Proceedings of SPIE 2016

Particle dynamics in discs with turbulence generated by the vertical shear instability

Authors: Moritz H. R. Stoll, Wilhelm Kley

Abstract: Among the candidates for generating turbulence in accretion discs in situations with low intrinsic ionization the vertical shear instability (VSI) has become an interesting candidate, as it relies purely on a vertical gradient in the angular velocity. Existing simulations have shown that $α$-values a few times $10^{-4}$ can be generated. The particle growth in the early planet formation phase is d… ▽ More Among the candidates for generating turbulence in accretion discs in situations with low intrinsic ionization the vertical shear instability (VSI) has become an interesting candidate, as it relies purely on a vertical gradient in the angular velocity. Existing simulations have shown that $α$-values a few times $10^{-4}$ can be generated. The particle growth in the early planet formation phase is determined by the dynamics of dust particles. Here, we address in particular the efficiency of VSI-turbulence in concentrating particles in order to generate overdensities and low collision velocities. We perform 3D numerical hydrodynamical simulations of accretion discs around young stars that include radiative transport and irradiation from the central star. The motion of particles within a size range of a fraction of mm up to several m is followed using standard drag formula. We confirm that under realistic conditions the VSI is able to generate turbulence in full 3D protoplanetary discs. The irradiated disc shows turbulence within 10 to 60au. The mean radial motion of the gas is such that it is directed inward near the midplane and outward in the surface layers. We find that large particles drift inward with the expected speed, while small particles can experience phases of outward drift. Additionally, the particles show bunching behaviour with overdensities reaching 5 times the average value, which is strongest for dimensionless stopping times around unity. Particles in a VSI-turbulent discs are concentrated in large scale turbulent eddies and show low relative speeds that allow for growing collisions. The reached overdensities will also allow for the onset streaming instabilities further enhancing particle growth. The outward drift for small particles at higher disk elevations allows for the transport of processed high temperature material in the Solar System to larger distances. △ Less

Submitted 8 July, 2016; originally announced July 2016.

Journal ref: A&A 594, A57 (2016)

A smooth particle hydrodynamics code to model collisions between solid, self-gravitating objects

Authors: Christoph M. Schäfer, Sven Riecker, Thomas I. Maindl, Roland Speith, Samuel Scherrer, Wilhelm Kley

Abstract: Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth… ▽ More Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request. △ Less

Submitted 12 April, 2016; originally announced April 2016.

Comments: 18 pages, 15 figures, accepted for publication in A&A

On the observability of bow shocks of Galactic runaway OB stars

Authors: D. M. -A. Meyer, A. -J. van Marle, R. Kuiper, W. Kley

Abstract: Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of dens… ▽ More Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of densities ranging from n_ISM=0.01 up to 10.0/cm3. Independently of their location in the Galaxy, we confirm that the infrared is the most appropriated waveband to search for bow shocks from massive stars. Their spectral energy distribution is the convenient tool to analyze them since their emission does not depend on the temporary effects which could affect unstable, thin-shelled bow shocks. Our numerical models of Galactic bow shocks generated by high-mass (~40 Mo) runaway stars yield H$α$ fluxes which could be observed by facilities such as the SuperCOSMOS H-Alpha Survey. The brightest bow shock nebulae are produced in the denser regions of the ISM. We predict that bow shocks in the field observed at Ha by means of Rayleigh-sensitive facilities are formed around stars of initial mass larger than about 20 Mo. Our models of bow shocks from OB stars have the emission maximum in the wavelength range 3 <= lambda <= 50 micrometer which can be up to several orders of magnitude brighter than the runaway stars themselves, particularly for stars of initial mass larger than 20 Mo. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 13 pages, 12 figures. Accepted to MNRAS (2016)

Dynamical friction for supersonic motion in a homogeneous gaseous medium

Authors: Daniel Thun, Rolf Kuiper, Franziska Schmidt, Wilhelm Kley

Abstract: The supersonic motion of gravitating objects through a gaseous medium constitutes a classical problem in theoretical astrophysics. Its application covers a broad range of objects and scales from planets up to galaxies. Especially the dynamical friction, caused by the forming wake behind the object, plays an important role for the dynamics of the system. To calculate the dynamical friction, standar… ▽ More The supersonic motion of gravitating objects through a gaseous medium constitutes a classical problem in theoretical astrophysics. Its application covers a broad range of objects and scales from planets up to galaxies. Especially the dynamical friction, caused by the forming wake behind the object, plays an important role for the dynamics of the system. To calculate the dynamical friction, standard formulae, based on linear theory are often used. It is our goal to check the general validity of these formulae and provide suitable expressions for the dynamical friction acting on the moving object, based on the basic physical parameters of the problem. We perform sequences of high resolution numerical studies of rigid bodies moving supersonically through a homogeneous medium, and calculate the total drag acting on the object, which is the sum of gravitational and hydro drag. We study cases without gravity with purely hydrodynamical drag, as well as gravitating objects. From the final equilibrium state of the simulations, we compute for gravitating objects the dynamical friction by direct numerical integration of the gravitational pull acting on the embedded object. The numerical experiments confirm the known scaling laws for the dependence of the dynamical friction on the basic physical parameters as derived in earlier semi-analytical studies. As a new important result we find that the shock's stand-off distance is revealed as the minimum spatial interaction scale of dynamical friction. Below this radius, the gas settles into a hydrostatic state, which causes no net gravitational pull onto the moving body. Finally, we derive an analytic estimate for the stand-off distance that can be used to calculate the dynamical friction force. △ Less

Submitted 28 January, 2016; originally announced January 2016.

Comments: 17 pages, 30 figures

Journal ref: A&A 589, A10 (2016)