Source author record

Julian H. Krolik

Julian H. Krolik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph gr-qc astro-ph.GA astro-ph.CO astro-ph.IM astro-ph.SR physics.comp-ph physics.plasm-ph

Catalog footprint

What is connected

56works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

preprint2022arXiv

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by an $m=1$ overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases, we find most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are $8$ times stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.

preprint2022arXiv

What powers the radio emission in TDE AT2019dsg: a long-lived jet or the disruption itself?

The tidal disruption event AT2019dsg was observed from radio to X-rays and was possibly accompanied by a high-energy neutrino. Previous interpretations have focused on continued injection by a central engine as the source of energy for radio emission. We show that continuous energy injection is unnecessary; the radio data can be explained by a single ejection of plasma that supplies all the energy needed. To support this assertion, we analyze the synchrotron self-absorbed spectra in terms of the equipartition model. Similar to previous analyses, we find that the energy in the radio-emitting region increases approximately $\propto t^{0.7}$ and the lengthscale of this region grows $\propto t$ at a rate $\simeq0.06c$. This event resembles the earliest stage of a supernova remnant: because the ejected mass is much greater than the shocked external mass, its velocity remains unchanged, while the energy in shocked gas grows with time. The radio-emitting material gains energy from the outflow, not continuing energy injection by the central object. Although energy injection from an accreting BH cannot be completely excluded, the energy injection rate is very different from the fallback luminosity, and maintaining constant outflow velocity requires fine-tuning demanding further physical explanation. If the neutrino association is real, the energy injection needed is much greater than for the radio emission, suggesting that the detected neutrino did not arise from the radio-emitting region.

preprint2021arXiv

High-energy emission from tidal disruption events in active galactic nuclei

Tidal disruption events (TDEs) taking place in active galactic nuclei (AGNs) are different from ordinary TDEs. In these events, the returning tidal debris stream drills through the pre-existing AGN accretion disk near the stream pericenter, destroying the inner disk in the process, and then intersects with the disk a second time at radii ranging from a few to hundreds of times the pericenter distance. The debris dynamics of such TDEs, and hence their appearance, are distinct from ordinary TDEs. Here we explore the observational signatures of this "second impact" of the stream with the disk. Strong shocks form as the dilute stream is stopped by the denser disk. Compton cooling of the shocked material produces hard X-rays, even soft gamma-rays, with most of the energy emitted between ~10 keV and 1 MeV. The luminosity follows the mass-return rate, peaking between ~$10^{42}$ and $10^{44}$ erg/s. The X-ray hardness and the smoothness of the light curve provide possible means for distinguishing the second impact from ordinary AGN flares, which exhibit softer spectra and more irregular light curves.

preprint2020arXiv

Future Simulations of Tidal Disruption Events

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

preprint2020arXiv

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprised of an arbitrary number of moving, local meshes, or patches, which are free to employ their own resolution, coordinate system/topology, physics equations, reference frame, and in our new approach, numerical method. Each local patch exchanges boundary data with a single global patch on which all other patches reside through a client-router-server parallelization model. In generalizing Patchwork to be compatible with arbitrary order time integration, PatchworkMHD and PatchworkWave have significantly improved the interpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in the original Patchwork code. Furthermore, we extend Patchwork to be multimethod by allowing multiple state vectors to be updated simultaneously, with each state vector providing its own interpatch interpolation and transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partial differential equations. We demonstrate our changes through the implementation of a scalar wave toy-model that is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

preprint2020arXiv

Radiative Emission Mechanisms of Tidal Disruption Events

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

preprint2020arXiv

Tidal disruption events in active galactic nuclei

A fraction of tidal disruption events (TDEs) occur in active galactic nuclei (AGNs) whose black holes possess accretion disks; these TDEs can be confused with common AGN flares. The disruption itself is unaffected by the disk, but the evolution of the bound debris stream is modified by its collision with the disk when it returns to pericenter. The outcome of the collision is largely determined by the ratio of the stream mass current to the azimuthal mass current of the disk rotating underneath the stream footprint, which in turn depends on the mass and luminosity of the AGN. To characterize TDEs in AGNs, we simulated a suite of stream--disk collisions with various mass current ratios. The collision excites shocks in the disk, leading to inflow and energy dissipation orders of magnitude above Eddington; however, much of the radiation is trapped in the inflow and advected into the black hole, so the actual bolometric luminosity may be closer to Eddington. The emergent spectrum may not be thermal, TDE-like, or AGN-like. The rapid inflow causes the disk interior to the impact point to be depleted within a fraction of the mass return time. If the stream is heavy enough to penetrate the disk, part of the outgoing material eventually hits the disk again, dissipating its kinetic energy in the second collision; another part becomes unbound, emitting synchrotron radiation as it shocks with surrounding gas.

preprint2018arXiv

Magnetorotational instability in eccentric disks

Eccentric disks arise in such astrophysical contexts as tidal disruption events, but it is unknown whether the magnetorotational instability (MRI), which powers accretion in circular disks, operates in eccentric disks as well. We examine the linear evolution of unstratified, incompressible MRI in an eccentric disk orbiting a point mass. We consider vertical modes of wavenumber $k$ on a background flow with uniform eccentricity $e$ and vertical Alfvén speed $v_\mathrm A$ along an orbit with mean motion $n$. We find two mode families, one with dominant magnetic components, the other with dominant velocity components; the former is unstable at $(1-e)^3f^2\lesssim3$, where $f\equiv kv_\mathrm A/n$, the latter at $e\gtrsim0.8$. For $f^2\lesssim3$, MRI behaves much like in circular disks, but the growth per orbit declines slowly with increasing $e$; for $f^2\gtrsim3$, modes grow by parametric amplification, which is resonant for $0<e\ll1$. MRI growth and the attendant angular momentum and energy transport happen chiefly near pericenter, where orbital shear dominates magnetic tension.

preprint2016arXiv

Discovery of transient infrared emission from dust heated by stellar tidal disruption flares

