Models of Compact Stars in the Bimetric Scalar-Tensor Theory of Gravitation

We investigate static spherically-symmetric configurations of gravitating masses in the bimetric scalar-tensor theory of gravitation. In the gravitational sector, the theory contains the metric tensor, a scalar field and a background metric as an absolute variable of the theory. The analysis is presented for the simplest version of the theory with a constant coupling function and a zero cosmological function. We show that, depending on the value of the theory parameter, the masses for superdense compact configurations can be essentially larger compared to the configurations in general relativity.


Introduction
The scalar-tensor theories are among the most popular alternatives of general relativity.Another class of modified gravity theories, the so-called f (R)-gravities, are also presented in the form of a scalar-tensor theory with the potential for a scalar field determined by the f (R) function in the gravitational Lagrangian density.The recent activity in scalar-tensor theories is closely related to unification schemes (models with extra dimensions, string theories) and to the construction of models for dark energy, dark matter and inflation.In scalar-tensor theories, the field added in the gravitational sector is a scalar one.In bimetric (or tensor-tensor) theories of gravity, the additional field is a second rank tensor (for recent reviews in various types of modified gravity theories, see, for instance, [1,2]).The latter can be either dynamical or non-dynamical (absolute).
Recent activity in considering theories of gravity with two metric tensors is partially related to non-linear extensions of Fierz-Pauli massive gravity (for a review, see [3]).In these models, the second tensor is required for the construction of non-linear generalizations of the Fierz-Pauli mass term.The existence of a particular set of ghost-free nonlinear interactions has been shown in [4][5][6].In the corresponding models, the second metric is non-dynamical and flat.Models with a general dynamical metric tensor have been discussed in [7][8][9], and the cosmological dynamics of these models were considered in [10][11][12][13] (for a more complete list of references on recent investigations within the framework of nonlinear massive gravity, see, for instance, [14]).More general F(R) massive theories of gravity have been studied in [15,16].Similar to the case of standard F(R)-gravities, these models can be presented in the form of a scalar-tensor theory.
In [17][18][19][20], we have suggested a variant of scalar-tensor theories involving a second non-dynamical metric tensor: the bimetric scalar-tensor theory (BSTT).The BSTT belongs to the class of metric theories of gravity and consequently obeys the Einstein equivalence principle (for a general review, see [21]).If the value of the coupling function of the theory at the recent stage of the universe's expansion is not too small, the post-Newtonian parameters of BSTT coincide with those in general relativity.This is in contrast to usual scalar-tensor theories (an example is the Brans-Dicke theory), where the comparison of the post-Newtonian approximation of the theory with the observational data imposes relatively strong restrictions on the parameters.Another difference is that, because of the presence of the second metric tensor, the gravitational field is characterized by the energy-momentum tensor [22,23].The gravitational waves in BSTT have been analyzed in [24,25].In this theory, the velocity of weak perturbations of the curved metric coincides with the speed of light in a vacuum.In a variant of the theory with zero cosmological function, the same is the case for the scalar field.Similar to general relativity, the BSTT is of class N 2 in E(2) classification of the gravitational wave polarization of metric theories of gravity.In particular, this means that in BSTT, the weak perturbations for the metric and scalar field propagate independently.In usual scalar-tensor theories (for example, in Brans-Dicke theory), the equation for weak perturbations of the metric tensor contain a contribution from the scalar field, as well (see, for example, [21]).As a consequence of that, the scalar and metric perturbations are mixed, and the theory belongs to the class N 3 .Related to the presence of a scalar degree of freedom, in addition to the standard quadrupole gravitational radiation, in BSTT, there is also dipole gravitational radiation.The corresponding Peters-Mathews parameters for the radiation from a gravitating system with small velocities and nonrelativistic internal structure have been determined in [24].The corresponding formalism for a system of gravitating bodies with small velocities, but with relativistic internal structure (modified Einstein-Infeld-Hoffmann formalism) has been considered in [25].
In the present paper, we consider a static, spherically-symmetric configuration of gravitating masses within the framework of BSTT.In a variant of the theory with a constant coupling function and a zero cosmological function, the corresponding solution outside the matter distribution (external solution) has been found in [18].The results for a numerical integration of the internal equations were presented in [18,26] for some special cases of the equation of state for gravitating masses.The paper is organized as follows.In Section 2, we present the action and the field equations in BSTT.The equations of the theory for a static, spherically-symmetric distribution of the matter are considered in Section 3. The solution outside the matter distribution is discussed.In Section 4, the results of the numerical investigations are presented for superdense stellar configurations corresponding to white dwarfs and neutron stars.The main results are summarized in Section 5.

