Electrically charged supermassive twin stars

By assuming that ultra dense hybrid neutron stars are endowed with a distribution of electric charge, we study the corresponding twin star solutions and their properties resulting from a sharp first order transition from confined hadronic to a deconfined quark phase. Two distinct quark matter equations of state with increasing stiffness are considered and the values for the maximum gravitational masses of the hadronic and hybrid twin configurations are obtained for different values of the total electric charge. Interestingly, our calculations indicate that sharp transitions make charged twin hybrid stars more massive than their neutral counterparts, and that the {$\SI{2}{M_{\odot}}$ constraint from PSR J0740+6620} is surpassed for standard values of electric charge and can be considered stable only satisfying $\partial M/\partial \epsilon_0>0$. In particular, our charged stellar models reach masses even higher than the unknown compact object measured in the GW190814 event.


I. INTRODUCTION
Compact stars are an unique laboratory to study the strong interactions theory in the ultra dense regime, beyond the nuclear regime, where a phase transition between the hadronic and deconfined quarks degrees of freedom is expected to occur (for a recent review see, e.g. Ref. [1]). If such a phase transition is sharp, a third family of compact stars can be present in the mass-radius diagram, with smaller radii but similar masses to those predicted for the disconnected second family branch [2][3][4][5][6][7][8][9]. The description of these twin star configurations has been a theme of debate in recent years [10][11][12][13][14][15][16][17][18][19][20][21], mainly motivated by the perspective that future X-ray astronomy will be able to measure the masses and radii of compact stars with a high precision. The studies performed, e.g., in Refs. [11,12,20,21] indicate that the properties of the twin stars is strongly sensitive to the treatment of the phase transition in ultra dense matter. As a consequence, a future discovery of twin stars will be important to improve our understanding about the equation of state (EoS) of compact objects.
One of the main constraints in the EoS comes from the mass of the most massive pulsars observed (≈ 2 M ), since a viable EoS has to be able to generate compact stars with a maximum mass higher than this constraint. In recent years, additional constraints have been obtained by the analysis of gravitational wave data for the coalescence of compact binary systems. In particular, such data has been used to obtain information about the sources, such as the progenitor masses (M prog1,2 ). For instance, the gravitational wave event GW170817 [22], apart from marking the beginning of the era of multimessenger astronomy [23,24], consistently allowed to identify the progenitors [22] as being neutron stars (NS) having M prog,1 ∈ [1.36, 1.60]M and M prog,2 ∈ [1.17, 1.36]M . The interpretation of these compact stars as being com-posed by hadronic and/or quark degrees of freedom is still a theme of debate. On the other hand, gravitational wave events without multimessenger counterpart and having one progenitor with 2.2 < M prog,1 /M < 5 are considerably more difficult to interpret 1 . For instance, the low-mass progenitor (M prog,1 ∈ [1.61, 2.52]M ) of GW190425 [26] is likely a NS, but not conclusively. Another example is the event GW190814 [27], since its lightest progenitor is an unknown object having M prog,1 ∈ [2.50, 2.67]M , possibly being a very light BH or a supermassive NS. Interestingly, if assumed to be a NS, it puts severe constraints on a large set of EoSs for dense matter 2 . Very recently, the analysis of these supermassive compact stars has been considered to constrain twin stars, with the results derived in Ref. [20] pointing out that a mass constraint of 2.5 M ruled out all twin star solutions, while Ref. [21] concluded that these ultra heavy NSs can be in a twin-mass configuration but the two stable branches must be connected.
In this paper, we extend these previous studies by considering the impact of the electric charge on the structure properties of hybrid compact stars, with particular emphasis on the ultra-dense branch, i.e. charged twin stars. We stress that although our study is strongly motivated by the higher masses reached by charged stellar configurations, specially in the view of the possibility that the unknown object in GW190814 to be a charged hybrid star, our main aim is to investigate the generic effects of electric charge on the structure and stability of hybrid stars with sharp phase transitions. For the partic-1 Notice that the recent events GW200105 and GW200115 [25] do not belong to this class since M prog,1 M prog,2 and are easily characterized as NS-BH mergers. 2 Many standard pure quark and nuclear-matter EoSs are able to reach this mass limit although not satisfying gravitational wave constraint for the radius of 1.4 M , see e.g. Refs. [28,29]. ular case of the GW190814 unknown object, the electric charge might be acquired during the coalescence process. In fact, several authors have discussed in recent years the possibility that during the gravitational collapse or the coalescence process of a compact binary system, matter can acquire a large amount of electric charge and, consequently, the structure equations must be solved taking into account an electric charge distribution inside the star. In particular, in Ref. [30], two of the authors have studied the properties of charged strange stars considering an EoS based on perturbative QCD and demonstrated that the presence of an electric charge distribution implies that the maximum mass is larger in comparison to the neutral counterpart, in agreement with the results obtained in Ref. [32] using a simplified approach (see also Refs. [31,[33][34][35]). Electrically charged hybrid stars were previously studied in Ref. [31]. However, these authors considered a smooth phase transition with a mixed phase, which does not generate twin stars. In contrast, in this work, we will consider the presence of an electric charge distribution within the compact hybrid star and assume that a sharp first order phase transition occurs.The impact of the electric charge will be investigated for two distinct phenomenological treatments of the phase transition, which give rise to twin star solutions, and the structural properties of the charged compact stars will be obtained. Besides, we will see that our findings indicate that the presence of an electric charge distribution implies that twin stars with masses larger than 2.6 M are stable, i.e. the associated configurations satisfy the classic criteria ∂M/∂ 0 > 0, which is the stability condition for sharp first order transitions characterized by a rapid conversion between phases [36]. However, notice that this can only be considered an application of our generic results but not the main finding which deals with different classifications for the transitions. This work is organized as follows. In Sec. II, we summarize the main aspects of the stellar structure equations for charged compact stars. Moreover, the EoSs considered in the description of the hybrid NS will be reviewed. In Sec. III, we present our results for the structural properties of the charged compact stars and a comparison with the neutral case will be performed. The twin star solutions will be presented for distinct values of the electric charge and different treatments of the phase transition. Finally, in Sec. IV we will summarize our main results and conclusions.

II. FRAMEWORK
In this section, we will present a brief review of the theoretical framework necessary to investigate electrically charged compact stars (see Refs. [30-32, 34, 35] for extensive details). Initially, we will present the electricallycharged structure equations and the model for the electric charge distribution used in this work. In the second part, the two treatments for the phase transition will be dis-cussed and the associated EoSs considered are presented. Finally, the predictions for the mass-radius diagram of electrically neutral stars will be shown for completeness.

A. Charged stellar structure equations
In the last years, several authors have discussed the basic structure equations needed to derive the properties of a charged compact star [30-32, 34, 35]. In our analysis, we will assume a static and spherically symmetric charged star. The stress-energy tensor will be expressed in terms of the contributions associated to the perfect fluid and to the electromagnetic field. For the case considered, this field is fully specified by the radial component of the electric field, which can be expressed in terms of the electric charge distribution. Following the approach presented in detail in, e.g. Ref. [30], one finds that the Einstein-Maxwell field equations imply that the charged stellar structure is determined by the following system of equations where ρ e (r), q(r) and m(r) are the electric charge density, charge and mass profiles, respectively, and Eq. (3) is the Tolman-Oppenheimer-Volkov (TOV) equation for a charged star. The metric potential e −2λ has the Reissner-Nordström form (see Refs. [30][31][32] for more details). The associated boundary conditions are established as follows. At the center of the star we have that the charge and mass are null q(0) = m(0) = 0 and we stipulate a central pressure p(0) = p 0 . Moreover, the stellar surface (r = R) is defined by the point where the pressure is null, i.e., p(R) = 0. At this point, we have the star's total gravitational mass m(R) = M and charge q(R) = Q, respectively. Moreover, in order to solve the structure equations, we must specify the charge distribution and the EoS that describes the matter inside the star, the latter being discussed in the next subsection. Regarding the charge distribution, as e.g. in Ref. [30], we will assume it is proportional to the energy density, i.e., ρ e = α , where, in geometric units, α is a dimensionless proportionality constant which can be considered a free parameter in our modeling. Such distribution will be denoted α-distribution in what follows. In this distribution, for a constant α one has that Q = αM . Requiring that the compact star does not collapse into a BH, i. e. that the exterior event horizon of the Reissner-Nordström spacetime has to be smaller than the stellar radius, which implies that Q < M , one should have that α < 1. As assumed in Refs. [30,31,37,38], we will also consider that the micro-physical effects of the electric charge on the EoS are negligible.

B. Equations of state
As demonstrated in previous studies, as e.g. in Refs. [11,12,20,21], twin star configurations arise when a sharp phase transition of the QCD matter occurs. In our analysis, we will assume that hadronic matter (HM) is described by the stiffest EoS of Hebeler et al. [39] which is obtained from effective chiral field theory, used between baryon densities of n B ∈ [0.6, 1.1]n 0 , and being supplemented by the BPS EoS [40] for the outer crust with n B < 0.6n 0 . For the quark phase, i.e. n B > 1.1n 0 , we consider two distinct approaches, which differ in the behaviour assumption for the speed of sound. In the first model, denoted CSS in what follows, we assume a constant speed of sound. In contrast, in second one, denoted MP, we consider that the speed of sound is dependent on the density and the associated equation of state can be described using a multipolytrope approach.
Following Refs. [10,12,20], we will consider that the entire EoS for the CSS model is given by where p t is the transitional pressure and HM (p t ) is the energy density for the point of transition. Moreover, ∆ is the discontinuity of the energy density at the transition and c QM is the speed of sound of quark matter (QM). In our analysis, we will assume c 2 QM = 1 in order to achieve the stiffest possible EoS. The resulting predictions are dependent on the values assumes for ∆ and p t , with the twin star solutions being generated for some combinations of these parameters. The studies performed in Refs. [12,20] have demonstrated that the maximum mass of the second branch is determined by ∆ , while the maximum mass of the first branch is influenced by p t . Moreover, these authors have classified the twin star solutions in four categories, depending on the masses of the twin stars, which are defined by: (a) Category I: the maximum masses of both stars are larger than 2 M ; (b) Category II: only the first maximum is larger than 2 M ; (c) Category III: the first maximum is in the range 1 ≤ M Had /M ≤ 2, while the second one is larger than 2 M ; and (d) Category IV: the first maximum is smaller than 1 M , while the second one is larger than 2 M . In Fig. 1 (left panel), we present the typical EoSs for these four categories, with the (∆ , p t ) values being given in the parenthesis. For completeness, we also present, in the right panel, the resulting predictions for the massradius relations derived by solving the TOV equation for a neutral compact star. The pure hadronic case (black thin curve) obtained from the Hebeler et al. stiffest EoS is presented for comparison. For the case where a phase transition is present, the maximum hadronic and hybrid masses are marked by filled circles. Moreover, the dotted lines indicate the unstable configurations. As it is clear from the figure, two separate branches, characteristic of twin stars, are present.
In our analysis, we also consider an EoS for QM with a speed of sound dependent on the density. Inspired by the study performed in Ref. [11], we will assume that high-density region can be described by a polytrope EoS given by P (n) = κn Γ , where n is the nuclear density and Γ is the polytrope parameter. As in Ref. [11], we will consider four different sets of parameters for κ and Γ, with the maximum values of Γ being determined by imposing that the speed of sound does not exceed the speed of light for the values of density reached in the center of the star. As in the CSS case, the hadronic phase will be described by the stiffest EoS of Hebeler et al. [39]. For all sets of parameters for the QM EoS, we assume that the critical pressure is p t = 63.177 MeV fm −3 , that the discontinuity of the energy density at the transition is ∆ = 253.89 MeV fm −3 and that the transitional baryon number density is 2.03 n 0 , in agreement with the previous estimation of these values 3 . The EoSs for the MP model are presented in Fig. 2 (left panel), while the corresponding predictions for the mass-radius relation are shown in the right panel, which were derived assuming a neutral compact star. One has that they share the same hadronic branch, and that the maximum mass of the second one increases for larger values of Γ. In particular, one has that for Set 1, the second branch is unstable and twin star solutions are not obtained. On the other hand, a third family of stable stars is present for the other sets, in agreement with the results obtained in Ref. [11].

III. RESULTS
In what follows, we will analyze the impact of the electric charge on the results derived in the previous section for the neutral twin stars predicted by the CSS and MP models. In general terms, one has for both models that the presence of an electric charge distribution implies stellar configurations that have a larger mass and radius compared to its neutral counterparts. However, depending on the amount of the electric charge, the new configurations become unstable and twin stars solutions are not present. Another important aspect for the case of stable charged twin stars configurations, is that the gap mass between the two branches is sensitive to the value of α, which establishes the relation between the electric charge distribution and the energy density (ρ e = α ).  [39] is presented. The maximum hadronic and hybrid masses are marked by filled circles. We also present the mass-radius constraints from the GW170817 event [22] considering Universal Relations (blue region) and a spectral EOS (yellow region which looks green when overlapped). Moreover, we add the NICER constraints [41,42] for PSR J0030+0451 (red regions), the updated constraint in mass [43,44] and radius [45,46] of the most massive pulsar PSR J0740+6620 (green region) and the mass of the low-mass companion of GW190814 [27] (brown horizontal band). Finally, the general relativity, finite pressure and causality limits [47] are also presented in different shades of grey. EoSs for the MP model considering four sets of values for the κ and Γ parameters (in parenthesis) that characterize quark matter in the high density regime. As in Ref. [11], for Set 4, we describe the high density regime in terms of two polytrope branches, depending on the nuclear density. Right panel: Corresponding predictions for the mass-radius relation of a neutral compact star. The maximum hadronic and hybrid masses are marked by filled circles. The same constraints and limits of Fig. 1 are presented.
able to predict stable charged twin stars with masses larger than 2.5 M in contrast to charged hybrid stars with Gibbs transitions like in Ref. [31]. Initially, in Fig. 3 we present our results for the distinct categories of the CSS model and different values of α. The neutral solution is also presented for comparison. The filled circles represent the maximum mass configurations for the hadronic and hybrid cases. We will focus on the values of α ≥ 0.25, where we expect a larger impact of the electric charge on the stellar configuration [30].
As a consequence, we will not show the results for twin stars configurations classified as Category II, since we have verified that they do not present a hybrid branch for α ≥ 0.2. The results presented in Fig. 3 demonstrate that the presence of charge increases the masses of the twin stars in both branches of all categories considered. However, we have found that Categories I, II and III do not present twin stable configurations for large values of α. In particular, for Category I, the hybrid branch is unstable for α > 0.25, even though the hadronic mass reaches the hypothetical constraint of GW190814 in the charged case. Moreover, for α = 0.25, the hybrid branch becomes smaller in comparison to the neutral counterpart. On the other hand, our results indicate that the hybrid branch becomes unstable for α > 0.8 for Category III, while Category IV has stable twin star solutions for the largest value of α allowed (close to unity), but that do not satisfy any of the previous mass-radius constraints.
A few comments are in order.
Regarding the GW170817 mass-radius constraint, it is only satisfied for Category IV twin stars. The remaining constraints are in general satisfied in every other category for larger values of alpha (when mentioned) for stable charged twin stars. In particular, the NICER constraints for PSR J0030+0451 at 68% confidence interval are satisfied by all of our 1.4 M charged stellar models which represent hadronic stars (except for Category IV, where 1.4 M stars are hybrid ones).
In Table I, we present our predictions for the maximum gravitational masses of the hadronic (Had) and hybrid (Hyb) configurations as well as the corresponding radii for the distinct categories of the CSS model. We have also included the radius of 1.4 M stars for comparison with the NICER results discussed above. In agreement with the previous discussions, the maximum masses and radii increase for larger values of α. In particular, for Categories III and IV, the hybrid branches for α = 0.5 clearly fulfill the hypothetical GW190814 mass constraint and is easily surpassed for larger values of α (not shown). However, the hadronic maximum mass of category IV barely surpasses 1 M even in the charged case. An interesting aspect is that the radius difference of the twin star solutions for category IV becomes larger than 4.5 km for α ≥ 0.5.
In order to verify the robustness of our findings, we also have studied the stellar configurations assuming that the quark EoSs are characterized by c 2 QM = 0.33 and 0.7. In both scenarios, we have found that the charged hybrid branch is completely unstable for Categories I and II for all values of α. Category III has unstable hybrid configurations for c 2 QM = 0.33, but has a stable one for c 2 QM = 0.7. In the latter, only the hadronic configurations surpass the 2.6 M limit. For Category IV, where the phase transition pressure is small, we have found a stable hybrid branch for both speeds of sound. Finally, M ≈ 2.6 M is reached for c 2 QM = 0.33 (0.7) when α ≥ 0.8 (0.6).
Now we pass to estimate the impact of the electric charge distribution on the predictions of the MP model for the different sets. In the previous section, we have shown that stable twin star solutions are not obtained for neutral stars when Set 1 is considered. One has verified that such conclusion is also valid for the charged case. Moreover, we have obtained that for Set 2 the stable twin star configurations only occur for small values of the electric charge (α ≤ 0.1). In contrast, such configurations are present in the Sets 3 and 4 for α ≤ 0.6 and 0.7, respectively. The corresponding predictions for the mass-radius relations of these sets are presented in Fig. 4.
The results for the neutral counterparts are shown for comparison. Our results indicate that for the Set 3, the GW190814 hypothetical mass constraint is satisfied by the hadronic branch for α = 0.5, with the associated hybrid branch barely reaching the lower limit. For α = 0.6, both branches satisfy the constraint. On the other hand, for larger values of α, all hybrid branches are unstable.
In the case of Set 4, we have verified that the GW190814 hypothetical constraint is fulfilled for α = 0.5, and stable charged twin star configurations, with larger masses and radii, are also present for α ≤ 0.7. In contrast, we have obtained that for α > 0.7 there are no stable hybrid branches.
In Table II, we present our predictions for the maximum gravitational masses of the hadronic (Had) and hybrid (Hyb) configurations as well as the corresponding radii for the Sets 3 and 4 of the MP model. We have not included the radius of 1.4 M stars in Table II because both Sets share the same hadronic configurations as the CSS model. So the R 1.4 are the same as the ones presented in Table I for categories I and III. Similarly to what was observed for the CSS model, the electric charge also implies in larger maximum masses and radii for larger values of α. An interesting aspect present in our predictions for Set 4 is that the impact of the electric charge on the hadronic and hybrid branches are distinct and dependent on the value of α, with increasing of the hadronic maximum mass being faster than the hybrid one. One has that for α < 0.5 the hybrid maximum mass is larger than the hadronic one, but for α ≥ 0.5 the hadronic maximum mass is already larger than the hybrid one. Finally, we have also calculated the baryon mass and compactness of charged twin stars considering the Set 4 and found that the charged stellar models also possess a negative binding energy, i.e. the baryon mass is larger than the gravitational one. Moreover, one has verified that the compactness of our charged twin star models are far below the 'charged Buchdahl limit' presented in Ref. [50] and, in fact, are even below 1/2 corresponding to Schwarzschild black holes. Both of these criteria serve to further verify the overall stability of our twin star configurations.

IV. SUMMARY
The description of compact stars have been subject of intense attention in recent years, mainly motivated by the beginning of the multi-messenger era in astrophysics, which has provided new constraints on the bulk properties of matter inside these objects and high precision data for their masses and radii is becoming available. Considering this perspective, the nuclear community is currently challenged to improve the modelling of complex processes, as for example, possible phase transitions in the core matter and the coalescence of binary compact objects. Our goal in this paper was to contribute to this effort by investigating, for the first time, the impact of an electric charge distribution on the stellar properties and stability of twin hybrid neutron stars. In particular, we have extended for the charged case the studies performed in Refs. [11,20], which have determined limi- tations on the masses of stable twin stars configurations assuming distinct assumptions for the phase transition which we emulate in order to establish our findings within a well-known framework. In particular, our calculations indicate that in general, the presence of large values of electric charge imply twin star configurations with larger masses than the maximum mass estimates made from gravitational-wave data [51,52] and charged hybrid stars with Gibbs transitions [31]. We have also verified that twin stars with masses of ≈ 2.6 M are stable in the lines of Ref. [36], which allow us to interpret the unknown compact object measured in the GW190814 event as a charged NS instead of a very light BH. However, such interpretation is still a theme of debate and should only be considered an interesting coincidence with our generic findings. Finally, our results motivate us to extend our analysis on the formation of stable charged twin hybrid stars which, according to Ref. [53], are difficult to occur in the zero-charge limit.