Our existence in the Universe resulted from a rare combination of circumstances. The same must hold for any highly developed extraterrestrial civilisation, and if they have ever existed in the Milky Way, they would likely be scattered over large distances in space and time. However, all technologically advanced species must be aware of the unique property of the galactic centre: it hosts Sagittarius A* (Sgr A*), the closest supermassive black hole to anyone in the Galaxy. A civilisation with sufficient technical know-how may have placed material in orbit around Sgr A* for research, energy extraction, and communication purposes. In either case, its orbital motion will necessarily be a source of gravitational waves. We show that a Jupiter-mass probe on the retrograde innermost stable circular orbit around Sgr A* emits, depending on the black hole spin, at a frequency of $f_{GW} = 0.63 - 1.07$ mHz and with a power of $P_{GW}=2.7 \times\, 10^{36} - 2.0 \times\, 10^{37}$ erg/s. We discuss that the energy output of a single star is sufficient to stabilise the location of an orbiting probe for a billion years against gravitational wave induced orbital decay. Placing and sustaining a device near Sgr A* is therefore astrophysically possible. Such a probe will emit an unambiguously artificial continuous gravitational wave signal that is observable with LISA-type detectors.

Context. The Event Horizon Telescope (EHT) collaboration recently obtained first images of the surroundings of the supermassive compact object M87* at the center of the galaxy M87. Aims. We want to develop a simple analytic disk model for the accretion flow of M87*. Compared to general-relativistic magnetohydrodynamic (GRMHD) models, it has the advantage of being independent of the turbulent character of the flow, and controlled by only few easy-to-interpret, physically meaningful parameters. We want to use this model to predict the image of M87* assuming that it is either a Kerr black hole, or an alternative compact object. Methods. We compute the synchrotron emission from the disk model and propagate the resulting light rays to the far-away observer by means of relativistic ray tracing. Such computations are performed assuming different spacetimes (Kerr, Minkowski, non-rotating ultracompact star, rotating boson star or Lamy spinning wormhole). We perform numerical fits of these models to the EHT data. Results. We discuss the highly-lensed features of Kerr images and show that they are intrinsically linked to the accretion-flow properties, and not only to gravitation. This fact is illustrated by the notion of secondary ring that we introduce. Our model of spinning Kerr black hole predicts mass and orientation consistent with the EHT interpretation. The non-Kerr images result in similar quality of the numerical fits and may appear very similar to Kerr images, once blurred to the EHT resolution. This implies that a strong test of the Kerr spacetime may be out of reach with the current data. We notice that future developments of the EHT could alter this situation. Conclusions. Our results show the importance of studying alternatives to the Kerr spacetime in order to be able to test the Kerr paradigm unambiguously.

We report on a new class of solutions of black hole accretion disks that we have found through three-dimensional, global, radiative magnetohydrodynamic simulations in general relativity. It combines features of the canonical thin, slim and thick disk models but differs in crucial respects from each of them. We expect these new solutions to provide a more realistic description of black hole disks than the slim disk model. We are presenting a disk solution for a non-spinning black hole at a sub-Eddington mass accretion rate, $\dot M=0.6\,\dot M_{\rm Edd}$. By the density scale-height measure the disk appears to be thin, having a high density core near the equatorial plane of height $h_{\rho} \sim 0.1 \,r$, but most of the inflow occurs through a highly advective, turbulent, optically thick, Keplerian region that sandwiches the core and has a substantial geometrical thickness comparable to the radius, $H \sim r$. The accreting fluid is supported above the midplane in large part by the magnetic field, with the gas and radiation to magnetic pressure ratio $\beta \sim 1$, this makes the disk thermally stable, even though the radiation pressure strongly dominates over gas pressure. A significant part of the radiation emerging from the disk is captured by the black hole, so the disk is less luminous than a thin disk would be at the same accretion rate.

Context. The properties of the accretion flow surrounding the supermassive central black hole of the Galaxy, Sgr A*, will be scrutinized by the new-generation instrument GRAVITY and the Event Horizon Telescope (EHT). Developing fast, robust, and simple models of such flows is thus important and very timely. Aims. We want to model the quiescent emission of Sgr A* from radio to mid-infrared by considering a magnetized compact torus and an extended jet. Results. We find perfect spectral fit both for face-on and edge-on views. These best fits give parameters values very close to that found by the most recent numerical simulations, which are much more complex than our model. The intrinsic radio size of Sgr A* is found to be in reasonable agreement with the centimetric observed constraints. Our best-fit infrared spectral index is in perfect agreement with the latest constraints. Our emission region at 1.3 mm, although larger than the Doeleman et al. (2008) Gaussian best-fit, does contain bright features at the 40 microarcsec scale. EHT-reconstructed images show that torus/jet-specific features persist after the reconstruction procedure, and that these features are sensitive to inclination. Conclusions. The main interest of our model is to give a simple and fast model of the quiescent state of Sgr A*, which gives extremely similar results as compared to state-of-the-art numerical simulations. Our model is easy to use and we publish all the material necessary to reproduce our spectra and images, so that anyone interested can use our results rather straightforwardly. We hope that such a public tool can be useful in the context of the recent and near-future GRAVITY and EHT results. Our model can in particular be easily used to test alternative compact objects models, or alternative gravity theories. The main limitation of our model is that we do not yet treat the X-ray emission.

A recent study of the X-ray source XMMUJ134736.6+173403 has revealed a strong quasi-periodic modulation of the X-ray flux. The observation of two quasiperiodic oscillations (QPOs) that occur on a daily timescale and exhibit a 3:1 frequency ratio strongly supports the evidence for the presence of an active galactic nucleus black hole (AGN BH). Assuming the orbital origin of QPOs we calculate the upper and lower limit on AGN BH mass M arriving at $M\approx 10^7-10^9M_{\odot}$. Comparing that to mass estimates of other sources, MMUJ134736.6+173403 appears to be the most massive source with commensurable QPO frequencies, and its mass represents the current observational upper limit on AGN BH mass based on QPO observations. We note it will be crucial for the falsification of particular resonance models of QPOs whether only single QPO with frequency that completes the harmonic sequence 3:2:1 is found in this source, or whether there is a new different pair of QPOs with frequencies in the 3:2 ratio. The former case would agree with the prediction of the 3:2 epicyclic resonance model and black hole mass $M\approx(5a^2 + 8a + 8)\times 10^{7}M_{\odot}$, where a is a dimensionless BH spin.

If primordial black holes with masses of $10^{25}\,\mbox{g}\gtrsim m \gtrsim 10^{17}\,\mbox{g}$ constitute a non-negligible fraction of the galactic dark-matter haloes, their existence should have observable consequences: they necessarily collide with galactic neutron stars, nest in their centers and accrete the dense matter, eventually converting them to neutron-star mass black holes while releasing the neutron-star magnetic field energy. Such processes may explain the fast radio bursts phenomenology, in particular their millisecond durations, large luminosities ${\sim}10^{43}$ erg/s, high rate of occurrence $\gtrsim 1000/\mbox{day}$, as well as high brightness temperatures, polarized emission and Faraday rotation. Longer than the dynamical timescale of the Bondi-like accretion for light primordial black holes allows for the repeating fast radio bursts. This explanation follows naturally from (assumed) existence of the dark matter primordial black holes and requires no additional unusual phenomena, in particular no unacceptably large magnetic fields of neutron stars. In our model, the observed rate of fast radio bursts throughout the Universe follows from the presently known number of neutron stars in the Galaxy.

The sum of squared epicyclic frequencies of nearly circular motion ($\omega_r^2+\omega_\theta^2$) in axially symmetric configurations of Newtonian gravity is known to depend both on the matter density and on the angular velocity profile of circular orbits. It was recently found that this sum goes to zero at the photon orbits of Schwarzschild and Kerr spacetimes. However, these are the only relativistic configurations for which such result exists in the literature. Here, we extend the above formalism in order to describe the analogous relation for geodesic motion in arbitrary static, axially symmetric, asymptotically flat solutions of general relativity. The sum of squared epicyclic frequencies is found to vanish at photon radii of vacuum solutions. In the presence of matter, we obtain that $\omega_r^2+\omega_\theta^2>0$ for perturbed timelike circular geodesics on the equatorial plane if the strong energy condition holds for the matter-energy fluid of spacetime; in vacuum, the allowed region for timelike circular geodesic motion is characterized by the inequality above. The results presented here may be of use to shed light on general issues concerning the stability of circular orbits once they approach photon radii, mainly the ones corresponding to stable photon motion.

We show that a recent assertion [arXiv:1608.01541] that gravitational wave emission can lead to a repulsive force explaining the accelerated expansion of the Universe is totally unfounded.

At super-Eddington rates accretion flows onto black holes have been described as slim (aspect ratio $H/R \lesssim 1$) or thick (H/R >1) discs, also known as tori or (Polish) doughnuts. The relation between the two descriptions has never been established, but it was commonly believed that at sufficiently high accretion rates slim discs inflate, becoming thick. We wish to establish under what conditions slim accretion flows become thick. We use analytical equations, numerical 1+1 schemes, and numerical radiative MHD codes to describe and compare various accretion flow models at very high accretion rates.We find that the dominant effect of advection at high accretion rates precludes slim discs becoming thick. At super-Eddington rates accretion flows around black holes can always be considered slim rather than thick.

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.

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.

