## Improving relativistic energy density functionals with

## tensor couplings

Stefan Typel 1,2,3* and Shalom Shlomo 4

1*

Texas A&amp;M University, Cyclotron Institute, College Station, 77843, Texas, USA.

Technische Universit¨t Darmstadt, Fachbereich Physik, Institut f¨r Kernphysik, a u Schlossgartenstraße 9, Darmstadt, 64289, Germany. 2 GSI Helmholtzzentrum f¨r Schwerionenforschung, Theorie, Planckstraße 1, u Darmstadt, 64291, Germany. 3 Helmholtz Forschungsakademie Hessen f¨r FAIR (HFHF), GSI Helmholtzzentrum u f¨r Schwerionenforschung, Campus Darmstadt, Schlossgartenstraße 2, Darmstadt, u 64289, Germany. 4

*Corresponding author(s). E-mail(s): stypel@ikp.tu-darmstadt.de, ORCID: 0000-0003-3238-9973; Contributing authors: s-shlomo@tamu.edu,

ORCID: 0000-0002-4049-7563;

## Abstract

Energy density functionals (EDFs) have been used extensively with great success to calculate properties of nuclei and to predict the equation of state (EOS) of dense nuclear matter. Besides non-relativistic EDFs, mostly of the Skyrme or Gogny type, relativistic EDFs of different types are in widespread use. In these latter approaches, the effective in-medium interaction is described by an exchange of mesons between nucleons. In most cases, only minimal meson-nucleon couplings are considered. The effects of additional tensor couplings were rarely investigated. In this work, a new relativistic EDF with tensor couplings and density dependent minimal meson-nucleon couplings will be presented. The parameters of the model are determined using a carefully selected set of experimental data with realistic uncertainties that are determined self-consistently. Predictions for various nuclear observables, the nuclear matter equation of state, and properties of neutron stars are discussed.

Keywords: Relativistic energy density functional, tensor couplings, parameterisation, nuclear properties, nuclear matter, equation of state, neutron stars

## 1 Introduction

For more than half a century, energy density functionals (EDFs) have been used with great success in the description of atomic nuclei across the nuclear chart and nuclear matter over a large range of densities and isospin asymmetries. Many improvements were introduced in these phenomenological

models leading to a more and more precise prediction of nuclear observables but, nevertheless, the quest for further improvements remains. There are different strategies to develop more accurate EDFs. They either follow a systematic approach, e.g., using expansions in power of small quantities, or try to incorporate known effects in a more heuristic way.

EDFs exist in different versions. Most well known are non-relativistic approaches that employ convenient representations of effective in-medium interactions between nucleons, e.g., the Skyrme or Gogny type models [1-4]. Relativistic or covariant EDFs start from a field-theoretic approach where the interaction is described by meson-exchange taking medium dependent effects into account in different ways [5-15]. Either self-interactions of the mesons are introduced in so-called nonlinear (NL) models or a density dependence (DD) of the the meson-nucleon couplings is utilized. Another class of relativistic EDFs uses only nucleons as degrees of freedom with point couplings and strengths that are constants or depend on the density, see, e.g., Refs. [16-19]. No explicit meson contributions appear in these formulations. These extensions of the original Walecka model [20] are essential to obtain reasonable properties of nuclei and nuclear matter. Applying certain approximations to the full many-body wave function of the system of interest, e.g., of Hartree, Hartree-Fock or Hartree-Fock-Bogoliubov type, a set of coupled equations is found which has to be solved self-consistently.

In the present work, new parameterisations of a relativistic EDF will be presented that include tensor couplings between nucleons and mesons in addition to the standard minimal couplings. Such an extension has been considered before in early applications of relativistic EDFs, see, e.g., Ref. [6], but it has been explored scarcely in depth. Effects of considering ω -tensor couplings in EDFs on the spin-orbit splittings and indirectly on the nucleon effective mass have been studied in Ref. [21]. A relativistic EDF with ω and ρ -tensor couplings in the framework of the NL coupling models has been developed in Ref. [22] taking into account also contributions from Dirac-sea nucleons. The parameters of this model were adjusted to the energies of a eight nuclei from 16 O to 208 Pb. Including ρ -tensor, but no ω -tensor, contributions in relativistic Hartree-Fock calculations showed that the single-nucleon spectrum could be improved [23]. Tensor couplings were also included in the framework of point couplings models in the Hartree-Fock approximation, see, e.g., Ref. [24]. The interest in relativistic EDFs with tensor couplings has been grown recently as it can affect the description of surface properties of nuclei and (indirectly) the equation of state of nuclear matter. New parametrisations were developed in Ref. [25] for relativistic models with meson-nucleon couplings that depend on vector or scalar densities of the medium in various combinations, but using a more limited number of nuclear observables in the parameter fit as compared to the present work. A new parameterisation of a covariant EDF with density dependent minimal couplings and density

independent tensor couplings in the Hartree-Bogoliubov approximation was developed recently [26]. The model parameters were determined using constraints from nuclear matter and nuclei with preset uncertainties of the fitted observables. The effects of tensor couplings on different properties of nuclei and, using the resulting equation of state (EOS) of dense matter, on neutron stars were explored in Ref. [27], however without a thorough determination of model parameters. Here, we shall embark on the development of a new relativistic EDF with minimal and tensor couplings using a careful selection of nuclear observables that are used in the fit of the model parameters. Particular attention is paid to the choice of the uncertainties of the considered observables. The medium dependence of the interaction is taken into account with a density dependence of the minimal couplings using the functional forms introduced in Ref. [11]. In addition to the description of nuclei, the implications of the new parameterisations on the equation of state of nuclear matter and the properties of neutron stars will be explored.

The content of this paper is structured as follows: The formalism of the relativistic energy density functional is presented in Section 2. The main quantities are introduced and the notation is established, in particular regarding the contribution with the tensor couplings. This section also provides the relevant equations for the description of spherical nuclei and nuclear matter in two subsections. In Section 3 the approach to determine the model parameters is outlined. The functional form of the density dependent couplings is specified and the selection of the observables and nuclei in the fit is presented. Details of the fitting strategy are also given. Results and their discussion follow in Section 4, where properties of nuclei, nuclear matter, and neutron stars are considered in individual subsections. Finally, we close with a summary and outlook in Section 5.

## 2 Relativistic energy density functional

In a relativistic theory of nuclei and nuclear matter, nucleons and mesons are assumed to be the relevant degrees of freedom. Protons and neutrons are represented as Dirac spinors Ψ η (with η = π for protons and η = ν for neutrons) following Fermi statistics, whereas mesons are described by bosonic fields. The strong interaction between nucleons can be represented effectively by an exchange of mesons of different type. They often share their properties with those mesons of the same name and quantum numbers that are observed in experiments but are not necessarily identical. Actually, it is possible to describe the main features of the strong interaction by a rather small set of effective mesons that reduces the complexity of the calculations and the number of model parameters.

In relativistic EDFs, it is common practise to introduce an iso-scalar, Lorentz-scalar σ meson to model the attraction and an iso-scalar, Lorentz-vector ω meson to model the repulsion between

nucleons. In order to take the isovector dependence of the interaction into account, an iso-vector, Lorentz-scalar δ meson and an iso-vector, Lorentz-vector ρ meson are also considered. The notation of these fields follows their properties, i.e., σ , ω µ , /vector δ , and /vector ρ µ with Lorentz indices µ and the vector symbol for the three isovector components. The pion as another important meson is not included as it does not appear in the final equations after applying the basic approximations of the approach. Besides the above mentioned mesons, the photon field, denoted by the vector A µ , completes the number of degrees of freedom to include also the electromagnetic interaction in the model.

The derivation of an EDF in a relativistic approach follows several steps that are detailed below. First, a Lagrangian density is formulated containing the contributions for the free propagation of the fields and their interaction. From it, the coupled field equations are derived that have to be solved self-consistently. This can be achieved in certain approximations and expressions for the energy of the system and other quantities of interest can be obtained. They are functions of the particle fields and their densities, so that the EDF is a generalized form of a density functional. It depends on a rather small number of parameters that need to be determined by appropriate procedures. In the formulas given below, we use a system of units where /planckover2pi1 = c = 1. The value /planckover2pi1 = 197 3269804 MeV fm/c is used for the conversion of units. .

## 2.1 Lagrangian density, field equations, and energy density

The development of the relativistic EDF starts from a covariant Lagrangian density

$$\mathcal { L } & = \sum _ { \eta = \pi, \nu } \overline { \Psi } _ { \eta } \left ( \gamma _ { \mu } i D ^ { \mu } - \sigma _ { \mu \lambda } \mathcal { T } ^ { \mu \lambda } - \mathcal { M } _ { \eta } \right ) \Psi _ { \eta } \\ & \quad + \frac { 1 } { 2 } \left ( \partial ^ { \mu } \sigma \partial _ { \mu } \sigma - M _ { \sigma } ^ { 2 } \, \sigma ^ { 2 } + \partial ^ { \mu } \delta \cdot \partial _ { \mu } \vec { \delta } - M _ { \delta } ^ { 2 } \, \vec { \delta } \cdot \vec { \delta } \\ & \quad - \frac { 1 } { 2 } G ^ { \mu \lambda } G _ { \mu \lambda } + M _ { \omega } ^ { 2 } \, \omega ^ { \mu } \omega _ { \mu } - \frac { 1 } { 2 } \, \vec { H } ^ { \mu \lambda } \cdot \vec { H } _ { \mu \lambda } + M _ { \rho } ^ { 2 } \, \vec { \rho } ^ { \mu } \cdot \vec { \rho } _ { \mu } - \frac { 1 } { 2 } F ^ { \mu \lambda } F _ { \mu \lambda } \right ) \,, \\ & \quad \dots \quad \dots \quad \dots$$

with several contributions. Here, the field tensors of the mesons and of the massless electromagnetic field are given by the usual expressions;

$$G ^ { \mu \lambda } = \partial ^ { \mu } \omega ^ { \lambda } - \partial ^ { \lambda } \omega ^ { \mu } \,, \quad \vec { H } ^ { \mu \lambda } = \partial ^ { \mu } \vec { \rho } ^ { \lambda } - \partial ^ { \lambda } \vec { \rho } ^ { \mu } \,, \quad F ^ { \mu \lambda } = \partial ^ { \mu } A ^ { \lambda } - \partial ^ { \lambda } A ^ { \mu } \,. \quad \quad ( 2 )$$

The minimal coupling to the nucleons is realized by the corresponding terms in the covariant derivative

$$i \mathcal { D } ^ { \mu } = i \partial ^ { \mu } - \tilde { \Gamma } _ { \omega } \omega ^ { \mu } - \tilde { \Gamma } _ { \rho } \vec { \rho } ^ { \mu } \cdot \vec { \tau } - \Gamma _ { \gamma } A ^ { \mu } \frac { 1 + \tau _ { 3 } } { 2 } \,,$$

for the Lorentz-vector mesons and the photon, and in the mass operator

$$\mathcal { M } _ { \eta } = M _ { \eta } - \tilde { \Gamma } _ { \sigma } \sigma - \tilde { \Gamma } _ { \delta } \vec { \delta } \cdot \vec { \tau } \,,$$

for the Lorentz-scalar mesons, respectively. The term in Eq. (1) with the quantity

$$\mathcal { T } ^ { \mu \lambda } = \frac { \Gamma _ { T \omega } } { 2 M _ { \pi } } G ^ { \mu \lambda } + \frac { \Gamma _ { T \rho } } { 2 M _ { \pi } } \vec { H } ^ { \mu \lambda } \cdot \vec { \tau } \,,$$

accounts for the additional tensor coupling of the nucleons to the mesons. Tensor couplings to the electromagnetic field are not included. The masses of the nucleons and the mesons are denoted by M η ( η = π, ν ) and M j ( j = σ, δ, ω, ρ ), respectively. The isospin dependence is expressed with help of the isopin matrices τ k ( k = 1 2 3) in analogy to the Pauli spin matrices , , σ k . The coupling factors ˜ Γ j in Eqs. (3), and (4) are assumed to be functionals of the nucleon fields Ψ η , and Ψ η = Ψ † η γ 0 . In contrast, the tensor couplings Γ Tω , and Γ Tρ , in Eq. (5) as well as the electromagnetic coupling Γ γ in Eq. (3) are simple constants. The devision by the proton mass in Eq. (5) is introduced to have dimensionless couplings Γ Tω and Γ Tρ .

The field equations of the nucleons, mesons, and the photon can be derived in the usual way from the Euler-Lagrange equations. However, they cannot be solved in the field-theoretical sense and approximations are needed to arrive at tractable equations. The fields of the mesons and the photons are not quantized but considered as classical fields and will be denoted by the same symbols as in the original Lagrangian density (1) in the following. For the nucleons, the Hartree approximation is applied, i.e., the many-nucleon wave function Ψ of the system of interest is written as a product of single-nucleon states Ψ ηk η m η ,

$$\Psi = \prod _ { \eta = \pi, \, \nu } \prod _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,,$$

with occupation factors w ηk η m η . Only values of w ηk η m η = 0 and 1 are allowed. This condition reflects the fermionic nature of the nucleons. The meaning of the quantum numbers k η , and m η , is explained in Subsection 2.2. Only single-nucleon states of positive-energy are considered here, corresponding to the no-sea approximation and only a finite number of states contributes in the many-nucleon wave function.

For the dependence of the couplings Γ ˜ j in Eqs. (3), and (4), a so-called vector density dependence is assumed, i.e., with the approximation mentioned above, they become functions Γ j of the

vector density, with the total current,

$$\varrho _ { v } & = \sqrt { j _ { \mu } j ^ { \mu } } \,,$$

$$j ^ { \mu } & = \sum _ { \eta = \pi, \nu } \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } \gamma ^ { \mu } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,, \\ \intertext {. The countries } \alpha \text{ is in fact $aT$ around $\alpha$ but has to be distinctinhead from the}$$

of the nucleons. The quantity /rho1 v is in fact a Lorentz scalar but has to be distinguished from the actual scalar density

$$\varrho _ { s } & = \sum _ { \eta = \pi, \nu } \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,, \\ \text{slims all the code} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \,\text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and}, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \, \text{end} \,$$

that is used for couplings with a scalar density dependence, c.f. Ref. [25].

In the actual application to nuclei and nuclear matter, a specific frame of reference is chosen and the covariance of the equations will be lost. In particular, a frame is chosen where the system is at rest, the couplings depend only on the zero-component of the current /rho1 v = j 0 , and /vector j = 0. Since only static solutions without a change of isospin are considered, only the zero-component of the Lorentz-vector fields and a single component of the isovector quantities have to be taken into account, since all other components vanish. Assuming spherical symmetry, the field equation of the nucleons assumes the form of a modified Dirac equation,

$$\hat { H } _ { \eta } \Psi _ { \eta k _ { \eta } m _ { \eta } } ( \vec { r } ) = \left [ \vec { \alpha } \cdot \hat { \vec { p } } + \beta \left ( M _ { \eta } - S _ { \eta } \right ) + V _ { \eta } + i \vec { \gamma } \cdot \frac { \vec { r } } { r } T _ { \eta } \right ] \Psi _ { \eta k _ { \eta } m _ { \eta } } ( \vec { r } ) = E _ { \eta k _ { \eta } } \Psi _ { \eta k _ { \eta } m _ { \eta } } ( \vec { r } ) \,, \quad ( 1 0 )$$

with /vector α = γ /vector 0 γ , and β = γ 0 , for single-nucleon states of energy E ηk η . The interaction with the meson and photon fields is expressed through the scalar potentials

$$S _ { \eta } ( \vec { r } ) = \Gamma _ { \sigma } \sigma + g _ { \eta } \Gamma _ { \delta } \delta \,,$$

vector potentials, and tensor potentials,

$$T _ { \eta } ( \vec { r } ) = - \frac { \Gamma _ { T \omega } } { M _ { \pi } } \frac { \vec { r } } { r } \cdot \vec { \nabla } \omega _ { 0 } - g _ { \eta } \frac { \Gamma _ { T \rho } } { M _ { \pi } } \frac { \vec { r } } { r } \cdot \vec { \nabla } \rho _ { 0 } \,,$$

with factors g π = 1, and g ν = -1. The vector density dependence of the couplings generates the rearrangement contribution

$$V ^ { ( R ) } ( \vec { r } ) = \frac { d \Gamma _ { \omega } } { d n ^ { ( v ) } } n _ { \omega } \omega _ { 0 } + \frac { d \Gamma _ { \rho } } { d n ^ { ( v ) } } n _ { \rho } \rho _ { 0 } - \frac { d \Gamma _ { \sigma } } { d n ^ { ( v ) } } n _ { \sigma } \sigma - \frac { d \Gamma _ { \delta } } { d n ^ { ( v ) } } n _ { \sigma } \delta \,,$$

$$V _ { \eta } ( \vec { r } ) = \Gamma _ { \omega } \omega _ { 0 } + g _ { \eta } \Gamma _ { \rho } \rho _ { 0 } + \frac { 1 + g _ { \eta } } { 2 } \Gamma _ { \gamma } A _ { 0 } + V ^ { ( R ) } \,,$$

