In the modern era of abundant X-ray detections and the increasing momentum of gravitational waves astronomy, tests of General Relativity in strong field regime become increasingly feasible and their importance for probing gravity cannot be understated. To this end, we study the characteristics of slowly rotating topological neutron stars in the tensor-multi-scalar theories of gravity following the static study of this new type of compact objects by two of the authors. We explore the moment of inertia and verify that universal relations known from General Relativity hold for this new class of compact objects. Furthermore, we study the properties of their innermost stable circular orbits (ISCO) and the epicyclic frequencies due to the latter's hinted link to observational quantities such as quasi-periodic X-ray spectrum features.

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.

Tensor-multi-scalar theories (TMST) are among the most natural generalizations of Einstein's theory, they are mathematically self-consistent and free from pathologies. They pass through all the known observations but contrary to standard scalar-tensor theories, TMST offer extremely rich spectrum of solutions and allow for large deviations from general relativity. One of the most interesting objects in these theories are the topological neutron stars. The main goal in the present paper is to study their radial stability since many branches of solutions can exist even for a fixed value of the topological charge. It turns out that only one of these branches is stable for each value of the topological charge while all the rest are unstable. The stable branch is exactly the one that spans from small to large neutron star masses having as well moderate values of the radii in agreement with the observations. As expected, it becomes unstable at the maximum of the mass. The frequencies of the radial modes are examined and it turns out that in most cases the mode frequencies are larger than the general relativistic ones and they increase with the increase of the topological charge. All this shows that the topological neutron stars are viable astrophysical objects that should be explored further in order to determine their observational manifestations.

In the present paper, we numerically construct new non-topological, spontaneously scalarized neutron stars in the tensor-multi-scalar theories of gravity 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 of a nontrivial map $\varphi : \text{\it spacetime} \to \text{\it target space}$. The theories of gravity admitting scalarization are characterized by the fact that the field equations always admit the general relativistic solution but for certain ranges of the parameters space it loses stability and nonlinear development of a scalar field is observed. Thus, in order to determine the values of the parameters where such scalarization is possible we studied the stability of the general relativistic solution within the framework of the considered tensor-multi-scalar theories. Based on these results we could obtain a family of scalarized branches characterized by the number of the scalar field nodes. These branches bifurcate from the general relativistic solution at the points where new unstable modes appear and they are energetically more favorable over the pure Einstein solutions. Interestingly, in certain parameter ranges, we could obtain non-uniqueness within a single branch of scalarized solutions.

In the present paper we investigate the structure of relativistic stars in 4D Einstein-Gauss-Bonnet gravity. The mass-radius relations are obtained for realistic hadronic and for strange quark star equations of state, and for a wide range of the Gauss-Bonnet coupling parameter $\alpha$. Even though the deviations from general relativity for nonzero values of $\alpha$ can be large enough, they are still comparable with the variations due to different modern realistic equations of state if we restrict ourselves to moderate value of $\alpha$. That is why the current observations of the neutron star masses and radii alone can not impose stringent constraints on the value of the parameter $\alpha$. Nevertheless some rough constraints on $\alpha$ can be put. The existence of stable stellar mass black holes imposes $\sqrt{\alpha}\lesssim 2.6 {\rm km}$ for $\alpha>0$ while the requirement that the maximum neutron star mass should be greater than two solar masses gives $\sqrt{|\alpha|}\lesssim 3.9 {\rm km}$ for $\alpha<0$. We also present an exact solution describing the structure of relativistic stars with uniform energy density in 4D Einstein-Gauss-Bonnet gravity.

In the present paper we demonstrate that there exists a new type of neutron stars in the tensor-multi-scalar theories of gravity. We call this new type of neutron stars topological neutron stars. In addition to the standard characteristics of the usual neutron stars the topological neutron stars are also characterized by a topological charge. We numerical construct explicit examples of topological neutron stars in tensor-multi-scalar theories whose target space is $\mathbb{S}^3$. Besides the topological charge the topological neutron stars also exhibit other attractive features which can place them among the realistic compact objects that could exist in Nature. For example they possess zero scalar charge and thus evades the strong binary pulsar constraints on the dipole scalar radiation.

In the present paper we consider special classes of massive tensor-multi-scalar theories of gravity whose target space metric admits Killing field(s) with a periodic flow. For such tensor-multi-scalar theories we show that there exist mixed configurations of tenor-multi-scalar solitons and relativistic (neutron) stars. The influence of the curvature of the target space on the structure of the mixed configurations is studied for two explicit scalar-tensor theories. The stability of the obtained compact objects is examined and it turns out that the stability region spans to larger central energy densities and central values of the scalar field compared to the pure neutron stars or the pure tensor-multi-scalar solitons.