Polish doughnuts (PDs) are geometrically thick disks that rotate with super-Keplerian velocities in their innermost parts, and whose long and narrow funnels along rotation axes collimate the emerging radiation into beams. In this paper we construct an extremal family of PDs that maximize both geometrical thickness and radiative efficiency. We then derive upper limits for these quantities and subsequently for the related ability to collimate radiation. PDs with such extreme properties may explain the observed properties of the ultraluminous X-ray sources without the need for the black hole masses to exceed ~ 10 solar masses. However, we show that strong advective cooling, which is expected to be one of the dominant cooling mechanisms in accretion flows with super-Eddington accretion rates, tends to reduce the geometrical thickness and luminosity of PDs substantially. We also show that the beamed radiation emerging from the PD funnels corresponds to "isotropic" luminosities that linearly scale with the mass accretion rate, and do not obey the familiar and well-known logarithmic relation.

In general relativity static gaseous atmospheres may be in hydrostatic balance in the absence of a supporting stellar surface, provided that the luminosity is close to the Eddington value. We construct analytic models of optically thin, spherically symmetric shells supported by the radiation pressure of a luminous central body in the Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. The inner parts of the atmospheres, where the luminosity locally has supercritical values, are characterized by a density and pressure inversion. The atmospheres are convectively and Rayleigh-Taylor stable, and there is no outflow of gas.

Possible astrophysical consequences of the Ho\v{r}ava quantum gravity theory have been recently studied by several authors. They usually employ the Kehagias-Sfetsos (KS) spacetime which is a spherically symmetric vacuum solution of a specific version of Ho\v{r}ava's gravity. The KS metric has several unusual geometrical properties that in the present article we examine by means of the often used technique of embedding diagrams. We pay particular attention to the transition between naked singularity and black-hole states, which is possible along some particular sequences of the KS metrics.

High frequency (millisecond) quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra of several microquasars and low mass X-ray binaries. Two distinct QPO peaks, so-called twin peak QPOs, are often detected simultaneously exhibiting their frequency ratio close or equal to 3/2. Following the analytic theory and previous studies of observable spectral signatures, we aim to model the twin peak QPOs as a spectral imprint of specific dual oscillation regime defined by a combination of the lowest radial and vertical oscillation mode of optically thick slender tori with constant specific angular momentum. We examined power spectra and fluorescent K$\alpha$ iron line profiles for two different simulation setups with the mode frequency relations corresponding to the epicyclic resonance HF QPOs model and modified relativistic precession QPOs model. We use relativistic ray-tracing implemented in parallel simulation code LSDplus. In the background of the Kerr spacetime geometry, we analyze the influence of the distant observer inclination and the spin of the central compact object. Relativistic optical projection of the oscillating slender torus is illustrated by images in false colours related to the frequency shift. We show that performed simulations yield power spectra with the pair of dominant peaks corresponding to the frequencies of radial and vertical oscillation modes with the proper ratio equal to 3/2 on a wide range of inclinations and spin values. We also discuss exceptional cases of a very small and very high inclination as well as unstable high spin relativistic precession-like configuration predicting constant frequency ratio equal to 1/2. We demonstrate signifiant dependency of broadened K$\alpha$ iron line profiles on the inclination of the distant observer.

Estimations of black hole spin in the three Galactic microquasars GRS 1915+105, GRO J1655-40, and XTE J1550-564 have been carried out based on spectral and timing X-ray measurements and various theoretical concepts. Among others, a non-linear resonance between axisymmetric epicyclic oscillation modes of an accretion disc around a Kerr black hole has been considered as a model for the observed high-frequency quasi-periodic oscillations (HF QPOs). Estimates of spin predicted by this model have been derived based on the geodesic approximation of the accreted fluid motion. Here we assume accretion flow described by the model of a pressure-supported torus and carry out related corrections to the mass-spin estimates. We find that for dimensionless black hole spin a<0.9, the resonant eigenfrequencies are very close to those calculated for the geodesic motion. Their values slightly grow with increasing torus thickness. These findings agree well with results of a previous study carried out in the pseudo-Newtonian approximation. The situation becomes different for a>0.9, in which case the resonant eigenfrequencies rapidly decrease as the torus thickness increases. We conclude that the assumed non-geodesic effects shift the lower limit of the spin, implied for the three microquasars by the epicyclic model and independently measured masses, from a~0.7 to a~0.6. Their consideration furthermore confirms compatibility of the model with the rapid spin of GRS 1915+105 and provides highly testable predictions of the QPO frequencies. Individual sources with a moderate spin (a<0.9) should exhibit a smaller spread of the measured 3:2 QPO frequencies than sources with a near-extreme spin (a~1). This should be further examined using the large amount of high-resolution data expected to become available with the next generation of X-ray instruments, such as the proposed Large Observatory for X-ray Timing (LOFT).

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

High frequency quasi-periodic oscillations (HF QPOs) are observed in the X-ray power-density spectra (PDS) of several microquasars and low mass X-ray binaries. Many proposed QPO models are based on oscillations of accretion toroidal fluid structures orbiting in the vicinity of a compact object. We study oscillating accretion tori orbiting in the vicinity of a Kerr black hole. We demonstrate that significant variation of the observed flux can be caused by the combination of radial and vertical oscillation modes of a slender, polytropic, perfect fluid, non-self-graviting torus with constant specific angular momentum. We investigate two combinations of the oscillating modes corresponding to the direct resonance QPO model and the modified relativistic precession QPO model.

Modified Newtonian potentials have been proposed for the description of relativistic effects acting on particles and fluids in permanent orbital motion around black holes. Here we further discuss spherically symmetric potentials like the one proposed by Artemova, Bjornson & Novikov (1996, Astrophysical Journal, 461, 565), and we illustrate their virtues by studying the acceleration along circular trajectories. We compare the results with exact expressions in the spacetime of a rotating (Kerr) black hole.

Context. The supermassive black hole, Sagittarius (Sgr) A*, in the centre of our Galaxy has the largest angular size in the sky among all astrophysical black holes. Its shadow, assuming no rotation, spans ~ 50 microarcsec. Resolving such dimensions has long been out of reach for astronomical instruments until a new generation of interferometers being operational during this decade. Of particular interest is the Event Horizon Telescope (EHT) with resolution ~ 20 microarcsec in the millimeter-wavelength range 0.87 mm - 1.3 mm. Aims. We investigate the ability of the fully general relativistic Komissarov (2006) analytical magnetized torus model to account for observable constraints at Sgr A* in the centimeter and millimeter domains. The impact of the magnetic field geometry on the observables is also studied. Methods. We calculate ray-traced centimeter- and millimeter-wavelength synchrotron spectra and images of a magnetized accretion torus surrounding the central black hole in Sgr A*. We assume stationarity, axial symmetry, constant specific angular momentum and polytropic equation of state. A hybrid population of thermal and non-thermal electrons is considered. Results. We show that the torus model is capable of reproducing spectral constraints in the millimeter domain, and in particular in the observable domain of the EHT. However, the torus model is not yet able to fit the centimeter spectrum. 1.3 mm images at high inclinations are in agreement with observable constraints. Conclusions. The ability of the torus model to account for observations of Sgr A* in the millimeter domain is interesting in the perspective of the future EHT. Such an analytical model allows very fast computations. It will thus be a suitable test bed for investigating large domains of physical parameters, as well as non-black-hole compact object candidates and alternative theories of gravity.

In the context of high-frequency quasi-periodic oscillation (HF QPOs) we further explore the appearance of an observable signal generated by hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. The aim of our investigation is to reveal whether observable characteristics of the Fourier power-spectral density can help us to distinguish between the two competing models, namely, the idea of bright spots orbiting on the surface of an accretion torus versus the scenario of intrinsic oscillations of the torus itself. We take the capabilities of the present observatories (represented by the Rossi X-ray Timing Explorer, RXTE) into account, and we also consider the proposed future instruments (represented here by the Large Observatory for X-ray Timing, LOFT).

In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{\mu \nu}$ and show that the necessary and sufficient condition for extraction of a black hole's rotational energy is analogous to that in the mechanical Penrose process: absorption of negative energy and negative angular momentum. We also show that a necessary condition for the Penrose process to occur is for the Noether current (the conserved energy-momentum density vector) to be spacelike or past directed (timelike or null) on some part of the horizon. In the particle case, our general criterion for the occurrence of a Penrose process reproduces the standard result. In the case of relativistic jet-producing "magnetically arrested disks" we show that the negative energy and angular-momentum absorption condition is obeyed when the Blandford-Znajek mechanism is at work, and hence the high energy extraction efficiency up to $\sim 300\%$ found in recent numerical simulations of such accretion flows results from tapping the black hole's rotational energy through the Penrose process. We show how black-hole rotational energy extraction works in this case by describing the Penrose process in terms of the Noether current.

Context. Some microquasars exhibit millisecond quasi-periodic oscillations (QPO) that are likely related to phenomena occuring in the immediate vicinity of the central black hole. Oscillations of accretion tori have been proposed to model these QPOs. Aims. Here, we aim at determining the observable spectral signature of slender accretion tori surrounding Kerr black holes. We analyze the impact of the inclination and spin parameters on the power spectra. Methods. Ray-traced power spectra of slender tori oscillation modes are computed in the Kerr metric. Results. We show that the power spectral densities of oscillating tori are very sensitive to the inclination and spin parameters. This strong dependency of the temporal spectra on inclination and spin may lead to observable constraints of these parameters. Conclusions. This work goes a step further in the analysis of the oscillating torus QPO model. It is part of a long-term study that will ultimately lead to comparison with observed data.