in the vector potential (12). It contains derivatives of the couplings with respect to the vector density n ( v ) = /rho1 v , as well as meson fields and source densities n j ( j = σ, δ, ω, ρ ) that also appear in the meson field equations;

$$- \Delta \sigma + M _ { \sigma } ^ { 2 } \sigma \, = \, \Gamma _ { \sigma } n _ { \sigma } \,,$$

$$- \Delta \delta + M _ { \delta } ^ { 2 } \delta = \Gamma _ { \delta } n _ { \delta } \,,$$

$$- \Delta \omega _ { 0 } + M _ { \omega } ^ { 2 } \omega _ { 0 } = \Gamma _ { \omega } n _ { \omega } + \frac { \Gamma _ { T \omega } } { m _ { n } } \vec { \nabla } \cdot \vec { \jmath } _ { \omega } ^ { ( t ) } \,,$$

$$\omega n _ { \omega } + \frac { \mathfrak { I } \, T \omega } { m _ { p } } \vec { \nabla } \cdot \vec { \mathfrak { J } } _ { \omega } ^ { ( t ) } \,,$$

$$- \Delta \rho _ { 0 } + M _ { \rho } ^ { 2 } \rho _ { 0 } = \Gamma _ { \rho } n _ { \rho } + \frac { \Gamma _ { T \rho } ^ { \cdot } } { m _ { p } } \vec { \nabla } \cdot \vec { \jmath } _ { \rho } ^ { ( t ) } \,,$$

for the non-vanishing components of the fields. They are given by

$$n _ { \sigma } = n ^ { ( s ) } = n _ { \pi } ^ { ( s ) } + n _ { \nu } ^ { ( s ) } \,, \quad n _ { \delta } = n _ { \pi } ^ { ( s ) } - n _ { \nu } ^ { ( s ) } \,,$$

$$n _ { \omega } = n ^ { ( v ) } = n _ { \pi } ^ { ( v ) } + n _ { \nu } ^ { ( v ) } \,, \quad n _ { \rho } = n _ { \pi } ^ { ( v ) } - n _ { \nu } ^ { ( v ) } \,,$$

with the vector and scalar densities;

$$n _ { \eta } ^ { ( v ) } ( \vec { r } ) = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( v ) } = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } \gamma ^ { 0 } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,,$$

and

$$n _ { \eta } ^ { ( s ) } ( \vec { r } ) & = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( s ) } = \sum _ { k _ { \eta } m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,, \\ \intertext { m o n c o n t i o l, t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e t h e n t i o n, e t h e n t i o n, e t h e n t i o n, e t h e n t i o n }$$

of nucleons, respectively. In Eqs. (17), and (18), also the divergences of the tensor currents

$$\vec { \lambda } _ { \omega } ^ { ( t ) } = \vec { \lambda } _ { \pi } ^ { ( t ) } + \vec { \lambda } _ { \nu } ^ { ( t ) } \,, \quad \vec { \lambda } _ { \rho } ^ { ( t ) } = \vec { \lambda } _ { \pi } ^ { ( t ) } - \vec { \lambda } _ { \nu } ^ { ( t ) } \,,$$

with,

$$\vec { j } _ { \eta } ^ { ( t ) } ( \vec { r } ) = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \vec { j } _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } i \vec { \alpha } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,, \\ \intertext { w i h t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t to t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o to t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t i t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t st t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t a t t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o ot o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t ot o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o o t t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t t o o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o o t to t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o o t o o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o to t o o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t ot ot o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t ot o o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o to o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o t o { \ } \, \quad \,$$

contribute to the source terms of the Lorentz-vector meson fields in Eqs. (17), and (18). Note that these tensor currents are not the vector part of the usual four-currents,

$$j _ { \eta } ^ { \mu } ( \vec { r } ) = \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } \overline { \Psi } _ { \eta k _ { \eta } m _ { \eta } } \gamma ^ { \mu } \Psi _ { \eta k _ { \eta } m _ { \eta } } \,.$$

Finally, the field equation for the electromagnetic field is given by the usual Poisson's equation

$$- \Delta A _ { 0 } = \Gamma _ { \gamma } n _ { \gamma } \,,$$

with the source density n γ = n ( v ) π .

In the application to nuclear matter, there is no spatial dependence of the source densities and meson fields. They become constants and there are no contributions from the tensor currents in Eqs. (17), and (18). The field equations are solved trivially and the tensor potential (13) in the Dirac equation (10) vanishes too. Nucleon states are then described by plane waves.

The energy density of the system can be obtained as a component of the energy-momentum tensor

$$T ^ { \mu \lambda } & = \sum _ { \varphi } \frac { \partial \mathcal { L } } { \partial ( \partial _ { \mu } \varphi ) } \partial ^ { \lambda } \varphi - g ^ { \mu \lambda } \mathcal { L } \,, \\ \intertext { v o r all f i d e t a } & = \mathfrak { W } \,. \quad \overline { \mathfrak { W } } \,. \quad \sigma \, \delta \,. \quad \sigma \, \delta \,. \quad \partial _ { \lambda } \, \text{ and the metric tensor} \, \text{is}$$

where the sum runs over all fields ϕ = Ψ ηk η m η , Ψ ηk η m η , σ, δ, ω 0 , ρ 0 , A 0 , and the metric tensor is g µλ = diag( -1 1 1 1). It can be expressed as a sum, , , ,

$$\varepsilon ( \vec { r } ) = T ^ { 0 0 } = \varepsilon _ { \text{mucleon} } ( \vec { r } ) + \varepsilon _ { \text{field} } ( \vec { r } ) \,,$$

with the contribution

$$\varepsilon _ { \text{nucleon} } ( \vec { r } ) = \sum _ { \eta = \pi, \nu } \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } E _ { \eta k _ { \eta } } n ^ { ( v ) } _ { \eta k _ { \eta } m _ { \eta } } \,, \\ \text{and to the field}$$

of the nucleons and that of the fields

$$\varepsilon _ { \text{field} } ( \vec { r } ) = \frac { 1 } { 2 } \left ( \vec { \nabla } \sigma \cdot \vec { \nabla } \sigma + M _ { \sigma } ^ { 2 } \sigma ^ { 2 } + \vec { \nabla } \delta \cdot \vec { \nabla } \delta + M _ { \delta } ^ { 2 } \delta ^ { 2 } \\ - \vec { \nabla } \omega _ { 0 } \cdot \vec { \nabla } \omega _ { 0 } - M _ { \omega } ^ { 2 } \omega _ { 0 } ^ { 2 } - \vec { \nabla } \rho _ { 0 } \cdot \vec { \nabla } \rho _ { 0 } - M _ { \rho } ^ { 2 } \rho _ { 0 } ^ { 2 } - \vec { \nabla } A _ { 0 } \cdot \vec { \nabla } A _ { 0 } \right ) - V ^ { ( R ) } n ^ { ( v ) } \quad ( 3 0 )$$

with an explicit rearrangement term.

## 2.2 Spherical nuclei

For describing nuclei, the coupled differential equations (10) for the single-nucleon wave functions, and (15) - (26), have to be solved self-consistently. In the case of spherical symmetry, it is convenient to represent the single-nucleon wave functions as

$$\Psi _ { \eta k _ { \eta } m _ { \eta } } ( \vec { r } ) = \frac { 1 } { r } \begin{pmatrix} F _ { \eta k _ { \eta } } ( r ) \mathcal { Y } _ { \kappa _ { \eta } m _ { \eta } } ( \hat { r } ) \\ i G _ { \eta k _ { \eta } } ( r ) \mathcal { Y } _ { - \kappa _ { \eta } m _ { \eta } } ( \hat { r } ) \end{pmatrix} \,,$$

with radial functions F ηk η , and G ηk η , and tensor spherical harmonics,

$$\mathcal { Y } _ { \kappa _ { \eta } m } ( \hat { r } ) = \sum _ { m _ { l } m _ { s } } ( l _ { \eta } \, m _ { l } \, s \, m _ { s } | j _ { \eta } \, m _ { \eta } ) Y _ { l _ { \eta } m _ { l } } ( \hat { r } ) \chi _ { s m _ { s } } \,,$$

in the upper and lower components of the Dirac spinor. Each state is identified by the index η = π, ν , and two quantum numbers n η , and κ η , that are collectively denoted as k η , and a total angular momentum projection m η . The first, n η = 0 1 2 , , , . . . , is the principal quantum number related to the number of nodes of the function F ηk η ; the second, κ η = ± 1 , ± 2 , . . . , defines the total angular momentum j η = | κ η | -1 / 2 and the orbital angular momentum l η = j η -κ / η (2 | κ η | ) of a state. Then the Dirac equation (10) is equivalent to the set of two coupled differential equations;

$$\left ( E _ { \eta k _ { \eta } } - V _ { \eta } - M _ { \eta } + S _ { \eta } \right ) F _ { \eta k _ { \eta } } + \left ( \frac { d } { d r } + \frac { \kappa _ { \eta } } { r } - T _ { \eta } \right ) G _ { \eta k _ { \eta } } = 0 \,,$$

and,

$$\left ( \frac { d } { d r } - \frac { \kappa _ { \eta } } { r } + T _ { \eta } \right ) F _ { \eta k _ { \eta } } - \left ( E _ { \eta k _ { \eta } } - V _ { \eta } + M _ { \eta } - S _ { \eta } \right ) G _ { \eta k _ { \eta } } = 0 \,, \\ \i e \, \text{radial functions that have to be solved with the boundary conditions } F _ { \eta k _ { \eta } } ( 0 ) = G _ { \eta k _ { \eta } } ( 0 ) = 0,$$

for the radial functions that have to be solved with the boundary conditions F ηk η (0) = G ηk η (0) = 0, and lim r → 0 F ηk η ( r ) = lim r → 0 G ηk η ( r ) = 0.

The representation (31) of the single-nucleon wavefunctions leads to the expressions

$$n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( s ) } ( r ) = \frac { 1 } { 4 \pi r ^ { 2 } } \left [ \left | F _ { \eta k _ { \eta } } ( r ) \right | ^ { 2 } - \left | G _ { \eta k _ { \eta } } ( r ) \right | ^ { 2 } \right ]$$

and

$$n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( v ) } ( r ) = \frac { 1 } { 4 \pi r ^ { 2 } } \left [ \left | F _ { \eta k _ { \eta } } ( r ) \right | ^ { 2 } + \left | G _ { \eta k _ { \eta } } ( r ) \right | ^ { 2 } \right ] & & ( 3 6 ) \\ \text{and vector single-nucleon densities in Eqs. (22), and (21), respectively, if an average}$$

for the scalar and vector single-nucleon densities in Eqs. (22), and (21), respectively, if an average over the quantum number m η is applied. This procedure is used to obtain spherical source densities in the field equations for the mesons and the photon in case that sub-shells of given j are not completely filled. The vector densities are normalized as,

$$\int d ^ { 3 } r \, n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( v ) } ( r ) = 1 \,,$$

but, there is no such condition for the scalar densities. The single-nucleon tensor currents have only a component in radial direction in a spherical system. Then one can define

$$n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } = \frac { 1 } { 4 \pi r ^ { 2 } } \left [ F _ { \eta k _ { \eta } } ^ { * } ( r ) G _ { \eta k _ { \eta } } ( r ) + F _ { \eta k _ { \eta } } ( r ) G _ { \eta k _ { \eta } } ^ { * } ( r ) \right ] \,,$$

so that,

$$\vec { \jmath } _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } = n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } \frac { \vec { r } } { r } \,, \quad \text{and} \quad \vec { \nabla } \cdot \vec { \jmath } _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } = \frac { 1 } { r ^ { 2 } } \frac { d } { d r } \left ( r ^ { 2 } n _ { \eta k _ { \eta } m _ { \eta } } ^ { ( t ) } \right ) \,, \\ \intertext { F a. ( 24 ) after anliving the averaging nrocedure. The occimation factors in the calculation of }$$

in Eq. (24) after applying the averaging procedure. The occupation factors in the calculation of the densities are determined such that,

$$\sum _ { k _ { \nu }, m _ { \pi } } w _ { \pi k _ { \pi } m _ { \pi } } = Z \,, \quad \text{and} \quad \sum _ { k _ { \nu } m _ { \nu } } w _ { \nu k _ { \nu } m _ { \nu } } = N \,,$$

with the proton and neutron numbers Z and N of the nucleus, respectively, with mass number A = Z + N .

The energy of the nucleus is finally obtained as

$$E _ { \text{nucleus} } = E _ { \text{sp} } + E _ { \text{field} } - E _ { \text{corr} } \,,$$

with the contributions of the single-particle states

$$E _ { \text{sp} } = \sum _ { \eta = \pi, \nu } \sum _ { k _ { \eta }, m _ { \eta } } w _ { \eta k _ { \eta } m _ { \eta } } E _ { \eta k _ { \eta } } \,,$$

and of the fields,

$$E _ { \text{field} } = \frac { 1 } { 2 } \int d ^ { 3 } r \left \{ \Gamma _ { \sigma } n _ { \sigma } \sigma + \Gamma _ { \delta } n _ { \delta } \delta - \Gamma _ { \gamma } n _ { \gamma } A _ { 0 } \\ & - \left [ \Gamma _ { \omega } n _ { \omega } + \frac { \Gamma _ { T \omega } } { M _ { p } } \vec { \nabla } \cdot \vec { \jmath } _ { T \omega } \right ] \omega _ { 0 } - \left [ \Gamma _ { \rho } n _ { \rho } + \frac { \Gamma _ { T \rho } } { M _ { p } } \vec { \nabla } \cdot \vec { \jmath } _ { T \rho } \right ] \rho _ { 0 } - V ^ { ( R ) } n ^ { ( v ) } \right \},$$

where the field equations were exploited after a partial integration using Eq. (30). Since the translational symmetry is broken, a correction for the center-of-mass motion has to be applied. It is calculated as the expectation value

$$E _ { \text{corr} } = \langle \Psi | \frac { \vec { P } ^ { 2 } } { 2 M ( Z, N ) } | \Psi \rangle \,,$$

of the total momentum

$$\vec { P } = \sum _ { j = 1, Z } \vec { p } _ { \pi j } + \sum _ { j = 1, N } \vec { p } _ { \nu j } \,, \\ \text{with to the total $\max_{\mathbb{ }CZ}$} \, M / Z \, N \, - \, Z \, M \, \perp \, N \, M \, \text{of the balance non-totime} }$$

in non-relativistic form with the total mass M Z,N ( ) = ZM π + NM ν of the nucleus, neglecting the binding energy.

For the numerical solution of the coupled differential equations (33), and (34), the Lagrangemesh technique is utilized, see Ref. [28] and references therein. This approach combines the virtues of an expansion in terms of basis functions and a discretization on a grid. At the same time it is numerically efficient and precise. The radial wave functions are represented with their values at grid points r i = hx ( α ) i , where h is a scale parameter, and x ( α ) i are the zeroes i = 1 , . . . , N zero of the Laguerre polynomials L ( α ) N ( x ) with N zero = 25, and α = 2 | κ | . The scale parameter h is chosen such that the largest zero of L (0) N is at a radius r = 0 95 . R , with R = r 0 A 1 / 3 +10 fm, and r 0 = 1 2 fm . depending on the mass number A of the nucleus. The field equations of the mesons, (15) - (18), and the photon, (26), are solved by expanding the fields and sources in terms of Riccati-Bessel functions. Special care has to be taken for the electromagnetic potential A 0 to obtain the correct asymtotic behavior for large radii. For more details on the numerical techniques, see Ref. [28].

## 2.3 Nuclear matter

The solution of the field equations (10), (15) - (18), and (26) simplifies considerably in homogeneous and isotropic nuclear matter. The source densities and meson fields become constants. The electromagnetic field has to be discarded and all contributions of the tensor terms vanish. Nucleons are described by plane waves of momentum k η following the dispersion relation

$$E _ { \eta k _ { \eta } } = \sqrt { k _ { \eta } ^ { 2 } + \left ( M _ { \eta } ^ { * } \right ) ^ { 2 } } + V _ { \eta } \,,$$

where the Dirac effective mass

$$M _ { \eta } ^ { * } = M _ { \eta } - S _ { \eta } \,,$$

and constant scalar and vector potentials (11), and (12), appear. At zero temperature, all relevant quantities can be calculated analytically. The scalar and vector densities assume the form;

$$n _ { \eta } ^ { ( s ) } = \gamma _ { \eta } \int _ { 0 } ^ { k _ { \eta } ^ { * } } \frac { d ^ { 3 } k _ { \eta } } { ( 2 \pi ) ^ { 3 } } \, \frac { M _ { \eta } ^ { * } } { \sqrt { ( k _ { \eta } ) ^ { 2 } + ( M _ { \eta } ^ { * } ) ^ { 2 } } } = \frac { \gamma _ { \eta } M _ { \eta } ^ { * } } { 4 \pi ^ { 2 } } \left [ k _ { \eta } ^ { * } \mu _ { \eta } ^ { * } - ( M _ { \eta } ^ { * } ) ^ { 2 } \ln \frac { k _ { \eta } ^ { * } + \mu _ { \eta } ^ { * } } { M _ { \eta } ^ { * } } \right ] \,, \quad \quad ( 4 8 )$$

