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.

We have investigated the toroidal analog of ellipsoidal shells of matter, which are of great significance in Astrophysics. The exact formula for the gravitational potential $\Psi(R,Z)$ of a shell with a circular section at the pole of toroidal coordinates is first established. It depends on the mass of the shell, its main radius and axis-ratio $e$ (i.e. core-to-main radius ratio), and involves the product of the complete elliptic integrals of the first and second kinds. Next, we show that successive partial derivatives $\partial^{n +m} \Psi/\partial_{R^n} \partial_{Z^m}$ are also accessible by analytical means at that singular point, thereby enabling the expansion of the interior potential as a bivariate series. Then, we have generated approximations at orders $0$, $1$, $2$ and $3$, corresponding to increasing accuracy. Numerical experiments confirm the great reliability of the approach, in particular for small-to-moderate axis ratios ($e^2 \lesssim 0.1$ typically). In contrast with the ellipsoidal case (Newton's theorem), the potential is not uniform inside the shell cavity as a consequence of the curvature. We explain how to construct the interior potential of toroidal shells with a thick edge (i.e. tubes), and how a core stratification can be accounted for. This is a new step towards the full description of the gravitating potential and forces of tori and rings. Applications also concern electrically-charged systems, and thus go beyond the context of gravitation.

We discuss the effects of electric charging on the equilibrium configurations of magnetized, rotating fluid tori around black holes of different mass. In the context of gaseous/dusty tori in galactic nuclei, the central black hole dominates the gravitational field and it remains electrically neutral, while the surrounding material acquires some electric charge and exhibits non-negligible self-gravitational effect on the torus structure. The structure of the torus is influenced by the balance between the gravitational and electromagnetic forces. A cusp may develop even in Newtonian tori due to the charge distribution.

Equilibrium configurations of electrically charged perfect fluid surrounding a central rotating black hole endowed with a test electric charge and embedded in a large-scale asymptotically uniform magnetic field are presented. Following our previous studies considering the central black hole non-rotating, we show that in the rotating case, conditions for the configurations existence change according to the spin of the black hole. We focus our attention on the charged fluid in rigid rotation which can form toroidal configurations centered in the equatorial plane or the ones hovering above the black hole, along the symmetry axis. We conclude that a non-zero value of spin changes the existence conditions and the morphology of the solutions significantly. In the case of fast rotation, the morphology of the structures is close to an oblate shape.

We study stationary, electrically charged fluid structures encircling a rotating compact object with a dipole magnetic field oriented along the rotation axis. This situation is described in an idealized way by the Kerr metric and a magnetic dipole "test" field, that does not affect the spacetime. The self-gravitational and self-electromagnetic field of the fluid are neglected and the fluid is assumed to be non conductive and in rigid motion. Our work generalizes a previous study by Kov\'a\v{r} et al.(2016) by taking into account the rotation of the central object. Therefore, we focus on the influence of the rotation onto the existence and position of bound fluid structures. Frame dragging effects allow the existence of polar clouds, which could not be found in non-rotating case. Furthermore counter-rotating equatorial tori become more preferred the faster the central object is spinning.

We present an analytical approach for the equilibrium of a self-gravitating charged fluid embedded in a spherical gravitational and dipolar magnetic fields produced by a central mass. Our scheme is proposed, as a toy-model, in the context of gaseous/dusty tori surrounding supermassive black holes in galactic nuclei. While the central black hole dominates the gravitational field and it remains electrically neutral, the surrounding material has a non-negligible self-gravitational effect on the torus structure. By charging mechanisms it also acquires non-zero electric charge density, so the two influences need to be taken into account to achieve a self-consistent picture. With our approach we discuss the impact of self-gravity, represented by the term dt (ratio of the torus total mass to the mass of the central body), on the conditions for existence of the equilibrium and the morphology and typology of the tori. By comparison with a previous work without self-gravity, we show that the conditions can be different. Although the main aim of the present paper is to discuss a framework for the classification of electrically charged, magnetized, self-gravitating tori, we also mention potential astrophysical applications to vertically stratified fluid configurations.

The softened point mass model is commonly used in simulations of gaseous discs including self-gravity while the value of associated length \lambda remains, to some degree, controversial. This ``parameter'' is however fully constrained when, in a discretized disc, all fluid cells are demanded to obey Newton's law. We examine the topology of solutions in this context, focusing on cylindrical cells more or less vertically elongated. We find that not only the nominal length depends critically on the cell's shape (curvature, radial extension, height), but it is either a real or an imaginary number. Setting \lambda as a fraction of the local disc thickness -- as usually done -- is indeed not the optimal choice. We then propose a novel prescription valid irrespective of the disc properties and grid spacings. The benefit, which amounts to 2-3 more digits typically, is illustrated in a few concrete cases. A detailed mathematical analysis is in progress.

We study charged-fluid toroidal structures surrounding a non-rotating charged black hole immersed in a large-scale, asymptotically uniform magnetic field. In continuation of our former study on electrically charged matter in approximation of zero conductivity, we demonstrate existence of orbiting structures in permanent rigid rotation in the equatorial plane, and charged clouds hovering near the symmetry axis. We constrain the range of parameters that allow stable configurations and derive the geometrical shape of equi-pressure surfaces. Our simplified analytical study suggests that these regions of stability may be relevant for trapping electrically charged particles and dust grains in some areas of the black hole magnetosphere, being thus important in some astrophysical situations.

The local character of self-gravity along with the number of spatial dimensions are critical issues when computing the potential and forces inside massive systems like stars and disks. This appears from the discretisation scale where each cell of the numerical grid is a self-interacting body in itself. There is apparently no closed-form expression yet giving the potential of a three-dimensional homogeneous cylindrical or spherical cell, in contrast with the Cartesian case. By using Green's theorem, we show that the potential integral for such polar-type 3D sectors -- initially, a volume integral with singular kernel -- can be converted into a regular line-integral running over the lateral contour, thereby generalising a formula already known under axial symmetry. It therefore is a step towards the obtention of another potential/density pair. The new kernel is a finite function of the cell's shape (with the simplest form in cylindrical geometry), and mixes incomplete elliptic integrals, inverse trigonometric and hyperbolic functions. The contour integral is easy to compute; it is valid in the whole physical space, exterior and interior to the sector itself and works in fact for a wide variety of shapes of astrophysical interest (e.g. sectors of tori or flared discs). This result is suited to easily providing reference solutions, and to reconstructing potential and forces in inhomogeneous systems by superposition. The contour integrals for the 3 components of the acceleration vector are explicitely given.

As hydrostatic equilibrium of gaseous discs is partly governed by the gravity field, we have estimated the component caused by a vertically homogeneous disc, with a special attention for the outer regions where self-gravity classically appears. The accuracy of the integral formula is better than 1%, whatever the disc thickness, radial extension and radial density profile. At order zero, the field is even algebraic for thin discs and writes $- 4 \pi G \Sigma(R) f_{edge} (R)$ at disc surface, thereby correcting Paczynski's formula by a multiplying factor $f_{edge} \gtrsim 1/2$, which depends on the relative distance to the edges and the local disc thickness. For very centrally condensed discs however, this local contribution can be surpassed by action of mass stored in the inner regions, possibly resulting in $f_{edge} \gg 1$. A criterion setting the limit between these two regimes is derived. These result are robust in the sense that the details of vertical stratification are not critical. We briefly discuss how hydrostatic equilibrium is impacted. In particular, the disc flaring should not reverse in the self-gravitating region, which contradicts what is usually obtained from Paczynski's formula. This suggests that i) these outer regions are probably not fully shadowed by the inner ones (important when illuminated by a central star), and ii) the flared shape of discs does not firmly prove the absence or weakness of self-gravity.

By using various properties of the complete elliptic integrals, we have derived an alternative expression for the gravitational potential of axially symmetric bodies, which is free of singular kernel in contrast with the classical form. This is mainly a radial integral of the local surface density weighted by a regular "mean Green function" which depends explicitly on the body's vertical thickness. Rigorously, this result stands for a wide variety of configurations, as soon as the density structure is vertically homogeneous. Nevertheless, the sensitivity to vertical stratification | the Gaussian profile has been considered | appears weak provided that the surface density is conserved. For bodies with small aspect ratio (i.e. geometrically thin discs), a first-order Taylor expansion furnishes an excellent approximation for this mean Green function, the absolute error being of the fourth order in the aspect ratio. This formula is therefore well suited to studying the structure of self-gravitating discs and rings in the spirit of the "standard model of thin discs" where the vertical structure is often ignored, but it remains accurate for discs and tori of finite thickness. This approximation which perfectly saves the properties of Newton's law everywhere (in particular at large separations), is also very useful for dynamical studies where the body is just a source of gravity acting on external test particles.