Next Article in Journal
A Light Mediator Relating Neutrino Reactions
Previous Article in Journal
Determination of the Magnetic Field Strength and Geometry in the Accretion Disks of AGNs by Optical Spectropolarimetry
Previous Article in Special Issue
Heating in Magnetar Crusts from Electron Captures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transport Coefficients of Hyperonic Neutron Star Cores

1
Ioffe Institute, 26 Politekhnicheskaya St., 194021 St. Petersburg, Russia
2
INFN Sezione di Catania, Via Santa Sofia 64, 95123 Catania, Italy
*
Author to whom correspondence should be addressed.
Universe 2021, 7(6), 203; https://doi.org/10.3390/universe7060203
Submission received: 22 April 2021 / Revised: 11 June 2021 / Accepted: 15 June 2021 / Published: 20 June 2021
(This article belongs to the Special Issue Fundamental Processes in Neutron Stars and Supernovae)

Abstract

:
We consider transport properties of the hypernuclear matter in neutron star cores. In particular, we calculate the thermal conductivity, the shear viscosity, and the momentum transfer rates for np Σ Λ e μ composition of dense matter in β –equilibrium for baryon number densities in the range 0.1–1 fm 3 . The calculations are based on baryon interactions treated within the framework of the non-relativistic Brueckner-Hartree-Fock theory. Bare nucleon-nucleon (NN) interactions are described by the Argonne v18 phenomenological potential supplemented with the Urbana IX three-nucleon force. Nucleon-hyperon (NY) and hyperon-hyperon (YY) interactions are based on the NSC97e and NSC97a models of the Nijmegen group. We find that the baryon contribution to transport coefficients is dominated by the neutron one as in the case of neutron star cores containing only nucleons. In particular, we find that neutrons dominate the total thermal conductivity over the whole range of densities explored and that, due to the onset of Σ which leads to the deleptonization of the neutron star core, they dominate also the shear viscosity in the high density region, in contrast with the pure nucleonic case where the lepton contribution is always the dominant one.

1. Introduction

The interior composition of neutron stars (NSs)—among the most dense material objects in the Universe—is not known yet and various phases of dense matter might be expected in their interiors, see, e.g., Ref. [1], for a recent review. With masses of the order of the Solar mass and ∼10 km radii, the density of matter in the innermost regions of NSs, their cores, can reach values a few times that of symmetric nuclear matter at saturation, n 0 = 0.16 fm 3 . In the simplest model, the matter of NS cores consists of uniform fluid of neutrons (n), protons (p), and leptons (electrons (e) and muons ( μ )) in equilibrium with respect to the weak interaction ( β -stable matter). However, at densities of about (2–3) n 0 simple energy arguments suggest that other baryonic degrees of freedom, such as, for instance, hyperons, can appear. This possibility was first proposed by Ambartsumyan and Saakyan [2] in 1960 and it has been later extensively studied in great detail by many authors using different phenomenological [3,4,5,6,7,8,9,10,11,12] or microscopical [13,14,15,16,17,18,19,20,21,22,23,24,25,26] approaches to the equation of state (EOS) of NS matter with hyperons. However, although the presence of hyperons in NSs seems to be energetically unavoidable, the softening that their presence induces on the EOS leads to NS maximum masses which are incompatible with the observation of the unusually high masses of the millisecond pulsars PSR J1903 + 0327 ( 1.667± 0.021 M ) [27], PSR J1614−2230 ( 1.928± 0.017 M ) [28], PSR J0348 + 0432 ( 2.01± 0.04 M ) [29] and the most recent one PSR J0740+6620 ( 2.14 0.09 + 0.10 M ) [30]. The solution of this problem, commonly known in the literature as the ‘hyperon puzzle’, is not easy and it has become a subject of very active research in recent years, see, e.g., Refs. [31,32] and references therein for a recent and comprehensive review. An alternative way to circumvent the hyperon puzzle is to invoke the appearance of other hadronic degrees of freedom that push the onset of hyperons to higher densities, such as for instance meson condensates. The possible existence of a Bose–Einstein condensate of negative kaons in the inner core of NSs has been also extensively considered in the literature (see e.g., Refs. [33,34,35,36,37] and references therein). When the K chemical potential, μ K , becomes smaller than the one of the electron, the process e K + ν e becomes energetically possible and kaons condensate in the NS interior. The onset density for this process was calculated to be in the range (2.5–5) n 0 [35,36]. However, the appearance of kaon condensation induces also a strong softening of the EOS and consequently leads to the reduction of the NS maximum mass, as in the case of hyperons, to values below the current observational limits. The role of kaons in NSs is out of the scope of the present paper and the interested reader is referred to the original works on this topic [33,34,35,36,37] for a comprehensive description of the implications of kaon condensation on the structure and evolution of NSs, and to Ref. [38] for a recent review of the present status of experimental and theoretical developments on kaon physics.
Neutron stars are evolving entities, where various dynamical processes are believed to occur. One can mention among others, for instance, NS cooling, NS oscillations due to undamped instabilities in rotating stars, or magnetic field dissipation. The theoretical modelling of these processes requires knowledge of the transport properties of dense NS matter [39]. In recent studies, the transport properties of ultradense nucleon or pure quark matter have received considerable attention, see e.g., Ref. [39] for a review. Conversely, the coverage of the strange sector has been much lower. Therefore, in this paper we make a step in the direction of filling this gap and consider transport coefficients of the hypernuclear NS core matter. Specifically, we focus on the thermal conductivity and shear viscosity as well as on the momentum transfer rates in binary collisions. The latter coefficients describe friction between the components of the matter in the diffusion process. If the diffusion is caused by the external electromagnetic field, these coefficients enter the generalized Ohm law and determine the electrical conductivity, which is important for studies of the magnetic field evolution of NSs, see, e.g., Refs. [39,40,41,42]. Our results are based on the classical formalism of the kinetic theory of Fermi systems [43] adapted for the NS context, see Refs. [39,44,45,46,47], and references therein.
The basic problem in the transport theory of strongly interacting systems lies in the poor knowledge of many-body interactions at high densities. Various methods are employed, often having poor control of the underlying assumptions and approximations in the high density region. Even within the same theoretical framework, differences in the microphysical input may lead to drastically different results. For instance, for nuclear matter it was found that various nuclear potential models, all behaving well at low densities, result in an order of magnitude difference of the transport properties at several nuclear densities [47]. In the case of hypernuclear matter the situation is even more tricky, because the nucleon-hyperon (NY) and hyperon-hyperon (YY) interactions are, unfortunately, still poorly constrained. Contrary to the nucleon-nucleon (NN) interaction, which is fairly well known due to the larger number of existing scattering data and measured properties of nuclei, experimental data in the hyperonic sector (see Refs. [31,32,38] and references therein) is not sufficient yet to constrain these interactions with a comparable precision.
In the present work we calculate the transport coefficients of np Σ Λ e μ matter. The number fractions of different matter constituents are found from the requirements of electric charge neutrality and β -stability. The EOS of the hyperonic matter is obtained within the non-relativistic Brueckner-Hartree-Fock (BHF) approach using realistic NN, NY, and YY interactions, see Ref. [48] for a general overview of this approach and Refs. [13,14,15,16,17,18] for details of its application to hypernuclear matter. As a first step, in this work, we use a couple of interaction models for the NY and YY interactions in order to explore the difference with the results for pure nucleonic compositions; a more systematic study of the dependence of the results on the choice of baryon interaction is left for future works. Within the nucleonic sector we use, in particular, the Argonne v18 (Av18 for short) potential [49] supplemented with an effective density-dependent two-body force derived from the three-nucleon force of the Urbana IX (UIX for short) model [50]. For the interactions involving hyperons we employ two versions of the Nijmegen Soft-Core 97 (NSC97 in the next) meson-exchange model [51,52], namely, the models NSC97e and NSC97a. Three-body interactions involving hyperons are, unfortunately, still quite uncertain and, therefore, are not included in this work. In the following, since the nucleonic part is described always with the Av18 NN potential plus the UIX NNN force, for simplicity we will use the labels NSC97e and NSC97a when referring to the two sets of models used.
The paper is organized as follows. In Section 2 we briefly review the basic expressions for the transport coefficients of Fermi systems. In Section 3 we shortly present the BHF approach of hypernuclear matter. We present the transport coefficients in Section 4, and conclude in Section 5. In what follows, we set = k B = c = 1 .

2. Transport Coefficients

