Next Article in Journal
Spin Susceptibility in Neutron Matter from Quantum Monte Carlo Calculations
Previous Article in Journal
Direct Photons from Hot Quark Matter in Renormalized Finite-Time-Path QED
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutron Star Cooling Within the Equation of State With Induced Surface Tension

CFisUC, Department of Physics, University of Coimbra, 3004-516 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 11 September 2020 / Revised: 14 November 2020 / Accepted: 19 November 2020 / Published: 21 November 2020

Abstract

:
We study the thermal evolution of neutron stars described within the equation of state with induced surface tension (IST) that reproduces properties of normal nuclear matter, fulfills the proton flow constraint, provides a high-quality description of hadron multiplicities created during the nuclear-nuclear collision experiments, and it is equally compatible with the constraints from astrophysical observations and the GW170817 event. The model features strong direct Urca processes for the stars above 1.91 M . The IST equation of state shows very good agreement with the available cooling data, even without introducing nuclear pairing. We also analysed the effect of the singlet proton/neutron and triplet neutron pairing on the cooling of neutron stars of different mass. We show that the description of the compact object in the center of the Cassiopeia A does not necessarily require an inclusion of neutron superfluidity and/or proton superconductivity. Our results indicate that data of Cassiopeia A can be adequately well reproduced by a 1.66 M star with an atmosphere of light elements. Moreover, the IST EoS reproduces each of the observational datasets for the surface temperature of Cassiopeia A either by a rapidly cooling ∼ 1.955 M star with paired and unpaired matter or by a 1.91 M star with the inclusion of neutron and proton pairings in the singlet channel.

1. Introduction

Born out of supernova explosions, neutron stars (NSs) are considered to start their life having very high internal temperatures and cool down through a combination of thermal radiation from their surface and neutrino emission from their interior. From the first day of their lives, when the temperature of their interior has already dropped from ∼ 10 11 K to ∼ 10 9 10 K making it transparent to neutrinos, up until the first million years of their existence, thermal energy is mainly carried away in the form of neutrino radiation. During this time, measurements of surface temperature and luminosity of the stars can provide significant information about the properties of matter in their depth, since the thermal evolution depends on factors, such as the internal composition and the thermodynamic properties of matter, are defined by the Equation of State (EoS), the chemical abundances of the envelope, and the type of pairing between the constituent particles. Therefore, it is necessary for any theoretical calculations aiming to describe the cooling process of NSs to consider the various combinations between those factors. The set of cooling curves extracted from these simulations can then be tested against the observable features of NSs in the X-ray part of the spectrum.
The particle composition of the core of NSs plays a prominent role in their cooling process, since it is the decisive factor of whether the direct Urca (DU) process of neutron β -decay and its inverse process are allowed to occur in the interior of the star [1]. These are the most efficient neutrino-emitting reactions that can happen in hadronic matter and, when permitted, they lead to a rapid cooling of the NS, substantially different from the case they are forbidden [2,3]. For these reactions to be possible, the Fermi momenta of the participating particles must satisfy the kinematic restriction of triangle inequality. While taking charge neutrality and the relation between the Fermi momenta and the number density of each particle into account, this constraint leads to the requirement that the proton fraction should be higher than ∼ 11 % of the total baryon density [2].
The composition and thermodynamic properties of compact stars were described by the EoS formulated within the framework of induced surface tension (IST). Originally, this approach was developed in order to correctly account for the short-range repulsion between the particles characterized by the different hard-core radii, which are known to be principal elements of hadron EoS at high temperatures typical for collisions of heavy ions [4]. One of the main advantages of the IST EoS compared to other approaches considering the hard-core repulsion is its ability to correctly account for the excluded volume of particles reflected in correct values of first four virial coefficients of hard spheres [5]. Thus, the IST EoS goes well beyond the usual van der Waals approximation, which increases the range of its causality and applicability. This allowed the present model to describe the nuclear matter phase diagram at low temperatures and intermediate densities, the region of the nuclear liquid-gas phase transition and its critical end-point [6]. In principle, the present model, in comparison to the Skyrme and Relativistic Mean Field models, can be extended for an arbitrarily large number of different particle species, as was done in Ref. [5]. This feature of the present approach also allowed it to successfully describe hadron yields that were measured in heavy-ion collision experiments, i.e., AGS (Alternating Gradient Synchrotron, Brookhaven National Laboratory), SPS (Super Proton Synchrotron, CERN), RHIC (Relativistic Heavy Ion Collider, Brookhaven National Laboratory), and LHC (Large Hadron Collider, CERN) [7,8]. Extensive analysis of nuclear-nuclear collision experimental data with the multicomponent IST EoS helped to constraint the repulsive interaction between the particles, described by the hard-core radii. The successful application of the IST approach in a wide range of temperatures and densities demonstrated its potential towards formulation of a unified EoS of hadron matter. Further development of this model led to incorporation of the long-range interaction within the mean-field framework [9], which allowed for the IST EoS to fulfill the flow constraint of Danielewicz et al. [10]. Recently, the present model was applied to the modeling of nuclear matter at very high density and vanishing temperature existing inside the compact astrophysical objects [11,12,13]. The present study of the thermal evolution of NSs is the next step towards constraining the IST EoS and developing a unified hadron EoS.
For a realistic description of the outer layers of the NS the IST EoS is supplemented by the Haensel-Zdunik (HZ) EoS for the outer crust and the Negele–Vautherin (NV) EoS for the inner crust [14,15].
An additional important factor worth considering is the superfluidity (superconductivity) of neutrons (protons) in the NS interior [16], since they are able to alter the thermal evolution of the star [17]. Despite suppressing the rates of neutrino emission from the rest of the processes, pairing between protons or neutrons introduces a new ν e emitting mechanism, due to the permanent formation and breaking of Cooper pairs, known as PBF [3]. This neutrino emitting process is activated in each baryon species as soon as the temperature reaches the respective critical temperatures T c , n , T c , p . Nevertheless, the details, such as the onset T, the form of the gap parameter, and the peak of the PBF emissivity, are strongly sensitive to the gap model that was adopted to describe the matter in this supercritical state. As a result, applying different models to study the NS cooling provides also the possibility of gaining a further insight on the behavior of paired matter.
One of the most intriguing and debatable object that serves as a test ground for many models of nucleon pairing and alternative approaches is Cassiopeia A (Cas A). This compact source being a very young (at an age of ≈340 yr) NS demonstrates a rapid cooling across a decade of observations. Thus, since the early observations by Heinke&Ho (2010), the reported surface temperature change is decreasing from 3.6 ± 0.6 % [18] to 3.5 ± 0.4 % [19], recently refined to < 2.4 % (< 3.3 % ) for constant (varying) hydrogen absorption [20], and 2.1 ± 0.2 % ( 2.7 ± 0.3 % ) for constant (varying) hydrogen absorption [21]. Moreover, as it was reported by Ho & Heinke (2009) [22], the Cas A X-ray spectra suggest a carbon composition of the atmosphere and a mass in a range between 1.4 M and 2.4 M for the NS in question. An interpretation of such a rapid cooling of the Cas A includes onset of superfluidity (superconductivity) in the core of the star [23,24,25,26], nuclear medium cooling scenario [27,28], phase transition to quark-gluon plasma [29,30], decay of magnetic field [31], heating by r-mode oscillations [32], etc. Unfortunately, for the time being, available observational data do not provide us with a possibility to differentiate between these scenarios.
Based on the IST EoS, in this work we study the thermal evolution of NSs, aiming to describe the available cooling data. Our investigation is focused on the implications of employing various gap models in order to describe neutron superfluidity and proton superconductivity. In addition, we examine the effect of envelope composition on the cooling process of compact stars. Modeling the thermal evolution of NSs was performed by using the NSCool code [2,33].
The paper is organized, as follows: in Section 2, we present a brief description of the model. In Section 3, we discuss the neutrino emission mechanisms taking place in a NS as well as the pairing models used in the simulations. Section 4 is dedicated to the results of the calculations, while Section 5 contains the summary and the conclusions.

