Next Article in Journal
On Cooling of Neutron Stars with a Stiff Equation of State Including Hyperons
Previous Article in Journal
The Black Hole Radiation in Massive Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

On Manifestation of In-Medium Effects in Neutron Stars and Heavy-Ion Collisions

by
Dmitry N. Voskresensky
1,2
1
National Research Nuclear University (MEPhI), Kashirskoe shosse 31, 115409 Moscow, Russia
2
Joint Institute for Nuclear Research, RU-141980 Dubna, Moscow, Russia
Universe 2018, 4(2), 28; https://doi.org/10.3390/universe4020028
Submission received: 27 November 2017 / Revised: 3 January 2018 / Accepted: 29 January 2018 / Published: 8 February 2018

Abstract

:
This review focuses on the demonstration of an interrelation between various in-medium effects, which are manifested in the phenomena occurring in neutron stars and heavy-ion collisions. More specifically, the equation of state of a baryon-rich cold hadron matter is considered. It is done within the relativistic mean-field approach with σ field-scaled hadron masses and couplings, and a cut-mechanism is discussed leading to an increase of the stiffness of the dense baryon matter. Then, I discuss the role of the viscosity and thermal conductivity in description of the first-order phase transitions occurring in heavy-ion collisions and neutron stars. Next, the p-wave polarization effects on pion and kaon spectra are studied beyond the mean-field level. In particular, the pion softening effect is detailed. Then, a role of in-medium effects in neutrino radiation of neutron stars is discussed and effects of the bulk and shear viscosities in the problem of r-mode damping in young rapidly rotating pulsars are considered.

1. Introduction

Often, researchers studying heavy-ion collisions, where one deals with a non-equilibrium and a quasi-equilibrium hot and dense hadron matter, do not care about possible application of their models in other branches of nuclear physics, e.g., in neutron star physics, where one deals with the nuclear matter but in cold state, and, vice versa, the researchers considering neutron-star phenomena often ignore information known from heavy-ion-collision physics. The aim of this contribution is to partially overcome this gap demonstrating an interrelation between description of various phenomena in the heavy-ion-collision physics and the neutron-star physics. In the phenomena, which will be reviewed, the key role is played by in-medium effects in strongly-interacting nuclear matter.

2. EoS of Baryon-Rich Dense Not Too Hot Matter

2.1. Existing Constraints on EoS

Equation of state (EoS) of the cold hadronic matter should:
  • satisfy experimental information on properties of dilute nuclear matter and not contradict the results of microscopical approaches like the Brueckner–Hartree–Fock (BHF) one;
  • place empirical constraints on global characteristics of atomic nuclei [1];
  • place constraints on the pressure of the nuclear mater from the description of particle transverse and elliptic flows and the K + production in heavy-ion collisions, cf. [2,3];
  • allow for the heaviest known compact stars [4];
  • allow for an adequate description of the compact star cooling in the absence of the direct Urca (DU) neutrino processes in the majority of the known pulsars detected in soft X rays [1];
  • and yield a mass–radius relation comparable with the empirical constraints including recent gravitation wave LIGO–Virgo detection [5];
  • being extended to non-zero temperature T (for T < T c , where T c is the critical temperature of the deconfinement), appropriately describe supernova explosions, proto-neutron stars, and heavy-ion collision data [6,7], etc.
For more details, see discussion in [1,6,7,8,9,10,11,12,13,14].
Most difficult is the ability to satisfy simultaneously the heavy-ion-collision flow and the maximum neutron-star mass constraints. Flow of matter in heavy-ion-collision is directed both forward and perpendicular to the beam axis. At high densities, spectator nucleons may shield the transversal flow into their direction and generate an inhomogeneous density and pressure profile in the transversal plane. This effect is usually referred to as elliptic flow and depends on the EoS. An analysis of the nucleon flow data, which depends essentially only on the isospin-independent part of the EoS, was carried out in [2]. The flow constraint [2] to be fulfilled requires rather soft EoS of isospin-symmetric matter (ISM), whereas the EoS of the beta-equilibrium matter (BEM) should be stiff so that maximum available mass of the neutron star would be higher than the measured mass M = 2.01 ± 0.04   M [4] of the pulsar PSR J0348+0432, being the heaviest at present. Various EoSs in comparison with the flow and neutron-star mass constraints are shown in Figure 1 taken from [1]. Models DD, D 3 C and DBHF (Bonn A) satisfy well the maximum neutron-star mass constraint but do not fulfill the flow constraint. On the contrary, the models NL ρ , NL ρ δ satisfy the flow constraint but do not fulfill the maximum neutron-star mass constraint, and only models with σ -field-dependent effective hadron masses and couplings like KVR, KVOR, and the DD-F relativistic mean-field (RMF) model with density-dependent couplings are able to satisfy the flow constraint and marginally the maximum neutron-star mass constraint. However, hyperons and Δ isobars are not incorporated within these models. Their inclusion additionally softens the EoSs [15,16,17] that may cause a problem with the fulfillment of the maximum neutron-star mass constraint. In conclusion, in-medium effects, e.g., those resulting in a modification of the hadron effective masses and couplings in matter, are required to fulfill both the flow and the maximum neutron-star mass constraints.

2.2. Cut-Mechanism for Stiffening of EoS

Ref. [12] studied various RMF models, where nucleons interact with scalar, vector and iso-vector mean fields and self- and cross- mean-field interaction terms might be also included. Usually, in such models, the magnitude of the scalar field increases monotonically with the nucleon density increase and the nucleon effective mass decreases. The latter quantity stops decreasing and the EoS stiffens, provided the mean-field self-interaction potential rises sharply in a narrow vicinity of the values of mean fields corresponding to nucleon densities n n > n 0 , where n 0 is the nuclear saturation density. Such a behaviour might be related to an excluded volume effect due to a strong repulsion at short distances between hadrons. As a result, the maximum neutron-star mass increases. Such a “cut”-procedure offered a simple way to stiffen the EoS at densities above n without altering it at n n 0 . The developed scheme allows an application to neutron stars of the RMF models, which are well fitted to finite nuclei but do not fulfill the experimental constraint on the maximum neutron-star mass. Then, the energy density is given by
E ( n p , n n , f ) = E kin ( n p , n n , f ) + E V ( n p , n n ) + E σ ( f ) , E kin ( n p , n n , f ) = ( 0 p F n + 0 p F p ) p 2 d p π 2 p 2 + m N 2 ( f ) , E V ( n p , n n ) = C ω 2 2 m N 2 ( n p + n n ) 2 + C ρ 2 8 m N 2 ( n p n n ) 2 , E σ ( f ) = m N 4 f 2 2 C σ 2 + U ( f ) , m N ( f ) = m N ( 1 f ) .
C j = g j N m N / m j , j = ( σ , ω , ρ ) are fitting parameters of the RMF model, g j N are meson-nucleon coupling constants, m j are masses of j-hadrons, m N is the nucleon mass, p F n and p F p are Fermi momenta of neutrons and protons, f = g σ N σ / m N , U is the effective potential; in ordinary nonlinear Walecka (NLW) models, one chooses U ( f ) = m N 4 ( b f 3 / 3 + c f 4 / 4 ) , where b and c are fitting parameters. We also exploit a modified potential U ˜ depending on a c σ parameter such that U ˜ undergoes a significant growth for f > f c u t (NLWcut model). The scalar field potential U ( f ) in the NLW model, and U ˜ ( f ) in the NLWcut model for various values c σ , as functions of the scalar field parameter f, are shown in Figure 2a. The right panel demonstrates the nucleon effective mass in the ISM as a function of the nucleon density for the same models.
In the left panel in Figure 3, we show the pressure P for the NLW and NLWcut models under consideration as a function of the nucleon density in the ISM for various values c σ , and the insertion demonstrates the binding energy per nucleon. The neutron-star mass and radius as functions of the central density for the same models are shown in the right panel. We see that the initial NLW model satisfies the flow constraint but does not fulfill the maximum neutron-star mass constraint, whereas two of four presented cut-versions of the model fulfill properly both constraints and one model fulfills the maximum neutron-star mass constraint marginally.

2.3. Hyperon and Δ -Isobar Puzzles and Their Resolution in RMF Models with Scaled Hadron Masses and Couplings

In the standard RMF models, hyperons and Δ -isobars may appear in neutron-star cores already for n (2–3) n 0 , which results in a decrease of the maximum neutron-star mass below the observed limit. The problems were named the hyperon [15,16] and the Δ [17,18] puzzles. These problems can be resolved in the RMF models with scaled hadron masses and couplings [9,10,11,14]. The model employed in [11,14] is a generalization of the NLW model with effective coupling constants g m b = g m b χ m b ( σ ) and hadron masses m i = m i Φ i ( σ ) , which are functions of the σ field [8], i = { b , m } , where m = { σ , ω , ρ , ϕ } lists the included mesonic fields, b = { N , H , Δ } indicates a baryon (nucleon N = p , n ; hyperon H = Λ , Σ , Ξ ; and Δ isobar). Quantities χ m b ( σ ) and Φ i ( σ ) are dimensionless scaling functions of the σ mean field. In the RMF approximation, the contribution of Δ s to the energy density has the same form as for spin- 1 / 2 fermions but with the spin degeneracy 4. With the standard solutions for the ω and ρ meson mean fields, being expressed through baryon densities, the energy density of the model takes the following form
E [ { n b } , { n l } , f ] = b E kin ( p F , b , m b Φ b ( f ) , s b ) + l = e , μ E kin ( p F , l , m l , s l ) + m N 4 f 2 2 C σ 2 η σ ( f ) + 1 2 m N 2 C ω 2 n V 2 η ω ( f ) + C ρ 2 n I 2 η ρ ( f ) + C ϕ 2 n S 2 η ϕ ( f ) , n V = b x ω b n b , n I = b x ρ b t 3 b n b , n S = H x ϕ H n H ,
where we introduced the dimensionless scalar field f = g σ N χ σ N ( σ ) σ / m N , the coupling constant ratios x M b = g M b / g M N , M = { σ , ω , ρ } , x ϕ H = g ϕ H / g ω N , g ϕ N = g ϕ Δ = 0 , and t 3 b is the isospin projection of baryon b; the Fermi momentum is related to the fermion density p F , j = ( 6 π 2 n j / ( 2 s j + 1 ) ) 1 / 3 with s j standing for spin of fermion j = ( b , l ) , l = ( e , μ ) ;
C M = g M N m N m M , C ϕ = g ω N m N m ϕ .
The scaling functions are
η ω ( f ) = Φ ω 2 ( f ) χ ω N 2 ( f ) , η ρ ( f ) = Φ ρ 2 ( f ) χ ρ N 2 ( f ) , η ϕ ( f ) = Φ ϕ 2 ( f ) χ ϕ H 2 ( f ) , η σ ( f ) = Φ σ 2 [ σ ( f ) ] χ σ N 2 [ σ ( f ) ] + 2 C σ 2 m N 4 f 2 U [ σ ( f ) ] .
As functions of f, they are constructed to fulfill all necessary constraints. The model uses universal mass scaling for nucleons and mesons, Φ N = Φ m = 1 f , and universal scaling for coupling constants χ ω b ( f ) = χ ω N ( f ) and χ ρ b ( f ) = χ ρ N ( f ) , the coupling constant for ϕ meson remains unscaled, χ ϕ H ( f ) = 1 , which results in the ϕ scaling function η ϕ = ( 1 f ) 2 . The scaling η σ ( f ) includes also the scalar meson self-interaction potential U ( σ ) . The vector-meson coupling ratios with baryons x ω B and x ρ B are chosen following the quark SU(6) symmetry. The baryon coupling ratios with the scalar field x σ B are deduced from the potentials U B ( n 0 ) = C ω 2 m N 2 x ω B n 0 x σ B ( m N + m N Φ N ( n 0 ) ) in ISM at the saturation baryon density n = n 0 . There are large uncertainties in fixing of the Δ potential U Δ ( n 0 ) U Δ in the literature. The value U Δ = 50 MeV, which results from calculations [19], is used as the most realistic choice but also a variation of U Δ is performed. The detailed description of parameters, which we use for all included baryon species, can be found in [11].
Refs. [9,10,11,14] use models of two kinds, the KVORcut-based models, where one exploits the cut in ω sector choosing accordingly the η ω ( f ) dependence, and the MKVORcut-based models, where one exploits the cut in ρ sector choosing accordingly the η ρ ( f ) dependence. In KVORcut models without inclusion of hyperons, the ρ condensate appears by the second-order phase transition and the influence of the ρ condensation on the EoS is minor. In MKVOR-based models, the ρ condensate (see in Section 2.4) appears by the first-order phase transition and the influence of the ρ condensation on the EoS is substantially stronger. In Figure 4, we show neutron-star masses as functions of the radii for MKVOR*-based models [14]. In MKVOR* models, we prevent the effective nucleon mass from vanishing at any density. In the MKVOR* Δ model, Δ isobars are included; in MKVOR*H Δ ϕ model hyperons, Δ -isobars and ϕ -mesons are included; in the model MKVOR* Δ ρ , Δ and ρ condensate are included; and in MKVOR*H Δ ϕ ρ model hyperons, Δ , ϕ and ρ are included, respectively. As we see, in the presence of Δ -isobars, the ρ condensate softens EoS substantially, but, even with its inclusion, we are able to satisfy the maximum neutron-star mass constraint.