Let us briefly review the formalism used to calculate the transport coefficients in the case of neutron star cores with hyperons. A detailed description of this formalism for the nuclear matter case can be found in Refs. [46,47,53].
The basic expressions for calculating the thermal conductivity κ c and the shear viscosity η c of particle species c, and the momentum transfer rates J c i in the binary collisions between particle species c and i are
κ c = π 2 3 T n c p F c λ c κ ,
η c = 1 5 n c p F c λ c η ,
J c i = n c p F c λ c i D 1 ,
where T is the temperature, n c and p F c are, respectively, the particle number densities and Fermi momenta, λ c α with α = κ , η , and λ c i D are the effective mean free paths which are, in general, different for the different transport coefficient in question. Having this in mind, below we drop the index α for brevity and reinstall it only when necessary. The effective mean free paths are obtained in the kinetic theory from the solution of the system of linearized transport equations for the distribution functions of the quasiparticles.
The simplest approximation, which however was proved to be very good in degenerate matter, is the lowest-order variational approximation, in which case the system of transport equations reduces to the system of linear equations for the effective mean free paths
i Λ c i λ i = 1 .
The matrix of this system is the transport matrix Λ c i , which is connected to the effective transport cross-sections σ c i , σ c i of binary collisions as
Λ c c = i n i σ c i + n c σ c c ,
Λ c i = n i σ c i , i c .
Effective transport cross-sections σ c i and σ c i in Equations (5) and (6) are obtained by averaging of a squared matrix element of the scattering matrix with certain angular factors. In our case of np Σ Λ e μ matter, particles interact via the electromagnetic and strong forces. Leptons participate only in the electromagnetic interactions, while the strong interaction is the dominant channel for the baryons (see, however, Section 4.1). The explicit integral expressions for the transport cross-sections convenient to use for the strong interactions can be found in Ref. [47] and the form of these quantities more convenient for electromagnetic interactions is given, e.g., in Ref. [54].
The momentum transfer rates in Equation (3) are usually treated in the lowest-order moment expansion of the kinetic theory which is practically equivalent to the variational method described above. The effective partial mean free paths for binary collisions are
λ c i D 1 = n i σ c i D ,
where σ c i D is the transport cross-section for the diffusion problem [47].
Notice, that the transport matrices in Equations (4)–(6) are not symmetric. Non-diagonal matrix elements of the transport matrices obey the relations n c p F c 1 Λ c i κ = n i p F i 1 Λ i c κ and n c p F c Λ c i η = n i p F i Λ i c η [47]. The momentum transfer rates in Equation (3) are symmetric, J c i = J i c .
In the quasiparticle approximation, only the collisions in the vicinity of the Fermi surface contribute to transport. If the matrix element does not depend on the energy transfer in collisions (which is ∼T), one finds the simple scaling rules valid for Fermi systems at low temperatures, σ c i m c 2 m i 2 T 2 , where m c and m i are the quasiparticles effective masses on the Fermi surface [43]. This is the case for collisions mediated by the nuclear forces. Transport cross-sections governed by the electromagnetic interactions in NS cores obey non-Fermi liquid temperature dependence as a result of the dynamical character of the plasma screening in the dominant transverse channel of the electromagnetic interaction [55,56]. The matrix elements of the ‘electromagnetic’ transport matrix do not follow ∝ T 2 scaling, but, generally, have weaker temperature dependence; the exact scaling is subject to the transport problem in question [39].

3. Brueckner–Hartree–Fock Approach of Hypernuclear Matter

Calculations of the transport cross-sections, and consequently of the transport coefficients, require as an input the squares of the different scattering transition matrix elements and the effective masses of the involved baryons which should be provided by a nuclear many-body theory. Various methods have been considered to solve the nuclear many-body problem: the variational approach [57], the correlated basis function (CBF) formalism [58], the self-consistent Green’s function (SCGF) technique [59,60], or the Brueckner–Bethe–Goldstone (BBG) [61] and the Dirac–Bruecker–Hartree–Fock (DBHF) theories [62,63,64]. Nevertheless, although all of them have been extensively applied to the study of nuclear matter, up to our knowledge, only the BBG theory in the BHF approximation [13,14,15,16,17,18,19,22,25,26], and very recently the DBHF theory [20], the V-low-k approach [21], and the quantum Monte Carlo method [23,24] have been extended to the hyperonic sector.
In this section we briefly review the non-relativistic Brueckner–Hartree–Fock approximation of hypernuclear matter. We would like to note that although this approach is non-relativistic we found that the speed of sound is, in all the range of densities considered here, always smaller than the speed of light. The reason is that the softening of the EOS induced by the appeareance of the hyperons guarantees causality. The BHF approximation starts with the construction of all baryon-baryon (NN, YN, and YY) G-matrices which describe the interaction between two baryons in the presence of a surrounding medium. The G-matrices are obtained by solving the well-known coupled-channel Bethe–Goldstone integral equation
k B 1 k B 2 | G ( ω ) | k B 3 k B 4 = k B 1 k B 2 | V | k B 3 k B 4 a a a a a a a a + B i B j k B 1 k B 2 | V | k B i k B j k B i k B j | Q | k B i k B j k B i k B j | G ( ω ) | k B 3 k B 4 ω E B i ( k B i ) E B j ( k B j ) + i η ,
where V is the bare baryon-baryon interaction, Q is the Pauli operator that prevents the intermediate baryons B i and B j from being scattered to states below their respective Fermi momenta, ω is the sum of the non-relativistic single-particle energies of the interacting baryons, and η is an infinitesimal positive quantity.
The single-particle energy of a baryon B i is given by
E B i ( k B i ) = m B i + k B i 2 2 m B i + Re [ U B i ( k B i ) ] ,
where m B i denotes the rest mass of the baryon, and the (complex) single-particle potential U B i ( k B i ) represents the average field ‘felt’ by the baryon owing to its interaction with other baryons. In the BHF approximation, U B i ( k B i ) is calculated through the ‘on-shell’ G-matrix, and is given by
U B i ( k B i ) = B j k B j n B j ( k B j ) k B i k B j G ω = E B i ( k B i ) + E B j ( k B j ) k B i k B j A .
Here n B j ( k B j ) is the occupation number of the baryon species B j , and the index A indicates that the matrix elements are properly antisymmetrized when baryons B i and B j belong to the same isomultiplet. We note here that the so-called continuous prescription has been adopted for the single-particle potentials when solving the Bethe–Goldstone equation, since, as it was shown in Refs. [65,66], the contribution to the energy per particle from three-hole line diagrams is minimized in this prescription. As we already said in the introduction all calculations have been carried out with the Av18 NN potential [49] supplemented with the UIX three-nucleon force [50], which, for the use in the BHF calculations, was reduced to a two-body density-dependent force by averaging over the spatial, spin, and isospin coordinates of the third nucleon in the medium [67,68,69]. This three-nucleon force contains two parameters that are fixed by requiring that the BHF calculation reproduces the energy and saturation density of the symmetric nuclear matter. For the NY and YY interactions, we employed the NSC97e and NSC97a models [51,52]. The reason behind the choice of these models is that both models result in the best predictions for the hypernuclear observables among the potentials constructed by the Nijmegen group. Three-body forces involving hyperons, i.e., forces of the type NNY, NYY, and YYY, have been ignored in the calculation due to the large uncertainties still existing about these kind of forces.
To solve Equations (8)–(10), one starts from a reasonable guess for the single-particle potentials and obtain initial value of the G-matrix from Equation (8), which is solved in a partial wave basis including contributions up to a total angular momentum J = 4 . Then the new single-particle potential is computed from Equation (10) and passed as an input to Equations (8)–(9). These iterations continue until the desired level of convergence is reached. Once a self-consistent solution of Equations (8)–(10) is obtained, the baryon contribution to the total energy density can be calculated simply as
ε B = n B B i k B i n ( k B i ) m B i + k B i 2 2 m B i + 1 2 Re [ U B i ( k B i ) ] ,
where n B is the baryon number density. Adding to ε B the contribution from the noninteracting leptons, ε L , the composition and the EOS of neutron star matter can then be obtained from the requirement of equilibrium under weak interaction processes, μ i = b i μ n q i μ e ( b i and q i denoting the baryon number and electric charge of the species i, respectively), and electric charge neutrality, i q i n i = 0 . The chemical potentials of various species and the pressure are computed from the usual thermodynamic relations, μ i = ε / n i and P = n B ε / n B ε , where ε = ε B + ε L is the total energy density. The EOS (pressure versus total energy density) and the particle fractions x i = n i / n B of the different species are shown respectively in panels (a) and (b) of Figure 1. Vertical dashed lines in Figure 1b show the threshold densities for appearance of muons ( 0.132 fm 3 ), Σ -hyperons ( 0.248 fm 3 for NSC97e and 0.278 fm 3 for NSC97a), and Λ -hyperons ( 0.584 fm 3 for NSC97e and 0.642 fm 3 for NSC97a).
The G-matrices needed in the calculation of the transport-cross sections are taken on-shell and on the corresponding Fermi surfaces of the involved baryons, i.e., the starting energy is taken as ω = E B i ( k F B i ) + E B j ( k F B j ) . The effective masses of the different baryons, required in the calculation of the transport coefficients, are calculated from their corresponding single-particle energies as
m B i ( k B i ) m B i = k B i m B i d E B i ( k B i ) d k B i 1 .
Figure 2 shows the effective masses of the different baryons evaluated at each density at their corresponding Fermi momenta. The calculation of the effective masses via Equation (12) requires numerical differentiation of the single-particle potentials resulting from the solution of the BBG equations, which is a delicate numerical task. The results frequently contain considerable numerical noise. It is customary to approximate the results with some smooth function of the baryon density [70]. Here we employ this approach and fit m c ( n B ) functions with fourth-order polynomials. The fit results are given in Appendix A.