2. Equation of State

In the present formulation, the IST EoS includes neutrons, protons, and electrons, assuming that an effect of heavy baryons and mesons, which can appear at high densities, is absorbed in the mean-field potential. Thus, the model accounts for the strong short-range repulsion and relatively weak long-range attraction between nucleons, while electrons are treated as an ideal Fermi gas. The former part of the nucleon-nucleon interaction is modelled with the hard-core radius, similarly to the famous Van der Waals EoS. The hard-core repulsion of nucleons leads to an appearance of an excluded volume. However, contrary to the Van der Waals approximation, the IST EoS instead of the constant excluded volume has a density-dependent one. The key element of the model is the IST coefficient. It accounts not only for the correct values of four virial coefficients of hard spheres, but it also extends the causality range of the model to the density range typical for the NSs [6,11]. Moreover, the long-range attraction and asymmetry between neutrons and protons are accounted via the mean-field potentials, the parameters of which were fitted to the properties of matter at saturation density, and, therefore, effectively accounts for many particles interaction of different species (for details see [9]). More detailed information regarding the model and its application to the NSs can be found in Refs. [6,11,12].
We adopt the softer parameterisation of the IST EoS that corresponds to the set B of Ref. [34] (for simplicity we keep referring to this parameterisation of the model as set B). The set B was chosen due to the better agreement with the astrophysical constraints and results coming from the GW170817 [35]. It gives the values of the symmetry energy E s y m = 30.0 MeV, symmetry energy slope L = 93.2 MeV and nuclear incompressibility factor K 0 = 201.0 MeV at normal nuclear density that is in full agreement with the present nuclear matter results [11,12]. As shown in Figure 1, the maximum mass M m a x = 2.08 M is consistent with the recent measurements of the most massive NSs, i.e., PSR J0348+0432 [36] and PSR J0740+6620 [37].
The relatively high slope of the symmetry energy of the present model requires a special comment. The analysis of correlation between this quantity and the radius of a 1.4 M NS performed within a wide set of Skyrme and Relativistic mean-field (RMF) models shows that L > 90 MeV is unlikely to be consistent with R 1.4 below 12 km [41,42,43,44]. However, as seen from Figure 1, the present model yields R 1.4 = 11.39 km. In order to explain this apparent discrepancy, we have to mention that a stellar characteristic, such as the NS radius, is defined by the total stiffness of an EoS. The considered model contains an isospin-independent term and a term related to the nuclear symmetry energy. While high values of L lead to stiff symmetry energy, this redundant stiffness can be moderated by the soft isospin-independent part of EoS, which is the case for the IST EoS. In order to demonstrate this, we compare EoS of symmetric nuclear matter calculated within the present model and the well known RMF model NL3 [45,46]. As shown on Figure 2, the symmetric IST EoS is indeed not only significantly softer than the symmetric NL3 model, but it also allows for the present model to fit the proton flow constraint [10] (the shaded gray area on Figure 2). Further, we compare the electrically neutral IST EoS with the three parameterizations of RMF model N L 3 ω ρ , which implies coupling of ω and ρ mesons leading to different slopes of the symmetry energy [47]. The values of L = 55 MeV, 88 MeV and 118 MeV (green, blue and magenta dashed curves, respectively) label the N L 3 ω ρ models. It is seen from Figure 2 that, despite its large L, the neutral IST EoS is actually rather soft, which makes it consistent with R 1.4 < 12 km.

