Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant applica
tion using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20%$, the large-scale electromagnetic jet efficiency is of order $10%$, and the total radiative efficiency that reaches large distances remains low at only order $1%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs. disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $zsim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.
The rotational period of isolated pulsars increases over time due to the extraction of angular momentum by electromagnetic torques. These torques also change the obliquity angle $alpha$ between the magnetic and rotational axes. Although actual pulsar
magnetospheres are plasma-filled, the time evolution of $alpha$ has mostly been studied for vacuum pulsar magnetospheres. In this work, we self-consistently account for the plasma effects for the first time by analysing the results of time-dependent 3D force-free and magnetohydrodynamic simulations of pulsar magnetospheres. We show that if a neutron star is spherically symmetric and is embedded with a dipolar magnetic moment, the pulsar evolves so as to minimise its spin-down luminosity: both vacuum and plasma-filled pulsars evolve toward the aligned configuration ($alpha=0$). However, they approach the alignment in qualitatively different ways. Vacuum pulsars come into alignment exponentially fast, with $alpha propto exp(-t/tau)$ and $tau sim$ spindown timescale. In contrast, we find that plasma-filled pulsars align much more slowly, with $alpha propto (t/tau)^{-1/2}$. We argue that the slow time evolution of obliquity of plasma-filled pulsars can potentially resolve several observational puzzles, including the origin of normal pulsars with periods of $sim1$ second, the evidence that oblique pulsars come into alignment over a timescale of $sim 10^7$ years, and the observed deficit, relative to an isotropic obliquity distribution, of pulsars showing interpulse emission.
Black hole (BH) accretion flows and jets are qualitatively affected by the presence of ordered magnetic fields. We study fully three-dimensional global general relativistic magnetohydrodynamic (MHD) simulations of radially extended and thick (height
$H$ to cylindrical radius $R$ ratio of $|H/R|sim 0.2--1$) accretion flows around BHs with various dimensionless spins ($a/M$, with BH mass $M$) and with initially toroidally-dominated ($phi$-directed) and poloidally-dominated ($R-z$ directed) magnetic fields. Firstly, for toroidal field models and BHs with high enough $|a/M|$, coherent large-scale (i.e. $gg H$) dipolar poloidal magnetic flux patches emerge, thread the BH, and generate transient relativistic jets. Secondly, for poloidal field models, poloidal magnetic flux readily accretes through the disk from large radii and builds-up to a natural saturation point near the BH. For sufficiently high $|a/M|$ or low $|H/R|$ the polar magnetic field compresses the inflow into a geometrically thin highly non-axisymmetric magnetically choked accretion flow (MCAF) within which the standard linear magneto-rotational instability is suppressed. The condition of a highly-magnetized state over most of the horizon is optimal for the Blandford-Znajek mechanism that generates persistent relativistic jets with $gtrsim 100$% efficiency for $|a/M|gtrsim 0.9$. A magnetic Rayleigh-Taylor and Kelvin-Helmholtz unstable magnetospheric interface forms between the compressed inflow and bulging jet magnetosphere, which drives a new jet-disk quasi-periodic oscillation (JD-QPO) mechanism. The high-frequency QPO has spherical harmonic $|m|=1$ mode period of $tausim 70GM/c^3$ for $a/Msim 0.9$ with coherence quality factors $Qgtrsim 10$. [abridged]
