by Hayasaki, Kimitake and Yagi, Kent and Tanaka, Takahiro and Mineshige, Shin
10 pages, 2 figures, submitted to ApJ

When binary black holes are embedded in a gaseous environment, a rotating disk surrounding them, the so-called circumbinary disk, will be formed. The binary exerts a gravitational torque on the circumbinary disk and thereby the orbital angular momentum is transferred to it, while the angular momentum of the circumbinary disk is transferred to the binary through the mass accretion. The binary undergoes an orbital decay due to both the gravitational wave emission and the binary-disk interaction. This causes the phase evolution of the gravitational wave signal. The precise measurement of the gravitational wave phase thus may provide information regarding the circumbinary disk. In this paper, we assess the detectability of the signature of the binary-disk interaction using the future space-borne gravitational wave detectors such as DECIGO and BBO by the standard matched filtering analysis. We find that the effect of the circumbinary disk around binary black holes in the mass range $latex 6M_sun\le{M}\lesssim3\times10^3M_sun$ is detectable at a statistically significant level in five year observation, provided that gas accretes onto the binary at a rate greater than $latex \dot{M}\sim1.4\times10^{17} [gs^{-1}] j^{-1}(M/10M_sun)^{33/23}$ with 10% mass-to-energy conversion efficiency, where j represents the efficiency of the angular momentum transfer from the binary to the circumbinary disk. We show that $latex O(0.1)$ coalescence events are expected to occur in sufficiently dense molecular clouds in five year observation. We also point out that the circumbinary disk is detectable, even if its mass at around the inner edge is by over 10 orders of magnitude less than the binary mass.

by Harko, Tiberiu and Mocanu, Gabriela
10 pages, 8 figures, accepted for publication in MNRAS

We analyze the general relativistic oscillations of thin accretion disks around compact astrophysical objects interacting with the surrounding medium through non-gravitational forces. The interaction with the external medium (a thermal bath) is modeled via a friction force, and a random force, respectively. The general equations describing the stochastically perturbed disks are derived by considering the perturbations of trajectories of the test particles in equatorial orbits, assumed to move along the geodesic lines. By taking into account the presence of a viscous dissipation and of a stochastic force we show that the dynamics of the stochastically perturbed disks can be formulated in terms of a general relativistic Langevin equation. The stochastic energy transport equation is also obtained. The vertical oscillations of the disks in the Schwarzschild and Kerr geometries are considered in detail, and they are analyzed by numerically integrating the corresponding Langevin equations. The vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases.

by Lousto, Carlos O. and Zlochower, Yosef and Dotti, Massimo and Volonteri, Marta
17 pages, 10 tables, 14 figures, revtex 4

We explore the newly discovered “hangup-kick” effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.