4. Results and Discussion

4.1. Transport Matrices and Mean Free Paths

We have found that the results for the NSC97e and NSC97a models are qualitatively similar. Therefore for brevity we mainly focus on the NSC97e model until indicated otherwise.
The transport matrix Λ c i in Equation (4) contains a diagonal part [Equation (5)] and non-diagonal elements [Equation (6)]. If the non-diagonal elements can be neglected (as it is frequently the case, see below), the equations for the different species decouple and the effective mean free paths are given by a simple formula λ c = Λ c c 1 . The diagonal elements of the transport matrix thus represent useful approximations for the inverse mean free paths. In Figure 3 we show the diagonal components of the transport matrix for the baryon sector for thermal conductivity [panel (a)] and shear viscosity cases [panel (b)]. In this figure only the strong interaction channels are included and temperature-independent combinations Λ c c T 8 2 are shown, with T 8 = T / ( 10 8 K ) . The behavior of Λ c c κ and Λ c c η is qualitatively similar. The lowest values of Λ c c are obtained for neutrons. Thus, the neutrons are expected to have the largest mean free paths and dominate the baryon transport. The protons and Σ hyperons scatter one-two orders of magnitude more effectively, and the values of Λ c c for the Λ hyperon lie in between. Notice almost identical results for protons and Σ hyperons in Figure 3 at n B 0.5 fm 3 .
These results can be understood by analyzing the partial contributions from the scattering on different species in Equation (5). In the matrix-shaped Figure 4 and Figure 5 with the solid lines we plot the relative contributions to Λ c c for thermal conductivity and shear viscosity cases, respectively. The rows in these figures are numerated with index c and columns with index i of Equation (5). Namely, we show in each panel with solid lines the ratio Υ ˜ c i = ( n i σ c i + δ c i n c σ c c ) Λ c c 1 (notice that the collisions of the like species contain additional n c σ c c contribution). It is clear that soon after the Σ appearance, their mutual collisions with protons become completely dominant scattering mechanisms of these species. Taking into account that at n B 0.5 fm 3 the proton and Σ number fractions are similar due to the requirement of charge neutrality in lack of charged leptons, see Figure 1b, the coincidence of their transport matrix elements in Figure 3 is natural. Because of the large values of the transport matrix elements, the proton and Σ mean free paths, λ p and λ Σ , are small and, therefore, these two baryons are not expected to influence the transport equations for other species. Having this in mind and taking into account the relatively small fraction of these particles (in comparison to neutrons) in β -stable matter, one expects that they do not give sizable contribution to overall transport coefficients.
To analyze the mutual impact of the mean free paths of different species, it is instructive to introduce the renormalized mean free paths defined as λ ˜ c = λ c Λ c c . Then the system (4) takes the form
1 = i Λ c i Λ i i λ ˜ i i Λ ˜ c i λ ˜ i ,
where the renormalized transport matrix elements Λ ˜ c i = Λ c i / Λ i i are introduced.
With red dashed (for positive values) and blue dash-dotted (for negative values) lines in Figure 4 and Figure 5 we show logarithms of the non-diagonal relative matrix elements log 10 | Λ ˜ c i | , c i . In the majority of cases, | Λ ˜ c i | 1 , which means that the influence of λ i on the equations for the c species in that cases is negligible. Please note that irrespectively of the values of | Λ ˜ c i | where one of the indices refers to protons or Σ hyperons, the corresponding terms can be safely neglected since λ p and λ Σ are small as discussed above. This means that protons and Σ hyperons can be treated as passive scatterers. The remaining Λ ˜ c i matrix elements correspond to the n Λ subsystem. One observes that Λ ˜ Λ n η and Λ ˜ n Λ η and also, for the thermal conductivity case and at the largest densities, Λ ˜ Λ n κ and Λ ˜ n Λ κ , are significant and may need to be taken into account.
Up to now, we discussed the collisions in the baryon subsystem mediated by the strong interaction. The charged particles, p Σ e μ in our case, also participate in the electromagnetic interactions1. We neglect the interference between the electromagnetic and strong amplitudes and consider these interaction channels separately. The important difference between the mean free paths governed by electromagnetic and strong interactions in NS cores is in the non-Fermi-liquid temperature dependence of the former. This stems from the dynamical character of the plasma screening in the dominant transverse channel of the electromagnetic interaction [55,56]. As already stated, the matrix elements of the ‘electromagnetic’ transport matrix do not follow the same T 2 scaling as the matrix elements for the strong interaction. For shear viscosity the contribution of the electromagnetic interactions of baryons is several orders of magnitude smaller than the contribution of strong interactions and, therefore, it can be always neglected, see e.g., Ref. [47]. For the thermal conductivity case at low temperatures, the diagonal elements of the transport matrix of charged particles are found to be given by a simple universal expression
Λ c c κ , e m = 6 ζ ( 3 ) π 2 α f T = 2.45× 10 6 T 8 cm 1 ,
irrespectively of the content of the charged particles in the matter. For protons and Σ hyperons at sufficiently low temperatures (e.g., T 10 7 K) it can become comparable with Λ c c κ resulting from strong interactions. This makes the mean free paths of p and Σ even smaller so that they can be obtained by treating the baryon subsystem separately. Furthermore, non-diagonal matrix elements of Λ c i e m are small both for viscosity and thermal conductivity [47]. As a consequence, lepton and baryon subsystems can indeed be treated separately. Moreover, in both subsystems p and Σ act as passive scatterers. The influence of their non-equilibrium distributions on transport properties of other constituents is negligible. The combined contribution of p and Σ to the total values of transport coefficients is found to be less than 10% for thermal conductivity and less then 5% for shear viscosity. In what follows, the transport coefficients of leptons are calculated according to Refs. [71,72].

4.2. Corrections to the Variational Solution

The variational solution gives a first approximation to the transport coefficients. In a multicomponent Fermi-liquid it is possible to obtain the exact solutions of the corresponding system of transport equations [44,45]. However, because of the non-Fermi-liquid behavior of the matrix elements of the electromagnetic interaction, obtaining the exact solution for the full n p Σ Λ e μ mixture is complicated and standard approaches [44,45] do not work. Nevertheless, as discussed above, it is sufficient to consider the n Λ and e μ subsystems separately, with protons and Σ particles acting as the passive scatterers. The corrections to the variational solution in the electromagnetic sector were considered in Refs. [71,72] and were found to be less than 5–10%.
The corrections to the variational solution in the n Λ subsystem can be calculated using the method described in Ref. [45], or by the numerical solution of the system of integral equations [53]. It is customary to express transport coefficients resulting from the exact solution, via the correction coefficients
κ ex = C ex κ κ var , η ex = C ex η η var .
In Figure 6 we show the correction factors C ex for the thermal conductivity (solid line) and shear viscosity (dashed line) for the n Λ subsystem, where protons and Λ -hyperons are treated as passive scatterers. The results are similar to the nuclear matter [46,47,53]. For the shear viscosity, the correction coefficient does not exceed 5%; this correction is truly minor having in mind other systematic uncertainties and it can be neglected. The correction coefficient for the thermal conductivity case is about 20% for the n Λ and can be included in calculations.
We do not discuss the corrections to the simplest approximation for the momentum transfer rates J c i . These quantities are required as an input to the studies of the magnetic field evolution in NS cores. In this case the exact solution of the transport equations in the presence of the magnetic field is required. Such solutions have never been constructed for multicomponent Fermi-systems to the best of our knowledge and require a separate study. Since similar transport cross-sections (with the same leading-order dependencies) enter expressions for the shear viscosity and the momentum transfer rates, one can expect that the corrections resulting from better approximations to the transport equations will not be large.

4.3. Transport Coefficients

