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

Wladimir Lyra

Wladimir Lyra 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.IMastro-ph.GAastro-ph.HEphysics.ao-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

23works
8topics
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

23 published item(s)

preprint2022arXiv

Exploring Jupiter's Polar Deformation Lengths with High Resolution Shallow Water Modeling

The polar regions of Jupiter host a myriad of dynamically interesting phenomena including vortex configurations, folded-filamentary regions (FFRs), and chaotic flows. Juno observations have provided unprecedented views of the high latitudes, allowing for more constraints to be placed upon the troposphere and the overall atmospheric energy cycle. Moist convective events are believed to be the primary drivers of energetic storm behavior as observed on the planet. Here, we introduce a novel single layer shallow water model to investigate the effects of polar moist convective events at high resolution, the presence of dynamical instabilities over long timescales, and the emergence of FFRs at high latitudes. We use a flexible, highly parallelizable, finite-difference hydrodynamic code to explore the parameter space set up by previous models. We study the long term effects of deformation length (Ld), injected pulse size, and injected geopotential. We find that models with Ld beyond 1500 km (planetary Burger number, Bu$=4.4\times10^{-4}$) tend to homogenize their potential vorticity (PV) in the form of dominant stable polar cyclones, while lower Ld cases tend to show less stability with regards to Arnol'd-type flows. We also find that large turbulent forcing scales consistently lead to the formation of high latitude FFRs. Our findings support the idea that moist convection, occurring at high latitudes, may be sufficient to produce the dynamical variety seen at the Jovian poles. Additionally, derived values of localized horizontal shear and Ld may constrain FFR formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae

Direct images of protoplanets embedded in disks around infant stars provide the key to understanding the formation of gas giant planets like Jupiter. Using the Subaru Telescope and Hubble Space Telescope, we find evidence for a jovian protoplanet around AB Aurigae orbiting at a wide projected separation (93 au), likely responsible for multiple planet-induced features in the disk. Its emission is reproducible as reprocessed radiation from an embedded protoplanet. We also identify two structures located at 430-580 au that are candidate sites of planet formation. These data reveal planet formation in the embedded phase and a protoplanet discovery at wide, > 50 au separations characteristic of most imaged exoplanets. With at least one clump-like protoplanet and multiple spiral arms, the AB Aur system may also provide the evidence for a long-considered alternative to the canonical model for Jupiter's formation: disk (gravitational) instability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Linking atmospheric chemistry of the hot Jupiter HD 209458b to its formation location through infrared transmission and emission spectra

The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D chemical kinetics model, treating both vertical mixing and photochemistry. We consider two atmospheric pressure-temperature profiles, one with and one without a thermal inversion. From each of the resulting 32 atmospheric composition profiles, we find that the molecules CH4, NH3, HCN, and C2H2 are more prominent in the atmospheres computed using a realistic non-inverted P-T profile in comparison to a prior equilibrium chemistry based work which used an analytical P-T profile. We also compute the synthetic transmission and emission spectra for these atmospheres and find that many spectral features vary with the location in the disc where the planet's envelope was accreted. By comparing with the species detected using the latest high-resolution ground-based observations, our model suggests HD 209458b could have accreted most of its gas between the CO2 and CH4 icelines with a super solar C/O ratio from its protostellar disc, which in turn directly inherited its chemical abundances from the protostellar cloud. Finally, we simulate observing the planet with the James Webb Space Telescope (JWST) and show that differences in spectral signatures of key species can be recognized. Our study demonstrates the enormous importance of JWST in providing new insights into hot Jupiter's formation environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gas and dust dynamics in starlight-heated protoplanetary disks

Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$μ$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $α\simeq 10^{-4}$. The value of $α$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $α$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Sublimative Evolution of (486958) Arrokoth

