BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shane W. Davis

Shane W. Davis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.IMastro-phastro-ph.COgr-qchep-phastro-ph.SRhep-thphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

31works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

31 published item(s)

preprint2022arXiv

CR Driven Multi-phase Gas Formed via Thermal Instability

Cosmic rays (CRs) are an important energy source in the circum-galactic medium (CGM) that impact the multi-phase gas structure and dynamics. We perform two-dimensional CR-magnetohydrodynamic simulations to investigate the role of CRs in accelerating multi-phase gas formed via thermal instability. We compare outflows driven by CRs to those driven by a hot wind with equivalent momentum. We find that CRs driven outflow produces lower density contrast between cold and hot gas due to non-thermal pressure support, and yields a more filamentary cloud morphology. While entrainment in a hot wind can lead to cold gas increasing due to efficient cooling, CRs tend to suppress cold gas growth. The mechanism of this suppression depends on magnetic field strength, with CRs either reducing cooling or shredding the clouds by differential acceleration. Despite the suppression of cold gas growth, CRs are able to launch the cold clouds to observed velocities without rapid destruction. The dynamical interaction between CRs ad multi-phase gas is also sensitive to the magnetic field strength. In relatively strong fields, the CRs are more important for direct momentum input to cold gas. In relatively weak fields, the CRs impact gas primarily by heating, which modifies gas pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ice Age : Chemo-dynamical modeling of Cha-MMS1 to predict new solid-phase species for detection with JWST

Chemical models and experiments indicate that interstellar dust grains and their ice mantles play an important role in the production of complex organic molecules (COMs). To date, the most complex solid-phase molecule detected with certainty in the ISM is methanol, but the James Webb Space Telescope (JWST) may be able to identify still larger organic species. In this study, we use a coupled chemo-dynamical model to predict new candidate species for JWST detection toward the young star-forming core Cha-MMS1, combining the gas-grain chemical kinetic code MAGICKAL with a 1-D radiative hydrodynamics simulation using Athena++. With this model, the relative abundances of the main ice constituents with respect to water toward the core center match well with typical observational values, providing a firm basis to explore the ice chemistry. Six oxygen-bearing COMs (ethanol, dimethyl ether, acetaldehyde, methyl formate, methoxy methanol, and acetic acid), as well as formic acid, show abundances as high as, or exceeding, 0.01% with respect to water ice. Based on the modeled ice composition, the infrared spectrum is synthesized to diagnose the detectability of the new ice species. The contribution of COMs to IR absorption bands is minor compared to the main ice constituents, and the identification of COM ice toward the core center of Cha-MMS1 with the JWST NIRCAM/Wide Field Slitless Spectroscopy (2.4-5.0 micron) may be unlikely. However, MIRI observations (5-28 micron) toward COM-rich environments where solid-phase COM abundances exceed 1% with respect to the water ice column density might reveal the distinctive ice features of COMs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Launching of Cosmic Ray Driven Outflows

Cosmic rays (CRs) are thought to be an important feedback mechanism in star-forming galaxies. They can provide an important source of pressure support and possibly drive outflows. We perform multidimensional CR-magnetohydrodynamic simulations including transport by streaming and diffusion to investigate wind launching from an initially hydrostatic atmosphere by CRs. We estimate a characteristic Eddington limit on the CR flux for which the CR force exceeds gravity and compare it to simulated systems. Scaling our results to conditions in star-forming galaxies, we find that CRs are likely to contribute to driving outflows for a broad range of star formation environments. We quantify the momentum and energy transfer between CRs and gas, along with the associated mass outflow rates under different assumptions about the relative importance of streaming and diffusion for transport. In simulations with streaming, we observe the growth and saturation of the CR acoustic instability, but the CRs and gas remain well coupled, with CR momentum transferred efficiently to the gas even when this instability is present. Higher CR fluxes transferr more energy to the gas and drive stronger outflows. When streaming is present, most of the transferred energy takes the form of Alfvén wave heating of the gas, raising its pressure and internal energy, with a lower fractional contribution to the kinetic energy of the outflow. We also consider runs with radiative cooling, which modifies gas temperature and pressure profiles but does not seem to have a large impact on the mass outflow for super-Eddington CR fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets

