Slowly Rotating Anisotropic Neutron Stars in General Relativity and Scalar-Tensor Theory
Some models (such as the Skyrme model, a low-energy effective field theory for QCD) suggest that the high-density matter prevailing in neutron star interiors may be significantly anisotropic. Anisotropy is known to affect the bulk properties of nonrotating neutron stars in General Relativity. In this paper we study the effects of anisotropy on slowly rotating stars in General Relativity. We also consider one of the most popular extensions of Einstein’s theory, namely scalar-tensor theories allowing for spontaneous scalarization (a phase transition similar to spontaneous magnetization in ferromagnetic materials). Anisotropy affects the moment of inertia of neutron stars (a quantity that could potentially be measured in binary pulsar systems) in both theories. We find that the effects of scalarization increase (decrease) when the tangential pressure is bigger (smaller) than the radial pressure, and we present a simple criterion to determine the onset of scalarization by linearizing the scalar-field equation. Our calculations suggest that binary pulsar observations may constrain the degree of anisotropy or even, more optimistically, provide evidence for anisotropy in neutron star cores.
pacs:04.40.Dg, 04.20.-q, 04.50.Kd, 21.65.Mn, 26.60.Kp
Most investigations of the structure of neutron stars (NSs) assume isotropic matter with a perfect-fluid equation of state (EoS) relating the pressure and density in the stellar interior. However, various physical effects can lead to local anisotropies (see  for a review). Anisotropy can occur for stars with a solid core  or strong magnetic fields [3, 4, 5]. Spaghetti- and lasagna-like structures would induce anisotropic elastic properties that could be important for NS quakes . Nuclear matter may be anisotropic at very high densities [7, 8], where the nuclear interactions must be treated relativistically and phase transitions (e.g. to pion condensates  or to a superfluid state ) may occur. For example, Nelmes and Piette  recently considered NS structure within the Skyrme model, a low energy, effective field theory for Quantum Chromodynamics (QCD), finding significant anisotropic strains for stars with mass (see also [12, 13]). From a mathematical point of view, two-fluid systems can be shown to be equivalent to a single anisotropic fluid . Anisotropy affects the bulk observable properties of NSs, such as the mass-radius relation and the surface redshift : it can increase the maximum NS mass for a given EoS [15, 16] and stabilize otherwise unstable stellar configurations . Incidentally, exotic objects such as gravastars  and boson stars [19, 20] are also equivalent to anisotropic fluids (i.e., they have anisotropic pressure).
It is known that rotation can induce anisotropy in the pressure due to anisotropic velocity distributions in low-density systems , but to the best of our knowledge – with the exception of some work by Bayin  – slowly rotating anisotropic stars have never been investigated in General Relativity (GR). The goal of this paper is to fill this gap using two different phenomenological models for anisotropy [15, 16], and to extend the analysis of slowly rotating anisotropic stars to scalar-tensor theories of gravity.
Scalar-tensor theories are among the simplest and best studied extensions of GR . In addition to the metric, in these theories gravity is also mediated by a scalar field. Scalar-tensor theories arise naturally from the dimensional reduction of higher-dimensional proposals to unify gravity with the Standard Model, and they encompass theories of gravity as special cases [23, 24]. The simplest variant of scalar-tensor theory, Brans-Dicke theory, is tightly constrained experimentally , but certain versions of the theory could in principle differ from GR by experimentally measurable amounts in the strong-field regime, as shown by Damour and Esposito-Farèse [26, 27].
From an astrophysical standpoint, compact objects such as black holes and NSs are the most plausible candidates to test strong-field gravity . Compared to black holes, NSs are a more promising strong-field laboratory to distinguish scalar-tensor gravity from GR, because a large class of scalar-tensor theories admits the same black-hole solutions as GR (see  and references therein), and the dynamics of black holes can differ from GR only if the black holes are surrounded by exotic forms of matter [30, 31, 32, 33] or if the asymptotic behavior of the scalar field is nontrivial [34, 35].
The study of NS structure in GR is textbook material [36, 37, 38, 39], and there is an extensive literature on stellar configurations in scalar-tensor theories as well (see e.g. [40, 41] and references therein). One of the most intriguing phenomena in this context is “spontaneous scalarization” , a phase transition analogous to the familiar spontaneous magnetization in solid state physics : in a certain range of central densities, asymptotically flat solutions with a nonzero scalar field are possible and energetically favored with respect to the corresponding GR solutions.
In the absence of anisotropy, the degree of scalarization depends on a certain (real) theory parameter , defined in Eq. (2) below. Theory predicts that scalarization cannot occur (in the absence of anisotropy) when . Present binary pulsar observations yield a rather tight experimental constraint: [44, 45]. One of our main findings is that the effects of scalarization, as well as the critical for spontaneous scalarization to occur, increase (decrease) for configurations in which the tangential pressure is bigger (smaller) than the radial pressure. Therefore binary pulsars can be used to constrain the degree of anisotropy at fixed , or to constrain for a given degree of anisotropy. This may open the door to experimental constraints on the Skyrme model via binary pulsar observations. Other notable findings of this study are (i) an investigation of the dependence of the stellar moment of inertia on the degree of anisotropy (more precisely, and , because we consider two different anisotropy models [15, 16]); and (ii) an investigation of the threshold for scalarization for different values of and in terms of a simple linear stability criterion, along the lines of recent work for black holes surrounded by matter [32, 33].
The plan of the paper is as follows. In Section 2 we introduce the equations of motion in scalar-tensor theory and the stress-energy tensor describing anisotropic fluids that will be used in the rest of the paper. In Section 3 we present the equations of structure for relativistic stars at first order in the slow-rotation expansion. The macroscopic properties of NSs obtained by integrating these equations for two different models of anisotropic stars are presented in Section 4. Section 4.3 shows that a linear approximation is sufficient to identify the threshold for spontaneous scalarization for different values of and . Section 5 summarizes our main conclusions and points out possible avenues for future work. Finally, in Appendix A we give a detailed derivation of an integral formula to compute the moment of inertia. Throughout this work, quantities associated with the Einstein (Jordan) frame will be labeled with an asterisk (tilde). We use geometrical units () unless stated otherwise and signature .
2 Anisotropic fluids in scalar-tensor theory of gravity
2.1 Overview of the theory
where is the bare gravitational constant, is the determinant of the Einstein-frame metric , is the Ricci curvature scalar of the metric , and is a massless scalar field. is the action of the matter fields, collectively represented by . Free particles follow geodesics of the Jordan-frame metric , where is a conformal factor. In this work we assume that has the form
The field equations of this theory, obtained by varying the action with respect to and , are given by
where is the Ricci tensor, (in the language of [27, 42]) is the “scalar-matter coupling function” and is the d’Alembertian operator associated to the metric . GR is obtained in the limit where the scalar field decouples from matter, i.e. . Under the particular choice of the conformal factor (2), this is equivalent to letting . In this paper, all equations will be derived within the context of scalar-tensor gravity.
Finally, is the energy-momentum tensor of the matter fields, defined as
and is its trace. The energy-momentum tensor in the Jordan frame , with trace , is defined in an analogous fashion:
The two energy-momentum tensors (and their traces) are related as follows:
The covariant divergence of the energy-momentum tensor satisfies
in the Einstein and Jordan frames, respectively.
2.2 Anisotropic fluids
where is the fluid four-velocity, is a unit radial vector () satisfying , and . is a projection operator onto a two-surface orthogonal to both and : indeed, defining a projected vector , one can easily verify that . At the center of symmetry of the fluid distribution the tangential pressure must vanish, since is not defined there . The trace of the Einstein-frame stress-energy tensor for an anisotropic fluid is
As emphasized by Bowers and Liang , and contain contributions from fluid pressures and other possible stresses inside the star, therefore they should not be confused with purely hydrostatic pressure. Additional stresses could be caused, for instance, by the presence of a solid core , strong magnetic fields  or a multi-fluid mixture . The derivation of a microphysical model for anisotropy is a delicate issue, so we will adopt a phenomenological approach. We will assume that is described by a barotropic EoS, i.e. . For brevity in this paper we focus on the APR EoS , but we have verified that our qualitative results do not depend on this choice. The APR EoS supports NS models with a maximum mass larger than , and therefore it is compatible with the recent observations of the pulsar PSR J1614-2230  and of the pulsar PSR J0348+0432 .
The functional form of the anisotropy [15, 46, 50] depends on microscopic relationships between , and , that unfortunately are not known. However we can introduce physically motivated functional relations for that allow for a smooth transition between the isotropic and anisotropic regimes. Many such functional forms have been studied in the literature. As an application of our general formalism we will consider two of these phenomenological relations, described below.
2.2.1 Quasi-local equation of state
Horvat et al.  proposed the following quasi-local equation for :
where . The “mass function” , defined in Eq. (15) below, is essentially the mass contained within the radius , so the quantity is a local measure of compactness, whereas is a free (constant) parameter that controls the degree of anisotropy.
The calculations of  show that, if anisotropy occurs due to pion condensation, , therefore could be of order unity . More recently, Nelmes and Piette  considered NS structure within a model consisting of a Skyrme crystal, which allows for the presence of anisotropic strains. They found that , as defined in Eq. (12), has a nearly constant value throughout the NS interior. The nonradial oscillations of anisotropic stars were studied in  using the model of Eq. . Following Doneva and Yazadjiev , we will consider values of in the range .
2.2.2 Bowers-Liang model
As a second possibility we will consider the functional form for proposed by Bowers and Liang , who suggested the relation111The factor of 1/3 in Eq. (13) is chosen for convenience. Also, there is a sign difference between our definition of and the one in . Our parameter is related with the Bowers-Liang (physically equivalent) parameter by .
The model is based on the following assumptions: (i) the anisotropy should vanish quadratically at the origin (the necessity for this requirement will become clear in Sec. 3), (ii) the anisotropy should depend nonlinearly on , and (iii) the anisotropy is (in part) gravitationally induced. The parameter controls the amount of anisotropy in the fluid.
This ansatz was used in  to obtain an exact solution for incompressible stars with constant. In their simple model, the requirement that equilibrium configurations should have finite central pressure implies that . The Newtonian limit of the Bowers-Liang ansatz was also considered in a recent study of the correspondence between superradiance and tidal friction . In our calculations we will assume that .
3 Stellar structure in the slow-rotation approximation
In this Section we approximate the metric of a slowly, rigidly rotating, anisotropic star following the seminal work by Hartle and Thorne [51, 52]. The idea is to consider the effects of rotation as perturbations of the spherically symmetric background spacetime of a static star. We generalize the results of [51, 52] (in GR) and  (in scalar-tensor theory) to account for anisotropic fluids up to first order in rotation, so we can study how anisotropy and scalarization affect the moment of inertia of the star and the dragging of inertial frames.
We remark that the moment of inertia , the star’s uniform angular velocity and the angular momentum are the same in the Jordan and Einstein frames (cf. [42, 53]). Therefore, to simplify the notation, we will drop asterisks and tildes on these quantities. Working at order , the line element of a stationary axisymmetric spacetime in the Jordan frame reads
is the mass function and is the angular velocity acquired by a particle falling from infinity as measured by a static asymptotic observer .
The four-velocity of the rotating fluid is such that , and it has components 
Using (14), at first order in the slow-rotation parameter we obtain:
where we defined , , and . The equations above reduce to the Tolman-Oppenheimer-Volkoff (TOV) equations for anisotropic stars in GR  when , to the results of  in the isotropic limit , and to the usual TOV equations when both quantities are equal to zero . In the GR limit, our frame-dragging equation (24) agrees with Bayin’s  result222In principle, as mentioned in the introduction, rotation may induce anisotropy. Therefore the Horvat et al. and Bowers-Liang models for should contain terms proportional to . However, Eq. (24) implies that such terms in would lead to corrections of second order in the angular velocity . These corrections are beyond the scope of the approximation considered in our work..
To obtain the interior solution we integrate the generalized TOV equations (20)-(24) from a point close to the stellar center outwards up to a point where the pressure vanishes, i.e. . This point specifies the Einstein-frame radius of the star. If , the Jordan-frame radius is
In practice, to improve numerical stability, given , , and (where the subscript c means that all quantities are evaluated at ) we use the following series expansions:
where and depend on the particular anisotropy model.
In the vacuum exterior we have . Eqs. (20)–(22) must be integrated outwards starting from the stellar radius to obtain the stellar mass, angular momentum and scalar charge. For large we can expand the relevant functions as follows:
where is the Arnowitt-Deser-Misner (ADM) mass of the NS, is the scalar charge, is the star’s angular momentum and is the (constant) cosmological value of the scalar field, here assumed to be zero. Under this assumption the mass is the same in the Jordan and Einstein frames . By matching the numerical solution integrated from the surface of the star with the asymptotic expansions (27)–(30) we can compute , and .
We compute the moment of inertia of the star in two equivalent ways. The first method consists of extracting the angular momentum as described above and using
(see Appendix A for details). As and we recover Hartle’s result , and in the isotropic limit we match the result of . The numerical values of obtained with (31) and (32) are in excellent agreement.
For each stellar model we also calculate the baryonic mass , defined as 
where g is the atomic mass unit and is the baryonic number density.
4 Numerical results
The tools developed so far allow us to investigate the effect of anisotropy on the bulk properties of rotating stars. In Section 4.1 we will focus on slowly rotating stars in GR. To the best of our knowledge – and to our surprise – rotating anisotropic stars have not been studied in the GR literature, with the only exception of a rather mathematical paper by Bayin . In Section 4.2 we extend our study to scalar-tensor theories. Our main motivation here is to understand whether anisotropy may increase the critical value above which spontaneous scalarization cannot happen, and therefore allow for observationally interesting modifications to the structure of NSs that would still be compatible with the stringent bounds from binary pulsars [44, 45].
4.1 The effect of anisotropy in GR
In the top panels of Figure 1 we show the mass-radius relation for anisotropic NS models in GR. All curves are truncated at the central density corresponding to the maximum NS mass, because models with larger central densities are unstable to radial perturbations [36, 37]. Solid lines correspond to , i.e. the isotropic fluid limit. The horizontal shaded band in the upper panels represents the largest measured NS mass (PSR J0348+0432: cf. ).
Recall that is proportional to and (with a positive proportionality constant) in both models, and that and represent the “radial” and “tangential” pressures, respectively. Therefore positive values of and mean that the radial pressure is larger than the tangential pressure (dashed lines); the opposite is true when the anisotropy parameters are negative (dotted lines).
The trend in the top panels of Figure 1 is clear: for both anisotropy models, positive (negative) anisotropy parameters yield smaller (larger) radii at fixed mass, and smaller masses at fixed radius. The lower panels of Figure 1 show that the stellar compactness decreases (for a given EoS and fixed central density) as the anisotropy degree increases. Nuclear matter EoSs are usually ordered in terms of a “stiffness” parameter, with stiffer EoSs corresponding to larger sound speeds (more incompressible matter) in the stellar interior, and larger values of the compactness . The qualitative effect of increasing anisotropy (with our sign conventions) is opposite (for a given EoS) to the qualitative effect of increasing stiffness.
Figure 2 is, to our knowledge, the first calculation of the effect of anisotropy on the moment of inertia . As in Figure 1, solid lines corresponds to the isotropic limit. In the right panel we use the quasi-local model of ; in the left panel, the Bowers-Liang model . Hypothetical future observations of the moment of inertia of star A, from the double pulsar PSR J0737-3039 [55, 56, 57], or preferably from large-mass NSs, may be used to constrain the degree of anisotropy under the assumptions that GR is valid and that the nuclear EoS is known.
4.2 The effect of anisotropy on spontaneous scalarization
In Figures 3 and 4 we display the properties of nonrotating, spontaneously scalarized stars within the anisotropy models of Horvat et al.  and Bowers-Liang , respectively. The main panel in each Figure shows the mass-radius relation as the anisotropy parameter increases (in increments of 1, and from top to bottom) in the range (Figure 3) or (Figure 4). Solid lines correspond to the GR limit; dotted, dashed and dash-dotted lines correspond to and , as indicated in the legend. The lower panels show the scalar charge as a function of the baryonic mass. In each of these panels we plot the scalar charge for a fixed value of and different anisotropy parameters.
For isotropic EoSs in GR, Harada  used catastrophe theory to show that scalarization is only possible when . We find that scalarization can occur for larger values of in the presence of anisotropy. For example, for a value of (compatible with the Skyrme model predictions of ) scalarization is possible when , and for scalarization produces rather large () deviations in the mass-radius relation. This qualitative conclusion applies to both anisotropy models considered by us. The lower panels show that: (i) for fixed (i.e., for a fixed theory) and for a fixed central density, the “strength” of scalarization – as measured by the scalar charge of the star – increases for large negative ’s, i.e. when the tangential pressure is significantly larger than the radial pressure, for both anisotropy models; (ii) scalarization occurs in a much wider range of baryonic masses, all of which are compatible with the range where anisotropy would be expected according to the Skyrme model predictions of . These calculations are of course preliminary and should be refined using microphysical EoS models. However, let us remark once again that the scalarization threshold in the absence of anisotropy is to a very good approximation EoS-independent, and stars only acquire significant scalar charge when (as shown in  and in Figure 6 below). In the admittedly unlikely event that binary pulsar observations were to hint at scalarization with , this would be strong evidence for the presence of anisotropy333An important caveat here is that fast rotation can also strengthen the effects of scalarization: according to , scalarization can occur for for NSs spinning at the mass-shedding limit. However the NSs found in binary pulsar systems are relatively old, as they are expected to be spinning well below the mass-shedding limit, where the slow-rotation approximation works very well . and even lead to experimental constraints on the Skyrme model and QCD.
As in Figure 2, in the left panel of Figure 5 we show the moment of inertia as a function of the stellar mass for the quasi-local model of , while the right panel refers to the Bowers-Liang model . Solid lines corresponds to the GR limit for different anisotropy parameters. Unsurprisingly, the largest modifications to the moment of inertia occur for large negative ’s, and they follow the same trends highlighted in our discussion of the mass-radius relation.
4.3 Critical scalarization point in the linearized approximation
The condition for spontaneous scalarization to occur can be found in a linearized approximation to the scalar-field equation of motion. The idea is that at the onset of scalarization the scalar field must be small, so we can neglect its backreaction on the geometry and look for bound states of the scalar field by dropping terms quadratic in the field [27, 42]. Here we study general conditions for the existence of bound states in the linearized regime, and we show that (as expected based on the previous argument) the linearized theory does indeed give results in excellent agreement with the full, nonlinear calculation.
Redefining the scalar field as and neglecting terms , Eq. (4) can be written as a Schrödinger-like equation:
where the tortoise radial coordinate is defined by . The effective potential is
where we have introduced an effective (position-dependent) mass
Eq. (34) with the potential (35) is a wave equation for a scalar field with effective mass . From Eq. (11) we see that anisotropy affects the effective mass (and therefore the scalarization threshold) because contains a term proportional to , that in turn is proportional to either or : cf. Eqs. (12) and (13). The case of spontaneous scalarization around black holes (studied in [32, 33]) can be recovered by setting .
The scalarization threshold can be analyzed by looking for the zero-energy () bound state solutions of Eq. (34). In this case, the scalar field satisfies the following boundary conditions:
and we impose , where the prime denotes derivative with respect to . To obtain the scalarization threshold we integrate Eq. (34) outwards, starting from , with the above boundary conditons. Since the equation is linear, is arbitrary. At infinity we impose that the first derivative of with respect to must be zero. This is a two-point boundary value problem that can be solved with a standard shooting method to find the critical value of the central density for which the above conditions are satisfied, given fixed values of and (or ). The solution is some
where is the smallest critical density at which scalarization can occur for the given . The largest critical density producing scalarization can be similarly obtained by looking for zero-energy bound state solutions to find some
It can be shown that in these two regimes (i.e., at the starting and ending points of the scalarization regime) the derivative of with respect to has opposite signs:
As a warm-up, in Figure 6 we compute the scalarization threshold for nonrotating isotropic stars with several nuclear-physics based EoSs. The original references for the subset of EoSs used here can be found in  (the one exception is SLy4: cf. ). The EoSs are sorted by stiffness, with APR EoS being the stiffest and G EoS the softest in our catalog. As a trend, for stiffer EoSs scalarization occurs at lower values of the central densities and at higher values of the compactness. The most remarkable fact is that the value above which scalarization cannot occur is very narrow: it ranges from for APR EoS to for F EoS . This is consistent with Harada’s study based on catastrophe theory, that predicts a threshold value (horizontal line in the figure) in the absence of anisotropy  (see also ).
In Table 1 we compare the values for and computed using (i) the linearized method described in this Section, and (ii) the full nonlinear set of equations for anisotropic models constructed using the APR EoS. The results agree remarkably well, showing that the onset of scalarization can be analyzed to an excellent degree of accuracy by neglecting the backreaction effects of the scalar field on the geometry. The last column of Table 1 lists , the value of above which scalarization cannot happen. We do not present results for because the resulting is already ruled out by binary pulsar observations .
In the left panels of Figure 7 we analyze the dependence of the critical on the central density, focusing on EoS APR and selecting different values of the anisotropy parameters (top) and (bottom). The shaded region at the top () is allowed by current binary pulsar observations [44, 45]. The horizontal line is the roughly EoS-independent threshold for isotropic stars. For a given theory, the starting and ending points of the scalarization regime are those for which a (horizontal) line crosses the curves. Anisotropic models have two distinctive features: (1) when the tangential pressure is larger than the radial pressure (dashed lines in Figure 7) scalarization can occur even for (for example, for the Horvat et al. model with we have , and for the Bowers-Liang model with we have ; cf. Table 1, Figure 4 and Figure 3); (2) when the tangential pressure is smaller than the radial pressure (dash-dotted lines in Figure 7) scalarized solutions may exist for a much wider range of .
In the right panels of Figure 7 we plot the critical as a function of the stellar compactness . For low compactness () all curves have the same behaviour regardless of or . This universality has two reasons: (1) all modern nuclear-physics based EoS have the same Newtonian limit (cf.  for an analytic treatment of this regime for constant density stars); (2) for any given EoS, the effects of anisotropy are suppressed in the Newtonian regime, where pressures and densities are low and the local compactness parameter is small: cf. Eqs. (12) and (13).
Binary pulsar observations require [44, 45], and even more stringent constraints are expected in the near future. As shown in Figure 6, most “ordinary” nuclear-physics based EoSs for nuclear matter predict that scalarization can only occur for . As binary pulsar observations get closer and closer to the limit , the spontaneous scalarization mechanism originally proposed by Damour and Esposito-Farèse [27, 42] looks more and more unlikely to be realized in Nature if neutron stars are isotropic.
The admittedly unlikely event of a binary-pulsar observation of scalarization with would be strong evidence for the presence of anisotropy, and it may even lead to experimental constraints on the Skyrme model and QCD. An important caveat here is that fast rotation can also strengthen the effects of scalarization: according to , scalarization can occur for when NSs spin at the mass-shedding limit. However the NSs found in binary pulsar systems are relatively old, are they are expected to spin well below the mass-shedding limit, where the slow-rotation approximation works very well .
Our work can be extended in several directions. An obvious extension is to consider the effects of anisotropy at second or higher order in the Hartle-Thorne expansion. This would allow us to assess whether the recently discovered “I-Love-Q” and “three-hair” universal relations between the multipole moments of the spacetime hold in the presence of anisotropy and scalarization [65, 66, 67, 68, 53]. A second obvious extension could consider fast rotating, anisotropic stars (cf. [58, 69]) and the orbital and epicyclic frequencies around these objects .
Anisotropy can lower the threshold for scalarization to occur, and this could be of interest to test scalar-tensor theories through gravitational-wave asteroseismology [71, 72, 73]. We also remark that our study used simplified, phenomenological models for anisotropy, when of course it would be desirable to study realistic microphysical models. Last but not least, our study should be extended to evaluate the stellar sensitivities [74, 75] and to identify exclusion regions in the parameter space using binary pulsar observations (cf. e.g. ).
We are grateful to K. Glampedakis, M. Horbatsch and G. Pappas for discussions, and to G. Pappas for providing EoS data tables. H. O. S. and E. B. were supported by NSF CAREER Grant No. PHY-1055103. C. F. B. M and L. C. B. C would like to thank the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), and Fundação Amazônia de Amparo a Estudos e Pesquisas do Pará (FAPESPA) for partial financial support. The authors also acknowledge support from the FP7-PEOPLE-2011-IRSES Grant No.295189. The figures were created using the python-based library matplotlib .
Appendix A Derivation of equation (32)
In this Appendix we present a derivation of the integral (32), used to compute the moment of inertia of slowly rotating stars in scalar-tensor theory. We begin by noting that
Multiplying the frame dragging equation (24) by and rearranging, we obtain
If we multiply by , integrate from to infinity and use the fact that
as , we finally get Eq. (32).
-  Herrera L and Santos N O 1997 Phys. Rep. 286 53–130
-  Kippenhahn R and Weigert A 1990 Stellar structure and evolution Astronomy and astrophysics library (Springer) ISBN 9783540502111
-  Yazadjiev S S 2012 Phys. Rev. D 85(4) 044030
-  Folomeev V and Dzhunushaliev V 2015 Phys.Rev. D91 044040 (Preprint 1501.06275)
-  Kamiab F, Broderick A E and Afshordi N 2015 (Preprint 1503.03898)
-  Heiselberg H and Hjorth-Jensen M 2000 Phys.Rept. 328 237–327 (Preprint nucl-th/9902033)
-  Ruderman M 1972 Ann.Rev.Astron.Astrophys. 10 427–476
-  Canuto V and Chitre S 1974 Phys.Rev. D9 1587–1613
-  Sawyer R 1972 Phys.Rev.Lett. 29 382–385
-  Carter B and Langlois D 1998 Nucl.Phys. B531 478–504 (Preprint gr-qc/9806024)
-  Nelmes S and Piette B 2012 Phys.Rev. D85 123004 (Preprint 1204.0910)
-  Adam C, Naya C, Sanchez-Guillen J, Vazquez R and Wereszczynski A 2015 Phys.Lett. B742 136–142 (Preprint 1407.3799)
-  Adam C, Naya C, Sanchez-Guillen J, Vazquez R and Wereszczynski A 2015 (Preprint 1503.03095)
-  Letelier P S 1980 Phys. Rev. D 22(4) 807–813
-  Bowers R L and Liang E P T 1974 ApJ 188 657
-  Horvat D, Ilijic S and Marunovic A 2011 Class.Quant.Grav. 28 025009 (Preprint 1010.0878)
-  Dev K and Gleiser M 2003 Gen.Rel.Grav. 35 1435–1457 (Preprint gr-qc/0303077)
-  Cattoen C, Faber T and Visser M 2005 Class.Quant.Grav. 22 4189–4202 (Preprint gr-qc/0505137)
-  Schunck F and Mielke E 2003 Class.Quant.Grav. 20 R301–R356 (Preprint 0801.0307)
-  Macedo C F, Pani P, Cardoso V and Crispino L C B 2013 Phys.Rev. D88 064046 (Preprint 1307.4812)
-  Bayin S S 1982 Phys.Rev. D26 1262
-  Fujii Y and Maeda K I 2003 The Scalar-Tensor Theory of Gravitation (Cambridge University Press)
-  Sotiriou T P and Faraoni V 2010 Rev. Mod. Phys. 82 451–497 arXiv:0805.1726 [gr-qc]
-  De Felice A and Tsujikawa S 2010 Living Rev.Rel. 13 3 (Preprint 1002.4928)
-  Will C M 2006 Living Reviews in Relativity 9 URL http://www.livingreviews.org/lrr-2006-3
-  Damour T and Esposito-Farèse G 1992 Class.Quant.Grav. 9 2093–2176
-  Damour T and Esposito-Farèse G 1993 Phys.Rev.Lett. 70 2220–2223
-  Berti E, Barausse E, Cardoso V, Gualtieri L, Pani P et al. 2015 (Preprint 1501.07274)
-  Sotiriou T P and Faraoni V 2012 Phys.Rev.Lett. 108 081103 (Preprint 1109.6324)
-  Stefanov I Z, Yazadjiev S S and Todorov M D 2008 Mod.Phys.Lett. A23 2915–2931 (Preprint 0708.4141)
-  Doneva D D, Yazadjiev S S, Kokkotas K D and Stefanov I Z 2010 Phys.Rev. D82 064030 (Preprint 1007.1767)
-  Cardoso V, Carucci I P, Pani P and Sotiriou T P 2013 Phys.Rev. D88 044056 (Preprint 1305.6936)
-  Cardoso V, Carucci I P, Pani P and Sotiriou T P 2013 Phys.Rev.Lett. 111 111101 (Preprint 1308.6587)
-  Horbatsch M and Burgess C 2012 JCAP 1205 010 (Preprint 1111.4009)
-  Berti E, Cardoso V, Gualtieri L, Horbatsch M and Sperhake U 2013 Phys.Rev. D87 124020 (Preprint 1304.2836)
-  Harrison B K, Thorne K S, Wakano M and Wheeler J A 1965 Gravitation Theory and Gravitational Collapse (Chicago: University of Chicago Press)
-  Misner C W, Thorne K S and Wheeler J A 1973 Gravitation (San Francisco: W.H. Freeman and Co.)
-  Shapiro S and Teukolsky S 1983 Black holes, white dwarfs, and neutron stars: The physics of compact objects (Wiley)
-  Friedman J L and Stergioulas N 2013 Rotating Relativistic Stars (Cambridge: Cambridge University Press)
-  Horbatsch M and Burgess C 2011 JCAP 1108 027 (Preprint 1006.4411)
-  Pani P, Berti E, Cardoso V and Read J 2011 Phys.Rev. D84 104035 (Preprint 1109.0928)
-  Damour T and Esposito-Farèse G 1996 Phys.Rev. D54 1474–1491 (Preprint gr-qc/9602056)
-  Harada T 1998 Phys.Rev. D57 4802–4811 (Preprint gr-qc/9801049)
-  Freire P C, Wex N, Esposito-Farese G, Verbiest J P, Bailes M et al. 2012 Mon. Not. Roy. Astron. Soc. 423 3328 (Preprint 1205.1450)
-  Wex N 2014 (Preprint 1402.5594)
-  Doneva D D and Yazadjiev S S 2012 Phys.Rev. D85 124023 (Preprint 1203.3963)
-  Akmal A, Pandharipande V and Ravenhall D 1998 Phys.Rev. C58 1804–1828 (Preprint nucl-th/9804027)
-  Demorest P B, Pennucci T, Ransom S M, Roberts M S E and Hessels J W T 2010 Nature 467 1081–1083 (Preprint 1010.5788)
-  Antoniadis J, Freire P C, Wex N, Tauris T M, Lynch R S et al. 2013 Science 340 6131 (Preprint 1304.6875)
-  Glampedakis K, Kapadia S J and Kennefick D 2014 Phys.Rev. D89 024007 (Preprint 1312.1912)
-  Hartle J B 1967 Astrophys.J. 150 1005–1029
-  Hartle J B and Thorne K S 1968 Astrophys.J. 153 807
-  Pani P and Berti E 2014 Phys.Rev. D90 024025 (Preprint 1405.4547)
-  Staykov K V, Doneva D D, Yazadjiev S S and Kokkotas K D 2014 J. Cosmology Astropart. Phys 10 006 (Preprint 1407.2180)
-  Lyne A, Burgay M, Kramer M, Possenti A, Manchester R et al. 2004 Science 303 1153–1157 (Preprint astro-ph/0401086)
-  Lattimer J M and Schutz B F 2005 Astrophys.J. 629 979–984 (Preprint astro-ph/0411470)
-  Kramer M and Wex N 2009 Class.Quant.Grav. 26 073001
-  Doneva D D, Yazadjiev S S, Stergioulas N and Kokkotas K D 2013 Phys.Rev. D88 084060 (Preprint 1309.0605)
-  Berti E, White F, Maniopoulou A and Bruni M 2005 Mon.Not.Roy.Astron.Soc. 358 923–938 (Preprint gr-qc/0405146)
-  Kokkotas K and Ruoff J 2001 Astron.Astrophys. 366 565 (Preprint gr-qc/0011093)
-  Douchin F and Haensel P 2001 Astron.Astrophys. 380 151 (Preprint astro-ph/0111092)
-  Arponen J 1972 Nucl.Phys. A191 257–282
-  Novak J 1998 Phys.Rev. D58 064019 (Preprint gr-qc/9806022)
-  Pani P, Cardoso V, Berti E, Read J and Salgado M 2011 Phys.Rev. D83 081501 (Preprint 1012.1343)
-  Yagi K and Yunes N 2013 Science 341 365–368 (Preprint 1302.4499)
-  Yagi K and Yunes N 2013 Phys.Rev. D88 023009 (Preprint 1303.1528)
-  Pappas G and Apostolatos T A 2014 Phys.Rev.Lett. 112 121101 (Preprint 1311.5508)
-  Yagi K, Kyutoku K, Pappas G, Yunes N and Apostolatos T A 2014 Phys.Rev. D89 124013 (Preprint 1403.6243)
-  Doneva D D, Yazadjiev S S, Staykov K V and Kokkotas K D 2014 Phys.Rev. D90 104021 (Preprint 1408.1641)
-  Doneva D D, Yazadjiev S S, Stergioulas N, Kokkotas K D and Athanasiadis T M 2014 Phys.Rev. D90 044004 (Preprint 1405.6976)
-  Sotani H and Kokkotas K D 2004 Phys.Rev. D70 084026 (Preprint gr-qc/0409066)
-  Sotani H and Kokkotas K D 2005 Phys.Rev. D71 124038 (Preprint gr-qc/0506060)
-  Silva H O, Sotani H, Berti E and Horbatsch M 2014 Phys.Rev. D90 124044 (Preprint 1410.2511)
-  Will C M and Zaglauer H W 1989 Astrophys.J. 346 366
-  Zaglauer H 1992 Astrophys.J. 393 685–696
-  Alsing J, Berti E, Will C M and Zaglauer H 2012 Phys.Rev. D85 064041 (Preprint 1112.4903)
-  Hunter J D 2007 Computing In Science & Engineering 9 90–95