2.4. Charged ρ Condensation

The non-Abelian ρ meson field is analogous to the massive gluon field. Ref. [20] demonstrated that, in a dense nucleon isospin-asymmetric matter at an appropriate decrease of the effective ρ meson mass, there may occur a charged ρ condensation. Then, the phenomenon of the charged ρ condensation was incorporated in the RMF models with scaled hadron masses and coupling constants in [8]. Similarly, the charged ρ condensation may appear in very strong magnetic fields [21,22].
The contribution of the charged ρ meson condensate to the energy density is then given by [14]
Δ E ch , ρ [ { n b } ; f ] = C ρ 2 2 m N 2 η ρ n I + n ρ 2 θ ( n I n ρ ) μ ch , ρ n ch , ρ ,
where θ ( n I n ρ ) = 1 for n I + n ρ < 0 and zero otherwise. The solution for n I + n ρ < 0 corresponds to ρ 0 ( 3 ) = μ ch , ρ m ρ Φ ρ g ρ χ ρ , | ρ c | 2 = ( n I n ρ ) θ ( n I n ρ ) 2 m ρ η ρ 1 / 2 χ ρ , where n ρ = a ( m ρ Φ ρ μ ch , ρ ) , a = m N 2 η ρ 1 / 2 Φ ρ C ρ 2 χ ρ . The ρ charged density is n ch , ρ = 2 m ρ Φ ρ ( f ) | ρ c | 2 < 0 . The scaling of the non-abelian ρ - ρ interaction χ ρ is put unity.

3. Viscosity and Thermal Conductivity Effects at First-Order Phase Transitions in Heavy-Ion Collisions and Neutron Stars

3.1. Pressure–Baryon Density Behavior at First-Order Phase Transitions

At various first-order phase transitions, the pressure–density behaviour gets a van der Waals form. It occurs for the nuclear liquid–vapor phase transition, whose signatures are manifested in low-energy heavy-ion collisions, and might be for quark–hadron phase transition in baryon-rich matter, for phase transitions to the pion- and kaon- and charged ρ -condensate states in rather massive neutron stars and might be for nucleon- Δ -resonance state transition in ISM, as follows in some models (cf. Figure 6, which will be demonstrated and discussed below).
In Figure 5, we show typical behaviour of the pressure isotherms as functions of the density at the first-order liquid–vapor-like phase transition. The horizontal dashed line shows the Maxwell construction (MC). At equilibrium, the chemical potentials μ A = μ D . The interval AB corresponds to a metastable supercooled vapor (SV) and the interval CD, to a metastable overheated liquid (OL). The interval BC corresponds to an unstable spinodal region.
Calculations [11] show that, with increase of the baryon density in the cold ISM, there may appear a phase transition to the Δ -rich matter. Pressures as functions of the density in the ISM for MKVOR* model and MKVOR* Δ model are shown in Figure 6. As is seen, for the values of the Δ optical potential U Δ ( n = n 0 ) > 56 MeV, the transition proves to be of the second-order and, for U Δ ( n = n 0 ) < 56 MeV, it is of the first-order. Although the most realistic value for the Δ optical potential is estimated as U Δ ( n = n 0 ) 50 MeV, a variation is allowed. In the case of the first-order phase transition, the pressure isotherms as functions of the baryon density demonstrate a back-bending. There are two narrow spinodal regions. A specifics of the back-bending is described in [11]. Note that the flow constraint [2] is better fulfilled in the presence of the first-order transition to the Δ -rich matter. A first-order phase transition into a Δ -rich state that we obtained within the MKVOR* Δ model for U Δ < 56 MeV could manifest itself in an increase of the pion yield at typical energies and momenta corresponding to the Δ π N decays in heavy-ion collision experiments.
In neutron stars, appearance of a strong phase transition may result in a second neutrino burst if the transition occurred right after supernova explosion and a hot neutron star formation. A possibility of the delayed second neutrino burst in SN 1987A event has been considered in [23,24,25]. It might be associated with a significant delay of the heat transport to the neutron-star surface if the system is close to the pion–condensate phase transition. Recently, new arguments have been expressed for that namely two neutrino bursts were measured during 1987A explosion, one delayed respectively the other one by 4.7 h, cf. [26]. The second burst could then be related to the phase transition of the neutron star to the pion condensate state. In addition, a phase transition to the pion- or to the kaon–condensate state could occur during 10-s-period of the neutron star formation. Moreover, the first-order phase transition, if it occurs, could result in a blowing off a star matter. In old neutron stars, the first-order phase transition, if it occurs, could result in a strong star-quake [25,27].

3.2. Hydrodynamical Description

