Relativistic Mean-Field Models with Different Parametrizations of Density Dependent Couplings

Relativistic mean-field models are successfully used for the description of finite nuclei and nuclear matter. Approaches with density-dependent meson-nucleon couplings assume specific functional forms and a dependence on vector densities in most cases. In this work, parametrizations with a larger sample of functions and dependencies on vector and scalar densities are investigated. They are obtained from fitting properties of finite nuclei. The quality of the description of nuclei and the obtained equations of state of symmetric nuclear matter and neutron matter below saturation are very similar. However, characteristic nuclear matter parameters, the equations of state and the symmetry energy at suprasaturation densities show some correlations with the choice of the density dependence and functional form of the couplings. Conditions are identified that can lead to problems for some of the parametrizations.


Introduction
A realistic description of dense matter is essential for the physics of compact stars and the simulation of core-collapse supernovae and neutron-star mergers; see, e.g., [1] for details.Depending on the thermodynamic conditions, a variety of particle species has to be considered in the construction of microscopic models; in particular, nucleons and electrons are indispensable in astrophysical applications.The thermodynamic properties of such matter are encoded in the equation of state (EoS).In general, the EoS depends on several variables, e.g., the baryon number density n B , the temperature T and the electron fraction Y e .They cover several orders of magnitude in the application, and the properties of matter change dramatically within the corresponding ranges.The main theoretical challenge is the description of the hadronic component, i.e., the subsystem composed of nucleons (and possibly other baryons such as hyperons at very high densities) that interact strongly.The leptonic component can be well treated in a simple Fermi gas model in most cases.
At baryon densities below the nuclear saturation density n sat ≈ 0.16 fm −3 and temperatures not exceeding approximately 1.7 × 10 11 K or k B T ≈ 15 MeV, matter is not homogeneous, and structures develop due to the competition of the short-range strong interaction, the long-range Coulomb interaction and the entropy.For instance, neutron-rich nuclear clusters can form that arrange on a lattice in the crust of neutron stars.At densities around and above n sat , e.g., in the core of compact stars, or at high temperatures, matter is expected to be homogeneous.The description of such uniform hadronic matter is the subject of this work.
There are many microscopic theoretical approaches to construct EoSs of strongly-interacting matter.One class of models is based on realistic interactions, which are fitted to nucleon-nucleon scattering data and properties of light nuclei, in combination with sophisticated techniques to solve the many-body problem.In contrast, in a second class of models, phenomenological approaches proceed more heuristically employing effective interactions with a number of parameters that can be determined from experimental data of different origin.These approaches are often derived from self-consistent mean-field (MF) models [2].They can be expressed as energy density functionals (EDF), which, in principle, are able to capture the essential features of the system.Amongst them, EDFs based on non-relativistic Hartree-Fock calculations with Skyrme or Gogny potentials are most familiar [3,4].Relativistic or covariant EDFs can be constructed from relativistic mean-field (RMF) models that describe the strong interaction by an exchange of mesons, which are are usually assumed to couple minimally to nucleons.
For a quantitative description of nuclei and nuclear matter, medium effects of the effective interaction have to be included.This can be achieved in different ways.In a large number of RMF models [5], nonlinear (NL) self-interactions of the mesons are incorporated in the Lagrangian density L that defines the model.Alternatively, the meson-nucleon couplings in L can contain an explicit dependence on the nucleon field operators that is mapped to a density dependence of the couplings in the derived EDF.The choice of the functional form of this dependence facilitates a very flexible variety of models.The parameters of the density dependence are customarily determined by applying the EDF to the description of finite nuclei and by fitting to a selected set of their properties.
A particular property of RMF models is the occurrence of two different particle number densities: vector densities and scalar densities.Their interplay is essential to describe the saturation of nuclear matter in the model.RMF models with density-dependent (DD) couplings were first applied to the description of finite nuclei by Fuchs and Lenske [6,7] considering dependencies on vector and scalar densities.The couplings were obtained from microscopic Dirac-Brueckner-Hartree-Fock calculations of nuclear matter.The first self-consistent RMF model with density-dependent couplings that were fitted to properties of finite nuclei used the so-called vector density in the couplings and specific forms for the functional dependence [8].Many subsequent models and parametrizations followed this approach, sometimes introducing different functions; see, e.g., [9][10][11][12][13][14][15].A dependence on other densities, in particular the scalar density, in the description of nuclei was not really explored.The aim of this work is a comparison of RMF models and the corresponding nuclear matter EoSs with DD couplings of different functional form and dependencies on vector and scalar densities that were fitted to the same set of nuclear observables.
The formalism of RMF models with DD couplings is presented in Section 2 with applications to homogeneous nuclear matter and finite nuclei.The parametrization of the meson-nucleon couplings is introduced in Section 3. The determination of the parameters and the form of the couplings is discussed in Section 4. Properties of nuclear matter and the EoS are summarized in Section 5 including the characteristic nuclear matter parameters and the density dependence of the symmetry energy.In Section 6, constraints on the possible density dependence of the couplings are considered.Conclusions are given in Section 7, and the actual values of the different parameter sets are collected in Appendix A. Natural units with h = c = k B = 1 are used throughout this paper.