Stars that pass within the Roche radius of a supermassive black hole will be tidally disrupted, yielding a sudden injection of gas close to the black hole horizon which produces an electromagnetic flare. A few dozen of these flares have been discovered in recent years, but current observations provide poor constraints on the bolometric luminosity and total accreted mass of these events. Using images from the Wide-field Infrared Survey Explorer (WISE), we have discovered transient 3.4 micron emission from several previously known tidal disruption flares. The observations can be explained by dust heated to its sublimation temperature due to the intense radiation of the tidal flare. From the break in the infrared light curve we infer that this hot dust is located ~0.1 pc from the supermassive black hole. Since the dust has been heated by absorbing UV and (potentially) soft X-ray photons of the flare, the reprocessing light curve yields an estimate of the bolometric flare luminosity. For the flare PTF-09ge, we infer that the most likely value of the luminosity integrated over frequencies at which dust can absorb photons is $8\times 10^{44}$ erg/s, with a factor of 3 uncertainty due to the unknown temperature of the dust. This bolometric luminosity is a factor ~10 larger than the observed black body luminosity. Our work is the first to probe dust in the nuclei of non-active galaxies on sub-parsec scales. The observed infrared luminosity implies a covering factor ~1% for the nuclear dust in the host galaxies.

preprint2016arXiv

Fe K$α$ Profiles from Simulations of Accreting Black Holes

We present first results from a new technique for the prediction of Fe K$α$ profiles directly from general relativistic magnetohydrodynamic (GRMHD) simulations. Data from a GRMHD simulation are processed by a Monte Carlo global radiation transport code, which determines the X-ray flux irradiating the disk surface and the coronal electron temperature self-consistently. With that irradiating flux and the disk's density structure drawn from the simulation, we determine the reprocessed Fe K$α$ emission from photoionization equilibrium and solution of the radiation transfer equation. We produce maps of the surface brightness of Fe K$α$ emission over the disk surface, which---for our example of a $10 M_\odot$, Schwarzschild black hole accreting at $1\%$ the Eddington value---rises steeply one gravitational radius outside the radius of the innermost stable circular orbit and then falls $\propto r^{-2}$ at larger radii. We explain these features of the Fe K$α$ radial surface brightness profile as consequences of the disk's ionization structure and an extended coronal geometry, respectively. We also present the corresponding Fe K$α$ line profiles as would be seen by distant observers at several inclinations. Both the shapes of the line profiles and the equivalent widths of our predicted K$α$ lines are qualitatively similar to those typically observed from accreting black holes. Most importantly, this work represents a direct link between theory and observation: in a fully self-consistent way, we produce observable results---iron fluorescence line profiles---from the theory of black hole accretion with almost no phenomenological assumptions.

preprint2016arXiv

How bright are the gaps in circumbinary disk systems?

When a circumbinary disk surrounds a binary whose secondary's mass is at least $\sim 10^{-2}\times$ the primary's mass, a nearly empty cavity with radius a few times the binary separation is carved out of the disk. Narrow streams of material pass from the inner edge of the circumbinary disk into the domain of the binary itself, where they eventually join onto the small disks orbiting the members of the binary. Using data from 3-d MHD simulations of this process, we determine the luminosity of these streams; it is mostly due to weak laminar shocks, and is in general only a few percent of the luminosity of adjacent regions of either the circumbinary disk or the "mini-disks". This luminosity therefore hardly affects the deficit in the thermal continuum predicted on the basis of a perfectly dark gap region.

preprint2016arXiv

Radiation-driven outflows from and radiative support in dusty tori of active galactic nuclei

Substantial evidence points to dusty, geometrically thick tori obscuring the central engines of active galactic nuclei (AGNs), but so far no mechanism satisfactorily explains why cool dust in the torus remains in a puffy geometry. Near-Eddington infrared (IR) and ultraviolet (UV) luminosities coupled with high dust opacities at these frequencies suggest that radiation pressure on dust can play a significant role in shaping the torus. To explore the possible effects of radiation pressure, we perform three-dimensional radiative hydrodynamics simulations of an initially smooth torus. Our code solves the hydrodynamics equations, the time-dependent multi-angle group IR radiative transfer (RT) equation, and the time-independent UV RT equation. We find a highly dynamic situation. IR radiation is anisotropic, leaving primarily through the central hole. The torus inner surface exhibits a break in axisymmetry under the influence of radiation and differential rotation; clumping follows. In addition, UV radiation pressure on dust launches a strong wind along the inner surface; when scaled to realistic AGN parameters, this outflow travels at $\sim 5000 (M/10^7 M_\odot)^{1/4} [L_\mathrm{UV}/(0.1 L_\mathrm E)]^{1/4} \mathrm{km}\,\mathrm s^{-1}$ and carries $\sim 0.1 (M/10^7 M_\odot)^{3/4} [L_\mathrm{UV}/(0.1 L_\mathrm E)]^{3/4} M_\odot\,\mathrm{yr}^{-1}$, where $M$, $L_\mathrm{UV}$, and $L_\mathrm E$ are the mass, UV luminosity, and Eddington luminosity of the central object respectively.

preprint2015arXiv

A Steady-State Alignment Front in an Accretion Disk Subjected to Lense-Thirring Torques

Using only physical mechanisms, i.e., 3D MHD with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian gravitomagnetism. In this simulation, we have shown how, in the presence of MHD turbulence, a time-steady transition can be achieved between an inner disk region aligned with the equatorial plane of the central mass's spin and an outer region orbiting in a different plane. The position of the equilibrium orientation transition is determined by a balance between gravitomagnetic torque and warp-induced inward mixing of misaligned angular momentum from the outer disk. If the mixing is interpreted in terms of diffusive transport, the implied diffusion coefficient is ~(0.6--0.8)c_s^2/Omega for sound speed c_s and orbital frequency Omega. This calibration permits estimation of the orientation transition's equilibrium location given the central mass, its spin parameter, and the disk's surface density and scaleheight profiles. However, the alignment front overshoots before settling into an equilibrium, signaling that a diffusive model does not fully represent the time-dependent properties of alignment fronts under these conditions. Because the precessional torque on the disk at the alignment front is always comparable to the rate at which misaligned angular momentum is brought inward to the front by warp-driven radial motions, no break forms between the inner and outer portions of the disk in our simulation. Our results also raise questions about the applicability to MHD warped disks of the traditional distinction between "bending wave" and "diffusive" regimes.

