We perform a bivariate Taylor expansion of the axisymmetric Green function in order to determine the exterior potential of a static thin toroidal shell having a circular section, as given by the Laplace equation. This expansion, performed at the centre of the section, consists in an infinite series in the powers of the minor-to-major radius ratio $e$ of the shell. It is appropriate for a solid, homogeneous torus, as well as for inhomogeneous bodies (the case of a core stratification is considered). We show that the leading term is identical to the potential of a loop having the same main radius and the same mass | this "similarity" is shown to hold in the ${\cal O}(e^2)$ order. The series converges very well, especially close to the surface of the toroid where the average relative precision is $\sim 10^{-3}$ for $e\! = \!0.1$ at order zero, and as low as a few $10^{-6}$ at second order. The Laplace equation is satisfied {\em exactly} in every order, so no extra density is induced by truncation. The gravitational acceleration, important in dynamical studies, is reproduced with the same accuracy. The technique also applies to the magnetic potential and field generated by azimuthal currents as met in terrestrial and astrophysical plasmas.

Complete integrability of geodesic motion, the well known feature of the fields of isolated stationary black holes, can easily be "spoilt" by the presence of some additional source (even if highly symmetric). In previous papers, we used various methods to show how free time-like motion becomes chaotic if the gravitational field of the Schwarzschild black hole is perturbed by that of a circular disc or ring, considering specifically the inverted first disc of the Morgan-Morgan counter-rotating family and the Bach-Weyl ring as the additional sources. The present paper focuses on two new points. First, since the Bach-Weyl thin ring is physically quite unsatisfactory, we now repeat some of the analysis for a different, Majumdar-Papapetrou--type (extremally charged) ring around an extreme Reissner-Nordstr\"om black hole, and compare the results with those obtained before. We also argue that such a system is in fact more relevant astrophysically than it may seem. Second, we check numerically, for the latter system as well as for the Schwarzschild black hole encircled by the inverted Morgan-Morgan disc, how indicative is the geometric (curvature) criterion for chaos suggested by Sota, Suzuki & Maeda (1996). We also add a review of the literature where the relevance of geometric criteria in general relativity (as well as elsewhere) has been discussed for decades.

In ApJ (1974), Will solved the perturbation of a Schwarzschild black hole due to a slowly rotating light concentric thin ring, using Green's functions expressed as infinite-sum expansions in multipoles and in the small mass and rotational parameters. In a previous paper (ApJS 2017, Paper I), we expressed the Green functions in closed form containing elliptic integrals, leaving just summation over the mass expansion. Such a form is more practical for numerical evaluation, but mainly for generalizing the problem to extended sources where the Green functions have to be integrated over the source. We exemplified the method by computing explicitly the first-order perturbation due to a slowly rotating thin disc lying between two finite radii. After finding basic parameters of the system -- mass and angular momentum of the black hole and of the disc -- we now add further properties, namely those which reveal how the disc gravity influences geometry of the black-hole horizon and those of circular equatorial geodesics (specifically, radii of the photon, marginally bound and marginally stable orbits). We also realize that, in the linear order, no ergosphere occurs and the central singularity remains point-like, and check the implications of natural physical requirements (energy conditions and subluminal restriction on orbital speed) for the single-stream as well as counter-rotating double-stream interpretations of the disc.

Will (1974) treated the perturbation of a Schwarzschild black hole due to a slowly rotating light concentric thin ring by solving the perturbation equations in terms of a multipole expansion of the mass-and-rotation perturbation series. In the Schwarzschild background, his approach can be generalized to the perturbation by a thin disc (which is more relevant astrophysically), but, due to a rather bad convergence properties, the resulting expansions are not suitable for specific (numerical) computations. However, we show that Green's functions represented by the Will's result can be expressed in a closed form (without multipole expansion) which is more useful. In particular, they can be integrated out over the source (thin disc in our case), to yield well converging series both for the gravitational potential and for the dragging angular velocity. The procedure is demonstrated, in the first perturbation order, on the simplest case of a constant-density disc, including physical interpretation of the results in terms of a one-component perfect fluid or a two-component dust on circular orbits about the central black hole. Free parameters are chosen in such a way that the resulting black hole has zero angular momentum but non-zero angular velocity, being just carried along by the dragging effect of the disc.

"Almost all" seems to be known about isolated stationary black holes in asymptotically flat space-times and about the behaviour of {\em test} matter and fields in their backgrounds. The black holes likely present in galactic nuclei and in some X-ray binaries are commonly being represented by the Kerr metric, but actually they are not isolated (they are detected only thanks to a strong interaction with the surroundings), they are not stationary (black-hole sources are rather strongly variable) and they also probably do not live in an asymptotically flat universe. Such "perturbations" may query the classical black-hole theorems (how robust are the latter against them?) and certainly affect particles and fields around, which can have observational consequences. In the present contribution we examine how the geodesic structure of the static and axially symmetric black-hole space-time responds to the presence of an additional matter in the form of a thin disc or ring. We use several different methods to show that geodesic motion may become chaotic, to reveal the strength and type of this irregularity and its dependence on parameters. The relevance of such an analysis for galactic nuclei is briefly commented on.

We continue the study of time-like geodesic dynamics in exact static, axially and reflection symmetric space-times describing the fields of a Schwarzschild black hole surrounded by thin discs or rings. In the first paper of this series, the rise (and decline) of geodesic chaos with ring/disc mass and position and with test particle energy was revealed on Poincar\'e sections and on time series of position or velocity and their power spectra. In the second paper we compared these results with those obtained by two recurrence methods, focusing on "sticky" orbits whose different parts show different degrees of chaoticity. Here we complement the analysis by using several Lyapunov-type coefficients which quantify the rate of orbital divergence. After comparing the results with those obtained by the previous methods, we specifically consider a system involving a black hole surrounded by a small thin disc or a large ring, having in mind the configuration which probably occurs in galactic nuclei. Within the range of parameters which roughly corresponds to our Galactic center, we found that the black-hole accretion disc does not have a significant gravitational effect on the dynamics of free motion at larger radii, while the inner circumnuclear molecular ring (concentrated above 1 parsec radius) can only induce some irregularity in motion of stars ("particles") on smaller radii if its mass reaches 10 to 30% of the central black hole (which is the upper estimate given in the literature), if it is sufficiently compact (which does not hold but maybe for its inner rim) and if the stars can get to its close vicinity. The outer dust ring between 60 and 100 parsecs appears to be less important for the geodesic dynamics in its interior.