We study the optical appearance and the apparent radiation flux of a thin accretion disk around the static Janis-Newman-Winicour naked singularity. We confine ourselves to the astrophysically most relevant case, when the solution possesses a photon sphere, assuming that the radiation emitted by the disk is described by the Novikov-Thorne model. The observable images resemble closely the visual appearance of the Schwarzschild black hole, as only quantitative differences are present. For the Janis-Newman-Winicour solution the accretion disk appears smaller, and its emission is characterized by a higher peak of the radiation flux. In addition, the most significant part of the radiation is concentrated in a closer neighbourhood of the flux maximum. The results are obtained independently by two alternative methods, consisting of a semi-analytical scheme using the spherical symmetry of the spacetime, and a fully numerical ray-tracing procedure valid for any stationary and axisymmetric spacetime.

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.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

In the present paper we consider special classes of tensor-multi-scalar theories of gravity whose target-space metric admits Killing field(s) with a periodic flow. For such tensor-multi-scalar theories we show that if the dynamics of the scalar fields is confined on the periodic orbits of the Killing field(s) then there exists a new type of compact objects - the tensor-multi-scalar solitons formed by a condensation of the gravitational scalars. The existence of the tensor-multi-scalar solitons is proven by solving the fully non-linear eigenvalue problem which follows from the dimensional reduction of the field equations of tensor-multi-scalar theories of gravity. The mass of the tensor-multi-scalar solitons can range at least from the mass of a neutron star to the mass of dark objects in the center of the galaxies in dependence of mass(es) of the gravitational scalars and which sector of massive scalars is excited. These facts show that the tensor-multi-scalar solitons could have important implications for the dark matter problem. The existence of the tensor-multi-scalar solitons points towards the possibility that the dark matter, or part of it, is made of condensed gravitational scalars.