There is a broad consensus that accretion onto supermassive black holes and consequent jet formation power the observed emission from active galactic nuclei (AGNs). However, there has been less agreement about how jets form in accretion flows, their possible relationship to black hole spin, and how they interact with the surrounding medium. There have also been theoretical concerns about instabilities in standard accretion disk models and lingering discrepancies with observational constraints. Despite seemingly successful applications to X-ray binaries, the standard accretion disk model faces a growing list of observational constraints that challenge its application to AGNs. Theoretical exploration of these questions has become increasingly reliant on numerical simulations owing to the dynamic nature of these flows and the complex interplay between hydrodynamics, magnetic fields, radiation transfer, and curved spacetime. We conclude the following: The advent of general relativistic magnetohydrodynamics (MHD) simulations has greatly improved our understanding of jet production and its dependence on black hole spin. Simulation results show both disks and jets are sensitive to the magnetic flux threading the accretion flow as well as possible misalignment between the angular momentum of the accretion flow and the black hole spin. Radiation MHD simulations are providing new insights into the stability of luminous accretion flows and highlighting the potential importance of radiation viscosity, UV opacity from atoms, and spiral density waves in AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dusty Cloud Acceleration with Multiband Radiation

We perform two-dimensional and three-dimensional simulations of cold, dense clouds, which are accelerated by radiation pressure on dust relative to a hot, diffuse background gas. We examine the relative effectiveness of acceleration by ultraviolet and infrared radiation fields, both independently and acting simultaneously on the same cloud. We study clouds that are optically thin to infrared emission but with varying ultraviolet optical depths. Consistent with previous work, we find relatively efficient acceleration and long cloud survival times when the infrared band flux dominates over the ultraviolet flux. However, when ultraviolet is dominant or even a modest percentage ($\sim 5-10$\%) of the infrared irradiating flux, it can act to compress the cloud, first crushing it and then disrupting the outer layers. This drives mixing of outer regions of the dusty gas with the hot diffuse background to the point where most dust is not likely to survive or stay coupled to the gas. Hence, the cold cloud is unable to survive for a long enough timescale to experience significant acceleration before disruption even though efficient infrared cooling keeps the majority of the gas close to radiative equilibrium temperature ($T \lesssim 100$K). We discuss implications for observed systems, concluding that radiation pressure driving is most effective when the light from star-forming regions is efficiently reprocessed into the infrared.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time Dependent Radiation Hydrodynamics on a Moving Mesh

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Covariant Radiative Transfer for Black Hole Spacetimes

It has now become possible to study directly, via numerical simulation, the evolution of relativistic, radiation-dominated flows around compact objects. With this in mind we set out explicitly covariant forms of the radiative transfer equation that are suitable for numerical integration in curved spacetime or flat spacetime in curvilinear coordinates. Our work builds on and summarizes in consistent form earlier work by Lindquist, Thorne, Morita and Kaneko, and others. We give explicitly the basic equations in spherical-polar coordinates for Minkowski space and the Kerr spacetime in Kerr-Schild coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Global Radiation Magneto-hydrodynamic Simulations of Sub-Eddington Accretion Disks around Supermassive Black Holes

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A broadband X-ray spectral study of the intermediate-mass black hole candidate M82 X-1 with NuSTAR, Chandra and Swift

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiation Magneto-hydrodynamic Simulations of the Formation of Hot Accretion Disk Coronae

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of warped accretion discs