3. Cooling Processes

The thermal evolution of NSs can be divided in two stages. During the first one, which is known as neutrino cooling era, the ν e emission generated from a plethora of emission mechanisms throughout the whole interior of the star dominates the cooling process [2,3]. On a timescale of a million years after the NS formation, when T c o r e has dropped below 10 8 K , neutrino emission from the core is less efficient due to the strong temperature dependence, and photon emission from the surface overtakes as the leading heat loss mechanism. Beyond this point, the neutrino cooling era is over and the photon cooling era begins. This shift is marked by an abrupt decline of the total luminosity and surface temperature of the star, with more massive stars exhibiting steeper drops of those quantities than less massive ones, since photon luminosity follows the typical blackbody radiation law.
During the neutrino cooling era, the leading neutrino-generating process varies over the regions of the star and over time. Overall, the main factors regulating the neutrino emissivity of each process are density, temperature, and the existing degree of Cooper pairing between particles. For example, in the envelope of the star a pair annihilation is the most productive mechanism, while in the outer crust neutrino emission is dominated by plasmon decay until its temperature reaches a few times 10 8 K , replaced by electron-ion bremsstrahlung beyond that point [3,17]. The latter remains the most efficient energy loss mechanism throughout the inner crust as well, as long as the neutrons of this region are paired [48]. Even in their supercritical state though, inner crust neutrons contribute to the total neutrino emission through the PBF process over a narrow range of densities and temperatures [3].
In the core of the star, apart from the bremsstrahlung processes between the free particles, the main neutrino-emitting processes that can occur are the DU process of the neutron β -decay and its inverse:
n p + e + ν ¯ e , p + e n + ν e .
The realization of these reactions depends on the proton fraction in the stellar interior, as mentioned above. For the IST EoS adopted in this work, these fast processes can only proceed in the core of stars with central densities n c higher than n D U = 0.862 f m 3 . For NSs with lower n c , where the Fermi momenta of the involved particles fail to satisfy the kinematic restriction of p F , n p F , p + p F , e , neutron β -decay and its inverse can still proceed with the help of a spectator proton or neutron that provides the extra p F , i needed for the conservation of momentum. In this case, though, the efficiency of these so-called modified Urca (MU) processes in removing heat from the star is lower than that of the DU ones, resulting in an overall slower cooling of the star. Nevertheless, neutrino emission from MU processes surpasses the rest of the emissivities in the core, unless n and p are in a paired state, which, in turn, results in a further suppression of the neutrino emission rates. However, pairing partially compensates for the delay in the cooling of the NS core as it introduces PBF, an additional efficient channel for carrying the heat away in the form of neutrino-antineutrino pairs [3].
Cooper pairing of core protons and crust neutrons sets in a few years after the NS birth, while core neutrons pairing is viable at a later time. As regards the types of pairing, all of them follow the standard BCS pattern, with free neutrons of the inner crust and protons of the core undergoing singlet-state pairing ( 1 S 0 ) , while neutrons of the core are expected to pair in the triplet-state ( 3 P 2 ) [2,49,50]. Although a complete, precise description of the effect for densities that are relevant to NS matter is still pending, it is widely applied for investigating its imprint on the thermal evolution of compact stars. The transition of neutrons (protons) to superfluid (superconducting) state also inflicts a suppression on their specific heat. Most studies agree that, once it switches on for a certain species, c υ , i is reduced by a factor R e Δ / T , where Δ is the gap parameter, linked to the respective critical temperature T c via the standard BCS pairing relation T c 0.57 Δ [2,49,50]. A further drop in temperature causes a further reduction of heat capacity to the point of being equal to that of leptons [3]. Concerning the neutrino emissivity of processes that involve paired baryons, suppression is also induced, since particles in such a supercritical state have to overcome the energy gap, in order to interact with another particle. In the range of temperatures T T c , this behavior is numerically illustrated while using a yet another set of control functions R χ T / T c , which differs between each particle species χ [17].
The onset temperatures T c , n ( T c , p ) of superfluidity (superconductivity), the associated Δ , and the profile of PBF neutrino emissivity vary over the several models developed to describe the paired matter. Therefore, when including nucleon pairing for simulating the thermal evolution of NSs, the effect is subject to the gap models employed. Our choice for the simulations were the following models: SFB [51] for the 1 S 0 superfluidity of neutrons, T72 [52] and AO [53] for their 3 P 2 superfluidity, and AO [54], CCDK [55] for the proton 1 S 0 channel in different combinations. The critical temperatures T c were calculated according to the phenomenological formula that was suggested by Kaminker et al. and the parameters used by Ho et al. [48,56]. It has to be to noted, though, that the vertex corrections were omitted in the study.
The relation between the critical temperature and the baryon density for all of the adopted gap models is shown on Figure 3. According to the SFB model, 1 S 0 pairing of crust neutrons appears first at densities around 0.15 n 0 , once the crust temperature has dropped below 5 · 10 9 K . As the crust keeps cooling, the region of superfluid crust neutrons expands both outwards and towards the crust-core interface. Regarding the superconductivity of protons in the core, the CCDK model suggests that the proton pairing occurs initially at a temperature not greater than 7 · 10 9 K in the layer that corresponds to 2 n 0 . On the contrary, the AO model proposes that the onset temperature is a few times lower, at 2 · 10 9 K , while the associated baryon density only slightly lower. Finally, as claimed by the T72 model, the threshold for the onset of 3 P 2 superfluidity of core neutrons is below 10 9 K and pairing occurs in a narrow, symmetric region centered around the normal nuclear density n 0 . The respective AO model, on the other hand, although it does not differ considerably in terms of T c , n , implies that neutron pairing can occur throughout the whole region of the core.
The composition of the envelope is another deciding factor for the surface photon luminosity of the star, which is ultimately the quantity of observational importance. Namely, heavier elements tend to delay heat transport from the outer crust to the surface, since, in this case, the electron thermal conductivity is reduced [57]. The envelope acts as a blanket that, due to its poor thermal conductivity, creates a temperature gradient between the surface and interior of the star. The standard approach of thermal evolution codes is to set an envelope model as an outer boundary condition that links the temperature at the bottom of the envelope ( T b ) to the temperature of the stellar surface ( T s ) . The fundamental assumptions behind this method are that the envelope has a thermal relaxation timescale that is much shorter than that of the crust, and that the neutrino emissivity in the envelope is negligible [2]. In this work, we used two distinct envelope models: one composed of heavy elements and a hydrogen-rich one that contains the fraction of light elements η = Δ M / M = 10 7 [58]. The presence of light elements (such as H) in the outer layers may suggest an accretion onto the NS after its formation, elements between He and C could occur as a result of diffusive nuclear burning on the surface of the star, while an existence of much heavier elements (up to Fe) give a hint of no accretion or the thermonuclear reactions after it [21].

