BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Min-Kai Lin

Min-Kai Lin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAMachine Learningphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

29 published item(s)

preprint2022arXiv

Revealing the drag instability in one-fluid nonideal MHD simulations of a 1D isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the postshock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Using Bayesian Deep Learning to infer Planet Mass from Gaps in Protoplanetary Disks

Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its predictions. In this paper, we introduce a Bayesian deep learning network "DPNNet-Bayesian" that can predict planet mass from disk gaps and provides uncertainties associated with the prediction. A unique feature of our approach is that it can distinguish between the uncertainty associated with the deep learning architecture and uncertainty inherent in the input data due to measurement noise. The model is trained on a data set generated from disk-planet simulations using the \textsc{fargo3d} hydrodynamics code with a newly implemented fixed grain size module and improved initial conditions. The Bayesian framework enables estimating a gauge/confidence interval over the validity of the prediction when applied to unknown observations. As a proof-of-concept, we apply DPNNet-Bayesian to dust gaps observed in HL Tau. The network predicts masses of $ 86.0 \pm 5.5 M_{\Earth} $, $ 43.8 \pm 3.3 M_{\Earth} $, and $ 92.2 \pm 5.1 M_{\Earth} $ respectively, which are comparable to other studies based on specialized simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Impact of Local Pressure Enhancements on Dust Concentration inTurbulent Protoplanetary Discs

We investigate the evolution of dust and gas in the vicinity of local pressure enhancements (pressure bumps) in a protoplanetary disc (PPD) with turbulence due to the Vertical Shear Instability (VSI). We perform global 2D axisymmetric and 3D simulations of dust and gas for a range of values for Z (ratio of dust-to-gas (d/g) surface mass densities or metallicity), particle Stokes numbers tau, and pressure bump amplitude A. Dust feedback onto the gas is included. For the first time we demonstrate in global 3D simulations the collection of dust in long-lived vortices induced by the VSI. Without a pressure bump and for Z~0.01 and tau~0.01 we find that such vortices reach d/g density ratios slightly below unity in the discs mid-plane, while for Z>0.05 long-lived vortices are largely absent. In presence of a pressure bump, for Z~0.01 & tau~0.01 a dusty vortex forms reaching d/g ratios of a few times unity, such that the SI is expected to develop, before it eventually shears out into a turbulent dust ring. For Z~0.03 this occurs for tau~0.005, with a weaker, more short-lived vortex, while for larger tau only a turbulent dust ring forms. For Z>0.03 we find that the dust ring becomes increasingly axisymmetric for increasing tau and d/g ratios reach ~1 for tau>0.005. Furthermore, the disc's vertical mass flow profile is strongly affected by dust for Z>0.03, such that gas is transported inward near the mid-plane and outward at larger heights, i.e. the reversed situation compared to simulations with zero or small amounts of dust. Viscous $α$-values drop moderately as 0.001-0.0001 for increasing Z=0-0.05. Our results suggest that the VSI can play an active role in planetesimal formation through the formation of vortices for plausible values of Z and tau. Also it may provide a natural explanation for the presence/absence of asymmetries of observed dust rings in PPDs, depending on the value of Z.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Streaming instabilities in accreting and magnetized laminar protoplanetary disks

The streaming instability is one of the most promising pathways to the formation of planetesimals from pebbles. Understanding how this instability operates under realistic conditions expected in protoplanetary disks is therefore crucial to assess the efficiency of planet formation. Contemporary models of protoplanetary disks show that magnetic fields are key to driving gas accretion through large-scale, laminar magnetic stresses. However, the effect of such magnetic fields on the streaming instability has not been examined in detail. To this end, we study the stability of dusty, magnetized gas in a protoplanetary disk. We find the streaming instability can be enhanced by passive magnetic torques and even persist in the absence of a global radial pressure gradient. In this case, instability is attributed to the azimuthal drift between dust and gas, unlike the classical streaming instability, which is driven by radial drift. This suggests that the streaming instability can remain effective inside dust-trapping pressure bumps in accreting disks. When a live vertical field is considered, we find the magneto-rotational instability can be damped by dust feedback, while the classic streaming instability can be stabilized by magnetic perturbations. We also find that Alfvén waves can be destabilized by dust-gas drift, but this instability requires nearly ideal conditions. We discuss the possible implications of these results for dust dynamics and planetesimal formation in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Machine Learning model to infer planet masses from gaps observed in protoplanetary disks

