Researcher profile

Willy Kley

Willy Kley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2013arXiv

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

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

preprint2012arXiv

Stellar irradiated discs and implications on migration of embedded planets I: equilibrium discs

The strength and direction of migration of low mass planets depends on the disc's thermodynamics. In discs where the viscous heating is balanced by radiative transport, the migration can be directed outwards, a process which extends the lifetime of growing planetary embryos. We investigate the influence of opacity and stellar irradiation on the disc thermodynamics. Utilizing the resulting disc structure, we determine the regions of outward migration. We perform two-dimensional numerical simulations of equilibrium discs with viscous heating, radiative cooling and stellar irradiation. We use the hydrodynamical code NIRVANA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration is studied by using torque formulae. In the constant opacity case, we reproduce the analytical results of a black-body disc: the stellar irradiation dominates in the outer regions -- leading to flaring -- while the viscous heating dominates close to the star. We find that the inner edge of the disc should not be significantly puffed-up by the stellar irradiation. If the opacity depends on the local density and temperature, the structure of the disc is different, and several bumps in the aspect ratio H/r appear, due to transitions between different opacity regimes. The bumps in the disc can shield the outer disc from stellar irradiation. Stellar irradiation is an important factor for determining the disc structure and has dramatic consequences for the migration of embedded planets. Compared to discs with only viscous heating, a stellar irradiated disc features a much smaller region of outward migration for a range of planetary masses. This suggests that the region where the formation of giant planet cores takes place is smaller, which in turn might lead to a shorter growth phase.

preprint2011arXiv

Evolution of inclined planets in three-dimensional radiative discs

While planets in the solar system only have a low inclination with respect to the ecliptic there is mounting evidence that in extrasolar systems the inclination can be very high, at least for close-in planets. One process to alter the inclination of a planet is through planet-disc interactions. Recent simulations considering radiative transport have shown that the evolution of migration and eccentricity can strongly depend on the thermodynamic state of the disc. We extend previous studies to investigate the planet-disc interactions of fixed and moving planets on inclined and eccentric orbits. We also analyse the effect of the disc&#39;s thermodynamic properties on the orbital evolution of embedded planets in detail. The protoplanetary disc is modelled as a viscous gas where the internally produced dissipation is transported by radiation. For locally isothermal discs, we confirm previous results and find inclination damping and inward migration for planetary cores. For low inclinations i < 2 H/r, the damping is exponential, while di/dt is proportional to i^-2 for larger i. For radiative discs, the planetary migration is very limited, as long as their inclination exceeds a certain threshold. If the inclination is damped below this threshold, planetary cores with a mass up to approximately 33 Earth masses start to migrate outwards, while larger cores migrate inwards right from the start. The inclination is damped for all analysed planet masses. In a viscous disc an initial inclination of embedded planets will be damped for all planet masses. This damping occurs on timescales that are shorter than the migration time. If the inclination lies beneath a certain threshold, the outward migration in radiative discs is not handicapped. Outward migration is strongest for circular and non-inclined orbits.

preprint2011arXiv

Orbital Migration Models under Test

Planet-disk interaction predicts a change in the orbital elements of an embedded planet. Through linear and fully hydrodynamical studies it has been found that migration is typically directed inwards. Hence, this migration process gives natural explanation for the presence of the &#39;hot&#39; planets orbiting close to the parent star, and it plays a mayor role in explaining the formation of resonant planetary systems. However, standard migration models for locally isothermal disks indicate a too rapid inward migration for small mass planets, and a large number of massive planets are found very far away from the star. Recent studies, including more complete disk physics, have opened up new paths to slow down or even reverse migration. The new findings on migration are discussed and connected to the observational properties of planetary systems.

preprint2011arXiv

Range of outward migration and influence of the disc&#39;s mass on the migration of giant planet cores