Here, we follow the analysis of works [28,29,30,31]. Assume that a nuclear system is rather close to a critical point of a liquid–vapor-type first-order phase transition. Then, the velocity of a fluctuation (seed), u , is much less than the mean thermal velocity and we may use equations of non-relativistic non-ideal hydrodynamics: the Navier–Stokes equation, the continuity equation, and equation for the heat transport:
m n t u i + ( u ) u i = i P + k η k u i + i u k 2 d δ i k div u + ζ δ i k div u ,
t n + div ( n u ) = 0 ,
T s t + div ( s u ) = div ( κ T ) + η k u i + i u k 2 d δ i k div u 2 + ζ ( div u ) 2 .
n is the density of the conserving charge (here, the baryon density), m is the effective baryon mass, P is the pressure, η and ζ are the shear and bulk viscosities, d shows the geometry of the fluctuation (droplets, rods, slabs), T is the temperature, s is the entropy density, and κ is the thermal conductivity.
As an example, we consider a van der Waals fluid. Then, adiabatic trajectories, s ˜ s / n const, for an expansion of a uniform fireball to vacuum, are shown in Figure 7. The SV and the OL regions are between the MC and the isothermal spinodal (ITS) curves, on the left and on the right, respectively. The adiabatic spinodal (AS) curve bounds the AS region from above. For s ˜ c r > s ˜ > s ˜ MC 2 , where s ˜ c r corresponds to the value of the specific entropy s ˜ at the critical point and the line with s ˜ MC 2 passes through the point n / n c r = 3 at T = 0 , the system traverses the OL state (the region OL in Figure 7), the ITS region (below the ITS line) and the AS region (below the AS line). For s ˜ > s ˜ c r , the system trajectory passes through the SV state (the region SV in Figure 7) and the ITS region.
Note that for the quark–hadron first-order phase transition the phase diagram looks a bit different since then T c increases with a decrease of the baryon density [32,33]. However presence of this peculiarity does not change our general analysis given here.
Note that, for the hadron–quark transition in heavy-ion collisions, besides a possibility of the first-order phase transition, one also discusses a possibility of the crossover. The latter transition is expected at collision energies achievable at RHIC (Brookhaven) and LHC (CERN). Lattice calculations performed for a baryon-less matter also indicate on the crossover transition. For lower collision energies relevant for NICA (Dubna) and FAIR (Darmstadt) facilities, one expects to find signatures of the first-order phase transition. Although it is not excluded that even in the NICA-FAIR range of the collision energies the hadron–quark transition is the crossover, in this review, I focus only on the description of the signatures of the first-order phase transition.
All thermodynamic quantities can be expanded near an arbitrary reference point ( n r , T r ) , or ( n r , s ˜ r ) , which we assume to be close to the critical point. The Landau free energy in the variables δ n = n n c r , δ T = T T c r , δ ( δ F L ) / δ ( δ n ) = P P f + P MC is as follows:
δ F L = d 3 x n c r c m [ ( δ n ) ] 2 2 + λ m 3 ( δ n ) 4 4 λ v 2 m ( δ n ) 2 2 ϵ δ n ,
where ϵ = P f P MC is expressed through the pressure at the MC. The maximum of the quantity ϵ is ϵ m = 4 λ v 3 / ( 3 3 ) . The first term in Equation (8) is due to the surface tension, δ F L , surf = σ S , S is the surface of the seed. For the van der Waals equation of state:
v 2 ( T ) = 4 δ T n c r 2 m 2 T c r , σ = σ 0 | δ T | 3 / 2 T c r 3 / 2 , σ 0 = 32 m n c r 2 T c r c .
Using Equation (6), and keeping only linear terms in u, since near the critical point processes develop slowly, v 2 δ T , we rewrite Equation (5) in the dimensionless variables δ ρ = v ψ , ξ i = x i / l , i = 1, ⋯, d, d = 3 for seeds of spherical geometry, τ = t / t 0 as
β 2 ψ τ 2 = Δ ξ Δ ξ ψ + 2 ψ ( 1 ψ 2 ) + ϵ ˜ ψ τ , l = 2 c λ v 2 1 / 2 , t 0 = 2 ( d ˜ η r + ζ r ) λ v 2 n c r m c r , ϵ ˜ = 2 ϵ λ v 3 , β = c n c r 2 m c r 2 ( d ˜ η r + ζ r ) 2 ,
and d ˜ = 2 ( d 1 ) / d . One often speculates about a manifestation of fluctuations as a measure of the closeness to the critical quark–hadron first-order phase transition point in the course of heavy-ion collisions. As we demonstrate, processes near the critical point are frozen (critical slowing down). Thereby, at least a region very close to the critical point can hardly be manifested via abnormal fluctuations in the course of heavy-ion collisions since these fluctuations do not have enough time to develop.
Equation (10) is supplemented by Equation (7) for the heat transport, which owing to Equation (6) after linearization simplifies to
T c r t δ s s c r ( n c r ) 1 t δ n = κ r Δ δ T .
The variation of the temperature is related to the variation of the entropy density s [ n , T ] by δ T T c r ( c V , r ) 1 δ s ( s / n ) T , c r δ n . Note that the reference point differs from the critical point, although we take the reference point in the vicinity of the critical point. This circumstance is important for the determination of the specific heat density c V , r and transport coefficients, whereas other quantities are already smooth functions of n , T .
Equation (10) differs from the standard Ginzburg–Landau equation exploited in condensed matter theory. The difference disappears, if one puts the bracket-term in the r.h.s. of Equation (10) to zero. This procedure is not legitimate for consideration of fluctuations at a short time. Thereby, there exists an initial stage of the dynamics of seeds ( t t init ), which is not described by the standard Ginzburg–Landau equation [28,29].
The time scale for the relaxation of the density following Equation (10) is t ρ R , where R is the size of a seed (after awhile, an overcritical seed begins to grow with a constant velocity already following the standard Ginzburg–Landau equation), and the time scale for the relaxation of the entropy/temperature, following Equation (11), is
t T = R 2 c V , r / κ r R 2 .
The evolution of a seed is governed by the slowest mode. Thus, for t T ( R ) < t ρ ( R ) , i.e., for R < R fog ( R fog is the typical size of the seed at which t ρ = t T ), dynamics of seeds is controlled by Equation (10) for the density mode. For seeds with sizes R > R fog , t T R 2 exceeds t ρ R and growth of seeds is slowed down. Thereby, the number of seeds with the size R R fog may increase with time. For the quark–hadron phase transition R fog 0.1–1 fm and for the nuclear liquid–vapor transition R fog 1–10 fm R fb ( t f . o . ) , where R fb ( t f . o . ) is the fireball size at the freeze out, t f . o . is the fireball evolution time till freeze out. Thus, thermal conductivity effects may manifest themselves in heavy-ion collision dynamics.
There are only two dimensionless parameters in Equation (10), ϵ ˜ and β . The parameter ϵ ˜ is responsible for a difference between the Landau free energies of the metastable and stable states. For t ρ t T (isothermal stage), ϵ ˜ c o n s t and dependence on this quantity disappears because of Δ ξ ϵ ˜ 0 . Then, dynamics is controlled by the parameter β , which characterizes an inertia. It is expressed in terms of the surface tension and the viscosity as
β = ( 32 T c r ) 1 [ d ˜ η r + ζ r ] 2 σ 0 2 m .
The larger is the viscosity and the smaller is the surface tension, the effectively more viscous (inertial) is the fluidity of seeds. For β 1 , one deals with the regime of effectively viscous fluid and at β 1 , with the regime of perfect fluid, for the nuclear liquid–vapor phase transition typically β 0.01 , and for the quark–hadron phase transition β 0.02–0.2, even for very low value of the ratio η / s 1 / ( 4 π ) . The latter quantity characterizes fluidity ultra-relativistic systems where viscosity proves to be very small but finite [34]. In the case of baryon-rich matter, one deals with effectively very viscous fluidity of density fluctuations in both nuclear liquid–vapor and quark–hadron phase transitions. An interplay between viscosity, surface tension, and thermal conductivity effects is responsible for the typical size of the fluctuations.
In neutron stars, an overcritical pion–condensate drop reaches a size R 0.1 km for t 10 3 s by the growth of the density mode. Then, it may reach R (1–10) km for t 10 s up to several hours (rather than for typical collapse time 10 3 s) with a delay owing to neutrino heat transport to the surface (effect of thermal conductivity) that strongly depends on the value of the pion softening (see in the next section), which is stronger for most massive neutron stars, cf. [23,24,25]. In addition, the dynamics of the pion–condensate transition is specific since it occurs in inhomogeneous state, k 0 . The seeds of the liquid-crystal-like state are elongated in the process of their growth [35], and such an effect is observed in liquid crystals.
It is well known that at, least in the mean-field approximation, the ITS line and AS line are different. From the exact thermodynamic relation, one gets
u s ˜ 2 = u T 2 + T n m c V P T n 2 .
The variable u s ˜ , u s ˜ 2 = m 1 ( P / n ) s ˜ has the meaning of the adiabatic sound velocity (at s ˜ = const) and u T , u T 2 = m 1 ( P / n ) T of the isothermal sound velocity (at T = const). The former quantity characterizes propagation of sound waves in ideal hydrodynamics. In non-ideal hydrodynamics at finite values of the thermal conductivity, κ , the propagation of sound waves is defined by the interplay between u T and u s ˜ . The conditions u T = 0 and u s ˜ = 0 define on a T ( n ) plane the ITS line and AS line. The maximum temperature points on these lines are the critical temperature T c r (on ITS line) and the adiabatic maximum temperature T P , m a x (on AS line). In the mean-field approximation, c V has finite non-negative values. Therefore, u T 2 0 on the AS line ( u s ˜ = 0 ). Note that, as we have mentioned, even in the region near the critical point, where thermodynamic fluctuations are strong (fluctuation region), a mean-field treatment can be applied in dynamics, since the system spends in this region a shorter time than the time needed for a development of the critical fluctuations. Calculations performed in mean-field models show that T c r is significantly higher than T P , m a x , e.g., within the Nambu-Jona-Lasinio (NJL) model, one obtains T P , m a x T c r / 2 45 MeV for the chiral transition. At such small temperatures as T P , m a x , the chiral and deconfinement transitions are unlikely in heavy-ion collisions.
To find onset of instabilities, consider dynamic variables a = ( n , s , T ) to be modulated as δ a = δ a 0 e x p [ γ t + i p r ] . For p p T , where p T is the mean thermal momentum, from the linearized equations of non-ideal hydrodynamics, one may find the increment γ ( p ) ,
γ 2 = p 2 u T 2 + ( 4 3 η + ζ ) γ m n + c p 2 + u s ˜ 2 u T 2 1 + κ p 2 ( c V γ ) .
The solutions are two density modes (existing even in the limit κ = 0 ) and one thermal mode (existing only for κ 0 ).
For the limiting case of zero shear and bulk viscosities and non-zero, but small thermal conductivity, for u T 2 1 , i.e., slightly below the ITS line, for the most rapidly growing mode γ = γ m , p = p m p T , one gets
γ m κ u T 4 4 c c V u s ˜ 2 , p m 2 u T 2 / ( 2 c ) .
For any small but finite value κ , the solution γ m results in instability already for u T 2 < 0 (i.e., below the ITS line). On the contrary, in case κ = η = ζ = 0 , i.e., within ideal hydrodynamics, the instability appears, when the system trajectory crosses the AS line rather than the ITS line. Since, in reality κ 0 , it is possible to conclude that the instability condition u s ˜ 2 < 0 should be replaced in favor of u T 2 < 0 . For large values of κ , the density mode proves to be the most rapidly growing one, but the condition for the onset of instability, u T 2 = 0 , is not changed. For β 1 (e.g. for a large viscosity), we get
γ m u T 2 m n c r 4 η / 3 + ζ , p m 2 u T 2 β / c .
Note that in condensed matter physics a transition from a liquid to a glass state can be interpreted as a first-order phase transition occurring within a spinodal region at a very high viscosity [35]. Then, there appears an order at several Å- scale, which transforms in a disorder at larger distances.
According to our findings, spinodal instabilities may manifest themselves in experiments with heavy ions in some collision energy interval that corresponds to the first-order phase transition region of the QCD phase diagram. One of the possible signatures is a manifestation of fluctuations with a typical size r 1 / p m in the rapidity spectra (see details in Ref. [30]). In conclusion, we note that viscosity and thermal conductivity are driving forces of the first-order liquid–vapor and quark-gluon phase transitions and the spinodal instability occurs for T below the ITS line.
The above considered are the first-order phase transitions of the liquid–vapor type when there exists only one conserved charge, as the baryon number in heavy-ion collisions. In BEM of neutron stars, there are two conserved charges, the baryon and the lepton ones. In this case, the equilibrium configuration is a mixed pasta phase [36], which appears provided a surface tension is below a critical value σ c , cf. [37]. Finite-size screening effects prove to be very important [38]. They were studied on examples of the quark–hadron first-order phase transition [38,39] and the kaon–condensate phase transition [40]. For the quark–hadron first-phase transition, σ c 60 MeV/fm 2 . It was shown that, in both cases, the P ( n ) behaviour is close to that given by the MC. In [41], the finite-size effects in the nuclear pasta phase in the inner neutron-star crust have been considered.

4. Pions and Kaons in Dense Baryon Matter