Observations of bright protoplanetary disks often show annular gaps in their dust emission. One interpretation of these gaps is disk-planet interaction. If so, fitting models of planetary gaps to observed protoplanetary disk gaps can reveal the presence of hidden planets. However, future surveys are expected to produce an ever-increasing number of protoplanetary disks with gaps. In this case, performing a customized fitting for each target becomes impractical owing to the complexity of disk-planet interaction. To this end, we introduce DPNNet (Disk Planet Neural Network), an efficient model of planetary gaps by exploiting the power of machine learning. We train a deep neural network with a large number of dusty disk-planet hydrodynamic simulations across a range of planet masses, disk temperatures, disk viscosities, disk surface density profiles, particle Stokes numbers, and dust abundances. The network can then be deployed to extract the planet mass for a given gap morphology. In this work, first in a series, we focus on the basic concepts of our machine learning framework. We demonstrate its utility by applying it to the dust gaps observed in the protoplanetary disk around HL Tau at $10$ au, $30$ au, and $80$ au. Our network predict planet masses of $80 \, M_{\rm Earth}$, $63 \, M_{\rm Earth}$, and $70 \, M_{\rm Earth}$, respectively, which are comparable to other studies based on specialized simulations. We discuss the key advantages of our DPNNet in its flexibility to incorporate new physics, any number of parameters and predictions, and its potential to ultimately replace hydrodynamical simulations for disk observers and modelers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How efficient is the streaming instability in viscous protoplanetary disks?

The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability with external turbulence modelled by gas viscosity and particle diffusion. We find the streaming instability is sensitive to turbulence, with growth rates becoming negligible for alpha-viscosity parameters $α\gtrsim \mathrm{St} ^{1.5}$, where $\mathrm{St}$ is the particle Stokes number. We explore the effect of non-linear drag laws, which may be applicable to porous dust particles, and find growth rates are modestly reduced. We also find that gas compressibility increase growth rates by reducing the effect of diffusion. We then apply linear theory to global models of viscous protoplanetary disks. For minimum-mass Solar nebula disk models, we find the streaming instability only grows within disk lifetimes beyond $\sim 10$s of AU, even for cm-sized particles and weak turbulence ($α\sim 10^{-4}$). Our results suggest it is rather difficult to trigger the streaming instability in non-laminar protoplanetary disks, especially for small particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Migrating Low-Mass Planets in Inviscid Dusty Protoplanetary Discs

Disc-driven planet migration is integral to the formation of planetary systems. In standard, gas-dominated protoplanetary discs, low-mass planets or planetary cores undergo rapid inwards migration and are lost to the central star. However, several recent studies indicate that the solid component in protoplanetary discs can have a significant dynamical effect on disc-planet interaction, especially when the solid-to-gas mass ratio approaches unity or larger and the dust-on-gas drag forces become significant. As there are several ways to raise the solid abundance in protoplanetary discs, for example through disc winds and dust-trapping in pressure bumps, it is important to understand how planets migrate through a dusty environment. To this end, we study planet migration in dust-rich discs via a systematic set of high-resolution, two-dimensional numerical simulations. We show that the inwards migration of low-mass planets can be slowed down by dusty dynamical corotation torques. We also identify a new regime of stochastic migration applicable to discs with dust-to-gas mass ratios $\gtrsim 0.3$ and particle Stokes numbers $\gtrsim 0.03$. In these cases, disc-planet interaction leads to the continuous development of small-scale, intense dust vortices that scatter the planet, which can potentially halt or even reverse the inwards planet migration. We briefly discuss the observational implications of our results and highlight directions for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stratified and vertically-shearing streaming instabilities in protoplanetary disks

