We integrate numerically axially symmetric stationary Einstein equations describing self-gravitating disks around spinless black holes. The numerical scheme is based on a method developed by Shibata, but contains important new ingredients. We derive a new general-relativistic Keplerian rotation law for self-gravitating disks around spinning black holes. Former results concerning rotation around spin-less black holes emerge in the limit of a vanishing spin parameter. These rotation curves might be used for the description of rotating stars, after appropriate modification around the symmetry axis. They can be applied to the description of compact torus--black hole configurations, including active galactic nuclei or products of coalescences of two neutron stars.

We formulate new general-relativistic extensions of Newtonian rotation laws for self-gravitating stationary fluids. They have been used to re-derive, in the first post-Newtonian approximation, the well known geometric dragging of frames. We derive two other general-relativistic weak-field effects within rotating tori: the recently discovered dynamic anti-dragging and a new effect that measures the deviation from the Keplerian motion and/or the contribution of the fluids selfgravity. One can use the rotation laws to study the uniqueness and the convergence of the post-Newtonian approximations, and the existence of the post-Newtonian limits.

In this paper we investigate a steady accretion within the Einstein-Straus vacuole, in the presence of the cosmological constant. The dark energy damps the mass accretion rate and --- above certain limit --- completely stops the steady accretion onto black holes, which in particular is prohibited in the inflation era and after (roughly) $10^{12}$ years from Big Bang (assuming the presently known value of the cosmological constant). Steady accretion would not exist in the late phases of the Penrose's scenario - known as the Weyl curvature hypothesis - of the evolution of the Universe.

We investigate accreting disk systems with polytropic gas in Keplerian motion. Numerical data and partial analytic results show that the self-gravitation of the disk speeds up its rotation -- its rotational frequency is larger than that given by the well known strictly Keplerian formula that takes into account the central mass only. Thus determination of central mass in systems with massive disks requires great care -- the strictly Keplerian formula yields only an upper bound. The effect of self-gravity depends on geometric aspects of disk configurations. Disk systems with a small (circa $10^{-4}$) ratio of the innermost radius to the outermost disk radius have the central mass close to the upper limit, but if this ratio is of the order of unity then the central mass can be smaller by many orders of magnitude from this bound.

The Keplerian motion of accretion disks in Active Galactic Nuclei (AGN) is usually believed to be generated by a heavy central mass. We investigate accreting disk systems with polytropic gas in Keplerian rotation and obtain a phenomenological formula that relates the Keplerian angular frequency to the ratio of disk and central masses. Central mass approaches the Keplerian value, if the inner boundary of a disk is close to the minimal stable orbit of a black hole. These results are applied to NGC 4258, the unique AGN with a finely measured Keplerian rotation curve of the central disk, with the conclusion that its rotation curve is, in fact, determined by the central black hole. The mass of the accretion disk exceeds 100 solar masses.

We studied a steadily accreting, geometrically thick disk model that selfconsistently takes into account selfgravitation of the polytropic gas, its interaction with the radiation and the mass accretion rate. The accreting mass is injected inward in the vicinity of the central $z=0$ plane, where also radiation is assumed to be created. The rest of the disk remains approximately stationary. Only conservation laws are employed and the gas-radiation interaction in the bulk of the disk is described in the thin-gas approximation. We demonstrate that this scheme is numerically viable and yields a structure of the bulk that is influenced by the radiation and (indirectly) by the prescribed mass accretion rate. The obtained disk configurations are typical for environments in Active Galactic Nuclei (AGN), with the central mass of the order of $10^7 M_{\astrosun}$ to $10^8 M_{\astrosun}$, quasi-Keplerian rotation curves, disk masses ranging from about $10^6 M_{\astrosun}$ to $10^7 M_{\astrosun}$, and the luminosity ranging from $10^6 L_{\astrosun}$ to $10^9 L_{\astrosun}$. These luminosities are much lower than the corresponding Eddington limit.

We estimate the minimal mass of selfgravitating polytropic disks using the famous Sobolev inequality. This bound resembles the well known mass formula for Lane-Emden stars. For ideal gas with the polytropic index n = 3 the minimal mass is not smaller than the Jeans mass. The accuracy of the estimate is verified in a number of numerical examples. The bound works well for heavy selfgravitating disks and is less useful for light disks.

We compare Newtonian and general relativistic descriptions of the stationary accretion of self-gravitating fluids onto compact bodies. Spherical symmetry and thin gas approximation are assumed. Luminosity depends, amongst other factors, on the temperature and the contribution of gas to the total mass, in both -- general relativistic ($L_{GR}$) and Newtonian ($L_N$) -- models. We discover a remarkable universal behaviour for transonic flows: the ratio of respective luminosities $L_{GR}/L_N$ is independent of the fractional mass of the gas and depends on asymptotic temperature. It is close to 1 in the regime of low asymptotic temperatures and can grow by one order of magnitude for high temperatures. These conclusions are valid for a wide range of polytropic equations of state.

We analyze stationary accretion of selfgravitating gas onto a compact center within general-relativistic radiation hydrodynamics. Spherical symmetry and thin gas approximation are assumed. Numerical investigation shows that transonic flows exist for small redshifts and they cease to exist for high redshifts and high luminosities. There exist two branches of flows (subsonic or supersonic) that originate at a bifurcation point and that embrace the set of subsonic solutions. The morphology of the set of subsonic solutions is essentially independent of redshifts and flows that belong to their boundary provide estimates of the gas abundance of subsonic solutions. It appears that prescribed boundary data guarantee uniqueness only of the bifurcation point, and that the latter has maximal luminosity.

