We develop a new method to measure neutron star parameters and derive constraints on the equation of state of dense matter by fitting the frequencies of simultaneous Quasi Periodic Oscillation modes observed in the X-ray flux of accreting neutron stars in low mass X-ray binaries. To this aim we calculate the fundamental frequencies of geodesic motion around rotating neutron stars based on an accurate general-relativistic approximation for their external spacetime. Once the fundamental frequencies are related to the observed frequencies through a QPO model, they can be fit to the data to obtain estimates of the three parameters describing the spacetime, namely the neutron star mass, angular momentum and quadrupole moment. From these parameters we derive information on the neutron star structure and equation of state. We present a proof of principle of our method applied to pairs of kHz QPO frequencies observed from three systems (4U1608-52, 4U0614+09 and 4U1728-34). We identify the kHz QPOs with the azimuthal and the periastron precession frequencies of matter orbiting the neutron star, and via our Bayesian inference technique we derive constraints on the neutrons stars' masses and radii. This method is applicable to other geodesic-frequency-based QPO models.

An important question when studying the light curves produced by hot spots on the surface of rotating neutron stars is, how might the shape of the light curve be affected by relaxing some of the simplifying assumptions that one would adopt on a first treatment of the problem, such as that of a point-like spot. In this work we explore the dependence of light curves on the size and shape of a single hot spot on the surface of slowly-rotating neutron stars. More specifically, we consider two different shapes for the hot spots (circular and annular) and examine the resulting light curves as functions of the opening angle of the hot spot (for both) and width (for the latter). We find that the point-like approximation can describe light curves reasonably well, if the opening angle of the hot spot is less than $\sim 5^\circ$. Furthermore, we find that light curves from annular spots are almost the same as those of the full circular spot if the opening angle is less than $\sim 35^\circ$ independently of the hot spot's width.

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

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.

One of the central goals of the newborn field of gravitational wave astronomy is to test gravity in the highly nonlinear, strong field regime characterizing the spacetime of black holes. In particular, "black hole spectroscopy" (the observation and identification of black hole quasinormal mode frequencies in the gravitational wave signal) is expected to become one of the main tools for probing the structure and dynamics of Kerr black holes. In this paper we take a significant step towards that goal by constructing a "post-Kerr" quasinormal mode formalism. The formalism incorporates a parametrized but general perturbative deviation from the Kerr metric and exploits the well-established connection between the properties of the spacetime's circular null geodesics and the fundamental quasinormal mode to provide approximate, eikonal limit formulae for the modes' complex frequencies. The resulting algebraic toolkit can be used in waveform templates for ringing black holes with the purpose of measuring deviations from the Kerr metric. As a first illustrative application of our framework, we consider the Johannsen-Psaltis deformed Kerr metric and compute the resulting deviation in the quasinormal mode frequency relative to the known Kerr result.

The problem of having an accurate description of the spacetime around neutron stars is of great astrophysical interest. For astrophysical applications, one needs to have a metric that captures all the properties of the spacetime around a neutron star. Furthermore, an accurate appropriately parameterised metric, i.e., a metric that is given in terms of parameters that are directly related to the physical structure of the neutron star, could be used to solve the inverse problem, which is to infer the properties of the structure of a neutron star from astrophysical observations. In this work we present such an approximate stationary and axisymmetric metric for the exterior of neutron stars, which is constructed using the Ernst formalism and is parameterised by the relativistic multipole moments of the central object. This metric is given in terms of an expansion on the Weyl-Papapetrou coordinates with the multipole moments as free parameters and is shown to be extremely accurate in capturing the physical properties of a neutron star spacetime as they are calculated numerically in general relativity. Because the metric is given in terms of an expansion, the expressions are much simpler and easier to implement, in contrast to previous approaches. For the parameterisation of the metric in general relativity, the recently discovered universal 3-hair relations are used to produce a 3-parameter metric. Finally, a straightforward extension of this metric is given for scalar-tensor theories with a massless scalar field, which also admit a formulation in terms of an Ernst potential.

The bulk properties of spherically symmetric stars in general relativity can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations. In previous work we developed a "post-TOV" formalism - inspired by parametrized post-Newtonian theory - which allows us to classify in a parametrized, phenomenological form all possible perturbative deviations from the structure of compact stars in general relativity that may be induced by modified gravity at second post-Newtonian order. In this paper we extend the formalism to deal with the stellar exterior, and we compute several potential astrophysical observables within the post-TOV formalism: the surface redshift $z_s$, the apparent radius $R_{\rm app}$, the Eddington luminosity at infinity $L_{\rm E}^\infty$ and the orbital frequencies. We show that, at leading order, all of these quantities depend on just two post-TOV parameters $\mu_1$ and $\chi$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