and,

$$n _ { \eta } ^ { ( v ) } = \gamma _ { \eta } \int _ { 0 } ^ { k _ { \eta } ^ { * } } \frac { d ^ { 3 } k _ { \eta } } { ( 2 \pi ) ^ { 3 } } = \frac { g _ { \eta } } { 6 \pi ^ { 2 } } \left ( k _ { \eta } ^ { * } \right ) ^ { 3 } \,,$$

in the continuum approximation, with the spin degeneracy factor γ η = 2. The Fermi momenta k ∗ η are related to the effective chemical potentials

$$\mu _ { \eta } ^ { * } = \mu _ { \eta } - V _ { \eta } = \sqrt { \left ( k _ { \eta } ^ { * } \right ) ^ { 2 } + \left ( M _ { \eta } ^ { * } \right ) ^ { 2 } } \,,$$

of the nucleons.

The energy density of cold nuclear matter can be written as,

$$\varepsilon = \sum _ { \eta = \pi, \nu } \frac { 1 } { 4 } \left [ 3 \mu _ { \eta } ^ { * } n _ { \eta } ^ { ( v ) } + M _ { \eta } ^ { * } n _ { \eta } ^ { ( s ) } \right ] + \frac { 1 } { 2 } \sum _ { j = \sigma, \delta, \omega, \rho } C _ { j } n _ { j } ^ { 2 } \,,$$

with the coupling-to-mass ratios

$$C _ { j } = \frac { \Gamma _ { j } ^ { 2 } } { M _ { j } ^ { 2 } } \,.$$

The pressure

$$P = \sum _ { \eta = \pi, \nu } \frac { 1 } { 4 } \left [ \mu _ { \eta } ^ { * } n _ { \eta } ^ { ( v ) } - M _ { \eta } ^ { * } n _ { \eta } ^ { ( s ) } \right ] - \frac { 1 } { 2 } \sum _ { j = \sigma, \delta } \left ( C _ { j } + n ^ { ( v ) } C _ { j } ^ { \prime } \right ) n _ { j } ^ { 2 } + \frac { 1 } { 2 } \sum _ { j = \omega, \rho } \left ( C _ { j } + n ^ { ( v ) } C _ { j } ^ { \prime } \right ) n _ { j } ^ { 2 },$$

also depends on the derivatives

$$C _ { j } ^ { \prime } & = \frac { d C _ { j } } { d n ^ { ( v ) } } \,,$$

and can be derived using the thermodynamic definition P = d ε/n ( ( v ) ) /dn ( v ) , or from the energymomentum tensor (27), as P = (1 / 3) ∑ 3 µ =1 T µµ . It is easily shown that energy density and pressure fulfill the Hugenberg-van-Hove relation [29]

$$\varepsilon + P = \sum _ { \eta = \pi, \nu } \mu _ { \eta } n _ { \eta } ^ { ( v ) } \,,$$

since the rearrangement terms in the pressure combine with those in the chemical potentials µ η appearing in Eq. (50).

The energy density (51), and pressure (53), are functions of two independent densities, the proton and neutron vector densities n ( v ) π , and n ( v ) ν . Alternatively, they can be regarded as functions of the baryon density

$$n _ { b } = \varrho _ { v } = n ^ { ( v ) } = n _ { \nu } ^ { ( v ) } + n _ { \pi } ^ { ( v ) } \,,$$

and the neutron-proton asymmetry

$$\mathfrak { m } _ { J } \quad & \alpha = \frac { n _ { \nu } ^ { ( v ) } - n _ { \pi } ^ { ( v ) } } { n _ { b } } \,,$$

that lies inside the interval [ -1 +1]. Then it is convenient to define the energy per nucleon ,

$$\frac { E } { A } ( n _ { b }, \alpha ) = \frac { \varepsilon } { n _ { b } } - M _ { \text{nuc} } \,,$$

substracting the average nucleon mass M nuc = ( M ν + M / π ) 2. A second derivative,

$$E _ { \text{sym} } ( n _ { b } ) = \frac { 1 } { 2 } \, \frac { d ^ { 2 } } { d \alpha ^ { 2 } } \frac { E } { A } \Big | _ { \alpha = 0 } \,, \\ \text{gy. It characterizes how the isospin dependence of the energy per nucleon}$$

gives the symmetry energy. It characterizes how the isospin dependence of the energy per nucleon changes with the baryon density.

The energy per nucleon of symmetric matter near saturation can be expanded as

$$E _ { 0 } ( n _ { b } ) = \frac { E } { A } ( n _ { b }, 0 ) = - B _ { \text{sat} } + \frac { 1 } { 2 } K x ^ { 2 } + \frac { 1 } { 6 } Q x ^ { 3 } + \dots \,,$$

with x = ( n b -n sat ) / (3 n sat ) without a term linear in x because the pressure at saturation is P n ( sat ) = 0. The incompressibility coefficient and skewness parameter

$$K = 9 n _ { \text{sat} } ^ { 2 } \, \frac { d ^ { 2 } } { d ( n _ { b } ) ^ { 2 } } E _ { 0 } \Big | _ { n _ { \text{sat} } } = 9 \, \frac { d P } { d n _ { b } } \Big | _ { n _ { \text{sat} } } \, \\ \text{ameter}$$

$$Q = 2 7 n _ { \text{sat} } ^ { 3 } \frac { d ^ { 3 } } { d ( n _ { b } ) ^ { 3 } } E _ { 0 } \Big | _ { n _ { \text{sat} } }, & & ( 6 2 ) \\ \text{third derivatives of } E _ { 0 }. \text{ Similarly, the symmetry energy can be written as}$$

∣ are found as second and third derivatives of E 0 . Similarly, the symmetry energy can be written as

$$E _ { \text{sym} } ( n _ { b } ) = J + L x + \frac { 1 } { 2 } K _ { \text{sym} } x ^ { 2 } + \dots \,,$$

with the symmetry energy at saturation J = E sym ( n sat ), the slope parameter and the symmetry incompressibility

$$L = 3 n _ { \text{sat} } \, \frac { d E _ { \text{sym} } } { d n _ { b } } \Big | _ { n _ { \text{sat} } } \,,$$

$$K _ { \text{sym} } = 9 n _ { \text{sat} } ^ { 2 } \, \frac { d ^ { 2 } E _ { \text{sym} } } { d ( n _ { b } ) ^ { 2 } } \Big | _ { n _ { \text{sat} } } \,.$$

Table 1 Masses of nucleons and Lorentz-vector mesons used in the model parameterisations of this work.

| particle    | p              | n             |   ω |   ρ |
|-------------|----------------|---------------|-----|-----|
| mass in MeV | 938 . 27208816 | 939 . 5654205 | 783 | 763 |

It is also useful to define the socalled volume part of the isospin compressibility

$$K _ { \tau, v } = K _ { \text{sym} } - 6 L - \frac { Q L } { K } \,,$$

that describes the change of the incompressibility with the variation in α [15, 30]. These quantities, together with the effective Dirac mass at saturation M ∗ sat = M ∗ nuc ( n sat ) = M nuc -C n σ ( s ) sat , are convenient to characterize the equation of state of nuclear matter around saturation.

## 3 Determination of model parameters

The relativistic EDF depends on several parameters: the masses of the particles and the coefficients that describe the strength of the couplings and their density dependence. The masses of the nucleons and the mesons, except that of the σ meson are usually chosen to be identical or close to the experimental values. In the present application of the relativistic EDF the δ meson is not included in the model. The remaining fixed mass parameters of the nucleons and Lorentz-vector mesons are given in Table 1 for the models in this work. The mass of the σ meson is used as an independent parameter because it is not well defined by experimental data and affects the surface properties of nuclei.

The actual values of the model parameters are determined by fitting selected observables of finite nuclei that are sensitive to different contributions in the EDF and are obtained more-or-less directly from experiment. More indirectly obtained quantities, e.g., characteristic nuclear matter parameters, that can depend strongly on the extraction method from actual experimental data, or other constraints, e.g., from ab-initio calculations of nuclear matter or from astrophysical observations, are not used in this procedure. A comparison to these data will be performed after the model parameters have been determined.

## 3.1 Functional form of density dependent couplings

The meson-nucleon couplings Γ j are assumed to be functions of the vector density /rho1 v = n ( v ) . In this work, the functional form of the dependence is chosen as introduced originally in [11] with

$$\Gamma _ { j } ( \varrho ) = \Gamma _ { j } ^ { ( 0 ) } f _ { j } ( x ) \,,$$

where Γ (0) j = Γ ( j /rho1 ref ) is the coupling at a reference density /rho1 ref and the functions f j depend on the parameter x = /rho1 v //rho1 ref ≥ 0. For the isoscalar mesons j = σ, ω a rational form

$$f _ { j } ( x ) = a _ { j } \frac { 1 + b _ { j } ( x + d _ { j } ) ^ { 2 } } { 1 + c _ { j } ( x + d _ { j } ) ^ { 2 } } \,,$$

with four coefficients a j , b j , c j , and d j is used. To avoid a pole in f j , the coefficients c j , and d j have to be non-negative. The constraints f j (1) = 1, and f ′′ j (0) = 0, reduce the number of independent coefficients to two with a j = [1 + c j (1 + d j ) 2 ] / [1 + b j (1 + d j ) 2 ], and c j = 1 (3 / d 2 j ). The coupling function for the isovector mesons j = ρ is chosen as an exponential,

$$f _ { \rho } ( x ) = \exp [ - a _ { \rho } ( x - 1 ) ] \,,$$

with only one coefficient a ρ . Obviouly, other functions can be and have been used to describe the density dependence of the couplings, see, e.g., Ref. [15]. Also a dependence on the scalar density n ( s ) or other densities can be explored, see, e.g., Ref. [25].

In the model parameterisations of this work, there are three independent parameters for the isoscalar mesons σ and for ω and two for the isovector meson ρ , thus eight parameters in total that need to be determined for the nuclear-matter calculations besides the mass of the σ meson. In the actual fitting procedure, the saturation density of symmetric nuclear matter n sat is used as reference density /rho1 ref , but the couplings Γ (0) j at /rho1 ref and the coefficients a j , b j , c j , and d j are not used directly as independent quantities that are varied. Instead, the following characteristic nuclear matter parameters are selected as variable quantities:

- 1. the effective Dirac mass at saturation M ∗ sat ,
- 2. the binding energy per nucleon at saturation B sat ,
- 3. the incompressibility coefficient at saturation K ,
- 4. the symmetry energy at saturation J ,
- 5. the slope parameter of the symmetry energy at saturation L .

Together with the condition of vanishing pressure at n sat , this gives already six parameters. The remaining two quantities that are used in the actual fit are the ratios

$$r _ { 1 0 } = \frac { f _ { \omega } ^ { \prime } ( 1 ) } { f _ { \omega } ( 1 ) } \,, \quad \text{and} \quad r _ { 2 1 } = \frac { f _ { \omega } ^ { \prime \prime } ( 1 ) } { f _ { \omega } ^ { \prime } ( 1 ) } \,,$$

with the density dependent function f ω ( x ) of the ω meson. It is assumed that the coupling functions are decreasing with the density, hence f ′ ω (1) has to be smaller than zero. Then the functional form of f ω leads to the constraint -3 &lt; r 21 &lt; 0 for the second ratio in Eq. (70). The actual procedure of the parameter conversion is given in detail in Appendix A.

In the description of finite nuclei, two additional tensor couplings Γ Tω and Γ Tρ appear if the tensor terms in the Lagrangian density (1) are considered. These constants are used directly as independent parameters in the fit. Hence there are in total 11 (9) independent parameters in models with (without) tensor couplings taking into account the mass of the σ meson M σ as another undetermined quantity.

## 3.2 Observables and nuclei in the fit

The choice of the observables that are included in fitting the model parameters will affect the final results. Only quantities should be selected that are well determined by experiments. Different strategies are followed by various authors. In the present work, the following nuclear observables of a nucleus are used:

- 1. the binding energy B ,
- 2. the charge radius r c ,
- 3. the diffraction radius r d ,
- 4. the surface thickness σ ,
- 5. the root-mean-square radius r n of a single valence neutron level above a shell closure,
- 6. the spin-orbit splitting ∆ E so for two levels of a nucleon η in states with the same quantum numbers n η and l η but different j η = l η ± 1 / 2,
- 7. the constraint isoscalar monopole giant resonance energy E mono .

They are calculated after the coupled equations for a nucleus ( Z, N ) are solved and can be compared with experimental data.

## 3.2.1 Values of observables from theory

The binding energy is derived as

$$B ( Z, N ) = M ( Z, N ) - E _ { \text{nucleus} } ( Z, N ) \,,$$

from the energy (41). The charge radius, the diffraction radius, and the surface thickness are obtained from the charge form factor of a spherical nucleus. It depends on a momentum /vector q and is

calculated as,

$$F _ { Z, N } ( q ) = \left [ F _ { Z } ( q ) F _ { \pi } ( q ) + F _ { N } ( q ) F _ { \nu } ( q ) \right ] F _ { \text{corr} } ( q ) \,,$$

with form factors of the proton and neutron distributions in the nucleus, single-nucleon form factors and a center-of-mass correction. The density distributions (36) of neutron and protons are used to calculate the form factors

$$F _ { Z } ( q ) = \int d ^ { 3 } r \sum _ { k _ { \pi }, m _ { \pi } } w _ { \pi k _ { \pi } m _ { \pi } } n _ { \pi k _ { \pi } m _ { \pi } } ( \vec { r } ) \, \exp \left ( i \vec { q } \cdot \vec { r } \right ) \,,$$

and,

$$F _ { N } ( q ) = \int d ^ { 3 } r \sum _ { k _ { \nu }, m _ { \nu } } w _ { \nu k _ { \nu } m _ { \nu } } n _ { \nu k _ { \nu } m _ { \nu } } ( \vec { r } ) \exp \left ( i \vec { q } \cdot \vec { r } \right ) \,, \\ \text{nherical} v \, \text{summetric in momentum since due to the spherical symmetric$v$ of the density$v}$$

that are spherically symmetric in momentum space due to the spherical symmetry of the density distributions. The single-nucleon form factors are parameterized [5, 31] as

$$F _ { \eta } ( q ) = \sum _ { k = 1 } ^ { 4 } \frac { A _ { \eta k } } { 1 + B _ { \eta k } q ^ { 2 } } \,,$$

with coefficients A ηk , and B ηk , that are given in Table 2. The center-of-mass correction

$$F _ { \text{corr} } ( q ) = \exp \left [ \frac { 3 q ^ { 2 } } { 1 6 E _ { \text{corr} } M ( Z, N ) } \right ] \,,$$

depends on the center-of-mass correction (44) of the energy. The charge density of the nucleus is finally obtained from

$$n _ { \text{charge} } ( r ) = \int \frac { d ^ { 3 } q } { ( 2 \pi ) ^ { 3 } } \, F _ { Z, N } ( q ) \, \exp \left ( - i \vec { q } \cdot \vec { r } \right ) \,. \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \end{cases}$$

It is used to calculate the root-mean-square charge radius

$$r _ { c } = \left [ \frac { 1 } { Z } \int d ^ { 3 } r \, r ^ { 2 } \, n _ { \text{charge} } ( r ) \right ] ^ { 1 / 2 } \,,$$

of the nucleus. The form factor (72) at zero momentum has the values F Z,N (0) = Z and decreases with q until it reaches a first zero at q 0 . This value is used to obtain the diffraction radius

$$r _ { d } = \frac { 4. 4 9 3 4 0 9 4 5 8 } { q _ { 0 } } \,,$$

Table 2 Coefficients [5, 31] in the parameterisation (75) of the nucleon form factors F η ( q ) ( η = π, ν ).

| k A πk      | B πk [fm 2 ]   | A νk   | B νk [fm 2 ]   |
|-------------|----------------|--------|----------------|
| 1 0 . 312   | 0 . 16667      | 1      | 0 . 04833      |
| 2 1 . 312   | 0 . 06658      | - 1    | 0 . 05833      |
| 3 - 0 . 709 | 0 . 02269      | 0      | 0              |
| 4 0 . 085   | 0 . 006485     | 0      | 0              |

where the numerical factor is the first zero of the Bessel function j 1 ( x ) = -j ′ 0 ( x ). Further increasing q leads to the first minimum of F Z,N ( q ) at q m that is used to calculate the surface thickness

$$\sigma ( Z, N ) = \frac { 1 } { q _ { m } } \left \{ 2 \ln \left [ \frac { Z F _ { \text{box} } ( q _ { m } ) } { F _ { Z, N } ( q _ { m } ) } \right ] \right \} ^ { 1 / 2 } \,,$$