4. Results

At the beginning, we focused on the thermal evolution of NSs without any sort of pairing between the nucleons. As you can see on Figure 4, all of the cooling curves that are depicted in color exhibit a slow cooling due to the domination of the MU processes, since the DU ones are not kinematically allowed. Finally, at central densities of over n D U = 0.862 f m 3 of beta-stable and charge neutral matter, which correspond to NSs with masses M 1.91 M , the DU processes are switched on in the core of the star, so that it undergoes enhanced cooling (see the black curves on Figure 4). In order to model the uncertainties of the heat-blanketing effect of the envelope, we compare the thermal evolution of NSs with a non-accreted envelope containing heavy elements (dashed curves on Figure 4) with the envelope containing light elements (solid curves on Figure 4). Remarkably, the obtained cooling curves for unpaired matter describe the experimental data very good.
Furthermore, we studied the effect of neutron superfluidity and proton superconductivity on the thermal evolution of NSs in two stages. First, we considered n 1 S 0 superfluidity with the SFB model [51] together with a p 1 S 0 superconductivity that was described by AO [54] and CCDK [55] models (see Figure 5). For both considered combinations of gaps, i.e., SFB + AO (see the left panel of Figure 5) and SFB + CCDK (see the right panel of Figure 5), the IST EoS is in very good agreement with the observational data. As in the non-superfluid case, the model predicts the surface temperature for the Cas A both with the fast cooling curve for M D U 2.0 M , as well as with the curves for low-mass stars with a heavy-elements envelope. For the SFB + CCDK gaps, the Cas A is only described with the curve for M = 1.8 M also with a heavy-elements envelope.
Finally, we studied the effect of the neutron triplet pairing on the cooling of NSs by considering a shallow gap with a maximum of the critical temperature at the saturation density n 0 (T72 model [52]) and an extended gap with a maximum of T c at 2 n 0 (AO model [54]). As it shown in Figure 6, the combination of models SFB (for n 1 S 0 ) + T72 (for n 3 P 2 )+CCDK (for p 1 S 0 ) (left panel) and SFB (for n 1 S 0 ) + AO (for n 3 P 2 ) + CCDK (for p 1 S 0 ) (right panel) give a qualitatively similar results. Adding the n 3 P 2 pairing in the core of NSs leads to more rapid cooling and makes it incompatible with most of the observational data. Therefore, we conclude that, within our model, neutron pairing in the triplet channel is inconsistent with the observational data.

Description of the Cas A Temperature Evolution

The fit results of the latest Cas A observational data within the IST EoS are presented on Figure 7. Namely, we find that the three data points (brown diamonds on Figure 7) reported by Posselt& Pavlov [20] can be described by a rapidly cooling 1.91 M star with light-elements atmosphere (blue dashed curve). Despite a slightly slower cooling rate, the medium-mass 1.66 M star with light-elements envelope (green dashed curve) indicates a good agreement with the data. Another set of 13 Chandra ACIS-S Graded data points for constant interstellar absorption column [21] is depicted as black circles on Figure 7. The red squares correspond to three points measured by Posselt& Pavlov [20], normalized by the observations of 2015 [60]. As far as one can see, these three points together with Chandra ACIS-S Graded data can be simultaneously fitted with a fast cooling massive star of 1.96 M with unpaired matter in its interior (red dashed curve), as well as with a 1.955 M star with inclusion of neutron and proton pairings in the singlet channel (purple solid curve). Both of these high-mass stars have switched-on DU processes that lead to such a fast temperature decrease. Therefore, we can conclude that the IST EoS is able to reproduce the apparent fast cooling of Cas A, even without the inclusion of neutron superfluidity and proton superconductivity.

5. Conclusions