We treated nuclear systems on the mean-field level above. Pions and kaons are pseudoscalar Goldstone-like particles and they do not contribute to EoS on the mean-field level. For their description, medium-polarization effects are important. Different processes involved in nuclear phenomena occurring on the hadron level are characterized by various energy-momentum and time-length scales. Of the order of the confinement radius r Λ are scales associated with masses of heavy mesons such as σ , ω and ρ ( r r σ 1 / m σ , r r ω 1 / m ω , r r ρ 1 / m ρ ) and with nucleons, r N 1 / m N . In description of in-medium pions and kaons at n several n 0 , T 100 MeV, we treat these scales as short-ranged. Typical low momenta are k m π and k p F 1 / ( 2 m π ) , where m π 140 MeV is the pion mass and p F is the Fermi momentum of the nucleon, which corresponds to the length scales r 1 / m π and r 1 / p F , respectively. These latter scales are treated as long-ranged microscopic scales. Time scales that correspond to low energies are 1 / ϵ F ( ϵ F is the Fermi energy of nucleons 40 MeV in ISM for n = n 0 ), ∼1/ m π and ∼1/ ω Δ ( ω Δ = m Δ m N 2.1 m π is the energy gap for the excitation of Δ isobars). Minimal energy scale of our interest is ω ˜ ( n ) , where ω ˜ m π is an effective pion gap (see Figure 9 below). This energy scale arises due to a collective pion softening effect appearing for nucleon densities n > n c ( 1 ) (0.5–0.8 ) n 0 , cf. [25,42].
Consider the Fermi-liquid approach with the explicit incorporation of the in-medium pion exchange. For zero temperature, it was formulated by A. B. Migdal in Ref. [43], and, then, the approach was generalized for finite temperatures and non-equilibrium systems, cf. [25,42,44,45]. Within this approach, the long-range processes are treated explicitly, whereas short-range processes are described by local quantities approximated by phenomenological, so-called Landau–Migdal, parameters. At low excitation energies, the retarded N N interaction amplitude is presented as follows:
Universe 04 00028 i001
where
Universe 04 00028 i002
The solid line stands for a nucleon, the double-line, for a Δ isobar. At a higher temperature, the Δ Δ 1 diagrams should be included [44]. The doubly-dashed line corresponds to the exchange of the free pion with inclusion of the contributions of the residual s-wave π N N interaction and π π scattering, i.e., the residual irreducible interaction to the nucleon particle-hole and delta-nucleon hole insertions. The empty block in Equation (15) is irreducible with respect to the particle-hole, delta-nucleon hole and pion states and is thereby essentially more local than contributions given by explicitly presented graphs. It is given by
Γ n α β , n γ δ ω = F n , n δ α β δ γ δ + G n , n σ α β σ γ δ ,
where σ is the Pauli matrix, F and G are functions dependent on the direction of the momenta of incoming and outgoing nucleon and hole at the Fermi surface. One introduces the dimensionless amplitudes f n , n = C 0 1 F n , n and g n , n = C 0 1 G n , n where C 0 1 = m N ( n 0 ) p F ( n 0 ) / π 2 is the density of states at the Fermi surface for n = n 0 . The amplitudes f p , p and g p , p are then expanded in the Legendre polynomials with the coefficients known as Landau–Migdal parameters. In most cases, it is sufficient to deal with zero and first harmonics. These harmonics can be extracted from comparison with the data or can be computed within some models.
The irreducible part of the interaction involving Δ isobar is constructed similarly to Equation (15). The main part of the N Δ interaction is due to the pion exchange. Although information on the local part of the N Δ interaction is scarce, one concludes [25] that the corresponding Landau–Migdal parameters are essentially smaller then those for N N interaction. Therefore, and also for the sake of simplicity, we neglect this term.
The spectrum of the particle excitations is determined by the spectral function given by the imaginary part of the retarded Green function A = 2 Im G R . Resummation of diagrams shown in Equation (14) yields the following Dyson equation for the retarded pion Green function G R and self-energy Σ R ,
Universe 04 00028 i003
Here, Σ r e s R is a residual retarded pion self-energy that includes the contribution of all diagrams, which are not presented explicitly in Equation (17), like the s-wave π N and π π scattering (included by doubly-dashed line in Equation (15)). The full vertex takes into account N N correlations
Universe 04 00028 i004
Resummation of Equation (14) yields in the spin-isospin channel
Universe 04 00028 i005
The total resummed N N interaction reads
Γ N N R = C 0 [ F + F τ 1 τ 2 + ( G + G τ 1 τ 2 ) σ 1 σ 2 ] + f π N N 2 T π σ 1 k · σ 2 k ,
and f π N N = 1 / m π is the pion–nucleon P-wave coupling. First, three square-bracketed terms of Equation (20) arise due to the loop-resummation of the empty block graphs in Equation (15), whereas the last term of Equation (20) is due to the contribution of the second diagram of Equation (15). Amplitudes F, F , G, G , and T π are expressed with the help of the dimensionless Landau–Migdal parameters f, f , g, g , related to the local N N -interaction in the particle-hole channel
F = f Γ ( f ) , F = f Γ ( f ) , G = g Γ ( g ) , G = g Γ ( g ) , T π = Γ 2 ( g ) G π R ( ω , k ) = Γ 2 ( g ) / [ ω 2 m π 2 k 2 Σ R ( ω , k ) ] .
The Landau–Migdal parameters f = ( f n n + f n p ) / 2 , f = ( f n n f n p ) / 2 in the scalar channel are expressed via the nucleon incompressibility and can be calculated already on the mean-field level, cf. [10,46]. Calculation of the spin-parameters g = ( g n n + g n p ) / 2 , g = ( g n n g n p ) / 2 needs to incorporate the contribution in the N N interaction beyond the mean-field Hartree-level.
For small energies ω m π and for momenta k p F , expression for Γ ( x ) , where x = f , f , g , g , simplifies, since at such energies the Δ – nucleon hole ( Δ N 1 ) contribution in Equation (17) is a smooth function of ω and k, and, doing particle–hole resummation, one can explicitly take into account only the N N 1 loop. Then,
Γ ( x ) 1 / 1 + 2 x p F / p F ( n 0 ) .
There were found two sets of the parameters fitted to atomic-nucleus experiments, cf. [47], being f 0.25 , f 1 , g 0.5 , g 1 and f 0 , f 0.5–0.6, g 0.05 ± 0.1 , g 1.1 ± 0.1 . Uncertainties in numerical values appear due to attempts to get the best fit to experimental data in each specific case modifying the parametrization used for the residual part of the N N interaction. All of these numerical values of the parameters relate to ISM and the density n n 0 , whereas there is no direct experimental information on their values for n > n 0 and for asymmetric matter. There exist various calculations of the Landau–Migdal parameters as functions of the density for symmetric nuclear matter and for the neutron matter.
To simplify, we may present the spectral function of pions in ISM and for π 0 at arbitrary isospin composition as
A π ( ω , k ) i 2 π δ ω ω i ( k ) 2 ω Re Σ π R ( ω , k ) ω ω i ( k ) + 2 β 0 k ω ω ˜ 4 ( k ) + β 0 2 k 2 ω 2 θ ( ω < p F k / m N ) ,
and the sum is over the quasiparticle-like branches. The second term is due to virtual particle-hole modes. For n = n 0 in ISM, one estimates β 0 0.7 .
The quantity
G π R 1 ( μ π , k ) = ω ˜ 2 ( k ) = k 2 + m π 2 μ π 2 + Re Σ π ( μ π , k )
has the meaning of the squared effective pion gap with a typical behaviour as shown in Figure 9. At n > n c ( 1 ) (in ISM n c ( 1 ) 0.5–0.8 n 0 ), the quantity ω ˜ 2 develops a roton-like minimum ω ˜ 2 ( k ) ω ˜ 2 ( k 0 ) + γ 0 ( k 2 k 0 2 ) 2 / 4 k 0 2 with some density dependent parameters γ 0 1 , k 0 p F and the effective pion gap ω ˜ ( k 0 ) ω ˜ proves to be < m π . The effective pion gap determines the degree of the pion-mode softening. The matter for n > n c ( 1 ) can be treated as a liquid phase of a quantum pion condensate [35,44,48]. At normal nuclear density in ISM, one estimates ω ˜ 2 ( n 0 ) 0.8 . For n < n c ( 1 ) , there is no minimum at a finite momentum k and there is no pion softening effect thereby. In the BEM, the pion chemical potentials ( μ π + μ π 0 , μ π 0 = 0 ) are determined from equilibrium conditions for the reactions involving pions. In the neutron star matter, μ π follows from the condition of the chemical quasi-equilibrium with respect to the reactions n p π and n p e ν ¯ : μ π = μ e = ε F n ε F p , where ϵ F n , ϵ F p are the Fermi energies of the neutron and proton. For small-size systems like atomic nucleus, one should put μ π + = μ π = μ π 0 = 0 .
The spectrum of pion quasiparticles possesses three branches for π ± , 0 in the ISM and for π 0 in asymmetric matter ( N Z ). The spectrum computed for n = n 0 in ISM is shown in Figure 8. In the region ω m π , there are two branches: the Δ branch and the pion branch. For ω < m π , there is the spin-sound branch (with ω 0 when k 0 ). The hatches show the regions on the ( ω , k ) plane with a non-vanishing pion width, calculated within the quasiparticle approximation for nucleons and Δ isobars. The pion spectral function calculated beyond the quasiparticle approximation for nucleons and Δ s is more involved even for ω m π . In the lower hatched region, at ω < k p F and k p F , there are no quasiparticle branches at all and the pion width cannot be neglected. This is the region of the Landau damping in the nucleon particle-hole channel. The pion spectral function is enhanced in this region of ω and k for n > n c ( 1 ) ( 0.5 0.8 ) n 0 . The lowest-energy state determining by the pole of the pion Green function is i β 0 ω ω ˜ 2 ( k 0 ) with β 0 > 0 appeared due to the Landau damping. Thus, for ω ˜ 2 ( k 0 ) > 0 , the pion excitations die out with time exponentially exp ( ω ˜ 2 ( k 0 ) t / β 0 ) . With an increase of n, the quantity ω ˜ 2 ( k 0 ) decreases and may become negative for n > n c , π .
The spectrum of π + and π in isospin asymmetric matter is also reconstructed. For purely neutron matter, it can be found in [25,43].
For ω ˜ 2 ( k 0 ) < 0 , the classical pion field begins to grow exponentially during a passage of a time as exp ( | ω ˜ 2 ( k 0 ) | t / β 0 ) . Thus, the change in the sign of ω ˜ 2 at n = n c π marks the critical point of a phase transition to a liquid-crystal-like state or solid-like state (depending on the parameters of the matter under consideration, pion species and other conditions) with a classical pion field (a pion condensate). The value of the critical density n c π depends on the values of the Landau–Migdal parameters, which are badly known especially for asymmetric matter and for densities significantly larger than n 0 . Nevertheless, some estimations can be given. Various experiments have shown that the pion condensation does not manifest itself in atomic nuclei as a volume effect, see Ref. [25]. Different, although model-dependent, estimations indicate that n c π ( 1.5 3 ) n 0 , depending on the pion species, the proton-to-neutron ratio and the model used. For example, variational calculations [51] yield n c π 2 n 0 for ISM and n c π 1.3 n 0 for π 0 mesons in the neutron matter.
Assuming that strong interactions manage to keep nuclear system in local quasi-equilibrium up to its breakup and considering the fireball break up as prompt, for the pion distribution of the given species at infinity one gets [25,42,52,53]
d N π d 3 r d 3 k = 2 m π 2 + k 2 0 d ω ( 2 π ) 4 A ( ω , k , T ( t 0 , r ) , n ( t 0 , r ) ) e ω / T ( t 0 , r ) 1 .
This expression relates the resulting particle yield to the in-medium spectral density, A, of excitations, calculated at a freeze out density n ( t 0 , r ) and the temperature T ( t 0 , r ) . It is assumed that the state has no time to be changed, provided a typical time scale characterizing the fireball expansion τ exp ( n ( t 0 , r ) ) ( u / u ˙ ) | t 0 < 1 / | ω m π 2 + k 2 | . This condition is definitely fulfilled for the pion branch and partially for the Δ branch for k 1.5 m π , but not for virtual particle-hole modes with pion quantum numbers, which can be treated in an opposite, adiabatic approximation. Thus, one can drop contribution of the second term in Equation (23) related to particle-hole excitations with the pion quantum numbers, which have time to die out at the break up stage. Note that the pion momentum is conserved in this sudden freeze out model, but the pion energy is not conserved. The energy of the purely pion sub-system is not conserved. The conservation of the total energy is recovered with taking into account of the energy of the nucleon sub-system.
Exploiting this model for the Au+Au collision at energy ∼1 GeV per nucleon using the in-medium pion spectrum at freeze-out parameters n b , T b from [42], we obtain the pion yield displayed on the right panel in Figure 8, cf. [49]. A good overall agreement with experimental data for pion momenta 300 MeV < k < 700 MeV is found. Note that the smaller the collision energy, the more pronounced should be the effect of the pion–nucleon interaction on the pion differential cross section. On the contrary, with the increase of the collision energy, the pion softening effects become less pronounced. However, for the baryon-less matter and for T m π as occurs at RHIC and LHC energies, the in-medium polarization effects become pronounced again [54,55].
Typical density behavior of ω ˜ 2 ( k 0 ) (for π ± , π 0 in ISM and for π 0 at N Z ) is demonstrated in Figure 9.
For n < n c ( 1 ) , ω ˜ 2 ( k ) has a minimum for k = k 0 = 0 , i.e., ω ˜ 2 ( k 0 = 0 ) = m π 2 . Thereby, one may expect a stiffening of the N N interaction for n < n c ( 1 ) and, oppositely, the softening for n > n c ( 1 ) . At the critical point of the pion condensation, n = n c π , the value ω ˜ 2 ( k 0 ) , at artificially neglected π π fluctuations, changes the sign (see the dashed line in Figure 9).1 Without taking into account pion fluctuation effects, the pion condensate appears via a second-order phase transition. In reality, the π π fluctuations are significant in the vicinity of the critical point and the phase transition is of the first order [44,48,56]. Therefore, we depict branches in Figure 9 (solid curves) with positive and negative values of ω ˜ 2 ( k 0 ) . Calculations in Ref. [48] performed in the Thomas–Fermi approximation demonstrated that, at n > n c π , the free energy of the state with ω ˜ 2 ( k 0 ) > 0 and without the pion mean field becomes larger than the free energy of the state with ω ˜ 2 ( k 0 ) < 0 and a finite pion mean field. Therefore, at n > n c π , the state with ω ˜ 2 ( k 0 ) > 0 is metastable and the state with ω ˜ 2 ( k 0 ) < 0 and the pion mean field φ π 0 becomes the ground state. Note that, although qualitative behaviours of the effective pion gap for π and π 0 mesons in BEM are similar, quantitatively the corresponding values ω ˜ are different. Only simplifying, one may exploit the very same dependence shown in Figure 9.
The quantity ω ˜ 2 ( k 0 ) demonstrates how much the virtual (particle-hole) mode with pion quantum numbers is softened at the given density. Typical momentum transferred in two-nucleon reaction is k p F . For the symmetric nuclear matter at n = n 0 , the ratio α = G π 0 , 1 ( 0 , p F ) / G π 1 ( 0 , p F ) 6 , where G π 0 , 1 ( 0 , p F ) m π 2 + p F 2 , and G π 1 ( 0 , p F ) = ω ˜ 2 ( p F ) . We note that this pion softening [25,43] does not significantly enhance the N N scattering cross section for n = n 0 because of the simultaneous essential suppression of the π N N vertex by nucleon–nucleon correlations. Indeed, the ratio of the N N cross sections calculated with the free one-pion exchange (FOPE) and medium one-pion exchange (MOPE) models for N N interaction equals
R = σ [ FOPE ] σ [ MOPE ] Γ 4 ( ω 0 , k p F ) ( m π 2 + p F 2 ) 2 ω ˜ 4 ( p F ) ,
where Γ is the vertex dressing factor determined by Equations (18) and (22), and Γ ( n 0 ) 0.4 . For n n 0 , one has R 1 , whereas already for n = 2 n 0 this estimate yields R 10 . Thus, following Equation (26), one can evaluate the N N interaction for n > n 0 with the help of the MOPE, i.e., one approximately has
Universe 04 00028 i006
Here, the bold wavy line relates to the in-medium pion. In the soft-pion approximation, the same one-pion exchange determines also interaction in the particle–particle channel
Universe 04 00028 i007
For example, this quantity determines the N N interaction entering neutrino emissivities of the two-nucleon processes (see the next section).
Kaons as pions are also modified in nuclear matter. Here, the most important are hyperon— nucleon-hole diagrams [57,58] and the s-wave scattering kaon–nucleon interaction being stronger in the kaon case. On the left panel of Figure 10 a spectral density of K excitations in ISM at saturation is shown. On the right panel, results are presented of calculations of the invariant differential cross section of K production in Ni+Ni collisions with energy 1.8 GeV per nucleon as the function of the kaon kinetic energy in the center-of-mass system [59] in comparison with the experimental data [60,61]. For these calculations, one used the spectral function of K shown on the left panel and the prompt freeze-out model (see Equation (25) above). As above, the fireball dynamics were considered within the expanding fireball model [42]. We see a good overall agreement of calculations taking into account of in-medium polarization effects with the data.

