In this paper we continue the first ever study of magnetized mini-disks coupled to circumbinary accretion in a supermassive binary black hole (SMBBH) approaching merger reported in Bowen et al. 2018. We extend this simulation from 3 to 12 binary orbital periods. We find that relativistic SMBBH accretion acts as a resonant cavity, where quasi-periodic oscillations tied to the the frequency at which the black hole's orbital phase matches a non-linear $m=1$ density feature, or ``lump'', in the circumbinary accretion disk permeate the system. The rate of mass accretion onto each of the mini-disks around the black holes is modulated at the beat frequency between the binary frequency and the lump's mean orbital frequency, i.e., $\Omega_{\rm beat} = \Omega_{\rm bin} - \bar{\Omega}_{\rm lump}$, while the total mass accretion rate of this equal-mass binary is modulated at two different frequencies, $\gtrsim \bar{\Omega}_{\rm lump}$ and $\approx 2 \Omega_{\rm beat}$. The instantaneous rotation rate of the lump itself is also modulated at two frequencies close to the modulation frequencies of the total accretion rate, $\bar{\Omega}_{\rm lump}$ and $2 \Omega_{\rm beat}$. Because of the compact nature of the mini-disks in SMBBHs approaching merger, the inflow times within the mini-disks are comparable to the period on which their mass-supply varies, so that their masses---and the accretion rates they supply to their black holes---are strongly modulated at the same frequency. In essence, the azimuthal symmetry of the circumbinary disk is broken by the dynamics of orbits near a binary, and this $m=1$ asymmetry then drives quasi-periodic variation throughout the system, including both accretion and disk-feeding. In SMBBHs approaching merger, such time variability could introduce distinctive, increasingly rapid, fluctuations in their electromagnetic emission.

The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a strong magnetic field is introduced to numerical simulations of thin disks, large deviations in efficiency are observed, in part due to mass and energy carried by jets and winds launched by the disk or BH spin. The total efficiency of accretion can be significantly enhanced beyond that predicted by NT but it has remained unclear how the radiative component is affected. In order to study the effect of a dynamically relevant large-scale magnetic field on radiatively efficient accretion, we have performed numerical 3D general relativistic - radiative - magnetohydroynamic (GRRMHD) simulations of a disk with scale height to radius ratio of $H/R~0.1$ around a moderately spinning BH (a=0.5) using the code HARMRAD. Our simulations are fully global and allow us to measure the jet, wind, and radiative properties of a magnetically arrested disk (MAD) that is kept thin via self-consistent transport of energy by radiation using the M1 closure scheme. Our fiducial disk is MAD out to a radius of ~16R_g and the majority of the total ~13\% efficiency of the accretion flow is carried by a magnetically driven wind. We find that the radiative efficiency is slightly suppressed compared to NT, contrary to prior MAD GRMHD simulations with an ad hoc cooling function, but it is unclear how much of the radiation and thermal energy trapped in the outflows could ultimately escape.

In accretion disks with large-scale ordered magnetic fields, the magnetorotational instability (MRI) is marginally suppressed, so other processes may drive angular momentum transport leading to accretion. Accretion could then be driven by large-scale magnetic fields via magnetic braking, but large-scale magnetic flux can build-up onto the black hole and within the disk leading to a magnetically-arrested disk (MAD). Such a MAD state is unstable to the magnetic Rayleigh-Taylor (RT) instability, which itself leads to vigorous turbulence and the emergence of low-density highly-magnetized bubbles. This instability was studied in a thin (ratio of half-height H to radius R, $H/R \approx 0.1$) MAD simulation, where it has a more dramatic effect on the dynamics of the disk than for thicker disks. We find that the low-density bubbles created by the magnetic RT instability decrease the stress (leading to angular momentum transport) in the disk rather than increasing magnetic torques. Indeed, we find that the dominant component of the stress is due to turbulent magnetic fields, despite the suppression of the axisymmetric MRI and the dominant presence of large-scale magnetic fields. This suggests that the magnetic RT instability plays a significant role in driving angular momentum transport in MADs.