with the form factor F box ( q ) = 3 j 1 ( qr d ) / qr ( d ) of a box of radius r d , c.f. Ref. [5], and Eq. (18a) of Ref. [6] with corrections. The root-mean-square radius of a single-neutron state is defined as

$$r _ { \nu } = \left [ \int d ^ { 3 } r \, r ^ { 2 } \, n _ { \nu k _ { \nu } m _ { \nu } } ( r ) \right ] ^ { 1 / 2 } \,,$$

similar as the charge radius (78). The spin-orbit splitting is given by the difference

$$\Delta E _ { s o } ( \eta n _ { \eta } l _ { \eta } ) = E _ { \eta k _ { \eta } ^ { - } } - E _ { \eta k _ { \eta } ^ { + } }$$

of the single-nucleon energies of the levels with identical principal quantum number n η , and k ± η = ( n , κ η ± η ), where κ ± η = ( j ± η -l η )(2 j ± η + 1) for j ± η = l η ± 1 / 2. The constraint isoscalar monopole giant resonance energy E mono is obtained from a constraint calculation of the nucleus as described in Ref. [1]. A contribution

$$V _ { \lambda } = - \lambda r ^ { 2 } \,,$$

with parameter λ , is added to the vector potential (12) of the nucleons and the coupled equations for the fields are solved self-consistently for given λ . Then the root-mean-square radius

$$\langle r ^ { 2 } \rangle _ { \lambda } = \int d ^ { 3 } r \, r ^ { 2 } \, n _ { \lambda } ^ { ( v ) } ( \vec { r } ) \,,$$

of the total vector density n ( v ) λ of the obtained configuration is calculated. The constraint isoscalar monopole giant resonance energy

$$E _ { \text{mono} } = \sqrt { \frac { m _ { 1 } } { m _ { - 1 } } } \,,$$

can be found from the ratio of moments

$$m _ { 1 } = 2 \frac { \hslash ^ { 2 } } { M _ { \text{nuc} } } \langle r ^ { 2 } \rangle _ { \lambda = 0 } \,,$$

and,

$$m _ { - 1 } = \frac { 1 } { 2 } \, \frac { d } { d \lambda } \langle r ^ { 2 } \rangle _ { \lambda } \Big | _ { \lambda = 0 } \. & & ( 8 7 ) \\ \text{five-point formula is used with } \lambda _ { n } = n \times 0. 1 2 5 \, \text{keV and } n = - 2, - 1, 0, 1, 2$$

∣ In the actual calculation, a five-point formula is used with λ n = n × 0 125 keV and . n = -2 , -1 0 1 2 , , , to calculate the derivative.

## 3.2.2 Selection of nuclei and experimental data

The description of nuclei with the relativistic EDF presented here can only be applied to nuclei with spherical symmetry. Also pairing effects are not taken into account. Thus a comparison of theory and experiment can be reasonably performed only for magic or semi-magic nuclei. This limits the amount of data that can be included in the fit of the model parameters. However, this restriction is not too severe since the essential features of strongly interacting systems are sufficiently well constrained by the selection of the fitted data in this approach. In order to fix the isospin dependence of the effective interaction well, nuclei with large | N -Z | should be included, preferably on both sides of the valley of stability. The nuclei actually considered can be read off from Tables 6, 7, and 9. There are 18 (semi)magic nuclei and 2 with one valence neutron above a doubly magic core. Not all data from the list in Subsection 3.2 are always available for all nuclei.

Binding energies B Z,N ( ) of nuclei are obtained from the Atomic Mass Evaluation 2020 [32] by multiplying the binding energy per nucleon given there with the mass number A = Z + N . These data are tabulated in Table 6. Experimental data for the charge radius, diffraction radius, surface thicknesses and neutron radii can be found in Table 7. The charge radii and quantities related to the charge form factor are extracted from [33-36] using the Fourier-Bessel parameterisation if available. The experimental valence neutron radii are taken from [37, 38]. Experimental spin-orbit splittings are given in Table 8. They are derived from the nuclear level schemes collected at the National Nucleat Data Center [39, 40]. Experimental constraint isoscalar monopole giant resonance energies [41] are listed in Table 9 for selected nuclei. The choice of the constraint isoscalar monopole giant resonance energy as an observable in the fit greatly helps to obtain a reasonable incompressibility coefficient K of nuclear matter, a quantity that is often set to a predefined value since it is not well constrained by other quantities of nuclei.

## 3.3 Objective function and uncertainties

In order to determine optimal parameter values of the relativistic EDF, an appropriate fitting procedure has to be applied. This is usually done by minimizing an objective function that quantifies the deviation of the theoretical predictions for certain observables from experimental data. A crucial question is the size of the uncertainties that are used for the different observables, e.g., in a χ 2 minimization. In the parameter fit of an EDF it is generally not reasonable to use the experimental uncertainties because they are for most quantities much smaller than the precision that can be achieved by an EDF. Instead, often some 'reasonable' values are chosen that may, however, not reflect the 'true' uncertainties that are expected for such an EDF. Thus, in the present work a fitting scheme is employed that also attempts to determine the uncertainties, c.f. Ref. [25].

In the present approach to fit the model parameters, N obs = 7 types of observables are considered as specified in Section 3.2 with N (obs) i data per observable ( i = 1 , . . . , N obs ). Thus an objective function,

$$\chi ^ { 2 } ( \{ p _ { k } \} ) & = \sum _ { i = 1 } ^ { N _ { o b s } } \chi _ { i } ^ { 2 } ( \{ p _ { k } \} ) \,, \\ \text{on } N _ { n a r } \text{ parameters } p _ { k } \, ( k = 1 \dots N _ { n a r } ) \text{ with } N _ { n a r } = 1 1 \, ( 9 ) \text{ for models}$$

is defined that depends on N par parameters p k ( k = 1 , . . . , N par ) with N par = 11 (9) for models with (without) tensor couplings. The contributions for each observable O i to (88) are given by the standard expression

$$\chi _ { i } ^ { 2 } ( \{ p _ { k } \} ) = \sum _ { n = 1 } ^ { N _ { i } ^ { ( \text{obs} ) } } \left [ \frac { \mathcal { O } _ { i } ^ { ( \text{model} ) } ( n, \{ p _ { k } \} ) - \mathcal { O } _ { i } ^ { ( \text{exp} ) } ( n ) } { \Delta \mathcal { O } _ { i } } \right ] ^ { 2 } \,,$$

where the n = 1 , . . . , N obs i data of an observable O i depend on the parameters { p k } . The total number of data is given by

$$N _ { \text{data} } = \sum _ { i = 1 } ^ { N _ { \text{obs} } } N _ { i } ^ { ( \text{obs} ) } \,, & & ( 9 0 ) \\ \colon & \text{ees of freedom is } N _ { \text{dof} } = N _ { \text{data}$$

and thus the number of degrees of freedom is N dof = N data -N par .

The uncertainties for each observable O i are denoted by ∆ O i . Since binding energies of stable nuclei are much more precisely known than those of exotic nuclei, their uncertainties vary over a large range. Thus the uncertainty of the observable O i = B Z,N ( ) in the function (89) is assumed to compose of two contributions as,

$$\Delta \mathcal { O } _ { i } = \left [ \left ( \Delta \mathcal { O } _ { i } ^ { ( \exp ) } \right ) ^ { 2 } + \left ( \Delta \mathcal { O } _ { i } ^ { ( \text{fit} ) } \right ) ^ { 2 } \right ] ^ { 1 / 2$$

with an experimental uncertainty ∆ O (exp) i and a fit uncertainty ∆ O (fit) i . For all other observables, the experimental uncertainty is neglected in the fits of this work. A reasonable estimate of the fit

uncertainties will be obtained if, and,

$$\frac { \chi _ { i } ^ { 2 } ( \{ p _ { k } \} ) } { N _ { i } ^ { ( \text{obs} ) } } = \frac { \chi ^ { 2 } ( \{ p _ { k } \} ) } { N _ { \text{data} } } \,, \quad \quad \quad$$

i.e., each observable contributes on the average equally to the total χ 2 . Theses conditions allow to determine the fit errors ∆ O (fit) i during the fit by a rescaling. Thus the general fit procedure proceeds in two steps: variation of the model parameter to minimize (88) until a certain convergence is reached, and a rescaling of the fit uncertainties. These two steps are executed repeatedly until the fit uncertainties and model parameters have converged to the desired precision. To find the optimal values of the model parameters in their space of dimension N par , a simplex method [42] was employed as in Ref. [25]. It only relies on values of the χ 2 function and does not need derivatives. Once the minimum of the total χ 2 function is found, the solution can be represented as a vector /vector p min = ( p min 1 , . . . , p min N par ) in the space of parameters. A true (local) minimum of χ 2 is reached when the matrix of second derivatives,

$$\text{aives,} \\ \mathcal { M } _ { i j } & = \frac { 1 } { 2 } \, \frac { \partial ^ { 2 } \chi ^ { 2 } } { \partial p _ { i } \partial p _ { j } } \Big | _ { \vec { p } _ { \min } } \,, \\ \text{$$

∣ has only positive eigenvalues. If this is not the case, the eigenvectors with negative eigenvalues indicate directions for a further reduction of χ 2 . This information is used every now and then to accelerate the minimization. However, the calculation of M ij is expensive as it needs 2 N 2 par + 1 function evaluations. Finally, it is possible to determine the uncertainty

$$\Delta \mathcal { O } = \sqrt { \overline { \left ( \Delta \mathcal { O } \right ) ^ { 2 } } } \,,$$

of any observable O and the correlation coefficient

$$c _ { \mathcal { O } _ { 1 } \mathcal { O } _ { 2 } } = \frac { \overline { \Delta \mathcal { O } _ { 1 } \Delta \mathcal { O } _ { 2 } } } { \sqrt { \overline { ( \Delta \mathcal { O } _ { 1 } ) ^ { 2 } } } \, \overline { ( \Delta \mathcal { O } _ { 2 } ) ^ { 2 } } } \,,$$

between two observables O 1 and O 2 from

$$\overline { \Delta \mathcal { O } _ { 1 } \Delta \mathcal { O } _ { 2 } } = \sum _ { i j } \frac { \partial \mathcal { O } _ { 1 } } { \partial p _ { i } } \left ( \mathcal { M } ^ { - 1 } \right ) _ { i j } \frac { \partial \mathcal { O } _ { 2 } } { \partial p _ { j } } \,,$$

that contains the inverse of the matrix M and first derivatives of the observables with respect to the parameters.

$$\chi ^ { 2 } ( \{ p _ { k } \} ) / N _ { \text{dof} } = 1 \,,$$

## 3.4 Types of energy density functionals

The main aim of the present work is the development of a new parameterisation of a relativistic EDF with tensor couplings, denoted by 'DDT' in the following, in addition to the usual minimal couplings of the mesons to the nucleons. In order to compare the main differences of this EDF to a model without tensor interactions, the results of the DD2 model [43] will be shown, that uses the same form of the relativistic EDF. Although a different fitting strategy with smaller numbers of nuclei and observables was used to determine the parameters of the DD2 model, a comparison of the results with those of models with tensor couplings will be instructive even though the differences will originate not completely from these additional contributions in the EDF.

A particular issue is the effect of exchange contributions that are neglected in the relativistic EDF, which is constructed in Hartree approximation from the Lagrangian density. This concerns particularly the Coulomb interaction. The source density of the electromagnetic field, that is found by solving Eq. (26), includes the charge distribution of all protons so that the Coulomb field in the vector potential (12), which a proton feels in the mean-field created by all particles, also contains a self-interaction contribution that would not arise with the proper use of an antisymmetrized manybody wave function. The exchange term is sometimes added in form of the Slater approximation [44]. However, this is not done in this work. As an alternative, a simple correction to reduce the self-interaction term is to multiply the field A 0 in Eq. (12) by the factor ( Z -1) /Z in the calculation of nuclei, as used in the DD2 parameterisation. Then the vector potential V π ( /vector r ) has the correct asymptotic form for large radii when a proton is separated from a nucleus. To study the effect of this modification, a variation with this reduction factor, DDTC, of the original DDT parameterisation is also developed.

A more detailed consideration of the contributions to the energy of a nucleus reveals, however, that it is admissible to neglect the exchange part in the Coulomb energy. In fact, there are two contributions that in total compensate for the exchange term. On the one hand, there is the Coulomb displacement energy of analog states in mirror nuclei that, in the mean-field approximation, e.g., for the 41 Sc 41 Ca pair is smaller by about 7%. This is known as the Nolen-Schiffer anomaly [45] which is solved by the taking into account the charge symmetry breaking (CSB) contribution in the nucleon-nucleon interaction, and long-range correlation (LRC) effects [46]. Ignoring the factor ( Z -1) /Z in the Coulomb direct term, which leads to including the Coulomb self-interaction terms, and ignoring the Coulomb exchange term corresponds to including the contributions of the CSB and LRC effects. Thus only the direct part of the Coulomb contribution to the energy should be included in the calculation, see also subparagraph IIA in Ref. [1].

Table 3 Parameters of the coupling functions, reference densities, and σ -meson masses for the different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| parameter                     | DDT                                                    | DDT                                                                                       |                                                                     | DDTC                                                             | DD2                                                     |
|-------------------------------|--------------------------------------------------------|-------------------------------------------------------------------------------------------|---------------------------------------------------------------------|------------------------------------------------------------------|---------------------------------------------------------|
| Γ (0) σ a σ b σ c σ d σ Γ (0) | 8 . 630588 1 . 142684 1 . 854752 2 . 193260 0 . 389847 | (0 . 289206) (0 . 009581) (0 . 627900) (0 . 741533) (0 . 065899) (0 . 29709) (1 . 032068) | 8 . 580686 1 . 311104 2 . 315517 3 . 213878 0 . 322051 11 . 09567 . | (0 . 846337) (0 . 048916) (2 . 062872) (2 . 826693) (0 . 146240) | 10 . 686681 1 . 357630 0 . 634442 1 . 005358 0 . 575810 |
| ω a ω b ω c ω d ω (0)         | 3 . 747188 0 . 684812                                  | (0 . (0 . 050891)                                                                         | 1 434242 9 . 292208 13 . 65208 0 . 156257                           | (0 . 086777) (0 . 054904)                                        | 3 . 626940 0 . 518903                                   |
| a ρ Γ                         | 3 . 200838                                             | (0 . 621091)                                                                              | 0 . 675146 3 . 489206                                               | (1 . 881621) (3 . 39598)                                         | 0 . 000000                                              |
| Tω Γ Tρ                       |                                                        |                                                                                           |                                                                     |                                                                  |                                                         |
|                               | 11 . 14407 2 . 891633                                  |                                                                                           |                                                                     | (0 . 86419) (0 . 204994)                                         | 13 . 342362 1 . 369718                                  |
| M σ                           | 511 . 9820                                             | (6 . 5387)                                                                                | 510 . 5184                                                          | (19 . 2068)                                                      | 546 . 2125                                              |
|                               | 1439 . 522                                             | (1290 . 251)                                                                              |                                                                     | (10 . 52430)                                                     | 0 . 496475                                              |
|                               | 4164 . 427                                             | (5086 . 006)                                                                              |                                                                     | (17 . 02927)                                                     | 0 . 817753                                              |
|                               | 0 . 008947                                             | (0 . 005460)                                                                              |                                                                     | (0 . 097456)                                                     | 0 . 638452                                              |
| Γ ρ                           |                                                        | 070045)                                                                                   | 3 . 693565                                                          |                                                                  |                                                         |
| -                             | .                                                      |                                                                                           | 0 .                                                                 |                                                                  |                                                         |
|                               | 18 . 67809                                             | (3 . 74586)                                                                               | 20 . 50426                                                          |                                                                  | 0 . 000000                                              |
| /rho1 ref [fm 3 ]             | 0 15493496                                             | (0 . 00075982)                                                                            | 15304168                                                            | (0 . 00094838)                                                   | 0 . 149065                                              |
| [MeV]                         |                                                        |                                                                                           |                                                                     |                                                                  |                                                         |

## 4 Results and discussion

After performing the fit of the parameters for the two models with tensor couplings, various quantities for nuclei, nuclear matter and neutron stars can be calculated and compared to experimental data. Before they are presented, in the following subsections, a discussion of the model parameters and of the obtained density dependence of the coupling functions is of interest since they are essential to obtaining a good quality fit.

The parameters of the three relativistic EDF models of the present work are specified in Table 3. For the DDT and DDTC parameterisations, uncertainties are given as determined with the method described in Subsection 3.3. For the DD2 model, these are not available because the parameter set does not correspond to a minimum of χ 2 with the present set of observables and model uncertainties. The most obvious differences are the large values of the tensor couplings Γ Tω , and Γ Tρ , and the small σ meson masses M σ as compared to the value of the DD2 model without tensor couplings. This is not surprising since Γ Tω and M σ directly affect the description of the surface properties of nuclei and thus there is a strong correlation between them. The DD2 model is also characterized by a rather small saturation density. The isoscalar couplings Γ (0) σ and Γ (0) ω at the reference density of the DDT and DDTC models are considerably smaller than those of the DD2 model. On the other hand, the couplings Γ (0) ρ and the coefficient a ρ are larger. Particularly striking are the large coefficients b ω and c ω with a small d ω for the DDT model. This will lead to an unusual density dependence of the coupling Γ ω , see below.

The correlation coefficients of the fit parameters in the DDT model are given in Figure 1 where boxes in red color indicate the largest absolute values whereas grey boxes correspond to the

smallest moduli of the coefficients. As expected, a first group with M σ , M ∗ sat , and Γ Tω shows strong correlations since they all affect the description of the nuclear surface. In contrast, Γ Tρ is hardly correlated to any other parameter of the fit. The other two parameters related to the isovector part of the interaction, J , and L , however, show a strong correlation. Apart from this dependence, J seems to be practically uncorrelated to all other fit parametes. Another pair of parameters, r 10 , and r 21 , determines the density dependence of the isoscalar couplings. They are also strongly correlated. A somewhat lesser strength of correlation is found for the incompressibility coefficient K , and r 10 with the first group of parameters.

The coupling functions Γ ( i n b ) generally decrease with increasing baryon density n b , see Figure 2. They approach a finite value in the high-density limit by construction, and go to zero for Γ ρ , in particular. The isoscalar couplings for the models with tensor couplings (DDT, and DDTC) are significantly smaller than those of the model without tensor couplings (DD2). This can also be observed for the absolute values of the couplings at the reference density in Table 3. The strong rise of the ω meson coupling in the DDT model towards very small densities is particularly noteworthy. This behavior is a consequence of the large values of the parameters b ω and c ω in combination with a small d ω . The isoscalar couplings, Γ ω and Γ σ , of the DDTC parameterisation have a much larger uncertainty as compared to those of the DDT model, but their shape is very similar to the

Fig. 1 Correlation matrix of the fit parameters in the DDT model.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-8ddcbdcd-67d9-4456-a89a-24148927f40a.jpg)

Fig. 2 Coupling functions Γ ( i n b ) of the (a) σ meson, (b) ω meson, and (c) ρ meson for the different parameterisations of the relativistic EDF with uncertainty bands for DDTC and DDT.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-6dcc4ef8-4409-495e-9927-8e0da49966f2.jpg)

Table 4 Fit uncertainties ∆ O (fit) i of the observables O i for DDT.

| observable O i              |   B [MeV] |   r c [fm] |   r d [fm] |   σ [fm] |   r n [fm] |   ∆ E so [MeV] |   E mono [MeV] |
|-----------------------------|-----------|------------|------------|----------|------------|----------------|----------------|
| fit uncertainty ∆ O (fit) i |  0.619311 |   0.013364 |   0.017155 | 0.026851 |   0.008249 |       0.240832 |       0.430714 |

corresponding couplings of the DD2 model. The density dependence of the ρ meson couplings is very similar for all models with slightly higher values for the models with tensor couplings.

## 4.1 Nuclei

Before analyzing the predictions of individual properties of finite nuclei, it is worthwhile to look at the fit uncertainties of the observables, see Table 4, and their contributions to the χ 2 function that are collected in Table 5 for the three parameterisations. The uncertainties of the energy quantities, i.e., the binding energy, the spin-orbit splitting, and the constraint isocalar monopole giant resonance energy, are all on a sub-MeV level, whereas the uncertainties of the size quantities, i.e., radii and skin thickness, are typically in the order of a few hundredths of a fm with the error in the root-mean-square (rms) valence neutron radius beeing the smallest. These results give an impression on the actual quality a description of nuclei with a relativistic EDF can achieve and how reasonable values have to be chosen in other attempts to find parameters of such models from fits to finite nuclei.

A comparison of the total χ 2 , see Table 5, and the contributions for the three models reveals interesting differences. First of all, it is evident that the DDT model gives the best description of the data with the lowest total χ 2 . By construction, its value is given by the number of degrees of freedom N dof = 39. The other parameterisations, in particular the one without tensor couplings,

are considerably worse. In case of the DDTC parameterisation this is mostly caused by the larger errors in the binding energies and the spin-orbit splittings. The values of χ 2 of the DD2 model is completely dominated by the difference of the valence neutron radii from the experimental values. Even neglecting their contribution to χ 2 , the descripton of experimental data with the DD2 model is almost always much worse than that of the DDT and DDTC models. All observables that are determined strongly by the surface description of a nucleus, viz., the radii, the skin thickness, and the spin-orbit splitting, are on the average better described in the EDFs with tensor couplings. Only the surface thickness σ is slightly better described by the DD2 model.

The first quantity of major interest is the binding energy B Z,N ( ) of a nucleus that is predicted by a given parameterisation of a relativistic EDF. The results of the binding energy calculations of the three models of the present work are listed in Table 6, for nuclei considered in the fit, in

Table 5 Contributions χ 2 i ( { p k } ) of the different observables and total χ 2 ( { p k } ) for the different parameterisations of the relativistic EDF. The bottom line gives the sum of the quantities in the same column above.

| observable   |   N (obs) i |   DDT |   DDTC |    DD2 |
|--------------|-------------|-------|--------|--------|
| B ( Z, N )   |          16 | 12.48 |  31.79 | 100.02 |
| r c          |           8 |  6.24 |  12.61 |  16.06 |
| r d          |           5 |  3.9  |   6.3  |  26.27 |
| σ            |           5 |  3.9  |   5.08 |   2.24 |
| r n          |           2 |  1.56 |   4.01 | 578.17 |
| ∆ E so       |          10 |  7.8  |   7.47 |  40.66 |
| E mono       |           4 |  3.12 |   2.86 |   3.61 |
|              |          50 | 39    |  70.11 | 767.01 |

Fig. 3 Binding energy differences B theory ( Z, N ) -B experiment ( Z, N ) of nuclei used in the parameter fits with theoretical uncertainties indicated by the error bars.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-9f6c7066-0b09-4c84-be58-2365e49981bf.jpg)

