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

Ramesh Narayan

Ramesh Narayan 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-phgr-qcastro-ph.GAastro-ph.COastro-ph.IMhep-thastro-ph.SRphysics.plasm-phhep-phnucl-exnucl-thphysics.flu-dynphysics.space-ph

Catalog footprint

What is connected

107works
14topics
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

107 published item(s)

preprint2022arXiv

Analytical Model of Disk Evaporation and State Transitions in Accreting Black Holes

State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process which is suitable for analytical and numerical studies. With radius $r$ scaled to Schwarzschild units and coronal mass accretion rate $\dot{m}_c$ to Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and AGN. The corona solution consists of two power-law segments separated at a break radius $r_b \sim10^3 \,(α/0.3)^{-2}$, where $α$ is the viscosity parameter. Gas evaporates from the disk to the corona for $r>r_b$, and condenses back for $r<r_b$. At $r_b$, $\dot{m}_c$ reaches its maximum, $\dot{m}_{c,{\rm max}} \approx 0.02\, (α/0.3)^3$. If at $r\gg r_b$ the thin disk accretes with $\dot{m}_d < \dot{m}_{c,{\rm max}} $, then the disk evaporates fully before reaching $r_b$, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model which includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, and $r_b$ remains unchanged. This model predicts strong coronal winds for $r>r_b$, and a $T\sim 5\times 10^8\,{\rm K}$ Compton-cooled corona for $r < r_b$. Two-temperature effects are ignored, but may be important at small radii.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density jump as a function of magnetic field for switch-on collisionless shocks in pair plasmas

The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. In \cite{BretJPP2018}, we devised a model for the kinetic history of the plasma through the shock front, allowing to self-consistently compute the downstream anisotropy, hence the density jump, in terms of the upstream parameters. This model dealt with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only 1 switch-on solution within a limited range of Alfvén Mach numbers, our model offers 2 kinds of solutions within a slightly different range of Alfvén Mach numbers. These 2 solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density jump for oblique collisionless shocks in pair plasmas: allowed solutions

Shockwaves in plasma are usually dealt with using Magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing to account for kinetic effects within an MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalog, and a density jump that is notably less in the quasi parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Detectability of Wandering Intermediate-Mass Black Holes in the Milky Way Galaxy from Radio to X-rays

Intermediate-mass black holes (IMBHs, $10^{3-6} \, \rm{M_\odot}$), are typically found at the center of dwarf galaxies and might be wandering, thus far undetected, in the Milky Way (MW). We use model spectra for advection-dominated accretion flows to compute the typical fluxes, in a range of frequencies spanning from radio to X-rays, emitted by a putative population of $10^5 \, \rm{M_\odot}$ IMBHs wandering in five realistic, volume-weighted, MW environments. We predict that $\sim 27\%$ of the wandering IMBHs can be detected in the X-ray with Chandra, $\sim 37\%$ in the near-infrared with the Roman Space Telescope, $\sim 49\%$ in the sub-mm with CMB-S4 and $\sim 57\%$ in the radio with ngVLA. We find that the brightest fluxes are emitted by IMBHs passing through molecular clouds or cold neutral medium, where they are always detectable. We propose criteria to facilitate the selection of candidates in multi-wavelength surveys. Specifically, we compute the X-ray to optical ratio ($α_{\rm ox}$) and the optical to sub-mm ratio, as a function of the accretion rate of the IMBH. We show that at low rates the sub-mm emission of IMBHs is significantly higher than the optical, UV and X-ray emission. Finally, we place upper limits on the number $N_\bullet$ of these objects in the MW: $N_\bullet<2000$ and $N_\bullet<100$, based on our detectability expectations and current lack of detections in molecular clouds and cold neutral medium, respectively. These predictions will guide future searches of IMBHs in the MW, which will be instrumental to understanding their demographics and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Jets in Magnetically Arrested Hot Accretion Flows: Geometry, Power and Black Hole Spindown

We present the results of nine simulations of radiatively-inefficient magnetically arrested disks (MADs) across different values of the black hole spin parameter $a_*$: $-0.9$, $-0.7$, $-0.5$, $-0.3$, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100,000\,GM/c^3$ to ensure disk inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD disks depends strongly on the black hole spin, confirming previous results. Prograde disks saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to $100\;GM/c^2$, the power-law index is $k \approx 0.27-0.42$. There is a strong correlation between the disk-jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disk variability: accretion rate variability increases with increasing spin for $a_*>0$ and remains almost constant for $a_*\lesssim 0$, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observational Constraints on Direct Electron Heating in the Hot Accretion Flows in Sgr A* and M87*

An important parameter in the theory of hot accretion flows around black holes is $δ$, which describes the fraction of ``viscously'' dissipated energy in the accretion flow that goes directly into heating electrons. For a given mass accretion rate, the radiative efficiency of a hot accretion flow is determined by $δ$. Unfortunately, the value of $δ$ is hard to determine from first principles. The recent Event Horizon Telescope Collaboration (EHTC) results on M87* and Sgr A* provide us with a different way of constraining $δ$. By combining the mass accretion rates in M87* and Sgr A* estimated by the EHTC with the measured bolometric luminosities of the two sources, we derive good constraints on the radiative efficiencies of the respective accretion flows. In parallel, we use a theoretical model of hot magnetically arrested disks (MAD) to calculate the expected radiative efficiency as a function of $δ$ (and accretion rate). By comparing the EHTC-derived radiative efficiencies with the theoretical results from MAD models, we find that Sgr A* requires $δ\ga 0.3$. %with the most likely value being $δ\sim 0.5$. A similar comparison in the case of M87* gives inconclusive results as there is still a large uncertainty in the accretion rate in this source.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observational properties of puffy disks: radiative GRMHD spectra of mildly sub-Eddington accretion

Numerical general relativistic radiative magnetohydrodynamic simulations of accretion disks around a stellar mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy disks. Above a dense and geometrically thin core of dimensionless thickness $h/r \sim 0.1$, crudely resembling a classic thin accretion disk, a puffed-up, geometrically thick layer of lower density and $h/r \sim 1.0$ is formed. We discuss the observational properties of puffy disks, in particular the geometrical obscuration of the inner disk by the elevated puffy region at higher observing inclinations, and collimation of the radiation along the accretion disk spin axis, which may explain the apparent super-Eddington luminosity of some X-ray objects. We also present synthetic spectra of puffy disks, and show that they are qualitatively similar to those of a Comptonized thin disk. We demonstrate that the existing xspec spectral fitting models provide good fits to synthetic observations of puffy disks, but cannot correctly recover the input black hole spin. The puffy region remains optically thick to scattering; in its spectral properties the puffy disk roughly resembles that of a warm corona sandwiching the disk core. We suggest that puffy disks may correspond to X-ray binary systems of luminosities above 0.3 of the Eddington luminosity in the intermediate spectral states.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Probing Plasma Physics with Spectral Index Maps of Accreting Black Holes on Event Horizon Scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. In this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximised near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disk (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Toward determining the number of observable supermassive black hole shadows

We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ${\sim}0.1$ $μ$as resolution and ${\sim}1$ $μ$Jy sensitivity could access ${>}10^6$ SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ${\lesssim}2$ $μ$as resolution and ${\lesssim}0.5$ mJy sensitivity, while doing the same for second-order photon rings requires ${\lesssim}0.1$ $μ$as resolution and ${\lesssim}5$ $μ$Jy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), while a next-generation EHT observing at 345 GHz should have access to ${\sim}$3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Brightness Asymmetry of Black Hole Images as a Probe of Observer Inclination

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Faraday depolarization and induced circular polarization by multi-path propagation with application to FRBs

We describe how the observed polarization properties of an astronomical object are related to its intrinsic polarization properties and the finite temporal and spectral resolutions of the observing device. Moreover, we discuss the effect that a scattering screen, with non-zero magnetic field, between the source and observer has on the observed polarization properties. We show that the polarization properties are determined by the ratio of observing bandwidth and coherence bandwidth of the scattering screen and the ratio of temporal resolution of the instrument and the variability time of screen, as long as the length over which the Faraday rotation induced by the screen changes by $\simπ$ is smaller than the size of the screen visible to the observer. We describe the conditions under which a source that is 100\% linearly polarized intrinsically might be observed as partially depolarized, and how the source's temporal variability can be distinguished from the temporal variability induced by the scattering screen. In general, linearly polarized waves passing through a magnetized scattering screen can develop a significant circular polarization. We apply the work to the observed polarization properties of a few fast radio bursts (FRBs), and outline potential applications to pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accretion disks around naked singularities

