We perform global three-dimensional MHD simulations of unstratified accretion disks in cataclysmic variables (CVs). By including mass inflow via an accretion stream, we are able to evolve the disk to a steady state. We investigate the relative importance of spiral shocks and the magnetorotational instability (MRI) in driving angular momentum transport and how each depend on the geometry and strength of seed magnetic field and the Mach number of the disk (where Mach number is ratio of the azimuthal velocity and the sound speed of gas). We use a locally isothermal equation of state and adopt temperature profiles that are consistent with CV disk observations. Our results indicate that the relative importance of spiral shocks and MRI in driving angular momentum transport is controlled by the gas Mach number and the seed magnetic field strength. MRI and spiral shocks provide comparable efficiency of angular momentum transport when the disk Mach number is around 10 and the seed magnetic field has plasma $\beta=400$ (where $\beta$ is ratio of gas pressure and magnetic pressure). The MRI dominates whenever the seed field strength, or the disk Mach number, is increased. Among all of our simulations, the effective viscosity parameter $\alpha_{eff} \sim 0.016-0.1$ after MRI saturates and the disk reaches steady state. Larger values of $\alpha_{eff}$ are favored when the seed magnetic field has vertical components or the flow has stronger magnetization ($1/\beta$). Our models all indicate that the role of MRI in driving angular momentum transport thus mass accretion in CV disks is indispensable, especially in cool disks with weak spiral shocks.

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two dimensional hydrodynamical simulations with radiative cooling to study CPDs, and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock driven accretion is, on the other hand, unsteady on a timescale of months/years due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet's Hill radius all the way to the planet surface, and the disk $\alpha$-coefficient characterizing angular momentum transport due to spiral shocks is $\sim$0.001-0.02. The disk surface density ranges from 10 to 1000 g cm$^{-2}$ in our simulations, which is at least 3 orders of magnitude smaller than the "minimum mass sub-nebula" model used to study satellite formation; instead it is more consistent with the "gas-starved" satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs.

Supermassive black hole binaries (SMBHs) are expected to result from galaxy mergers, and thus are natural byproducts (and probes) of hierarchical structure formation in the Universe. They are also the primary expected source of low-frequency gravitational wave emission. We search for binary BHs using time-variable velocity shifts in broad Mg II emission lines of quasars with multi-epoch observations. First, we inspect velocity shifts of the binary SMBH candidates identified in Ju et al. (2013), using SDSS spectra with an additional epoch of data that lengthens the typical baseline to ~10 yr. We find variations in the line-of-sight velocity shifts over 10 years that are comparable to the shifts observed over 1-2 years, ruling out the binary model for the bulk of our candidates. We then analyze 1438 objects with 8 yr median time baselines, from which we would expect to see velocity shifts >1000 km/s from sub-pc binaries. We find only one object with an outlying velocity of 448 km/s, indicating, based on our modeling, that ~< 1 per cent (the value varies with different assumptions) of SMBHs that are active as quasars reside in binaries with ~0.1 pc separations. Binaries either sweep through these small separations rapidly or stall at larger radii.

We present results from the first global 3D MHD simulations of accretion disks in Cataclysmic Variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. 3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate alpha_{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. 4) Using the mass accretion rate measured in our simulations to directly measure alpha defined in standard thin-disk theory, we find 0.02 < alpha_{eff} < 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae (DNe). In this regime the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime, and therefore future models with more realistic thermodynamics and non-ideal MHD are warranted.

Supermassive black hole (SMBH) binaries are expected in a Lambda CDM cosmology given that most (if not all) massive galaxies contain a massive black hole at their center. So far, however, direct evidence for such binaries has been elusive. We use cross-correlation to search for temporal velocity shifts in the MgII broad emission lines of 0.36 < z < 2 quasars with multiple observations in the Sloan Digital Sky Survey. For ~ 10^9 Msun BHs in SMBH binaries, we are sensitive to velocity drifts for binary separations of ~ 0.1 pc with orbital periods of ~100 years. We find seven candidate sub-pc--scale binaries with velocity shifts > 3.4 sigma ~ 280 km/s, where sigma is our systematic error. Comparing the detectability of SMBH binaries with the number of candidates (N < 7), we can rule out that most 10^9 Msun BHs exist in ~ 0.03-0.2 pc scale binaries, in a scenario where binaries stall at sub-pc scales for a Hubble time. We further constrain that < one-third of quasars host SMBH binaries after considering gas-assisted sub-pc evolution of SMBH binaries, although this result is very sensitive to the assumed size of the broad line region. We estimate the detectability of SMBH binaries with ongoing or next-generation surveys (e.g., BOSS, Subaru Prime Focus Spectrograph), taking into account the evolution of the sub-parsec binary in circumbinary gas disks. These future observations will provide longer time baselines for searches similar to ours and may in turn constrain the evolutionary scenarios of SMBH binaries.