BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Aleksander Sadowski

Aleksander Sadowski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.GAastro-ph.COastro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

28 published item(s)

preprint2016arXiv

Levitating atmospheres of Eddington-luminosity neutron stars

We construct models of static, spherically symmetric shells supported by the radiation flux of a luminous neutron star in the Schwarzschild metric. The atmospheres are disconnected from the star and levitate above its surface. Gas pressure and density inversion appear in the inner region of these atmospheres, which is a purely relativistic phenomenon. We account for the scattering opacity dependence on temperature and utilize the relativistic M1 closure scheme for the radiation tensor, hence allowing for a fully GR-consistent treatment of the photon flux and radiation tensor anisotropy. In this way we are able to address atmospheres of both large and moderate/low optical depths with the same set of equations. We discuss properties of the levitating atmospheres and find that they may indeed be optically thick, with the distance between star surface and the photosphere expanding as luminosity increases. These results may be relevant for the photosphereric radius expansion X-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Variability and mm/IR flares in GRMHD Models of Sgr A* from Strong-Field Gravitational Lensing

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HERO: A 3D General Relativistic Radiative Postprocessor for Accretion Discs around Black Holes

HERO (Hybrid Evaluator for Radiative Objects) is a 3D general relativistic radiative transfer code which has been tailored to the problem of analyzing radiation from simulations of relativistic accretion discs around black holes. HERO is designed to be used as a postprocessor. Given some fixed fluid structure for the disc (i.e. density and velocity as a function of position from a hydrodynamics or magnetohydrodynamics simulation), the code obtains a self-consistent solution for the radiation field and for the gas temperatures using the condition of radiative equilibrium. The novel aspect of HERO is that it combines two techniques: 1) a short characteristics (SC) solver that quickly converges to a self consistent disc temperature and radiation field, with 2) a long characteristics (LC) solver that provides a more accurate solution for the radiation near the photosphere and in the optically thin regions. By combining these two techniques, we gain both the computational speed of SC and the high accuracy of LC. We present tests of HERO on a range of 1D, 2D and 3D problems in flat space and show that the results agree well with both analytical and benchmark solutions. We also test the ability of the code to handle relativistic problems in curved space. Finally, we discuss the important topic of ray-defects, a major limitation of the SC method, and describe our strategy for minimizing the induced error.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HEROIC: 3D General Relativistic Radiative Postprocessor with Comptonization for Black Hole Accretion Discs

We describe HEROIC, an upgraded version of the relativistic radiative post-processor code HERO described in a previous paper, but which now Includes Comptonization. HEROIC models Comptonization via the Kompaneets equation, using a quadratic approximation for the source function in the short characteristics radiation solver. It employs a simple form of accelerated lambda iteration to handle regions of high scattering opacity. In addition to solving for the radiation field, HEROIC also solves for the gas temperature by applying the condition of radiative equilibrium. We present benchmarks and tests of the Comptonization module in HEROIC with simple 1D and 3D scattering problems. We also test the ability of the code to handle various relativistic effects using model atmospheres and accretion flows in a black hole space-time. We present two applications of HEROIC to general relativistic MHD simulations of accretion discs. One application is to a thin accretion disc around a black hole. We find that the gas below the photosphere in the multi-dimensional HEROIC solution is nearly isothermal, quite different from previous solutions based on 1D plane parallel atmospheres. The second application is to a geometrically thick radiation-dominated accretion disc accreting at 11 times the Eddington rate. The multi-dimensional HEROIC solution shows that, for observers who are on axis and look down the polar funnel, the isotropic equivalent luminosity could be more than ten times the Eddington limit, even though the spectrum might still look thermal and show no signs of relativistic beaming.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High frequency oscillations in outbursts of Kerr-metric slim disks

We numerically investigate the thermally unstable accretion disks around black holes. We adopt an evolutionary viscous stress equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find a kind of interesting oscillations on some running models in limit-cycle outburst state. The oscillations arise near the inner boundary and propagate radially outwards. We deem that they are the trapped $p$-mode oscillations excited by sonic-point instability. We directly integrate the local radiation cooling fluxes to construct the mimic bolometric light-curve. We find a series of overtones beside the fundamental harmonic on the power spectra of mimic light-curves. The frequency of the fundamental harmonic is very close to the maximum epicyclic frequency of the disk and the frequency ratio of the fundamental harmonic and overtones is a regular integer series. We suggest that the code for ray-tracing calculation must be time-dependent in virtual observation and point out the robustness of the black hole spin measurement with high frequency QPOs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Jet and disk luminosities in tidal disruption events