Table 6 Binding energies B Z,N ( ) of nuclei in MeV from experiment [32] in comparison to theoretical values from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets. The bottom line gives the root-mean-square deviation of the theoretical predictions from the experiment in MeV.

| nucleus   | experiment   |      DDT | DDT     |     DDTC | DDTC    |      DD2 |
|-----------|--------------|----------|---------|----------|---------|----------|
| 16 O      | 127.619      |  128.193 | (0.262) |  128.469 | (0.285) |  127.947 |
| 24 O      | 168.952      |  168.368 | (0.327) |  167.89  | (0.440) |  169.13  |
| 28 O      | 167.664      |  167.888 | (0.701) |  167.672 | (0.839) |  173.005 |
| 34 Si     | 283.464      |  284.471 | (0.300) |  283.732 | (0.300) |  284.044 |
| 34 Ca     | 243.882      |  244.06  | (0.234) |  245.259 | (0.259) |  246.245 |
| 40 Ca     | 342.052      |  342.161 | (0.362) |  342.366 | (0.444) |  342.203 |
| 48 Ca     | 416.001      |  416.545 | (0.224) |  416.225 | (0.265) |  414.929 |
| 48 Ni     | 347.328      |  346.74  | (0.290) |  348.87  | (0.347) |  350.092 |
| 56 Ni     | 483.996      |  483.435 | (0.380) |  483.523 | (0.412) |  482.346 |
| 68 Ni     | 590.408      |  589.675 | (0.374) |  588.546 | (0.367) |  590.595 |
| 78 Ni     | 642.564      |  642.398 | (0.470) |  642.46  | (0.464) |  642.048 |
| 90 Zr     | 783.897      |  782.939 | (0.321) |  782.332 | (0.376) |  780.727 |
| 100 Sn    | 825.163      |  825.984 | (0.490) |  825.978 | (0.551) |  826.231 |
| 132 Sn    | 1102.843     | 1103.17  | (0.394) | 1103.26  | (0.456) | 1103.39  |
| 140 Ce    | 1172.683     | 1172.89  | (0.423) | 1173.1   | (0.625) | 1172.81  |
| 208 Pb    | 1636.430     | 1636.35  | (0.531) | 1636.48  | (0.530) | 1635.92  |
|           |              |    0.562 |         |    0.919 |         |    1.907 |

comparison to experimental values from the Atomic Mass Evaluation 2020 [32]. In order to have a better impression of the quality of the models, the differences between the experimental and theoretical data are depicted in Figure 3. It is seen that the models with tensor coupling show less fluctuations and deviate on the average less strongly from the experimental data than the models without tensor couplings. This is also reflected in the root-mean-square deviation of the theoretical values from the experimental data in the bottom line of Table 6. The trend of these numbers agrees with the variation of the binding energy contribution to the total χ 2 shown in Table 5.

Also interesting observations for individual sets of nuclei can be made. In constrast to the DD2 parameterisation, the binding energy of 28 O is smaller than that of 24 O in the DDT and DDTC parameterisations, a result that is difficult to achieve in most other relativistic EDFs. Similar patterns can be seen in Figure 3 when comparing different parameterisations, e.g., for the heaviest nuclei. There are also striking effects for the set of Ni nuclei with different trends when comparing the models. It is seen in Table 6 that the calculated binding energies of 48 Ni obtained for the DDT model without the ( Z -1) /Z factor is significantly closer to the experimental value than that of the DDTC and DD2 models with the ( Z -1) /Z factor.

Theoretical values of the quantities related to the charge form factor are given in Table 7 in comparison to experimental data for the nuclei in the fit. The theoretical uncertainties are generally in the order of 0.005 fm. Differences between theory and experiment are depicted in Figure 4 with three panels for the charge radius, the diffraction radius, and the surface thickness. There is a good agreement with typical differences between theory and experiment of less than 0.02 fm, and a clear trend in the deviations is not seen. Only the surface thickness of 90 Zr is overestimated by

Table 7 Experimental data of charge radii, diffraction radii, surface thicknesses [33-36] and valence neutron radii [37, 38] in fm of nuclei used in the fit, in comparison to theoretical values from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

|                            | nucleus   |   experiment |   DDT | DDT     |   DDTC | DDTC    |   DD2 |
|----------------------------|-----------|--------------|-------|---------|--------|---------|-------|
| charge radius r c          | 16 O      |        2.699 | 2.705 | (0.005) |  2.689 | (0.011) | 2.691 |
|                            | 40 Ca     |        3.478 | 3.453 | (0.003) |  3.44  | (0.004) | 3.463 |
|                            | 48 Ca     |        3.477 | 3.489 | (0.005) |  3.482 | (0.007) | 3.478 |
|                            | 68 Ni     |        3.887 | 3.871 | (0.004) |  3.864 | (0.004) | 3.887 |
|                            | 90 Zr     |        4.269 | 4.269 | (0.004) |  4.264 | (0.004) | 4.282 |
|                            | 132 Sn    |        4.709 | 4.707 | (0.005) |  4.708 | (0.005) | 4.732 |
|                            | 140 Ce    |        4.877 | 4.886 | (0.005) |  4.888 | (0.005) | 4.897 |
|                            | 208 Pb    |        5.501 | 5.504 | (0.006) |  5.507 | (0.008) | 5.54  |
| diffraction radius r d     | 16 O      |        2.764 | 2.767 | (0.006) |  2.756 | (0.011) | 2.707 |
|                            | 40 Ca     |        3.849 | 3.832 | (0.004) |  3.825 | (0.005) | 3.816 |
|                            | 48 Ca     |        3.963 | 3.973 | (0.005) |  3.971 | (0.010) | 3.953 |
|                            | 90 Zr     |        5.04  | 5.015 | (0.005) |  5.015 | (0.006) | 5.019 |
|                            | 208 Pb    |        6.776 | 6.787 | (0.008) |  6.799 | (0.009) | 6.829 |
| surface thickness σ        | 16 O      |        0.851 | 0.828 | (0.005) |  0.81  | (0.008) | 0.837 |
|                            | 40 Ca     |        0.968 | 0.964 | (0.005) |  0.946 | (0.009) | 0.968 |
|                            | 48 Ca     |        0.89  | 0.909 | (0.005) |  0.899 | (0.006) | 0.897 |
|                            | 90 Zr     |        0.957 | 1.001 | (0.005) |  0.993 | (0.007) | 0.994 |
|                            | 208 Pb    |        0.919 | 0.914 | (0.004) |  0.906 | (0.006) | 0.912 |
| valence neutron radius r n | 17 O      |        3.36  | 3.356 | (0.008) |  3.377 | (0.008) | 3.514 |
|                            | 41 Ca     |        3.99  | 3.999 | (0.006) |  3.99  | (0.011) | 4.116 |

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-56320887-d243-4980-8a34-0bc43b11b0e7.jpg)

Sn Ce Pb

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-573ba033-f619-411e-bf13-780eff72c217.jpg)

Pb

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-3c72813a-61d7-4637-8c7f-2da3198ffe67.jpg)

Pb

Fig. 4 Differences between theoretical and experimental charge radii r c , diffraction radii r d and surface thicknesses σ , respetively, of nuclei used in the fit of the model parameters. Theoretical uncertainties are indicated by the error bars.

more than 0.04 fm. The values of the valence neutron radii of 17 O and 41 Ca can also be found in Table 7. The deviations between theory and experiment are slightly larger than those of the quantities related to the charge form factor, but the agreement with the experimental data is still satisfactory, given the large theoretical uncertainties.

The differences between experimental and theoretical values of the spin-orbit splittings can be explored by studying Table 8. Considering the theoretical uncertainties of somewhat less than 0.1 MeV, there is a fair agreement of these data. The difference between the experimental spinorbit splittings of the ν 0 p and π p 0 levels in 16 O ( -0 140 MeV) is not reproduced by the theoretical . models as their differences are positive, 0 182 MeV and 0 170 MeV for the DDT and DDTC models, . . respectively. The spin-orbit splitting of the ν 0 p levels is larger than that of the π p 0 levels in contrast

to the experiment, in particular for the DD2 parameterisation. The spin-orbits splittings of the ν 1 p levels in 48 Ca are underestimated in the theoretical models in comparison to the experimental value. Larger differences between theory and experiment are also observed for the π d 1 levels in 132 Sn, and 208 Pb. In general, the results of the DD2 model show a larger deviation from the experiment as compared to the models with tensor couplings.

The contributions with tensor couplings in the EDF modify the strength of the spin-orbit interaction as compared to models without these couplings. Thus a modification in the singleparticle level structure can be expected. This could also affect the appearance of shell gaps. As an example, the semi-magic proton shell gaps in 140 Ce and 218 U at Z = 58 and Z = 92, respectively, are examined here. It has been observed that many relativistic EDFs predict large spurious shell gaps for these nuclei, see, e.g., Refs. [24, 47, 48]. In the DD2 model the π g 0 7 / 2 -π d 1 5 / 2 splitting in 140 Ce amounts to 3 054 MeV and the . π h 0 9 / 2 -π f 1 7 / 2 splitting in 218 U is 3 248 MeV. In fact, . reduced values of 2 288 MeV (2 543 MeV) and 2 665 MeV (2 999 MeV) are found for the . . . . Z = 58 and Z = 92 proton shell gaps in the DDT (DDTC) model.

The energies of the isoscalar monopole giant resonance are crucial, in the fit of the model parameters, to obtain a reasonable value of the incompressibility coefficient K of nuclear matter. Experimental data and theoretical values are given in Table 9 with theoretical uncertainties in the order of 0.2 MeV. In general, there is a good reproduction of the experimental values by theory, by the models with and without tensor couplings. Only the prediction for the resonance energy of 208 Pb is on the lower side as compared to the experiment.

In the determination of the model parameters, several exotic nuclei on both sides of the valley of stability were included in order to better constrain the isovector part of the effective in-medium interaction. In contrast to most earlier relativistic EDFs with density dependent couplings, three

Table 8 Experimental data of spin-orbit splittings [39, 40] in MeV of single-particle levels in nuclei used in the fit in comparison to theoretical values from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| nucleus   | level             | experiment        | DDT               | DDT                     | DDTC              | DDTC                    | DD2               |
|-----------|-------------------|-------------------|-------------------|-------------------------|-------------------|-------------------------|-------------------|
| 16 O      | ν 0 p π 0 p       | 6.180 6.320       | 6.127 5.945       | (0.095) (0.094)         | 6.231 6.061       | (0.214) (0.209)         | 6.867 6.819       |
| 48 Ca     | ν 1 p             | 2.020             | 1.841             | (0.068)                 | 1.849             | (0.075)                 | 1.510             |
| 90 Zr     | π 1 p             | 1.510             | 1.598             | (0.049)                 | 1.649             | (0.121)                 | 1.701             |
| 132 Sn    | π 1 d π 0 g       | 1.480 6.140       | 1.749 6.254       | (0.080) (0.097)         | 1.825 6.083       | (0.202) (0.169)         | 2.014 6.720       |
| 140 Ce    | ν 2 p             | 0.475             | 0.560             | (0.068)                 | 0.622             | (0.049)                 | 0.517             |
| 208 Pb    | ν 2 p π 1 d π 0 h | 0.898 1.330 5.550 | 1.066 1.680 5.712 | (0.065) (0.037) (0.088) | 0.993 1.732 5.527 | (0.037) (0.064) (0.085) | 0.930 1.868 6.203 |

