The search of new solutions for the Einstein-Maxwell field equations is an important area of research because it allows describe compact objects with strong gravitational fields as neutron stars, white dwarfs and quark stars [ 1,2]. Within this context it is appropriate to mention the findings of Delgaty and Lake [ 3] who constructed several analytic solutions that can describe realistic stellar configurations and satisfy all the necessary conditions to be physically acceptable. These exact solutions have also made it possible the way to study cosmic censorship and analyze the formation of naked singularities [ 4].

In the development of the first stellar models it is important to mention the pioneering research of Schwarzschild [ 5], Tolman [ 6], Oppenheimer and Volkoff [ 7] and Chandrasekhar [ 8]. Schwarzschild [ 5] obtained interior solutions that allows describing a star with uniform density, Tolman [ 6] generated new solutions for static spheres of fluid, Oppenheimer, and Volkoff [ 7] studied the gravitational equilibrium of neutron stars using Tolman’s solutions and Chandrasekhar [ 8] produced new models of white dwarfs in presence of relativistic effects.

A great number of exact models from the Einstein-Maxwell field equations have been generated by Gupta and Maurya [ 9], Kiess [ 10], Mafa Takisa and Maharaj [ 11], Malaver and Kasmaei [ 12], Malaver [ 13,14], Ivanov [ 15] and Sunzu et al [ 16]. For the construction of these models, several forms of equations of state can be considered [ 17]. Komathiraj and Maharaj [ 18], Malaver [ 19], Bombaci [ 20], Thirukkanesh and Maharaj [ 21], Dey et al. [ 22] and Usov [ 23] assume linear equation of state for quark stars. Feroze and Siddiqui [ 24] considered a quadratic equation of state for the matter distribution and specified particular forms for the gravitational potential and electric field intensity. Mafa Takisa and Maharaj [ 11] obtained new exact solutions to the Einstein-Maxwell system of equations with a polytropic equation of state. Thirukkanesh and Ragel [ 25] have obtained particular models of anisotropic fluids with polytropic equation of state which are consistent with the reported experimental observations. Malaver [ 26] generated new exact solutions to the Einstein-Maxwell system considering Van der Waals modified equation of state with polytropic exponent. Malaver and Kasmaei proposed a new model of compact star with charged anisotropic matter using a cosmological Chaplygin fluid [ 27]. Tello-Ortiz et al. [ 28] found an anisotropic fluid sphere solution of the Einstein-Maxwell field equations with a modified Chaplygin equation of state. More recently, Malaver et. al [ 29] obtained new solutions of Einstein’s field equations in a Buchdahl spacetime considering a nonlinear electromagnetic field.

The analysis of compact objects with anisotropic matter distribution is very important, because that the anisotropy plays a significant role in the studies of relativistic spheres of fluid [ 30,31,32,33,34,35,36,37,38,39,40,41,42]. Anisotropy is defined as $\Delta =p_t-p_r$ where$p_r$ is the radial pressure and $p_t$is the tangential pressure. The existence of solid core, presence of type 3A superfluid [ 43], magnetic field, phase transitions, a pion condensation and electric field [ 23] are most important reasonable facts that explain the presence of tangential pressures within a star. Many astrophysical objects as X-ray pulsar, Her X-1, 4U1820-30 and SAXJ1804.4-3658 have anisotropic pressures. Bowers and Liang [ 42] include in the equation of hydrostatic equilibrium the case of local anisotropy. Bhar et al. [ 44] have studied the behavior of relativistic objects with locally anisotropic matter distribution considering the Tolman VII form for the gravitational potential with a linear relation between the energy density and the radial pressure. Malaver [ 45,46], Feroze and Siddiqui [ 24,47] and Sunzu et al.[ 16] obtained solutions of the Einstein-Maxwell field equations for charged spherically symmetric spacetime by assuming anisotropic pressure.

In this paper we generated new classes of exact solutions for anisotropic charged distribution with a consistent with quark matter. New models have been obtained by specifying a particular form for one of the metric potentials and for the electric field intensity. The paper has been organized as follows: In section 2, we present the Einstein-Maxwell field equations. In section 3, we have chosen a particular form for the metric potential and for the electric field intensity. In Section.4, physical requirements for the new models are described. In section 5 we present the physical analysis of the new models. In section 6, we conclude.

