Nuclear Equation of State in the Relativistic Point-Coupling Model Constrained by Excitations in Finite Nuclei

: Nuclear equation of state is often described in the framework of energy density functional. However, the isovector channel in most functionals has been poorly constrained, mainly due to rather limited available experimental data to probe it. Only recently, the relativistic nuclear energy density functional with an effective point-coupling interaction was constrained by supplementing the ground-state properties of nuclei with the experimental data on dipole polarizability and isoscalar monopole resonance energy in 208 Pb, resulting in DD-PCX parameterization. In this work, we pursue a complementary approach by introducing a family of 8 relativistic point-coupling functionals that reproduce the same nuclear ground-state properties, including binding energies and charge radii, but in addition have a constrained value of symmetry energy at saturation density in the range J = 29, 30, ..., 36 MeV. In the next step, this family of functionals is employed in studies of excitation properties such as dipole polarizability and magnetic dipole transitions, and the respective experimental data are used to validate the optimal choice of functional as well as to assess reliable values of the symmetry energy and slope of the symmetry energy at saturation. functional. The results show that the approach based on construction of a family of point-coupling interactions spanning the range of values of J and its restriction by the α D data provides a similar constraint on the symmetry energy as that of the DD-PCX functional obtained by including α D value for 208 Pb directly in minimizing the χ 2 function to determine the model parameters. Both strategies result in successful optimization of the isovector channel of the relativistic EDF.