The migration of planets plays an important role in the early planet-formation process. An important problem has been that standard migration theories predict very rapid inward migration, which poses problems for population synthesis models. However, it has been shown recently that low-mass planets (20-30 Earth Masses) that are still embedded in the protoplanetary disc can migrate outwards under certain conditions. Simulations have been performed mostly for planets at given radii for a particular disc model. Here, we plan to extend previous work and consider different masses of the disc to quantify the influence of the physical disc conditions on planetary migration. We perform three-dimensional (3D) radiation hydrodynamical simulations of embedded planets in protoplanteary discs. For planets on circular orbits at various locations we measure the radial dependece of the torques. For all considered planet masses (20-30 Earth masses) in this study we find outward migration within a limited radial range of the disc, typically from about 0.5 up to 1.5-2.5 a_Jup. Inside and outside this intervall, migration is inward and given by the Lindblad value for large radii. Because outward migration stops at a certain location in the disc, there exists a zero-torque distance for planetary embryos, where they can continue to grow without moving too fast. For higher disc masses (M_disc > 0.02 M_Sol) convection ensues, which changes the structure of the disc and therefore the torque on the planet as well. Outward migration stops at different points in the disc for different planetary masses, resulting in a quite extended region where the formation of larger cores might be easier. In higher mass discs, convection changes the disc&#39;s structure resulting in fluctuations in the surface density, which influence the torque acting on the planet, and therefore its migration rate.

preprint2010arXiv

Orbital evolution of eccentric planets in radiative discs

With an average eccentricity of about 0.29, the eccentricity distribution of extrasolar planets is markedly different from the solar system. Among other scenarios considered, it has been proposed that eccentricity may grow through planet-disc interaction. Recently, it has been noticed that the thermodynamical state of the disc can significantly influence the migration properties of growing protoplanets. However, the evolution of planetary eccentricity in radiative discs has not been considered yet. In this paper we study the evolution of planets on eccentric orbits that are embedded in a three-dimensional viscous disc and analyse the disc&#39;s effect on the orbital evolution of the planet. We use the three-dimensional hydrodynamical code NIRVANA that includes full tensor viscosity and implicit radiation transport in the flux-limited diffusion approximation. The code uses the FARGO-algorithm to speed up the simulations. First we measure the torque and power exerted on the planet by the disc for fixed orbits, and then we let the planet start with initial eccentricity and evolve it in the disc. For locally isothermal we confirm previous results and find eccentricity damping and inward migration for planetary cores. In the case of radiative discs, the planets experience an inward migration as long as its eccentricity lies above a certain threshold. After the damping of eccentricity cores with masses below 33 Earthmasses begin to migrate outward in radiative discs, while higher mass cores always migrate inward. For all planetary masses studied (up to 200 Earthmasses) we find eccentricity damping. In viscous discs the orbital eccentricity of embedded planets is damped during the evolution independent of the mass. Hence, planet-disc interaction does not seem to be a viable mechanism to explain the observed high eccentricity of exoplanets.

preprint2010arXiv

Variable accretion as a mechanism for brightness variations in T Tau S

(Note: this is a shortened version of the original A&A-style structured abstract). The physical nature of the strong photometric variability of T Tau Sa, the more massive member of the Southern &#34;infrared companion&#34; to T Tau, has long been debated. Intrinsic luminosity variations due to variable accretion were originally proposed but later challenged in favor of apparent fluctuations due to time-variable foreground extinction. In this paper we use the timescale of the variability as a diagnostic for the underlying physical mechanism. Because the IR emission emerging from Sa is dominantly thermal emission from circumstellar dust at <=1500K, we can derive a minimum size of the region responsible for the time-variable emission. In the context of the variable foreground extinction scenario, this region must be (un-) covered within the variability timescale, which implies a minimum velocity for the obscuring foreground material. If this velocity supercedes the local Kepler velocity we can reject foreground extinction as a valid variability mechanism. The variable accretion scenario allows for shorter variability timescales since the variations in luminosity occur on much smaller scales, essentially at the surface of the star, and the disk surface can react almost instantly on the changing irradiation with a higher or lower dust temperature and according brightness. We have detected substantial variations at long wavelengths in T Tau S: +26% within four days at 12.8 micron. We show that this short-term variability cannot be due to variable extinction and instead must be due to variable accretion. Using a radiative transfer model of the Sa disk we show that variable accretion can in principle also account for the much larger (several magnitude) variations observed on timescales of several years. For the long-term variability, however, also variable foreground extinction is a viable mechanism.

preprint2009arXiv

Formation of massive planets in binary star systems