BSTT Action and the Field Equations
The BSTT belongs to the class of metric theories of gravity with a preferred geometry.In addition to the curved metric tensor g ik , it contains a dynamical scalar field ϕ(x) and a non-dynamical metric γ ik .The latter is the absolute variable of the theory.The action of the theory in its general formulation is given by the expression (in units c = 1): where ζ(ϕ) is a dimensionless coupling function, Λ(ϕ) is the cosmological function, L m is the Lagrangian density for nongravitational fields q a and ∇ l is the covariant derivative operator with respect to the metric g ik .The function Λ g is given by the expression: where Γl ik = Γ l ik − Γl ik is the affine deformation tensor and Γ l ik and Γl ik are the Christoffel symbols for the metrics g ik and γ ik , respectively.Note that the affine deformation tensor can be written in the form: with ∇i being the covariant derivative with respect to the metric γ ik .The action (1) contains the derivatives of the fields less than the second order and two functions: the coupling function ζ(ϕ) and the cosmological function Λ(ϕ).The simplest variant of the theory corresponds to a constant coupling function ζ(ϕ) = ζ = const and to the zero cosmological function Λ(ϕ) = 0.In the nongravitational part of the Lagrangian density, the gravitational field enters through the metric tensor g ik only, and hence, the theory obeys the Einstein equivalence principle.In usual scalar-tensor theories, instead of Λ g in ( 1), the scalar curvature R for the metric tensor g ik appears.These quantities differ by a total divergence: Because of the spacetime dependence of ϕ, the field equations following from (1) differ from those for scalar-tensor theories.
In the bimetric formulation of general relativity, the scalar Λ g appears as the Lagrangian density for the gravitational field instead of the scalar curvature R. Since these two scalars differ by the total divergence ∇ l wl , the corresponding field equations coincide.In particular, the second metric does not appear in these equations.However, the presence of the second metric tensor allows one to formulate the energy-momentum characteristics of the gravitational field in terms of tensorial objects (see, for example, [27]).
The field equations for the metric tensor and the scalar field, obtained from (1), have the form: where )/δg ik is the metric energy-momentum tensor for the nongravitational matter, T = g mn T mn ; the prime means the derivative with respect to the field ϕ, and the brackets in the index expression of the first equation mean the symmetrization with respect to the indices included.By using the first equation in (5), the equation for the scalar field can also be written in the form: From the equation for the nongravitational matter, δL m /δq a = 0, it follows that ∇ k T k i = 0. Note that, unlike the usual scalar-tensor theories, in BSTT, the latter equation does not follow from the field equations.This is a consequence of the presence of non-dynamical metric γ ik .

