In the present paper we consider multi-scalar extension of Einstein-Gauss-Bonnet gravity. We focus on multi-scalar Einstein-Gauss-Bonnet models whose target space is a three-dimensional maximally symmetric space, namely either $\mathbb{S}^3$, $\mathbb{H}^3$ or $\mathbb{R}^3$, and in the case when the map $\text{\it spacetime} \to \text{\it target space}$ is nontrivial. We prove numerically the existence of black holes in this class of models for several Gauss-Bonnet coupling functions, including the case of scalarization. We also perform systematic study of a variety of black hole characteristics and the space-time around them, such as the area of the horizon, the entropy and the radius of the photon sphere. One of the most important properties of the obtained solutions is that the scalar charge is zero and thus the scalar dipole radiation is suppressed which leads to much weaker observational constraints compared to the majority of modified theories possessing a scalar degree of freedom. For one of the coupling functions we could find branches of scalarized black holes which have a nontrivial structure -- there is non-uniqueness of the scalarized solutions belonging to a single branch and there is a region of the parameter space where most probably stable scalarized black holes coexist with the stable Schwarzschild black holes. Such a phenomena can have a clear observational signature.

In the present paper we consider the extended scalar-tensor-Gauss-Bonnet gravity with a massive scalar field. We prove numerically the existence of Gauss-Bonnet black holes for three different forms of the coupling function including the case of spontaneous scalarization. We have performed a systematic study of the black hole characteristics such as the area of the horizon, the entropy and the temperature for these coupling functions and compared them to the Schwarzschild solutions. The introduction of scalar field mass leads to a suppression of the scalar field and the increase of this mass brings the black holes closer to the Schwarzschild case. For linear and exponential coupling, a nonzero scalar field mass expands the domain of existence of black holes solutions. Larger deviations from the GR case are observed only for small masses and these differences decrease with the increase of the scalar field mass. In the case of a coupling function which leads to scalarization the scalar field mass has a significant influence on the bifurcation points where the scalarized black holes branch out of the Schwarzschild solution. The largest deviation from the case with a massless scalar field are observed for black hole masses close to the bifurcation point.

Testing modified theories of gravity with direct observations of the parameters of a neutron star is not the optimal way of testing gravitational theories. However, observing electromagnetic signals originating from the close vicinity of the compact object my turn out an excellent way of probing spacetime in strong field regime. A promising candidate for doing so are the so-called quasi-periodic oscillations, observed in the X-ray light curves of some pulsars. Although the origin of those oscillations is unknown, one thing most of the models describing them have in common is that in one way or another they incorporate the radius of the innermost stable circular obit, and the orbital and the epicyclic frequencies of particles moving around the compact object. In this paper we study the aforementioned quantities in the context of massive scalar-tensor theory and massive scalar-tensor theory with self-interaction, both of which in strong regime allow for significant deviations from General relativity for values for the free parameters of the theory in correlation with the observations.

We are investigating universal relations between different normalisations of the moment of inertia and the compactness of neutron stars in slow rotation approximation. We study the relations in particular class of massive scalar-sensor theories with self-interaction, for which significant deviations from General Relativity are allowed for values of the parameters that are in agreement with the observations. Moment of inertia-compactness relations are examined for different normalisation of the moment of inertia. It is shown that for all studied cases the deviations from EOS universality are small for the examined equations of state. On the other hand the scalarization can lead to large deviations from the general relativistic universal relations for values of the parameters that are in agreement with the current observations that can be potentially used to set further test the scalar-tensor theories.

Binary pulsar observations and gravitational wave detections seriously constrained scalar-tensor theories with massless scalar field allowing only small deviations from general relativity. If we consider a nonzero mass of the scalar field, though, significant deviations from general relativity are allowed for values of the parameters that are in agreement with the observations. In the present paper we extend this idea and we study scalar-tensor theory with massive field with self-interaction term in the potential. The additional term suppresses the scalar field in the neutron star models in addition to the effect of the mass of the scalar field but still, large deviations from pure GR can be observed for values of the parameters that are in agreement with the observations.

In the present paper the axial quasi-normal modes of neutron stars in $f(R)$ gravity are examined using a large set of equations of state. The numerical calculations are made using two different approaches -- performing time evolution of the perturbation equations and solving the time-independent representation of the equations as a boundary value problem. According to the results the mode frequencies and the damping times decrease with the increase of the free parameter of the theory in comparison to the pure general relativistic case. While the frequencies deviate significantly from Einstein's theory for all realistic neutron star masses (say above $1M_\odot$), the damping times reach non-negligible differences only for the more massive models. We have constructed as well universal (equation of state independent) gravitational wave asteroseismology relations involving the frequencies and the damping times. It turns out that the equation of state independence is preserved using the same normalization as in pure general relativity and the qualitative differences of the phenomenological relations with respect to Einstein's theory of gravity can be large for large values of the free parameter in $f(R)$ gravity.