Tidal disruption events (TDE) in which a star is devoured by a massive black hole at a galac- tic center pose a challenge to our understanding of accretion processes. Within a month the accretion rate reaches super-Eddington levels. It then drops gradually over a time scale of a year to sub-Eddington regimes. The initially geometrically thick disk becomes a thin one and eventually an ADAF at very low accretion rates. As such, TDEs explore the whole range of accretion rates and configurations. A challenging question is what the corresponding light curves of these events are. We explore numerically the disk luminosity and the conditions within the inner region of the disk using a fully general relativistic slim disk model. Those conditions determine the magnitude of the magnetic field that engulfs the black hole and this, in turn, determines the Blandford-Znajek jet power. We estimate this power in two different ways and show that they are self-consistent. We find, as expected earlier from analytic argu- ments (Krolik & Piran 2012), that neither the disk luminosity nor the jet power follows the accretion rate throughout the disruption event. The disk luminosity varies only logarithmi- cally with the accretion rate at super-Eddington luminosities. The jet power follows initially the accretion rate but remains a constant after the transition from super- to sub- Eddington. At lower accretion rates at the end of the MAD phase the disk becomes thin and the jet may stop altogether. These new estimates of the jet power and disk luminosity that do not simply follow the mass fallback rate should be taken into account when searching for TDEs and analysing light curves of TDE candidates. Identification of some of the above mentioned transitions may enable us to estimate better TDE parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nature of the Extreme Ultraluminous X-ray Sources

In this proof-of-concept study we demonstrate that in a binary system mass can be transferred toward an accreting compact object at extremely high rate. If the transferred mass is efficiently converted to X-ray luminosity (with disregard of the classical Eddington limit) or if the X-rays are focused into a narrow beam then binaries can form extreme ULX sources with the X-ray luminosity of Lx>10^42 erg/s. For example, Lasota & King argued that the brightest known ULX (HLX-1) is a regular binary system with a rather low-mass compact object (a stellar-origin black hole or a neutron star). The predicted formation efficiencies and lifetimes of binaries with the very high mass transfer rates are large enough to explain all observed systems with extreme X-ray luminosities. These systems are not only limited to binaries with stellar-origin black hole accretors. Noteworthy, we have also identified such objects with neutron stars. Typically, a 10 Msun black hole is fed by a massive (10 Msun) Hertzsprung gap donor with Roche lobe overflow rate of 10^-3 Msun/yr (2600 MEdd). For neutron star systems the typical donors are evolved low-mass (2 Msun) helium stars with Roche lobe overflow rate of 10^-2 Msun/yr. Our study does not prove that any particular extreme ULX is a regular binary system, but it demonstrates that any ULX, including the most luminous ones, may potentially be a short-lived phase in the life of a binary star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Simulation of Hot Accretion Flows (III): Revisiting wind properties using trajectory approach

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $ε_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Photon-conserving Comptonization in simulations of accretion disks around black holes

We introduce a new method for treating Comptonization in computational fluid dynamics. By construction, this method conserves the number of photons. Whereas the traditional &#34;blackbody Comptonization&#34; approach assumes that the radiation is locally a perfect blackbody and therefore uses a single parameter, the radiation temperature, to describe the radiation, the new &#34;photon-conserving Comptonization&#34; approach treats the photon gas as a Bose-Einstein fluid and keeps track of both the radiation temperature and the photon number density. We have implemented photon-conserving Comptonization in the general relativistic radiation magnetohydrodynamical code KORAL and we describe its impact on simulations of mildly super-critical black hole accretion disks. We find that blackbody Comptonization underestimates the gas and radiation temperature by up to a factor of two compared to photon-conserving Comptonization. This discrepancy could be serious when computing spectra. The photon-conserving simulation indicates that the spectral color correction factor of the escaping radiation in the funnel region of the disk could be as large as 5.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Powerful radiative jets in super-critical accretion disks around non-spinning black holes