Under the right conditions, the streaming instability between imperfectly coupled dust and gas is a powerful mechanism for planetesimal formation as it can concentrate dust grains to the point of gravitational collapse. In its simplest form, the streaming instability can be captured by analyzing the linear stability of unstratified disk models, which represent the midplane of protoplanetary disks. We extend such studies by carrying out vertically-global linear stability analyses of dust layers in protoplanetary disks. We find the dominant form of instability in stratified dust layers is one driven by the vertical gradient in the rotation velocity of the dust-gas mixture, but also requires partial dust-gas coupling. These vertically-shearing streaming instabilities grow on orbital timescales and occur on radial length scales $\sim10^{-3}H_\mathrm{g}$, where $H_\mathrm{g}$ is the local pressure scale height. The classic streaming instability, associated with the relative radial drift between dust and gas, occur on radial length scales $\sim10^{-2}H_\mathrm{g}$, but have much smaller growth rates than vertically-shearing streaming instabilities. Including gas viscosity is strongly stabilizing and leads to vertically-elongated disturbances. We briefly discuss the potential effects of vertically-shearing streaming instabilities on planetesimal formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the gravitational stability of gravito-turbulent accretion disks

Low mass, self-gravitating accretion disks admit quasi-steady,`gravito-turbulent' states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: 1) cooling, which reduces pressure support; and/or 2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond $\sim 60$AU dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cooling Requirements for the Vertical Shear Instability in Protoplanetary Disks

The vertical shear instability (VSI) offers a potential hydrodynamic mechanism for angular momentum transport in protoplanetary disks (PPDs). The VSI is driven by a weak vertical gradient in the disk's orbital motion, but must overcome vertical buoyancy, a strongly stabilizing influence in cold disks, where heating is dominated by external irradiation. Rapid radiative cooling reduces the effective buoyancy and allows the VSI to operate. We quantify the cooling timescale $t_c$ needed for efficient VSI growth, through a linear analysis of the VSI with cooling in vertically global, radially local disk models. We find the VSI is most vigorous for rapid cooling with $t_c<Ω_\mathrm{K}^{-1}h|q|/(γ-1)$ in terms of the Keplerian orbital frequency, $Ω_\mathrm{K}$; the disk's aspect-ratio, $h\ll1$; the radial power-law temperature gradient, $q$; and the adiabatic index, $γ$. For longer $t_c$, the VSI is much less effective because growth slows and shifts to smaller length scales, which are more prone to viscous or turbulent decay. We apply our results to PPD models where $t_c$ is determined by the opacity of dust grains. We find that the VSI is most effective at intermediate radii, from $\sim5$AU to $\sim50$AU with a characteristic growth time of $\sim30$ local orbital periods. Growth is suppressed by long cooling times both in the opaque inner disk and the optically thin outer disk. Reducing the dust opacity by a factor of 10 increases cooling times enough to quench the VSI at all disk radii. Thus the formation of solid protoplanets, a sink for dust grains, can impede the VSI.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gap formation and stability in non-isothermal protoplanetary discs

Several observations of transition discs show lopsided dust-distributions. A potential explanation is the formation of a large-scale vortex acting as a dust-trap at the edge of a gap opened by a giant planet. Numerical models of gap-edge vortices have thus far employed locally isothermal discs, but the theory of this vortex-forming or `Rossby wave' instability was originally developed for adiabatic discs. We generalise the study of planetary gap stability to non-isothermal discs using customised numerical simulations of disc-planet systems where the planet opens an unstable gap. We include in the energy equation a simple cooling function with cooling timescale $t_c=βΩ_k^{-1}$, where $Ω_k$ is the Keplerian frequency, and examine the effect of $β$ on the stability of gap edges and vortex lifetimes. We find increasing $β$ lowers the growth rate of non-axisymmetric perturbations, and the dominant azimuthal wavenumber $m$ decreases. We find a quasi-steady state consisting of one large-scale, over-dense vortex circulating the outer gap edge, typically lasting $O(10^3)$ orbits. Vortex lifetimes were found to generally increase with cooling times up to an optimal value, beyond which vortex lifetimes decrease. This non-monotonic dependence is qualitatively consistent with recent studies using strictly isothermal discs that vary the disc aspect ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One-armed spirals in locally isothermal, radially structured self-gravitating discs

