Showing 1–50 of 52 results for author: Umurhan, O M

Rapid protoplanet formation in vortices: three-dimensional local simulations with selfgravity

Authors: Wladimir Lyra, Chao-Chin Yang, Jacob B. Simon, Orkan M. Umurhan, Andrew N. Youdin

Abstract: Disk vortices, seen in numerical simulations of protoplanetary disks and found observationally in ALMA and VLA images of these objects, are promising sites for planet formation given their pebble trapping abilities. Previous works have shown strong concentration of pebbles in vortices, but gravitational collapse has only been shown in low-resolution, two-dimensional, global models. In this letter,… ▽ More Disk vortices, seen in numerical simulations of protoplanetary disks and found observationally in ALMA and VLA images of these objects, are promising sites for planet formation given their pebble trapping abilities. Previous works have shown strong concentration of pebbles in vortices, but gravitational collapse has only been shown in low-resolution, two-dimensional, global models. In this letter, we aim to study the pebble concentration and gravitational collapse of pebble clouds in vortices via high-resolution, three-dimensional, local models. We performed simulations of the dynamics of gas and solids in a local shearing box where the gas is subject to convective overstability, generating a persistent giant vortex. We find that the vortex produces objects of Moon and Mars mass, with mass function of power law $d\ln N/d\ln M=-1.6\pm 0.3$. The protoplanets grow rapidly, doubling in mass in about 5 orbits, following pebble accretion rates. The mass range and mass doubling rate are in broad agreement with previous low resolution global models. We conclude that Mars-mass planetary embryos are the natural outcome of planet formation inside the disk vortices seen in millimeter and radio images of protoplanetary disks. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 9 pages, 5 figures, accepted for publication in ApJ letters

Length and Velocity Scales in Protoplanetary Disk Turbulence

Authors: Debanjan Sengupta, Jeffrey N. Cuzzi, Orkan M. Umurhan, Wladimir Lyra

Abstract: In the theory of protoplanetary disk turbulence, a widely adopted \emph{ansatz}, or assumption, is that the turnover frequency of the largest turbulent eddy, $Ω_L$, is the local Keplerian frequency $Ω_K$. In terms of the standard dimensionless Shakura-Sunyaev $α$ parameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given r… ▽ More In the theory of protoplanetary disk turbulence, a widely adopted \emph{ansatz}, or assumption, is that the turnover frequency of the largest turbulent eddy, $Ω_L$, is the local Keplerian frequency $Ω_K$. In terms of the standard dimensionless Shakura-Sunyaev $α$ parameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrtαH$ and $\sqrtαc$, in which $H$ and $c$ are the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as Vertical Shear Instability. Here we explore the more general case where $Ω_L\geΩ_K$ and show that under these conditions, the characteristic length and velocity scales are respectively $\sqrt{α/R'}H$ and $\sqrt{αR'}c$, where $R'\equiv Ω_L/Ω_K$ is twice the Rossby number. It follows that $α=\alphat/R'$, where $\sqrt{\alphat} c$ is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R'>1$ the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

A solution for the density dichotomy problem of Kuiper Belt objects with multi-species streaming instability and pebble accretion

Authors: Manuel H. Cañas, Wladimir Lyra, Daniel Carrera, Leonardo Krapp, Debanjan Sengupta, Jacob B. Simon, Orkan M. Umurhan, Chao-Chin Yang, Andrew Youdin

Abstract: Kuiper belt objects show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of $^{26}$Al.… ▽ More Kuiper belt objects show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of $^{26}$Al. We aim to explain the density trend in the context of streaming instability and pebble accretion. Small pebbles experience lofting into the atmosphere of the disk, being exposed to UV and partially losing their ice via desorption. Conversely, larger pebbles are shielded and remain more icy. We use a shearing box model including gas and solids, the latter split into ices and silicate pebbles. Self-gravity is included, allowing dense clumps to collapse into planetesimals. We find that the streaming instability leads to the formation of mostly icy planetesimals, albeit with an unexpected trend that the lighter ones are more silicate-rich than the heavier ones. We feed the resulting planetesimals into a pebble accretion integrator with a continuous size distribution, finding that they undergo drastic changes in composition as they preferentially accrete silicate pebbles. The density and masses of large KBOs are best reproduced if they form between 15 and 22\,AU. Our solution avoids the timing problem because the first planetesimals are primarily icy, and $^{26}$Al is mostly incorporated in the slow phase of silicate pebble accretion. Our results lend further credibility to the streaming instability and pebble accretion as formation and growth mechanisms. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 24 pages, 13 figures, accepted to The Planetary Science Journal

Retention of CO Ice and Gas Within 486958 Arrokoth

Authors: Samuel P. D. Birch, Orkan M. Umurhan