RMF Model with Density-Dependent Couplings
The theoretical description of nuclei and nuclear matter in the present RMF approach proceeds in the usual way as presented, e.g., in [8].The strong interaction is described by an exchange of isoscalar σ and ω mesons and isovector ρ and δ mesons.They carry the same quantum numbers as the corresponding experimentally-known mesons, but they are not necessarily the same.However, in this way, it is possible to capture the main features of the effective in-medium nuclear interaction.The mesons couple minimally to the nucleons i = n, p that are represented by four-spinor operators Ψ i forming an isospin doublet.The coupling with a meson j = σ, ω, ρ or δ is denoted by Γj , a functional that can depend on suitable combinations of Ψ i and Ψ i = Ψ † i γ 0 , where γ µ is a Dirac matrix.In the present work, a dependence on the vector density: with the nucleon current: or the scalar density: is considered.Despite the name, the vector density ( 1) is a Lorentz scalar, as well as the scalar density (3).This guarantees the Lorentz covariance of the approach.In addition to the meson fields, a coupling to the electromagnetic field, represented by the vector field A µ , with constant coupling strength Γ γ is taken into account.

Lagrangian Density and Energy Density Functional
The starting point of the formalism is the Lagrangian density: with three contributions.The first: contains the covariant derivative: where τ k (k = 1, 2, 3) are isospin matrices in analogy with the Pauli matrices σ k , but acting in isospin space.The effective mass operator has the form: with the nucleon masses m i in vacuum.The ρ and δ terms carry an arrow since they have three components in isospin space.In contrast, the σ and ω fields are isoscalar quantities.The meson term in Equation ( 4) has the standard form: with meson masses m j and field tensors: for the vector mesons.The last contribution in Equation ( 4) is the Lagrangian density of the electromagnetic/photon field: with It is only relevant for finite nuclei, and not for homogeneous nuclear matter.The field equations for all degrees of freedom are derived with the help of the Euler-Lagrange equations in the standard mean-field and no-sea approximation.Mesons and photons are treated as classical fields.This leads to the Dirac equation for the nucleons, Klein-Gordon and Proca equations for scalar and vector mesons, respectively, and the Maxwell equation for the electromagnetic field.Using symmetries, e.g., by a restriction to a stationary system and choosing a particular frame of reference, i.e., breaking Lorentz covariance, the equations simplify further.
In the MF approximation of the present approach, the coupling functionals Γj become simple functions Γ j (n (v) , n (s) ) of the total vector density n p with: and the total scalar density n (s) = n p with: where w ik denotes the occupation factor of the single-particle state k of nucleon i.The nucleon wave functions ψ ik are solutions of the time-independent Dirac equation: with scalar and vector potentials S i and V i , respectively.They are given by: with factors g nσ = g pσ = g nω = g pω = −g nρ = g pρ = −g nδ = g pδ = g pγ = 1 and g nγ = 0 that reflect the different coupling of neutrons and protons to the fields.Due to symmetries, only a single component of the Lorentz vector and isospin vector fields remains; hence, the notation is simplified to δ, ω 0 , ρ 0 and A 0 without an additional index for the isospin in the following.The potentials (13) and ( 14) contain the rearrangement terms: that arise due to the density dependence of the couplings Γ j .These terms are identical for all particles since the couplings depend only on the sum of the scalar and vector densities, n (s) and n (s) in the present work.The source densities in Equations ( 15) and ( 16) are found from: for scalar mesons j = σ, δ and: for vector mesons j = ω, ρ by simple summations over the nucleons.The energy density: is obtained from the energy-momentum tensor T µν where the brackets denote the summation over all occupied states in the system.The contribution of the nucleons has the form: and the contribution of the fields reads as: The energy density ( 19) is a functional of the nucleons (ψ ik , ψ ik ), the fields (σ, ω 0 , δ, ρ, A 0 ) and their derivatives ( ∇σ, ∇ω 0 , ∇δ, ∇ρ 0 , ∇A 0 ).Similar to the Lagrangian density (4), it can be used to derive the field equations: with the source density: in the Poisson Equation (26) for the electromagnetic potential A 0 .

