Compact stars in the non-minimally coupled electromagnetic fields to gravity

We investigate the gravitational models with the non-minimal Y(R)F2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y(R)F^2$$\end{document} coupled electromagnetic fields to gravity, in order to describe charged compact stars, where Y(R) denotes a function of the Ricci curvature scalar R and F2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F^2$$\end{document} denotes the Maxwell invariant term. We determine two parameter family of exact spherically symmetric static solutions and the corresponding non-minimal model without assuming any relation between energy density of matter and pressure. We give the mass-radius, electric charge-radius ratios and surface gravitational redshift which are obtained by the boundary conditions. We reach a wide range of possibilities for the parameters k and α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} in these solutions. Lastly we show that the models can describe the compact stars even in the more simple case α=3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =3$$\end{document}.


Introduction
Spherically symmetric solutions in gravity are fundamental tools in order to describe the structure and physical properties of compact stars. There is a large number of interior exact spherically symmetric solutions of Einstein's theory of gravitation (for reviews see [1,2]). But, very few of them satisfy the necessary physical and continuity conditions for a compact fluid. Some of them were given by Mak and Harko [3][4][5] for an isotropic neutral spherically symmetric matter distribution.
A charged compact star may be more stable [6] and prevent the gravitationally collapse [7,8], therefore it is interesting to consider the case with charge. The charged solutions of Einstein-Maxwell field equations which describe a strange quark star were found by Mak and Harko [9] considering a symmetry of conformal motions with MIT bag model.
One can consider that a charged astrophysical object can described by the minimal coupling between the gravitational and electromagnetic fields known as Einstein-Maxwell theory. However, the above problems of Einstein's gravity at large scales can also lead to investigate the Y (R)F 2 type modification of Einstein-Maxwell theory [20][21][22][23][24][25][26][27][28][29][30]. Such nonminimal modifications can also be found in  to obtain more information on the interaction between electromagnetic and gravitational fields and all other energy forms. The nonminimal couplings also can arise in such compact objects as black holes, quark stars and neutron stars which have very high density gravitational and electromagnetic fields [42]. If the extreme situations are disappeared, that is far from the compact stars the model turns out to be the Einstein-Maxwell case.
Here we consider the non-minimal Y (R)F 2 type modification to the Einstein-Maxwell theory and generalize the exact solutions for radiation fluid case k = 1 in [42] to the cases with k = 1, inspired by the study [9]. We obtain the inner region solutions and construct the corresponding model which turns to the Einstein-Maxwell theory in the outer region. We note that the inner solutions recover the solution obtained by Misner and Zapolsky [43] for charge-less case and b = 0. We find the surface gravitational redshift, matter mass, total mass and charge in terms of boundary radius and the parameters k and α via the continuity and boundary conditions.
We organize the paper as follows: In Sect. 2, we give the non-minimal Y (R)F 2 gravity model and field equations in order to describe a compact fluid. In Sect. 3, we obtain exact static, spherically symmetric solutions of the model under the conformal symmetry and the corresponding Y (R) function. In Sect. 4, we determine the total mass, charge and gravitational redshift of the star in terms of boundary radius and the parameters of the model k and α. We summarize the results in the last section.