Table 9 Experimental data of constraint isoscalar monopole giant resonance energies [41] in MeV of nuclei used in the fit in comparison to theoretical values from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| nucleus   | experiment   | DDT      | DDT     |   DDTC | DDTC    |    DD2 |
|-----------|--------------|----------|---------|--------|---------|--------|
| 90 Zr     | 17 . 810     | 17 . 597 | (0.214) | 17.714 | (0.214) | 17.584 |
| 116 Sn    | 15 . 900     | 16 . 112 | (0.193) | 16.189 | (0.196) | 16.085 |
| 144 Sm    | 15 . 250     | 15 . 158 | (0.212) | 15.22  | (0.318) | 15.205 |
| 208 Pb    | 14 . 180     | 13 . 489 | (0.188) | 13.519 | (0.232) | 13.417 |

independent isovector parameters are considered in the parameterisations that include tensor couplings. It is interesting to study, how this affects the predictions of the neutron skin thickness of nuclei, ∆ r skin = r ν -r π , that is obtained by the difference between the neutron and proton rms radii. The theoretical values for the nuclei in the fit are collected in Table 10. Nuclei with N = Z exhibit a small negative value of ∆ r skin that becomes even more negative for neutron deficient nuclei like 34 Ca and 48 Ni. On the other hand, neutron-rich nuclei can develop a substantial neutron skin. This is most clearly seen for 24 O, 28 O, 78 Ni, and 132 Sn. The theoretical uncertainties, however, become much larger for the most exotic nuclei as compared to stable nuclei. The neutron skin thicknesses of more exotic nuclei are generally larger for the DD2 models as compared to the DDT and DDTC models. Also stable nuclei such as 48 Ca and 208 Pb show a sizeable neutron skin thickness. It can be determined experimentally by studying the parity-violation asymmetry in high-energy electron scattering experiments. A comparison of the theoretical predictions in Table 10 with the experimental data of the CREX experiment [49] with ∆ r skin ( 48 Ca) = 0 121 . ± 0 025 (exp.) . ± 0 024 (model) fm . and those of the combined PREX experiments [50] with ∆ r skin ( 208 Pb) = 0 283 . ± 0 071 fm shows . that the theory can nicely reproduce the size of the neutron skin in 48 Ca within the uncertainties

Table 10 Theoretical neutron skin thickness ∆ r skin in fm of selected nuclei used in the fit from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| nucleus   | DDT       | DDT     | DDTC      | DDTC    | DD2       |
|-----------|-----------|---------|-----------|---------|-----------|
| 16 O      | - 0 . 031 | (0.001) | - 0 . 029 | (0.001) | - 0 . 029 |
| 24 O      | 0 . 450   | (0.014) | 0 . 480   | (0.018) | 0 . 595   |
| 28 O      | 0 . 622   | (0.019) | 0 . 650   | (0.042) | 0 . 769   |
| 34 Si     | 0 . 181   | (0.007) | 0 . 183   | (0.007) | 0 . 203   |
| 34 Ca     | - 0 . 299 | (0.007) | - 0 . 302 | (0.008) | - 0 . 323 |
| 40 Ca     | - 0 . 053 | (0.001) | - 0 . 052 | (0.001) | - 0 . 053 |
| 48 Ca     | 0 . 125   | (0.008) | 0 . 127   | (0.008) | 0 . 191   |
| 48 Ni     | - 0 . 251 | (0.007) | - 0 . 251 | (0.006) | - 0 . 316 |
| 56 Ni     | - 0 . 057 | (0.001) | - 0 . 057 | (0.001) | - 0 . 054 |
| 68 Ni     | 0 . 166   | (0.008) | 0 . 167   | (0.009  | 0 . 200   |
| 78 Ni     | 0 . 251   | (0.012) | 0 . 249   | (0.013) | 0 . 329   |
| 90 Zr     | 0 . 040   | (0.007) | 0 . 039   | (0.007) | 0 . 083   |
| 100 Sn    | - 0 . 083 | (0.001) | - 0 . 083 | (0.001) | - 0 . 082 |
| 132 Sn    | 0 . 191   | (0.011) | 0 . 189   | (0.013) | 0 . 261   |
| 140 Ce    | 0 . 086   | (0.009) | 0 . 083   | (0.010) | 0 . 150   |
| 208 Pb    | 0 . 134   | (0.010) | 0 . 130   | (0.013) | 0 . 198   |

Table 11 Experimental binding energy differences in MeV of mirror nuclei in comparison to theoretical values from different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| pair of mirror nuclei   |   experiment |    DDT | DDT     |   DDTC | DDTC    |    DD2 |
|-------------------------|--------------|--------|---------|--------|---------|--------|
| 41 Ca - 41 Sc           |        7.278 |  7.652 | (0.031) |  7.35  | (0.030) |  7.041 |
| 34 Si - 34 Ca           |       39.528 | 40.411 | (0.170) | 38.473 | (0.168) | 37.79  |
| 48 Ca - 48 Ni           |       68.673 | 69.805 | (0.272) | 67.355 | (0.244) | 64.837 |

but it predicts a considerably smaller neutron skin thickness in 208 Pb, an issue that is not solved so far. The neutron skin thickness can also be determined using hadronic probes using precisely known proton radii from electron scattering. A former analysis of 1 GeV proton scattering on 48 Ca [51] gave a neutron skin thickness of ∆ r skin = 0 10 . ± 0 03 fm only slightly smaller but consistent within . the errors with the CREX result. A more recent study of neutron radii using proton scattering [52] gave ranges of 0 056 fm . &lt; ∆ r skin ( 48 Ca) &lt; . 0 102 fm and 0 083 fm . &lt; ∆ r skin ( 208 Pb) &lt; . 0 111 fm that point to rather small neutron thicknesses of these nuclei.

The inclusion of pairs of mirror nuclei in the fit also allows to better constrain the mirror displacement energies. The binding energy differences of three pairs of mirror nuclei are given in Table 11. There is a good agreement between experiment and theory with theoretical results of the DDT model that are somewhat larger than the experimental values. The differences between experimental and theoretical values are 0 374, 0 883, and 1 132 MeV for the pairs with mass number . . . A = 41, 34, and 48 for the DDT model, respectively, thus overestimating the binding energy differences. In contrast, they are systematically smaller with 0 072, . -1 055, and . -1 318 MeV for . the DDTC model and -0 237, . -1 738, and . -3 836 MeV for the DD2 model, respectively, where . the Z Z ( -1) factor in the Coulomb field was introduced. In addition, the neutron-deficient nuclei were not considered in the fit of the parameters in the latter parameterisation. Thus the parameter set DDT gives the best description on the average. The theoretical uncertainties increase with the isospin difference of the pair of nuclei.

## 4.2 Nuclear matter

The energy per nucleon E/A , Eq. (58), in homogeneneous nuclear matter is a function of the baryon density n b = n ν + n π and the isospin asymmetry α = ( n ν -n π ) /n b . For α = 0, the energy per nucleon of symmetric nuclear matter E 0 = E/A n , ( b 0) is obtained. It is depicted in Figure 5 for the parameterisations DDT, DDTC, and DD2. The minimum of E 0 for these equations of state defines the saturation density n sat , see Table 12. The region of negative slope, i.e., negative pressure, signals an instability of homogeneous matter that is related to the occurrence of the liquid-gas phase transition. At high densities, the predictions for the energy per nucleon become

less certain as the widening of the uncertainty bands of the DDT and, significantly larger, DDTC models shows. The energy per nucleon of pure neutron matter, i.e., E/A for α = 1, is also shown in Figure 5. It increases monotonically with the baryon density.

The curves of the DDT and DD2 parameterisations in Fig. 5 look very similar at first sight. But there are differences in detail, in particular at lower densities. The curves of the parameterisation DDTC follow along the lower limit of the uncertainty bands of the DDT model, hence signaling a softer equation of state. The neutron matter equation of state of the DDT model shows a peculiar feature at the lowest densities. This can be seen in the inset of Figure 5. The energy per nucleon does not approach zero as smoothly as in the DD2 and DDTC models. This is related to the specific density dependence of the ω -meson coupling Γ ω , see panel (b) of Figure 2. The increase of the coupling at very small densities causes an additional repulsion, which also applies for symmetric nuclear matter at very low densities (not shown in the inset).

Assuming a quadratic dependence of E/A on α for given baryon density n b , the symmetry energy (59) can be obtained as a difference E sym ( n b ) = E/A n , ( b 1) -E/A n , ( b 0). This quantity is shown in Figure 6 as function of n b . Particularly noteworthy is the change of the slope above saturation density, for the DDT and DDTC models, with a widening of the uncertainty bands.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-e4d74f66-3807-4a2d-bc62-952d66718cdc.jpg)

b

Fig. 5 Equation of state of symmetric nuclear matter and pure neutron matter, respectively, with theoretical uncertainty bands for the DDT and DDTC models. The inset shows the low-energy equation of state of pure neutron matter.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-4dec1b99-09e4-4337-8bd6-8a6e457d5be7.jpg)

b

Fig. 6 Symmetry energy E sym as function of the baryon density n b with theoretical uncertainty bands for the DDT and DDTC models. The inset shows the low density part with a comparison to experimental data, see text for details.

At densities below saturation, symmetry energy is well constrained by experiments and ab-initio calculations of pure neutron matter, see, e.g., Refs. [53, 54] and references therein. At densities above saturation, the predictions for the equation of state of nuclear matter are less precise as indicated by the increasing uncertainty bands.

The E sym curves of the DDT, DDTC, and DD2 model cross at a density below n sat . It has been observed before that fitting model parameters of EDFs to properties of nuclei effectively determines the symmetry energy around this subsaturation density, see, e.g., Ref. [55]. A comparison of the lowdensity symmetry energy of the DDT, DDTC, and DD2 models with experimental data is depicted in the inset of Figure 6. These data come from the determination of the dipole polarizability of 208 Pb (violet), the analysis of EDFs (green, light blue), and isobaric analog states (orange), see Ref. [56] and references therein. The symmetry energy of the parameterisations DDT and DDTC lie slightly above these points, corresponding to the small slope parameter L , whereas the DD2 curve matches the data somewhat better.

Characteristic nuclear matter parameters are given in Table 12 for the models DDT, DDTC, and DD2. Most remarkable are the increased effective mass and the reduced incompressibility coefficient of the models with tensor couplings (DDT and DDTC) in comparison to that without

Table 12 Characteristic nuclear matter parameters for different parameterisations of the relativistic EDF. Uncertainties of the quantities are given in brackets.

| quantity   | unit   | DDT        | DDT         | DDTC       | DDTC        | DD2        |
|------------|--------|------------|-------------|------------|-------------|------------|
| n sat      | fm - 3 | 0 . 15493  | (0 . 00076) | 0 . 15304  | (0 . 00095) | 0 . 14907  |
| M ∗ sat    | M nuc  | 0 . 65729  | (0 . 01315) | 0 . 66305  | (0 . 03803) | 0 . 56252  |
| B sat      | MeV    | 16 . 226   | (0 . 025)   | 16 . 226   | (0 . 040)   | 16 . 022   |
| K          | MeV    | 229 . 20   | (7 . 99)    | 215 . 05   | (7 . 71)    | 242 . 72   |
| Q          | MeV    | 88 . 565   | (230 . 061) | - 291 . 86 | (763 . 87)  | 168 . 77   |
| J          | MeV    | 31 . 209   | (0 . 498)   | 30 . 379   | (0 . 568)   | 31 . 670   |
| L          | MeV    | 34 . 236   | (4 . 332)   | 33 . 317   | (7 . 847)   | 55 . 033   |
| K sym      | MeV    | - 69 . 538 | (15 . 454)  | - 75 . 562 | (26 . 416)  | - 93 . 226 |
| K τ,v      | MeV    | - 288 . 18 | (54 . 09)   | - 230 . 25 | (131 . 21)  | - 461 . 69 |

these couplings (DD2). Most existing relativistic models show values of K that are rather high when compared with non-relativistic EDFs, cf. Refs. [3, 15]. The tensor models are an exception and they are also characterized by rather large saturation densities. The higher-order coefficients Q and K sym show considerable differences and are not very well determined by the parameter fits as they have large uncertainties. The smaller values of K and Q of the DDTC model indicate a softer equation of state as compared to the DDT model as seen already in Figure 5. The DDT and DDTC models also have symmetry energy slope parameters that are smaller than the value of L of the DD2 model. This difference correlates with the smaller neutron skin thicknesses for exotic nuclei as discussed in Subsection 4.1.

## 4.3 Neutron stars

The properties of a non-rotating, spherically symmetric neutron star, are obtained by solving the first-order Tolman-Oppenheimer-Volkoff equation [64, 65]

$$\frac { d P } { d r } = - G \, \frac { M ( r ) \varepsilon ( r ) } { r ^ { 2 } } \left [ 1 + \frac { P ( r ) } { \varepsilon ( r ) } \right ] \left [ 1 + \frac { 4 \pi r ^ { 3 } P ( r ) } { M ( r ) } \right ] \left [ 1 - \frac { 2 G M ( r ) } { r } \right ] ^ { - 1 } \,,$$

with the mass

$$M ( r ) = 4 \pi \int _ { 0 } ^ { r } d r ^ { \prime } \left ( r ^ { \prime } \right ) ^ { 2 } \varepsilon ( r ^ { \prime } ) \,, \\ \dim \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m n } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m } \, \mathbb { m m }$$

inside the volume of radius r . The factors in square brackets are the general relativistic corrections to the classical expression for the change of the pressure with the radius. Assuming a particular density n central in the center of the star, the first order differential equation (98) can be solved by an outward integration until the pressure P r ( ) becomes sufficiently small. This defines the stellar radius R and the total mass M R ( ). The essential ingredient in the calculation is the equation of state of charge-neutral matter in β equilibrium. It is obtained by adding the leptonic contributions (electrons and muons) to the nuclear equation of state and determining the isospin asymmetry α for given baryon density n b , by requiring vanishing lepton chemical potentials µ e + µ q , and µ µ + µ q ,

Fig. 7 Mass-radius relation of neutron stars for the parameterisations DDTC (blue), DDT (red), with uncertainty bands, and DD2 (black) of the relativistic EDF. Open circles indicate neutron stars where the central density is an integer multiple n central = kn sat of the saturation density starting with k = 1 at the lowest masses. (The point with k = 1 of the DDTC parameterisation is at a radius R = 21 62 fm outside the figure.) Yellow lines give the masses . of heavy neutron stars [57-59] with errors indicated on the left. Masses and radii of neutron stars from NICER observations [60-63] are shown in green color.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00425v2-figure-aa27b909-da21-4bc4-be45-f8deb01163de.jpg)

with charge chemical potential µ q = µ π -µ ν . At densities below nuclear saturation, an appropriate crust equation of state has to be used as the matter is no longer homogeneous. Here, the crust description of the generalized relativistic density functional with the DD2 parameterisation of the couplings [43] is used where the masses of the nuclei in the crust are taken from the Duflo-Zuker model DZ31 [66], if not available from the atomic mass evaluation AME2020 [32].

The mass-radius relations of neutron stars for the different parameterisations are depicted in Figure 7. Open circles on the curves denote the properties of neutron stars with central densities that are an integer multiple of the saturation density n sat . The curves show the well-known inverse S-shaped rise of the mass with increasing central density and the corresponding variation of the radius. Ordinary neutron stars with masses of 1.4 times the solar mass M sol have central densities that are typically between twice and three times the saturation density. At certain central densities n (max) central , the most massive neutron star configurations are obtained. They have central densities of five to six times the saturation density. Neutron stars with higher central density are not stable any more and this last branch of the curve will not be populated. The maximum masses of the shown models are always above the largest masses of observed neutron stars [57-59] since their equation of state is sufficiently stiff. Masses and radii are avalaible from NICER observations [6063], but these data have large uncertainties. All models are consistent with these observations. Characteristic properties of neutron stars of the models in this work are collected in Table 13.

Table 13 Properties of neutron stars for the different parameterisations of the relativistic EDF.

| quantity        | unit   |      DDT |     DDTC |      DD2 |
|-----------------|--------|----------|----------|----------|
| R 1 . 4         | km     | 12.744   | 12.566   | 13.172   |
| M max           | M sol  |  2.43    |  2.371   |  2.417   |
| R max           | km     | 11.739   | 11.504   | 11.869   |
| n (max) central | fm - 3 |  0.85785 |  0.89663 |  0.85114 |

The properties of the maximum mass stars are very similar in all models but the radii of the DD2 model are systematically larger that those of the parameterisation DDT and DDTC. The latter model predicts the lowest maximum mass, a higher central density at this mass, and smaller radii as compared to the other models, corresponding to the softer equation of state.

## 5 Summary and Conclusions