We describe a new mechanism that leads to the destabilisation of non-axisymmetric waves in astrophysical discs with an imposed radial temperature gradient. This might apply, for example, to the outer parts of protoplanetary discs. We use linear density wave theory to show that non-axisymmetric perturbations generally do not conserve their angular momentum in the presence of a forced temperature gradient. This implies an exchange of angular momentum between linear perturbations and the background disc. In particular, when the disturbance is a low-frequency trailing wave and the disc temperature decreases outwards, this interaction is unstable and leads to the growth of the wave. We demonstrate this phenomenon through numerical hydrodynamic simulations of locally isothermal discs in 2D using the FARGO code and in 3D with the ZEUS-MP and PLUTO codes. We consider radially structured discs with a self-gravitating region which remains stable in the absence of a temperature gradient. However, when a temperature gradient is imposed we observe exponential growth of a one-armed spiral mode (azimuthal wavenumber $m=1$) with co-rotation radius outside the bulk of the spiral arm, resulting in a nearly-stationary one-armed spiral pattern. The development of this one-armed spiral does not require the movement of the central star, as found in previous studies. Because destabilisation by a forced temperature gradient does not explicitly require disc self-gravity, we suggest this mechanism may also affect low-frequency one-armed oscillations in non-self-gravitating discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational instability of planetary gaps and its effect on orbital migration

Gap formation by giant planets in self-gravitating disks may lead to a gravitational edge instability (GEI). We demonstrate this GEI with global 3D and 2D self-gravitating disk-planet simulations using the ZEUS, PLUTO and FARGO hydrodynamic codes. High resolution 2D simulations show that an unstable outer gap edge can lead to outwards orbital migration. Our results have important implications for theories of giant planet formation in massive disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear stability of magnetized massive protoplanetary disks

Magneto-rotational instability (MRI) and gravitational instability (GI) are the two principle routes to turbulent angular momentum transport in accretion disks. Protoplanetary disks may develop both. This paper aims to reinvigorate interest in the study of magnetized massive protoplanetary disks, starting from the basic issue of stability. The local linear stability of a self-gravitating, uniformly magnetized, differentially rotating, three-dimensional stratified disk subject to axisymmetric perturbations is calculated numerically. The formulation includes resistivity. It is found that the reduction in the disk thickness by self-gravity can decrease MRI growth rates; the MRI becomes global in the vertical direction, and MRI modes with small radial length scales are stabilized. The maximum vertical field strength that permits the MRI in a strongly self-gravitating polytropic disk with polytropic index $Γ=1$ is estimated to be $B_{z,\mathrm{max}} \simeq c_{s0}Ω\sqrt{μ_0/16πG} $, where $c_{s0}$ is the midplane sound speed and $Ω$ is the angular velocity. In massive disks with layered resistivity, the MRI is not well-localized to regions where the Elsasser number exceeds unity. For MRI modes with radial length scales on the order of the disk thickness, self-gravity can enhance density perturbations, an effect that becomes significant in the presence of a strong toroidal field, and which depends on the symmetry of the underlying MRI mode. In gravitationally unstable disks where GI and MRI growth rates are comparable, the character of unstable modes can transition smoothly between MRI and GI. Implications for non-linear simulations are discussed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Instabilities at planetary gap edges in 3D self-gravitating disks

Numerical simulations are presented to study the stability of gaps opened by giant planets in 3D self-gravitating disks. In weakly self-gravitating disks, a few vortices develop at the gap edge and merge on orbital time-scales. The result is one large but weak vortex with Rossby number -0.01. In moderately self-gravitating disks, more vortices develop and their merging is resisted on dynamical time-scales. Self-gravity can sustain multi-vortex configurations, with Rossby number -0.2 to -0.1, over a time-scale of order 100 orbits. Self-gravity also enhances the vortex vertical density stratification, even in disks with initial Toomre parameter of order 10. However, vortex formation is suppressed in strongly self-gravitating disks and replaced by a global spiral instability associated with the gap edge which develops during gap formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-barotropic linear Rossby wave instability in three-dimensional disks

