Showing 1–50 of 88 results for author: Blaes, O

Dynamical Effects of Magnetic Opacity in Neutron Star Accretion Columns

Authors: Xin Sheng, Lizhong Zhang, Omer Blaes, Yan-Fei Jiang

Abstract: We present relativistic, radiation magnetohydrodynamic simulations of supercritical neutron star accretion columns in Cartesian geometry, including temperature-dependent, polarization-averaged Rosseland mean opacities accounting for classical electron scattering in a magnetic field. Just as in our previous pure Thomson scattering simulations, vertical oscillations of the accretion shock and horizo… ▽ More We present relativistic, radiation magnetohydrodynamic simulations of supercritical neutron star accretion columns in Cartesian geometry, including temperature-dependent, polarization-averaged Rosseland mean opacities accounting for classical electron scattering in a magnetic field. Just as in our previous pure Thomson scattering simulations, vertical oscillations of the accretion shock and horizontally propagating entropy waves (photon bubbles) are present in all our simulations. However, at high magnetic fields $\gtrsim10^{12}$~G, the magnetic opacities produce significant differences in the overall structure and dynamics of the column. At fixed accretion rate, increasing the magnetic field strength results in a shorter accretion column, despite the fact that the overall opacity within the column is larger. Moreover, the vertical oscillation amplitude of the column is reduced. Increasing the accretion rate at high magnetic fields restores the height of the column. However, a new, slower instability takes place at these field strengths because they are in a regime where the opacity increases with temperature. This instability causes both the average height of the column and the oscillation amplitude to substantially increase on a time scale of $\sim10$~ms. We provide physical explanations for these results, and discuss their implications for the observed properties of these columns, including mixed fan-beam/pencil-beam emission patterns caused by the oscillations. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS

Dynamics of Neutron Star Accretion Columns in Split-Monopole Magnetic Fields

Authors: Lizhong Zhang, Omer Blaes, Yan-Fei Jiang

Abstract: We perform 2D axisymmetric radiative relativistic MHD simulations of radiation pressure supported neutron star accretion columns in split-monopole magnetic fields. The accretion columns exhibit quasi-periodic oscillations, which manifest in the luminosity power spectrum as 2-10 kHz peaks, together with broader extensions to somewhat higher frequencies. The peak frequency decreases for wider column… ▽ More We perform 2D axisymmetric radiative relativistic MHD simulations of radiation pressure supported neutron star accretion columns in split-monopole magnetic fields. The accretion columns exhibit quasi-periodic oscillations, which manifest in the luminosity power spectrum as 2-10 kHz peaks, together with broader extensions to somewhat higher frequencies. The peak frequency decreases for wider columns or higher mass accretion rates. In contrast to the case of shorter columns in uniform magnetic fields, pdV work contributes substantially to maintaining the radiation pressure inside the column against sideways radiative cooling. This is in part due to the compression associated with accretion along the converging magnetic field lines towards the stellar surface. Propagating entropy waves which are associated with the slow-diffusion photon bubble instability form in all our simulations. Radial advection of radiation from the oscillation itself as well as the entropy waves is also important in maintaining radiation pressure inside the column. The time-averaged profile of our fiducial simulation accretion is approximately consistent with the classical 1D stationary model provided one incorporates the correct column shape. We also quantify the porosity in all our accretion column simulations so that this may also in principle be used to improve 1D models. △ Less

Submitted 3 January, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: 18 pages, 17 figures, and 3 table, accepted by MNRAS

Radiative Relativistic Magnetohydrodynamic Simulations of Neutron Star Column Accretion in Cartesian Geometry

Authors: Lizhong Zhang, Omer Blaes, Yan-Fei Jiang

Abstract: High luminosity accretion onto a strongly magnetized neutron star results in a radiation pressure dominated, magnetically confined accretion column. We investigate the dynamics of these columns using two-dimensional radiative relativistic magnetohydrodynamic simulations, restricting consideration to modest accretion rates where the height of the column is low enough that Cartesian geometry can be… ▽ More High luminosity accretion onto a strongly magnetized neutron star results in a radiation pressure dominated, magnetically confined accretion column. We investigate the dynamics of these columns using two-dimensional radiative relativistic magnetohydrodynamic simulations, restricting consideration to modest accretion rates where the height of the column is low enough that Cartesian geometry can be employed. The column structure is dynamically maintained through high-frequency oscillations of the accretion shock at $\simeq 10-25$~kHz. These oscillations arise because it is necessary to redistribute the power released at the accretion shock through bulk vertical motions, both to balance the cooling and to provide vertical pressure support against gravity. Sideways cooling always dominates the loss of internal energy. In addition to the vertical oscillations, photon bubbles form in our simulations and add additional spatial complexity to the column structure. They are not themselves responsible for the oscillations, and they do not appear to affect the oscillation period. However, they enhance the vertical transport of radiation and increase the oscillation amplitude in luminosity. The time-averaged column structure in our simulations resembles the trends in standard 1D stationary models, the main difference being that the time-averaged height of the shock front is lower because of the higher cooling efficiency of the 2D column shape. △ Less

Submitted 13 August, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 21 pages, 19 figures, and 1 table, published in MNRAS

Radiative MHD Simulations of Photon Bubbles in Radiation-Supported Magnetized Atmospheres of Neutron Stars with Isotropic Thomson Scattering

Authors: Lizhong Zhang, Omer Blaes, Yan-Fei Jiang

Abstract: A major uncertainty in the structure and dynamics of magnetized, radiation pressure dominated neutron star accretion columns in X-ray pulsars and pulsating ultraluminous X-ray sources is that they are thought to be subject to the photon bubble instability. We present the results of two dimensional radiation relativistic magnetohydrodynamic simulations of a non-accreting, static atmosphere to study… ▽ More A major uncertainty in the structure and dynamics of magnetized, radiation pressure dominated neutron star accretion columns in X-ray pulsars and pulsating ultraluminous X-ray sources is that they are thought to be subject to the photon bubble instability. We present the results of two dimensional radiation relativistic magnetohydrodynamic simulations of a non-accreting, static atmosphere to study the development of this instability assuming isotropic Thomson scattering in the slow diffusion regime that is relevant to neutron star accretion columns. Photon bubbles generally grow faster toward shorter wavelengths, until a maximum growth rate is achieved at the radiation viscosity length scale, which is generally quite small and requires high numerical resolution to simulate. We confirm the consistency between our simulation results and linear theory in detail, and show that the nonlinear evolution inevitably leads to collapse of the atmosphere with the higher resolution simulation collapsing faster due to the presence of shorter length scale nonlinear structures. At least in static atmospheres with horizontally periodic boundary conditions, this resolution dependence may make simulations of the nonlinear dynamics of photon bubble instability in neutron star accretion columns challenging. It remains to be seen whether these difficulties will persist upon inclusion of an accretion flow through the top and magnetically-confined horizontal boundaries through which photons can escape. Our results here provide a foundation for such future work. △ Less