Figure 7 and Figure 8 show, respectively, the partial contributions of the matter constituents to the thermal conductivity and shear viscosity coefficients. These contributions are obtained via the simplest variational solution. The results are shown for the temperature-independent (in a normal Fermi-liquid) combinations κ T 8 and η T 8 2 . Due to the non-Fermi liquid behavior of the lepton transport coefficients, we show results for three temperature values, T = 10 7 , 10 8 , and 10 9 K. Thick solid lines in Figure 7 and Figure 8 show total thermal conductivity and shear viscosity coefficients, respectively, including the corrections described in Section 4.2.
Among the baryon sector, neutrons always dominate the transport, such as in the case of nuclear matter [46,47]. In the case of the thermal transport, neutrons dominate also the total value of the thermal conductivity except for T = 10 9 K where the lepton contribution becomes comparable to the nucleon one, see Figure 7c. This behavior is similar to that found in the nucleonic NS cores [46,47].
A different situation is observed in the case of the shear viscosity, see Figure 8. Here the leptonic contribution is dominant at lower densities (the upper boundary of this region depends on temperature), while neutrons give the main contribution at higher densities due to decrease in the lepton number fractions induced by the onset of the Σ hyperon. This is in contrast with the results of the nucleonic matter where the lepton contribution was found to be always dominant in the non-superfluid matter [46,47]. The total baryon contribution to η is similar for both nucleon and hyperon NS core compositions. However, as a result of the suppression of lepton contribution, the total shear viscosity of the hypernuclear matter can be several orders of magnitude smaller than those in the nuclear matter.
Let us investigate the model dependence of the obtained results. Comparing the results for the NSC9e and NSC97a potentials we find that although the partial transport coefficients exhibit some variations, the total values dominated by neutrons and leptons do not change much. This is illustrated in Figure 9, where we show thermal conductivity (panel (a)) and shear viscosity (panel (b)) for the NSC97e (solid lines) and NSC97a (dashed lines) for three values of temperature as indicated in the plot. Notice that here the absolute values of κ and η (without temperature rescaling) are plotted for convenience. The difference between the two models is not large. The reason is that the NSC97e and NSC97a models for the NY and YY interactions are quite similar and differ only on the values of a few meson-baryon coupling constants, cut-off values of the vertex form factors, and some other parameters. Both describe, as mentioned before, with the same accuracy the existing hypernuclear data. In some sense this resembles the results for the nuclear matter of Ref. [47], where it was found that the main source for difference between the various models originates from the model for the three-body interactions which we do not include here for hyperons and do not alter for nucleons. Neutrons, as the most significant particles in terms of transport properties, are not affected by the change of model. The situation is different for other particles, see below.
Finally, in Figure 10 we show the momentum transfer rates in the baryon collisions. Panel (a) corresponds to the NSC97e model and panel (b) to the NSC97a model. Here, we plot the temperature-independent combination J c i T 8 2 . The contribution of the electromagnetic interactions to J c i is negligible like in the case of the shear viscosity above. The momentum transfer rates mediated by the electromagnetic interactions can be calculated following the expressions in Refs. [42,73]. They are several orders of magnitude smaller than J c i mediated by the strong interaction, and we do not show them here. The momentum transfer rate between the nucleons, J n p , is in an order-of-magnitude agreement with the results of the pure nucleonic matter [47]. A large value of J p Σ draws attention in Figure 10. This is due to the particularly strong interaction in this channel predicted by the NSC97 models, see Refs. [51,52]. Notice that for this transport coefficient the difference between the NSC97e and NSC97a model becomes apparent.
To facilitate the application of our study, we give the numerical results for the total thermal conductivity and shear viscosity coefficients, as well as the momentum transfer rates in Table A3 and Table A4 in Appendix B.

5. Conclusions

In this work, we calculated the transport coefficients of the hypernuclear matter in the NS cores. In particular, we calculated the thermal conductivity, the shear viscosity, and the momentum transfer rates for np Σ Λ e μ matter in β -equilibrium for baryon number densities between 0.1 and 1 fm 3 . The hypernuclear matter EOS was described within the non-relativistic BHF approach using realistic NN, NY, and YY interactions. Our results have shown that the proton, Σ , and Λ do not give sizable contributions to the overall transport coefficients, being the total baryon contribution dominated by the neutron one as in the case of neutron star cores with only nucleons. However, these species are important in mediating the neutron mean free path. In particular, we found that neutrons dominate also the total thermal conductivity over the whole range of densities explored and that, due to the onset of Σ which leads to the deleptonization of the neutron star core, they dominate also the shear viscosity on the high density region, in contrast with the pure nucleonic case where the lepton contribution is always the dominant one.
Our paper makes only a first step in the studies of the transport properties of hypernuclear NS cores with realistic baryon interactions. The present study can be improved and extended in various aspects.
In particular, we employed particular potentials both for nucleon and strange sectors and a specific many-body theory. Based on the experience for the nuclear sector [47], we expect that the variations of baryon interaction can lead to a substantial modifications of the values of the transport coefficients at large (e.g., ≳ 2 n 0 ) baryon densities, especially when the three-body interactions are altered. Nevertheless, the qualitative conclusions presented above should be stable over the choice of a particular strong interaction model, which is seen comparing NSC97e and NSC97a results.
Next, the neutron star cores can contain strong magnetic fields. In this case the transport becomes anisotropic due to the Larmor rotation of the charged particles. The inclusion of these effects in not too strong magnetic fields, when the transport cross-sections are not modified, is relatively straightforward. For instance, the momentum transfer rates as calculated here constitute the microphysical input for the magnetic field evolution studies, see, e.g., Ref. [41,42]. Large fields, e.g., ≳ 10 14 G, may affect the polarization tensor of the plasma making the scattering of the charged particles due to electromagnetic interactions anisotropic. This case is more complicated and requires a separate study. In even stronger fields, e.g., ≳ 10 18 G, the particle motion becomes strongly quantized and the description of the transport based on the semiclassical distribution functions becomes invalid; however the presence of such ultra-high magnetic fields can hardly be considered to be realistic.
Finally, due to the presence of the attractive character of some partial wave channels of the baryon interaction, the baryon can be subject to pairing instabilities, see, e.g., Ref. [74], for a review. Naturally, the strength and the character of the particular partial wave interaction is density-dependent. It is expected that the protons in NS cores pair in the 1 S 0 channel, while the neutrons should pair in the 3 P 2 3 F 2 channel. Hyperons may also form superfluids if their interactions are attractive enough. However, a quantitative estimation of the hyperon pairing, contrary to the nucleonic one, has not received so much attention and just a few calculations exist in the literature [75,76,77,78,79,80,81]. In the presence of the paired species, the dynamical processes in NSs are described by a complicated multifluid hydrodynamics of the superfluid/suprerconducting mixtures e.g., [82], and references therein. Baryon pairing affects the transport coefficients of the NS core matter in several ways (see, e.g., [39], for review). We briefly outline the main effects expected in the presence of pairing. In this case the transport of the normal component of matter is governed by the Bogoliubov quasiparticles, which spectra have a gap at the Fermi surface. This results in exponential suppression of the scattering rates of quasiparticles at low temperatures. In NS core context, this effect was considered in Ref. [83] for neutron thermal conductivity and in Refs. [71,72] for lepton transport coefficients in the presence of the proton pairing. Moreover, the collision integral for the Bogouluibov quasiparticles modifies and should include processes with non-conserved number of quasiparticles, e.g., [84]. In addition, if charged species are paired, the electromagnetic polarization of this superconducting matter is modified which affects the scattering rate of all charged quasiparticles (even non-superconducting ones, such as leptons) [39,47]. At low temperatures, the gapless Goldstone modes, e.g., superfluid phonons, can contribute to the transport properties of matter, see, e.g., [85], for a review.
Even more complicated is the case of the magnetized superfluid/superconducting matter, where—if the type II superconductivity is realized—the topological structures called Abrikosov vortices appear. The study of all these effects is much desirable but is beyond the scope of the present paper.

Author Contributions

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

Funding

A part of this work was done during the PHAROS COST STSM #CA16214-44662. The work was partially supported by the Foundation for the Advancement of Theoretical Physics and Mathematics “BASIS”, grant 17-13-205-1.

Data Availability Statement

The data presented in this study are available on request from the corresponding author

Acknowledgments

We thank the anonymous referees for suggestions which helped us to improve the paper. P.S. Thanks the INFN Sez. di Catania for hospitality.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NSNeutron star
EOSEquation of state
BHFBrueckner-Hartree-Fock
Av18Argonne v18
UIXUrbana IX
NSC97Nijmegen Soft-Core 97

Appendix A

We fit the dependence of effective masses on baryon density via the following expression
m c m c = k = 0 4 a c k n B 1 fm 3 k
for c = n , p , Σ , and Λ . The fit coefficients a c k are given in Table A1 for NSC97e model and in Table A2 for NSC97a model. Notice that in the latter case we used third-order polynomials.
Table A1. Coefficients of the polynomial fit (A1) for NSC97e model. Last column gives the density range where the expression (A1) is valid.
Table A1. Coefficients of the polynomial fit (A1) for NSC97e model. Last column gives the density range where the expression (A1) is valid.
Species a 0 a 1 a 2 a 3 a 4 n B Range
n 0.940 0.668 0.678 0.029 0.121 0.10–1 fm 3
p 0.933 2.343 6.504 6.472 2.316 0.10–1 fm 3
Σ 0.853 0.0418 0.590 0.695 0.239 0.25–1 fm 3
Λ 4.256 19.554 28.78 19.27 5.077 0.60–1 fm 3
Table A2. Coefficients of the polynomial fit (A1) for NSC97a model. Last column gives the density range where the expression (A1) is valid. Notice that the a 4 coefficient is absent.
Table A2. Coefficients of the polynomial fit (A1) for NSC97a model. Last column gives the density range where the expression (A1) is valid. Notice that the a 4 coefficient is absent.
Species a 0 a 1 a 2 a 3 n B Range
n 0.934 0.591 0.453 0.029 0.10–1 fm 3
p 0.826 1.092 2.171 1.000 0.10–1 fm 3
Σ 0.905 0.459 0.401 0.078 0.30–1 fm 3
Λ 0.892 3.916 8.742 4.844 0.65–1 fm 3