preprint2015arXiv

Disk Emission from Magneto-hydrodynamic Simulations of Spinning Black Holes

We present the results of a new series of global 3D relativistic magneto-hydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of $H/R\sim 0.05$ and spin parameters $a/M=0, 0.5, 0.9$, and $0.99$. Using the ray-tracing code Pandurata, we generate broad-band thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Lastly, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well-known to the continuum fitting method of measuring black hole spin.

preprint2015arXiv

Evolution of a Binary Black Hole with a Retrograde Circumbinary Accretion Disk

We consider the evolution of a supermassive black hole binary (SMBHB) surrounded by a retrograde accretion disk. Assuming the disk is exactly in the binary plane and transfers energy and angular momentum to the binary via direct gas accretion, we calculate the time evolution of the binary's semi-major axis $a$ and eccentricity $e$. Because the gas is predominantly transferred when the binary is at apocenter, we find the eccentricity grows rapidly while maintaining constant $a(1+e)$. After accreting only a fraction of the secondary's mass, the eccentricity grows to nearly unity; from then on, gravitational wave emission dominates the evolution, preserving constant $a(1-e)$. The high-eccentricity waveforms redistribute the peak gravitational wave power from the nHz to $μ$Hz bands, substantially affecting the signal that might be detected with pulsar timing arrays. We also estimate the torque coupling binaries of arbitrary eccentricity with obliquely aligned circumbinary disks. If the outer edge of the disk is not an extremely large multiple of the binary separation, retrograde accretion can drive the binary into the gravitational wave-dominated state before these torques align the binary with the angular momentum of the mass supply.

preprint2015arXiv

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

preprint2015arXiv

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

preprint2015arXiv

Three Dimensional MHD Simulation of Circumbinary Accretion Disks -2. Net Accretion Rate

When an accretion disk surrounds a binary rotating in the same sense, the binary exerts strong torques on the gas. Analytic work in the 1D approximation indicated that these torques sharply diminish or even eliminate accretion from the disk onto the binary. However, recent 2D and 3D simulational work has shown at most modest diminution. We present new MHD simulations demonstrating that for binaries with mass ratios of 1 and 0.1 there is essentially no difference between the accretion rate at large radius in the disk and the accretion rate onto the binary. To resolve the discrepancy with earlier analytic estimates, we identify the small subset of gas trajectories traveling from the inner edge of the disk to the binary and show how the full accretion rate is concentrated onto them.

preprint2014arXiv

Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $α$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $α$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

preprint2014arXiv

Observational Signatures of Binary Supermassive Black Holes

Observations indicate that most massive galaxies contain a supermassive black hole, and theoretical studies suggest that when such galaxies have a major merger, the central black holes will form a binary and eventually coalesce. Here we discuss two spectral signatures of such binaries that may help distinguish them from ordinary AGN. These signatures are expected when the mass ratio between the holes is not extreme and the system is fed by a circumbinary disk. One such signature is a notch in the thermal continuum that has been predicted by other authors; we point out that it should be accompanied by a spectral revival at shorter wavelengths and also discuss its dependence on binary properties such as mass, mass ratio, and separation. In particular, we note that the wavelength $λ_n$ at which the notch occurs depends on these three parameters in such a way as to make the number of systems displaying these notches $\propto λ_n^{16/3}$; longer wavelength searches are therefore strongly favored. A second signature, first discussed here, is hard X-ray emission with a Wien-like spectrum at a characteristic temperature $\sim 100$ keV produced by Compton cooling of the shock generated when streams from the circumbinary disk hit the accretion disks around the individual black holes. We investigate the observability of both signatures. The hard X-ray signal may be particularly valuable as it can provide an indicator of black hole merger a few decades in advance of the event.

preprint2013arXiv

A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes

We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption terms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code, and also present results from a number of numerical tests to assess its accuracy and convergence.

preprint2013arXiv

Alignment of supermassive black hole binary orbits and spins

Recent studies of accretion onto supermassive black hole binaries suggest that much, perhaps most, of the matter eventually accretes onto one hole or the other. If so, then for binaries whose inspiral from ~1 pc to 0.001 - 0.01 pc is driven by interaction with external gas, both the binary orbital axis and the individual black hole spins can be reoriented by angular momentum exchange with this gas. Here we show that, unless the binary mass ratio is far from unity, the spins of the individual holes align with the binary orbital axis in a time few-100 times shorter than the binary orbital axis aligns with the angular momentum direction of the incoming circumbinary gas; the spin of the secondary aligns more rapidly than that of the primary by a factor ~(m_1/m_2)^{1/2}>1. Thus the binary acts as a stabilizing agent, so that for gas-driven systems, the black hole spins are highly likely to be aligned (or counteraligned if retrograde accretion is common) with each other and with the binary orbital axis. This alignment can significantly reduce the recoil speed resulting from subsequent black hole merger.

preprint2013arXiv

MHD Simulation of a Disk Subjected to Lense-Thirring Precession

When matter orbits around a central mass obliquely with respect to the mass's spin axis, the Lense-Thirring effect causes it to precess at a rate declining sharply with radius. Ever since the work of Bardeen & Petterson (1975), it has been expected that when a fluid fills an orbiting disk, the orbital angular momentum at small radii should then align with the mass's spin. Nearly all previous work has studied this alignment under the assumption that a phenomenological "viscosity" isotropically degrades fluid shears in accretion disks, even though it is now understood that internal stress in flat disks is due to anisotropic MHD turbulence. In this paper we report a pair of matched simulations, one in MHD and one in pure (non-viscous) HD in order to clarify the specific mechanisms of alignment. As in the previous work, we find that disk warps induce radial flows that mix angular momentum of different orientation; however, we also show that the speeds of these flows are generically transonic and are only very weakly influenced by internal stresses other than pressure. In particular, MHD turbulence does not act in a manner consistent with an isotropic viscosity. When MHD effects are present, the disk aligns, first at small radii and then at large; alignment is only partial in the HD case. We identify the specific angular momentum transport mechanisms causing alignment and show how MHD effects permit them to operate more efficiently. Lastly, we relate the speed at which an alignment front propagates outward (in the MHD case) to the rate at which Lense-Thirring torques deliver angular momentum at smaller radii.