Submitted 6 October, 2021; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 21 pages, 18 figures and 3 tables, Published in MNRAS

A characteristic optical variability timescale in astrophysical accretion disks

Authors: Colin J. Burke, Yue Shen, Omer Blaes, Charles F. Gammie, Keith Horne, Yan-Fei Jiang, Xin Liu, Ian M. McHardy, Christopher W. Morgan, Simone Scaringi, Qian Yang

Abstract: Accretion disks around supermassive black holes in active galactic nuclei produce continuum radiation at ultraviolet and optical wavelengths. Physical processes in the accretion flow lead to stochastic variability of this emission on a wide range of timescales. We measure the optical continuum variability observed in 67 active galactic nuclei and the characteristic timescale at which the variabili… ▽ More Accretion disks around supermassive black holes in active galactic nuclei produce continuum radiation at ultraviolet and optical wavelengths. Physical processes in the accretion flow lead to stochastic variability of this emission on a wide range of timescales. We measure the optical continuum variability observed in 67 active galactic nuclei and the characteristic timescale at which the variability power spectrum flattens. We find a correlation between this timescale and the black hole mass, extending over the entire mass range of supermassive black holes. This timescale is consistent with the expected thermal timescale at the ultraviolet-emitting radius in standard accretion disk theory. Accreting white dwarfs lie close to this correlation, suggesting a common process for all accretion disks. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Author's version; 49 pages; Science (published Aug 13)

Journal ref: Science, Vol. 373, Issue 6556, pp. 789-792 (2021)

Investigating lack of accretion disk eccentricity growth in a global 3D MHD simulation of a superhump system

Authors: Bryance Oyang, Yan-Fei Jiang, Omer Blaes

Abstract: We present the results of a 3D global magnetohydrodynamic (MHD) simulation of an AM CVn system that was aimed at exploring eccentricity growth in the accretion disc self-consistently from a first principles treatment of the MHD turbulence. No significant eccentricity growth occurs in the simulation. In order to investigate the reasons why, we ran 2D alpha disc simulations with alpha values of 0.01… ▽ More We present the results of a 3D global magnetohydrodynamic (MHD) simulation of an AM CVn system that was aimed at exploring eccentricity growth in the accretion disc self-consistently from a first principles treatment of the MHD turbulence. No significant eccentricity growth occurs in the simulation. In order to investigate the reasons why, we ran 2D alpha disc simulations with alpha values of 0.01, 0.1, and 0.2, and found that only the latter two exhibit significant eccentricity growth. We present an equation expressing global eccentricity evolution in terms of contributing forces and use it to analyze the simulations. As expected, we find that the dominant term contributing to the growth of eccentricity is the tidal gravity of the companion star. In the 2D simulations, the alpha viscosity directly contributes to eccentricity growth. In contrast, the overall magnetic forces in the 3D simulation damp eccentricity. We also analyzed the mode-coupling mechanism of Lubow, and confirmed that the spiral wave excited by the 3:1 resonance was the dominant contributor to eccentricity growth in the 2D $α=0.1$ simulations, but other waves also contribute significantly. We found that the $α=0.1$ and 0.2 simulations had more relative mass at larger radii compared to the $α=0.01$ and 3D MHD simulation, which also had an effective $α$ of 0.01. This suggests that in 3D MHD simulations without sufficient poloidal magnetic flux, MRI turbulence does not saturate at a high enough $α$ to spread the disc to large enough radii to reproduce the superhumps observed in real systems. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 17 pages, 17 figures (accepted by MNRAS)

Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

Authors: Sean M. Ressler, Eliot Quataert, Christopher J. White, Omer Blaes