As of today over 40 planetary systems have been discovered in binary star systems. In all cases the configuration appears to be circumstellar, where the planets orbit around one of the stars, the secondary acting as a perturber. The formation of planets in binary star systems is more difficult than around single stars due to the gravitational action of the companion on the dynamics of the protoplanetary disk. In this contribution we first briefly present the relevant observational evidence for planets in binary systems. Then the dynamical influence that a secondary companion has on a circumstellar disk will be analyzed through fully hydrodynamical simulations. We demonstrate that the disk becomes eccentric and shows a coherent precession around the primary star. Finally, fully hydrodynamical simulations of evolving protoplanets embedded in disks in binary star systems are presented. We investigate how the orbital evolution of protoplanetary embryos and their mass growth from cores to massive planets might be affected in this very dynamical environment. We consider, in particular, the planet orbiting the primary in the system Gamma Cephei.

preprint2009arXiv

Planet migration in three-dimensional radiative discs

The migration of growing protoplanets depends on the thermodynamics of the ambient disc. Standard modelling, using locally isothermal discs, indicate in the low planet mass regime an inward (type-I) migration. Taking into account non-isothermal effects, recent studies have shown that the direction of the type-I migration can change from inward to outward. In this paper we extend previous two-dimensional studies, and investigate the planet-disc interaction in viscous, radiative discs using fully three-dimensional radiation hydrodynamical simulations of protoplanetary accretion discs with embedded planets, for a range of planetary masses. We use an explicit three-dimensional (3D) hydrodynamical code NIRVANA that includes full tensor viscosity. We have added implicit radiation transport in the flux-limited diffusion approximation, and to speed up the simulations significantly we have newly adapted and implemented the FARGO-algorithm in a 3D context. First, we present results of test simulations that demonstrate the accuracy of the newly implemented FARGO-method in 3D. For a planet mass of 20 M_earth we then show that the inclusion of radiative effects yields a torque reversal also in full 3D. For the same opacity law used the effect is even stronger in 3D than in the corresponding 2D simulations, due to a slightly thinner disc. Finally, we demonstrate the extent of the torque reversal by calculating a sequence of planet masses. Through full 3D simulations of embedded planets in viscous, radiative discs we confirm that the migration can be directed outwards up to planet masses of about 33 M_earth. Hence, the effect may help to resolve the problem of too rapid inward migration of planets during their type-I phase.

preprint2003arXiv

Orbital Migration and Mass Accretion of Protoplanets in 3D Global Computations with Nested Grids

We investigate the evolution of protoplanets with different masses embedded in an accretion disk, via global fully three-dimensional hydrodynamical simulations. We consider a range of planetary masses extending from one and a half Earth&#39;s masses up to one Jupiter&#39;s mass, and we take into account physically realistic gravitational potentials of forming planets. In order to calculate accurately the gravitational torques exerted by disk material and to investigate the accretion process onto the planet, the flow dynamics has to be thoroughly resolved on long as well as short length scales. We achieve this strict resolution requirement by applying a nested-grid refinement technique which allows to greatly enhance the local resolution. Our results from altogether 51 simulations show that for large planetary masses, approximately above a tenth of the Jupiter&#39;s mass, migration rates are relatively constant, as expected in type II migration regime and in good agreement with previous two-dimensional calculations. In a range between seven and fifteen Earth&#39;s masses, we find a dependency of the migration speed on the planetary mass that yields time scales considerably longer than those predicted by linear analytical theories. This property may be important in determining the overall orbital evolution of protoplanets. The growth time scale is minimum around twenty Earth-masses, but it rapidly increases for both smaller and larger mass values. Significant differences between two- and three-dimensional calculations are found in particular for objects with masses smaller than ten Earth-masses. We also derive an analytical approximation for the numerically computed mass growth rates.

preprint2003arXiv

Thermo-Hydrodynamics of Circumstellar Disks with High-mass Planets