High frequency quasi-periodic oscillations (HF QPOs) appear in the X-ray variability of several accreting low-mass binaries. In a series of works it was suggested that these QPOs may have connection to inhomogeneities orbiting close to an inner edge of the accretion disc. In this paper we explore the appearance of an observable signal generated by small radiating circular hot spots moving along quasi-elliptic trajectories close to the innermost stable circular orbit in the Schwarzschild spacetime. Our consideration takes into account the capabilities of observatories that have been operating in the past two decades represented by the Rossi X-ray Timing Explorer (RXTE) and the proposed future instruments represented by the Large Observatory for X-ray Timing (LOFT). For these purposes we choose such model parameters that lead to lightcurves comparable to those observed in Galactic black hole sources, in particular the microquasar GRS 1915+105. We find that when a weak signal corresponding to the hot-spot Keplerian frequency is around the limits of the RXTE detectability, the LOFT observations can clearly reveal its first and second harmonics. Moreover, in some specific situations the radial epicyclic frequency of the spot can be detected as well. Finally, we also compare the signal produced by the spots to the signal produced by axisymmetric epicyclic disc-oscillation modes and discuss the key differences that could be identified via the proposed future technology. We conclude that the ability to recognize the harmonic content of the signal can help to distinguish between the different proposed physical models.

We discuss photon and test-particle orbits in the Kehagias-Sfetsos (KS) metric. For any value of the Ho\v{r}ava parameter $\omega$, there are values of the gravitational mass $M$ for which the metric describes a naked singularity, and this is always accompanied by a vacuum "antigravity sphere" on whose surface a test particle can remain at rest (in a zero angular momentum geodesic), and inside which no circular geodesics exist. The observational appearance of an accreting KS naked singularity in a binary system would be that of a quasi-static spherical fluid shell surrounded by an accretion disk, whose properties depend on the value of $M$, but are always very different from accretion disks familiar from the Kerr-metric solutions. The properties of the corresponding circular orbits are qualitatively similar to those of the Reissner-Nordstr\"om naked singularities. When event horizons are present, the orbits outside the Kehagias-Sfetsos black hole are qualitatively similar to those of the Schwarzschild metric.

Context: An impulsive increase in luminosity by one half or more of the Eddington value will lead to ejection of all optically thin plasma from Keplerian orbits around the radiating star, if gravity is Newtonian and the Poynting-Robertson drag is neglected. Radiation drag may bring some particles down to the stellar surface. On the other hand, general relativistic calculations show that gravity may be balanced by a sufficiently intense radiation field at a certain distance from the star. Aims: We investigate the motion of test particles around highly luminous stars to determine conditions under which plasma may be ejected from the system. Results: In Einstein's gravity, if the outburst is close to the Eddington luminosity, all test particles orbiting outside an "escape sphere" will be ejected from the system, while all others will be captured from their orbits onto the surface of another sphere, which is well above the stellar surface, and may even be outside the escape sphere, depending on the value of luminosity. Radiation drag will bring all the captured particles to rest on this "Eddington capture sphere," where they will remain suspended in an equilibrium state as long as the local flux of radiation does not change and remains at the effective Eddington value.

Context. One of the often discussed models for X-ray binaries high-frequency quasi-periodic oscillations is the oscillating torus model that considers oscillation modes of slender accretion tori. Aims. Here, we aim at developing this model by considering the observable signature of an optically thick slender accretion torus subject to simple periodic deformations. Methods. We compute light curves and power spectra of a slender accretion torus subject to simple periodic deformations: vertical or radial translation, rotation, expansion and shear. Results. We show that different types of deformations lead to very different Fourier power spectra and therefore could be observationally distinguished. Conclusions. This work is a first step in a longer-term study of the observable characteristics of the oscillating torus model. It gives promising perspectives on the possibility to constrain this model by studying the observed power spectra of quasi-periodic oscillations.

We show that results of a simple dynamical gedanken experiment interpreted according to standard Newton's gravitational theory, may reveal that three-dimensional space is curved. The experiment may be used to reconstruct the curved geometry of space, i.e. its non-Euclidean metric. The perihelion of Mercury advance and the light bending calculated from the Poisson equation and the equation of motion in the curved geometry have the correct (observed) values. Independently, we also show that Newtonian gravity theory may be enhanced to incorporate the curvature of three dimensional space by adding an extra equation which links the Ricci scalar with the density of matter. Like in Einstein's general relativity, matter is the source of curvature. In the spherically symmetric (vacuum) case, the metric of space 3gik that follows from this extra equation agrees, to the expected accuracy, with the metric measured by the Newtonian gedanken experiment mentioned above.

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).

Three radii are associated with a circle: the "geodesic radius" R_1 which is the distance from circle's center to its perimeter, the "circumferential radius" R_2 which is the length of the perimeter divided by 2 pi and the "curvature radius" R_3 which is circle's curvature radius in the Frenet sense. In the flat Euclidean geometry it is R_1 = R_2 = R_3, but in a curved space these three radii are different. I show that although Newton's dynamics uses Euclidean geometry, its equations that describe circular motion in spherical gravity always unambiguously refer to one particular radius of the three --- geodesic, circumferential, or curvature. For example, the gravitational force is given by F = -GMm/(R_2)^2, and the centrifugal force by mv^2/R_3. Building on this, I derive a Newtonian formula for the perihelion of Mercury advance.

A dormant Swift source J 164449.3+573451 (Sw 164449+57)recently experienced a powerful outburst, caused most probably by a tidal disruption of a star by the super-massive black hole at the center of the source. During the outburst, a quasi periodic oscillation (QPO) was detected in the observed X-ray flux from Sw 164449+57. We show that if the observed QPO belongs to a "3:2 twin peak QPO" (with the second frequency not observed), the mass of the black hole in Sw 164449+57 is rather low, M ~ 10^5 M_sun, and the source belongs to a class of intermediate mass black holes. The low mass of the source has been pointed out previously by several authors.

The Galactic center is an ideal laboratory to study strong-field general relativistic phenomena, as the supermassive black hole Sgr A* has the biggest angular Schwarzschild radius among all black holes. This article presents three different ways of using the immediate surroundings of Sgr A* as a probe in order to either constrain its spin, or even test the very nature of this compact object.

We present a toy model of mildly super-Eddington, optically thin accretion onto a compact star in the Schwarzschild metric, which predicts periodic variations of luminosity when matter is supplied to the system at a constant accretion rate. These are related to the periodic appearance and disappearance of the Eddington Capture Sphere. In the model the frequency is found to vary inversely with the luminosity. If the input accretion rate varies (strictly) periodically, the luminosity variation is quasi-periodic, and the quality factor is inversely proportional to the relative amplitude of mass accretion fluctuations, with its largest value approximately Q= 1/(10 |delta Mdot/Mdot|) attained in oscillations at about 1 to 2 kHz frequencies for a 2 solar mass star.

Collisions of particles in black holes' ergospheres may result in an arbitrarily large center of mass energy. This led recently to the suggestion (Banados et al., 2009) that black holes can act as ultimate particle accelerators. If the energy of an outgoing particle is larger than the total energy of the infalling particles the energy excess must come from the rotational energy of the black hole and hence this must involve a Penrose process. However, while the center of mass energy diverges the position of the collision makes it impossible for energetic particles to escape to infinity. Following an earlier work on collisional Penrose processes (Piran & Shaham 1977) we show that even under the most favorable idealized conditions the maximal energy of an escaping particle is only a modest factor above the total initial energy of the colliding particles. This implies that one shouldn't expect collisions around a black hole to act as spectacular cosmic accelerators.

Test particles infalling from infinity onto a compact spherical star with a mildly super-Eddington luminosity at its surface are typically trapped on the "Eddington Capture Sphere" and do not reach the surface of the star. The presence of a sphere on which radiation pressure balances gravity for static particles was first discovered some twenty five years ago. Subsequently, it was shown to be a capture sphere for particles in radial motion, and more recently also for particles in non-radial motion, in which the Poynting-Robertson radiation drag efficiently removes the orbital angular momentum of the particles, reducing it to zero. Here we develop this idea further, showing that "levitation" on the Eddington sphere (above the stellar surface) is a state of stable equilibrium, and discuss its implications for Hoyle-Lyttleton accretion onto a luminous star. When the Eddington sphere is present, the cross-section of a compact star for actual accretion is typically less than the geometrical cross-section (pi Rsquared), direct infall onto the stellar surface only being possible for relativistic particles, with the required minimum particle velocity at infinity typically ~1/2 the speed of light. We further show that particles on typical trajectories in the vicinity of the stellar surface will also be trapped on the Eddington Capture Sphere.

We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.

The two high frequency quasi periodic oscillations (HFQPOs) recently reported in the black hole candidate IGR J17091-3624 by Altamirano and Belloni (2012) are in a 5:2 frequency ratio (164 Hz to 66 Hz). This ratio is strongly suggestive of period doubling and nonlinear resonance analogous to phenomena known in RV Tauri-type pulsating stars (and recently discovered also in oscillations of RR Lyrae-type and of BL Herculis-type variables). An interpretation of the frequency ratio in terms of nonlinear interactions and a comparison with the HFQPOs reported in GRS 1915+105 may imply a mass of about 6 solar masses for the black hole in IGR J17091- 3624.

We explore effects of the Shakura-Sunyaev alpha-viscosity on the dynamics and oscillations of slender tori. We start with a slow secular evolution of the torus. We show that the angular-momentum profile approaches the Keplerian one on the timescale longer than a dynamical one by a factor of the order of 1/\alpha. Then we focus our attention on the oscillations of the torus. We discuss effects of various angular momentum distributions. Using a perturbation theory, we have found a rather general result that the high-order acoustic modes are damped by the viscosity, while the high-order inertial modes are enhanced. We calculate a viscous growth rates for the lowest-order modes and show that already lowest-order inertial mode is unstable for less steep angular momentum profiles or very close to the central gravitating object.

The frequencies of kHz QPOs are widely interpreted as being indicative of the values of characteristic frequencies related to orbital motion around neutron stars, e.g., the radial epicyclic frequency. In regions directly exposed to the radiation from the luminous neutron star these frequencies change with the luminosity. Including radiative corrections will change the neutron star mass value inferred from the QPO frequencies. Radiative forces may also be behind the puzzling phenomenon of parallel tracks.

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.

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.

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.

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

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

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.