by King, Andrew
to appear in proceedings of `The Golden Age of Cataclysmic Variables’, Memorie Societa’ Astronomica Italiana, 2012 (F. Giovannelli and L. Sabau-Graziati eds.)

I briefly review the progress of accretion disc theory since the seminal paper of Shakura and Sunyaev.

by Zubovas, Kastytis and King, Andrew R.
6 pages, 1 figure, 2 tables; accepted for publication in ApJ Letters

It is widely suspected that AGN activity ultimately sweeps galaxies clear of their gas. We work out the observable properties required to achieve this. Large-scale AGN-driven outflows should have kinetic luminosities $latex \sim \eta\le/2 \sim 0.05\le$ and momentum rates $latex \sim 20\le/c$, where $latex \le$ is the Eddington luminosity of the central black hole and $latex \eta\sim 0.1$ its radiative accretion efficiency. This creates an expanding two-phase medium in which molecular species coexist with hot gas, which can persist after the central AGN has switched off. This picture predicts outflow velocities $latex \sim 1000 – 1500$ km\,s$latex ^{-1}$ and mass outflow rates up to $latex 4000 \msun\,{\rm yr}^{-1}$ on kpc scales, fixed mainly by the host galaxy velocity dispersion (or equivalently black hole mass). All these features agree with those of outflows observed in galaxies such as Mrk231. This strongly suggests that AGN activity is what sweeps galaxies clear of their gas on a dynamical timescale and makes them red and dead. We suggest future observational tests of this picture.

by Burkert, Andreas and Schartmann, Mark and Alig, Christian and Gillessen, Stefan and Genzel, Reinhard and Fritz, Tobias and Eisenhauer, Frank
22 pages, 13 figures, submitted to ApJ

The origin, structure and evolution of the small gas cloud, G2, is investigated, that is on an orbit almost straight into the Galactic central supermassive black hole (SMBH). G2 is a sensitive probe of the hot accretion zone of Sgr A*, requiring gas temperatures and densities that agree well with models of captured shock-heated stellar winds. Its mass is equal to the critical mass below which cold clumps would be destroyed quickly by evaporation. Its mass is also constrained by the fact that at apocenter its sound crossing timescale was equal to its orbital timescale. Our numerical simulations show that the observed structure and evolution of G2 can be well reproduced if it formed in pressure equilibrium with the surrounding in 1995 at a distance from the SMBH of 7.6e16 cm. If the cloud would have formed at apocenter in the ‘clockwise’ stellar disk as expected from its orbit, it would be torn into a very elongated spaghetti-like filament by 2011 which is not observed. This problem can be solved if G2 is the head of a larger, shell-like structure that formed at apocenter. Our numerical simulations show that this scenario explains not only G2’s observed kinematical and geometrical properties but also the Br_gamma observations of a low surface brightness gas tail that trails the cloud. In 2013, while passing the SMBH G2 will break up into a string of droplets that within the next 30 years mix with the surrounding hot gas and trigger cycles of AGN activity.

by Nixon, Chris and King, Andrew
8 pages, 9 figures. Accepted for publication in MNRAS

We simulate the viscous evolution of an accretion disc around a spinning black hole. In general any such disc is misaligned, and warped by the Lense-Thirring effect. Unlike previous studies we use effective viscosities constrained to be consistent with the internal fluid dynamics of the disc. We find that nonlinear fluid effects, which reduce the effective viscosities in warped regions, can promote the breaking of the disc into two distinct planes. This occurs when the Shakura & Sunyaev dimensionless viscosity parameter alpha is ~ 45 degrees. The break can be a long-lived feature, propagating outwards in the disc on the usual alignment timescale, after which the disc is fully co- or counter-aligned with the hole. Such a break in the disc may be significant in systems where we know the inclination of the outer accretion disc to the line of sight, such as some X-ray binaries: the inner disc, and so any jets, may be noticeably misaligned with respect to the orbital plane.

by Barai, Paramita and Proga, Daniel and Nagamine, Kentaro
23 pages, 11 figures, submitted. Uploaded version contains low-resolution color figures. Version with high-resolution figures can be found at: http://www.physics.unlv.edu/~barai/AllPages/Images-Movies/BHaccr_MultiPhase.pdf

(Abridged) We investigate non-spherical behavior of gas accreting onto a central supermassive black hole performing simulations using the SPH code GADGET-3 including radiative cooling and heating by the central X-ray source. As found in earlier 1D studies, our 3D simulations show that the accretion mode depends on the X-ray luminosity (L_X) for a fixed density at infinity and accretion efficiency. In the low L_X limit, gas accretes in a stable, spherically symmetric fashion. In the high L_X limit, the inner gas is significantly heated up and expands, reducing the central mass inflow rate. The expanding gas can turn into a strong enough outflow capable of expelling most of the gas at larger radii. For some intermediate L_X, the accretion flow becomes unstable developing prominent non-spherical features, the key reason for which is thermal instability (TI) as shown by our analyses. Small perturbations of the initially spherically symmetric accretion flow that is heated by the intermediate L_X quickly grow to form cold and dense clumps surrounded by overheated low density regions. The cold clumps continue their inward motion forming filamentary structures; while the hot infalling gas slows down because of buoyancy and can even start outflowing through the channels in between the filaments. We found that the ratio between the mass inflow rates of the cold and hot gas is a dynamical quantity depending on several factors: time, spatial location, and L_X; and ranges between 0 and 4. We briefly discuss astrophysical implications of such TI-driven fragmentation of accreting gas on the formation of clouds in narrow and broad line regions of AGN, the formation of stars, and the observed variability of the AGN luminiosity.

by Taylor, Paul A. and Miller, John C.
14 pages, 9 figures, submitted to MNRAS

A commonly cited drawback of SPH is the introduction of spurious shear viscosity by the artificial viscosity term in situations involving rotation. Existing approaches for understanding its effect include approximative analytic formulae and disc-averaged behaviour in specific ring-spreading simulations, based on the kinematic contribution of the artificial viscosity. In this work, we have developed a simple, general technique for evaluating the local effect of artificial viscosity directly from the entropic function of each SPH particle. This is simple and quick to implement, and it allows a detailed characterization of its effects as a function of position. Several advantages of this local method are discussed, including its ease in evaluation, its greater accuracy and its broad applicability to arbitary flow geometries and equations of state. Here, we apply this approach to various disc flows, including simulations which implement the commonly-used Balsara switch. Comparisons with existing analytic estimates are made, and examples of quantifying explicit dependencies of the effective viscosity in terms of SPH and flow parameters are given. Additionally, a method for the initial placement of SPH particles which reduces numerical fluctuations is discussed and utilised.

by Shankar, Francesco and Weinberg, David H. and Miralda-Escude’, Jordi
27 pages, 16 Figures. Submitted to MNRAS

We develop semi-empirical models of the supermassive black hole and active galactic nucleus (AGN) populations, which incorporate the black hole growth implied by the observed AGN luminosity function assuming a radiative efficiency \epsilon, and a distribution of Eddington ratios \lambda. By generalizing these continuity-equation models to allow a distribution P(\lambda|mbh,z) we are able to draw on constraints from observationally estimated P(\lambda) distributions and active galaxy fractions while accounting for the luminosity thresholds of observational samples. We consider models with a Gaussian distribution of log \lambda, and Gaussians augmented with a power-law tail to low \lambda. Within our framework, reproducing the high observed AGN fractions at low redshift requires a characteristic Eddington ratio \lambda_c that declines at late times, and matching observed Eddington ratio distributions requires P(\lambda) that broadens at low redshift. To reproduce the observed increase of AGN fraction with black hole or galaxy mass, we also require a \lambda_c that decreases with increasing black hole mass, reducing the AGN luminosity associated with the most massive black holes. Finally, achieving a good match to the high mass end of the local black hole mass function requires an increased radiative efficiency at high black hole mass. We discuss the potential impact of black hole mergers or a \lambda-dependent bolometric correction, and we compute evolutionary predictions for black hole and galaxy specific accretion rates. Despite the flexibility of our framework, no one model provides a good fit to all the data we consider.

by Garcia-Senz, Domingo and Cabezon, Ruben M. and Escartin, Jose Antonio
15 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

In this paper we develop and check a fully conservative SPH scheme based on a tensor formulation which can be applied to simulate astrophysical systems. In the proposed scheme derivatives are calculated from an integral expression which leads to a tensor, rather than vectorial, estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the interpolation of physical magnitudes, leading to a set of conservative equations which looks similar to those of standard SPH. The resulting scheme was checked using a variety of well known tests, all of them simulated in two dimensions. An application of the proposed tensor method to astrophysics was also discussed by simulating the stability of a sun-like polytrope calculated in three dimensions.

by Rice, W. K. M. and Forgan, D. H. and Armitage, P. J.
9 pages, 11 figures, MNRAS in press

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

by Laibe, Guillaume and Price, Daniel J.
Accepted for publication in MNRAS

In a companion paper (Laibe & Price 2011b), we have presented an algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momentum in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non- linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including i) the analytic solutions of the dustybox problem and ii) solutions of the dustywave, dustyshock, dustysedov and dustydisc obtained with explicit timestepping. We find that the implicit method is 1- 10 times faster than the explicit temporal integration when the ratio r between the the timestep and the drag stopping time is 1 < r < 1000.

by Laibe, Guillaume and Price, Daniel J.
Accepted for publication in MNRAS

We present a new algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH), systematically addressing a number of key issues including the generalised SPH density estimate in multi-fluid systems, the consistent treatment of variable smoothing length terms, finite particle size, time step stability, thermal coupling terms and the choice of kernel and smoothing length used in the drag operator. We find that using double-hump shaped kernels improves the accuracy of the drag interpolation by a factor of several hundred compared to the use of standard SPH bell-shaped kernels, at no additional computational expense. In order to benchmark our algorithm, we have developed a comprehensive suite of standardised, simple test problems for gas and dust mixtures: dustybox, dustywave, dustyshock, dustysedov and dustydisc, the first three of which have known analytic solutions. We present the validation of our algorithm against all of these tests. In doing so, we show that the spatial resolution criterion \Delta < cs ts is a necessary condition in all gas+dust codes that becomes critical at high drag (i.e. small stopping time ts) in order to correctly predict the dynamics. Implicit timestepping and the implementation of realistic astrophysical drag regimes are addressed in a companion paper.

by Krolik, Julian H. and Piran, Tsvi

Tidal disruption of main sequence stars by black holes has generally been thought to lead to a signal dominated by UV emission. If, however, the black hole spins rapidly and the poloidal magnetic field intensity on the black hole horizon is comparable to the inner accretion disk pressure, a powerful jet may form whose luminosity can easily exceed the thermal UV luminosity. When the jet beam points at Earth, its non-thermal luminosity can dominate the emitted spectrum. The thermal and non-thermal components decay differently with time. In particular, the thermal emission should remain roughly constant for a significant time after the period of maximum accretion, beginning to diminish only after a delay, whereas after the peak accretion rate, the non-thermal jet emission decays, but then reaches a plateau. When the newly-found flare source Swift J2058 is analyzed in terms of this model, it is found to be consistent with an event in which a main sequence solar-type star is disrupted by a black hole of mass at least $latex \sim 10^7 M_{\odot}$. Swift may have already observed the beginning of the flat phase in the non-thermal emission from this source. Optical photometry over the first $latex \simeq 40$ d of this flare is also consistent with this picture, but there is a large uncertainty in the bolometric correction. We suggest that future searches for main sequence tidal disruptions use methods sensitive to jet radiation as well as to thermal UV radiation.

by Willott, Chris J.
4 pages, 2 figures, ApJ Letters, in press

It has been claimed that there is a large population of obscured, accreting black holes at high-redshift and that the integrated black hole density at z=6 as inferred from X-ray observations is ~100 times greater than inferred from optical quasars. I have performed a stacking analysis of very deep Chandra X-ray data at the positions of photometrically-selected z=6 galaxy candidates. It is found that there is no evidence for a stacked X-ray signal in either the soft (0.5-2 keV) or hard (2-8 keV) X-ray bands. Previous work which reported a significant signal is affected by an incorrect method of background subtraction which underestimates the true background within the target aperture. The puzzle remains of why the z=6 black hole mass function has such a flat slope and a low normalization compared to the stellar mass function.

by Chen, Songbai and Jing, Jiliang
13 pages, 5 figures. References added, Expanded discussion of the marginally stable orbit and its consequence. arXiv admin note: substantial text overlap with arXiv:1106.5183

We study the accretion process in the thin disk around a rotating non-Kerr black hole with a deformed parameter and an unbound rotation parameter. Our results show that the presence of the deformed parameter $latex \epsilon$ modifies the standard properties of the disk. For the case in which the black hole is more oblate than a Kerr black hole, the larger deviation leads to the smaller energy flux, the lower radiation temperature and the fainter spectra luminosity in the disk. For the black hole with positive deformed parameter, we find that the effect of the deformed parameter on the disk becomes more complicated. It depends not only on the rotation direction of the black hole and the orbit particles, but also on the sign of the difference between the deformed parameter $latex \epsilon$ and a certain critical value $latex \epsilon_{c}$. These significant features in the mass accretion process may provide a possibility to test gravity in the strong field regime in future astronomical observations.

by Park, KwangHo and Ricotti, Massimo
27 pages, 10 figures, submitted to ApJ, for associated mpeg files, see http://www.astro.umd.edu/~kpark/research.html

In this paper, the second of a series on radiation-regulated accretion onto black holes(BHs) from galactic scales, we focus on the effects that radiation pressure and angular momentum of the gas have on the periodic and short-lived luminosity bursts found when thermal pressure of the ionized sphere around the BH regulates the accretion rate. Our simulations focus on intermediate-mass BH, but we derive general scaling relationships that are solutions of the classic Bondi problem when radiation feedback is considered. We find that for ambient gas densities(n) exceeding a critical value n (5×10^6 cm^{-3})/M_2, where M_2 is the mass of the BH in units of 100 solar masses, the period of the oscillations decreases rapidly and the duty cycle increases from 6% to 50%. However, the maximum and mean accretion rates become Eddington limited only if n>n_Edd n_cr/T_4 where T_4 is the ambient gas temperature in units of 10^4 K. In the sub-Eddington regime, the mean accretion rate onto BH is about 1% T_4^{2.5} of the Bondi rate, thus is proportional to the thermal pressure of the ambient medium. The period of the oscillations coincides with depletion time scale of the gas inside the hot ionized bubble surrounding the BH. For n~n_cr accretion onto the BH becomes the dominant gas depletion mechanism, explaining the rapid decrease of the period of the oscillations and increasing the duty cycle. Angular momentum of the accreting gas produces a time delay between the accretion rate near the sonic point and the luminosity output. Generally, assuming reasonable values of the time delay calculated using an alpha-model for a thin disk, angular momentum does not affect significantly the accretion rate and period of the oscillations.

by Dotti, M. and Sesana, A. and Decarli, R.
4 Figures. Accepted for publication in Advances in Astronomy

The study of the dynamical evolution of massive black hole pairs in mergers is crucial in the context of a hierarchical galaxy formation scenario. The timescales for the formation and the coalescence of black hole binaries are still poorly constrained, resulting in large uncertainties in the expected rate of massive black hole binaries detectable in the electromagnetic and gravitational wave spectra. Here we review the current theoretical understanding of the black hole pairing in galaxy mergers, with a particular attention to recent developments and open issues. We conclude with a review of the expected observational signatures of massive binaries, and of the candidates discussed in literature to date.

by Begelman, Mitchell C.
13 pages, 8 figures, to appear in Monthly Notices of the Royal Astronomical Society

We reformulate the adiabatic inflow-outflow (ADIOS) model for radiatively inefficient accretion flows, treating the inflow and outflow zones on an equal footing. For purely adiabatic flows (i.e., with no radiative losses), we show that the mass flux in each zone must satisfy Mdot ~ R^n with n=1, in contrast to previous work in which 0<n< 1 is a free parameter but in rough agreement with numerical simulations. We also demonstrate that the resulting two-zone ADIOS models are not dynamically self-consistent without the introduction of an energy source close in to the central regions of the flow; we identify this with the energy liberated by accretion. We explore the parameter space of non-radiative models and show that both powerful winds and gentle breezes are possible. When small radiative losses (with fixed efficiency) are included, any centrally injected energy flux is radiated away and the system reverts to a power-law behavior with n < 1, where n falls in a small range determined by the fractional level of radiative losses. We also present an ADIOS model for hypercritical (super-Eddington) disk accretion, in which the radiative losses are closely related to the flow geometry. We suggest that hyperaccretion can lead to either winds or breezes.

by Wang, Yan and Li, Xiang-Dong
22 pages, 8 figures. Accepted for publication in ApJ

If an accretion disk around a black hole is illuminated by hard X-rays from non-thermal coronae, fluorescent iron lines will be emitted from the inner region of the accretion disk. The emission line profiles will show a variety of strong field effects, which may be used as a probe of the spin parameter of the black hole and the structure of the accretion disk. In this paper we generalize the previous relativistic line profile models by including both the black hole spinning effects and the non-axisymmetries of warped accretion disks. Our results show different features from the conventional calculations for either a flat disk around a Kerr black hole or a warped disk around a Schwarzschild black hole by presenting, at the same time, multiple peaks, rather long red tails and time variations of line profiles with the precession of the disk. We show disk images as seen by a distant observer, which are distorted by the strong gravity. Although we are primarily concerned with the iron K-shell lines in this paper, the calculation is general and is valid for any emission lines produced from a warped accretion disk around a black hole.

by Valtonen, Mauri J. and Villforth, Carolin and Wiik, Kaj
to appear in Monthly Notices of Royal Astronomical Society

We study the variation of the optical polarization angle in the blazar OJ287 and compare it with the precessing binary black hole model with a ‘live’ accretion disk. First, a model of the variation of the jet direction is calculated, and the main parameters of the model are fixed by the long term optical brightness evolution. Then this model is compared with the variation of the parsec scale radio jet position angle in the sky. Finally, the variation of the polarization angle is calculated using the same model, but using a magnetic field configuration which is at a constant angle relative to the optical jet. It is found that the model fits the data reasonably well if the field is almost parallel to the jet axis. This may imply a steady magnetic field geometry, such as a large-scale helical field.

by Blandford, Roger D.
10 pages, 10 figures, to be submitted to MNRAS

In this paper we consider Roche accretion in an Extreme Mass-Ratio Inspiral (EMRI) binary system formed by a star orbiting a massive black hole. The ultimate goal is to detect the mass and spin of the black hole and provide a test of general relativity in the strong-field regime from the resultant quasi-periodic signals. Before accretion starts, the stellar orbit is presumed to be circular and equatorial, and shrinks due to gravitational radiation. New fitting formulae are presented for the inspiral time and the radiation-reaction torque in the relativistic regime. If the inspiralling star fills its Roche lobe outside the Innermost Stable Circular Orbit (ISCO) of the hole, gas will flow through the inner Lagrange point (L1) to the hole. We give new relativistic interpolation formulae for the volume enclosed by the Roche lobe. If this mass-transfer happens on a time scale faster than the thermal time scale but slower than the dynamical time scale, the star will evolve adiabatically, and, in most cases, will recede from the hole while filling its Roche lobe. We calculate how the stellar orbital period and mass-transfer rate will change through the “Roche evolution” for various types of stars in the relativistic regime. We envisage that the mass stream eventually hits the accretion disc, where it forms a hot spot orbiting the hole and may ultimately modulate the luminosity with the stellar orbital frequency. The observability of such a modulation is discussed along with a possible interpretation of an intermittent 1 hour period in the X-ray emission of RE J1034+396.

by Alexander, R. D. and Smedley, S. L. and Nayakshin, S. and King, A. R.
7 pages, 3 figures. Accepted for publication in MNRAS

We study the dynamical evolution of stars and gas close to the centre of the Milky Way. Any plausible means of forming the young stars observed at the Galactic Centre leaves behind a residual gas disc at ~0.01pc radii. We show that the combined effects of viscous accretion and gravitational interactions with stars do not remove the residual gas efficiently, and that a substantial gas disc, interior to the stellar disc, persists for >10Myr after the stars form. Since no such disc is currently seen at the Galactic Centre we argue that it has been accreted by the super-massive black hole. This scenario offers an attractive connection between nuclear star formation and black hole feeding, and we suggest that the “missing” gas may have been used to power Sgr A*.

by Li, Yuexing
8 pages, 4 figures. Submitted to ApJ

The origin of the supermassive black holes that power the most distant quasars observed is largely unknown. One hypothesis is that they grew rapidly from intermediate-mass seeds (~100 M_sun) left by the first stars. However, some previous studies argued that accretion onto these black holes was too low to build up the mass due to strong suppression by radiative feedback. Here, we re-exam the accretion process of such a black hole embedded in a primordial gas cloud, by considering a wide range of physical and numerical parameters not explored before. We find that, while radiative heating and pressure indeed suppress accretion effectively, self-gravity of the gas eventually overcomes the feedback effects and boosts the accretion to the Eddington rate after one free-fall timescale of the cloud. Moreover, for a given black hole mass, there exists a critical density above which the accretion can reach Eddington limit. Furthermore, we find a universal correlation between black hole accretion rate and ambient gas density, which may serve as a realistic recipe for black hole growth in simulations.

by Gil-Merino, Rodrigo and Goicoechea, Luis J. and Shalyapin, Vyacheslav N. and Braga, Vittorio F.
21 pages, 7 text pages, 13 figures, 2 tables, accepted by The Astrophysical Journal

We present the results of our X-ray, UV and optical monitoring campaign of the first gravitationally lensed AGN from late 2009 to mid 2010. The trailing (B) image of the AGN 0957+561 shows the intrinsic continuum variations that were predicted in advance based on observations of the leading (A) image in the gr optical bands. This multiwavelength variability of the B image allows us to carry out a reverberation mapping analysis in the radio-loud AGN 0957+561 at redshift z = 1.41. We find that the U-band and r-band light curves are highly correlated with the g-band record, leading and trailing it by 3 +/- 1 days (U band) and 4 +/- 1 days (r band). These 1-sigma measurements are consistent with a scenario in which flares originated in the immediate vicinity of the supermassive black hole are thermally reprocessed in a standard accretion disk at about 10-20 Schwarzschild radii from the central dark object. We also report that the light curve for the X-ray emission with power-law spectrum is delayed with respect to those in the Ugr bands by about 32 days. Hence, the central driving source can not be a standard corona emitting the observed power-law X-rays. This result is also supported by X-ray reprocessing simulations and the absence of X-ray reflection features in the spectrum of 0957+561. We plausibly interpret the lack of reflection and the 32-day delay as evidence for a power-law X-ray source in the base of the jet at a typical height of about 200 Schwarzschild radii. A central EUV source would drive the variability of 0957+561.

by Springel, Volker
35 pages, 10 figures; invited review for the volume “Tessellations in the Sciences: Virtues, Techniques and Applications of Geometric Tilings”, eds. R. van de Weijgaert, G. Vegter, J. Ritzerveld and V. Icke, Springer (accepted)

At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be adaptively refined during a calculation. Particle-based methods on the other hand discretize the mass instead of the volume, yielding an approximately Lagrangian approach. It is also possible to achieve Lagrangian behavior in mesh-based methods if the mesh is allowed to move with the flow. However, such approaches have often been fraught with substantial problems related to the development of irregularity in the mesh topology. Here we describe a novel scheme that eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order Godunov scheme with an exact Riemann solver. A particularly powerful feature of the approach is that the mesh-generating points can in principle be moved arbitrarily. If they are given the velocity of the local flow, a highly accurate Lagrangian formulation of continuum hydrodynamics is obtained that is free of mesh distortion problems, while it is at the same time fully Galilean-invariant, unlike ordinary Eulerian codes. We describe the formulation and implementation of our new Voronoi-based hydrodynamics, and we discuss a number of illustrative test problems that highlight its performance in practical applications.

by Springel, Volker
43 pages, 11 figures, invited review that appeared in ARA&A

This review discusses Smoothed Particle Hydrodynamics (SPH) in the astrophysical context, with a focus on inviscid gas dynamics. The particle-based SPH technique allows an intuitive and simple formulation of hydrodynamics that has excellent conservation properties and can be coupled to self-gravity easily and highly accurately. The Lagrangian character of SPH allows it to automatically adjust its resolution to the clumping of matter, a property that makes the scheme ideal for many applications in astrophysics, where often a large dynamic range in density is encountered. We discuss the derivation of the basic SPH equations in their modern formulation, and give an overview about extensions of SPH developed to treat physics such as radiative transfer, thermal conduction, relativistic dynamics or magnetic fields. We also briefly describe some of the most important applications areas of SPH in astrophysical research. Finally, we provide a critical discussion of the accuracy of SPH for different hydrodynamical problems, including measurements of its convergence rate for important classes of problems.

by Valtonen, M. J. and Mikkola, S. and Lehto, H. J. and Gopakumar, A. and Hudec, R. and Polednikova, J.
14 pages, 14 figures

We examine the ability to test the black hole no-hair theorem at the 10% level in this decade using the binary black hole in OJ287. In the test we constrain the value of the dimensionless parameter q that relates the scaled quadrupole moment and spin of the primary black hole: q2 = -q 2 . At the present we can say that q = 1 \pm 0.3 (one), in agreement with General Relativity and the no-hair theorems. We demonstrate that this result can be improved if more observational data is found in historical plate archives for the 1959 and 1971 outbursts. We also show that the predicted 2015 and 2019 outbursts will be crucial in improving the accuracy of the test. Space-based photometry is required in 2019 July due the proximity of OJ287 to the Sun at the time of the outburst. The best situation would be to carry out the photometry far from the Earth, from quite a different vantage point, in order to avoid the influence of the nearby Sun. We have considered in particular the STEREO space mission which would be ideal if it has a continuation in 2019 or LORRI on board the New Horizons mission to Pluto.

by Zanotti, Olindo
6 pages, 4 figures. Accepted by New Astronomy

We show the results of two dimensional general relativistic inviscid and isothermal hydrodynamical simulations comparing the behavior of co-rotating (with respect to the black hole rotation) and counter-rotating circumbinary quasi-Keplerian discs in the post merger phase of a supermassive binary black hole system. While confirming the spiral shock generation within the disc due to the combined effects of mass loss and recoil velocity of the black hole, we find that the maximum luminosity of counter-rotating discs is a factor ~(2-12) higher than in the co-rotating case, depending on the spin of the black hole. On the other hand, the luminosity peak happens ~10 days later with respect to the co-rotating case, for a binary with a total mass M~10^6 M_\odot. Although the global dynamics of counter-rotating discs in the post merger phase of a merging event is very similar to that for co-rotating discs, an important difference has been found. In fact, increasing the spin of the central black hole produces more luminous co-rotating discs while less luminous counter-rotating ones.

by Tchekhovskoy, Alexander and Narayan, Ramesh and McKinney, Jonathan C.
5 pages, 2 figures, MNRAS, submitted

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 BH can swallow. 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. Thus, the gravitational mass of the BH decreases with time. This 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 AGN with apparent eta ~ few x 100%.

by Sadowski, A.
Ph.D. thesis written under the supervision of Prof. Marek Abramowicz and defended on July 1, 2011 at the Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences; 206 pages

In this thesis, I study hydrodynamical models of slim accretion disks — advective, optically thick disks which generalize the standard models of radiatively efficient thin disks to all accretion rates. I start with a general introduction to the theory of accretion onto compact objects. It is followed by a derivation of the commonly-used standard models of thin disks. In the subsequent section I introduce the equations describing slim disks, explain the numerical methods I used to solve them and discuss properties of such solutions. I also give a general derivation of non-stationary equations and present the time evolution of thermally unstable accretion disks. I introduce a state-of-the-art approach coupling the radial and vertical structures of an advective accretion disk and discuss the improvements it brings to vertically-averaged solutions. I also present a numerical model of self-illuminated slim accretion disks. Finally, I present and discuss applications of slim accretion disks: estimating of spin of the central black hole in LMC X-3 through X-ray continuum fitting basing on high-luminosity data, spinning-up of black holes by super-critical accretion flows and normalizing of magnetohydrodynamical global simulations.

by Rice, W. K. M. and Armitage, P. J. and Mamatsashvili, G. R. and Lodato, G. and Clarke, C. J.
MNRAS, in press

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress – which can be represented as an effective viscous alpha – out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

by Tanaka, Takamitsu and Haiman, Zoltán and Menou, Kristen
16 pages with 5 figures, submitted to MNRAS

Pulsar timing arrays (PTAs) are expected to detect gravitational waves (GWs) from individual low-redshift (z10^9 Msun) black hole (SMBH) binaries with orbital periods of approx. 0.1 – 10 yrs. Identifying the electromagnetic (EM) counterparts of these sources would provide confirmation of putative direct detections of GWs, present a rare opportunity to study the environments of compact SMBH binaries, and could enable the use of these sources as standard sirens for cosmology. Here we consider the feasibility of such an EM identification. We show that because the host galaxies of resolved PTA sources are expected to be exceptionally massive and rare, it should be possible to find unique hosts of resolved sources out to redshift z=0.2. At higher redshifts, the PTA error boxes are larger, and may contain as many as 100 massive-galaxy interlopers. The number of candidates, however, remains tractable for follow-up searches in upcoming wide-field EM surveys. We develop a toy model to characterize the dynamics and the thermal emission from a geometrically thin, gaseous disc accreting onto a PTA-source SMBH binary. Our model predicts that at optical and infrared frequencies, the source should appear similar to a typical luminous active galactic nucleus (AGN). However, owing to the evacuation of the accretion flow by the binary’s tidal torques, the source might have an unusually low soft X-ray luminosity and weak UV and broad optical emission lines, as compared to an AGN powered by a single SMBH with the same total mass. For sources near z=1, the decrement in the rest-frame UV should be observable as an extremely red optical color. These properties would make the PTA sources stand out among optically luminous AGN, and could allow their unique identification.

by Sesana, A. and Roedig, C. and Reynolds, M. T. and Dotti, M.
19 pages, 11 figures, submitted to MNRAS

We demonstrate that very massive (>10^8\msun), cosmologically nearby (z10^-13 erg s^-1 cm^-2 will be in the reach of upcoming X-ray observatories. Double relativistic K\alpha lines may be observable in a handful of low redshift (z<0.3) sources by proposed deep X-ray probes, such as Athena. (Abridged)

by Ponce, Marcelo and Faber, Joshua A. and Lombardi, James C.
Submitted to ApJ. Movies available as ‘ancillary files’

Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion disk. We study here the effect of an impulsive kick delivered to a black hole on the dynamical evolution of its accretion disk using a smoothed particle hydrodynamics code, focusing attention on the role played by the kick angle with respect to the orbital angular momentum vector of the pre-kicked disk. We find that for more vertical kicks, for which the angle between the kick and the normal vector to the disk $latex \theta\lesssim 30^\circ$, a gap remains present in the inner disk, in accordance with the prediction from an analytic collisionless Keplerian disk model, while for more oblique kicks with $latex \theta\gtrsim 45^\circ$, matter rapidly accretes toward the black hole. There is a systematic trend for higher potential luminosities for more oblique kick angles for a given black hole mass, disk mass and kick velocity, and we find large amplitude oscillations in time in the case of a kick oriented $latex 60^\circ$ from the vertical.

by Zhang, Xue-Guang
17 pages, 11 figures, MNRAS accepted

In this manuscript, the structure of broad emission line regions (BLRs) of well-mapping double-peaked emitter (AGN with broad double-peaked low-ionization emission lines) 3C390.3 is studied. Besides the best fitted results for double-peaked broad optical balmer lines of 3C390.3 by theoretical disk model, we try to find another way to further confirm the origination of double-peaked line from accretion disk. Based on the long-period observed spectra in optical band around 1995 collected from AGN WATCH project, the theoretical disk parameters of disk-like BLRs supposed by elliptical accretion disk model (Eracleous et al. 1995) have been well determined. Through the theoretical disk-like BLRs, characters of observed light-curves of broad double-peaked H$latex \alpha$ of 3C390.3 can be well reproduced based on the reverberation mapping technique. Thus the accretion disk model is preferred as one better model for BLRs of 3C390.3. Furthermore, we can find that different disk parameters should lead to some different results about size of BLRs of 3C390.3 from the one measured through observational data, which indicates the measured disk parameters are significantly valid for 3C390.3. After that, the precession of theoretical elliptical disk-like BLRs being considered, we can find that the expected line profile in 2000 by theoretical model is consistent with the observed line profile by HST around 2000. Based on the results, we can further believe that the origination of broad double-peaked balmer emission lines of 3C390.3 are from accretion disk around central black hole.

by Wong, Ka-Wah and Irwin, Jimmy A. and Yukita, Mihoko and Million, Evan T. and Mathews, William G. and Bregman, Joel N.
5 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are resolving, for the first time, the accretion flow within the Bondi radius of a SMBH. We show that the temperature is rising toward the galaxy center as expected in all accretion models in which the black hole is gravitationally capturing the ambient gas. There is no hard central point source that could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 4 arcsec-5 arcsec (188-235 pc), suggesting a SMBH of 2 x 10^9 M_sun that is consistent with the upper end of the optical results. The density profile within the Bondi radius has a power law index of 1.03^{+0.23}_{-0.21} which is consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate at the Bondi radius is determined to be {\dot M}_B = 2.2 x 10^{-2} M_sun yr^{-1}. Thus, the accretion luminosity with 10% radiative efficiency at the Bondi radius (10^{44} ergs s^{-1}) is about six orders of magnitude higher than the upper limit of the X-ray luminosity of the nucleus.

by Kiuchi, Kenta and Shibata, Masaru and Montero, Pedro J. and Font, José A.
4 pages, 4 figure, to be published in PRL

Black hole (BH)–torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $latex m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observable by forthcoming ground-based and spacecraft detectors.

by Farris, Brian D. and Liu, Yuk Tung and Shapiro, Stuart L.
submitted to PRD

Simultaneous gravitational and electromagnetic wave observations of merging black hole binaries (BHBHs) can provide unique opportunities to study gravitation physics, accretion and cosmology. Here we perform fully general relativistic, hydrodynamic simulations of equal-mass, nonspinning BHBHs coalescing in a circumbinary disk. We evolve the metric using the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of Einstein’s field equations with standard moving puncture gauge conditions. We handle the hydrodynamics via a high-resolution shock-capturing (HRSC) scheme. We track the inspiral starting from a binary separation of 10M, where M is the total binary mass. We take the disks to have an inner radius at R_in~15M to account for the hollow created by the binary torques. Our disks extend to R=65M and have an initial scale height of H/R=0.03-0.11. The gas is governed by a Gamma-law EOS, with Gamma equal to 5/3, 4/3, and 1.1. Disks are allowed to relax in the “early inspiral” epoch to provide quasistationary realistic initial data. We then evolve the metric and matter during the “late inspiral and merger” epoch. The later simulations are designed to track BHBH inspiral following disk-binary decoupling, through merger and ringdown, terminating before viscosity has time to fill the hollow about the remnant. We compute the gas flow and accretion rate and estimate the electromagnetic luminosity due to bremsstrahlung and synchrotron emission as a perturbation for optically thin disks. The synchrotron component of the luminosity peaks in the infrared band and should be detectable by WFIRST and possibly the LSST for a 10^8 M_sun binary embedded in a disk with a density n~10^12/cm^3 at z=1, beginning with a maximum value of $L~10^46 n_12^2 M_8^3 erg/s at decoupling, and decreasing steadily over a timescale of ~100 M_8 hours to a value of L~10^45 n_12^2 M_8^3 erg/s at merger.

by Noble, Scott C. and Krolik, Julian H. and Schnittman, Jeremy D. and Hawley, John F.
Submitted to ApJ, 26 pages, 12 figures (some in color), AASTEX

Recent general relativistic magneto-hydrodynamic (MHD) simulations of accretion onto black holes have shown that, contrary to the basic assumptions of the Novikov-Thorne model, there can be substantial magnetic stress throughout the plunging region. Additional dissipation and radiation can therefore be expected. We use data from a particularly well-resolved simulation of accretion onto a non-spinning black hole to compute both the radiative efficiency of such a flow and its spectrum if all emitted light is radiated with a thermal spectrum whose temperature matches the local effective temperature. This disk is geometrically thin enough (H/r ~= 0.06) that little heat is retained in the flow. In terms of light reaching infinity (i.e., after allowance for all relativistic effects and for photon capture by the black hole), we find that the radiative efficiency is at least ~=6-10% greater than predicted by the Novikov-Thorne model (complete radiation of all heat might yield another ~6%). We also find that the spectrum more closely resembles the Novikov-Thorne prediction for a/M ~= 0.2–0.3 than for the correct value, a/M=0. As a result, if the spin of a non-spinning black hole is inferred by model-fitting to a Novikov-Thorne model with known black hole mass, distance, and inclination, the inferred a/M is too large by ~= 0.2–0.3.

by Armijo, Matias Montesinos and Pacheco, José A. de Freitas
14 pages, 11 figures, Submitted to the ApJ

The evolution of an accretion disk, formed as a consequence of the disruption of a star by a black hole, is followed by solving numerically the hydrodynamic equations. The present investigation aims to study the dependence of resulting light curves on dynamical and physical properties of such a transient disk during its existence. One of main results derived from our simulations is that black body fits of X-ray data tend to overestimate the true mean disk temperature. The temperature derived from black body fits should be identified with the color X-ray temperature rather than the average value derived from the true temperature distribution along the disk. The time interval between the beginning of the circularization of the bound debris and the beginning of the accretion process by the black hole is determined by the viscous timescale, which fixes also the raising part of the resulting light curve. The luminosity peak coincides with the beginning of matter accretion by the black hole and the late evolution of the light curve depends on the evolution of the debris fallback rate. Peak bolometric luminosities are in the range 10^45-10^46 erg s^-1 whereas peak luminosities in soft X-rays (0.2-2.0 keV) are typically one order of magnitude lower. The timescale derived from our preferred models for the flare luminosity to decay by two orders of magnitude is about 3-4 years. Predicted soft X-ray light curves were fitted to data on galaxies in which a variable X-ray emission, related to tidal events, was detected.

by Abramowicz, Marek A. and Fragile, P. Chris
58 pages, 22 figures, submitted to Living Reviews in Relativity

In this review, we summarize the present status of knowledge of black hole accretion disks by discussing in some details the fundamental concepts on which the theory rests, and by describing the most often used analytic and semi-analytic models of them. We also describe numerical simulations of black hole accretion.

by Kocsis, Bence and Yunes, Nicolas and Loeb, Abraham
42 pages, 8 figures, 3 tables, submitted to Phys. Rev. D

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

by Guzman, F. S. and Lora-Clavijo, F. D.
9 pages, 24 eps figures, 2 tables. Accepted for publication in MNRAS

Based on the numerical solution of the time-dependent relativistic Euler equations onto a fixed Schwarzschild background space-time, we estimate the accretion rate of radial flow toward the horizon of a test perfect fluid obeying an ideal gas equation of state. We explore the accretion rate in terms of the initial density of the fluid for various values of the inflow velocity in order to investigate whether or not sufficiently arbitrary initial conditions allow a steady state accretion process depending on the values of the pressure. We extrapolate our results to the case where the fluid corresponds to dark matter and the black hole is a supermassive black hole seed. Then we estimate the equation of state parameters that provide a steady state accretion process. We found that when the pressure of the dark matter is zero, the black hole’s mass grows up to values that are orders of magnitude above $latex 10^{9}M_{\odot}$ during a lapse of 10Gyr, whereas in the case of the accretion of the ideal gas dark matter with non zero pressure the accreted mass can be of the order of $latex \sim 1M_{\odot}/10Gyr$ for black holes of $latex 10^{6}M_{\odot}$. This would imply that if dark matter near a supermassive black hole acquires an equation of state with non trivial pressure, the contribution of accreted dark matter to the supermassive black hole growth could be small, even though only radial accretion is considered.

by Zhuravlev, V. and Ivanov, P.
accepted for publication in MNRAS

(abbreviated) In this paper we derive equations describing dynamics and stationary configurations of a twisted fully relativistic thin accretion disc around a slowly rotating black hole. We find that the disc dynamics and stationary shapes are determined by a pair of equations for two complex variables describing orientation of the disc rings and velocity perturbations in the disc. We analyse shapes of stationary twisted configurations. It is shown that the stationary configurations depend on two parameters – the $latex \alpha $parameter and $latex \tilde \delta = \delta_{*}/\sqrt a$, where $latex \delta_{*}\sim h/r$ is the disc opening angle (h is the disc halfthickness) and $latex a$ is the black hole rotational parameter. When $latex a > 0$ and $latex \tilde \delta \ll 1$ the shapes depend drastically on value of $latex \alpha$. When $latex \alpha $ is small the disc inclination angle oscillates with radius with amplitude and radial frequency of the oscillations dramatically increasing towards the last stable orbit. For moderate values of $latex \alpha $ the oscillations do not take place but the disc does not align with the equatorial plane at small radii. Its inclination angle is either increasing towards $latex R_{ms}$ or exhibits a non-monotonic dependence on the radial coordinate. Finally, when $latex \alpha $ is sufficiently large the disc aligns with the equatorial plane at small radii. When $latex a < 0$ the disc aligns with the equatorial plane for all values of $latex \alpha $. The results reported here may have implications for determining structure and variability of accretion discs close to $latex R_{ms}$ as well as for modelling of emission spectra coming from different sources, which are supposed to contain black holes.

by Inogamov, N. A. and Sunyaev, R. A.
Full text of the paper with two Appendices is published in Astronomy Letters, Vol. 36, p. 835 (2010)

The accretion of hot slowly rotating gas onto a supermassive black hole is considered. Rotation velocities at the Bondi radius r_B are small in comparison with speed of sound c_s. The centrifugal barrier at a depth r_c = l^2/G M_BH <r_c) and inner (r<r_c) disks are formed. The gas enters the zone of the internal ADAF flow along the accretion disk (r<r_c). Since the angular momentum is conserved, the outer disk removes outward an excess of angular momentum along with part of the matter falling into the torus. Such outer Keplerian disk was observed by Hubble Space Telescope around the nucleus of the M87 galaxy in the optical emission lines. Turbulence causes rotation. We discuss the characteristic times during which the turbulence should lead to the changes in the orientation of the torus, accretion disk and, possibly, of the jet.

by Zamaninasab, M. and Eckart, A. and Dovciak, M. and Karas, V. and Schoedel, R. and Witzel, G. and Sabha, N. and Garcia-Marin, M. and Kunneriath, D. and Muzic, K. and Straubmeier, C. and Valencia-S, M. and Zensus, J. A.
12 pages, 13 figures, accepted for publication in MNRAS

Several massive black holes exhibit flux variability on time scales that correspond to source sizes of the order of few Schwarzschild radii. We survey the potential of near-infrared and X-ray polarimetry to constrain physical properties of such black hole systems, namely their spin and inclination. We have focused on a model where an orbiting hot spot is embedded in an accretion disk. A new method of searching for the time-lag between orthogonal polarization channels is developed and applied to an ensemble of hot spot models that samples a wide range of parameter space. We found that the hot spot model predicts signatures in polarized light which are in the range to be measured directly in the near future. However, our estimations are predicted upon the assumption of a Keplerian velocity distribution inside the flow where the dominant part of the magnetic field is toroidal. We also found that if the right model of the accretion flow can be chosen for each source (e.g. on the basis of magnetohydrodynamics simulations) then the black hole spin and inclination can be constrained to a small two-dimensional area in the spin-inclination space. The results of the application of the method to the available near-infrared polarimetric data of Sagittarius A* (Sgr A*) is presented. It is shown that even with the currently available data the spin and inclination of Sgr A* can be constrained. Next generations of near-infrared and X-ray polarimeters should be able to exploit this tool.

by Lodato, Giuseppe and Clarke, Cathie C.
7 pages, 2 figures. Accepted by MNRAS

Stimulated by recent results by Meru and Bate (2010a,b), we revisit the issue of resolution requirements for simulating self-gravitating accretion discs with Smoothed Particle Hydrodynamics (SPH). We show that the results by Meru and Bate (2010a) are consistent with those of Meru and Bate (2010b) if they are both interpreted as driven by resolution effects, therefore implying that the resolution criterion for cooling gaseous discs is a function of the imposed cooling rate. We discuss two possible numerical origins of such dependence, which are both consistent with the limited number of available data. Our results tentatively indicate that convergence for current simulations is being reached for a number of SPH particles approaching 10 millions (for a disc mass of order 10 per cent of the central object mass), which would set the critical cooling time for fragmentation at about $latex 15\Omega^{-1}$, roughly a factor two larger than previously thought. More in general, we discuss the extent to which the large number of recent numerical results are reliable or not. We argue that those results that pertain to the dynamics associated with gravitational instabilities (such as the locality of angular momentum transport, and the relationship between density perturbation and induced stress) are robust, while those pertaining to the thermodynamics of the system (such as the determination of the critical cooling time for fragmentation) can be affected by poor resolution.

by Pang, Bijia and Pen, Ue-Li and Matzner, Christopher D. and Green, Stephen R. and Liebendörfer, Matthias

We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars.

Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow.

We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.