We consider the history of New Horizons target (486958) Arrokoth in the context of its sublimative evolution. Shortly after the Sun's protoplanetary disk (PPD) cleared, the newly intense sunlight sparked a sublimative period in Arrokoth's early history that lasted for ~10-100 Myr. Although this sublimation was too weak to significantly alter Arrokoth's spin state, it could drive mass transport around the surface significant enough to erase topographic features on length scales of ~10-100 m. This includes craters up to ~50-500 m in diameter, which suggests that the majority of Arrokoth's craters may not be primordial (dating from the merger of Arrokoth's lobes), but rather could date from after the end of this sublimative period. Thereafter, Arrokoth entered a Quiescent Period (which lasts to the present day), in which volatile production rates are at least 13 orders of magnitude less than the ~10^24 molecules/s detection limit of the New Horizons spacecraft (Lisse et al. 2020). This is insufficient to drive either mass transport or sublimative torques. These results suggest that the observed surface of Arrokoth is not primordial, but rather dates from the Quiescent Period. By contrast, the inability of sublimative torques to meaningfully alter Arrokoth's rotation state suggests that its shape is indeed primordial, and its observed rotation is representative of its spin state after formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Disentangling Planets from Photoelectric Instability in Gas-Rich Optically Thin Dusty Disks

Structures in circumstellar disks such as gaps and rings are often attributed to planets. This connection has been difficult to show unequivocally, as other processes may also produce these features. Particularly, a photoelectric instability (PEI) has been proposed, operating in gas-rich optically thin disks, that generates structures predicted by planet--disk interactions. We examine the question of how to disentangle planetary effects on disk structure from the effects of the PEI. We use the Pencil Code to perform 2D global hydrodynamical models of the dynamics of gas and dust in a thin disk, with and without planetary perturbers. Photoelectric heating is modeled with an equation of state where pressure is proportional to dust surface density. The drag force on grains and its backreaction on the gas are included. Analyzing the situation without PEI, we find that gas--dust interactions alter the shape of the planetary gap from the dust-free case when the local dust-to-gas ratio $\varepsilon$ approaches unity. This result applies also to primordial disks, because dust drifting inwards accumulates at the edge of the planetary gap, and any initial dust-to-gas ratio eventually achieves $\varepsilon=1$ if the dust reservoir is sufficient. We find a result particular to high dust-to-gas ratio disks as well: as dust drifts inwards, the dust front becomes a sharp transition, and the backreaction triggers the Rossby wave instability. When PEI is included, we find that it obscures structures induced by planets unless the planet's mass is sufficiently large to carve a noticeable gap. Specifically, the instability generates arcs and rings of regular spacing: a planet is discernible when it carves a dust gap wider than the wavelength of the PEI.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

CSI 2264: Characterizing Young Stars in NGC 2264 with Short-Duration, Periodic Flux Dips in their Light Curves

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On shocks driven by high-mass planets in radiatively inefficient disks. I. Two-dimensional global disk simulations

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet-disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1-10 M_jup for a 1 M_sun star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M>5 M_jup) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M_jup planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet-disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On shocks driven by high-mass planets in radiatively inefficient disks. II. Three-dimensional global disk simulations

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $α$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle Trapping and Streaming Instability in Vortices

We analyse the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, i.e. baroclinic vortex growth. Besides the gas drag acting on particles we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our 2D, local, shearing sheet simulations we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, ${\rm St}$. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for ${\rm St}=1$. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity-assisted planetesimal formation even for small particles ($\rm{St}=0.01$) and low initial dust-to-gas ratios (1:$10^4$)

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convective overstability in accretion disks: 3D linear analysis and nonlinear saturation

Recently, Klahr & Hubbard (2014) claimed that a hydrodynamical linear overstability exists in protoplanetary disks, powered by buoyancy in the presence of thermal relaxation. We analyse this claim, confirming it through rigorous compressible linear analysis. We model the system numerically, reproducing the linear growth rate for all cases studied. We also study the saturated properties of the overstability in the shearing box, finding that the saturated state produces finite amplitude fluctuations strong enough to trigger the subcritical baroclinic instability. Saturation leads to a fast burst of enstrophy in the box, and a large-scale vortex develops in the course of the next $\approx$100 orbits. The amount of angular momentum transport achieved is of the order of $α\approx 10^{-3}$, as in compressible SBI models. For the first time, a self-sustained 3D vortex is produced from linear amplitude perturbation of a quiescent base state.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radial Stresses and Energy Transport in Accretion Disks

Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sharp eccentric rings in planetless hydrodynamical models of debris disks

Debris disks should not be completely gas-free, since there is second generation gas from outgassing of planetesimals and dust grains via sublimation, photodesorption, or collisions, generating a system of dust-to-gas ratio close to unity, where hydrodynamics cannot be ignored. A clumping instability exists in this configuration, that has been hitherto explored only in one-dimensional, incompressible models. We performed 2D numerical compressible models of a disk with comparable amounts of gas and dust to study the growth and development of this instability. Our model solves the momentum equation for the gas and dust, together with energy and continuity equations. We uncover that the backreaction of the drag force from the gas onto the dust shepherds rings, similar to those observed in debris disks and usually attributed to the presence of hypothetical undetected planets. We also uncover that the eccentricity of these rings, usually presented as convincing evidence for the presence of a planet, can actually be simply explained by a standing wave propagating along the ring. The rings support a spectrum of oscillations, with one particular mode representing epicyclic motion. The apparent eccentricity matches the eccentricity in observed systems. This suggests that the planet possibility, though thrilling, is not necessarily required to explain these systems.

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
preprint2012arXiv

A Well-Posed Kelvin-Helmholtz Instability Test and Comparison

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two dimensional Kelvin-Helmholtz instability. We ran the problem in the Pencil Code, Athena, Enzo, NDSPHMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of Smoothed Particle Hydrodynamics. We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in Smoothed Particle Hydrodynamics interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the connection between the magneto-elliptic and magneto-rotational instabilities

