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.

The accuracy of the numerical findings of Lukes-Gerakopoulos PRD 86, 044013 (2012), regarding the existence of additional first integrals in the Zipoy-Voorhees space-time, was recently questioned by Maciejewski et al. PRD 88, 064003 (2013). In this comment, it is shown that the discrepancy between the results of Lukes-Gerakopoulos (2012) and Maciejewski et al. (2013) is not due to issues related to numerical accuracy, as claimed in Maciejewski et al (2013), but due to different choice of coordinates used in Maciejewski et al.(2013).

In this article we discuss a possible way of testing the Kerr black hole hypothesis by taking advantage of phenomena correlated with chaotic motion in the final stages of an accretion disk around a bumpy black hole. We anticipate that these phenomena should have an imprint in the electromagnetic spectrum coming from the accretion disk.

Black hole candidates in X-ray binary systems and at the centers of galaxies are expected to be the Kerr black holes of General Relativity, but the actual nature of these objects has yet to be verified. In this paper, we consider the possibility that they are exotic compact objects and we describe their exterior gravitational field with a subclass of the Manko-Novikov metrics, which are exact solutions of the vacuum Einstein's equations and can describe the spacetime geometry around bodies with arbitrary mass-multipole moments. We point out that around a Manko-Novikov object there may exist many disconnected non-plunging regions at small radii, with no counterpart in the Kerr background, and that their existence may be tested. For instance, in the presence of an accretion disk, they may be filled by the accreting gas, forming a ring structure that might remind the rings of Saturn. We suggest that the existence of these regions may have a clear observational signature in the waveform of the gravitational radiation emitted by an EMRI: in the last stage of the inspiral, the waveform would be the combination of "regular chirps", produced when the small object orbits in one of the non-plunging regions, and "bursts", released when the small object jumps from a non-plunging region to another one at smaller radii. Our conclusions are supported by some numerical calculations of trajectories in the geodesic approximation, in which a particle plunges from the ISCO and then seems to get trapped in the potential well at smaller radii.

We study the dynamics of the scalar field FLRW flat cosmological models within the framework of the Unified Dark Matter (UDM) scenario. In this model we find that the main cosmological functions such as the scale factor of the Universe, the scalar field, the Hubble flow and the equation of state parameter are defined in terms of hyperbolic functions. These analytical solutions can accommodate an accelerated expansion, equivalent to either the dark energy or the standard $\Lambda$ models. Performing a joint likelihood analysis of the recent supernovae type Ia data and the Baryonic Acoustic Oscillations traced by the SDSS galaxies, we place tight constraints on the main cosmological parameters of the UDM cosmological scenario. Finally, we compare the UDM scenario with various dark energy models namely $\Lambda$ cosmology, parametric dark energy model and variable Chaplygin gas. We find that the UDM scalar field model provides a large and small scale dynamics which are in fair agreement with the predictions by the above dark energy models although there are some differences especially at high redshifts.

We study the dynamics of the closed scalar field FRW cosmological models in the framework of the so called Unified Dark Matter (UDM) scenario. Performing a theoretical as well as a numerical analysis we find that there is a strong indication of chaos in agreement with previous studies. We find that a positive value of the spatial curvature is essential for the appearance of chaoticity, though the Lyapunov number seems to be independent of the curvature value. Models that are close to flat exhibit a chaotic behavior after a long time while pure flat models do not exhibit any chaos. Moreover, we find that some of the semiflat models in the UDM scenario exhibit similar dynamical behavior with the Lambda cosmology despite their chaoticity. Finally, we compare the measured evolution of the Hubble parameter derived from the differential ages of passively evolving galaxies with that expected in the semiflat unified scalar field cosmology. Based on a specific set of initial conditions we find that the UDM scalar field model matches well the observational data.