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.

GX~339$-$4 is a black hole X-ray binary that is a key focus of accretion studies since it goes into outburst roughly every two-to-three years. Tracking of its radio, IR and X-ray flux during multiple outbursts reveals tight broadband correlations. The radio emission originates in a compact, self-absorbed jet, however the origin of the X-ray emission is still debated: jet base or corona? We fit 20 quasi-simultaneous radio, IR, optical and X-ray observations of GX~339$-$4 covering three separate outbursts in 2005, 2007, 2010--2011, with a composite corona + jet model, where inverse Compton emission from both regions contributes to the X-ray emission. Using a recently-proposed identifier of the X-ray variability properties known as power-spectral hue, we attempt to explain both the spectral and evolving timing characteristics, with the model. We find the X-ray spectra are best fit by inverse Compton scattering in a dominant hot corona ($kT_{\rm e}\sim$ hundreds of keV). However, radio and IR-optical constraints imply a non-negligible contribution from inverse Compton scattering off hotter electrons ($kT_{\rm e} \ge 511$~keV) in the base of the jets, ranging from a few up to $\sim50$\% of the integrated 3--100~keV flux. We also find that the physical properties of the jet show interesting correlations with the shape of the broadband X-ray variability of the source, posing intriguing suggestions for the connection between the jet and corona.