Abstract: We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulenc… ▽ More We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: Accepted by MNRAS. For animations see [https://smressle.bitbucket.io/animations.html]

Periodicities in the K2 lightcurve of HP Librae

Authors: Siddhant Solanki, Thomas Kupfer, Omer Blaes, Elmé Breedt, Simone Scaringi

Abstract: We analyse Kepler/K2 lightcurve data of the AM CVn system HP Librae. We detect with confidence four photometric periodicities in the system: the orbital frequency, both positive and negative superhumps, and the positive apsidal precession frequency of the accretion disc. This is only the second time that the apsidal precession frequency has ever been directly detected in the photometry of a helium… ▽ More We analyse Kepler/K2 lightcurve data of the AM CVn system HP Librae. We detect with confidence four photometric periodicities in the system: the orbital frequency, both positive and negative superhumps, and the positive apsidal precession frequency of the accretion disc. This is only the second time that the apsidal precession frequency has ever been directly detected in the photometry of a helium accreting system, after SDSS J135154.46-064309.0. We present phase-folded lightcurves and sliding power spectra of each of the four periodicities. We measure rates of change of the positive superhump period of $10^{-7}$ days/day. We also redetect a QPO at 300 cyc/day, a feature that has been stable over decades, and show that it is harmonically related to two other QPOs, the lowest of which is centered on the superhump/orbital frequency. The continuum power spectrum is consistent with a single power law with no evidence of any breaks within our observed frequency range. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted to MNRAS. The movie of the lightcurve is available at https://www.youtube.com/watch?v=fseDT2Mr-bM

Opacity Driven Convection and Variability in Accretion Disks around Supermassive Black Holes

Authors: Yan-Fei Jiang, Omer Blaes

Abstract: We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes th… ▽ More We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes the disk to be convective unstable. Turbulence generated by convection puffs up the disk due to additional turbulent pressure support and enhances the local angular momentum transport. This also results in strong fluctuations in surface density and heating of the disk. The opacity drops with increasing temperature and convection is suppressed. The disk cools down and the whole process repeats again. This causes strong oscillations of the disk scale height and luminosity variations by more than a factor of $\approx 3-6$ over a few years' timescale. Since the iron opacity bump will move to different locations of the disk for black holes with different masses and accretion rates, we suggest that this is a physical mechanism that can explain the variability of AGN with a wide range of amplitudes over a time scale of years to decades. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 15 pages, 14 figures, resubmitted after including the referee's comments

Journal ref: ApJ, 900, 25, 2020

The Effects of Tilt on the Images of Black Hole Accretion Flows

Authors: Christopher J. White, Jason Dexter, Omer Blaes, Eliot Quataert

Abstract: We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed… ▽ More We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets. △ Less

Submitted 30 March, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

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

The Comptonising medium of the neutron star in 4U 1636-53 through its lower kilohertz Quasi-Periodic Oscillations

Authors: Konstantinos Karpouzas, Mariano Méndez, Evandro M. Ribeiro, Diego Altamirano, Omer Blaes, Federico García

Abstract: Inverse Compton scattering dominates the high energy part of the spectra of neutron star (NS) low mass X-ray binaries (LMXBs). It has been proposed that inverse Compton scattering also drives the radiative properties of kilohertz quasi periodic oscillations (kHz QPOs). In this work, we construct a model that predicts the energy dependence of the rms amplitude and time lag of the kHz QPOs. Using th… ▽ More Inverse Compton scattering dominates the high energy part of the spectra of neutron star (NS) low mass X-ray binaries (LMXBs). It has been proposed that inverse Compton scattering also drives the radiative properties of kilohertz quasi periodic oscillations (kHz QPOs). In this work, we construct a model that predicts the energy dependence of the rms amplitude and time lag of the kHz QPOs. Using this model, we fit the rms amplitude and time lag energy spectra of the lower kHz QPO in the NS LMXB 4U 1636-53 over 11 frequency intervals of the QPO and report three important findings: (i) A medium that extends 1-8 km above the NS surface is required to fit the data; this medium can be sustained by the balance between gravity and radiation pressure, without forcing any equilibrium condition. (ii) We predict a time delay between the oscillating NS temperature, due to feedback, and the oscillating electron temperature of the medium which, with the help of phase resolved spectroscopy, can be used as a probe of the geometry and the feedback mechanism. (iii) We show that the observed variability as a function of QPO frequency is mainly driven by the oscillating electron temperature of the medium. This provides strong evidence that the Comptonising medium in LMXBs significantly affects, if not completely drives, the radiative properties of the lower kHz QPOs regardless of the nature of the dynamical mechanism that produces the QPO frequencies. △ Less

Submitted 14 December, 2019; v1 submitted 11 December, 2019; originally announced December 2019.

Comments: 20 pages, 9 figures

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

Tilted Disks around Black Holes: A Numerical Parameter Survey for Spin and Inclination Angle

Authors: Christopher J. White, Eliot Quataert, Omer Blaes

Abstract: We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quick… ▽ More We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quickly reach a warped and twisted shape that rigidly precesses about the black hole spin axis with deformations in shape large enough to hinder the application of linear bending wave theory. Magnetized polar outflows form, oriented along the disk rotation axes. At sufficiently high inclinations a pair of standing shocks develops in the disks. These shocks dramatically affect the flow at small radii, driving angular momentum transport. At high spins they redirect material more effectively than they heat it, reducing the dissipation rate relative to the mass accretion rate and lowering the radiative efficiency of the flow. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 20 pages, 14 figures, to be published in ApJ

Journal ref: ApJ 878 (2019) 51

Three Dimensional Radiation Hydrodynamic Simulations of Massive Star Envelopes

Authors: Yan-Fei Jiang, Matteo Cantiello, Lars Bildsten, Eliot Quataert, Omer Blaes, James Stone

Abstract: (Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some… ▽ More (Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: This is the authors' original version. The paper is published in the September 27th issue of Nature: https://rdcu.be/7PgW

Journal ref: Nature 2018, 561, 498

A simple framework for modelling the dependence of bulk Comptonization by turbulence on accretion disc parameters

Authors: Jason Kaufman, Omer M. Blaes, Shigenobu Hirose

Abstract: Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated di… ▽ More Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated discs, in these systems bulk Comptonization by turbulence may provide a physical basis in the disc itself for warm Comptonization models. We model the dependence of bulk Comptonization on fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and $α$. In addition to constraining warm Comptonization models, our model can help distinguish contributions from bulk Comptonization to the soft X-ray excess from those due to other physical mechanisms, such as absorption and reflection. By linking the time variability of bulk Comptonization to fluctuations in the disc vertical structure due to magnetorotational instability (MRI) turbulence, our results show that observations of the soft X-ray excess can be used to study disc turbulence in the radiation pressure dominated regime. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes. △ Less

Submitted 3 April, 2018; v1 submitted 30 March, 2018; originally announced April 2018.

Comments: 35 pages, 50 figures, accepted to MNRAS

Convection Enhances Magnetic Turbulence in AM CVn Accretion Disks

Authors: Matthew S. B. Coleman, Omer Blaes, Shigenobu Hirose, Peter H. Hauschildt

Abstract: We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state ta… ▽ More We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $α$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $α$ near the tip of the upper branch caused by convection; this increase in $α$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight. △ Less

Submitted 12 March, 2018; originally announced March 2018.

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

Convective Quenching of Field Reversals in Accretion Disc Dynamos

Authors: Matthew S. B. Coleman, Evan Yerger, Omer Blaes, Greg Salvesen, Shigenobu Hirose

Abstract: Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antis… ▽ More Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antisymmetric about the disc midplane. Using data from the simulations of Hirose et al. (2014), we demonstrate that the lack of field reversals in the presence of convection is due to hydrodynamic mixing of magnetic field from the more strongly magnetized upper layers into the midplane, which then annihilate field reversals that are starting there. Our convective simulations differ in several respects from those reported in previous work by others, in which stronger magnetization likely plays a more important role than convection. △ Less

Submitted 27 January, 2017; originally announced January 2017.

Comments: 12 pages, 10 figures, 1 table, accepted for publication in MNRAS

The contribution of bulk Comptonization to the soft X-ray excess in AGN

Authors: Jason Kaufman, Omer Blaes, Shigenobu Hirose

Abstract: Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynami… ▽ More Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynamic (MHD) shearing box simulations. We characterize our results with temperatures and optical depths to make contact with other warm Comptonization models of the soft excess. We chose our fiducial mass, $M = 2 \times 10^6 M_{\odot}$, and accretion rate, $L/L_{\rm Edd} = 2.5$, to correspond to those fit to the super-Eddington narrow line Seyfert 1 (NLS1) RE1034+396. The temperatures, optical depths, and Compton $y$ parameters we find broadly agree with those fit to RE1034+396. The effect of bulk Comptonization is to shift the Wien tail to higher energy and lower the gas temperature, broadening the spectrum. Observations of the soft excess in NLS1s can constrain the properties of disc turbulence if the bulk Comptonization contribution can be separated out from contributions from other physical effects, such as reflection and absorption. △ Less

Submitted 23 January, 2017; originally announced January 2017.

The Effects of Magnetic Fields on the Dynamics of Radiation Pressure Dominated Massive Star Envelopes

Authors: Yan-Fei Jiang, Matteo Cantiello, Lars Bildsten, Eliot Quataert, Omer Blaes

Abstract: We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy… ▽ More We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfvén velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars. △ Less

Submitted 19 December, 2016; originally announced December 2016.

Comments: 12 pages, 10 figures, 2 tables, movies available at https://goo.gl/3kYbtg

Dwarf Nova Outbursts with Magnetorotational Turbulence

Authors: M. S. B. Coleman, I. Kotko, O. Blaes, J. -P. Lasota, S. Hirose

Abstract: The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $α$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dw… ▽ More The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $α$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $α$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $α$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $α$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior. △ Less

Submitted 3 August, 2016; originally announced August 2016.

Comments: 19 pages, 14 figures, 2 tables, accepted for publication in MNRAS

Bulk Comptonization by Turbulence in Accretion Disks

Authors: J. Kaufman, O. M. Blaes

Abstract: Radiation pressure dominated accretion discs around compact objects may have turbulent velocities that greatly exceed the electron thermal velocities within the disc. Bulk Comptonization by the turbulence may therefore dominate over thermal Comptonization in determining the emergent spectrum. Bulk Comptonization by divergenceless turbulence is due to radiation viscous dissipation only. It can be t… ▽ More Radiation pressure dominated accretion discs around compact objects may have turbulent velocities that greatly exceed the electron thermal velocities within the disc. Bulk Comptonization by the turbulence may therefore dominate over thermal Comptonization in determining the emergent spectrum. Bulk Comptonization by divergenceless turbulence is due to radiation viscous dissipation only. It can be treated as thermal Comptonization by solving the Kompaneets equation with an equivalent "wave" temperature, which is a weighted sum over the power present at each scale in the turbulent cascade. Bulk Comptonization by turbulence with non-zero divergence is due to both pressure work and radiation viscous dissipation. Pressure work has negligible effect on photon spectra in the limit of optically thin turbulence, and in this limit radiation viscous dissipation alone can be treated as thermal Comptonization with a temperature equivalent to the full turbulent power. In the limit of extremely optically thick turbulence, radiation viscous dissipation is suppressed, and the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids. We discuss the consequences of these effects for self-consistently resolving and interpreting turbulent Comptonization in spectral calculations in radiation MHD simulations of high luminosity accretion flows. △ Less

Submitted 4 April, 2016; v1 submitted 31 March, 2016; originally announced April 2016.

Comments: Accepted for publication in MNRAS, corrected typos

High-Frequency and Type-C QPOs from Oscillating, Precessing Hot, Thick Flow

Authors: P. Chris Fragile, Odele Straub, Omer Blaes

Abstract: Motivated by recent studies showing an apparent correlation between the high-frequency quasi-periodic oscillations (QPOs) and the low-frequency, type-C QPO in low-mass, black hole X-ray binaries (LMXBs), we explore a model that explains all three QPOs in terms of an oscillating, precessing hot flow in the truncated-disk geometry. Our model favors attributing the two high-frequency QPOs, often occu… ▽ More Motivated by recent studies showing an apparent correlation between the high-frequency quasi-periodic oscillations (QPOs) and the low-frequency, type-C QPO in low-mass, black hole X-ray binaries (LMXBs), we explore a model that explains all three QPOs in terms of an oscillating, precessing hot flow in the truncated-disk geometry. Our model favors attributing the two high-frequency QPOs, often occurring in a near 3:2 frequency ratio, to the breathing and vertical epicyclic frequency modes of the hot, thick flow, although we can not rule out the Keplerian and m=-1 radial epicyclic modes. In either case, the type-C QPO is attributed to precession. The correlation of the QPOs comes from the fact that all three frequencies are associated with the same geometrical structure. While the exact QPO frequencies are sensitive to the black hole mass and spin, their evolution over the course of an outburst is mainly tied to the truncation radius between the geometrically thin, optically thick disk and the inner, hot flow. We show that, in the case of the LMXB GRO J1655-40, this model can explain the one simultaneous observation of all three QPOs and that an extrapolation of the model appears to match lower frequency observations where only two of the three components are seen. Thus, this model may be able to unify multiple QPO observations using the properties of a single, simple, geometrical model. △ Less

Submitted 17 June, 2016; v1 submitted 25 February, 2016; originally announced February 2016.

Comments: 7 pages, 4 figures, accepted to MNRAS

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

Authors: Yan-Fei Jiang, Matteo Cantiello, Lars Bildsten, Eliot Quataert, Omer Blaes

Abstract: We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on th… ▽ More We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $τ_0$, to $τ_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $τ_0 \gg τ_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $τ_0 \ll τ_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius. △ Less

Submitted 17 September, 2015; originally announced September 2015.

Comments: 20 pages, 21 figures. submitted to ApJ

Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst

Authors: Shigenobu Hirose, Omer Blaes, Julian H. Krolik, Matthew S. B. Coleman, Takayoshi Sano

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

Submitted 11 April, 2014; v1 submitted 12 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ

Journal ref: 2014 ApJ 787 1

A model of the steep power law spectra and high-frequency quasi-periodic oscillations in luminous black hole X-ray binaries

Authors: Jason Dexter, Omer Blaes

Abstract: We propose a new model of the steep power law state of luminous black hole X-ray binaries. The model uses the fact that at high luminosities, the inner radii of radiation pressure dominated accretion discs are expected to i) become effectively optically thin and ii) produce significant luminosities. The gas temperature therefore rises sharply inwards, producing local saturated Compton spectra with… ▽ More We propose a new model of the steep power law state of luminous black hole X-ray binaries. The model uses the fact that at high luminosities, the inner radii of radiation pressure dominated accretion discs are expected to i) become effectively optically thin and ii) produce significant luminosities. The gas temperature therefore rises sharply inwards, producing local saturated Compton spectra with rapidly increasing peak energy. These spectra sum together to form a steep power law tail to the spectrum. A given photon energy on this tail corresponds to a narrow range in radius, so that local vertical oscillations of the disc naturally produce high quality high-frequency quasi-periodic oscillations (HFQPOs) in the hard X-ray band. The two lowest order modes have a robust frequency ratio of ~1.53. This model explains the appearance of steep power law spectra and HFQPOs at high luminosity, the 3:2 HFQPO frequency ratios, and their association with the power law spectral component. We predict an increase in QPO quality factor when the power spectrum is restricted to a narrower photon energy band, and an increase in HFQPO frequency at higher X-ray energies or lower luminosities. Future X-ray telescopes could detect additional HFQPOs from higher order modes. We demonstrate how this model could be used to measure black hole spin from HFQPOs, and qualitatively estimate the spin of GRO J1655-40 as a/M ~ 0.4 - 0.7. △ Less