We investigate here the thermal properties of accretion disks in a spacetime for some galactic density profiles in spherical symmetry. The matter distributions have a finite outer radius with a naked central singularity. The luminosities of the accretion disks for some density profile models are found to be higher than those for a Schwarzschild black hole of the same mass. The slopes for the luminosity distributions with respect to frequencies are significantly different, especially at higher frequencies, from that in the Schwarzschild black hole case. Such features may be used to distinguish black holes from naked singularities. The efficiencies for the conversion of the mass energy of the accreting gas into radiation and the strength of naked singularities are analyzed. The novel feature that we find is, the strength of the singularity is different depending on the profiles considered, and the stronger the singularity is, the higher is the efficiency for the accretion disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reconnection-driven particle acceleration in relativistic shear flows

Particle energization in shear flows is invoked to explain non-thermal emission from the boundaries of relativistic astrophysical jets. Yet, the physics of particle injection, i.e., the mechanism that allows thermal particles to participate in shear-driven acceleration, remains unknown. With particle-in-cell simulations, we study the development of Kelvin-Helmholtz (KH) instabilities seeded by the velocity shear between a relativistic magnetically-dominated electron-positron jet and a weakly magnetized electron-ion ambient plasma. We show that, in their nonlinear stages, KH vortices generate kinetic-scale reconnection layers, which efficiently energize the jet particles, thus providing a first-principles mechanism for particle injection into shear-driven acceleration. Our work lends support to spine-sheath models of jet emission - with a fast core/spine surrounded by a slower sheath - and can explain the origin of radio-emitting electrons at the boundaries of relativistic jets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal Interferometric Signatures of a Black Hole's Photon Ring

The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Black Hole Paradoxes

We propose here that the well-known black hole paradoxes such as the information loss and teleological nature of the event horizon are restricted to a particular idealized case, which is the homogeneous dust collapse model. In this case, the event horizon, which defines the boundary of the black hole, forms initially, and the singularity in the interior of the black hole at a later time. We show that, in contrast, gravitational collapse from physically more realistic initial conditions typically leads to the scenario in which the event horizon and space-time singularity form simultaneously. We point out that this apparently simple modification can mitigate the causality and teleological paradoxes, and also lends support to two recently suggested solutions to the information paradox, namely, the `firewall' and `classical chaos' proposals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy flows in thick accretion disks and their consequences for black hole feedback

We study energy flows in geometrically thick accretion discs, both optically thick and thin, using general relativistic, three-dimensional simulations of black hole accretion flows. We find that for non-rotating black holes the efficiency of the total feedback from thick accretion discs is $3\%$ - roughly half of the thin disc efficiency. This amount of energy is ultimately distributed between outflow and radiation, the latter scaling weakly with the accretion rate for super-critical accretion rates, and returned to the interstellar medium. Accretion on to rotating black holes is more efficient because of the additional extraction of rotational energy. However, the jet component is collimated and likely to interact only weakly with the environment, whereas the outflow and radiation components cover a wide solid angle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extreme Brightness Temperatures and Refractive Substructure in 3C273 with RadioAstron

Earth-space interferometry with RadioAstron provides the highest direct angular resolution ever achieved in astronomy at any wavelength. RadioAstron detections of the classic quasar 3C273 on interferometric baselines up to 171,000 km suggest brightness temperatures exceeding expected limits from the "inverse-Compton catastrophe" by two orders of magnitude. We show that at 18 cm, these estimates most probably arise from refractive substructure introduced by scattering in the interstellar medium. We use the scattering properties to estimate an intrinsic brightness temperature of 7*10^12 K, which is consistent with expected theoretical limits, but which is ~15 times lower than estimates that neglect substructure. At 6.2 cm, the substructure influences the measured values appreciably but gives an estimated brightness temperature that is comparable to models that do not account for the substructure. At 1.35 cm, the substructure does not affect the extremely high inferred brightness temperatures, in excess of 10^13 K. We also demonstrate that for a source having a Gaussian surface brightness profile, a single long-baseline estimate of refractive substructure determines an absolute minimum brightness temperature, if the scattering properties along a given line of sight are known, and that this minimum accurately approximates the apparent brightness temperature over a wide range of total flux densities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High Resolution Linear Polarimetric Imaging for the Event Horizon Telescope

Images of the linear polarization of synchrotron radiation around Active Galactic Nuclei (AGN) identify their projected magnetic field lines and provide key data for understanding the physics of accretion and outflow from supermassive black holes. The highest resolution polarimetric images of AGN are produced with Very Long Baseline Interferometry (VLBI). Because VLBI incompletely samples the Fourier transform of the source image, any image reconstruction that fills in unmeasured spatial frequencies will not be unique and reconstruction algorithms are required. In this paper, we explore extensions of the Maximum Entropy Method (MEM) to linear polarimetric VLBI imaging. In contrast to previous work, our polarimetric MEM algorithm combines a Stokes I imager that uses only bispectrum measurements that are immune to atmospheric phase corruption with a joint Stokes Q and U imager that operates on robust polarimetric ratios. We demonstrate the effectiveness of our technique on 7- and 3-mm wavelength quasar observations from the VLBA and simulated 1.3-mm Event Horizon Telescope observations of Sgr A* and M87. Consistent with past studies, we find that polarimetric MEM can produce superior resolution compared to the standard CLEAN algorithm when imaging smooth and compact source distributions. As an imaging framework, MEM is highly adaptable, allowing a range of constraints on polarization structure. Polarimetric MEM is thus an attractive choice for image reconstruction with the EHT.

0 citations0 reviewsNo code linkedOpen access
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
preprint2016arXiv

The Optics of Refractive Substructure

Newly recognized effects of refractive scattering in the ionized interstellar medium have broad implications for very long baseline interferometry (VLBI) at extreme angular resolutions. Building upon work by Blandford & Narayan (1985), we present a simplified, geometrical optics framework, which enables rapid, semi-analytic estimates of refractive scattering effects. We show that these estimates exactly reproduce previous results based on a more rigorous statistical formulation. We then derive new expressions for the scattering-induced fluctuations of VLBI observables such as closure phase, and we demonstrate how to calculate the fluctuations for arbitrary quantities of interest using a Monte Carlo technique.

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

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 "blackbody Comptonization" 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 "photon-conserving Comptonization" 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

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Event Horizon of M87

The 6 billion solar mass supermassive black hole at the center of the giant elliptical galaxy M87 powers a relativistic jet. Observations at millimeter wavelengths with the Event Horizon Telescope have localized the emission from the base of this jet to angular scales comparable to the putative black hole horizon. The jet might be powered directly by an accretion disk or by electromagnetic extraction of the rotational energy of the black hole. However, even the latter mechanism requires a confining thick accretion disk to maintain the required magnetic flux near the black hole. Therefore, regardless of the jet mechanism, the observed jet power in M87 implies a certain minimum mass accretion rate. If the central compact object in M87 were not a black hole but had a surface, this accretion would result in considerable thermal near-infrared and optical emission from the surface. Current flux limits on the nucleus of M87 strongly constrain any such surface emission. This rules out the presence of a surface and thereby provides indirect evidence for an event horizon.

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
preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A High-Frequency Doppler Feature in the Power Spectra of Simulated GRMHD Black Hole Accretion Disks

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collisionless Weibel shocks: full formation mechanism and timing

Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Confirmation Via the Continuum-Fitting Method that the Spin of the Black Hole in Cygnus X-1 is Extreme

In Gou et al. (2011), we reported that the black hole primary in the X-ray binary Cygnus X-1 is a near-extreme Kerr black hole with a spin parameter a*>0.95(3σ). We confirm this result while setting a new and more stringent limit: a*>0.983 at the 3σ(99.7%) level of confidence. The earlier work, which was based on an analysis of all three useful spectra that were then available, was possibly biased by the presence in these spectra of a relatively strong Compton power-law component: The fraction of the thermal seed photons scattered into the power law was f_s=23-31%, while the upper limit for reliable application of the continuum-fitting method is f_s<25%. We have subsequently obtained six additional spectra of Cygnus X-1 suitable for the measurement of spin. Five of these spectra are of high quality with f_s in the range 10% to 19%, a regime where the continuum-fitting method has been shown to deliver reliable results. Individually, the six spectra give lower limits on the spin parameter that range from a*>0.95 to a*>0.98, allowing us to conservatively conclude that the spin of the black hole is a*>0.983 (3σ).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distinguishing black holes from naked singularities through their accretion disk properties