Abstract: Kuiper Belt Objects (KBOs) represent some of the most ancient remnants of our solar system, having evaded significant thermal or evolutionary processing. This makes them important targets for exploration as they offer a unique opportunity to scrutinize materials that are remnants of the epoch of planet formation. Moreover, with recent and upcoming observations of KBOs, there is a growing interest… ▽ More Kuiper Belt Objects (KBOs) represent some of the most ancient remnants of our solar system, having evaded significant thermal or evolutionary processing. This makes them important targets for exploration as they offer a unique opportunity to scrutinize materials that are remnants of the epoch of planet formation. Moreover, with recent and upcoming observations of KBOs, there is a growing interest in understanding the extent to which these objects can preserve their most primitive, hypervolatile ices. Here, we present a theoretical framework that revisits this issue for small, cold classical KBOs like Arrokoth. Our analytical approach is consistent with prior studies but assumes an extreme cold end-member thermophysical regime for Arrokoth, enabling us to capture the essential physics without computationally expensive simulations. Under reasonable assumptions for interior temperatures, thermal conductivities, and permeabilities, we demonstrate that Arrokoth can retain its original CO stock for Gyrs if it was assembled long after the decay of radionuclides. The sublimation of CO ice generates an effective CO `atmosphere' within Arrokoth's porous matrix, which remains in near vapor-pressure equilibrium with the ice layer just below, thereby limiting CO loss. According to our findings, Arrokoth expels no more than $\approx 10^{22}$ particles s$^{-1}$, in agreement with upper limits inferred from \textit{New Horizons}' 2019 flyby observations. While our framework challenges recent predictions, it can serve as a benchmark for existing numerical models and be applied to future KBO observations from next-generation telescopes. △ Less

Submitted 1 March, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: Accepted on April 29, 2024 for publication in ICARUS special conference issue: ACM 14

Formation of the First Planetesimals via the Streaming Instability in Globally Turbulent Protoplanetary Disks?

Authors: Paul R. Estrada, Orkan M. Umurhan

Abstract: Using self-consistent models of turbulent particle growth in an evolving protoplanetary nebula of solar composition we find that recently proposed local metallicity and Stokes number criteria necessary for the streaming instability to generate gravitationally bound particle overdensities are generally not approached anywhere in the disk during the first million years, an epoch in which meteoritic… ▽ More Using self-consistent models of turbulent particle growth in an evolving protoplanetary nebula of solar composition we find that recently proposed local metallicity and Stokes number criteria necessary for the streaming instability to generate gravitationally bound particle overdensities are generally not approached anywhere in the disk during the first million years, an epoch in which meteoritic and observational evidence strongly suggests that the formation of the first planetesimals and perhaps giant planet core accretion is already occurring. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 14 pages, 4 figures, 1 appendix. Accepted to ApJ

Turbulence in particle laden midplane layers of planet forming disks

Authors: Debanjan Sengupta, Orkan M. Umurhan

Abstract: We examine the settled particle layers of planet forming disks in which the streaming instability (SI) is thought to be either weak or inactive. A suite of low-to-moderate resolution three-dimensional simulations in a $0.2H$ sized box, where $H$ is the pressure scale height, are performed using PENCIL for two Stokes numbers, \St$=0.04$ and $0.2$, at 1\% disk metallicity. We find a complex of Ekman… ▽ More We examine the settled particle layers of planet forming disks in which the streaming instability (SI) is thought to be either weak or inactive. A suite of low-to-moderate resolution three-dimensional simulations in a $0.2H$ sized box, where $H$ is the pressure scale height, are performed using PENCIL for two Stokes numbers, \St$=0.04$ and $0.2$, at 1\% disk metallicity. We find a complex of Ekman-layer jet-flows emerge subject to three co-acting linearly growing processes: (1) the Kelvin-Helmholtz instability (KHI), (2) the planet-forming disk analog of the baroclinic Symmetric Instability (SymI), and (3) a later-time weakly acting secondary transition process, possibly a manifestation of the SI, producing a radially propagating pattern state. For \St$=0.2$, KHI is dominant and manifests as off-midplane axisymmetric rolls, while for \St$=0.04$ the axisymmetric SymI mainly drives turbulence. SymI is analytically developed in a model disk flow, predicting that it becomes strongly active when the Richardson number (Ri) of the particle-gas midplane layer transitions below 1, exhibiting growth rates $\le\sqrt{2/\Ri - 2}\cdotΩ$, where $Ω$ is local disk rotation rate. For fairly general situations absent external sources of turbulence it is conjectured that the SI, when and if initiated, emerges out of a turbulent state primarily driven and shaped by at least SymI and/or KHI. We also find that turbulence produced in $256^3$ resolution simulations are not statistically converged and that corresponding $512^3$ simulations may be converged for \St$=0.2$. Furthermore, we report that our numerical simulations significantly dissipate turbulent kinetic energy on scales less than 6-8 grid points. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 55 pages, 27 figures, accepted for publication in ApJ

Global Modeling of Nebulae With Particle Growth, Drift, and Evaporation Fronts. II. The Influence of Porosity on Solids Evolution

Authors: Paul R. Estrada, Jeffrey N. Cuzzi, Orkan M. Umurhan

Abstract: Incremental particle growth in turbulent protoplanetary nebulae is limited by a combination of barriers that can slow or stall growth. Moreover, particles that grow massive enough to decouple from the gas are subject to inward radial drift which could lead to the depletion of most disk solids before planetesimals can form. Compact particle growth is probably not realistic. Rather, it is more likel… ▽ More Incremental particle growth in turbulent protoplanetary nebulae is limited by a combination of barriers that can slow or stall growth. Moreover, particles that grow massive enough to decouple from the gas are subject to inward radial drift which could lead to the depletion of most disk solids before planetesimals can form. Compact particle growth is probably not realistic. Rather, it is more likely that grains grow as fractal aggregates which may overcome this so-called radial drift barrier because they remain more coupled to the gas than compact particles of equal mass. We model fractal aggregate growth and compaction in a viscously evolving solar-like nebula for a range of turbulent intensities $α_{\rm{t}} = 10^{-5}-10^{-2}$. We do find that radial drift is less influential for porous aggregates over much of their growth phase; however, outside the water snowline fractal aggregates can grow to much larger masses with larger Stokes numbers more quickly than compact particles, leading to rapid inward radial drift. As a result, disk solids outside the snowline out to $\sim 10-20$ AU are depleted earlier than in compact growth models, but outside $\sim 20$ AU material is retained much longer because aggregate Stokes numbers there remain lower initially. Nevertheless, we conclude even fractal models will lose most disk solids without the intervention of some leap-frog planetesimal forming mechanism such as the Streaming Instability (SI), though conditions for the SI are generally never satisfied, except for a brief period %for a brief stage around $\sim 0.2$ Myr at the snowline for $α_{\rm{t}}=10^{-5}$. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 45 pages, 19 figures, 2 tables

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

Authors: G. Lesur, B. Ercolano, M. Flock, M. -K. Lin, C. -C. Yang, J. A. Barranco, P. Benitez-Llambay, J. Goodman, A. Johansen, H. Klahr, G. Laibe, W. Lyra, P. Marcus, R. P. Nelson, J. Squire, J. B. Simon, N. Turner, O. M. Umurhan, A. N. Youdin

Abstract: The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t… ▽ More The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

A Near Surface Temperature Model of Arrokoth

Authors: O. M. Umurhan, W. M. Grundy, M. K. Bird, R. Beyer, J. T. Keane, I. R. Linscott, S. Birch, C. Bierson, L. A. Young, S. A. Stern, C. M. Lisse, C. J. A. Howett, S. Protopapa, J. R. Spencer, R. P. Binzel, W. B. Mckinnon, T. R. Lauer, H. A. Weaver, C. B. Olkin, K. N. Singer, A. J. Verbiscer, A. H. Parker

Abstract: A near surface thermal model for Arrokoth is developed based on the recently released $10^5$ facet model of the body. This thermal solution takes into account Arrokoth's surface re-radiation back onto itself. The solution method exploits Arrokoth's periodic orbital character to develop a thermal response using a time-asymptotic solution method, which involves a Fourier transform solution of the he… ▽ More A near surface thermal model for Arrokoth is developed based on the recently released $10^5$ facet model of the body. This thermal solution takes into account Arrokoth's surface re-radiation back onto itself. The solution method exploits Arrokoth's periodic orbital character to develop a thermal response using a time-asymptotic solution method, which involves a Fourier transform solution of the heat equation, an approach recently used by others. We display detailed thermal solutions assuming that Arrokoth's near surface material's thermal inertia ${\cal I} = $ 2.5 W/m$^{-2}$K$^{-1}$s$^{1/2}$. We predict that at New Horizons' encounter with Arrokoth its encounter hemisphere surface temperatures were $\sim$ 57-59 K in its polar regions, 30-40 K on its equatorial zones, and 11-13 K for its winter hemisphere. Arrokoth's orbitally averaged temperatures are around 30-35 K in its polar regions, and closer to 40 K near its equatorial zones. Thermal reradiation from the surrounding surface amounts to less than 5\% of the total energy budget, while the total energy ensconced into and exhumed out Arrokoth's interior via thermal conduction over one orbit is about 0.5\% of the total energy budget. As a generalized application of this thermal modeling together with other KBO origins considerations, we favor the interpretation that New Horizons' REX instrument's $29 \pm 5$K brightness temperature measurement is consistent with Arrokoth's near surface material's being made of sub-to-few mm sized tholin-coated amorphous \water ice grains with 1 W/m$^{-2}$K$^{-1}$s$^{1/2}$ $< {\cal I} < $10-20 W/m$^{-2}$K$^{-1}$s$^{1/2}$, and which are characterized by an X-band emissivity in the range 0.9 and 1. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: Accepted on Feb. 21, 2022 for publication in Planetary Science Journal (PSJ)

Hypotheses for Triton's Plumes: New Analyses and Future Remote Sensing Tests

Authors: Jason D. Hofgartner, Samuel P. D. Birch, Julie Castillo, Will M. Grundy, Candice J. Hansen, Alexander G. Hayes, Carly J. A. Howett, Terry A. Hurford, Emily S. Martin, Karl L. Mitchell, Tom A. Nordheim, Michael J. Poston, Louise M. Prockter, Lynnae C. Quick, Paul Schenk, Rebecca N. Schindhelm, Orkan M. Umurhan

Abstract: At least two active plumes were observed on Neptune's moon Triton during the Voyager 2 flyby in 1989. Models for Triton's plumes have previously been grouped into five hypotheses, two of which are primarily atmospheric phenomena and are generally considered unlikely, and three of which include eruptive processes and are plausible. These hypotheses are compared, including new arguments, such as com… ▽ More At least two active plumes were observed on Neptune's moon Triton during the Voyager 2 flyby in 1989. Models for Triton's plumes have previously been grouped into five hypotheses, two of which are primarily atmospheric phenomena and are generally considered unlikely, and three of which include eruptive processes and are plausible. These hypotheses are compared, including new arguments, such as comparisons based on current understanding of Mars, Enceladus, and Pluto. An eruption model based on a solar-powered, solid-state greenhouse effect was previously considered the leading hypothesis for Triton's plumes, in part due to the proximity of the plumes to the subsolar latitude during the Voyager 2 flyby and the distribution of Triton's fans that are putatively deposits from former plumes. The other two eruption hypotheses are powered by internal heat, not solar insolation. Based on new analyses of the ostensible relation between the latitude of the subsolar point on Triton and the geographic locations of the plumes and fans, we argue that neither the locations of the plumes nor fans are strong evidence in favor of the solar-powered hypothesis. We conclude that all three eruption hypotheses should be considered further. Five tests are presented that could be implemented with remote sensing observations from future spacecraft to confidently distinguish among the eruption hypotheses for Triton's plumes. The five tests are based on the: (1) composition and thickness of Triton's southern hemisphere terrains, (2) composition of fan deposits, (3) distribution of active plumes, (4) distribution of fans, and (5) surface temperature at the locations of plumes and/or fans. The tests are independent, but complementary, and implementable with a single flyby mission such as the Trident mission concept. We note that, in the case of the solar-driven hypothesis, the 2030s and 2040s may be the last ... △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Icarus

A minimal model for vertical shear instability in protoplanetary accretion disks

Authors: Ron Yellin-Bergovoy, Eyal Heifetz, Orkan M. Umurhan

Abstract: The Vertical Shear Instability is an axisymmetric effect suggested to drive turbulence in the magnetically inactive zones of protoplanetary accretion disks. Here we examine its physical mechanism in analytically tractable ``minimal models" in three settings that include a uniform density fluid, a stratified atmosphere, and a shearing-box section of a protoplanetary disk. Each of these analyses sho… ▽ More The Vertical Shear Instability is an axisymmetric effect suggested to drive turbulence in the magnetically inactive zones of protoplanetary accretion disks. Here we examine its physical mechanism in analytically tractable ``minimal models" in three settings that include a uniform density fluid, a stratified atmosphere, and a shearing-box section of a protoplanetary disk. Each of these analyses show that the vertical shear instability's essence is similar to the slantwise convective symmetric instability in the mid-latitude Earth atmosphere, in the presence of vertical shear of the baroclinic jet stream, as well as mixing in the top layers of the Gulf Stream. We show that in order to obtain instability the fluid parcels' slope should exceed the slope of the mean absolute momentum in the disk radial-vertical plane. We provide a detailed and mutually self-consistent physical explanation from three perspectives: in terms of angular momentum conservation, as a dynamical interplay between a fluid's radial and azimuthal vorticity components, and from an energy perspective involving a generalized Solberg-Høiland Rayleigh condition. Furthermore, we explain why anelastic dynamics yield oscillatory unstable modes and isolate the oscillation mechanism from the instability one. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: accepted for publication

Initial results from the New Horizons exploration of 2014 MU69, a small Kuiper Belt Object

Authors: S. A. Stern, H. A. Weaver, J. R. Spencer, C. B. Olkin, G. R. Gladstone, W. M. Grundy, J. M. Moore, D. P. Cruikshank, H. A. Elliott, W. B. McKinnon, J. Wm. Parker, A. J. Verbiscer, L. A. Young, D. A. Aguilar, J. M. Albers, T. Andert, J. P. Andrews, F. Bagenal, M. E. Banks, B. A. Bauer, J. A. Bauman, K. E. Bechtold, C. B. Beddingfield, N. Behrooz, K. B. Beisser , et al. (180 additional authors not shown)

Abstract: The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a fl… ▽ More The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 43 pages, 8 figure

Journal ref: Science 364, eaaw9771 (2019)

The Geology and Geophysics of Kuiper Belt Object (486958) Arrokoth

Authors: J. R. Spencer, S. A. Stern, J. M. Moore, H. A. Weaver, K. N. Singer, C. B. Olkin, A. J. Verbiscer, W. B. McKinnon, J. Wm. Parker, R. A. Beyer, J. T. Keane, T. R. Lauer, S. B. Porter, O. L. White, B. J. Buratti, M. R. El-Maarry, C. M. Lisse, A. H. Parker, H. B. Throop, S. J. Robbins, O. M. Umurhan, R. P. Binzel, D. T. Britt, M. W. Buie, A. F. Cheng , et al. (53 additional authors not shown)

Abstract: The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger t… ▽ More The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Journal ref: Science, 367, aay3999 (2020)

The solar nebula origin of (486958) Arrokoth, a primordial contact binary in the Kuiper belt

Authors: W. B. McKinnon, D. C. Richardson, J. C. Marohnic, J. T. Keane, W. M. Grundy, D. P. Hamilton, D. Nesvorny, O. M. Umurhan, T. R. Lauer, K. N. Singer, S. A. Stern, H. A. Weaver, J. R. Spencer, M. W. Buie, J. M. Moore, J. J. Kavelaars, C. M. Lisse, X. Mao, A. H. Parker, S. B. Porter, M. R. Showalter, C. B. Olkin, D. P. Cruikshank, H. A. Elliott, G. R. Gladstone , et al. (4 additional authors not shown)

Abstract: The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collaps… ▽ More The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: Published in Science 28 Feb 2020 (First release 13 Feb 2020)

Journal ref: Science 367, eaay6620 (2020)

Color, Composition, and Thermal Environment of Kuiper Belt Object (486958) Arrokoth

Authors: W. M. Grundy, M. K. Bird, D. T. Britt, J. C. Cook, D. P. Cruikshank, C. J. A. Howett, S. Krijt, I. R. Linscott, C. B. Olkin, A. H. Parker, S. Protopapa, M. Ruaud, O. M. Umurhan, L. A. Young, C. M. Dalle Ore, J. J. Kavelaars, J. T. Keane, Y. J. Pendleton, S. B. Porter, F. Scipioni, J. R. Spencer, S. A. Stern, A. J. Verbiscer, H. A. Weaver, R. P. Binzel , et al. (24 additional authors not shown)

Abstract: The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrog… ▽ More The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: 31 pages, 8 figures

Streaming Instability in Turbulent Protoplanetary Disks

Authors: Orkan. M. Umurhan, Paul. R. Estrada, Jeffrey N. Cuzzi

Abstract: The streaming instability for solid particles in protoplanetary disks is re-examined assuming the familiar alpha ($α$) model for isotropic turbulence. Turbulence always reduces the growth rates of the streaming instability relative to values calculated for globally laminar disks. While for small values of the turbulence parameter, $α< 10^{-5}$, the wavelengths of the fastest-growing disturbances a… ▽ More The streaming instability for solid particles in protoplanetary disks is re-examined assuming the familiar alpha ($α$) model for isotropic turbulence. Turbulence always reduces the growth rates of the streaming instability relative to values calculated for globally laminar disks. While for small values of the turbulence parameter, $α< 10^{-5}$, the wavelengths of the fastest-growing disturbances are small fractions of the local gas vertical scale height $H$, we find that for moderate values of the turbulence parameter, i.e., $α\sim 10^{-5}-10^{-3}$, the lengthscales of maximally growing disturbances shift toward larger scales, approaching $H$. At these moderate turbulent intensities and for local particle to gas mass density ratios $ε< 0.5$, the vertical scales of the most unstable modes begin to exceed the corresponding radial scales so that the instability appears in the form of vertically oriented sheets extending well beyond the particle scale height. We find that for hydrodynamical turbulent disk models reported in the literature, with $α= 4\times 10^{-5} - 5\times 10^{-4}$, together with state of the art global evolution models of particle growth, the streaming instability is predicted to be viable within a narrow triangular patch of $α$--$τ_s$ parameter space centered on Stokes numbers, $τ_s \sim 0.01$ and $α\sim 4\times 10^{-5}$ and, further, exhibits growth rates on the order of several hundred to thousands of orbit times for disks with 1 percent ($Z= 0.01$) cosmic solids abundance or metallicity. Our results are consistent with, and place in context, published numerical studies of streaming instabilities. △ Less

Submitted 24 May, 2020; v1 submitted 12 June, 2019; originally announced June 2019.

Comments: Accepted to ApJ, April 14, 2020

Reorientation of Sputnik Planitia implies a Subsurface Ocean on Pluto

Authors: F. Nimmo, D. P. Hamilton, W. B. McKinnon P. M. Schenk, R. P. Binzel, C. J. Bierson, R. A. Beyer, J. M. Moore, S. A. Stern, H. A. Weaver, C. Olkin, L. A. Young, K. E. Smith, J. R. Spencer, M. Buie, B. Buratti, A. Cheng, D. Cruikshank, C. Dalle Ore, A. Earle, R. Gladstone, W. Grundy, A. D. Howard, T. Lauer, I. Linscott, J. Parker , et al. (38 additional authors not shown)

Abstract: The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topogra… ▽ More The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7]. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Journal ref: Pluto, Nature, 540, 94-96 (2016)

Convection in a volatile nitrogen-ice-rich layer drives Pluto's geological vigor

Authors: William B. McKinnon, Francis Nimmo, Teresa Wong, Paul M. Schenk, Oliver L. White, J. H. Roberts, J. M. Moore, J. R. Spencer, A. D. Howard, O. M. Umurhan, S. A. Stern, H. A. Weaver, C. B. Olkin, L. A. Young, K. E. Smith, R. Beyer, R. P. Binzel, M. Buie, B. Buratti, A. Cheng, D. Cruikshank, C. Dalle Ore, A. Earle, R. Gladstone, W. Grundy , et al. (39 additional authors not shown)

Abstract: The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on availa… ▽ More The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Journal ref: Nature, 534, 82-85 (2016)

Formation of Charon's Red Poles From Seasonally Cold-Trapped Volatiles

Authors: W. M. Grundy, D. P. Cruikshank, G. R. Gladstone, C. J. A. Howett, T. R. Lauer, J. R. Spencer, M. E. Summers, M. W. Buie, A. M. Earle, K. Ennico, J. Wm. Parker, S. B. Porter, K. N. Singer, S. A. Stern, A. J. Verbiscer, R. A. Beyer, R. P. Binzel, B. J. Buratti, J. C. Cook, C. M. Dalle Ore, C. B. Olkin, A. H. Parker, S. Protopapa, E. Quirico, K. D. Retherford , et al. (16 additional authors not shown)

Abstract: A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphe… ▽ More A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphere provides a potential feed stock for production of complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase showing the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon, the supply and temporary cold-trapping of material escaping from Pluto, and, while cold-trapped, its photolytic processing into more complex and less volatile molecules. The model results are consistent with the proposed mechanism producing the observed colour pattern on Charon. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Journal ref: Published in 2016 in Nature 539, 65-68

The CH4 cycles on Pluto over seasonal and astronomical timescales

Authors: T. Bertrand, F. Forget, O. M. Umurhan, J. M. Moore, L. A. Young, S. Protopapa, W. M. Grundy, B. Schmitt, R. D. Dhingra, R. P. Binzel, A. M. Earle, D. P. Cruikshank, S. A. Stern, H. A. Weaver, K. Ennico, C. B. Olkin

Abstract: New Horizons observations suggest that CH4 on Pluto has a complex history, involving reservoirs of different composition, thickness and stability controlled by volatile processes occurring on different timescales. In order to interpret these observations, we use a Pluto volatile transport model able to simulate the cycles of N2 and CH4 ices over millions of years. By assuming fixed solid mixing ra… ▽ More New Horizons observations suggest that CH4 on Pluto has a complex history, involving reservoirs of different composition, thickness and stability controlled by volatile processes occurring on different timescales. In order to interpret these observations, we use a Pluto volatile transport model able to simulate the cycles of N2 and CH4 ices over millions of years. By assuming fixed solid mixing ratios, we explore how changes in surface albedos, emissivities and thermal inertias impact volatile transport. This work is therefore a direct and natural continuation of the work by Bertrand et al. (2018), which only explored the N2 cycles. Results show that bright CH4 deposits can create cold traps for N2 ice outside Sputnik Planitia, leading to a strong coupling between the N2 and CH4 cycles. Depending on the assumed albedo for CH4 ice, the model predicts CH4 ice accumulation (1) at the same equatorial latitudes where the Bladed Terrain Deposits are observed, supporting the idea that these CH4-rich deposits are massive and perennial, or (2) at mid-latitudes (25°N-70°N), forming a thick mantle which is consistent with New Horizons observations. In our simulations, both CH4 ice reservoirs are not in an equilibrium state and either one can dominate the other over long timescales, depending on the assumptions made for the CH4 albedo. This suggests that long-term volatile transport exists between the observed reservoirs. The model also reproduces the formation of N2 deposits at mid-latitudes and in the equatorial depressions surrounding the Bladed Terrain, as observed by New Horizons. At the poles, only seasonal CH4 and N2 deposits are obtained in Pluto's current orbital configuration. Finally, we show that Pluto's atmosphere always contained, over the last astronomical cycles, enough gaseous CH4 to absorb most of the incoming Lyman-flux. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: Accepted in Icarus

Great Expectations: Plans and Predictions for New Horizons Encounter with Kuiper Belt Object 2014 MU69 ('Ultima Thule')

Authors: Jeffrey M. Moore, William B. McKinnon, Dale P. Cruikshank, G. Randall Gladstone, John R. Spencer, S. Alan Stern, Harold A. Weaver, Kelsi N. Singer, Mark R. Showalter, William M. Grundy, Ross A. Beyer, Oliver L. White, Richard P. Binzel, Marc W. Buie, Bonnie J. Buratti, Andrew F. Cheng, Carly Howett, Cathy B. Olkin, Alex H. Parker, Simon B. Porter, Paul M. Schenk, Henry B. Throop, Anne J. Verbiscer, Leslie A. Young, Susan D. Benecchi , et al. (9 additional authors not shown)

Abstract: The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far,… ▽ More The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity. △ Less

Submitted 6 August, 2018; originally announced August 2018.

The nitrogen cycles on Pluto over seasonal and astronomical timescales

Authors: T. Bertrand, F. Forget, O. M. Umurhan, W. M. Grundy, B. Schmitt, S. Protopapa, A. M. Zangari, O. L. White, P. M. Schenk, K. N. Singer, A. Stern, H. A. Weaver, L. A. Young, K. Ennico, C. B. Olkin, the New Horizons Science Team

Abstract: Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. The… ▽ More Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. These observations suggest a complex history involving volatile and glacial processes on different timescales. We present numerical simulations of volatile transport on Pluto performed with a model able to simulate the N2 cycle over millions of years (Myrs), taking into account the changes of obliquity and orbital parameters as experienced by Pluto. Results show that over one obliquity cycle, the latitudes of SP between 25°S-30°N are dominated by N2 condensation, while the latitudes between 30-50°N are dominated by N2 sublimation. We find that a net amount of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs. By comparing these results with the observed geology of SP, we can relate the formation of the pits and the brightness of the ice to the ice flux occurring at the annual timescale, while the glacial flows at its eastern edge and the erosion of the water ice mountains all around the ice sheet are related to the astronomical timescale. We also perform simulations with a glacial flow scheme which shows that SP is currently at its minimum extent. We also explore the stability of N2 ice outside SP. Results show that it is not stable at the poles but rather in the equatorial regions, in particular in depressions, where thick deposits may persist over tens of Myrs, before being trapped in SP. Finally, another key result is that the minimum and maximum surface pressures obtained over the simulated Myrs remain in the range of mm-Pa and Pa, respectively. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 49 pages, 16 figures

Modeling glacial flow on and onto Pluto's Sputnik Planitia

Authors: O. M. Umurhan, A. D. Howard, J. M. Moore, A. M. Earle, R. P. Binzel, S. A. Stern, P. M. Schenk, R. A. Beyer, O. L. White, F. NImmo, W. B. McKinnon, K. Ennico, C. B. Olkin, H. A. Weaver, L. A. Young

Abstract: Observations of Pluto's surface made by the New Horizons spacecraft indicates present-day nitrogen ice glaciation in and around the basin known as Sputnik Planum. Motivated by these observations, we have developed an evolutionary glacial flow model of solid nitrogen ice taking into account its published thermophysical and rheologies properties. This model assumes that glacial ice layers flow lamin… ▽ More Observations of Pluto's surface made by the New Horizons spacecraft indicates present-day nitrogen ice glaciation in and around the basin known as Sputnik Planum. Motivated by these observations, we have developed an evolutionary glacial flow model of solid nitrogen ice taking into account its published thermophysical and rheologies properties. This model assumes that glacial ice layers flow laminarly and have low aspect ratios which permits a vertically integrated mathematical formulation. We assess the conditions for the validity of laminar nitrogen ice motion by revisiting the problem of the onset of solid-state buoyant convection of nitrogen ice for a variety of bottom thermal boundary conditions. Subject to uncertainties in nitrogen ice rheology, nitrogen ice layers are estimated to flow laminarly for thicknesses less than 400-1000 meters. The resulting mass-flux formulation for when the nitrogen ice flows as a laminar dry glacier is characterized by an Arrhenius-Glen functional form. The flow model developed is used here to qualitatively answer some questions motivated by observed glacial flow features found on Sputnik Planum. We find that the wavy transverse dark features found along the northern shoreline of Sputnik Planum may be a transitory imprint of shallow topography just beneath the ice surface suggesting the possibility that a major shoreward flow event happened relatively recently within the last few hundred years. Model results also support the interpretation that the prominent darkened features resembling flow lobes observed along the eastern shoreline of the Sputnik Planum basin may be a result of wet nitrogen glacial ice flowing into the basin from the pitted highlands of eastern Tombaugh Regio. △ Less

Submitted 26 July, 2018; v1 submitted 18 June, 2016; originally announced June 2016.

Comments: This definitive version has several typos corrected including equation expressions in Appendix C.2 and C.3. Updated version with associated errata in process at Icarus (July 26, 2018). All previous results appearing in journal version --especially those pertaining to glacial flow and their related timescales (Sections 4 & 5) -- are unchanged

The Geology of Pluto and Charon Through the Eyes of New Horizons

Authors: Jeffrey M. Moore, William B. McKinnon, John R. Spencer, Alan D. Howard, Paul M. Schenk, Ross A. Beyer, Francis Nimmo, Kelsi N. Singer, Orkan M. Umurhan, Oliver L. White, S. Alan Stern, Kimberly Ennico, Cathy B. Olkin, Harold A. Weaver, Leslie A. Young, Richard P. Binzel, Marc W. Buie, Bonnie J. Buratti, Andrew F. Cheng, Dale P. Cruikshank, Will M. Grundy, Ivan R. Linscott, Harold J. Reitsema, Dennis C. Reuter, Mark R. Showalter , et al. (16 additional authors not shown)

Abstract: NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent tran… ▽ More NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt. △ Less

Submitted 19 April, 2016; originally announced April 2016.

Journal ref: Science 351, 1284 (2016)

Mean radius and shape of Pluto and Charon from New Horizons images

Authors: Francis Nimmo, Orkan M Umurhan, Carey M Lisse, Carver J Bierson, Tod R Lauer, Marc W Buie, Henry B Throop, Josh A Kammer, James H Roberts, William B McKinnon, Amanda M Zangari, Jeffrey M Moore, S Alan Stern, Leslie A Young, Harold A Weaver, Cathy B Olkin, Kim Ennico, the New Horizons GGI team

Abstract: Approach images taken by the LORRI imaging system during the New Horizons spacecraft encounter have been used to determine the mean radii and shapes of Pluto and Charon. The primary observations are limb locations derived using three independent approaches. The resulting mean radii of Pluto and Charon are 1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The corresponding densities a… ▽ More Approach images taken by the LORRI imaging system during the New Horizons spacecraft encounter have been used to determine the mean radii and shapes of Pluto and Charon. The primary observations are limb locations derived using three independent approaches. The resulting mean radii of Pluto and Charon are 1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The corresponding densities are 1854 +/- 11 kg m-3 and 1701 +/- 33 kg m-3 (2-sigma). The Charon radius value is consistent with previous Earth-based occultation estimates. The Pluto radius estimate is consistent with solar occultation measurements performed by the ALICE and Fine Sun Sensor instruments on New Horizons. Neither Pluto nor Charon show any evidence for tidal/rotational distortions; upper bounds on the oblateness are <0.6% and <0.5%, respectively. △ Less

Submitted 22 June, 2016; v1 submitted 2 March, 2016; originally announced March 2016.

Comments: revised version submitted to Icarus (Pluto special issue)

Vortex formation in protoplanetary discs induced by the vertical shear instability

Authors: Samuel Richard, Richard P. Nelson, Orkan M. Umurhan

Abstract: We present the results of 2D and 3D hydrodynamic simulations of idealized protoplanetary discs that examine the formation and evolution of vortices by the vertical shear instability (VSI). In agreement with recent work, we find that discs with radially decreasing temperature profiles and short thermal relaxation time-scales, are subject to the axisymmetric VSI. In three dimensions, the resulting v… ▽ More We present the results of 2D and 3D hydrodynamic simulations of idealized protoplanetary discs that examine the formation and evolution of vortices by the vertical shear instability (VSI). In agreement with recent work, we find that discs with radially decreasing temperature profiles and short thermal relaxation time-scales, are subject to the axisymmetric VSI. In three dimensions, the resulting velocity perturbations give rise to quasi-axisymmetric potential vorticity perturbations that break-up into discrete vortices, in a manner that is reminiscent of the Rossby wave instability. Discs with very short thermal evolution time-scales (i.e. τ<0.1 local orbit periods) develop strong vorticity perturbations that roll up into vortices that have small aspect ratios (χ<2) and short lifetimes (~ a few orbits). Longer thermal time-scales give rise to vortices with larger aspect ratios (6<χ<10), and lifetimes that depend on the entropy gradient. A steeply decreasing entropy profile leads to vortex lifetimes that exceed the simulation run times of hundreds of orbital periods. Vortex lifetimes in discs with positive or weakly decreasing entropy profiles are much shorter, being 10s of orbits at most, suggesting that the subcritical baroclinic instability plays an important role in sustaining vortices against destruction through the elliptical instability. Applied to the outer regions of protoplanetary discs, where the VSI is most likely to occur, our results suggest that vortices formed by the VSI are likely to be short lived structures. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Critical layers and protoplanetary disk turbulence

Authors: Orkan M. Umurhan, Karim Shariff, Jeffrey N. Cuzzi

Abstract: A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to… ▽ More A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability. △ Less

Submitted 9 July, 2016; v1 submitted 4 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ (25 pages). Results of previously appearing versions unchanged. Changes and additions include: updated references and nomenclature, viscous solutions added, higher resolution solutions displayed, expanded discussion, and shortened abstract

The Pluto system: Initial results from its exploration by New Horizons

Authors: S. A. Stern, F. Bagenal, K. Ennico, G. R. Gladstone, W. M. Grundy, W. B. McKinnon, J. M. Moore, C. B. Olkin, J. R. Spencer, H. A. Weaver, L. A. Young, T. Andert, J. Andrews, M. Banks, B. Bauer, J. Bauman, O. S. Barnouin, P. Bedini, K. Beisser, R. A. Beyer, S. Bhaskaran, R. P. Binzel, E. Birath, M. Bird, D. J. Bogan , et al. (126 additional authors not shown)

Abstract: The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly ext… ▽ More The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 8 pages - Initial Science paper from NASA's New Horizons Pluto Encounter

Journal ref: Science, vol 350, 292, 2015

Linear analysis of the vertical shear instability: outstanding issues and improved solutions (Research Note)

Authors: O. M. Umurhan, R. P. Nelson, O. Gressel

Abstract: The Vertical Shear Instability is one of two known mechanisms potentially active in the so-called dead zones of protoplanetary accretion disks. A recent analysis indicates that a subset of unstable modes shows unbounded growth - both as resolution is increased and when the nominal lid of the atmosphere is extended, possibly indicating ill-posedness in previous attempts of linear analysis. The redu… ▽ More The Vertical Shear Instability is one of two known mechanisms potentially active in the so-called dead zones of protoplanetary accretion disks. A recent analysis indicates that a subset of unstable modes shows unbounded growth - both as resolution is increased and when the nominal lid of the atmosphere is extended, possibly indicating ill-posedness in previous attempts of linear analysis. The reduced equations governing the instability are revisited and the generated solutions are examined using both the previously assumed separable forms and an improved non-separable solution form that is herewith introduced. Analyzing the reduced equations using the separable form shows that, while the low-order body modes have converged eigenvalues and eigenfunctions (as both the vertical boundaries of the atmosphere are extended and with increased radial resolution), it is also confirmed that the corresponding high-order body modes and the surface modes do indeed show unbounded growth rates. However, the energy contained in both the higher-order body modes and surface modes diminishes precipitously due to the disk's Gaussian density profile. Most of the energy of the instability is contained in the low-order modes. An inseparable solution form is introduced which filters out the inconsequential surface modes leaving only body modes (both low and high-order ones). The analysis predicts a fastest growing mode with a specific radial length scale. The growth rates associated with the fundamental corrugation and breathing modes matches the growth and length scales observed in previous nonlinear studies of the instability. △ Less

Submitted 28 May, 2015; v1 submitted 7 May, 2015; originally announced May 2015.

Comments: 8 pages, 5 figures, submitted to A+A as a Research Note (Version 2 reflects typo corrections to some formulae - results unaltered from previous version - still under review)

On the mechanism of self gravitating Rossby interfacial waves in proto-stellar accretion discs

Authors: Ron Yellin-Bergovoy, Eyal Heifetz, Orkan M. Umurhan

Abstract: The dynamical response of edge waves under the influence of self-gravity is examined in an idealized two-dimensional model of a proto-stellar disc, characterized in steady state as a rotating vertically infinite cylinder of fluid with constant density except for a single density interface at some radius r0. The fluid in basic state is prescribed to rotate with a Keplerian profile… ▽ More The dynamical response of edge waves under the influence of self-gravity is examined in an idealized two-dimensional model of a proto-stellar disc, characterized in steady state as a rotating vertically infinite cylinder of fluid with constant density except for a single density interface at some radius r0. The fluid in basic state is prescribed to rotate with a Keplerian profile $Ω_k(r)\sim r^{-3/2}$ modified by some additional azimuthal sheared flow. A linear analysis shows that there are two azimuthally propagating edge waves, kin to the familiar Rossby waves and surface gravity waves in terrestrial studies, which move opposite to one another with respect to the local basic state rotation rate at the interface. Instability only occurs if the radial pressure gradient is opposite to that of the density jump (unstably stratified) where self-gravity acts as a wave stabilizer irrespective of the stratification of the system. The propagation properties of the waves are discussed in detail in the language of vorticity edge waves. The roles of both Boussinesq and non-Boussinesq effects upon the stability and propagation of these waves with and without the inclusion of self-gravity are then quantified. The dynamics involved with self-gravity non- Boussinesq effect is shown to be a source of vorticity production where there is a jump in the basic state density, in addition, self-gravity also alters the dynamics via the radial main pressure gradient, which is a Boussinesq effect . Further applications of these mechanical insights are presented in the conclusion including the ways in which multiple density jumps or gaps may or may not be stable. △ Less

Submitted 12 January, 2016; v1 submitted 29 March, 2015; originally announced March 2015.

Comments: 20 pages, 5 figures and in GAFD article sytle

Intercomparison of General Circulation Models for Hot Extrasolar Planets

Authors: Inna Polichtchouk, James Y-K. Cho, Chris Watkins, Heidar Thor Thrastarson, Orkan M. Umurhan, Manuel de la Torre Juarez

Abstract: We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equa… ▽ More We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ different numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and `cubed-sphere' grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably different. When subject to a typical `hot-Jupiter'-like forcing, all five models show quantitatively different behavior -- although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst. △ Less

