Variability in GRMHD simulations of Sgr A∗^*: Implications for EHT closure phase observations

The observable quantities that carry the most information regarding the structures of the images of black holes in the interferometric observations with the Event Horizon Telescope are the closure phases along different baseline triangles. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A$^*$ to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of a small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the jet-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the $u-v$ space that show a high level of variability will need to be identified and excised.Comment: submitted to apj. 12 pages, 12 figure

Studies of Thermally Unstable Accretion Disks around Black Holes with Adaptive Pseudospectral Domain Decomposition Method. II. Limit-Cycle Behavior in accretion disks around Kerr black holes

For the first time ever, we derive equations governing the time-evolution of fully relativistic slim accretion disks in the Kerr metric, and numerically construct their detailed non-stationary models. We discuss applications of these general results to a possible limit-cycle behavior of thermally unstable disks. Our equations and numerical method are applicable in a wide class of possible viscosity prescriptions, but in this paper we use a diffusive form of the "standard alpha prescription" that assumes the viscous torque is proportional to the total pressure. In this particular case, we find that the parameters which dominate the limit-cycle properties are the mass-supply rate and the value of the alpha-viscosity parameter. Although the duration of the cycle (or the outburst) does not exhibit any clear dependence on the black hole spin, the maximal outburst luminosity (in the Eddington units) is positively correlated with the spin value. We suggest a simple method for a rough estimate of the black hole spin based on the maximal luminosity and the ratio of outburst to cycle durations. We also discuss a temperature-luminosity relation for the Kerr black hole accretion discs limit-cycle. Based on these results we discuss the limit-cycle behavior observed in microquasar GRS 1915+105. We also extend this study to several non-standard viscosity prescriptions, including a "delayed heating" prescription recently stimulated by the recent MHD simulations of accretion disks.Comment: 36 pages, 6 figures, 1 table. Accepted by ApJ

Slim disks around Kerr black holes revisited

We investigate stationary slim accretion disks around Kerr black holes. We construct a new numerical method based on the relaxation technique. We systematically cover the whole parameter space relevant to stellar mass X-ray binaries. We also notice some non-monotonic features in the disk structure, overlooked in previous studies.Comment: 12 pages, 11 figures. The relativistic slim accretion disk solutions have been published online at http://users.camk.edu.pl/as/slimdisk. An appendix containing the manual for the data available online has been added to the previous version of the pape

Vertical dissipation profiles and the photosphere location in thin and slim accretion disks

We calculate optically thick but geometrically thin (and slim) accretion disk models and perform a ray-tracing of photons (in the Kerr geometry) to calculate the observed disk spectra. Previously, it was a common practice to ray-trace photons assuming that they are emitted from the Kerr geometry equatorial plane, z = 0. We show that the spectra calculated with this assumption differ from these calculated under the assumption that photons are emitted from the actual surface of the disc, z = H(r). This implies that a knowledge of the location of the thin disks effective photosphere is relevant for calculating the spectra. In this paper we investigate, in terms of a simple toy model, a possible influence of the (unknown, and therefore ad hoc assumed) vertical dissipation profiles on the vertical structure of the disk and thus on the location of the effective photosphere, and on the observed spectra. For disks with moderate and high mass accretion rates (\dot m>0.01\dot m_C) we find that the photosphere location in the inner, radiation pressure dominated, disk region (where most of the radiation comes from) does not depend on the dissipation profile and therefore emerging disk spectra are insensitive to the choice of the dissipation function. For lower accretion rates the photosphere location depends on the assumed vertical dissipation profile down to the disk inner edge, but the dependence is very weak and thus of minor importance. We conclude that the spectra of optically thick accretion disks around black holes should be calculated with the ray-tracing from the effective photosphere and that, fortunately, the choice of a particular vertical dissipation profile does not substantially influence the calculated spectrum.Comment: 7 pages, 7 figure

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

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

Launching of jets by cold, magnetized disks in Kerr Metric

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

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

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

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

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