The model for a compact star
Compact stars have very intense energy density, pressure and gravitational fields. They also can have very high electric fields in order to balance the huge gravitational pulling [6][7][8][44][45][46][47][48][49]. Even if they collapse, very high electric fields necessary to explain the formation of electrically charged black holes [50,51]. Moreover they can have very strong magnetic fields [52]. Therefore the new non-minimal interactions between these electromagnetic and gravitational fields with Y (R)F 2 type can arise under the extreme situations. Thus we consider the following non-minimal model for compact stars, which involves the Lagrangian of the electromagnetic source A ∧ J and the matter part L mat in the action [42] I (e a , ω a b , Here e a is the orthonormal co-frame 1-form, ω a b is the Levi-Civita connection 1-form obtaining by the relation T a = de a + ω a b ∧ e b = 0, F is the electromagnetic tensor 2-form which is derived from the electromagnetic potential A via the exterior derivative, that is F = d A; λ a is the Lagrange multiplier which gives torsion-less space-times, T a = 0; R is the Ricci scalar and J is the electromagnetic current density in the star. The co-frame variation of the action (1) is given by the following gravitational field equation after eliminating the connection variation, where we have used the velocity 1-form u = u a e a for time-like inertial observer and dY d R = Y R . Furthermore we have considered the matter energy momentum tensor has the energy density ρ and the pressure p as diagonal elements for the isotropic matter in the star. It is worth to note that the total energy-momentum tensor which is right hand side of the gravitational field Eq. (2) satisfies the conservation relation [42]. The electromagnetic potential variation of the action gives the following modified Maxwell field equation We have also the identity d F = d(d A) = 0. We will find solutions to the field equations (2) and (3) under the condition which eliminates the instabilities of the higher order derivatives in the theory. Here we note that the constant k 1 determines the strength of the non-minimal coupling between gravitational and electromagnetic fields. The case with k = 0 leads to Y (R) = constant and this case corresponds to minimal Einstein-Maxwell theory which can be considered as the exterior vacuum Reissner-Nordström solution with R = 0. The additional features of the constraint (4) can be found in [42]. On the other hand, the trace of the non-minimally coupled gravitational field equation (2) gives Since the case with k = 1 or ρ = 3 p is investigated as the radiation fluid stars in [42], we concentrate on the case with k = 1 or ρ = 3 p in this study.

Spherically symmetric solutions under conformal symmetry
We take the following static, spherically symmetric metric and Maxwell 2-form with the electric component E which has only radial dependence Then the charge in the star can be obtained from the integral of the current density 3-form J over the volume V with radius r using the Maxwell equation (3) The Ricci scalar for the metric (6) is calculated as For k = 1 or ρ = 3 p the gravitational field equation (2) gives under the condition (4) The conservation of the total energy-momentum tensor for the gravitational field equation (2) requires that Assuming the metric (6) has the conformal symmetry L ξ g ab = φ(r )g ab which describe the interior gravitational field of stars [9,[53][54][55], the metric functions f 2 (r ) and g 2 (r ) were found in [53] as Here φ(r ) is an arbitrary function of r , L ξ is Lie derivative of the metric tensor along the vector field ξ , X is a new function, φ 0 and a are arbitrary constants. Then we obtain the following differential equation system from (10)- (14) under the symmetry, with the condition (13). Since we aim to extend the solution given in [42] to compact stars without introducing an equation of state, we take the following metric function with the real numbers α > 2 and b = 0 inspired by [9], We see that the metric function is regular at center of coordinate r = 0 and leads to the following regular Ricci scalar Then we found the following class of solutions to the system of equation (16)-(20) where C is an integration constant. We see that the solutions are dependent on the parameter k introduced by (13) together with the parameter α in the metric function (20). The charge of the star (8) inside a spherical volume with radius r can be calculated by the charge-radius equality It can be seen that the charge is regular at r = 0 for α > 2. By solving r in terms of R as the inverse function r (R) from (21), we re-express the non-minimal coupling function (25) as We can consider that the exterior vacuum region is described by Reissner-Nordström metric with R = 0. Then the nonminimal function becomes Y (R) = C and our model involves the Einstein-Maxwell theory as a minimal case at the exterior region with C = 1. Thus we have the more general Lagrangian of the non-minimal theory which describe the compact stars for α > 2 The field equations of the Lagrangian (28) accept the solutions with the electric field (22), pressure (23), energy density (24), electric charge (26) and the following metric tensor in the star We will determine the power α and k in the model from observational data and the constant b from the matching and continuity conditions (33). In the absence of the electromagnetic source and matter, the Lagrangian of the non-minimally with the exterior electric field E(r ) = Q r 2 , where Q is total electric charge of the star.

Continuity conditions
The continuity of the interior (29) and exterior Reissner-Nordström metric at the boundary of the charged star r = r b leads to At the boundary of the star the pressure (23) should be zero the condition determines the constant b in the non-minimal theory (28) as Using this constant, variation of the pressure and energy density as a function of the radial distance r inside the star is given in Fig. 1. The interior of the star is considered as the specific fluid which has very high gravitational fields, electromagnetic fields and matter. Then the electromagnetic fields obey the modified Maxwell field equation d * Y F = J in matter. The integral of the modified Maxwell equation (3) gives the charge in volume with radius r in the star, q(r ) = Y Er 2 (26). On the other hand, the Ricci scalar is zero at the exterior of the star and then, using (27), we can take the non-minimal function as Y = 1 and obtain the Maxwell field equation d * F = 0 which leads to the total charge Q = Er 2 at the exterior region. Here the displacement field Y E inside of the star turns out to be the electric field E outside. Setting by r = r b in Eq. (26) the total charge of the star Q = q(r b ) is found as Then the outside electric field is given by E = Q/r 2 . The total charge-boundary radius ratio obtained from (35) is plotted by Fig. 2a dependent on the parameter α for some different k values.
The figure of the matter mass M m is shown in Fig. 2b for increasing α and some different k values. We see that the total charge and matter mass of the star increase with the increasing k values. When we compare the constant b given by the Eqs. (34) and (32), we find the following mass-charge relation for the model By substituting the total charge (35) in (37), the total massradius ratio M/r b can be found in terms of the parameters α and k of the model . (38) The gravitational redshift z at the boundary is obtained from By taking the limit α, k → ∞, the maximum redshift is found as z = √ 3 − 1 ≈ 0.732 same with k = 1 in [42]. The upper redshift bound is smaller than the Buchdahl bound z = 2 and the bound given in [56].
The mass-radius relation and the gravitational redshiftradius relation are shown in Fig. 3a, b, respectively, depending on the parameter α for some k values. We see that as the k value increases, the mass and redshift increases.
To obtain an interval for the parameter k we consider the energy density condition inside the star At the center of the star the condition gives k ≥ −1. On the other hand, at the boundary r = r b the condition (40) turns out to be The solution of the inequality (42) is Then we can choose k ≥ 1 without loss of generality. The derivative of the pressure p(r ) according to the energy density ρ(r ) is calculated as The phase speed of the sound waves in the star is defined by ( dp dρ ) 1/2 and the speed satisfies the causality condition ( dp dρ ) 1/2 < 1 for k ≥ 1 and α > 2 values (see Fig. 6a for some values), where the speed of light c = 1. Thus each possible values of the parameters k > 1 and α > 2 in the modified model give a mass, charge-radius ratio and redshift z.

The simple model with α = 3
We consider the simple case setting by α = 3 in the nonminimally coupled model (28) Here we note that the non-minimal coupling function Y (R) can be expanded as the binomial series in power of R 3 for The non-minimal model (45) accept the interior metric as a solution ds 2 = −a 2 r 2 dt 2 + 3 1 + br 3 dr 2 + r 2 (dθ 2 + sin 2 θ dφ 2 ) together with the following electric field, pressure, energy density and electric charge from (22)-(26) inside of the star Then under the boundary and matching conditions the total charge, mass and gravitational redshift can be expressed by the followings Fig. 4 The dimensionless quantities which are related with the matter mass (a) and total electric charge (b) versus the dimensionless parameter k for α = 3 from (35)- (39). We note that the redshift has the upper bound z ≈ 0.633 for α = 3, which is smaller than the general Fig. 6 Variation of the phase speed v = dp dρ as a function of radial distance r using the boundary radius r b = 10 km and some k values for α = 3 (a) and α = 5 (b) Table 1 The dimensionless parameter k, the charge-radius ratio κ 2 Q 2 r 2 b and the surface redshift z obtained by using the observational mass M and the radius r b for some neutron stars redshift bound of the model which is √ 3 − 1 ≈ 0.732. As the simple model with α = 3, we depict the related physical quantities in Figs. 4, 5 and 6. We give the corresponding k values by taking the observed mass-radius ratios of some known stars and calculate the other quantities in Table 1.
Here we note that each star have its own k value and the each k value can be determined by the observational mass and radius of the star.
In order to obtain an approximate equation of state of the matter inside the star from Eqs. (49) and (50), we fit the parametric p − ρ curve with the equation where we have used B 0 = 1 km −2 as a dimension-full constant. We note that if k = 1 we obtain the radiation fluid case p = c 2 ρ/3 [42]. The curve fitting is shown in Fig. 7.
It is interesting to see that the fitting equation of state corresponds to the MIT bag model [61,62] which is given by p =  Table 1 for realistic compact stars. In the case α = 3, the non-minimal Einstein-Maxwell model with k = 1.015 gives the gravitational mass M = 0.65M using the radius r b = 9.46 km. We see that this mass is less than the mass obtained for conformal symmetric charged stellar models in the minimal Einstein-Maxwell theory [9]. We note that the both models have the conformal symmetry and very similar interior metric solutions. Although this modified Einstein-Maxwell model has free parameters which can be set to be consistent with various compact star observations, the model in [9] determines a unique charged configuration of quark matter in terms of the Bag constant. Alternatively, the mass in our model increases to 2.09M for α = 4, and 2.64M for α = 5. That is, we can confine k as k = 1.015 which gives the Bag constant and consider α is a free parameter, in order to describe compact stars.

Conclusion
We have extended the solutions of the previous the nonminimally coupled Y (R)F 2 theory [42] to the case with ρ = 3 p under the symmetry of conformal motions. In the case without any assumption of the equation of state, we have acquired one more parameter k in the solutions and the corresponding model. The pressure and energy density in the solutions are decreasing function with r in the interior of the star. By matching the interior solution with the exterior Reissner-Nordström solution and applying the zero pressure condition at the boundary radius of star r = r b , we determine some physical properties of the star such as the ratio of the total mass and charge to boundary radius r b and gravitational redshift depending on the parameters k and α.
We note that the total mass and charge increase with increasing k values and we have not reach an upper bound for k in the extended non-minimal model. But the increasing k values give an upper bound for the gravitational redshift, z ≈ 0.732, which is smaller than the more general restriction found in [56] for compact charged objects. It is interesting to note that each α and k value in this model (28) determines a different non-minimally coupled theory and each theory with the different parameters gives different mass-radius, chargeradius ratios and gravitational redshift configuration. We calculated k values and the corresponding other quantities of the compact stars with the simple case α = 3 via the some observed mass-radius values in Table 1. In this case, we also obtained the approximate equation of state (55) by fitting the p − ρ curve of the model. By comparing the fitting equation of state with the MIT bag model for α = 3, we found the gravitational mass M = 0.65M which is smaller than the mass obtained for conformal symmetric quark stars [9] in the Einstein-Maxwell theory. However the mass in our model increases as α increases. The non-minimally coupled model has the arbitrary parameters α and k which can be set in order to be consistent with compact star observations. Even for α = 3, each k value can describe a charged compact star.