We studied the cooling of NSs with the novel IST EoS that was previously applied to the analysis of nuclear matter properties, heavy-ion collision experimental data, and was recently generalized for the description of matter inside NSs. The considered model parameterisation is in good agreement with the GW170817 and the recent measurements of the most massive NSs. For the stars with M 1.91 M the model allows for the occurrence of DU processes that leads to much faster cooling. An effect of neutron pairing in 1 S 0 and 3 P 2 channels, as well as proton pairing in 1 S 0 channel, is thoroughly analysed with different superfluidity/superconductivity models. In addition, the effects of different types of envelopes on the thermal evolution of NSs were also compared. We show that, for unpaired matter, the obtained cooling curves are fully consistent with the currently existing observational data.
We report that an explanation of the rapid temperature drop of Cas A does not necessarily require an inclusion of neutron superfluidity and/or proton conductivity in the model. The surface temperature of Cas A reported in Refs. [21,60] is reproduced, within the IST EoS, by the cooling curves for medium and high-mass stars. Thus, a medium-mass 1.66 M NS with an envelope of light elements offers a reasonably good description and, at the same time, is the most congruous with the observational results. Alternatively, the data can be explained by the unpaired matter with a rapidly cooling 1.96 M star, as well as with a 1.955 M star (light-elements envelope) with included proton and neutron pairing in the singlet channel. Another set of measurements of the surface temperature of Cas A [20] can be explained by a heavy-mass 1.91 M star with an allowed DU process.
The two considered models of the triplet neutron pairing in the core of the star (a shallow and an extended neutron superfluidity) lead to too rapid cooling of all NSs, such that the old stars ( t > 10 5 yr) cannot be reproduced. Thus, we can conclude that the calculations favour the vanishing neutron triplet pairing gap.
As a future step, we plan to extend the list of particles by including heavy baryons and mesons. Despite their small effect on EoS, these particles do modify the thermal evolution of NSs and, therefore, have to be considered in the cooling simulations.

Author Contributions

Both authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

V.S. acknowledges the support from the Fundação para a Ciência e Tecnologia (FCT) within the project UID/04564/2020.

Acknowledgments