The radiative and jet efficiencies of thin magnetized accretion disks around black holes (BHs) are affected by BH spin and the presence of a magnetic field that, when strong, could lead to large deviations from Novikov-Thorne (NT) thin disk theory. To seek the maximum deviations, we perform general relativistic magnetohydrodynamic (GRMHD) simulations of radiatively efficient thin (half-height $H$ to radius $R$ of $H/R\approx 0.10$) disks around moderately rotating BHs with $a/M=0.5$. First, our simulations, each evolved for more than $70,000r_g/c$ (gravitational radius $r_g$ and speed of light $c$), show that large-scale magnetic field readily accretes inward even through our thin disk and builds-up to the magnetically-arrested disk (MAD) state. Second, our simulations of thin MADs show the disk achieves a radiative efficiency of $\eta_{\rm r}\approx 15\%$ (after estimating photon capture), which is about twice the NT value of $\eta_{\rm r}\sim 8\%$ for $a/M=0.5$ and gives the same luminosity as a NT disk with $a/M\approx 0.9$. Compared to prior simulations with $\lesssim 10\%$ deviations, our result of an $\approx 80\%$ deviation sets a new benchmark. Building on prior work, we are now able to complete an important scaling law which suggest that observed jet quenching in the high-soft state in BH X-ray binaries is consistent with an ever-present MAD state with a weak yet sustained jet.

The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) simulation of a spinning BH (spin $a/M=0.8$) accreting at $\sim 50$ times Eddington shows a total efficiency $\sim 50\%$ when time-averaged and total efficiency $\gtrsim 100\%$ in moments. Magnetic compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation escapes via advection by a magnetized wind and via transport through a low-density channel created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial loading of field lines by optically thick radiation, leading to BZ efficiency $\sim 40\%$ on the horizon and BZ efficiency $\sim 5\%$ by $r\sim 400r_g$ (gravitational radii) via absorption by the wind. Importantly, radiation escapes at $r\sim 400r_g$ with efficiency $\eta\approx 15\%$ (luminosity $L\sim 50L_{\rm Edd}$), similar to $\eta\approx 12\%$ for a Novikov-Thorne thin disk and beyond $\eta\lesssim 1\%$ seen in prior GRRMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion.

The role played by magnetic fields in the intracluster medium (ICM) of galaxy clusters is complex. The weakly collisional nature of the ICM leads to thermal conduction that is channelled along field lines. This anisotropic heat conduction profoundly changes the stability of the ICM atmosphere, with convective stabilities being driven by temperature gradients of either sign. Here, we employ the Athena magnetohydrodynamic code to investigate the local non-linear behavior of the heat-flux driven buoyancy instability (HBI), relevant in the cores of cooling-core clusters where the temperature increases with radius. We study a grid of 2-d simulations that span a large range of initial magnetic field strengths and numerical resolutions. For very weak initial fields, we recover the previously known result that the HBI wraps the field in the horizontal direction thereby shutting off the heat flux. However, we find that simulations which begin with intermediate initial field strengths have a qualitatively different behavior, forming HBI-stable filaments that resist field-line wrapping and enable sustained vertical conductive heat flux at a level of 10--25% of the Spitzer value. While astrophysical conclusions regarding the role of conduction in cooling cores require detailed global models, our local study proves that systems dominated by HBI do not necessarily quench the conductive heat flux.

Chandra X-ray observations of the high redshift (z =1.532) radio-loud quasar 3C270.1 in 2008 February show the nucleus to have a power-law spectrum, Gamma = 1.66 +/- 0.08, typical of a radio-loud quasar, and a marginally-detected Fe Kalpha emission line. The data also reveal extended X-ray emission, about half of which is associated with the radio emission from this source. The southern emission is co-spatial with the radio lobe and peaks at the position of the double radio hotspot. Modeling this hotspot including Spitzer upper limits rules out synchrotron emission from a single power-law population of electrons, favoring inverse-Compton emission with a field of ~11nT, roughly a third of the equipartition value. The northern emission is concentrated close to the location of a 40 deg. bend where the radio jet is presumed to encounter external material. It can be explained by inverse Compton emission involving Cosmic Microwave Background photons with a field of ~3nT, roughly a factor of nine below the equipartition value. The remaining, more diffuse X-ray emission is harder (HR=-0.09 +/- 0.22). With only 22.8+/-5.6 counts, the spectral form cannot be constrained. Assuming thermal emission with a temperature of 4 keV yields an estimate for the luminosity of 1.8E44 erg/s, consistent with the luminosity-temperature relation of lower-redshift clusters. However deeper Chandra X-ray observations are required to delineate the spatial distribution, and better constrain the spectrum of the diffuse emission to verify that we have detected X-ray emission from a high-redshift cluster.