Introduction
The construction of a universal energy density functional (EDF) that can make precise predictions on the properties of finite nuclei and can provide a reasonable description of nuclear matter properties is one of the major challenges in nuclear physics, and it is also relevant for nuclear astrophysics applications. The equation of state (EoS) of nuclear matter represents an essential ingredient for understanding many astrophysical phenomena, such as supernovae and neutron stars, as well as for the description of various nuclear processes in a stellar environment (for more details, see reviews [1,2]). In order to provide improved and universal nuclear energy density functionals as well as to improve their predictive power, many studies have been performed, e.g., see References [3][4][5][6][7][8][9]. The EoS around saturation density is often constrained by the ground-state properties of finite nuclei, and the nuclear excitation properties provide valuable information to optimize the EDFs used in the description of the EoS. More details about constraining the EoS from ground and collective excited states are given in an extensive review in Reference [10].
The general strategy in optimization of the functionals is to constrain their parameters in χ 2 minimization using nuclear ground-state properties such as binding energies, charge radii, spin-orbit splittings, etc. The properties of symmetric nuclear matter around saturation are rather well-constrained by fitting the model parameters to the ground-state observables of nuclei. However, these observables are not able to constrain well the isovector channel of the functionals, which in turn leads to poorly constrained symmetry energy parameters of the EoS, and additional constraints on the pseudo-observables on the nuclear matter are introduced to cure this deficiency. Constraining the isovector channel of the EDFs with appropriate tools or observables is one of the principal aims of the nuclear physics community to obtain reliable information about the symmetry energy parameters around the saturation density as well as the properties of neutron-rich nuclei. In recent years, many works have been devoted to determining reliable values for the symmetry energy J, the slope of the symmetry energy L, and incompressibility of nuclear matter K 0 around the saturation density (see Reference [10] and the references therein) as well as finding possible observables to constrain their values [11][12][13][14][15]. These studies have shown that the symmetry energy parameter J is expected between 31 and 32 MeV while the slope of the symmetry energy L is obtained between 30 and 90 MeV, and it is subject to larger uncertainties than the J value. Moreover, the incompressibility parameter K 0 , which is known as the curvature of the nuclear matter EoS, is predicted to be between 210 and 250 MeV [10,15]. There are also studies indicating larger incompressibilities if the ratio of the surface and volume coefficients c = K sur f /K vol from leptodermous expansion in powers of A −1/3 of the incompressibility of finite nuclei K A is allowed to vary [16]. The decription of the EoS at higher densities represents a great challenge for ongoing research. Heavy-ion collisions can compress nuclear matter to densities present in neutron stars and core-collapse supernovae and, thus, can provide additional constraints of the EoS at higher densities [17].
Recent studies have shown that the isovector channel of the EDFs is strongly correlated with the neutron skin thickness, neutron star mass-radii, and isovector dipole excitation properties of nuclei and that these observables can be used to constrain the isovector channel of the functionals (see References [10,12,18] and the references therein). There are recent experimental constraints on the neutron skin thickness, obtained from the parity violating elastic electron scattering on 208 Pb (PREX) [19]; however, the experimental uncertainties despite recent improvements remain rather large [20]. Among relevant observables, the isovector dipole excitations of nuclei have been known for a long time and represent a very useful constraint for optimization of the functionals. The dipole polarizability of nuclei, which is calculated from the isovector dipole excitation spectra, can also be used as a reliable observable to constrain the isovector channel and symmetry energy parameters of the nuclear EoS around the saturation density [21]. Charge-exchange transitions in nuclei also provide useful constraints for the isovector channel and the symmetry energy of the EDFs, in particular isobaric analog states [22], Gamow-Teller resonances [23], and anti-analog giant dipole resonances [24,25].
Recently, a new parameterization for the relativistic point-coupling interaction, denoted as DD-PCX [26], was optimized directly in χ 2 minimization using the isovector dipole and isoscalar monopole excitation properties of 208 Pb along with the nuclear groundstate properties of a set of nuclei. To this aim, a unified framework of the relativistic Hartree-Bogoliubov (RHB) model and self-consistent relativistic quasiparticle random phase approximation (R(Q)RPA) was used in the optimization procedure and further analysis [26,27]. Using the collective excitation properties of 208 Pb, both the isoscalar and isovector channels of the EDF were constrained in a unique way, in which the symmetry energy and the slope of the symmetry energy around the saturation density were obtained as J = 31.12 ± 0.32 MeV and L = 46.32 ± 1.68 MeV, respectively. These values are also consistent with the model predictions from various approaches given in References [10,18]. Including the constrained isoscalar giant monopole resonance (ISGMR) energy in the optimization of the DD-PCX functional, the incompressibility of nuclear matter was also obtained as K 0 = 213.03 ± 3.54 MeV, which is also consistent with the predictions of previous studies [10,15].
The aim of this work is to introduce a complementary approach to constrain the symmetry energy parameters using relativistic EDFs with point-coupling interactions. First, a family of 8 functionals was established, which reproduced the same nuclear groundstate properties for a given set of nuclei, whereas the symmetry energy at the saturation density of the EoS was constrained to the values in the range of J = 29, 30, . . . , 36 MeV in the χ 2 minimization procedure. Furthermore, this family of functionals was used in studies of excitation properties such as dipole polarizability and magnetic dipole transitions, and by using respective experimental data, the optimal choice of the functional is discussed. In this way, insights into the values of the symmetry energy and slope of the symmetry energy at saturation are obtained. The results of the present study are also analyzed in comparison to the existing point-coupling functionals, DD-PC1 and DD-PCX, as well as with previous studies based on different EDFs [21,28]. Various families of the functionals have already been established both in the nonrelativistic and relativistic framework. [10,21,28]. Since the respective results are subject to considerable model dependence, it is important to provide several independent families of the EDFs. The isovector properties of the EDFs are also closely related to the size of the neutron skin thickness ∆R np . Thus, the relation of the excitation properties with the ∆R np values is also explored in the present study.
This paper is organized in the following way. In Section 2, the basic formalism of the effective relativistic point-coupling interaction is introduced. Section 3 includes the results of the present analysis of the neutron skin thickness as well as the electric and magnetic dipole transitions and their relation to the symmetry energy properties. Finally, in Section 4, the conclusions of the present study are given.