Submitted 20 November, 2013; originally announced November 2013.

Comments: 29 pages, 11 figures: Accepted for publication in Icarus

The Equations of Magnetoquasigeostrophy

Authors: O. M. Umurhan

Abstract: The dynamics contained in magnetized layers of exoplanet atmospheres are important to understand in order to characterize what observational signatures they may provide for future observations. It is important to develop a framework to begin studying and learning the physical processes possible under those conditions and what, if any, features contained in them may be observed in future observatio… ▽ More The dynamics contained in magnetized layers of exoplanet atmospheres are important to understand in order to characterize what observational signatures they may provide for future observations. It is important to develop a framework to begin studying and learning the physical processes possible under those conditions and what, if any, features contained in them may be observed in future observation missions. The aims of this study is to formally derive, from scaling arguments, a manageable reduced set of equations for analysis, i.e. a magnetic formulation of the equations of quasigeostrophy appropriate for a multi-layer atmosphere. The main goal is to provide a simpler theoretical platform to explore the dynamics possible within confined magnetized layers of exoplanet atmospheres. We primarily use scaling arguments to derive the reduced equations of "magnetoquasigeostrophy" which assumes dynamics to take place in an atmospheric layer which is vertically thin compared to its horizontal scales. The derived equation set retains features existing in standard shallow-water magnetohydrodynamic equations but are absent in more classical derivations of the quasi-geostrophic limit, namely, the non-divergence of the in-plane components of the magnetic field. We liken this non-divergence of the in-plane magnetic fields as indicative of a quantity whose behaviour mimics a two-dimensional "pseudo"-magnetic monopole source. We also find, using the same scaling argument procedures, appropriate limits of the fundamental parameters of the system which yield reduced equations describing the flow dynamics primarily characterized by magnetostrophic balance. The standard scaling arguments employed here show how traditional magnetized quasigeostrophic equations connect to their magnetized shallow water forms. The equations derived are amenable to analysis using well-known techniques. △ Less

