## Impact of dark energy on the structure of neutron stars: The vacuum case

L. F. Ara´ ujo, 1, ∗ J. A. S. Lima, 1, † and G. Lugones 2, ‡

1

Departamento de Astronomia, Universidade de S˜ ao Paulo Rua do Mat˜ ao, 1226 05508-900, S˜ ao Paulo, SP, Brazil 2 Universidade Federal do ABC, Centro de Ciˆ encias Naturais e Humanas, Avenida dos Estados 5001- Bang´, CEP 09210-580, Santo Andr´ e, SP, Brazil. u

The potential role of a cosmic vacuum dark component in the properties of neutron stars is investigated. It is assumed that the static, spherically symmetric distribution of matter within neutron stars is supported by two distinct components: ordinary matter and a vacuum fluid. For normal matter we use a set of state-of-the-art nuclear matter equations of state, each grounded in nuclear physics experiments. The vacuum energy component is inhomogeneously distributed within the star and obeys the standard equation of state ( p v = -ϵ v ). This is characterized by an energy density fraction y = ϵ m / ϵ ( m + ϵ v ), which we model as either a constant or radius-dependent. Our findings reveal that the inclusion of vacuum energy significantly affects the mass-radius relationships in neutron stars, influencing both the maximum achievable masses and the qualitative form of these relationships. Some constraints from current multimessenger observational data limiting the amount of vacuum energy within neutron stars are also discussed.

Keywords: compact objects, vacuum energy, dark energy, mass-radius relation

## I. INTRODUCTION

Numerous astrophysical and cosmological studies suggest that most of the observed Universe is made up of dark matter and dark energy, both of which have unknown origins and compositions. The clustering of dark matter helps explain the observed rotation of galaxies and the Universe's large-scale structure, while a uniform dark energy component is required to account for the Universe's late-time acceleration. The standard explanation for dark matter assumes the existence of yet unknown particles weakly coupled to the Standard Model. Meanwhile, it remains uncertain whether dark energy represents the manifestation of a novel field, a challenge to General Relativity thereby leading to alternative gravity theories, or an entirely different phenomenon.

While the role of dark energy at large scales has been extensively discussed, our understanding of its effects at more local scales, particularly at the level of individual astronomical objects such as neutron stars, remains much less explored. Over the years, some models have been proposed to understand the presence of dark energy in compact objects [1]. In the 1960s, Gliner put forward the idea of objects with a p = -ϵ equation of state (EOS) at their centers [2]. The Bardeen spacetime, proposed as the first concrete solution, represented a nonsingular, asymptotically flat, spherically symmetric spacetime with zero, one, or two event horizons, depending on a specific parameter [3-6]. In the 1980s, the formation and gravitational effects of false vacuum bubbles in true vacuum (and vice versa) were studied, with some analysis suggesting their potential role in wormhole formation and localized inflation [7-10]. The 2000s witnessed a resurgence of interest in the p = -ϵ concept. Chapline et al. proposed 'dark energy stars', gravitationally stable compact objects with a dark energy core and a microscopic quantum critical layer in place of an event horizon [11]. Another model, the stiff shell gravastar, was introduced by Mazur and Mottola, characterized by a positive pressure stiff matter surface layer connected to the dark energy core and external vacuum by junction layers [12, 13]. Research has extensively examined these models for their potential behaviors and stability [14-17]. Additional models, including non-singular black holes [18], have also been explored to understand dark energy in compact objects, as well as in the gravitational collapse in its variant associated with the vacuum energy density, minimally-coupled scalar fields, and in anisotropic regimes [19-23].

∗ loreanyfa@usp.br

† jas.lima@iag.usp.br

‡ german.lugones@ufabc.edu.br

In contrast to the more theoretical exploration of dark energy compact objects, there exists a significant body of work in the literature detailing the possible role of dark matter within realistic astrophysical objects. Numerous theoretical investigations have proposed that dark matter can become part of neutron stars either through direct gravitational capture or through non-gravitational interactions, resulting in so-called 'admixed' stars that contain a mixture of standard model matter and dark matter [24-32]. In these admixed neutron stars, it is suggested that dark matter may exist as a small core [31], or as a halo around the star that does not affect its visible radius [32]. In some extreme cases, neutron stars might be almost entirely composed of dark matter, creating highly compact dark objects with very low masses [33]. Several studies have also contemplated the existence of neutron stars composed of particles from a hidden dark-QCD sector, asymmetric dark matter fermions, or asymmetric bosons [34-39]. Many observational consequences have been hypothesized as a result of the existence of dark

matter within neutron stars. It has been suggested that dark matter could influence the mass-radius relation [40], affect the tidal deformability [41], and potentially modify the gravitational wave signal emitted during a binary neutron star merger [42].

In parallel with the gravitational capture of dark matter by neutron stars, the potential influence of dark energy in compact objects, given their strong gravitational fields, is an intriguing hypothesis that merits further exploration. While prior studies have essentially addressed gravitationally stable compact entities composed partially or entirely of pure dark energy, the coexistence of dark energy with standard model matter in such contexts remains unexplored. Our investigation adopts a perspective akin to 'admixed' dark matter stars, proposing the interweaving of dark energy and ordinary matter within the star. This 'admixed' framework gains enhanced astrophysical plausibility, representing a more feasible scenario.

Among several possibilities for describing the dark energy candidates, the simplest and most compelling one is the cosmological constant, Λ. In the context of the Einstein field equations, the Λ-term behaves like an additional energy component with a constant energy density, ϵ v = Λ 8 / πG , and pressure, p v = -ϵ v . When combined with the cold dark matter (CDM) component, the resulting model is commonly referred to as the cosmic concordance cosmology (ΛCDM). Over the last two decades, it has become the standard cosmological model. However, the ΛCDM cosmology is currently facing observational challenges, particularly from the socalled Supernovae-CMB tension [43]. Briefly, estimates of the Hubble constant based on supernovae observations from the SHOES collaboration, measurements at intermediate redshifts from various tests [44], and observations of strong lensing time-delays suggest a best-fit value of H 0 = 74km s Mpc, whereas the latest Planck / / CMB data favors a value closer to 68 km / / s Mpc. This indicates a discrepancy of around 5 6 . σ (see, e.g., [45] for a recent review discussing the main problems with the ΛCDM model).

The cosmological constant and coincidence problems have motivated the idea that Λ in the homogeneous and isotropic case would be a time-dependent quantity, i.e., Λ = Λ( ) due to its interaction with all the matter fields t in the Universe. In particular, it has been shown that ˙ Λ( ) t &lt; 0, suggesting that Λ is small today because the Universe is too old [46]. This result is also required by the second law of thermodynamics because entropy growth ( ˙ S &gt; 0) occurs only when the vacuum energy density decreases during expansion [47, 48]. In principle, this explains the current low value of the vacuum energy density and, potentially, the coincidence problem [49]. Naturally, in an inhomogeneous static case, such as in compact objects, the cosmological constant would be a function only of position, Λ = Λ( r ).

In this context, prompted by cosmology, it is of interest to investigate the implications of spatially varying vac- uum energy density on compact objects such as neutron stars. In cosmology, we have seen that the strength of Λ is greater at early times when the energy density of the Universe is very high. Similarly, one may conjecture that for neutron stars, the vacuum energy density diminishes from their core to their surface, as occurs for the pressure within the star itself.

On the other hand, recent developments in multimessenger astronomy have significantly advanced our understanding of neutron star properties. Pulsars such as PSR J1614-2230 [50, 51] and PSR J0348+0432 [52], each with masses around 2 M ⊙ , have provided crucial data for studying neutron star matter. The subsequent discovery of PSR J0740+6620 [53, 54], another high-mass pulsar, reinforces the need for equations of state to be capable of accommodating neutron stars with masses exceeding 2 M ⊙ . Collaborative efforts between the NICER telescope and XMM-Newton have yielded valuable constraints on the mass and radius of pulsars like PSR J0740+6620 [55, 56] and PSR J0030+0451 [57, 58]. Additionally, the LIGO-Virgo gravitational wave observatory has contributed crucial data through neutron star merger events, such as GW170817 [59, 60], providing insights into neutron star masses, radii, and tidal deformabilities. These combined observations offer opportunities to constrain the presence of vacuum energy within compact objects, contributing to our understanding of fundamental physics in astrophysical contexts.