Relativistic Point-Coupling Functionals for Studies of the Symmetry Energy
Among the variety of EDFs currently used in the description of static and dynamic properties of finite nuclei along the nuclide map as well as the nuclear EoS [10], the relativistic EDF represents one of successful frameworks [29]. It has already been employed in the description of various nuclear excitation properties [30][31][32][33][34][35][36] and astrophysically relevant weak interaction processes [27,[37][38][39][40][41]. The theoretical framework used in this work is based on a self-consistent relativistic mean-field model with density dependent point-coupling interaction. It is derived from the effective Lagrangian density including the isoscalar-scalar, isoscalar-vector, and isovector-vector four-fermion interactions [4,29]: In addition to the free-nucleon terms, the Lagrangian density includes point-coupling interaction terms and coupling of protons to the electromagnetic field. The derivative term accounting for the leading effects of finite-range interactions is also necessary for a quantitative description of nuclear density distribution and radii. The couplings of the interaction terms a s , a v , and a TV are functionals of the vector density ρ v = j µ j µ with the nucleon four-current j µ =ψγ µ ψ [29]. The Dirac equation for the nucleons is obtained from the variation of the Lagrangian density (1) with respect toψ [29]: with the nucleon self-energies defined by the following: where Σ µ , Σ S , and Σ µ TV denote the isoscalar-vector, isoscalar-scalar, and isovector-vector self-energies, respectively. The isoscalar-vector self-energy Σ µ also includes the "rear-rangement" terms Σ µ R , which are essential for energy-momentum conservation and the thermodynamical consistency of the model [42]. The currents and densities are given in more details in Reference [29].
For the functional form of the couplings, the following practical ansatz is chosen [4,29]: with x = ρ/ρ sat and ρ sat denoting the nucleon density at saturation in symmetric nuclear matter. The point-coupling model has ten parameters: and δ S . Following Reference [29], the parameters c V , a TV , and c TV are set to zero. In order to describe the properties of open-shell nuclei, the relativistic Hartree-Bogoliubov (RHB) model [29] is used. The pairing field is formulated using separable pairing force, which also includes two parameters for the pairing strength (G p and G n ) [43].
In addition to the nuclear ground-state and excitation properties, the relativistic EDF can also be employed in the description of nuclear matter properties. In the translationally invariant infinite nuclear matter, the derivative term in Equation (1) drop out and the electromagnetic field is neglected [44]. More details about the implementation of the density dependent point-coupling interaction in nuclear matter description are given in Reference [44]. The binding energy per particle around the saturation density can be expanded as power series around δ → 0 [10], where the isospin asymmetry term is δ = (ρ n − ρ p )/(ρ n + ρ p ). The first term of Equation (8) is the EoS of symmetric nuclear matter, given by where K 0 represents the incompressibility of symmetric nuclear matter and is given by The second term of Equation (8) corresponds to the symmetry energy and can be expanded in terms of the density, where J ≡ S(ρ 0 ) corresponds to the symmetry energy at saturation density, L ≡ 3ρ 0

S(ρ)
∂ρ | ρ=ρ 0 is the slope of the symmetry energy at saturation, and K sym ≡ 9ρ 2 0 S(ρ) ∂ρ 2 | ρ=ρ 0 is the incompressibility of the symmetry energy at saturation. The symmetry energy properties will be discussed in the following section in relation to electric and magnetic excitations in finite nuclei.
In Reference [26] the relativistic point-coupling interaction DD-PCX has been constrained using the properties of isoscalar giant monopole resonance (ISGMR) energy and dipole polarizability of 208 Pb along with the nuclear ground-state properties: binding energies (34 nuclei), charge radii (26 nuclei), and mean pairing gaps (15 nuclei) of the selected open-shell nuclei. In addition to the ground-state properties of nuclei, the excitation data have been used as an additional constraint for the relevant properties of the nuclear matter, incompressibility parameter (K 0 ), symmetry energy (J), and slope of the symmetry energy (L) at saturation. In this work, a complementary approach is used, i.e., a set of 8 density-dependent point-coupling functionals is established by χ 2 minimization, using not only the given ground-state properties of nuclei but also an additional constraint on the symmetry energy at the saturation in the range of values, J = 29-36 MeV. In optimization of the interaction, the average neutron and proton pairing gap values are calculated using [45] where v and u are occupation factors of the single-particle states and ∆ k is the state dependent single-particle gap. On the other side, the pairing gaps were extracted from the experimental data on binding energies using a five-point formula [28]. In constraining the model parameters of the DD-PC family of functionals, the following experimental data on nuclear ground-state properties from light to heavy nuclei were used: binding energies (72 nuclei), charge radii (36 nuclei), diffraction radii (22 nuclei), surface thickness (18 nuclei), neutron pairing gaps (15 nuclei), and proton pairing gaps (14 nuclei) [28]. In Table 1, the selected nuclei are presented alongside the relevant data used in the fitting of the functionals. A similar choice of observables and the set of experimental data were adopted in adjusting the model parameters of Skyrme functionals in Reference [28]. In the present calculation, the adopted errors for the binding energies, charge radii, diffraction radii, surface thickness, and pairing gaps were taken as 1.0 MeV, 0.02 fm, 0.05 fm, 0.05 fm, and 0.15 MeV, respectively. In addition, nuclear matter properties were constrained with the given adopted errors: saturation density ρ sat = 0.152 fm −3 (10%), binding energy at saturation E/A = −16.14 MeV (10%), effective mass m * /m = 0.58 (10%), and incompressibility parameter K 0 = 230 MeV (1%), as used in Reference [4].   [46,47] and charge radii r c (36 nuclei) [48] are used alongside with the neutron (15 nuclei) and proton (14 nuclei) mean gap values ∆ n(p) , diffraction radius r di f f r (22 nuclei), and surface thickness σ (18 nuclei) values for the selected nuclei [28] in the fitting of the DD-PC (J = 29, 30, . . . , 36 MeV) functionals. In the second line, we display the adopted errors for each observable. The adopted error is also multiplied by a further integer weight factor that is given in the parentheses next to each observable. Although the same ground-state data were used systematically in the fitting protocol, the DD-PC functionals can accommodate a broad range of J values for the symmetry energy while providing similar quality in the description of binding energies and radii. Table 2 shows the resulting parameter sets of 8 functionals from the DD-PC family with J = 29, 30, . . . , 36 MeV. In the next step, this set of functionals will be examined by using nuclear excitation data and the possibility to constrain the optimal range of the symmetry energy parameters will be investigated.

Results
In this section, we present the results for the nuclear matter properties in relation to the electric and magnetic dipole excitation properties using the new family of pointcoupling functionals spanning the range of values for the symmetry energy at saturation J = 29, 30, . . . , 36 MeV. Thus, we can analyze the effects of varying the symmetry energy in a systematic way.

Nuclear Matter Properties and Equation of State
In Table 3, the nuclear matter properties at saturation density are listed for the family of point-coupling functionals introduced in this work, DD-PC-J29, . . . , DD-PC-J36, including energy per nucleon E/A; incompressibility modulus K 0 , spanning the range of values of the symmetry energy J = 29, 30, . . . , 36 MeV; and the slope of the symmetry energy at saturation density L. In addition, the respective values are also shown for the DD-PC1 [4] and DD-PCX [26] functionals. Except for the DD-PCX functional, the incompressibility parameters K 0 is limited to the same value K 0 = 230 MeV. As mentioned above, the DD-PCX functional was optimized using the isoscalar monopole resonance energy of 208 Pb to constrain the value of K 0 , which in turn leads to a lower incompressibility parameter around the saturation densities. Because of the strong correlation between the J and L values [49], it is also seen that the L values increase with increasing J values for the DD-PC family of interactions. In Figure 1, the EoS of symmetric nuclear matter (panels (a) and (c)) and pure neutron matter (panels (b) and (d)) are displayed for the introduced family of point-coupling functionals together with the ones for the DD-PC1 and DD-PCX interactions. Considering Figure 1a, the results are the same for the densities below 0.3 fm −3 for symmetric nuclear matter. Only when moving toward higher densities, differences are obtained in the values of the energy per nucleon. The divergence of the nuclear matter EoS at high-density regions, obtained when using different functionals, is a well-known problem in nuclear physics. Although some probes have been found to constrain this region [1,17,50], the model predictions are quite different and the high-density behavior of the EoS still remains unknown. It is also seen that the DD-PC functionals, which are obtained by varying the symmetry energy parameters, give the same results for the symmetric nuclear matter at each density. In other words, the EoS of symmetric nuclear matter is not sensitive to the changes only in the symmetry energy parameters of the DD-PC functionals. In Figure 1c, the results are presented up to twice the nuclear saturation density 0.34 fm −3 , and we also compare our results with the order-by-order chiral effective-field theory (χEFT) predictions, with NN and 3N interactions up to N 3 LO (see Reference [51] for more details). Considering the LO and NLO predictions for the symmetric nuclear matter, binding energy per nucleon E/A decreases smoothly with increasing density, i.e., their behavior is not compatible with all other model predictions. However, the behavior of the EoS with the higher-order predictions of the χEFT theory, i.e, N 2 LO and especially N 3 LO, appears in qualitative agreement with the DD-PC results.
The results for the energy per nucleon of pure neutron matter are presented in Figure 1b,d. In contrast to the symmetric nuclear matter, neutron matter appears rather sensitive to the choice of symmetry energy parameters for the DD-PC family of functionals. As it is known, the symmetry energy is stiffer (softer) for functionals with larger (lower) values of J and L. Increasing the value of the symmetry energy at saturation J, the energy per nucleon value increases at a faster rate at higher densities, as expected. In Figure 1d, we also compare our results for the pure neutron matter with those from the χEFT theory. The respective EoS of the DD-PC functionals are found to be in agreement with the N 2 LO and N 3 LO results. In particular, the DD-PCX functional appears in remarkable agreement with the N 3 LO results. Among the family of DD-PC functionals, the DD-PC-J31 fits best to the N 3 LO neutron matter EoS.  The dependence of the symmetry energy on the density is presented in Figure 2a, for the same family of point-coupling functionals with J = 29, 30, . . . , 36 MeV alongside the DD-PCX and DD-PC1 functionals. Below the saturation density, the symmetry energy values do not change around ρ 0 ≈ 0.1 fm −3 ; the results are the same for all the functionals. By increasing the density, the symmetry energy values also start to change in accordance with their different values for J and L. The symmetry energy as a function of density is higher (lower) for the functionals with larger (lower) values of the symmetry energy parameters, as expected [10]. In Figure 2b, we also display the symmetry energy of other models at lower densities to compare our results. Apart from the χEFT theory results (LO, NLO, N 2 LO, and N 3 LO), we also present the symmetry energy limits obtained using the Isobaric analog state (IAS) and neutron skin (∆R np ) constraints (see Reference [22]). It is known that the IAS can be used to constrain the isovector channel of the nuclear energy density functionals. In this work, we find that almost all models predict increasing symmetry energy with increasing density, except for N 2 LO, which decreases slightly after ρ > 0.25 fm −3 . We find that the predictions and behavior of the symmetry energy using the point-coupling DD-PCX and DD-PC-J31 functionals are in good agreement with the results of the NLO and N 3 LO. Using the χEFT theory, the symmetry energy and its slope was constrained to J = 31.7 ± 1.1 MeV and L = 59.8 ± 4.1 MeV around the nuclear saturation density [51]. Moreover, the symmetry energy was constrained between 0.04 ≤ ρ ≤ 0.13 fm −3 using the IAS + ∆R np information, and the J and L values were obtained between 30.2-33.7 MeV and 35-70 MeV around the saturation densities [22], respectively. Considering the predictions in the low density region and around the saturation densities, we find that point-coupling DD-PCX (J = 31.12 ± 0.32 MeV) and DD-PC-J31 functionals appear very compatible with the results from the χEFT theory and the IAS + ∆R np constraints on the symmetry energy. In the following section the family of DD-PC functionals will also be tested with the dipole polarizability data on finite nuclei.   [51]). The isobaric analog state (IAS) and neutron skin (∆R np ) constraints on the symmetry energy (see Reference [22] for more details) are also presented.