5. Pion Softening and Neutrino Emission from Neutron Stars

After first tens of seconds (at most hours), the typical temperature of a neutron star decreases below the neutrino-opacity temperature T opac MeV. For T < T opac , the neutrino mean-free path becomes longer than the star radius and neutrinos are radiated directly from the star interior without subsequent rescattering [62,63,64,65,66]. Hence, the star can be considered for neutrinos as a warm “white” body. In this sense, the problem is somewhat similar to the problem of the di-lepton radiation from the nuclear fireball prepared in heavy-ion collisions. Typically averaged lepton energy (∼several T) is much larger than the nucleon particle width Γ N T 2 / ε F . Therefore, nucleons can be treated within the quasiparticle approximation. Then, the processes can be separated according to their phase space: the direct Urca (DU) one-nucleon processes, n p e ν ¯ (if they are not forbidden by the energy-momentum conservation) have the largest emissivity, ϵ ν T 6 for non-superfluid systems, the modified Urca (MU) two-nucleon processes n n n p e ν ¯ , n p p p e ν ¯ and nucleon bremstrahlung (NB) processes n n n n ν ν ¯ , n p n p ν ν ¯ , have a lesser emissivity, ϵ ν T 8 , and so on. The DU emissivity proves to be so high that if DU reactions occurred in a majority of neutron stars, they could not be visible in soft X-rays in disagreement with observations. The requirement of the absence of the DU processes in neutron stars with M < 1.5   M puts a constraint [1,8,67] on the density dependence of the symmetry energy, which is also an important quantity for the description of heavy-ion reactions. Within EoSs that we use (KVORcut-based and MKVOR-based), the one-nucleon DU processes are forbidden by energy-momentum conservation at least for light ( M 1   M ) and medium-heavy 1.5   M neutron stars and the most efficient reactions become processes involving two nucleons.
In the so-called standard scenario of the neutron star cooling [68,69], the processes were calculated without taking into account in-medium effects. The MU processes were considered as the most important channel for temperatures up to T 10 8 10 9 K. Then, based on the standard scenario, the so-called minimal cooling paradigm was proposed [70], where medium effects were assumed to play only a minor role. For evaluation of the MU and the NB neutrino emissivity within this paradigm, one uses the FOPE model of Friman and Maxwell [71]. The density dependence of the reaction rates calculated with the FOPE model is very weak and thereby the neutrino radiation from a neutron star depends very weakly on the star mass. It proves to be that the existing data on the time dependence of the surface temperatures of pulsars are then hardly explained. The agreement can be achieved only with inclusion of the efficient DU reaction in a majority of neutron stars and at the price that neutron stars measured in soft X-rays should then have approximately the same masses. The latter assumption disagrees with a broad distribution of neutron stars over the masses, as follows from the population synthesis modeling. Based on the assumption that the mass distribution of those objects for which surface temperatures are measured is similar to the one extracted, e.g., from a population synthesis analysis and from supernova simulations, the very efficient DU reaction should be forbidden in the majority of the former neutron stars. We add here that measurements of the high masses of the pulsars PSR J1614-2230 [72,73] and PSR J0348-0432 [4] and of the low masses for PSR J0737-3039B [74] for the companion of PSR J1756-2251 [75,76] and the companion of PSR J0453+1559 have provided the proof for the existence of neutron stars with masses varying at least from 1.2 to 1.97   M .
Refs. [25,62,63,64,65,66] calculated the rate of the two-nucleon processes taking into account in-medium effects. As we have argued, for n n 0 , the N N interaction amplitude might be mainly controlled by the soft pion exchange. Thus, instead of the FOPE-based MU diagram
Universe 04 00028 i008
(zig-zag line corresponds to the free pion, dot, to the free vertices), one has as a generalization a medium one-pion exchange (MOPE) for the medium-modified Urca (MMU) processes:
Universe 04 00028 i009
where the bold-wavy line corresponds to the dressed pion and the hatched vertex takes into account the N N correlations. The first diagram naturally generalizes the corresponding MU (FOPE) contribution (29). Evaluations have shown that it gives a smaller contribution in the emissivity for n n 0 than what comes from the processes occurring through intermediate reaction states, as shown by the second and third diagrams (30). Note that these diagrams are absent, if one approximates the nucleon–nucleon interaction by a two-body potential.
At T < T c (0.1–1), MeV nucleons form Cooper pairs. The emissivities of the processes with participation of nucleons are then suppressed by the corresponding phase-space suppression factors. However, a new nucleon pair-breaking-formation (PBF) reaction channel N N ν ν ¯ is opened up for T < T c , cf. [64,65,77]. This process was also included in calculations of the neutrino emissivity within the minimal cooling paradigm [70,78].
On the other hand, the nuclear medium cooling scenario was developed, where besides the MMU processes and the medium-modified nucleon bremstrahlung (MNB) processes [63], the PBF processes [64,65,79,80,81] and other ones were incorporated. The results of calculations were confronted to the data [67,82,83,84,85,86] demonstrating an overall agreement of calculations with the data already without necessity to include the DU channel. An example of such calculations is presented in Figure 11, where the cooling curves are shown being computed with MKVOR EoS [87]. The effective pion gap was used following the solid curve 1a + 1b in Figure 9. The pairing 1 S 0 gaps correspond to calculations with the EEHO model, cf. similar results of calculations with DD2 EoS in Figure 3 of [86].

6. Viscosity of Neutron-Star Matter and r-Modes in Rotating Pulsars

Neutron stars are typically formed rapidly rotating with a Kepler rotation frequency ν i n 1 kHz. However, the majority of young ( 10 5 y age) pulsars have rotation frequencies less than 10 Hz, and the fastest young pulsar PSR J0537-6910 (of age ∼5000 y) rotates with the frequency ν = ν m a x young = 62 Hz. An efficient mechanism allowing to decelerate rotation of a neutron star already at an early stage of its evolution is the r-mode instability [88,89]. The r-mode oscillations, related to Rossby waves in Earth atmosphere and oceans, lead to emission of gravitational waves, which carry away most of the initial angular momentum of a star and the star rotation decelerates. The r-modes would grow exponentially for any rotation frequency, if not a damping because of a viscosity of a warm neutron star matter [90]. The typical temperature in the interior of a pulsar of age ∼5000 y is ∼ 2 × 10 7 5 × 10 8 K, depending on the mass of the star and the cooling scenario, whereas the r-mode instability proves to be the strongest at higher temperatures T 10 9 K. Hence, young pulsars have passed through an instability phase during their early history and thus one should demonstrate that, after the instability phase, the pulsar frequencies remain larger than ν m a x young . Different mechanisms for suppression of the instability were studied. Most attempts were spent to find appropriate arguments to increase the values of shear and bulk viscosities, increasing thereby damping of modes. However, it was concluded, e.g., see [91], that the minimum value of the frequency at the r-mode stability boundary ν c ( T ) proves to be less than ν m a x young , if one uses standard dissipation mechanisms only not including the DU reactions.
Refs. [92,93] studied how the viscosity of the neutron-star matter influences the r-mode instability in rotating neutron stars provided a pion softening effect is included. For n > n c π , two possibilities were studied: one corresponding to saturation of the pion softening (see lines 1b and 1c in Figure 9) and the other including pion condensation (lines 3 and 2). In the latter case, the pion Urca (PU) processes, like p π c p e ν ¯ , n π c 0 p e ν ¯ , were included. The neutron–neutron and proton–proton Cooper pairing effects were also incorporated. In the shear viscosity, the lepton contribution was calculated with taking into account of the Landau damping in the photon exchange as was done in [94], the nucleon contribution described by the MOPE, and some other terms, such as the phonon contribution in the 1S 0 superfluid neutron phase, and the neutrino term in the neutrino opacity region of temperatures. It proves to be that the term in the bulk viscosity from the two-nucleon MMU reactions possesses the strongest density dependence (rising by several orders of magnitude for massive stars) because of the pion softening. In addition, contributions to the bulk viscosity arising from other reactions induced by charged weak currents, e.g., in the PU processes and DU processes, were included. The radiative bulk viscosity induced by charged and neutral weak currents in the region of the neutrino transparency of the star was also calculated while taking into account in-medium effects.
The characteristic time scale of the evolution of the r-mode amplitude is given by
τ 1 = τ G 1 + τ η 1 + τ ζ 1 ,
where τ G is the typical time of the gravitational radiation, τ η is the relaxation time due to the shear viscosity, and τ ζ is the relaxation time due to the bulk viscosity. The r-modes become stable, provided τ 1 > 0 .
We exploited the EoS similar to the Akmal–Pandharipande–Ravenhall EoS up to 4 n 0 , but stiffer at higher densities, producing the maximum neutron star mass compatible with observations. Within this EoS, the DU processes do not appear up to 5 n 0 (corresponding to the star mass M 1.9   M ). To compute MMU processes, we use curves 1a + 1b in Figure 9 and, for MMU + PU, ones we exploit curves 1a + 2 for MMU and 3 for PU.
Equation τ 1 = 0 was solved with all contributions to the shear and bulk viscosities taken into account. The solutions are presented in Figure 12 as functions of the temperature for neutron star masses from M < M < 2.05   M . Within the minimal cooling scenario, when in-medium effects are not included, ν c , min exceeds ν max young only for the masses M > 2.03   M (when DU reaction is already efficient), i.e., very close to the maximum mass, 2.05   M . However, if initially the star passes through the instability region, the developed r-modes may blow off some part of the star matter. Thus, its final mass can hardly be very close to the maximum mass. Alternatively, the experimental value of the frequency of the pulsar PSR J0537-6910 could be explained within the minimal cooling scenario, if one exploited the EoS that allows for a lower DU threshold density. Within the nuclear medium cooling scenario, the DU processes are not needed to explain the stability of PSR J0537-6910. We may explain it for M > 1.804   M , if calculations of MMU processes are done using curves 1a + 2, and 3 for PU, and for M > 1.776 provided we use curves 1a + 1c to calculate MMU, without PU. For the MMU case (for curves 1a + 1b for MMU without PU), PSR J0537-6910 mass should be M > 1.84   M . In the figure legend symbol, MMU + DU means that MMU and DU (for n > n c DU ) reactions are included and MMU + PU + DU means that PU processes (for n > n c π ) are included as well.

7. Conclusions