Submitted 3 December, 2013; originally announced December 2013.

Comments: 5 pages, 4 figures, re-submitted to MNRAS Letters after minor revisions

Physical Properties of the Inner Shocks in Hot, Tilted Black Hole Accretion Flows

Authors: Aleksey Generozov, Omer Blaes, P. Chris Fragile, Ken B. Henisey

Abstract: Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories… ▽ More Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories in the simulations through the shocks, and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction ($\simeq3-12$ percent) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted. △ Less

Submitted 21 November, 2013; originally announced November 2013.

Comments: Accepted to The Astrophysical Journal

Effects of local dissipation profiles on magnetized accretion disk spectra

Authors: Ted Tao, Omer Blaes

Abstract: We present spectral calculations of non-LTE accretion disk models appropriate for high luminosity stellar mass black hole X-ray binary systems. We first use a dissipation profile based on scaling the results of shearing box simulations of Hirose et al. (2009) to a range of annuli parameters. We simultaneously scale the effective temperature, orbital frequency and surface density with luminosity an… ▽ More We present spectral calculations of non-LTE accretion disk models appropriate for high luminosity stellar mass black hole X-ray binary systems. We first use a dissipation profile based on scaling the results of shearing box simulations of Hirose et al. (2009) to a range of annuli parameters. We simultaneously scale the effective temperature, orbital frequency and surface density with luminosity and radius according to the standard α-model (Shakura & Sunyaev, 1973). This naturally brings increased dissipation to the disk surface layers (around the photospheres) at small radii and high luminosities. We find that the local spectrum transitions directly from a modified black body to a saturated Compton scattering spectrum as we increase the effective temperature and orbital frequency while decreasing midplane surface density. Next, we construct annuli models based on the parameters of a L/L_Edd=0.8 disk orbiting a 6.62 solar mass black hole using two modified dissipation profiles that explicitly put more dissipation per unit mass near the disk surface. The new dissipation profiles are qualitatively similar to the one found by Hirose et al. (2009), but produce strong near power-law spectral tails. Our models also include physically motivated magnetic acceleration support based once again on scaling the Hirose et al. (2009) results. We present three full-disk spectra each based on one of the dissipation prescriptions. Our most aggressive dissipation profile results in a disk spectrum that is in approximate quantitative agreement with certain observations of the steep power law (SPL) spectral states from some black hole X-ray binaries. △ Less