In this work, an extension of a relativistic EDF is presented that includes tensor couplings between the Lorentz vector mesons, ω and ρ , and nucleons in addition to the standard minimal mesonnucleon couplings. The latter are assumed to depend on the baryon density. The field equation are derived in the mean-field approximation from a covariant Lagrangian density and solved selfconsistently for nuclei and nuclear matter.

The parameters of this phenomenological approach are found from a fit to properties of 20 selected nuclei that comprise binding energies, charge and diffraction radii, surface thicknesses, spin-orbit splittings, valence neutron radii, and constraint isoscalar monopole giant resonance energies. The model uncertainties of these quantities are determined self-consistently during the fitting procedure. The model predictions are compared to experimental data and given with their theoretical uncertainties, which were, to our knowledge, not specified in earlier parameterisations of relativistic EDFs. Two new sets of parameters, DDT and DDTC, are developed that differ in the treatment of the Coulomb interaction.

The inclusion of tensor couplings in the relativistic EDF leads to a considerable improvement in the description of nuclear properties, in particular their binding energies. Quantities related to the nuclear charge form factor show a similar quality as those of an earlier parameterisation (DD2) without tensor couplings. The binding energy differences of mirror nuclei are also closer to the experimental data than in the DD2 model, that did not consider exotic, neutron deficient nuclei in the determination of the parameters. Predictions for the neutron skin thickess of bound and unstable nuclei are made.

The equations of state for symmetric nuclear matter and pure neutron matter are calculated with the new parameters. They look similar as those of the DD2 model but show some specific

differences. These are most apparent by comparing characteristic nuclear matter parameters. The new parameterisations exhibit larger saturation densities and effective nucleon masses, and smaller incompressibility coefficients and symmetry energy slope parameters.

The description of charge-neutral nuclear matter in β equilibrium with the extended EDFs allows to calculate properties of neutron stars and, in particular, their mass-radius relation. The results are similar as for the DD2 model but with slightly smaller radii. Due to the stiffness of the EOS at high densities, the maximum masses are predicted to lie above the data of observed heavy neutron stars. The obtained mass-radius relations are consistent with results from NICER measurements that still have substantial uncertainties, in particular for the radii.

EDFs with the new parameterisations can be applied in further calculations to explore their quality in the description of finite nuclei. The study of deformed nuclei is of particular interest as it allows to investigate additional observables, e.g., higher multipole moments of the density distributions, or the variation of energies of single-particle states in isotopic and isotonic chains of nuclei. The tensor couplings could affect the position of the levels in the spectrum and give rise to the appearance of new shell gaps or their suppression. Effects of pairing have to be included in the model too when open-shell nuclei are considered. However, additional work with proper numerical methods is required to include deformation and pairing in the calculations. This is beyond the scope of the present work. The new models can also be used to prepare EOS tables of stellar matter at finite temperature for larger ranges in isospin asymmetry and density, in order to be used in simulations of core-collapse supernovae and neutron-star mergers [67]. For this, the calculations have to be extended by including additional degrees of freedom, e.g., light clusters at subsaturation densities [43, 68], or hyperons at high densities, see, e.g., Ref. [69].

## Acknowledgments

S. Typel acknowledges the support by the Cyclotron Institute, Texas A&amp;M University. S. Shlomo is partially supported by the US Department of Energy under Grant no. DE-FG03-93ER-40773.

## Statements and Declarations

Open Access funding enabled and organized by Projekt DEAL. The authors have no relevant financial or non-financial interests to disclose. Data sets generated during the current study are available from the corresponding author on reasonable request.

## Appendix A Transformation between parameters of nuclear matter and of coupling functions

The parameters of the coupling functions (67) can be derived from the characteristic nuclear matter parameters. To simplify this conversion, the average nucleon mass M nuc = ( M π + M / ν ) 2 is used instead of proton and neutron masses in the full calculation. However, this does not pose any problem.

- 1. The reference density /rho1 ref is chosen as the baryon density at saturation n sat and defines the

Fermi momentum

$$k _ { \text{sat} } ^ { * } & = \left ( \frac { 3 \pi ^ { 2 } } { 2 } n _ { \text{sat} } \right ) ^ { 1 / 3 } \,. & \quad \quad \quad$$

- 2. The Dirac effective mass at saturation M ∗ sat gives the effective chemical potential

$$\mu _ { \text{sat} } ^ { * } & = \sqrt { ( k _ { \text{sat} } ^ { * } ) ^ { 2 } + ( M _ { \text{sat} } ^ { * } ) ^ { 2 } }$$

and the scalar density at saturation