Astrophysical disks with localized radial structure, such as protoplanetary disks containing dead zones or gaps due to disk-planet interaction, may be subject to the non-axisymmetric Rossby wave instability (RWI) that lead to vortex-formation. The linear instability has recently been demonstrated in three-dimensional (3D) barotropic disks. It is the purpose of this study to generalize the 3D linear problem to include an energy equation, thereby accounting for baroclinity in three-dimensions. Linear stability calculations are presented for radially structured, vertically stratified, geometrically-thin disks with non-uniform entropy distribution in both directions. Polytropic equilibria are considered but adiabatic perturbations assumed. The unperturbed disk has a localized radial density bump making it susceptible to the RWI. The linearized fluid equations are solved numerically as a partial differential equation eigenvalue problem. Emphasis on the ease of method implementation is given. It is found that when the polytropic index is fixed and adiabatic index increased, non-uniform entropy has negligible effect on the RWI growth rate, but pressure and density perturbation magnitudes near a pressure enhancement increases away from the midplane. The associated meridional flow is also qualitatively changed from homentropic calculations. Meridional vortical motion is identified in the nonhomentropic linear solution, as well as in a nonlinear global hydrodynamic simulation of the RWI in an initially isothermal disk evolved adiabatically. Numerical results suggest buoyancy forces play an important role in the internal flow of Rossby vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orbital migration of giant planets induced by gravitationally unstable gaps: the effect of planet mass

It has been suggested that long-period giant planets, such as HD 95086b and HR 8799bcde, may have formed through gravitational instability of protoplanetary discs. However, self-gravitating disc-satellite interaction can lead to the formation of a gravitationally unstable gap. Such an instability significantly affects the orbital migration of gap-opening perturbers in massive discs. We use 2D hydrodynamical simulations to examine the role of planet mass on the gravitational stability of gaps and its impact on orbital migration. We consider giant planets with planet-to-star mass ratio q=0.0003 to q=0.003, in a self-gravitating disc with disc-to-star mass ratio M_d/M_*=0.08, aspect ratio h=0.05, and Keplerian Toomre parameter Q = 1.5 at 2.5 times the planet's initial orbital radius. Fixed-orbit simulations show that all planet masses we consider open gravitationally unstable gaps, but the instability is stronger and develops sooner with increasing planet mass. The disc-on-planet torques typically become more positive with increasing planet mass. In freely-migrating simulations, we observe faster outward migration with increasing planet mass, but only for planet masses capable of opening unstable gaps early on. For q=0.0003, the planet undergoes rapid inward type III migration before it can open a gap. For q=0.0013 we find it is possible to balance the tendency for inward migration by the positive torques due to an unstable gap, but only for a few 10's of orbital periods. We find the unstable outer gap edge can trigger outward type III migration, sending the planet to twice its initial orbital radius on dynamical timescales. This triggered type III migration may be a way of bringing giant planets to large orbital distances. On the other hand, our results also imply gap-opening is not a sufficient condition for the in situ formation of giant planets on wide orbits through disc fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Steady state of dust distributions in disk vortices: Observational predictions and applications to transitional disks

The Atacama Large Millimeter Array (ALMA) has been returning images of transitional disks in which large asymmetries are seen in the distribution of mm-sized dust in the outer disk. The explanation in vogue borrows from the vortex literature by suggesting that these asymmetries are the result of dust trapping in giant vortices, excited via Rossby wave instability (RWI) at planetary gap edges. Due to the drag force, dust trapped in vortices will accumulate in the center, and diffusion is needed to maintain a steady state over the lifetime of the disk. While previous work derived semi-analytical models of the process, in this paper we provide analytical steady-state solutions. Exact solutions exist for certain vortex models. The solution is determined by the vortex rotation profile, the gas scale height, the vortex aspect ratio, and the ratio of dust diffusion to gas-dust friction. In principle, all these quantities can be derived from observations, which would give validation of the model, also giving constrains on the strength of the turbulence inside the vortex core. Based on our solution, we derive quantities such as the gas-dust contrast, the trapped dust mass, and the dust contrast at the same orbital location. We apply our model to the recently imaged Oph IRS 48 system, finding values within the range of the observational uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing large-scale vortex formation against viscous layers in three-dimensional discs

