Black hole accretion is one of nature's most efficient energy extraction processes. When gas falls in, a significant fraction of its gravitational binding energy is either converted into radiation or flows outwards in the form of black hole-driven jets and disk-driven winds. Recently, the Event Horizon Telescope (EHT), an Earth-size sub-millimetre radio interferometer, captured the first images of M87's black hole (or M87*). These images were analysed and interpreted using general-relativistic magnetohydrodynamics (GRMHD) models of accretion disks with rotation axes aligned with the black hole spin axis. However, since infalling gas is often insensitive to the black hole spin direction, misalignment between accretion disk and black hole spin may be a common occurrence in nature. In this work, we use the general-relativistic radiative transfer (GRRT) code \texttt{BHOSS} to calculate the first synthetic radio images of (highly) tilted disk/jet models generated by our GPU-accelerated GRMHD code \texttt{HAMR}. While the tilt does not have a noticeable effect on the system dynamics beyond a few tens of gravitational radii from the black hole, the warping of the disk and jet can imprint observable signatures in EHT images on smaller scales. Comparing the images from our GRMHD models to the 43 GHz and 230 GHz EHT images of M87, we find that M87 may feature a tilted disk/jet system. Further, tilted disks and jets display significant time variability in the 230 GHz flux that can be further tested by longer-duration EHT observations of M87.

General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black-hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multi-faceted optimizations that, collectively, accelerate computation by 2-5 orders of magnitude for a wide range of applications. Firstly, it involves a novel implementation of a spherical-polar grid with 3D adaptive mesh refinement that operates in each of the 3 dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-res computational performance. Secondly, we demonstrate that local adaptive time-stepping (LAT) on a logarithmic spherical-polar grid accelerates computation by a factor of $\lesssim10$ compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of $\sim10^9$ zone-cycles-per-second-per-node on 5,400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We demonstrate its computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of $13$,$440\times4$,$608\times8$,$092$ cells, and a total of $\lesssim22$ billion cells and $\sim0.65\times10^8$ timesteps, it is among the largest astrophysical simulations ever performed. We find that frame-dragging by the black hole tears up the disk into two independently precessing sub-disks. The innermost sub-disk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

Accreting black holes produce collimated outflows, or jets, that traverse many orders of magnitude in distance, accelerate to relativistic velocities, and collimate into tight opening angles. Of these, perhaps the least understood is jet collimation due to the interaction with the ambient medium. In order to investigate this interaction, we carried out axisymmetric general relativistic magnetohydrodynamic simulations of jets produced by a large accretion disc, spanning over 5 orders of magnitude in time and distance, at an unprecedented resolution. Supported by such a disc, the jet attains a parabolic shape, similar to the M87 galaxy jet, and the product of the Lorentz factor and the jet half-opening angle, $\gamma\theta\ll 1$, similar to values found from very long baseline interferometry (VLBI) observations of active galactic nuclei (AGN) jets; this suggests extended discs in AGN. We find that the interaction between the jet and the ambient medium leads to the development of pinch instabilities, which produce significant radial and lateral variability across the jet by converting magnetic and kinetic energy into heat. Thus pinched regions in the jet can be detectable as radiating hotspots and may provide an ideal site for particle acceleration. Pinching also causes gas from the ambient medium to become squeezed between magnetic field lines in the jet, leading to enhanced mass-loading of the jet and potentially contributing to the spine-sheath structure observed in AGN outflows.

Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) tend to be surrounded by geometrically thin, radiatively cooled accretion discs. According to both theory and observations, these are -- in many cases -- highly misaligned with the black hole spin axis. In this work we present the first general relativistic magnetohydrodynamic simulations of very thin ($h/r \sim 0.015-0.05$) accretion discs around rapidly spinning ($a \sim 0.9$) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with $h/r \lesssim 0.03$ align with the black hole equator, though at smaller radii than predicted by theoretical work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame-dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that at all tilt values the system produces a pair of relativistic jets. At small distances the jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.