preprint2013arXiv

Relaxation of Warped Disks: the Case of Pure Hydrodynamics

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e. with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square-root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

preprint2013arXiv

Testing Convergence for Global Accretion Disks

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

preprint2013arXiv

X-ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with strong hard X-ray emission: a power-law extending to high energies, a Compton reflection hump, and a broad iron line. Varying only the mass accretion rate, we are able to reproduce a wide range of X-ray states seen in most galactic black hole sources. The temperature in the corona is T_e ~ 10 keV in a boundary layer near the disk and rises smoothly to T_e >~ 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to ~ 6M as the accretion rate decreases, we find that the shape of the Fe Kαline is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

preprint2012arXiv

Circumbinary MHD Accretion into Inspiraling Binary Black Holes

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

preprint2012arXiv

Jet Luminosity From Neutrino-Dominated Accretion Flows in Gamma-Ray Bursts

A hyperaccretion disk formed around a stellar mass black hole is a plausible model for the central engine that powers gamma-ray bursts (GRBs). If the central black hole rotates and a poloidal magnetic field threads its horizon, a powerful relativistic jet may be driven by a process resembling the Blandford-Znajek mechanism. We estimate the luminosity of such a jet assuming that the poloidal magnetic field strength is comparable to the inner accretion disk pressure. We show that the jet efficiency attains its maximal value when the accretion flow is cooled via optically-thin neutrino emission. The jet luminosity is much larger than the energy deposition through neutrino-antineutrino annihilation provided that the black hole is spinning rapidly enough. When the accretion rate onto a rapidly spinning black hole is large enough (> 0.003-0.01M_sun/sec), the predicted jet luminosity is sufficient to drive a GRB.

preprint2012arXiv

Jets from Tidal Disruptions of Stars by Black Holes

Tidal disruption of main sequence stars by black holes has generally been thought to lead to a signal dominated by UV emission. If, however, the black hole spins rapidly and the poloidal magnetic field intensity on the black hole horizon is comparable to the inner accretion disk pressure, a powerful jet may form whose luminosity can easily exceed the thermal UV luminosity. When the jet beam points at Earth, its non-thermal luminosity can dominate the emitted spectrum. The thermal and non-thermal components decay differently with time. In particular, the thermal emission should remain roughly constant for a significant time after the period of maximum accretion, beginning to diminish only after a delay, whereas after the peak accretion rate, the non-thermal jet emission decays, but then reaches a plateau. Both transitions are tied to a characteristic timescale $t_{\rm Edd}$ at which the accretion rate falls below Eddington. Making use of this timescale in a new parameter-inference formalism for tidal disruption events with significant emission from a jet, we analyze the recent flare source Swift J2058. It is consistent with an event in which a main sequence solar-type staris disrupted by a black hole of mass $\sim 4 \times 10^7 M_{\odot}$. The beginning of the flat phase in the non-thermal emission from this source can possibly be seen in the late-time lightcurve. Optical photometry over the first $\simeq 40$ d of this flare is also consistent with this picture, but is only weakly constraining because the bolometric correction is very uncertain. We suggest that future searches for main sequence tidal disruptions use methods sensitive to jet radiation as well as to thermal UV radiation.

preprint2012arXiv

Three Dimensional MHD Simulation of Circumbinary Accretion Disks: Disk Structures and Angular Momentum Transport

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations, strong torques by the binary can maintain a gap of radius $\simeq 2a$. Streams curve inward from $r \simeq 2a$ toward the binary; some of their mass passes through the inner boundary, while the remainder swings back out to the disk. However, we also find that near its inner edge the disk develops both a strong $m=1$ asymmetry and growing orbital eccentricity. Because the MHD stresses introduce more matter into the gap, the total torque per unit disk mass is $\simeq 14$ times larger than found previously. The inner boundary accretion rate per unit disk mass is $\simeq 40$ times greater than found from previous hydrodynamical calculations. The implied binary shrinkage rate is determined by a balance between the rate at which the binary gains angular momentum by accretion and loses it by gravitational torque. The large accretion rate brings these two rates nearly into balance, but in net, we find that $\dot a/a < 0$, and its magnitude is about 2.7 times larger than predicted by the earlier hydrodynamic simulations. If the binary comprises two massive black holes, the accretion rate may be great enough for one or both to be AGN, with consequences for the physical state of the gas both in the disk body and in its inner gap.

preprint2011arXiv

Assessing Quantitative Results in Accretion Simulations: From Local to Global

Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial gridscale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has not often been extensively examined in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly under-resolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

preprint2011arXiv

Dissipation and Vertical Energy Transport in Radiation-Dominated Accretion Disks

Standard models of radiation supported accretion disks generally assume that diffusive radiation flux is solely responsible for vertical heat transport. This requires that heat must be generated at a critical rate per unit volume if the disk is to be in hydrostatic and thermal equilibrium. This raises the question of how heat is generated and how energy is transported in MHD turbulence. By analysis of a number of radiation/MHD stratified shearing-box simulations, we show that the divergence of the diffusive radiation flux is indeed capped at the critical rate, but deep inside the disk, substantial vertical energy flux is also carried by advection of radiation. Work done by radiation pressure is a significant part of the energy budget, and much of this work is dissipated later through damping by radiative diffusion. We show how this damping can be measured in the simulations, and identify its physical origins. Radiative damping accounts for as much as tens of percent of the total dissipation, and is the only realistic physical mechanism for dissipation of turbulence that can actually be resolved in numerical simulations of accretion disks. Buoyancy associated with dynamo-driven, highly magnetized, nearly-isobaric nonlinear slow magnetosonic fluctuations is responsible for the radiation advection flux, and also explains the persistent periodic magnetic upwelling seen at all values of the radiation to gas pressure ratio. The intimate connection between radiation advection and magnetic buoyancy is the first example we know of in astrophysics in which a dynamo has direct impact on the global energetics of a system.

preprint2011arXiv

Radiative efficiency and thermal spectrum of accretion onto Schwarzschild black holes