Submitted 2 January, 2013; originally announced January 2013.

Comments: 13 pages. Under consideration for publication in Astronomy and Astrophysics

Non-Dissipative Saturation of the Magnetorotational Instability in Thin Disks

Authors: Edward Liverts, Yuri Shtemler, Michael Mond, Orkan M. Umurhan, Dmitry V. Bisikalo

Abstract: A new non-dissipative mechanism is proposed for the saturation of the axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that are subject to an axial magnetic field. That mechanism relies on the energy transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction is enabled due to the vertical stratification of the disk that results in the discretization of i… ▽ More A new non-dissipative mechanism is proposed for the saturation of the axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that are subject to an axial magnetic field. That mechanism relies on the energy transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction is enabled due to the vertical stratification of the disk that results in the discretization of its MRI spectrum, as well as by applying the appropriate boundary conditions. A second order Duffing-like amplitude equation for the initially unstable MRI modes is derived. The solutions of that equation exhibit bursty nonlinear oscillations with a constant amplitude that signifies the saturation level of the MRI. Those results are verified by a direct numerical solution of the full nonlinear reduced set of thin disk magnetohydrodynamics equations. △ Less

Submitted 19 October, 2012; originally announced October 2012.

Comments: 4 pages, 3 figures; Phys. Rev. Lett. accepted

