-
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.
-
BOWIE-ALIGN: How formation and migration histories of giant planets impact atmospheric compositions
Authors:
Anna B. T. Penzlin,
Richard A. Booth,
James Kirk,
James E. Owen,
Eva-Maria Ahrer,
Duncan A. Christie,
Alastair B. Claringbold,
Emma Esparza-Borges,
M. López-Morales,
N. J. Mayne,
Mason McCormack,
Annabella Meech,
Vatsal Panwar,
Diana Powell,
Denis E. Sergeev,
Jake Taylor,
Peter J. Wheatley,
Maria Zamyatina
Abstract:
Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the composition…
▽ More
Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the compositions of 8 hot Jupiters around F stars, 4 with orbits aligned with the stellar rotation axis and 4 misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to create a full, unbiased model sample from which we can estimate the range of final atmospheric compositions. For high metallicity atmospheres (O/H > 10 times solar), the C/O ratios of aligned and misaligned planets diverge, with aligned planets having lower C/O (< 0.25) due to the accretion of oxygen-rich silicates from the inner disc. However, silicates may rain out instead of releasing their oxygen into the atmosphere. This would significantly increase the C/O of aligned planets (C/O > 0.6), inverting the trend between the aligned and misaligned planets. Nevertheless, by comparing statistically significant samples of aligned and misaligned planets, we expect atmospheric composition to constrain how planets form.
△ Less
Submitted 4 July, 2024; v1 submitted 3 July, 2024;
originally announced July 2024.
-
BOWIE-ALIGN: A JWST comparative survey of aligned vs misaligned hot Jupiters to test the dependence of atmospheric composition on migration history
Authors:
James Kirk,
Eva-Maria Ahrer,
Anna B. T. Penzlin,
James E. Owen,
Richard A. Booth,
Lili Alderson,
Duncan A. Christie,
Alastair B. Claringbold,
Emma Esparza-Borges,
Chloe E. Fisher,
Mercedes López-Morales,
N. J. Mayne,
Mason McCormack,
Annabella Meech,
Vatsal Panwar,
Diana Powell,
Jake Taylor,
Denis E. Sergeev,
Daniel Valentine,
Hannah R. Wakeford,
Peter J. Wheatley,
Maria Zamyatina
Abstract:
A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft…
▽ More
A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft break with different orbital alignments. Under the assumption that aligned planets migrate through the inner disc, while misaligned planets migrate after disc dispersal, the act of migrating through the inner disc should lead to a measurable difference in the C/O between aligned and misaligned planets. We expect the amplitude and sign of this difference to depend on the amount of planetesimal accretion and whether silicates accreted from the inner disc release their oxygen. Here, we identify all known exoplanets that are suitable for testing this hypothesis, describe our JWST survey, and use noise simulations and atmospheric retrievals to estimate our survey's sensitivity. With the selected sample of four aligned and four misaligned hot Jupiters, we will be sensitive to the predicted differences in C/O between aligned and misaligned hot Jupiters for a wide range of model scenarios.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
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.
-
The Santa Barbara Binary-Disk Code Comparison
Authors:
Paul C. Duffell,
Alexander J. Dittmann,
Daniel J. D'Orazio,
Alessia Franchini,
Kaitlin M. Kratter,
Anna B. T. Penzlin,
Enrico Ragusa,
Magdalena Siwek,
Christopher Tiede,
Haiyang Wang,
Jonathan Zrake,
Adam M. Dempsey,
Zoltan Haiman,
Alessandro Lupi,
Michal Pirog,
Geoffrey Ryan
Abstract:
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed…
▽ More
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.
△ Less
Submitted 20 June, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
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.
-
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.
-
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.
-
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.
-
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.
-
On moving shadows and pressure bumps in HD 169142
Authors:
Gesa H. -M. Bertrang,
Mario Flock,
Miriam Keppler,
Trifon Trifonov,
Anna B. T. Penzlin,
Henning Avenhaus,
Thomas Henning,
Matias Montesinos
Abstract:
The search for young planets had its first breakthrough with the detection of the accreting planet PDS70b. In this study, we aim to broaden our understanding towards the formation of multi-planet systems such as HR8799 or the Solar System. Our previous study on HD169142, one of the closest Herbig stars, points towards a shadow-casting protoplanetary candidate. Here, we present follow-up observatio…
▽ More
The search for young planets had its first breakthrough with the detection of the accreting planet PDS70b. In this study, we aim to broaden our understanding towards the formation of multi-planet systems such as HR8799 or the Solar System. Our previous study on HD169142, one of the closest Herbig stars, points towards a shadow-casting protoplanetary candidate. Here, we present follow-up observations to test our previously proposed hypothesis. We set our new data into context with previous observations to follow structural changes in the disk over the course of 6 years. We find spatially resolved systematic changes in the position of the previously described surface brightness dip in the inner ring. We further find changes in the brightness structure in azimuthal direction along the ring. And finally, a comparison of our SPHERE data with recent ALMA observations reveals a wavelength dependent radial profile of the bright ring. The time-scale on which the changes in the ring's surface brightness occur suggest that they are caused by a shadow cast by a 1-10Mj planet surrounded by dust, an orbit comparable to those of the giant planets in our own Solar System. Additionally, we find the first indications for temperature-induced instabilities in the ring. And finally, we trace a pressure maxima, for the first time spatially resolved, with a width of 4.5au. The density distribution of the ring at mm wavelengths around the pressure maxima could further indicate effects from snow lines or even the dynamics and feedback of the larger grains.
△ Less
Submitted 29 April, 2021; v1 submitted 22 July, 2020;
originally announced July 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.
-
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.
-
Dust-driven viscous ring-instability in protoplanetary disks
Authors:
C. P. Dullemond,
A. B. T. Penzlin
Abstract:
Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. Th…
▽ More
Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. The presence of copious amounts of dust tends to strongly reduce the conductivity of the gas, thereby inhibiting the magneto-rotational instability, and thus reducing the turbulence in the disk. From viscous disk theory it is known that a disk tends to increase its surface density in regions where the viscosity (i.e. turbulence) is low. Local maxima in the gas pressure tend to attract dust through radial drift, increasing the dust content even more. We investigate mathematically if this could potentially lead to a feedback loop in which a perturbation in the dust surface density could perturb the gas surface density, leading to increased dust drift and thus amplification of the dust perturbation and, as a consequence, the gas perturbation. We find that this is indeed possible, even for moderately small dust grain sizes, which drift less efficiently, but which are more likely to affect the gas ionization degree. We speculate that this instability could be triggered by the small dust population initially, and when the local pressure maxima are strong enough, the larger dust grains get trapped and lead to the familiar ring-like shapes. We also discuss the many uncertainties and limitations of this model.
△ Less
Submitted 27 November, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.