Recent general relativistic magneto-hydrodynamic (MHD) simulations of accretion onto black holes have shown that, contrary to the basic assumptions of the Novikov-Thorne model, there can be substantial magnetic stress throughout the plunging region. Additional dissipation and radiation can therefore be expected. We use data from a particularly well-resolved simulation of accretion onto a non-spinning black hole to compute both the radiative efficiency of such a flow and its spectrum if all emitted light is radiated with a thermal spectrum whose temperature matches the local effective temperature. This disk is geometrically thin enough (H/r ~= 0.06) that little heat is retained in the flow. In terms of light reaching infinity (i.e., after allowance for all relativistic effects and for photon capture by the black hole), we find that the radiative efficiency is at least ~=6-10% greater than predicted by the Novikov-Thorne model (complete radiation of all heat might yield another ~6%). We also find that the spectrum more closely resembles the Novikov-Thorne prediction for a/M ~= 0.2--0.3 than for the correct value, a/M=0. As a result, if the spin of a non-spinning black hole is inferred by model-fitting to a Novikov-Thorne model with known black hole mass, distance, and inclination, the inferred a/M is too large by ~= 0.2--0.3.

preprint2011arXiv

Swift J1644+57: A White Dwarf Tidally Disrupted by a 10^4 M_{odot} Black Hole?

We propose that the remarkable object Swift J1644+57, in which multiple recurring hard X-ray flares were seen over a span of several days, is a system in which a white dwarf was tidally disrupted by an intermediate mass black hole. Disruption of a white dwarf rather than a main sequence star offers a number of advantages in understanding the multiple, and short, timescales seen in the light curve of this system. In particular, the short internal dynamical timescale of a white dwarf offers a more natural way of understanding the short rise times ($\sim 100$ s) observed. The relatively long intervals between flares ($\sim 5 \times 10^4$ s) may also be readily understood as the period between successive pericenter passages of the remnant white dwarf. In addition, the expected jet power is larger when a white dwarf is disrupted. If this model is correct, the black hole responsible must have mass $\lesssim 10^5 M_{\odot}$.

preprint2010arXiv

Dependence of inner accretion disk stress on parameters: the Schwarzschild case

We explore the parameter dependence of inner disk stress in black hole accretion by contrasting the results of a number of simulations, all employing 3-d general relativistic MHD in a Schwarzschild spacetime. Five of these simulations were performed with the intrinsically conservative code HARM3D, which allows careful regulation of the disk aspect ratio, H/R; our simulations span a range in H/R from 0.06 to 0.17. We contrast these simulations with two previously reported simulations in a Schwarzschild spacetime in order to investigate possible dependence of the inner disk stress on magnetic topology. In all cases, much care was devoted to technical issues: ensuring adequate resolution and azimuthal extent, and averaging only over those time-periods when the accretion flow is in approximate inflow equilibrium. We find that the time-averaged radial-dependence of fluid-frame electromagnetic stress is almost completely independent of both disk thickness and poloidal magnetic topology. It rises smoothly inward at all radii (exhibiting no feature associated with the ISCO) until just outside the event horizon, where the stress plummets to zero. Reynolds stress can also be significant near the ISCO and in the plunging region; the magnitude of this stress, however, depends on both disk thickness and magnetic topology. The two stresses combine to make the net angular momentum accreted per unit rest-mass 7-15% less than the angular momentum of the ISCO.

preprint2010arXiv

Space density of optically-selected type 2 quasars

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5008 luminosity function for 10^8.3 Lsun < L[OIII] < 10^10 Lsun (corresponding to intrinsic luminosities up to M[2400A]-28 mag or bolometric luminosities up to 4x10^47 erg/sec). This luminosity function provides strong lower limits to the actual space density of obscured quasars, due to our selection criteria, the details of the spectroscopic target selection, as well as other effects. We derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars; then, we determine the ratio of type 2/type 1 quasar number densities. Our best data constrain this ratio to be at least 1.5:1 for 10^8.3 Lsun < L[OIII] < 10^9.5 Lsun at z < 0.3, and at least 1.2:1 for L[OIII]=10^10 Lsun at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars in the relatively nearby Universe (z < 0.8) for the highest luminosities.

preprint2009arXiv

Estimating the Prompt Electromagnetic Luminosity of a Black Hole Merger

Although recent work in numerical relativity has made tremendous strides in quantifying the gravitational wave luminosity of black hole mergers, very little is known about the electromagnetic luminosity that might occur in immediate conjunction with these events. We show that whenever the heat deposited in the gas near a pair of merging black holes is proportional to its total mass, and the surface density of the gas in the immediate vicinity is greater than the (quite small) amount necessary to make it optically thick, the characteristic scale of the luminosity emitted in direct association with the merger is the Eddington luminosity independent of the gas mass. The duration of the photon signal is proportional to the gas mass, and is generally rather longer than the merger event. At somewhat larger distances, dissipation associated with realigning the gas orbits to the new spin orientation of the black hole can supplement dissipation of the energy gained from orbital adjustment to the mass lost in gravitational radiation; these two heat sources can combine to augment the electromagnetic radiation over longer timescales.

preprint2009arXiv

Host Galaxies of Luminous Type 2 Quasars at z ~ 0.5

We present deep Gemini GMOS optical spectroscopy of nine luminous quasars at redshifts z ~ 0.5, drawn from the SDSS type 2 quasar sample. Our targets were selected to have high intrinsic luminosities (M_V < -26 mag) as indicated by the [O III] 5007 A emission-line luminosity (L_[O III]). Our sample has a median black hole mass of ~ 10^8.8 M_sun inferred assuming the local M_BH-sigma_* relation and a median Eddington ratio of ~ 0.7, using stellar velocity dispersions sigma_* measured from the G band. We estimate the contamination of the stellar continuum from scattered quasar light based on the strength of broad H-beta, and provide an empirical calibration of the contamination as a function of L_[O III]; the scattered light fraction is ~ 30% of L_5100 for objects with L_[O III] = 10^9.5 L_sun. Population synthesis indicates that young post-starburst populations (< 0.1 Gyr) are prevalent in luminous type 2 quasars, in addition to a relatively old population (> 1 Gyr) which dominates the stellar mass. Broad emission complexes around He II 4686 A with luminosities up to 10^8.3 L_sun are unambiguously detected in three out of the nine targets, indicative of Wolf-Rayet populations. Population synthesis shows that ~ 5-Myr post-starburst populations contribute substantially to the luminosities (> 50% of L_5100) of all three objects with Wolf-Rayet detections. We find two objects with double cores and four with close companions. Our results may suggest that luminous type 2 quasars trace an early stage of galaxy interaction, perhaps responsible for both the quasar and the starburst activity.