Journal ref: Phys. Rev. Lett. 109, 224501, 2012

Linear and nonlinear evolution of the vertical shear instability in accretion discs

Authors: Richard P. Nelson, Oliver Gressel, Orkan M. Umurhan

Abstract: (Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as… ▽ More (Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear' linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance. △ Less

Submitted 12 September, 2012; originally announced September 2012.

Comments: 19 pages, 15 figures, submitted to MNRAS for publication. A version with full resolution figures is available at http://www.maths.qmul.ac.uk/~rpn/Corrugation/NGU.pdf

Potential vorticity dynamics in the framework of disk shallow-water theory: II. Mixed Barotropic-Baroclinic Instability

Authors: O. M. Umurhan

Abstract: We extend exploration of potential vorticity instabilities in cold astrophysical disks whose mean states are baroclinic. In particular, we seek to demonstrate the potential existence of traditional baroclinic instabilities of meteorological studies in a simplified two-layer Philips Disk Model. Each disk layer is of constant but differing densities. The resulting mean azimuthal velocity profile sho… ▽ More We extend exploration of potential vorticity instabilities in cold astrophysical disks whose mean states are baroclinic. In particular, we seek to demonstrate the potential existence of traditional baroclinic instabilities of meteorological studies in a simplified two-layer Philips Disk Model. Each disk layer is of constant but differing densities. The resulting mean azimuthal velocity profile shows a variation in the vertical direction implying that the system is baroclinic in the mean state. The stability of the system is treated in the context of disk shallow water theory wherein azimuthal disturbances are much longer than the corresponding radial or vertical scales. The normal-mode problem is solved numerically using two different methods. The results of a symmetric single layer barotropic model is considered and it is found that instability persists for models in which the potential vorticity profiles are not symmetric, consistent with previous results. The instaiblity is interpreted in terms of interacting Rossby waves. For a two layer system in which the flow is fundamentally baroclinic we report here that instability takes on the form of mixed barotropic-baroclinic type: instability occurs but it qualitatively follows the pattern of instability found in the barotropic models. Instability arises because of the phase locking and interaction of the Rossby waves between the two layers. The strength of the instability weakens as the density contrast between layers increases. (For full abstract see article.) △ Less

Submitted 7 March, 2012; originally announced March 2012.

Comments: Under consideration for publication in Astronomy and Astrophysics. Comments to the author welcome

Potential vorticity dynamics in the framework of disk shallow-water theory: I. The Rossby wave instability

Authors: O. M. Umurhan

Abstract: The Rossby wave instability in astrophysical disks is as a potentially important mechanism for driving angular momentum transport in disks. We aim to understand this instability in an approximate three-dimensional disk model environment which we assume to be a single homentropic annular layer we analyze using disk shallow-water theory. We consider the normal mode stability analysis of two kinds of… ▽ More The Rossby wave instability in astrophysical disks is as a potentially important mechanism for driving angular momentum transport in disks. We aim to understand this instability in an approximate three-dimensional disk model environment which we assume to be a single homentropic annular layer we analyze using disk shallow-water theory. We consider the normal mode stability analysis of two kinds of radial profiles of the mean potential vorticity: The first type is a single step and the second kind is a symmetrical step of finite width describing either a localized depression or peak of the mean potential vorticity. For single potential vorticity steps we find there is no instability. There is no instability when the symmetric step is a localized peak. However, the Rossby wave instability occurs when the symmetrical step profile is a depression, which, in turn, corresponds to localized peaks in the mean enthalpy profile. This is in qualitative agreement with previous two-dimensional investigations of the instability. For all potential vorticity depressions, instability occurs for regions narrower than some maximum radial length scale. We interpret the instability as resulting from the interaction of at least two Rossby edgewaves. We identify the Rossby wave instability in the restricted three-dimensional framework of disk shallow water theory. Additional examinations of generalized barotropic flows are needed. Viewing disk vortical instabilities from the conceptual perspective of interacting edgewaves can be useful. △ Less

Submitted 12 August, 2010; originally announced August 2010.

Comments: To appear in Astronomy and Astrophysics

Non-exponential hydrodynamical growth in density-stratified thin Keplerian discs

Authors: Yu. M. Shtemler, M. Mond, G. Rudiger, O. Regev, O. M. Umurhan

Abstract: The short time evolution of three dimensional small perturbations is studied. Exhibiting spectral asymptotic stability, thin discs are nonetheless shown to host intensive hydrodynamical activity in the shape of non modal growth of initial small perturbations. Two mechanisms that lead to such behavior are identified and studied, namely, non-resonant excitation of vertically confined sound waves b… ▽ More The short time evolution of three dimensional small perturbations is studied. Exhibiting spectral asymptotic stability, thin discs are nonetheless shown to host intensive hydrodynamical activity in the shape of non modal growth of initial small perturbations. Two mechanisms that lead to such behavior are identified and studied, namely, non-resonant excitation of vertically confined sound waves by stable planar inertia-coriolis modes that results in linear growth with time, as well as resonant coupling of those two modes that leads to a quadratic growth of the initial perturbations. It is further speculated that the non modal growth can give rise to secondary strato-rotational instabilities and thus lead to a new route to turbulence generation in thin discs. △ Less

Submitted 12 March, 2010; v1 submitted 8 February, 2010; originally announced February 2010.

Low magnetic-Prandtl number flow configurations for cold astrophysical disk models: speculation and analysis

Authors: O. M. Umurhan

Abstract: Simulations of astrophysical disks in the shearing box that are subject to the magnetorotational instability (MRI) show that activity appears to be reduced as the magnetic Prandtl number (Pm) is lowered. On the other hand, calculations for laboratory experiments show that saturation is achieved through modification of the background shear for Pm << 1. Guided by the results of calculations approp… ▽ More Simulations of astrophysical disks in the shearing box that are subject to the magnetorotational instability (MRI) show that activity appears to be reduced as the magnetic Prandtl number (Pm) is lowered. On the other hand, calculations for laboratory experiments show that saturation is achieved through modification of the background shear for Pm << 1. Guided by the results of calculations appropriate for laboratory experiments when Pm is very low, the axisymmetric stability of inviscid disturbances in a shearing box model immersed in a constant vertical background magnetic field is considered under a variety of shear profiles and boundary conditions in order to evaluate the hypothesis that modifications of the shear bring about saturation of the instability. It is found that the emergence and stability of the MRI is sensitive to the boundary conditions adopted. Channel modes do not appear to be stabilized through any modification of the background shear whose average remains Keplerian. However, systems that have non-penetrative boundaries can saturate the MRI through modification of the background shear. Conceptually equating the qualitative results from laboratory experiments to the conditions in a disk may therefore be misleading. △ Less

Submitted 2 February, 2010; originally announced February 2010.

Comments: To Appear in Astronomy and Astrophysics

Global transient dynamics of three-dimensional hydrodynamical disturbances in a thin viscous accretion disk

Authors: Paola Rebusco, Orkan M. Umurhan, Wlodek Kluzniak, Oded Regev

Abstract: Thin viscous Keplerian accretion disks are considered asymptotically stable, even though they can show significant dynamic activity on short timescales. In this paper the dynamics of non-axisymmetric hydrodynamical disturbances of disks are investigated analytically building upon the steady state three-dimensional structure and evolution of axisymmetric perturbations explored in previous work. A… ▽ More Thin viscous Keplerian accretion disks are considered asymptotically stable, even though they can show significant dynamic activity on short timescales. In this paper the dynamics of non-axisymmetric hydrodynamical disturbances of disks are investigated analytically building upon the steady state three-dimensional structure and evolution of axisymmetric perturbations explored in previous work. Assuming a polytropic equation of state solutions are found by means of an asymptotic expansion in the small parameter measuring the ratio of the disk thickness to characteristic radius. In-depth analysis shows that every perturbation that disturbs the radial velocity induces significant transient growth in the (acoustic) energy of the evolving disturbance. This effect is most evident in the density and vertical velocity. The transient growth observed is tied to the non-separable nature of the solutions where, in particular, pattern evolution is controlled by a similarity variable composed of the radial coordinate and time. This leads to growing winding perturbations that display successive radial peaks and troughs. We argue that these transient non-axisymmetric structures may precipitate secondary instabilities which, consequently, may be a critical element for a new alternative picture of turbulence arousal in non-magnetized astrophysical disks. △ Less

Submitted 11 June, 2009; v1 submitted 1 June, 2009; originally announced June 2009.

Comments: Accepted in Physics of Fluids

Linear dynamics of weakly viscous accretion disks: A disk analog of Tollmien-Schlichting waves

Authors: O. M. Umurhan, G. Shaviv

Abstract: This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is much larger than the radial and vertical scales. We examine the sta… ▽ More This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is much larger than the radial and vertical scales. We examine the stability of a non-axisymmetric potential vorticity mode, i.e. a PV-anomaly. in a configuration in which buoyant convection and the strato-rotational instability do not to operate. We consider a series of boundary conditions which show the PV-anomaly to be unstable both on a finite and semi-infinite radial domains. We find these conditions leading to an instability which is the disk analog of Tollmien-Schlichting waves. When the viscosity is weak, evidence of the instability is most pronounced by the emergence of a vortex sheet at the critical layer located away from the boundary where the instability is generated. For some boundary conditions a necessary criterion for the onset of instability for vertical wavelengths that are a sizable fraction of the layer's thickness and when the viscosity is small is that the appropriate Froude number of the flow be greater than one. This instability persists if more realistic boundary conditions are applied, although the criterion on the Froude number is more complicated. The unstable waves studied here share qualitative features to the instability seen in rotating Blasius boundary layers. The implications of these results are discussed. An overall new strategy for exploring and interpreting disk instability mechanisms is also suggested. △ Less

Submitted 4 November, 2008; v1 submitted 30 October, 2008; originally announced October 2008.

Comments: Accepted for publication in Astronomy and Astrophysics. 18 pages. This version 3 with corrected style file

A shallow-water theory for annular sections of Keplerian Disks

Authors: O. M. Umurhan

Abstract: A scaling argument is presented that leads to a shallow water theory of non-axisymmetric disturbances in annular sections of thin Keplerian disks. To develop a theoretical construction that will aid in physically understanding the relationship of known two-dimensional vortex dynamics to their three-dimensional counterparts in Keplerian disks. Using asymptotic scaling arguments varicose disturban… ▽ More A scaling argument is presented that leads to a shallow water theory of non-axisymmetric disturbances in annular sections of thin Keplerian disks. To develop a theoretical construction that will aid in physically understanding the relationship of known two-dimensional vortex dynamics to their three-dimensional counterparts in Keplerian disks. Using asymptotic scaling arguments varicose disturbances of a Keplerian disk are considered on radial and vertical scales consistent with the height of the disk while the azimuthal scales are the full $2π$ angular extent of the disk. The scalings lead to dynamics which are radially geostrophic and vertically hydrostatic. It follows that a potential vorticity quantity emerges and is shown to be conserved in a Lagrangian sense. Uniform potential vorticity linear solutions are explored and the theory is shown to contain an incarnation of the strato-rotational instability under channel flow conditions. Linearized solutions of a single defect on an infinite domain is developed and is shown to support a propagating Rossby edgewave. Linear non-uniform potential vorticity solutions are also developed and are shown to be similar in some respects to the dynamics of strictly two-dimensional inviscid flows. Based on the framework of this theory, arguments based on geophysical notions are presented to support the assertion that the strato-rotational instability is in a generic class of barotropic/baroclinic potential vorticity instabilities. Extensions of this formalism are also proposed. The shallow water formulation achieved by the asymptotic theory developed here opens a new approach to studying disk dynamics. △ Less

Submitted 22 September, 2008; v1 submitted 23 February, 2008; originally announced February 2008.

Comments: Accepted (July 21, 2008), now in final form

The instability of counter-propagating kernel gravity waves in a constant shear flow

Authors: O. M. Umurhan, E. Heifetz, N. Harnik, F. Lott

Abstract: The mechanism describing the recently developed notion of kernel gravity waves (KGWs) is reviewed and such structures are employed to interpret the unstable dynamics of an example stratified plane parallel shear flow. This flow has constant vertical shear, is infinite in the vertical extent, and characterized by two density jumps of equal magnitude each decreasing successively with height, in wh… ▽ More The mechanism describing the recently developed notion of kernel gravity waves (KGWs) is reviewed and such structures are employed to interpret the unstable dynamics of an example stratified plane parallel shear flow. This flow has constant vertical shear, is infinite in the vertical extent, and characterized by two density jumps of equal magnitude each decreasing successively with height, in which the jumps are located symmetrically away from the midplane of the system. We find that for a suitably defined bulk-Richardson number there exists a band of horizontal wavenumbers which exhibits normal-mode instability. The instability mechanism closely parallels the mechanism responsible for the instability seen in the problem of counter-propagating Rossby waves. In this problem the instability arises out of the interaction of counter-propagating gravity waves. We argue that the instability meets the Hayashi-Young criterion for wave instability. We also argue that the instability is the simplest one that can arise in a stratified atmosphere with constant shear flow. The counter propagating gravity waves mechanism detailed here explains why the Rayleigh criteria for shear flow instability in the unstratified case does not need to be satisfied in the stratified case. This illustrates how the Miles-Howard theorem may support destabilization through stratification. A normal mode analysis of a foamy layer consisting of two density jumps of unequal magnitude is also analyzed. The results are considered in terms of observations made of sea-hurricane interfaces. △ Less

Submitted 29 November, 2007; originally announced November 2007.

Comments: Submitted to Physics of Fluids, 13 pages

On the viability of the shearing box approximation for numerical studies of MHD turbulence in accretion disks

Authors: O. Regev, O. M. Umurhan

Abstract: Most of our knowledge on the nonlinear development of the magneto-rotational instability (MRI) relies on the results of numerical simulations employing the shearing box (SB) approximation. A number of difficulties arising from this approach have recently been pointed out in the literature. We thoroughly examine the effects of the assumptions made and numerical techniques employed in SB simulatio… ▽ More Most of our knowledge on the nonlinear development of the magneto-rotational instability (MRI) relies on the results of numerical simulations employing the shearing box (SB) approximation. A number of difficulties arising from this approach have recently been pointed out in the literature. We thoroughly examine the effects of the assumptions made and numerical techniques employed in SB simulations. This is done in order to clarify and gain better understanding of those difficulties as well as of a number of additional serious problems, raised here for the first time, and of their impact on the results. Analytical derivations and estimates as well as comparative analysis to methods used in the numerical study of turbulence are used. Numerical experiments are performed to support some of our claims and conjectures. The following problems, arising from the (virtually exclusive) use of the SB simulations as a tool for the understanding and quantification of the nonlinear MRI development in disks, are analyzed and discussed: (i) inconsistencies in the application of the SB approximation itself; (ii) the limited spatial scale of the SB; (iii) the lack of convergence of most ideal MHD simulations; (iv) side-effects of the SB symmetry and the non-trivial nature of the linear MRI; (v) physical artifacts arising on the too small box scale due to periodic boundary conditions. The computational and theoretical challenge posed by the MHD turbulence problem in accretion disks cannot be met by the SB approximation, as it has been used to date. A new strategy to confront this challenge is proposed, based on techniques widely used in numerical studies of turbulent flows - developing (e.g., with the help of local numerical studies) a sub-grid turbulence model and implementing it in global calculations. △ Less

Submitted 27 February, 2008; v1 submitted 5 November, 2007; originally announced November 2007.

Comments: Accepted for publication in Astronomy and Astrophysics

Hydrodynamic response of rotationally supported flows in the Small Shearing Box model

Authors: A. Sternberg, O. M. Umurhan, Y. Gil, O. Regev

Abstract: The hydrodynamic response of the inviscid small shearing box model of a midplane section of a rotationally supported astrophysical disk is examined. An energy functional ${\cal E}$ is formulated for the general nonlinear problem. It is found that the fate of disturbances is related to the conservation of this quantity which, in turn, depends on the boundary conditions utilized: ${\cal E}$ is con… ▽ More The hydrodynamic response of the inviscid small shearing box model of a midplane section of a rotationally supported astrophysical disk is examined. An energy functional ${\cal E}$ is formulated for the general nonlinear problem. It is found that the fate of disturbances is related to the conservation of this quantity which, in turn, depends on the boundary conditions utilized: ${\cal E}$ is conserved for channel boundary conditions while it is not conserved in general for shearing box conditions. Linearized disturbances subject to channel boundary conditions have normal-modes described by Bessel Functions and are qualitatively governed by a quantity $Σ$ which is a measure of the ratio between the azimuthal and vertical wavelengths. Inertial oscillations ensue if $Σ>1$ - otherwise disturbances must in general be treated as an initial value problem. We reflect upon these results and offer a speculation. △ Less

Submitted 4 March, 2008; v1 submitted 14 March, 2007; originally announced March 2007.

Comments: 6 pages, resubmitted to Astronomy and Astrophysics, shortened with references added

On the nonlinear saturation of the magnetorotational instability near threshold in a thin-gap Taylor-Couette setup

Authors: O. M. Umurhan, O. Regev, K. Menou

Abstract: We study the saturation near threshold of the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible fluid in a thin-gap Taylor-Couette configuration. A vertical magnetic field, Keplerian shear and no-slip, conducting radial boundary conditions are adopted. The weakly non-linear theory leads to a real Ginzburg-Landau equation for the disturbance amplitude, like… ▽ More We study the saturation near threshold of the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible fluid in a thin-gap Taylor-Couette configuration. A vertical magnetic field, Keplerian shear and no-slip, conducting radial boundary conditions are adopted. The weakly non-linear theory leads to a real Ginzburg-Landau equation for the disturbance amplitude, like in our previous idealized analysis. For small magnetic Prandtl number (P \ll 1), the saturation amplitude scales as P^{2/3} while the magnitude of angular momentum transport scales as P^{4/3}. The difference with the previous scalings (~P^{1/2} and P respectively) is attributed to the emergence of radial boundary layers. Away from those, steady-state non-linear saturation is achieved through a modest reduction in the destabilizing shear. These results will be useful to understand MRI laboratory experiments and associated numerical simulations. △ Less

Submitted 1 July, 2007; v1 submitted 28 February, 2007; originally announced February 2007.

Comments: Accepted with minor changes to Physical Review E

An idealized model for dust-gas interaction in a rotating channel

Authors: O. M. Umurhan

Abstract: A 2D model representing the dynamical interaction of dust and gas in a planetary channel is explored. The two components are treated as interpenetrating fluids in which the gas is treated as a Boussinesq fluid while the dust is treated as pressureless. The only coupling between both fluid states is kinematic drag. The channel gas experiences a temperature gradient in the spanwise direction and i… ▽ More A 2D model representing the dynamical interaction of dust and gas in a planetary channel is explored. The two components are treated as interpenetrating fluids in which the gas is treated as a Boussinesq fluid while the dust is treated as pressureless. The only coupling between both fluid states is kinematic drag. The channel gas experiences a temperature gradient in the spanwise direction and it is adverse the constant force of gravity. The latter effects only the gas and not the dust component which is considered to free float in the fluid. The channel is also considered on an f-plane so that the background vorticity gradient can cause any emerging vortex structure to drift like a Rossby wave. A linear theory analysis is explored and a nonlinear amplitude theory is developed for disturbances of this arrangement. It is found that the presence of the dust can help generate and shape emerging convection patterns and dynamics in the gas so long as the state of the gas exceeds a suitably defined Rayleigh number appropriate for describing drag effects. In the linear stage the dust particles collect quickly onto sites in the gas where the vorticity is minimal, i.e. where the disturbance vorticity is anticylonic which is consistent with previous studies. The nonlinear theory shows that, in turn, the local enhancement of dust concentration in the gas effects the vigor of the emerging convective roll by modifying the effective local Rayleigh number of the fluid. It is also found that without the f-plane approximation built into the model the dynamics there is an algebraic runaway caused by unrestrained growth in the dust concentration. The background vorticity gradient forces the convective roll to drift like a Rossby wave and this causes the dust concentration enhancements to not runaway. △ Less

Submitted 11 January, 2007; v1 submitted 9 January, 2007; originally announced January 2007.

Comments: Updated with more current references

A weakly nonlinear analysis of the magnetorotational instability in a model channel flow

Authors: O. M. Umurhan, K. Menou, O. Regev

Abstract: We show by means of a perturbative weakly nonlinear analysis that the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible rotating shear flow in a thin channel gives rise to a real Ginzburg-Landau equation for the disturbance amplitude. For small magnetic Prandtl number (${\cal P}_{\rm m}$), the saturation amplitude is $\propto \sqrt{{\cal P}_{\rm m}}$ and th… ▽ More We show by means of a perturbative weakly nonlinear analysis that the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible rotating shear flow in a thin channel gives rise to a real Ginzburg-Landau equation for the disturbance amplitude. For small magnetic Prandtl number (${\cal P}_{\rm m}$), the saturation amplitude is $\propto \sqrt{{\cal P}_{\rm m}}$ and the resulting momentum transport scales as ${\cal R}^{-1}$, where $\cal R$ is the {\em hydrodynamic} Reynolds number. Simplifying assumptions, such as linear shear base flow, mathematically expedient boundary conditions and continuous spectrum of the vertical linear modes, are used to facilitate this analysis. The asymptotic results are shown to comply with numerical calculations using a spectral code. They suggest that the transport due to the nonlinearly developed MRI may be very small in experimental setups with ${\cal P}_{\rm m} \ll 1$. △ Less

Submitted 8 January, 2007; v1 submitted 17 August, 2006; originally announced August 2006.

Comments: Accepted to Physical Review Letters - Nov. 30, 2006. In final form

Journal ref: Phys.Rev.Lett.98:034501,2007

On the stratorotational instability in the quasi-hydrostatic semi-geostrophic limit

Authors: O. M. Umurhan

Abstract: The linear normal-mode stratorotational instability (SRI) is analytically reexamined in the inviscid limit where the length scales of horizontal disturbances are large compared their vertical and radial counterparts. Boundary conditions different than channel walls are also considered. This quasi-hydrostatic, semi-geostrophic (QHSG) approximation allows one to examine the effect of a vertically… ▽ More The linear normal-mode stratorotational instability (SRI) is analytically reexamined in the inviscid limit where the length scales of horizontal disturbances are large compared their vertical and radial counterparts. Boundary conditions different than channel walls are also considered. This quasi-hydrostatic, semi-geostrophic (QHSG) approximation allows one to examine the effect of a vertically varying Brunt-Vaisaila frequency, $N^2$. It is found that the normal-mode instability persists when $N^2$ increases quadratically with respect to the disc vertical coordinate. However we also find that the SRI seems to exist in this inviscid QHSG extreme only for channel wall conditions: when one or both of the reflecting walls are removed there is no instability in the asymptotic limit explored here. It is also found that only exponential-type SRI modes (as defined by Dubrulle et al. 2005) exist under these conditions. These equations also admit non-normal mode behaviour. Fixed Lagrangian pressure conditions on both radial boundaries predicts there to be no normal mode behaviour in the QHSG limit. The mathematical relationship between the results obtained here and that of the classic Eady (1949) problem for baroclinic instability is drawn. We conjecture as to the mathematical/physical nature of the SRI. The general linear problem, analyzed without approximation in the context of the Boussinesq equations, admits a potential vorticity-like quantity that is advectively conserved by the shear. Its existence means that a continuous spectrum \emph{is a generic feature of this system}. It also implies that in places where the Brunt-Vaisaila frequency becomes dominant the linearized flow may two-dimensionalize by advectively conserving its vertical vorticity. △ Less

Submitted 13 September, 2005; v1 submitted 1 June, 2005; originally announced June 2005.

Comments: 16 pages. Accepted in MNRAS (09/05). New sections added and abstract changed

Journal ref: Mon.Not.Roy.Astron.Soc. 365 (2006) 85-100

Hydrodynamical stability of thin accretion discs: transient growth of global axisymmetric perturbations

Authors: O. M. Umurhan, A. Nemirovsky, O. Regev, G. Shaviv

Abstract: The purpose of this paper is to explore how accretion discs manifest the phenomenon of transient growth on a global scale. We investigate analytically the time response of a thin accretion disc to particular axisymmetric perturbations. To facilitate an analytical treatment we replace the energy equation with a general polytropic assumption. The asymptotic expansion of Kluźniak & Kita (2000), whi… ▽ More The purpose of this paper is to explore how accretion discs manifest the phenomenon of transient growth on a global scale. We investigate analytically the time response of a thin accretion disc to particular axisymmetric perturbations. To facilitate an analytical treatment we replace the energy equation with a general polytropic assumption. The asymptotic expansion of Kluźniak & Kita (2000), which extended the method of Regev (1983) to a full steady polytropic disc (with $n=3/2$), is further developed and implemented for both the steady (for any polytropic index) and time-dependent problems. The spatial form and temporal behaviour of selected dynamical disturbances are studied in detail. We identify the perturbation space which leads to transient growth and provide analytical solutions which manifest this expected transient growth behaviour. Three terms (physical causes) responsible for the appearance of transient growth are identified. Two depend explicitly on the viscosity while the third one is relevant also for inviscid discs. The main conclusion we draw is that the phenomenon of transient growth exists in discs on a global scale. △ Less

Submitted 20 February, 2005; originally announced February 2005.

Comments: 20 pages, 4 figues (last figure in jpg format). Submitted to MNRAS

On the nature of the hydrodynamic stability of accretion disks

Authors: O. M. Umurhan, G. Shaviv

Abstract: The linear stability of accretion disks is revisited. The governing equations are expanded asymptotically and solved to first order in the expansion parameter $ε$ defined by the ratio of the disk's vertical thickness to its radial extent. Algebraically growing solutions are found for global perturbations on the radial accretion flow of thin inviscid compressible Keplerian disks. The algebraic te… ▽ More The linear stability of accretion disks is revisited. The governing equations are expanded asymptotically and solved to first order in the expansion parameter $ε$ defined by the ratio of the disk's vertical thickness to its radial extent. Algebraically growing solutions are found for global perturbations on the radial accretion flow of thin inviscid compressible Keplerian disks. The algebraic temporal behavior is exhibited in the vertical velocities and the thermodynamic variables and has the form $t\sinΩ_0 t$ locally in the disk where $Ω_0$ is the Keplerian rotation rate. The physical implications and relations to the Solberg-Hoiland stability criteria are discussed. △ Less

Submitted 30 January, 2005; v1 submitted 15 November, 2004; originally announced November 2004.

Comments: Submitted to A+A Letters on November 15, 2004. 4 pages, no figures. Accepted January 20, 2005. Section 3 rewritten and references updated