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.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

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.

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.

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.

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.

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.

This is a review with the ambitious goal of covering the recent progress in: 1) universal relations (in general relativity and alternative theories of gravity), and 2) neutron star models in alternative theories. We also aim to be complementary to recent reviews in the literature that are covering the topic of universal relations of neutron stars.

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.

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.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

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 in detail the f-mode secular instability for rapidly rotating neutron stars, putting emphasis on supermassive models which do not have a stable nonrotating counterpart. Such neutron stars are thought to be the generic outcome of the merger of two standard mass neutron stars. In addition we take into account the effects of strong magnetic field and r-mode instability, that can drain a substantial amount of angular momentum. We find that the gravitational wave signal emitted by supramassive neutron stars can reach above the Advance LIGO sensitivity at distance of about 20Mpc and the detectability is substantially enhanced for the Einstein Telescope. The event rate will be of the same order as the merging rates, while the analysis of the signal will carry information for the equation of state of the post-merging neutron stars and the strength of the magnetic fields.

In the present paper we examine gravitational wave asteroseismology relations for f-modes of rapidly rotating neutron stars. An approach different to the previous studies is employed - first, the moment of inertia is used instead of the stellar radius, and second, the normalization of the oscillation frequencies and damping times is different. It was shown that in the non-rotating case this can lead to a much stronger equation of state independence and our goal is to generalize the static relations to the rapidly rotating case and values of the spherical mode number $l\ge2$. We employ realistic equations of state that cover a very large range of stiffness in order to check better the universality of the relations. At the end we explore the inverse problem, i.e. obtain the neutron star parameters from the observed gravitational frequencies and damping times. It turns out that with this new set of relations we can solve the inverse problem with a very good accuracy using three frequencies that was not possible in the previous studies where one needs also the damping times. The asteroseismology relations are also particularly good for the massive rapidly rotating models that are subject to secular instabilities.

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.

$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.

In the present paper we study the oscillations of fast rotating neutron stars with realistic equations of state (EoS) within the Cowling approximation. We derive improved empirical relations for gravitational wave asteroseismology with f-modes and for the first time we consider not only quadrupolar oscillations but also modes with higher spherical order (l=|m|=3,4). After performing a systematic comparison with polytropic EoS, it is shown that the empirical relations found in this case approximately also hold for realistic EoS. Even more, we show that these relations will not change significantly even if the Cowling approximation is dropped and the full general relativistic case is considered, although the normalization used here(frequencies and damping times in the nonrotating limit) could differ considerably. We also address the inverse problem, i.e. we investigate in detail what kind of observational data is required in order to determine characteristical neutron star parameters. It is shown that masses, radii and rotation rates can be estimated quite accurately using the derived asteroseismology relations. We also compute the instability window for certain models, i.e. the limiting curve in a T-\Omega-plane where the secular Chandrasekhar-Friedman-Schutz(CFS) instability overcomes dissipative effects, and show that some of the modern realistic EoS will lead to a larger instability window compared to all of the polytropic ones presented so far in the literature. Additionally, we calculate the r-mode instability window and compare it with the f-mode-case. The overall results for the instability window suggest that it is vital to take into account oscillations with l=3,4 when considering gravitational wave asteroseismology using the f-mode in rapidly rotating neutron stars, as these modes can become CFS unstable for a much larger range of parameters than pure quadrupolar oscillations.

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.

The imminent detection of gravitational waves will trigger precision tests of gravity through observations of quasinormal ringing of black holes. While General Relativity predicts just two polarizations of gravitational waves, the so-called plus and cross polarizations, numerous alternative theories of gravity predict up to six different polarizations which will potentially be observed in current and future generations of gravitational wave detectors. Bekenstein's Tensor-Vector-Scalar (TeVeS) theory and its generalization fall into one such class of theory that predict the full gamut of six polarizations of gravitational waves. In this paper we begin the study of quasinormal modes (QNMs) in TeVeS by studying perturbations of the scalar field in a spherically symmetric background. We show that, at least in the case where superluminal propagation of perturbations is not present, black holes are generically stable to this kind of perturbation. We also make a unique prediction that, as the limit of the various coupling parameters of the theory tend to zero, the QNM spectrum tends to $1/\sqrt{2}$ times the QNM spectrum induced by scalar perturbations of a Schwarzschild black hole in General Relativity due to the intrinsic presence of the background vector field. We further show that the QNM spectrum does not vary significantly from this value for small values of the theory's coupling parameters, however can vary by as much as a few percent for larger, but still physically relevant parameters.