Appendix B

Table A3 and Table A4 contain the numerical values of the transport coefficients for the hypernuclear NS core composition discussed in the main text for NSC97e and NSC97a models, respectively. We show κ T 8 (columns 2–4) and η T 8 2 (columns 5–7) for three values of the temperature. Total (lepton plus baryon) values of the transport coefficients are given. In the baryon subsystem the corrections for the variational solution for n and Λ are included. We also tabulate the momentum transfer rates J c i T 8 2 in collisions of six possible distinguishable baryon pairs (columns 7–13).
Table A3. Transport coefficients of the hypernuclear matter as a function of baryon density for the NSC97e model. Columns 24: total thermal conductivity for three values of temperature. Columns 57: total shear viscosity for three values of temperature. For η and κ exact solutions of transport equations in the sense described in the text are shown. Columns 814: momentum transfer rates J c i .
Table A3. Transport coefficients of the hypernuclear matter as a function of baryon density for the NSC97e model. Columns 24: total thermal conductivity for three values of temperature. Columns 57: total shear viscosity for three values of temperature. For η and κ exact solutions of transport equations in the sense described in the text are shown. Columns 814: momentum transfer rates J c i .
n B log 10 κ T 8 22 log 10 η T 8 2 18 log 10 J ci T 8 2 29
fm 3 erg cm 1 s 1  K 1 g cm 1 s 1 g cm 3 s 1
10 7  K 10 8  K 10 9  K 10 7  K 10 8  K 10 9  K np n Σ n Λ p Σ p Λ Σ Λ
0.100−0.0320.0380.384−0.201−0.0780.0120.469
0.1250.1750.2300.5370.0510.1920.3020.449
0.1500.3400.3920.6780.240.4010.5320.433
0.1750.4720.5230.8070.3920.5700.7190.423
0.2000.5780.6270.9120.5280.7190.8830.418
0.2250.6570.7070.9950.6510.8541.0320.431
0.2500.6960.7481.0350.7020.8711.0040.455−1.407−0.869
0.2750.7140.7631.0560.6680.8450.9930.510−0.4610.292
0.3000.7500.7951.0710.6390.8170.9740.567−0.2190.611
0.3250.7870.8271.0850.6110.7890.9550.586−0.0730.933
0.3500.7960.8341.0840.5800.7550.9240.6570.0361.067
0.3750.8030.8391.0800.5500.7200.8910.7220.1261.246
0.4000.8060.8411.0740.5210.6840.8550.7870.2101.360
0.4250.8070.8401.0650.4930.6480.8170.8490.2881.466
0.4500.8060.8371.0550.4650.6120.7770.9100.3591.565
0.4750.8030.8321.0430.4370.5750.7360.9700.4251.661
0.5000.7980.8261.0290.4100.5390.6941.0280.4831.749
0.5250.7920.8191.0130.3840.5030.6511.0830.5371.845
0.5500.7850.8100.9960.3590.4680.6081.1380.5861.918
0.5750.7760.8000.9770.3350.4330.5651.1910.6312.006
0.6000.7640.7860.9550.3110.3990.5221.2410.671−1.3982.086−2.312−1.210
0.6250.7520.7730.9330.2870.3660.4801.2870.707−0.7462.172−1.470−0.641
0.6500.7380.7570.9080.2620.3330.4371.3310.741−0.3782.257−0.992−0.328
0.6750.7220.7400.8840.2380.3000.3951.3740.773−0.1182.343−0.655−0.103
0.7000.7070.7240.8640.2130.2670.3521.4170.8040.0782.428−0.4000.071
0.7250.6880.7050.8420.1880.2340.3101.4560.8350.2622.522−0.1690.216
0.7500.6720.6880.8200.1640.2030.2691.4950.8650.4002.6120.0100.334
0.7750.6480.6640.7930.1380.1710.2271.5330.8930.5482.7070.2230.459
0.8000.6320.6470.7700.1150.1410.1881.5700.9200.6552.7970.3570.540
0.8250.6100.6250.7440.0900.1110.1501.6050.9480.7822.9070.5080.612
0.8500.5880.6010.7150.0650.0820.1131.6400.9750.8973.0190.6540.680
0.8750.5670.5790.6870.0410.0540.0781.6751.0011.0003.1370.7760.734
0.9000.5420.5540.6560.0150.0250.0431.7091.0261.1183.2690.8950.774
0.9250.5080.5200.616−0.015−0.0080.0051.7421.0501.2533.4221.0450.823
0.9500.4480.4590.554−0.062−0.057−0.0471.7761.0731.4493.6391.2000.857
0.9750.3870.3980.491−0.109−0.106−0.0991.8091.0961.6453.8551.3550.890
1.0000.2510.2640.368−0.190−0.187−0.1831.8431.1181.9544.1921.4180.892
Table A4. Same as Table A3 but for the NSC97a model.
Table A4. Same as Table A3 but for the NSC97a model.
n B log 10 κ T 8 22 log 10 η T 8 2 18 log 10 J ci T 8 2 29
fm 3 erg cm 1 s 1 K 1 g cm 1 s 1 g cm 3 s 1
10 7 K 10 8 K 10 9 K 10 7 K 10 8 K 10 9 K np n Σ n Λ p Σ p Λ Σ Λ
0.100 0.015 0.0520.391 0.193 0.07 0.0220.446
0.1250.1820.2370.5400.0530.1950.3060.438
0.1500.3380.390.6770.2390.4010.5320.435
0.1750.4650.5170.8040.390.5670.7160.432
0.2000.5680.6190.9070.5250.7160.880.432
0.2250.6480.6980.9910.6480.8511.0280.447
0.2500.7150.7661.0620.7640.9771.1640.467
0.2750.7700.8211.1230.8711.0941.2920.496
0.3000.7590.8131.1220.8030.9961.1580.559 0.513 0.215
0.3250.8060.8541.1430.7720.9711.1460.587 0.235 0.673
0.3500.8120.8581.1430.7350.9371.1210.664 0.077 0.886
0.3750.8160.8601.1400.6990.9001.0900.7330.0421.082
0.4000.8180.8601.1350.6620.8611.0560.7980.1441.208
0.4250.8170.8581.1270.6260.8211.0180.8600.2361.318
0.4500.8140.8541.1180.5890.7800.9780.9180.3181.417
0.4750.8090.8481.1070.5530.7380.9360.9730.3901.509
0.5000.8030.8411.0940.5160.6960.8911.0270.4541.597
0.5250.7950.8321.0800.4800.6520.8451.0790.5101.681
0.5500.7870.8221.0640.4440.6080.7971.1300.5061.764
0.5750.7770.8111.0470.4080.5630.7481.1790.6061.846
0.6000.7650.7981.0280.3720.5180.6971.2270.6481.925
0.6250.7530.7841.0070.3360.4720.6451.2750.6872.003
0.6500.7290.7600.9790.3010.4270.5911.3200.723 1.405 2.079 1.646 1.392
0.6750.7070.7370.9510.2640.3810.5381.3620.754 0.551 2.153 0.797 0.563
0.7000.6850.7140.9220.2270.3350.4841.4020.784 0.161 2.227 0.400 0.202
0.7250.6600.6880.8890.1900.2880.4281.4420.8120.1152.300 0.114 0.045
0.7500.6340.6600.8540.1510.2400.3711.4800.8410.3322.3750.1120.237
0.7750.6050.6300.8160.1110.1910.3121.5180.8680.5202.4470.3050.394
0.8000.5740.5980.7810.0700.1400.2511.5550.8950.6842.5230.4720.526
0.8250.5280.5520.7390.0200.0830.1841.5910.9210.8872.6290.6750.654
0.8500.5050.5290.710 0.017 0.0350.1241.6270.9480.9792.6860.7700.745
0.8750.4690.4920.672 0.062 0.018 0.0591.6620.9741.1102.7730.8880.828
0.9000.4340.4570.633 0.105 0.069 0.004 1.6961.0011.2212.8570.9930.908
0.9250.4100.4320.600 0.142 0.113 0.061 1.7301.0271.3072.9251.0540.949
0.9500.3740.3950.557 0.184 0.162 0.121 1.7621.0531.4163.0191.1350.997
0.9750.3520.3720.523 0.215 0.199 0.169 1.7941.0801.4853.0971.1791.032
1.0000.3110.3300.473 0.257 0.246 0.224 1.8241.1061.6033.2131.2441.060

Note

1
Neutrons and Λ hyperons also couple to the electromagnetic fields due to their magnetic moments; however, this contribution is negligible.