It was suggested by several authors that hypothetical primordial black holes (PBHs) may contribute to the dark matter in our Galaxy. There are strong constraints based on the Hawking evaporation that practically exclude PBHs with masses m~1e15-1e16g and smaller as significant contributors to the Galactic dark matter. Similarly, PBHs with masses greater than about 1e26g are practically excluded by the gravitational lensing observation. The mass range between 10e16g<m<10e26g is unconstrained. In this paper, we examine possible observational signatures in the unexplored mass range, investigating hypothetical collisions of PBHs with main sequence stars, red giants, white dwarfs, and neutron stars in our Galaxy. This has previously been discussed as possibly leading to an observable photon eruption due to shock production during the encounter. We find that such collisions are either too rare to be observed (if the PBH masses are typically larger than about 1e20g), or produce too little power to be detected (if the masses are smaller than about 1e20g).

In the mass range of 1e15 g up to 1e26 g, primordial black holes (PBHs) as a possible contribution to the dark matter are still unexplored. In this contribution, we investigate the possibility of an electromagnetic signal from PBH interactions with astrophysical objects in the Galaxy. We find that a signal from passages cannot be observed, since, depending on the mass, either the interaction probability or the energy loss is too small. Further, we discuss possible effects from high-mass PBHs at masses >1e26 g, where PBHs can still contribute to the dark matter at the order of ~10%. Here, we find that a significant fraction of PBHs can be captured in the Hubble time. These captures could therefore lead to detectable effects.

Paczynski realized that a properly chosen gravitational potential may accurately model (in a "pseudo Newtonian" theory) general relativistic effects that determine motion of matter near a non-rotating black hole. Paczynski's choice, known today as the "Paczynski-Wiita potential", proved to be very practical. It was used by numerous researchers in the black hole accretion theory, and became a standard tool in relativistic astrophysics. The model is an example of Paczynski's admired ability to invent "out of nowhere" simple ideas that were brilliant, deep and useful. Paczynski has guessed intuitively the form of the potential. However, it could be also derived by a a step-by-step formal procedure. I show the derivation here. My derivation is based on a standard definition of the relativistic "effective potential" in the Schwarzschild spacetime. The relativistic effective potential may be uniquely divided into its "gravitational" and "centrifugal" part. The gravitational part differs from the Paczynski-Wiita potential only by a constant.

We construct a new family of analytic models of black hole accretion disks in dynamical equilibria. Our construction is based on assuming distributions of angular momentum and entropy. For a particular choice of the distribution of angular momentum, we calculate the shapes of equipressure surfaces. The equipressure surfaces we find are similar to those in thick, slim and thin disks, and to those in ADAFs.

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.

We explore some properties of twin kilohertz quasiperiodic oscillations (QPOs) in a simple toy-model consisting of two oscillation modes coupled by a general nonlinear force. We examine resonant effects by slowly varying the values of the tunable, and nearly commensurable, eigenfrequencies. The behavior of the actual oscillation frequencies and amplitudes during a slow transition through the 3:2 resonance is examined in detail and it is shown that both are significantly affected by the nonlinearities in the governing equations. In particular, the amplitudes of oscillations reflect a resonant exchange of energy between the modes, as a result the initially weaker mode may become dominant after the transition. We note that a qualitatively similar behavior has been recently reported in several neutron star sources by Torok (2008, arXiv:0812.4751), who found that the difference of amplitudes in neutron star twin peak QPOs changes sign as the observed frequency ratio of the QPOs passes through the value 3:2.

I review five of Bohdan Paczynski's ideas on black hole accretion disk theory. They formed my understanding of the subject and often guided intuition in my research. They are fundamentally profound, rich in physical consequences, mathematically elegant and clever, and in addition are useful in several technically difficult practical applications.

Recently [arXiv:astro-ph/0612155] we presented a formal mathematical proof that, contrary to a widespread misconception, cosmological expansion cannot be understood as the motion of galaxies in non-expanding space. We showed that the cosmological redshift must be physically interpreted as the expansion of space. Although our proof was generally accepted, a few authors disagreed. We rebut their criticism in this Note.

A long discussion has been devoted to the issue of clustering of the kHz QPO frequency ratios in the neutron star sources. While the distribution of ratios inferred from an occurrence of a single QPO seems to be consistent with a random walk, the distribution based on simultaneous detections of both peaks indicates a preference of ratios of small integers. Based on the public RXTE data we further investigate this issue for the source 4U 1636-53. Quality factors and rms amplitudes of both the QPOs nearly equal to the points where the frequencies are commensurable, and where the twin QPO detections cluster. We discuss a connection of the clustering with the varying properties of the two QPO modes. Assuming approximative relations for the observed correlations of the QPO properties, we attempt to reproduce the frequency and ratio distributions using a simple model of a random-walk evolution along the observed frequency-frequency correlation. We obtain results which are in qualitative agreement with the observed distributions.

A recently published study on long term evolution of the frequencies of the kilohertz quasi-periodic oscillations (QPOs) in the atoll source 4U 1636-53 concluded that there is no preferred frequency ratio in a distribution of twin QPOs that was inferred from the distribution of a single frequency alone. However, we find that the distribution of the ratio of actually observed pairs of kHz QPO frequencies is peaked close to the 3/2 value, and possibly also close to the 5/4 ratio. To resolve the apparent contradiction between the two studies, we examine in detail the frequency distributions of the lower kHz QPO and the upper kHz QPO detected in our data set. We demonstrate that for each of the two kHz QPOs (the lower or the upper), the frequency distribution in all detections of a QPO differs from the distribution of frequency of the same QPO in the subset of observations where both the kHz QPOs are detected. We conclude that detections of individual QPOs alone should not be used for calculation of the distribution of the frequency ratios.

Context. The high-frequency quasi-periodic oscillations (HF QPOs) in neutron star and stellar-mass black hole X-ray binaries may be the result of a resonance between the radial and vertical epicyclic oscillations in strong gravity. Aims. In this paper we investigate the resonant coupling between the epicyclic modes in a torus in a strong gravitational field. Methods. We perform numerical simulations of axisymmetric constant angular momentum tori in the pseudo-Newtonian potential. The epicyclic motion is excited by adding a constant radial velocity to the torus. Results. We verify that slender tori perform epicyclic motions at the frequencies of free particles, but the epicyclic frequencies decrease as the tori grow thicker. More importantly, and in contrast to previous numerical studies, we do not find a coupling between the radial and vertical epicyclic motions. The appearance of other modes than the radial epicyclic motion in our simulations is rather due to small numerical deviations from exact equilibrium in the initial state of our torus. Conclusions. We find that there is no pressure coupling between the two axisymmetric epicyclic modes as long as the torus is symmetric with respect to the equatorial plane. However we also find that there are other modes in the disc that may be more attractive for explaining the HF QPOs.

We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which is important for the theory of high-frequency quasi-periodic oscillations (QPOs), occurs when the fluid tori epicyclic frequencies nu_r and nu_z are about 15% higher than the ones corresponding to free particles. Our results therefore suggest that previous estimates of black hole spins from QPOs have produced values that are too high.

The rather wide-spread belief that cosmological expansion of a flat 3D--space (with spatial curvature k=0) cannot be observationally distinguished from a kinematics of galaxies moving in a flat and non-expanding space is erroneous. We suggest that the error may have its source in a non relativistic intuition that imagines the Universe not as a spacetime but separates space from time and pictures the cosmological expansion as space evolving in time. The physical reality, however, is fundamentally different -- the expanding Universe is necessarily a curved spacetime. We show here that the fact that the spacetime is curved implies that the interpretation of the observed cosmological redshift as being due to the expansion of the cosmological 3D--space is observationally verifiable. Thus it is impossible to mimic the true cosmological redshift by a Doppler effect caused by motion of galaxies in a non-expanding 3D-space, flat or curved. We summarize our points in simple spacetime diagrams that illustrate a gedanken experiment distinguishing between expansion of space and pure kinematics. We also provide all relevant mathematics. None of the previously published discussions of the issue, including a recent popular Scientific American article, offered a similarly clear way out of the confusion.

We consider a new version of the twin paradox. The twins move along the same circular free photon path around the Schwarzschild center. In this case, despite their different velocities, all twins have the same non-zero acceleration. On the circular photon path, the symmetry between the twins situations is broken not by acceleration (as it is in the case of the classic twin paradox), but by the existence of an absolute standard of rest (timelike Killing vector). The twin with the higher velocity with respect to the standard of rest is younger on reunion. This closely resembles the case of periodic motions in compact (non-trivial topology) 3-D space recently considered in the context of the twin paradox by Barrow and Levin, except that there accelerations of all twins were equal to zero, and that in the case considered here, the 3-D space has trivial topology.

In the Galactic microquasars with double peak kHz quasi-periodic oscillations (QPOs), the ratio of the two frequencies is 3:2. This supports the suggestion that double peak kHz QPOs are due to a non-linear resonance between two modes of accretion disk oscillations. For the microquasars with known mass, we briefly compare the black hole spin estimates based on the orbital resonance model with the recently reported spin predictions obtained by fitting the spectral continua. Results of these two approaches are not in good agreement. We stress that if the spectral fit estimates are accurate and can be taken as referential (which is still questionable), the disagreement between the predicted and referential values would represent a rather generic problem for any relativistic QPO model, as no spin influence would appear in the observed 1/M scaling of the QPO frequencies. The epicyclic frequencies relevant in these models are often considered to be equal to those of a test particle motion. However modifications of the frequencies due to the disc pressure or other non-geodesic effects may play an important role, and the inaccuracy introduced in the spin estimates by the test particle approximation could be crucial.

