Collisionless Accretion onto Black Holes: Dynamics and Flares.

We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic and fluid simulations are identified. Magnetic reconnection is more efficient, and rapidly accelerates a nonthermal particle population, in our kinetic approach. In addition, the plasma in the kinetic simulations develops significant departures from thermal equilibrium, including pressure anisotropy that excites kinetic-scale instabilities, and a large field-aligned heat flux near the horizon that approaches the free-streaming value. We discuss the implications of our results for modeling event-horizon scale observations of Sgr A* and M87 by GRAVITY and the Event Horizon Telescope.

Black holes up close.

Recent developments have ushered in a new era in the field of black-hole astrophysics, providing a direct view of the remarkable environment near black-hole event horizons. These observations have enabled astronomers to confirm long-standing ideas on the physics of gas flowing into black holes with temperatures that are hundreds of times greater than at the centre of the Sun. At the same time, the observations have conclusively shown that light rays near a black hole experience large deflections that cause a dark shadow in the centre of the image, an effect predicted by Einstein's theory of general relativity. With further investment, this field is poised to deliver decades of advances in our understanding of gravity and black holes through stringent tests of general relativity, as well as insights into the role of black holes as the central engines powering a wide range of astronomical phenomena.

Shearing-box simulations of MRI-driven turbulence in weakly collisional accretion discs.

We present a systematic shearing-box investigation of MRI-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure anisotropy (Δp) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by considering different Reynolds numbers and magnetic-field configurations, including net vertical flux, net toroidal-vertical flux and zero net flux. Remarkably, we find that the level of turbulence and angular-momentum transport are not greatly affected by large anisotropic viscosities: the Maxwell and Reynolds stresses do not differ much from the MHD result. Angular-momentum transport in Braginskii MHD still depends strongly on isotropic dissipation, e.g., the isotropic magnetic Prandtl number, even when the anisotropic viscosity is orders of magnitude larger than the isotropic diffusivities. Braginskii viscosity nevertheless changes the flow structure, rearranging the turbulence to largely counter the parallel rate of strain from the background shear. We also show that the volume-averaged pressure anisotropy and anisotropic viscous transport decrease with increasing isotropic Reynolds number (Re); e.g., in simulations with net vertical field, the ratio of anisotropic to Maxwell stress (α A/α M) decreases from ~ 0.5 to ~ 0.1 as we move from Re ~ 103 to Re ~ 104, while 〈4πΔp/B 2〉 → 0. Anisotropic transport may thus become negligible at high Re. Anisotropic viscosity nevertheless becomes the dominant source of heating at large Re, accounting for ≳50% of the plasma heating. We conclude by briefly discussing the implications of our results for RIAFs onto black holes.

Where are the most ancient stars in the Milky Way?

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disc galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars (z form ≳ 5) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at z = 0 than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed in situ but are accreted during hierarchical assembly. These ex situ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed in situ in the solar neighbourhood, and in some simulations, even in the galactic centre. Secondly, old stars formed in situ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at z ≳ 2, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines away from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before z = 5. We predict most of these targets to have higher metallicity (-3 < [Fe/H] < -2) than the most extreme metal-poor stars.

Outbursts of luminous blue variable stars from variations in the helium opacity.

Luminous blue variables are massive, evolved stars that exhibit large variations in luminosity and size on timescales from months to years, with high associated rates of mass loss1-5. In addition to this on-going variability, these stars exhibit outburst phases, during which their size increases and as a result their effective temperature decreases, typically to about 9,000 kelvin3,6. Outbursts are believed to be caused by the radiation force on the cooler, more opaque, outer layers of the star balancing or even exceeding the force of gravity, although the exact mechanisms are unknown and cannot be determined using one-dimensional, spherically symmetric models of stars because such models cannot determine the physical processes that occur in this regime7. Here we report three-dimensional simulations of massive, radiation-dominated stars, which show that helium opacity has an important role in triggering outbursts and setting the observed effective temperature during outbursts of about 9,000 kelvin. It probably also triggers the episodic mass loss at rates of 10-7 to 10-5 solar masses per year. The peak in helium opacity is evident in our three-dimensional simulations only because the density and temperature of the stellar envelope (the outer part of the star near the photosphere) need to be determined self-consistently with convection, which cannot be done in one-dimensional models that assume spherical symmetry. The simulations reproduce observations of long-timescale variability, and predict that convection causes irregular oscillations in the radii of the stars and variations in brightness of 10-30 per cent on a typical timescale of a few days. The amplitudes of these short-timescale variations are predicted to be even larger for cooler stars (in the outburst phase). This short-timescale variability should be observable with high-cadence observations.