We examine the effects of higher-order multipole contributions of rotating neutron star (NS) spacetimes on the propagation of corrugation (c-)modes within a thin accretion disk. We find that the Lense-Thirring precession frequency, which determines the propagation region of the low-frequency fundamental corrugation modes, can experience a turnover allowing for c-modes to become self-trapped for sufficiently high dimensionless spin $j$ and quadrupole rotational deformability $\alpha$. If such self-trapping c-modes can be detected, e.g. through phase-resolved spectroscopy of the iron line for a high-spin low-mass accreting neutron star, this could potentially constrain the spin-induced NS quadrupole and the NS equation of state.

In recent years, a lot of work was done that has revealed some very interesting properties of neutron stars. One can relate the first few multipole moments of a neutron star, or quantities that can be derived from them, with relations that are independent of the equation of state (EoS). This is a very significant result that has great implications for the description of neutron stars and in particular for the description of the spacetime around them. Additionally, it was recently shown that there is a four parameter analytic spacetime, known as the two-soliton spacetime, which can accurately capture the properties of the geometry around neutron stars. This allows for the possibility of describing in a unified formalism the astrophysically relevant properties of the spacetime around a neutron star independently of the particulars of the EoS for the matter of the star. More precisely, the description of these astrophysical properties is done using an EoS omniscient spacetime that can describe the exterior of any neutron star. In the present work we investigate properties such as the location of the innermost stable circular orbit $R_{ISCO}$ (or the surface of the star when the latter overcomes the former), the various frequencies of perturbed circular equatorial geodesics, the efficiency of an accretion disc, its temperature distribution, and other properties associated with the emitted radiation from the disc, in a way that holds for all possible choices of a realistic EoS for the neutron star. Furthermore, we provide proof of principle that if one were to measure the right combinations of pairs of these properties, with the additional knowledge of the mass of the neutron star, one could determine the EoS of the star.

Besides their astrophysical interest, compact stars also provide an arena for understanding the properties of theories of gravity that differ from Einstein's general relativity. Numerous studies have shown that different modified theories of gravity can modify the bulk properties (such as mass and radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. Motivated by this observation, in this paper we take an alternative route and construct a stellar structure formalism which, without adhering to any particular theory of gravity, describes in a simple parametrized form the departure from compact stars in general relativity. This "post-Tolman-Oppenheimer-Volkoff (TOV)" formalism for spherical static stars is inspired by the well-known parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable correction terms to the fully relativistic TOV equations. We show how neutron star properties are modified within our formalism, paying special attention to the effect of each correction term. We also show that the formalism is equivalent to general relativity with an "effective" (gravity-modified) equation of state.

The formalism for describing a metric and the corresponding scalar in terms of multipole moments has recently been developed for scalar-tensor theories. We take advantage of this formalism in order to obtain expressions for the observables that characterise geodesics in terms of the moments. These expressions provide some insight into how the structure of a scalarized compact object affects observables. They can also be used to understand how deviations from general relativity are imprinted on the observables.

Stationary, asymptotically flat spacetimes in general relativity can be characterized by their multipole moments. The moments have proved to be very useful tools for extracting information about the spacetime from various observables and, more recently, for establishing universalities in the structure of neutron stars. As a first step toward extending these methods beyond general relativity, we develop the formalism that allows one to define and calculate the multipole moments in scalar-tensor theories of gravity.

Black holes are said to have no hair because all of their multipole moments can be expressed in terms of just their mass, charge and spin angular momentum. The recent discovery of approximately equation-of-state-independent relations among certain multipole moments in neutron stars suggests that they are also approximately bald. We here explore the yet unknown origin for this universality. First, we investigate which region of the neutron star's interior and of the equation of state is most responsible for the universality. We find that the universal relation between the moment of inertia and the quadrupole moment is dominated by the star's outer-core, a shell of width (50-95)% of the total radius, which corresponds to the density range $(10^{14}$-$10^{15})$g/cm$^3$. Second, we study the impact on the universality of approximating stellar isodensity contours as self-similar ellipsoids. An analytical calculation in the Newtonian limit reveals that the shape of the ellipsoids does not affect the universal relations, but relaxing the self-similarity assumption can completely destroy it. Third, we investigate the eccentricity profiles of rotating relativistic stars and find that the ellipticity is roughly constant, with variations of roughly (20-30)% in the region that matters to the universal relations. Fourth, we repeat the above analysis for non-compact, regular stars and find that the ellipticity is not constant, with variations that easily exceed 100% and universality is lost. These findings suggest that universality arises as an emergent approximate symmetry: as one flows in the stellar-structure phase space from non-compact star region to the relativistic star region, the eccentricity variation inside stars decreases, leading to approximate self-similarity in their isodensity contours, which then leads to the universal behavior observed in their exterior multipole moments.

