Showing 1–50 of 61 results for author: Davis, S W

Pre-peak Emission in Tidal Disruption Events

Authors: Xiaoshan Huang, Shane W. Davis, Yan-fei Jiang

Abstract: The rising part of a tidal disruption event light curve provides unique insight into early emission and the onset of accretion. Various mechanisms are proposed to explain the pre-peak emission, including shocks from debris interaction and reprocessing of disk emission. We study the pre-peak emission and its influence on the gas circularization by a series of gray radiation hydrodynamic simulations… ▽ More The rising part of a tidal disruption event light curve provides unique insight into early emission and the onset of accretion. Various mechanisms are proposed to explain the pre-peak emission, including shocks from debris interaction and reprocessing of disk emission. We study the pre-peak emission and its influence on the gas circularization by a series of gray radiation hydrodynamic simulations with varying black hole mass. We find that given a super-Eddington fallback rate of 10\dot{M}_{Edd}, the stream-stream collision can occur multiple times and drive strong outflows of up to 9\dot{M}_{Edd}. By dispersing gas to \gtrsim 100rs, the outflow can delay gas circularization and leads to sub-Eddington accretion rates during the first few stream-stream collisions. The stream-stream collision shock and circularization shock can sustain a luminosity of ~10^{44}erg/s for days. The luminosity is generally sub-Eddington and shows a weak correlation with accretion rate at early time. The outflow is optically thick, yielding a reprocessing layer with a size of ~10^{14} cm and photospheric temperature of ~4\times10^{4}K. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 21 Pages, 19 figures, submitted to ApJ, comments welcome

Line Driven Instabilities due to Continuum Radiation Transport in Stellar Winds

Authors: Sergei Dyda, Shane W. Davis

Abstract: We study line driven stellar winds using time-dependent radiation hydrodynamics where the continuum radiation couples to the gas via either a scattering or absorption opacity and there is an additional radiation force due to spectral lines that we model in the Sobolev approximation. We find that in winds with scattering opacities, instabilties tend to be suppressed and the wind reaches a steady st… ▽ More We study line driven stellar winds using time-dependent radiation hydrodynamics where the continuum radiation couples to the gas via either a scattering or absorption opacity and there is an additional radiation force due to spectral lines that we model in the Sobolev approximation. We find that in winds with scattering opacities, instabilties tend to be suppressed and the wind reaches a steady state. Winds with absorption opacities are unstable and remain clumpy at late times. Clumps persist because they are continually regenerated in the subcritical part of the flow. Azimuthal gradients in the radial velocity distribution cause a drop in the radial radiation force and provide a mechanism for generating clumps. These clumps form on super-Sobolev scales, but at late times become Sobolev-length sized indicating that our radiation transfer model is breaking down. Inferring the clump distribution at late times therefore requires radiation-hydrodynamic modeling below the Sobolev scale. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 10 pages, 8 figures

Time-Dependent AGN Disc Winds I -- X-ray Irradiation

Authors: Sergei Dyda, Shane W. Davis, Daniel Proga

Abstract: We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing rad… ▽ More We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing radiation is affected by scattering, absorption and re-emission by the intervening gas. In a pure attenuation model, the disc launches an episodic wind, as previous studies have shown. Including scattering or re-emission tends to weaken the wind, lowering the mass flux and outflow velocity and if sufficiently dominant, suppressing the outflow entirely. However, the exponential nature of radiative attenuation means only a modest, factor of a few, increase in the absorption cross section can overcome the wind suppression due to scattering and re-emission. We find mass outflow rates of $\sim 20\%$ or more of the assumed inflow rate through the disk, indicating that radiation driven winds may significantly alter the structure of the accretion flow. The winds also supply a large, time-varying column of material above the nominal constant disk scale height, which will determine the geometry of reprocessed emission from the central source. Our results suggest the need for accurate photoionization modeling, radiation transport as well as accretion disc physics, to study their effects on the AGN disc winds △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 13 pages, 4 figures

Spin measurement of 4U 1543-47 with Insight-HXMT and NICER from its 2021 outburst: A test of accretion disk models at high luminosities

Authors: E. S. Yorgancioglu, Q. C. Bu, A. Santangelo, L. Tao, S. W. Davis, A. Vahdat, L. D. Kong, S. Piraino, M. Zhou, S. N. Zhang