We show that, in principle, a slowly evolving gravitationally collapsing perfect fluid cloud can asymptotically settle to a static spherically symmetric equilibrium configuration with a naked singularity at the center. We consider one such asymptotic final configuration with a finite outer radius, and construct a toy model in which it is matched to a Schwarzschild exterior geometry. We examine the properties of circular orbits in this model. We then investigate observational signatures of a thermal accretion disk in this spacetime, comparing them with the signatures expected for a disk around a black hole of the same mass. Several notable differences emerge. A disk around the naked singularity is much more luminous than one around an equivalent black hole. Also, the disk around the naked singularity has a spectrum with a high frequency power law segment that carries a major fraction of the total luminosity. Thus, at least some naked singularities can, in principle, be distinguished observationally from black holes of the same mass. We discuss possible implications of these results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. I. Compression-Driven Instabilities and the Electron Heating Mechanism

In systems accreting well below the Eddington rate, the plasma in the innermost regions of the disk is collisionless and two-temperature, with the ions hotter than the electrons. Yet, whether a collisionless faster-than-Coulomb energy transfer mechanism exists in two-temperature accretion flows is still an open question. We study the physics of electron heating during the growth of ion velocity-space instabilities, by means of multi-dimensional particle-in-cell (PIC) simulations. A large-scale compression - embedded in a novel form of the PIC equations - continuously amplifies the field. This constantly drives a pressure anisotropy P_perp > P_parallel, due to the adiabatic invariance of the particle magnetic moments. We find that, for ion plasma beta values beta_i ~ 5-30 appropriate for the midplane of low-luminosity accretion flows, mirror modes dominate if the electron-to-proton temperature ratio is > 0.2, whereas if it is < 0.2 the ion cyclotron instability triggers the growth of strong Alfven-like waves, that pitch-angle scatter the ions to maintain marginal stability. We develop an analytical model of electron heating during the growth of the ion cyclotron instability, which we validate with PIC simulations. We find that for cold electrons (beta_e < m_e/m_i), the electron energy gain is controlled by the magnitude of the E-cross-B velocity induced by the ion cyclotron waves. This term is independent of the initial electron temperature, so it provides a solid energy floor even for electrons starting with extremely low temperatures. On the other hand, the electron energy gain for beta_e > m_e/m_i - governed by the conservation of the magnetic moment in the growing fields of the instability - is proportional to the initial electron temperature. Our results have implications for two-temperature accretion flows as well as the solar wind and intracluster plasmas. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Event-Horizon-Telescope Evidence for Alignment of the Black Hole in the Center of the Milky Way with the Inner Stellar Disk

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

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

Globular Clusters and Dark Satellite Galaxies through the Stream Velocity

The formation of purely baryonic globular clusters with no gravitationally bound dark matter is still a theoretical challenge. We show that these objects might form naturally whenever there is a relative stream velocity between baryons and dark matter. The stream velocity causes a phase shift between linear modes of baryonic and dark matter perturbations, which translates to a spatial offset between the two components when they collapse. For a 2sigma (3sigma) density fluctuation, baryonic clumps with masses in the range 1e5 - 2.5e6 Msun (1e5 - 4e6 Msun) collapse outside the virial radii of their counterpart dark matter halos. These objects could survive as long-lived dark matter-free objects and might conceivably become globular clusters. In addition, their dark matter counterparts, which were deprived of gas, might become dark satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hot Accretion Flows Around Black Holes

Black hole accretion flows can be divided into two broad classes: cold and hot. Cold accretion flows, which consist of cool optically thick gas, are found at relatively high mass accretion rates. Prominent examples are the standard thin disk, which occurs at a fraction of the Eddington mass accretion rate, and the slim disk at super-Eddington rates. These accretion flows are responsible for luminous systems such as active galactic nuclei radiating at or close to the Eddington luminosity and black hole X-ray binaries in the soft state. Hot accretion flows, the topic of this review, are virially hot and optically thin. They occur at lower mass accretion rates, and are described by models such as the advection-dominated accretion flow and luminous hot accretion flow. Because of energy advection, the radiative efficiency of these flows is in general lower than that of a standard thin accretion disk. Moreover, the efficiency decreases with decreasing mass accretion rate. Observations show that hot accretion flows are associated with jets. In addition, theoretical arguments suggest that hot flows should produce strong winds. Hot accretion flows are believed to be present in low-luminosity active galactic nuclei and in black hole X-ray binaries in the hard and quiescent states. The prototype is Sgr A*, the ultra-low-luminosity supermassive black hole at our Galactic center. The jet, wind and radiation from a supermassive black hole with a hot accretion flow can interact with the external interstellar medium and modify the evolution of the host galaxy. Details of this "maintenance-mode feedback" could, in principle, be worked out through theoretical studies and numerical simulations of hot accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Imaging an Event Horizon: Mitigation of Scattering Toward Sagittarius A*

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Thermal Electron Acceleration in Low Mach Number Collisionless Shocks. I. Particle Energy Spectra and Acceleration Mechanism