Homogeneous Nuclear Matter
In this case, the theoretical description simplifies further.The solutions of the Dirac Equation ( 12) are plane waves with momentum p and energy: for particles (η = +1) and antiparticles (η = −1), which have to be included when finite temperatures are considered.The meson fields and densities are constants in space.There is no contribution from the electromagnetic field because the source density n γ is constant, and, thus A 0 = 0.The energy density can be written as: with the quantities: and: and the kinetic contribution: in the continuum limit containing the degeneracy factor g i = 2 for the spin 1/2 nucleons.Rearrangement contributions in the energy density (29) appear only if the couplings depend on the scalar density n (s) , i.e., D = 0.The Fermi-Dirac distribution functions: in Equation (32) depend on the energy E (η) i , the chemical potential µ i and the temperature T. These functions also appear in the vector and scalar densities: with the Dirac effective mass: Additional thermodynamic quantities are easily calculated.The entropy density assumes the standard form: The pressure: can be obtained from the energy-momentum tensor.The result is identical to that of the thermodynamic definition: as a derivative of the energy density with respect to the total baryon density n B = n p for constant temperature T and isospin asymmetry δ = (n The rearrangement terms are essential for the thermodynamic consistency.
For vanishing temperature, as considered in most cases of this work, there are no antiparticle contributions, and analytical results are available for the vector density: with the Fermi momentum p * i and the scalar density: with the effective chemical potential: The expression (32) reduces to: and the pressure can be written as: with: and the kinetic contribution: Obviously, the entropy density (37) vanishes for T = 0.

Finite Nuclei
The energy of a nucleus with N neutrons and Z protons in the RMF calculation is found as: with single-particle energies E ik and corresponding occupation numbers w ik of the single-particle states k of nucleon i.They satisfy the normalization conditions: Using the field equations, the contribution of the meson and electromagnetic fields in Equation (47) can be written as: with factors: that contain derivatives of the coupling functions.
In the present work, only spherical nuclei are considered.In this case, the single-particle wave functions can be written as: with real radial wave functions F iκ k and G iκ k and spin-spherical harmonics Y κ k m k .They depend on the quantum number κ k = (j k − l k )(2j k + 1) = ±1, ±2, . . ., which determines the total angular momentum the state, and the projection m k = −j k , . . ., j k , i.e., there is a 2j k + 1 = |κ k |-fold degeneracy of the level with energy E ik .Assuming equal occupation of the sub-states, the vector and scalar single-particle densities are given by: with the normalization: This assures the sphericity of the source densities and potentials.
For a comparison with experimental data, the MF energy (47) and quantities related to the density distributions have to be corrected.The Coulomb potential A 0 in the calculation of the single-particle wave functions of nucleus with Z protons is multiplied with a factor (Z − 1)/Z to have the correct asymptotic dependence of the field when a proton is separated from the nucleus.This is necessary because in the MF approximation, exchange terms are missing that would correct this error.A similar correction for the meson fields is not applied because they are of short range.
Since the nuclear wave functions are fixed in the calculation to the origin of a spherical coordinate system, the translation symmetry is broken and a center-of-mass correction has to be applied.The energy (47) contains a 'localization' contribution E (cm) N,Z that has to be subtracted.Here, in a non-relativistic approximation, the expectation value: for the nucleus with A = N + Z nucleons with the total momentum ˆ P = ∑ A n=1 ˆ p n is used.In the same spirit, a correction for the density distributions is implemented.In a first step, the point particle distributions are converted to form factors by Fourier transformations.These are multiplied by the correction factor: with M N,Z = Nm n + Zm p for momentum q and form factors for the charge distributions of neutrons and protons if required.An inverse Fourier transformation yields the corrected distributions that can be used, e.g., in the calculation of root-mean-square radii.Further corrections, e.g., from pairing or particle-vibration couplings, are not taken into account.They can be considered in future extensions of this work.

Parametrization of Couplings
The density-dependent couplings Γ j for j = σ, δ, ω, ρ are the central quantities that determine the quality of the relativistic density functional.They are usually written in the form: with a constant coupling Γ (0) j = Γ j (n ref ) at a reference density n ref and an arbitrary function f j that depends on the ratio x = n/n ref .The density n can be any density that is formed as a Lorentz scalar from the single-particle wave functions ψ i and ψ j .The most frequent choice is a dependence on the total vector density that reduces to the sum: with the total neutron and proton vector densities n (n) n and n (n) p , respectively, in a system at rest without nucleon currents.In this work, also a dependence on the scalar density: is explored that has been suggested already when the RMF model with DD couplings was developed, but the quality of describing nuclear matter or finite nuclei was not examined in detail.
For a dependence on n (v) , the reference density is chosen as the vector density at saturation n sat , whereas the scalar density at saturation n sat is used for a dependence on n (s) .The most widely-used form for f j is a rational function: with four parameters a j , b j , c j and d j .It was introduced in [8] for the couplings of the σ and ω meson because such a function could describe very well effective density-dependent couplings that were extracted from self-energies in Dirac-Brueckner calculations of nuclear matter.For the ρ meson, a simple exponential form: with only one parameter a i was used in [8].Subsequently, also other functions were devised.
For instance, a generalization of (60) as: was used in [12] with also four parameters as in Equation (59).A specific feature of the functions (59) and ( 60) is that they are well-behaved for x → ∞ approaching a constant a j b j /c j or zero, respectively.In contrast, the function (61) diverges for c j = 0.A still different form was introduced in [13] as: for all mesons with an additional modification for the ω meson close to the saturation density.The latter two functional forms will not be used in this work.
In order to reduce the number of free parameters in the rational function (59), several conditions are demanded.First of all, it is required that f j (1) = 1.This fixes the first parameter in (59) as: and leads to the fact that the prefactor in Equation (56) can really be identified with the value of the coupling at the reference density n ref .
The function (60) automatically conforms to this condition.In [8], a second condition was introduced, requiring that the curvature of the function (59) vanishes at x = 0, i.e., the derivative is f j (0) = 0.This is met if the relation 1 = 3c j d 2 j holds.In [8], only the solution with d j > 0 was explored, whereas the second possibility d j < 0 will be investigated here, as well.Instead of introducing a condition for the second derivative of the function f j at x = 0, the first derivative at this point can be set to a specific value as a third option.In this work, the choice f j (0) = 0 is examined, corresponding to d j = 0.With these conditions on the function at x = 0 and x = 1, the number of independent parameters reduces to two, one more than for the exponential function (60).
Several combinations for the choice of the density in the argument x, the functional form of the function f j and the conditions on f j are explored in this work.In order to summarize these options in a concise form, a three-letter abbreviation is introduced for the identification.There are three cases for the first letter: ref and f σ and ref .
The second letter indicates the condition on the rational function for the couplings of the ω and σ mesons: The last letter denotes the coupling of the ρ and δ mesons: • 'E': exponential form (60) of the function f j , • 'R': rational form (59) of the function f j with f j (0) = 0, d j = 0 and f j (1)/ f j (1) = f j (1)/ f j (1).
The last condition of the case 'R' is motivated by the density dependence of the exponential function close to x = 1 and reduces the number of independent parameters by one.In total, there are 3 × 3 × 2 = 18 different combinations of functions tested in this work.

Determination of Parameters and Couplings
As every phenomenological model, the RMF approach with DD couplings depends on a number of parameters that need to be determined by a comparison of model predictions to experimental data.The resulting values will depend on the chosen observables and the method of the fitting procedure.In the present model, the parameters comprise the masses of the nucleons m i (i = n, p) and mesons m j (j = σ, ω, δ, ρ), the couplings Γ j (n ref ) of the mesons at the reference density, the parameters of the functions f j and the electromagnetic coupling constant Γ γ .Not all of them are assumed to be free parameters that can be varied more-or-less arbitrarily.
The masses of the nucleons and of the mesons, except the σ meson, are taken as preassigned with m n = 939.565413MeV, m p = 938.272081MeV, m ω = 783 MeV, m ρ = 763 MeV and m δ = 980 MeV.The coupling Γ γ is given by the experimental value.This leaves one mass (m σ ), four couplings at the reference density (Γ j (n ref )) and six parameters of the functions f j (two for two isoscalar and one for two isovector mesons) as free parameters of the model.Thus, in total, there are at most 11 parameters that have to be determined.In this work, a contribution of the δ meson will not be considered, hence there are N para = 9 free parameters left, which are denoted by p k in the following.In the actual fitting procedure, the six parameters for the couplings of the isoscalar mesons are not used directly because they are highly correlated.Instead, they are replaced by the characteristic values of nuclear matter parameters, i.e., the saturation density, the binding energy per nucleon and the effective mass at saturation, the incompressibility and two parameters related to the ratios of the coupling function The conversion between these parameter sets is analytic and easily implemented.
The actual set of parameter values {p k } is obtained from a least-squares fit by minimizing the function: and weighted by the inverse of an assigned uncertainty ∆O n .The latter are not the experimental errors, but values that reflect the hopefully achievable uncertainties of the model.The relative size for different observables also determines their relative importance in the fit.The selection of the observables and their uncertainties directly influences the results.There are different approaches to set up the fitting protocol.In this work, only observables of finite nuclei are considered, but no data that are derived more indirectly like nuclear matter parameters.This choice is motivated by the perception that data should be included that are close to the experimental determination and that show a different sensitivity to the various model parameters.In the parameter fits presented in this work, binding energies of nuclei, data related to the charge form factor and spin-orbit splittings are used as observables as far as available for a set of 12 magic and semi-magic nuclei.See Table 1 for the selected nuclei, observables and the assumed uncertainties.In principle, a larger set of nuclei could be included in the fit, but the calculation time would increase substantially.However, the selected set already contains sufficient information to determine the parameters and to allow a comparison of different functional forms of the density-dependent couplings.The binding energies are taken from the 2016 atomic mass evaluation (AME2016) [16].Information on the size of a nucleus and the density distribution is encoded in the charge form factor F(q).It determines the charge radius r charge , the diffraction radius r diff and the surface thickness σ surf .These three quantities are related to the curvature of F at momentum transfer q = 0, the position of the first zero and the height of the second extremum.For details, see [17].The most recent updated values for charge radii are taken from [18].Diffraction radii and surface thicknesses are extracted from charge form factors calculated with charge distributions of nuclei in [19].Spin-orbit splittings are deduced from the level spectra of the nuclei included in the fit and their neighbors; see [20].There is, however, sometimes an ambiguity in extracting the level energies due to their uncertain identification and possible level splittings.In all cases, absolute errors for the observables are employed, in particular for the total binding energy, because a percentage error, as used, e.g., in [9] or [11], would determine the energy of a light nucleus like 16 O much more precisely than that of a heavy nucleus like 208 Pb.
The parameters of the obtained best fits are given in Appendix A in Tables A1-A3 for the vector, scalar and mixed functional dependencies of the couplings on the densities.The quality of the fits can be assessed from the quantity χ 2 /N dof with the number of degrees of freedom N dof = N data − N para = 35.Explicit values are given in Table 2.They are much larger than one, indicating that the assumed errors (see Table 1) are estimated too small.This is particularly true for the binding energies BE and can be seen from their root-mean-square errors ∆BE that are also presented in Table 2.The results are very similar for all parametrizations without a strong preference for a particular functional form of the density dependence of the couplings.Perhaps the pure vector density dependence has a small advantage.Larger differences between the parametrizations are seen when the density dependencies of the couplings are compared for a wide range of densities.The couplings of the ω, σ, and ρ meson are depicted in Figure 1.The parametrizations of the ω meson couplings with d ω > 0 and d ω = 0 (full and dashed lines in Panel (a) of Figure 1) show a smooth decreasing trend with increasing density without a strong variation.In contrast, the functional form of Γ ω for parametrizations with d ω < 0 (dotted lines in Panel (a) of Figure 1) is very different with a minimum at densities slightly above the nuclear saturation density and a strong increase for larger densities, in particular for the vector (light and dark blue lines) and scalar (light and dark green lines) dependencies.A less strong increase is observed for a mixed dependency (orange and red lines).There is no strong influence on the ω meson coupling from the choice of the functional form of the ρ meson coupling, whether exponential or rational.The density dependencies of the σ meson couplings, depicted in Panel (b) of Figure 1, show a very similar pattern as the ω meson couplings.The main differences are the somewhat smaller absolute values.Again, the parametrizations with d σ < 0 stick out with a rising trend of the couplings at high densities.
The density dependence of the ρ meson couplings, shown in Panel (c) of Figure 1, is almost the same for all parametrizations.There is a slightly larger spread at densities above saturation, but in all cases, a decrease of the coupling with increasing density is obtained.The difference between the exponential and rational form of the functions can be recognized in the region close to zero density.

Properties of Nuclear Matter and Equation of State
Studying the EoS allows a further comparison of the different parametrizations.The nuclear matter parameters that characterize the EoS close to the saturation point and the EoS for symmetric matter and neutron matter can be examined.The density dependence of the symmetry energy is of particular interest for astrophysical applications.

Nuclear Matter Parameters
The energy per nucleon can be written as a power series: in squares of the isospin asymmetry: which depends on the difference between the neutron and proton vector densities.Here, the neutron-proton mass difference is neglected in the expansion.The first contribution in Equation ( 65) is the energy per nucleon in symmetric nuclear matter: that only depends on the total baryon density: It can be expanded close to the saturation point in powers of: measuring the deviation from the saturation density n sat .Similarly, the symmetry energy can be expanded as: with explicit terms up to second order in x.The coefficients in Equation (67) are the average nucleon mass, m nuc = (m n + m p )/2, the binding energy per nucleon at saturation, B sat , the incompressibility, K, and the skewness, Q.There is no term linear in x because the expansion is around the minimum of the energy in symmetric matter where the pressure vanishes.The coefficients in (70) are the symmetry energy at saturation, J, the slope parameter, L, and the symmetry incompressibility, K sym .All coefficients can be obtained from appropriate derivatives of the energy per nucleon E(n B , δ) with respect to n B and δ at n p = n sat /2.The six characteristic nuclear parameters together with the saturation density, n sat , are presented in Table 3 for all 18 models of the present study.In addition, the average Dirac effective mass 36), in symmetric nuclear matter and the average Landau effective mass m L av = (m L n + m L p )/2, which is related to the density of states at the Fermi energy, are given in units of the average nuclear mass in a vacuum, m nuc .Here, the definition: of the Landau effective mass with the Fermi momentum p * i is used.The data are compared to average values, including uncertainty ranges, of more than 200 existing parametrizations of the Skyrme-Hartree-Fock (SHF) and RMF models collected in [3] and [5], respectively.The scattering of the n sat and B sat values in the present RMF parametrizations is rather small, and the obtained data are consistent with the expectations from previous RMF models, but lower than the average value of SHF models.A larger spread is obtained for the incompressibilities, K, but their values are within the error bands of SHF and RMF models.A precise determination of K from fits to properties of nuclei seems to be difficult.A recent comprehensive analysis of experimental information on giant monopole resonances in [21] indicated an acceptable range of 250 MeV ≤ K ≤ 315 MeV, which is on the high side of the values from the present fits.The values for the skewness Q span a wide range with a clear correlation with the sign of the d j parameters for the ω and σ meson, i.e., the constraint on the function (59) at zero density.This fact becomes even more evident when the correlation of the incompressibility K with the skewness Q is investigated as shown in Panel (a) of Figure 2.For d ω > 0 and d σ > 0, the values of Q are the lowest of all parametrizations (diamonds), but for d ω < 0 and d σ < 0, they are the highest (squares).The sets with d ω = 0 and d σ = 0 (circles) are in between.At the same time, there is a systematic trend of larger K values with smaller d j parameters.The large positive values for Q indicate that the EoS of symmetric matter will be rather stiff at high baryon densities.From the inspection of Tables A1-A3, also a clear correlation of the mass of the σ meson with the sign of the d j parameter for j = ω, σ is found.Negative values of d j prefer to be associated with the largest σ meson masses.
The different parametrizations predict symmetry energies at saturation J within a narrow range, similar as for the B sat values.The data are close to the values expected from SHF and RMF models and inside the uncertainty band.The extracted slope parameters, L, cover a somewhat larger range that is more consistent with SHF parametrizations than the average of the RMF models.Certain trends can also be seen in the L-J correlation plot in Panel (b) of Figure 2. If models with exponential (E) and rational (R) density dependence of the ρ meson are compared separately, there is an indication that values for L are systematically larger for models with pure scalar density dependence of the couplings and smaller for models with pure vector density dependence.Models with a mixed density dependence lie in between.On the other hand, there is the trend that 'E' models have on average a larger slope parameter L than 'R' models.A large spread of the symmetry incompressibilities K sym is seen in Table 3, and even the sign of K sym cannot be determined unambiguously.Again, a systematic variation is observed as for the skewness Q or the slope parameter L. The obtained values for K sym are consistent on the whole with those of the SHF and RMF models in the compilations.The Dirac and Landau effective masses at saturation, m D av and m L av , are systematically lower as compared to the averages of the RMF and SHF models.This observation is correlated with spin-orbit splittings that are predicted on average somewhat larger than in the experiment.