Dipole Polarizability and Neutron Skin
The dipole polarizability α D was recently established as a useful additional constraint for the isovector channel in EDFs [21]. It is defined as where E is the excitation energy and m −1 (E1) is the inverse energy weight sum of the isovector dipole transition strength. The strength function is given by where |Φ denotes the ground-state of nucleus andQ (IV−E1) is the isovector transition dipole operator [27]. The transition matrix elements can be calculated using the quasiparticle random phase approximation [27,52].
In numerous studies both relativistic and nonrelativistic EDFs have been tested using the experimental data on α D for 48 Ca, 68 Ni, 120 Sn, and 208 Pb (see [10] and the references therein). Another useful isovector quantity to be considered in probing the isovector sector of the EDFs is the neutron-skin thickness [49]. However, its data are still rather limited, often model dependent, or have large uncertainties such as in the case of parity violating electron scattering experiment (PREX) [19]. Therefore, for any novel functional, it is interesting to explore both the dipole polarizability and neutron-skin thickness, and when possible, to compare the results of model calculations with the experimental data.
We start our analysis with the neutron-skin thickness ∆R np for the point-coupling functionals DD-PC1 and DD-PCX. The neutron skin thickness of a nucleus is defined as ∆R np = r 2 n 1/2 − r 2 p 1/2 , and r 2 n 1/2 and r 2 p 1/2 represent the neutron and proton root-mean-square radii, respectively. In Figure 3, we display the available results from antiprotonic atom experiments [53,54] along with the DD-PC1 and DD-PCX results. We also display the experimental and theoretical results for the Sn chain in Figure 4. As can be seen in Figure 3, the experimental data have large error bars for most of the nuclei. Comparing the experimental results with the DD-PC1 and DD-PCX ones, it is seen that both functionals provide reasonable results and that the predictions are generally within the experimental error bars. It is known that the symmetry energy (J) and slope of the symmetry energy (L) at saturation are strongly correlated with the neutron skin thickness in finite nuclei [49]. Since the J and L values are lower for the DD-PCX functional than it is for DD-PC1, in the former case, smaller neutron skin radii are predicted for all nuclei, as expected.
As discussed above, the symmetry energy parameters are strongly correlated with the dipole polarizability as well as with the neutron skin thickness of nuclei. In Figure 5, we display the relationship between the dipole polarizability and neutron skin thickness results using the previously introduced family of point-coupling functionals obtained for 48 Ca, 68 Ni, 120 Sn and 208 Pb. An increase in the symmetry energy parameters leads to an increase in the neutron skin thickness of nuclei. As shown in figure, by increasing the J value of the functional, both the neutron thickness and dipole polarizability predictions increase linearly. The horizontal violet bands represent the measured dipole polarizability values of nuclei with error bars from References [55][56][57][58][59][60]. The vertical green bands denote the neutron skin thickness values obtained from model-averaged results (see Reference [60]). The overlap between the two bands defines an area denoted by red rectangle, which restricts the original family of DD-PC interactions to those that can simultaneously describe the dipole polarizability and neutron-skin thickness predictions. For the nuclei considered in the analysis, the results of the DD-PC functional with J = 31 MeV are in good agreement with all the predictions for α D and ∆R np . In addition, we can observe that the DD-PCX interaction for all nuclei fits into the restricted range of α D and ∆R np values. As can be seen from Figure 5, the DD-PC1 functional predicts higher values for dipole polarizability and neutron skin thickness compared to DD-PCX. In addition, the DD-PC1 results are outside of the ranges given by the experiment and model calculations.   [53,54] in which the proton distributions were obtained from electron scattering data for Sn nuclei [61] and from muonic atom data for others [62][63][64].  . Dipole polarizability as a function of neutron skin thickness for 48 Ca, 68 Ni, 120 Sn and 208 Pb using the relativistic point-coupling functionals. The experimental data for the dipole polarizability values are taken from References [55][56][57][58][59][60]. The recent experimental data for 120 Sn is taken from Reference [58] and displayed (α D = 8.08 ± 0.6 fm 3 ) with the violet shaded sphere. The band for the neutron radii represents the model-averaged results using many functionals (see Reference [60] for more information).
Since the DD-PCX is optimized using the α D data for 208 Pb, the results for this nucleus are found in the middle of the intersection point of the experimental data for dipole polarizability and model-averaged results for neutron skin thickness. Although we only used the dipole polarizability of 208 Pb in the optimization of the DD-PCX functional, it is clear that it can successfully predict the dipole polarizabilities and neutron-skin thicknesses of other nuclei. The predictions of the DD-PCX for α D and ∆R np also indicate that the isovector channel of the respective relativistic EDF is constrained successfully.
In Figure 6, the dipole polarizability results of the 208 Pb is presented along with the results of the 48 Ca, 68 Ni, and 120 Sn. We present our calculations using the DD-PC family of functionals varying the J value as well as the DD-PC1 and DD-PCX ones, and the theoretical results using the other functionals were taken from Reference [60]. The ranges of experimental data are also denoted with vertical and horizontal bands for the considered nuclei. The band overlap area, denoted by the red line, reduces a large set of nonrelativistic and relativistic functionals to only a few that successfully reproduce the experimental data on α D for all nuclei under consideration. The most reliable functionals that can predict the dipole polarizability of all nuclei are the DD-PC functional with J = 31 MeV and the DD-PCX functional. The results show that the approach based on construction of a family of point-coupling interactions spanning the range of values of J and its restriction by the α D data provides a similar constraint on the symmetry energy as that of the DD-PCX functional obtained by including α D value for 208 Pb directly in minimizing the χ 2 function to determine the model parameters. Both strategies result in successful optimization of the isovector channel of the relativistic EDF.  Figure 6. Comparison of the theoretical results for the dipole polarizability with the experimental data for 48 Ca-208 Pb, 68 Ni-208 Pb, and 120 Sn-208 Pb. The theoretical results for the KDE0-J, DD-ME, Skyrme, SAMI-J, NL3Λ, and FSU functionals are taken from Reference [60]. The experimental data with the error bars are shown as the shaded regions over the graph. The experimental data for 48 Ca (α D = 2.07 ± 0.22 fm 3 ), 68 Ni (α D = 3.88 ± 0.31 fm 3 ), and 208 Pb (α D = 19.6 ± 0.6 fm 3 ) are taken from References [55], [56,60], and [59,60], respectively. For 120 Sn, the experimental value from References [57,60] (α D = 8.59 ± 0.37 fm 3 ) is also shown using the shaded regions, and the recent data from Reference [58] (α D = 8.08 ± 0.6 fm 3 ) is displayed with a green shaded sphere.