Accretion discs are present around both stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. A wide variety of circumstantial evidence implies that many of these discs are warped. The standard Bardeen--Petterson model attributes the shape of the warp to the competition between Lense--Thirring torque from the central black hole and viscous angular-momentum transport within the disc. We show that this description is incomplete, and that torques from the companion star (for X-ray binaries) or the self-gravity of the disc (for active galactic nuclei) can play a major role in determining the properties of the warped disc. Including these effects leads to a rich set of new phenomena. For example, (i) when a companion star is present and the warp arises from a misalignment between the companion's orbital axis and the black hole's spin axis, there is no steady-state solution of the Pringle--Ogilvie equations for a thin warped disc when the viscosity falls below a critical value; (ii) in AGN accretion discs, the warp can excite short-wavelength bending waves that propagate inward with growing amplitude until they are damped by the disc viscosity. We show that both phenomena can occur for plausible values of the black hole and disc parameters, and briefly discuss their observational implications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Line driven winds and the UV turnover in AGN accretion discs

AGN SEDs generally show a turnover at lambda 1000A, implying a maximal accretion disc (AD) temperature of T_max~50,000K. Massive O stars display a similar T_max, associated with a sharp rise in a line driven mass loss Mdot_wind with increasing surface temperature. AGN AD are also characterized by similar surface gravity to massive O stars. The Mdot_wind of O stars reaches ~10^-5 Msun/year. Since the surface area of AGN AD can be 10^6 larger, the implied Mdot_wind in AGN AD can reach the accretion rate Mdot. A rise to Mdot_wind Mdot towards the AD center may therefore set a similar cap of T_max~50,000K. To explore this idea, we solve the radial structure of an AD with a mass loss term, and calculate the implied AD emission using the mass loss term derived from observations of O stars. We find that Mdot_wind becomes comparable to Mdot typically at a few 10s of GM/c^2. Thus, the standard thin AD solution is effectively truncated well outside the innermost stable orbit. The calculated AD SED shows the observed turnover at lambda~1000A, which is weakly dependent on the AGN luminosity and black hole mass. The AD SED is generally independent of the black hole spin, due to the large truncation radius. However, a cold AD (low Mdot, high black hole mass) is predicted to be windless, and thus its SED should be sensitive to the black hole spin. The accreted gas may form a hot thick disc with a low radiative efficiency inside the truncation radius, or a strong line driven outflow, depending on its ionization state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Saturation of the MRI in Strongly Radiation Dominated Accretion Disks

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Radiation Transfer Solver for Athena using Short Characteristics

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-linear Evolution of Rayleigh-Taylor Instability in a Radiation Supported Atmosphere

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Eye of the Storm: Light from the Inner Plunging Region of Black Hole Accretion Discs