preprint2009arXiv

Transport of Large Scale Poloidal Flux in Black Hole Accretion

We report on a global, three-dimensional GRMHD simulation of an accretion torus embedded in a large scale vertical magnetic field orbiting a Schwarzschild black hole. This simulation investigates how a large scale vertical field evolves within a turbulent accretion disk and whether global magnetic field configurations suitable for launching jets and winds can develop. We find that a "coronal mechanism" of magnetic flux motion, which operates largely outside the disk body, dominates global flux evolution. In this mechanism, magnetic stresses driven by orbital shear create large-scale half-loops of magnetic field that stretch radially inward and then reconnect, leading to discontinuous jumps in the location of magnetic flux. In contrast, little or no flux is brought in directly by accretion within the disk itself. The coronal mechanism establishes a dipole magnetic field in the evacuated funnel around the orbital axis with a field intensity regulated by a combination of the magnetic and gas pressures in the inner disk. These results prompt a reevaluation of previous descriptions of magnetic flux motion associated with accretion. Local pictures are undercut by the intrinsically global character of magnetic flux. Formulations in terms of an "effective viscosity" competing with an "effective resistivity" are undermined by the nonlinearity of of the magnetic dynamics and the fact that the same turbulence driving mass motion (traditionally identified as "viscosity") can alter magnetic topology.

preprint2009arXiv

Turbulent Stresses in Local Simulations of Radiation-Dominated Accretion Disks, and the Possibility of the LIghtman-Eardley Instability

We present the results of a series of radiation-MHD simulations of a local patch of an accretion disk, with fixed vertical gravity profile but with different surface mass densities and a broad range of radiation to gas pressure ratios. Each simulation achieves a thermal equilibrium that lasts for many cooling times. After averaging over times long compared to a cooling time, we find that the vertically integrated stress is approximately proportional to the vertically-averaged total thermal (gas plus radiation) pressure. We map out--for the first time on the basis of explicit physics--the thermal equilibrium relation between stress and surface density: the stress decreases (increases) with increasing surface mass density when the simulation is radiation (gas) pressure dominated. The dependence of stress on surface mass density in the radiation pressure dominated regime suggests the possibility of a Lightman-Eardley inflow instability, but global simulations or shearing box simulations with much wider radial boxes will be necessary to confirm this and determine its nonlinear behavior.

preprint2009arXiv

What Is The Numerically Converged Amplitude of MHD Turbulence in Stratified Shearing Boxes?

We study the properties of the turbulence driven by the magnetorotational instability (MRI) in a stratified shearing box with outflow boundary conditions and an equation of state determined by self-consistent dissipation and radiation losses. A series of simulations with increasing resolution are performed within a fixed computational box. We achieve numerical convergence with respect to radial and azimuthal resolution. As vertical resolution is improved, the ratio of stress to pressure increases only slowly, but the absolute levels of both the stress and the pressure increase noticeably. The vertical correlation length of the magnetic field within the core of the disk is highly time-dependent, but averaged over time, it is $\simeq 3$ times larger than found in previous unstratifed simulations. This correlation length decreases slowly as vertical resolution increases. We suggest that this greater degree of correlation in stratified than unstratified disks, even near the midplane, is due to field buoyancy. We further show that the undulatory Parker instability drives the magnetic field upwelling at several scaleheights from the midplane that is characteristic of stratified MHD-turbulent disks.

preprint2009arXiv

X-ray Polarization from Accreting Black Holes: Coronal Emission

We present new calculations of X-ray polarization from accreting black holes (BHs), using a Monte-Carlo ray-tracing code in full general relativity. In our model, an optically thick disk in the BH equatorial plane produces thermal seed photons with polarization oriented parallel to the disk surface. These seed photons are then inverse-Compton scattered through a hot (but thermal) corona, producing a hard X-ray power-law spectrum. We consider three different models for the corona geometry: a wedge ``sandwich'' with aspect ratio H/R and vertically-integrated optical depth tau_0 constant throughout the disk; an inhomogeneous ``clumpy'' corona with a finite number of hot clouds distributed randomly above the disk within a wedge geometry; and a spherical corona of uniform density, centered on the BH and surrounded by a truncated thermal disk with inner radius R_edge. In all cases we find a characteristic transition from horizontal polarization at low energies to vertical polarization above the thermal peak; the vertical direction is defined as the projection of the BH spin axis on the plane of the sky. We show how the details of the spectropolarization signal can be used to distinguish between these models and infer various properties of the corona and BH. Although the bulk of this paper focuses on stellar-mass BHs, we also consider the effects of coronal scattering on the X-ray polarization signal from supermassive BHs in active galactic nuclei.

preprint2009arXiv

X-ray Polarization from Accreting Black Holes: II. The Thermal State

We present new calculations of X-ray polarization from black hole (BH) accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. For moderately large spin parameters (a/M >~ 0.9), this scattered returning radiation dominates the polarization signal at energies above the thermal peak, giving a net rotation in the polarization angle of 90 deg. We show how these new features of the polarization spectra from BHs in the thermal state may be developed into a powerful tool for measuring BH spin and probing the gas flow in the innermost disk. In addition to determining the emission profile, polarization observations can be used to constrain other properties of the system such as BH mass, inclination, and distance. New instruments currently under development should be able to exploit this tool in the near future.

preprint2009arXiv

X-ray Polarization from Black Holes in the Thermal State

We present new calculations of X-ray polarization from black hole accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. We show how these new features of the polarization spectra may be developed into a powerful tool for measuring black hole spin and probing the gas flow in the innermost disk.

preprint2008arXiv

Direct calculation of the radiative efficiency of an accretion disk around a black hole

