Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

The spin-curvature coupling as captured by the so-called Mathisson-Papapetrou-Dixon (MPD) equations is the leading order effect of the finite size of a rapidly rotating compact astrophysical object moving in a curved background. It is also a next-to-leading order effect in the phase of gravitational waves emitted by extreme-mass-ratio inspirals (EMRIs), which are expected to become observable by the LISA space mission. Additionally, exploring the Hamiltonian formalism for spinning bodies is important for the construction of the so-called Effective-One-Body waveform models that should eventually cover all mass ratios. The MPD equations require supplementary conditions determining the frame in which the moments of the body are computed. We review various choices of these supplementary spin conditions and their properties. Then, we give Hamiltonians either in proper-time or coordinate-time parametrization for the Tulczyjew-Dixon, Mathisson-Pirani, and Kyrian-Semer\'ak conditions. Finally, we also give canonical phase-space coordinates parametrizing the spin tensor. We demonstrate the usefulness of the canonical coordinates for symplectic integration by constructing Poincar\'e surfaces of section for spinning bodies moving in the equatorial plane in Schwarzschild space-time. We observe the motion to be essentially regular for EMRI-ranges of the spin, but for larger values the Poincar\'e surfaces of section exhibit the typical structure of a weakly chaotic system. A possible future application of the numerical integration method is the inclusion of spin effects in EMRIs at the precision requirements of LISA.

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.

The Kerr black hole is stationary and axisymmetric, which leads to conservation of energy and azimuthal angular momentum along the orbits of free test particles in its vicinity, but also to conservation laws for the evolution of continuum matter fields. However, the Kerr space-time possesses an additional "hidden symmetry" which exhibits itself in an unexpected conserved quantity along geodesics known as the Carter constant. We investigate the possibility of using this hidden symmetry to obtain conservation laws and other identities which could be used to test astrophysical simulations of the evolution of matter fields near spinning black holes. After deriving such identities, we set up a simple numerical toy model on which we demonstrate how they can detect the violations of evolution equations in a numerical simulation. Even though one of the expressions we derive is in the form of a conservation law, we end up recommending an equivalent but simpler expression that is not in the form of a conservation law for practical implementation.

Pseudo-Newtonian potentials are a tool often used in theoretical astrophysics to capture some key features of a black-hole space-time in a Newtonian framework. As a result, one can use Newtonian numerical codes, and Newtonian formalism in general, in an effective description of important astrophysical processes such as accretion onto black holes. In this paper we develop a general pseudo-Newtonian formalism which pertains to the motion of particles, light, and fluids in stationary space-times. In return, we are able to assess the applicability of the pseudo-Newtonian scheme. The simplest and most elegant formulas are obtained in space-times without gravitomagnetic effects, such as the Schwarzschild rather than the Kerr space-time; the quantitative errors are smallest for motion with low binding energy. Included is a ready-to-use set of fluid equations in Schwarzschild space-time in Cartesian and radial coordinates.

Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydrodynamics. When Noether symmetries associated with Killing vectors or tensors are present in geodesic flows, they lead to constants of motion polynomial in the momenta. We generalize these concepts to hydrodynamic flows. Moreover, the Hamiltonian descriptions of ideal magnetohydrodynamics allow one to cast the evolution equations into a hyperbolic form useful for evolving rotating or binary compact objects with magnetic fields in numerical general relativity. Conserved circulation laws, such as those of Kelvin, Alfv\'en and Bekenstein-Oron, emerge simply as special cases of the Poincar\'e-Cartan integral invariant of Hamiltonian systems. We use this approach to obtain an extension of Kelvin's theorem to baroclinic (non-isentropic) fluids, based on a temperature-dependent time parameter. We further extend this result to perfectly or poorly conducting baroclinic magnetoflows. Finally, in the barotropic case, such magnetoflows are shown to also be geodesic, albeit in a Finsler (rather than Riemann) space.

The dynamical system studied in previous papers of this series, namely a bound time-like geodesic motion in the exact static and axially symmetric space-time of an (originally) Schwarzschild black hole surrounded by a thin disc or ring, is considered to test whether the often employed "pseudo-Newtonian" approach (resorting to Newtonian dynamics in gravitational potentials modified to mimic the black-hole field) can reproduce phase-space properties observed in the relativistic treatment. By plotting Poincar\'e surfaces of section and using two recurrence methods for similar situations as in the relativistic case, we find similar tendencies in the evolution of the phase portrait with parameters (mainly with mass of the disc/ring and with energy of the orbiters), namely those characteristic to weakly non-integrable systems. More specifically, this is true for the Paczy\'nski--Wiita and a newly suggested logarithmic potential, whereas the Nowak--Wagoner potential leads to a different picture. The potentials and the exact relativistic system clearly differ in delimitation of the phase-space domain accessible to a given set of particles, though this mainly affects the chaotic sea whereas not so much the occurrence and succession of discrete dynamical features (resonances). In the pseudo-Newtonian systems, the particular dynamical features generally occur for slightly smaller values of the perturbation parameters than in the relativistic system, so one may say that the pseudo-Newtonian systems are slightly more prone to instability. We also add remarks on numerics (a different code is used than in previous papers), on the resemblance of dependence of the dynamics on perturbing mass and on orbital energy, on the difference between the Newtonian and relativistic Bach--Weyl rings, and on the relation between Poincar\'e sections and orbital shapes within the meridional plane.