It is generally thought that the light coming from the inner plunging region of black hole accretion discs contributes negligibly to the disc's overall spectrum, i.e. the plunging fluid is swallowed by the black hole before it has time to radiate. In the standard disc model used to fit X-ray observations of accretion discs, the plunging region is assumed to be perfectly dark. However, numerical simulations that include the full physics of the magnetized flow predict that a small fraction of the disc's total luminosity emanates from the plunging region. We investigate the observational consequences of this neglected inner light. We compute radiative transfer based disc spectra that correspond to 3D general relativistic magnetohydrodynamic simulated discs (which produce light inside their plunging regions). In the context of black hole spin estimation, we find that the neglected inner light only has a modest effect (this bias is less than typical observational systematic errors). For rapidly spinning black holes, we find that the combined emission from the plunging region produces a weak power-law tail at high energies. This indicates that infalling matter is the origin for some of the `coronal' emission observed in the thermal dominant and steep power-law states of X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cold Accretion Disks and Lineless Quasars

The optical-UV continuum of quasars is broadly consistent with the emission from a geometrically thin optically thick accretion disk (AD). The AD produces the ionizing continuum which powers the broad and narrow emission lines. The maximum AD effective temperature is given by Teff=fmax(Mdot/M^2)^1/4, where M is the black hole mass, Mdot the accretion rate, and fmax is set by the black hole spin a_*. For a low enough value of Mdot/M^2 the AD may become too cold to produce ionizing photons. Such an object will form a lineless quasar. This occurs for a local blackbody (BB) AD with a luminosity Lopt=10^46 erg/s for M>3.6E9 Msun, when a_*=0, and for M>1.4E10 Msun, when a_*=0.998. Using the AD based Mdot, derived from M and Lopt, and the reverberation based M, derived from Lopt and the Hbeta FWHM, v, gives Teff \propto Lopt^-0.13v^-1.45. Thus, Teff is mostly set by v. Quasars with a local BB AD become lineless for v> 8,000 km/s, when a_*=0, and for v> 16,000 km/s, when a_*=0.998. Higher values of v are required if the AD is hotter than a local BB. The AD becoming non-ionizing may explain why line emitting quasars with v>10,000 km/s are rare. Weak low ionization lines may still be present if the X-ray continuum is luminous enough, and such objects may form a population of weak emission line quasars (WLQ). If correct, such WLQ should show a steeply falling SED at lambda<1000A. Such an SED was observed by Hryniewicz et al. in SDSS J094533.99+100950.1, a WLQ observed down to 570A, which is well modeled by a rather cold AD SED. UV spectroscopy of z~1-2 quasars is required to eliminate potential intervening Lyman limit absorption by the intergalactic medium (IGM), and to explore if the SEDs of lineless quasars and some additional WLQ are also well fit by a cold AD SED.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring Black Hole Spin by the Continuum-Fitting Method: Effect of Deviations from the Novikov-Thorne Disc Model

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring the Spins of Accreting Black Holes

A typical galaxy is thought to contain tens of millions of stellar-mass black holes, the collapsed remnants of once massive stars, and a single nuclear supermassive black hole. Both classes of black holes accrete gas from their environments. The accreting gas forms a flattened orbiting structure known as an accretion disk. During the past several years, it has become possible to obtain measurements of the spins of the two classes of black holes by modeling the X-ray emission from their accretion disks. Two methods are employed, both of which depend upon identifying the inner radius of the accretion disk with the innermost stable circular orbit (ISCO), whose radius depends only on the mass and spin of the black hole. In the Fe K method, which applies to both classes of black holes, one models the profile of the relativistically-broadened iron line with a special focus on the gravitationally redshifted red wing of the line. In the continuum-fitting method, which has so far only been applied to stellar-mass black holes, one models the thermal X-ray continuum spectrum of the accretion disk. We discuss both methods, with a strong emphasis on the continuum-fitting method and its application to stellar-mass black holes. Spin results for eight stellar-mass black holes are summarized. These data are used to argue that the high spins of at least some of these black holes are natal, and that the presence or absence of relativistic jets in accreting black holes is not entirely determined by the spin of the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Cool Accretion Disk in ESO 243-49 HLX-1: Further Evidence of an Intermediate Mass Black Hole

With an inferred bolometric luminosity exceeding 10^42 erg/s, HLX-1 in ESO 243-49 is the most luminous of ultraluminous X-ray sources and provides one of the strongest cases for the existence of intermediate mass black holes. We obtain good fits to disk-dominated observations of the source with BHSPEC, a fully relativistic black hole accretion disk spectral model. Due to degeneracies in the model arising from the lack of independent constraints on inclination and black hole spin, there is a factor of 100 uncertainty in the best-fit black hole mass M. Nevertheless, spectral fitting of XMM-Newton observations provides robust lower and upper limits with 3000 Msun < M < 3 x 10^5 Msun, at 90% confidence, placing HLX-1 firmly in the intermediate-mass regime. The lower bound on M is entirely determined by matching the shape and peak energy of the thermal component in the spectrum. This bound is consistent with (but independent of) arguments based solely on the Eddington limit. Joint spectral modelling of the XMM-Newton data with more luminous Swift and Chandra observations increases the lower bound to 6000 Msun, but this tighter constraint is not independent of the Eddington limit. The upper bound on M is sensitive to the maximum allowed inclination i, and is reduced to M < 10^5 Msun if we limit i < 75 deg.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a/M>0.95 (3σ). For a less probable (synchronous) dynamical model, we find a/M>0.92 (3σ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The X-Ray Polarization Signature of Quiescent Magnetars: Effect of Magnetospheric Scattering and Vacuum Polarization

In the magnetar model, the quiescent non-thermal soft X-ray emission from Anomalous X-ray Pulsars and Soft-Gamma Repeaters is thought to arise from resonant comptonization of thermal photons by charges moving in a twisted magnetosphere. Robust inference of physical quantities from observations is difficult, because the process depends strongly on geometry and current understanding of the magnetosphere is not very deep. The polarization of soft X-ray photons is an independent source of information, and its magnetospheric imprint remains only partially explored. In this paper we calculate how resonant cyclotron scattering would modify the observed polarization signal relative to the surface emission, using a multidimensional Monte Carlo radiative transfer code that accounts for the gradual coupling of polarization eigenmodes as photons leave the magnetosphere. We employ a globally-twisted, self-similar, force-free magnetosphere with a power-law momentum distribution, assume a blackbody spectrum for the seed photons, account for general relativistic light deflection close to the star, and assume that vacuum polarization dominates the dielectric properties of the magnetosphere. The latter is a good approximation if the pair multiplicity is not much larger than unity. Phase-averaged polarimetry is able to provide a clear signature of the magnetospheric reprocessing of thermal photons and to constrain mechanisms generating the thermal emission. Phase-resolved polarimetry, in addition, can characterize the spatial extent and magnitude of the magnetospheric twist angle at ~100 stellar radii, and discern between uni- or bidirectional particle energy distributions, almost independently of every other parameter in the system. We discuss prospects for detectability with GEMS.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precise Measurement of the Spin Parameter of the Stellar-Mass Black Hole M33 X-7

In prior work, {\it Chandra} and Gemini-North observations of the eclipsing X-ray binary M33 X-7 have yielded measurements of the mass of its black hole primary and the system's orbital inclination angle of unprecedented accuracy. Likewise, the distance to the binary is known to a few percent. In an analysis based on these precise results, fifteen {\it Chandra} and {\it XMM-Newton} X-ray spectra, and our fully relativistic accretion disk model, we find that the dimensionless spin parameter of the black hole primary is $a_* = 0.77 \pm 0.05$. The quoted 1-$σ$ error includes all sources of observational uncertainty. Four {\it Chandra} spectra of the highest quality, which were obtained over a span of several years, all lead to the same estimate of spin to within statistical errors (2%), and this estimate is confirmed by 11 spectra of lower quality. There are two remaining uncertainties: (1) the validity of the relativistic model used to analyze the observations, which is being addressed in ongoing theoretical work; and (2) our assumption that the black hole spin is approximately aligned with the angular momentum vector of the binary, which can be addressed by a future X-ray polarimetry mission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing Accretion Disk Structure with Suzaku data of LMC X-3

The Suzaku observation of LMC X-3 gives the best data to date on the shape of the accretion disk spectrum. This is due to the combination of very low absorbing column density along this line of sight which allows the shape of the disk emisison to be constrained at low energies by the CCD's, while the tail can be simultaneously determined up to 30 keV by the high energy detectors. These data clearly demonstrate that the observed disk spectrum is broader than a simple 'sum of blackbodies', and relativistic smearing of the emission is strongly required. However, the intrinsic emission should be more complex than a (color-corrected) sum of blackbodies as it should also contain photo-electric absorption edges from the partially ionised disk photosphere. These are broadened by the relativistic smearing, but the models predict ~ 3-5 per cent deviations for 1/3- 1 solar abundance around the edge energies, significantly stronger than observed. This indicate that the models need to include more physical processes such as self-irradiation, bound-bound (line) absorption, and/or emission from recombination continuua and/or lines. Alternatively, if none of these match the data, it may instead require that the accretion disk density and/or emissivity profile with height is different to that assumed. Thus these data demonstrate the feasibility of observational tests of our fundamental understanding of the vertical structure of accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Radiative Efficiency of Accretion Flows in Individual AGN

The radiative efficiency of AGN is commonly estimated based on the total mass accreted and the total AGN light emitted per unit volume in the universe integrated over time (the Soltan argument). In individual AGN, thin accretion disk model spectral fits can be used to deduce the absolute accretion rate Mdot, if the black hole mass M is known. The radiative efficiency η is then set by the ratio of the bolometric luminosity L_bol to Mdot c^2. We apply this method to determine η in a sample of 80 PG quasars with well determined L_bol, where Mdot is set by thin accretion disk model fits to the optical luminosity density, and the M determination based on the bulge stellar velocity dispersion (13 objects) or the broad line region (BLR). For the BLR-based masses, we derive a mean log η = -1.05 +/- 0.52 consistent with the Soltan argument based estimates. We find a strong correlation of η with M, rising from η ~ 0.03 at M = 10^7 M{\odot} and L/L_Edd ~ 1 to η ~ 0.4 at M = 10^9 M{\odot} and L/L_Edd ~ 0.3. This trend is related to the overall uniformity of L_opt/L_bol in our sample, particularly the lack of the expected increase in L_opt/L_bol with increasing M (and decreasing L/L_Edd), which is a generic property of thermal disk emission at fixed η. The significant uncertainty in the M determination is not large enough to remove the correlation. The rising η with M may imply a rise in the black hole spin with M, as proposed based on other indirect arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Eddington Limit in Cosmic Rays: An Explanation for the Observed Faintness of Starbursting Galaxies

We show that the luminosity of a star forming galaxy is capped by the production and subsequent expulsion of cosmic rays from its interstellar medium. By defining an Eddington luminosity in cosmic rays, we show that the star formation rate of a given galaxy is limited by its mass content and the cosmic ray mean free path. When the cosmic ray luminosity and pressure reaches a critical value as a result of vigorous star formation, hydrostatic balance is lost, a cosmic ray-driven wind develops, and star formation is choked off. Cosmic ray pressure-driven winds are likely to produce wind velocities significantly in excess of the galactic escape velocity. It is possible that cosmic ray feedback results in the Faber-Jackson relation for a plausible set of input parameters that describe cosmic ray production and transport, which are calibrated by observations of the Milky Way's interstellar cosmic rays as well as other galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Spin of the Near-Extreme Kerr Black Hole GRS 1915+105

Based on a spectral analysis of the X-ray continuum that employs a fully relativistic accretion-disk model, we conclude that the compact primary of the binary X-ray source GRS 1915+105 is a rapidly-rotating Kerr black hole. We find a lower limit on the dimensionless spin parameter of a* greater than 0.98. Our result is robust in the sense that it is independent of the details of the data analysis and insensitive to the uncertainties in the mass and distance of the black hole. Furthermore, our accretion-disk model includes an advanced treatment of spectral hardening. Our data selection relies on a rigorous and quantitative definition of the thermal state of black hole binaries, which we used to screen all of the available RXTE and ASCA data for the thermal state of GRS 1915+105. In addition, we focus on those data for which the accretion disk luminosity is less than 30% of the Eddington luminosity. We argue that these low-luminosity data are most appropriate for the thin alpha-disk model that we employ. We assume that there is zero torque at the inner edge of the disk, as is likely when the disk is thin, although we show that the presence of a significant torque does not affect our results. Our model and the model of the relativistic jets observed for this source constrain the distance and black hole mass and could thus be tested by determining a VLBA parallax distance and improving the measurement of the mass function. Finally, we comment on the significance of our results for relativistic-jet and core-collapse models, and for the detection of gravitational waves.

0 citations0 reviewsNo code linkedOpen access