Equation of State and Symmetry Energy
The different functional forms of the couplings, as depicted in Figure 1, also influence the equation of state.Here, we consider two cases for T = 0.The EoS of symmetric nuclear matter is shown in Panel (a) of Figure 3.For nucleon densities below approximately 0.2 fm −3 , all parametrizations predict energies per baryon that are practically indistinguishable.Only at higher densities can different trends be seen.The most prominent feature is the large stiffness of parametrizations with negative d ω and d σ and the strong increase of the energy per nucleon with increasing density.This behavior is expected because of the high values of the skewness parameter Q, cf.Table 3.A similar observation is made for the case of pure neutron matter, shown in Panel (b) of Figure 3.The curves for parametrizations with negative d j again stick out because they are the stiffest.All other lines are within a band that corresponds to a softer neutron matter EoS.At densities below saturation, the curves are almost identical.
The density dependence of the symmetry energy E sym is also easily extracted from the general EoS.It is depicted in Figure 4 for all parametrizations of the present work.There are only small differences between the curves for densities below approximately 0.12 fm −3 .Larger variations are found at higher densities, but a different systematics is observed as compared to the EoS of symmetric matter or of neutron matter.Parametrizations with a pure dependence of the couplings on scalar densities predict the stiffest symmetry energy on average, whereas models with a vector densities dependence display the softest symmetry energies.This observation is consistent with the ordering of the models in Panel (b) of Figure 2. Parametrizations with larger values of J and L give stiffer symmetry energies; however, this is partly compensated by smaller values of the symmetry incompressibility K sym .The results clearly show that a fit of parametrizations to properties of finite nuclei fix the EoS at sub-saturation densities fairly well, but the extrapolation to higher densities depends strongly on the functional form of the couplings and the choice of the argument, i.e., whether a scalar or vector density dependence is used.