In the given contribution, I considered some specific but common problems in neutron-star and heavy-ion-collision physics, where in-medium effects play an important role. First, I discussed an effect of a scaling of the effective hadron masses and couplings on the equation of state of the baryon matter. Then, an important role was clarified, played by the shear and bulk viscosities and the thermal conductivity as driving forces of the first-order phase transitions. Next, the baryon-medium polarization effects on the pion and kaon propagation in dense matter were reviewed. Special attention was attended to the pion softening effect, and its particular role was indicated in the problem of neutrino radiation from neutron stars. Finally, viscosity effects, by taking into account the pion softening on the stability of r-modes in young pulsars, were analysed.

Acknowledgments

I thank David Blaschke, Hovik Grigorian, Evgeni E. Kolomeitsev and Konstantin A. Maslov for the help and advice. The research was supported by the Ministry of Education and Science of the Russian Federation within the state assignment, project No. 3.6062.2017/BY, by the Russian Foundation for Basic Research (RFBR) according to the research project No. 16-02-00023-A (Section 2.3, Section 2.4 and Section 4), and by the Russian Science Foundation, Grant No. 17-12-01427 (Section 3.2).

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

All necessary abbreviations are introduced in the paper body. For convenience some of them are repeated below:
EoSEquation of state
BEMBeta-equilibrium matter
ISMIsospin-symmetric matter
RMF modelRelativistic mean-field model
NLW modelNonlinear Walecka model
DU processDirect Urca process
MUmodified Urca
MMUmedium modified Urca
PUpion (condensate) Urca
NBnucleon bremstrahlung
MNBmedium-modified nucleon bremstrahlung
PBFPair-breaking-formation
FOPEFree one-pion exchange
MOPEMedium one-pion exchange