We consider a spherically symmetric, static and homogeneous spacetime. In Schwarzschild coordinates the metric is given by:

where $\nu (r)$and a$\lambda (r)$are two arbitrary functions.

Using the transformations,$x=C{r}^2$, $Z(x)=e^{-\text{2λ}(r)}$ and $A^{\text{2}}{y}^{\text{2}}(x)=e^{\text{2ν}(r)}$ with arbitrary constants A and c>0, suggested by Durgapal and Bannerji [ 48], the Einstein-Maxwell field equations shown in [ 21] can be written as:

$\rho$ is the energy density, $p_r$ is the radial pressure, $E$ is electric field intensity,$p_t$ is the tangential pressure, σ is the charge density, $\Delta =p_t-p_r$ is the measure of anisotropy and dots denote differentiations with respect to x.

With the transformations of [ 48], the mass within a radius r of the sphere take the form

In this paper, we asume the following linear equation of state in the bag model

where B is the bag constant. We can write the Einstein-Maxwell field equations with the eq. (8) in the following form

The equations (9), (10), (11), (12), (13) governs the gravitational behavior of an anisotropic charged quark star.

In this research, we have chosen the Thirukanesh-Ragel-Malaver ansatz [ 25,37,49] as metric potential which has the form$Z(x)={\left(1-ax\right)}^n$, where a is a real constant and n is an adjustable parameter. This potential is regular at the stellar center and well behaved in the interior of the sphere. We have considered the particular cases n=1, 2. For the electric field intensity we make the choice

This electric field is finite at the centre of the star and remains continuous in the interior. For the case n=1, substituting Z(x) and eq. (14) in eq. (2) we obtain

Using eq. (15) in eq. (8), the radial pressure can be written in the form

and for the mass function we obtain

Substituting (14) and Z(x) in eq. (12) we have

Integrating eq. (18)

where $c_1$ is the constant within integration procedures.

For convenience we have let

For the metric functions $e^{2\lambda }$, $e^{2\nu }$

With eq. (14) and Z(x) in eq. (13), the charge density is

and the anisotropy F044 can be written as

With n=2, the expression for the energy density is

replacing eq. (25) in eq. (8), we have for the radial pressure

and the mass function is

Substituting eq. (14) and Z(x) in eq. (12) we obtain

Integrating eq. (28)

where $c_2$is the constant of integration

Again for convenience

The charge density can be written as

and for the metric functions$e^{2\lambda }$, $e^{2\nu }$ and anisotropy F044 we have

For a model to be physically acceptable, the following conditions should be satisfied [ 3,4,25]:

The metric potentials $e^{2\lambda }$and $e^{2\nu }$assume finite values throughout the stellar interior and are singularity-free at the center r=0.

The energy density ρ should be positive and a decreasing function inside the star.

The radial pressure also should be positive and a decreasing function of radial parameter.

The radial pressure and density gradients $\raisebox{1ex}{$d{p}_r$}\!\left/ \!\raisebox{-1ex}{$dr$}\right.\le 0$ and $\raisebox{1ex}{$d\rho $}\!\left/ \!\raisebox{-1ex}{$dr$}\right.\le 0$for $0\le r\le R$ .

The anisotropy is zero at the center r=0, i.e. Δ(r=0) =0.

The interior solution should match with the exterior of the Reissner-Nordstrom spacetime, for which the metric is given by

through the boundary r=R where M and Q are the total mass and the total charge of the star, respectively.

The conditions (ii) and (iv) imply that the energy density must reach a maximum at the centre and decreasing towards the surface of the sphere.