Following our previous work of Lachowicz & Czerny (2005), we explore further the application of the continuous wavelet transform to X-ray astronomical signals. Using the public archive of the XMM-Newton satellite, we analyze all available EPIC-pn observations for nearby Seyfert 1 galaxies MCG-6-30-15 and NGC 4051. We confine our analysis to 0.002-0.007 Hz frequency band in which, on the way of theoretically motivated premises, some quasi-periodic oscillations (QPOs) are expected to be found. We find that indeed wavelet power histogram analysis reveals such QPOs centered at two frequencies of ~2.5E-3 Hz and 4-6E-3 Hz, respectively. We show that these quasi-periodic features can be disentangled from the Poisson noise contamination level what is hardly to achieve with the standard Fourier analysis. Interestingly, we find some of them to be in 2:1 or 3:2 ratio. If real, our finding may be considered as a link between QPOs observed in AGN and kHz QPOs seen in X-ray binary systems.

Kilohertz Quasi-Periodic Oscillations (QPOs) have been detected in many accreting X-ray binaries. It has been suggested that the highest QPO frequencies observed in the modulation of the X-ray flux reflect a non-linear resonance between two modes of accreting disk oscillation. This hypothesis implies certain very general predictions, several of which have been borne out by observations. Some of these follow from properties of non-linear oscillators, while the others are specific to oscillations of fluid in strong gravity. A 3:2 resonant ratio of frequencies can be clearly recognized in the black-hole as well as in the neutron-star QPO data.

There is strong evidence that the observed kHz Quasi Periodic Oscillations (QPOs) in the X-ray flux of neutron star and black hole sources in LMXRBs are linked to Einstein's General Relativity. Abramowicz&Klu\'zniak (2001) suggested a non-linear resonance model to explain the QPOs origin: here we summarize their idea and the development of a mathematical toy-model which begins to throw light on the nature of Einstein's gravity non-linear oscillations.

Quasi periodic oscillations (QPO) seen in the X-ray fluxes of individual neutron stars and black hole sources are one of most intriguing phenomena in today's astrophysics. The QPO nature is visibly determined by super-strong Einstein's gravity. I argue here that it also profoundly depends on the MRI turbulence in accretion flows. Understanding the QPO physics may therefore guide accretion theory out of its present state of confusion.

Fluids in external gravity may oscillate with frequencies characteristic of the epicyclic motions of test particles. We explicitly demonstrate that global oscillations of a slender, perfect fluid torus around a Kerr black hole admit incompressible vertical and radial epicyclic modes. Our results may be directly relevant to one of the most puzzling astrophysical phenomena -- high (hundreds of hertz) frequency quasiperiodic oscillations (QPOs) detected in X-ray fluxes from several black hole sources. Such QPOs are pairs of stable frequencies in the 3/2 ratio. It seems that they originate a few gravitational radii away from the black hole and thus observations of them have the potential to become an accurate probe of super-strong gravity.

Double peak kHz QPO frequencies in neutron star sources varies in time by a factor of hundreds Hz while in microquasar sources the frequencies are fixed and located at the line \nu_2 = 1.5 \nu_1 in the frequency-frequency plot. The crucial question in the theory of twin HFQPOs is whether or not those observed in neutron-star systems are essentially different from those observed in black holes. In black hole systems the twin HFQPOs are known to be in a 3:2 ratio for each source. At first sight, this seems not to be the case for neutron stars. For each individual neutron star, the upper and lower kHz QPO frequencies, \nu_2 and \nu_1, are linearly correlated, \nu_2=A \nu_1 + B, with the slope A < 1.5, i.e., the frequencies definitely are not in a 1.5 ratio. In this contribution we show that when considered jointly on a frequency-frequency plot, the data for the twin kHz QPO frequencies in several (as opposed to one) neutron stars uniquely pick out a certain preferred frequency ratio that is equal to 1.5 for the six sources examined so far.

Recent observations, reported by Warner and Woudt, of Dwarf Nova Oscillations (DNOs) exhibiting frequency drift, period doubling, and 1:2:3 harmonic structure, can be understood as disc oscillations that are excited by perturbations at the spin frequency of the white dwarf or of its equatorial layers. Similar quasi-periodic disc oscillations in black hole low-mass X-ray binary (LMXB) transients in a 2:3 frequency ratio show no evidence of frequency drift and correspond to two separate modes of disc oscillation excited by an internal resonance. Just as no effects of general relativity play a role in white dwarf DNOs, no stellar surface or magnetic field effects need be invoked to explain the black hole QPOs.

We present and discuss a short and simple derivation of orbital epicyclic frequencies for circular geodesic orbits in stationary and axially symmetric spacetimes. Such spacetimes include as special cases analytically known black hole Kerr and Schwarzschild spacetimes, as well as the analytic Hartle-Thorne spacetime and all numerically constructed spacetimes relevant for rotating neutron stars. Our derivation follows directly from energy and angular momentum conservation and it uses the concept of the effective potential. It has never been published, except for a few special cases, but it has already become a part of the common knowledge in the field.

The familiar Newtonian version of the Rayleigh criterion demands that for dynamical stability the specific angular momentum should increase with the increasing circumferential radius of circular trajectories of matter. However, sufficiently close to a black hole the Rayleigh criterion reverses: in stable fluids, the angular momentum must decrease with the increasing circumferential radius. The geometrical reason for this reversal is that the space is so very strongly warped that it turns inside-out.

The theory of highly super-Eddington black hole accretion was developed in the 1980s in Warsaw by Paczynski and his collaborators in terms of "Polish doughnuts", i.e. low viscosity, rotating accretion flows that are optically thick, radiation pressure supported, cooled by advection, and radiatively very inefficient. Polish doughnuts resemble fat tori with two narrow funnels along rotation axis. The funnels collimate radiation into beams with highly super-Eddington luminosities. "Slim disks" introduced later by Abramowicz, Lasota and collaborators, have only moderately super-Eddington accretion rates, rather disk-like shapes, and almost thermal spectra.

The origin of quasi periodic modulations of flux in the kilohertz range (kHz QPOs), observed in low-mass X-ray binaries, is usually assumed to be physically distinct from that of the ``normal branch oscillations'' (NBOs) in the Z-sources. We show that a low-frequency modulation of the kHz QPOs is a natural consequence of the non-linear relativistic resonance suggested previously to explain the properties of the high-frequency twin peaks. The theoretical results discussed here are reminiscent of the 6 Hz variations of frequency and amplitude of the kHz QPOs reported by Yu, van der Klis and Jonker (2001).

We show that a luminous torus in the Schwarzschild metric oscillating along its own axis gives rise to a periodically varying flux of radiation, even though the source of radiation is steady and perfectly axisymmetric. This implies that the simplest oscillation mode in an accretion flow, axisymmetric up-and-down motion at the meridional epicyclic frequency, may be directly observable when it occurs in the inner parts of accretion flow around neutron stars and black holes. The high-frequency modulations of the X-ray flux observed in low-mass X-ray binaries at two frequencies (twin kHz QPOs) could then be a signature of strong gravity both because radial and meridional oscillations have different frequencies in non-Newtonian gravity, and because strong gravitational deflection of light rays causes the flux of radiation to be modulated at the higher frequency.

We have performed numerical simulations of a radially perturbed "accretion" torus around a black hole or neutron star and find that the torus performs radial and vertical motions at the appropriate epicyclic frequencies. We find clear evidence that vertical motions are excited in a non-linear resonance when the applied perturbation is periodic in time. The strongest resonant response occurs when the frequency difference of the two oscillations is equal to one-half the forcing frequency, precisely as recently observed in the accreting pulsar, SAX J1808.4-3658, where the observed kHz QPO peak separation is half the spin frequency of 401 Hz.

Recently, twin-peak QPOs have been observed in a 3:2 ratio for three Galactic black-hole microquasars with frequencies that have been shown to scale as 1/M, as expected for general relativisitic motion near a black hole. It may be possible to extend this result to distinguish between the following two disparate models that have been proposed for the puzzling ultraluminous X-ray sources (ULXs): (1) an intermediate-mass black hole M ~1000 solar mass) radiating very near the Eddington limit and (2) a conventional black hole (M ~ 10 solar mass) accreting at a highly super-Eddington rate with its emission beamed along the rotation axis. We suggest that it may be possible to distinguish between these models by detecting the counterpart of a Galactic twin-peak QPO in a ULX: the expected frequency for the intermediate-mass black hole model is only about 1 Hz, whereas, for the conventional black hole model the expected frequency would be the ~100 Hz value observed for the Galactic microquasars.

The presence of a kHz frequency in LMXBs has been expected from scaling laws, by analogy with the QPO phenomenon in HMXB X-ray pulsars. Interpretation of the two kHz frequencies, observed in accreting neutron stars, in terms of non-linear resonance in strong-field gravity led to the prediction of twin QPOs in black hole systems, in a definite frequency ratio (such as 2/3). The imprint of a subharmonic of the 401 Hz rotation rate in the frequencies of the QPOs detected in the accreting millisecond pulsar is at once a signature of non-linear resonance and of coupling between accretion disk modes and the neutron star spin.

Many Galactic black hole and neutron star sources in low X-ray mass binaries show QPOs (quasi periodic oscillations) in their observed X-ray fluxes, i.e. peaks in the Fourier variability power spectra. Pairs of twin peaks are observed, and in black-hole systems their frequencies (upp), (down) are in rational ratios. For example, in all four microquasars with twin peaks observed, (upp):(down) = 3:2. The rational ratios have been postulated in a model that explained twin peak QPOs as a non-linear resonance between modes of accretion disk oscillations. For those microquasars where the mass of the X-ray source is known, we determine the black-hole spin, following from the observed QPO frequencies within various models of resonance.

In all the microquasars with two hHz QPOs, the ratio of the frequencies is 3:2, supporting our suggestion that a non-linear resonance between two modes of oscillation in the accretion disk plays a role in exciting the observed modulations of the X-ray flux. We discuss the evidence in favor of this interpretation, and we relate the black hole spin to the frequencies expected for various types of resonances that may occur in nearly Keplerian disks in strong gravity. For those microquasars where the mass of the central X-ray source is known, the black hole spin can be deduced from a comparison of the observed and expected frequencies.

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.

Twin quasi-periodic millisecond modulations of the X-ray flux (kHz QPOs) have recently been reported from an accreting 2.5 ms X-ray pulsar. We identify modes of disk oscillations whose frequencies are in agreement with the observed ones when the rotating neutron star is modeled with realistic equations of state. The frequency difference of the twin QPOs, equal to about one half of the neutron-star spin rate, clearly indicates that resonant oscillations of the accretion disk have been observed. Similar non-linear resonances may also be spontaneously excited in the accretion disk. The two QPO frequencies in the pulsar system are close to a 5:7 ratio and this suggests a link with the QPOs in black hole systems, where frequency ratios of 2:3 and 3:5 have been reported.

Millisecond variability detected in low-mass X-ray binaries (LMXBs) probably reflects motions of accreting fluid in the strong-field gravity of neutron stars or neutron stars. Parametric resonance between two oscillators with unequal eigenfrequencies provides a natural explanation for the 2:3 frequency ratio observed in some black-hole systems, and probably also in the brightest neutron stars. The two oscillators likely correspond to the meridional and radial epicyclic motions.

We present three-dimensional MHD simulations of rotating radiatively inefficient accretion flows onto black holes. In the simulations, we continuously inject magnetized matter into the computational domain near the outer boundary, and we run the calculations long enough for the resulting accretion flow to reach a quasi-steady state. We have studied two limiting cases for the geometry of the injected magnetic field: pure toroidal field and pure poloidal field. In the case of toroidal field injection, the accreting matter forms a nearly axisymmetric, geometrically-thick, turbulent accretion disk. The disk resembles in many respects the convection-dominated accretion flows found in previous numerical and analytical investigations of viscous hydrodynamic flows. Models with poloidal field injection evolve through two distinct phases. In an initial transient phase, the flow forms a relatively flattened, quasi-Keplerian disk with a hot corona and a bipolar outflow. However, when the flow later achieves steady state, it changes in character completely. The magnetized accreting gas becomes two-phase, with most of the volume being dominated by a strong dipolar magnetic field from which a thermal low-density wind flows out. Accretion occurs mainly via narrow slowly-rotating radial streams which `diffuse' through the magnetic field with the help of magnetic reconnection events.

We have explored the ideas that parametric resonance affects nearly geodesic motion around a black hole or a neutron star, and that it may be relevant to the high frequency (twin) quasi-periodic oscillations occurring in some low-mass X-ray binaries. We have assumed the particles or fluid elements of an accretion disc to be subject to an isotropic perturbation of a hypothetical but rather general form. We find that the parametric resonance is indeed excited close to the radius where epicyclic frequencies of radial and meridional oscillations are in a 2:3 ratio. The location and frequencies of the highest amplitude excitation vary with the strength of the perturbation. These results agree with actual frequency ratios of twin kHz QPOs that have been reported in some black hole candidates, and they may be consistent also with correlation of the twin peaks in Sco X-1.

Recently, several ways of obtaining observational proof of the existence of black-hole horizons have been proposed. We argue here that such proof is fundamentally impossible: observations can provide arguments, sometimes very strong ones, in favour of the existence of the event horizon, but they cannot prove it. This applies also to future observations, which will trace very accurately the details of the spacetime metric of a body suspected of being a black hole.

We find evidence that the two high frequency QPOs in Sco X-1 are, more often than not, approximately in the 2:3 frequency ratio familiar from studies of black hole candidates (e.g., XTE J1550-564, Remillard et al. 2002). This implies that the double kHz QPO phenomenon in neutron stars has its origin in properties of strong-field gravity and has little to do with the rotation of a stellar surface or any magnetic field structure anchored in the star.

We apply Feynman's principle, ``The same equations have the same solutions'', to Kepler's problem and show that Newton's dynamics in a properly curved 3-D space is identical with that described by Einstein's theory in the 3-D optical geometry of Schwarzschild's spacetime. For this reason, rather unexpectedly, Newton's formulae for Kepler's problem, in the case of nearly circular motion in a static, spherically spherical gravitational potential accurately describe strong field general relativistic effects, in particular vanishing of the radial epicyclic frequency at the marginally stable orbit.

It is shown that optical geometry of the Reissner-Nordstrom exterior metric can be embedded in a hyperbolic space all the way down to its outer horizon. The adopted embedding procedure removes a breakdown of flat-space embeddings which occurs outside the horizon, at and below the Buchdahl-Bondi limit (R/M=9/4 in the Schwarzschild case). In particular, the horizon can be captured in the optical geometry embedding diagram. Moreover, by using the compact Poincare ball representation of the hyperbolic space, the embedding diagram can cover the whole extent of radius from spatial infinity down to the horizon. Attention is drawn to advantages of such embeddings in an appropriately curved space: this approach gives compact embeddings and it distinguishes clearly the case of an extremal black hole from a non-extremal one in terms of topology of the embedded horizon.

Radiatively inefficient accretion flows onto black holes are unstable due to both an outwardly decreasing entropy (`convection') and an outwardly decreasing rotation rate (the `magnetorotational instability'; MRI). Using a linear magnetohydrodynamic stability analysis, we show that long-wavelength modes are primarily destabilized by the entropy gradient and that such `convective' modes transport angular momentum inwards. Moreover, the stability criteria for the convective modes are the standard Hoiland criteria of hydrodynamics. By contrast, shorter wavelength modes are primarily destabilized by magnetic tension and differential rotation. These `MRI' modes transport angular momentum outwards. The convection-dominated accretion flow (CDAF) model, which has been proposed for radiatively inefficient accretion onto a black hole, posits that inward angular momentum transport and outward energy transport by long-wavelength convective fluctuations are crucial for determining the structure of the accretion flow. Our analysis suggests that the CDAF model is applicable to a magnetohydrodynamic accretion flow provided the magnetic field saturates at a sufficiently sub-equipartition value (plasma beta >> 1), so that long-wavelength convective fluctuations can fit inside the accretion disk. Numerical magnetohydrodynamic simulations are required to determine whether such a sub-equipartition field is in fact obtained.

Non-linear acoustic coupling of modes in accretion disks allows a parametric resonance between epicyclic motions. In black hole disks, such a resonance first occurs, and is strongest, when the radial and vertical epicyclic frequencies are approximately in a 2:3 ratio, in agreement with the 300 and 450 Hz frequencies reported in GRO J1655-40 and the 184 and 276 Hz frequencies reported in XTE J1550-564. The first overtone in the same resonance is in a 5:3 ratio with the fundamental, in agreement with the 69.2 and 41.5 Hz frequencies reported in GRS 1915+105. The narrow width and the frequency stability of the corresponding QPOs follow from general properties of parametric resonance.

Quantized orbital structures are typical for many aspects of classical gravity (Newton's as well as Einstein's). The astronomical phenomenon of orbital resonances is a well-known example. Recently, Rothman, Ellis and Murugan (2001) discussed quantized orbital structures in the novel context of a holonomy invariance of parallel transport in Schwarzschild geometry. We present here yet another example of quantization of orbits, reflecting both orbital resonances and holonomy invariance. This strong-gravity effect may already have been directly observed as the puzzling kilohertz quasi-periodic oscillations (QPOs) in the X-ray emission from a few accreting galactic black holes and several neutron stars.

We show that, at first order in the angular velocity, the general relativistic description of Rossby-Haurwitz waves (the analogues of r-waves on a thin shell) can be obtained from the corresponding Newtonian one after a coordinate transformation. As an application, we show that the results recently obtained by Rezzolla and Yoshida (2001) in the analysis of Newtonian Rossby-Haurwitz waves of a slowly and differentially rotating, fluid shell apply also in General Relativity, at first order in the angular velocity.

We consider the radial structure of radiatively inefficient hydrodynamic accretion flows around black holes. We show that low-viscosity flows consist of two zones: an outer convection-dominated zone and an inner advection-dominated zone. The transition between these two zones occurs at ~50 Schwarzschild radii.

Heyl (2000) made an interesting suggestion that the observed shifts in QPO frequency in type I X-ray bursts could be influenced by the same geometrical effect of strong gravity as the one that causes centrifugal force reversal discovered by Abramowicz and Lasota (1974). However, his main result contains a sign error. Here we derive the correct formula and conclude that constraints on the M(R) relation for neutron stars deduced from the rotational-modulation model of QPO frequency shifts are of no practical interest because the correct formula implies a weak condition R* > 1.3 Rs, where Rs is the Schwarzschild radius. We also argue against the relevance of the rotational-modulation model to the observed frequency modulations.

We note that the recently discovered 450 Hz frequency in the X-ray flux of the black hole candidate GRO J1655-40 is in a 3:2 ratio to the previously known 300 Hz frequency of quasi-periodic oscillations (QPO) in the same source. If the origin of high frequency QPOs in black hole systems is a resonance between orbital and epicyclic motion of accreting matter, as suggested previously, the angular momentum of the black hole can be accurately determined, given its mass. We find that the dimensionless angular momentum is in the range $0.2<j<0.65$ if the mass is in the (corresponding) range of 5.5 to 7.9 solar masses.