The glitch activity of young pulsars arises from the exchange of angular momentum between the crust and the interior of the star. Recently, it was inferred that the moment of inertia of the crust of a neutron star is not sufficient to explain the observed glitches. Such estimates are presumed in the Einstein's general relativity in describing the hydrostatic equilibrium of neutron stars. The crust of the neutron star has a space-time curvature of 14 orders of magnitude larger than that probed in solar system tests. This makes gravity the weakest constrained physics input in the crust related processes. We calculate the ratio of the crustal to the total moment of inertia of neutron stars in the scalar-tensor theory of gravity and the non-perturbative $f({\cal R})={\cal R}+ a {\cal R}^2$ gravity. We find for the former that the crust to core ratio of the moment of inertia does not change significantly from what is inferred in general relativity. For the latter we find that the ratio increases significantly from what is inferred in general relativity in the case of high mass objects. Our results suggest that the glitch activity of pulsars may be used to probe gravity models, although the gravity models explored in this work are not appropriate candidates.

We study the electromagnetic spectrum of accretion disks around neutron and strange stars in $\mathcal{R}^2$ gravity. Both static and rapidly rotating models are investigated. The results are compared with the General Relativistic results. We found difference between the results in both theories of about 50\% for the electromagnetic flux and about 20\% in the luminosity for models with equal mass and angular velocity in both theories. The observed differences are much lower for models rotating with Kelperian velocity and with equal masses.

We are investigating universal relations between different normalisations of the moment of inertia and the compactness of neutron and strange stars. Slowly rotating as well as rapidly rotating models are studied in General Relativity, $\mathcal{R}^2$ gravity and scalar--tensor theories of gravity. Moment of inertia -- compactness relations are examined for different normalisations of the moment of inertia. It is shown that for all studied cases the deviations from EOS universality are small for the examined equations of state. It turns out that in some of the cases the examined relations are also theory independent to a good extent. Universality in relations between the maximum mass and the moment of inertia for some unstable models is also investigated.

According to various models, the orbital and the epicyclic frequencies of particles moving on a circular orbit around compact objects are related to the quasi-periodic oscillations observed in the X-ray flux of some pulsars or black hole candidates. It is expected that they originate from the inner edge of the accretion discs, deep into the gravitational field of the compact objects. Considering the planned new generation X-ray timing observatories with large collective areas, the quasi-periodic oscillations might be an excellent tool for testing gravity in strong field regime and respectively alternative gravitational theories. We examine the orbital and the epicyclic frequencies of a particle moving on a circular orbit around neutron or strange stars in $R^2$ gravity. The case of slow rotation is considered too. The $R^2$ gravity results are compared to the General Relativistic case. We comment the deviations from General Relativity, as well as the deviations due to rotation in both theories.

We study how rapid rotation influences the relation between the normalized moment of inertia $\bar{I}$ and quadrupole moment $\bar{Q}$ for scalarized neutron stars. The questions one has to answer are whether the EOS universality is preserved in this regime and what are the deviations from general relativity. Our results show that the $\bar{I}-\bar{Q}$ relation is nearly EOS independent for scalarized rapidly rotating stars, but the differences with pure Einstein's theory increase compared to the slowly rotating case. In general, smaller negative values of the scalar field coupling parameters $\beta$ lead to larger deviations, but these deviations are below the expected accuracy of the future astrophysical observations if one considers values of $\beta$ in agreement with the current observational constraint. An important remark is that although the normalized $\bar{I}-\bar{Q}$ relation is quite similar for scalar-tensor theories and general relativity, the unnormalized moment of inertia and quadrupole moment can be very different in the two theories. This demonstrates that although the $\bar{I}-\bar{Q}$ relations are potentially very useful for some purposes, they might not serve us well when trying to distinguish between different theories of gravity.

In the present paper we investigate self-consistently slowly rotating neutron and strange stars in R-squared gravity. For this purpose we first derive the equations describing the structure of the slowly rotating compact stars in $f(R)$-gravity and then simultaneously solve the exterior and the interior problem. The structure of the slowly rotating neutron stars is studied for two different hadronic equations of state and a strange matter equation of state. The moment of inertia and its dependence on the stellar mass and the $R$-squared gravity parameter $a$ is also examined in details. We find that the neutron star moment of inertia for large values of the parameter $a$ can be up to $30\%$ larger compared to the corresponding general relativistic models. This is much higher than the change in the maximum mass induced by $R$-squared gravity and is beyond the EOS uncertainty. In this way the future observations of the moment of inertia of compact stars could allow us to distinguish between general relativity and $f(R)$ gravity, and more generally to test the strong field regime of gravity.

In the present paper we investigate non-perturbatively and self-consistently the structure of neutron stars in $R$-squared gravity by simultaneously solving the interior and exterior problem. The mass-radius relations are obtained for several equations of state and for wide range of the $R$-squared gravity parameter $a$. Even though the deviation from general relativity for nonzero values of $a$ can be large, they are still comparable with the variations due to different modern realistic equations of state. That is why the current observations of the neutron star masses and radii alone can not put constraints on the value of the parameter $a$. We also compare our results with those obtained within the perturbative method and we discuss the differences between them.