References

  1. Blaschke, D.; Chamel, N. Phases of Dense Matter in Compact Stars. In Astrophysics and Space Science Library; Rezzolla, L., Pizzochero, P., Jones, D.I., Rea, N., Vidaña, I., Eds.; Springer: Cham, Switzerland, 2018; Volume 457, p. 337. [Google Scholar] [CrossRef] [Green Version]
  2. Ambartsumyan, V.A.; Saakyan, G.S. The Degenerate Superdense Gas of Elementary Particles. Sov. Astron. 1960, 4, 187. [Google Scholar]
  3. Balberg, S.; Gal, A. An effective equation of state for dense matter with strangeness. Nucl. Phys. A 1997, 625, 535. [Google Scholar] [CrossRef] [Green Version]
  4. Balberg, S.; Lichtenstadt, I.; Cook, G.B. Roles of hyperons in neutron stars. Astrophys. J. Suppl. Ser. 1999, 121, 515. [Google Scholar] [CrossRef]
  5. Glendenning, N.K. The hyperon composition of neutron stars. Phys. Lett. B 1982, 114, 392. [Google Scholar] [CrossRef] [Green Version]
  6. Glendenning, N.K. Neutron stars are giant hypernuclei? Astrophys. J. 1985, 293, 470. [Google Scholar] [CrossRef]
  7. Glendenning, N.K. Hyperons in neutron stars. Z. Phys. A 1987, 326, 57. [Google Scholar] [CrossRef]
  8. Glendenning, N.K.; Moszkowski, S.A. Reconciliation of neutron-star masses and binding of the Λ hypernuclei. Phys. Rev. Lett. 1991, 67, 2414. [Google Scholar] [CrossRef]
  9. Weber, F.; Weigel, M.K. Baryon Composition and Macroscopic Properties of Neutron Stars. Nucl. Phys. A 1989, 505, 779. [Google Scholar] [CrossRef]
  10. Knorren, R.; Prakash, M.; Ellis, P.J. Strangeness in hadronic stellar matter. Phys. Rev. C 1995, 52, 3470. [Google Scholar] [CrossRef] [Green Version]
  11. Schaffner, J.; Mishustin, I. Hyperon-rich matter in neutron stars. Phys. Rev. C 1996, 53, 1416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Huber, H.; Weber, F.; Weigel, M.K.; Schaab, C. Neutron star properties with relativistic equations of state. Int. J. Mod. Phys. E 1998, 7, 310. [Google Scholar] [CrossRef] [Green Version]
  13. Schulze, H.J.; Baldo, M.; Lombardo, U.; Cugnon, J.; Lejeune, A. Hypernuclear matter in the Brueckner–Hartree–Fock approximation. Phys. Lett. B 1995, 355, 21. [Google Scholar] [CrossRef] [Green Version]
  14. Schulze, H.J.; Baldo, M.; Lombardo, U.; Cugnon, J.; Lejeune, A. Hyperonic nuclear matter in Brueckner theory. Phys. Rev. C 1998, 57, 704. [Google Scholar] [CrossRef] [Green Version]
  15. Baldo, M.; Burgio, G.F.; Schulze, H.J. Onset of hyperon formation in neutron star matter from Brueckner theory. Phys. Rev. C 1998, 58, 3688. [Google Scholar] [CrossRef]
  16. Vidaña, I.; Polls, A.; Ramos, A.; Hjorth–Jensen, M.; Stoks, V.G.J. Strange nuclear matter within Brueckner–Hartree–Fock theory. Phys. Rev. C 2000, 61, 025802. [Google Scholar] [CrossRef] [Green Version]
  17. Baldo, M.; Burgio, G.F.; Schulze, H.J. Hyperon stars in the Brueckner–Bethe–Goldstone theory. Phys. Rev. C 1998, 61, 055801. [Google Scholar] [CrossRef] [Green Version]
  18. Vidaña, I.; Polls, A.; Ramos, A.; Engvik, L.; Hjorth–Jensen, M. Hyperon-hyperon interactions and properties of neutron star matter. Phys. Rev. C 2000, 62, 035801. [Google Scholar] [CrossRef] [Green Version]
  19. Schulze, H.J.; Polls, A.; Ramos, A.; Vidaña, I. Maximum mass of neutron stars. Phys. Rev. C 2006, 73, 058801. [Google Scholar] [CrossRef] [Green Version]
  20. Sammarruca, F. Effect of Λ hyperons on the nuclear equation of state in a Dirac–Brueckner–Hartree–Fock model. Phys. Rev. C 2009, 79, 034301. [Google Scholar] [CrossRef] [Green Version]
  21. Dapo, H.; Schaeffer, B.J.; Wambach, J. Appearance of hyperons in neutron stars. Phys. Rev. C 2010, 81, 035803. [Google Scholar] [CrossRef] [Green Version]
  22. Schulze, H.J.; Rijken, T. Maximum mass of hyperon stars with the Nijmegen ESC08 model. Phys. Rev. C 2011, 84, 035801. [Google Scholar] [CrossRef] [Green Version]
  23. Lonardoni, D.; Pederiva, F.; Gandolfi, S. Accurate determination of the interaction between Λ hyperons and nucleons from auxiliary field diffusion Monte Carlo calculations. Phys. Rev. C 2014, 89, 014314. [Google Scholar] [CrossRef] [Green Version]
  24. Lonardoni, D.; Lovato, A.; Gandolfi, S.; Pederiva, F. Hyperon Puzzle: Hints from Quantum Monte Carlo Calculations. Phys. Rev. Lett. 2015, 114, 092301. [Google Scholar] [CrossRef] [Green Version]
  25. Petschauer, S.; Haidenbauer, J.; Kaiser, N.; Meissner, U.G.; Weise, W. Hyperons in nuclear matter from SU(3) chiral effective field theory. Eur. Phys. J. A 2016, 52, 15. [Google Scholar] [CrossRef] [Green Version]
  26. Logoteta, D.; Vidaña, I.; Bombaci, I. Impact of chiral hyperonic three-body forces on neutron stars. Eur. Phys. J. A 2019, 57, 207. [Google Scholar] [CrossRef]
  27. Champion, D.J.; Ransom, S.M.; Lazarus, P.; Camilo, F.; Bassa, C.; Kaspi, V.M.; Nice, D.J.; Freire, P.C.C.; Stairs, I.H.; van Leeuwen, J.; et al. An eccentric binary millisecond pulsar in the galactic plane. Science 2008, 320, 1309. [Google Scholar] [CrossRef] [Green Version]
  28. Demorest, P.B.; Pennucci, T.; Ransom, S.M.; Roberts, M.S.E.; Hessels, J.W.T. A two-solar-mass neutron star measured using Shapiro delay. Nature 2010, 467, 1081. [Google Scholar] [CrossRef]
  29. 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, 1233232. [Google Scholar] [CrossRef] [Green Version]
  30. 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. Astron. 2020, 4, 72. [Google Scholar] [CrossRef] [Green Version]
  31. Chatterjee, D.; Vidaña, I. Do hyperons exists in the interior of neutron stars? Eur. Phys. J. A 2016, 52, 29. [Google Scholar] [CrossRef] [Green Version]
  32. Vidaña, I. Hyperons: the strange ingredients of the nuclear equation of state. Proc. R. Soc. Lond. Ser. A 2018, 474, 20180145. [Google Scholar] [CrossRef]
  33. Kaplan, D.B.; Nelson, A.E. Strange goings on in dense nucleonic matter. Phys. Lett. B 1986, 175, 57–63, Erratum in 1986, 179, 409–410. [Google Scholar] [CrossRef]
  34. Brown, G.E.; Lee, C.H.; Rho, M.; Thorsson, V. From kaon-nuclear interactions to kaon condensation. Nucl. Phys. A 1994, 567, 937–956. [Google Scholar] [CrossRef] [Green Version]
  35. Thorsson, V.; Prakash, M.; Lattimer, J.M. Composition, structure and evolution of neutron stars with kaon condensates. Nucl. Phys. A 1994, 572, 693–731. [Google Scholar] [CrossRef] [Green Version]
  36. Lee, C.H. Kaon condensation in dense stellar matter. Phys. Rep. 1996, 275, 255–341. [Google Scholar] [CrossRef] [Green Version]
  37. Glendenning, N.K.; Schaffner-Bielich, J. Kaon Condensation and Dynamical Nucleons in Neutron Stars. Phys. Rev. Lett. 1998, 81, 4564–4567. [Google Scholar] [CrossRef] [Green Version]
  38. Tolos, L.; Fabbietti, L. Strangeness in nuclei and neutron stars. Prog. Part. Nucl. Phys. 2020, 112, 103770. [Google Scholar] [CrossRef] [Green Version]
  39. Schmitt, A.; Shternin, P. Reaction Rates and Transport in Neutron Stars. In The Physics and Astrophysics of Neutron Stars; Rezzolla, L., Pizzochero, P., Jones, D.I., Rea, N., Vidaña, I., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 455–574. [Google Scholar] [CrossRef] [Green Version]
  40. Yakovlev, D.G.; Shalybkov, D.A. Electrical conductivity of neutron star cores in the presence of a magnetic field—I. General solution for a multicomponent Fermi liquid—II.—A free particle model of npeΣ matter. Astroph. Space Sci. 1991, 176, 171–189. [Google Scholar] [CrossRef]
  41. Goldreich, P.; Reisenegger, A. Magnetic field decay in isolated neutron stars. Astrophys. J. 1992, 395, 250–258. [Google Scholar] [CrossRef] [Green Version]
  42. Dommes, V.A.; Gusakov, M.E.; Shternin, P.S. Dissipative relativistic magnetohydrodynamics of a multicomponent mixture and its application to neutron stars. Phys. Rev. D 2020, 101, 103020. [Google Scholar] [CrossRef]
  43. Baym, G.; Pethick, C. Landau Fermi-Liquid Theory: Concepts and Applications; John Wiley & Sons, Inc.: New York, NY, USA; Chichester, UK; Brisbane, Australia; Toronto, ON, Canada; Singapore, 1991. [Google Scholar] [CrossRef]
  44. Flowers, E.; Itoh, N. Transport properties of dense matter. II. Astrophys. J. 1979, 230, 847–858. [Google Scholar] [CrossRef]
  45. Anderson, R.H.; Pethick, C.J.; Quader, K.F. Transport properties of a multicomponent Fermi liquid. Phys. Rev. B 1987, 35, 1620–1629. [Google Scholar] [CrossRef]
  46. Shternin, P.S.; Baldo, M.; Haensel, P. Transport coefficients of nuclear matter in neutron star cores. Phys. Rev. C 2013, 88, 065803. [Google Scholar] [CrossRef] [Green Version]
  47. Shternin, P.S.; Baldo, M. Transport coefficients of nucleon neutron star cores for various nuclear interactions within the Brueckner-Hartree-Fock approach. Phys. Rev. D 2020, 102, 063010. [Google Scholar] [CrossRef]
  48. Baldo, M. (Ed.) Nuclear Methods and the Nuclear Equation of State; Volume 8, International Review of Nuclear Physics; World Scientific: Singapore, 1999. [Google Scholar] [CrossRef]
  49. Wiringa, R.B.; Stoks, V.G.J.; Schiavilla, R. Accurate nucleon-nucleon potential with charge-independence breaking. Phys. Rev. C 1995, 51, 38–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Carlson, J.; Pandharipande, V.R.; Wiringa, R.B. Three-nucleon interaction in 3-, 4- and N-body systems. Nucl. Phys. A 1983, 401, 59–85. [Google Scholar] [CrossRef]
  51. Stoks, V.G.J.; Rijken, T.A. Soft-core baryon-baryon potentials for the complete baryon octet. Phys. Rev. C 1999, 59, 3009–3020. [Google Scholar] [CrossRef] [Green Version]
  52. Rijken, T.A.; Stoks, V.G.J.; Yamamoto, Y. Soft-core hyperon-nucleon potentials. Phys. Rev. C 1999, 59, 21–40. [Google Scholar] [CrossRef] [Green Version]
  53. Shternin, P.; Baldo, M.; Schulze, H. Transport coefficients in neutron star cores in BHF approach. Comparison of different nucleon potentials. J. Phys. Conf. Ser. 2017, 932, 012042. [Google Scholar] [CrossRef]
  54. Shternin, P. Transport coefficients of leptons in superconducting neutron star cores. Phys. Rev. D 2018, 98, 063015. [Google Scholar] [CrossRef] [Green Version]
  55. Heiselberg, H.; Pethick, C.J. Transport and relaxation in degenerate quark plasmas. Phys. Rev. D 1993, 48, 2916–2928. [Google Scholar] [CrossRef] [PubMed]
  56. Heiselberg, H.; Baym, G.; Pethick, C.J.; Popp, J. Transport coefficients of relativistic plasmas. Nucl. Phys. A 1992, 544, 569–572. [Google Scholar] [CrossRef]
  57. Akmal, A.; Pandharipande, V.R.; Ravenhall, D.G. Equation of state of nucleon matter and neutron star structure. Phys. Rev. C 1998, 58, 1804–1828. [Google Scholar] [CrossRef] [Green Version]
  58. Fabrocini, F.; Fantoni, S. Correlated basis function results for the Argonne models of nuclear matter. Phys. Lett. B 1993, 298, 263–266. [Google Scholar] [CrossRef]
  59. Kadanoff, L.P.; Baym, G. Quantum Statistical Mechanics; Benjamin: New York, NY, USA, 1962. [Google Scholar]
  60. Kraeft, W.D.; Kremp, W.; Ebeling, W.; Röpke, G. Quantum Statistical of Charged Paerticle Systems; Akademie–Verlag: Berlin, Germany, 1986. [Google Scholar]
  61. Day, B.D. Elements of the Brueckner–Goldstone Theory of Nuclear Matter. Rev. Mod. Phys. 1967, 39, 719–744. [Google Scholar] [CrossRef]
  62. ter Haar, B.; Malfliet, R. Nucleons, mesons and deltas in nuclear matter a relativistic Dirac–Brueckner approach. Phys. Rep. 1987, 149, 207–286. [Google Scholar] [CrossRef]
  63. ter Haar, B.; Malfliet, R. Pion production, pion absorption, and nucleon properties in dense nuclear matter: Relativistic Dirac–Brueckner approach at intermediate and high energies. Phys. Rev. C 1987, 36, 1611–1620. [Google Scholar] [CrossRef] [PubMed]
  64. Brockman, R.; Machleidt, R. Relativistic nuclear structure. I. Nuclear matter. Phys. Rev. C 1990, 42, 1965–1980. [Google Scholar] [CrossRef]
  65. Song, H.Q.; Baldo, M.; Giansiracura, G.; Lombardo, U. Bethe–Brueckner–Goldstone expansion in nuclear matter. Phys. Rev. Lett. 1998, 81, 1584. [Google Scholar] [CrossRef]
  66. Baldo, M.; Giansiracusa, G.; Lombardo, U.; Song, H.Q. Bethe–Brueckner–Goldstone expansion in neutron matter. Phys. Lett. B 2000, 473, 1. [Google Scholar] [CrossRef]
  67. Loiseau, A.; Nogami, Y.; Ross, C.K. Nucleon-nucleon correlation and two-pion-exchange three-body force in nuclear matter. Nucl. Phys. A 1971, 165, 601. [Google Scholar] [CrossRef]
  68. Grangé, P.; Martzolff, M.; Nogami, Y.; Sprung, D.W.L.; Ross, C.K. Three-body force in nuclear matter. Phys. Lett. B 1976, 60, 237. [Google Scholar] [CrossRef]
  69. Baldo, M.; Ferreira, L. Nuclear liquid-gas phase transition. Phys. Rev. C 1999, 59, 682. [Google Scholar] [CrossRef]
  70. Baldo, M.; Burgio, G.F.; Schulze, H.J.; Taranto, G. Nucleon effective masses within the Brueckner-Hartree-Fock theory: Impact on stellar neutrino emission. Phys. Rev. C 2014, 89, 048801. [Google Scholar] [CrossRef] [Green Version]
  71. Shternin, P.S.; Yakovlev, D.G. Electron-muon heat conduction in neutron star cores via the exchange of transverse plasmons. Phys. Rev. D 2007, 75, 103004. [Google Scholar] [CrossRef] [Green Version]
  72. Shternin, P.S.; Yakovlev, D.G. Shear viscosity in neutron star cores. Phys. Rev. D 2008, 78, 063006. [Google Scholar] [CrossRef] [Green Version]
  73. Shternin, P.S. Exchange of transverse plasmons and electrical conductivity of neutron star cores. JETP 2008, 107, 212–227. [Google Scholar] [CrossRef]
  74. Sedrakian, A.; Clark, J.W. Superfluidity in nuclear systems and neutron stars. Eur. Phys. J. A 2019, 55, 167. [Google Scholar] [CrossRef] [Green Version]
  75. Balberg, S.; Barnea, N. S-wave pairing of Λ hyperons in dense matter. Phys. Rev. C 1998, 57, 409–416. [Google Scholar] [CrossRef] [Green Version]
  76. Takatsuka, T.; Tamagaki, R. Superfluidity of Λ hyperons admixed in neutron star cores. Prog. Theor. Phys. 1999, 102, 1043–1048. [Google Scholar] [CrossRef] [Green Version]
  77. Takatsuka, T.; Nishizaki, S.; Yamamoto, Y.; Tamagaki, R. Possibility of hyperon superfluidity in neutron star cores. Prog. Theor. Phys. 2000, 107, 179–184. [Google Scholar] [CrossRef] [Green Version]
  78. Takatsuka, T.; Nishizaki, S.; Yamamoto, Y.; Tamagaki, R. Superfluidity of hyperon-mixed neutron stars. Prog. Theor. Phys. Suppl. 2002, 146, 279–288. [Google Scholar] [CrossRef] [Green Version]
  79. Vidaña, I.; Tolós, L. Superfluidity of Σ hyperons in β-stable neutron star matter. Phys. Rev. C 2004, 70, 02882. [Google Scholar] [CrossRef] [Green Version]
  80. Zhou, X.R.; Schulze, H.J.; Pan, E.; Drayer, J.P. Strong hyperon-nucleon pairing in neutron stars. Phys. Rev. Lett. 2005, 95, 051101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Wang, Y.N.; Shen, H. Superfluidity of Λ hyperons in neutron stars. Phys. Rev. C 2010, 81, 025801. [Google Scholar] [CrossRef] [Green Version]
  82. Andersson, N. A Superfluid Perspective on Neutron Star Dynamics. Universe 2021, 7, 17. [Google Scholar] [CrossRef]
  83. Baiko, D.A.; Haensel, P.; Yakovlev, D.G. Thermal conductivity of neutrons in neutron star cores. Astron. Astrophys. 2001, 374, 151–163. [Google Scholar] [CrossRef] [Green Version]
  84. Vollhardt, D.; Wölfle, P. The Superfluid Phases of Helium 3; Taylor & Francis: Bristol, UK, 1990. [Google Scholar] [CrossRef]
  85. Manuel, C.; Tolos, L. Transport Properties of Superfluid Phonons in Neutron Stars. Universe 2021, 7, 59. [Google Scholar] [CrossRef]