We explain the origin of the puzzling high frequency peaks (QPOs) in the variability power spectra of accreting neutron stars and black holes as a non-linear 1:2 or 1:3 resonance between orbital and radial epicyclic motion. These resonances are present because the gravitational field deviates strongly from a Newtonian 1/r potential. Our theory agrees with the recently reported observations of two QPOs, at 300 Hz and 450 Hz, in the black hole candidate GRO J1655-40.

Recent hydrodynamical simulations of radiatively inefficient black hole accretion flows with low viscosity have demonstrated that these flows differ significantly from those described by an advection-dominated model. The black hole flows are advection-dominated only in their inner parts, but convectively dominated at radii R>100R_g. In such flows, the radiative output comes mostly from the convection part, and the radiative efficiency is independent of accretion rate and equals ~0.001. This value gives a limit for how dim an accreting black hole could be. It agrees with recent Chandra observations which indicate that accreting black holes in low-mass X-ray binaries are by factor about 100 dimmer that neutron stars accreting with the same accretion rates.

Recent numerical simulations of radiatively inefficient accretion flows onto compact objects have shown that convection is a general feature in such flows. Dissipation of rotational and gravitational energies in the accretion flows results in inward increase of entropy and development of efficient convective motions. Convection-dominated accretion flows (CDAFs) have a structure that is modified significantly in comparison with the canonical advection-dominated and Bondi-like accretion flows. The flows are characterized by the flattened radial density profiles, ~R^{-1/2}, and have reduced mass accretion rates. Convection transports outward a significant amount of the released binding energy of the accretion flow. We discuss basic dynamical and observational properties of ADAFs using numerical models and self-similar analytical solutions.

The origin and the profile of the iron K-alpha line from active galactic nuclei is being studied. The model with a quasi-spherical distribution of clouds around a compact galactic nucleus is explored in which the intrinsic spectral feature around 6-7 keV is substantially affected by Comptonization of primary X-rays. The predicted spectral profiles are influenced mainly by motion of the clouds in strong gravitational field of the centre, and by orientation of the source with respect to the observer. The case of the Seyfert galaxy MCG-6-30-15 is examined in more detail and free parameters of the model are determined. A robust result is obtained: acceptable fits to the observed profile correspond to the spherical distribution of the clouds with strong self-obscuration effects and with rather high value of ionization parameter. This is a different situation from that considered in the model of a cold illuminated disc. Nevertheless, geometrically flat (disc-type) configurations are not rejected.

We analyse the radiation pressure instability of the standard, optically thisc disc and propose it as an explanation of the obseved variability of the superluminal microquasar GRS 1915+105. The quantitave agreement with the observed outburst amplitudes is obtained under the assumption that substantial part of the energy is carried away by the jet.

We argue that the proper way to treat thin-thick accretion-disk transitions should take into account the 2-D nature of the problem. We illustrate the physical inconsistency of the 1-D vertically integrated approach by discussing a particular example of the convective transport of energy.

Advection-dominated accretion flows (ADAFs) are known to be convectively unstable for low values of the viscosity parameter alpha. Two-dimensional axisymmetric hydrodynamic simulations of such flows reveal a radial density profile which is significantly flatter than the rho~r^{-3/2} expected for a canonical ADAF. The modified density profile is the result of inward transport of angular momentum by axisymmetric convective eddies. We present three-dimensional hydrodynamic simulations of convective ADAFs which are free of the assumption of axisymmetry. We find that the results are qualitatively and quantitatively similar to those obtained with two-dimensional simulations. In particular, the convective eddies are nearly axisymmetric and transport angular momentum inward.

We present a systematic numerical study of two-dimensional axisymmetric accretion flows around black holes. The flows have no radiative cooling and are treated in the framework of the hydrodynamical approximation. The models calculated in this study cover the large range of the relevant parameter space. There are four types of flows, determined by the values of the viscosity parameter $\alpha$ and the adiabatic index $\gamma$: convective flows, large-scale circulations, pure inflows and bipolar outflows. Thermal conduction introduces significant changes to the solutions, but does not create a new flow type. Convective accretion flows and flows with large-scale circulations have significant outward-directed energy fluxes, which have important implications for the spectra and luminosities of accreting black holes.

We consider height-integrated equations of an advection-dominated accretion flow (ADAF), assuming that there is no mass outflow. We include convection through a mixing length formalism. We seek self-similar solutions in which the rotational velocity and sound speed scale as R^{-1/2}, where R is the radius, and consider two limiting prescriptions for the transport of angular momentum by convection. In one limit, the transport occurs down the angular velocity gradient, so convection moves angular momentum outward. In the other, the transport is down the specific angular momentum gradient, so convection moves angular momentum inward. We also consider general prescriptions which lie in between the two limits. When convection moves angular momentum outward, we recover the usual self-similar solution for ADAFs in which the mass density scales as rho ~ R^{-3/2}. When convection moves angular momentum inward, the result depends on the viscosity coefficient alpha. If alpha>alpha_{crit1} ~ 0.05, we once again find the standard ADAF solution. For alpha<alpha_{crit}, however, we find a non-accreting solution in which rho ~ R^{-1/2}. We refer to this as a "convective envelope" solution or a "convection-dominated accretion flow". Two-dimensional numerical simulations of ADAFs with values of alpha<0.03 have been reported by several authors. The simulated ADAFs exhibit convection. By virtue of their axisymmetry, convection in these simulations moves angular momentum inward, as we confirm by computing the Reynolds stress. The simulations give rho ~ R^{-1/2}, in good agreement with the convective envelope solution. The R^{-1/2} density profile is not a consequence of mass outflow.