$$n _ { \text{sat} } ^ { ( s ) } = \frac { M _ { \text{sat} } ^ { * } } { \pi ^ { 2 } } \left [ k _ { \text{sat} } ^ { * } \mu _ { \$$

This determines the scalar coupling at saturation

$$C _ { \sigma } & = \frac { M _ { \text{nuc} } - M _ { \text{sat} } ^ { * } } { n _ { \text{sat} } ^ { ( s ) } } \,,$$

and the auxiliary quantity

$$\i v y & & f _ { \text{sat} } = 3 \left ( \frac { n _ { \text{sat} } ^ { ( s ) } } { M _ { \text{sat} } ^ { * } } - \frac { n _ { \$$

that will be used below.

- 3. The energy per nucleon at saturation ( E/A ) sat gives the energy density at saturation

$$\varepsilon _ { \text{sat} } = \left [ \left ( E / A \right ) _ { \text{sat} } + M _ { \text{nuc} } \right ] n _ { \text{sat} } \,,$$

and the vector coupling at saturation

$$C _ { \omega } = \frac { 2 } { n _ { \text{sat} } ^ { 2 } } \left [ \varepsilon _ { \text{sat} } - \frac { 1 } { 4 } \left ( 3 \mu _ { \text{sat} } ^ { * } n _ { \text{sat} } + M _ { \text{sat} } ^ { * } n _ { \text{sat} } ^ { ( s ) } \right ) - \frac { 1 } { 2 } C _ { \sigma } \left ( n _ { \text{sat} } ^ { ( s ) } \right ) ^ { 2 } \right ] \,.$$

Furthermore, the chemical potential

$$\mu _ { \text{sat} } = \frac { \varepsilon _ { \text{sat} } } { n _ { \text{sat} } } = ( E / A ) _ { \text{sat} } + M _ { \text{sat} } \,,$$

and the rearrangement potential

$$V ^ { ( R ) } = \mu _ { \text{sat} } - \mu _ { \text{sat} } ^ { * } - C _ { \omega } n _ { \text{sat} } \,,$$

are found.

## 4. The ratio r 10 determines the derivatives

$$C _ { \omega } ^ { \prime } & = 2 r _ { 1 0 } \frac { C _ { \omega } } { n _ { \text{sat} } } \,,$$

and

$$C _ { \sigma } ^ { \prime } & = \frac { 1 } { \left ( n _ { \text{sat} } ^ { ( s ) } \right ) ^ { 2 } } \left ( C _ { \omega } ^ { \prime } n _ { \text{sat} } ^ { 2 } - 2 V ^ { ( R ) } \right ) \,, & ( A 1 1 ) \\ \intertext { n, the second and third auxiliary variables }$$

at saturation. Then, the second and third auxiliary variables

$$g _ { \text{sat} } = \frac { \frac { M _ { \text{sat} } ^ { * } } { \mu _ { \text{sat} } ^ { * } } - C _ { \sigma } ^ { \prime } n _ { \text{sat} } ^ { ( s ) } f _ { \text{sat} } } { 1 + C _ { \sigma } f _ { \text{sat} } } \,,$$

$$h _ { \text{sat} } = - \frac { C _ { \sigma } ^ { \prime } n _ { \text{sat} } ^ { ( s ) } + C _ { \sigma } \frac { M _ { \text{sat} } ^ { * } } { \mu _ { \text{sat} } ^ { * } } } { 1 + C _ { \sigma } f _ { \text{sat} } } \,,$$

can be calculated.

- 5. The ratio r 21 is used to obtain the second derivative

$$C _ { \omega } ^ { \prime \prime } & = 2 r _ { 1 0 } ( r _ { 1 0 } + r _ { 2 1 } ) \frac { C _ { \omega } } { n _ { \text{sat} } ^ { 2 } } \,,$$

for the ω meson at saturation.

- 6. The incompressibility coefficient K gives the second derivative,

$$C _ { \sigma } ^ { \prime \prime } & = \frac { 2 } { n _ { \text{sat} } \left ( n _ { \text{sat} } ^ { ( s ) } \right ) ^ { 2 } } \left [ \frac { 1 } { 3 } \left ( \mu _ { \text{sat} } ^ { * } - M _ { \text{sat} } ^ { * } g _ { \text{sat} } \right ) - \frac { K } { 9 } + \left ( C _ { \omega } + 2 C _ { \omega } ^ { \prime } n _ { \text{sat} } + \frac { 1 } { 2 } C _ { \omega } ^ { \prime \prime } n _ { \text{sat} } ^ { 2 } \right ) n _ { \text{sat} } \\ & \quad - \left ( C _ { \sigma } + C _ { \sigma } ^ { \prime } n _ { \text{sat} } \right ) n _ { \text{sat} } ^ { ( s ) } g _ { \text{sat} } - C _ { \sigma } ^ { \prime } \left ( n _ { \text{sat} } ^ { ( s ) } \right ) ^ { 2 } \right ] \,,$$

and,

for the σ meson.

- 7. The symmetry energy at saturation J is used to calculate the coupling of the ρ meson at saturation,

where

$$\Gamma _ { \rho } ^ { ( 0 ) } & = \sqrt { C _ { \rho } } M _ { \rho } \,,$$

$$C _ { \rho } = \frac { 2 } { n _ { \text{sat} } } \left [ J - \frac { ( k _ { \text{sat} } ^ { * } ) ^ { 2 } } { 6 \mu _ { \text{sat} } ^ { * } } \right ] \,.$$

8. The slope parameter L of the symmetry energy at saturation gives the derivative

$$C _ { \rho } ^ { \prime } = \frac { 2 } { 3 n _ { \text{sat} } ^ { 2 } } \left \{ L - 3 J - \frac { ( k _ { \text{sat} } ^ { * } ) ^ { 2 } } { 3 \mu _ { \text{sat} } ^ { * } } \left [ \frac { 1 } { 2 } \left ( \frac { M _ { \text{sat} } ^ { * } } { \mu _ { \text{sat} } ^ { * } } \right ) ^ { 2 } \left ( 1 - 3 \frac { n _ { \text{sat} } } { M _ { \text{sat} } ^ { * } } h _ { \text{sat} } \right ) - 1 \right ] \right \} \,, \quad ( A 1 8 )$$

for the ρ meson at saturation. Then the parameter

$$a _ { \rho } = - \frac { 1 } { 2 } \frac { C _ { \rho } ^ { \prime } } { C _ { \rho } } n _ { \text{sat} } \,,$$

in the function f ρ ( x ) is found.

- 9. The values of C j and their derivatives determine the couplings and first derivatives

$$\Gamma _ { j } ^ { ( 0 ) } = \sqrt { C _ { j } } M _ { j } \,, \quad \text{and} \quad \Gamma _ { j } ^ { \prime } = \frac { 1 } { 2 } \Gamma _ { j } ^ { ( 0 ) } \frac { C _ { j } ^ { \prime } } { C _ { j } } \,,$$

as well as the second derivatives

$$\Gamma _ { j } ^ { \prime \prime } = \frac { 1 } { 2 } \Gamma _ { j } ^ { ( 0 ) } \left [ \frac { C _ { j } ^ { \prime \prime } } { C _ { j } } - \frac { 1 } { 2 } \left ( \frac { C _ { j } ^ { \prime } } { C _ { j } } \right ) ^ { 2 } \right ] \,,$$

for j = σ and ω . Then also the derivatives

$$f _ { j } ^ { \prime } ( 1 ) = \frac { \Gamma _ { j } ^ { \prime } } { \Gamma _ { j } ^ { ( 0 ) } } \,, \quad \text{and} \quad f _ { j } ^ { \prime \prime } ( 1 ) = \frac { \Gamma _ { j } ^ { \prime \prime } } { \Gamma _ { j } ^ { ( 0 ) } } \,,$$

of the functions f j ( x ) are known.

10. Finally, the parameters of the rational functions f j ( x ) for the isoscalar mesons can be found with help of the explicit expressions;

$$f _ { j } ^ { \prime } ( x ) & = 2 a _ { j } ( b _ { j } - c _ { j } ) \frac { x + d _ { j } } { [ 1 + c _ { j } ( x + d _ { j } ) ^ { 2 } ] ^ { 2 } } \,, & ( A 2 3 )$$

$$f _ { j } ^ { \prime \prime } ( x ) \, = \, 2 a _ { j } ( b _ { j } - c _ { j } ) \frac { \dot { \lambda } ^ { \lambda } \dot { \lambda } ^ { \lambda } } { [ 1 + c _ { j } ( x + d _ { j } ) ^ { 2 } ] ^ { 3 } } \,,$$

for the first and second derivatives of the functions f j ( x ). The ratio,

$$r _ { 2 1 } = \frac { f _ { \omega } ^ { \prime \prime } ( 1 ) } { f ( 1 ) } = \frac { d _ { j } ^ { 2 } - ( 1 + d _ { j } ) ^ { 2 } } { [ d _ { j } ^ { 2 } + ( 1 + d _ { j } ) ^ { 2 } / 3 ] ( 1 + d _ { j } ) } \,,$$

is a function of d j only because of the constraint 1 = 3 c d j 2 j that follows from f ′′ j (0) = 0. Thus the value of d j can be obtained by a simple numerical search or iteration. Then c j = 1 (3 / d 2 j ) is known as well and b j is calculated from

$$b _ { j } = \frac { c _ { j } z + r _ { 1 0 } } { z - ( 1 + d _ { j } ) ^ { 2 } r _ { 1 0 } } \,,$$

with the auxiliary quantity

$$z = \frac { 2 ( 1 + d _ { j } ) } { 1 + c _ { j } ( 1 + d _ { j } ) ^ { 2 } } \,.$$

Finally, the normalization condition f j (1) = 1 gives

$$a _ { j } = \frac { 1 + c _ { j } ( 1 + d _ { j } ) ^ { 2 } } { 1 + b _ { j } ( 1 + d _ { j } ) ^ { 2 } } \,.$$

## References

- [1] B.K. Agrawal, S. Shlomo, V.K. Au, Determination of the parameters of a Skyrme type effective interaction using the simulated annealing approach. Phys. Rev. C 72 , 014310 (2005). https://doi.org/10.1103/PhysRevC.72.014310. arXiv:nucl-th/0505071
- [2] J.R. Stone, P.G. Reinhard, The Skyrme Interaction in finite nuclei and nuclear matter. Prog. Part. Nucl. Phys. 58 , 587-657 (2007). https://doi.org/10.1016/j.ppnp.2006.07.001. arXiv:nucl-th/0607002
- [3] M. Dutra, O. Lourenco, J.S. Sa Martins, A. Delfino, J.R. Stone, P.D. Stevenson, Skyrme Interaction and Nuclear Matter Constraints. Phys. Rev. C 85 , 035201 (2012). https://doi.org/10.1103/PhysRevC.85.035201. arXiv:1202.3902 [nucl-th]
- [4] R. Sellahewa, A. Rios, Isovector properties of the Gogny interaction. Phys. Rev. C 90 (5), 054327 (2014). https://doi.org/10.1103/PhysRevC.90.054327. arXiv:1407.8138 [nucl-th]
- [5] P.G. Reinhard, M. Rufa, J. Maruhn, W. Greiner, J. Friedrich, Nuclear Ground State Properties in a Relativistic Meson Field Theory. Z. Phys. A 323 , 13-25 (1986)

- [6] M. Rufa, P.G. Reinhard, J.A. Maruhn, W. Greiner, M.R. Strayer, Optimal parametrization for the relativistic mean-field model of the nucleus. Phys. Rev. C 38 , 390-409 (1988). https://doi.org/10.1103/PhysRevC.38.390
- [7] P.G. Reinhard, The Relativistic Mean Field Description of Nuclei and Nuclear Dynamics. Rept. Prog. Phys. 52 , 439 (1989). https://doi.org/10.1088/0034-4885/52/4/002
- [8] B.D. Serot, Quantum hadrodynamics. Rept. Prog. Phys. 55 , 1855-1946 (1992). https://doi.org/10.1088/0034-4885/55/11/001
- [9] C. Fuchs, H. Lenske, H.H. Wolter, Density dependent hadron field theory. Phys. Rev. C 52 , 3043-3060 (1995). https://doi.org/10.1103/PhysRevC.52.3043. arXiv:nucl-th/9507044
- [10] P. Ring, Relativistic mean field in finite nuclei. Prog. Part. Nucl. Phys. 37 , 193-263 (1996). https://doi.org/10.1016/0146-6410(96)00054-3
- [11] S. Typel, H.H. Wolter, Relativistic mean field calculations with density dependent meson nucleon coupling. Nucl. Phys. A 656 , 331-364 (1999). https://doi.org/10.1016/S0375-9474(99)00310-3
- [12] D. Vretenar, A.V. Afanasjev, G.A. Lalazissis, P. Ring, Relativistic Hartree Bogoliubov theory: static and dynamic aspects of exotic nuclear structure. Phys. Rept. 409 , 101-259 (2005). https://doi.org/10.1016/j.physrep.2004.10.001
- [13] J. Meng, H. Toki, S.G. Zhou, S.Q. Zhang, W.H. Long, L.S. Geng, Relativistic Continuum Hartree Bogoliubov theory for ground state properties of exotic nuclei. Prog. Part. Nucl. Phys. 57 , 470-563 (2006). https://doi.org/10.1016/j.ppnp.2005.06.001. arXiv:nucl-th/0508020
- [14] T. Niksic, D. Vretenar, P. Ring, Relativistic Nuclear Energy Density Functionals: Mean-Field and Beyond. Prog. Part. Nucl. Phys. 66 , 519-548 (2011). https://doi.org/10.1016/j.ppnp.2011.01.055. arXiv:1102.4193 [nucl-th]
- [15] M. Dutra, O. Louren¸ co, S.S. Avancini, B.V. Carlson, A. Delfino, D.P. Menezes, C. Providˆncia, e S. Typel, J.R. Stone, Relativistic Mean-Field Hadronic Models under Nuclear Matter Constraints. Phys. Rev. C 90 (5), 055203 (2014). https://doi.org/10.1103/PhysRevC.90.055203. arXiv:1405.3633 [nucl-th]
- [16] B.A. Nikolaus, T. Hoch, D.G. Madland, Nuclear ground state properties in a relativistic point coupling model. Phys. Rev. C 46 , 1757-1781 (1992).

## https://doi.org/10.1103/PhysRevC.46.1757

- [17] B¨ urvenich, T. and Madland, D. G. and Maruhn, J. A. and Reinhard, P. G., Nuclear ground state observables and QCD scaling in a refined relativistic point coupling model. Phys. Rev. C 65 , 044308 (2002). https://doi.org/10.1103/PhysRevC.65.044308. arXiv:nucl-th/0111012
- [18] T. Niksic, D. Vretenar, P. Ring, Relativistic Nuclear Energy Density Functionals: Adjusting parameters to binding energies. Phys. Rev. C 78 , 034318 (2008). https://doi.org/10.1103/PhysRevC.78.034318. arXiv:0809.1375 [nucl-th]
- [19] P.W. Zhao, Z.P. Li, J.M. Yao, J. Meng, New parametrization for the nuclear covariant energy density functional with point-coupling interaction. Phys. Rev. C 82 , 054319 (2010). https://doi.org/10.1103/PhysRevC.82.054319. arXiv:1002.1789 [nucl-th]
- [20] J.D. Walecka, A Theory of highly condensed matter. Annals Phys. 83 , 491-529 (1974). https://doi.org/10.1016/0003-4916(74)90208-5
- [21] R.J. Furnstahl, J.J. Rusnak, B.D. Serot, The Nuclear spin orbit force in chiral effective field theories. Nucl. Phys. A 632 , 607-623 (1998). https://doi.org/10.1016/S0375-9474(98)00004-9. arXiv:nucl-th/9709064
- [22] G.j. Mao, Effect of tensor couplings in a relativistic Hartree approach for finite nuclei. Phys. Rev. C 67 , 044318 (2003). https://doi.org/10.1103/PhysRevC.67.044318. arXiv:nucl-th/0211034
- [23] W.H. Long, H. Sagawa, N. Van Giai, J. Meng, Shell Structure and rho-Tensor Correlations in Density-Dependent Relativistic Hartree-Fock theory. Phys. Rev. C 76 , 034314 (2007). https://doi.org/10.1103/PhysRevC.76.034314. arXiv:0706.3497 [nucl-th]
- [24] Q. Zhao, Z. Ren, P. Zhao, J. Meng, Covariant density functional theory with localized exchange terms. Phys. Rev. C 106 (3), 034315 (2022). https://doi.org/10.1103/PhysRevC.106.034315. arXiv:2207.01764 [nucl-th]
- [25] S. Typel, D. Alvear Terrero, Parametrisations of relativistic energy density functionals with tensor couplings. Eur. Phys. J. A 56 (6), 160 (2020). https://doi.org/10.1140/epja/s10050-020-00172-2. arXiv:2003.02085 [nucl-th]
- [26] F. Mercier, J.P. Ebran, E. Khan, Covariant energy density functionals with and without tensor couplings at the Hartree-Bogoliubov level. Phys. Rev. C 107 (3), 034309 (2023).

https://doi.org/10.1103/PhysRevC.107.034309. arXiv:2210.11142 [nucl-th]

- [27] M. Salinas, J. Piekarewicz, Impact of tensor couplings with scalar mixing on covariant energy density functionals. Phys. Rev. C 109 (4), 045807 (2024). https://doi.org/10.1103/PhysRevC.109.045807. arXiv:2312.13474 [nucl-th]
- [28] S. Typel, Lagrange-Mesh Method for Deformed Nuclei With Relativistic Energy Density Functionals. Front. in Phys. 6 , 73 (2018). https://doi.org/10.3389/fphy.2018.00073
- [29] N.M. Hugenholtz, L. van Hove, A theorem on the single particle energy in a Fermi gas with interaction. Physica 24 , 363-376 (1958). https://doi.org/10.1016/S0031-8914(58)95281-9
- [30] I. Vidana, C. Providencia, A. Polls, A. Rios, Density dependence of the nuclear symmetry energy: A Microscopic perspective. Phys. Rev. C 80 , 045806 (2009). https://doi.org/10.1103/PhysRevC.80.045806. arXiv:0907.1165 [nucl-th]
- [31] G.G. Simon, C. Schmitt, F. Borkowski, V.H. Walther, Absolute electron Proton Cross-Sections at Low Momentum Transfer Measured with a High Pressure Gas Target System. Nucl. Phys. A 333 , 381-391 (1980). https://doi.org/10.1016/0375-9474(80)90104-9
- [32] M. Wang, W.J. Huang, F.G. Kondev, G. Audi, S. Naimi, The AME 2020 atomic mass evaluation (II). Tables, graphs and references. Chin. Phys. C 45 (3), 030003 (2021). https://doi.org/10.1088/1674-1137/abddaf
- [33] H. De Vries, C.W. De Jager, C. De Vries, Nuclear charge and magnetization density distribution parameters from elastic electron scattering. Atom. Data Nucl. Data Tabl. 36 , 495-536 (1987). https://doi.org/10.1016/0092-640X(87)90013-1
- [34] G. Fricke, C. Bernhardt, K. Heilig, L.A. Schaller, L. Schellenberg, E.B. Shera, C.W. de Jager, Nuclear Ground State Charge Radii from Electromagnetic Interactions. Atom. Data Nucl. Data Tabl. 60 , 177-285 (1995). https://doi.org/10.1006/adnd.1995.1007
- [35] E.G. Nadjakov, K.P. Marinova, Y.P. Gangrsky, Systematics of Nuclear Charge Radii. Atom. Data Nucl. Data Tabl. 56 , 133-157 (1994). https://doi.org/10.1006/adnd.1994.1004
- [36] I. Angeli, K.P. Marinova, Table of experimental nuclear ground state charge radii: An update. Atom. Data Nucl. Data Tabl. 99 (1), 69-95 (2013). https://doi.org/10.1016/j.adt.2011.12.006

- [37] N. Kalantar-Nayestanaki, et al., Magnetic structure of 17 O at high momentum. Phys. Rev. Lett. 60 , 1707-1710 (1988). https://doi.org/10.1103/PhysRevLett.60.1707
- [38] S. Platchkov, A. Amroun, P. Bricault, J.M. Cavedon, P.K.A. de Witt Huberts, P. Dreux, B. Frois, C.D. Goodman, D. Goutte, J. Martino, V. Meot, G.A. Peterson, X.H. Phan, S. Raman, I. Sick, Measurement of the 1 f 7 2 -neutron-orbit radius in 41 Ca. Phys. Rev. Lett. 61 , 1465-1468 (1988). https://doi.org/10.1103/PhysRevLett.61.1465. URL https://link.aps.org/doi/10.1103/PhysRevLett.61.1465
- [39] B.N.L. National Nuclear Data Center. Evaluated nuclear structure data files. URL https://www.nndc.bnl.gov/ensdf/
- [40] S. Typel, Relativistic model for nuclear matter and atomic nuclei with momentum-dependent self-energies. Phys. Rev. C 71 , 064301 (2005). https://doi.org/10.1103/PhysRevC.71.064301. arXiv:nucl-th/0501056
- [41] D.H. Youngblood, H.L. Clark, Y.W. Lui, Incompressibility of Nuclear Matter from the Giant Monopole Resonance. Phys. Rev. Lett. 82 , 691-694 (1999). https://doi.org/10.1103/PhysRevLett.82.691
- [42] W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, Numerical recipes - the art of scientific computing (Cambridge University Press, Cambridge, 1986)
- [43] S. Typel, G. R¨ opke, T. Kl¨ ahn, D. Blaschke, H.H. Wolter, Composition and thermodynamics of nuclear matter with light clusters. Phys. Rev. C 81 , 015803 (2010). https://doi.org/10.1103/PhysRevC.81.015803. arXiv:0908.2344 [nucl-th]
- [44] J.C. Slater, A Simplification of the Hartree-Fock Method. Phys. Rev. 81 , 385 (1951). https://doi.org/https://doi.org/10.1103/PhysRev.81.385
- [45] J.A. Nolen, Jr., J.P. Schiffer, Coulomb energies. Ann. Rev. Nucl. Part. Sci. 19 , 471-526 (1969). https://doi.org/10.1146/annurev.ns.19.120169.002351
- [46] S. Shlomo, W.G. Love, Core Polarization and Coulomb Displacement Energies. Phys. Scr. 26 , 280 (1982). https://doi.org/10.1088/0031-8949/26/4/005
- [47] J. Geng, J.J. Li, W.H. Long, Y.F. Niu, S.Y. Chang, Pseudospin symmetry restoration and the in-medium balance between nuclear attractive and repulsive interactions. Phys. Rev. C 100 (5), 051301 (2019). https://doi.org/10.1103/PhysRevC.100.051301. arXiv:1907.00597 [nucl-th]

- [48] B. Wei, Q. Zhao, Z.H. Wang, J. Geng, B.Y. Sun, Y.F. Niu, W.H. Long, Novel relativistic mean field Lagrangian guided by pseudo-spin symmetry restoration. Chin. Phys. C 44 (7), 074107 (2020). https://doi.org/10.1088/1674-1137/44/7/074107
- [49] D. Adhikari, et al., Precision Determination of the Neutral Weak Form Factor of 48 Ca. Phys. Rev. Lett. 129 (4), 042501 (2022). https://doi.org/10.1103/PhysRevLett.129.042501. arXiv:2205.11593 [nucl-ex]
- [50] D. Adhikari, et al., Accurate Determination of the Neutron Skin Thickness of 208 Pb through Parity-Violation in Electron Scattering. Phys. Rev. Lett. 126 (17), 172502 (2021). https://doi.org/10.1103/PhysRevLett.126.172502. arXiv:2102.10767 [nucl-ex]
- [51] S. Shlomo, R. Schaeffer, The difference between neutron and proton radii in the Ca isotopes. Phys. Lett. B 83 , 5-8 (1979). https://doi.org/10.1016/0370-2693(79)90875-X
- [52] B.C. Clark, L.J. Kerr, S. Hama, Neutron densities from a global analysis of medium-energy proton nucleus elastic scattering. Phys. Rev. C 67 , 054605 (2003). https://doi.org/10.1103/PhysRevC.67.054605. arXiv:nucl-th/0209052
- [53] K. Chatziioannou, H.T. Cromartie, S. Gandolfi, I. Tews, D. Radice, A.W. Steiner, A.L. Watts, Neutron stars and the dense matter equation of state: from microscopic theory to macroscopic observations (2024). arXiv:2407.11153 [nucl-th]
- [54] I. Tews, R. Somasundaram, D. Lonardoni, H. G¨ ottling, R. Seutin, J. Carlson, S. Gandolfi, K. Hebeler, A. Schwenk, Neutron matter from local chiral EFT interactions at large cutoffs (2024). arXiv:2407.08979 [nucl-th]
- [55] W.G. Lynch, M.B. Tsang, Decoding the density dependence of the nuclear symmetry energy. Phys. Lett. B 830 , 137098 (2022). https://doi.org/10.1016/j.physletb.2022.137098. arXiv:2106.10119 [nucl-th]
- [56] C.Y. Tsang, M.B. Tsang, W.G. Lynch, R. Kumar, C.J. Horowitz, Determination of the equation of state from nuclear experiments and neutron star observations. Nature Astron. 8 (3), 328-336 (2024). https://doi.org/10.1038/s41550-023-02161-z. arXiv:2310.11588 [nucl-th]
- [57] M. Falanga, E. Bozzo, A. Lutovinov, J.M. Bonnet-Bidaud, Y. Fetisova, J. Puls, Ephemeris, orbital decay, and masses of ten eclipsing high-mass X-ray binaries. Astron. Astrophys. 577 , A130 (2015). https://doi.org/10.1051/0004-6361/201425191. arXiv:1502.07126 [astro-ph.HE]

- [58] J. Antoniadis, et al., A Massive Pulsar in a Compact Relativistic Binary. Science 340 , 6131 (2013). https://doi.org/10.1126/science.1233232. arXiv:1304.6875 [astro-ph.HE]
- [59] E. Fonseca, et al., The NANOGrav Nine-year Data Set: Mass and Geometric Measurements of Binary Millisecond Pulsars. Astrophys. J. 832 (2), 167 (2016). https://doi.org/10.3847/0004-637X/832/2/167. arXiv:1603.00545 [astro-ph.HE]
- [60] M.C. Miller, et al., PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter. Astrophys. J. Lett. 887 (1), L24 (2019). https://doi.org/10.3847/2041-8213/ab50c5. arXiv:1912.05705 [astro-ph.HE]
- [61] T.E. Riley, et al., A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation. Astrophys. J. Lett. 887 (1), L21 (2019). https://doi.org/10.3847/2041-8213/ab481c. arXiv:1912.05702 [astro-ph.HE]
- [62] M.C. Miller, et al., The Radius of PSR J0740+6620 from NICER and XMM-Newton Data. Astrophys. J. Lett. 918 (2), L28 (2021). https://doi.org/10.3847/2041-8213/ac089b. arXiv:2105.06979 [astro-ph.HE]
- [63] T.E. Riley, et al., A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy. Astrophys. J. Lett. 918 (2), L27 (2021). https://doi.org/10.3847/2041-8213/ac0a81. arXiv:2105.06980 [astro-ph.HE]
- [64] R.C. Tolman, Static solutions of Einstein's field equations for spheres of fluid. Phys. Rev. 55 , 364-373 (1939). https://doi.org/10.1103/PhysRev.55.364
- [65] J.R. Oppenheimer, G.M. Volkoff, On massive neutron cores. Phys. Rev. 55 , 374-381 (1939). https://doi.org/10.1103/PhysRev.55.374
- [66] J. Duflo, A.P. Zuker, Microscopic mass formulae. Phys. Rev. C 52 , R23 (1995). https://doi.org/10.1103/PhysRevC.52.R23. arXiv:nucl-th/9505011
- [67] M. Oertel, M. Hempel, T. Kl¨ ahn, S. Typel, Equations of state for supernovae and compact stars. Rev. Mod. Phys. 89 (1), 015007 (2017). https://doi.org/10.1103/RevModPhys.89.015007. arXiv:1610.03361 [astro-ph.HE]
- [68] S. Typel, Equations of state for astrophysical simulations from generalized relativistic density functionals. J. Phys. G 45 (11), 114001 (2018). https://doi.org/10.1088/1361-6471/aadea5

[69] J.R. Stone, V. Dexheimer, P.A.M. Guichon, A.W. Thomas, S. Typel, Equation of state of hot dense hyperonic matter in the Quark-Meson-Coupling (QMC-A) model. Mon. Not. Roy. Astron. Soc. 502 (3), 3476-3490 (2021). https://doi.org/10.1093/mnras/staa4006. arXiv:1906.11100 [nucl-th]
Stefan Typel Shalom Shlomo
Improving relativistic energy density functionals with tensor couplings
Energy density functionals (EDFs) have been used extensively with great
success to calculate properties of nuclei and to predict the equation of state
(EOS) of dense nuclear matter. Besides non-relativistic EDFs, mostly of the
Skyrme or Gogny type, relativistic EDFs of different types are in widespread
use. In these latter approaches, the effective in-medium interaction is
described by an exchange of mesons between nucleons. In most cases, only
minimal meson-nucleon couplings are considered. The effects of additional
tensor couplings were rarely investigated. In this work, a new relativistic EDF
with tensor couplings and density dependent minimal meson-nucleon couplings
will be presented. The parameters of the model are determined using a carefully
selected set of experimental data with realistic uncertainties that are
determined self-consistently. Predictions for various nuclear observables, the
nuclear matter equation of state, and properties of neutron stars are
discussed.