Sensitivity of Magnetic-Dipole Excitation to Symmetry Energy
In the following, we investigate the sensitivity of magnetic dipole (M1) excitation properties on the symmetry energy parameters by using the same family of the relativistic point-coupling functionals as discussed in previous sections. Thus, the RRPA calculations for the isovector (IV) M1 excitations in 48 Ca and 208 Pb nuclei are performed by using the DD-PC functionals with J = 29-36 MeV. In References [36,70], the relativistic EDF with point-coupling interaction has been for the first time employed in the description of IV-M1 excitations by employing the RHB + R(Q)RPA model. The same framework and model parameters from Reference [36] were used in the present study, except that its DD-PC parameterization was replaced to account for different values of the symmetry energy at the saturation density, J ≡ S(ρ = ρ 0 ). In the case of 48 Ca and 208 Pb, our quasiparticle random phase approximation (QRPA) method reduces to the normal Random Phase Approximation (RPA).
Up to the one-body-RPA level, the IV-M1 operator, which depends on the spinŝ ν and orbital angular momentuml ν , is given aŝ ll ν (k) + g where µ N is the nuclear magneton,τ 3 = 1 (−1) for protons (neutrons), and g coefficients are given as g (IV) l = 0.5 and g (IV) s = 4.706 [71,72]. For simplicity, we neglect the quenching effect in g coefficients [73][74][75]. Then, the M1-excitation strength is computed as For demonstration purposes, instead of δ(E −hω i ), we smear the strength with a Cauchy-Lorenz profile, for which the full width at half maximum (FWHM) is 1.0 MeV.
In Figure 7, the transition strength distributions for M1 excitations are presented for the family of DD-PC interactions spanning the symmetry energy in the range J = 29-36 MeV. One can observe a systematic dependence of the M1 excitation energy on the J value of the effective interaction used in the calculation. Namely, the small (large) symmetry energy coincides with the high (low) excitation energy obtained for the M1 transition. However, when considering the corresponding transition strengths, no remarkable change is obtained among J = 29-36 results. Next, we consider the moments of M1 transitions, that is, where K = 0 and K = 1 indicate the non-energy-weighted and energy-weighted sum of the B (M1) transition strength, respectively. Similar to that in Section 3.2, here, we consider the relationship between the M1 excitation properties and the neutron-skin thickness for the family of DD-PC interactions. In Figure 8, the energy-weighted sum of B (M1) strength for 48 Ca is presented as a function of the neutron-skin thickness ∆R np , obtained by using the DD-PC functionals spanning the range of J values as denoted in the figure. In the lower panel, the dependence of the M1 transition centroid energy on ∆R np is shown. Figure 9 shows the same analysis but for 208 Pb. Both for the 48 Ca and 208 Pb nuclei, the m 1 (M1) value becomes reduced (enhanced) for the larger (smaller) neutron-skin thickness ∆R np or equivalently larger (smaller) J value. This is consistent with the sensitivity found in Figure 7. Thus, the M1 transitions may also be used as an additional constraint for the symmetry energy properties. Notice also that, for the 208 Pb nucleus, since the dB M1 /dE distribution includes two major peaks, its m 1 (M1) value does not simply indicate the major excitation energy, in contrast to the 48 Ca nucleus. For the non-energy-weighted sum of the M1 strength m 0 (M1), no significant sensitivity to J is obtained. Therefore, as shown in Figure 8, the centroid energy,Ē ≡ m 1 (M1)/m 0 (M1), displays the same tendency as the m 1 (M1) value. This constant m 0 (M1) value is in agreement with previous studies [70,76], where the non-energy-weighted sum rule does not change for the M1 mode in considering the case when the pairing correlation is negligible. We note that, in this work, we performed the first investigation of the relation between M1 excitation and symmetry energy parameters. On the other hand, for further applications, there remain several tasks. Especially in theoretical calculations, the isovectorpseudovector (IV-PV) interaction plays an important role in the description of M1-excitation energy [36,70]. Its coupling parameter, however, still has a finite ambiguity [36,70]. Because of this ambiguity, the present RPA calculations have not completely reproduced the experimental data on M1 centroid energies, that is,Ē exp. = 10.2 MeV for 48 Ca [77][78][79][80] andĒ exp. = 7.3 MeV for 208 Pb [81]. We are planning more systematic optimization of theoretical parameters to remedy this gap. Nevertheless, we emphasize that the results of the present analysis display the sensitivity of M1 excitations to symmetry energy at saturation density.