Spherically-Symmetric Static Configurations
In this section, we consider a static spherically-symmetric distribution of gravitating masses (for an investigation of superdense stellar configurations of degenerate gaseous masses in general relativity and in scalar-tensor theories, see, for example, [28]).In spherical spacetime coordinates (t, r, θ, φ) with the origin at the center of the configuration, we will take the metric tensors in the form: g ik = diag(e ν , −e λ , −r 2 e µ , −r 2 e µ sin 2 θ), (7) where ν, λ, µ are functions of the radial coordinate tending to zero at large distances.The energy-momentum tensor will be taken in the form: with ρ and P being the energy density and the pressure of the matter and u i is the corresponding four-velocity.We will consider the variant of the theory with ζ(ϕ) = ζ = const and Λ(ϕ) = 0.
For the metric tensors from (7), one has: where the prime stands for the derivative with respect to r.From the Equation (5) for the scalar field, we get: for 0 r < ∞.After integration, this gives: Assuming that ϕ (0) is finite, from here, we find µ(r) = λ(r).
With this relation, the remaining equations are presented in the form: where z = (ν + λ)/2 and: In addition to (12), the equation of state ρ = ρ(P) should be given.For nongravitational sources obeying the strong energy condition, one has ρ + 3P 0, and the function ν(r) is monotonically increasing.
Outside the mass distribution with the radius r 1 , r > r 1 , one has ρ = 0 = P.In this region, M(r) = M(r 1 ) ≡ M, β(r) = β(r 1 ) ≡ β and ϕ(r) → ϕ 0 = 1/(8πG) for r → ∞, with G being the Newton gravitational constant.The latter condition is required to ensure the limiting transition to the Newtonian gravity for weak fields.In the case ζ = 2, the corresponding solution for the function ν = ν(r) is given by: where l = r/r g , r g = 2GM and: For the remaining functions, one has: Note that the Tolmen formula for the mass M of the configuration, obtained in general relativity, is valid in BSTT, as well (see [29]).The external solution is determined by the mass M and by the parameter α.In the limit ζ → ∞, the solution given by ( 13) and ( 15) reduces to the Schwarzschild solution in general relativity described in isotropic coordinates: The horizon corresponds to the spherical surface r = r g /4.In the special case ζ = 2, the solution is specified to: This solution can also be obtained from ( 13) and ( 15) in the limit ζ → 2. Under the condition −ζ/2 α 2 < 1 − ζ/2, the external solution is regular for all r 1 .For the values of the parameter α outside this region, the external solution is singular at: Near the point l = l c , the leading terms in the metric functions and the scalar field behave as e ν ∝ (l At large distances from the spherically-symmetric configuration, one has: Similar to the post-Newtonian considerations, in the expansion for the g 00 component of the metric tensor, we have also kept the next-to-leading correction.As seen, the theory parameter enters into the corrections in the form of the ratio α/ζ.For α ζ (this is the case for nonrelativistic sources and for ζ 1), the variation in the scalar field is much smaller than that for the metric tensor components.Note that in the solar system, α 5 × 10 −6 .For this type of matter source, one has l c ≈ (1/4) 1 − 2/ζ, and the spherical surface l = l c is inside the Schwarzschild horizon in general relativity.
Similar to the case of Brans-Dicke theory, for large values of the parameter ζ, the results obtained within the framework of BSTT tend to those for general relativity.However, the constraints on the parameters of these theories are essentially different.For the Brans-Dicke theory, the observational data within the framework of the solar system give rise to the strong constraint (mainly from Cassini measurements of the Shapiro time delay) on the theory parameter, namely ζ > 4 × 10 4 .In BSTT, the restriction is much weaker, ζ v 2 , where v 2 ∼ P/ρ (P/ρ ∼ 10 −5 in the solar system) is the velocity in a post-Newtonian system with the pressure P and the energy density ρ.Under this condition, the post-Newtonian parameters of BSTT coincide with those in general relativity, whereas in the Brans-Dicke theory for the post-Newtonian parameter γ, one has γ = (1 + ζ)/(2 + ζ).This difference between the two scalar-tensor theories is related to the fact that in the Brans-Dicke theory the scalar field is sourced by the scalar curvature R (in the equation for the scalar field, that is the analog of the second equation in (5), R stands instead of Λ g ), which is of the order v 2 and, hence, the same for the variations of scalar field.In BSTT, the scalar field is sourced by Λ g , which has the order v 4 , and related to that, the variations of scalar field in post-Newtonian systems are of the order v 4 if the theory parameter is not too small (in the variations, the parameter ζ enters in the form 1/ζ (see, for example, (18))).

Models of Superdense Stellar Configurations
In this section, we consider models of superdense configurations based on the numerical integration of Equation (12).In superdense celestial bodies, except for their relatively thin surface layer, the real temperatures are always considerably lower than the degeneracy temperatures of the particles inside the configuration.Consequently, the temperature corrections to the equation of state are negligible.In the construction of models for degenerate stellar configurations, it is convenient to use the pressure as an independent variable and to consider the equation of state of the form ρ = ρ(P).
For a given value of the pressure at the center of the configuration, P(0), the quantities ν(0), λ(0), ϕ(0) and the mass M are determined by the matching conditions for the external and internal solutions on the surface of the gravitating body r = r 1 , where P(r 1 ) = 0.In Figure 1, we have plotted the dependence of the mass (in units of the solar mass M ) of a superdense configuration as a function of the central pressure for separate values of the theory parameter ζ (numbers near the curves).In the numerical integration, we have used the equation of state from [30,31].It is obtained from the relations (6.1) in [30] substituting the pressure and the chemical potential from the formulas (6.6), (6.7), (6.9) and (6.11) from that reference.The curves with ζ = ∞ correspond to configurations in general relativity.In general relativity, the configurations corresponding to the monotonically-increasing segment on the left of the first local maximum correspond to white dwarfs.The monotonically-increasing segment between the first and second local maxima corresponds to configurations that are unstable against radial perturbations.The monotonically-increasing segment between the second and third local maxima corresponds to neutron stars.Note that in BSTT, the issue of the stability of the static configurations requires an additional investigation.As seen from the graphs, depending on the value of the theory parameter ζ, in BSTT, the masses of the corresponding configurations for a fixed value of the pressure at the center can be essentially larger compared to the corresponding values in general relativity.For large values of the parameter, ζ 1, the curves in BSTT and in general relativity are practically indistinguishable.Note that the study of neutron star masses presently has attracted considerable interest both in the fields of theory and observations (see, for instance, [32,33]).This, in particular, was motivated by the observation of neutron stars with masses near 2M .
In addition to the mass of the configuration, the external geometry is determined by: In terms of this parameter, one has α = 3β/M.In Figure 2, we have plotted the ratio β/M as a function of the central pressure for the values of the theory parameter ζ = 0.1, 0.4, ∞.As we could expect for configurations corresponding to white dwarfs, β/M 1.Note that for the examples corresponding to Figure 2, one has a > 0, and the external solution is given by the first formula in (13).The numerical analysis shows that the difference of that solution from the solution on general relativity with the same value of the mass is rather small.The numerical integration shows that for ζ > 2, the parameters of the configurations in BSTT are very close to those for general relativity.Let us consider examples of configurations in the variant of BSTT with ζ = 0.4, corresponding to a neutron star with the mass M = 1.5M .For the equation of state we have used, one has r 1 ≈ 36 km and α ≈ 0.044.The numerical analysis shows that the difference of the corresponding external solution, given by the first formula in (13), from the solution of general relativity with the same value of the mass is rather small.In particular, the function ϕ(r) is a monotonically-increasing function outside the configuration and 1 − ϕ(r)/ϕ 0 < 0.001.This means that with the approach to the surface of the star in the exterior region, the effective gravitational constant increases.For this example, the external gravitational field is close to the one for a neutron star with the same mass in general relativity.The numerical integration inside the configuration shows that the function ϕ(r) monotonically increases in the interior region, as well, and the corresponding variation from the center to the surface is small, [ϕ(0) − ϕ(r 1 )] /ϕ 0 ≈ 0.0012.Hence, the effective gravitational constant takes its maximum value at the center of the configuration.The variations of the scalar field are relatively larger for larger values of the mass.For example, for the configuration with M = 3.8M , one has ϕ(0)/ϕ 0 ≈ 0.705 and [ϕ(0) − ϕ(r 1 )] /ϕ 0 ≈ 0.0087.As in the previous case, the function ϕ(r) is monotonically increasing for all values of the radial coordinate.For this configuration, the value of the parameter α is larger, α ≈ 0.135, and in accordance with the analysis presented before, the variations of the scalar field are more significant.

Conclusions
In the present paper, we have investigated static spherically-symmetric configurations of gravitating masses in the simplest version of BSTT with a constant coupling function and with the zero cosmological function.The constraints on the coupling constant, obtained from the observations within the framework of the post-Newtonian approximation, are essentially weaker than those in usual scalar-tensor theories (Brans-Dicke theory).The corresponding external solution is given by ( 13) and (15).In addition to the mass, the solution depends on the parameter α, which is determined in terms of the integral involving the pressure of the gravitating matter.For the integration of the field equations inside the configuration, we have used the equation of state from [30,31].The corresponding numerical data are presented in Figures 1-3.These results show that, for strong gravitational fields, depending on the value of the BSTT parameter ζ, the characteristics of the superdense configurations may essentially differ from those in general relativity.Note that we have considered spherically-symmetric solutions with a variable gravitational scalar.The possibility for solutions with a constant scalar field has been discussed in [34,35].

Figure 1 .
Figure 1.The mass of a superdense star as a function of the central pressure P(0).The numbers near the curves correspond to the values of the bimetric scalar-tensor theory (BSTT) parameter ζ and P 33 = P/10 33 erg • cm −3 .The curve ζ = ∞ corresponds to configurations in general relativity.

Figure 2 .
Figure 2. The dependence of the ratio β/M on central pressure.The numbers near the curves correspond to the values of the BSTT parameter ζ.

Figure 3
Figure 3 presents the radius of the configuration versus the central pressure for different values of the parameter ζ (numbers near the curves).The graphs show that, depending on the value of ζ, the radii of the configurations may notably exceed the corresponding values in general relativity.

Figure 3 .
Figure 3.The radius of a superdense star versus the central pressure.The numbers near the curves correspond to the values of the parameter ζ.The curve ζ = ∞ corresponds to the results in general relativity.