Numerical simulation of magnetohydrodynamic (MHD) turbulence makes it possible to study accretion dynamics in detail. However, special effort is required to connect inflow dynamics (dependent largely on angular momentum transport) to radiation (dependent largely on thermodynamics and photon diffusion). To this end we extend the flux-conservative, general relativistic MHD code HARM from axisymmetry to full 3D. The use of an energy conserving algorithm allows the energy dissipated in the course of relativistic accretion to be captured as heat. The inclusion of a simple optically thin cooling function permits explicit control of the simulated disk's geometric thickness as well as a direct calculation of both the amplitude and location of the radiative cooling associated with the accretion stresses. Fully relativistic ray-tracing is used to compute the luminosity received by distant observers. For a disk with aspect ratio H/r ~ 0.1 accreting onto a black hole with spin parameter a/M = 0.9, we find that there is significant dissipation beyond that predicted by the classical Novikov-Thorne model. However, much of it occurs deep in the potential, where photon capture and gravitational redshifting can strongly limit the net photon energy escaping to infinity. In addition, with these parameters and this radiation model, significant thermal and magnetic energy remains with the gas and is accreted by the black hole. In our model, the net luminosi ty reaching infinity is 6% greater than the Novikov-Thorne prediction. If the accreted thermal energy were wholly radiated, the total luminosity of the accretion flow would be ~20% greater than the Novikov-Thorne value.

preprint2008arXiv

Mid-infrared spectra of optically selected type 2 quasars

Type 2 quasars are luminous Active Galactic Nuclei whose central engines are seen through large amounts of gas and dust. We present Spitzer spectra of twelve type 2 quasars selected on the basis of their optical emission line properties. Within this sample, we find a surprising diversity of spectra, from those that are featureless to those showing strong PAH emission, deep silicate absorption at 10 micron, hydrocarbon absorption, high-ionization emission lines and H_2 rotational emission lines. About half of the objects in the sample are likely Compton-thick, including the two with the deepest Si absorption. The median star-formation luminosity of the objects in our sample measured from the strength of the PAH features is 5x10^11 L_sun, much higher than for field galaxies or for any other AGN sample, but similar to other samples of type 2 quasars. This suggests an evolutionary link between obscured quasars and peak star formation activity in the host galaxy. Despite the high level of star formation, the bolometric output is dominated by the quasar in all cases. For a given strength of 10 micron Si absorption, ULIRGs are significantly colder than are type 2 quasars (their F_nu[14.5 micron]/F_nu[27.5 micron] ratio is 0.5 dex lower), perhaps reflecting different obscuration geometries in these sources. We find that the appearance of the 10 micron feature (i.e., whether it shows in emission or in absorption) is well-correlated with the optical classification in type 1 and type 2 quasars, contrary to some models of clumpy obscuration. Furthermore, this correlation is significantly stronger in quasars (L_bol>10^45 erg/s) than it is in Seyfert galaxies (L_bol<<10^45 erg/s).

preprint2008arXiv

Radiation-Dominated Disks Are Thermally Stable

When the accretion rate is more than a small fraction of Eddington, the inner regions of accretion disks around black holes are expected to be radiation-dominated. However, in the alpha-model, these regions are also expected to be thermally unstable. In this paper, we report two 3-d radiation MHD simulations of a vertically-stratified shearing box in which the ratio of radiation to gas pressure is ~ 10, and yet no thermal runaway occurs over a timespan ~ 40 cooling times. Where the time-averaged dissipation rate is greater than the critical dissipation rate that creates hydrostatic equilibrium by diffusive radiation flux, the time-averaged radiation flux is held to the critical value, with the excess dissipated energy transported by radiative advection. Although the stress and total pressure are well-correlated as predicted by the alpha-model, we show that stress fluctuations precede pressure fluctuations, contrary to the usual supposition that the pressure controls the saturation level of the magnetic energy. This fact explains the thermal stability. Using a simple toy-model, we show that independently-generated magnetic fluctuations can drive radiation pressure fluctuations, creating a correlation between the two while maintaining thermal stability.

preprint2008arXiv

SDSSJ092712.65+294344.0: NGC 1275 at z=0.7?

SDSSJ092712.65+294344.0 was identified by the SDSS as a quasar, but has the unusual property of having two emission line systems offset by 2650 km/s. One of these contains the usual combination of broad and narrow lines, the other only narrow lines. In the first paper commenting on this system (Komossa et al. 2008), it was interpreted as a galaxy in which a pair of black holes had merged, imparting a several thousand km/s recoil to the new, larger black hole. In two other papers (Bogdanovic, Eracleous & Sigurdsson 2008; Dotti et al. 2008), it was interpreted as a small-separation binary black hole. We propose a new interpretation: that this system is a more distant analog of NGC1275, a large and a small galaxy interacting near the center of a rich cluster.

preprint2008arXiv

The Infrared Afterglow of Supermassive Black Hole Mergers

We model the spectra and light curves of circumbinary accretion disks during the time after the central black holes merge. The most immediate effect of this merger is the dissipation of energy in the outer regions of the disk due to the gravitational wave energy and linear momentum flux released at merger. This has the effect of perturbing the gas in the disk, which then radiates the dissipated energy over a cooling timescale, giving a characteristic infrared signal for tens of thousands of years when the total black hole mass is M~10^8 M_sun. On the basis of a simple cosmological merger model in which a typical supermassive black hole undergoes a few major mergers during its lifetime, we predict that ~10^4-10^5 of these IR sources should be observable today and discuss the possibility of identifying them with multi-wavelength surveys such as SWIRE/XMM-LSS/XBootes and COSMOS.

preprint2007arXiv

AGN Obscuring Tori Supported by Infrared Radiation Pressure

Explicit 2-d axisymmetric solutions are found to the hydrostatic equilibrium, energy balance, and photon diffusion equations within obscuring tori around active galactic nuclei. These solutions demonstrate that infrared radiation pressure can support geometrically thick structures in AGN environments subject to certain constraints: the bolometric luminosity must be roughly 0.03--1 times the Eddington luminosity; and the Compton optical depth of matter in the equatorial plane should be order unity, with a tolerance of about an order of magnitude up or down. Both of these constraints are at least roughly consistent with observations. In addition, angular momentum must be redistributed so that the fractional rotational support against gravity rises from the inner edge of the torus to the outer in a manner specific to the detailed shape of the gravitational potential. This model also predicts that the column densities observed in obscured AGN should range from about 10^{22} to about 10^{24} cm^{-2}.