We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates $\gtrsim 10 \dot M_{\rm Edd}$, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity $\approx 0.3c$. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as $10^{48}\,\rm erg\,s^{-1}$ for a $10^7 M_\odot$ BH accreting at $10^3$ Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor $γ> 2$ inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly super-Eddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-dimensional simulations of super-critical black hole accretion disks --- luminosities, photon trapping and variability

We present a set of four three-dimensional, general relativistic, radiation MHD simulations of black hole accretion at super-critical mass accretion rates, $\dot{M} > \dot{M}_{\rm Edd}$. We use these simulations to study how disk properties are modified when we vary the black hole mass, the black hole spin, or the mass accretion rate. In the case of a non-rotating black hole, we find that the total efficiency is of order $3\%\dot M c^2$, approximately a factor of two less than the efficiency of a standard thin accretion disk. The radiation flux in the funnel along the axis is highly super-Eddington, but only a small fraction of the energy released by accretion escapes in this region. The bulk of the $3\%\dot M c^2$ of energy emerges farther out in the disk, either in the form of photospheric emission or as a wind. In the case of a black hole with a spin parameter of 0.7, we find a larger efficiency of about $8\%\dot M c^2$. By comparing the relative importance of advective and diffusive radiation transport, we show that photon trapping is effective near the equatorial plane. However, near the disk surface, vertical transport of radiation by diffusion dominates. We compare the properties of our fiducial three-dimensional run with those of an equivalent two-dimensional axisymmetric model with a mean-field dynamo. The latter simulation runs nearly 100 times faster than the three-dimensional simulation, and gives very similar results for time-averaged properties of the accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $η\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the effect of the Sgr A* accretion flow on the appearance of G2 after pericenter

We study the dynamical interaction of the G2 cloud with the accretion flow around Sgr A* by means of three-dimensional, hydrodynamic simulations. We show the effects of the rotating accretion flow on the evolution of G2 by projecting the cloud density onto the plane of the sky, and extracting position-velocity diagrams. We study a number of possible orientations of the cloud orbit with respect to the disk. We find that once the center of mass of the cloud has crossed the pericenter, the differences between models becomes significant. Models with the cloud counter-rotating with respect to the disk are expected to reach higher blue-shifted line of sight velocities. The spatial extent of the emission depends strongly on the cloud-to-disk inclination angle. Future imaging and spectroscopy of G2 emission will shed light both on the structure of the Sgr A* disk and on the properties of the cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Power of Imaging: Constraining the Plasma Properties of GRMHD Simulations using EHT Observations of Sgr A*

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black-holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast GPU-based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all twelve models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black-holes and their thermodynamic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General Relativistic Magnetohydrodynamic Simulations of Blandford-Znajek Jets and the Membrane Paradigm

Recently it has been observed that the scaling of jet power with black hole spin in galactic X-ray binaries is consistent with the predictions of the Blandford-Znajek (BZ) jet model. These observations motivate us to revisit the BZ model using general relativistic magnetohydrodynamic simulations of magnetized jets from accreting (h/r ~ 0.3), spinning (0 < a_* < 0.98) black holes. We have three main results. First, we quantify the discrepancies between the BZ jet power and our simulations: assuming maximum efficiency and uniform fields on the horizon leads to a ~10% overestimate of jet power, while ignoring the accretion disk leads to a further ~50% overestimate. Simply reducing the standard BZ jet power prediction by 60% gives a good fit to our simulation data. Our second result is to show that the membrane formulation of the BZ model correctly describes the physics underlying simulated jets: torques, dissipation, and electromagnetic fields on the horizon. This provides intuitive yet rigorous pictures for the black hole energy extraction process. Third, we compute the effective resistance of the load region and show that the load and the black hole achieve near perfect impedance matching. Taken together, these results increase our confidence in the BZ model as the correct description of jets observed from astrophysical black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes

A numerical scheme is described for including radiation in multi-dimensional general-relativistic conservative fluid dynamics codes. In this method, a covariant form of the M1 closure scheme is used to close the radiation moments, and the radiative source terms are treated semi-implicitly in order to handle both optically thin and optically thick regimes. The scheme has been implemented in a conservative general relativistic radiation hydrodynamics code KORAL. The robustness of the code is demonstrated on a number of test problems, including radiative relativistic shock tubes, static radiation pressure supported atmosphere, shadows, beams of light in curved spacetime, and radiative Bondi accretion. The advantages of M1 closure relative to other approaches such as Eddington closure and flux-limited diffusion are discussed, and its limitations are also highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure

Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20\%$, the large-scale electromagnetic jet efficiency is of order $10\%$, and the total radiative efficiency that reaches large distances remains low at only order $1\%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs.\ disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

GRMHD Simulations of Magnetized Advection Dominated Accretion on a Non-Spinning Black Hole: Role of Outflows

We present results from two long-duration GRMHD simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200,000GM/c^3 and achieved inflow equilibrium out to a radius \sim90GM/c^2. Even at this relatively large radius, the mass outflow rate \dot{M}_{out} is found to be only 60% of the net mass inflow rate \dot{M}_{BH} into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100,000GM/c^3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius \sim170GM/c^2. Here, \dot{M}_{out} becomes equal to \dot{M}_{BH} at r\sim 160GM/c^2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Shakura-Sunyaev Viscosity Prescription with Variable alpha(r)

Almost all hydrodynamic accretion disk models parametrize viscosity with the dimensionless parameter alpha. There is no detailed model for alpha, so it is usually taken to be a constant. However, global simulations of magnetohydrodynamic disks find that alpha varies with distance from the central object. Also, Newtonian simulations tend to find smaller alpha&#39;s than general relativistic simulations. We seek a one-dimensional model for alpha that can reproduce these two observations. We are guided by data from six general relativistic magnetohydrodynamic accretion disk simulations. The variation of alpha in the inner, laminar regions of the flow results from stretching of mean magnetic field lines by the flow. The variation of alpha in the outer, turbulent regions results from the dependence of the magnetorotational instability on the dimensionless shear rate. We give a one-dimensional prescription for alpha(r) that captures these two effects and reproduces the radial variation of alpha observed in the simulations. For thin disks, the prescription simplifies to the formula alpha(r)=0.025[q(r)/1.5]^6, where the shear parameter, q(r), is an analytical function of radius in the Kerr metric. The coefficient and exponent are inferred from our simulations and will change as better simulation data becomes available. We conclude that the alpha-viscosity prescription can be extended to the radially varying alpha&#39;s observed in simulations. It is possible that Newtonian simulations find smaller alpha&#39;s than general relativistic simulations because the shear parameter is lower in Newtonian flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraints on the black hole spin in the quasar SDSS J094533.99+100950.1

The spin of the black hole is an important parameter which may be responsible for the properties of the inflow and outflow of the material surrounding a black hole. Broad band IR/optical/UV spectrum of the quasar SDSS J094533.99+100950.1 is clearly disk-dominated, with the spectrum peaking up in the observed frequency range. Therefore, disk fitting method usually used for Galactic black holes can be used in this object to determine the black hole spin. We develop the numerical code for computing disk properties, including radius-dependent hardening factor, and we apply the ray-tracing method to incorporate all general relativity effects in light propagation. We show that the simple multicolor disk model gives a good fit, without any other component required, and the disk extends down to the marginally stable orbit. The best fit accretion rate is 0.13, well below the Eddington limit, and the black hole spin is moderate, 0.3. The contour error for the fit combined with the constraints for the black hole mass and the disk inclination gives a constraint that the spin is lower than 0.8. We discuss the sources of possible systematic errors in the parameter determinations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Relativistic slim disks with vertical structure

We report on a scheme for incorporating vertical radiative energy transport into a fully relativistic, Kerr-metric model of optically thick, advective, transonic alpha disks. Our code couples the radial and vertical equations of the accretion disk. The flux was computed in the diffusion approximation, and convection is included in the mixing-length approximation. We present the detailed structure of this &#34;two-dimensional&#34; slim-disk model for alpha=0.01. We then calculated the emergent spectra integrated over the disk surface. The values of surface density, radial velocity, and the photospheric height for these models differ by 20%-30% from those obtained in the polytropic, height-averaged slim disk model considered previously. However, the emission profiles and the resulting spectra are quite similar for both types of models. The effective optical depth of the slim disk becomes lower than unity for high values of the alpha parameter and for high accretion rates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spinning up black holes with super-critical accretion flows

