Charged anisotropic strange stars in Finslerian geometry

TL;DR

This paper models charged anisotropic strange stars within Finslerian geometry, deriving their structural properties, stability conditions, and maximum mass limits using a simplified approach with constant flag curvature and MIT bag model EOS.

Contribution

It introduces a novel Finslerian geometric framework for charged strange stars, incorporating anisotropy and charge distribution, and explores their equilibrium and stability characteristics.

Findings

Maximum charge of about 10^{20} C

Electric field around 10^{22} V/m

Mass-radius relations and stability conditions derived

Abstract

We investigate a simplified model for the strange stars in the framework of Finslerian spacetime geometry, composed of charged fluid. It is considered that the fluid consisting of three flavor quarks including a small amount of non-interacting electrons to maintain the chemical equilibrium and assumed that the fluid is compressible by nature. To obtain the simplified form of charged strange star we considered constant flag curvature. Based on geometry, we have developed the field equations within the localized charge distribution. We considered that the strange quarks distributed within the stellar system are compiled with the MIT bag model type of equation of state (EOS) and the charge distribution within the system follows a power law. We represent the exterior spacetime by the Finslerian Ressiner-Nordstr{\"o}m space-time. The maximum anisotropic stress is obtained at the surface of the system. Whether the system is in equilibrium or not, has been examined with respect to the Tolman-Oppenheimer-Volkoff (TOV) equation, Herrera cracking concept, different energy conditions and adiabatic index. We obtain that the total charge is of the order of 10$^{20}$ C and the corresponding electric field is of around 10$^{22}$ V/m. The central density and central pressure vary inversely with the charge. Varying the free parameter (charge constant) of the model, we find the generalized mass-radius variation of strange stars and determine the maximum limited mass with the corresponding radius. Furthermore, we also considered the variation of mass and radius against central density respectively.

Figures

Tables

Table 1. Table 1 : The physical parameters are predicted from the proposed model where 1 ​ M ⊙ = 1.475 ​ k ​ m 1 subscript 𝑀 direct-product 1.475 𝑘 𝑚 1~{}{{M}_{\odot}}=1.475~{}km for G = c = 1 𝐺 𝑐 1 G=c=1 . The set is valid for B g = 83 Mev/fm 3 , R ​ i ​ c ¯ ¯ 𝑅 𝑖 𝑐 \overline{Ric} = 1.2 and α 𝛼 \alpha = 0.0010 k ​ m − 2 𝑘 superscript 𝑚 2 km^{-2}

Strange	Observed	Predicted	${\rho }_c$	$P_c$	Q	E	Surface
Stars	Mass	Radius	($gm/c{m}^3$)	($dyne/c{m}^2$)	($Coulomb$)	(V/m)	redshift
	$\left(M_{\odot }\right)$	$km$	$\times {10}^{14}$	$\times {10}^{34}$	$\times {10}^{20}$	$\times {10}^{22}$
$PSRJ1614-2230$	$1.97\pm 0.04$[54]	${11.08}_{-0.07}^{+0.06}$	$7.703$	$5.314$	$1.578$	$1.156$	0.315
$CygX-2$	$1.78\pm 0.23$[55]	${10.75}_{-0.43}^{+0.39}$	$7.555$	$4.873$	$1.441$	$122$	0.284
$PSRB1913+16$	$1.44\pm 0.002$[56]	${10.09}_{-0.01}^{+0.01}$	$7.254$	$3.976$	$1.192$	$1.053$	0.231
$CenX-3$	$1.09\pm 0.08$[57]	${9.27}_{-0.22}^{+0.20}$	$6.933$	$3.022$	$0.924$	$0.967$	0.180
$LMCX-4$	$1.29\pm 0.05$[58]	${9.76}_{-0.12}^{+0.11}$	$7.116$	$3.565$	$1.079$	$1.018$	0.209
$4U1820-30$	$1.58\pm 0.06$[59]	${10.38}_{-0.13}^{+0.12}$	$7.535$	$4.276$	$1.298$	$1.083$	0.253
$4U1608-52$	$1.74\pm 0.14$[60]	${10.68}_{-0.26}^{+0.25}$	$7.509$	$4.736$	$1.413$	$1.115$	0.278
$EXO1785-248$	$1.30\pm 0.20$[61]	${9.78}_{-0.49}^{+0.44}$	$7.140$	$3.637$	$1.085$	$1.021$	0.253
$VelaX-1$	$1.77\pm 0.16$[62]	${10.73}_{-0.29}^{+0.28}$	$7.535$	$4.276$	$1.433$	$1.120$	0.211

Table 2. Table 2 : Numerical values of physical parameters for different charge constant α 𝛼 \alpha for the strange star L ​ M ​ C ​ X − 4 𝐿 𝑀 𝐶 𝑋 4 LMC~{}X-4 of mass 1.29 M ⊙ subscript 𝑀 direct-product M_{\odot} ( 1 ​ M ⊙ = 1.475 ​ k ​ m 1 subscript 𝑀 direct-product 1.475 𝑘 𝑚 1~{}{{M}_{\odot}}=1.475~{}km ) with B g = 83 Mev/fm 3 .

Value of $\alpha$	$\alpha$ = 0.0000	$\alpha$ = 0.0005	$\alpha$ = 0.0010	$\alpha$ = 0.0015
	$k{m}^{-2}$	$k{m}^{-2}$	$k{m}^{-2}$	$k{m}^{-2}$
Predicted Radius ($km$)	9.79	9.78	9.76	9.72
Central density $\left(gm/c{m}^3\right)$	$7.421\times {10}^{14}$	$7.359\times {10}^{14}$	$7.116\times {10}^{14}$	$6.754\times {10}^{14}$
Central pressure $\left(dyne/c{m}^2\right)$	$4.484\times {10}^{34}$	$4.290\times {10}^{34}$	$3.565\times {10}^{34}$	$2.488\times {10}^{34}$
Charge $\left(Coulomb\right)$	$0.000$	$0.543\times {10}^{20}$	$1.079\times {10}^{20}$	$1.598\times {10}^{20}$
Electric field ($V/m$)	$0.000$	$0.510\times {10}^{22}$	$1.016\times {10}^{22}$	$1.521\times {10}^{22}$
Surface redshift ($Z_s$)	0.216	0.214	0.209	0.201

Table 3. Table 3 : Numerical values of physical parameters for different R ​ i ​ c ¯ ¯ 𝑅 𝑖 𝑐 \overline{Ric} for the strange star L ​ M ​ C ​ X − 4 𝐿 𝑀 𝐶 𝑋 4 LMC~{}X-4 of mass 1.29 M ⊙ subscript 𝑀 direct-product M_{\odot} ( 1 ​ M ⊙ = 1.475 ​ k ​ m 1 subscript 𝑀 direct-product 1.475 𝑘 𝑚 1~{}{{M}_{\odot}}=1.475~{}km ) with B g = 83 Mev/fm 3 and α 𝛼 \alpha = 0.0010 k ​ m − 2 𝑘 superscript 𝑚 2 km^{-2} .

Value of $\overline{Ric}$	$\overline{Ric}$ = 1	$\overline{Ric}$ = 1.1	$\overline{Ric}$ = 1.2
Predicted Radius ($km$)	9.68	9.73	9.76
Central density $\left(gm/c{m}^3\right)$	$7.531\times {10}^{14}$	$7.372\times {10}^{14}$	$7.116\times {10}^{14}$
Central pressure $\left(dyne/c{m}^2\right)$	$4.812\times {10}^{34}$	$4.028\times {10}^{34}$	$3.565\times {10}^{34}$
Charge $\left(Coulomb\right)$	$1.052\times {10}^{20}$	$1.068\times {10}^{20}$	$1.079\times {10}^{20}$
Electric field ($V/m$)	$1.010\times {10}^{22}$	$1.015\times {10}^{22}$	$1.018\times {10}^{22}$
Surface redshif	0.274	0.237	0.209

Table 4. Table 4 : Numerical values of physical parameters for different Bag values for the strange star L ​ M ​ C ​ X − 4 𝐿 𝑀 𝐶 𝑋 4 LMC~{}X-4 of mass 1.29 M ⊙ subscript 𝑀 direct-product M_{\odot} ( 1 ​ M ⊙ = 1.475 ​ k ​ m 1 subscript 𝑀 direct-product 1.475 𝑘 𝑚 1~{}{{M}_{\odot}}=1.475~{}km ) with R ​ i ​ c ¯ ¯ 𝑅 𝑖 𝑐 \overline{Ric} = 1.2.

Value of Bg	Bg = 70 Mev/fm3	Bg = 80 Mev/fm3	Bg = 90 Mev/fm3
Predicted Radius($km$)	10.33	9.88	9.50
Central density $\left(gm/c{m}^3\right)$	$5.875\times {10}^{14}$	$6.833\times {10}^{14}$	$7.779\times {10}^{14}$
Central pressure $\left(dyne/c{m}^2\right)$	$2.631\times {10}^{34}$	$3.361\times {10}^{34}$	$4.055\times {10}^{34}$
Charge $\left(Coulomb\right)$	$1.279\times {10}^{20}$	$1.119\times {10}^{20}$	$\mathrm{1.0.995}\times {10}^{20}$
Electric field ($V/m$)	$1.078\times {10}^{22}$	$1.031\times {10}^{22}$	$0.992\times {10}^{22}$
Surface redshift	0.192	0.205	0.217

Equations

$$\frac{d^2{x}^{\mu }}{d{\tau }^2}+2G^{\mu }(x,\dot{x})=0,,$$

$$G^{\mu }=\frac{1}{4}{g}^{\mu \nu }\left(\frac{{\partial }^2{F}^2}{\partial x^{\lambda }\partial y^{\nu }}{y}^{\lambda }-\frac{\partial F^2}{\partial x^{\nu }}\right).$$

$$g_{\mu \nu }=\frac{\partial }{\partial y^{\mu }}\frac{\partial }{\partial y^{\nu }}\left(\frac{1}{2}{F}^2\right),$$

$$F^2=e^{\lambda (r)}{y}^t{y}^t-e^{\nu (r)}{y}^r{y}^r-r^2{\overline{F}}^2(\theta ,\varphi ,y^{\theta },y^{\varphi }).$$

$$g_{\mu \nu }=diag(e^{\lambda (r)},-e^{\nu (r)},-r^2{\overline{g}}_{ij}).$$

$$G^t=\frac{1}{2}{\lambda }^{\prime }{y}^t{y}^r,$$

$$G^r=\frac{1}{4}\left({\nu }^{\prime }{y}^r{y}^r+{\lambda }^{\prime }{e}^{\lambda -\nu }{y}^t{y}^t-2r{e}^{-\nu }{\overline{F}}^2\right),$$

$$G^{\theta }=\frac{1}{r}{y}^{\theta }{y}^r+{\overline{G}}^{\theta },$$

$$G^{\varphi }=\frac{1}{r}{y}^{\theta }{y}^r+{\overline{G}}^{\varphi }.$$

$$Ri{c}_{\mu \nu }=\frac{{\partial }^2}{\partial y^{\mu }\partial y^{\nu }}\left(\frac{1}{2}{F}^{2}Ric\right).$$

$$Ric=R_{\mu }^{\mu }=\frac{1}{F^2}\left[2\frac{\partial G^{\mu }}{\partial x^{\mu }}-y^{\lambda }\frac{{\partial }^2{G}^{\mu }}{\partial x^{\lambda }\partial y^{\mu }}+2G^{\lambda }\frac{{\partial }^2{G}^{\mu }}{\partial y^{\lambda }\partial y^{\mu }}-\frac{\partial G^{\mu }}{\partial y^{\lambda }}\frac{\partial G^{\lambda }}{\partial y^{\mu }}\right],$$

$$F^{2}Ric=\left[\frac{{\lambda }^{\prime \prime }}{2}+\frac{{\lambda }^{\prime 2}}{4}-\frac{{\lambda }^{\prime }{\nu }^{\prime }}{4}+\frac{{\lambda }^{\prime }}{r}\right]e^{(\lambda -\nu )}{y}^t{y}^t+\left[-\frac{{\lambda }^{\prime \prime }}{2}-\frac{{\lambda }^{\prime 2}}{4}+\frac{{\lambda }^{\prime }{\nu }^{\prime }}{4}+\frac{{\nu }^{\prime }}{r}\right]y^r{y}^r$$

$$+\left(\overline{Ric}-e^{-\nu }+\frac{r{\nu }^{\prime }{e}^{-\nu }}{2}-\frac{r{\lambda }^{\prime }{e}^{-\nu }}{2}\right).$$

$$G_{\nu }^{\mu }=g^{\mu \nu }Ri{c}_{\mu \nu }-\frac{1}{2}S,$$

$$S=e^{-\nu }\left[{\lambda }^{\prime \prime }+\frac{{\lambda }^{\prime 2}}{2}-\frac{{\lambda }^{\prime }{\nu }^{\prime }}{2}+\frac{{\lambda }^{\prime }-{\nu }^{\prime }}{r}\right]-\frac{2}{r^2}\left[\overline{Ric}-e^{-\nu }+\frac{r{e}^{-\nu }}{2}({\nu }^{\prime }-{\lambda }^{\prime })\right].$$

$$T_j^i=(\rho +p_t)u^i{u}_j-p_t{\delta }_j^i-(p_t-p_r)v^i{v}_j,$$

$$E_j^i=\frac{1}{4\pi }(-{\mathcal{F}}^{im}{\mathcal{F}}_{jm}+{\delta }_j^i{\mathcal{F}}^{mn}{\mathcal{F}}_{mn}),$$

$$G_{\nu }^{\mu }=8{\pi }_F(T_j^i+E_j^i)=8{\pi }_F{T}_{\nu }^{\mu },$$

$$T_{\nu ;\mu }^{\mu }=0.$$

$$\displaystyle T^{\mu}_{\nu}=\left(\begin{tabular}[]{cccc}$\rho+\frac{q^{2}}{8\pi r^{4}}$&0&0&0\\ 0&$-(p_{r}-\frac{q^{2}}{8\pi r^{4}})$&0&0\\ 0&0&$-(p_{t}+\frac{q^{2}}{8\pi r^{4}})$&0\\ 0&0&0&$-(p_{t}+\frac{q^{2}}{8\pi r^{4}})$\\ \end{tabular}\right),$$

$$\frac{{\nu }^{\prime }{e}^{-\nu }}{r}-\frac{e^{-\nu }}{r^2}+\frac{\overline{Ric}}{r^2}$$

$$\frac{{\lambda }^{\prime }{e}^{-\nu }}{r}+\frac{e^{-\nu }}{r^2}-\frac{\overline{Ric}}{r^2}$$

$$e^{-\nu }\left[\frac{{\lambda }^{\prime \prime }}{2}+\frac{{\lambda }^{\prime 2}}{4}-\frac{{\lambda }^{\prime }{\nu }^{\prime }}{4}+\frac{{\lambda }^{\prime }-{\nu }^{\prime }}{2r}\right]$$

$$\frac{dm}{dr}=4\pi r^2\rho +\frac{q}{r}\frac{dq}{dr}.$$

$$m(0)=0,q(0)=0,\rho (0)={\rho }_{c}and{p}_r(0)=p_c.$$

$$F^2=f_{RN}{y}^t{y}^t-f_{RN}^{-1}{y}^r{y}^r-r^2{\overline{F}}^2(\theta ,\varphi ,y^{\theta },y^{\varphi }),$$

$$q(r)=Q{\left(\frac{r}{R}\right)}^3=\alpha r^3,$$

$$\rho (r)={\rho }_c\left[1-\left(1-\frac{{\rho }_0}{{\rho }_c}\right)\frac{r^2}{R^2}\right],$$

$$p_r=\sum _{f=u,d,s}{p}^f-B,$$

$$\rho =\sum _{f=u,d,s}{\rho }^f+B,$$

$$p_r=\frac{1}{3}(\rho -4B).$$

$$p_r=\frac{1}{3}(\rho -{\rho }_0).$$

$$\nu \left(r\right)=-\ln \left(\frac{{\nu }_1{r}^4+{\nu }_2{r}^2+A{R}^5}{R^5}\right),$$

$$\displaystyle\lambda(r)=-\frac{32}{\lambda_{{1}}\nu_{1}}\Bigg{[}\lambda_{1}\,{\rm arctanh}\left({\frac{{R}^{5}{\alpha}^{2}+16\,B\pi\,{R}^{3}-M}{\lambda_{{2}}}}\right)+\lambda_{{2}}\left(\lambda_{3}-\lambda_{4}\right)\ln\lambda_{5}$$

Full text

1 Introduction

The theoretical study of the nature and properties of the strange quark stars is an attractive topic of research, not only due to a distinct branch of compact stars, but also the strange matter EOS is worthy of explaining a few astrophysical compact objects. The possibility of strange quark matter (SQM) [1, 2, 3], made up of equal unconfined up, down and strange quarks may be counted as the basic state of the strong interaction, proposed by Bodmer [4] and later on by Witten [5]. Theoretical presence of strange quark star was proposed by Itoh [6].

To provide the chemical equilibrate of the strange stars, a small number of electrons should be included in SQM. This small number of electrons (not bounded by any strong interaction) plays a significant role in the constitution of electric dipole layer at the surface. Therefore, a non-negligible electrical energy density comparable to radial pressure has developed a local non-neutral strong electric field within the star, i.e. between electron layer and positively charged core. Strange quark stars in the presence of the strong electric field can be modeled using the Maxwell-Einstein field equations. Strong impacts of this electric filed on the gravitationally bounded system are already shown by the authors [7, 8, 9]. The electric field on the surface region is in the order of $10^{19}\sim 10^{20}$ V/m. This surface electric field would be more extreme if the strange matter is made up of color superconducting strange matter [10]. The strange star demands more charge to be in stable equilibrium in a strong gravitational field. A significant amount of charge can induce an acute electric field. Non-zero charge modulate the structure of the strange star in a different manner: (i) Curvature of space-time, i.e. the metric, (ii) Energy density affiliated to the electric field enriches the total mass of the system, (iii) Coulombian interaction has a finite contribution to the hydrostatic equilibrium of the system, (iv) Charge also contributes to the anisotropic stress. It reduces the amount of anisotropic stress of the system. The charge and mass density contribute finitely to form an equilibrium configuration of the charged fluid.

Within last two decades, various charge distribution has been frequently applied to describe the effects of charge on the interior of strange stars [11, 12, 13]. Polytropic stars with charge density related to the energy density have been studied in [14, 15]. In their pioneer work, Negreiros et al. [11] investigated the strange star with the Gaussian charge distribution. They found that the charge gradient ($dQ/dr$) is not dependent on the width of charge distribution. This is true for any relativistic stellar object, provided that the distribution is narrowly spread over the system. Several authors worked with the charge distribution in the form of power law $q(r)=Q(r/R)^{n}$, where $Q$ and $R$ are total charge and radius of the system respectively. Felice et al. [16, 17], in his works studied the system for n$\geq$ 3 to make sure that the charge density does not diverge at the origin. Arbañil and Malheiro, in their work [13] studied for $n=3$ for perfect fluid system and later on, Deb et al. [18] studied the strange star for the same. Anninos and Rothman [19] studied the stellar system with a complex type of charge distribution.

Considering the homogeneous distribution of matter and constant surface charge, the stability of the charged star investigated in [20] and found that the incompressible star with no charge is less stable than a star with small surface charge and constant energy density. Glazer developed Chandrasekhar’s pulsation equation for the charged fluid system [21]. Stability of an incompressible fluid star can be increased by introducing charge [22]. The hydrostatic equilibrium and collapse of a charged fluid system studied in [23]. An extensive study of the equilibrium dependency and stability on the charge distribution is provided in [13]. They found that the stability of the strange star inversely varies with the total charge due to a certain range of the central energy density, whereas for a range of total mass stability increases with the increment of charge. Also, studied the hydrostatic equilibrium and stability with the radial perturbation for the varying central density, charge and charge-radius ratio.

In this article, we have studied the generalized structure of strange stars in the presence of electric charge, in the framework of Finsler geometry, and study the stability of the system. The reason for choosing Finsler geometry is that its length elements are not bounded by any quadratic restriction and the geometry depends on dynamics along with the position of the system [24]. Extensions of Einstein gravity in Fineslerain geometry is studied in several literature [25, 26, 27]. Nowadays several authors are using Finsler geometry to describe the violation of Lorentzian invariance and anisotropy of the Universe [28, 29, 30, 31, 32, 33, 34, 35]. Pfeifer and Wohlfarth [36] studied the casual structure and the generalized theory of electrodynamics in Finslerian space-time. However, the solutions to the field equations are not discussed. Li [37] has obtained a Finslerian Reissner-Nordström solution for vacuum space. He also found the eigenfunction of the Finslerian Laplacian operator. The generalized form of the Maxwell-Einstein Field equation in the following geometry is obtained from the geodesic equation of motion in the geometry. We have considered simplified MIT bag model EOS and the charge distribution in the form as assumed in [13] (power law). For simplicity, we constrained ourselves in constant flag curvature. Akbar-Zahed [38] already discussed the helpfulness of considering constant flag curvature and generality remain same. We obtained the exact solution of the Maxwell-Einstein Field equation. As a consequence, obtained maximum enclosed charge and field on the surface. The generalized variation of mass-radius relation for the strange stars due to a definite value of bag has been enumerated. We found a range of total mass, radius and charge respective of central density along with the theoretical bounds. The explicit study of the system are shown graphically and in tabular format for $\overline{Ric}=1.2$ and bag value $83~{}Mev/fm^{3}$, which is within the well accepted range.

Our paper is structured as follows: Definition and formation of field equation of Finsler geometry are outlined in section 2. $Adhoc$ relations are stated in section 3. We provided the formalism of basic stellar equations in section 4, we have presented the solution of the Maxwell-Einstein field equations. Physical acceptability and stability of the stellar system are verified in section 5 by studying mass-radius relation, energy conditions verification, the stability of the stellar model and compactification. Finally, the conclusion of our study with a discussion is provided in section 6.

2 BASIC STELLAR EQUATIONS

We briefly present fundamental geometrical concepts from the theory of Finsler spaces and generates the respective field equations, as well as discuss about the EOS of the stellar system.

2.1 Basic Formalism

We consider on a manifold $\mathpzc{M}$, the Finsler metric is $F$. In standard coordinate notation, $F=F(x,\dot{x})$ is the function of $(x^{\mu},y^{\mu})$ in $\in T\mathpzc{M}$.

The equation of Geodesic of the Finsler metric $(F)$ can be written as follows:

[TABLE]

where the geodesic spray is given by,

[TABLE]

The metric structure coefficient is given by

[TABLE]

with $(g^{\mu\nu})$ = $(g_{\mu\nu})^{-1}$.

Let, the Finsler structure be

[TABLE]

Using eq. (2.1), the Finsler metric potential can be defined as

[TABLE]

The respective geodesics sprays of the system are as follows:

[TABLE]

In Finsler geometry, Ricci tensor is introduced by Akbar-Zadeh [38], given by

[TABLE]

Ricci scalar in the Finsler geometry is,

[TABLE]

where, $R^{\mu}_{\mu}$ is insensitive to connections, it only depends on Finsler structure.

With the help of eqs. (2.2), (2.3) and (2.10), we obtain,

[TABLE]

The scalar curvature can define as $S=g^{\mu\nu}Ric_{\mu\nu}$.

Therefore, in Finsler geometry, the modified form of Einstein tensor reads,

[TABLE]

with the consideration of the flag curvature reads $\overline{F}^{2}=y^{\theta}y^{\theta}+f(\theta)y^{\phi}y^{\phi}$.

The explicit form of $S$ for the Finsler structure eq. 2.3 as follows:

[TABLE]

The energy momentum tensor for the anisotropic fluid distribution and electromagnetic field within the system, has the following forms respectively,

[TABLE]

where $\rho$, $p_{r}$ and $p_{t}$ represent the energy density, radial and tangential pressures, respectively.Here, $u_{i}$ and $v_{i}$ represent four-velocity and radial four-vector, respectively. $\mathcal{F}^{i}_{j}$ is the anti-symmetric electromagnetic field tensor.

Hence, in the Finsler geometry the Einstein-Maxwell field equations is

[TABLE]

where, we consider the geometrized unit, i.e., $G=1=c$.

The covariant divergence of the stress-energy tensor is

[TABLE]

The effective energy-momentum tensor for the locally anisotropic charged fluid distribution can be written as,

[TABLE]

where, the electric charge ($\mathbf{q}$) and the corresponding field ($\mathbf{E}$) is related as $\frac{E^{2}}{8\,\pi}=\frac{q^{2}}{8\,\pi\,r^{4}}$.

So, the components of the Einstein-Maxwell field equations are given by

[TABLE]

The effected gravitational mass equation due the spherically symmetric charged compact stellar object is defined by

[TABLE]

2.2 Junction Condition

The stellar structure extended from centre towards the surface of the system within the following constrains

[TABLE]

The interior spacetime of the stellar system should be matched smoothly at the boundary with the exterior spacetime. The Finslerian Ressiner-Nordström metric to represent the exterior spacetime of the following form [37]

[TABLE]

where $f_{RN}=\left(C-\frac{2GM}{R}+\frac{GQ^{2}}{r^{2}}\right)$ with $M$ and $Q$ being the total mass and charge of the system, respectively, and $C$ is a constant.

3 Charge distribution, Density profile and Equation of state

To study the effects of charge on the anisotropic strange stars, the charge distribution definitely should have a form. Considered the charge distribution is in the form of power law following Felice et al. [16] as $q(r)=Q(r/R)^{n}$. For the simplicity, we assume that $n=3$ as follows,

[TABLE]

where $Q$ is the total charge and $R$ is the radius of the system, respectively. $\alpha\left(=\frac{Q}{R^{3}}\right)$ is a constant (charge constant).

We consider that the energy density of the fluid inside the strange stars maintaining the form [39]

[TABLE]

where $\rho_{c}$ and $\rho_{0}$ are the central and surface density, respectively.

The SQM within the system is described by the phenomenal MIT bag model. We assumed that the quarks are massless and non-interacting (Here we considered up($u~{}$), down($d~{}$) and strange ($s~{}$) quarks). The corrected form of pressure can be defined as after introduction of $adhoc$ bag function (B) as follows,

[TABLE]

where $p^{f}$ is the pressure of each type of quarks, i.e. u, d and s. The corresponding corrected energy density is as follow,

[TABLE]

Now, substituting the relation between The pressure and energy density due to the each quark flavor given by $p^{f}=\frac{1}{3}\rho^{f}$ and eq. (3.4) into eq. (3.3) we have the final form of the MIT bag EOS as follows

[TABLE]

The radial pressure ($p_{r}$) vanishes on the surface. Therefore, we can consider the surface density ($\rho_{0}$) as 4B. Following that the eqs. (3.5) can be rewritten as,

[TABLE]

4 Solution of the Einstein-Maxwell Field Equation

We obtained the following expression for the gravitational potentials ($\nu$ and $\lambda$), density, radial and tangential pressures, respectively, by solving the eqs. (2.23)-(2.25) with the help of eqs. (3.1)-(3.5), given by

[TABLE]

[TABLE]

where $\nu_{1},\nu_{2},\lambda_{1},\lambda_{2},\lambda_{3},\lambda_{4}$ and $\lambda_{5}$ are constants and their expressions are provided at the Appendix.

The variation of the gravitational potentials, viz. $\textrm{e}^{\nu}$ and $\textrm{e}^{\lambda}$ as a function of the fractional radial coordinates (r/R) at the interior of the stellar system are exhibited in the left and right panel of figure 1, respectively. The variations of the physical quantities like $\rho$, $p_{r}$ and $p_{t}$ are shown in the left, middle and right panel of figure 2, respectively.

The anisotropic stress ($\Delta$) of a stellar system can be expressed as additional tangential pressure over radial direction, i.e. $p_{t}-p_{r}$, is given by

[TABLE]

The variation of the anisotropic stress of the stellar system as a function of fractional radial coordinates ($r/R$) is shown in figure 3.

The variation of the electrical charge distribution (q(r)) and respective electrical energy density ($E^{2}/8\,\pi$) are exhibited in figure 4 in the left and right panel respectively as a function of fractional radial coordinate. It is clear from the plots that there are no charge density and corresponding field at the centre of the stellar distribution.

5 Physical features of the Stellar System

In this section, we are going to verify whether our model is in the stable equilibrium and physically valid.

5.1 physically acceptable behaviour

5.1.1 Adiabatic Index

The stiffness of the EOS of a stellar system for a defined density variation can be characterized from the adiabatic index ($\Gamma$). For an infinitesimal radial perturbation, the dynamical stability has been studied by several authors [40, 41, 42, 43, 44]. In the following work [45] it is shown that the adiabatic index must be greater than $4/3$ which can be defined as

[TABLE]

The variation of the adiabatic index as a function of the fractional radial coordinate is shown in figure 5.

5.1.2 Energy Conditions

Energy conditions depict the observer’s measurement of the matter distribution in the space-time. The conditions are always positive, which define that the flow of matter should be along null or time-like world line. Stavrinos and Alexiou [46] provided the required energy conditions for the Finslerian system. A stellar system is said to be a physically valid system if the following inequalities are simultaneously satisfied:

[TABLE]

Here NEC, WEC, SEC and DEC denote the null energy condition, weak energy condition, strong energy condition and dominant energy condition, respectively. The variation of the different energy conditions with the fractional radial coordinate due to different values of $\alpha$ is shown in figure 6.

5.1.3 Mass-Radius Relation

According to [47], the mass-radius ratio for a stable stellar system should maintain $2m(r)/r\leq 1$ throughout the region. Andréasson [48] generalized the maximum mass-radius ratio for a charged stellar system in the following form

[TABLE]

For the present stellar system the mass variation can be written as

[TABLE]

The variation of total mass respect to total radius is exhibited in figure 7. The variation is drawn for different charge constant. As a result from variation, we obtain that the mass increases with the charge constant. Stellar mass is normalized in solar mass ($M_{\odot}$).

5.1.4 Compactification Factor and Red Shift

The compactification factor of a stellar system can express as,

[TABLE]

The redshift function ($Z_{s}$) of a stellar system is defined as,

[TABLE]

The variation of $Z_{s}$ is shown in the figure 8 as function of radial coordinate.

5.2 Verification of the Stability

To justify the equilibrium configuration, we considered (i) TOV equation and (ii) Herrera cracking condition.

5.2.1 TOV Equation

A stellar system is said to be in equilibrium, if the resultant of forces on a system is nil. Imbalance in force drives the system to an unstable configuration. For the stellar system these counterbalancing forces are given by Tolman [49] and Oppenheimer-Volkoff [50].

According to TOV equation

[TABLE]

Here, the terms are defined as hydrostatic force ($F_{h}$), gravitational force ($F_{g}$), electrostatic force ($F_{e}$) and the last term defines the anisotropic force ($F_{a}$) respectively of the eqs. 5.5. To maintain the equilibrium, the outward forces $F_{e}$, $F_{a}$ and $F_{h}$ must balance by the attractive pull $F_{g}$.

The variation of the forces for different charge distributions are exhibit in figure 9.

5.2.2 Herrera cracking Condition

The physical stability of a stellar system can be verified with respect to the causality condition also. According to the condition, the square of the radial ($v_{sr}^{2}=\frac{dp_{r}}{d\rho}$) and the tangential ($v_{st}^{2}=\frac{dp_{t}}{d\rho}$) sound speed lies between 0 $\rightarrow$ 1, i.e. 0 $\leq v_{si}^{2}\leq 1$ (where, $i=r,t$). The region is called potentially stable if the radial sound speed is greater than the tangential sound speed.

According to the following works [51, 52] the difference of the squares of tangential and radial sound speeds must hold the sign inside the stellar system. Hence, according to Herrera’s condition $\mid v_{st}^{2}-v_{sr}^{2}\mid\leq$ 1.

In figure 10 the variation of the radial and tangential sound speeds and their differences are shown in the left panel and right panel, respectively.

6 Discussion and conclusion

In the paper at hand, we have performed a detail investigation of the physical impacts of the charge distribution on the structure and behaviour of the strange stars. We have generalized the description of the charged strange stars in Finslerian background by studying the modified form of the Maxwell-Einstein field equation. We derived the gravitational potentials from the family of field equations and smoothly matched them with exterior Finslerian Ressiner-Nordström solution. Figure 1 indicates the variation of the metric potentials viz., $e^{\nu(r)}$ and $e^{\lambda(r)}$, which are monotonically increasing and geometrically non-singular by nature. Interior quarks and electrons distribution profile of the system is related to the radial pressure of the system by the MIT bag equation. The density profile of the system decreases from the central density $\rho_{c}$ (A.8) holding the same sign of the slope. The radius of the stellar body is predicted when, there is no overlying matter against the gravitational attractive force, which endorse zero radial pressure. Due to the reason, the surface density is constant, this is one of the limitations of the approach. The radial and tangential pressure also decreases monotonically with the radial variation. The expression and variation of density, radial and tangential pressure are provided in eqs. (4.3)-(4.5) and in figure 2, respectively. As a consequence of using aniotropic fluid, our model shows that the anisotropic stress is maximum at surface area and there is no anisotropy at the centre of the system, as predicted by Deb et al. [53] in Riemannian frame. Due to the introduction of charge, the anisotropy reduces with increment of the charge constant, i.e., charge can remove the anisotropy of a system. The expression is provided in eq. (4.5) and the figure 3. The variation of charge distribution and the corresponding field is portrayed in figure 6. With the considered assumption, the charge and the following electric field reaches its maximum at the surface. The acceptability of the system is examined on the basis of the energy conditions, Herrera cracking condition, TOV equation and mass-radius relation.

The anisotropic flow along with the pressure gradient and coulomb repulsion is supported by gravitational pull of the matter (mass), inwards to it. It defines the overlying matter density reduces radially. The variation of the forces endorses that our stellar structure is non-varying in terms of the equilibrium of forces.

Cracking concept describes the nature of the deviation of the system from equilibrium. The notion depends on the theory of gravitation not on the geometry. We found that our stellar structure is consistent with both conditions (i) the causality relation and (ii) the Herrera cracking concept.

The red shift is monotonically decreasing with radial variation toward the surface. Variation of the redshift as a function of fractional radial coordinate ($r/R$) for the strange star $LMC~{}X-4$ is portrait in figure 8, where the compactification factor defines the mass boundness of the structure.

We have predicted the radii of few candidates of strange stars by using their observed masses. In addition we studied different physical parameters related to the structure for the predicted radius. All studies are provided in the tabular form, Ref. Table 1.

The maximum limiting mass and correlated radius gradually increase with the increasing values of $\alpha$ ref. figure 7 (the solid circles in the variation of the total mass M of the the stellar structure corresponds to maximum mass limit of the model for a given value of $\overline{Ric}=1.2$ and $B_{g}=83\,MeV/fm{{}^{3}}$). We found that the maximum mass ($M_{max}$) for $\alpha=0.0015~{}km^{-2}$ is higher than $M_{max}$ for zero charge by $4.77\%$ and the respective radius $R_{max}$ is increased by 0.98%. A stellar configuration is a stable or unstable equilibrium, can be categorized from the relation $\frac{dM}{d\rho_{c}}$. Positive slope indicates stable configuration. Maximum mass point for $\alpha=0.0015~{}km^{-2}$ is increased by $7.6\%$ than the maximum mass point for zero charge respective of $\rho_{c}$ variation. The variations of $M$ (normalized in $M_{\odot}$) and $R$ with respect to $\rho_{c}$ are shown in the figures 11 and 12, respectively. In addition, the variation of total charge ($Q$) with respect to central density is shown in figure 13, which features that with the increasing values of $Q$ the central density of the strange stars decreases gradually. Again, for the charge constant $\alpha=0.0015~{}km^{-2}$ the maximum total charge is obtained at $\rho_{c}=963~{}Mev/fm^{3}$, whereas the maximum mass $M_{max}=5.03M_{\odot}$ is obtained for $\rho_{c}=965~{}Mev/fm^{3}$. From the above numerical study, it is interesting to mention that, in the presence of charge, the amount of total mass increases more intensely than the radius, which indicates that the more compact systems can be justified in the framework.

To have a better understanding of the present model in Tables 2 and 3 we present a comparative study of the different physical parameters, viz., radius, central density, central pressure, surface redshift, total charge of the structure and the corresponding field, etc., due to the chosen parametric values of $\alpha$ and $\overline{Ric}$, respectively. Further, we have predicted values of the above mentioned physical parameters due to the chosen values of bag constant, viz. $B=70,80$ and $90~{}MeV/fm^{3}$ in Table 4. Interestingly, the numerical analysis reveals that with the increasing values of $\overline{Ric}$ it is possible to pack more mass in the stellar system and the total charge of the system also increases gradually.

There are also few limitations in our approach, due to which simulation is not completely analogues with the real life system, formally applicable to define a static stellar structure. Within all constraints, it is worthy of noting that in Finsler geometry the ultra-high dense strange stars can be successfully represented, which leaves the very interesting area for future research.

Acknowledgments

SR and FR are thankful to the Inter-University Centre for Astronomy and Astrophysics (IUCAA), Pune, India for providing Visiting Associateship under which a part of this work was carried out. SR is also thankful to the authority of The Institute of Mathematical Sciences, Chennai, India and the Centre for Theoretical Studies, IIT Kharagpur, India for providing short term visits under which a part of this work was carried out. FR is also grateful to DST-SERB (EMR/2016/000193), Government of India for providing financial support. A part of this work was completed while SRC and DD were visiting IUCAA and the authors gratefully acknowledge the warm hospitality and facilities there.

Appendix A Appendix

Values of the different constants viz., $\nu_{1},\nu_{2},\lambda_{1},\lambda_{2},\lambda_{3},\lambda_{4}$ and $\lambda_{5}$ are given by

[TABLE]