Submitted 14 May, 2013; originally announced May 2013.

Comments: Accepted by ApJ

Journal ref: ApJ 2013, 770, 55

General Overview of Black Hole Accretion Theory

Authors: Omer Blaes

Abstract: I provide a broad overview of the basic theoretical paradigms of black hole accretion flows. Models that make contact with observations continue to be mostly based on the four decade old alpha stress prescription of Shakura & Sunyaev (1973), and I discuss the properties of both radiatively efficient and inefficient models, including their local properties, their expected stability to secular pertu… ▽ More I provide a broad overview of the basic theoretical paradigms of black hole accretion flows. Models that make contact with observations continue to be mostly based on the four decade old alpha stress prescription of Shakura & Sunyaev (1973), and I discuss the properties of both radiatively efficient and inefficient models, including their local properties, their expected stability to secular perturbations, and how they might be tied together in global flow geometries. The alpha stress is a prescription for turbulence, for which the only existing plausible candidate is that which develops from the magnetorotational instability (MRI). I therefore also review what is currently known about the local properties of such turbulence, and the physical issues that have been elucidated and that remain uncertain that are relevant for the various alpha-based black hole accretion flow models. △ Less

Submitted 17 April, 2013; originally announced April 2013.

Comments: To be published in Space Science Reviews and as hard cover in the Space Sciences Series of ISSI: The Physics of Accretion on to Black Holes (Springer Publisher)

Variability from Nonaxisymmetric Fluctuations Interacting with Standing Shocks in Tilted Black Hole Accretion Disks

Authors: Ken B. Henisey, Omer M. Blaes, P. Chris Fragile

Abstract: We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is incon… ▽ More We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is inconsistent with a modal interpretation. We find instead that previously identified, transient, over-dense clumps orbiting on roughly Keplerian trajectories appear generically in our global simulations, independent of tilt. Associated with these fluctuations are acoustic spiral waves interior to the orbits of the clumps. We show that the two nonaxisymmetric standing shock structures that exist in the inner regions of these tilted flows effectively amplify the variability caused by these spiral waves to markedly higher levels than in untilted flows, which lack standing shocks. Our identification of clumps, spirals, and spiral-shock interactions in these fully general relativistic, magnetohydrodynamical simulations suggests that these features may be important dynamical elements in models which incorporate tilt as a way to explain the observed variability in black hole accretion flows. △ Less

Submitted 14 November, 2012; v1 submitted 9 November, 2012; originally announced November 2012.

Comments: accepted in ApJ October 1, 2012

Local Radiation MHD Instabilities in Magnetically Stratified Media

Authors: Ted Tao, Omer Blaes

Abstract: We study local radiation magnetohydrodynamic instabilities in static, optically thick, vertically stratified media with constant flux mean opacity. We include the effects of vertical gradients in a horizontal background magnetic field. Assuming rapid radiative diffusion, we use the zero gas pressure limit as an entry point for investigating the coupling between the photon bubble instability and th… ▽ More We study local radiation magnetohydrodynamic instabilities in static, optically thick, vertically stratified media with constant flux mean opacity. We include the effects of vertical gradients in a horizontal background magnetic field. Assuming rapid radiative diffusion, we use the zero gas pressure limit as an entry point for investigating the coupling between the photon bubble instability and the Parker instability. Apart from factors that depend on wavenumber orientation, the Parker instability exists for wavelengths longer than a characteristic wavelength lambda_{tran}, while photon bubbles exist for wavelengths shorter than lambda_{tran}. The growth rate in the Parker regime is independent of the orientation of the horizontal component of the wavenumber when radiative diffusion is rapid, but the range of Parker-like wavenumbers is extended if there exists strong horizontal shear between field lines (i.e. horizontal wavenumber perpendicular to the magnetic field). Finite gas pressure introduces an additional short wavelength limit to the Parker-like behavior, and also limits the growth rate of the photon bubble instability to a constant value at short wavelengths. We also consider the effects of differential rotation with accretion disk applications in mind. Our results may explain why photon bubbles have not yet been observed in recent stratified shearing box accretion disk simulations. Photon bubbles may physically exist in simulations with high radiation to gas pressure ratios, but higher spatial resolution will be needed to resolve the asymptotically growing unstable wavelengths. △ Less

Submitted 3 October, 2011; originally announced October 2011.

Comments: The Astrophysical Journal, in press

Intrinsic disc emission and the Soft X-ray Excess in AGN

Authors: Chris Done, Shane Davis, Chichuan Jin, Omer Blaes, Martin Ward