Constraints on the Density Dependence of the Couplings
The selection of a particular density as the argument of the couplings Γ j has also consequences for the EoS under specific conditions.This can serve as a criterion to exclude some parametrizations.Two particular cases can be distinguished: a dependence of the vector meson couplings on the scalar densities and a dependence of the scalar meson couplings on the vector densities.In the following, only symmetric matter is considered to keep the discussion simple.
Usually, the scalar densities n i vanish at some finite total vector density n (v) = n ω .The couplings and thus the functions (30) in general have a smooth dependence on the density.The scalar potential (13) at zero scalar densities n (s) i is given by the rearrangement term only as: with the derivative (31) of the ω meson coupling function since n (s) = n σ = 0 in this case.If D (s) ω (0) < 0, the Dirac effective mass (36) becomes zero for m i = S i , and there is a limiting (vector) density: up to which the EoS can be calculated in the model for a particle i.This density is usually much larger than typical vector densities considered in applications of the model.However, in order to avoid the problem of a collapsing effective mass, the coupling of the ω meson should have a non-negative derivative at density zero, i.e., D ω (0) ≥ 0, if it depends on the scalar density.This condition would exclude parametrizations SPE, SNE, SPR and SNR.
Antiparticles contribute differently than particles to the vector and scalar densities (34) and (35), respectively.In the former case, their densities partly cancel, and in the latter case, they add to the total contribution of a particle species i (sum over η).For vanishing baryon densities, i.e., n (v) i = 0 for all nucleons i, the effective chemical potentials µ * i = µ i − V i have to vanish and the vector potentials ω 0 and ρ 0 in Equation ( 14) are zero, as well.However, the scalar densities n (s) i are positive at finite temperature and rise with increasing T. Hence, the scalar potential σ is finite.If the couplings of the scalar mesons depend on the total vector density n (v) , the vector potentials have a contribution from the rearrangement terms.For symmetric nuclear matter, they are given by: with the derivative (45) of the σ meson coupling function.If the derivative D (v) σ (0) is unequal zero, the chemical potential µ i = V i is finite and not zero as expected from the condition of vanishing vector densities.The mismatch between the conditions n (v) i = 0 and µ i = 0 becomes more severe with increasing temperature.This feature can be avoided by excluding the parametrizations VPE, VNE, VPR, VNR and VZE to guarantee that the derivative of the vector meson couplings are zero at vanishing total baryon density.The only the parametrization VZR is admitted.
The two problems above do not appear for the M parametrizations where the couplings of vector (scalar) mesons depend only on vector (scalar) densities.This form of modeling the effective in-medium interaction closely corresponds to the structure of most of the earlier RMF models with NL self-interactions of the mesons.In these approaches, there are no cross-terms of σ and ω contributions in the Lagrangian density.Only self-couplings of the form: and: were considered additionally in Equation ( 4), cf.[5].This leads to field equations: for the σ and ω mesons without a cross-coupling.Couplings of the form ω µ ω µ ρ ν • ρ ν , which were introduced later in some models to modify the density dependence of the isovector part of the interaction, also do not violate the vector-scalar separation.
Taking the above considerations into account, only the SZE, SZR, VZR and all M parametrizations are really viable since they do not show the identified problems.

