Orkan M. Umurhan

Linear and non-linear evolution of the vertical shear instability in accretion discs

We analyse the stability and non-linear dynamics of power-law accretion disc models. These have mid-plane densities that follow radial power laws and have either temperature or entropy distributions that are strict power-law functions of cylindrical radius, R. We employ two different hydrodynamic codes to perform high-resolution 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 disc scaleheight, the thermal relaxation time of the fluid and the disc viscosity. We present an accompanying stability analysis of the problem, based on asymptotic methods, that we use to guide our interpretation of the simulation results. We find that axisymmetric disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of perturbations whose most obvious character is modes with horizontal and vertical wavenumbers that satisfy vertical bar k(R)/k(Z)vertical bar 1. Instability occurs only 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 appear to exhibit the Goldreich-Schubert-Fricke or vertical shear linear instability. Closer inspection of the simulation results uncovers the growth of two distinct modes. The first are characterized by very short radial wavelength perturbations that grow rapidly at high latitudes in the disc, and descend down towards the mid-plane on longer time-scales. We refer to these as finger modes because they display k(R)/k(Z) 1. The second appear at slightly later times in the main body of the disc, including near the mid-plane. These body modes have somewhat longer radial wavelengths. Early on they manifest themselves as fundamental breathing modes, but quickly become corrugation modes as symmetry about the mid-plane is broken. The corrugation modes are a prominent feature of the non-linear saturated state, leading to strong vertical oscillation of the disc mid-plane. In a viscous disc with aspect ratio H/r = 0.05, instability is found to operate when the viscosity parameter alpha < 4 x 10(-4). In three dimensions the instability generates a quasi-turbulent flow, and the associated Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha similar to 10(-3) by the end of the simulation. The evolution and saturation of the vertical shear instability in astrophysical disc models which include realistic treatments of the thermal physics has yet to be examined. Should it occur on either global or local scales, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties and observational appearance.

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

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Plutos vigorous geological activity. Composed of molecular nitrogen, methane, and carbon monoxide ices, but dominated by nitrogen ice, this layer is organized into cells or polygons, typically about 10 to 40 kilometres across, that resemble the surface manifestation of solid-state convection. Here we report, on the basis of available rheological measurements, that solid layers of nitrogen ice with a thickness in excess of about one kilometre should undergo convection for estimated present-day heat-flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-kilometre-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of nitrogen-ice viscosity implies that the ice layer convects in the so-called sluggish lid regime, a unique convective mode not previously definitively observed in the Solar System. Average surface horizontal velocities of a few centimetres a year imply surface transport or renewal times of about 500,000 years, well under the ten-million-year upper-limit crater retention age for Sputnik Planum. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help to explain the high albedos shown by some of these bodies.

The formation of Charon’s red poles from seasonally cold-trapped volatiles

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 tholin-like organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location on Charon implicates the temperature extremes that result from Charon’s high obliquity and long seasons in the production of this material. The escape of Pluto’s atmosphere provides a potential feedstock for a complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon’s winter pole was proposed as an explanation for the dark coloration 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 that show the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon and the supply and temporary cold-trapping of material escaping from Pluto, as well as the photolytic processing of this material into more complex and less volatile molecules while cold-trapped. The model results are consistent with the proposed mechanism for producing the observed colour pattern on Charon.

Vortex formation in protoplanetary discs induced by the vertical shear instability

This work used the DiRAC Complexity system, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk). The equipment is funded by BIS National E-Infrastructure capital grant ST/K000373/1 and STFC Operations grant ST/K0003259/1. DiRAC is part of the national E-Infrastructure.

Interacting vorticity waves as an instability mechanism for magnetohydrodynamic shear instabilities

The interacting vorticity wave formalism for shear flow instabilities is extended here to the magnetohydrodynamic (MHD) setting, to provide a mechanistic description for the stabilising and destabilising of shear instabilities by the presence of a background magnetic field. The interpretation relies on local vorticity anomalies inducing a non-local velocity field, resulting in action-at-a-distance. It is shown here that the waves supported by the system are able to propagate vorticity via the Lorentz force, and waves may interact; existence of instability then rests upon whether the choice of basic state allows for phase-locking and constructive interference of the vorticity waves via mutual interaction. To substantiate this claim, we solve the instability problem of two representative basic states, one where a background magnetic field stabilises an unstable flow and the other where the field destabilises a stable flow, and perform relevant analyses to show how this mechanism operates in MHD.

Linear analysis of the vertical shear instability: outstanding issues and improved solutions

This equipment is funded by BIS National E-Infrastructure capital grant ST/K000373/1 and STFC DiRAC Operations grant ST/K0003259/1. DiRAC is part of the National E-Infrastructure.

Great Expectations: Plans and Predictions for New Horizons Encounter With Kuiper Belt Object 2014 MU69 (“Ultima Thule”)

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 Ultimas composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

Corrigendum: Convection in a volatile nitrogen-ice-rich layer drives Pluto's geological vigour.

This corrects the article DOI: 10.1038/nature18289