We are thankful for the fruitful discussions and suggestions to M. Fortin, O. Ivanytskyi and C. Providência.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lattimer, J.M.; Pethick, C.J.; Prakash, M.; Haensel, P. Direct URCA process in neutron stars. Phys. Pev. Lett. 1991, 66, 2701. [Google Scholar] [CrossRef]
  2. Page, D.; Geppert, U.; Weber, F. The cooling of compact stars. Nucl. Phys. A 2006, 777, 497. [Google Scholar] [CrossRef]
  3. Potekhin, A.Y.; Pons, J.A.; Page, D. Neutron stars-cooling and transport. Sp. Sci. Rev. 2015, 191, 239. [Google Scholar] [CrossRef] [Green Version]
  4. Andronic, A.; Braun-Munzinger, P.; Stachel, J. Hadron production in central nucleus nucleus collisions at chemical freeze-out. Nucl. Phys. A 2006, 772, 167. [Google Scholar] [CrossRef] [Green Version]
  5. Bugaev, K.A.; Sagun, V.V.; Ivanytskyi, A.I.; Yakimenko, I.P.; Nikonov, E.G.; Taranenko, A.V.; Zinovjev, G.M. Going beyond the second virial coefficient in the hadron resonance gas model. Nucl. Phys. A 2018, 970, 133. [Google Scholar] [CrossRef]
  6. Sagun, V.V.; Ivanytskyi, A.I.; Bugaev, K.A.; Mishustin, I.N. The statistical multifragmentation model for liquid-gas phase transition with a compressible nuclear liquid. Nucl. Phys. A 2014, 924, 24. [Google Scholar] [CrossRef] [Green Version]
  7. Sagun, V.V.; Bugaev, K.A.; Ivanytskyi, A.I.; Yakimenko, I.P.; Nikonov, E.G.; Taranenko, A.V.; Greiner, C.; Blaschke, D.B.; Zinovjev, G.M. Hadron resonance gas model with induced surface tension. Eur. Phys. J. A 2018, 54, 100. [Google Scholar] [CrossRef] [Green Version]
  8. Sagun, V.V.; Bugaev, K.A.; Ivanytskyi, A.I.; Oliinychenko, D.R.; Mishustin, I.N. Effects of induced surface tension in nuclear and hadron matter. Eur. Phys. J. Web Conf. 2017, 137, 09007. [Google Scholar] [CrossRef] [Green Version]
  9. Ivanytskyi, A.I.; Bugaev, K.A.; Sagun, V.V.; Bravina, L.V.; Zabrodin, E.E. Influence of flow constraints on the properties of the critical endpoint of symmetric nuclear matter. Phys. Rev. C 2018, 97, 064905. [Google Scholar] [CrossRef] [Green Version]
  10. Danielewicz, P.; Lacey, R.; Lynch, W.G. Determination of the equation of state of dense matter. Science 2002, 198, 1592. [Google Scholar] [CrossRef] [Green Version]
  11. Sagun, V.V.; Lopes, I.; Ivanytskyi, A.I. The Induced surface tension contribution for the equation of state of neutron stars. Astrophys. J. 2019, 871, 157. [Google Scholar] [CrossRef] [Green Version]
  12. Sagun, V.; Lopes, I.; Ivanytskyi, A. Neutron stars meet constraints from high and low energy nuclear physics. Nucl. Phys. A 2019, 982, 883. [Google Scholar] [CrossRef]
  13. Sagun, V.V.; Lopes, I. Neutron stars: A novel equation of state with induced surface tension. Astrophys. J. 2017, 850, 75. [Google Scholar] [CrossRef] [Green Version]
  14. Haensel, P.; Zdunik, J.L. Non-equilibrium processes in the crust of an accreting neutron star. Astron. Astrophys. 1990, 229, 117. [Google Scholar]
  15. Negele, J.W.; Vautherin, D. Neutron star matter at sub-nuclear densities. Nucl. Phys. A 1973, 207, 298. [Google Scholar] [CrossRef]
  16. Flowers, E.; Ruderman, M.; Sutherl, P. Neutrino pair emission from finite-temperature neutron superfluid and the cooling of young neutron stars. Astrophys. J. 1976, 205, 541. [Google Scholar] [CrossRef]
  17. Yakovlev, D.G.; Kaminker, A.D.; Gnedin, O.Y.; Haensel, P. Neutrino emission from neutron stars. Phys. Rep. 2001, 354, 1. [Google Scholar] [CrossRef] [Green Version]
  18. Heinke, C.O.; Ho, W.C.G. Direct observation of the cooling of the Cassiopeia A neutron star. Astrophys. J. 2010, 719, L167. [Google Scholar] [CrossRef] [Green Version]
  19. Elshamouty, K.G.; Heinke, C.O.; Sivakoff, G.R.; Ho, W.C.G.; Shternin, P.S.; Yakovlev, D.G.; Patnaude, D.J.; David, L. Measuring the Cooling of the Neutron Star in Cassiopeia A with all Chandra X-Ray Observatory Detectors. Astrophys. J. 2013, 777, 22. [Google Scholar] [CrossRef] [Green Version]
  20. Posselt, B.; Pavlov, G.G. Upper limits on the rapid cooling of the central compact object in Cas A. Astrophys. J. 2018, 864, 135. [Google Scholar] [CrossRef]
  21. Wijngaarden, M.J.P.; Ho, W.C.G.; Chang, P.; Heinke, C.O.; Page, D.; Beznogov, M.; Patnaude, D.J. Diffusive nuclear burning in cooling simulations and application to new temperature data of the Cassiopeia A neutron star. Mon. Not. R. Astron. Soc. 2019, 484, 974. [Google Scholar] [CrossRef] [Green Version]
  22. Ho, W.C.G.; Heinke, O.C. A neutron star with a carbon atmosphere in the Cassiopeia A supernova remnant. Nature 2009, 462, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yakovlev, D.G.; Ho, W.C.G.; Shternin, P.S.; Heinke, C.O.; Potekhin, A.Y. Cooling rates of neutron stars and the young neutron star in the Cassiopeia A supernova remnant. Mon. Not. R. Astron. Soc. 2011, 411, 1977. [Google Scholar] [CrossRef] [Green Version]
  24. Shternin, P.S.; Yakovlev, D.G.; Heinke, C.O.; Ho, W.C.G.; Patnaude, D.J. Cooling neutron star in the Cassiopeia A supernova remnant: Evidence for superfluidity in the core. Mon. Not. R. Astron. Soc. 2011, 412, L108. [Google Scholar] [CrossRef] [Green Version]
  25. Page, D.; Prakash, M.; Lattimer, J.M.; Steiner, A.W. Rapid Cooling of the Neutron Star in Cassiopeia A Triggered by Neutron Superfluidity in Dense Matter. Phys. Rev. Lett. 2011, 106, 081101. [Google Scholar] [CrossRef] [Green Version]
  26. Taranto, G.; Burgio, F.G.; Schulze, H.-J. Cassiopeia A and direct Urca cooling. Mon. Not. R. Astron. Soc. 2016, 456, 1451–1458. [Google Scholar] [CrossRef] [Green Version]
  27. Blaschke, D.; Grigorian, H.; Voskresensky, D.N.; Weber, F. Cooling of the neutron star in Cassiopeia A. Phys. Rev. C 2012, 85, 022802. [Google Scholar] [CrossRef] [Green Version]
  28. Blaschke, D.; Grigorian, H.; Voskresensky, D.N. Nuclear medium cooling scenario in light of new Cas A cooling data and the 2.0 M pulsar mass measurements. Phys. Rev. C 2013, 88, 065805. [Google Scholar] [CrossRef] [Green Version]
  29. Sedrakian, A. Rapid cooling of Cassiopeia A as a phase transition in dense QCD. Astron. Astrophys. 2013, 555, L10. [Google Scholar]
  30. Noda, T.; Hashimoto, M.; Yasutake, N.; Maruyama, T.; Tatsumi, T.; Fujimoto, M. Cooling of compact stars with color superconducting phase in quark-hadron mixed phase. Astrophys. J. 2013, 765, 1. [Google Scholar] [CrossRef] [Green Version]
  31. Bonanno, A.; Baldo, M.; Burgio, G.F.; Urpin, V. The neutron star in Cassiopeia A: Equation of state, superfluidity, and Joule heating. Astron. Astrophys. 2013, 561, L5. [Google Scholar] [CrossRef] [Green Version]
  32. Yang, S.-H.; Pi, C.-M.; Zheng, X.-P. Rapid Cooling of the Neutron Star in Cassiopeia A and r-mode Damping in the Core. Astrophys. J. Lett. 2011, 735, L29. [Google Scholar] [CrossRef] [Green Version]
  33. Page, D.; Reddy, S. Dense matter in compact stars: Theoretical developments and observational constraints. Annu. Rev. Nucl. Part Sci. 2006, 56, 327. [Google Scholar] [CrossRef] [Green Version]
  34. Sagun, V.; Panotopoulos, G.; Lopes, I. Asteroseismology: Radial oscillations of neutron stars with realistic equation of state. Phys. Rev. D 2020, 101, 063025. [Google Scholar] [CrossRef] [Green Version]
  35. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. (LIGO and Virgo Collab.) GW170817: Observation of gravitational waves from a binary neutron star inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [Green Version]
  36. Antoniadis, J.; Freire, P.C.C.; Wex, N.; Tauris, T.M.; Lynch, R.S.; van Kerkwijk, M.H.; Kramer, M.; Bassa, C.; Dhillon, V.S.; Driebe, T.; et al. A massive pulsar in a compact relativistic binary. Science 2013, 340, 448. [Google Scholar] [CrossRef] [Green Version]
  37. Cromartie, H.T.; Fonseca, E.; Ransom, S.M.; Demorest, P.B.; Arzoumanian, Z.; Blumer, H.; Brook, P.R.; DeCesar, M.E.; Dolch, T.; Ellis, J.A.; et al. Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar. Nat. Astr. 2020, 4, 72. [Google Scholar] [CrossRef] [Green Version]
  38. Steiner, A.W.; Lattimer, J.M.; Brown, E.F. The equation of state from observed masses and radii of neutron stars. Astrophys. J. 2010, 722, 33. [Google Scholar] [CrossRef]
  39. Steiner, A.W.; Lattimer, J.M.; Brown, E.F. The neutron star mass-radius relation and the equation of state of dense matter. Astrophys. J. 2013, 765, L5. [Google Scholar] [CrossRef]
  40. Özel, F.; Freire, P. Masses, radii, and the equation of state of neutron stars. An. Rev. Astr. Astroph. 2016, 54, 401. [Google Scholar] [CrossRef] [Green Version]
  41. Alam, N.; Agrawal, B.K.; Fortin, M.; Pais, H.; Providência, C.; Raduta, A.R.; Sulaksono, A. Strong correlations of neutron star radii with the slopes of nuclear matter incompressibility and symmetry energy at saturation. Phys. Rev. C 2016, 94, 052801. [Google Scholar] [CrossRef]
  42. Li, B.-A.; Magno, M. Curvature-slope correlation of nuclear symmetry energy and its imprints on the crust-core transition, radius and tidal deformability of canonical neutron stars. Phys. Rev. C 2020, 102, 045807. [Google Scholar] [CrossRef]
  43. Li, B.-A.; Han, X. Constraining the neutron-proton effective mass splitting using empirical constraints on the density dependence of nuclear symmetry energy around normal density. Phys. Rev. B 2013, 727, 276. [Google Scholar] [CrossRef] [Green Version]
  44. Oertel, M.; Hempel, M.; Klähn, T.; Typel, S. Equations of state for supernovae and compact stars. Rev. Mod. Phys. 2017, 89, 015007. [Google Scholar] [CrossRef]
  45. Lalazissis, G.A.; König, J.; Ring, P. New parametrization for the Lagrangian density of relativistic mean field theory. Phys. Rev. C 1997, 55, 540. [Google Scholar] [CrossRef] [Green Version]
  46. Carriere, J.; Horowitz, C.J.; Piekarewicz, J. Low-mass neutron stars and the equation of state of dense matter. Astrophys. J. 2003, 593, 463. [Google Scholar] [CrossRef] [Green Version]
  47. Pais, H.; Providência, C. Vlasov formalism for extended relativistic mean field models: The crust-core transition and the stellar matter equation of state. Phys. Rev. C 2016, 94, 015808. [Google Scholar] [CrossRef]
  48. Kaminker, A.D.; Haensel, P.; Yakovlev, D.G. Nucleon superfluidity vs. observations of cooling neutron stars. Astron. Astrophys. 2001, 373, 2. [Google Scholar]
  49. Bardeen, J.; Cooper, L.N.; Schrieffer, J.R. Theory of superconductivity. Phys. Rev. 1957, 108, 1175. [Google Scholar] [CrossRef] [Green Version]
  50. Becker, W. Neutron Stars and Pulsars; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  51. Schwenk, A.; Friman, B.; Brown, G.E. Renormalization group approach to neutron matter: Quasiparticle interactions, superfluid gaps and the equation of state. Nucl. Phys. A 2003, 713, 191. [Google Scholar] [CrossRef] [Green Version]
  52. Takatsuka, T. Energy gap in neutron-star matter. Prog. Theor. Phys. 1972, 48, 1517. [Google Scholar] [CrossRef] [Green Version]
  53. Amundsen, L.; Østgaard, E. Superfluidity of neutron matter: (I). Singlet pairing. Nucl. Phys. A 1985, 437, 487. [Google Scholar] [CrossRef]
  54. Amundsen, L.; Østgaard, E. Superfluidity of neutron matter: (II). Triplet pairing. Nucl. Phys. A 1985, 442, 163. [Google Scholar] [CrossRef]
  55. Chen, J.M.C.; Clark, J.W.; Davé, R.D.; Khodel, V.V. Pairing gaps in nucleonic superfluids. Nucl. Phys. A 1993, 555, 59. [Google Scholar] [CrossRef]
  56. Ho, W.C.G.; Elshamouty, K.G.; Heinke, C.O.; Potekhin, A.Y. Tests of the nuclear equation of state and superfluid and superconducting gaps using the Cassiopeia A neutron star. Phys. Rev. C 2015, 91, 015806. [Google Scholar] [CrossRef]
  57. Yakovlev, D.G.; Urpin, V.A. Thermal and electrical conductivity in white dwarfs and neutron stars. Soviet Astron. 1980, 24, 303. [Google Scholar]
  58. Potekhin, A.Y.; Yakovlev, D.G.; Chabrier, G.; Gnedin, O.Y. Thermal structure and cooling of superfluid neutron stars with accreted magnetized envelopes. Astrophys. J. 2003, 594, 404. [Google Scholar] [CrossRef] [Green Version]
  59. Beznogov, M.V.; Yakovlev, D.G. Statistical theory of thermal evolution of neutron stars. Mon. Not. R. Astron. Soc. 2015, 447, 1598–1609. [Google Scholar] [CrossRef] [Green Version]
  60. Ho, W.C.G.; Wijngaarden, M.J.P.; Chang, P.; Heinke, C.O.; Page, D.; Beznogov, M.; Patnaude, D.J. Cooling of the Cassiopeia A neutron star and the effect of diffusive nuclear burning. AIP Conf. Proc. 2019, 2127, 020007. [Google Scholar]