Conclusions and Outlook
Relativistic mean-field models are widely-used phenomenological approaches to describe the properties of nuclear matter and finite nuclei.Within a subclass of these models, i.e., those with density-dependent meson-nucleon couplings, the effects of choosing different functional forms of the dependence on vector or scalar particle densities were studied going beyond the standard choice of functions and assumed vector density dependence.The model parameters were obtained in all cases by a fit to observables of a selected set of spherical nuclei.
Despite the differences of the obtained energy density functionals, the description of nuclei has practically the same quality for all parametrizations, and the equations of state of symmetric nuclear matter and pure nuclear matter below the nuclear saturation density look very similar.In contrast, differences in some of the characteristic nuclear matter parameters and the equations of state above saturation are found.This is most evident for the incompressibility K and the skewness Q that correlate with the d j parameter of the rational function used for the density dependence of the isoscalar mesons.The differences in the nuclear matter parameters are reflected in the high-density behavior of the equations of state.Similarly, there is a connection between the symmetry energy at saturation J, the slope parameter L and the choice of the argument, i.e., scalar or vector density.This also affects the stiffness of the symmetry energy.Robust constraints at densities above nuclear saturation are clearly needed to select proper parametrizations for further applications of the model.Some of the parametrizations studied in this work can lead to problems, e.g., the breakdown of the description of nuclear matter at a finite baryon density or the non-vanishing of the baryon chemical potential at finite temperature and zero baryon density.As a consequence, certain combinations of functional forms and arguments for the density dependence of the couplings have to be rejected, in particular those where the couplings of the vector mesons depend on the scalar density with negative derivative or couplings of scalar mesons depending on vector densities.
In the present study, only models with σ, ω and ρ mesons were considered.In a next step, also the δ meson should be included, which could affect particularly the density dependence of the symmetry energy at high baryon densities.Furthermore, the effect of tensor couplings of the vector mesons with the nucleons could be investigated.In this study, no mechanism for taking care of pairing effects was included in the description.All these future extensions of the model will increase the number of independent parameters, and a more extensive fitting procedure is required.Furthermore, the selection of observables and the size of their uncertainties can be reconsidered and will affect the final predictions of the models.the three types of density dependence considered in this work.The parameters a j of the rational function (59) are not specified since they are determined by Relation (63).
64) in the multidimensional parameter space.The χ 2 function is a summation of contributions from N data observables, comparing experimental data O (exp) n with the parameter-dependent model values O (model) n