Abstract: (Abridged) Narrow Line Seyfert 1 (NLS1) galaxies have low mass black holes and mass accretion rates close to (or exceeding) Eddington, so a standard blackbody accretion disc should peak in the EUV. However, the lack of true absorption opacity in the disc means that the emission is better approximated by a colour temperature corrected blackbody, and this colour temperature correction is large enoug… ▽ More (Abridged) Narrow Line Seyfert 1 (NLS1) galaxies have low mass black holes and mass accretion rates close to (or exceeding) Eddington, so a standard blackbody accretion disc should peak in the EUV. However, the lack of true absorption opacity in the disc means that the emission is better approximated by a colour temperature corrected blackbody, and this colour temperature correction is large enough ($\sim 2.4$) that the bare disc emission from a zero spin black hole can extend into the soft X-ray bandpass. Part of the soft X-ray excess seen in these objects must be intrinsic emission from the disc unless the vertical structure is very different to that predicted. However, the soft excess is much broader than predicted by a bare disc spectrum, indicating some Compton upscattering by cool, optically thick material. We associate this with the disc itself, so it must ultimately be powered by mass accretion. We build an energetically self consistent model assuming that the emission thermalises at large radii, but that at smaller radii the gravitational energy is split between powering optically thick Comptonised disc emission (forming the soft X-ray excess) and an optically thin corona above the disc (forming the tail to higher energies). We show examples of this model fit to the extreme NLS1 REJ1034+396, and to the much lower Eddington fraction Broad Line Seyfert 1 PG1048+231. We use these to guide our fits and interpretations of three template spectra made from co-adding multiple sources to track out a sequence of AGN spectra as a function of $L/L_{Edd}$. The new model is publically available within the {\sc xspec} spectral fitting package. △ Less

Submitted 6 September, 2011; v1 submitted 27 July, 2011; originally announced July 2011.

Comments: 13 pages, 9 figures, accepted in MNRAS

Dissipation and Vertical Energy Transport in Radiation-Dominated Accretion Disks

Authors: Omer Blaes, Julian H. Krolik, Shigenobu Hirose, Natalia Shabaltas

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

Submitted 25 March, 2011; originally announced March 2011.

Comments: accepted for publication in The Astrophysical Journal

Excitation of Trapped Waves in Simulations of Tilted Black Hole Accretion Disks with Magnetorotational Turbulence

Authors: Ken B. Henisey, Omer M. Blaes, P. Chris Fragile, Bárbara T. Ferreira

Abstract: We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic simulations of accretion flows onto a black hole with dimensionless spin parameter a/M=0.9. We consider both the case where the angular momentum of the accretion material is aligned with the black hole spin axis (an untilted flow) and where it is misaligned by 15 degrees (a tilted flow). In comparison… ▽ More We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic simulations of accretion flows onto a black hole with dimensionless spin parameter a/M=0.9. We consider both the case where the angular momentum of the accretion material is aligned with the black hole spin axis (an untilted flow) and where it is misaligned by 15 degrees (a tilted flow). In comparison to the untilted simulation, the tilted simulation exhibits a clear excess of inertial variability, that is, variability at frequencies below the local radial epicyclic frequency. We further study the radial structure of this inertial-like power by focusing on a radially extended band at 118 (M/10Msol)^-1 Hz found in each of the three analyzed fluid variables. The three dimensional density structure at this frequency suggests that the power is a composite oscillation whose dominant components are an over dense clump corotating with the background flow, a low order inertial wave, and a low order inertial-acoustic wave. Our results provide preliminary confirmation of earlier suggestions that disk tilt can be an important excitation mechanism for inertial waves. △ Less

Submitted 9 October, 2009; originally announced October 2009.

Comments: 8 Pages, 6 Figures, accepted for publication in ApJ

Journal ref: Astrophys.J.706:705-711,2009

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

Authors: Shigenobu Hirose, Omer Blaes, Julian H. Krolik

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

Submitted 7 August, 2009; originally announced August 2009.

Comments: accepted for publication in The Astrophysical Journal

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

AGN accretion disks as spatially resolved by polarimetry

Authors: Makoto Kishimoto, Robert Antonucci, Omer Blaes, Andy Lawrence, Catherine Boisson, Marcus Albrecht, Christian Leipski

Abstract: A crucial difficulty in understanding the nature of the putative accretion disk in AGNs is that some of its key intrinsic spectral signatures cannot be observed directly. The strong emissions from the broad-line region (BLR) and the obscuring torus, which are generally yet to be spatially resolved, essentially 'bury' such signatures. Here we argue that we can actually isolate the disk emission s… ▽ More A crucial difficulty in understanding the nature of the putative accretion disk in AGNs is that some of its key intrinsic spectral signatures cannot be observed directly. The strong emissions from the broad-line region (BLR) and the obscuring torus, which are generally yet to be spatially resolved, essentially 'bury' such signatures. Here we argue that we can actually isolate the disk emission spectrum by using optical and near-infrared polarization of quasars and uncover the important spectral signatures. In these quasars, the polarization is considered to originate from electron scattering interior to the BLR, so that the polarized flux shows the disk spectrum with all the emissions from the BLR and torus eliminated. The polarized flux observations have now revealed a Balmer edge feature in absorption and a blue near-infrared spectral shape consistent with a specific and robust theoretical prediction. These results critically verify the long-standing picture of an optically thick and locally heated disk in AGNs. △ Less

Submitted 8 October, 2008; originally announced October 2008.

Comments: Proceedings for "The Universe under the Microscope" (AHAR 2008), held in Bad Honnef (Germany) in April 2008, to be published in Journal of Physics: Conference Series by Institute of Physics Publishing, R. Schoedel, A. Eckart, S. Pfalzner, and E. Ros (eds.)

Journal ref: J.Phys.Conf.Ser.131:012039,2008

Radiation-Dominated Disks Are Thermally Stable

Authors: Shigenobu Hirose, Julian H. Krolik, Omer Blaes

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

Submitted 10 September, 2008; originally announced September 2008.

Comments: Accepted by ApJ

Journal ref: Astrophys.J.691:16-31,2009

The characteristic blue spectra of accretion disks in quasars as uncovered in the infrared

Authors: Makoto Kishimoto, Robert Antonucci, Omer Blaes, Andy Lawrence, Catherine Boisson, Marcus Albrecht, Christian Leipski