Conclusions
Constraining the isovector channel of the EDFs, which is of paramount importance for description of the symmetry energy of the nuclear EoS, represents a continuous challenge for the nuclear physics and astrophysics community. In this work, we investigated the possible constraints to symmetry energy by exploiting the properties of dipole electric and magnetic excitations in finite nuclei. For this purpose, a novel family of 8 relativistic densitydependent point-coupling interactions is established, which varies the symmetry energy at saturation in the range of values J = 29-36 MeV. We exploited this family of functionals in the analysis of possible isovector constraints, in particular, the dipole polarizability and magnetic dipole excitation energy in nuclei. It is shown that, in this way, the symmetry energy can be further restricted by using available experimental data on electric and magnetic excitation properties. Through this analysis, we also validated the isovector properties of the DD-PCX functional, which was recently established not only using the nuclear ground-state properties of nuclei but also using additional constraints from the experimental data on dipole polarizability and isoscalar monopole resonance energy in 208 Pb directly in χ 2 optimization. The DD-PCX functional, together with the family of DD-PC functionals introduced in this work that cover a range of values for the symmetry energy at the saturation density, represent useful tools for further studies necessary to constrain the properties of nuclear EoS. The new set of functionals also provides an interesting possibility to establish EoS for supernovae and compact stars [2] in order to explore the sensitivity of their properties on variations in the EoS in a systematic way. Available information on the neutron star mass-radius relationship can also provide useful constraints on the DD-PC functionals.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

EDF
Energy density functional RHB Relativistic Hartree-Bogoliubov QRPA Quasiparticle random phase approximation EoS Equation of state