In recent years there has been a surge of interest in what has come to be known as the `universal relations' between various global properties of neutron stars. These universal relations are equation of state independent relations between quantities such as the moment of inertia $I$, the tidal deformability or Love number $\lambda$, and the quadrupole $Q$ (I-Love-Q relations), or the relativistic multipole moments (3-hair relations). While I-Love-Q relations have been studied extensively in both general relativity and various alternatives, 3-hair relations have been studied only in general relativity. Recent progress on the definition of the multipole moments of a compact object in the case of scalar-tensor theories allows for the study of 3-hair relations in modified theories of gravity. Specifically, the aim of this work is to study them for scalarized stars in scalar-tensor theories with a massless scalar field that admit spontaneous scalarization. We find that the 3-hair relations between the mass and angular momentum moments that hold in general relativity hold for scalarized stars as well. The scalar moments also exhibit a universal behaviour, which is equation of state independent within one specific theory, but differs between different theories. Combining astrophysical observations one can in principle measure the different properties of scalarized neutron star and tell different theories apart.

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.

Recently new scalarized black hole solutions were constructed in the extended scalar-tensor-Gauss-Bonnet gravity, where the scalar field is sourced by the curvature of the spacetime via the Gauss-Bonnet invariant. A natural extension of these results is to consider the case of nonzero black hole charge. In addition we have explored a large set of coupling functions between the Gauss-Bonnet invariant and the scalar field, that was not done until now even in the uncharged case, in order to understand better the behavior of the solutions and the deviations from pure general relativity. The results show that in the case of nonzero black hole charge two bifurcation points can exist - one at larger masses where the scalarized solutions bifurcated from the Reissner-Nordstrom one, and one at smaller masses where the scalar charge of the solutions decreases again to zero and the branch merges again with the GR one. All of the constructed scalarized branches do not reach an extremal limit. We have examined the entropy of the black holes with nontrivial scalar field and it turns out, that similar to the uncharged case, the fundamental branch which possesses scalar field with no nodes is thermodynamically favorable over the Reissner-Nordstrom one for the considered coupling functions, while the rest of the branches possessing scalar field with one or more zeros have lower entropy compared to the GR case and they are supposed to be unstable.

We present the first numerical models of differentially rotating stars in alternative theories of gravity. We chose a particular class of scalar-tensor theories of gravity that is indistinguishable from GR in the weak field regime but can lead to significant deviations when strong fields are considered. We show that the maximum mass that a differentially rotating neutron star can sustain increases significantly for scalarized solutions and such stars can reach larger angular momenta. In addition, the presence of a nontrivial scalar field has the effect of increasing the required axis ratio for reaching a given value of angular momentum, when compared to a corresponding model of same rest mass in general relativity. We find that the scalar field also makes rapidly rotating models less quasi-toroidal than their general-relativistic counterparts. For large values of the angular momentum and values of the coupling parameter that are in agreement with the observations, we find a second turning point for scalarized models along constant angular momentum sequences, which could have interesting implications for the stability of remnants created in a binary neutron star merger.

We revisit the shadow of rotating traversable wormholes discussing the role of the wormhole throat in the shadow formation. For certain classes of wormholes the throat serves as a potential barrier for light rays with particular impact parameters, thus modifying the shadow shape. We consider a couple of wormhole solutions and examine the structure of their shadow images, and the intrinsic mechanisms for their formation. Some of the shadows possess cuspy edges, which arise due to the interplay of two distinct families of unstable spherical orbits. These solutions provide examples, in which the explicit mechanism for cusp formation can be uncovered.

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 recent detection of gravitational waves from a neutron star merger was a significant step towards constraining the nuclear matter equation of state by using the tidal Love numbers (TLNs) of the merging neutron stars. Measuring or constraining the neutron star TLNs allows us in principle to exclude or constraint many equations of state. This approach, however, has the drawback that many modified theories of gravity could produce deviations from General Relativity similar to the deviations coming from the uncertainties in the equation of state. The first and the most natural step in resolving the mentioned problem is to quantify the effects on the TLNs from the modifications of General Relativity. With this motivation in mind, in the present paper we calculate the TLNs of (non-rotating) neutron stars in $f(R)$ gravity. For this purpose, we first derived the equations describing both the polar and the axial stationary perturbations of neutron stars in a particular class of $f(R)$ gravity, the so-called $R^2$-gravity. Then, by solving numerically the perturbation equations, we calculate explicitly the polar and the axial $l=2$ TLNs of the neutron stars in $R^2$-gravity for three characteristic realistic equations of state. Our results show that while the polar TLNs are slightly influenced by the $R^2$ modification of General Relativity, the axial TLNs can be several times larger (in terms of the absolute value) compared to the general relativistic case.

We study the local shadow of the Schwarzschild black hole with a quadrupole distortion and the influence of the external gravitational field on the photon dynamics. The external matter sources modify the light ring structure and lead to the appearance of multiple shadow images. In the case of negative quadrupole moments we identify the most prominent mechanism causing multiple shadow formation. Furthermore, we obtain a condition under which this mechanism can be realized. This condition depends on the quadrupole moment, but also on the position of the observer and the celestial sphere.

In the present paper we study models of neutron stars in a class of extended scalar-tensor Gauss-Bonnet (ESTGB) theories for which the scalar degree of freedom is exited only in the strong curvature regime. We show that in the framework of the ESTGB theories under consideration there exist new neutron star solutions which are formed via spontaneous scalarization of the general relativistic neutron stars. In contrast to the spontaneous scalarization in the standard scalar-tensor theories which is induced by the presence of matter, in our case the scalarization is induced by the spacetime curvature.

In the present paper we consider a class of extended scalar-tensor-Gauss-Bonnet (ESTGB) theories for which the scalar degree of freedom is excited only in the extreme curvature regime. We show that in the mentioned class of ESTGB theories there exist new black hole solutions which are formed by spontaneous scalarization of the Schwarzaschild balck holes in the extreme curvature regime. In this regime, below certain mass, the Schwarzschild solution becomes unstable and new branch of solutions with nontrivial scalar field bifurcate from the Schwarzschild one. As a matter of fact, more than one branches with nontrivial scalar field can bifurcate at different masses but only the first one is supposed to be stable. This effect is quite similar to the spontaneous scalarization of neutron stars. In contrast with the standard spontaneous scalarization of neutron stars which is induced by the presence of matter, in our case the scalarization is induced by the curvature of the spacetime.

We perform the first study of the oscillation frequencies of rapidly rotating neutron stars in alternative theories of gravity, focusing mainly on the fundamental $f$-modes. We concentrated on a particular class of alternative theories - the (massive) scalar-tensor theories. The generalization to rapid rotation is important because on one hand the rapid rotation can magnify the deviations from general relativity compared to the static case and on the other hand some of the most efficient emitters of gravitational radiation, such as the binary neutron star merger remnants, are supposed to be rotating close to their Kepler (mass-shedding) limits shortly after their formation. We have constructed several sequences of models starting from the nonrotating case and reaching up to the Kepler limit, with different values of the scalar-tensor theory coupling constant and the scalar field mass. The results show that the deviations from pure Einstein's theory can be significant especially in the case of nonzero scalar field mass. An important property of the oscillation modes of rapidly rotating stars is that they can become secularly unstable due to the emission of gravitational radiation, that is so-called Chandrasekhar-Friedman-Schutz instability. Such unstable modes are efficient emitters of gravitational radiation. Our studies show that the inclusion of nonzero scalar field would decrease the threshold value of the normalized angular momentum where this instability stars to operate, but the growth time of the instability seems to be increased compared to pure general relativity.

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 construct rapidly rotating neutron star models in scalar-tensor theories with a massive scalar field. The fact that the scalar field has nonzero mass leads to very interesting results since the allowed range of values of the coupling parameters is significantly broadened. These deviations from pure general relativity can be very large for values of the parameters that are in agreement with the observations. The rapid rotation can magnify the differences several times compared to the static case. The universal relations between the normalized moment of inertia and quadrupole moment are also investigated both for the slowly and rapidly rotating cases. The results show that these relations are still EOS independent up to a large extend and the deviations from pure general relativity can be large. This places the massive scalar-tensor theories amongst the few alternative theories of gravity that can be tested via the universal $I$-Love-$Q$ relations.

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.

In the scalar-tensor theories with a massive scalar field the coupling constants, and the coupling functions in general, which are observationally allowed, can differ significantly from those in the massless case. This fact naturally implies that the scalar-tensor neutron stars with a massive scalar field can have rather different structure and properties in comparison with their counterparts in the massless case and in general relativity. In the present paper we study slowly rotating neutron stars in scalar-tensor theories with a massive gravitational scalar. Two examples of scalar-tensor theories are examined - the first example is the massive Brans-Dicke theory and the second one is a massive scalar-tensor theory indistinguishable from general relativity in the weak field limit. In the later case we study the effect of the scalar field mass on the spontaneous scalarization of neutron stars. Our numerical results show that the inclusion of a mass term for the scalar field indeed changes the picture drastically compared to the massless case. It turns out that mass, radius and moment of inertia for neutron stars in massive scalar-tensor theories can differ drastically from the pure general relativistic solutions if sufficiently large masses of the scalar field are considered.

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 consider the problem for the classification of static and asymptotically flat Einstein-Maxwell-dilaton spacetimes with a photon sphere. It is first proven that the photon spheres in Einstein-Maxwell-dilaton gravity have constant mean and constant scalar curvature. Then we derive some relations between the mean curvature and the physical characteristics of the photon spheres. Using further the symmetries of the dimensionally reduced Einstein-Maxwell-dilaton field equations we show that the lapse function, the electrostatic potential and the dilaton field are functionally dependent in the presence of a photon sphere. Using all this we prove the main classification theorem by explicitly constructing all Einstein-Maxwell-dilaton solutions possessing a non-extremal photon sphere.

In the present paper we study the behavior of the normalized $I$-$Q$ relation for neutron stars in a particular class of $f(R)$ theories of gravity, namely the $R^2$ gravity that is one of the most natural and simplest extensions of general relativity in the strong field regime. We study both the slowly and rapidly rotating cases. The results show that the $I$-$Q$ relation remain nearly equation of state independent for fixed values of the normalized rotational parameter, but the deviations from universality can be a little bit larger compared to the general relativistic case. What is the most interesting in our studies, is that the differences with the pure Einstein's theory can be large reaching above 20\%. This is qualitative different from the majority of alternative theories of gravity, where the normalized $I$-$Q$ relations are almost indistinguishable from the general relativistic case, and can lead to observational constraints on the $f(R)$ theories in the future.

In the present paper we prove a uniqueness theorem for the static and asymptotically flat solutions to the Einstein-scalar field equations which possess a photon sphere. We show that such solutions are uniquely specified by their mass $M$ and scalar charge $q$ and that they are isometric to the Janis-Newman-Winicour solution with the same mass and scalar charge subject to the inequality $\frac{q^2}{M^2}<3$.

$f(R)$ theories of gravity are one of the most popular alternative explanations for dark energy and therefore studying the possible astrophysical implications of these theories is an important task. In the present paper we make a substantial advance in this direction by considering rapidly rotating neutron stars in $R^2$ gravity. The results are obtained numerically and the method we use is non-perturbative and self-consistent. The neutron star properties, such as mass, radius and moment of inertia, are studied in detail and the results show that rotation magnifies the deviations from general relativity and the maximum mass and moment of inertia can reach very high values. This observation is similar to previous studies of rapidly rotating neutron stars in other alternative theories of gravity, such as the scalar-tensor theories, and it can potentially lead to strong astrophysical manifestations.

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.

We study the orbital and epicyclic frequencies of particles orbiting around rapidly rotating neutron stars and strange stars in a particular scalar-tensor theory of gravity. We find very large deviations of these frequencies, when compared to their corresponding values in general relativity, for the maximum-mass rotating models. In contrast, for models rotating with spin frequency of 700Hz (approximately the largest known rotation rate of neutron stars), the deviations are generally small. Nevertheless, for a very stiff equation of state and a high mass the deviation of one of the epicyclic frequencies from its GR value is appreciable even at a spin frequency of 700Hz. In principle, such a deviation could become important in models of quasi-periodic oscillations in low-mass x-ray binaries and could serve as a test of strong gravity (if other parameters are well constraint). Even though the present paper is concentrated mainly on orbital and epicyclic frequencies, we present here for the first time rapidly rotating, scalarized equilibrium compact stars with realistic hadronic equations of state and strange matter equation of state. We also provide analytical expressions for the exterior spacetime of scalarized neutron stars and their epicyclic frequencies in the nonrotating limit.

We present the field equations governing the equilibrium of rapidly rotating neutron stars in scalar-tensor theories of gravity, as well as representative numerical solutions. The conditions for the presence of a nontrivial scalar field and the deviations from the general relativistic solutions are studied. Two examples of scalar-tensor theories are examined - one case that is equivalent to the Brans-Dicke theory and a second case, that is perturbatively equivalent to Einstein's General Relativity in the weak field regime, but can differ significantly for strong fields. Our numerical results show that rapidly rotating neutron star models with a nontrivial scalar field exist in both cases and that the effect of the scalar field is stronger for rapid rotation. If we consider values of the coupling parameters in accordance with current observations, only the second example of scalar-tensor theories has significant influence on the neutron star structure. We show that scalarized, rapidly rotating neutron stars exist for a larger range of the parameters than in the static case, since a nontrivial scalar field is present even for values of the coupling constant $\beta>-4.35$, and that these solutions are energetically more favorable than the general relativistic ones. In addition, the deviations of the rapidly rotating scalar-tensor neutron stars from the general-relativistic solutions is significantly larger than in the static case.

It was shown recently, that normalized relations between the moment of inertia (I), the quadrupole moment (Q) and the tidal deformability (Love number) exist and for slowly-rotating neutron stars they are almost independent of the equation of state (EOS). We extend the computation of the I-Q relation to models rotating up to the mass-shedding limit and show that the universality of the relations is lost. With increasing rotation rate, the normalized I-Q relation departs significantly from its slow-rotation limit deviating up to 40% for neutron stars and up to 75% for strange stars. The deviation is also EOS dependent and for a broad set of hadronic and strange matter EOS the spread due to rotation is comparable to the spread due to the EOS, if one considers sequences with fixed rotational frequency. Still, for a restricted sample of modern realistic EOS one can parameterize the deviations from universality as a function of rotation only. The previously proposed I-Love-Q relations should thus be used with care, because they lose their universality in astrophysical situations involving compact objects rotating faster than a few hundred Hz.

We explore the shadow of certain class of rotating traversable wormholes within classical general relativity. The images depend on the angular momentum of the wormhole, and the inclination angle of the observer. We compare the results with the case of the Kerr black hole. For small angular momenta the shadows for the two solutions are nearly identical, however with the increasing of the angular momentum they start to deviate considerably.

We model the spacetime of low-mass X-ray binaries with the Tomimatsu-Sato {\delta} = 2 (TS2) metric and study the properties of the orbital and the epicyclic frequencies. The numerical analysis shows that the properties of the characteristic frequencies of oscillation do not differ qualitatively from those of the Kerr black hole. Estimates for the angular momenta of the three stellar mass black hole candidates GRO 1655-40, XTE 1550-564 and GRS 1915+105 are made with the application of the nonlinear resonance model. We find agreement between the predictions based on the 3 : 2 nonlinear resonance model for a TS2 background and the current estimates found in the literature.

One of the most common assumptions in the study of neutron star models and their oscillations is that the pressure is isotopic, however there are arguments that this may not be correct. Thus in the present paper we make a first step towards studying the nonradial oscillations of neutron stars with an anisotropic pressure. We adopt the so-called Cowling approximation where the spacetime metric is kept fixed and the oscillation spectrum for the first few fluid modes is obtained. The effect of the anisotropy on the frequencies is apparent, although with the present results it might be hard to distinguish it from the changes in the frequencies caused by different equations of state.

In the present paper we study the oscillation spectrum of neutron stars containing both ordinary matter and dark energy in different proportions. Within the model we consider, the equilibrium configurations are numerically constructed and the results show that the properties of the mixed neuron-dark-energy star can differ significantly when the amount of dark energy in the stars is varied. The oscillations of the mixed neuron-dark-energy stars are studied in the Cowling approximation. As a result we find that the frequencies of the fundamental mode and the higher overtones are strongly affected by the dark energy content. This can be used in the future to detect the presence of dark energy in the neutron stars and to constrain the dark-energy models.

In the present paper we focus on building simple nonperturbative analytical relativistic models of magnetars. With this purpose in mind we first develop a method for generating exact interior solutions to the static and axisymmetric Einstein-Maxwell-hydrodynamic equations with anisotropic perfect fluid and with pure poloidal magnetic field. Then using an explicit exact solution we present a simple magnetar model and calculate some physically interesting quantities as the surface elipticity and the total energy of the magnetized star.

Adopting the phantom (ghost) scalar field description of dark energy, we construct a general class of exact interior solutions describing mixed relativistic stars containing both ordinary matter and dark energy in different proportions. The exterior solution that continuously matches the interior solutions is found. Exact solutions describing extremal configuration with zero ordinary matter pressure are also constructed.

The purpose of the current Letter is to give some relations between gravitational lensing in the strong-deflection limit and the frequencies of the quasinormal modes of spherically symmetric, asymptotically flat black holes. On the one side, the obtained relations can give a physical interpretation of the strong-deflection limit parameters. On the other side, they also give an alternative method for the measurement of the frequencies of the quasinormal modes of spherically symmetric, asymptotically flat black holes. They could be applied to the localization of the sources of gravitational waves and could tell us what frequencies of the gravitational waves we could expect from a black hole acting simultaneously as a gravitational lens and a source of gravitational waves.

In the present work we study numerically quasi-equatorial lensing by the charged, stationary, axially-symmetric Kerr-Sen dilaton-axion black hole in the strong deflection limit. In this approximation we compute the magnification and the positions of the relativistic images. The most outstanding effect is that the Kerr-Sen black hole caustics drift away from the optical axis and shift in clockwise direction with respect to the Kerr caustics. The intersections of the critical curves on the equatorial plane as a function of the black hole angular momentum are found, and it is shown that they decrease with the increase of the parameter $Q^{2}/M$. All of the lensing quantities are compared to particular cases as Schwarzschild, Kerr and Gibbons-Maeda black holes.

New nondiagonal $G_{2}$ inhomogeneous cosmological solutions are presented in a wide range of scalar-tensor theories with a stiff perfect fluid as a matter source. The solutions have no big-bang singularity or any other curvature singularities. The dilaton field and the fluid energy density are everywhere regular, too. The geodesic completeness of the solutions is investigated.

We construct exact static, axisymmetric solutions of Einstein-Maxwell-dilaton gravity presenting distorted charged dilaton black holes. The thermodynamics of such distorted black holes is also discussed.

We numerically construct static and spherically symmetric electrically charged black hole solutions in Einstein-Born-Infeld gravity with massive dilaton. The numerical solutions show that the dilaton potential allows many more black hole causal structures than the massless dilaton. We find that depending on the black hole mass and charge and the dilaton mass the black holes can have either one, two, or three horizons. The extremal solutions are also found out. As an interesting peculiarity we note that there are extremal black holes with an inner horizon and with triply degenerated horizon.

We investigate numerically models of the static spherically symmetric boson-fermion stars in the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of such stars is interpreted as a nonlinear two-parametric eigenvalue problem with unknown internal boundary. To solve this problem the Continuous Analogue of Newton Method is used.

We develop the general theory of stars in Saa's model of gravity with propagating torsion and study the basic stationary state of neutron star. Our numerical results show that the torsion force decreases the role of the gravity in the star configuration leading to significant changes in the neutron star masses depending on the equation of state of star matter. The inconsistency of the Saa's model with Roll-Krotkov-Dicke and Braginsky-Panov experiments is discussed.

It is shown that the recently proposed interpretation of the transposed equi-affine theory of gravity as a theory with variable Plank "constant" is inconsistent with basic solar system gravitational experiments.