In the course of its evolution, a black hole (BH) accretes gas from a wide range of directions. Given a random accretion event, the typical angular momentum of an accretion disc would be tilted by $\sim$60$^\circ$ relative to the BH spin. Misalignment causes the disc to precess at a rate that increases with BH spin and tilt angle. We present the first general-relativistic magnetohydrodynamic (GRMHD) simulations spanning a full precession period of highly tilted (60$^\circ$), moderately thin ($h/r=0.1$) accretion discs around a rapidly spinning ($a\simeq0.9$) BH. While the disc and jets precess in phase, we find that the corona, sandwiched between the two, lags behind by $\gtrsim 10^{\circ}$. For spectral models of BH accretion, the implication is that hard non-thermal (corona) emission lags behind the softer (disc) emission, thus potentially explaining some properties of the hard energy lags seen in Type-C low frequency quasi-periodic oscillations in X-Ray binaries. While strong jets are unaffected by this disc-corona lag, weak jets stall when encountering the lagging corona at distances $r \sim 100$ black hole radii. This interaction may quench large-scale jet formation.

Prevalent around luminous accreting black holes, thin discs are challenging to resolve in numerical simulations. When the disc and black hole angular momentum vectors are misaligned, the challenge becomes extreme, requiring adaptive meshes to follow the disc proper as it moves through the computational grid. With our new high-performance general relativistic magnetohydrodynamic (GRMHD) code H-AMR we have simulated the thinnest accretion disc to date, of aspect ratio $H/R\approx0.03\approx1.7^\circ$ , around a rapidly spinning ($a=0.9375$) black hole, using a cooling function. Initially tilted at $10^\circ$, the disc warps inside $\sim 5r_g$ into alignment with the black hole, where $r_g$ is the gravitational radius. This is the first demonstration of Bardeen-Petterson alignment in MHD with viscosity self-consistently generated by magnetized turbulence. The disc develops a low-density high-viscosity ($\alpha_{\rm eff}\sim1.0$) magnetic-pressure--dominated inner region at $r\lesssim25r_g$ that rapidly empties itself into the black hole. This inner region may in reality, due to thermal decoupling of ions and electrons, evaporate into a radiatively inefficient accretion flow if, as we propose, the cooling time exceeds the accretion time set by the order unity effective viscosity. We furthermore find the unexpected result that even our very thin disc can sustain large-scale vertical magnetic flux on the black hole, which launches powerful relativistic jets that carry $20{-}50\%$ of the accretion power along the angular momentum vector of the outer tilted disc, providing a potential explanation for the origin of jets in radio-loud quasars.

Accreting black holes launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code HAMR to study global 3D black hole accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that accretion disc turbulence in a radially-extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the black hole till it becomes dynamically-important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10,000. The jets do not show significant kink or pinch instabilities, accelerate to $\gamma \sim 10$ over 3 decades in distance, and follow a collimation profile similar to the observed M87 jet.

Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of tilted thick accretion discs around rapidly spinning BHs, at the highest resolution to date. We explored the limit where disc thermal pressure dominates magnetic pressure, and showed for the first time that, for different magnetic field strengths on the BH, these flows launch magnetized relativistic jets propagating along the rotation axis of the tilted disc (rather than of the BH). If strong large-scale magnetic flux reaches the BH, it bends the inner few gravitational radii of the disc and jets into partial alignment with the BH spin. On longer time scales, the simulated disc-jet system as a whole undergoes Lense-Thirring precession and approaches alignment, demonstrating for the first time that jets can be used as probes of disc precession. When the disc turbulence is well-resolved, our isolated discs spread out, causing both the alignment and precession to slow down.

Recent Herschel/SPIRE maps of the Small and Large Magellanic Clouds (SMC, LMC) exhibit in each thousands of clouds. Observed at 250 microns, they must be cold, T ~ 15 K, hence the name "Herschel cold clouds" (HCCs). From the observed rotational velocity profile and the assumption of spherical symmetry, the Galactic mass density is modeled in a form close to that of an isothermal sphere. If the HCCs constitute a certain fraction of it, their angular size distribution has a specified shape. A fit to the data deduced from the SMC/LMC maps supports this and yields for their radius 2.5 pc, with a small change when allowing for a spread in HCC radii. There are so many HCCs that they will make up all the missing Halo mass density if there is spherical symmetry and their average mass is of order 15,000 Mo. This compares well with the Jeans mass of circa 40,000 Mo and puts forward that the HCCs are in fact Jeans clusters, constituting all the Galactic dark matter and much of its missing baryons, a conclusion deduced before from a different field of the sky (Nieuwenhuizen, Schild and Gibson 2011). A preliminary analysis of the intensities yields that the Jeans clusters themselves may consist of some billion MACHOs of a few dozen Earth masses. With a size of dozens of solar radii, they would mostly obscure stars in the LMC, SMC and towards the Galactic center, and may thus have been overlooked in microlensing.