With n=1, the metric potentials $e^{2\lambda }$and $e^{2\nu }$have finite values and remain positive throughout the stellar interior. At the center $e^{2\lambda \left(0\right)}=1$ and $e^{2\nu \left(0\right)}=A^2{c}_1^2{\left(-1\right)}^{2D}$. We show that in r=0 ${\left(e^{2\lambda (r)}\right)}^{\prime }{}_{r=0}={\left(e^{2\nu (r)}\right)}^{\prime }{}_{r=0}=0$ and this makes is possible to verify that the gravitational potentials are regular at the center. The energy density and radial pressure are positive and well behaved between the center and the surface of the star. In the center $\rho (r=0)=6aC$ and$p_r(r=0)=aC-\frac{4B}{3}$, therefore the energy density will be non-negative in r=0 and $p_r(r=0)$ > 0. In the surface of the star $r=R$ and we have $p_r(r=R)=0$ and $R=\sqrt{\frac{4B-3aC}{\left((3a^2-a)C-4Ba\right)C}}$. For the radial pressure of density gradients we obtain

The energy density and radial pressure decrease from the centre to the surface of the star. From eq.(17), the mass function can be written as

and the total mass of the star is

On the boundary r=R, the solution must match the Reissner–Nordström exterior space–time as

and therefore, the continuity of $e^{2\lambda }$and $e^{2\nu }$across the boundary r=R is

Then for the matching conditions, we obtain:

For the case n=2, we have for the metric potentials $e^{2\lambda \left(0\right)}=1$, $e^{2\nu \left(0\right)}=A^2{c}_1^2{\left(-1\right)}^{2E}{e}^{-2F}$and ${\left(e^{2\lambda (r)}\right)}^{\prime }{}_{r=0}={\left(e^{2\nu (r)}\right)}^{\prime }{}_{r=0}=0$ at the centre r=0. Again the gravitational potentials are regular in the origin. The energy density and radial pressure also are positive and well behaved in the stellar interior. In the center $\rho (r=0)=2aC$ and $p_r(r=0)=2aC-\frac{4B}{3}$, therefore the energy density will be non-negative in r=0 and $p_r(r=0)$ > 0 . In the surface of the star $r=R$ and we have $p_r(r=R)=0$ and $R=\frac{\sqrt{30aC-20B}}{5aC}$ . For the radial pressure of density gradients we obtain

The energy density and radial pressure decrease from the centre to the surface of the star. For the mass function we have

The total mass of the star is

On the boundary r=R, the solution must match the Reissner–Nordström exterior space–time and therefore for the matching conditions, we obtain:

The figures 1, 2, 3, 4, 5, 6 and 7 present the dependence of $\rho$ , $\frac{d\rho }{dr}$ , $p_r$, $\frac{d{p}_r}{dr}$, M, Δ and ${\sigma }^2$with the radial coordinate respectively with a=0.2, B=0.05, C=1 for n=1 and n=2. We considered r=1.8 Km .

In theFigure 1 is shown that the energy density remains positive, continuous and is monotonically decreasing function throughout the stellar interior for all values of a. In theFigure 2 it is noted that for the radial variation of energy density gradient $\frac{d\rho }{dr}$< 0 in the two cases studied. The radial pressure showed the same behavior by the energy density, that is, it is growing within the star and vanishes at a greater radial distance and its results are shown inFigure 3. Again, according toFigure 4, the profile of $\frac{d{p}_r}{dr}$ shows that radial pressure gradient is negative inside the star for n=1 and n=2. InFigure 5, the mass function is continuous, strictly increasing and well behaved for all the cases. The anisotropic factor is plotted inFigure 6 and it shows that vanishes at the centre of the star, i.e. &#x00394;(r=0) =0 [ 30,37].Figure 7 shows that the charge density is regular at the centre, non-negative and grows with the radial parameter.

In this paper we have generated new models of anisotropic stars considering the Thirukkanesh-Ragel-Malaver ansatz for the gravitational potential and the MIT bag model equation of state. These models may be used in the description of compact objects in absence of charge and in the study of internal structure of strange quark stars. We show that the developed configuration obeys the physical conditions required for the physical viability of the stellar model. A graphical analysis shows that the radial pressure, energy density, mass function and anisotropy are regular at the origin and well behaved in the interior. The new solutions match smoothly with the Reissner&#x02013;Nordstr&#x000f6;m exterior metric at the boundary r=R because matter variables and the gravitational potentials of this work are consistent with the physical analysis of these stars. It is expected that the results of this research can contribute to modeling of relativistic compact objects and configurations with anisotropic matter distribution.