preprint2006arXiv

The Weak Field Limit of the Magnetorotational Instability

We investigate the behavior of the magneto-rotational instability in the limit of extremely weak magnetic field, i.e., as the ratio of ion cyclotron frequency to orbital frequency (X) becomes small. Considered only in terms of cold two-fluid theory, instability persists to arbitrarily small values of X, and the maximum growth rate is of order the orbital frequency except for the range m_e/m_i < |X| < 1, where it can be rather smaller. In this range, field aligned with rotation (X > 0) produces slower growth than anti-aligned field (X < 0). The maximum growth rate is generally achieved at smaller and smaller wavelengths as |X| diminishes. When |X| < m_e/m_i, new unstable "electromagnetic-rotational" modes appear that do not depend on the equilibrium magnetic field. Because the most rapidly-growing modes have extremely short wavelengths when |X| is small, they are often subject to viscous or resistive damping, which can result in suppressing all but the longest wavelengths, for which growth is much slower. We find that this sort of damping is likely to curtail severely the frequently-invoked mechanism for cosmological magnetic field growth in which a magnetic field seeded by the Biermann battery is then amplified by the magneto-rotational instability. On the other hand, the small |X| case may introduce interesting effects in weakly-ionized disks in which dust grains carry most of the electric charge.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$α$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$λ5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$α$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

preprint1996arXiv

The Vertical Structure and Ultraviolet Spectrum of X-Ray Irradiated Accretion Disks in Active Galactic Nuclei

Motivated by recent work indicating that the UV continuum in AGN may be produced by reradiation of energy absorbed from X-rays irradiating an accretion disk, we present a calculation of the vertical structures and ultraviolet spectra of X-ray irradiated accretion disks around massive non-rotating black holes. After finding the radial dependence of vertically-integrated quantities for these disks, we solve the equations of hydrostatic equilibrium, energy balance, and frequency-dependent radiation transfer as functions of altitude. To solve the last set of equations, we use a variable Eddington factor method. We include electron scattering, free-free, and HI, HeI, and HeII bound-free opacities and the corresponding continuum cooling processes. While the overall spectral shape predicted by X-ray irradiation may be compatible with observations, the Lyman edge emission feature it predicts is not. This finding raises questions for many otherwise plausible models in which X-ray irradiation plays a major role.

Julian H. Krolik

What is connected

Connect this record

See the researcher in context

Building this map preview

56 published item(s)

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

What powers the radio emission in TDE AT2019dsg: a long-lived jet or the disruption itself?

High-energy emission from tidal disruption events in active galactic nuclei

Future Simulations of Tidal Disruption Events

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

Radiative Emission Mechanisms of Tidal Disruption Events

Tidal disruption events in active galactic nuclei

Magnetorotational instability in eccentric disks

Discovery of transient infrared emission from dust heated by stellar tidal disruption flares

Fe K$α$ Profiles from Simulations of Accreting Black Holes

How bright are the gaps in circumbinary disk systems?

Radiation-driven outflows from and radiative support in dusty tori of active galactic nuclei

A Steady-State Alignment Front in an Accretion Disk Subjected to Lense-Thirring Torques

Disk Emission from Magneto-hydrodynamic Simulations of Spinning Black Holes

Evolution of a Binary Black Hole with a Retrograde Circumbinary Accretion Disk

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

Three Dimensional MHD Simulation of Circumbinary Accretion Disks -2. Net Accretion Rate

Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst

Observational Signatures of Binary Supermassive Black Holes

A Monte Carlo Code for Relativistic Radiation Transport Around Kerr Black Holes

Alignment of supermassive black hole binary orbits and spins

MHD Simulation of a Disk Subjected to Lense-Thirring Precession

Relaxation of Warped Disks: the Case of Pure Hydrodynamics

Testing Convergence for Global Accretion Disks

X-ray Spectra from MHD Simulations of Accreting Black Holes

Circumbinary MHD Accretion into Inspiraling Binary Black Holes

Jet Luminosity From Neutrino-Dominated Accretion Flows in Gamma-Ray Bursts

Jets from Tidal Disruptions of Stars by Black Holes

Three Dimensional MHD Simulation of Circumbinary Accretion Disks: Disk Structures and Angular Momentum Transport

Assessing Quantitative Results in Accretion Simulations: From Local to Global

Dissipation and Vertical Energy Transport in Radiation-Dominated Accretion Disks

Radiative efficiency and thermal spectrum of accretion onto Schwarzschild black holes

Swift J1644+57: A White Dwarf Tidally Disrupted by a 10^4 M_{odot} Black Hole?

Dependence of inner accretion disk stress on parameters: the Schwarzschild case

Space density of optically-selected type 2 quasars

Estimating the Prompt Electromagnetic Luminosity of a Black Hole Merger

Host Galaxies of Luminous Type 2 Quasars at z ~ 0.5

Transport of Large Scale Poloidal Flux in Black Hole Accretion

Turbulent Stresses in Local Simulations of Radiation-Dominated Accretion Disks, and the Possibility of the LIghtman-Eardley Instability

What Is The Numerically Converged Amplitude of MHD Turbulence in Stratified Shearing Boxes?

X-ray Polarization from Accreting Black Holes: Coronal Emission

X-ray Polarization from Accreting Black Holes: II. The Thermal State

X-ray Polarization from Black Holes in the Thermal State

Direct calculation of the radiative efficiency of an accretion disk around a black hole

Mid-infrared spectra of optically selected type 2 quasars

Radiation-Dominated Disks Are Thermally Stable

SDSSJ092712.65+294344.0: NGC 1275 at z=0.7?

The Infrared Afterglow of Supermassive Black Hole Mergers

AGN Obscuring Tori Supported by Infrared Radiation Pressure

The Weak Field Limit of the Magnetorotational Instability

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

The Vertical Structure and Ultraviolet Spectrum of X-Ray Irradiated Accretion Disks in Active Galactic Nuclei