This paper is structured as follows. Sec. II details the derivation of the equations governing stellar structure, specifically incorporating the effects of vacuum energy into these equations. Our study examines two different models. In Sec. II A, we consider a simpler model where the ratio of the energy densities of ordinary matter to vacuum energy is treated as a constant, acting as a free parameter. Sec. II B extends this approach to include the derivation of structure equations that allow for densitydependent variations in this ratio. In Section III, we describe a set of state-of-the-art nuclear EOS used to represent normal matter in neutron stars. Sec. IV presents the outcomes of our calculations, with a particular focus on the mass-radius relationship and tidal deformabilities of neutron stars under these models. Finally, in Sec. V, we summarize our findings and discuss their implications for understanding the potential role of dark energy in neutron star physics.

## II. STELLAR STRUCTURE EQUATIONS WITH A VACUUM COMPONENT

In this section, we will incorporate the vacuum energy into the composition of neutron stars, assuming that it is admixed with ordinary matter. To this end, we will consider these compact stars to be spherically symmetric, static, and isotropic. Consequently, we can express the metric of these objects as:

$$\L a c { - \quad d s ^ { 2 } } = e ^ { 2 \nu ( r ) } d t ^ { 2 } - e ^ { 2 \mu ( r ) } d r ^ { 2 } - r ^ { 2 } ( d \theta ^ { 2 } + \sin ^ { 2 } \theta d \phi ^ { 2 } ). \ \ ( 1 )$$

It is common practice to introduce a new metric function m r ( ) defined by e 2 µ r ( ) = 1 -2 m r ( ) r . The quantity m r ( ) carries the interpretation of 'mass enclosed within a radius r ', and, from now on, units such that c = G = 1 are adopted.

It will be assumed here that the vacuum component ( v ) and ordinary matter ( m ) are characterized as separated perfect fluid components interacting only gravitationally. The total energy-momentum tensor, which encompasses both vacuum energy and ordinary matter, can be expressed as:

$$T _ { \mu \nu } ^ { ( t ) } = T _ { \mu \nu } ^ { ( m ) } + T _ { \mu \nu } ^ { ( v ) }, \quad \quad ( 2 ) \quad \text{wh}$$

where the individual components have the following form:

$$T _ { \mu \nu } = ( p + \epsilon ) u _ { \mu } u _ { \nu } - p g _ { \mu \nu }, \text{ \quad \ \ } ( 3 )$$

being p and ϵ the pressure and energy density, respectively.

As one may check, for a comoving observer ( u µ = ( √ g 00 ) -1 δ µ 0 , u µ = √ g 00 δ µ 0 ), the Tolman-OppenheimerVolkoff (TOV) equations describing the structure of stellar objects take the following form:

