arxiv: v2 [nucl-th] 21 May 2012

Transcription

1 Constraining mean-fied modes of the nucear matter equation of state at ow densities arxiv: v2 [nuc-th] 21 May 2012 M.D. Voskresenskaya, S. ype GSI Hemhotzzentrum für Schwerionenforschung GmbH, heorie, Panckstraße 1, D Darmstadt, Germany Abstract An extension of the generaized reativistic mean-fied (grmf mode with density dependent coupings is introduced in order to describe thermodynamica properties and the composition of dense nucear matter for astrophysica appications. Bound states of ight nucei and two-nuceon scattering correations are considered as expicit degrees of freedom in the thermodynamica potentia. hey are represented by quasipartices with medium-dependent properties. he mode describes the correct ow-density imit given by the viria equation of state (VEoS and reproduces RMF resuts around nucear saturation density where custers are dissoved. A comparison between the fugacity expansions of the VEoS and the grmf mode provides consistency reations between the quasipartices properties, the nuceon-nuceon scattering phase shifts and the meson-nuceon coupings of the grmf mode at zero density. Reativistic effects are found to be important at temperatures that are typica in astrophysica appications. Neutron matter and symmetric nucear matter are studied in detai. Keywords: mean-fied modes, nucear matter, viria equation of state, correations, custer formation, custer dissoution, neutron matter Preprint submitted to Nucear Physics A December 21, 2013

2 1. Introduction Correations are an essentia feature of an interacting many-body system such as nucear matter. hey have an impact on the thermodynamica properties and the composition of the matter since the constituent partices can form new nucear species as bound states. his aspect is particuary important in the description of matter beow nucear saturation density that is empoyed in theoretica simuations of astrophysica phenomena such as corecoapse supernovae. here, the knowedge of the equation of state (EoS of neutron-proton asymmetric matter is required in a wide range of densities and temperatures [1, 2]. he description of dense matter is important for the investigation of various stages in supernova exposions and the structure of neutron stars [3]. It is known to affect the effectiveness of the neutrino reheating of the shock wave [4]. he sensitivity of the coapse dynamics on the properties of matter and the composition in this density regime coud strongy infuence the structure of the proto-neutron star [5]. Athough there are many approaches to describe dense matter, in particuar at zero temperature, ony few modes cover the fu parameter space needed in most astrophysica appications. Often, they do not suppy sufficient information on the thermodynamica and compositiona detais. hus, for many years, a very sma number of EoS tabes was avaiabe which have been used in simuations of dynamica astrophysica processes [6 8]. However, in recent years, the interest in deveoping EoS has surged and severa new EoS tabes were provided [9 16]. his deveopment was triggered by the needs of astrophysica modeers, the progress in the theoretica description of nucear matter and an increase of computationa power. Nevertheess, approximations and simpifications are sti needed for practica purposes. here are two major paths to buid an EoS of warm and dense nucear matter for practica appications. One approach starts with an idea mixture of nuceons and nucei eading to a nucear statistica equiibrium (NSE description [11]. he effect of interactions between a constituents can be incorporated with the hep of viria corrections. However, at present, ony nuceons and ight nucei are considered in practice in such a viria equation of state (VEoS [17 19] that provides the correct finite-temperature EoS in the imit of very ow densities. he resuts are mode independent since they depend ony on experimentay determined data, i.e. binding energies of nucei and scattering phase shifts. Unfortunatey, the appication of this approach is imited to rather ow densities. he dissoution of nucei and the transition to 2

3 uniform neutron-proton matter with increasing density cannot be described propery. In order to simuate such an effect, the heuristic excuded voume mechanism was empoyed frequenty [11]. A second cass of EoS modes for astrophysica appications is based on sef-consistent mean-fied methods with neutrons and protons as fundamenta constituents. hey are considered as quasipartices with sef-energies that contain the information on the interaction which is usuay modeed in an effective way and not taken from a reaistic nuceon-nuceon (NN interaction. hese mean-fied modes, both non-reativistic and reativistic, can be very successfu in describing finite nucei and nucear matter around saturation density n sat 0.16 fm [20 23]. But at ow densities n n sat they fai to incude correations propery that give rise to inhomogeneities and the formation of many-body bound states, i.e. nucei. here were attempts to modify the ow-density behavior of meanfied modes guided by microscopic cacuations, see, e.g., Ref. [24] for the zero temperature case. However, it is sti chaenging to incude the formation of custers propery. he deficiencies of the various modes ead to the strategy of patching up or merging different approaches providing a uniform description of warm dense matter. A recent pubication[25] considered a quantum statistica approach to nucear matter and devised a generaized reativistic mean-fied (grmf mode, which is an extension of the reativistic mean-fied mode with density dependent coupings (DD-RMF from Ref. [26]. he grmf mode incudes, besides nuceons, ight nucei ( 2 H= d, 3 H= t, 3 He= h, 4 He= α as hadronic degrees of freedom in the Lagrangian. he dissoution of these custers is modeed by a medium-dependent shift of their binding energies originating mainy from the action of the Paui principe. A hadronic constituents of the mode can be viewed as quasipartices that interact via the exchange of effective mesons. he grmf mode achieves bridging the two views on nucear matter, i.e. a mixture of nuceon and nucei at ow densities and nuceons as quasipartices moving in mean fieds at high densities. Even though the mode contains the same reevant partices at ow densities as the VEoS, it does not reproduce exacty the thermodynamica properties of the VEoS as wi be shown beow. In fact, a presenty empoyed EoS for astrophysica appications that are based on mean-fied concepts share the probem of reproducing the VEoS in the ow-density imit. hus the question arises, how the modes can be modified in order to approach the correct ow-density imit. Is a modification of the effective interaction sufficient to meet this aim or are there more severe changes of the modes required? 3

4 In this work we propose an approach that introduces bound and scattering states in an effective way as additiona degrees of freedom in the thermodynamica potentia. A reevant quantities are derived in a thermodynamicay consistent way. hese many-body bound and scattering states are represented by quasipartices with medium-dependent properties and temperature dependent energies. he ow-density behavior of nucear matter at finite temperatures is considered in detai by comparing the VEoS with the grmf approach by means of a series expansion of the grand canonica potentia in powers of the nuceon fugacities. From the comparison of the VEoS and grmf expansions we wi derive consistency reations that connect quasipartice properties with the meson-nuceon coupings in the vacuum and the phase shifts or effective-range parameters of nuceon-nuceon scattering. We investigate different choices of the meson-nuceon coupings and quasipartice properties set by these consistency reations. In addition, reativistic corrections to the traditiona VEoS are obtained from these reations that become arger than the effects of partice correations in the ow-density imit. Various parametrizations of the quasipartice medium dependence are discussed. Studying thermodynamica properties of the matter we consider temperatures higher than the critica temperatures for pairing and Bose condensation in order to avoid these more compex effects. hus our consideration of the ow and zero temperature imits shoud be corrected in the future to incorporate condensation phenomena. Our approach is not imited to the chosen particuar grmf mode with density-dependent mesonnuceon coupings. It can be appied easiy to other mean-fied approaches as we. his paper is organized as foows. In Section 2 we present a genera formuation of the VEoS that serves as one starting point of the discussion. Basic quantities are defined and the notation is estabished. We show that correations in continuum states can be represented effectivey by temperature dependent resonances. Anaytic expressions for the second viria coefficients are found with the hep of the effective-range expansion for the s-wave phase shifts. he connection of the VEoS and NSE modes is estabished. he reation between the VEoS and the generaized Beth-Uhenbeck approach, which takes into account modifications due to medium effects, is discussed. he grmf mode with density-dependent meson-nuceon coupings is sketched in Section 3 and a series expansion of the grand canonica potentia Ω in powers of fugacities is derived for ow densities. he power series for Ω in the VEoS and grmf are compared up to second-order eading 4

5 to consistency reations between the modes that are studied in various imits. In Sect. 4 we consider an extension of the grmf mode: in order to satisfy consistency reations we introduce temperature dependent resonance energies of the scattering states and modified degeneracy factors that aow to make a smooth interpoation between VEoS and grmf modes. We find that reativistic corrections to the standard VEoS are important in these consistency reations aready in the first order of the expansion. he particuar exampe of neutron matter is presented in Section 5 and symmetric nucear matter is discussed in Section 6. he foowing section deas with the transition from ow to high densities and the occurring probems. Concuding remarks and an outook are given in Section 8. he comparison of our notations with the ones given in Ref. [18] can be found in Appendix A. An expansion of the energy per partice in neutron matter at zero temperature in powers of the Fermi momentum is considered in Appendix B. hroughout this work we use natura units such that = c = Equation of state in the viria imit 2.1. Genera formaism he VEoS represents a mode-independent approach in the cacuation of thermodynamica properties of ow-density matter at finite temperatures. In the non-reativistic imit, which is usuay considered, it describes a system of interacting partices i,j,... with non-reativistic chemica potentias µ i in a voume V, provided the fugacities z i = exp(µ i / are sma (z i 1. he range of densities where the viria expansion is vaid can be estimated by the reation n i λ 3 i 1 where n i is the number density of partice i with mass m i and λ i = 2π/(m i is the therma waveength. he foowing presentation gives a generaized and more symmetric formuation of the viria approach as compared to [17 19]. he S-matrix formuation by Ref. [27] that was recenty appied to nucear matter in Ref. [28] gives essentiay identica resuts. Under the presupposed conditions the grand canonica partition function Q can be expanded in powers of fugacities as Q(,V,µ i = 1+ Q i z i + 1 Q ij z i z j + 1 Q ijk z i z j z k +... (1 2 6 i ij with one-, two-, three-,...many-body canonica partition functions Q i, Q ij, Q ijk,... In cassica non-reativistic mechanics, the singe-partice canonica 5 ijk

6 partition function is Q i = g i (2π 3 d 3 r i d 3 p i exp( βh i = g i V (2 λ 3 i with the spin degeneracy factor g i (= 2 for each nuceon, β = 1/ and the singe-partice Hamitonian H i = p 2 i/(2m i that ony contains a kinetic contribution with momentum p i. For the two-body canonica partition function we have Q ij = g ig j d 3 r (2π 6 i d 3 p i d 3 r j d 3 p j exp( βh ij (3 withthetwo-bodyhamitonianh ij = p 2 i /(2m i+p 2 j /(2m j+v ij thatincudes a two-body potentia V ij responsibe for the correations. hen, the grand canonica potentia Ω(,V,µ i = nq(,v,µ i = pv, (4 that is directy reated to the pressure p, can be written in the form Ω(,V,µ i = (5 ( b i V z λ 3 i + b ij z i i ij λ 3/2 i λ 3/2 i z j + b ijk z i z j z k +... λ j i λ j λ k ijk with the (symmetrized dimensioness custer (or viria coefficients b i ( = λ3 i V Q i, (6 b ij ( = λ3/2 i λ 3/2 j (Q ij Q i Q j, (7 2V b ijk ( = λ iλ j λ k 6V (Q ijk Q i Q jk Q j Q ik Q k Q ij +2Q i Q j Q k (8 that depend on the temperature. Without interaction, the two-, three-,...many-body partition functions factorize, e.g. Q ij = Q i Q j, and the second, third,...custer coefficients b ij, b ijk,...vanish. Individua partice number densities are found from the reation n i = 1 V Ω z i µ i = b i +2,V,µj i λ 3 i j b ij 6 z i z j λ 3/2 i λ 3/2 j +3 jk b ijk z i z j z k λ i λ j λ k +... (9

7 with contributions from free partices (first term in the second equation and correated pairs, tripes, etc. he entropy is S = Ω (10 V,µi ( = 5Ω 2 +V c i z λ 3 i + c ij z i i ij λ 3/2 i λ 3/2 i z j + c ijk z i z j z k +... λ j i λ j λ k ijk with coefficients c i = db i d µ i, (11 c ij = db ij d µ i +µ j, (12 c ijk = db ijk d µ i +µ j +µ k. (13 Other reevant thermodynamica quantities such as the free energy and the interna energy F = Ω+V µ i n i (14 i and the interna energy E = F +S (15 can be obtained immediatey. heformuasgivenabovecanbegeneraizedtoincudefour-,five-,...manybody correations, but aready contributions of the three-body term are hardy ever considered in practice. For a given temperature, two-body correations wi aways dominate higher-order correations when the density decreases. Hence, in the foowing, we wi truncate the expansion at second order in the fugacities of the basic constituents. he quantum mechanica generaization of the viria expansion up to second order was given by Beth and Uhenbeck [29, 30]. In cassica mechanics the interaction potentia between the partices is the reevant quantity that appears in the cacuation of the viria coefficients. In quantum mechanics, the Schrödinger equation has to be soved with this potentia and the obtained eigenstates comprise bound and scattering states. Integrations over phase space are repaced by sums over a eigenstates of the one- and twobody system and it is obvious that the density of states becomes the reevant 7

8 quantity. heresut forthefirst viriacoefficient b i isidentica tothecassica vaue (6 in the continuum approximation, corresponding to the repacement of the sum over discrete momentum states by an integra. he second viria coefficient can be expressed as an integra over center-of-mass energies E b ij ( = 1+δ ij 2 λ 3/2 i λ 3/2 j λ 3 ij with λ ij = 2π/[(m i +m j ] and the quantity de exp( βed ij (E±δ ij g i 2 5/2 (16 D ij (E = k g (ij k δ(e E (ij k + g (ij π dδ (ij de (17 that measures the difference between the density of states of an interacting and free two-partice system. he ast term in equation (16 is a quantum statistica correction with the positive (negative sign if i = j are identica bosons (fermions. he first contribution to D ij is a sum over a two-body bound states k of thesystem ij withdegeneracy factorsg (ij k andenergiese (ij k < 0. hesecond term takes into account continuum correations in the two-partice system in via the energy-dependent scatter-. In genera, the index stands for the tota and anguar orbita momenta of a particuar channe. If the phase shift is dominated by a narrow resonance, it jumps by π in a very sma energy interva around the resonance energy and the contribution to the viria coefficient resembes that of a bound state. With experimentay known bound state energies and phase shifts, the second viria coefficient b ij and thus the ow-density EoS can be cacuated in a mode-independent way. he quantum mechanica generaization of the third viria coefficient was discussed in Ref. [31] but actua cacuations turn out to be difficut in practice and, thus, they are not considered here. he derivation of the cassica VEoS is based on non-reativistic kinematics. In order to incude reativistic effects, at east the reativistic dispersion reation of the partices has to be used. he resuting modifications of the viria coefficients wi be discussed when the VEoS is compared to the grmf approach in Section 3.4. a channes with degeneracy factors g (ij ing phase shifts δ (ij 2.2. Appication to nucear matter with arbitrary neutron to proton ratio In ow-density nucear matter at not too high temperatures, neutrons and protons are the reevant fundamenta constituents. Hence, the grand 8

9 canonica potentia up to second order in the fugacities becomes Ω(,V,µ n,µ p (18 ( b n = V z λ 3 n + b p z n λ 3 p + b nn z 2 p λ 3 n + b pp z 2 n λ 3 p +2 b np z p λ 3/2 n λ 3/2 n z p p where the symmetry b np = b pn is used. he Couomb interaction is usuay negected in the VEoS and b pp is repaced by b nn in (18, cf. Ref. [18]. he second viria coefficient b nn receives contributions ony fromscattering states. It is convenient to spit b np = b np0 + b np1 into two contributions with tota isospin = 0 and = 1. he ony bound two-body state in the two-nuceon system, the deuteron, appears in the isospin zero channe with g d = 3 and E (np d = B d = MeV [32]. According to this observation we can write b nn ( = λ3 n λ 3 nn b pp ( = λ3 p λ 3 pp b np1 ( = 1 2 b np0 ( = 1 2 with the viria integra λ 3/2 n λ 3/2 p λ 3 np λ 3/2 n λ 3/2 p λ 3 np I (ij ( = g (nn I (nn g n 2 5/2, (19 g (pp I (pp g p 2 5/2, (20 0 [ g (np1 I (np1, (21 ( Bd g d exp + ] g (np0 I (np0 (22 de π dδ (ij ( de exp E. (23 Formay, we can write the sum of the scattering contributions in Eqs. (19- (22 in the form of a singe bound state contribution ( g (ij I (ij de dδ ij = π de exp E [ = ĝ ij exp E ] ij( (24 0 with a temperature dependent effective resonance energy E ij ( in each ij channe by summing the contributions of a partia waves in an effective tota phase shift g (ij δ (ij. (25 δ ij (E = 9

10 he effective degeneracy factor ĝ ij = ±g (ij 0 is chosen such that it is identica to the degeneracy factor in the s-wave. he sign is determined by that of the integra. In order to iustrate the behavior of the phase shifts and effective resonance energies, neutron-neutron scattering is considered in the foowing. In this exampe, the sum of the experimenta phase shifts receives contributions ony from certain partia waves with tota anguar momentum J and orbita anguar momenta due to the Paui principe. aking into account a possibe partia waves with = 0, 1 and 2, i.e. the channes 1 S 0, 3 P 0, 3 P 1, 3 P 2 and 1 D 2 in spectroscopic notation 2S+1 L J, one has δ nn (E = δ (nn 00 +δ (nn 01 +3δ (nn 11 +5δ (nn 21 +3δ (nn 12 (26 using doube indices J in δ (ij J to indicate the partia wave. In Fig. 1 the effective phase shift for the sum of a partia waves with 2 is depicted in comparison with the pure s-wave phase shift. Since there is no experimenta data on the nn phase shifts, we use the np data taken from the Nijmegen partia wave anaysis [33] to iustrate the comparison between the experimentay known phase shifts and the resuts using the effective range approximation. Note, that the aboratory energy E ab = 2E of a neutron scattered on the target at rest is used as the argument. At very ow energies, both curves are rapidy rising. At energies above 10 MeV, contributions of higher partia waves become important. At ow temperatures, the viria integra is dominated by the ow-energy phaseshifts δ (ij and ony the s-wave contributes significanty due to the arge derivative with respect to the energy. he effective-range approximation [34] for the s-wave phase shift kcotδ (ij 0 = 1 a ij r ijk 2 (27 with momentum k = 2µ ij E and effective mass µ ij = m i m j /(m i +m j can be used to specify the energy dependence of δ (ij 0 with hep of the scattering ength a ij and the effective range r ij. Fig. 1 shows the resut for the s-wave phase shift using ony the scattering ength or both the scattering ength and the effective range parameter. he experimenta phase shift is nicey reproduced by the effective range expansion at very ow energies. he contribution of the k 2 term with r ij is essentia to obtain a decrease of the s-wave phase shift at higher energies. In this case, 10

11 δ nn (E ab [deg] S exp S exp +P exp +D exp a nn a nn +r nn E ab [MeV] Figure 1: Effective phase shift δ nn (E ab for neutron-neutron scattering versus aboratory energy E ab = 2E. Resuts are shown using experimenta data (incuding s, p and d-waves or s-waves ony and the effective-range approximation with and without the contribution depending on the effective range parameter, respectivey. the resut foows more cosey the experimenta data, however, differences with respect to the tota effective phase shift δ nn (E remain. he effective resonance energy E nn for nn scattering derived from the phase shifts in Eq. (24 is compared for different approximations in Fig. 2. In genera, E nn rises smoothy with increasing temperature indicating that it is not dominated by a particuar resonance. Contributions of higher partia waves reduce the effective resonance energy as compared to the pure s-wave resut. he s-wave viria integra can be obtained anayticay in the effectiverange approximation as I (ij 0 ( = 1 2 [ ] B (η ij 1 exp erfc A (η 2µ ij A (η ij ij η=±1 1 2µ ij A (η ij (28 for 2r ij /a ij 1. It depends on the coefficients ( A (η ij = a2 ij 1 r ij +η 1 2 r ij 2 a ij a ij (29 11

12 E nn [MeV] S exp S exp +P exp +D exp a nn a nn +r nn unitary imit modified a nn +r nn [MeV] Figure 2: emperature dependence of the effective resonance energy E nn for neutronneutron scattering in different approximations. See caption of Figure 1 and text for detais. and For r ij = 0 the integra reduces to B (η ij = a ( ij 1+η 1 2 r ij 2 a ij I (ij 0 ( = a ( ij 2 a ij exp 1 2µ ij a 2 ij. (30 erfc 1. (31 2µ ij a 2 ij Using these resuts, the effective resonance energy can be expressed expicity as a function of the temperature and the effective-range parameters. We use a nn = fm and r nn = fm from Ref. [35]. As shown in Figure 2, E nn in the effective range approximation rises more strongy with temperature when the contribution of the r nn term is taken into account. At high temperatures there remains a difference to the exact resut that takes higher partia waves into account. However, by a readjustment of the scatteringengthtoa nn = 11.31fmandtheeffectiverangetor nn = 1.06fm, the correct effective resonance energy obtained with a partia waves is very we reproduced for the reevant temperatures with the modified parameters 12

13 (denoted by modified a nn + r nn in Figure 2. Of course, these modified vaues for a nn and r nn are not physica and do not describe the scattering in a particuar partia wave. hey are ony an effective means to describe the temperature dependence of the effective resonance energy. here are severa interesting imits for the viria integra. For ow temperatures A (η ij /(2µ ij one finds I (ij 0 ( a ij µij 2π +... (32 (independent of the effective range r ij. In the unitary imit a ij and r ij = 0, that describes the situation where a bound state/resonance is just at the continuum threshod, I (ij 0 ( = 1/2 and the effective resonance energy is given by E ij = n2. his resut is aso depicted in Fig. 2. For positive scattering ength a ij, e.g. in the deuteron channe, the integra is negative and compensates part of the correation strength that is ocated in the bound state of the same channe. In Ref. [18] aso the α-partice was considered as a fundamenta constituent in the VEoS and the additiona viria coefficients b α, b nα, b pα and b αα appeared in the formaism. (For a comparison of the viria coefficients in the two formuations see Appendix A. In our approach, the corresponding contributions to Ω represent four-, five- and eight-body correations that wi become irreevant in the ow-density imit as compared to two-nuceon correations. Since chemica equiibrium requires µ α = 2µ n +2µ p +B α with the α-partice binding energy B α = MeV [32], the α-partice fugacity z α = exp(µ α /isnotindependent oftheneutronandprotonfugacities. he VEoS with neutrons, protons and α-partices was further extended in Ref. [19] to incude 3 H and 3 He as constituent partices corresponding to threebody correations. Care has to be taken in order to avoid a doube counting of states since, e.g., the α-partice states can be obtained in the 3 H-p and 3 He-nchannes(orthenegected 2 H- 2 Hchannes thatarepresent four-body correations when nuceons are assumed to be the ony fundamenta partices Reation to the nucear statistica equiibrium approach In principe, the viria expansion coud be extended to third, fourth,...order to incude more and more many-body correations. he second, third,... viria coefficients wi contain bound state contributions that cor- 13

14 respond to ground and excited states of heavy nucei. Negecting the manybody continuum correations, the grand canonica potentia becomes a sum Ω = V g ( A,Z µa,z exp (33 λ 3 A,Z (A,Z over a nucei with mass number A, charge number Z, chemica potentia µ A,Z = (A Zµ n + Zµ p and therma waveength λ A,Z = 2π/(m A,Z. he effective degeneracy factor g A,Z ( = g (0 A,Z + k ( g (k A,Z exp E k (34 is a sum over the ground state 0 and a excited states k of the nuceus A Z with excitation energies E k, tota anguar momenta J k and g (k A,Z = (2J k+1. he summation over individua states is often repaced by an integra over energy with an appropriate eve density, see, e.g. Ref. [11]. he temperature dependence of the degeneracy factor is a simpe and efficient means to incude excited states of nucei. It wi be used beow in the formuation of the grmf approach Generaized Beth-Uhenbeck approach he description of matter in the conventiona VEoS is imited to rather ow densities. However, the approach can be generaized in order to reach higher densities and to incude the effects of dissoution of composite partices, i.e. the Mott effect. his generaized Beth-Uhenbeck (gbu approach is based on a quantum statistica description with thermodynamic Greens functions. For detais see Ref. [36]. he grand canonica potentia Ω in the gbu approach assumes a simiar form as in the VEoS. he constituents are considered as quasipartices with the correct statistics (Fermi-Dirac or Bose-Einstein and with sef-energies that contain aready some effect of the mutua interaction. he remaining two-body correations (between quasipartices! are contained in a modified second viria coefficient that is given by b ij ( = 1+δ ij λ 3/2 i λ 3/2 j de f ij (E +E cont F ij (E (35 2 with the correct distribution function f ij (E of the composite system, i.e. Bose-Einstein for the two-nuceon-states, and the quantity 14

15 F ij (E = k + g (ij d 3 P k P>P Mott (2π g (ij π d 3 P 3 δ(e E(ij (2π 3 2sin2 δ (ij k (36 dδ (ij de that is reated to the in-medium density of states. he main differences as compared to Eqs. (16 and (17 are caused by medium effects. he properties of a two-body system ij depend on the tota c.m. momentum P of the state with respect to the medium. Hence, there is an additiona summation (integration over P. Bound states ony appear for P arger than the Mott momentum P Mott because at ower tota momenta their formation is suppressed by the action of the Paui principe. he phase shifts δ (ij (P,,µ i,µ j are determined by the in-medium -matrix describing the scattering of two quasipartices. Simiary E (ij k (P,,µ i,µ j is the mediumdependent energy of a bound state. Energies in Eq. (36 are measured with respect to the continuum edge E cont that corresponds to the energy in the frame of the medium where the reative momentum of the scattering partices is zero. he 2sin 2 δ (ij factor takes into account that part of interaction effects is aready contained in the sef-energies of the quasipartices. he second term in Eq. (36 is reated to the number of correated quasipartices that is different to the number of correated partices as cacuated by the usua Beth-Uhenbeck approach in Eq. ( Generaized reativistic mean-fied mode 3.1. Genera formaism In Ref. [25] a RMF mode with density-dependent meson-nuceon coupings was generaized to incude ight custers (bound states of 2 H, 3 H, 3 He, 4 He as additiona degrees of freedom with medium-dependent binding energies. hese custers coupe minimay to the meson fieds with a strength that was assumed to be a mutipe of the nuceon couping strength. Scattering states were not considered in the RMF mode of in Ref. [25]. he approach can easiy be extended when more than the above mentioned nucei need to be incuded. Such a mode with nucei as expicit degrees of freedom we ca the grmf mode in this paper. 15

16 he grand canonica potentia of the mode can be formuated as a function Ω(,V, µ i,ω,ρ,σ,δ depending on the temperature, voume V, reativistic chemica potentias µ i = m i + µ i of the partices and the spatiay constantmesonfiedsω, ρ, σ andδ underconsideration. Here, afourmesons are incuded to cover the four possibe combinations of scaars and vectors in Lorentz and isospin space. We negect antipartices in the foowing because they give sizabe contributions ony at temperatures much higher than those reevant in the considered astrophysica appications. In the genera case of inhomogeneous matter, not discussed here, Ω becomes a functiona of the spatiay varying meson fieds and their gradients. For uniform nucear matter, the grand canonica potentia can be written as Ω = [ 1 ( Ω i V m 2 2 ω ω 2 +m 2 ρ ρ2 m 2 σ σ2 m 2 δ δ2 (37 i + ( Γ ωωn ω +Γ ρρn ρ Γ σσn σ Γ δδn δ (nn +n p ] where Γ i are the derivatives of the density dependent coupings which determine the nuceon-meson couping strength. Contributions from individua partices (nuceons and custers are [ ( d 3 k Ω i = g i V (2π n 1±exp E ] i µ i (38 3 where the upper (ower signs refer to fermions (bosons. In the case of Maxwe-Botzmann statistics, Eq. (38 reduces to ( d 3 k Ω i = g i V (2π exp E i µ i. (39 3 We assume that the degeneracy factors g i can depend in genera on the temperature as in the case of NSE modes, cf. Section 2.3. he densities n ω, n ρ, n σ and n δ appearing in (37 are themseves functions of the temperature, chemica potentias and meson fieds (see beow. he mass of a custer i = d,t,h,α,... with N i neutrons and Z i protons is written as m i = N i m n +Z i m p (1 δ in (1 δ ip B (vac i (40 with neutron and proton rest masses m n and m p, respectivey, and the vacuum binding energy B (vac i > 0 for composite partices. For neutrons and 16

17 protons the binding energy does not appear. Chemica equiibrium eads to the constraint µ i = N i µ n +Z i µ p (41 for the reativistic chemica potentias of the custers. With every partice i = n,p,d,t,h,α,... of momentum k, an energy E i ( k = V i + k 2 +(m i S i 2 (42 is associated that contains the vector and scaar sef-energies V i and S i, respectivey. hese potentias are given by V i = Γ iω ω +Γ iρ ρ+(δ in +δ ip V (r, (43 S i = Γ iσ σ +Γ iδ δ (1 δ in (1 δ ip B i (44 with meson fieds m = ω,ρ,σ,δ and their coupings Γ im to the partice i. We assume that the couping strength of the custers is a mutipe of that of the nuceons, i.e. Γ im = g im Γ m ( (45 with g iω = g iσ = N i +Z i, (46 g iρ = g iδ = N i Z i. (47 Other choices are possibe and shoud be expored in the future. he mesonnuceon coupings Γ m ( depend on the tota nuceon density = n n + n p ony. he functiona form and the parameters of the density dependent coupings Γ m ( used in the grmf mode were chosen in such a way to describe the properties of finite nucei and nucear matter parameters at saturation density. he vector sef-energy contains the rearrangement term V (r = Γ ωn ω ω +Γ ρn ρ ρ Γ σn σ σ Γ δn δ δ (48 that is ony reevant for nuceons. In the medium, the actua binding energy of a custer B i = B (vac i B i wi be shifted with respect to the vacuum vaue B (vac i by a quantity B i. his medium-dependent binding energy shift B i for composite partices is parametrized as a function of temperature and the Lorentz vector meson fieds ω and σ, see, e.g. Ref. [25]. he specific form of these shifts is not important in the present discussion of the ow-density 17

18 behavior of the EoS but it affects the transition to high densities, see Section 7. he source densities in Eqs. (37 and (48 n ω = i n σ = i g iω n i, n ρ = i g iσ n (s i, n δ = i g iρ n i, (49 g iδ n (s i (50 depend on the vector and scaar partice number densities n i and n (s i, respectivey, that are defined as n i = Ω d 3 k µ i = g i,v, µj i (2π f 3 i, (51 d 3 k = g i (2π f m i S i 3 i (52 k2 (m i S i 2 n (s i with the Fermi-Dirac (Bose-Einstein distribution function f i ( k = [ exp ( 1 Ei µ i ±1] (53 or the Maxwe-Botzmann distribution function ( f i ( k = exp E i µ i (54 depending on the partice statistics. he fied equations for the mesons are found from the functiona Ω with the hep of the Euer-Lagrange equations. hey have the form m 2 ωω = Γ ω n ω + i (1 δ in (1 δ ip n (s B i i ω, (55 m 2 ρρ = Γ ρ n ρ + (1 δ in (1 δ ip n (s B i i ρ, i (56 m 2 σ σ = Γ σn σ, (57 m 2 δ δ = Γ δn δ (58 18

19 with additiona contributions to the source terms due to the dependence of the binding energies B i on the vector meson fieds. Finay, the entropy is obtained as S = V d 3 k g i (2π [f inf 3 i ±(1 f i n(1 f i ] (59 i V i (1 δ in (1 δ ip n (s B i i i Ω i dg i g i d with the usua singe-partice contribution, a term due to the temperature dependence of the custer binding energy shifts B i and a contribution caused by interna excitations of the custers, i.e. the introduction of the temperature dependence of the degeneracy factors, described in Sect. 4. In case of the Maxwe-Botzmann statistics for a partice i, the integrand contains ony the term f i nf i Scheme of the fugacity expansion Introducing the modified fugacity ( ( µi m i +S i V i Si V i z i = exp = exp z i (60 we write the distribution function as f i = [ z 1 i exp ( ei ±1 ] 1 z i exp ( e i = (61 1± z i exp ( e i with the kinetic energy e i (k = k 2 +(m i S i 2 (m i S i 0. For z i 1 the distribution function can be expanded in a power series in z i exp( e i /. he appearing momentum space integras can be evauated expicity [37] with the integra representation of the modified Besse functions K ν (x [38]. For the contributions (38 of the individua partices to the grand canonica potentia we find Ω i = V g i n=0 λ 3 i ( 3/2 mi S i (62 m i ( 1 n (n+1 5/2k 2 [ (n+1 m ] [ i S i exp (n+1 S ] i V i z n+1 i 19

20 with the non-reativistic fugacity z i and the functions 2x k ν (x = π exp(xk ν(x. (63 Simiary, the vector density is obtained as n i = g i λ 3 i ( 3/2 mi S i (64 n=0 m i ( 1 n (n+1 3/2k 2 [ (n+1 m ] [ i S i exp (n+1 S ] i V i zi n+1. For the scaar density n (s i, k 2 has to be repaced with k 1. In case of Maxwe- Botzmann statistics, ony the n = 0 term remains. Without interaction, i.e. S i = V i = 0, the resuts for a reativistic Fermi or Bose gas are recovered. For 0, i.e. x = (n+1(m i S i /, one finds from the asymptotic expansion [38] k ν (x = 1+ µ 1 8x + (µ 1(µ 9 + (µ 1(µ 9(µ (65 2!(8x 2 3!(8x 3 with µ = 4ν 2 the common expressions and Ω i = V g i λ 3 i n i = g i λ 3 i m=0 m=0 ( 1 m (j +1 5/2zm+1 m, (66 ( 1 m (m+1 3/2zm+1 m (67 for non-reativistic partices. In this imit, n (s i = n i, because the scaar density differs from the baryon density ony in the reativistic description. Since the scaar and vector sef-energies are themseves functions of the densities, a second expansion of the series (62 and (64 is required Fugacity expansion of the grand canonica potentia up to second order In order to compare the true series expansion of the grand canonica potentia Ω of the grmf mode with the form (18 in the VEoS approach at ow densities, we expand (37 up to second order in the fugacities of neutrons 20

21 and protons. In the foowing we wi ony consider neutrons, protons and deuterons in the density expansion. Nucei with mass numbers A 3 do not contribute in the second order of expansion since we ony consider nuceons as basic constituents. he expicit contribution of the deuteron ground state appears in the two-nuceon correation term with ( Bd z d = z n z p exp. (68 For sufficienty ow nuceon densities, the sef-energies of the nuceons are approximatey inear in the nuceon densities. Hence, the contribution of the individua nuceon can be approximated as Ω i V g i λ 3 i 1 2 5/2k 2 { ( mi [ ( k 2 1+ S i ( 2 m } i zi 2. Appying the recursion reation k ν (x = ( 1 1 k 2 ( mi ( k mi m i V i ] z i (69 2ν 1 k ν k ν 1 (x, (70 2x the expression (69 for the nuceons reduces to [( Ω i V 1 V i x i + S i y i g ( i k 2 5/2 λ 3 2 i 2 m i ] zi 2 (71 with the abbreviations x i = g i λ 3 i y i = g i λ 3 i ( mi k 2 k 1 ( mi z i, (72 z i. (73 For the deuteron contribution we have Ω d = V g ( d md k λ 3 2 z n z p exp d ( Bd (74 21

22 without sef-energy terms which contribute ony in higher order of the fugacities. Contrary the nuceon sef-energies contribute aready in the owest order of the fugacities with coefficients V n C ω (x n +x p +C ρ (x n x p, (75 V p C ω (x n +x p C ρ (x n x p, (76 S n C σ (y n +y p +C δ (y n y p, (77 S p C σ (y n +y p C δ (y n y p (78 C m = Γ2 m(0 m 2 m (79 for m = ω,ρ,σ,δ that depend on the density dependent meson coupings Γ m ( taken in the imit 0. he rearrangement term in the vector sefenergy does not contribute at this eve because it is at east quadratic in the densities. With the hep of the fied equations (55-(58, the mesonic contributions in (37 can be expressed as quadratic forms of the nuceon densities. Finay, after performing a necessary expansions up to second order in the nuceon fugacities, we obtain the resuting grand canonica potentia Ω(,V,µ n,µ p = (80 V g [ ( n mn k λ 3 2 z n 1 ( n 2 5/2k 2 2 m ] n zn 2 V g [ ( p mp k λ 3 2 z p 1 ( p 2 5/2k2 2 2 m ] p zp 2 V g ( d md ( Bd k λ 3 2 exp z n z p d + V g 2 [ ( n (C 2 λ 6 ω +C ρ k2 2 mn ( (C σ +C δ k 2 mn ] 1 zn 2 n + V g 2 [ ( p (C 2 λ 6 ω +C ρ k2 2 mp ( (C σ +C δ k 2 mp ] 1 zp 2 p +V g [ ( ng p mn ( mp (C λ 3 nλ 3 ω C ρ k 2 k 2 p ( mn ( mp ] (C σ C δ k 1 k 1 z n z p with contributions from free nuceons and their correation due to statistics, the deuteron bound state and the interaction. 22

23 3.4. Comparison of fugacity expansions A comparison of Eq. (18 with the corresponding expansion (80 permits to extract the viria coefficients b n, b p, b nn, b pp and b np of the grmf mode. he first-order coefficients receive a reativistic correction with ( mn b n = g n k 2, (81 ( mp b p = g p k 2. (82 heydependonthetemperaturenow. Inthenon-reativisticimit /m i 0 the correction wi vanish due to the asymptotic expansion of the function k 2, cf. Eq. (65. Beow we wi show that these corrections are important and shoud be taken into account for an improved description of ow-density nucear matter. Simiary, there is a reativistic modification to the statistica corrections in b nn and b pp, i.e. Eqs. (19 and (20, respectivey, and to the deuteron contribution to (22. From the comparison of (18 and (80, three independent reations for the channes nn, pp and np remain ( ] λ 3 n = 1 2 λ 3 p = 1 2 [ [ b nn + g n 2 5/2k 2 gn 2 λ 6 n 2 m n [ (C ω +C ρ k2 2 [ b pp + g p 2 5/2k 2 2 = = g ng p λ 3 n λ3 p ( 2 m p = λ 3 nn λ 3 pp k 2 ( 2mn ( mn ( mn ] (C σ +C δ k1 2, ] = ( 2mp k 2 ( mp ( mp g (nn I (nn (83 g (pp I (pp (84 gp 2 [ ] (C λ 6 ω +C ρ k2 2 (C σ +C δ k1 2, p b np g ( d md ( ] Bd k λn 3/2 λp 3/2 λ 3 2 exp (85 d ( mn +m p [ ] k λ 3 2 g (np1 I (np1 +g (np0 I (np0 np [ ( mn ( mp (C ω C ρ k 2 k 2 ( mn ( mp ] (C σ C δ k 1 k 1 that connect the viria integras with the strengths C m of the zero-density meson-nuceon coupingsγ m (0ofthegRMFmode. heexpected correction 23

24 factors k 2 (m i / due to the reativistic effects were added to the VEoS part (expressed through viria integras I ij of reations (83 to ( emperature independent imit of consistency conditions he consistency reations (83 to (85 were derived from a comparison of two fugacity expansions that are vaid at ow densities and finite temperatures, i.e. in the viria imit. hese conditions arise from the contributions in the fugacity expansion that are quadratic in the fugacities. We perform an expansion of Eqs. (83 to (85 in powers of sma /m i and keep the owest-order terms which are independent of. his is equivaent of taking the imit 0 in (83 to (85. In this imit there are no reativistic corrections, ony s-wave contributions are reevant and the imit (32 of the integra with the scattering engths can be used. In the grmf mode the nuceonnuceon interaction is identica in the nn and pp systems. hus, the first two equations, (83 and (84 cannot be considered independent but shoud be combined. Denoting the scattering engths in the s-wave channes expicity with a (nn 1 S 0, a (pp 1 S 0, a (np 1 S 0 and a (np 3 S 1 and considering the degeneracy factors g n = g p = 2, g (nn 1 S 0 = g (pp 1 S 0 = g (np 1 S 0 = 1 and g (np 3 S 1 = 3, the two consistency conditions C ω C σ = π { [ (nn 1 a 1 S 0 + a 2 m n (pp 1 S 0 m p ] + m (np n +m p a m n m p 1 S 0 +3a (np 3 S 1 4 } (86 and C ρ C δ = π { [ (nn 1 a 1 S 0 + a 2 m n (pp 1 S 0 m p ] m (np n +m p a m n m p 1 S 0 +3a (np 3 S 1 4 } (87 for the isovector and isoscaar meson coupings are obtained. he eft side of Eqs. (86 to (87 corresponds to the grmf expansion whie the right side represents the viria expansion in the effective-range approximation. With the couping constants at zero density Γ m (0 and the meson masses m m the couping coefficients C m can be cacuated for different RMF modes. In abe 1 the resuts are shown for the grmf mode with parametrization DD2 without δ meson foowing Ref. [25] and for a new density dependent RMF parametrization DD-MEδ from [39] based on ab-initio cacuations in nucear matter incuding a Lorentz scaar isovector δ meson. In the DD2 parametrization the differences of the couping coefficients are 24

25 abe 1: Couping coefficients (79 in the DD2 parametrization of the grmf mode [25] and the DD-MEδ parametrization of Ref. [39]. Mode DD2 DD-MEδ meson m C m [fm 2 ] C m [fm 2 ] ω σ ρ δ C ω C σ = 5.39 fm 2 and C ρ C δ = 2.48 fm 2, respectivey. In contrast, the cacuation with the scattering engths from Ref. [35] gives C ω C σ = fm 2 and C ρ C δ = 9.61 fm 2. It is obvious that the couping coefficients of the grmf mode do not obey the consistency reations to guarantee the agreement with the VEoS at ow densities. For the recenty deveoped parametrization DD-MEδ [39], the differences ofthecouping coefficients turnouttobec ω C σ = 4.90 fm 2 andc ρ C δ = 2.55 fm 2. hese numbers are cose to those of the DD2 parametrization. here are sti differences to the vaues required by the consistency conditions (86 and (87 and a more drastic modification of the ow-density coupings seems to be necessary. In the DD2 parameter fit of the grmf approach a particuar form of the density dependence of the couping functions Γ m ( was assumed. he choice was motivated by resuts from Dirac-Brueckner cacuations of nucear matter [26]. However, the vaues at zero density are found from a not we constrained extrapoation to sma densities. Ony near the saturation density of symmetric nucear matter they are we determined from mode fits to the properties of atomic nucei. hus, in genera, one coud modify the density dependence of Γ m ( for a mesons m at ow densities in order to satisfy the constraints (86 and (87. Since C ρ C δ has to be negative, it is cear that the standard ω, σ and ρ mesons are not sufficient and a δ meson is needed in this case. Since Eqs. (86 and (87 give ony two reations for four couping coefficients C m, obtained from the consistency with the viria expansion, two further constraints are required to fix a coupings unambiguousy. Hence, in addition to the viria constraints, we tried to fit the vacuum coupings 25

26 coefficients Γ m (0 by describing ow-energy nuceon-nuceon scattering, i.e. experimentay determined NN phase shifts [33], simutaneousy satisfying equations (86 and (87. he corresponding Schrödinger equation in coordinate space was soved using σ, ω, ρ and δ mesons with a one-boson exchange potentia [40] in order to fit the couping constants of the nuceon-meson interaction and to describe the experimentay known phase shifts and effective range parameters given in [35], [41]. he gradient terms in the potentia ead to unphysica singuarities at r 0, thus to avoid these singuarities a certain cutoff was introduced as in [42]. Couping constants derived in this way woud serve as an additiona constraint at zero densities, improving the density dependence of the couping functions Γ m (. However, even when the contribution of the pion was incuded, that does not appear in the mean-fied approximation, it was impossibe to find a reasonabe parametrization. Quite drastic modifications of the ow-density meson-nuceon coupings and meson masses were required. his faiure is not very surprising because the number of meson fieds in the RMF mode is much smaer than that in reaistic NN potentias resuting in a smaer number of degrees of freedom. In addition, the meson fieds in the RMF mode do not represent rea physica mesons and their masses, that are used as adjustabe parameters, e.g. in case of the σ meson, are not identica to the true meson masses. Hence, we do not foow this approach but extend the grmf mode further. 4. Extension of the grmf approach he faiure of the grmf mode to describe the viria imit suggests that a further extension of the approach is inevitabe. In addition to the bound state two-body correation, the deuteron, we introduce additiona degrees of freedom in the grand canonica potentia that represent the possibe twonuceon states, nn, pp and np, in the continuum. In the foowing, we wi consider the isospin 0 and 1 channes for the np states separatey. Each of the two-body channes is represented by a temperature dependent resonance energy (Ẽnn, Ẽpp, Ẽnp0, Ẽnp1 andacorresponding effective degeneracy factor ( g nn, g pp, g np0, g np1. In genera, the resonance energies and degeneracy factors do not necessariy have to be the same as the resonance energies of Eq. (24 or the standard degeneracy factors in Eqs. (19 to (22. wo-body correations in the continuum are considered as two-body custers with mass m ij = N ij m n +Z ij m p +Ẽij (88 26

27 with a temperature dependent resonance energy Ẽij. We assume that the energy shift of a custer B ij that appears in the scaar sef-energy (44 is identica to the binding energy shift of the deuteron B d for a continuum channes. Of course, this most simpe choice can be substituted by improved descriptions in the future. he specific functiona form of the medium-dependent energy shift wi affect the form of the transition from very ow densities, where the VEoS is appicabe, to the high-density regime. However, for the ow-density matching of the VEoS and grmf mode it is irreevant. he tota rest mass of the custer in the medium is given by M ij = m ij + B ij. (89 In case of the viria approach, the appearing masses are given by m ij = N ij m n + Z ij m p +E ij and M ij = m ij, respectivey. We incude the effect of the resonance energies in the reativistic correction through the k 2 function and the definition of the therma wave engths λ ij = [2π/(M ij ] 1/2 and simiary for λ ij. With these modifications we can rederive consistency reations between this extended grmf mode and the VEoS. he viria integras (23 in the viria coefficients are represented by an exponentia of the effective resonance energy as defined in (24. he two np channes are separated by requiring for the np isospin 1 channe a simiar form as for the nn and pp channes. Comparing the fugacity expansion of the grmf mode with 2-body correations with the VEoS at ow densities, we finay obtain an extension of the reations (83 to (85, where an additiona term with temperature dependent resonance energy Ẽij( appears in the right-hand side of the Eqs. (90 to (91, i.e. ( Mij k λ 3 2 ij ĝ ij exp ( E ij ( ( = λ3 Mij k 2 g (eff ij exp ij Ẽij + 1 g i g [ ( j mi ( mj (C 2λ 3 ω +C ρ k 2 k 2 i λ3 j ( mi ( mj ] (C σ +C δ k 1 k 1 (90 27

28 for the channes ij = nn, pp, np1 and ( Mnp0 k λ 3 2 np0 ĝ np0 exp ( = Mnp0 k λ 3 2 np0 + 1 g n g [ p 2λ 3 nλ 3 p (C σ 3C δ k 1 ( mn ( E np0 g (eff np0 exp ( Ẽnp0 ( mn ( mp (C ω 3C ρ k 2 k 2 ( mp ] k 1 (91 for the channe np0. hese consistency reations eave sufficient freedom to reproduce the exact VEoS at ow densities in the grmf mode. here are different possibiities to achieve this goa depending on the choice of the zero-density coupings and on the choice of the effective resonance energies and degeneracy factors that represent the continuum contributions in the extended grmf mode. he structure on the eft-hand side shares some simiarities with the generaized Beth-Uhenbeck approach in Section 2.4. One part of the effects of the two-body correations is incuded in the expicit continuum contributions and the remaining part is covered by the mean-fied terms. here is the freedom to shift the correation strength between these terms rather arbitrariy. In the foowing, a few choices wi be discussed. In Section 3.5 the case without introducing expicit continuum channes in the grmf mode, i.e. shifting the fu correation strength into a redefinition of the vacuum coupings, turned out to be impossibe. he other extreme is to represent the correations by the introduction of the effective continuum states described by temperature dependent resonance energies Ẽij with medium corrections B ij as described above. In this case the contribution of the viria integra in the eft-hand side of the Eqs. (90 and (91 is fuy represented by the the resonance energies Ẽij in the first term in the right-hand side of these equations. herefore the mean-fied part has to vanish in the non-reativistic imit, eading to C ω C σ = 0 and C ρ C δ = 0. he couping coefficients C m in reations (86 to (87 were expressed through the viria integra taken in the effective range approximation. In the new reations, given above, the viria integra and thus the scattering engths are impemented in the definition of the resonance energies Ẽij(. herefore, negecting the term proportiona to g (eff ij exp ( Ẽij/ 28 in the reations (90 and (91, we