When a black hole is accreting well below the Eddington rate, a geometrically thick, radiatively inefficient state of the accretion disk is established. There is a limited number of closed-form physical solutions for geometrically thick (non-selfgravitating) toroidal equilibria of perfect fluids orbiting a spinning black hole, and these are predominantly used as initial conditions for simulations of accretion in the aforementioned mode. However, different initial configurations might lead to different results and thus observational predictions drawn from such simulations. We expand the known equilibria by a number of closed multiparametric solutions with various possibilities of rotation curves and geometric shapes. Then, we ask whether choosing these as initial conditions influences the onset of accretion and the asymptotic state of the disk. We investigate a set of examples from the derived solutions in detail; we analytically estimate the growth of the magneto-rotational instability (MRI) from their rotation curves and evolve the analytically obtained tori using the 2D magneto-hydrodynamical code HARM. Properties of the evolutions are then studied through the mass, energy, and angular-momentum accretion rates. The rotation curve has a decisive role in the numerical onset of accretion in accordance with our analytical MRI estimates: In the first few orbital periods, the average accretion rate is linearly proportional to the initial MRI rate in the toroids. The final state obtained from any initial condition within the studied class after an evolution of $\gtrsim 10$ orbital periods is mostly qualitatively identical and the quantitative properties vary within a single order of magnitude. The average values of the energy of the accreted fluid have an irregular dependency on initial data, and in some cases fluid with energies many times its rest mass is systematically accreted.

We consider circular motion in the NUT (Newman-Unti-Tamburino) space-time. Among other things, we determine the location of circular time-like geodesic orbits, in particular of the innermost stable circular orbit (ISCO) and of the marginally bound circular orbit. Moreover, we discuss the von Zeipel cylinders with respect to the stationary observers and with respect to the Zero Angular Momentum Observers (ZAMOs). We also investigate the relation of von Zeipel cylinders to inertial forces, in particular in the ultra-relativistic limit. Finally, we generalise the construction of thick accretion tori ("Polish doughnuts") which are well known on the Schwarzschild or Kerr background to the case of the NUT metric. We argue that, in principle, a NUT source could be distinguished from a Schwarzschild or Kerr source by observing the features of circular matter flows in its neighbourhood.

We consider the motion of classical spinning test particles in Schwarzschild and Kerr metrics and investigate innermost stable circular orbits (ISCO). The main goal of this work is to find analytically the small-spin corrections for the parameters of ISCO (radius, total angular momentum, energy, orbital angular frequency) of spinning test particles in the case of vectors of black hole spin, particle spin and orbital angular momentum being collinear to each other. We analytically derive the small-spin linear corrections for arbitrary Kerr parameter $a$. The cases of Schwarzschild, slowly rotating and extreme Kerr black hole are considered in details. For a slowly rotating black hole the ISCO parameters are obtained up to quadratic in $a$ and particle's spin $s$ terms. From the formulae obtained it is seen that the spin-orbital coupling has attractive character when spin and angular momentum are parallel and repulsive when they are antiparallel. For the case of the extreme Kerr black hole with co-rotating particle we succeed to find the exact analytical solution for the limiting ISCO parameters for arbitrary spin. It has been shown that the limiting values of ISCO radius and frequency do not depend on the particle's spin while values of energy and total angular momentum depend on it. We have also considered circular orbits of arbitrary radius and have found small-spin linear corrections for the total angular momentum and energy at given radius. System of equations for numerical calculation of ISCO parameters for arbitrary $a$ and $s$ is also explicitly written.