We show that recently published assertions that advection dominated accretion flows (ADAFs) require the presence of strong winds are unfounded because they assume that low radiative efficiency in flows accreting at low rates onto black holes implies vanishing radial energy and angular momentum fluxes through the flow (which is also formulated in terms of the `Bernoulli function' being positive). This, however, is a property only of self-similar solutions which are an inadequate representation of global accretion flows. We recall general properties of accretion flows onto black holes and show that such, necessarily transonic, flows may have either positive or negative Bernoulli function depending on the flow viscosity. Flows with low viscosities (alpha<0.1 in the alpha-viscosity model) have a negative Bernoulli function. Without exception, all 2-D and 1-D numerical models of low viscosity flows constructed to date experience no significant outflows. At high viscosities the presence of outflows depends on the assumed viscosity, the equation of state and on the outer boundary condition. The positive sign of the Bernoulli function invoked in this context is irrelevant to the presence of outflows. As an illustration, we recall 2-D numerical models with moderate viscosity that have positive values of the Bernoulli function and experience no outflows. ADAFs, therefore, do not differ from this point of view from thin Keplerian discs: they may have, but they do not have to have strong winds.

Observed spectra of Active Galactic Nuclei (AGN) and luminous X-ray binaries in our Galaxy suggest that both hot (~10^9 K) and cold (~10^6 K) plasma components exist close to the central accreting black hole. Hard X-ray component of the spectra is usually explained by Compton upscattering of optical/UV photons from optically thick cold plasma by hot electrons. Observations also indicate that some of these objects are quite efficient in converting gravitational energy of accretion matter into radiation. Existing theoretical models have difficulties in explaining the two plasma components and high intensity of hard X-rays. Most of the models assume that the hot component emerges from the cold one due to some kind of instability, but no one offers a satisfactory physical explanation for this. Here we propose a solution to these difficulties that reverses what was imagined previously: in our model the hot component forms first and afterward it cools down to form the cold component. In our model, accretion flow has initially a small angular momentum, and thus it has a quasi-spherical geometry at large radii. Close to the black hole, the accreting matter is heated up in shocks that form due to the action of the centrifugal force. The hot post-shock matter is very efficiently cooled down by Comptonization of low energy photons and condensates into a thin and cold accretion disk. The thin disk emits the low energy photons which cool the hot component.

We further investigate the runaway instability of relativistic tori around black holes assuming a power law distribution of the specific angular momentum. We neglect the self-gravity of the torus; thus the gravitational field of the system is determined by the central rotating black hole alone (the Kerr geometry). We find that for any positive (i.e. consistent with the local Rayleigh stability condition) power law index of the angular momentum distribution, the tori are always runaway stable for a sufficiently large spin of the black hole. However, before concluding that some (astrophysically realistic) tori could indeed be runaway stable, one should include, in full general relativity, the destabilizing effect of self-gravity.

We numerically construct slim, global, vertically integrated models of optically thin, transonic accretion discs around black holes, assuming a regularity condition at the sonic radius and boundary conditions at the outer radius of the disc and near the black hole. In agreement with several previous studies, we find two branches of shock-free solutions, in which the cooling is dominated either by advection, or by local radiation. We also confirm that the part of the accretion flow where advection dominates is in some circumstances limited in size: it does not extend beyond a certain outer limiting radius. New results found in our paper concern the location of the limiting radius and properties of the flow near to it. In particular, we find that beyond the limiting radius, the advective dominated solutions match on to Shapiro, Lightman & Eardley (SLE) discs through a smooth transition region. Therefore, the full global solutions are shock-free and unlimited in size. There is no need for postulating an extra physical effect (e.g. evaporation) for triggering the ADAF-SLE transition. It occurs due to standard accretion processes described by the classic slim disc equations.

In the innermost region of the advective accretion disk orbiting a black hole of high spin, the inertia of heat stored in the accreting gas is comparable to that of the gas rest mass itself. Accounting for this effect, we derive additional terms in the disk structure equations, and show that the heat inertia plays a significant role in the global energy conservation and dynamics of accretion in the relativistic advective disks.

The morphology of adiabatic accretion flows with small non-zero specific angular momentum has been investigated in the axisymmetric and nonviscous limit. For an initial state characterized by a Bondi flow with the specific angular momentum distributed with respect to polar angle, a travelling shock wave forms which propagates more rapidly in the equatorial plane than in the plane perpendicular to it resulting in the formation of a hot torus. In cases where the incoming flow is restricted to lie near the equatorial plane, a strong wind forms directed away from this plane with the tendency for the formation of a nonsteady shock structure. As the height of the incoming flow is increased it is found that the resultant wind weakens. The parameter regime which delineates accretion flows characterized by a travelling shock and a nearly standing shock structure is presented.

We re-derive the equation for vertical hydrodynamical equilibrium of stationary thin and slim accretion disks in the Kerr spacetime. All the previous derivations have been unsatisfactory, yielding unphysical singularities. Our equation is non-singular, more general, and simpler than those previously derived.

We obtained global transonic solutions representing optically thin advection-dominated accretion flow by solving the full set of differential equations describing such systems. We found that far from the sonic point self-similar solutions are an excellent approximation of the global flow structure if the accretion rate is well below the maximum value above which no optically thin solutions exist.

We derive all relevant equations needed for constructing a global general relativistic model of advectively cooled, very hot, optically thin accretion disks around black holes and present solutions which describe advection dominated flows in the gravitational field of a Kerr black hole.

A simple re-scaling of velocities calculated from the Paczy{\'n}ski \& Wiita (1980) pseudo-Newtonian potential makes them consistent with special relativity and greatly improves the agreement with the exact relativistic calculations. The improvement is relevant in calculations of the Doppler effect, spectra and radiation transfer which all involve effects of special relativity.

The mass of the central black hole in the active galaxy NGC 4258 (M106) has been measured to be $M=3.6\times10^7\Msun$ (Miyoshi et al. 1995). The Eddington luminosity corresponding to this mass is $L_E=4.5\times10^{45}$ erg s$^{-1}$. By contrast the X-ray luminosity of the nucleus of NGC 4258 between $2-10$ keV is $(4\pm 1)\times10^{40}~{\rm erg\,s^{-1}}$ while the optical/UV luminosity is less than $1.5\times10^{42} ~{\rm erg\,s^{-1}}$. The luminosity of NGC 4258 is therefore extremely sub-Eddington, $L\sim10^{-5}L_E$ in X-rays and $L\sim3\times10^{-4} L_E$ even if we take the maximum optical/UV luminosity. Assuming the usual accretion efficiency of 0.1 would imply accretion rates orders of magnitude lower than in Seyfert galaxies and quasars. We show that the properties of the AGN in NGC 4258 can be explained by an accretion flow in the form of a very hot, optically-thin plasma which advects most of the viscously generated thermal energy into the central black hole and radiates only a small fraction of the energy. In this case the accretion rate in Eddington units could be as high as $\sim 0.16\alpha$, where $\alpha$ is the standard viscosity parameter; and the size of the hot disk should be larger than $\sim 10$ times the Schwarzschild radius. We compare the predictions of this model with observations and discuss its application to other low luminosity AGN.

Thermal instability is examined for advection-dominated one-temperature accretion disks. We consider axisymmetric perturbations with short wavelength in the radial direction. The viscosity is assumed to be sufficiently small for the vertical hydrostatic balance to hold in perturbed states. The type of viscosity is given either by the $\alpha$-viscosity or by a diffusion-type stress tensor. Optically thick disks are found to be in general more unstable than optically thin ones. When the thermal diffusion is present, the optically thin disks become stable, but the optically thick disks are still unstable. The instability of the advection-dominated disks is different from that of the geometrically thin disks without advection. In the case of no advection, the thermal mode behaves under no appreciable surface density change. In the case of advection-dominated disks, however, the thermal mode occurs with no appreciable pressure change (compared with the density change), when local perturbations are considered. The variations of angular momentum and of surface density associated with the perturbations lead to a thermal instability. The astrophysical implications of this instability are briefly discussed.

Two-dimensional accretion flows near black holes have been investigated by time-dependent hydrodynamical calculations. We assume that the flow is axisymmetric and that radiative losses of internal energy are negligible, so that the disc is geometrically thick and hot. Accretion occurs due to the overflow of the effective potential barrier near the black hole, similar to the case of the Roche lobe overflowing star in a binary system. We make no pre-assumptions on the properties of the flow, instead our models evolve self-consistently from an initially non-accreting state. The viscosity is due to the the small-scale turbulence and it is described by the $\alpha$-viscosity prescription. We confirm earlier suggestions that viscous accretion flows are convectively unstable. We found that the instability produces transient eddies of various length-scales. The eddies contribute to the strength of the viscosity in the flow by redistributing the angular momentum. They also introduce low amplitude oscillatory variations which have a typical frequency about $100\; \msun/M$ Hz for a system of mass $M$. This may be relevant to the high frequency ($\sim 4-10$ Hz) quasi-periodic oscillations observed in the Galactic black hole candidate X-ray sources.

We reexamine the hypothesis that the optical/UV/soft X-ray continuum of Active Galactic Nuclei is thermal emission from an accretion disk. Previous studies have shown that fitting the spectra with the standard, optically thick and geometrically thin accretion disk models often led to luminosities which contradict the basic assumptions adopted in the standard model. There is no known reason why the accretion rates in AGN should not be larger than the thin disk limit. In fact, more general, slim accretion disk models are self-consistent even for moderately super-Eddington luminosities. We calculate here spectra from a set of thin and slim, optically thick accretion disks. We discuss the differences between the thin and slim disk models, stressing the implications of these differences for the interpretation of the observed properties of AGN. We found that the spectra can be fitted not only by models with a high mass and a low accretion rate (as in the case of thin disk fitting) but also by models with a low mass and a high accretion rate. In the first case fitting the observed spectra in various redshift categories gives black hole masses around 10^9 solar masses for a wide range of redshifts, and for accretion rates ranging from 0.4 to 8 solar masses/year. In the second case the accretion rate is around 10^2 solar masses/year for all AGN and the mass ranges from 3*10^6 to 10^8 solar masses. Unlike the disks with a low accretion rate, the spectra of the high-accretion-rate disks extend into the soft X-rays. A comparison with observations shows that such disks could produce the soft X-ray excesses claimed in some AGNs. We show also that the sequence of our models with fixed mass and different accretion rates can explain the time evolution of the observed spectra in Fairall 9.

We show that for realistic equations of state of dense matter, the universal proportionality factor relating the maximum rotation rate of neutron stars due to mass-shedding limit to the mass and radius of maximum allowable mass configuration of non-rotating models results from a universal proportionality between masses and radii of static maximum-mass neutron stars and those of maximally rotating configurations. These empirical relations cannot be obtained in the slow rotation approximation.

The global nonlinear time dependent evolution of accretion disk-corona systems in X-ray binary sources has been investigated to provide an understanding of the low frequency ($\sim 0.04$~Hz) quasi-periodic oscillations (QPOs) observed recently in the Rapid Burster MXB 1730-335 and in some black hole candidates sources (Cyg X--1 and GRO~J0422+32). We consider $\alpha$-viscosity models in which the viscous stress is proportional to the total pressure. In contrast to previous time dependent studies it is assumed that all mass accretion and angular momentum transport take place in an optically thick disk, but that a fraction of the gravitational energy that is released is dissipated in a corona. It is found that the coronal energy dissipation can effectively reduce the amplitudes of the mass flow variations generated from the thermal and viscous instabilities (in comparison to models without a corona) to a state in which the disk exhibits mild oscillatory nonsteady behavior. These oscillations are globally coherent in the unstable regions of the disk. A model for the high and low states of black hole candidate systems is also proposed. It is suggested that the low state, which is characterized by a hard X-ray spectrum, corresponds to a disk configuration in which the inner disk is in an advection dominated hot optically thin state whereas the high state corresponds to a configuration in which the inner disk is in an optically thick state surrounded by a corona. In this model, the mass accretion rate in the system is higher in the low state than in the high state. The hard X-ray spectrum of QPOs observed in the low state can be naturally explained by such a model.

We provide a unified description of thermal equilibria of black hole accretion disks, including the newly-discovered advection-dominated solutions. We classify the solutions on the basis of optical depth and importance of advection cooling. We demonstrate that only four physically distinct topological types of equilibria exist. Two of the types correspond to optically thin and optically thick equilibria, while the other two types are distinguished by whether advection is negligible or dominant. A stable Shakura-Sunyaev disk exists only for accretion rates $\dot M$ below a certain maximum. However, there is a critical viscosity parameter $\alpha_{\rm crit}$, which is a function of radius, such that for $\alpha>\alpha_{\rm crit}$ advection-dominated solutions exist for all $\dot M$. Even when $\alpha< \alpha_{\rm crit}$, the advection-dominated solutions are available for a wide range of $\dot M$ except for a gap around the Eddington rate. We therefore suggest that advection-dominated flows may be more common than standard thin disks in black hole systems. For certain ranges of radii and $\dot M$, no stable steady state solution is possible. In these cases, we suggest that limit cycle behavior may occur, leading to variability.

We show that most of hot, optically thin accretion disk models which ignore advective cooling are not self-consistent. We have found new types of optically thin disk solutions where cooling is dominated by radial advection of heat. These new solutions are thermally and viscously stable.

The Seyfert galaxy NGC 6814 is known to show periodic variation of its X-ray luminosity. We found that the sequences of peaks (variability patterns) in the folded X-ray light curves constructed from the {\it Exosat} and {\it Ginga} data are remarkably similar when one ignores amplitudes of the peaks and considers only their phases. The stable pattern consists of five peaks which are present in the both curves. The phases of the corresponding peaks coincide with an accuracy of about 10 degrees. The probability that this coincidence occurs by chance is less than about 1\% according to the most conservative estimate. The observed stable pattern of peaks may be produced by a stable distribution of ``bright spots'' on the accretion disk surface, {\it e.g.} by strong vortices or magnetic flux tubes.