It has been recently suggested that the magneto-rotational instability (MRI) is a limiting case of the magneto-elliptic instability (MEI). This limit is obtained for horizontal modes in the presence of rotation and an external vertical magnetic field, when the aspect ratio of the elliptic streamlines tends to infinite. In this paper we unveil the link between these previously unconnected mechanisms, explaining both the MEI and the MRI as different manifestations of the same Magneto-Elliptic-Rotational Instability (MERI). The growth rates are found and the influence of the magnetic and rotational effects is explained, in particular the effect of the magnetic field on the range of negative Rossby numbers at which the horizontal instability is excited. Furthermore, we show how the horizontal rotational MEI in the rotating shear flow limit links to the MRI by the use of the local shearing box model, typically used in the study of accretion discs. In such limit the growth rates of the two instability types coincide for any power-type background angular velocity radial profile with negative exponent corresponding to the value of the Rossby number of the rotating shear flow. The MRI requirement for instability is that the background angular velocity profile is a decreasing function of the distance from the centre of the disk which corresponds to the horizontal rotational MEI requirement of negative Rossby numbers. Finally a physical interpretation of the horizontal instability, based on a balance between the strain, the Lorentz force and the Coriolis force is given.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Orbital migration of interacting low-mass planets in evolutionary radiative turbulent models

The torques exerted by a locally isothermal disk on an embedded planet lead to rapid inward migration. Recent work has shown that modeling the thermodynamics without the assumption of local isothermality reveals regions where the net torque on an embedded planet is positive, leading to outward migration of the planet. When a region with negative torque lies directly exterior to this, planets in the inner region migrate outwards and planets in the outer region migrate inwards, converging where the torque is zero. We incorporate the torques from an evolving non-isothermal disk into an N-body simulation to examine the behavior of planets or planetary embryos interacting in the convergence zone. We find that mutual interactions do not eject objects from the convergence zone. Small numbers of objects in a laminar disk settle into near resonant orbits that remain stable over the 10 Myr periods that we examine. However, either or both increasing the number of planets or including a correlated, stochastic force to represent turbulence drives orbit crossings and mergers in the convergence zone. These processes can build gas giant cores with masses of order ten Earth masses from sub-Earth mass embryos in 2-3 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rossby wave instability at dead zone boundaries in 3D resistive magnetohydrodynamical global models of protoplanetary disks

It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parametrized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in 3D resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the Pencil Code to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the $α\approx 10^{-2}$ level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped, or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Formation of planetary cores at Type I migration traps

One of the longstanding unsolved problems of planet formation is how solid bodies of a few decimeters in size can "stick" to form large planetesimals. This is known as the "meter size barrier". In recent years it has become increasingly clear that some form of "particle trapping" must have played a role in overcoming the meter size barrier. Particles can be trapped in long-lived local pressure maxima, such as those in anticyclonic vortices, zonal flows or those believed to occur near ice lines or at dead zone boundaries. Such pressure traps are the ideal sites for the formation of planetesimals and small planetary embryos. Moreover, they likely produce large quantities of such bodies in a small region, making it likely that subsequent N-body evolution may lead to even larger planetary embryos. The goal of this Letter is to show that this indeed happens, and to study how efficient it is. In particular, we wish to find out if rocky/icy bodies as large as 10 M_earth can form within 1 Myr, since such bodies are the precursors of gas giant planets in the core accretion scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Meridional circulation in turbulent protoplanetary disks

Based on the viscous disk theory, a number of recent studies have suggested there is large scale meridional circulation in protoplanetary disks. Such a flow could account for the presence of crystalline silicates, including calcium- and aluminum-rich inclusions (CAIs), at large distances from the sun. This paper aims at examining whether such large-scale flows exist in turbulent protoplanetary disks. High-resolution global hydrodynamical and magnetohydrodynamical (MHD) numerical simulations of turbulent protoplanetary disks were used to infer the properties of the flow in such disks. By performing hydrodynamic simulations using explicit viscosity, we demonstrate that our numerical setup does not suffer from any numerical artifact. The aforementioned meridional circulation is easily recovered in viscous and laminar disks and is quickly established. In MHD simulations, the magnetorotational instability drives turbulence in the disks. Averaging out the turbulent fluctuations on a long timescale, the results fail to show any large-scale meridional circulation. A detailed analysis of the simulations show that this lack of meridional circulation is due to the turbulent stress tensor having a vertical profile different from the viscous stress tensor. A simple model is provided that successfully accounts for the structure of the flow in the bulk of the disk. In addition to those results, possible deviations from standard vertically averaged alpha disk models are suggested by the simulations and should be the focus of future work. Global MHD numerical simulations of fully ionized and turbulent protoplanetary disks are not consistent with the existence of a large-scale meridional flow. As a consequence, the presence of crystalline silicates at large distance for the central star cannot be accounted for by that process as suggested by recent models based on viscous disk theory.

0 citations0 reviewsNo code linkedOpen access