Figure 1 .
Figure 1.Density dependence of the meson-nucleon coupling Γ j on the vector density n (v) or scalar density n (s) for the ω meson (a), σ meson (b) and ρ meson (c).The coding of the lines is given in the legend on the right.

Figure 2 .
Figure 2. Correlation of the incompressibility K with the skewness Q in (a) and of the symmetry energy J with the slope parameter L in (b).

Figure 3 .Figure 4 .
Figure 3. Equation of state of symmetric nuclear matter in (a) and of pure neutron matter in (b) for T = 0.The coding of the lines is the same as in Figure 1.

i
and the Dirac effective masses m D i are larger than zero.However, if the couplings depend on the total scalar density n (s) = n

Table 1 .
Selected nuclei and values of experimental observables used in the fitting procedure: binding energies per nucleon BE/A, charge radii r charge , diffraction radii r diff , surface thicknesses σ surf and spin-orbit splittings ∆ (nl) for neutron and proton levels, respectively, with principal quantum number n and orbital angular momentum l.The last line gives the assumed uncertainties.

Table 2 .
Quality of the parametrizations measured with the quantities χ 2 per number of degrees of freedom N dof = N data − N para and the root-mean-square error ∆BE of the binding energy.

Table 3 .
Nuclear matter parameters of the DD-RMF parametrizations determined in Section 4 in comparison with averages of Skyrme-Hartree-Fock (SHF) and RMF models.

Table A1 .
Parameter sets from the fitting of DD-RMF models with different vector density dependencies of the couplings to observables of nuclei.

Table A2 .
Parameter sets from the fitting of DD-RMF models with different scalar density dependencies of the couplings to observables of nuclei.

Table A3 .
Parameter sets from the fitting of DD-RMF models with different mixed density dependencies of the couplings to observables of nuclei.