Figure 1. Mass-radius relation for non-rotating neutron stars (NSs) calculated for set B of the induced surface tension Equation of State (IST EoS) [34]. Horizontal bands correspond to the two most massive NSs, e.g., PSR J0348+0432 [36] (magenta band) and PSR J0740+6620 [37] (blue band). The shaded grey area represents the M-R constraint taken from Refs. [38,39], while the constraint depicted as a cyan area was taken from Ref. [40]. The red line represents the allowed range of NS radius, according to GW170817 event.
Figure 1. Mass-radius relation for non-rotating neutron stars (NSs) calculated for set B of the induced surface tension Equation of State (IST EoS) [34]. Horizontal bands correspond to the two most massive NSs, e.g., PSR J0348+0432 [36] (magenta band) and PSR J0740+6620 [37] (blue band). The shaded grey area represents the M-R constraint taken from Refs. [38,39], while the constraint depicted as a cyan area was taken from Ref. [40]. The red line represents the allowed range of NS radius, according to GW170817 event.
Particles 03 00045 g001
Figure 2. Pressure as a function of baryonic density for symmetric and β -equilibrated nuclear matter for the IST EoS (solid and dashed black curves, correspondingly) and the N L 3 ω ρ EoS. The solid red curve corresponds to the N L 3 ω ρ EoS for symmetric matters (or NL3, for simplicity). The N L 3 ω ρ EoS for the β -equilibrated nuclear matter is shown for L = 55 MeV (dashed green curve), L = 88 MeV (dashed blue curve) and L=118 MeV (dashed magenta curve). The shaded gray area corresponds to the flow constraint for symmetric matter [10].
Figure 2. Pressure as a function of baryonic density for symmetric and β -equilibrated nuclear matter for the IST EoS (solid and dashed black curves, correspondingly) and the N L 3 ω ρ EoS. The solid red curve corresponds to the N L 3 ω ρ EoS for symmetric matters (or NL3, for simplicity). The N L 3 ω ρ EoS for the β -equilibrated nuclear matter is shown for L = 55 MeV (dashed green curve), L = 88 MeV (dashed blue curve) and L=118 MeV (dashed magenta curve). The shaded gray area corresponds to the flow constraint for symmetric matter [10].
Particles 03 00045 g002
Figure 3. Density dependence of the critical temperature for the considered singlet and triplet neutron gaps (left panel) and proton singlet pairing gaps (right panel).
Figure 3. Density dependence of the critical temperature for the considered singlet and triplet neutron gaps (left panel) and proton singlet pairing gaps (right panel).
Particles 03 00045 g003
Figure 4. Cooling curves for stars of different mass M / M = 1.2 , 1.4 , 1.6 , 1.8 , 2.0 for the case of unpaired matter. T S denotes the surface temperature at infinity. The solid curves correspond to the light-element envelope ( η = 10 7 ), while the dashed curves were obtained for the heavy-element envelope. The data points are taken from [59].
Figure 4. Cooling curves for stars of different mass M / M = 1.2 , 1.4 , 1.6 , 1.8 , 2.0 for the case of unpaired matter. T S denotes the surface temperature at infinity. The solid curves correspond to the light-element envelope ( η = 10 7 ), while the dashed curves were obtained for the heavy-element envelope. The data points are taken from [59].
Particles 03 00045 g004
Figure 5. The same as Figure 4, but considering the effect of neutron superfluidity in the 1 S 0 channel via the SFB model [51] and proton superconductivity in the 1 S 0 channel with the AO model [53] (left panel) and CCDK model [55] (right panel).
Figure 5. The same as Figure 4, but considering the effect of neutron superfluidity in the 1 S 0 channel via the SFB model [51] and proton superconductivity in the 1 S 0 channel with the AO model [53] (left panel) and CCDK model [55] (right panel).
Particles 03 00045 g005
Figure 6. The same as lower panel of Figure 4, but supplemented with the triplet neutron pairing in the core of the star described by the T72 model [52] (left panel) and the AO model [54] (right panel).
Figure 6. The same as lower panel of Figure 4, but supplemented with the triplet neutron pairing in the core of the star described by the T72 model [52] (left panel) and the AO model [54] (right panel).
Particles 03 00045 g006
Figure 7. Temperature decrease rate for the NS in Cas A. The brown diamonds ( N H of 2006 [20]) and red squares ( N H of 2015 [60]) correspond to Chandra ACIS subarray mode data. Black circles indicate temperature measured using Chandra ACIS-S Graded spectra [21]. Fit results are shown in blue, green, purple, and red.
Figure 7. Temperature decrease rate for the NS in Cas A. The brown diamonds ( N H of 2006 [20]) and red squares ( N H of 2015 [60]) correspond to Chandra ACIS subarray mode data. Black circles indicate temperature measured using Chandra ACIS-S Graded spectra [21]. Fit results are shown in blue, green, purple, and red.
Particles 03 00045 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsiopelas, S.; Sagun, V. Neutron Star Cooling Within the Equation of State With Induced Surface Tension. Particles 2020, 3, 693-705. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3040045

AMA Style

Tsiopelas S, Sagun V. Neutron Star Cooling Within the Equation of State With Induced Surface Tension. Particles. 2020; 3(4):693-705. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3040045

Chicago/Turabian Style

Tsiopelas, Stefanos, and Violetta Sagun. 2020. "Neutron Star Cooling Within the Equation of State With Induced Surface Tension" Particles 3, no. 4: 693-705. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3040045

Article Metrics

Back to TopTop