The origin of the diverse morphologies and kinematics of Milky Way-mass galaxies in the FIRE-2 simulations.

We use hydrodynamic cosmological zoom-in simulations from the Feedback in Realistic Environments project to explore the morphologies and kinematics of 15 Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90 per cent of their stellar mass within 2.5 kpc to well-ordered discs that reach ≳15 kpc. The gas in our galaxies always forms a thin, rotation-supported disc at z = 0, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disc-like orbits and with the radial extent of the stellar disc. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disc after z ~ 1 develop well-ordered stellar discs. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-z mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact discs in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in z = 0 MW-mass galaxies (even z = 0 bulge stars) form in a disc: ≳60-90 per cent of stars begin their lives rotationally supported.

Gas kinematics in FIRE simulated galaxies compared to spatially unresolved HI observations.

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotationally supported gas discs produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing V/σ and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of ~2000 low-redshift, HI-detected galaxies with Mstar = 107-11 M☉ and compare to the simulated galaxies. At Mstar ≳ 1010 M☉, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotationally supported discs. Both the observed and simulated line profiles become more Gaussian like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotational to dispersion support more strongly: at Mstar 108-10 M, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion supported. Most of the lower-mass-simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotational versus dispersion support.

Origin of the heavy elements in binary neutron-star mergers from a gravitational-wave event.

The cosmic origin of elements heavier than iron has long been uncertain. Theoretical modelling shows that the matter that is expelled in the violent merger of two neutron stars can assemble into heavy elements such as gold and platinum in a process known as rapid neutron capture (r-process) nucleosynthesis. The radioactive decay of isotopes of the heavy elements is predicted to power a distinctive thermal glow (a 'kilonova'). The discovery of an electromagnetic counterpart to the gravitational-wave source GW170817 represents the first opportunity to detect and scrutinize a sample of freshly synthesized r-process elements. Here we report models that predict the electromagnetic emission of kilonovae in detail and enable the mass, velocity and composition of ejecta to be derived from observations. We compare the models to the optical and infrared radiation associated with the GW170817 event to argue that the observed source is a kilonova. We infer the presence of two distinct components of ejecta, one composed primarily of light (atomic mass number less than 140) and one of heavy (atomic mass number greater than 140) r-process elements. The ejected mass and a merger rate inferred from GW170817 imply that such mergers are a dominant mode of r-process production in the Universe.

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma.

We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disk. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatiotemporally variable. Energy spectra suggest an Alfvén-wave cascade at large scales and a kinetic-Alfvén-wave cascade at small scales, with strong small-scale density fluctuations and weak nonaxisymmetric density waves. Ions undergo nonthermal particle acceleration, their distribution accurately described by a κ distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.

Numerical simulations of internal wave generation by convection in water.

Water's density maximum at 4°C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4°C, but stably stratified above 4°C. We present numerical simulations of a waterlike fluid near its density maximum in a two-dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. To isolate the physical mechanism exciting internal waves, we use the spectral code dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal-wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations overexcite high-frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find that the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal-wave generation by convection.

A possible relativistic jetted outburst from a massive black hole fed by a tidally disrupted star.

Gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, but most MBHs are considered dormant. Occasionally, a star passing too near an MBH is torn apart by gravitational forces, leading to a bright tidal disruption flare (TDF). Although the high-energy transient Sw 1644+57 initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that observations suggest a sudden accretion event onto a central MBH of mass about 10(6) to 10(7) solar masses. There is evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 to a temporary smaller-scale blazar.

Black hole accretion.

Black holes are most often detected by the radiation produced when they gravitationally pull in surrounding gas, in a process called accretion. The efficiency with which the hot gas radiates its thermal energy strongly influences the geometry and dynamics of the accretion flow. Both radiatively efficient thin disks and radiatively inefficient thick disks are observed. When the accreting gas gets close to the central black hole, the radiation it produces becomes sensitive to the spin of the hole and the presence of an event horizon. Analysis of the luminosities and spectra of accreting black holes has yielded tantalizing evidence for both rotating holes and event horizons. Numerical simulations imply that the relativistic jets often seen from accreting black holes may be powered in part by the spin of the hole.