Vortex formation through the Rossby wave instability (RWI) in protoplanetary discs has been invoked to play a role in planet formation theory, and suggested to explain the observation of large dust asymmetries in several transitional discs. However, whether or not vortex formation operates in layered accretion discs, i.e. models of protoplanetary discs including dead zones near the disc midplane --- regions that are magnetically inactive and the effective viscosity greatly reduced --- has not been verified. As a first step toward testing the robustness of vortex formation in layered discs, we present non-linear hydrodynamical simulations of global 3D protoplanetary discs with an imposed kinematic viscosity that increases away from the disc midplane. Two sets of numerical experiments are performed: (i) non-axisymmetric instability of artificial radial density bumps in viscous discs; (ii) vortex-formation at planetary gap edges in layered discs. Experiment (i) shows that the linear instability is largely unaffected by viscosity and remains dynamical. The disc-planet simulations also show the initial development of vortices at gap edges, but in layered discs the vortices are transient structures which disappear well into the non-linear regime. We suggest that the long term survival of columnar vortices, such as those formed via the RWI, requires low effective viscosity throughout the vertical extent of the disc, so such vortices do not persist in layered discs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamical instabilities in disc-planet interactions

Protoplanetary discs may become dynamically unstable due to structure induced by an embedded giant planet. In this thesis, I discuss the stability of such systems and explore the consequence of instability on planetary migration. I begin with non-self-gravitating, low viscosity discs and show that giant planets induce shocks inside its co-orbital region, leading to a profile unstable to vortex formation around a potential vorticity minimum. This instability is commonly known as the vortex or Rossby wave instability. Vortex-planet interaction lead to episodic phases of migration, which can be understood in the framework of type III migration. I then examine the effect of disc self-gravity on gap stability. The linear theory of the Rossby wave instability is extended to include disc gravity, which shows that self-gravity is effective at stabilising the vortex instability at small azimuthal wavenumber. This is consistent with the observation that more vortices develop with increasing disc mass in hydrodynamic simulations. Vortices in self-gravitating discs also resist merging, and is most simply understood as pair-vortices undergoing mutual horsehoe turns upon encounter. I show that in sufficiently massive discs vortex modes are suppressed. Instead, global spiral instabilities develop which are associated with a potential vorticity maximum at the gap edge. These edge modes can be physically understood as a result of unstable interaction between the gap edge and the exterior disc through gravity. I show the spiral arms can provide a positive torque on the planet, leading to fast migration outwards. I confirm the above results, obtained from razor-thin disc models, persist in three-dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of upper disc boundary conditions on the linear Rossby wave instability

The linear Rossby wave instability (RWI) in global, 3D polytropic discs is revisited with a much simpler numerical method than that previously employed by the author. The governing partial differential equation is solved with finite differences in the radial direction and spectral collocation in the vertical direction. RWI modes are calculated subject to different upper disc boundary conditions. These include free surface, solid boundaries and variable vertical domain size. Boundary conditions that oppose vertical motion increase the instability growth rate by a few per cent. The magnitude of vertical flow throughout the fluid column can be affected but the overall flow pattern is qualitatively unchanged. Numerical results support the notion that the RWI is intrinsically two dimensional. This implies that inconsistent upper disc boundary conditions, such as vanishing enthalpy perturbation, may inhibit the RWI in 3D.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rossby wave instability in locally isothermal and polytropic disks: three-dimensional linear calculations