References

  1. Klähn, T.; Blaschke, D.; Typel, S.; Van Dalen, E.N.; Faessler, A.; Fuchs, C.; Gaitanos, T.; Grigorian, H.; Ho, A.; Kolomeitsev, E.E.; et al. Constraints on the high-density nuclear equation of state from the phenomenology of compact stars and heavy-ion collisions. Phys. Rev. C 2006, 74, 035802. [Google Scholar] [CrossRef]
  2. Danielewicz, P.; Lacey, R.; Lynch, W.G. Determination of the equation of state of dense matter. Science 2002, 298, 1592–1596. [Google Scholar] [CrossRef] [PubMed]
  3. Lynch, W.G.; Tsang, M.B.; Zhang, Y.; Danielewicz, P.; Famiano, M.; Li, Z.; Steiner, A.W. Probing the Symmetry Energy with Heavy Ions. Prog. Part. Nucl. Phys. 2009, 62, 427–432. [Google Scholar] [CrossRef]
  4. Antoniadis, J.; Freire, P.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] [PubMed]
  5. 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. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [PubMed]
  6. Khvorostukhin, A.S.; Toneev, V.D.; Voskresensky, D.N. Equation of State for Hot and Dense Matter: σ-ω-ρ Model with Scaled Hadron Masses and Couplings. Nucl. Phys. A 2007, 791, 180–221. [Google Scholar] [CrossRef]
  7. Khvorostukhin, A.S.; Toneev, V.D.; Voskresensky, D.N. Relativistic Mean-Field Model with Scaled Hadron Masses and Couplings. Nucl. Phys. A 2008, 813, 313–346. [Google Scholar] [CrossRef]
  8. Kolomeitsev, E.E.; Voskresensky, D.N. Relativistic mean-field models with effective hadron masses and coupling constants, and rho- condensation. Nucl. Phys. A 2005, 759, 373–413. [Google Scholar] [CrossRef]
  9. Maslov, K.A.; Kolomeitsev, E.E.; Voskresensky, D.N. Solution of the Hyperon Puzzle within a Relativistic Mean-Field Model. Phys. Lett. B 2015, 748, 369–375. [Google Scholar] [CrossRef] [Green Version]
  10. Maslov, K.A.; Kolomeitsev, E.E.; Voskresensky, D.N. Relativistic Mean-Field Models with Scaled Hadron Masses and Couplings: Hyperons and Maximum Neutron Star Mass. Nucl. Phys. A 2016, 950, 64–109. [Google Scholar] [CrossRef]
  11. Kolomeitsev, E.E.; Maslov, K.A.; Voskresensky, D.N. Delta isobars in relativistic mean-field models with σ-scaled hadron masses and couplings. Nucl. Phys. A 2016, 961, 106–141. [Google Scholar] [CrossRef]
  12. Maslov, K.A.; Kolomeitsev, E.E.; Voskresensky, D.N. Making a soft relativistic mean-field equation of state stiffer at high density. Phys. Rev. C 2016, 92, 052801. [Google Scholar] [CrossRef]
  13. Voskresensky, D.N. Comments on manifestation of in-medium effects in heavy-ion collisions. Eur. Phys. J. A 2016, 52, 223. [Google Scholar] [CrossRef]
  14. Kolomeitsev, E.E.; Maslov, K.A.; Voskresensky, D.N. Charged ρ-meson condensation in neutron stars. Nucl. Phys. A 2018, 970, 291–315. [Google Scholar] [CrossRef]
  15. Schaffner-Bielich, J. Hypernuclear Physics for Neutron Stars. Nucl. Phys. A 2008, 804, 309–321. [Google Scholar] [CrossRef]
  16. Djapo, H.; Schaefer, B.J.; Wambach, J. On the appearance of hyperons in neutron stars. Phys. Rev. C 2010, 81, 035803. [Google Scholar]
  17. Drago, A.; Lavagno, A.; Pagliara, G.; Pigato, D. Early appearance of Δ isobars in neutron stars. Phys. Rev. C 2014, 90, 065809. [Google Scholar] [CrossRef]
  18. Drago, A.; Lavagno, A.; Pagliara, G.; Pigato, D. The scenario of two families of compact stars: Part 1. Equations of state, mass–radius relations and binary systems. Eur. Phys. J. A 2016, 52, 41. [Google Scholar] [CrossRef]
  19. Riek, F.; Lutz, M.F.M.; Korpa, C.L. Photoabsorption off nuclei with self consistent vertex corrections. Phys. Rev. C 2008, 80, 024902. [Google Scholar] [CrossRef]
  20. Voskresensky, D.N. On the possibility of the condensation of the charged rho meson field in dense isospin asymmetric baryon matter. Phys. Lett. B 1997, 392, 262–266. [Google Scholar] [CrossRef]
  21. Chernodub, M.N. Electromagnetic superconductivity of vacuum induced by strong magnetic field. Lect. Notes Phys. 2013, 871, 143–180. [Google Scholar]
  22. Mallick, R.; Schramm, S.; Dexheimer, V.; Bhattacharyya, A. On the possibility of rho-meson condensation in neutron stars. Mon. Not. R. Astron. Soc. 2015, 449, 1347–1351. [Google Scholar] [CrossRef]
  23. Voskresensky, D.N.; Senatorov, A.V.; Kampfer, B.; Haubold, H.J. A possible explanation of the second neutrino burst in SN1987A. Astrophys. Space Sci. 1987, 138, 421–424. [Google Scholar] [CrossRef]
  24. Haubold, H.J.; Kampfer, B.; Senatorov, A.V.; Voskresensky, D.N. A Tentative approach to the second neutrino burst in SN1987A. Astron. Astrophys. 1988, 191, L22–L24. [Google Scholar]
  25. Migdal, A.B.; Saperstein, E.E.; Troitsky, M.A.; Voskresensky, D.N. Pion degrees of freedom in nuclear matter. Phys. Rept. 1990, 192, 179–437. [Google Scholar] [CrossRef]
  26. Galeotti, P.; Pizzella, G. New analysis for the correlation between gravitational wave and neutrino detectors during SN1987A. Eur. Phys. J. C 2016, 76, 426. [Google Scholar] [CrossRef]
  27. Prakash, M.; Bombaci, I.; Prakash, M.; Ellis, P.J.; Lattimer, J.M.; Knorren, R. Composition and structure of protoneutron stars. Phys. Rept. 1997, 280, 1–77. [Google Scholar] [CrossRef]
  28. Skokov, V.V.; Voskresensky, D.N. Hydrodynamical description of a hadron–quark first-order phase transition. JETP Lett. 2009, 90, 223–227. [Google Scholar] [CrossRef]
  29. Skokov, V.V.; Voskresensky, D.N. Hydrodynamical description of first-order phase transitions: Analytical treatment and numerical modeling. Nucl. Phys. A 2009, 828, 401–438. [Google Scholar] [CrossRef]
  30. Skokov, V.V.; Voskresensky, D.N. Thermal conductivity in dynamics of first-order phase transition. Nucl. Phys. A 2010, 847, 253–267. [Google Scholar] [CrossRef]
  31. Skokov, V.V.; Voskresensky, D.N. Viscosity and thermal conductivity effects at first-order phase transitions in heavy-ion collisions. Phys. Atom. Nucl. 2012, 75, 770–775. [Google Scholar]
  32. Randrup, J. Phase transition dynamics for baryon-dense matter. Phys. Rev. C 2009, 79, 054911. [Google Scholar] [CrossRef]
  33. Steinheimer, J.; Randrup, J. Spinodal density enhancements in simulations of relativistic nuclear collisions. Phys. Rev. C 2013, 87, 054903. [Google Scholar] [CrossRef]
  34. Romatschke, P.; Romatschke, U. Viscosity Information from Relativistic Nuclear Collisions: How Perfect is the Fluid Observed at RHIC? Phys. Rev. Lett. 2007, 99, 172301. [Google Scholar] [CrossRef] [PubMed]
  35. Voskresensky, D.N. Quasiclassical description of condensed systems by a complex order parameter. Phys. Scr. 1993, 47, 333–354. [Google Scholar] [CrossRef]
  36. Glendenning, N.K. First order phase transitions with more than one conserved charge: Consequences for neutron stars. Phys. Rev. D 1992, 46, 1274–1287. [Google Scholar] [CrossRef]
  37. Heiselberg, H.; Pethick, C.J.; Staubo, E.F. Quark matter droplets in neutron stars. Phys. Rev. Lett. 1993, 70, 1355–1359. [Google Scholar] [CrossRef] [PubMed]
  38. Voskresensky, D.N.; Yasuhira, M.; Tatsumi, T. Charge screening at first order phase transitions. Phys. Lett. B 2002, 541, 93–100. [Google Scholar] [CrossRef]
  39. Voskresensky, D.N.; Yasuhira, M.; Tatsumi, T. Charge screening at first order phase transitions and hadron quark mixed phase. Nucl. Phys. A 2003, 723, 291–339. [Google Scholar] [CrossRef]
  40. Maruyama, T.; Tatsumi, T.; Voskresensky, D.N.; Tanigawa, T.; Endo, T.; Chiba, S. Finite size effects on kaonic pasta structures. Phys. Rev. C 2006, 73, 035802. [Google Scholar] [CrossRef]
  41. Maruyama, T.; Tatsumi, T.; Voskresensky, D.N.; Tanigawa, T.; Chiba, S. Nuclear pasta structures and the charge screening effect. Phys. Rev. C 2005, 72, 015802. [Google Scholar] [CrossRef]
  42. Voskresensky, D.N. Many particle effects in nucleus-nucleus collisions. Nucl. Phys. A 1993, 555, 293–328. [Google Scholar] [CrossRef]
  43. Migdal, A.B. Pion Fields in Nuclear Matter. Rev. Mod. Phys. 1978, 50, 107–172. [Google Scholar] [CrossRef]
  44. Voskresensky, D.N.; Mishustin, I.N. Polarization operator of pions at finite temperatures. Yad. Fiz. 1982, 35, 1139–1156. [Google Scholar]
  45. Kolomeitsev, E.E.; Voskresensky, D.N. Superfluid nucleon matter in and out of equilibrium and weak interactions. Phys. Atom. Nucl. 2011, 74, 1316–1363. [Google Scholar] [CrossRef]
  46. Roepke, G.; Voskresensky, D.; Kryukov, I.; Blaschke, D. Fermi liquid, clustering, and structure factor in dilute warm nuclear matter. Nucl. Phys. A 2018, 970, 224–258. [Google Scholar] [CrossRef]
  47. Migdal, A.B. Theory of Finite Fermi Systems and Properties of Atomic Nuclei, Wiley and Sons: New York, NY, USA, 1967, 2nd ed.; Nauka: Moscow, Russia, 1983. [Google Scholar]
  48. Dyugaev, A.M. Crystalline And Liquid Phases Of A Pion Condensate. Pisma Zh. Eksp. Teor. Fiz. 1982, 35, 341–344. [Google Scholar]
  49. Kolomeitsev, E.E.; Voskresensky, D.N. Meson particle hole dynamics. arXiv, 2000; arXiv:nucl-th/0001062. [Google Scholar]
  50. Wagner, A.; Muntz, C.; Oeschler, H.; Sturm, C.; Barth, R.; Cieslak, M.; Debowski, M.; Grosse, E.; Koczon, P.; Mang, M.; et al. Evidence for different freezeout radii of high-energy and low-energy pions emitted in Au + Au collisions at 1-A/GeV. Phys. Lett. B 1998, 420, 20–24. [Google Scholar] [CrossRef]
  51. Akmal, A.; Pandharipande, V.R.; Ravenhall, D.G. The Equation of state of nucleon matter and neutron star structure. Phys. Rev. C 1998, 58, 1804–1828. [Google Scholar] [CrossRef]
  52. Senatorov, A.V.; Voskresensky, D.N. Pion Dynamics in Heavy Ion Collisions. Phys. Lett. B 1989, 219, 31–34. [Google Scholar] [CrossRef]
  53. Voskresensky, D.N.; Blaschke, D.; Roepke, G.; Schulz, H. Nonequilibrium approach to dense hadronic matter. Int. J. Mod. Phys. E 1995, 4, 1–45. [Google Scholar] [CrossRef]
  54. Voskresensky, D.N. Hadron liquid with a small baryon chemical potential at finite temperature. Nucl. Phys. A 2004, 744, 378–444. [Google Scholar] [CrossRef]
  55. Voskresensky, D.N. Thermodynamics of resonances and blurred particles. Nucl. Phys. A 2008, 812, 158–185. [Google Scholar] [CrossRef]
  56. Dyugaev, A.M. Nature of Phase Transition in Case of pi- Condensation. Pisma Zh. Eksp. Teor. Fiz. 1975, 22, 181–185. [Google Scholar]
  57. Kolomeitsev, E.E.; Kampfer, B.; Voskresensky, D.N. Kaon polarization in nuclear matter. Nucl. Phys. A 1995, 588, 889–917. [Google Scholar] [CrossRef]
  58. Kolomeitsev, E.E.; Voskresensky, D.N. Negative kaons in dense baryonic matter. Phys. Rev. C 2003, 68, 015803. [Google Scholar] [CrossRef]
  59. Kolomeitsev, E.E.; Kampfer, B.; Voskresensky, D.N. The impact of kaon polarization in nuclear matter on the K- production in heavy ion collisions. Int. J. Mod. Phys. E 1996, 5, 313–328. [Google Scholar] [CrossRef]
  60. Barth, R.; Senger, P.; Ahner, W.; Beckerle, P.; Bormann, C.; Brill, D.; Cieślak, M.; Grosse, E.; Koczoń, P.; Kohlmeyer, B.; et al. Subthreshold production of kaons and anti-kaons in nucleus-nucleus collisions at equivalent beam energies. Phys. Rev. Lett. 1997, 78, 4007–4010. [Google Scholar] [CrossRef]
  61. Schröter, A.; Berdermann, E.; Geissel, H.; Gillitzer, A.; Homolka, J.; Kienle, P.; Koenig, W.; Povh, B.; Schumacher, F.; Ströher, H. Subthreshold anti-proton and K- production in heavy ion collisions. Z. Phys. A Hadron. Nucl. 1994, 350, 101–113. [Google Scholar] [CrossRef]
  62. Voskresensky, D.N.; Senatorov, A.V. Pion Excitations In A Nucleonic Medium May Be Pertinent To The Luminosity Of Neutron Stars. Pisma Zh. Eksp. Teor. Fiz. 1984, 40, 395–398. [Google Scholar]
  63. Voskresensky, D.N.; Senatorov, A.V. Emission of Neutrinos by Neutron Stars. Zh. Eksp. Teor. Fiz. 1986, 90, 1505–1526. [Google Scholar]
  64. Senatorov, A.V.; Voskresensky, D.N. Collective Excitations in Nucleonic Matter and the Problem of Cooling of Neutron Stars. Phys. Lett. B 1987, 184, 119–124. [Google Scholar] [CrossRef]
  65. Voskresensky, D.N.; Senatorov, A.V. Description of Nuclear Interaction in Keldysh’s Diagram Technique and Neutrino Luminosity of Neutron Stars. Yad. Fiz. 1987, 45, 657–669. (In Russian) [Google Scholar]
  66. Voskresensky, D.N. Neutrino cooling of neutron stars: Medium effects. In Physics of Neutron Star Interiors; Springer: Berlin/Heidelberg, Germany, 2001; pp. 467–503. [Google Scholar]
  67. Blaschke, D.; Grigorian, H.; Voskresensky, D.N. Cooling of neutron stars: Hadronic model. Astron. Astrophys. 2004, 424, 979–992. [Google Scholar] [CrossRef]
  68. Tsuruta, S. Thermal properties and detectability of neutron stars. I. Cooling and heating of neutron stars. Phys. Rep. 1979, 56, 237–277. [Google Scholar] [CrossRef]
  69. Shapiro, S.; Teukolsky, S.A. Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects; Wiley: New York, NY, USA, 1983; Chapter 11. [Google Scholar]
  70. Page, D.; Lattimer, J.M.; Prakash, M.; Steiner, A.W. Minimal cooling of neutron stars: A New paradigm. Astrophys. J. Suppl. Ser. 2004, 155, 623–650. [Google Scholar] [CrossRef]
  71. Friman, B.L.; Maxwell, O.V. Neutron Star Neutrino Emissivities. Astrophys. J. 1979, 232, 541–557. [Google Scholar] [CrossRef]
  72. Demorest, P.; Pennucci, T.; Ransom, S.; Roberts, M.; Hessels, J. Shapiro Delay Measurement of A Two Solar Mass Neutron Star. Nature 2010, 467, 1081–1083. [Google Scholar] [CrossRef] [PubMed]
  73. Fonseca, E.; Pennucci, T.T.; Ellis, J.A.; Stairs, I.H.; Nice, D.J.; Ransom, S.M.; Demorest, P.B.; Arzoumanian, Z.; Crowter, K.; Dolch, T.; et al. The NANOGrav Nine-year Data Set: Mass and Geometric Measurements of Binary Millisecond Pulsars. Astrophys. J. 2016, 832, 167. [Google Scholar] [CrossRef]
  74. Kramer, M.; Stairs, I.H.; Manchester, R.N.; McLaughlin, M.A.; Lyne, A.G.; Ferdman, R.D.; Burgay, M.; Lorimer, D.R.; Possenti, A.; D’amico, N.; et al. Tests of general relativity from timing the double pulsar. Science 2006, 314, 97–102. [Google Scholar] [CrossRef] [PubMed]
  75. Faulkner, A.J.; Kramer, M.; Lyne, A.G.; Manchester, R.N.; McLaughlin, M.A.; Stairs, I.H.; Hobbs, G.; Possenti, A.; Lorimer, D.R.; D’Amico, N.; et al. PSR J1756-2251: A New relativistic double neutron star system. Astrophys. J. 2004, 618, L119–L122. [Google Scholar] [CrossRef]
  76. Ferdman, R.D.; Stairs, I.H.; Kramer, M.; Janssen, G.H.; Bassa, C.G.; Stappers, B.W.; Demorest, P.B.; Cognard, I.; Desvignes, G.; Theureau, G.; et al. PSR J1756-2251: A pulsar with a low-mass neutron star companion. Mon. Not. R. Astron. Soc. 2014, 443, 2183–2196. [Google Scholar] [CrossRef] [Green Version]
  77. Flowers, E.; Ruderman, M.; Sutherland, P. Neutrino pair emission from finite-temperature neutron superfluid and the cooling of young neutron stars. Astrophys. J. 1976, 205, 541–544. [Google Scholar] [CrossRef]
  78. Yakovlev, D.G.; Kaminker, A.D.; Gnedin, O.Y.; Haensel, P. Neutrino emission from neutron stars. Phys. Rep. 2001, 354, 1–55. [Google Scholar] [CrossRef]
  79. Kolomeitsev, E.E.; Voskresensky, D.N. Neutrino emission due to Cooper-pair recombination in neutron stars revisited. Phys. Rev. C 2008, 77, 065808. [Google Scholar] [CrossRef]
  80. Kolomeitsev, E.E.; Voskresensky, D.N. Neutral weak currents in nucleon superfluid Fermi liquids: Larkin-Migdal and Leggett approaches. Phys. Rev. C 2010, 81, 065801. [Google Scholar] [CrossRef]
  81. Kolomeitsev, E.E.; Voskresensky, D.N. Spin excitonic and diffusive modes in superfluid Fermi liquids. Phys. Rev. C 2011, 84, 068801. [Google Scholar] [CrossRef]
  82. Schaab, C.; Voskresensky, D.; Sedrakian, A.D.; Weber, F.; Weigel, M.K. Impact of medium effects on the cooling of nonsuperfluid and superfluid neutron stars. Astron. Astrophys. 1997, 321, 591–604. [Google Scholar]
  83. Grigorian, H.; Voskresensky, D.N. Medium effects in cooling of neutron stars and 3P(2) neutron gap. Astron. Astrophys. 2005, 444, 913–929. [Google Scholar] [CrossRef]
  84. Blaschke, D.; Grigorian, H.; Voskresensky, D.N.; Weber, F. On the Cooling of the Neutron Star in Cassiopeia A. Phys. Rev. C 2012, 85, 022802. [Google Scholar] [CrossRef]
  85. Blaschke, D.; Grigorian, H.; Voskresensky, D.N. Nuclear medium cooling scenario in the light of new Cas A cooling data and the 2M pulsar mass measurements. Phys. Rev. C 2013, 88, 065805. [Google Scholar]
  86. Grigorian, H.; Voskresensky, D.N.; Blaschke, D. Influence of the stiffness of the equation of state and in-medium effects on the cooling of compact stars. Eur. Phys. J. A 2016, 52, 1–19. [Google Scholar] [CrossRef]
  87. Grigorian, H.; Kolomeitsev, E.E.; Maslov, K.A.; Voskresensky, D.N. On cooling of neutron stars with stiff equation of state including hyperons. arXiv, 2017; arXiv:1801.00040. [Google Scholar]
  88. Andersson, N. A New class of unstable modes of rotating relativistic stars. Astrophys. J. 1998, 502, 708–713. [Google Scholar] [CrossRef]
  89. Friedman, J.L.; Morsink, S.M. Axial instability of rotating relativistic stars. Astrophys. J. 1998, 502, 714–720. [Google Scholar] [CrossRef]
  90. Lindblom, L.; Owen, B.J.; Morsink, S.M. Gravitational radiation instability in hot young neutron stars. Phys. Rev. Lett. 1998, 80, 4843–4846. [Google Scholar] [CrossRef]
  91. Andersson, N.; Kokkotas, K.D. The R mode instability in rotating neutron stars. Int. J. Mod. Phys. D 2001, 10, 381–441. [Google Scholar] [CrossRef] [Green Version]
  92. Kolomeitsev, E.E.; Voskresensky, D.N. Mechanism of r-mode stability in young rapidly rotating pulsars. Eur. Phys. J. A 2014, 50, 180. [Google Scholar] [CrossRef]
  93. Kolomeitsev, E.E.; Voskresensky, D.N. Viscosity of neutron star matter and r-modes in rotating pulsars. Phys. Rev. C 2015, 91, 025805. [Google Scholar] [CrossRef]
  94. 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]
