Abstract
Introduction:
This paper investigates the impact of differential rotation on the bulk properties and onset of rotational instabilities in neutron stars at finite temperatures up to 50 MeV.
Methods:
Utilizing the relativistic Brueckner-Hartree-Fock (RBHF) formalism in full Dirac space, the study constructs equation of state (EOS) models for hot neutron star matter, including conditions relevant for high temperatures. These finite-temperature EOS models are applied to compute the bulk properties of differentially rotating neutron stars with varying structural deformations.
Results:
The findings demonstrate that the stability of these stars against bar-mode deformation, a key rotational instability, is only weakly dependent on temperature. Differential rotation significantly affects the maximum mass and radius of neutron stars, and the threshold for the onset of bar-mode instability shows minimal sensitivity to temperature changes within the examined range.
Discussion:
These findings are crucial for interpreting observational data from neutron star mergers and other high-energy astrophysical events. The research underscores the necessity of incorporating differential rotation and finite temperature effects in neutron star models to predict their properties and stability accurately.
1 Introduction
Neutron stars provide a unique, naturally occurring laboratory for studying matter at extreme pressures and densities not reproducible by experiments in terrestrial laboratories (see, for instance, [1–3]). The cold, highly isospin asymmetric matter within the core of a massive neutron star can reach densities up to an order of magnitude higher than nuclear saturation density. During a binary neutron star merger event, the resulting matter may promptly collapse into a black hole or form a remnant neutron star. If formed, the remnant star is characterized by a high mass, extreme temperatures on the order of 50–100 MeV, and rapid differential rotation [4–7]. These massive, differentially rotating remnant stars may also deviate from spherical or axial symmetry by exhibiting extreme triaxial deformations. The structural deformation, thermal pressure, and differential rotation allow the remnant to remain stable on short, dynamical timescales, for masses that would be otherwise unstable in the static and uniform-rotation cases.
Differential rotation in neutron stars has been explored in the literature through numerical simulations [8–11], with more recent studies incorporating finite temperature equation of state (EOS) models [12, 13]. The inclusion of temperature when modeling the EOS of neutron star matter, however, is a formidable task. Theoretical modeling of neutron star matter as a dense, many-body system can be done in a phenomenological or ab initio framework [3, 14–16, 16–18, 20]. Phenomenological models employ density functional theories with effective nucleon-nucleon (NN) interactions to reproduce the empirical saturation properties of symmetric nuclear matter while adhering to constraints extracted from nuclear physics and astrophysics [19, 21–23]. In contrast, ab initio methods use realistic NN interactions determined by nucleon-nucleon scattering data and the properties of the deuteron. Relativistic ab initio methods, such as the relativistic Brueckner-Hartree-Fock (RBHF) approximation, closely reproduce the saturation properties of empirical data [24–29]. The RBHF approximation couples the propagation of baryons to the many-body background and encapsulates dynamical correlations between baryons, computed using a relativistic scattering matrix. Importantly, the RBHF method does not involve adjustable parameters, distinguishing it from many phenomenological models.
While most studies using the RBHF approximation are conducted at zero temperature, our previous work [30] extended the approximation to model nuclear matter at finite temperatures, where EOS models were derived in different temperature regimes. The nuclear EOS models incorporate finite temperatures in a self-consistent manner, unlike previous studies that added thermal effects to models of cold nuclear matter [31–33]. This framework ensures a more comprehensive and accurate representation of the thermodynamic properties of dense nuclear matter. In previous work, these EOS models were used to determine bulk properties of non-rotating, uniformly rotating, and differentially rotating neutron stars.
In this paper, we expand on previous work by exploring more deformed stars and dynamical rotation instabilities. This paper is organized as follows: Section 2 describes the theoretical framework for deriving EOS models at finite temperatures using the RBHF approximation and for constructing equilibrium models of differentially rotating stars. Section 3 presents the calculated results, where Section 3.1 shows stellar sequences over ranges of structural deformation, Section 3.2 discusses the stability of calculated models to dynamical bar mode excitation, Section 3.3 presents density and frequency profiles for stars with high degrees of differential rotation and structural deformation, and Section 3.4 discusses how various approximations introduced in the numerical calculations may influence the presented results. Section 4 gives a summary of the work presented.
2 Theoretical framework
This section discusses the theoretical framework for constructing equation of state (EOS) models for neutron star matter at finite temperatures using the relativistic Brueckner-Hartree-Fock (RBHF) theory. The EOS models are used as input to construct equilibrium models of differentially rotating objects, for which the theory is described below.
2.1 Relativistic Brueckner-Hartree-Fock theory at finite temperatures
The essential structure of modeling nuclear matter using RBHF theory is outlined in this section. A more detailed explanation of the approach can be found in Poschenrieder and Weigel [24]; Weber [14], with finite temperature extensions given in our previous work [30, 34]. Nuclear matter at supranuclear densities can be described as a complex, many-body system whose dynamics are governed by the Lagrangian density:where denotes the Lagrangian of non-interacting nucleons, is the Lagrangian density of different free meson fields, and describes the interaction between nucleons and mesons. In the relativistic framework, nucleons are treated as effective Dirac particles, which are described by the relativistic Dirac equation. The equations of motion for the various particle fields within the many-body system are derived from the Euler-Lagrange equation and solved using the Martin-Schwinger hierarchy of coupled Green’s functions [35].
The formal structure of the RBHF approach is to solve a system of highly nonlinear, coupled equations, which include the Dyson equation for the two-body Green’s function , a Bethe-Salpeter type integral equation for the scattering matrix , and the equation for the self-energy . In this formalism, the Dyson equation is used to determine the two-body Green’s function for all nucleons as:where and and are the Dirac matrices. Following the determination of the , the in-medium scattering matrix is determined with Bethe-Salpeter type integral equation:where is the intermediate nucleon-nucleon propagator and represents the repeated sums of two-particle interactions given by a one-boson-exchange (OBE) potential, which describes the interaction among two nucleons in terms of the exchange of scalar, pseudo-scalar, and vector mesons. In this work, we use the so-called Bonn B potential [36] to describe the OBE interaction, which employs a pseudoscalar type of pion-nucleon coupling. This potential is widely used in high-density nuclear matter studies due to its reliability and numerical stability, particularly in RBHF theory. While newer relativistic OBE potentials have been proposed (e.g., in [37]), the Bonn-B potential has been extensively tested in dense nuclear environments, making it a robust and reliable choice for modeling the equation of state under extreme conditions.
The final coupled equation in the formal scheme is for the self-energy , written in momentum space as:The self-consistent calculations are carried out using a complete basis of particles and antiparticles to decouple the integral equations and make the two-body propagator diagonal, where are the helicity eigenvalues.
The particle propagator in Equation 1 takes the form of the Brueckner propagator, which is defined as the following at finite temperatures:where signifies the Fermi-Dirac distribution functions, given by:where “1” indicates the positive energy states and “2” indicates the thermally-excited negative energy states. We recall here that at finite temperatures, the behavior of nuclear matter undergoes important modification, attributed to thermal baryonic excitations surpassing the Fermi surface. As , the Fermi-Dirac distribution for positive energy states becomes:and for negative energy states, .
An elegant technique used to make the many-body equations numerically tractable and to calculate the key quantities of many-body systems is to utilize the spectral representation of the function [24]. can then be defined in Fourier space at finite temperatures as [14, 38]:where represents the spectral function which is dependent on the single-particle energy , is the nucleon chemical potential, and is used to circumvent a singularity occurring as integrals are carried out in the complex plane. The temperature inclusion arises through , where is the Stefan-Boltzmann constant and is the temperature.
Once a self-consistent solution to the coupled system of equations is found, the self-energy and spectral function are used to determine the EOS. The number density of the system follows from:where is the time-like component of the spectral function and denotes the Fermi-Dirac distributions. The pressure of the system at finite temperatures is determined from the free energy per nucleon, denoted as , which is defined as:where is the internal energy, is the temperature, and is the entropy ter Haar and Malfliet [39]. In this approach, is held constant and is determined by:Both particles and antiparticles contribute to , but the antiparticle contribution is very small . Therefore, Equation 2 can be approximated with only the particle contribution, . Once the entropy , and subsequently the free energy , are calculated, the pressure is derived as:
Using the outlined theory, two models for the EOS of neutron star matter are constructed at temperatures T = 10 and 50 MeV, shown visually in Figure 1. As shown in previous work Farrell and Weber [30], the maximum mass of each EOS for non-rotating and uniformly rotating stellar sequences at their mass-shedding limit is over 2 , as required by observational constraints. These models will be used as input to determine equilibrium models for differentially rotating objects, of which the formalism is discussed below.
FIGURE 1
2.2 Differential rotation
The theoretical framework for modeling differential rotation in neutron stars described in this work follows from the framework laid out by Komatsu et al. [40], which was then modified in Cook et al. [8] (referred to as CST throughout the text). The equations shown in this section directly follow the modifications introduced in CST.
To model differentially rotating neutron stars, we begin with the definition of the line element [8]:where the metric potentials , , , and are dependent on both the radial and polar coordinates. For both uniform and differential rotation, Equation 3 models neutron stars as stationary, axisymmetric configurations of a (self-gravitating) perfect fluid [8]. Under the assumption of neutron star matter as a perfect fluid, sources of non-isotropic stresses such as magnetic fields or heat transport are ignored [41]. This assumption also allows neutron star matter to be described by the energy-momentum (or stress-energy) tensor given by:where is the fluid’s 4-velocity, and are indices ranging from 0 to 3, and and are given by the underlying EOS. Equilibrium models for neutron stars must obey the equation of hydrostatic equilibrium as Einstein’s field equation, given as:where is the Ricci tensor, is the curvature scalar, and is the metric tensor.
Equilibrium models for neutron stars must obey the equation of hydrostatic equilibrium, which has the form:where is enthalpy as a function of the pressure, is the time-like component of the 4-velocity , and the subscripts and denote the variable’s value at the pole or center, respectively. The integrand of the final integral term, , is the function that defines the rotation law of the matter in the case of differential rotation. Following CST, we define as a linear rotation law:where is a parameter that dictates the degree of differential rotation within the star. In the case of uniform rotation, Equation 4 disappears as the value for the frequency at the center of the star is constant throughout (i.e., ). It is important to note that the choice of rotation law directly impacts the maximum mass for a given EOS. While many implementations of CST use the linear rotation law in Equation 4 (see, for example, [42, 43]), other studies of differential rotation have explored the impact of using either modified versions of the linear law or non-linear rotation laws, as shown in Galeazzi et al. [10]; Hanauske et al. [11]; Zhou et al. [44].
Using the linear rotation law in Equation 4, the equation of hydrostatic equilibrium can be integrated to give:where is the rotation parameter scaled as . The matrix of angular frequency can be derived using the following equation:As the rotation parameter appears throughout the numerical scheme scaled and inverted, we follow the lead of previous work which parameterized calculated sequences of differentially rotating stars by values of = 0.3, 0.5, 0.7, and 1.0 [8, 10]. Uniform rotation is obtained in the limit , and an upper bound of the scaled rotation parameter in this study is .
In the numerical scheme, the metric potentials (, , , and ) are used to determine a value for the radius at the equator, of the star:which is equivalent to Equation 5 evaluated at the location of the maximum (denoted by subscript ) density of the star. In CST, this location is assumed to be at the star’s center. The equation for updating changes when the maximum density within the star is not in the center, which is the case in very deformed configurations. We modify the original CST algorithm by instead fixing the maximum interior density, thus redefining Equation 6 as:where the subscript denotes the quantity’s value at the coordinates of the maximum density.
3 Results
Equilibrium models of differentially rotating stars at finite temperatures are computed using the theoretical formalism described in the sections above. For the two EOS models at and 50 MeV, both stellar sequences and individual stellar models are constructed to explore the impact varying degrees of differential rotation has on bulk properties and structural deformation. Computed models are tested for stability against the dynamical bar-mode instability. It is important to note that while the core is treated at finite temperatures, the crust of the neutron star is nevertheless treated at zero temperature, which may impact the properties of neutron stars at the high-temperature end of this study.Additionally, approximations are introduced both in the determination of the finite temperature EOS models using RBHF theory and when numerically modeling differential rotation in neutron stars. These approximations, which may introduce some level of uncertainty in the results presented in this section, are discussed in detail below in Section 3.4.
3.1 Stellar sequences: Varying the ratio of polar to equatorial radius
In this section, stellar sequences are constructed over a range of constant central densities for the two EOS models at temperatures of 10 and 50 MeV. These sequences are calculated with a fixed value for the rotation parameter set the degree of differential rotation; as mentioned in Section 2.2, sequences are parameterized by fixing to be 0.3, 0.5, 0.7, or 1.0. When assuming the maximum density of the star is no longer in the center, more extreme configurations with lower values of can be calculated. Sequences presented in this section are also calculated over a range of values to demonstrate the parameter’s impact on the maximum mass for a given EOS.
For each EOS, sequences over a range of constant central densities are computed for varying values of . These sequences are shown for the EOS model with temperature MeV in Figure 2 and for MeV in Figure 3. In both figures, a sequence of stars rotating uniformly at their mass-shedding limit is plotted for comparison. As echoed in previous work, more extreme differential rotation paired with more extreme structural deformation results in higher masses at lower central densities, which then taper off to the Kepler limit at higher central densities. This trend is also encountered in Morrison et al. [43], who make similar modifications to the CST algorithm.
FIGURE 2
FIGURE 3
3.2 Dynamical bar-mode instability
Rotating neutron stars formed from a core-collapse suprenova (CCSN) or binary stellar mergers may experience nonaxisymmetric instabilities that directly impact their rotation rates and overall stability. Previous studies (see [45, 46]) in Newtonian gravity have shown rotational instabilities arise from non-radial toroidal modes, i.e., , which result in the stability parameter, or the ratio of rotational to gravitational energy , exceeding some critical value . These rotational instabilities are likely to impact a star’s gravitational radiation signal, making the study of such instabilities an important topic in the wake of new gravitational wave detectors. In this section, we will focus on determining if calculated stellar models are subject to the so-called bar-mode instability, where , which is expected to be the fastest-growing mode and the subject of many instability studies for both uniformly and differentially rotating neutron stars.
Two mechanisms cause rotating stars to be unstable to bar-mode deformation: secular and dynamical instabilities. In Newtonian theory, uniformly rotating incompressible neutron stars become secularly unstable to bar-mode deformation at a critical value of with similar findings in studies of post-Newtonian theories [47]. In general, the secular instability grows only in the presence of a small dissipative mechanism like viscosity or gravitational radiation at lower rotation rates [48]. The secular instability usually has a longer growth time when compared to the dynamic timescale of the system. A similar critical value of has been observed in numerical studies of relativistic stars but is also dependent on the compaction (M/R) of the star and the dissipative mechanism. For example, viscosity-driven secular instability has been shown to occur at in more compact configurations with higher rotation rates, but gravitational radiation-driven instabilities occur at at lower rotation rates. A more in-depth review of this topic is given by Paschalidis and Stergioulas [49].
The dynamical bar-mode instability occurs independent of any dissipative mechanism and with a growth rate determined by the dynamical timescale of the system, which is generally faster than the timescale of growth for secular instabilities. Therefore, numerical simulations of hydrodynamical equations are necessary to determine the onset threshold of the dynamical bar-mode instability. Many simulations have been carried out in Newtonian theory, the consensus of which gives the critical value [4]. Simulations of the dynamical bar-mode instability in general relativity are less common, as solving the nonlinear hydrodynamical equations in full relativity is more complex. However, there have been reliable studies carried out by Shibata and Uryu (2000); Saijo et al. [50] for uniformly rotating stars and by Bodenheimer and Ostriker [51]; Shibata et al. [52, 53]; Camarda et al. [54]; Di Giovanni et al. [55] for differentially rotating stars.
The relativistic simulation of differentially rotating stars carried out by Shibata et al. [52] finds that the critical value of the stability parameter for the dynamic bar-mode instability is , slightly lower than the Newtonian limit. These simulations were carried out using the same linear rotation law as in this chapter (and the proceeding one) for similar values of . Therefore, we adhere to their specified threshold of to determine whether differentially rotating stars are stable or unstable to the dynamical bar-mode deformation. Specifically, we use the upper limit for the critical threshold, so stars deemed unstable in this section will have a stability parameter .
A visual representation of stable and unstable models using this criterion at MeV for two degrees of differential rotation is shown in Figure 4. The two panels in each row show the same information, mass vs. central density for a range of values, where the figures in the right column categorize each stellar model as stable with a blue square or unstable with a red dot based on the calculated stability parameter for each stellar model. The top row presents a lesser degree of differential rotation , where 81.6% of calculated stellar models are considered stable against dynamical bar-mode deformation. In contrast, the bottom row has the highest degree of differential rotation , and the increase in differential rotation decreases the percentage of stable models to 67.2%. For both cases of differential rotation, the majority of the unstable configurations have small values for , a trend that will be echoed in the next section.
FIGURE 4
Table 1 gives the percentage of unstable models and the average of those unstable models for four degrees of differential rotation for both temperatures of 10 and 50 MeV. The general trend shows the percentage of unstable models decreases as temperature increases. For both temperatures, generally tends to fall between 0.4 and 0.5. Stars with are extremely deformed and thus likely vulnerable to dynamical instabilities. These findings imply temperature plays less of a role when compared to the star’s deformation on its stability against dynamical bar-mode excitation.
TABLE 1
| EOS | % unstable | ||
|---|---|---|---|
| T = 10 MeV | 0.3 | 23.2 | 0.47 |
| 0.5 | 27.6 | 0.46 | |
| 0.7 | 25.5 | 0.51 | |
| 1.0 | 27.0 | 0.52 | |
| T = 50 MeV | 0.3 | 18.4 | 0.43 |
| 0.5 | 21.4 | 0.43 | |
| 0.7 | 21.7 | 0.49 | |
| 1.0 | 32.8 | 0.47 |
The percentage of unstable models and the average , , of unstable models for four degrees of differential rotation EOS models at 10 and 50 MeV.
3.3 Structural deformation
In this section, we examine density and frequency maps of individual stellar models. For each temperature, two models are computed: the first with a lesser degree of differential rotation ( = 0.3) and the second with a larger degree of differential rotation ( = 1.0). All four models are computed with an = 0.4, a significant degree of structural deformation. Information on the stars’ masses, radii, and stability parameters are given in Table 2.
TABLE 2
| EOS | () | () | (km) | (km) | ||||
|---|---|---|---|---|---|---|---|---|
| MeV | 0.3 | 0.4 | 1.69 | 1.84 | 20.41 | 8.16 | 1.13 | 0.01 |
| 1.0 | 0.4 | 1.91 | 2.61 | 12.94 | 5.18 | 5.88 | 0.63 | |
| MeV | 0.3 | 0.4 | 1.78 | 1.93 | 22.25 | 8.90 | 1.13 | 0.01 |
| 1.0 | 0.4 | 2.08 | 2.84 | 13.65 | 5.19 | 5.95 | 0.75 |
Bulk properties of highly deformed neutron stars, all with = 0.4, at two temperatures: 10 MeV and 50 MeV. For both temperatures, results are presented for a small degree of differential rotation and a large degree of differential rotation . Calculated properties include the gravitational mass , baryonic mass , equatorial radius , polar radius , the ratio of central to equatorial frequency , and stability parameter . The MeV entries correspond to Figure 5 and the MeV entries correspond to Figure 6; see text for more details.
The density and frequency maps for MeV are shown in Figure 5, where the top row shows results for = 0.3 and the bottom for . The lesser degree of differential rotation is shown very clearly in the frequency map, where the difference of frequency values in the star only spans 320 Hz. In contrast, the star with a higher degree of differential rotation sees a span of 5,600 Hz from the highest to lowest frequency values. The density maps show an interesting depiction of the overall structural deformation. For , the star takes an ellipsoid-like shape, but for , the star instead takes a quasi-toroidal shape. While each star was initialized with the same central density of 400 MeV/fm3, the star with experiences much higher densities outside of the center of the star. A similar story is seen for MeV in Figure 6 but with even higher densities observed at the higher temperature.
FIGURE 5
FIGURE 6
As shown in Section 3.2, the stability of the star depends not only on the deformation characterized by but also on the degree of differential rotation. In Table 2, all four models have the same value of . For both temperatures, the stars computed with have stability parameters well below the critical limit of 0.25. However, the stars with are well above the critical limit for , where the stability parameter is 0.64 for MeV and 0.75 for MeV. These stars are well beyond the threshold for the dynamical bar-mode instability and thus likely unphysical.
The timescale over which the dynamical bar-mode instability develops, also known as the dynamical timescale, is proportional to . While this is not a definite indication of the growth time of the instability, which would instead require a full simulation in both time and space (see Ref. [52] for a good example), for extremely unstable configurations the dynamical timescale can provide some idea of how long these stars may exist with the bar-mode instability excited before collapse. For the unstable configuration at MeV, the dynamical timescale is 35.8 s, and for the unstable configuration at MeV, the dynamical timescale is 34.9 s.
3.4 Key approximations
The preceding sections present results dependent on the the underlying theoretical frameworks discussed in Section 2, which employ important approximations, discussed below.
The finite temperature EOS models at and 50 MeV are derived from calculations utilizing the RBHF approximation, described in Section 2.1. While this method offers significant improvements over other approaches to modeling relativistic numerical matter, it still involves key approximations [36]. For instance, nucleon-nucleon interactions are modeled using OBE potentials, with the Bonn-B potential selected in this paper for stability at high densities. Although RBHF theory does not require adjustable parameters, which is an advantage, it omits higher-order quantum corrections [56, 57]. One of the most important higher-order effects missing in RBHF is the inclusion of three-nucleon forces, which are known to significantly influence high-density matter [58]. Other potentially influential effects excluded in the RBHF approximation include higher-order relativistic many-body corrections, ring diagrams, vertex corrections, self-energy insertions, and medium polarization effects; a more complete discussion of these effects can be found in [14, 59, 60].
Variations in the choice of OBE potential and omission of quantum corrections, such as three-body forces, can introduce errors in the RBHF EOS of up to 10%, as demonstrated in Brockmann and Machleidt [36]. Using chiral effective field theory, Hebeler and Schwenk [61] and Tews et al. [62] concluded that including three-body forces affects neutron star radius and mass predictions by 5%–10%.
In addition, the EOS models include finite temperatures (10 MeV and 50 MeV), which are essential for modeling neutron stars formed in extreme events like supernovae or binary neutron star mergers Steiner et al. [63]. However, as mentioned before, the neutron star crust is assumed to be at zero temperature, while the core is modeled at finite temperatures. This approximation may introduce inconsistencies in the determination of properties like neutron star mass and radius, especially in higher temperature regimes. In particular, the radius increases due to thermal pressure in the crust, potentially making the star slightly larger.
When modeling differential rotation in neutron stars, the rotation profile, given in Equation 4, is parameterized by , representing the degree of differential rotation. This profile, while consistent with the original CST algorithm discussed in Section 2.2, is chosen for mathematical convenience rather than based on astrophysical observations, potentially limiting its accuracy for real neutron stars. Uncertainties in the choice of the rotation profile may cause errors of 5%–15% in star mass and radius predictions; for more on uncertainties related to differential rotation profiles and their impact on neutron star stability and mass limits, see Baumgarte and Shapiro [64]. Additionally, the assumption of axial symmetry is a simplification. In reality, neutron stars may exhibit more complex geometries, especially under differential rotation and dynamical instabilities.
For a comprehensive review discussing EOS uncertainties, neutron star mass-radius relationships, and the challenges of matching theoretical models with astrophysical data, see Lattimer and Prakash [65].
4 Discussion and conclusion
In this paper, we present a comprehensive investigation into the properties of differentially rotating neutron stars at finite temperatures up to 50 MeV. In Section 2.1, we detailed the process of constructing models for the equation of state of neutron star matter at two temperatures, 10 and 50 MeV, using the relativistic Brueckner-Hartree-Fock (RBHF) formalism modified to include thermal effects, utilizing the Bonn-B potential for the one-boson exchange (OBE) interaction. The inclusion of temperature is essential for the realistic modeling of extreme astrophysical events, such as binary neutron star (BNS) mergers or core-collapse supernovae, where differential rotation is prevalent.
The RBHF approach goes well beyond standard relativistic mean-field (RMF) calculations and relativistic Hartree-Fock (RHF) methods. The RBHF formalism includes the relativistic scattering matrix, which accounts for dynamical correlations among nucleons that are ignored in RMF and RHF. Furthermore, in sharp contrast to RMF and RHF, there are no adjustable parameters in the RBHF approach, making it a more fundamental and predictive theory compared to RMF and RHF models.
The two EOS models were used as input to the numerical scheme to determine the bulk properties of differentially rotating compact objects, focusing on heavily deformed objects characterized by the ratio of their polar to equatorial radii, . In Section 3.1, we constructed stellar sequences with varying values (0.8–0.3) over a range of central densities for all four finite temperature EOS models. Our results are consistent with previous literature at zero temperature, demonstrating that a higher degree of differential rotation results in higher mass stars at lower central densities compared to lower values of , especially as decreases.
In Section 3.2, we explored the stability of the calculated stellar models against rotational instabilities, specifically the dynamical bar-mode instability. At a temperature of 10 MeV, the average value for stars unstable against dynamical bar-mode excitement is 0.49. For the higher temperature of 50 MeV, the average for unstable stars is 0.455. These results highlight that deformation, characterized by , has a more significant impact on the star’s stability against dynamical instabilities compared to temperature.
In Section 3.3, we presented individual stellar maps of density and frequency distributions for two degrees of differential rotation for each EOS model. As expected, the frequency range for higher degrees of differential rotation (and higher values of was much wider. Higher values of also resulted in more pronounced structural deformations, indicating a strong dependence of the rotational profile on the degree of differential rotation.Finally, in Section 3.4, key approximations and their influence on calculated results are discussed. Variations in the choice of OBE potential and omission of quantum corrections like three-body forces can introduce errors in the RBHF EOS of up to 10%, which can, in turn, affect neutron star radius and mass predictions by 5%–10%. Additional uncertainty may be introduced by the choice of a linear rotational profile described in Section 2.2, on the order of 5%–15% in star mass and radius predictions.In conclusion, our study underscores the necessity of incorporating differential rotation and finite temperature effects in neutron star models to predict their properties and stability accurately. The stability of differentially rotating neutron stars against bar-mode deformation, a key rotational instability, is only weakly dependent on temperature. Differential rotation significantly affects the maximum mass and radius of neutron stars, and the threshold for the onset of bar-mode instability shows minimal sensitivity to temperature changes within the examined range. These findings are crucial for interpreting observational data from neutron star mergers and other high-energy astrophysical events (see, for instance, [66–74]).
The insights gained from this research are particularly relevant in the context of observations from current and future gravitational wave detectors such as LIGO [75], Virgo [76], KAGRA [77], and the upcoming Einstein Telescope [78]. Moreover, X-ray and radio telescopes, including the Chandra X-ray Observatory [79], XMM-Newton [80], the Very Large Array (VLA) [81], the Square Kilometre Array (SKA) [82], and the Five-hundred-meter Aperture Spherical Radio Telescope (FAST) [83] provide crucial observational data that can further constrain the models presented here.
Future work should extend these models to include additional physical effects, such as magnetic fields and more sophisticated treatments of thermal transport processes, to provide an even more comprehensive understanding of neutron star dynamics and stability under extreme conditions.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
DF: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. FW: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Supervision, Validation, Writing–original draft, Writing–review and editing. RN: Conceptualization, Funding acquisition, Methodology, Software, Validation, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. DF and FW are supported by the National Science Foundation (USA) under Grant No. PHY-2012152.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
BeckerW. Neutron stars and pulsars. Astrophysics and space science library. Berlin Heidelberg: Springer (2009).
2.
ReddyS. Neutron stars as astrophysical laboratories. Rev Mod Phys (2012) 84:85–120. 10.1103/RevModPhys.84.85
3.
BurgioGSchulzeH-JVidañaIWeiJ-B. Neutron stars and the nuclear equation of state. Prog Part Nucl Phys (2021) 120:103879. 10.1016/j.ppnp.2021.103879
4.
ShibataMUryuK. Simulation of merging binary neutron stars in full general relativity: γ= 2 case. Phys Rev D (2000) 61:064001. 10.1103/physrevd.61.064001
5.
ShibataMKiuchiK. Gravitational waves from remnant massive neutron stars of binary neutron star merger: viscous hydrodynamics effects. Phys Rev D (2017) 95:123003. 10.1103/physrevd.95.123003
6.
FujibayashiSKiuchiKNishimuraNSekiguchiYShibataM. Mass ejection from the remnant of a binary neutron star merger: viscous-radiation hydrodynamics study. Astrophysical J (2018) 860:64. 10.3847/1538-4357/aabafd
7.
PeregoABernuzziSRadiceD. Thermodynamics conditions of matter in neutron star mergers. The Eur Phys J A (2019) 55:124. 10.1140/epja/i2019-12810-7
8.
CookGBShapiroSLTeukolskySA. Rapidly rotating neutron stars in general relativity: realistic equations of state. Astrophysical J (1994) 424:823–45. 10.1086/173934
9.
SchaabCWeigelMK. Differential rotation in neutron stars. Phys Lett B (1999) 450:187–94. 10.1016/S0370-2693(99)00195-3
10.
GaleazziFYoshidaSEriguchiY. Differentially-rotating neutron star models with a parametrized rotation profile. Astron & Astrophysics (2012) 541:A156. 10.1051/0004-6361/201016316
11.
HanauskeMTakamiKBovardLRezzollaLFontJAGaleazziFet alRotational properties of hypermassive neutron stars from binary mergers. Phys Rev D (2017) 96:043004. 10.1103/physrevd.96.043004
12.
MarquesMOertelMHempelMNovakJ. New temperature dependent hyperonic equation of state: application to rotating neutron star models and i-q relations. Phys Rev C (2017) 96:045806. 10.1103/physrevc.96.045806
13.
CamelioGDietrichTMarquesMRosswogS. Rotating neutron stars with nonbarotropic thermal profile. Phys Rev D (2019) 100:123001. 10.1103/physrevd.100.123001
14.
WeberF. Pulsars as astrophysical laboratories for nuclear and particle physics series in high energy physics, cosmology and gravitation. FL, USA: CRC Press (1999). 10.1201/9780203741719
15.
GlendenningN. Compact stars: nuclear physics particle physics and general relativity. Astronomy and astrophysics library (2012).
16.
SpinellaWMWeberF. Dense baryonic matter in the cores of neutron stars. In: Topics on strong gravity: a modern view on theories and experiments (world scientific) (2020). p. 85–152.
17.
WeberFFarrellDSpinellaWMMalfattiGOrsariaMGContreraGAet alPhases of hadron-quark matter in (proto) neutron stars. Universe (2019) 5:169. 10.3390/universe5070169
18.
MalfattiGOrsariaMGContreraGAWeberFRanea-SandovalIF. Hot quark matter and (proto-) neutron stars. Phys Rev C (2019) 100:015803. 10.1103/PhysRevC.100.015803
19.
MalfattiGOrsariaMGRanea-SandovalIFContreraGAWeberF. Delta baryons and diquark formation in the cores of neutron stars. Phys Rev D (2020) 102:063008. 10.1103/PhysRevD.102.063008
20.
SedrakianALiJJWeberF. Heavy baryons in compact stars. Prog Part Nucl Phys (2023) 475:104041. 10.1016/j.ppnp.2023.104041
21.
LattimerJMSteinerAW. Constraints on the symmetry energy using the mass-radius relation of neutron stars. Eur Phys J A (2014) 50:40–24. 10.1140/epja/i2014-14040-y
22.
LattimerJM. Neutron star mass and radius measurements. Universe (2019) 5:159. 10.3390/universe5070159
23.
StoneJR. Nuclear physics and astrophysics constraints on the high density matter equation of state. Universe (2021) 7:257. 10.3390/universe7080257
24.
PoschenriederPWeigelM. Nuclear matter problem in the relativistic green’s function approach. Phys Rev C (1988) 38:471–86. 10.1103/physrevc.38.471
25.
HorowitzCSerotBD. The relativistic two-nucleon problem in nuclear matter. Nucl Phys A (1987) 464:613–99. 10.1016/0375-9474(87)90370-8
26.
ZuoWBombaciILombardoU. Asymmetric nuclear matter from an extended brueckner-Hartree-Fock approach. Phys Rev C (1999) 60:024605. 10.1103/PhysRevC.60.024605
27.
WangSZhaoQRingPMengJ. Nuclear matter in relativistic brueckner-Hartree-Fock theory with bonn potential in the full Dirac space. Phys Rev C (2021) 103:054319. 10.1103/physrevc.103.054319
28.
LovatoABombaciILogotetaDPiarulliMWiringaRB. Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials. Phys Rev C (2022) 105:055808. 10.1103/PhysRevC.105.055808
29.
TongHWangCWangS. Nuclear matter and neutron stars from relativistic brueckner–Hartree–Fock theory. Astrophysical J (2022) 930:137. 10.3847/1538-4357/ac65fc
30.
FarrellDWeberF. Relativistic brueckner–Hartree–Fock calculations for cold and hot neutron stars. Astrophysical J (2024) 969:49. 10.3847/1538-4357/ad47b8
31.
FaberJARasioFA. Binary neutron star mergers. Living Rev Relativity (2012) 15:8–83. 10.12942/lrr-2012-8
32.
BaiottiLRezzollaL. Binary neutron star mergers: a review of einstein’s richest laboratory. Rep Prog Phys (2017) 80:096901. 10.1088/1361-6633/aa67bb
33.
RaithelCAÖzelFPsaltisD. Finite-temperature extension for cold neutron star equations of state. Astrophysical J (2019) 875:12. 10.3847/1538-4357/ab08ea
34.
FarrellDAlpAWeberFSpinellaWMalfattiGOrsariaMGet alHot neutron star matter and proto-neutron stars. World Scientific (2023) 199–259. 10.1142/9789811220913_0005
35.
MartinPCSchwingerJ. Theory of many-particle systems. i. Phys Rev (1959) 115:1342–73. 10.1103/physrev.115.1342
36.
BrockmannRMachleidtR. Relativistic nuclear structure. i. nuclear matter. Phys Rev C (1990) 42:1965–80. 10.1103/physrevc.42.1965
37.
WangCWangSTongHHuJYaoJ. Low-momentum relativistic nucleon-nucleon potentials: nuclear matter. Phys Rev C (2024) 109:034002. 10.1103/physrevc.109.034002
38.
DolanLJackiwR. Symmetry behavior at finite temperature. Phys Rev D (1974) 9:3320–41. 10.1103/PhysRevD.9.3320
39.
ter HaarBMalflietR. Equation of state of nuclear matter in the relativistic Dirac-brueckner approach. Phys Rev Lett (1986) 56:1237–40. 10.1103/physrevlett.56.1237
40.
KomatsuHEriguchiYHachisuI. Rapidly rotating general relativistic stars–i. numerical method and its application to uniformly rotating polytropes. Monthly Notices R Astronomical Soc (1989) 237:355–79. 10.1093/mnras/237.2.355
41.
StergioulasN. Rotating stars in relativity. Living Rev Relativity (2003) 6:3–109. 10.12942/lrr-2003-3
42.
ShapiroSL. Differential rotation in neutron stars: magnetic braking and viscous damping. Astrophysical J (2000) 544:397–408. 10.1086/317209
43.
MorrisonIABaumgarteTWShapiroSL. Effect of differential rotation on the maximum mass of neutron stars: realistic nuclear equations of state. Astrophysical J (2004) 610:941–7. 10.1086/421897
44.
ZhouETsokarosAUryuKXuRShibataM. Differentially rotating strange star in general relativity. Phys Rev D (2019) 100:043015. 10.1103/physrevd.100.043015
45.
ChandrasekharS. Ellipsoidal figures of equilibrium. New York: Dover (1987).
46.
ShapiroSLTeukolskySA. Black holes, white dwarfs, and neutron stars: the physics of compact objects. John Wiley & Sons (2008).
47.
CutlerCLindblomL. Post-Newtonian frequencies for the pulsations of rapidly rotating neutron stars. Astrophysical J (1992) 385:630–41. 10.1086/170968
48.
HunterC. On secular stability, secular instability, and points of bifurcation of rotating gaseous masses. Astrophysical J (1977) 213:497–517. 10.1086/155181
49.
PaschalidisVStergioulasN. Rotating stars in relativity. Living Rev Relativity (2017) 20:7–169. 10.1007/s41114-017-0008-x
50.
SaijoMShibataMBaumgarteTWShapiroSL. Dynamical bar instability in rotating stars: effect of general relativity. Astrophysical J (2001) 548:919–31. 10.1086/319016
51.
BodenheimerPOstrikerJP. Rapidly rotating stars. viii. zero-viscosity polytropic sequences. Astrophysical J (1973) 180:159–70. 10.1086/151951
52.
ShibataMBaumgarteTWShapiroSL. The bar-mode instability in differentially rotating neutron stars: simulations in full general relativity. Astrophysical J (2000) 542:453–63. 10.1086/309525
53.
ShibataMKarinoSEriguchiY. Dynamical instability of differentially rotating stars. Monthly Notices R Astronomical Soc (2002) 334:L27–L31. 10.1046/j.1365-8711.2002.05724.x
54.
CamardaKDAnninosPFragilePCFontJA. Dynamical bar-mode instability in differentially rotating magnetized neutron stars. Astrophysical J (2009) 707:1610–22. 10.1088/0004-637x/707/2/1610
55.
Di GiovanniFSanchis-GualNCerdá-DuránPZilhaoMHerdeiroCFontJAet alDynamical bar-mode instability in spinning bosonic stars. Phys Rev D (2020) 102:124009. 10.1103/physrevd.102.124009
56.
RingP. Relativistic mean field theory in finite nuclei. Prog Part Nucl Phys (1996) 37:193–263. 10.1016/0146-6410(96)00054-3
57.
AkmalAPandharipandeVRavenhallD. Equation of state of nucleon matter and neutron star structure. Phys Rev C (1998) 58:1804–28. 10.1103/physrevc.58.1804
58.
EpelbaumEHammerH-WMeissnerU-G. Modern theory of nuclear forces. Rev Mod Phys (2009) 81:1773–825. 10.1103/RevModPhys.81.1773
59.
SammarrucaFChenBCoraggioLItacoNMachleidtR. Dirac-brueckner-Hartree-Fock versus chiral effective field theory. Phys Rev C—Nuclear Phys (2012) 86:054317. 10.1103/physrevc.86.054317
60.
ShenSLiangHMengJRingPZhangS. Fully self-consistent relativistic brueckner-Hartree-Fock theory for finite nuclei. Phys Rev C (2017) 96:014316. 10.1103/physrevc.96.014316
61.
HebelerKSchwenkA. Chiral three-nucleon forces and neutron matter. Phys Rev C (2010) 82:014314. 10.1103/physrevc.82.014314
62.
TewsIKrügerTHebelerKSchwenkA. Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory. Phys Rev Lett (2013) 110:032504. 10.1103/physrevlett.110.032504
63.
SteinerAHempelMFischerT. Core-collapse supernova equations of state. Astrophysical J (2013) 774:17. 10.1088/0004-637x/774/1/17
64.
BaumgarteTShapiroS. On the maximum mass of differentially rotating neutron stars. Astrophysical J (2003) 585:921–8. 10.1086/312425
65.
LattimerJPrakashM. Neutron star structure and the equation of state. Astrophysical J (2001) 550:426–42. 10.1086/319702
66.
LattimerJMPrakashM. Neutron star observations: progress and problems. Phys Rep (2007) 442:109–65. 10.1016/j.physrep.2007.02.003
67.
ChirentiCLRezzollaL. The merger of two neutron stars: gravitational waves and electromagnetic signals. Astrophysical J (2012) 760:150. 10.1088/0004-637X/760/2/150
68.
MetzgerBDBergerE. Kilonovae. Annu Rev Astron Astrophysics (2012) 50:503–46. 10.1146/annurev-astro-081811-125339
69.
ZhangB. The origin of long-duration gamma-ray bursts. Nat Phys (2016) 12:737–41. 10.1038/nphys3805
70.
AbbottBPAbbottRAbbottTAcerneseFAckleyKAdamsCet alGw170817: observation of gravitational waves from a binary neutron star inspiral. Phys Rev Lett (2017) 119:161101. 10.1103/PhysRevLett.119.161101
71.
KasenDMetzgerBBarnesJQuataertERamirez-RuizEKilonova from the ejecta of a neutron star merger. Nature (2017) 551:80–84. 10.1038/nature24063
72.
MargalitBMetzgerBD. Constraining the maximum mass of neutron stars from multi-messenger observations of GW170817. Astrophysical J Lett (2017) 850:L19. 10.3847/2041-8213/aa991c
73.
PiroAL. The role of neutron star merger kinematics in explaining the diversity of short gamma-ray bursts. Astrophysical J (2019) 885:10. 10.3847/1538-4357/ab3f99
74.
OrsariaMGMalfattiGMarianiMRanea-SandovalIFGarcíaFSpinellaWMet alPhase transitions in neutron stars and their links to gravitational waves. J Phys G: Nucl Part Phys (2019) 46:073002. 10.1088/1361-6471/ab1d81
75.
AasiJAbbottBAbbottRAbernathyMRAckleyKAdamsCet alAdvanced ligo. Class Quan Gravity (2015) 32, 074001, 10.1088/0264-9381/32/7/074001
76.
AcerneseFAgathosMAgatsumaKAisaDAllemandouNAlloccaAet alAdvanced virgo: a second-generation interferometric gravitational wave detector. Classical Quan Gravity (2015) 32:024001. 10.1088/0264-9381/32/2/024001
77.
AsoYMichimuraYSomiyaKAndoMMioNSekiguchiTet alKagra: 2.5 generation interferometric gravitational wave detector. Phys Rev D (2020) 88:043007. 10.1103/physrevd.88.043007
78.
PunturoMAbernathyMAcerneseFAllenBAnderssonNArunKet alThe einstein telescope: a third-generation gravitational wave observatory. Classical Quan Gravity (2010) 27:194002. 10.1088/0264-9381/27/19/194002
79.
WeisskopfMTananbaumHVan SpeybroeckLO’DellSThe chandra x-ray observatory: exploring the high-energy universe, 112. WA, USA: Publications of the Astronomical Society of the Pacific (2000). 586.
80.
JansenFLumbDAltieriBClavelJEhleMErdCet alThe xmm-Newton observatory. Astron & Astrophysics (2001) 365. 10.1051/0004-6361:20000036
81.
PerleyRChandlerCButlerBWrobelJ. The expanded very large array: a new telescope for new science. Astrophysical J Lett (2011) 739:L1. 10.1088/2041-8205/739/1/l1
82.
BraunRBourkeTGreenJAKeaneEWaggJThe square kilometre array. Proc Sci AASKA14 (2015) 174.
83.
NanRLiDJinCWangQZhuYZhuWet alThe five-hundred-meter aperture spherical radio telescope (fast) project. Int J Mod Phys D (2011) 20:989–1024. 10.1142/s0218271811019335
Summary
Keywords
neutron star, differential rotation, equation of state, brueckner-Hartree-Fock, finite temperature field theory, bar mode instability
Citation
Farrell D, Weber F and Negreiros R (2024) Differential rotation in neutron stars at finite temperatures. Front. Phys. 12:1474615. doi: 10.3389/fphy.2024.1474615
Received
01 August 2024
Accepted
30 September 2024
Published
22 October 2024
Volume
12 - 2024
Edited by
James Lattimer, Stony Brook University, United States
Reviewed by
Jacobo Ruiz De Elvira, Complutense University of Madrid, Spain
Jiangming Yao, Zhuhai Campus, China
Updates
Copyright
© 2024 Farrell, Weber and Negreiros.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Delaney Farrell, dfarrell@sdsu.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.