Abstract: Quasars are thought to be powered by supermassive black holes accreting surrounding gas. Central to this picture is a putative accretion disk which is believed to be the source of the majority of the radiative output. It is well known, however, that the most extensively studied disk model -- an optically thick disk which is heated locally by the dissipation of gravitational binding energy -- is… ▽ More Quasars are thought to be powered by supermassive black holes accreting surrounding gas. Central to this picture is a putative accretion disk which is believed to be the source of the majority of the radiative output. It is well known, however, that the most extensively studied disk model -- an optically thick disk which is heated locally by the dissipation of gravitational binding energy -- is apparently contradicted by observations in a few major respects. In particular, the model predicts a specific blue spectral shape asymptotically from the visible to the near-infrared, but this is not generally seen in the visible wavelength region where the disk spectrum is observable. A crucial difficulty was that, toward the infrared, the disk spectrum starts to be hidden under strong hot dust emission from much larger but hitherto unresolved scales, and thus has essentially been impossible to observe. Here we report observations of polarized light interior to the dust-emiting region that enable us to uncover this near-infrared disk spectrum in several quasars. The revealed spectra show that the near-infrared disk spectrum is indeed as blue as predicted. This indicates that, at least for the outer near-infrared-emitting radii, the standard picture of the locally heated disk is approximately correct. The model problems at shorter wavelengths should then be directed toward a better understanding of the inner parts of the revealed disk. The newly uncovered disk emission at large radii, with more future measurements, will also shed totally new light on the unanswered critical question of how and where the disk ends. △ Less

Submitted 23 July, 2008; originally announced July 2008.

Comments: published in Nature, 24 July 2008 issue. Supplementary Information can be found at http://www.mpifr-bonn.mpg.de/div/ir-interferometry/suppl_info.pdf Published version can be accessed from http://www.nature.com/nature/journal/v454/n7203/pdf/nature07114.pdf

Epicyclic Motions and Standing Shocks in Numerically Simulated Tilted Black-Hole Accretion Disks

Authors: P. Chris Fragile, Omer M. Blaes

Abstract: This work presents a detailed analysis of the overall flow structure and unique features of the inner region of the tilted disk simulations described in Fragile et al. (2007). The primary new feature identified in the main disk body is a latitude-dependent radial epicyclic motion driven by pressure gradients attributable to the gravitomagnetic warping of the disk. The induced motion of the gas i… ▽ More This work presents a detailed analysis of the overall flow structure and unique features of the inner region of the tilted disk simulations described in Fragile et al. (2007). The primary new feature identified in the main disk body is a latitude-dependent radial epicyclic motion driven by pressure gradients attributable to the gravitomagnetic warping of the disk. The induced motion of the gas is coherent over the scale of the entire disk and is fast enough that it could be observable in features such as relativistic iron lines. The eccentricity of the associated fluid element trajectories increases with decreasing radius, leading to a crowding of orbit trajectories near their apocenters. This results in a local density enhancement akin to a compression. These compressions are sufficiently strong to produce a pair of weak shocks in the vicinity of the black hole. These shocks are roughly aligned with the line-of-nodes between the black-hole symmetry plane and disk midplane, with one shock above the line-of-nodes on one side of the black hole and the other below on the opposite side. These shocks enhance angular momentum transport and energy dissipation near the hole, forcing some material to plunge toward the black hole from well outside the innermost stable circular orbit. A new, extended simulation, which was evolved for more than a full disk precession period, allows us to confirm that these shocks and the previously identified ``plunging streams'' precess with the disk in such a way as to remain aligned relative to the line-of-nodes, as expected based on our physical understanding of these phenomena. Such a precessing structure would likely present a strong quasi-periodic signal. △ Less

Submitted 15 July, 2008; originally announced July 2008.

Comments: 10 pages, 14 figures. 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

Epicyclic Oscillations of Fluid Bodies: Newtonian Nonslender Torus

Authors: Omer M. Blaes, Eva Sramkova, Marek A. Abramowicz, Wlodek Kluzniak, Ulf Torkelsson

Abstract: We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which… ▽ More We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which is important for the theory of high-frequency quasi-periodic oscillations (QPOs), occurs when the fluid tori epicyclic frequencies nu_r and nu_z are about 15% higher than the ones corresponding to free particles. Our results therefore suggest that previous estimates of black hole spins from QPOs have produced values that are too high. △ Less

Submitted 29 June, 2007; originally announced June 2007.

Comments: ApJ, in press

Journal ref: Astrophys.J.665:642-653,2007

Global General Relativistic MHD Simulation of a Tilted Black-Hole Accretion Disk

Authors: P. Chris Fragile, Omer M. Blaes, Peter Anninos, Jay D. Salmonson

Abstract: This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the… ▽ More This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range. △ Less

Submitted 28 June, 2007; originally announced June 2007.

Comments: Accepted to ApJ. 32 pages, 16 fig. Please see http://www.cofc.edu/~fragilep/tilted.html for movies and high-resolution figures

Surface Structure in an Accretion Disk Annulus with Comparable Radiation and Gas Pressure

Authors: Omer Blaes, Shigenobu Hirose, Julian H. Krolik

Abstract: We have employed a 3-d energy-conserving radiation MHD code to simulate the vertical structure and thermodynamics of a shearing box whose parameters were chosen so that the radiation and gas pressures would be comparable. The upper layers of this disk segment are magnetically-dominated, creating conditions appropriate for both photon bubble and Parker instabilities. We find little evidence for p… ▽ More We have employed a 3-d energy-conserving radiation MHD code to simulate the vertical structure and thermodynamics of a shearing box whose parameters were chosen so that the radiation and gas pressures would be comparable. The upper layers of this disk segment are magnetically-dominated, creating conditions appropriate for both photon bubble and Parker instabilities. We find little evidence for photon bubbles, even though the simulation has enough spatial resolution to see them and their predicted growth rates are high. On the other hand, there is strong evidence for Parker instabilities, and they appear to dominate the evolution of the magnetically supported surface layers. The disk photosphere is complex, with large density inhomogeneities at both the scattering and effective (thermalization) photospheres of the evolving horizontally-averaged structure. Both the dominant magnetic support and the inhomogeneities are likely to have strong effects on the spectrum and polarization of thermal photons emerging from the disk atmosphere. The inhomogeneities are also large enough to affect models of reflection spectra from the atmospheres of accretion disks. △ Less

Submitted 2 May, 2007; originally announced May 2007.

Comments: ApJ, in press

Journal ref: Astrophys.J.664:1057-1071,2007

Thermodynamics of an Accretion Disk Annulus with Comparable Radiation and Gas Pressure

Authors: Julian H. Krolik, Shigenobu Hirose, Omer Blaes

