-
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.
-
On the effects of photoionization feedback on second-generation star formation in globular clusters of different masses
Authors:
Asiyeh Yaghoobi,
Joakim Rosdahl,
Francesco Calura,
Sareh Ataiee
Abstract:
We simulate the formation of second-generation stars in young clusters with masses of $10^{5}$ and $10^6\ {\rm M_\odot}$ within $30-100\ {\rm Myr}$ after the formation of clusters. We assume the clusters move through a uniform interstellar medium with gas densities of $10^{-24}$ and $10^{-23}\ {\rm gcm}^{-3} $ and consider the stellar winds from asymptotic giant branch (AGB) stars, gas accretion o…
▽ More
We simulate the formation of second-generation stars in young clusters with masses of $10^{5}$ and $10^6\ {\rm M_\odot}$ within $30-100\ {\rm Myr}$ after the formation of clusters. We assume the clusters move through a uniform interstellar medium with gas densities of $10^{-24}$ and $10^{-23}\ {\rm gcm}^{-3} $ and consider the stellar winds from asymptotic giant branch (AGB) stars, gas accretion onto the cluster, ram pressure, star formation, and photoionization feedback of our stellar systems including binary stars. We find that second-generation (SG) stars can be formed only within the $10^6\ {\rm M_\odot}$ cluster in the high-density simulation, where the cluster can accrete sufficient pristine gas from their surrounding medium, leading to efficient cooling required for the ignition of SG formation and sufficient dilution of the AGB ejecta. Hence, our results indicate that a denser environment is another requirement for the AGB scenario to explain the presence of multiple populations in globular clusters. On the other hand, the ionizing feedback becomes effective in heating the gas in our low-density simulations. As a result, the clusters cannot accumulate a considerable amount of pristine gas at their center. The gas mass within the clusters in these simulations is similar to that in young massive clusters (YMCs). Hence, our studies can provide a possible reason for the lack of gas, star formation, and SG stars in YMCs. Our results indicate that the ionizing stellar feedback is not a severe problem for SG formation; rather, it can help the AGB scenario to account for some observables.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
TOI-1130: A photodynamical analysis of a hot Jupiter in resonance with an inner low-mass planet
Authors:
J. Korth,
D. Gandolfi,
J. Šubjak,
S. Howard,
S. Ataiee,
K. A. Collins,
S. N. Quinn,
A. J. Mustill,
T. Guillot,
N. Lodieu,
A. M. S. Smith,
M. Esposito,
F. Rodler,
A. Muresan,
L. Abe,
S. H. Albrecht,
A. Alqasim,
K. Barkaoui,
P. G. Beck,
C. J. Burke,
R. P. Butler,
D. M. Conti,
K. I. Collins,
J. D. Crane,
F. Dai
, et al. (37 additional authors not shown)
Abstract:
The TOI-1130 is a known planetary system around a K-dwarf consisting of a gas giant planet, TOI-1130 c, on an 8.4-day orbit, accompanied by an inner Neptune-sized planet, TOI-1130 b, with an orbital period of 4.1 days. We collected precise radial velocity (RV) measurements of TOI-1130 with the HARPS and PFS spectrographs as part of our ongoing RV follow-up program. We perform a photodynamical mode…
▽ More
The TOI-1130 is a known planetary system around a K-dwarf consisting of a gas giant planet, TOI-1130 c, on an 8.4-day orbit, accompanied by an inner Neptune-sized planet, TOI-1130 b, with an orbital period of 4.1 days. We collected precise radial velocity (RV) measurements of TOI-1130 with the HARPS and PFS spectrographs as part of our ongoing RV follow-up program. We perform a photodynamical modeling of the HARPS and PFS RVs, and transit photometry from the Transiting Exoplanet Survey Satellite (TESS) and the TESS Follow-up Observing Program. We determine the planet masses and radii of TOI-1130 b and TOI-1130 c to be Mb = 19.28 $\pm$ 0.97 M$_\oplus$ and Rb = 3.56 $\pm$ 0.13 R$_\oplus$, and Mc = 325.59 $\pm$ 5.59 M$_\oplus$ and Rc = 13.32+1.55-1.41 R$_\oplus$, respectively. We spectroscopically confirm TOI-1130 b that was previously only validated. We find that the two planets orbit with small eccentricities in a 2:1 resonant configuration. This is the first known system with a hot Jupiter and an inner lower mass planet locked in a mean-motion resonance. TOI-1130 belongs to the small yet increasing population of hot Jupiters with an inner low-mass planet that challenges the pathway for hot Jupiter formation. We also detect a linear RV trend possibly due to the presence of an outer massive companion.
△ Less
Submitted 24 May, 2023;
originally announced May 2023.
-
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.
-
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.
-
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.
-
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.
-
Stability of the co-orbital resonance under dissipation: Application to its evolution in protoplanetary discs
Authors:
Adrien Leleu,
Gavin Coleman,
Sareh Ataiee
Abstract:
Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. In this paper we investigate how a pair of co-orbital exoplanets would fare during their migration in a protoplanetary disc. To this end, we computed a stability cr…
▽ More
Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. In this paper we investigate how a pair of co-orbital exoplanets would fare during their migration in a protoplanetary disc. To this end, we computed a stability criterion of the Lagrangian equilibria L4 and L5 under generic dissipation and slow mass evolution. Depending on the strength and shape of these perturbations, the system can either evolve towards the Lagrangian equilibrium, or tend to increase its amplitude of libration, possibly all the way to horseshoe orbits or even exiting the resonance. We estimated the various terms of our criterion using a set of hydrodynamical simulations, and show that the dynamical coupling between the disc perturbations and both planets have a significant impact on the stability: the structures induced by each planet in the disc perturb the dissipative forces applied on the other planets over each libration cycle. Amongst our results on the stability of co-orbitals, several are of interest to constrain the observability of such configurations: long-distance inward migration and smaller leading planets tend to increase the libration amplitude around the Lagrangian equilibria, while leading massive planets and belonging to a resonant chain tend to stabilise it. We also show that, depending on the strength of the dissipative forces, both the inclination and the eccentricity of the smaller of the two co-orbitals can be significantly increased during the inward migration of the co-orbital pair, which can have a significant impact on the detectability by transit of such configurations.
△ Less
Submitted 19 November, 2019; v1 submitted 22 January, 2019;
originally announced January 2019.
-
The formation of Jupiter by hybrid pebble-planetesimal accretion
Authors:
Yann Alibert,
Julia Venturini,
Ravit Helled,
Sareh Ataiee,
Remo Burn,
Luc Senecal,
Willy Benz,
Lucio Mayer,
Christoph Mordasini,
Sascha P. Quanz,
Maria Schönbächler
Abstract:
The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in me…
▽ More
The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in meteorites provided evidence for the existence of two reservoirs in the early Solar System. These reservoirs remained separated from ~1 until ~ 3 Myr after the beginning of the Solar System's formation. This separation is interpreted as resulting from Jupiter growing and becoming a barrier for material transport. In this framework, Jupiter reached ~20 Earth masses within ~1 Myr and slowly grew to ~50 Earth masses in the subsequent 2 Myr before reaching its present-day mass. The evidence that Jupiter slowed down its growth after reaching 20 Earth masses for at least 2 Myr is puzzling because a planet of this mass is expected to trigger fast runaway gas accretion. Here, we use theoretical models to describe the conditions allowing for such a slow accretion and show that Jupiter grew in three distinct phases. First, rapid pebble accretion brought the major part of Jupiter's core mass. Second, slow planetesimal accretion provided the energy required to hinder runaway gas accretion during 2 Myr. Third, runaway gas accretion proceeded. Both pebbles and planetesimals therefore have an important role in Jupiter's formation.
△ Less
Submitted 14 September, 2018;
originally announced September 2018.
-
Metallicity effect and planet mass function in pebble-based planet formation models
Authors:
Natacha Brügger,
Yann Alibert,
Sareh Ataiee,
Willy Benz
Abstract:
One of the main scenarios of planet formation is the core accretion model where a massive core forms first and then accretes a gaseous envelope. This core forms by accreting solids, either planetesimals, or pebbles. A key constraint in this model is that the accretion of gas must proceed before the dissipation of the gas disc. Classical planetesimal accretion scenario predicts that the time needed…
▽ More
One of the main scenarios of planet formation is the core accretion model where a massive core forms first and then accretes a gaseous envelope. This core forms by accreting solids, either planetesimals, or pebbles. A key constraint in this model is that the accretion of gas must proceed before the dissipation of the gas disc. Classical planetesimal accretion scenario predicts that the time needed to form a giant planets core is much longer than the time needed to dissipate the disc. This difficulty led to the development of another accretion scenario, in which cores grow by accretion of pebbles, which are much smaller and thus more easily accreted, leading to a more rapid formation. The aim of this paper is to compare our updated pebble-based planet formation model with observations, in particular the well studied metallicity effect. We adopt the Bitsch et al. 2015a disc model and the Bitsch et al. 2015b pebble model and use a population synthesis approach to compare the formed planets with observations. We find that keeping the same parameters as in Bitsch et al. 2015b leads to no planet growth due to a computation mistake in the pebble flux (Bitsch et al. 2017). Indeed a large fraction of the heavy elements should be put into pebbles (Zpeb/Ztot = 0.9) in order to form massive planets using this approach. The resulting mass functions show a huge amount of giants and a lack of Neptune mass planets, which are abundant according to observations. To overcome this issue we include the computation of the internal structure for the planetary atmosphere to our model. This leads to the formation of Neptune mass planets but no observable giants. Reducing the opacity of the planetary envelope finally matches observations better. We conclude that modeling the internal structure for the planetary atmosphere is necessary to reproduce observations.
△ Less
Submitted 31 August, 2018;
originally announced August 2018.
-
How much does turbulence change the pebble isolation mass for planet formation?
Authors:
S. Ataiee,
C. Baruteau,
Y. Alibert,
W. Benz
Abstract:
When a planet becomes massive enough, it gradually carves a partial gap around its orbit in the protoplanetary disk. A pressure maximum can be formed outside the gap where solids that are loosely coupled to the gas, typically in the pebble size range, can be trapped. The minimum planet mass for building such a trap, which is called the pebble isolation mass (PIM), is important for two reasons: it…
▽ More
When a planet becomes massive enough, it gradually carves a partial gap around its orbit in the protoplanetary disk. A pressure maximum can be formed outside the gap where solids that are loosely coupled to the gas, typically in the pebble size range, can be trapped. The minimum planet mass for building such a trap, which is called the pebble isolation mass (PIM), is important for two reasons: it marks the end of planetary growth by pebble accretion, and the trapped dust forms a ring that may be observed with millimetre observations. We study the effect of disk turbulence on the pebble isolation mass and find its dependence on the gas turbulent viscosity, aspect ratio, and particles Stokes number. By means of 2D gas hydrodynamical simulations, we found the minimum planet mass to form a radial pressure maximum beyond the orbit of the planet, which is the necessary condition to trap pebbles. We then carried out 2D gas plus dust hydrodynamical simulations to examine how dust turbulent diffusion impacts particles trapping at the pressure maximum. We finally provide a semi-analytical calculation of the PIM based on comparing the radial drift velocity of solids and the root mean square turbulent velocity fluctuations around the pressure maximum. From our results of gas simulations, we provide an expression for the PIM versus disk aspect ratio and turbulent viscosity. Our gas plus dust simulations show that the effective PIM can be nearly an order of magnitude larger in high-viscosity disks because turbulence diffuse particles out of the pressure maximum. This is quantified by our semi-analytical calculation, which gives an explicit dependence of the PIM with Stokes number of particles. We conclude than disk turbulence can significantly alter the PIM, depending on the level of turbulence in regions of planet formation.
△ Less
Submitted 5 April, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
A Multi-Wavelength Analysis of Dust and Gas in the SR 24S Transition Disk
Authors:
P. Pinilla,
L. M. Pérez,
S. Andrews,
N. van der Marel,
E. F. van Dishoeck,
S. Ataiee,
M. Benisty,
T. Birnstiel,
A. Juhász,
A. Natta,
L. Ricci,
L. Testi
Abstract:
We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-l…
▽ More
We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the $^{13}$CO and C$^{18}$O (J=2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is $\lesssim 0.12\,M_{\bigoplus}$, which gives a disk mass ratio in dust between the two components of $M_{\mathrm{dust, SR\,24S}}/M_{\mathrm{dust, SR\,24N}}\gtrsim840$. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to mm-sizes is inhibited there and that only warm gas, as seen by ro-vibrational CO emission inside the truncation radii of the binary, is present.
△ Less
Submitted 27 March, 2017;
originally announced March 2017.
-
Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet
Authors:
A. Pohl,
P. Pinilla,
M. Benisty,
S. Ataiee,
A. Juhasz,
C. P. Dullemond,
R. Van Boekel,
T. Henning
Abstract:
Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isoth…
▽ More
Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isothermal hydrodynamical simulations including viscous heating and a cooling prescription are combined with three-dimensional dust continuum radiative transfer models for direct comparison to observations. We find that the resulting contrast between the spirals and the surrounding disc in scattered light is by far higher for pressure scale height variations, i.e. thermal perturbations, than for pure surface density variations. Self-gravity effects suppress any vortex modes and tend to reduce the opening angle of planet-induced spirals, making them more tightly wound. If the disc is only marginally gravitationally stable with a Toomre parameter around unity, an embedded massive planet (planet-to-star mass ratio of $10^{-2}$) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.
△ Less
Submitted 18 August, 2015;
originally announced August 2015.
-
Gas and dust structures in protoplanetary disks hosting multiple planets
Authors:
P. Pinilla,
M. de Juan Ovelar,
S. Ataiee,
M. Benisty,
T. Birnstiel,
E. F. van Dishoeck,
M. Min
Abstract:
Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of par…
▽ More
Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. We study the gas and dust distributions of a disk hosting two massive planets as function of different disk and dust parameters. Observational signatures, such as the spectral energy distribution, sub-millimetre, and polarised images are simulated for the various parameters. We confirm that planets can lead to particle trapping, although for a disk with high viscosity ($α_{\rm{turb}}=10^{-2}$), the planet should be more massive than $5 M_{\rm{Jup}}$ and dust fragmentation should occur with low efficiency ($v_{f}\sim30\rm{m s}^{-1}$). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity ($α_{\rm{turb}}=10^{-2}$) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity ($α_{\rm{turb}}=10^{-3}$), the planet mass is constrained to be less than $5 M_{\rm{Jup}}$. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.
△ Less
Submitted 18 November, 2014; v1 submitted 22 October, 2014;
originally announced October 2014.
-
Planet-vortex interaction:How a vortex can shepherd a planetary embryo
Authors:
S. Ataiee,
C. P. Dullemond,
W. Kley,
Zs. Regaly,
H. Meheut
Abstract:
Context: Anticyclonic vortices are considered as a favourable places for trapping dust and forming planetary embryos. On the other hand, they are massive blobs that can interact gravitationally with the planets in the disc. Aims: We aim to study how a vortex interacts gravitationally with a planet which migrates toward it or a planet which is created inside the vortex. Methods: We performed hydrod…
▽ More
Context: Anticyclonic vortices are considered as a favourable places for trapping dust and forming planetary embryos. On the other hand, they are massive blobs that can interact gravitationally with the planets in the disc. Aims: We aim to study how a vortex interacts gravitationally with a planet which migrates toward it or a planet which is created inside the vortex. Methods: We performed hydrodynamical simulations of a viscous locally isothermal disc using GFARGO and FARGO-ADSG. We set a stationary Gaussian pressure bump in the disc in a way that RWI is triggered. After a large vortex is established, we implanted a low mass planet in the outer disc or inside the vortex and allowed it to migrate. We also examined the effect of vortex strength on the planet migration and checked the validity of the final result in the presence of self-gravity. Results: We noticed regardless of the planet's initial position, the planet is finally locked to the vortex or its migration is stopped in a farther orbital distance in case of a stronger vortex. For the model with the weaker vortex, we studied the effect of different parameters such as background viscosity, background surface density, mass of the planet and different planet positions. In these models, while the trapping time and locking angle of the planet vary for different parameters, the main result, which is the planet-vortex locking, remains valid. We discovered that even a planet with a mass less than 5 * 10^{-7} M_{\star} comes out from the vortex and is locked to it at the same orbital distance. For a stronger vortex, both in non-self-gravitated and self-gravitating models, the planet migration is stopped far away from the radial position of the vortex. This effect can make the vortices a suitable place for continual planet formation under the condition that they save their shape during the planetary growth.
△ Less
Submitted 1 October, 2014;
originally announced October 2014.
-
Trapping of giant-planet cores - I. Vortex aided trapping at the outer dead zone edge
Authors:
Zs. Regaly,
Zs. Sandor,
P. Csomos,
S. Ataiee
Abstract:
In this paper the migration of a 10 Earth-mass planetary core is investigated at the outer boundary of the dead zone of a protoplanetary disc by means of 2D hydrodynamic simulations done with the graphics processor unit version of the FARGO code. In the dead zone, the effective viscosity is greatly reduced due to the disc self-shielding against stellar UV radiation, X-rays from the stellar magneto…
▽ More
In this paper the migration of a 10 Earth-mass planetary core is investigated at the outer boundary of the dead zone of a protoplanetary disc by means of 2D hydrodynamic simulations done with the graphics processor unit version of the FARGO code. In the dead zone, the effective viscosity is greatly reduced due to the disc self-shielding against stellar UV radiation, X-rays from the stellar magnetosphere and interstellar cosmic rays. As a consequence, mass accumulation occurs near the outer dead zone edge, which is assumed to trap planetary cores enhancing the efficiency of the core-accretion scenario to form giant planets. Contrary to the perfect trapping of planetary cores in 1D models, our 2D numerical simulations show that the trapping effect is greatly dependent on the width of the region where viscosity reduction is taking place. Planet trapping happens exclusively if the viscosity reduction is sharp enough to allow the development of large-scale vortices due to the Rossby wave instability. The trapping is only temporarily, and its duration is inversely proportional to the width of the viscosity transition. However, if the Rossby wave instability is not excited, a ring-like axisymmetric density jump forms, which cannot trap the 10 Earth-mass planetary cores. We revealed that the stellar torque exerted on the planet plays an important role in the migration history as the barycentre of the system significantly shifts away from the star due to highly non-axisymmetric density distribution of the disc. Our results still support the idea of planet formation at density/pressure maximum, since the migration of cores is considerably slowed down enabling them further growth and runaway gas accretion in the vicinity of an overdense region.
△ Less
Submitted 28 June, 2013; v1 submitted 24 May, 2013;
originally announced May 2013.
-
Asymmetric transition disks: Vorticity or eccentricity?
Authors:
S. Ataiee,
P. Pinilla,
A. Zsom,
C. P. Dullemond,
C. Dominik,
J. Ghanbari
Abstract:
Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet formation. It has…
▽ More
Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet formation. It has been shown that an elongated giant vortex produced in a disk with a strong viscosity jump strikingly resembles the observed asymmetric rings.
Aims. We aim to study a similar behavior for a disk in which a giant planet is embedded. However, a giant planet can induce two kinds of asymmetries: (1) a giant vortex, and (2) an eccentric disk. We studied under which conditions each of these can appear, and how one can observationally distinguish between them. This is important because only a vortex can trap particles both radially and azimuthally, while the eccentric ring can only trap particles in radial direction.
Methods. We used the FARGO code to conduct the hydro-simulations. We set up a disk with an embedded giant planet and took a radial grid spanning from 0.1 to 7 times the planet semi-major axis. We ran the simulations with various viscosity values and planet masses for 1000 planet orbits to allow a fully developed vortex or disk eccentricity. Afterwards, we compared the dust distribution in a vortex-holding disk with an eccentric disk using dust simulations.
Results. We find that vorticity and eccentricity are distinguishable by looking at the azimuthal contrast of the dust density. While vortices, as particle traps, produce very pronounced azimuthal asymmetries, eccentric features are not able to accumulate millimeter dust particles in azimuthal direction, and therefore the asymmetries are expected to be modest.
△ Less
Submitted 5 April, 2013;
originally announced April 2013.