Quasi-stationary flows of gas accreting onto a compact center are analyzed in the framework of general-relativistic radiation hydrodynamics, under assumptions of spherical symmetry and thin gas approximation. Numerical investigation shows that luminosity, redshift and gas abundance are correlated. The gas can constitute up to one third of the total mass of brightest low-redshift sources, but its abundance goes down to 1/30 for sources with luminosities close to the Eddington limit.

We investigate the newtonian stationary accretion of a polytropic perfect fluid onto a central body with a hard surface. The selfgravitation of the fluid and its interaction with luminosity is included in the model. We find that for a given luminosity, asymptotic mass and temperature of the fluid there exist two solutions with different cores.

We investigate newtonian description of accreting compact bodies with hard surfaces, including luminosity and selfgravitation of polytropic perfect fluids. This nonlinear integro-differential problem reduces, under appropriate boundary conditions, to an algebraic relation between luminosity and the gas abundance in stationary spherically symmetric flows. There exist, for a given luminosity, asymptotic mass and the asymptotic temperature, two sub-critical solutions that bifurcate from a critical point. They differ by the fluid content and the mass of the compact centre.

We investigate the accretion onto luminous bodies with hard surfaces within the framework of newtonian theory. The accreting gases are assumed to be polytropes and their selfgravitation is included. A remarkable feature of the model is that under proper boundary conditions some parameters of the sonic point are the same as in the Bondi model and the relation between luminosity and the gas abundance reduces to an algebraic relation. All that holds under assumptions of stationarity and spherical symmetric. Assuming data that are required for the complete specification of the system, one finds that generically for a given luminosity there exist two solutions with different compact cores.

We study the propagation of axial gravitational waves in Friedman universes. The evolution equation is obtained in the Regge-Wheeler gauge. The gravitational waves obey the Huygens principle in the radiation dominated era, but in the matter dominated universe their propagation depends on their wavelengths, with the scale fixed essentially by the Hubble radius. Short waves practically satisfy the Huygens principle while long waves can backscatter off the curvature of a spacetime.

Corrected misprints and omissions.

The gravitational redshift formula is usually derived in the geometric optics approximation. In this note we consider an exact formulation of the problem in the Schwarzschild space-time, with the intention to clarify under what conditions this redshift law is valid. It is shown that in the case of shocks the radial component of the Poynting vector can scale according to the redshift formula, under a suitable condition. If that condition is not satisfied, then the effect of the backscattering can lead to significant modifications. The obtained results imply that the energy flux of the short wavelength radiation obeys the standard gravitational redshift formula while the energy flux of long waves can scale differently, with redshifts being dependent on the frequency.

Gravitational waves in the linear approximation propagate in the Schwarzschild spacetime similarly as electromagnetic waves. A fraction of the radiation scatters off the curvature of the geometry. The energy of the backscattered part of an initially outgoing pulse of the quadrupole gravitational radiation is estimated by compact formulas depending on the initial energy, the Schwarzschild radius, and the location and width of the pulse. The backscatter becomes negligible in the short wavelength regime.

Electromagnetic waves propagate in the Schwarzschild spacetime like in a nonuniform medium with a varying refraction index. A fraction of the radiation scatters off the curvature of the geometry. The energy of the backscattered part of an initially outgoing pulse of electromagnetic radiation can be estimated, in the case of dipole radiation, by a compact formula depending on the initial energy, the Schwarzschild radius and the pulse location. The magnitude of the backscattered energy depends on the frequency spectrum of the initial configuration. This effect becomes negligible in the short wave limit, but it can be significant in the long wave regime. Similar results hold for the massless scalar fields and are expected to hold also for weak gravitational waves.

We describe general-relativistically a spherically symmetric stationary fluid accretion onto a black hole. Relativistic effects enhance mass accretion, in comparison to the Bondi model predictions, in the case when backreaction is neglected. That enhancement depends on the adiabatic index and the asymptotic gas temperature and it can magnify accretion by one order in the ultrarelativistic regime.

This article investigates the interaction of a spherically symmetric massless scalar field with a strong gravitational field. It focuses on the propagation of waves in regions outside any horizons. The two factors acting on the waves can be identified as a redshift and a backscattering. The influence of backscattering on the intensity of the outgoing radiation is studied and rigorous quantitative upper bounds obtained. These show that the total flux may be decreased if the sources are placed in a region adjoining an apparent horizon. Backscattering can be neglected in the case $2m_0 /R<< 1$, that is when the emitter is located at a distance from a black hole much larger than the Schwarzschild radius. This backscattering may have noticeable astrophysical consequences.

The Cosmological Principle states that the universe is both homogeneous and isotropic. This, alone, is not enough to specify the global geometry of the spacetime. If we were able to measure both the Hubble constant and the energy density we could determine whether the universe is open or closed. Unfortunately, while some agreement exists on the value of the Hubble constant, the question of the energy density seems quite intractable. This Letter describes a possible way of avoiding this difficulty and shows that if one could measure the rate at which light-rays emerging from a surface expand, one might well be able to deduce whether the universe is closed.

This article gives necessary and sufficient conditions for the formation of trapped surfaces in spherically symmetric initial data defined on a closed manifold. Such trapped surfaces surround a region in which there occurs an enhancement of matter over the average. The conditions are posed directly in terms of physical variables and show that what one needs is a relatively large amount of excess matter confined to a small volume. The expansion of the universe and an outward flow of matter oppose the formation of trapped surfaces; an inward flow of matter helps. The model can be regarded as a Friedmann-Lema\^\i tre-Walker cosmology with localized spherical inhomogeneities. We show that the total excess mass cannot be too large.