Numerical calculations of the linear Rossby wave instability (RWI) in global three-dimensional (3D) disks are presented. The linearized fluid equations are solved for vertically stratified, radially structured disks with either a locally isothermal or polytropic equation of state, by decomposing the vertical dependence of the perturbed hydrodynamic quantities into Hermite and Gegenbauer polynomials, respectively. It is confirmed that the RWI operates in 3D. For perturbations with vertical dependence assumed above, there is little difference in growth rates between 3D and two-dimensional (2D) calculations. Comparison between 2D and 3D solutions of this type suggest the RWI is predominantly a 2D instability and that three-dimensional effects, such as vertical motion, to be interpreted as a perturbative consequence of the dominant 2D flow. The vertical flow around co-rotation, where vortex-formation is expected, is examined. In locally isothermal disks the expected vortex center remains in approximate vertical hydrostatic equilibrium. For polytropic disks the vortex center has positive vertical velocity, whose magnitude increases with decreasing polytropic index $n$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vortex and spiral instabilities at gap edges in three-dimensional self-gravitating disc-satellite simulations

Numerical simulations of global three-dimensional (3D), self-gravitating discs with a gap opened by an embedded planet are presented. The simulations are customised to examine planetary gap stability. Previous results, obtained by Lin & Papaloizou from two-dimensional (2D) disc models, are reproduced in 3D. These include (i) the development of vortices associated with local vortensity minima at gap edges and their merging on dynamical timescales in weakly self-gravitating discs, (ii) the increased number of vortices as the strength of self-gravity is increased and their resisted merging, and (iii) suppression of the vortex instability and development of global spiral arms associated with local vortensity maxima in massive discs. The vertical structure of these disturbances are examined. In terms of the relative density perturbation, the vortex disturbance has weak vertical dependence when self-gravity is neglected. Vortices become more vertically stratified with increasing self-gravity. This effect is seen even when the unperturbed region around the planet's orbital radius has a Toomre stability parameter ~10. The spiral modes display significant vertical structure at the gap edge, with the midplane density enhancement being several times larger than that near the upper disc boundary. However, for both instabilities the vertical Mach number is typically a few per cent,and on average vertical motions near the gap edge do not dominate horizontal motions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Edge modes in self-gravitating disc-planet interactions

We study the stability of gaps opened by a giant planet in a self-gravitating protoplanetary disc. We find a linear instability associated with both the self-gravity of the disc and local vortensity maxima which coincide with gap edges. For our models, these edge modes develop and extend to twice the orbital radius of a Saturn mass planet in discs with disc-to-star mass ratio >0.06, corresponding to a Toomre Q < 1.5 at the outer disc boundary. Unlike the local vortex-forming instabilities associated with gap edges in weakly or non-self-gravitating low viscosity discs, the edge modes are global and exist only in sufficiently massive discs, but for the typical viscosity values adopted for protoplanetary discs. Analytic modelling and linear calculations show edge modes may be interpreted as a localised disturbance associated with a gap edge inducing activity in the extended disc, through the launching of density waves excited at Lindblad resonances. Nonlinear hydrodynamic simulations are performed to investigate the evolution of edge modes in disc-planet systems. The form and growth rates of unstable modes are consistent with linear theory. Their dependence on viscosity and gravitational softening is also explored. We also performed a first study of the effect of edge modes on planetary migration. We found that if edge modes develop, then the average disc-on-planet torque becomes more positive with increasing disc mass. In simulations where the planet was allowed to migrate, although a fast type III migration could be seen that was similar to that seen in non-self-gravitating discs, we found that it was possible for the planet to interact gravitationally with the spiral arms associated with an edge mode and that this could result in the planet being scattered outwards. Thus orbital migration is likely to be complex and non monotonic in massive discs of the type we consider.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Outward migration of a giant planet with a gravitationally unstable gap edge

We present numerical simulations of disc-planet interactions where the planet opens a gravitationally unstable gap in an otherwise gravitationally stable disc. In our disc models, where the outer gap edge can be unstable to global spiral modes, we find that as we increase the surface density scale the gap becomes more unstable and the planet migrates outwards more rapidly. We show that the positive torque is provided by material brought into the planet's coorbital region by the spiral arms. This material is expected to execute horseshoe turns upon approaching the planet and hence torque it. Our results suggest that standard type II migration, applicable to giant planets in non-self-gravitating viscous discs, is likely to be significantly modified in massive discs when gravitational instabilities associated with the gap occur.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spin down of protostars through gravitational torques