We study the process of spinning up black holes by accretion from slim disks in a wide range of accretion rates. We show that for super-Eddington accretion rates and low values of the viscosity parameter alpha (\lesssim 0.01) the limiting value of the dimensionless spin parameter a_* can reach values higher than a_*=0.9978 inferred by Thorne (1974) in his seminal study. For \mdot=10\mdot_{Edd} and alpha=0.01 spin equilibrium is reached at a_*=0.9994. We show that the equilibrium spin value depends strongly on the assumed value of alpha. We also prove that for high accretion rates the impact of captured radiation on spin evolution is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stability of radiation-pressure dominated disks. I. The dispersion relation for a delayed heating alpha-viscosity prescription

We derive and investigate the dispersion relation for accretion disks with retarded or advanced heating. We follow the alpha-prescription but allow for a time offset (τ) between heating and pressure perturbations, as well as for a diminished response of heating to pressure variations. We study in detail solutions of the dispersion relation for disks with radiation-pressure fraction 1 - β. For τ<0 (delayed heating) the number and sign of real solutions for the growth rate depend on the values of the time lag and the ratio of heating response to pressure perturbations, ξ. If the delay is larger than a critical value (e.g., if Ωτ<-125 for α=0.1, β=0 and ξ=1) two real solutions exist, which are both negative. These results imply that retarded heating may stabilize radiation-pressure dominated accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing slim-disk models on the thermal spectra of LMC X-3

Slim-disk models describe accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3. LMC X-3 has the virtues of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity (thin-disk) and high luminosity (slim-disk) regimes. We also show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thin Disk Theory with a Non-Zero Torque Boundary Condition and Comparisons with Simulations

We present an analytical solution for thin disk accretion onto a Kerr black hole that extends the standard Novikov-Thorne alpha-disk in three ways: (i) it incorporates nonzero stresses at the inner edge of the disk, (ii) it extends into the plunging region, and (iii) it uses a corrected vertical gravity formula. The free parameters of the model are unchanged. Nonzero boundary stresses are included by replacing the Novikov-Thorne no torque boundary condition with the less strict requirement that the fluid velocity at the innermost stable circular orbit is the sound speed, which numerical models show to be the correct behavior for luminosities below ~30% Eddington. We assume the disk is thin so we can ignore advection. Boundary stresses scale as alpha*h and advection terms scale as h^2 (where h is the disk opening angle (h=H/r)), so the model is self-consistent when h < alpha. We compare our solution with slim disk models and general relativistic magnetohydrodynamic disk simulations. The model may improve the accuracy of black hole spin measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Launching of jets by cold, magnetized disks in Kerr Metric

We confirm recent discovery by Cao that in the vicinity of fast rotating black holes jets can be launched centrifugally by cold, magnetized disks even for nearly vertically shaped magnetic flux surfaces. Outflows produced under such extreme conditions are investigated via studying kinematics of test particles in the force-free magnetosphere approximation. Implications of a possibility of magneto-centrifugal launching of very well collimated central outflows around the fast rotating black holes are discussed in the general context of the jet formation scenarios in AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Leaving the ISCO: the inner edge of a black-hole accretion disk at various luminosities

The &#34;radiation inner edge&#34; of an accretion disk is defined as the inner boundary of the region from which most of the luminosity emerges. Similarly, the &#34;reflection edge&#34; is the smallest radius capable of producing a significant X-ray reflection of the fluorescent iron line. For black hole accretion disks with very sub-Eddington luminosities these and all other &#34;inner edges&#34; locate at ISCO. Thus, in this case, one may rightly consider ISCO as the unique inner edge of the black hole accretion disk. However, even for moderate luminosities, there is no such unique inner edge as differently defined edges locate at different places. Several of them are significantly closer to the black hole than ISCO. The differences grow with the increasing luminosity. For nearly Eddington luminosities, they are so huge that the notion of the inner edge losses all practical significance.

0 citations0 reviewsNo code linkedOpen access