Abstract: We explore the thermodynamic and global structural properties of a local patch of an accretion disk whose parameters were chosen so that radiation pressure and gas pressure would be comparable in magnitude. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and… ▽ More We explore the thermodynamic and global structural properties of a local patch of an accretion disk whose parameters were chosen so that radiation pressure and gas pressure would be comparable in magnitude. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure and energy balance of this disk segment over a span of more than forty cooling times. As is also true when gas pressure dominates, the disk's upper atmosphere is magnetically-supported. However, unlike the gas-dominated case, no steady-state is reached; instead, the total (i.e., radiation plus gas) energy content fluctuates by factors of 3--4 over timescales of several tens of orbits, with no secular trend. Because the radiation pressure varies much more than the gas pressure, the ratio of radiation pressure to gas pressure varies over the approximate range 0.5--2. The volume-integrated dissipation rate generally increases with increasing total energy, but the mean trend is somewhat slower than linear, and the instantaneous dissipation rate is often a factor of two larger or smaller than the mean for that total energy level. Locally, the dissipation rate per unit volume scales approximately in proportion to the current density; the time-average dissipation rate per unit mass is proportional to m^{-1/2}, where m is the horizontally-averaged mass column density to the nearer of the top or bottom surface. As in our earlier study of a gas-dominated shearing-box, we find that energy transport is completely dominated by radiative diffusion, with Poynting flux carrying less than 1% of the energy lost from the box. △ Less

Submitted 2 May, 2007; originally announced May 2007.

Comments: ApJ, in press

Journal ref: Astrophys.J.664:1045-1056,2007

Accretion Disks in AGNs

Authors: Omer Blaes

Abstract: I briefly review the theoretical models of radiatively efficient, geometrically thin and optically thick accretion disk spectra that currently exist for AGN. I then discuss three recent observational developments that have real potential to teach us about the physics of these flows. Finally, I present results on the most recent, thermodynamically consistent simulations of magnetorotational turbu… ▽ More I briefly review the theoretical models of radiatively efficient, geometrically thin and optically thick accretion disk spectra that currently exist for AGN. I then discuss three recent observational developments that have real potential to teach us about the physics of these flows. Finally, I present results on the most recent, thermodynamically consistent simulations of magnetorotational turbulence and discuss what these simulations are suggesting about the vertical structure of accretion disks. △ Less

Submitted 22 March, 2007; originally announced March 2007.

Comments: Review Talk, to appear in "The Central Engine of Active Galactic Nuclei", ed. L. C. Ho and J.-M. Wang, (San Francisco: ASP)

3D Relativistic MHD Simulation of a Tilted Accretion Disk Around a Rapidly Rotating Black Hole

Authors: P. Chris Fragile, Peter Anninos, Omer M. Blaes, Jay D. Salmonson

Abstract: We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk… ▽ More We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole. △ Less

Submitted 9 January, 2007; originally announced January 2007.

Comments: 3 pages, 3 figures, to appear in "Proceedings of the 11th Marcel Grossmann Meeting on General Relativity", eds H. Kleinert, R. T. Jantzen, R. Ruffini

The near-IR shape of the big blue bump emission from quasars: under the hot dust emission

Authors: Makoto Kishimoto, Robert Antonucci, Omer Blaes

Abstract: One primary difficulty in understanding the nature of the putative accretion disk in the central engine of AGNs is that some of its key intrinsic spectral signatures are buried under the emissions from the surrounding regions, i.e. the broad line region (BLR) and the obscuring torus. We argue here that these signatures can be revealed by using optical and near-IR polarization. At least in some q… ▽ More One primary difficulty in understanding the nature of the putative accretion disk in the central engine of AGNs is that some of its key intrinsic spectral signatures are buried under the emissions from the surrounding regions, i.e. the broad line region (BLR) and the obscuring torus. We argue here that these signatures can be revealed by using optical and near-IR polarization. At least in some quasars, the polarization is seen only in the continuum and is not shared by emission lines. In this case, the polarized flux is considered to show the intrinsic spectrum interior to the BLR, removing off the emissions from the BLR and torus. We have used this polarization to reveal the Balmer-edge feature and near-IR spectral shape of the central engine, both of which are important for testing the fundamental aspects of the models. △ Less

Submitted 20 December, 2006; originally announced December 2006.

Comments: 4 pages, 2 figures. Conference proceedings to appear in "The Central Engine of Active Galactic Nuclei", ed. L. C. Ho and J.-M. Wang (San Francisco: ASP)

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

Quasi-Periodic Oscillations from Magnetorotational Turbulence

Authors: Phil Arras, Omer Blaes, Neal J. Turner

Abstract: Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole binaries have been widely interpreted as being due to standing wave modes in accretion disks. These disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric so… ▽ More Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole binaries have been widely interpreted as being due to standing wave modes in accretion disks. These disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric sound waves and radial epicyclic motions driven by MRI turbulence give rise to narrow, distinct peaks in the temporal power spectrum. Inertial waves, on the other hand, do not give rise to distinct peaks which rise significantly above the continuum noise spectrum set by MRI turbulence, even when the fluid motions are projected onto the eigenfunctions of the modes. This is a serious problem for QPO models based on inertial waves. △ Less

Submitted 13 February, 2006; originally announced February 2006.

Comments: 4 pages, 2 figures. submitted to apj

Journal ref: Astrophys.J. 645 (2006) L65-L68

Testing Accretion Disk Theory in Black Hole X-ray Binaries

Authors: Shane W. Davis, Chris Done, Omer M. Blaes

Abstract: We present results from spectral modeling of three black hole X-ray binaries: LMC X-3, GRO J1655-40, and XTE J1550-564. Using a sample of disk dominated observations, we fit the data with a range of spectral models that includes a simple multitemperature blackbody (DISKBB), a relativistic accretion disk model based on color-corrected blackbodies (KERRBB), and a relativistic model based on non-LT… ▽ More We present results from spectral modeling of three black hole X-ray binaries: LMC X-3, GRO J1655-40, and XTE J1550-564. Using a sample of disk dominated observations, we fit the data with a range of spectral models that includes a simple multitemperature blackbody (DISKBB), a relativistic accretion disk model based on color-corrected blackbodies (KERRBB), and a relativistic model based on non-LTE atmosphere models within an alpha prescription (BHSPEC). BHSPEC provides the best fit for a BeppoSAX observation of LMC X-3, which has the broadest energy coverage of our sample. It also provides the best fit for multiple epochs of Rossi X-ray Timing Explorer (RXTE) data in this source, except at the very highest luminosity (L/L_{Edd} >~ 0.7), where additional physics must be coming into play. BHSPEC is also the best-fit model for multi-epoch RXTE observations of GRO J1655-40 and XTE J1550-564, although the best-fit inclination of the inner disk differs from the binary inclination. All our fits prefer alpha=0.01 to alpha=0.1, in apparent disagreement with the large stresses inferred from the rapid rise times observed in outbursts of these two sources. In all three sources our fits imply moderate black hole spins (a/M ~ 0.1 - 0.8), but this is sensitive to the reliability of independent measurements of these system parameters and to the physical assumptions which underly our spectral models. △ Less

Submitted 7 May, 2006; v1 submitted 10 February, 2006; originally announced February 2006.

Comments: 15 pages, emulateapj, ApJ, in press

Journal ref: Astrophys.J.647:525-538,2006