With a series of numerical simulations, we analyze the thermo-hydrodynamical evolution of circumstellar disks containing Jupiter-size protoplanets. In the framework of the two-dimensional approximation, we consider an energy equation that includes viscous heating and radiative effects in a simplified, yet consistent form. Multiple nested grids are used in order to study both global and local features around the planet. By means of different viscosity prescriptions, we investigate various temperature regimes. A planetary mass range from 0.1 to 1 Mj is examined. Computations show that gap formation is a general property which affects density, pressure, temperature, optical thickness, and radiated flux distributions. However, it remains a prominent feature only when the kinematic viscosity is on the order of 10^(15) cm^2/s or lower. Around accreting planets, a circumplanetary disk forms that has a surface density profile decaying exponentially with the distance and whose mass is 5-6 orders of magnitudes smaller than Jupiter&#39;s mass. Circumplanetary disk temperature profiles decline roughly as the inverse of the distance from the planet. Temperatures range from some 10 to ~1000 K. Planetary accretion and migration rates depend on the viscosity regime, with discrepancies within an order of magnitude. Estimates of growth and migration time scales inferred by these models are on the same orders of magnitude as those previously obtained with locally isothermal simulations.

preprint1999arXiv

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter sized planets is studied numerically. The disc is assumed to be flat and non-self gravitating, which is modeled by the planar (two-dimensional) Navier-Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass $M_{Jup}$ each and the radial distances from the star are one and two semi-major axis of Jupiter, respectively. During the evolution both planets increase their mass due to accretion of gas from the disc; after about 2500 orbital periods of the inner planet it has reached a mass of 2.3 and the outer planet of 3.2~$M_{Jup}$. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc, while the semi-major axis of the inner planet remains constant. When the distance of close approach eventually becomes smaller than their mutual Hill radius the eccentricities increase strongly and the system may turn unstable. The implications for the origin of the solar system and the extrasolar planets are discussed.

preprint1999arXiv

The Migration and Growth of Protoplanets in Protostellar Discs

We investigate the gravitational interaction of a Jovian mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions have been investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Nonlinear hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near circular orbit after a characteristic viscous time scale of approximately 10,000 initial orbital periods. This was found to be independent of whether the protoplanet was allowed to accrete mass or not. Inward migration is helped through the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion onto the central star.Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is approximately 5 Jupiter masses, because of the reduced accretion rates obtained for planets of increasing mass. Assuming some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets (Marcy, Cochran, & Mayor 1999; Marcy & Butler 1998) as well as their inward orbital migration can be accounted for by consideration of disc--protoplanet interactions during the late stages of giant planet formation. Maximally accreting protoplanets reached about four Jovian masses on reaching the neighbourhood of the central star.

preprint1998arXiv

Mass Flow and Accretion through gaps in Accretion Discs

We study the structure and dynamics of the gap created by a protoplanet in an accretion disc. The hydrodynamic equations for a flat, two-dimensional, non-selfgravitating protostellar accretion disc with an embedded, Jupiter sized protoplanet on a circular orbit are solved. To simulate possible accretion of mass onto the protoplanet we continually remove mass from the interior of the planet&#39;s Roche lobe which is monitored. Firstly, it is shown that consistent results independent on numerical issues (such as boundary or initial conditions, artificial viscosity or resolution) can be obtained. Then, a detailed parameter study delineates the influence of the disc viscosity and pressure on the magnitude of the accretion rate. We find that, even after the formation of a gap in the disc, the planet is still able to accrete more mass from the disc. This accretion occurs from regions of the disc which are radially exterior and interior to the planet&#39;s orbital radius. The rate depends on the magnitude of the viscosity and vertical thickness of the disc. For a disc viscosity alpha=10^{-3} and vertical thickness H/r=0.05 we estimate the time scale for the accumulation of one Jupiter mass to be of order hundred thousand years. For a larger(smaller) viscosity and disc thickness this accretion rate is increasing(decreasing). For a very small viscosity (alpha < 5 10^{-4}) the mass accretion rate through the gap onto the planet is markedly reduced, and the corresponding accretion time scale becomes larger than the viscous evolution time of the disc.

preprint1997arXiv

On relativistic discs and rings

Sequences of infinitesimally thin, uniformly rotating, self-gravitating relativistic discs with internal two-dimensional pressure have been constructed. It is shown that in weaker relativistic configurations the sequences undergo a continuous bifurcation from a disc to a ring structure, while in stronger relativistic cases the sequences terminate at the mass-shed limit where gravitational forces are exactly balanced by centrifugal forces.