Using hydrodynamic simulations and photoionization calculations, we demonstrate that quasar emission line spectra contain information on the driving mechanism of galaxy-scale outflows. Outflows driven by a hot shocked bubble are expected to exhibit LINER-like optical line ratios, while outflows driven by radiation pressure are expected to exhibit Seyfert-like line ratios. Driving by radiation pressure also has a distinct signature in the narrow UV lines, which is detected in an HST-COS spectrum of a nearby quasar hosting a large-scale wind.

We study disks and jets in various accretion states (SANE and MAD) using novel, GPU-accelerated general-relativistic magneto-hydrodynamic (GR-MHD) code which we developed, based on HARM. This code, written in CUDA-c and uses OpenMP to parallelize multi-GPU setups, allows high resolution simulations of accretion disks and the formation and structure of jets without the need of multi-node supercomputer infrastructure. A 2563 simulation is well within the reach of an Nvidia DGX-V100 server, with the computation being a factor about 100 times faster if only the CPU was used.

We use this code to examine several disk structures, wind and jet properties in the MAD and SANE states. In the MAD state, we find that the magnetic flux threading the horizon mostly depends on the spin of the BH. This implies that the jet structure and power are strong functions of the spin, with non-spinning BHs have the widest jets.