1.
In order to avoid misunderstanding, by the pion condensation we understand the liquid-crystal-like or solid-like phases (of an inhomogeneous k 0 0 condensate) rather than a liquid phase mentioned above.
Figure 1. Figure is taken from [1]. (a) pressure region consistent with experimental flow data in ISM (dark shaded region). The light shaded region extrapolates this region to higher densities within an upper (UB) and lower (LB) borders; (b) mass–radius constraints from thermal radiation of the isolated neutron star RX J1856.5-3754 (grey hatched region) and from quasi-periodic oscillations (QPOs) in the Low-Mass X-Ray Binary (LMXB) 4U 0614+09 (green hatched area). The constraint from 4U 1636-536 (orange hatched region) 2.0 ± 0.1   M was withdrawn and should be replaced by a similar constraint 2.01 ± 0.04   M obtained later in [4] for PSR J0348+0432.
Figure 1. Figure is taken from [1]. (a) pressure region consistent with experimental flow data in ISM (dark shaded region). The light shaded region extrapolates this region to higher densities within an upper (UB) and lower (LB) borders; (b) mass–radius constraints from thermal radiation of the isolated neutron star RX J1856.5-3754 (grey hatched region) and from quasi-periodic oscillations (QPOs) in the Low-Mass X-Ray Binary (LMXB) 4U 0614+09 (green hatched area). The constraint from 4U 1636-536 (orange hatched region) 2.0 ± 0.1   M was withdrawn and should be replaced by a similar constraint 2.01 ± 0.04   M obtained later in [4] for PSR J0348+0432.
Universe 04 00028 g001
Figure 2. Figure from [12]. (a) the scalar field potential U ( f ) in the NLW model and U ˜ ( f ) in the NLWcut model for various values c σ , as functions of the scalar field parameter f. The symbols indicate the values of [ f ( n ) ; U ( f ( n ) ) , U ˜ ( f ( n ) ) ] for the density n of the ISM varying from n 0 to 10 n 0 with the 1 n 0 step; (b) the nucleon effective mass in the ISM as a function of the nucleon density for the same models.
Figure 2. Figure from [12]. (a) the scalar field potential U ( f ) in the NLW model and U ˜ ( f ) in the NLWcut model for various values c σ , as functions of the scalar field parameter f. The symbols indicate the values of [ f ( n ) ; U ( f ( n ) ) , U ˜ ( f ( n ) ) ] for the density n of the ISM varying from n 0 to 10 n 0 with the 1 n 0 step; (b) the nucleon effective mass in the ISM as a function of the nucleon density for the same models.
Universe 04 00028 g002
Figure 3. Figure from [12]. (a) total pressure P for the NLW and NLWcut models as a function of the nucleon density in the ISM for various values c σ ; the insertion demonstrates the binding energy per nucleon. Two shadowed areas show the constraints on the pressure from the particle flow and the kaon production data in heavy-ion collisions extracted in [2,3]. Bold dots show the extrapolation of the pressure consistent with the giant monopole resonance (GMR) description [3]; (b) the neutron-star mass and radius as functions of the central density for models under consideration. The hatched band denotes the uncertainty in the value of the maximum measured neutron-star mass 2.01 ± 0.04   M .
Figure 3. Figure from [12]. (a) total pressure P for the NLW and NLWcut models as a function of the nucleon density in the ISM for various values c σ ; the insertion demonstrates the binding energy per nucleon. Two shadowed areas show the constraints on the pressure from the particle flow and the kaon production data in heavy-ion collisions extracted in [2,3]. Bold dots show the extrapolation of the pressure consistent with the giant monopole resonance (GMR) description [3]; (b) the neutron-star mass and radius as functions of the central density for models under consideration. The hatched band denotes the uncertainty in the value of the maximum measured neutron-star mass 2.01 ± 0.04   M .
Universe 04 00028 g003
Figure 4. Neutron star masses versus radii for MKVOR* Δ , MKVOR*H Δ ϕ , MKVOR* Δ ρ , MKVOR*H Δ ϕ ρ models.
Figure 4. Neutron star masses versus radii for MKVOR* Δ , MKVOR*H Δ ϕ , MKVOR* Δ ρ , MKVOR*H Δ ϕ ρ models.
Universe 04 00028 g004
Figure 5. Schematic pressure isotherms as functions of the density at a liquid–vapor-like phase transition. P c , n c and T c are the pressure, density and temperature at the critical point.
Figure 5. Schematic pressure isotherms as functions of the density at a liquid–vapor-like phase transition. P c , n c and T c are the pressure, density and temperature at the critical point.
Universe 04 00028 g005
Figure 6. Figure from [11]. Pressure as a function of the density in the ISM for MKVOR* model (dashed line) and MKVOR* Δ model (dotted and solid lines). The hatched region indicates the nucleon flow constraint [2] in heavy-ion collisions. MC lines indicate the Maxwell construction.
Figure 6. Figure from [11]. Pressure as a function of the density in the ISM for MKVOR* model (dashed line) and MKVOR* Δ model (dotted and solid lines). The hatched region indicates the nucleon flow constraint [2] in heavy-ion collisions. MC lines indicate the Maxwell construction.
Universe 04 00028 g006
Figure 7. Figure from [31]. The phase diagram of the van der Waals EoS, T ( n ) -plane. The bold solid, dashed and dash-dotted curves demonstrate the boundaries of the MC, the spinodal region at T = const and s ˜ = const, respectively. The short dashed lines show adiabatic trajectories of the system evolution: the curve labeled s c r passes through the critical point; s m , through the maximum pressure point P ( n P , m a x ) on the P ( n ) plane.
Figure 7. Figure from [31]. The phase diagram of the van der Waals EoS, T ( n ) -plane. The bold solid, dashed and dash-dotted curves demonstrate the boundaries of the MC, the spinodal region at T = const and s ˜ = const, respectively. The short dashed lines show adiabatic trajectories of the system evolution: the curve labeled s c r passes through the critical point; s m , through the maximum pressure point P ( n P , m a x ) on the P ( n ) plane.
Universe 04 00028 g007
Figure 8. Figure from [49]. (a) pion spectrum in ISM at the nuclear saturation density. Dash-dotted lines are the vacuum spectra of pions and Δ . Solid lines show quasiparticle branches of pion excitations in medium provided one neglects effects of the particle widths. The contour plot depicts the spectral density of virtual pions; (b) invariant differential cross section of pion production in Au+Au collisions with energy 1 GeV per nucleon as a function of the pion kinetic energy in the center-of-mass system in comparison with experimental data [50] (squares stand for π + mesons and circles for π ).
Figure 8. Figure from [49]. (a) pion spectrum in ISM at the nuclear saturation density. Dash-dotted lines are the vacuum spectra of pions and Δ . Solid lines show quasiparticle branches of pion excitations in medium provided one neglects effects of the particle widths. The contour plot depicts the spectral density of virtual pions; (b) invariant differential cross section of pion production in Au+Au collisions with energy 1 GeV per nucleon as a function of the pion kinetic energy in the center-of-mass system in comparison with experimental data [50] (squares stand for π + mesons and circles for π ).
Universe 04 00028 g008
Figure 9. Typical density dependence of the effective pion gap squared. The line 1a shows ω ˜ 2 ( k 0 , n ) for n c ( 1 ) < n < n c π . For n > n c π , the stable state is shown by the line 3. Line 2 demonstrates the quantity ω ˜ 2 ( k 0 , n ) that characterizes a softening of the pion spectrum in the presence of the pion condensate. Line 1c shows a metastable phase in the system where the ground state contains the pion condensate and line 1b demonstrates a possible saturation of the pion softening effect that could be, if the Landau–Migdal parameter g , or g, may increase with an increase of the density.
Figure 9. Typical density dependence of the effective pion gap squared. The line 1a shows ω ˜ 2 ( k 0 , n ) for n c ( 1 ) < n < n c π . For n > n c π , the stable state is shown by the line 3. Line 2 demonstrates the quantity ω ˜ 2 ( k 0 , n ) that characterizes a softening of the pion spectrum in the presence of the pion condensate. Line 1c shows a metastable phase in the system where the ground state contains the pion condensate and line 1b demonstrates a possible saturation of the pion softening effect that could be, if the Landau–Migdal parameter g , or g, may increase with an increase of the density.
Universe 04 00028 g009
Figure 10. Figure from [49]. (a) spectral density of K excitations in ISM at saturation. The upper curve shows the position of the quasiparticle kaon branch. Dashed curves border the ( Λ p 1 ) continuum. Thin lines between them depict the ascending levels of kaon spectral density; (b) invariant differential cross section of K production in Ni+Ni collisions with energy 1.8 GeV per nucleon as the function of the kaon kinetic energy in the center-of-mass system in comparison with experimental data [60,61]. The solid line depicts calculations with the in-medium spectral density, and the dashed line shows the results for the free kaon spectrum.
Figure 10. Figure from [49]. (a) spectral density of K excitations in ISM at saturation. The upper curve shows the position of the quasiparticle kaon branch. Dashed curves border the ( Λ p 1 ) continuum. Thin lines between them depict the ascending levels of kaon spectral density; (b) invariant differential cross section of K production in Ni+Ni collisions with energy 1.8 GeV per nucleon as the function of the kaon kinetic energy in the center-of-mass system in comparison with experimental data [60,61]. The solid line depicts calculations with the in-medium spectral density, and the dashed line shows the results for the free kaon spectrum.
Universe 04 00028 g010
Figure 11. Figure from H. Grigorian calculations with MKVOR EoS [87]. T s is the redshifted surface temperature, t is the neutron star age. The effective pion gap is given by the solid curve 1a + 1b in Figure 9. The pairing 1 S 0 gaps correspond to calculations with the EEHO model. The mass range is shown in the legend.
Figure 11. Figure from H. Grigorian calculations with MKVOR EoS [87]. T s is the redshifted surface temperature, t is the neutron star age. The effective pion gap is given by the solid curve 1a + 1b in Figure 9. The pairing 1 S 0 gaps correspond to calculations with the EEHO model. The mass range is shown in the legend.
Universe 04 00028 g011
Figure 12. Figure from [93]. Critical frequency of the r-mode instability as a function of the temperature for a band of neutron star masses. Dashed line shows calculation with MU + DU for M = 2.05   M . The solid line demonstrates calculation done for the threshold mass of the DU processes within MMU + PU + DU scenario. The dashed-dotted line shows calculation with MMU computed with 1a + 1c curve in Figure 9 for M = 2.05   M .
Figure 12. Figure from [93]. Critical frequency of the r-mode instability as a function of the temperature for a band of neutron star masses. Dashed line shows calculation with MU + DU for M = 2.05   M . The solid line demonstrates calculation done for the threshold mass of the DU processes within MMU + PU + DU scenario. The dashed-dotted line shows calculation with MMU computed with 1a + 1c curve in Figure 9 for M = 2.05   M .
Universe 04 00028 g012

Share and Cite

MDPI and ACS Style

Voskresensky, D.N. On Manifestation of In-Medium Effects in Neutron Stars and Heavy-Ion Collisions. Universe 2018, 4, 28. https://doi.org/10.3390/universe4020028

AMA Style

Voskresensky DN. On Manifestation of In-Medium Effects in Neutron Stars and Heavy-Ion Collisions. Universe. 2018; 4(2):28. https://doi.org/10.3390/universe4020028

Chicago/Turabian Style

Voskresensky, Dmitry N. 2018. "On Manifestation of In-Medium Effects in Neutron Stars and Heavy-Ion Collisions" Universe 4, no. 2: 28. https://doi.org/10.3390/universe4020028

APA Style

Voskresensky, D. N. (2018). On Manifestation of In-Medium Effects in Neutron Stars and Heavy-Ion Collisions. Universe, 4(2), 28. https://doi.org/10.3390/universe4020028

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