Figure 1. Properties of β -stable hypernuclear matter. Panel (a): Pressure as function of a total energy density for the NSC97e (solid lines) and NSC97a (dashed lines) models. Panel (b): Particle fractions as a function of total baryon number density for the NSC97e (solid lines) and NSC97a (dashed lines) models. Vertical dashed lines indicate μ , Σ , and Λ thresholds.
Figure 1. Properties of β -stable hypernuclear matter. Panel (a): Pressure as function of a total energy density for the NSC97e (solid lines) and NSC97a (dashed lines) models. Panel (b): Particle fractions as a function of total baryon number density for the NSC97e (solid lines) and NSC97a (dashed lines) models. Vertical dashed lines indicate μ , Σ , and Λ thresholds.
Universe 07 00203 g001
Figure 2. Effective baryon masses on the Fermi surface as a function of the total baryon number density for the NSC97e (solid lines) and NSC97a (dashed lines) models.
Figure 2. Effective baryon masses on the Fermi surface as a function of the total baryon number density for the NSC97e (solid lines) and NSC97a (dashed lines) models.
Universe 07 00203 g002
Figure 3. Diagonal elements of the transport matrix for baryons as a function of the total baryon number density. Different line types correspond to different particles, as indicated in the legend. Only the strong interaction is taken into account in this plot and temperature-independent ( Λ c c T 8 2 ) combinations are shown. Results are shown for the NSC97e model, respectively, for the thermal conductivity (panel (a)) and the shear viscosity cases (panel (b)).
Figure 3. Diagonal elements of the transport matrix for baryons as a function of the total baryon number density. Different line types correspond to different particles, as indicated in the legend. Only the strong interaction is taken into account in this plot and temperature-independent ( Λ c c T 8 2 ) combinations are shown. Results are shown for the NSC97e model, respectively, for the thermal conductivity (panel (a)) and the shear viscosity cases (panel (b)).
Universe 07 00203 g003
Figure 4. Normalized elements of the transport matrix for baryons for the thermal conductivity case and the NSC97e model. Solid lines in each row (c = n, p, Λ , Σ ) show relative contributions Υ ˜ c i κ to diagonal transport matrix elements from collisions with different baryon species. Other lines show the non-diagonal relative matrix elements Λ ˜ c i κ , where red dashed lines correspond to positive values and the blue dash-dotted ones to negative values, respectively. Only strong forces are included, see text for details.
Figure 4. Normalized elements of the transport matrix for baryons for the thermal conductivity case and the NSC97e model. Solid lines in each row (c = n, p, Λ , Σ ) show relative contributions Υ ˜ c i κ to diagonal transport matrix elements from collisions with different baryon species. Other lines show the non-diagonal relative matrix elements Λ ˜ c i κ , where red dashed lines correspond to positive values and the blue dash-dotted ones to negative values, respectively. Only strong forces are included, see text for details.
Universe 07 00203 g004
Figure 5. Normalized elements of the transport matrix for baryons for the shear viscosity case and NSC97e model. Solid lines in each row (c = n, p, Λ , Σ ) show relative contributions Υ ˜ c i η to diagonal transport matrix elements. Other lines show the non-diagonal relative matrix elements Λ ˜ c i η , where red dashed lines correspond to positive values and the blue dash-dotted ones to negative values, respectively. Only strong forces are included, see text for details.
Figure 5. Normalized elements of the transport matrix for baryons for the shear viscosity case and NSC97e model. Solid lines in each row (c = n, p, Λ , Σ ) show relative contributions Υ ˜ c i η to diagonal transport matrix elements. Other lines show the non-diagonal relative matrix elements Λ ˜ c i η , where red dashed lines correspond to positive values and the blue dash-dotted ones to negative values, respectively. Only strong forces are included, see text for details.
Universe 07 00203 g005
Figure 6. Corrections to the variational solution for the n Λ subsystem (obtained with the NSC97e model). Solid and dashed lines show correction coefficients for the thermal conductivity and shear viscosity cases, respectively.
Figure 6. Corrections to the variational solution for the n Λ subsystem (obtained with the NSC97e model). Solid and dashed lines show correction coefficients for the thermal conductivity and shear viscosity cases, respectively.
Universe 07 00203 g006
Figure 7. Partial contributions to the thermal conductivity coefficient (variational solution) for the β -stable matter in NS cores. Results for the NSC97e model are shown. Particle species labels are shown near the corresponding curves. The total thermal conductivity including corrections to the variational solution is shown with the thick solid line. Results for three values of the temperature are shown: T = 10 7 K (a), T = 10 8 K (b) and T = 10 9 K (c).
Figure 7. Partial contributions to the thermal conductivity coefficient (variational solution) for the β -stable matter in NS cores. Results for the NSC97e model are shown. Particle species labels are shown near the corresponding curves. The total thermal conductivity including corrections to the variational solution is shown with the thick solid line. Results for three values of the temperature are shown: T = 10 7 K (a), T = 10 8 K (b) and T = 10 9 K (c).
Universe 07 00203 g007
Figure 8. Partial contributions to the shear viscosity coefficient (variational solution) for the beta-stable mater in NS cores. Results for the NSC97e model are shown. Particle species labels are shown near the corresponding curves. The total shear viscosity including corrections to the variational solution is shown with the thick solid line. Results for three values of the temperature are shown: T = 10 7 K (a), T = 10 8 K (b) and T = 10 9 K (c).
Figure 8. Partial contributions to the shear viscosity coefficient (variational solution) for the beta-stable mater in NS cores. Results for the NSC97e model are shown. Particle species labels are shown near the corresponding curves. The total shear viscosity including corrections to the variational solution is shown with the thick solid line. Results for three values of the temperature are shown: T = 10 7 K (a), T = 10 8 K (b) and T = 10 9 K (c).
Universe 07 00203 g008
Figure 9. Total thermal conductivity (a) and shear viscosity (b) in hypernuclear NS cores for the NSC97e (solid lines) and NSC97 (dashed lines) models. Here we do not perform temperature rescaling of the transport coefficients. The results are shown for three values of temperature, T = 10 7 , 10 8 , and 10 9 K as indicated in the plot.
Figure 9. Total thermal conductivity (a) and shear viscosity (b) in hypernuclear NS cores for the NSC97e (solid lines) and NSC97 (dashed lines) models. Here we do not perform temperature rescaling of the transport coefficients. The results are shown for three values of temperature, T = 10 7 , 10 8 , and 10 9 K as indicated in the plot.
Universe 07 00203 g009
Figure 10. Momentum transfer rates J c i mediated by the strong interactions for the NSC97e (panel a) and NSC97a (panel b) models. The temperature-independent combination J c i T 8 2 is shown. Different lines correspond to different baryon pairs as shown in the legend. Vertical dashed lines indicate Σ and Λ thresholds.
Figure 10. Momentum transfer rates J c i mediated by the strong interactions for the NSC97e (panel a) and NSC97a (panel b) models. The temperature-independent combination J c i T 8 2 is shown. Different lines correspond to different baryon pairs as shown in the legend. Vertical dashed lines indicate Σ and Λ thresholds.
Universe 07 00203 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shternin, P.; Vidaña, I. Transport Coefficients of Hyperonic Neutron Star Cores. Universe 2021, 7, 203. https://doi.org/10.3390/universe7060203

AMA Style

Shternin P, Vidaña I. Transport Coefficients of Hyperonic Neutron Star Cores. Universe. 2021; 7(6):203. https://doi.org/10.3390/universe7060203

Chicago/Turabian Style

Shternin, Peter, and Isaac Vidaña. 2021. "Transport Coefficients of Hyperonic Neutron Star Cores" Universe 7, no. 6: 203. https://doi.org/10.3390/universe7060203

APA Style

Shternin, P., & Vidaña, I. (2021). Transport Coefficients of Hyperonic Neutron Star Cores. Universe, 7(6), 203. https://doi.org/10.3390/universe7060203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop