BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Omer Blaes

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

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.SRastro-ph.CO

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convection Causes Enhanced Magnetic Turbulence in Accretion Disks in Outburst

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of local dissipation profiles on magnetized accretion disk spectra

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General Overview of Black Hole Accretion Theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dissipation and Vertical Energy Transport in Radiation-Dominated Accretion Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

(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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Local Radiation MHD Instabilities in Magnetically Stratified Media

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

AGN accretion disks as spatially resolved by polarimetry

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Radiation-Dominated Disks Are Thermally Stable

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

0 citations0 reviewsNo code linkedOpen access