Young protostars embedded in circumstellar discs accrete from an angular momentum-rich mass reservoir. Without some braking mechanism, all stars should be spinning at or near break-up velocity. In this paper, we perform simulations of the self-gravitational collapse of an isothermal cloud using the ORION adaptive mesh refinement code and investigate the role that gravitational torques might play in the spin-down of the dense central object. While magnetic effects likely dominate for low mass stars, high mass and Population III stars might be less well magnetised. We find that gravitational torques alone prevent the central object from spinning up to more than half of its breakup velocity, because higher rotation rates lead to bar-like deformations that enable efficient angular momentum transfer to the surrounding medium. We also find that the long-term spin evolution of the central object is dictated by the properties of the surrounding disc. In particular, spiral modes with azimuthal wavenumber $m=2$ couple more effectively to its spin than the lopsided $m=1$ mode, which was found to inhibit spin evolution. We suggest that even in the absence of magnetic fields, gravitational torques may provide an upper limit on stellar spin, and that moderately massive circumstellar discs can cause long-term spin down.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effect of self-gravity on vortex instabilities in disc-planet interactions

We study the effect of disc self-gravity on vortex-forming instabilities associated with gaps opened by a Saturn mass planet in a protoplanetary disc. It is shown analytically and confirmed through linear calculations that vortex modes with low azimuthal mode number $,m,$ are stabilised by increasing self-gravity if the basic state is fixed. However, linear calculations show that the combined effect of self-gravity through the background and through the linear response shifts the most unstable vortex mode to higher $m.$ Nonlinear hydrodynamic simulations of planetary gaps show more vortices develop with increasing strength of self-gravity. For sufficiently large disc mass the vortex instability is suppressed and replaced by a new global instability, consistent with analytical expectations. In the nonlinear regime, vortex merging is increasingly delayed as the disc mass increases and multiple vortices may persist until the end of simulations. With self-gravity, the post-merger vortex is localised in azimuth and has similar structure to a Kida-like vortex. This is unlike the case without self-gravity where vortices merge to form a single vortex extended in azimuth. We also performed a series of supplementary simulations of co-orbital Kida-like vortices and found that self-gravity enables such vortices to execute horseshoe turns upon encountering each other. As a result vortex merging is avoided on time-scales where it would occur without self-gravity. Thus we suggest that mutual repulsion of self-gravitating vortices in a rotating flow is responsible for the delayed vortex merging above. The effect of self-gravity on vortex-induced migration is briefly discussed. We found that when self-gravity is included, the vortex-induced type III migration of Lin & Papaloizou (2010) is delayed but the extent of migration is unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Type III migration in a low viscosity disc

We study the type III migration of a Saturn mass planet in low viscosity discs. The planet is found to experience cyclic episodes of rapid decay in orbital radius, each amounting to a few Hill radii. We find this to be due to the scattering of large- scale vortices present in the disc. The origin and role of vortices in the context of type III migration is explored. It is shown through numerical simulations and semi- analytical modelling that spiral shocks induced by a sufficiently massive planet will extend close to the planet orbital radius. The production of vortensity across shock tips results in thin high vortensity rings with a characteristic width of the local scale height. For planets with masses equal to and above that of Saturn, the rings are co-orbital features extending the entire azimuth. Linear stability analysis show there exists unstable modes that are localised about local vortensity minima which coincide with gap edges. Simulations show that vortices are non-linear a outcome. We used hydrodynamic simulations to examine vortex-planet interactions. Their effect is present in discs with kinematic viscosity less than about an order of magnitude smaller than the typically adopted value of ν= 10^{-5}Ω_pr_p(0)^2, where r_p(0) and Ω_p are the initial orbital radius and angular velocity of the planet respectively. We find that the magnitude of viscosity affects the nature of type III migration but not the extent of the orbital decay. The role of vortices as a function of initial disc mass is also explored and it is found that the amount of orbital decay during one episode of vortex-planet interaction is independent of initial disc mass. We incorporate the concept of the co-orbital mass deficit in the analysis of our results and link it to the presence of vortices at gap edges.

0 citations0 reviewsNo code linkedOpen access