Electron acceleration to non-thermal energies in low Mach number (M<5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M=3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15 percent of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p~2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration. The accelerated electrons are then reflected back upstream, where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock, for further energization via shock drift acceleration. In summary, the self-generated waves allow for repeated cycles of shock drift acceleration, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Thermal Electron Acceleration in Low Mach Number Collisionless Shocks. II. Firehose-Mediated Fermi Acceleration and its Dependence on Pre-Shock Conditions

Electron acceleration to non-thermal energies is known to occur in low Mach number (M<5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M=3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves, that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (>20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Observational Evidence for Black Holes

Astronomers have discovered two populations of black holes: (i) stellar-mass black holes with masses in the range 5 to 30 solar masses, millions of which are present in each galaxy in the universe, and (ii) supermassive black holes with masses in the range 10^6 to 10^{10} solar masses, one each in the nucleus of every galaxy. There is strong circumstantial evidence that all these objects are true black holes with event horizons. The measured masses of supermassive black hole are strongly correlated with properties of their host galaxies, suggesting that these black holes, although extremely small in size, have a strong influence on the formation and evolution of entire galaxies. Spin parameters have recently been measured for a handful of black holes. Based on the data, there is an indication that the kinetic power of at least one class of relativistic jet ejected from accreting black holes may be correlated with black hole spin. If verified, it would suggest that these jets are powered by a generalized Penrose process mediated by magnetic fields.

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
preprint2014arXiv

Tidal Disruption and Magnetic Flux Capture: Powering a Jet from a Quiescent Black Hole

The transient Swift J1644+57 is believed to have been produced by an unlucky star wandering too close to a supermassive black hole (BH) leading to a tidal disruption event. This unusual flare displayed highly super-Eddington X-ray emission which likely originated in a relativistic, collimated jet. This presents challenges to modern accretion and jet theory as upper limits of prior BH activity, which we obtain from the radio afterglow of this event, imply that both the pre-disruption BH and stellar magnetic fluxes fall many orders of magnitude short of what is required to power the observed X-ray luminosity. We argue that a pre-existing, "fossil" accretion disc can contain a sufficient reservoir of magnetic flux and that the stellar debris stream is capable of dragging this flux into the BH. To demonstrate this, we perform local, 3D magnetohydrodynamic simulations of the disc--stream interaction and demonstrate that the interface between the two is unstable to mixing. This mixing entrains a sufficient amount of fossil disc magnetic flux into the infalling stellar debris to power the jet. We argue that the interaction with the fossil disc can have a pronounced effect on the structure and dynamics of mass fallback and likely the resulting transient. Finally, we describe possible ramifications of these interactions on unresolved problems in tidal disruption dynamics, in particular, the efficiency of debris circularization, and effects of the disruption on the preexisting black hole system. Animations online: http://goo.gl/T84tLs

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new equilibrium torus solution and GRMHD initial conditions

General relativistic magnetohydrodynamic (GRMHD) simulations are providing influential models for black hole spin measurements, gamma ray bursts, and supermassive black hole feedback. Many of these simulations use the same initial condition: a rotating torus of fluid in hydrostatic equilibrium. A persistent concern is that simulation results sometimes depend on arbitrary features of the initial torus. For example, the Bernoulli parameter (which is related to outflows), appears to be controlled by the Bernoulli parameter of the initial torus. In this paper, we give a new equilibrium torus solution and describe two applications for the future. First, it can be used as a more physical initial condition for GRMHD simulations than earlier torus solutions. Second, it can be used in conjunction with earlier torus solutions to isolate the simulation results that depend on initial conditions. We assume axisymmetry, an ideal gas equation of state, constant entropy, and ignore self-gravity. We fix an angular momentum distribution and solve the relativistic Euler equations in the Kerr metric. The Bernoulli parameter, rotation rate, and geometrical thickness of the torus can be adjusted independently. Our torus tends to be more bound and have a larger radial extent than earlier torus solutions. While this paper was in preparation, several GRMHD simulations appeared based on our equilibrium torus. We believe it will continue to provide a more realistic starting point for future simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black Hole Spin via Continuum Fitting and the Role of Spin in Powering Transient Jets

The spins of ten stellar black holes have been measured using the continuum-fitting method. These black holes are located in two distinct classes of X-ray binary systems, one that is persistently X-ray bright and another that is transient. Both the persistent and transient black holes remain for long periods in a state where their spectra are dominated by a thermal accretion disk component. The spin of a black hole of known mass and distance can be measured by fitting this thermal continuum spectrum to the thin-disk model of Novikov and Thorne; the key fit parameter is the radius of the inner edge of the black hole's accretion disk. Strong observational and theoretical evidence links the inner-disk radius to the radius of the innermost stable circular orbit, which is trivially related to the dimensionless spin parameter a_* of the black hole (|a_*| < 1). The ten spins that have so far been measured by this continuum-fitting method range widely from a_* \approx 0 to a_* > 0.95. The robustness of the method is demonstrated by the dozens or hundreds of independent and consistent measurements of spin that have been obtained for several black holes, and through careful consideration of many sources of systematic error. Among the results discussed is a dichotomy between the transient and persistent black holes; the latter have higher spins and larger masses. Also discussed is recently discovered evidence in the transient sources for a correlation between the power of ballistic jets and black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy Extraction from Spinning Black Holes via Relativistic Jets

It has for long been an article of faith among astrophysicists that black hole spin energy is responsible for powering the relativistic jets seen in accreting black holes. Two recent advances have strengthened the case. First, numerical general relativistic magnetohydrodynamic simulations of accreting spinning black holes show that relativistic jets form spontaneously. In at least some cases, there is unambiguous evidence that much of the jet energy comes from the black hole, not the disk. Second, spin parameters of a number of accreting stellar-mass black holes have been measured. For ballistic jets from these systems, it is found that the radio luminosity of the jet correlates with the spin of the black hole. This suggests a causal relationship between black hole spin and jet power, presumably due to a generalized Penrose process.

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

Generation of Primordial Magnetic Fields on Linear Over-density Scales

Magnetic fields appear to be present in all galaxies and galaxy clusters. Recent measurements indicate that a weak magnetic field may be present even in the smooth low density intergalactic medium. One explanation for these observations is that a seed magnetic field was generated by some unknown mechanism early in the life of the Universe, and was later amplified by various dynamos in nonlinear objects like galaxies and clusters. We show that a primordial magnetic field is expected to be generated in the early Universe on purely linear scales through vorticity induced by scale-dependent temperature fluctuations or equivalently, a spatially varying speed of sound of the gas. Residual free electrons left over after recombination tap into this vorticity to generate magnetic field via the Biermann battery process. Although the battery operates even in the absence of any relative velocity between dark matter and gas at the time of recombination, the presence of such a relative velocity modifies the predicted spatial power spectrum of the magnetic field. At redshifts of order a few tens, we estimate a root mean square field strength of order 1e-25 -1e-24 G on comoving scales~10 kpc. This field, which is generated purely from linear perturbations, is expected to be amplified significantly after reionization, and to be further boosted by dynamo processes during nonlinear structure formation.

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

Thermal Stability in Turbulently Mixed Accretion Discs

The standard thin accretion disc model predicts that discs around stellar mass black holes become radiation pressure dominated and thermally unstable once their luminosity exceeds L>0.02 L_Edd. Observationally, discs in the high/soft state of X-ray binaries show little variability in the range 0.01 L_Edd < L < 0.5 L_Edd, implying that these discs in nature are in fact quite stable. In an attempt to reconcile this conflict, we investigate one-zone disc models including turbulent and convective modes of vertical energy transport. We find both mixing mechanisms to have a stabilizing effect, leading to an increase in the L threshold up to which the disc is thermally stable. In the case of stellar mass black hole systems, convection alone leads to only a minor increase in this threshold, up to ~5 per cent of Eddington. However turbulent mixing has a much greater effect -- the threshold rises up to ~20 per cent Eddington under reasonable assumptions. In optimistic models with superefficient turbulent mixing, we even find solutions that are completely thermally stable for all accretion rates. Similar results are obtained for supermassive black holes, except that all critical accretion rates are a factor ~10 lower in Eddington ratio.

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

Double-peaked Narrow-Line Signatures of Dual Supermassive Black Holes in Galaxy Merger Simulations

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN should often be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

General Relativistic Modeling of Magnetized Jets from Accreting Black Holes

Recent advances in general relativistic magnetohydrodynamic modeling of jets offer unprecedented insights into the inner workings of accreting black holes that power the jets in active galactic nuclei (AGN) and other accretion systems. I will present the results of recent studies that determine spin-dependence of jet power and discuss the implications for the AGN radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Global Structure of Optically Thin, Magnetically Supported, Two-Temperature, Black Hole Accretion Disks

We present global solutions of optically thin, two-temperature black hole accretion disks incorporating magnetic fields. We assume that the πϕ-component of the Maxwell stress is proportional to the total pressure, and prescribe the radial dependence of the magnetic flux advection rate in order to complete the set of basic equations. We obtained magnetically supported (low-β) disk solutions, whose luminosity exceeds the maximum luminosity for an advection-dominated accretion flow (ADAF), L > 0.4 α^2 L_Edd, where L_Edd is the Eddington luminosity. The accretion flow is composed of the outer ADAF, a luminous hot accretion flow (LHAF) inside the transition layer from the outer ADAF to the low-β disk, the low-β disk, and the inner ADAF. The low-β disk region becomes wider as the mass-accretion rate increases further. In the low-β disk, the magnetic heating balances the radiative cooling, and the electron temperature decreases from ~ 10^9.5 K to ~ 10^8 K as the luminosity increases. These results are consistent with the anti-correlation between the energy cutoff in X-ray spectra (hence the electron temperature) and the luminosity when L > 0.1 L_Edd, observed in the bright/hard state during the bright hard-to-soft transitions of transient outbursts in galactic black hole candidates.

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

Jet Power and Black Hole Spin: Testing an Empirical Relationship and Using it to Predict the Spins of Six Black Holes

Using 5 GHz radio luminosity at light-curve maximum as a proxy for jet power and black-hole spin measurements obtained via the continuum-fitting method, Narayan & McClintock (2012) presented the first direct evidence for a relationship between jet power and black hole spin for four transient black-hole binaries. We test and confirm their empirical relationship using a fifth source, H1743-322, whose spin was recently measured. We show that this relationship is consistent with Fe-line spin measurements provided that the black hole spin axis is assumed to be aligned with the binary angular momentum axis. We also show that, during a major outburst of a black hole transient, the system reasonably approximates an X-ray standard candle. We further show, using the standard synchrotron bubble model, that the radio luminosity at light-curve maximum is a good proxy for jet kinetic energy. Thus, the observed tight correlation between radio power and black hole spin indicates a strong underlying link between mechanical jet power and black hole spin. Using the fitted correlation between radio power and spin for the above five calibration sources, we predict the spins of six other black holes in X-ray/radio transient systems with low-mass companions. Remarkably, these predicted spins are all relatively low, especially when compared to the high measured spins of black holes in persistent, wind-fed systems with massive companions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Mass Distribution and Birth Masses of Neutron Stars

We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radio Synchrotron Emission from a Bow Shock around the Gas Cloud G2 Heading toward the Galactic Center

A dense ionized cloud of gas has been recently discovered to be moving directly toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 Schwarzschild radii from the black hole and will move supersonically through the ambient hot gas with a velocity of v_p ~ 5400 km/s. A bow shock is likely to form in front of the cloud and could accelerate electrons to relativistic energies. We estimate via particle-in-cell simulations the energy distribution of the accelerated electrons and show that the non-thermal synchrotron emission from these electrons might exceed the quiescent radio emission from Sgr A* by a factor of several. The enhanced radio emission should be detectable at GHz and higher frequencies around the time of pericentric passage and in the following months. The bow shock emission is expected to be displaced from the quiescent radio emission of Sgr A* by ~33 mas. Interferometric observations could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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'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's observed in simulations. It is possible that Newtonian simulations find smaller alpha's than general relativistic simulations because the shear parameter is lower in Newtonian flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Angular momentum transport in convectively unstable shear flows

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bondi flow from a slowly rotating hot atmosphere

A supermassive black hole in the nucleus of an elliptical galaxy at the centre of a cool-core group or cluster of galaxies is immersed in hot gas. Bondi accretion should occur at a rate determined by the properties of the gas at the Bondi radius and the mass of the black hole. X-ray observations of massive nearby elliptical galaxies, including M87 in the Virgo cluster, indicate a Bondi accretion rate Mdot which roughly matches the total kinetic power of the jets, suggesting that there is a tight coupling between the jet power and the mass accretion rate. While the Bondi model considers non-rotating gas, it is likely that the external gas has some angular momentum, which previous studies have shown could decrease the accretion rate drastically. We investigate here the possibility that viscosity acts at all radii to transport angular momentum outward so that the accretion inflow proceeds rapidly and steadily. The situation corresponds to a giant Advection Dominated Accretion Flow (ADAF) which extends from beyond the Bondi radius down to the black hole. We find solutions of the ADAF equations in which the gas accretes at just a factor of a few less than Mdot. These solutions assume that the atmosphere beyond the Bondi radius rotates with a sub-Keplerian velocity and that the viscosity parameter is large, alpha~0.1. The infall time of the ADAF solutions is no more than a few times the free-fall time. Thus the accretion rate at the black hole is closely coupled to the surrounding gas, enabling tight feedback to occur. We show that jet powers of a few per cent of Mdot c^2 are expected if either a fraction of the accretion power is channeled into the jet or the black hole spin energy is tapped by a strong magnetic field pressed against the black hole by the pressure of the accretion flow.(Truncated)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraints on Cold Magnetized Shocks in Gamma-Ray Bursts

We consider a model in which the ultra-relativistic jet in a gamma-ray burst (GRB) is cold and magnetically accelerated. We assume that the energy flux in the outflowing material is partially thermalized via internal shocks or a reverse shock, and we estimate the maximum amount of radiation that could be produced in such magnetized shocks. We compare this estimate with the available observational data on prompt gamma-ray emission in GRBs. We find that, even with highly optimistic assumptions, the magnetized jet model is radiatively too inefficient to be consistent with observations. One way out is to assume that much of the magnetic energy in the post-shock, or even pre-shock, jet material is converted to particle thermal energy by some unspecified process, and then radiated. This can increase the radiative efficiency sufficiently to fit observations. Alternatively, jet acceleration may be driven by thermal pressure rather than magnetic fields. In this case, which corresponds to the traditional fireball model, sufficient prompt GRB emission could be produced either from shocks at a large radius or from the jet photosphere closer to the center.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Efficient Generation of Jets from Magnetically Arrested Accretion on a Rapidly Spinning Black Hole

We describe global, 3D, time-dependent, non-radiative, general-relativistic, magnetohydrodynamic simulations of accreting black holes (BHs). The simulations are designed to transport a large amount of magnetic flux to the center, more than the accreting gas can force into the BH. The excess magnetic flux remains outside the BH, impedes accretion, and leads to a magnetically arrested disc. We find powerful outflows. For a BH with spin parameter a = 0.5, the efficiency with which the accretion system generates outflowing energy in jets and winds is eta ~ 30%. For a = 0.99, we find eta ~ 140%, which means that more energy flows out of the BH than flows in. The only way this can happen is by extracting spin energy from the BH. Thus the a = 0.99 simulation represents an unambiguous demonstration, within an astrophysically plausible scenario, of the extraction of net energy from a spinning BH via the Penrose-Blandford-Znajek mechanism. We suggest that magnetically arrested accretion might explain observations of active galactic nuclei with apparent eta ~ few x 100%.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Equilibrium configurations from gravitational collapse

We develop here a new procedure within Einstein's theory of gravity to generate equilibrium configurations that result as the final state of gravitational collapse from regular initial conditions. As a simplification, we assume that the collapsing fluid is supported only by tangential pressure. We show that the equilibrium geometries generated by this method form a subset of static solutions to the Einstein equations, and that they can either be regular or develop a naked singularity at the center. When a singularity is present, there are key differences in the properties of stable circular orbits relative to those around a Schwarzschild black hole with the same mass. Therefore, if an accretion disk is present around such a naked singularity it could be observationally distinguished from a disk around a black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring the Spins of Accreting Black Holes

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observational Evidence for a Correlation Between Jet Power and Black Hole Spin

We show that the 5-GHz radio flux of transient ballistic jets in black hole binaries correlates with the dimensionless black hole spin parameter a* estimated via the continuum-fitting method. The data suggest that jet power scales either as the square of a* or the square of the angular velocity of the horizon. This is the first direct evidence that jets may be powered by black hole spin energy. The observed correlation validates the continuum-fitting method of measuring spin. In addition, for those black holes that have well-sampled radio observations of ballistic jets, the correlation may be used to obtain rough estimates of their spins.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Mass of the Black Hole in Cygnus X-1

Cygnus X-1 is a binary star system that is comprised of a black hole and a massive giant companion star in a tight orbit. Building on our accurate distance measurement reported in the preceding paper, we first determine the radius of the companion star, thereby constraining the scale of the binary system. To obtain a full dynamical model of the binary, we use an extensive collection of optical photometric and spectroscopic data taken from the literature. By using all of the available observational constraints, we show that the orbit is slightly eccentric (both the radial velocity and photometric data independently confirm this result) and that the companion star rotates roughly 1.4 times its pseudosynchronous value. We find a black hole mass of M =14.8\pm1.0 M_{\sun}, a companion mass of M_{opt}=19.2\pm1.9 M_{\sun}, and the angle of inclination of the orbital plane to our line of sight of i=27.1\pm0.8 deg.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Trigonometric Parallax of Cygnus X-1

We report a direct and accurate measurement of the distance to the X-ray binary Cygnus X-1, which contains the first black hole to be discovered. The distance of 1.86 (-0.11,+0.12) kpc was obtained from a trigonometric parallax measurement using the Very Long Baseline Array. The position measurements are also sensitive to the 5.6 d binary orbit and we determine the orbit to be clockwise on the sky. We also measured the proper motion of Cygnus X-1 which, when coupled to the distance and Doppler shift, gives the three-dimensional space motion of the system. When corrected for differential Galactic rotation, the non-circular (peculiar) motion of the binary is only about 21 km/s, indicating that the binary did not experience a large "kick" at formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variability in Blazars: Clues from PKS 2155-304

Rapid variability on a time scale much faster than the light-crossing time of the central supermassive black hole has been seen in TeV emission from the blazar PKS 2155-304. The most plausible explanation of this puzzling observation is that the radiating fluid in the relativistic jet is divided into a large number of sub-regions which move in random directions with relativistic speeds. We consider two versions of this "jets in a jet" model. In the first, the "subjets" model, stationary regions in the mean jet frame emit relativistic subjets that produce the observed radiation. The variability time scale is determined by the size of the sub-regions in the mean jet frame. This model, which is motivated by magnetic reconnection, has great difficulty explaining the observations in PKS 2155-304. In the alternate "turbulence" model, various sub-regions move relativistically in random directions and the variability time scale is determined by the size of these regions in their own comoving frames. This model fits the data much more comfortably. We consider collisions between TeV photons emitted from different sub-regions and find that, in both the subjets and turbulence models, the mean bulk Lorentz factor of the jet needs to be greater than 25 to avoid the pair catastrophe.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Black Hole Spin and the Radio Loud/Quiet Dichotomy of Active Galactic Nuclei

Radio loud active galactic nuclei (AGN) are on average 1000 times brighter in the radio band compared to radio quiet AGN. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin of the central black holes that power the radio-emitting jets. Using general relativistic magnetohydrodynamic simulations, we construct steady state axisymmetric numerical models for a wide range of black hole spins (dimensionless spin parameter 0.1 <= a <= 0.9999 and a variety of jet geometries. We assume that the total magnetic flux through the black hole horizon at radius r_H(a) is held constant. If the black hole is surrounded by a thin accretion disk, we find that the total black hole power output depends approximately quadratically on the angular frequency of the hole, P \propto Ω_H^2 \propto (a/r_H)^2. We conclude that, in this scenario, differences in the black hole spin can produce power variations of only a few tens at most. However, if the disk is thick such that the jet subtends a narrow solid angle around the polar axis, then the power dependence becomes much steeper, P \propto Ω_H^4 or even \propto Ω_H^6. Power variations of 1000 are then possible for realistic black hole spin distributions. We derive an analytic solution that accurately reproduces the steeper scaling of jet power with Ω_H, and we provide a numerical fitting formula that reproduces all our simulation results. We discuss other physical effects that might contribute to the observed radio loud/quiet dichotomy of AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetohydrodynamic Simulations of Gamma-Ray Burst Jets: Beyond the Progenitor Star

Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet's Lorentz factor γand opening angle Θ_j satisfies (γΘ_j) >> 1 at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give (γΘ_j) <~ 1, suggesting that MHD models may be inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet-envelope interface as a collimating rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to 10^4.5 times that of the central compact object, the jet achieves an asymptotic Lorentz factor γ~ 500 far outside the star and an asymptotic opening angle Θ_j ~ 0.04 rad ~ 2 deg, giving (γΘ_j) ~ 20. These values are consistent with observations of typical long-duration GRBs, and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations of Magnetized Disks Around Black Holes: Effects of Black Hole Spin, Disk Thickness, and Magnetic Field Geometry

The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Black Hole Mass Distribution in the Galaxy

We use dynamical mass measurements of 16 black holes in transient low-mass X-ray binaries to infer the stellar black hole mass distribution in the parent population. We find that the observations are best described by a narrow mass distribution at 7.8 +/- 1.2 Msolar. We identify a selection effect related to the choice of targets for optical follow-ups that results in a flux-limited sample. We demonstrate, however, that this selection effect does not introduce a bias in the observed distribution and cannot explain the absence of black holes in the 2-5 solar mass range. On the high mass end, we argue that the rapid decline in the inferred distribution may be the result of the particular evolutionary channel followed by low-mass X-ray binaries. This is consistent with the presence of high-mass black holes in the persistent, high-mass X-ray binary sources. If the paucity of low-mass black holes is caused by a sudden decrease of the supernova explosion energy with increasing progenitor mass, this would have observable implications for ongoing transient surveys that target core-collapse supernovae. Our results also have significant implications for the calculation of event rates from the coalescence of black hole binaries for gravitational wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Constant Inner-Disk Radius of LMC X-3: A Basis for Measuring Black Hole Spin

The black-hole binary system LMC X-3 has been observed by virtually every X-ray mission since the inception of X-ray astronomy. Among the persistent sources, LMC X-3 is uniquely both habitually soft and highly variable. Using a fully relativistic accretion-disk model, we analyze hundreds of spectra collected during eight X-ray missions that span 26 years. For a selected sample of 391 RXTE spectra we find that to within ~2 percent the inner radius of the accretion disk is constant over time and unaffected by source variability. Even considering an ensemble of eight X-ray missions, we find consistent values of the radius to within ~4-6 percent. Our results provide strong evidence for the existence of a fixed inner-disk radius. The only reasonable inference is that this radius is closely associated with the general relativistic innermost stable circular orbit (ISCO). Our findings establish a firm foundation for the measurement of black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Spin of the Black Hole in the Soft X-ray Transient A0620--00

During its yearlong outburst in 1975--76, the transient source A0620--00 reached an intensity of 50 Crab, an all-time record for any X-ray binary. The source has been quiescent since. We recently determined accurate values for the black hole mass, orbital inclination angle and distance. Building on these results, we have measured the radius of the inner edge of the accretion disk around the black hole primary by fitting its thermal continuum spectrum to our version of the relativistic Novikov-Thorne thin-disk model. We have thereby estimated the spin of the black hole. Although our spin estimate depends on a single high quality spectrum, which was obtained in 1975 by OSO-8, we are confident of our result because of the consistent values of the inner-disk radius that we have obtained for hundreds of observations of other sources: H1743-322, XTE J1550-564, and notably LMC X-3. We have determined the dimensionless spin parameter of the black hole to be a*=0.12+/-0.19, with a*<0.49 and a*>-0.59 at the 3 sigma level of confidence. This result takes into account all sources of observational and model-parameter uncertainties. Despite the low spin, the intensity and properties of the radio counterpart, both in outburst and quiescence, attest to the presence of a strong jet. If jets are driven by black hole spin, then current models indicate that jet power should be a steeply increasing function of a*. Consequently, the low spin of A0620--00 suggests that its jet may be disk-driven.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Spin of the Black Hole Microquasar XTE J1550-564 via the Continuum-Fitting and Fe-Line Methods

We measure the spin of XTE J1550-564 in two ways: by modelling the thermal continuum spectrum of the accretion disc, and independently by modeling the broad red wing of the reflection fluorescence Fe-K line. We find that the spin measurements conducted independently using both leading methods are in agreement with one another. For the continuum-fitting analysis, we use a data sample consisting of several dozen RXTE spectra, and for the Fe-K analysis, we use a pair of ASCA spectra from a single epoch. Our spin estimate for the black hole primary using the continuum-fitting method is -0.11 < a* < 0.71 (90 per cent confidence), with a most likely spin of a* = 0.34. In obtaining this result, we have thoroughly explored the dependence of the spin value on a wide range of model-dependent systematic errors and observational errors; our precision is limited by uncertainties in the distance and orbital inclination of the system. For the Fe-line method, our estimate of spin is a* = 0.55(+0.15,-0.22). Combining these results, we conclude that the spin of this black hole is moderate, a* = 0.49(+0.13,-0.20), which suggests that the jet of this microquasar is powered largely by its accretion disc rather than by the spin energy of the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A New Dynamical Model for the Black Hole Binary LMC X-1

We present a dynamical model of the high mass X-ray binary LMC X-1 based on high-resolution optical spectroscopy and extensive optical and near-infrared photometry. From our new optical data we find an orbital period of P=3.90917 +/- 0.00005 days. We present a refined analysis of the All Sky Monitor data from RXTE and find an X-ray period of P=3.9094 +/- 0.0008 days, which is consistent with the optical period. A simple model of Thomson scattering in the stellar wind can account for the modulation seen in the X-ray light curves. The V-K color of the star (1.17 +/- 0.05) implies A_V = 2.28 +/- 0.06, which is much larger than previously assumed. For the secondary star, we measure a radius of R_2 = 17.0 +/- 0.8 solar radii and a projected rotational velocity of V_rot*sin(i) = 129.9 +/- 2.2 km/s. Using these measured properties to constrain the dynamical model, we find an inclination of i = 36.38 +/- 1.92 deg, a secondary star mass of M_2 = 31.79 +/- 3.48 solar masses, and a black hole mass of 10.91 +/- 1.41 solar masses. The present location of the secondary star in a temperature-luminosity diagram is consistent with that of a star with an initial mass of 35 solar masses that is 5 Myr past the zero-age main sequence. The star nearly fills its Roche lobe (~90% or more), and owing to the rapid change in radius with time in its present evolutionary state, it will encounter its Roche lobe and begin rapid and possibly unstable mass transfer on a timescale of a few hundred thousand years.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Simple Comptonization Model

We present an empirical model of Comptonization for fitting the spectra of X-ray binaries. This model, "simpl", has been developed as a package implemented in XSPEC. With only two free parameters, simpl is competitive as the simplest empirical model of Compton scattering. Unlike other empirical models, such as the standard power-law model, simpl incorporates the basic physics of Compton scattering of soft photons by energetic coronal electrons. Using a simulated spectrum, we demonstrate that simpl closely matches the behavior of physical Comptonization models which consider the effects of optical depth, coronal electron temperature, and geometry. We present fits to RXTE spectra of the black-hole transient H1743-322 and a BeppoSAX spectrum of LMC X-3 using both simpl and the standard power-law model. A comparison of the results shows that simpl gives equally good fits and a comparable spectral index, while eliminating the troublesome divergence of the standard power-law model at low energies. Importantly, simpl is completely flexible and can be used self-consistently with any seed spectrum of photons. We show that simpl - unlike the standard power law - teamed up with diskbb (the standard model of disk accretion) gives results for the inner-disk radius that are unaffected by strong Comptonization, a result of great importance for the determination of black hole spin via the continuum-fitting method.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Host galaxies, clustering, Eddington ratios, and evolution of radio, X-ray, and infrared-selected AGNs

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Measuring Black Hole Spin via the X-ray Continuum Fitting Method: Beyond the Thermal Dominant State

All prior work on measuring the spins of stellar-mass black holes via the X-ray continuum-fitting method has relied on the use of weakly-Comptonized spectra obtained in the thermal dominant state. Using a self-consistent Comptonization model, we show that one can analyze spectra that exhibit strong power-law components and obtain values of the inner disk radius, and hence spin, that are consistent with those obtained in the thermal dominant state. Specifically, we analyze many RXTE spectra of two black hole transients, H1743-322 and XTE J1550-564, and we demonstrate that the radius of the inner edge of the accretion disk remains constant to within a few percent as the strength of the Comptonized component increases by an order of magnitude, i.e., as the fraction of the thermal seed photons that are scattered approaches 25%. We conclude that the continuum-fitting method can be applied to a much wider body of data than previously thought possible, and to sources that have never been observed to enter the thermal dominant state (e.g., Cyg X-1).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Event Horizon of Sagittarius A*

Black hole event horizons, causally separating the external universe from compact regions of spacetime, are one of the most exotic predictions of General Relativity (GR). Until recently, their compact size has prevented efforts to study them directly. Here we show that recent millimeter and infrared observations of Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, all but requires the existence of a horizon. Specifically, we show that these observations limit the luminosity of any putative visible compact emitting region to below 0.4% of Sgr A*'s accretion luminosity. Equivalently, this requires the efficiency of converting the gravitational binding energy liberated during accretion into radiation and kinetic outflows to be greater than 99.6%, considerably larger than those implicated in Sgr A*, and therefore inconsistent with the existence of such a visible region. Finally, since we are able to frame this argument entirely in terms of observable quantities, our results apply to all geometric theories of gravity that admit stationary solutions, including the commonly discussed f(R) class of theories.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Inferring the Inclination of a Black Hole Accretion Disk from Observations of its Polarized Continuum Radiation

Spin parameters of stellar-mass black holes in X-ray binaries are currently being estimated by fitting the X-ray continuum spectra of their accretion disk emission. For this method, it is necessary to know the inclination of the X-ray-producing inner region of the disk. Since the inner disk is expected to be oriented perpendicular to the spin axis of the hole, the usual practice is to assume that the black hole spin is aligned with the orbital angular momentum vector of the binary, and to estimate the inclination of the latter from ellipsoidal modulations in the light curve of the secondary star. We show that the inclination of the disk can be inferred directly if we have both spectral and polarization information on the disk radiation. The predicted degree of polarization varies from 0% to 5% as the disk inclination changes from face-on to edge-on. With current X-ray polarimetric techniques the polarization degree of a typical bright X-ray binary could be measured to an accuracy of 0.1% by observing the source for about 10 days. Such a measurement would constrain the disk inclination to within a degree or two and would significantly improve the reliability of black hole spin estimates. In addition, it would provide new information on the tilt between the black hole spin axis and the orbital rotation axis of the binary, which would constrain any velocity kicks experienced by stellar-mass black holes during their formation.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Generation of Type I X-ray Burst Oscillations by Unstable Surface Modes

The Rossi X-ray Timing Explorer has detected nearly coherent oscillations in the tails of type I X-ray bursts from 17 low-mass X-ray binaries. The oscillations are thought to be generated by brightness fluctuations associated with a surface mode on the rotating neutron star. The mechanism that drives the modes is, however, not understood, since the burning layer is stable to thermal perturbations. We show here via a linear perturbation analysis that, even under conditions when pure thermal perturbations are stable, nonradial surface modes may still be unstable by the epsilon mechanism. Specifically, we find that, if helium-burning reactions supply a reasonable fraction of the outgoing flux during burst decay, nonradial surface modes will grow in time. On the other hand, the same modes are likely to be stable in the presence of hydrogen burning via the rp-process. The results naturally explain why oscillations in the decay phase of type I X-ray bursts are detected only from short-duration bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Hydrogen-Triggered Type I X-ray Bursts in a Two-Zone Model

We use the two-zone model of Cooper & Narayan to study the onset and time evolution of hydrogen-triggered type I X-ray bursts on accreting neutron stars. At the lowest accretion rates, thermally unstable hydrogen burning ignites helium as well and produces a mixed hydrogen and helium burst. For somewhat higher accretion rates, thermally unstable hydrogen burning does not ignite helium and thus triggers only a weak hydrogen flash. The peak luminosities of weak hydrogen flashes are typically much lower than the accretion luminosity. These results are in accord with previous theoretical work. We find that a series of weak hydrogen flashes generates a massive layer of helium that eventually ignites in an energetic pure helium flash. Although previously conjectured, this is the first time such bursting behavior has been actually demonstrated in a theoretical model. For yet higher accretion rates, hydrogen burning is thermally stable and thus steadily generates a layer of helium that ultimately ignites in a pure helium flash. We find that, for a narrow range of accretion rates between the mixed hydrogen and helium burst and weak hydrogen flash regimes, unstable hydrogen burning ignites helium only after a short series of weak hydrogen flashes has generated a sufficiently deep layer of helium. These bursts have fluences that are intermediate between those of normal mixed hydrogen and helium bursts and energetic pure helium flashes.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A Two-Zone Model for Type I X-ray Bursts on Accreting Neutron Stars

We construct a two-zone model to describe H and He burning on the surface of an accreting neutron star and use it to study the triggering of type I X-ray bursts. Although highly simplified, the model reproduces all of the bursting regimes seen in the more complete global linear stability analysis of Narayan & Heyl (2003), including the regime of delayed mixed bursts. The results are also consistent with observations of type I X-ray bursts. At accretion rates Mdot < 0.1 Mdot_Edd, thermonuclear He burning via the well-known thin-shell thermal instability triggers bursts. As Mdot increases, however, the trigger mechanism evolves from the fast thermal instability to a slowly growing overstability involving both H and He burning. The competition between nuclear heating via the beta-limited CNO cycle and the triple-alpha process on the one hand, and radiative cooling via photon diffusion and emission on the other hand, drives oscillations with a period approximately equal to the H-burning timescale. If these oscillations grow, the gradually rising temperature at the base of the helium layer eventually provokes a thin-shell thermal instability and hence a delayed mixed burst. For Mdot > 0.25 Mdot_Edd, there is no instability or overstability, and there are no bursts. Nearly all other theoretical models predict that bursts should occur for all Mdot < Mdot_Edd, in conflict with both our results and observations. We suggest that this discrepancy arises from the assumed strength of the hot CNO cycle breakout reaction 15O(alpha,gamma)19Ne in these other models. That observations agree much better with the results of Narayan & Heyl and our two-zone model, both of which neglect breakout reactions, may imply that the true 15O(alpha,gamma)19Ne cross section is much smaller than assumed in previous investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

On the Physics of Type I X-ray Bursts on Accreting Neutron Stars at High Accretion Rates

We investigate the effect of the hot CNO cycle breakout reaction 15O(alpha,gamma)19Ne on the occurrence of type I X-ray bursts on accreting neutron stars. For f_rp <~ 0.1, where f_rp is a dimensionless factor by which we multiply the 15O(alpha,gamma)19Ne reaction rate of Caughlan & Fowler (1988), our model predicts that bursts should occur only for accretion rates below a critical value of approximately 0.3 times the Eddington limit. This agrees with observations. For larger values of f_rp, including the standard choice f_rp = 1, the model switches to a new regime in which bursts occur all the way up to roughly the Eddington limit. Since the latter regime disagrees with observations, we suggest that the true 15O(alpha,gamma)19Ne reaction rate is lower than usually assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Black Holes in Astrophysics

This article reviews the current status of black hole astrophysics, focusing on topics of interest to a physics audience. Astronomers have discovered dozens of compact objects with masses greater than 3 solar masses, the likely maximum mass of a neutron star. These objects are identified as black hole candidates. Some of the candidates have masses of 5 to 20 solar masses and are found in X-ray binaries, while the rest have masses from a million to a billion solar masses and are found in galactic nuclei. A variety of methods are being tried to estimate the spin parameters of the candidate black holes. There is strong circumstantial evidence that many of the objects have event horizons. Recent MHD simulations of magnetized plasma accreting on rotating black holes seem to hint that relativistic jets may be produced by a magnetic analog of the Penrose process.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Magnetically Arrested Disk: An Energetically Efficient Accretion Flow

We consider an accretion flow model originally proposed by Bisnovatyi-Kogan & Ruzmaikin (1974), which has been confirmed in recent 3D MHD simulations. In the model, the accreting gas drags in a strong poloidal magnetic field to the center such that the accumulated field disrupts the axisymmetric accretion flow at a relatively large radius. Inside the disruption radius, the gas accretes as discrete blobs or streams with a velocity much less than the free-fall velocity. Almost the entire rest mass energy of the gas is released as heat, radiation and mechanical/magnetic energy. Even for a non-rotating black hole, the efficiency of converting mass to energy is of order 50% or higher. The model is thus a practical analog of an idealized engine proposed by Geroch and Bekenstein.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Models of Galaxy Clusters with Thermal Conduction

We present a simple model of hot gas in galaxy clusters, assuming hydrostatic equilibrium and energy balance between radiative cooling and thermal conduction. For five clusters, A1795, A1835, A2199, A2390 and RXJ1347.5-1145, the model gives a good description of the observed radial profiles of electron density and temperature, provided we take the thermal conductivity $κ$ to be about 30% of the Spitzer conductivity. Since the required $κ$ is consistent with the recent theoretical estimate of Narayan & Medvedev (2001) for a turbulent magnetized plasma, we consider a conduction-based equilibrium model to be viable for these clusters. We further show that the hot gas is thermally stable because of the presence of conduction. For five other clusters, A2052, A2597, Hydra A, Ser 159-03 and 3C295, the model requires unphysically large values of $κ$ to fit the data. These clusters must have some additional source of heat, most likely an active galactic nucleus since all the clusters have strong radio galaxies at their centers. We suggest that thermal conduction, though not dominant in these clusters, may nevertheless play a significant role by preventing the gas from becoming thermally unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Neutrino trapping and accretion models for Gamma-Ray Bursts

Many models of Gamma Ray Bursts invoke a central engine consisting of a black hole of a few solar masses accreting matter from a disk at a rate of a fraction to a few solar masses per second. Popham et al. and Narayan et al. have shown that, for Mdot >~ 0.1 Msun/s, accretion proceeds via neutrino cooling and neutrinos can carry away a significant amount of energy from the inner regions of the disks. We improve on these calculations by including a simple prescription for neutrino transfer and neutrino opacities in such regions. We find that the flows become optically thick to neutrinos inside a radius R~6-40R_s for Mdot in the range of 0.1 -10 Msun/s, where R_s is the black hole Schwarzchild radius. Most of the neutrino emission comes from outside this region and, the neutrino luminosity stays roughly constant at a value L_ν ~ 10^{53} erg/s. We show that, for Mdot > 1 Msun/s, neutrinos are sufficiently trapped that energy advection becomes the dominant cooling mechanism in the flow. These results imply that neutrino annihilation in hyperaccreting black holes is an inefficient mechanism for liberating large amounts of energy. Extraction of rotational energy by magnetic processes remains the most viable mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Observational Prospects for Afterglows of Short Duration Gamma-ray Bursts

If the efficiency for producing $γ$-rays is the same in short duration ($\siml 2$ s) Gamma-Ray Bursts (GRBs) as in long duration GRBs, then the average kinetic energy of short GRBs must be $\sim 20$ times less than that of long GRBs. Assuming further that the relativistic shocks in short and long duration GRBs have similar parameters, we show that the afterglows of short GRBs will be on average 10--40 times dimmer than those of long GRBs. We find that the afterglow of a typical short GRB will be below the detection limit ($\siml 10 \microJy$) of searches at radio frequencies. The afterglow would be difficult to observe also in the optical, where we predict $R \simg 23$ a few hours after the burst. The radio and optical afterglow would be even fainter if short GRBs occur in a low-density medium, as expected in NS-NS and NS-BH merger models. The best prospects for detecting short-GRB afterglows are with early ($\siml 1$ day) observations in X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Deep radio observation of the gravitational lens candidate QSO2345+007

The double QSO2345+007 comprises two optical components separated by 7.1\ arcseconds and is the most prominent `dark matter' gravitational lens candidate. Despite being known for more than a decade, optical spectroscopy and imaging have been unable to determine whether this double QSO is a binary QSO or a gravitational lens system. In this note we report a deep VLA observation of this system, yielding a map with a noise level of 8.5\ts $μ$Jy per beam. We have a $4σ$ detection of a radio source within one arcsecond of the optical position of the brighter A-component of the QSO, but no significant detection of any radio counterpart of the B component. Given that the flux ratio in the optical waveband is $\sim$3--4, the gravitational lens hypothesis for this system would predict a radio flux of image B of $\lesssim 10$\ts $μ$Jy. The non-detection of the B component is thus consistent with, but does not prove, the lens interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Advection-Dominated Accretion Model of Sagittarius A$^*$ and Other Accreting Black Holes

Viscous rotating accretion flows around black holes become advection-dominated when the accretion rate $\dot M$ is sufficiently low. Most of the accretion energy in such flows is stored within the gas and advected radially inward. The temperature is therefore very high, and much of the radiation comes out in hard X-rays and $γ$-rays. We have constructed an advection-dominated accretion flow model for the Galactic Center source Sagittarius A$^*$. The model consists of a $7\times10^5M_\odot$ black hole accreting at $\dot M=1.2\times10^{-5}α\, M_\odot{\rm yr^{-1}}$, where $α$ is the usual viscosity parameter. The model spectrum fits the observations from radio to $γ$-rays quite well and explains the unusually low luminosity of the source. Since the model explicitly makes use of a horizon at the inner edge to swallow the advected energy, the success of the model strongly suggests that the central object in Sgr A$^*$ is a black hole. We further show that, if $α$ is not much smaller than unity, then advection-dominated models can be applied even to higher luminosity black holes. The existence of Low and High States in black hole X-ray binaries, and the abrupt transition between the two states, find a natural explanation. The models also explain the close similarity in the hard X-ray/$γ$-ray spectra of black hole X-ray binaries and active galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Models of Quiescent Black Hole Neutron Star Soft X-Ray Transients

When the mass accretion rate onto a black hole (BH) falls below a critical rate, ${\dot M}_{crit}\sim α^2 {\dot M}_{Edd}$, accretion can occur via a hot optically thin flow where most of the dissipated energy is advected inward. We present such an advection-dominated model for the soft X-ray transient (SXT) A0620-00. This source has a puzzlingly low X-ray luminosity in quiescence, $\sim 6\times 10^{30} ~{\rm erg\,s^{-1}}$, despite a relatively high mass accretion rate $\sim10^{-10} ~M_{\odot}{\rm yr^{-1}}$ deduced from its optical flux. The accreting gas makes a transition from a standard thin disc at large radii to an advection-dominated flow at small radii. The transition occurs when the effective temperature of the thin disc is $\sim 10^4$ K. Because of the very low accretion efficiency, $\sim 10^{-3}-10^{-4}$, in the inner flow, the model fits both the optical and X-ray data. We also present models for V404 Cyg and Nova Mus 1991 in quiescence. Quiescent neutron star (NS) transients are expected to appear very different from BH systems because the advected energy is re-radiated from the NS surface whereas a BH swallows the advected energy. We discuss models for NS SXTs.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Do Gamma-Ray Burst Sources Repeat

Following the discovery by Quashnock and Lamb (1993) of an apparent excess of $γ$-ray burst pairs with small angular separations, we reanalyze the angular distribution of the bursts in the BATSE catalogue. We find that in addition to an excess of close pairs, there is also a comparable excess of antipodal bursts, i.e pairs of bursts separated by about 180 degrees in the sky. Both excesses have only modest statistical significance. We reject the hypothesis put forward by Quashnock and Lamb that burst sources are repeaters, since it is obvious that this hypothesis does not predict an excess of antipodal coincidences. Lacking any physical model of bursts that can explain the antipodal pairs, we suggest that the two excesses seen in the data are either due to an unusual statistical fluctuation or caused by some unknown selection effect.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Hydrodynamics of Relativistic Fireballs

We examine analytically and numerically the evolution of a relativistic fireball. We show that, after an early rearrangement phase, most of the matter and energy in the fireball is concentrated within a narrow shell. The shell propagates at nearly the speed of light, with a frozen radial profile and according to a simple set of scaling laws. The spectrum of the escaping radiation is harder at early times and softer later on. The results are applicable to models of $γ$-ray bursts. Submitted to Monthly Notices of the RAS Figures available on request from the authors

0 citations0 reviewsNo code linkedOpen access