$$\frac { d m } { d r } \, = \, 4 \pi \epsilon r ^ { 2 }, \, \quad \ \ ( 4 ) \quad T \, \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\frac { \ddots } { d r } \ = \ - \frac { ( p + \epsilon ) \left [ m + 4 \pi r ^ { 3 } p \right ] } { r ^ { 2 } - 2 m r }, \quad \ \ ( 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\frac { d \nu } { d r } \ = \ \frac { m + 4 \pi r ^ { 3 } p } { r ^ { 2 } - 2 m r },$$

where m p , and ϵ are the sum of contributions from ordinary matter and vacuum:

$$m \, = \, m _ { m } + m _ { v }, \quad \quad \ ( 7 ) \quad ( 5 )$$

$$p \, = \, p _ { m } + p _ { v },$$

$$\epsilon \, = \, \epsilon _ { m } + \epsilon _ { v }. \quad \quad \quad ( 9 )$$

To solve the TOV equations, it is necessary to supplement them with equations of state (EOSs) that describe the relationship between pressure and energy density within the stellar interior. In the current scenario involving both ordinary matter and vacuum energy, we will specify separate EOSs for each component. For vacuum energy, we will adopt the standard relation p v = -ϵ v . Regarding ordinary matter, we will use a set of state-ofthe-art EOSs for hadronic matter. These EOSs will be thoroughly described in Section III. All of these hadronic EOSs are in agreement with current astrophysical constraints on neutron stars and meet experimental restrictions at the nuclear saturation density, as detailed in Sec. III.

Once the EOSs for ordinary matter and vacuum energy are specified, the three TOV equations will contain four unknown functions that need to be determined: m r ( ), p m ( r ), p v ( r ), and ν r ( ). To close the system and ensure a complete solution, an additional condition is required. In order to simplify the analysis, we adopt the following approach: we assume that ordinary matter constitutes a fraction y of the total energy density, while vacuum energy constitutes a fraction 1 -y . The fraction y will be specified in different ways to investigate the impact of varying vacuum energy contents on the stellar structure. In Sec. II A, we will assign y a set of constant values. Alternatively, in Sec. II B, we will allow y to vary across different layers of the star. Mathematically, y can be expressed as:

$$y = \frac { \epsilon _ { m } } { \epsilon _ { v } + \epsilon _ { m } }, \quad \quad ( 1 0 )$$

while 1 -y reads:

$$1 - y = \frac { \epsilon _ { v } } { \epsilon _ { v } + \epsilon _ { m } }. \quad \quad \quad ( 1 1 )$$

When y = 1, the situation corresponds to the absence of vacuum energy, while y = 0 represents a scenario in which the fluid is entirely composed of vacuum energy.

## A. Constant ratio of ordinary matter and vacuum

The TOV equations can be reformulated using only variables associated with ordinary matter and the y fraction. This approach allows us to rewrite Eq. (4) as follows:

$$\frac { d m } { d r } = 4 \pi r ^ { 2 } ( \epsilon _ { m } + \epsilon _ { v } ) = 4 \pi r ^ { 2 } \frac { \epsilon _ { m } } { y }. \quad \ \ ( 1 2 )$$

When y remains constant, the vacuum energy density is a fixed fraction of ϵ m , and, as such, ϵ v = ϵ v ( ϵ m ), and Eq. (5) can be rewritten using the following expressions:

$$\frac { d p } { d r } \, = \, \frac { d p _ { m } } { d r } - \frac { d \epsilon _ { v } } { d r } = \frac { d p _ { m } } { d r } - \frac { 1 - y } { y } \frac { d \epsilon _ { m } } { d r }, \quad ( 1 3 )$$

$$\frac { d \epsilon _ { m } } { d r } \, = \, \frac { d \epsilon _ { m } } { d p _ { m } } \frac { d p _ { m } } { d r }.$$

Substituting these formulas into Eq. (5), we obtain:

$$\ u \lim _ { \cdot \epsilon _ { v }. } \quad \frac { d p _ { m } } { d r } = & - \frac { \left ( p _ { m } + \epsilon _ { m } \right ) \left [ m + 4 \pi r ^ { 3 } \left ( p _ { m } - \frac { 1 - y } { y } \epsilon _ { m } \right ) \right ] } { r ^ { 2 } - 2 m r } \\ \intertext { n i c } \quad \times \frac { 1 } { \left ( 1 - \frac { 1 - y } { y } \frac { d \epsilon _ { m } } { d p _ { m } } \right ) }.$$

Note that the derivative dϵ m dp m is now required in order to solve the complete set of equations.

The denominator of Eq. (15) must always be positive to ensure a negative pressure gradient within the star. Consequently, we require y to fulfill the condition:

$$\left ( 1 - \frac { 1 - y } { y } \frac { d \epsilon _ { m } } { d p _ { m } } \right ) > 0. \text{ \quad \ \ } \text{(16)}$$

Recalling that the speed of sound in matter is given by c 2 s = dp m dϵ m (in units of the speed of light), we derive a lower limit for y :

$$y > \frac { 1 } { c _ { s } ^ { 2 } + 1 }. \quad \quad \quad ( 1 7 )$$

Due to causality, c s &lt; 1, setting an absolute lower limit of y &gt; 0 5. . However, in practical terms, many realistic EOSs for nuclear matter display a speed of sound below the conformal limit of 1 / √ 3, suggesting y &gt; 0 75. . Consequently, in the scenario of constant y , a maximum fraction of 25% of vacuum energy would be feasible at each layer of the star.

## B. Variable fractions of ordinary matter and vacuum

We will now shift our focus to the broader scenario where the fraction y of ordinary matter exhibits nontrivial variations throughout the star's interior. In this context, the derivative of the total pressure can be expressed as:

$$\frac { d p } { d r } = \frac { d p _ { m } } { d r } - \frac { d \epsilon _ { v } } { d r } = \frac { d p _ { m } } { d r } - \frac { 1 - y } { y } \frac { d \epsilon _ { m } } { d r } + \frac { \epsilon _ { m } } { y ^ { 2 } } \frac { d y } { d r }. \ \ ( 1 8 ) \quad \text{the a} _ { \text{gover} }.$$

Assuming, for the sake of simplicity, that the fraction y depends solely on the energy density of ordinary matter ( y = y ϵ ( m )), we obtain:

$$\frac { d p } { d r } = \frac { d p _ { m } } { d r } \left [ 1 + \frac { d \epsilon _ { m } } { d p _ { m } } \left ( - \frac { 1 - y } { y } + \frac { \epsilon _ { m } } { y ^ { 2 } } \frac { d y } { d \epsilon _ { m } } \right ) \right ]. \quad ( 1 9 ) \quad \text{ly} \ \varsigma$$

By replacing the previous expression in Eq. (5), one obtains:

$$\frac { d p _ { m } } { d r } \left [ 1 + \frac { d \epsilon _ { m } } { d p _ { m } } \left ( - \frac { 1 - y } { y } + \frac { \epsilon _ { m } } { y ^ { 2 } } \frac { d y } { d \epsilon _ { m } } \right ) \right ] = \\ = - \frac { ( p _ { m } + \epsilon _ { m } ) \left [ m + 4 \pi r ^ { 3 } \left ( p _ { m } - \frac { 1 - y } { y } \epsilon _ { m } \right ) \right ] } { r ^ { 2 } - 2 m r }. \\ \widehat { \mu } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

To describe the variation of y with the energy density of ordinary matter, we propose the following ansatz:

$$y ( \epsilon _ { m } ) = 1 + \frac { \beta - 1 } { 1 + \left ( \frac { \epsilon _ { * } } { \epsilon _ { m } } \right ) ^ { 4 } }, \quad \quad ( 2 1 ) \quad \text{pro$$

where β and ϵ ⋆ are free parameters. The behavior of the function y ϵ ( m ) is illustrated in Fig. 1. When β = 1, the function simplifies to a constant, y ϵ ( m ) = 1, representing the absence of dark energy throughout the stellar interior. For β values between 0 and 1, the function approaches a limit of 1 as ϵ m tends towards zero, indicating that vacuum energy vanishes near the stellar surface. As ϵ m approaches infinity, the function tends towards β , representing a scenario where the fraction of vacuum energy density attains the maximum possible value allowed by

FIG. 1. Fractions of ordinary matter, y , and vacuum energy, (1 -y ), as functions of ϵ m /ϵ ⋆ for β = 0 8. . Horizontal lines indicate the asymptotic limits, β for y and 1 -β for (1 -y ). At low densities, y tends toward 1 and (1 -y ) toward zero, evidencing the complete disappearance of vacuum energy. At high densities, the vacuum energy fraction asymptotically approaches 1 -β , representing the maximum fraction allowable under this ansatz.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-82431b6a-6ac4-4b39-ba50-ff4c05f0fa73.jpg)

the ansatz. The parameter ϵ ⋆ acts as a scaling factor that governs the range of ϵ m over which y differs significantly from the limiting values β and 1. Note that for ϵ ⋆ = 0, the value of y equals β at any energy density. In other words, the ansatz of Eq. (21) simplifies to the scenario of constant y as described in Sec. II A.

By using Eq. (21), we can easily determine the derivative dy/dϵ m , which is used in Eq. (20):

$$\frac { d y } { d \epsilon _ { m } } = \frac { 4 ( \beta - 1 ) ( \epsilon _ { * } / \epsilon _ { m } ) ^ { 4 } } { \epsilon _ { m } [ 1 + ( \epsilon _ { * }$$

## III. EQUATIONS OF STATE FOR ORDINARY MATTER

In this work, we employ a piecewise polytropic EOS to effectively represent complex, state-of-the-art nuclear EOSs. This approach simplifies intricate EOSs by approximating them with segmented polytropic relations, offering a versatile method for modeling dense matter within neutron stars. The number of segments and the polytropic indices in each segment can be adjusted to best fit different EOSs, retaining sufficient accuracy and enabling the representation of diverse microphysical approaches in a practical and flexible manner. Here, we adopt the Generalized Piecewise Polytropic (GPP) model as detailed in Ref. [61], which uses an ansatz that imposes continuity not only in pressure and energy density but also in the speed of sound.

The pressure and the energy density of each segment are written as functions of the rest mass density ρ as

follows:

$$p _ { m } ( \epsilon _ { m } ) = K \rho _ { m } ^ { \Gamma } + A,$$

$$\epsilon _ { m } ( \rho _ { m } ) = \frac { K } { \Gamma - 1 } \rho _ { m } ^ { \Gamma } + ( 1 + a ) \rho _ { m } - A, \quad ( 2 4 ) \quad \quad \quad$$

from which we derive:

$$\frac { d \epsilon _ { m } } { d p _ { m } } = \frac { 1 } { \Gamma - 1 } + \frac { ( 1 + a ) } { \Gamma K \rho _ { m } ^ { \Gamma - 1 } }. \quad \quad ( 2 5$$

The coefficients of adjacent segments are related by imposing continuity and differentiability of p and ϵ at the dividing density ρ i :

$$K _ { i + 1 } = K _ { i } \frac { \Gamma _ { i } } { \Gamma _ { i + 1 } } \rho _ { i } ^ { \Gamma _ { i } - \Gamma _ { i + i } }, \\.$$

$$A _ { i + 1 } = A _ { i } + \left ( 1 - \frac { \Gamma _ { i } } { \Gamma _ { i + 1 } } \right ) K _ { i } \rho _ { i } ^ { \Gamma _ { i } }, \quad \quad ( 2 7 )$$

$$a _ { i + 1 } = a _ { i } + \Gamma _ { i } \frac { \Gamma _ { i + 1 } - \Gamma _ { i } } { \left ( \Gamma _ { i + 1 } - 1 \right ) \left ( \Gamma _ { i } - 1 \right ) } K _ { i } \rho _ { i } ^ { \Gamma _ { i } - 1 }. \quad ( 2 8 ) \quad \vdots$$

From the extensive set of EOSs reproduced in Ref. [61], we have selected only those that match the predictions of chiral effective field theory (EFT) interactions. Chiral EFT provides a framework for a systematic expansion of nuclear forces at low momenta, explaining the hierarchy of two-, three-, and weaker higher-body forces [62-67]. Currently, microscopic calculations based on chiral EFT interactions allow for a reliable determination of the properties of neutron star matter up to densities somewhat above the nuclear saturation density n 0 . Specifically, the pressure is currently known to be roughly within ± 20% accuracy at n 0 [65, 67].

The selected models are presented in Fig. 2(a), along with the region delineating the constraints imposed by chiral EFT. Each model reflects different assumptions regarding the nature of the forces between nucleons and the potential emergence of exotic particles such as hyperons at extreme densities. The APR EOS is based on variational calculations using potentials resulting from fits to nucleon-nucleon scattering and properties of light nuclei [70]. The KDE0V, KDE0V1, SKI4, SKI6, SLY2, SLY230A, SLY4, and SLY9 models are part of the Skyrme interaction family [71]. In Skyrme-type models, the effective mass and nuclear potential are parametrized to reproduce properties of bulk nuclear matter and/or finite nuclei. Each variant of the Skyrme interaction represents a different parameter set, fitted to specific nuclear properties and experimental data. The MPA1 EOS is based on relativistic Brueckner-Hartree-Fock calculations, with a particular focus on contributions resulting from tensor interactions due to the exchange of π and ρ mesons and the dependence on neutron-proton asymmetry [72].

In Fig. 2(b), the mass-radius relationship for neutron stars is shown using the EOSs presented in Fig. 2(a), without considering any vacuum energy effects. The

FIG. 2. (a) Selected EOSs for ordinary matter used in this work. The gray shaded area in the lower-left corner represents the constraint imposed by chiral effective field theory up to 1 1 . n 0 [68, 69], a constraint satisfied by all of our hadronic EOSs. (b) Mass-radius relationship for the selected EOSs (without any vacuum contribution), illustrating the differences in predicted neutron star sizes and masses. Overlaying the theoretical curves are shaded regions, corresponding to observational constraints from ∼ 2 M ⊙ pulsars, NICER observations, the GW170817 event detected by LIGO/Virgo, and the very low mass object HESS J1731.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-598d17d2-47fd-4c94-a5dc-57eda6553ee2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-bb3aea6f-d6fe-48bc-bb33-0afda3010984.jpg)

KDE0V and KDE0V1 models fail to meet the observational constraints imposed by PSR J0740+6620 and 2 M ⊙ pulsar measurements. However, the remaining EOSs generally comply with all astrophysical constraints, despite some satisfying the HESS J1731 constraint only marginally. Subsequent figures will particularly emphasize the APR and MPA1 models. The APR model is somewhat less stiff than the MPA1 model, resulting in a lower predicted maximum mass. Furthermore, stellar radii within the APR model are considerably smaller than those predicted by the MPA1 model. This choice will provide deeper insight into the role of vacuum energy in EOSs with different degrees of stiffness.

FIG. 3. Effect of vacuum energy on APR and MPA1 EOSs. The dashed and dotted curves incorporate the effect of vacuum energy for different β (or y in the constant case). Note that there are significant changes in the stiffness of the EOS, especially at low-density regimes for ϵ ⋆ = 0 ϵ 0 (constant y ). This effect does not contradict nuclear physics, as the modified EOS only becomes relevant in strong gravitational fields. In general, the EOS exhibits only slight softening at high densities.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-d3f8752a-b7bf-4897-80f6-bee3b79a31b8.jpg)

## IV. RESULTS

contribution at low, intermediate, or high densities, respectively.

In this section, we will solve the stellar structure equations described in the previous sections, using the set of EOSs outlined in Sec. III. This will involve adopting different proportions and distributions of vacuum energy within the interior of the star, as explored in Sec. IVA. In Sec. IVB, our focus will be on the mass-radius relationship, while in Sec. IVC, we will analyze tidal deformabilities.

## A. Vacuum energy effects in the equation of state

Let us analyze in detail how two of the previously selected nuclear matter EOSs, namely APR and MPA1, are altered when vacuum energy is incorporated according to the prescriptions discussed earlier. As pointed out at the end of Sec. II B, the scenario with a constant y corresponds to the case of variable y with ϵ ⋆ = 0, and this notation will henceforth denote the constant y model. This scenario is depicted in Fig. 3(a) for the APR EOS and in Fig. 3(e) for the MPA1 EOS. The remaining panels of Fig. 3 show the impact of vacuum energy within the variable y scenario, employing different values for ϵ ⋆ , specifically 2 ϵ 0 , 4 ϵ 0 , and 8 ϵ 0 , with ϵ 0 = 2 7 . × 10 14 g / cm 3 representing the nuclear saturation density. These values were chosen to set the onset of the vacuum energy

As expected, in all the cases presented in Fig. 3, the EOS softens upon the addition of vacuum energy due to its negative pressure contribution. It is also observed that the smaller the parameter β , the more pronounced the softening of the EOS. This trend suggests that increasing the contribution of vacuum energy will progressively decrease the maximum mass and shift the mass-radius curves toward smaller radii, as will indeed be demonstrated in Sec. IVB. An interesting aspect of the models with constant y is that, as shown in Figs. 3(a) and 3(e), the total pressure vanishes at quite high densities, and is negative below 10 14 -10 15 g / cm . 3 This occurs because, in the constant y scenario, the vacuum contribution to the pressure, represented by ϵ v = (1 -y y ) -1 ϵ m , remains substantial even at low pressures and eventually surpasses the matter pressure at a relatively high density. Conversely, as shown in Fig. 3, in the variable y scenario, the effect of vacuum energy is mitigated at densities lower than approximately ϵ ⋆ (also illustrated in Fig. 1), preventing the total pressure from becoming negative in most cases.

FIG. 4. Mass-Radius relation for APR and MPA1 with the inclusion of vacuum energy related with the ordinary matter in the star by a function y , and observational constraints (GW170817 from [73], GW190425 from [74], PSR J0030+0451 from [75], PSR J0348+0432 from [76], PSR J1614+2230 from [77] and PSR J0740+6620 from [78]).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-5196790f-c4ed-418a-a01d-f5a269032461.jpg)

## B. The influence of vacuum energy on stellar structure

We have calculated the mass-radius relationship for the set of nuclear matter EOSs selected previously. For brevity and focus, we will only show plots for two of these EOSs, specifically APR (Figs. 4(a)-4(d)) and MPA1 (Figs. 4(e)-4(h)). Each panel corresponds to a different scenario regarding the proportion of vacuum energy. Observational constraints are overlaid on the curves, shown by the shaded regions and labeled with the names of the observed objects. Figs. 4(a) and 4(e) illustrate the massradius relationship for neutron stars under the assumption of a constant y ( ϵ ⋆ = 0). Figs. 4(b)-4(d) and 4(f)4(h) display the results with density-dependent values of y under varying thresholds of ϵ ⋆ , specifically at 2 ϵ 0 , 4 ϵ 0 , and 8 ϵ 0 , respectively. Within each panel various selections for the parameter β are explored.

larger, simply because the presence of vacuum energy is restricted to higher density regimes, thus exerting a lesser impact on the star's structure (cf. Fig. 1). For instance, with the APR EOS, when ϵ ⋆ = 4 ϵ 0 , the results for β ≲ 0 96 become incompatible with the 2 . M ⊙ constraint. However, for ϵ ⋆ = 8 ϵ 0 , a significantly broader range of β values (1 &lt; β ≲ 0 87) yield . M -R relationships that are consistent with the 2 M ⊙ constraint.

Across all examined scenarios, a consistent trend emerges: both radii and masses tend to be smaller as the amount of vacuum energy is increased. This effect is particularly pronounced in the scenario where the parameter y is held constant (Figs. 4(a) and 4(e)) because vacuum energy exerts a substantial influence on the EOS throughout the entire density spectrum. On the other hand, in the scenario with density dependent y the effect is substantially reduced as the parameter ϵ ⋆ becomes

Another significant effect of vacuum energy is related to the qualitative shape of the mass-radius curves. In the case of constant y , even a small fraction of vacuum energy leads to a qualitative change in the curves, making them resemble those typical of self-bound stars, such as strange quark stars [79-82]. In fact, unlike standard neutron stars, which exhibit larger radii at lower masses, strange quark stars show a decrease in radius with decreasing mass. This effect is observed in Figs. 4(a) and 4(e) with vacuum energy. For density-dependent y , the sensitivity to vacuum energy strongly depends on the value of ϵ ⋆ . When vacuum appears at high densities, its impact is similar to what is observed in hybrid stars, which feature a quark matter core surrounded by nuclear matter [83, 84]: There is a decrease in M max and a shift of the M -R curves towards smaller radii, particularly noticeable in Figs. 4(c), 4(d), 4(g), and 4(h). This shift affects objects with higher masses and leaves low mass neutron stars unaffected, ensuring the qualitative shape of the curve remains unchanged. On the other

hand, when sufficiently large values of β along with sufficiently small values of ϵ ⋆ are adopted, it is still possible to achieve curves with the same shape as self-bound stars (see Figs. 4(b), 4(f) with β = 0 76 0 87). . , . However, in all the cases we analyzed, these curves are incompatible with the 2 M ⊙ constraint.

The behavior described above can be understood as follows: The vacuum energy contributes additively to density but negatively to pressure: ϵ = ϵ m + ϵ v and p = p m -ϵ v . Therefore, vacuum energy invariably softens the total EOS. When vacuum energy is restricted to the stellar core, this softening leads to the well-known effects of reducing M max and shifting the mass-radius curves towards smaller radii [83, 85]. However, as the influence of vacuum energy extends to lower densities, it begins to significantly bind the normal matter at the star's surface. Specifically, the condition p R ( ) = 0 (which implies p m ( R ) = ϵ v ( R )) allows for configurations where the matter pressure at the surface is non-vanishing, while the total pressure remains zero. This effect 'truncates' the star compared to models without vacuum energy, resulting in a much smaller radius. This mechanism is similar to that which shapes strange quark stars, the only difference being that the vacuum energy in strange quark stars originates from strong interactions [80, 81, 86].

Observed masses and radii of neutron stars can provide valuable constraints on the amount of vacuum energy they contain. These constraints, however, depend on the specific EOS employed. EOSs that are close to exceeding the limits set by astrophysical observations will become incompatible with these observations, even with minor additions of vacuum energy to the matter EOS. In the case of constant y , radii are strongly reduced even for tiny additions of vacuum energy. As a consequence, the radius limits set by NICER are the first observational restrictions to be violated by models with vacuum energy. For example, with y = 0 995, the APR curves are . pushed beyond the accepted radius limits of PSR J0740 and PSR J0030. For the MPA1 EOS, the impact is similarly substantial, but because the curve without vacuum energy is positioned much further to the right, reflecting its stiffer nature, the β value at which it fails to meet NICER's constraints is lower: y ≈ 0 92. . In the case of density-dependent y , the decrease in radius tends to be less pronounced, particularly in models where vacuum energy is confined to the stellar core. Consequently, the first astrophysical constraint to be infringed by models with vacuum energy can vary; sometimes it is the NICER constraint and other times it is the two solar mass limit.

All the features previously described for APR and MPA1 are also corroborated by the results we obtained for the other EOSs included in our analysis. To summarize the findings across these EOSs, we will focus primarily on the 2 M ⊙ constraint, as it is the most robust among the astrophysical restrictions. In Table I, we present the minimum β required for compliance with the 2 M ⊙ constraint for all EOSs and for all selected vacuum energy density thresholds ϵ ⋆ . As anticipated, EOSs that are al-

TABLE I. For each EOS, this table presents the minimum values of the parameter β required for agreement with the 2 M ⊙ constraint across different vacuum energy density thresholds: ϵ ⋆ = 0 2 , ϵ 0 , 4 ϵ 0 , 8 ϵ 0 . The column M 0 max indicates the maximum mass predicted by each EOS in the absence of vacuum energy.

| EOS     | M 0 max   | β min   | β min       | β min       | β min       |
|---------|-----------|---------|-------------|-------------|-------------|
|         |           | ϵ ⋆ = 0 | ϵ ⋆ = 2 ϵ 0 | ϵ ⋆ = 4 ϵ 0 | ϵ ⋆ = 8 ϵ 0 |
| APR     | 2.170     | 0.970   | 0.970       | 0.965       | 0.870       |
| MPA1    | 2.465     | 0.870   | 0.910       | 0.870       | 0.760 ∗     |
| SKI4    | 2.170     | 0.970   | 0.970       | 0.960       | 0.870       |
| SKI6    | 2.191     | 0.970   | 0.970       | 0.950       | 0.860       |
| SLY2    | 2.055     | 0.990   | 0.990       | 0.985       | 0.980       |
| SLY230A | 2.100     | 0.980   | 0.980       | 0.975       | 0.930       |
| SLY4    | 2.053     | 0.990   | 0.990       | 0.985       | 0.980       |
| SLY9    | 2.157     | 0.980   | 0.980       | 0.970       | 0.870       |

ready near the 2 M ⊙ maximum mass threshold in the absence of vacuum energy present a higher β min .

## C. Tidal deformability

Neutron stars experience significant tidal deformation during the inspiral phase of a neutron star-neutron star merger leaving a detectable imprint on the observed gravitational waveform of the merger. The extent of this effect is quantifiable in terms of the star's dimensionless tidal deformability, which characterizes the induced mass quadrupole moment in response to an external tidal field. Specifically, the induced mass-quadrupole moment Q ij and the applied tidal field ϵ ij are related to this parameter by Q ij = -Λ M ϵ 5 ij .

The dimensionless tidal deformability parameter can be expressed as:

$$\Lambda = \frac { 2 } { 3 } k _ { 2 } C ^ { - 5 },$$

where C ≡ M/R is the compactness of the star and k 2 is the tidal Love number, which, in turn, is given by

$$\i u n \\ \text{the} \\ \text{dels} \\ \mathbb { E } \, & \, = \frac { 8 C ^ { 5 } } { 5 } \left ( 1 - 2 C \right ) ^ { 2 } \left [ 2 + 2 C \left ( y _ { R } - 1 \right ) - y _ { R } \right ] \\ \text{mit} \, & \, \times \left [ 4 C ^ { 3 } \left [ 1 3 - 1 1 y _ { R } + C \left ( 3 y _ { R } - 2 \right ) + 2 C ^ { 2 } \left ( 1 + y _ { R } \right ) \right ] \\ \text{and} \\ \text{in} \\ \text{mal} \, & \, + 2 C \left [ 6 - 3 y _ { R } + 3 C \left ( 5 y _ { R } - 8 \right ) \right ] \right ] ^ { - 1 }. \\ \text{mar} \, & \, + 2 C \left [ 6 - 3 y _ { R } + 3 C \left ( 5 y _ { R } - 8 \right ) \right ] \right ] ^ { - 1 }. \\ \text{long} \, & \, \text{Here, $y_{R}\equiv y_{R}(\text{), being $y_{r}(\text{)}$ the solution of the following}$$

Here, y R ≡ y R ( ), being y r ( ) the solution of the following first-order differential equation [87]:

$$r \frac { d y ( r ) } { d r } + y ( r ) ^ { 2 } + y ( r ) F ( r ) + r ^ { 2 } Q ( r ) = 0. \quad ( 3 1 )$$

FIG. 5. Tidal deformability for the APR and MPA1 EOSs with vacuum energy. The color coding and organization are the same as in Figs. 3 and 4. The constraint on Λ 1 4 . from the merger event GW170817 [60] is shown as vertical error bars in each panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01006v1-figure-2412f7f3-44b3-43c2-a9cf-ce0a86478f91.jpg)

In this equation, the coefficients are:

$$F ( r ) \ = \ & \left [ 1 - 4 \pi r ^ { 2 } ( \epsilon - p ) \right ] \left [ 1 - \frac { 2 m } { r } \right ] ^ { - 1 }, \quad & ( 3 2 ) \quad & \text{tidal} \quad & \text{to th} \\ Q ( r ) \ = \ & 4 \pi \left [ 5 \epsilon + 9 p + \frac { \epsilon + p } { c _ { s } ^ { 2 } } - \frac { 6 } { 4 \pi r ^ { 2 } } \right ] \left [ 1 - \frac { 2 m } { r } \right ] ^ { - 1 } \quad & \text{as th} \\ \quad & - \frac { 4 m ^ { 2 } } { r ^ { 4 } } \left [ 1 + \frac { 4 \pi r ^ { 3 } p } { m } \right ] ^ { 2 } \left [ 1 - \frac { 2 m } { r } \right ] ^ { - 2 }, \quad & ( 3 3 ) \quad & \text{energy} \\ \quad \text{--} \cdots \ 2 \colon + \colon + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \delta + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dts + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dists + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dclots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \ dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \ Dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \Dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \sigma + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dstots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \frown + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dsts + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \d.: + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \d: + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \d:[ + 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 + 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 2 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 -1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1$$

where c 2 s is the squared speed of sound, which within the present model can be expressed as:

$$c _ { s } ^ { 2 } \equiv \frac { d p } { d \epsilon } = \frac { \left ( \frac { \epsilon _ { m } } { y ^ { 2 } } \frac { d y } { d \epsilon _ { m } } - \frac { 1 } { y } + 1 + \frac { d p _ { m } } { d \epsilon _ { m } } \right ) } { \left ( - \frac { \epsilon _ { m } } { y ^ { 2 } } \frac { d y } { d \epsilon _ { m } } + \frac { 1 } { y } \right ) }. \quad ( 3 4 ) \quad \text{MP} ^ { \prime } \text{inv} \\ \text{per} }$$

As in previous sections, m = m m + m v , p = p m + p v , and ϵ = ϵ m + ϵ v . The boundary condition for Eq. (31) at r = 0 is given by y (0) = 2. In summary, the tidal Love number can be obtained once an EOS is supplied and the TOV equations together with Eq. (31) are integrated.

In Fig. 5 we show Λ as a function of M for the EOSs of Fig. 3. The results indicate that the highest values of tidal deformability for a given mass always correspond to the curves with β = 1 (no vacuum energy effect). As the amount of vacuum energy increases (indicated by lower β values), tidal deformability decreases for a given mass. Also, as the threshold for the onset of vacuum energy effects increases ( ϵ ∗ goes from 0 to 8 ϵ 0 ), the impact of vacuum energy on tidal deformability becomes less pronounced for both EOS, with the curves beginning to converge, especially for low M . For the APR EOS, the effect of vacuum energy is more noticeable at lower ϵ ∗ thresholds, as evidenced by the significant spread of the Λ curves seen in Figs. 5(a) through 5(c). For the MPA1 EOS, a similar pattern is observed; however, the influence of vacuum energy appears to be less severe compared to the APR EOS. At higher vacuum energy density thresholds, as seen in Figs. 5(d) and 5(h), the curves for both EOSs behave more similarly. Overall, our results show a clear trend of decreasing tidal deformability with increasing contributions of vacuum energy, particularly at lower energy density thresholds.

The primary constraint on the value of Λ currently stems from the detection of gravitational waves from the binary neutron star merger event GW170817 by LIGO/Virgo [59, 60, 88]. Assuming that both objects have spins within the range observed in Galactic binary neutron stars and that they share the same EOS, the tidal deformability of a 1 4 . M ⊙ neutron star can be estimated to be Λ 1 4 . = 190 +390 -120 at the 90% level [60].

We now examine how these calculated values of Λ compare with the restriction imposed by the GW170817 merger event. For both EOS, models without vacuum energy ( β = 1) are already consistent with the error bars of GW170817. As the amount of vacuum energy increases, especially in cases with lower ϵ ∗ thresholds, the curves diverge substantially from the error bars. For the APR EOS, the broad spread of theoretical curves at

lower ϵ ∗ thresholds means that only curves with higher β values, closer to 1, are consistent with the observational constraints from GW170817. Stiffer models such as MPA1 might accommodate a broader range of vacuum energy contributions while still being consistent with GW170817. At higher vacuum energy density thresholds, with ϵ ∗ increasing to 8 ϵ 0 , the theoretical curves for different β tend to coincide (see Figs. 5(d) and 5(h)). This convergence results in all models becoming compatible with the Λ 1 4 . constraints from GW170817, since the variations in the curves caused by the addition of vacuum energy occur at masses significantly larger than 1 4 . M ⊙ .

## V. SUMMARY AND CONCLUSIONS

In this work, we explored the role of vacuum energy within neutron stars, assuming it is admixed with ordinary matter. To describe ordinary matter we employed a polytropic representation of a set of state-of-the-art nuclear EOSs, while for vacuum energy, we adopted the standard relation p v = -ϵ v . To complete the system of equations, we assumed that, at each point of the star, ordinary matter constitutes a fraction y of the total energy density, and vacuum energy constitutes a fraction 1 -y . The fraction y was specified according to different models. Initially, we assumed that y remains constant throughout the stellar interior. We then adopted a density-dependent approach as defined in Eq. (21), where at low densities y → 1, indicating the complete disappearance of vacuum energy, and at high densities, y → β , where β represents the maximum possible fraction of vacuum energy present in the system.

Our selection of ordinary matter EOSs was specifically aimed to align with the predictions of chiral EFT interactions near the nuclear saturation density. Detailed analysis in Sec. IVA shows that incorporating vacuum energy modifies these EOSs, typically resulting in their softening due to the negative pressure contribution of vacuum energy. The degree of softening intensifies with higher values of the parameter β . In models maintaining a constant fraction y , total pressure vanishes at relatively high densities, as the substantial contribution from vacuum energy surpasses the matter pressure. In contrast, for models with a density-dependent y , where vacuum energy's influence disappears at densities below approximately ϵ ⋆ , the total pressure remains positive in all relevant cases.

We subsequently examined the mass-radius relationship for a set of nuclear matter EOSs, with a particular focus on the APR EOS and the significantly stiffer MPA1 EOS. Across all scenarios, we observed reductions in both stellar radii and masses as the contributions of vacuum energy increased. In models with a constant y , even minimal increases in vacuum energy profoundly affect the qualitative shape of the mass-radius curves, rendering them similar to those observed in strange quark stars, where R → 0 as M → 0. In contrast, in models with a density-dependent y , the changes are more quantitative in nature. Indeed, the impact of vacuum energy mirrors that seen in quark-hadron hybrid stars, i.e. it reduces the maximum masses and shifts the mass-radius curves towards smaller radii. These changes predominantly affect more massive neutron stars while leaving those of lower mass unaffected, thereby preserving the overall shape of the curves.

We evaluated our theoretical curves against observed masses and radii of neutron stars. For models with constant y , even minimal increases in vacuum energy significantly reduce radii, leading these models to first infringe the radius limits set by NICER. In models with densitydependent y , the radius reduction is generally milder, especially when vacuum energy is concentrated in the stellar core. Therefore, the first astrophysical constraints to be challenged by these models vary, sometimes being the NICER constraints and at other times the two solar mass limit. Given the robustness of the 2 M ⊙ constraint, we detailed in Table I the minimum β values needed to satisfy this limit for all selected EOSs across various vacuum energy density thresholds ϵ ⋆ . Predictably, EOSs close to the 2 M ⊙ maximum mass threshold without vacuum energy require a higher β min .

Finally, we showed that, at fixed stellar mass, the tidal deformability Λ diminishes as the amount of vacuum energy increases. Moreover, as the threshold for vacuum energy effects rises (from 0 to 8 ϵ 0 ), the influence of vacuum energy on tidal deformability becomes less significant, leading to a convergence of curves, particularly at lower masses. We then juxtaposed the theoretical values of Λ against the constraint from the GW170817 event. Our EOSs without vacuum energy match the observational error bars from GW170817. However, as the vacuum energy fraction is increased, particularly at lower ϵ ∗ thresholds, the curves significantly deviate from these error bars. At higher ϵ ∗ , the theoretical curves for different β values tend to coincide, especially at low masses, and all of them align with the Λ 1 4 . constraint from GW170817.

To conclude, we mention some challenges of our approach and the broader implications of our findings. Notably, the levels of vacuum energy we investigated exceed, by more than 40 orders of magnitude, those derived from cosmological estimates, where vacuum energy density approximates ϵ ∼ 10 -29 g / cm . A significant the3 oretical challenge lies in elucidating the mechanisms by which dark energy could be significantly concentrated within neutron stars via gravitational interactions. Furthermore, even assuming such a substantial accumulation of dark energy is possible, determining its presence within neutron stars observationally remains problematic. The influence of vacuum energy could be indistinguishable from effects attributed to other complex stellar phenomena, such as those associated with strange quark stars or hybrid stars, both of which are well-documented in existing astrophysical literature. This overlap complicates the interpretation of observational data and highlights the need for innovative approaches to distinguish between

these various influences on neutron star behavior.

## ACKNOWLEDGMENTS

LFA was supported by Coordena¸ c˜o a de Aperfei¸oamento c de Pessoal de N´ ıvel Superior (CAPES),

- [1] P. Beltracchi and P. Gondolo, Formation of dark energy stars, Phys. Rev. D 99 , 044037 (2019), arXiv:1810.12400 [gr-qc].
- [2] E. Gliner, Algebraic properties of the energy-momentum tensor and vacuum-like states of matter, Soviet Journal of Experimental and Theoretical Physics 22 , 378 (1966).
- [3] J. M. Bardeen, Non-singular general-relativistic gravitational collapse, in Proceedings of the International Conference GR5 (Tbilisi, USSR, 1968) pp. 174-175.
- [4] A. Borde, Open and closed universes, initial singularities and inflation, Phys. Rev. D 50 , 3692 (1994), arXiv:grqc/9403049.
- [5] A. Borde, Regular black holes and topology change, Phys. Rev. D 55 , 7615 (1997), arXiv:gr-qc/9612057.
- [6] S. Zhou, J. Chen, and Y. Wang, Geodesic Structure of Test Particle in Bardeen Spacetime, Int. J. Mod. Phys. D 21 , 1250077 (2012), arXiv:1112.5909 [gr-qc].
- [7] S. R. Coleman and F. De Luccia, Gravitational Effects on and of Vacuum Decay, Phys. Rev. D 21 , 3305 (1980).
- [8] K. Sato, M. Sasaki, H. Kodama, and K.-i. Maeda, Creation of Wormholes by First Order Phase Transition of a Vacuum in the Early Universe, Prog. Theor. Phys. 65 , 1443 (1981).
- [9] S. K. Blau, E. I. Guendelman, and A. H. Guth, The Dynamics of False Vacuum Bubbles, Phys. Rev. D 35 , 1747 (1987).
- [10] E. Farhi and A. H. Guth, An Obstacle to Creating a Universe in the Laboratory, Phys. Lett. B 183 , 149 (1987).
- [11] G. Chapline, Dark energy stars, eConf C041213 , 0205 (2004), arXiv:astro-ph/0503200.
- [12] P. O. Mazur and E. Mottola, Gravitational Condensate Stars: An Alternative to Black Holes, Universe 9 , 88 (2023), arXiv:gr-qc/0109035.
- [13] P. O. Mazur and E. Mottola, Gravitational vacuum condensate stars, Proc. Nat. Acad. Sci. 101 , 9545 (2004), arXiv:gr-qc/0407075.
- [14] M. Visser and D. L. Wiltshire, Stable gravastars: An Alternative to black holes?, Class. Quant. Grav. 21 , 1135 (2004), arXiv:gr-qc/0310107.
- [15] R. Chan, M. F. A. da Silva, and P. Rocha, Gravastars and Black Holes of Anisotropic Dark Energy, Gen. Rel. Grav. 43 , 2223 (2011), arXiv:1009.4403 [gr-qc].
- [16] A. DeBenedictis, D. Horvat, S. Ilijic, S. Kloster, and K. S. Viswanathan, Gravastar solutions with continuous pressures and equation of state, Class. Quant. Grav. 23 , 2303 (2006), arXiv:gr-qc/0511097.
- [17] B. M. N. Carter, Stable gravastars with generalised exteriors, Class. Quant. Grav. 22 , 4551 (2005), arXiv:grqc/0509087.
- [18] I. Dymnikova, Vacuum nonsingular black hole, Gen. Rel. Grav. 24 , 235 (1992).

Brazil. JASL is partially supported by the Brazilian agencies, CNPq under grant 310038/2019-7, CAPES (88881.068485/2014) and FAPESP (LLAMA Project no. 11/51676-9). GL acknowledges the partial financial support from CNPq (grant 316844/2021-7) and FAPESP (grant 2022/02341-9).

- [19] M. Campos and J. A. S. Lima, Black hole formation with an interacting vacuum energy density, Physical Review D 86 , 10.1103/physrevd.86.043012 (2012).
- [20] S. Smerechynskyi, M. Tsizh, and B. Novosyadlyj, Impact of dynamical dark energy on the neutron star equilibrium, Journal of Cosmology and Astroparticle Physics 2021 (02), 045.
- [21] K. G. Sagar, B. Pandey, and N. Pant, Hybrid star model in tolman-buchdahl metric potentials with coupled dark energy and baryonic matter, Astrophysics and Space Science 367 , 72 (2022).
- [22] P. Rej and A. Karmakar, Charged strange star coupled to anisotropic dark energy in tolman-kuchowicz spacetime, The European Physical Journal C 83 , 699 (2023).
- [23] K. P. Das, U. Debnath, and S. Ray, Dark energy star: Physical constraints on the bounds, Fortschritte der Physik 71 , 2200148 (2023).
- [24] X. Li, F. Wang, and K. S. Cheng, Gravitational effects of condensate dark matter on compact stellar objects, JCAP 10 , 031, arXiv:1210.1748 [astro-ph.CO].
- [25] P. Mukhopadhyay and J. Schaffner-Bielich, Quark stars admixed with dark matter, Phys. Rev. D 93 , 083009 (2016), arXiv:1511.00238 [astro-ph.HE].
- [26] G. Baym, D. H. Beck, P. Geltenbort, and J. Shelton, Testing dark decays of baryons in neutron stars, Phys. Rev. Lett. 121 , 061801 (2018), arXiv:1802.08282 [hepph].
- [27] J. Ellis, G. H¨tsi, K. Kannike, L. Marzola, M. Raidal, and u V. Vaskonen, Dark Matter Effects On Neutron Star Properties, Phys. Rev. D 97 , 123007 (2018), arXiv:1804.01418 [astro-ph.CO].
- [28] Y. Dengler, J. Schaffner-Bielich, and L. Tolos, Second Love number of dark compact planets and neutron stars with dark matter, Phys. Rev. D 105 , 043013 (2022), arXiv:2111.06197 [astro-ph.HE].
- [29] B. Kain, Dark matter admixed neutron stars, Phys. Rev. D 103 , 043009 (2021), arXiv:2102.08257 [gr-qc].
- [30] D. Sen and A. Guha, Implications of feebly interacting dark sector on neutron star properties and constraints from GW170817, Mon. Not. Roy. Astron. Soc. 504 , 3354 (2021), arXiv:2104.06141 [hep-ph].
- [31] J. C. Jim´ enez and E. S. Fraga, Radial Oscillations of Quark Stars Admixed with Dark Matter, Universe 8 , 34 (2022), arXiv:2111.00091 [hep-ph].
- [32] Z. Miao, Y. Zhu, A. Li, and F. Huang, Dark Matter Admixed Neutron Star Properties in the Light of X-Ray Pulse Profile Observations, Astrophys. J. 936 , 69 (2022), arXiv:2204.05560 [astro-ph.HE].
- [33] C. J. Horowitz and S. Reddy, Gravitational Waves from Compact Dark Objects in Neutron Stars, Phys. Rev. Lett. 122 , 071102 (2019), arXiv:1902.04597 [astroph.HE].

- [34] S. Maedan, Influence of the equation of state on a compact star made of hidden-sector nucleons, PTEP 2020 , 033B07 (2020), arXiv:1908.00711 [hep-ph].
- [35] G. Narain, J. Schaffner-Bielich, and I. N. Mishustin, Compact stars made of fermionic dark matter, Phys. Rev. D 74 , 063003 (2006), arXiv:astro-ph/0605724.
- [36] C. Kouvaris and N. G. Nielsen, Asymmetric Dark Matter Stars, Phys. Rev. D 92 , 063526 (2015), arXiv:1507.00959 [hep-ph].
- [37] M. I. Gresham and K. M. Zurek, Asymmetric Dark Stars and Neutron Star Stability, Phys. Rev. D 99 , 083008 (2019), arXiv:1809.08254 [astro-ph.CO].
- [38] K.-L. Leung, M.-c. Chu, and L.-M. Lin, Tidal deformability of dark matter admixed neutron stars, Phys. Rev. D 105 , 123010 (2022), arXiv:2207.02433 [astro-ph.HE].
- [39] N. Rutherford, G. Raaijmakers, C. Prescod-Weinstein, and A. Watts, Constraining bosonic asymmetric dark matter with neutron star mass-radius measurements, Phys. Rev. D 107 , 103051 (2023), arXiv:2208.03282 [astro-ph.HE].
- [40] G. Panotopoulos and I. Lopes, Dark matter effect on realistic equation of state in neutron stars, Phys. Rev. D 96 , 083004 (2017), arXiv:1709.06312 [hep-ph].
- [41] A. Nelson, S. Reddy, and D. Zhou, Dark halos around neutron stars and gravitational waves, JCAP 07 , 012, arXiv:1803.03266 [hep-ph].
- [42] J. Ellis, A. Hektor, G. H¨tsi, K. Kannike, L. Marzola, u M. Raidal, and V. Vaskonen, Search for Dark Matter Effects on Gravitational Signals from Neutron Star Mergers, Phys. Lett. B 781 , 607 (2018), arXiv:1710.05540 [astro-ph.CO].
- [43] A. G. Riess, S. Casertano, W. Yuan, L. M. Macri, and D. Scolnic, Large magellanic cloud cepheid standards provide a 1% foundation for the determination of the hubble constant and stronger evidence for physics beyond λ cdm, The astrophysical Journal 876 , 85 (2019).
- [44] J. V. Cunha, L. Marassi, and J. A. S. Lima, Constraining H from the Sunyaev-Zel'dovich effect, galaxy cluster X0 ray data and baryon oscillations, Mon. Not. R. Astron. Soc. 379 , L1 (2007).
- [45] E. Abdalla, Abell´n, G. Franco, A. Aboubrahim, A. Aga nello, ¨ . O Akarsu, Y. Akrami, G. Alestas, D. Aloni, L. Amendola, L. A. Anchordoqui, et al. , Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies, Journal of High Energy Astrophysics 34 , 49 (2022).
- [46] J. C. Carvalho, J. A. S. Lima, and I. Waga, Cosmological consequences of a time-dependent λ -term, Phys. Rev. D 46 , 7 (1992).
- [47] J. A. S. Lima, Thermodynamics of decaying vacuum cosmologies, Physical Review D 54 , 54 2571 (1996).
- [48] J. A. S. Lima, Alternative dark energy models: an overview, Brazilian J. Physics 34 , 194 (2004).
- [49] J. A. S. Lima, S. Basilakos, and J. Sol`, Thermodynama ical aspects of running vacuum models, The European Physical Journal C 76 , 10.1140/epjc/s10052-016-4060-6 (2016).
- [50] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, A two-solar-mass neutron star measured using Shapiro delay, Nature 467 , 1081 (2010), arXiv:1010.5788 [astro-ph.HE].
- [51] Z. Arzoumanian et al. , The NANOGrav 11-year data set: High-precision timing of 45 millisecond pulsars, The Astrophysical Journal Supplement Series 235 , 37 (2018).
- [52] J. Antoniadis et al. , A Massive Pulsar in a Compact Relativistic Binary, Science 340 , 448 (2013), arXiv:1304.6875 [astro-ph.HE].
- [53] H. T. Cromartie et al. , Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar, Nature Astronomy 4 , 72 (2020), arXiv:1904.06759 [astroph.HE].
- [54] E. Fonseca et al. , Refined Mass and Geometric Measurements of the High-mass PSR J0740+6620, ApJL 915 , L12 (2021), arXiv:2104.00880 [astro-ph.HE].
- [55] M. C. Miller et al. , The Radius of PSR J0740+6620 from NICER and XMM-Newton Data, ApJ 918 , L28 (2021), arXiv:2105.06979 [astro-ph.HE].
- [56] T. E. Riley et al. , A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy, ApJ 918 , L27 (2021), arXiv:2105.06980 [astro-ph.HE].
- [57] M. C. Miller and et al., PSR j0030+0451 mass and radius from NICER data and implications for the properties of neutron star matter, Astrophys. J. Lett. 887 , L24 (2019).
- [58] T. E. Riley and et al., A NICER view of PSR j0030+0451: Millisecond pulsar parameter estimation, Astrophys. J. Lett. 887 , L21 (2019).
- [59] B. P. Abbott et al. (LIGO Scientific, Virgo), GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral, Phys. Rev. Lett. 119 , 161101 (2017), arXiv:1710.05832 [gr-qc].
- [60] B. Abbott et al. (LIGO Scientific, Virgo), GW170817: Measurements of neutron star radii and equation of state, Phys. Rev. Lett. 121 , 161101 (2018), arXiv:1805.11581 [gr-qc].
- [61] M. F. O'Boyle, C. Markakis, N. Stergioulas, and J. S. Read, Parametrized equation of state for neutron star matter with continuous sound speed, Phys. Rev. D 102 , 083027 (2020).
- [62] L. Coraggio, J. W. Holt, N. Itaco, R. Machleidt, and F. Sammarruca, Reduced regulator dependence of neutron-matter predictions with perturbative chiral interactions, Phys. Rev. C 87 , 014322 (2013), arXiv:1209.5537 [nucl-th].
- [63] S. Gandolfi, J. Carlson, and S. Reddy, The maximum mass and radius of neutron stars and the nuclear symmetry energy, Phys. Rev. C 85 , 032801 (2012), arXiv:1101.1921 [nucl-th].
- [64] J. W. Holt, N. Kaiser, and W. Weise, Chiral Fermi liquid approach to neutron matter, Phys. Rev. C 87 , 014338 (2013), arXiv:1209.5296 [nucl-th].
- [65] K. Hebeler and A. Schwenk, Chiral three-nucleon forces and neutron matter, Phys. Rev. C 82 , 014314 (2010), arXiv:0911.0483 [nucl-th].
- [66] F. Sammarruca, B. Chen, L. Coraggio, N. Itaco, and R. Machleidt, Dirac-Brueckner-Hartree-Fock versus chiral effective field theory, Phys. Rev. C 86 , 054317 (2012), arXiv:1209.5001 [nucl-th].
- [67] I. Tews, T. Kr¨ger, K. Hebeler, and A. Schwenk, Neuu tron matter at next-to-next-to-next-to-leading order in chiral effective field theory, Phys. Rev. Lett. 110 , 032504 (2013), arXiv:1206.0025 [nucl-th].
- [68] K. Hebeler, J. M. Lattimer, C. J. Pethick, and A. Schwenk, Equation of state and neutron star properties constrained by nuclear physics and observation,

- Astrophys. J. 773 , 11 (2013), arXiv:1303.4662 [astroph.SR].
- [69] E. Annala, T. Gorda, A. Kurkela, J. N¨ attil¨, a and A. Vuorinen, Evidence for quark-matter cores in massive neutron stars, Nature Physics 16 , 907 (2020), arXiv:1903.09121 [astro-ph.HE].
- [70] A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, The Equation of state of nucleon matter and neutron star structure, Phys. Rev. C 58 , 1804 (1998), arXiv:nuclth/9804027.
- [71] M. Dutra, O. Lourenco, J. S. Sa Martins, A. Delfino, J. R. Stone, and P. D. Stevenson, Skyrme Interaction and Nuclear Matter Constraints, Phys. Rev. C 85 , 035201 (2012), arXiv:1202.3902 [nucl-th].
- [72] H. M¨ uther, M. Prakash, and T. L. Ainsworth, The nuclear symmetry energy in relativistic Brueckner-HartreeFock calculations, Phys. Lett. B 199 , 469 (1987).
- [73] C. A. Raithel, F. Ozel, and D. Psaltis, Tidal deforma-¨ bility from GW170817 as a direct probe of the neutron star radius, The Astrophysical Journal Letters 857 , L23 (2018).
- [74] B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, et al. , GW190425: Observation of a compact binary coalescence with total mass ∼ 3 4 . M ⊙ , The Astrophysical Journal Letters 892 , L3 (2020).
- [75] M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, A. K. Harding, W. C. G. Ho, et al. , PSR J0030+0451 Mass and Radius from NICER data and implications for the properties of neutron star matter, The Astrophysical Journal Letters 887 , L24 (2019).
- [76] J. Antoniadis, P. C. C. Freire, N. Wex, T. M. Tauris, R. S. Lynch, M. H. van Kerkwijk, et al. , A Massive Pulsar in a Compact Relativistic Binary, Science 340 , 448 (2013), arXiv:1304.6875 [astro-ph.HE].
- [77] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, A two-solar-mass neutron star measured using Shapiro delay, Nature 467 , 1081

(2010), arXiv:1010.5788 [astro-ph.HE].

- [78] H. T. Cromartie, E. Fonseca, S. M. Ransom, P. B. Demorest, Z. Arzoumanian, H. Blumer, P. R. Brook, M. E. DeCesar, T. Dolch, et al. , Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar, Nature Astronomy 4 , 72 (2020), arXiv:1904.06759 [astroph.HE].
- [79] C. Alcock, E. Farhi, and A. Olinto, Strange stars, Astrophys. J. 310 , 261 (1986).
- [80] G. Lugones and J. E. Horvath, High-density QCD pairing in compact star structure, Astron. Astrophys. 403 , 173 (2003), arXiv:astro-ph/0211638.
- [81] C. V´ asquez Flores and G. Lugones, Constraining color flavor locked strange stars in the gravitational wave era, Phys. Rev. C 95 , 025808 (2017), arXiv:1702.02081 [astroph.HE].
- [82] G. Lugones and A. G. Grunfeld, Excluded volume effects in the quark-mass density-dependent model: implications for the equation of state and compact star structure, Phys. Rev. D 109 , 063025 (2023), arXiv:2312.16095 [nucl-th].
- [83] C. H. Lenzi and G. Lugones, Hybrid stars in the light of the massive pulsar PSR J1614-2230, Astrophys. J. 759 , 57 (2012), arXiv:1206.4108 [astro-ph.SR].
- [84] M. Mariani and G. Lugones, Quark-hadron pasta phase in neutron stars: The role of medium-dependent surface and curvature tensions, Phys. Rev. D 109 , 063022 (2024), arXiv:2308.13973 [nucl-th].
- [85] G. Lugones, M. Mariani, and I. F. Ranea-Sandoval, A model-agnostic analysis of hybrid stars with reactive interfaces, JCAP 03 , 028, arXiv:2106.10380 [nucl-th].
- [86] G. Lugones and J. E. Horvath, Quark - diquark equation of state and compact star structure, Int. J. Mod. Phys. D 12 , 495 (2003), arXiv:astro-ph/0203069.
- [87] S. Postnikov, M. Prakash, and J. M. Lattimer, Tidal Love Numbers of Neutron and Self-Bound Quark Stars, Phys. Rev. D 82 , 024016 (2010), arXiv:1004.5098 [astro-ph.SR].
- [88] B. P. Abbott et al. (LIGO Scientific, Virgo), Properties of the binary neutron star merger GW170817, Phys. Rev. X 9 , 011001 (2019), arXiv:1805.11579 [gr-qc].
L. F. Araujo J. A. S. Lima G. Lugones
Impact of dark energy on the structure of neutron stars: The vacuum case
The potential role of a cosmic vacuum dark component in the properties of
neutron stars is investigated. It is assumed that the static, spherically
symmetric distribution of matter within neutron stars is supported by two
distinct components: ordinary matter and a vacuum fluid. For normal matter we
use a set of state-of-the-art nuclear matter equations of state, each grounded
in nuclear physics experiments. The vacuum energy component is inhomogeneously
distributed within the star and obeys the standard equation of state ($p_v =
-\epsilon_v$). This is characterized by an energy density fraction $y =
\epsilon_m/ (\epsilon_m + \epsilon_v)$, which we model as either a constant or
radius-dependent. Our findings reveal that the inclusion of vacuum energy
significantly affects the mass-radius relationships in neutron stars,
influencing both the maximum achievable masses and the qualitative form of
these relationships. Some constraints from current multimessenger observational
data limiting the amount of vacuum energy within neutron stars are also
discussed.