Abstract: 4U 1543--47 is one of a handful of known black hole candidates located in the Milky Way Galaxy, and has undergone a very bright outburst in 2021, reaching a total of $\sim$9 Crab, as observed by the Monitor of All-sky Image (MAXI), and exceeding twice its Eddington luminosity. The unprecedented bright outburst of 4U 1543--47 provides a unique opportunity to test the behavior of accretion disk mode… ▽ More 4U 1543--47 is one of a handful of known black hole candidates located in the Milky Way Galaxy, and has undergone a very bright outburst in 2021, reaching a total of $\sim$9 Crab, as observed by the Monitor of All-sky Image (MAXI), and exceeding twice its Eddington luminosity. The unprecedented bright outburst of 4U 1543--47 provides a unique opportunity to test the behavior of accretion disk models at high luminosities and accretion rates. In addition, we explore the possibility of constraining the spin of the source at high accretion rates, given that previous spin measurements of 4U 1543--47 have been largely inconsistent with each other. We measure the spectral evolution of the source throughout its outburst as observed by Insight-HXMT, and compare the behavior of both the thin disk model kerrbb2, as well as the slim disk model slimbh up to the Eddington limit for two different values of disk $α$-viscosity. In addition, given the behavior of these two models, we identify two `golden' epochs for which it is most suitable to measure the spin with continuum fitting. △ Less

Submitted 21 July, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 10 pages, 6 figures

Journal ref: A&A 677, A79 (2023)

Spectral calculations of 3D RMHD simulations of super-Eddington accretion onto a stellar-mass black hole

Authors: Brianna S. Mills, Shane W. Davis, Yan-Fei Jiang, Matthew J. Middleton

Abstract: We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations prod… ▽ More We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations produced moderately super-Eddington accretion rates onto a 6.62 $M_\odot$ black hole. Since we only achieve inflow equilibrium out to 20-25 gravitational radii, we focus on the hard X-ray emission. We provide a comparison between the MC and RMHD simulations showing that the treatment of Compton scattering in the gray RMHD simulations underestimates the gas temperature in the regions above and below the accretion disk. In contrast, the restarted multi-group snapshots provides a treatment for the radiation field that is more consistent with the MC calculations, and result in post-processed spectra with harder X-ray emission compared to their gray snapshot counterparts. We characterize these MC post-processed spectra using commonly employed phenomenological models used for spectral fitting. We also attempt to fit our MC spectra directly to observations of the ultraluminous X-ray source (ULX) NGC 1313 X-1, finding best fit values that are competitive to phenomenological model fits, indicating that first principle models of super-Eddington accretion may adequately explain the observed hard X-ray spectra in some ULX sources. △ Less

Submitted 16 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ; 20 pages, 15 figures,

A Bright First Day for Tidal Disruption Event

Authors: Xiaoshan Huang, Shane W. Davis, Yan-fei Jiang

Abstract: Stream-stream collision may be an important pre-peak energy dissipation mechanism in tidal disruption events (TDEs). We perform local three-dimensional radiation hydrodynamic simulations in a wedge geometry including the gravity to study stream self-crossing, with emphasis on resolving the collision and following the subsequent outflow. We find that the collision can contribute to pre-peak optical… ▽ More Stream-stream collision may be an important pre-peak energy dissipation mechanism in tidal disruption events (TDEs). We perform local three-dimensional radiation hydrodynamic simulations in a wedge geometry including the gravity to study stream self-crossing, with emphasis on resolving the collision and following the subsequent outflow. We find that the collision can contribute to pre-peak optical emissions by converting $\gtrsim5\%$ of stream kinetic energy to radiation, yielding prompt emission of $\sim10^{42-44}\rm erg~s^{-1}$. The radiative efficiency is sensitive to stream mass fallback rates, and strongly depends on the downstream gas optical depth. Even for a sub-Eddington ($10\%$) mass fallback rate, the strong radiation pressure produced in the collision can form a local super-Eddington region near the collision site, where a fast, aspherical outflow is launched. Higher mass fallback rate usually leads to more optically-thick outflow and lower net radiative efficiency. For $\dot{M}\gtrsim0.1\dot{M}_{\rm Edd}$, the estimated photosphere size of the outflow can expand by one to two orders of magnitudes reaching $\sim10^{14}\rm cm$. The average gas temperature at this photospheric surface is a few $\times10^{4}$K, roughly consistent with inferred pre-peak photosphere properties for some optical TDEs. We find that the dynamics is sensitive to collision angle and collision radius, but the radiative efficiency or outflow properties show more complex dependency than is often assumed in ballistic models. △ Less

Submitted 6 July, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: 18 pages, 18 figures, Submitted to ApJ. Comments are welcomed and appreciated!

An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization

Authors: Christopher J. White, Patrick D. Mullen, Yan-Fei Jiang, Shane W. Davis, James M. Stone, Viktoriya Morozova, Lizhong Zhang

Abstract: We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical… ▽ More We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems. △ Less

Submitted 28 March, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: To be published in ApJS

Report number: LA-UR-23-21167

Global Three-Dimensional Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar Mass Black Hole at Sub- and Near-critical Accretion Rates

Authors: Jiahui Huang, Yan-Fei Jiang, Hua Feng, Shane W. Davis, James M. Stone, Matthew J. Middleton

Abstract: We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundr… ▽ More We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 $r_{\rm g}$. The energy dissipation happens close to the mid-plane in the near-critical runs and near the disk surface in the low accretion rate run. The total radiative luminosity inside $\sim$20 $r_{\rm g}$ is about 1% to 30% the Eddington limit, with a radiative efficiency of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability (MRI) leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel with a terminal velocity of $\sim$0.1$c$ are seen only in the near-critical runs. We conclude that these magnetic pressure dominated disks are thermally stable and thicker than the $α$ disk, and the effective temperature profiles are much flatter than that in the $α$ disks. The magnetic pressure of these disks are comparable within an order of magnitude with the previous analytical magnetic pressure dominated disk model. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 17 pages, 13 figures, 3 tables, accepted for publication in ApJ

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

Authors: Mihwa Jin, Ka Ho Lam, Melissa K. McClure, Jeroen Terwisscha van Scheltinga, Zhi-Yun Li, Adwin Boogert, Eric Herbst, Shane W. Davis, Robin T. Garrod

Abstract: 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-… ▽ More 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. △ Less

Submitted 9 July, 2022; originally announced July 2022.

Comments: Accepted to ApJ (28 Pages, 11 figures)

CR Driven Multi-phase Gas Formed via Thermal Instability

Authors: Xiaoshan Huang, Yan-fei Jiang, Shane W. Davis

Abstract: 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… ▽ More 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. △ Less

Submitted 11 June, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 22 pages, 18 figures

Journal ref: 2022 ApJ 931 140

The Launching of Cosmic Ray Driven Outflows

Authors: Xiaoshan Huang, Shane W. Davis

Abstract: 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 Ed… ▽ More 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. △ Less

Submitted 25 February, 2022; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: 18 pages, 19 figures, submitted to MNRAS

The black hole spin in GRS 1915+105, revisited

Authors: Brianna S. Mills, Shane W. Davis, Matthew J. Middleton

Abstract: We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion disk-dominated spectral states as described in McClintock et al. (2006) and Middleton e… ▽ More We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion disk-dominated spectral states as described in McClintock et al. (2006) and Middleton et al. (2006), which previously gave discrepant spin estimates with this method. We find that, using the new system parameters, the spin in both datasets increased, providing a best-fit spin of $a_*=0.86$ for the Middleton et al. data and a poor fit for the McClintock et al. dataset, which becomes pegged at the BHSPEC model limit of $a_*=0.99$. We explore the impact of the uncertainties in the system parameters, showing that the best-fit spin ranges from $a_*= 0.4$ to 0.99 for the Middleton et al. dataset and allows reasonable fits to the McClintock et al. dataset with near maximal spin for system distances greater than $\sim 10$ kpc. We discuss the uncertainties and implications of these estimates. △ Less

Submitted 11 July, 2024; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 14 pages, 5 figures, Accepted for publication in the Astrophysical Journal

Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets

Authors: Shane W. Davis, Alexander Tchekhovskoy

Abstract: 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 instabili… ▽ More 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. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 35 pages, 8 figures, Invited review published in Annual Reviews of Astronomy and Astrophysics

Time Dependent Radiation Hydrodynamics on a Moving Mesh

Authors: Philip Chang, Shane W. Davis, Yan-Fei Jiang

Abstract: 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… ▽ More 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. △ Less

Submitted 25 February, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: 13 pages, 8 figures, accepted to MNRAS

Covariant Radiative Transfer for Black Hole Spacetimes

Authors: Shane W. Davis, Charles F. Gammie

Abstract: 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 consisten… ▽ More 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. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 29 pages, 0 figures

Dusty Cloud Acceleration with Multiband Radiation

Authors: Xiaoshan Huang, Shane W. Davis, Dong Zhang

Abstract: 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 emiss… ▽ More 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. △ Less

Submitted 16 April, 2020; v1 submitted 5 August, 2019; originally announced August 2019.

Comments: 15 Pages, 14 Figures

Journal ref: ApJ 893 (2020) 50

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

Authors: Yan-Fei Jiang, Omer Blaes, James Stone, Shane W. Davis

Abstract: 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 dis… ▽ More 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. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 13 pages, 12 figures, submitted to ApJ

Spectral Hardening in Black Hole Accretion: Giving Spectral Modelers an f

Authors: Shane W. Davis, Samer El-Abd

Abstract: By fitting synthetic spectral models computed via the TLUSTY code, we examine how the spectra from thin accretion disks are expected to vary in accreting black hole systems. We fit color-corrected blackbody models to our synthetic spectra to estimate the spectral hardening factor f, which parameterizes the departure from blackbody and is commonly used to help interpret multitemperature blackbody f… ▽ More By fitting synthetic spectral models computed via the TLUSTY code, we examine how the spectra from thin accretion disks are expected to vary in accreting black hole systems. We fit color-corrected blackbody models to our synthetic spectra to estimate the spectral hardening factor f, which parameterizes the departure from blackbody and is commonly used to help interpret multitemperature blackbody fitting results. We find we can define a reasonably robust f value to spectra when the effects of Compton scattering dominate radiation transfer. We examine the evolution of f with black hole mass and accretion rate, typically finding a moderate variation (f ~ 1.4-2) for accretion rates between 1% and 100% of the Eddington rate. Consistent with most previous work, we find f tends to increase with accretion rate, but we also infer a weaker correlation of f with black holes mass. We find that f is rarely much larger than 2 unless the disk becomes photon starved, in contention with some previous calculations. Significant spectral hardening (f > 2) is only found when the disk mass surface density is lower than expected for alpha-disk models unless alpha is near unity or larger. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 9 pages, submitted to ApJ

Super-Eddington Accretion Disks around Supermassive black Holes

Authors: Yan-Fei Jiang, James Stone, Shane W. Davis

Abstract: We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases wi… ▽ More We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $\sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $\sim 0.1-0.4c$. When the accretion rate is smaller than $\sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $\sim 5\%-7\%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1\%$. We always find the kinetic energy luminosity associated with the outflow is only $\sim 15\%-30\%$ of the radiative luminosity. The mass flux lost in the outflow is $\sim 15\%-50\%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks. △ Less

Submitted 8 September, 2017; originally announced September 2017.

Comments: 21 pages, 17 figures, submitted to ApJ, movies available here https://goo.gl/5yR3ig

Dusty Cloud Acceleration by Radiation Pressure in Rapidly Star-Forming Galaxies

Authors: Dong Zhang, Shane W. Davis, Yan-Fei Jiang, James M. Stone

Abstract: We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable o… ▽ More We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable of accelerating the clouds to hundreds of kilometers per second while remaining dense and cold, consistent with observations. We compare these results to simulations where acceleration is provided by entrainment in a hot wind, where the momentum injection of the hot flow is comparable to the momentum in the radiation field. We find that the survival time of the cloud accelerated by the radiation field is significantly longer than that of a cloud entrained in a hot outflow. We show that the dynamics of the irradiated cloud depends on the initial optical depth, temperature of the cloud, and the intensity of the flux. Additionally, gas pressure from the background may limit cloud acceleration if the density ratio between the cloud and background is $\lesssim 10^{2}$. In general, a 10 pc-scale optically thin cloud forms a pancake structure elongated perpendicular to the direction of motion, while optically thick clouds form a filamentary structure elongated parallel to the direction of motion. The details of accelerated cloud morphology and geometry can also be affected by other factors, such as the cloud lengthscale, the reduced speed of light approximation, spatial resolution, initial cloud structure, and the dimensionality of the run, but these have relatively little affect on the cloud velocity or survival time. △ Less

Submitted 16 January, 2018; v1 submitted 9 August, 2017; originally announced August 2017.

Comments: 22 pages, 21 figures, 2 table, Accepted for publication in ApJ

Radiation Hydrodynamic Simulations of Dust-Driven Winds

Authors: Dong Zhang, Shane W. Davis

Abstract: We study dusty winds driven by radiation pressure in the atmosphere of a rapidly star-forming environment. We apply the variable Eddington tensor algorithm to re-examine the two-dimensional radiation hydrodynamic problem of a column of gas that is accelerated by a constant infrared radiation flux. In the absence of gravity, the system is primarily characterized by the initial optical depth of the… ▽ More We study dusty winds driven by radiation pressure in the atmosphere of a rapidly star-forming environment. We apply the variable Eddington tensor algorithm to re-examine the two-dimensional radiation hydrodynamic problem of a column of gas that is accelerated by a constant infrared radiation flux. In the absence of gravity, the system is primarily characterized by the initial optical depth of the gas. We perform several runs with different initial optical depth and resolution. We find that the gas spreads out along the vertical direction, as its mean velocity and velocity dispersion increase. In contrast to previous work using flux-limited diffusion algorithm, we find little evolution in the trapping factor. The momentum coupling between radiation and gas in the absence of gravity is similar to that with gravity. For Eddington ratio increasing with the height in the system, the momentum transfer from the radiation to the gas is not merely $\sim L/c$, but amplified by a factor of $1+ητ_{\rm IR}$, where $τ_{\rm IR}$ is the integrated infrared optical depth through the system, and $η\sim0.5-0.9$, decreasing with the optical depth. We apply our results to the atmosphere of galaxies and conclude that radiation pressure may be an important mechanism for driving winds in the most rapidly star-forming galaxies and starbursts. △ Less

Submitted 23 March, 2017; v1 submitted 30 November, 2016; originally announced December 2016.

Comments: 12 pages, 10 figures, 1 table, Accepted for publication in ApJ

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

Authors: Murray Brightman, Fiona A. Harrison, Didier Barret, Shane W. Davis, Felix Fürst, Kristin K. Madsen, Matthew Middleton, Jon M. Miller, Daniel Stern, Lian Tao, Dominic J. Walton

Abstract: 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 p… ▽ More 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. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: Accepted for publication in ApJ

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

Authors: Yan-Fei Jiang, James M. Stone, Shane W. Davis

Abstract: 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 tra… ▽ More 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. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: 15 pages, 16 figures, 1 table, accepted to ApJ

Journal ref: 2014, ApJ, 796, 106

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

Authors: Yan-Fei Jiang, James M. Stone, Shane W. Davis

Abstract: (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… ▽ More (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. △ Less

Submitted 24 March, 2014; originally announced March 2014.

Comments: 15 pages, 16 figures, submitted to ApJS

Journal ref: ApJS, 213, 7, 2014

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

Authors: Shane W. Davis, Yan-Fei Jiang, James M. Stone, Norman Murray

Abstract: 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… ▽ More 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. △ Less

Submitted 7 March, 2014; originally announced March 2014.

Comments: 16 pages, 11 figures, submitted to ApJ

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

Authors: Yan-Fei Jiang, James M. Stone, Shane W. Davis

Abstract: 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 magn… ▽ More 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. △ Less

Submitted 12 February, 2014; originally announced February 2014.

Comments: 7 pages, 5 figures, accepted by ApJ

Line driven winds and the UV turnover in AGN accretion discs

Authors: Ari Laor, Shane W. Davis

Abstract: 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 t… ▽ More 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. △ Less

Submitted 29 December, 2013; v1 submitted 12 December, 2013; originally announced December 2013.

Comments: Accepted for publication in MNRAS. Minor typos corrected

Dynamics of warped accretion discs

Authors: Scott Tremaine, Shane W. Davis

Abstract: 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-momen… ▽ More 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. △ Less

Submitted 8 August, 2013; originally announced August 2013.

Comments: 28 pages, 11 figures

Saturation of the MRI in Strongly Radiation Dominated Accretion Disks

Authors: Yan-Fei Jiang, James M. Stone, Shane W. Davis

Abstract: 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 ver… ▽ More 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. △ Less

Submitted 7 March, 2013; originally announced March 2013.

Comments: 15 pages, 8 figures, 4 tables, accepted by ApJ

Journal ref: 2013, ApJ 767 148

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

Authors: Yan-Fei Jiang, Shane W. Davis, James Stone

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2012; originally announced December 2012.

Comments: 15 pages, 1 table, 13 figures, accepted by ApJ

Journal ref: 2013, ApJ 763, 102

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

Authors: Yucong Zhu, Shane W. Davis, Ramesh Narayan, Akshay K. Kulkarni, Robert F. Penna, Jeffrey E. McClintock

Abstract: 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 si… ▽ More 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. △ Less

Submitted 25 April, 2012; v1 submitted 7 February, 2012; originally announced February 2012.

Comments: 21 pages, 15 figures, accepted for publication in MNRAS

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

Authors: Yan-Fei Jiang, James M. Stone, Shane W. Davis

Abstract: 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 e… ▽ More 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. △ Less

Submitted 10 January, 2012; originally announced January 2012.

Comments: 42 pages, 22 figures, 2 tables, accepted by ApJS

Journal ref: ApJS, 2012, 199, 14

A Radiation Transfer Solver for Athena using Short Characteristics

Authors: Shane W. Davis, James M. Stone, Yan-Fei Jiang

Abstract: 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… ▽ More 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) △ Less

Submitted 10 January, 2012; originally announced January 2012.

Comments: 20 pages, 13 figures, accepted for publication in ApJ Supplement Series

Cold Accretion Disks and Lineless Quasars

Authors: Ari Laor, Shane W. Davis

Abstract: 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_… ▽ More 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. △ Less

Submitted 24 June, 2011; originally announced June 2011.

Comments: Accepted for publication in MNRAS

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

Authors: Lijun Gou, Jeffrey E. McClintock, Mark J. Reid, Jerome A. Orosz, James F. Steiner, Ramesh Narayan, Jingen Xiang, Ronald A. Remillard, Keith A. Arnaud, Shane W. Davis

Abstract: 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'… ▽ More 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. △ Less

Submitted 30 September, 2011; v1 submitted 18 June, 2011; originally announced June 2011.

Comments: Paper III of three papers on Cygnus X-1; 21 pages including 5 figures and 12 tables, ApJ in press. The paper is significantly restructured; two further tests of the robustness of our spin measurement are presented, and our error analysis has been substantially improved; the conclusions are unchanged

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

Authors: Shane W. Davis, Ramesh Narayan, Yucong Zhu, Didier Barret, Sean A. Farrell, Olivier Godet, Mathieu Servillat, Natalie A. Webb

Abstract: 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… ▽ More 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. △ Less

Submitted 13 April, 2011; originally announced April 2011.

Comments: 10 pages, 7 figures, accepted for publication in ApJ

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

Authors: Akshay K. Kulkarni, Robert F. Penna, Roman V. Shcherbakov, James F. Steiner, Ramesh Narayan, Aleksander Sadowski, Yucong Zhu, Jeffrey E. McClintock, Shane W. Davis, Jonathan C. McKinney

Abstract: 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… ▽ More 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. △ Less

Submitted 2 February, 2011; v1 submitted 31 January, 2011; originally announced February 2011.

Comments: 13 pages, 7 figures, accepted for publication in MNRAS. v2: fixed typo in author name, updated acknowledgments

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

Authors: Rodrigo Fernández, Shane W. Davis

Abstract: 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 magnetosp… ▽ More 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. △ Less

Submitted 3 February, 2011; v1 submitted 4 January, 2011; originally announced January 2011.

Comments: Accepted for publication in ApJ. Added figure (number 1) showing the geometry of the problem. Minor changes otherwise

Measuring the Spins of Accreting Black Holes

Authors: Jeffrey E. McClintock, Ramesh Narayan, Shane W. Davis, Lijun Gou, Akshay Kulkarni, Jerome A. Orosz, Robert F. Penna, Ronald A. Remillard, James F. Steiner

Abstract: 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 measurement… ▽ More 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. △ Less

Submitted 29 January, 2011; v1 submitted 4 January, 2011; originally announced January 2011.

Comments: To appear in Classical and Quantum Gravity; Special volume for GR19, eds. D. Marolf and D. Sudarsky; 28 pages, 6 figures; Includes corrections made to proofs, which are significant only for Section 5

The Radiative Efficiency of Accretion Flows in Individual AGN

Authors: Shane W. Davis, Ari Laor

Abstract: 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… ▽ More 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. △ Less

Submitted 14 December, 2010; originally announced December 2010.

Comments: 20 pages, accepted for publication in The Astrophysical Journal

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

Authors: Aya Kubota, Chris Done, Shane W. Davis, Tadayasu Dotani, Tsunefumi Mizuno, Yoshihiro Ueda

Abstract: 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 clearl… ▽ More 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. △ Less

Submitted 17 March, 2010; originally announced March 2010.

Comments: 12pages, 7 figures, Accepted by ApJ

Sustained Magnetorotational Turbulence in Local Simulations of Stratified Disks with Zero Net Magnetic Flux

Authors: S. W. Davis, J. M. Stone, M. E. Pessah

Abstract: We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, cont… ▽ More We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles. △ Less

Submitted 8 September, 2009; originally announced September 2009.

Comments: 13 pages, 18 figures, submitted to ApJ

Journal ref: Astrophys.J.713:52-65,2010

The Effects of Magnetic Fields and Inhomogeneities on Accretion Disk Spectra and Polarization

Authors: S. W. Davis, O. M. Blaes, S. Hirose, J. H. Krolik

Abstract: We present the results of one and three-dimensional radiative transfer calculations of polarized spectra emerging from snapshots of radiation magnetohydrodynamical simulations of the local vertical structure of black hole accretion disks. The simulations cover a wide range of physical regimes relevant for the high/soft state of black hole X-ray binaries. We constrain the uncertainties in theoret… ▽ More We present the results of one and three-dimensional radiative transfer calculations of polarized spectra emerging from snapshots of radiation magnetohydrodynamical simulations of the local vertical structure of black hole accretion disks. The simulations cover a wide range of physical regimes relevant for the high/soft state of black hole X-ray binaries. We constrain the uncertainties in theoretical spectral color correction factors due to the presence of magnetic support of the disk surface layers and strong density inhomogeneities. For the radiation dominated simulation, magnetic support increases the color correction factor by about ten percent, but this is largely compensated by a ten percent softening due to inhomogeneities. We also compute the effects of inhomogeneities and Faraday rotation on the resulting polarization. Magnetic fields in the simulations are just strong enough to produce significant Faraday depolarization near the spectral peak of the radiation field. X-ray polarimetry may therefore be a valuable diagnostic of accretion disk magnetic fields, being able to directly test simulations of magnetorotational turbulence. △ Less

Submitted 4 August, 2009; originally announced August 2009.

Comments: 18 pages, accepted for publication in ApJ

A Determination of the Spin of the Black Hole Primary in LMC X-1

Authors: Lijun Gou, Jeffrey E. McClintock, Jifeng Liu, Ramesh Narayan, James F. Steiner, Ronald A. Remillard, Jerome A. Orosz, Shane W. Davis, Ken Ebisawa, Eric M. Schlegel

Abstract: The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In the 1980s, its compact primary was established as the fourth dynamical black-hole candidate. Recently, we published accurate values for the mass of the black hole and the orbital inclination angle of the binary system. Building on these results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a selected sampl… ▽ More The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In the 1980s, its compact primary was established as the fourth dynamical black-hole candidate. Recently, we published accurate values for the mass of the black hole and the orbital inclination angle of the binary system. Building on these results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a selected sample of 18 of these spectra, we have determined the dimensionless spin parameter of the black hole to be a* = 0.92(-0.07,+0.05). This result takes into account all sources of observational and model-parameter uncertainties. The standard deviation around the mean value of a* for these 18 X-ray spectra, which were obtained over a span of several years, is only 0.02. When we consider our complete sample of 53 RXTE spectra, we find a somewhat higher value of the spin parameter and a larger standard deviation. Finally, we show that our results based on RXTE data are confirmed by our analyses of selected X-ray spectra obtained by the XMM-Newton, BeppoSAX and Ginga missions. △ Less

Submitted 29 June, 2009; v1 submitted 8 January, 2009; originally announced January 2009.

Comments: 40 pages including 9 figures and 4 tables, ApJ in press

Journal ref: Astrophys.J.701:1076-1090,2009

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

Authors: Jifeng Liu, Jeffery E. McClintock, Ramesh Narayan, Shane W. Davis, Jerome A. Orosz

Abstract: 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… ▽ More 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. △ Less

Submitted 11 May, 2010; v1 submitted 12 March, 2008; originally announced March 2008.

Comments: 14 pages, 3 figures, 1 table, published in ApJ Letters; as explained in the erratum at the end of the text, the spin parameter has been corrected upward from a*=0.77 to a*=0.84. Apart from the addition of the erratum, the paper is unchanged.

Journal ref: The Astrophysical Journal, 679, L37-L40, 2008 May 20

Angular Momentum Transport in Accretion Disks and its Implications for Spin Estimates in Black Hole Binaries

Authors: Chris Done, Shane W. Davis

Abstract: The accretion flow in the disk dominated state of black hole binaries has peak temperature and luminosity which vary together in such a way as to indicate an approximately constant emitting area. The association of this with the last stable orbit gives one of the few ways to estimate spin when the mass of the black hole is known. However, deriving this radius requires knowledge of how the disk s… ▽ More The accretion flow in the disk dominated state of black hole binaries has peak temperature and luminosity which vary together in such a way as to indicate an approximately constant emitting area. The association of this with the last stable orbit gives one of the few ways to estimate spin when the mass of the black hole is known. However, deriving this radius requires knowledge of how the disk spectrum is modified by radiative transfer through the vertical structure of the disk, as well as special and general relativistic effects on the propagation of this radiation. Here we investigate the extent to which differences in vertical structure change the derived disk spectra by calculating these for a range of different stress prescriptions. We find that at a given mass accretion rate the spectra are almost identical for accretion rates of L/L_Edd <~ 0.1. The spectra are remarkably similar even up to the highest luminosities considered (L/L_Edd ~ 0.6) as long as the stresses do not dissipate more than about 10 per cent of the gravitational energy above the effective photosphere. This is exceeded only by classic alpha disks with alpha >~ 0.1, but these models give spectral variation which is incompatible with existing data. Therefore, we conclude that disk spectral modelling can place interesting constraints on angular momentum transport, but still provide a robust estimate of the spin of the black hole. △ Less

Submitted 6 June, 2008; v1 submitted 4 March, 2008; originally announced March 2008.

Comments: accepted by ApJ

The UV Continuum of Quasars: Models and SDSS Spectral Slopes

Authors: Shane W. Davis, Jong-Hak Woo, Omer M. Blaes

Abstract: We measure long (2200-4000 ang) and short (1450-2200 ang) wavelength spectral slopes α(F_νproportional to ν^α) for quasar spectra from the Sloan Digital Sky Survey. The long and short wavelength slopes are computed from 3646 and 2706 quasars with redshifts in the z=0.76-1.26 and z=1.67-2.07 ranges, respectively. We calculate mean slopes after binning the data by monochromatic luminosity at 2200… ▽ More We measure long (2200-4000 ang) and short (1450-2200 ang) wavelength spectral slopes α(F_νproportional to ν^α) for quasar spectra from the Sloan Digital Sky Survey. The long and short wavelength slopes are computed from 3646 and 2706 quasars with redshifts in the z=0.76-1.26 and z=1.67-2.07 ranges, respectively. We calculate mean slopes after binning the data by monochromatic luminosity at 2200 ang and virial mass estimates based on measurements of the MgII line width and 3000 ang continuum luminosity. We find little evidence for mass dependent variations in the mean slopes, but a significant luminosity dependent trend in the near UV spectral slopes is observed with larger (bluer) slopes at higher luminosities. The far UV slopes show no clear variation with luminosity and are generally lower (redder) than the near UV slopes at comparable luminosities, suggesting a slightly concave quasar continuum shape. We compare these results with Monte Carlo distributions of slopes computed from models of thin accretion disks, accounting for uncertainties in the mass estimates. The model slopes produce mass dependent trends which are larger than observed, though this conclusion is sensitive to the assumed uncertainties in the mass estimates. The model slopes are also generally bluer than observed, and we argue that reddening by dust intrinsic to the source or host galaxy may account for much of the discrepancy. △ Less

Submitted 10 July, 2007; originally announced July 2007.

Comments: To be published in ApJ, 18 pages, 10 figures

Turbulent Comptonization in Relativistic Accretion Disks

Authors: Aristotle Socrates, Shane W. Davis, Omer Blaes

Abstract: Turbulent Comptonization, a potentially important damping and radiation mechanism in relativistic accretion flows, is discussed. Particular emphasis is placed on the physical basis, relative importance, and thermodynamics of turbulent Comptonization. The effects of metal-absorption opacity on the spectral component resulting from turbulent Comptonization is considered as well. Turbulent Comptonization, a potentially important damping and radiation mechanism in relativistic accretion flows, is discussed. Particular emphasis is placed on the physical basis, relative importance, and thermodynamics of turbulent Comptonization. The effects of metal-absorption opacity on the spectral component resulting from turbulent Comptonization is considered as well. △ Less

Submitted 28 September, 2006; originally announced September 2006.

Comments: 4 page COSPAR conference proceeding. New discussion on photo-ionization balance included

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

Authors: Aristotle Socrates, Shane W. Davis, Enrico Ramirez-Ruiz

Abstract: 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… ▽ More 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. △ Less

Submitted 22 January, 2008; v1 submitted 28 September, 2006; originally announced September 2006.

Comments: 15 pages, accepted for publication in the ApJ

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

Authors: Jeffrey E. McClintock, Rebecca Shafee, Ramesh Narayan, Ronald A. Remillard, Shane W. Davis, Li-Xin Li

Abstract: 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 analysi… ▽ More 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. △ Less

Submitted 16 August, 2006; v1 submitted 4 June, 2006; originally announced June 2006.

Comments: 58 pages, 18 figures. Accepted for publication in ApJ. New in this version is a proposed observational test of our spin model and the kinematic model of the radio jets

Journal ref: Astrophys.J.652:518-539,2006