Astrophysical charge-free black holes are known to satisfy no-hair relations through which all multipole moments can be specified in terms of just their mass and spin angular momentum. We here investigate the possible existence of no-hair-like relations among multipole moments for neutron stars and quark stars that are independent of their equation of state. We calculate the multipole moments of these stars up to hexadecapole order by constructing uniformly-rotating and unmagnetized stellar solutions to the Einstein equations. For slowly-rotating stars, we construct stellar solutions to quartic order in spin in a slow-rotation expansion, while for rapidly-rotating stars, we solve the Einstein equations numerically with the LORENE and RNS codes. We find that the multipole moments extracted from these numerical solutions are consistent with each other. We confirm that the current-dipole is related to the mass-quadrupole in an approximately equation of state independent fashion, which does not break for rapidly rotating neutron stars or quark stars. We further find that the current-octupole and the mass-hexadecapole moments are related to the mass-quadrupole in an approximately equation of state independent way to $\sim 10%$, worsening in the hexadecapole case. All of our findings are in good agreement with previous work that considered stellar solutions to leading-order in a weak-field expansion. The quartic in spin, slowly-rotating solutions found here allow us to estimate the systematic errors in the measurement of the neutron star's mass and radius with future X-ray observations, such as NICER and LOFT. We find that the effect of these quartic-in-spin terms on the quadrupole and hexadecapole moments and stellar eccentricity may dominate the error budget for very rapidly-rotating neutron stars. The new universal relations found here should help to reduce such systematic errors.

Recently it was shown that slowly rotating neutron stars exhibit an interesting correlation between their moment of inertia $I$, their quadrupole moment $Q$, and their tidal deformation Love number $\lambda$ (the I-Love-Q relations), independently of the equation of state of the compact object. In the present work a similar, more general, universality is shown to hold true for all rotating neutron stars within General Relativity; the first four multipole moments of the neutron star are related in a way independent of the nuclear matter equation of state we assume. By exploiting this relation, we can describe quite accurately the geometry around a neutron star with fewer parameters, even if we don't know precisely the equation of state. Furthermore, this universal behavior displayed by neutron stars, could promote them to a more promising class of candidates (next to black holes) for testing theories of gravity.

In this article we present some recent results on identifying correctly the relativistic multipole moments of numerically constructed spacetimes, and the consequences that this correction has on searching for appropriate analytic spacetimes that can approximate well the previously mentioned numerical spacetimes. We also present expressions that give the quadrupole and the spin octupole as functions of the spin parameter of a neutron star for various equations of state and in a range of masses for every equation of state used. These results are relevant for describing the exterior spacetime of rotating neutron stars that are made up of matter obeying realistic equations of state.

We have tested the appropriateness of two-soliton analytic metric to describe the exterior of all types of neutron stars, no matter what their equation of state or rotation rate is. The particular analytic solution of the vaccuum Einstein equations proved quite adjustable to mimic the metric functions of all numerically constructed neutron-star models that we used as a testbed. The neutron-star models covered a wide range of stiffness, with regard to the equation of state of their interior, and all rotation rates up to the maximum possible rotation rate allowed for each such star. Apart of the metric functions themselves, we have compared the radius of the innermost stable circular orbit $R_{\rm{ISCO}}$, the orbital frequency $\Omega\equiv\frac{d\phi}{dt}$ of circular geodesics, and their epicyclic frequencies $\Omega_{\rho}, \Omega_z$, as well as the change of the energy of circular orbits per logarithmic change of orbital frequency $\Delta\tilde{E}$. All these quantities, calculated by means of the two-soliton analytic metric, fitted with good accuracy the corresponding numerical ones as in previous analogous comparisons (although previous attempts were restricted to neutron star models with either high or low rotation rates). We believe that this particular analytic solution could be considered as an analytic faithful representation of the gravitation field of any rotating neutron star with such accuracy, that one could explore the interior structure of a neutron star by using this space-time to interpret observations of astrophysical processes that take place around it.

We show how one can estimate the multipole moments of the space-time, assuming that the quasi-periodic modulations of the X-ray flux (quasi-periodic oscillations), observed from accreting neutron stars or black holes, are due to orbital and precession frequencies (relativistic precession model). The precession frequencies $\Omega_{\rho}$ and $\Omega_z$ can be expressed as expansions on the orbital frequency $\Omega$, in which the moments enter the coefficients in a prescribed form. Thus, observations can be fitted to these expressions in order to evaluate the moments. If the compact object is a neutron star, constrains can be imposed on the equation of state. The same analysis can be used for black holes as a test for the validity of the no-hair theorem. Alternatively, instead of fitting for the moments, observations can be matched to frequencies calculated from analytic models that are produced so as to correspond to realistic neutron stars described by various equations of state. Observations can thus be used to constrain the equation of state and possibly other physical parameters (mass, rotation, quadrupole, etc.) Some distinctive features of the frequencies, which become evident by using the analytic models, are discussed.

Identifying the relativistic multipole moments of a spacetime of an astrophysical object that has been constructed numerically is of major interest, both because the multipole moments are intimately related to the internal structure of the object, and because the construction of a suitable analytic metric that mimics a numerical metric should be based on the multipole moments of the latter one, in order to yield a reliable representation. In this note we show that there has been a widespread delusion in the way the multipole moments of a numerical metric are read from the asymptotic expansion of the metric functions. We show how one should read correctly the first few multipole moments (starting from the quadrupole mass-moment), and how these corrected moments improve the efficiency of describing the metric functions with analytic metrics that have already been used in the literature, as well as other consequences of using the correct moments.