Abstract
We review the experimental knowledge on the dipole polarizability (DP) of nuclei and its relation to the neutron skin thickness and properties of the neutron-rich matter equation of state (EOS). The discussion focuses on recent experiments using relativistic Coulomb excitation in inelastic proton scattering at extreme forward angles covering a mass range from 40Ca to 208Pb. Constraints on the neutron skins and the density dependence of the symmetry energy are derived from a systematic comparison to calculations based on density functional theory (DFT) and ab initio methods utilizing interactions derived from chiral effective field theory (EFT). The results consistently favor a soft EOS around or slightly below the saturation point. An outlook is provided on possible improvements in the precision achievable in stable nuclei and studies of exotic neutron-rich unstable nuclei with upcoming experimental facilities.
1 Introduction
The nuclear equation of state (EOS) describes the energy per nucleon of nuclear matter as a function of proton and neutron densities [1]. It governs the properties of nuclei and neutron stars [, ] as well as the dynamics of core-collapse supernovae [4] and neutron star mergers [5]. As an example, Figure 1A illustrates the bounds of the mass–radius dependence of neutron stars predicted by different EOS models. A systematic description of the EOS from nuclear densities to those in neutron stars is a central goal of current physics []. Despite a wealth of new data at high densities from observations on the properties of neutron stars and neutron star mergers [7] and information on the intermediate density regime from central heavy ion collisions [, 9], experimental constraints on the EOS around the saturation density of nuclear matter are still insufficient.
FIGURE 1
The nuclear matter EOS can be approximately written as a sum of the energy per nucleon of symmetric matter and an asymmetry termwhere the nucleon density and the asymmetry parameter are defined by the neutron and proton density as
The symmetry energy factor in Equation 1 can be expanded around the saturation density as
Here, is the slope parameter at density .
The first term in Equation 1 representing symmetric nuclear matter is fairly well constrained by the compressibility derived from systematic measurements of the isoscalar giant monopole resonance (ISGMR) in nuclei [1]. Figures 1B,C [, ] illustrate the variety of experimental and theoretical constraints on (also called in the literature) and defined in Equation 3. While these confine possible values of to a range of approximately 30–35 MeV, the uncertainties of are much larger.
As detailed below, all relevant theoretical models predict a strong correlation between and two experimentally accessible quantities, viz., the neutron skin thickness and the dipole polarizability. The connection is illustrated in Figure 2. The density distributions of neutrons and protons in the ground state can be determined from the condition of minimum energy. They approximately have the shape of Fermi distributions, as illustrated in Figure 2, for a nucleus with . The mean square radius of neutrons, , is slightly larger than that of protons . The difference between the two, , is defined as the neutron skin thickness.
FIGURE 2
The neutron skin thickness is sensitive to the value due to the following reason. As discussed above, the symmetry energy of nuclear matter at a given nucleon density depends on the square of the asymmetry parameter , defined in Equation 2. The first-order density dependence of the symmetry energy is represented by the slope parameter . Suppose that the density distributions in Figure 2 were determined for an value to have the minimum energy. There are density differences between the neutrons and protons in the inner part (higher nucleon density) and at the surface part (lower nucleon density). For a larger value, the density distributions change to have less density difference in the inner part, thereby reducing the symmetry energy in the higher density part. Consequently, the neutron skin thickness and the symmetry energy at the surface become larger for a conserved number of neutrons and protons.
The neutron skin thicknesses of medium-mass and heavy nuclei have been extracted from experiments studying elastic proton scattering [11], coherent production [12], antiprotonic atoms [, 14], and the isovector (IV) spin–dipole resonance []. Of particular importance are experiments using parity-violating polarized elastic electron scattering [16]. The parity-violating part of the reaction is mediated by the weak interaction and, due to the dominance of the neutron form factor, allows for extracting the neutron density distribution in an almost model-independentway []. Such experiments have been performed for 208Pb (lead radius experiment, or PREX) [1] and 48Ca (calcium radius experiment, or CREX) []. Neutron skins were determined by comparison with the well-known charge radii.
The dipole polarizability (DP) of nuclei can be obtained from measurements of the photoabsorption cross-sections. A connection between DP, neutron skin thickness, and parameters of the symmetry energy can only be made through models. Such calculations are presently based either on density functional theory (DFT) [] or ab initio coupled-cluster calculations [] using interactions derived from EFT []. Both types of models predict a strong correlation between the magnitudes of dipole polarizability, , and . The considerable experimental challenges of direct measurements of the neutron skin and the model dependencies of methods that extract the neutron skin from the difference of mass and charge radius [23] call for an alternative experimental observable. Because properties of the symmetry energy cannot be extracted directly from experiments but require theory input, measurement of the dipole polarizability provides independent constraints.
While several experimental techniques to measure the DP are discussed, the present review mainly focuses on recent progress using relativistic Coulomb excitation in forward-angle proton scattering at energies of several hundred MeV [24]. One advantage of this method is consistent results across the neutron separation energy, while many of the other experimental techniques are limited to either the energy region below or above. Even more important, measurements of the strength with relativistic Coulomb excitation can be extended to exotic nuclei at rare isotope beam facilities like RIKEN, FRIB, and the GSI Facility for Antiproton and Ion Research (FAIR). Such experiments are performed in inverse kinematics, where the virtual photon flux can be boosted by using a high- target and efficient setups with almost solid angle coverage for detection of neutron emission above [] and emission below the neutron threshold [26, 27]. In combination with the large cross sections, this will permit access to nuclei with extremely large neutron excess, much closer to the properties of neutron-rich matter relevant to the physics of neutron stars. In addition to the antiProton Unstable Matter Annihilation (PUMA) project [] aiming at the neutron skin thickness in unstable nuclei using antiproton annihilation, dipole polarizability measurements with relativistic Coulomb excitation are probably the only experimental probe promising insight into properties of the symmetry energy over a wide range of neutron-to-proton ratios.
The paper is organized as follows. Section 2 discusses how information on the neutron skin thickness and symmetry energy can be inferred from model calculations based on DFT (Section 2.1) and ab initio methods (Section 2.2). Section 3 is devoted to experimental issues. A short discussion of the available techniques in Section 3.1 is followed by a description of methods to disentangle electric and magnetic contributions to the DP in Section 3.2. The relevance of experimental information in the energy region of the isovector giant dipole resonance (IVGDR), as well as below the neutron threshold and above the IVGDR, is compared in Section 3.3, Section 3.4, and Section 3.5. The comparison of experimental and theoretical results (Section 4) for a range of nuclei from 40Ca to 208Pb and constraints on neutron skin thickness and the parameters of the symmetry energy extracted thereof are presented in Section 4.1 for DFT and Section 4.2 for ab initio approaches. Section 4.3 focuses on the difficulties of simultaneously describing the results of parity-violating elastic electron scattering and DP experiments with present-day models. Finally, Section 4.4 discusses the systematics of the DP and the role of volume and surface contributions to the symmetry energy. A summary and an outlook are given in Section 5.
2 Relation between dipole polarizability, neutron skin thickness, and symmetry energy
In this section, we discuss how information on the neutron skin thickness and parameters of the symmetry energy can be inferred from the comparison of the experimental dipole polarizability to theoretical predictions. At the moment, there are two classes of models, based on either DFT or an ab initio coupled-cluster approach. Because isovector observables are not well constrained in DFT, quantitative predictions of the DP can vary considerably. However, one can establish a robust correlation between the parameters and of the symmetry energy through . With ab initio-based models, one aims at an absolute prediction of , and the underlying symmetry energy parameters of the interaction can be used to calculate the EOS.
The dipole polarizability is related to the reduced transition strengths and the photoabsorption cross sections by
While the integral runs to infinity in principle, because of the inverse energy weighting a measurement of the strength up to excitation energies of about 60 MeV in light [] or 30 MeV in heavy nuclei [] is sufficient to achieve saturation. Thus, is dominated by the isovector giant dipole resonance (IVGDR), but contributions from the energy regions below and above are non-negligible, as discussed in Section 3.
2.1 Connections in density functional theory
An approximately linear correlation between and was demonstrated in Hartee–Fock calculations of 208Pb with relativistic [31] and Skyrme [] density functionals, as illustrated in Figure 3A. A comprehensive investigation of correlations between IV experimental observables and the bulk parameters of DFT models [33] demonstrates that these two quantities are also correlated with in heavy nuclei. Figure 3B shows, as an example, the correlations with the neutron form factor of 208Pb, which can be derived from a parity-violating elastic electron scattering experiment.
FIGURE 3
While this type of correlation is observed for all interactions, absolute values show large differences. In general, the magnitude of IV quantities like is not well constrained in DFT models because the model parameters are typically fitted to binding energies and charge radii of selected nuclei, which show little sensitivity to the IV parts of the nuclear interaction. A study of the relation between and the neutron skin in 208Pb for a large number of interactions illustrates the problem [34]. Figure 4A shows that the predictions for the neutron skin vary from 0.12 fm to 0.32 fm, and for a given value of , predictions for scatter wildly. However, the product of plotted versus (or ) shows a linear dependence with a high correlation coefficient [34], cf.Figure 4B. This relation can be understood within the droplet model [35] and provides a correlated range of values, as indicated for the case of 208Pb [36] in Figure 1C.
FIGURE 4
2.2 Connections in ab initio models
Ab initio calculations based on interactions derived from EFT play an important role in the attempt to systematically describe the EOS of neutron-rich matter at all densities []. Figure 5A displays examples of next-to-next-to-next-to-leading order predictions of the density behavior in the nuclear regime []. The upper and lower parts present the neutron and symmetric matter results, respectively, for two different families of interactions with somewhat different symmetry energy parameters shown in the left and right columns. The gray boxes indicate the value of the saturation density. The colored curves correspond to different cutoff parameters of the model space; for details, see [].
FIGURE 5
Predictions of and correlations with proton and neutron radii based on EFT interactions have been obtained from calculations based on a coupled-cluster expansion of the wave functions [] combined with the Lorentz-integral-transform approach to extract the strength []. An example of such calculations for 48Ca [] is presented in Figure 5B, where the correlation of with , and is displayed together with representative examples of DFT predictions. While the DFT results predict neutron skin values ranging from 0.16 fm to 0.22 fm, the ab initio results based on a set of interactions from [, ] consistently favor rather small values varying from 0.12 fm to 0.15 fm.
A major difference between the two theoretical approaches lies in the predicted relation between the proton and neutron radii. The DFT predictions of are approximately constant, most likely because the charge radius of 48Ca is in all cases part of the data set used to fix the model parameters. The ab initio calculations, on the other hand, predict a linear correlation, leading to the approximate constancy of the neutron skin. The absolute value of in the ab initio models shows a larger variation than the DFT calculations but can be well described by a linear correlation similar to . As discussed in the following, these correlations allow extracting constraints on the range of symmetry energy parameters based on the successful description of experimentally measured polarizabilities and charge radii. This type of calculation has been limited so far to closed-(sub)shell nuclei. For recent attempts of an extension to open-shell nuclei, see [, ].
3 Dipole polarizability from experiment
In this section, we discuss the experimental methods to extract the distribution in nuclei and the DP. It is technically difficult to directly measure the DP of nuclei as the response to a static electric field, although there exist exceptional cases of very light nuclei; see, for example, the works studying the deviation of elastic scattering cross section from Rutherford scattering [44, 45]. Instead, or distributions are measured and integrated to determine the DP by Equation 4. Some of the experimental methods discussed below are restricted in the accessible excitation energy range; that is, the techniques are applicable below or above the neutron emission threshold only. Thus, the role of contributions to the DP below from the IVGDR and from the energy region above the IVGDR is discussed in more detail. Both and transitions are excited, and possible ways of their distinction are briefly presented.
3.1 Experimental methods
3.1.1 Photoneutron measurement
The photoexcitation of nuclei above the neutron separation energy was intensively studied by using photoneutron measurements. The photoneutron reaction is conventionally written as , where stands for the number of emitted neutrons after photoexcitation. The cross-section is the sum of the , , …+ above the respective thresholds. From the 1960s to the 1980s, positron annihilation in flight was used for producing a quasi-monoenergetic -ray beam at Lawrence Livermore National Laboratory (LLNL) and at Saclay. Neutrons emitted after interaction with a target were thermalized and detected. For details, see [, ].
Neutron emission is the dominant decay process after photoexcitation for a nucleus as heavy as 208 because charged-particle decays are strongly suppressed by the Coulomb barrier and the decay branch is as low as 1%–2% []. Thus, the cross-sections in heavy nuclei can be compared with total photoabsorption cross-sections. The 208 cross-sections are plotted in Figure 6. The data were taken at LLNL (open light blue circles [] and solid blue circles []) and at Saclay (half-filled red circles [51]). The results from the two laboratories show clear discrepancies, which is also true for some other nuclei. Kawano et al. [] reported that “in general, the Saclay cross sections are larger than the Livermore data, whereas the Saclay cross-sections are smaller than the corresponding Livermore data.”
FIGURE 6
Later, quasi-monoenergetic photon beams produced by laser Compton backscattering (LCBS) became available at the National Institute of Advanced Industrial Science and Technology (AIST) [54], the High Intensity -ray Source (HIS) facility at the Triangle University National Laboratory [55] and the NewSUBARU facility [, ]. An electron beam in a storage ring is irradiated by laser photons to produce high-energy photons by head-on collisions []. The scattered photons are collimated to have a narrow energy distribution. The photon energy is variable either by changing the electron beam energy or the laser frequency. The 208 cross-section data measured at NewSUBARU [] are plotted as solid black circles in Figure 6.
3.1.2 Total photoabsorption
Total photon absorption was studied by applying transmission measurements. In this method, the attenuation of photons in a thick target was measured as a function of the photon energy for extraction of the photoabsorption cross sections. At the Mainz electron accelerator, a narrow photon beam was produced by the bremsstrahlung of an electron beam. The average photon flux was photons/MeV at 20 MeV. Two identical Compton spectrometers monitored the photon flux before and after a natural abundance target with a thickness of 40–200 cm []. The dominant atomic photoabsorption cross sections needed to be subtracted. A high-resolution transmission measurement at AIST was reported for 28 using a HPGe detector []. Recently, an experimental setup for photon transmission measurements has been in operation at the photon tagger NEPTUN [60] at the S-DALINAC accelerator in Darmstadt.
3.1.3 Compton scattering
Compton scattering from 208 was measured at Mainz using quasi-monoenergetic photons produced by positron annihilation in flight [61] up to a photon energy of 143 MeV. The flux of the photon beam was monitored with a Compton spectrometer. Elastically scattered photons were detected by large-volume NaI scintillation counters. The multipolarity-dependent cross sections were analyzed using the angular distributions. The imaginary part of the scattering cross sections at zero degrees is related to the total photon cross section.
3.1.4 Bremsstrahlung excitation functions
Photonuclear cross sections have been extracted from the radioactive decay of residual nuclei populated in particle emission after irradiation with thick-target bremsstrahlung. The excitation energy dependence can be determined by variation of the bremsstrahlung endpoint energy with an unfolding procedure [62]. However, this requires precise knowledge of the bremsstrahlung spectra, which is experimentally not available. While such spectra can be reliably calculated [63] with present-day Monte Carlo codes such as GEANT4 [], older versions contained poor approximations [65]. Results deduced from phenomenological approximations or using the analytical description of thin-target bremsstrahlung have potentially very large systematic uncertainties, typically not included in the quoted errors.
3.1.5 Relativistic Coulomb excitation
Relativistic Coulomb excitation is an important experimental tool to study the electric dipole response at radioactive ion beam (RIB) facilities. At beam energies of several hundred MeV/nucleon, cross sections are large and cover an excitation energy range including the IVGDR. The small number of beam particles can be compensated for neutron-rich nuclei by placing a neutron detector under because, at highly relativistic velocities, a small angular opening is sufficient to cover the full 4 solid angle range in the center-of-mass system. The method has been applied to study, for example, halo nuclei [66] and neutron-rich oxygen isotopes [67]. DP measurements in heavier nuclei have been performed for 68Ni [26, 27] and 130,132Sn [].
The method has also been developed to study stable nuclei using inelastic proton scattering under extreme forward angles, including . Such experiments require advanced methods to remove the background from the beam halo and atomic small-angle scattering in the target. Zero-degree setups have been realized at the Research Center for Nuclear Physics (RCNP), Osaka, Japan, for proton energies up to 400 MeV [68] and at the iThemba Laboratory for Accelerator-Based Science, Faure, South Africa, for 200 MeV [69]. An overview of experiments, data analysis, and physics problems addressed is provided in [24].
3.1.6 Nuclear resonance fluorescence
Nuclear resonance fluorescence (NRF) or experiments study the emission after resonant absorption of a photon. The reaction selectively excites states with large ground-state branching ratios. The cross-section contributions due to the decay to excited states can be estimated in spherical and vibrational nuclei from the population of the lowest excited states. The experiments can be performed with Ge detectors and thus offer a unique energy resolution. The measured quantities depend on the product of photoabsorption cross sections and ground-state branching ratios; thus, the method is limited to excitation energies below the neutron threshold because of the dominance of particle decay widths in the continuum. Experimental methods, physics, and applications are discussed in a recent review [70].
3.2 Decomposition of and contributions
A general problem of all experimental methods discussed above is the removal of magnetic contributions to the photoabsorption cross sections and the derived DP. Overall, contributions of strength to the DP are small except for very light nuclei []. However, they become relevant in the excitation energy region of the spinflip resonance []. Therefore, techniques to decompose and contributions are important.
No decomposition can be performed for the photoneutron and total photoabsorption experiments. In the excitation energy regime relevant to determining the DP, they can be distinguished in Compton scattering by combining measurements at forward and backward angles. The multipolarity can also be determined in NRF experiments using transversely polarized photons [73]. Measurements of the response relative to the polarization plane permit a unique assignment of the electric or magnetic character of the emitted radiation. A polarized beam can be extracted from off-axis bremsstrahlung or LCBS. The latter method is particularly efficient because the polarization of the laser light is fully transferred to the photon beam [70].
In relativistic Coulomb excitation, the virtual photon spectrum in the forward direction is dominated by . However, in the proton scattering experiments close to , one must consider contributions to the cross sections due to the nuclear excitation of the spinflip strength. Two independent methods have been applied to separate and cross section parts based either on the total spin transfer derived from the combined information of polarization-transfer observables or from a multipole decomposition analysis (MDA) of the cross section angular distributions [24]. Figure 7 presents some illustrative examples. The lower panel of Figure 7A displays the total spin transfer at for the nucleus 120Sn [] derived from measurements of the polarization-transfer observables and [74]. The upper panel presents the differential cross sections at and their decomposition in non-spinflip ( from Coulomb excitation) and spinflip ( from nuclear excitation) parts.
FIGURE 7
An example of the MDA analysis is presented in Figure 7B for 40Ca []. Spectra at different scattering angles are displayed in the upper panel, demonstrating strongly forward-peaked cross sections in the energy region of the IVGDR expected for Coulomb excitation. The lower panel shows the partial contributions to the cross sections at the most forward angle measured resulting from the MDA: (red), multipoles (orange), and nuclear background from quasifree scattering (green). Note that strength was neglected for 40Ca because it is concentrated in a single state []. A comparison of the two independent methods was made in studies of 96Mo [76], 120Sn [], and 208Pb [36], and good correspondence of the resulting and cross sections was found. The results shown in Figure 7A have also been confirmed in an MDA analysis []. Because the polarization-transfer measurements require secondary scattering, statistics are limited. Thus, in most experiments, the decomposition was restricted to MDA.
3.3 Contributions from the IVGDR
The largest contribution to the DP stems from the IVGDR, whose energy centroids lie well above . It is experimentally accessible with different techniques, and the comparison of results for the same nucleus provides an estimate of the typical accuracy of the DP. One can also average over results obtained with independent methods, thereby reducing the error bars. Some illustrative examples are presented in Figures 8A,B for 48Ca and 116Sn, respectively.
FIGURE 8
The strength distributions in 48Ca obtained from [] (blue circles) and [81] (green triangles) agree well. Results derived from the bremsstrahlung-induced activity of 47Ca [82] agree on the low-energy wing of the resonance but are significantly larger than the other data on the high-energy side. This can probably be traced back to the problems discussed in Sec. 3.1.4. The second example (B) compares photoabsorption cross sections for 116Sn from relativistic Coulomb excitation [] (blue squares) with data [, 84] (green left arrows and red right arrows). Reasonable agreement is observed in the maximum region of the IVGDR, but one finds significant differences on the low-energy flank. Such deviations are systematically observed in the stable Sn isotope chain, and for some isotopes also at the high-energy flank [].
In general, studies of the and reactions with LCB beams at NewSUBARU agree well with the results from the RCNP; see, for example, [85] for a study of Sn isotopes (black upward arrows in Figure 8B) or for 208Pb [] in Figure 6. However, a puzzling result reported for 209Bi is shown in Figure 8C. Although it differs from 208Pb by only one extra neutron, additional strength is seen on the high-energy side of the IVGDR, leading to a difference in not predicted by any model. This particular case certainly needs further investigation.
3.4 Contributions from the PDR
All particle-emission coincidence experiments accessing the strength are limited to the excitation region above the lowest particle separation threshold. Experimental evidence has accumulated that in nuclei with significant neutron excess strength—often concentrated in a resonance-like structure commonly termed pygmy dipole resonance (PDR)—can be found below [, 87]. Low-energy strength is also found in lighter nuclei with . Its contribution to the DP can be significant because of the inverse energy weighting, cf.Equation 4. As examples, they amount to approximately 10% in 58Ni [] and 8%–13% in the stable Sn isotopes [].
Most data on low-energy strength stem from experiments [70]. They suffer from the problem that branching ratios to excited states are typically unknown, and the extracted strength based on the g.s. transitions represents a lower limit only. Taking 120Sn as an example, the resulting strength distribution [88] reasonably agrees with a experiment [] measuring the total excitation strength up to approximately 6.5 MeV but totally underestimates the strength at higher excitation energies, cf.Figure 9A. Attempts have been made to model the inelastic contributions assuming statistical decay (see, e.g., [90]) but tend to overestimate contributions at low excitation energies. However, progress has been made recently by analyzing the cumulative population of the first state in even-even nuclei [70]. For the quoted example 120Sn, good agreement between the two experimental methods is achieved [88].
FIGURE 9
The origin of the low-energy strength in nuclei with neutron excess is a topic of current debate. It has been suggested to arise from an oscillation of the excess neutrons forming a skin against the (approximately) isospin-saturated core [, 93]. If true, its strength should be related to the neutron skin thickness and, in turn, to the parameters of the symmetry energy [–]. However, a recent study of the Sn isotope chain for mass numbers 111–124 casts doubts on such a picture [91]. The correlation between neutron excess and neutron skin thickness in Sn isotopes has been experimentally demonstrated with different methods [], but based on combined data from Oslo [98–100] and experiments [], only a minor fraction of the photoabsorption cross section (expressed as fraction of the Thomas-Reiche-Kuhn (TRK) sum rule) can be related to the PDR [101]. A decomposition into the tail of the IVGDR and two resonance-like structures is shown in Figure 9B [91]. The contribution interpreted as PDR is much smaller than those of the IVGDR and the prominent structure approximately 8 MeV. These findings rather point to an interpretation of the PDR as a low-energy part of a toroidal mode [102, 103]. At present, understanding the nature of the PDR remains an open problem. It is clear, however, that DFT predictions restricted to 1p-1h excitations cannot reliably estimate the low-energy strength distribution for cases where data are unavailable [91].
3.5 Contributions from high excitation energies
At excitation energies beyond the giant resonance region, photonuclear cross sections typically contribute a few percent only to the DP. However, for precision results, they must be considered. Data up to the pion threshold have been measured for a few cases, viz., natCa [], natSn [104], and 208Pb [61, 51]. They show approximately constant cross sections as a function of excitation energy and were considered for the extraction of the DP from experiments [, 36], neglecting an isotopic dependence. The dominant excitation mechanism in this energy regime is the quasi-deuteron effect [105]. It has been pointed out by Roca-Maza et al. [106] that these contributions are not included in model calculations based on DFT and should thus be removed compared to theoretical predictions. For heavy nuclei, they can be estimated using [], while in light nuclei, they are negligible in the energy range covered by the models [, , ].
The ratio of Coulomb excitation to quasifree cross sections in the experiments [24] drops with decreasing mass number limiting, in some cases, the excitation energy range accessible with an MDA for the extraction of cross sections. In such cases, model-dependent corrections must be applied. In the study of the Sn isotopic chain [], these were based on quasiparticle random phase approximation (QRPA) calculations folded with a Lorentzian to reproduce the experimentally measured width of the IVGDR. A particularly promising approach is discussed in [] for the example of 58Ni. An extension of the QRPA calculations to include quasiparticle vibration coupling has been successful in describing the width of the ISGMR and curing a longstanding discrepancy between the compressibility values extracted from 208Pb and lighter nuclei [108, 109]. The application to 58Ni demonstrates that the predicted high-energy tail of the IVGDR is largely independent of the chosen interaction []. This can be understood to result from the dominance of stochastic coupling [110]; that is, the strength distribution is mainly determined by the density of states and an average coupling matrix element between the 1p-1h and more complex states.
4 Extracting neutron skin thickness and symmetry energy properties from dipole polarizability data
In this section, we discuss constraints on the neutron skin thickness and symmetry energy properties derived from the comparison between model predictions and experimental studies of the DP. These refer to specific nuclei like 40Ca, 48Ca, and 208Pb but also systematic isotopic trends or a global mass dependence. The difficulties that presently available models have in simultaneously accounting for measured polarizabilities and asymmetries in parity-violating elastic electron scattering are illuminated.
4.1 Constraints based on density functional theory
The DPs of 40Ca and 48Ca have been studied in [, ], respectively. Figure 10A depicts their correlation and a comparison to selected DFT results. The four functionals are representative of widely used forms: non-relativistic Skyrme functionals SV [] and RD [] with different forms of density dependence, and relativistic functionals DD [113] with finite-range meson-exchange coupling and PC [114] with point coupling. All four have been calibrated to the same set of ground-state data to determine the model parameters.
FIGURE 10
The predictions are displayed as filled ellipses that represent the error as defined in [116]. The DD functional performs rather well. The other models tend to slightly overestimate the experimental mean values of both 40Ca and 48Ca, but their error ellipses do overlap with the experimental bands, except for PC. In all cases, the values for both nuclei are highly correlated. The dashed ellipses show the effect of additionally including the experimental value of 208Pb [36] in the fit, yielding functionals denoted “-alpha.” This improves the agreement with the experiment and shrinks the error ellipsoids. The models incorporate a span of symmetry energy parameters MeV and MeV for the calculations excluding (including) the 208Pb data point.
The strength distribution of the unstable neutron-rich nucleus 68Ni determined in an experiment measuring Coulomb excitation in inverse kinematics [26] is displayed in Figure 10B. The DP was extracted from a comparison to the model of [115]. The model results show a sensitivity to the assumed neutron skin thickness, as illustrated by the colored curves. A value of 0.17 (2) fm was extracted for the neutron skin thickness from the correlation between the two quantities.
A study of the DP in a long isotopic chain is particularly suited to investigate the connection with the neutron skin thickness. This can be best done in the Sn isotopes with neutron numbers between 50 and 82, where the proton shell closure stabilizes the g.s. deformation. There are many stable isotopes, and a study of the systematics of the DP was presented in []. The results are summarized in Figure 10C, which shows the evolution of between mass numbers 112 and 124. All DFT calculations predict an approximately linear increase as a function of neutron excess with roughly the same slope. The experimental results indicate a saturation between mass numbers 120 and 124, but the uncertainties (blue band) do not exclude a mass dependence similar to the theoretical results. The rightmost part of Figure 10C shows the predictions of the different models for the 208Pb DP after subtraction of the quasi-deuteron part (see the next paragraph). The models closest in absolute magnitude to the data tend to underpredict (208Pb), while those reproducing it overshoot the absolute values in the Sn chain, indicating that the functionals cannot yet fully describe the mass dependence of the DP. We note that strength distributions have also been measured for the unstable neutron-rich isotopes 130,132Sn [] but the extracted values of cannot be compared directly to the results of [] because the experiment only provided data above neutron threshold.
Roca-Maza et al. [106] combined the experimental DP data for 68Ni [26], 120Sn [], and 208Pb [36] to test a large variety of density functionals. Because the DFT calculations do not include contributions from the quasi-deuteron process dominating the photoabsorption cross sections above the energy region of the IVGDR, these had to be removed for a comparison [106]. Figure 10D presents correlation plots between the experimental results and theoretical predictions from a wide range of DFT interactions. Only a handful (marked in red) are capable of simultaneously describing all three data points. Based on this reduced set, systematic predictions of for other masses, and the symmetry energy parameters could be derived. The experimental results for 40,48Ca discussed above are fairly well described by these predictions.
4.2 Constraints based on ab initio models
An experimental study of the DP in 48Ca [] is of particular interest because it is accessible for both DFT and ab initio calculations, and a measurement of the neutron skin with parity-violating electron scattering is available []. The comparison is summarized in Figure 11A, where the blue band describes the experimental uncertainty. Ab initio results for the set of interactions from [, ] are displayed as green triangles, and a prediction from [] based on a normalization to the 48Ca charge radius is displayed as a green bar. Results from the set of density functionals described in [] are shown as red squares with some representative error bars, and the prediction from the analysis of [106] discussed above is shown as a black bar.
FIGURE 11
The DFT results tend to be somewhat high compared to the experiment. The ab initio results show a significant dependence on the chosen interaction, but it can be well approximated by a linear dependence. In principle, this allows for the derivation of boundaries on the neutron skin thickness and the symmetry energy. However, while the ab initio results shown were truncated in the coupled-cluster expansion at the 2p-2h level, subsequent work [118] demonstrated that inclusion of 3p-3h correlations lowers the values by %. The refined results in 48Ca are plotted in Figure 11C against corresponding calculations for 40Ca []. A high correlation similar to the DFT results shown in Figure 10A is observed. The purple uncertainty band from the ab initio results overlaps with the crossing of the experimental error bands. In particular, the NNLOsat interaction [] accurately describing binding energies and radii of nuclei up to 40Ca and the saturation point of symmetric nuclear matter now reproduces both DP values. The importance of including 3p-3h correlations has also been demonstrated in a recent measurement of the 68Ni charge radius []. Figure 11B illustrates the improvement in reproducing the correlation between the charge radius and [26] when going from the 2p-2h level (light blue band) to the inclusion of 3p-3h correlations (dark blue band).
As noted in Section 2.2, independent of the chosen interaction, a neutron skin thickness of approximately 0.14 fm is predicted for 48Ca, consistent with the value deduced from the measurement of the weak form factor []. The simultaneous description of the data in 40,48Ca and 68Ni implies that the underlying symmetry energy parameters are correct. A conservative estimate is provided by taking the full range of values from the set of ab initio interactions, viz., MeV and MeV.
Recent work has, for the first time, been able to extend the range of ab initio DP calculations based on EFT interactions to 208Pb []. A different technique was used to construct the interactions by history matching [120] using selected experimental observables in light nuclei. Moreover, low-energy nucleon–nucleon scattering phase shifts were additionally considered. The latter are responsible for tight constraints to rather small values of the resulting neutron skin thickness ( fm for 208Pb). The variation of the density dependence of the symmetry energy in these calculations is MeV.
4.3 Tension between polarizability and parity-violating elastic electron scattering in 208Pb
While in 48Ca there is fair agreement between the neutron skin thickness and symmetry energy properties derived from the different experiments, the parity-violating elastic electron scattering experiment on 208Pb [] finds a much larger neutron skin fm than most other work. Accordingly, an extraction of symmetry energy parameters based on the correlations established in DFT (see Section 2.1) leads to large symmetry energy values of MeV and MeV in contradiction to limits derived from astrophysical observations of neutron star radii and masses as well as the tidal deformability of neutron star mergers [7]. All astrophysical constraints point toward a softer EOS. This has led to speculations about a phase transition at intermediate densities [].
Because of the strong correlation between and for a given nucleus and values of different nuclei in DFT models, Reinhard et al. [116, 122] investigated whether it is possible to construct a DFT interaction capable of simultaneously describing the data for 48Ca and 208Pb. The analysis was based on representative families of non-relativistic and relativistic functionals. The isovector properties of EDFs are typically not well constrained by the input data used to fit the model parameters. As illustrated in Figure 12A for the case of 208Pb, it is possible to vary the symmetry energy parameters—and thereby the predicted and —over a fairly large range maintaining comparable description of ground-state properties [116]. Figure 12B [122] demonstrates that the polarizabilities and the neutron skin thickness of 48Ca could be consistently described, but it was impossible to construct an EDF simultaneously accounting for the neutron skin thickness of 208Pb extracted from the PREX experiment []. Similar conclusions were drawn in [123, 124]. In another recent attempt [125], a DFT interaction reasonably accounting for the measured parity-violating asymmetries in both the PREX and CREX experiments was constructed, but at the expense of unusual properties of the symmetry energy curvature and a very strong isovector coupling leading to density fluctuations in the nuclear interior.
FIGURE 12
4.4 Volume and surface contributions to the symmetry energy
Another way of extracting properties of the symmetry energy is a study of the mass dependence of the DP. A simple power law based on a model of two interpenetrating fluids has been given by Migdal, where denotes the second inverse moment of the photoabsorption cross sections and in units of mb/MeV ([126] and Refs. therein). A proportionality constant has been determined by Orce [127] from a fit to data over a wide mass range. Figure 13 [128] shows a comparison with a combined data set of measurements in light nuclei [] with the then-available (2016) data from relativistic Coulomb excitation for heavier nuclei as a green short-dashed line. Note that results for from [] are neglected because the hydrodynamical picture is highly questionable and corrections due to the magnetic polarizability are large [] for these very light nuclei. Results are severely underestimated in lighter nuclei where charged-particle decay dominates. The mass dependence is reasonably described for larger masses, but the proportionality coefficient of [127] is too low because additional contributions from the strength below the neutron threshold, as discussed in Sec. 3.3, must be considered.
FIGURE 13
For masses , surface contributions must be considered, modifying the volume term of the symmetry energy dominating for heavy nuclei. These can be parameterized as [128]
Here , and and denote the surface and volume coefficients of the symmetry energy, respectively. The numerical coefficient in Equation 5 is obtained from Migdal’s approach. A fit with parameters from binding energies of isobaric nuclei [130] shown in Figure 13 as a long-dashed blue line still underestimates the lower-mass data. Parameters of the study of [129] provide a better description (dotted red line). Results of a free fit of Equation 5 crucially depend on the inclusion (dotted-dashed black line) or exclusion (solid black line) of the 12C data point. The latter provides a better fit with MeV, [128] close to [129]. can be interpreted as , but measured at about of the saturation density [34, ].
5 Conclusion and outlook
We present a review of methods to measure the isovector response in nuclei and the extraction of the dipole polarizability from these data. The discussion focuses on recent results obtained with inelastic proton scattering under extreme forward angles at RCNP. At energies of a few hundred MeV, relativistic Coulomb excitation dominates the cross sections in these kinematics. The method combines certain advantages compared to other experimental techniques: 1) it measures the absorption and is thus independent of the knowledge of branching ratios; 2) a separation of and contributions to the cross sections can be achieved with different independent approaches; 3) the relevant excitation energy region from well below the neutron threshold across the IVGDR can be covered in a single experiment.
Constraints on the neutron skin thickness of nuclei and the parameters of the symmetry energy can be extracted from the strong correlations between these three quantities seen in all microscopic models. Results from nuclei covering a mass range between 40Ca and 208Pb consistently favor small neutron skins and a soft density dependence of the EOS around saturation density. In 208Pb serving as a benchmark for theory, this finding is at variance with the PREX results, while a similar study of 48Ca by the CREX collaboration conforms. The PREX result, hard to interpret in the framework of present theory, has led to an initiative (called Mainz radius experiment, or MREX) for a study with improved statistical and systematic errors at the new high-current Mainz energy-recovering superconducting accelerator (MESA) [132].
While the mass dependence of the DP is reasonably well-covered by the available data, future work should explore other degrees of freedom, such as the variation of neutron excess along isotopic chains and the role of deformation. The experimental uncertainties of the DP for key nuclei can be improved by the availability of independent measurements, as illustrated in Figure 6. New high-brilliance LCBS photon beam facilities are under construction at the Extreme Light Infrastructure–Nuclear Physics (ELI-NP) in Bucharest [, 134] and the Shanghai Laser Electron Gamma Source (SLEGS) at the Shanghai Synchrotron Radiation Facility [135]. Combined with advanced techniques for neutron detection [], these facilities promise a new quality of precision for experiments.
Major steps can be expected in the future at radioactive ion beam facilities, providing access to cases with much larger neutron excess than achievable for stable nuclei. Experimental tools for measuring relativistic Coulomb excitation in reverse kinematics are available, and pioneering studies of the dipole response in unstable nuclei have been performed at GSI [, 26, 27]. First results for the neutron-rich isotope 52Ca investigated at RIKEN have been reported [137]. Because of the high energy/nucleon availability, future experiments at FAIR are particularly promising for research on the dipole polarizability of exotic neutron-rich nuclei [].
Statements
Author contributions
PN-C: writing – original draft and writing – review and editing. AT: writing – original draft and writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Contract No. SFB 1245 (Project ID No. 79384907), by the Research Council of Norway through its grant to the Norwegian Nuclear Research Centre (Project No. 341985), by the JSPS KAKENHI Grant Number 25H00641, and by the Japan-South Africa Bilateral Funding Grant Number JPJSBP 120246502.
Acknowledgments
PvNC thanks the nuclear physics group at the University of Oslo for their kind hospitality during a stay where major parts of this work were done.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
Roca-MazaXPaarN. Nuclear equation of state from ground and collective excited state properties of nuclei. Prog Part Nucl Phys (2018) 101:96–176. 10.1016/j.ppnp.2018.04.001
2.
HuthSPangPTHTewsIDietrichTLe FèvreASchwenkAet alConstraining neutron-star matter with microscopic and macroscopic collisions. Nature (2022) 606:276–80. 10.1038/s41586-022-04750-w
3.
LattimerJ. Neutron stars and the nuclear matter equation of state. Annu Rev Nucl Part Sci (2021) 71:433–64. 10.1146/annurev-nucl-102419-124827
4.
YasinHSchäferSArconesASchwenkA. Equation of state effects in core-collapse supernovae. Phys Rev Lett (2020) 124:092701. 10.1103/PhysRevLett.124.092701
5.
RaaijmakersGGreifSKHebelerKHindererTNissankeSSchwenkAet alConstraints on the dense matter equation of state and neutron star properties from NICER’s mass–radius estimate of PSR J0740+6620 and multimessenger observations. Astrophys J Lett (2021) 918:L29. 10.3847/2041-8213/ac089a
6.
KoehnHRoseHPangPTHSomasundaramRReedBTTewsIet alFrom existing and new nuclear and astrophysical constraints to stringent limits on the equation of state of neutron-rich dense matter. Phys Rev X (2025) 15:021014. 10.1103/PhysRevX.15.021014
7.
LattimerJ. Constraints on nuclear symmetry energy parameters. Particles (2023) 6:30–56. 10.3390/particles6010003
8.
EsteeJLynchWGTsangCYBarneyJJhangGTsangMBet alProbing the symmetry energy with the spectral pion ratio. Phys Rev Lett (2021) 126:162701. 10.1103/PhysRevLett.126.162701
9.
TsangMBStoneJRCameraFDanielewiczPGandolfiSHebelerKet alConstraints on the symmetry energy and neutron skins from experiments and theory. Phys Rev C (2012) 86:015803. 10.1103/PhysRevC.86.015803
10.
DrischlerCFurnstahlRJMelendezJAPhillipsDR. How well do we know the neutron-matter equation of state at the densities inside neutron stars? A Bayesian approach with correlated uncertainties. Phys Rev Lett (2020) 125:202702. 10.1103/PhysRevLett.125.202702
11.
ZenihiroJSakaguchiHMurakamiTYosoiMYasudaYTerashimaSet alNeutron density distributions of 204,206,208Pb deduced via proton elastic scattering at Ep = 295 MeV. Phys Rev C (2010) 82:044611. 10.1103/physrevc.82.044611
12.
TarbertCMWattsDPGlazierDIAguarPAhrensJAnnandJRMet alNeutron skin of 208Pb from coherent pion photoproduction. Phys Rev Lett (2014) 112:242502. 10.1103/PhysRevLett.112.242502
13.
KłosBTrzcińskaAJastrzebskiJCzosnykaTKisielińskiMLubińskiPet alNeutron density distributions from antiprotonic 208Pb and 209Bi atoms. Phys Rev C (2007) 76:014311. 10.1103/PhysRevC.76.014311
14.
TrzcińskaAJastrzȩbskiJLubińskiPHartmannFJSchmidtRvon EgidyTet alNeutron density distributions deduced from antiprotonic atoms. Phys Rev Lett (2001) 87:082501. 10.1103/PhysRevLett.87.082501
15.
KrasznahorkayAFujiwaraMvan AarlePAkimuneHDaitoIFujimuraHet alExcitation of isovector spin-dipole resonances and neutron skin of nuclei. Phys Rev Lett (1999) 82:3216–9. 10.1103/PhysRevLett.82.3216
16.
MammeiJMHorowitzCJPiekarewiczJReedBTSfientiC. Neutron skins: weak elastic scattering and neutron stars. Annu Rev Nucl Part Sci (2024) 74:321–42. 10.1146/annurev-nucl-102122-024207
17.
HorowitzCJPollockSJSouderPAMichaelsR. Parity violating measurements of neutron densities. Phys Rev C (2001) 63:025501. 10.1103/PhysRevC.63.025501
18.
AdhikariDAlbatainehHAndroicDAniolKArmstrongDSAverettTet alAccurate determination of the neutron skin thickness of 208Pb through parity-violation in electron scattering. Phys Rev Lett (2021) 126:172502. 10.1103/PhysRevLett.126.172502
19.
AdhikariDAlbatainehHAndroicDAniolKAArmstrongDSAverettTet alPrecision determination of the neutral weak form factor of 48Ca. Phys Rev Lett (2022) 129:042501. 10.1103/PhysRevLett.129.042501
20.
BenderMHeenenP-HReinhardP-G. Self-consistent mean-field models for nuclear structure. Rev Mod Phys (2003) 75:121–80. 10.1103/RevModPhys.75.121
21.
HagenGPapenbrockTHjorth-JensenMDeanD. Coupled-cluster computations of atomic nuclei. Rep Prog Phys (2014) 77:096302. 10.1088/0034-4885/77/9/096302
22.
EpelbaumEHammerH-WMeißnerU-G. Modern theory of nuclear forces. Rev Mod Phys (2009) 81:1773–825. 10.1103/RevModPhys.81.1773
23.
ThielMSfientiCPiekarewiczJHorowitzCJVanderhaeghenM. Neutron skins of atomic nuclei: per aspera ad astra. J Phys G (2019) 46:093003. 10.1088/1361-6471/ab2c6d
24.
von Neumann-CoselPTamiiA (2019). Electric and magnetic dipole modes in high-resolution inelastic proton scattering at 0°. Eur Phys J A55, 110. 10.1140/epja/i2019-12781-7
25.
AdrichPKlimkiewiczAFallotMBoretzkyKAumannTCortina-GilDet alEvidence for pygmy and giant dipole resonances in 130Sn and 132Sn. Phys Rev Lett (2005) 95:132501. 10.1103/PhysRevLett.95.132501
26.
RossiDMAdrichPAksouhFAlvarez-PolHAumannTBenlliureJet alMeasurement of the dipole polarizability of the unstable neutron-rich nucleus 68Ni. Phys Rev Lett (2013) 111:242503. 10.1103/PhysRevLett.111.242503
27.
WielandOBraccoACameraFBenzoniGBlasiNBrambillaSet alSearch for the pygmy dipole resonance in 68Ni at 600 MeV/nucleon. Phys Rev Lett (2009) 102:092502. 10.1103/PhysRevLett.102.092502
28.
AumannTBartmannWBoine-FrankenheimOBouvardABrocheAButinFet alPUMA, antiproton unstable matter annihilation. Eur Phys J A (2022) 58:88. 10.1140/epja/s10050-022-00713-x
29.
FearickRWvon Neumann-CoselPBaccaSBirkhanJBonaitiFBrandhermIet alElectric dipole polarizability of 40Ca. Phys Rev Res (2023) 5:L022044. 10.1103/PhysRevResearch.5.L022044
30.
HashimotoTKrumbholzAMReinhardP-GTamiiAvon Neumann-CoselPAdachiTet alDipole polarizability of 120Sn and nuclear energy density functionals. Phys Rev C (2015) 92:031305. 10.1103/PhysRevC.92.031305
31.
TypelSBrownBA. Neutron radii and the neutron equation of state in relativistic models. Phys Rev C (2001) 64:027302. 10.1103/PhysRevC.64.027302
32.
BrownBA. Neutron radii in nuclei and the neutron equation of state. Phys Rev Lett (2000) 85:5296–9. 10.1103/PhysRevLett.85.5296
33.
ReinhardP-GNazarewiczW. Information content of a new observable: the case of the nuclear neutron skin. Phys Rev C (2010) 81:051303. 10.1103/PhysRevC.81.051303
34.
Roca-MazaXBrennaMColòGCentellesMViñasXAgrawalBKet alElectric dipole polarizability in 208Pb: insights from the droplet model. Phys Rev C (2013) 88:024316. 10.1103/PhysRevC.88.024316
35.
MyersWSwiateckiW. The nuclear droplet model for arbitrary shapes. Ann Phys (N.Y.) (1974) 84:186–210. 10.1016/0003-4916(74)90299-1
36.
TamiiAPoltoratskaIvon Neumann-CoselPFujitaYAdachiTBertulaniCAet alComplete electric dipole response and the neutron skin in 208Pb. Phys Rev Lett (2011) 107:062502. 10.1103/PhysRevLett.107.062502
37.
DrischlerCHebelerKSchwenkA. Chiral interactions up to next-to-next-to-next-to-leading order and nuclear saturation. Phys Rev Lett (2019) 122:042501. 10.1103/PhysRevLett.122.042501
38.
EkströmAJansenGRWendtKAHagenGPapenbrockTCarlssonBDet alAccurate nuclear radii and binding energies from a chiral interaction. Phys Rev C (2015) 91:051301. 10.1103/PhysRevC.91.051301
39.
HebelerKBognerSKFurnstahlRJNoggaASchwenkA. Improved nuclear matter calculations from chiral low-momentum interactions. Phys Rev C (2011) 83(R):031301. 10.1103/PhysRevC.83.031301
40.
HagenGEkströmAForssénCJansenGRNazarewiczWPapenbrockTet alNeutron and weak-charge distributions of the 48Ca nucleus. Nat Phys (2016) 12:186. 10.1038/nphys3529
41.
BaccaSPastoreS. Electromagnetic reactions on light nuclei. J Phys G (2014) 41:123002. 10.1088/0954-3899/41/12/123002
42.
BonaitiFBaccaSHagenGJansenGR. Electromagnetic observables of open-shell nuclei from coupled-cluster theory. Phys Rev C (2024) 110:044306. 10.1103/PhysRevC.110.044306
43.
BrandhermIBonaitiFvon Neumann-CoselPBaccaSColòGJansenGRet alElectric dipole polarizability of 58Ni. Phys Rev C (2025) 111:024312. 10.1103/PhysRevC.111.024312
44.
LynchWGTsangMBBhangHCCramerJGPuighRJ. Relativity, nuclear polarizability, and screening in sub-Coulomb elastic scattering. Phys Rev Lett (1982) 48:979–82. 10.1103/PhysRevLett.48.979
45.
RodningNLKnutsonLDLynchWGTsangMB. Measurement of the electric polarizability of the deuteron. Phys Rev Lett (1982) 49:909–12. 10.1103/PhysRevLett.49.909
46.
BermanBLFultzSC. Measurements of the giant dipole resonance with monoenergetic photons. Rev Mod Phys (1975) 47:713–61. 10.1103/RevModPhys.47.713
47.
DietrichSSBermanBL. Atlas of photoneutron cross sections obtained with monoenergetic photons. Data Nucl Data Tables (1988) 38:199–338. 10.1016/0092-640X(88)90033-2
48.
BeeneJRBertrandFEHorenDJAubleRLBurksBLGomez del CampoJet alHeavy-ion Coulomb excitation and photon decay of the giant dipole resonance in 208Pb. Phys Rev C (1990) 41:920–32. 10.1103/PhysRevC.41.920
49.
HarveyRRCaldwellJTBramblettRLFultzSC. Photoneutron cross sections of Pb206, Pb207, Pb208, and Bi209. Phys Rev (1964) 136:B126–31. 10.1103/PhysRev.136.B126
50.
BermanBLPywellREDietrichSSThompsonMNMcNeillKGJuryJW. Absolute photoneutron cross sections for Zr, I, Pr, Au, and Pb. Phys Rev C (1987) 36:1286–92. 10.1103/PhysRevC.36.1286
51.
VeyssiereABeilHBergèreRCarlosPLeprêtreA. Photoneutron cross sections of 208Pb and 197Au. Nucl Phys A (1970) 159:561–76. 10.1016/0375-9474(70)90727-X
52.
KawanoTChoYDimitriouPFilipescuDIwamotoNPlujkoVet alIAEA photonuclear data library 2019. Nucl Data Sheets (2020) 163:109–62. 10.1016/j.nds.2019.12.002
53.
GheorgheIGorielySWagnerNAumannTBaumannMvan BeekPet alPhotoneutron cross section measurements on 208Pb in the giant dipole resonance region. Phys Rev C (2024) 110:014619. 10.1103/PhysRevC.110.014619
54.
ToyokawaHGokoSHoharaSKaihoriTKanekoFKurodaRet alRecent progress in generation and application of AIST laser-Compton gamma-ray beam. Nucl Instrum Methods A (2009) 608:S41–3. 10.1016/j.nima.2009.05.062
55.
WellerHRAhmedMWGaoHTornowWWuYKGaiMet alResearch opportunities at the upgraded HIγS facility. Prog Part Nucl Phys (2009) 62:257–303. 10.1016/j.ppnp.2008.07.001
56.
AmanoSHorikawaKIshiharaKMiyamotoSHayakawaTShizumaTet alSeveral-MeV γ-ray generation at NewSUBARU by laser Compton backscattering. Nucl Instrum Methods A (2009) 602:337–41. 10.1016/j.nima.2009.01.010
57.
HorikawaKMiyamotoSAmanoSMochizukiT. Measurements for the energy and flux of laser compton scattering γ-ray photons generated in an electron storage ring: NewSUBARU. Nucl Instrum Methods A (2010) 618:209–15. 10.1016/j.nima.2010.02.259
58.
AhrensJBorchertHCzockKEpplerHGimmHGundrumHet alTotal nuclear photon absorption cross sections for some light elements. Nucl Phys A (1975) 251:479–92. 10.1016/0375-9474(75)90543-6
59.
HaradaHFurutakaKOhgakiHToyokawaH. Fine structure of giant resonance in the 28Si(γ, abs) reaction. J Nucl Sci Technol (2001) 38:465–9. 10.1080/18811248.2001.9715055
60.
SavranDLindenbergKGloriusJLöherBMüllerSPietrallaNet alThe low-energy photon tagger NEPTUN. Nucl Instrum Methods A (2010) 613:232–9. 10.1016/j.nima.2009.11.038
61.
SchelhaasKHennebergJSanzone-ArenhövelMWieloch-LaufenbergNZurmühlUZieglerBet alNuclear photon scattering by 208Pb. Nucl Phys A (1988) 489:189–224. 10.1016/0375-9474(88)90149-2
62.
PenfoldASLeissJE. Analysis of photonuclear cross sections. Phys Rev (1959) 114:1332–7. 10.1103/PhysRev.114.1332
63.
RasulovaFAAksenovNVAlekseevSIBelyshevSSChuprakovIFursovaNYet alPhotonuclear reactions on stable isotopes of molybdenum at bremsstrahlung endpoint energies of 10–23 MeV. Phys Rev C (2025) 111:024604. 10.1103/PhysRevC.111.024604
64.
AgostinelliSAllisonJAmakoKApostolakisJAraujoHArcePet alGEANT4—a simulation toolkit. Nucl Instrum Methods A (2003) 506:250–303. 10.1016/S0168-9002(03)01368-8
65.
von Neumann-CoselPHuxelNRichterASpielerCCarrollJCollinsC. Absolute calibration of low energy, thick target bremsstrahlung. Nucl Instrum Methods A (1994) 338:425. 10.1016/0168-9002(94)91326-9
66.
NakamuraTVinodkumarAMSugimotoTAoiNBabaHBazinDet alObservation of strong low-lying E1 strength in the two-neutron halo nucleus 11Li. Phys Rev Lett (2006) 96:252502. 10.1103/PhysRevLett.96.252502
67.
LeistenschneiderAAumannTBoretzkyKCortinaDCubJPramanikUDet alPhotoneutron cross sections for unstable neutron-rich oxygen isotopes. Phys Rev Lett (2001) 86:5442–5. 10.1103/PhysRevLett.86.5442
68.
TamiiAFujitaYMatsubaraHAdachiTCarterJDozonoMet alMeasurement of high energy resolution inelastic proton scattering at and close to zero degrees. Nucl Instrum Methods A (2009) 605:326–38. 10.1016/j.nima.2009.03.248
69.
NevelingRFujitaHSmitFAdachiTBergGButheleziEet alHigh energy-resolution zero-degree facility for light-ion scattering and reactions at iThemba LABS. Nucl Instrum Methods A (2011) 654:29–39. 10.1016/j.nima.2011.06.077
70.
ZilgesABalabanskiDIsaakJPietrallaN. Photonuclear reactions – from basic research to applications. Prog Part Nucl Phys (2022) 122:103903. 10.1016/j.ppnp.2021.103903
71.
KnüpferWRichterA. Effect of the nuclear magnetic susceptibility on the evaluation of the electric dipole polarizability. Phys Lett B (1981) 107:325–6. 10.1016/0370-2693(81)90333-6
72.
HeydeKvon Neumann-CoselPRichterA. Magnetic dipole excitations in nuclei: elementary modes of nucleonic motion. Rev Mod Phys (2010) 82:2365–419. 10.1103/RevModPhys.82.2365
73.
PietrallaN. Photons and nuclear structure. Eur Phys J A (2024) 60:108. 10.1140/epja/s10050-024-01285-8
74.
SuzukiT. Polarization transfer invariants at forward angle and spin-dependent excitations in nuclei. Prog Theo Phys (2000) 103:859–64. 10.1143/PTP.103.859
75.
GrossWMeuerDRichterASpamerETitzeOKnüpferW. A strong magnetic dipole excitation in 40Ca observed in high-resolution inelastic electron scattering and coherent spin-flip transitions due to ground-state correlations. Phys Lett B (1979) 84:296–300. 10.1016/0370-2693(79)90043-1
76.
MartinDvon Neumann-CoselPTamiiAAoiNBassauerSBertulaniCAet alTest of the Brink-Axel hypothesis for the pygmy dipole resonance. Phys Rev Lett (2017) 119:182503. 10.1103/PhysRevLett.119.182503
77.
BassauerSvon Neumann-CoselPReinhardP-GTamiiAAdachiSBertulaniCAet alElectric and magnetic dipole strength in 112,114,116,118,120,124Sn. Phys Rev C (2020) 102:034327. 10.1103/PhysRevC.102.034327
78.
BirkhanJMiorelliMBaccaSBassauerSBertulaniCAHagenGet alElectric dipole polarizability of 48Ca and implications for the neutron skin. Phys Rev Lett (2017) 118:252501. 10.1103/PhysRevLett.118.252501
79.
BassauerSvon Neumann-CoselPReinhardPGTamiiAAdachiSBertulaniCAet alEvolution of the dipole polarizability in the stable tin isotope chain. Phys Lett B (2020) 810:135804. 10.1016/j.physletb.2020.135804
80.
GorielySPéruSColòGRoca-MazaXGheorgheIFilipescuDet alE1 moments from a coherent set of measured photoneutron cross sections. Phys Rev C (2020) 102:064309. 10.1103/PhysRevC.102.064309
81.
StrauchSvon Neumann-CoselPRangacharyuluCRichterASchriederGSchwedaKet alGiant resonances in the doubly magic nucleus 48Ca from the (e, e′n) reaction. Phys Rev Lett (2000) 85:2913–6. 10.1103/PhysRevLett.85.2913
82.
O’KeefeGThompsonMAssafiriYPywellRShodaK. The photonuclear cross sections of 48Ca. Nucl Phys A (1987) 469:239–52. 10.1016/0375-9474(87)90108-4
83.
FultzSCBermanBLCaldwellJTBramblettRLKellyMA. Photoneutron cross sections for Sn116, Sn117, Sn118, Sn119, Sn120, Sn124, and indium. Phys Rev (1969) 186:1255–70. 10.1103/PhysRev.186.1255
84.
LeprêtreABeilHBergèreRCarlosPDe MiniacAVeyssièreAet alA study of the giant dipole resonance of vibrational nuclei in the 103 ≤ A ≤ 133 mass region. Nucl Phys A (1974) 219:39–60. 10.1016/0375-9474(74)90081-5
85.
UtsunomiyaHGorielySKamataMKondoTItohOAkimuneHet alγ-ray strength function for 116,117Sn with the pygmy dipole resonance balanced in the photoneutron and neutron capture channels. Phys Rev C (2009) 80:055806. 10.1103/PhysRevC.80.055806
86.
BraccoALanzaETamiiA. Isoscalar and isovector dipole excitations: nuclear properties from low-lying states and from the isovector giant dipole resonance. Prog Part Nucl Phys (2019) 106:360–433. 10.1016/j.ppnp.2019.02.001
87.
SavranDAumannTZilgesA. Experimental studies of the pygmy dipole resonance. Prog Part Nucl Phys (2013) 70:210–45. 10.1016/j.ppnp.2013.02.003
88.
MüscherMWilhelmyJMassarczykRSchwengnerRGriegerMIsaakJet alHigh-sensitivity investigation of low-lying dipole strengths in 120Sn. Phys Rev C (2020) 102:014317. 10.1103/PhysRevC.102.014317
89.
KrumbholzAvon Neumann-CoselPHashimotoTTamiiAAdachiTBertulaniCet alLow-energy electric dipole response in 120Sn. Phys Lett B (2015) 744:7. 10.1016/j.physletb.2015.03.023
90.
RusevGSchwengnerRBeyerRErhardMGrosseEJunghansARet alEnhanced electric dipole strength below particle-threshold as a consequence of nuclear deformation. Phys Rev C (2009) 79:061302. 10.1103/PhysRevC.79.061302
91.
MarkovaMvon Neumann-CoselPLitvinovaE. Systematics of the low-energy electric dipole strength in the Sn isotopic chain. Phys Lett B (2025) 860:139216. 10.1016/j.physletb.2024.139216
92.
LanzaEPellegriLVitturiAAndrésM. Theoretical studies of pygmy resonances. Prog Part Nucl Phys (2023) 129:104006. 10.1016/j.ppnp.2022.104006
93.
PaarNVretenarDKhanEColòG. Exotic modes of excitation in atomic nuclei far from stability. Rep Prog Phys (2007) 70:R02. 10.1088/0034-4885/70/5/R02
94.
BertulaniCA. Pygmy resonances and symmetry energy. Eur Phys J A (2019) 55:240. 10.1140/epja/i2019-12883-2
95.
CarboneAColòGBraccoACaoL-GBortignonPFCameraFet alConstraints on the symmetry energy and neutron skins from pygmy resonances in 68Ni and 132Sn. Phys Rev C (2010) 81:041301. 10.1103/PhysRevC.81.041301
96.
KlimkiewiczAPaarNAdrichPFallotMBoretzkyKAumannTet alNuclear symmetry energy and neutron skins derived from pygmy dipole resonances. Phys Rev C (2007) 76:051603(R). 10.1103/PhysRevC.76.051603
97.
KrasznahorkayAAkimuneHvan den BergABlasiNBrandenburgSCsatlósMet alNeutron-skin thickness in neutron-rich isotopes. Nucl Phys A (2004) 731:224–34. 10.1016/s0375-94740490024-3
98.
MarkovaMLarsenACTvetenGMvon Neumann-CoselPEriksenTKBello GarroteFLet alNuclear level densities and γ-ray strength functions of 111,112,113Sn isotopes studied with the Oslo method. Phys Rev C (2023) 108:014315. 10.1103/PhysRevC.108.014315
99.
MarkovaMLarsenACvon Neumann-CoselPBassauerSGörgenAGuttormsenMet alNuclear level densities and γ-ray strength functions in 120,124Sn isotopes: impact of Porter-Thomas fluctuations. Phys Rev C (2022) 106:034322. 10.1103/PhysRevC.106.034322
100.
MarkovaMvon Neumann-CoselPLarsenACBassauerSGörgenAGuttormsenMet alComprehensive test of the Brink-Axel hypothesis in the energy region of the pygmy dipole resonance. Phys Rev Lett (2021) 127:182501. 10.1103/PhysRevLett.127.182501
101.
MarkovaMLarsenACvon Neumann-CoselPLitvinovaEChoplinAGorielySet alSystematic study of the low-lying electric dipole strength in Sn isotopes and its astrophysical implications. Phys Rev C (2024) 109:054311. 10.1103/PhysRevC.109.054311
102.
RepkoANesterenkoVOKvasilJReinhardP-G. Systematics of toroidal dipole modes in Ca, Ni, Zr, and Sn isotopes. Eur Phys J A (2019) 55:242. 10.1140/epja/i2019-12770-x
103.
von Neumann-CoselPNesterenkoVOBrandhermIVishnevskiyPIReinhardP-GKvasilJet alCandidate toroidal electric dipole mode in the spherical nucleus 58Ni. Phys Rev Lett (2024) 133:232502. 10.1103/PhysRevLett.133.232502
104.
LeprêtreABeilHBergèreRCarlosPFagotJDe MiniacAet alMeasurements of the total photonuclear cross sections from 30 MeV to 140 Mev for Sn, Ce, Ta, Pb and U nuclei. Nucl Phys A (1981) 367:237–68. 10.1016/0375-9474(81)90516-9
105.
LevingerJS. The high energy nuclear photoeffect. Phys Rev (1951) 84:43–51. 10.1103/PhysRev.84.43
106.
Roca-MazaXViñasXCentellesMAgrawalBKColòGPaarNet alNeutron skin thickness from the measured electric dipole polarizability in 68Ni, 120Sn, and 208Pb. Phys Rev C (2015) 92:064304. 10.1103/PhysRevC.92.064304
107.
ChadwickMBObloinskýPHodgsonPEReffoG. Pauli-blocking in the quasideuteron model of photoabsorption. Phys Rev C (1991) 44:814–23. 10.1103/PhysRevC.44.814
108.
LiZZNiuYFColòG. Toward a unified description of isoscalar giant monopole resonances in a self-consistent quasiparticle-vibration coupling approach. Phys Rev Lett (2023) 131:082501. 10.1103/PhysRevLett.131.082501
109.
LitvinovaE. Relativistic approach to the nuclear breathing mode. Phys Rev C (2023) 107:L041302. 10.1103/PhysRevC.107.L041302
110.
von Neumann-CoselPPonomarevVYRichterAWambachJ. Gross, intermediate and fine structure of nuclear giant resonances: Evidence for doorway states. Eur Phys J A (2019) 55:224. 10.1140/epja/i2019-12795-1
111.
KlüpfelPReinhardP-GBürvenichTJMaruhnJA. Variations on a theme by Skyrme: a systematic study of adjustments of model parameters. Phys Rev C (2009) 79:034310. 10.1103/PhysRevC.79.034310
112.
ErlerJKlüpfelPReinhardP-G. Misfits in Skyrme-Hartree-Fock. J Phys G (2010) 37:064001. 10.1088/0954-3899/37/6/064001
113.
NikšićTVretenarDFinelliPRingP. Relativistic Hartree-Bogoliubov model with density-dependent meson-nucleon couplings. Phys Rev C (2002) 66:024306. 10.1103/PhysRevC.66.024306
114.
NikšićTVretenarDRingP. Relativistic nuclear energy density functionals: adjusting parameters to binding energies. Phys Rev C (2008) 78:034318. 10.1103/PhysRevC.78.034318
115.
PiekarewiczJ. Pygmy resonances and neutron skins. Phys Rev C (2011) 83:034319. 10.1103/PhysRevC.83.034319
116.
ReinhardP-GRoca-MazaXNazarewiczW. Information content of the parity-violating asymmetry in 208Pb. Phys Rev Lett (2021) 127:232501. 10.1103/PhysRevLett.127.232501
117.
KaufmannSSimonisJBaccaSBillowesJBissellMLBlaumKet alCharge radius of the short-lived 68Ni and correlation with the dipole polarizability. Phys Rev Lett (2020) 124:132502. 10.1103/PhysRevLett.124.132502
118.
MiorelliMBaccaSHagenGPapenbrockT. Computing the dipole polarizability of 48Ca with increased precision. Phys Rev C (2018) 98:014324. 10.1103/PhysRevC.98.014324
119.
HuBJiangWMiyagiTSunZEkströmAForssénCet alAb initio predictions link the neutron skin of 208Pb to nuclear forces. Nat Phys (2022) 18:1196–200. 10.1038/s41567-022-01715-8
120.
VernonIGoldsteinMBowerR. Galaxy formation: Bayesian history matching for the observable universe. Statist Sci (2014) 29:81. 10.1214/12-STS412
121.
HanSSteinerAW. Tidal deformability with sharp phase transitions in binary neutron stars. Phys Rev D (2019) 99:083014. 10.1103/PhysRevD.99.083014
122.
ReinhardP-GRoca-MazaXNazarewiczW. Combined theoretical analysis of the parity-violating asymmetry for 48Ca and 208Pb. Phys Rev Lett (2022) 129:232501. 10.1103/PhysRevLett.129.232501
123.
PiekarewiczJ. Implications of PREX-2 on the electric dipole polarizability of neutron-rich nuclei. Phys Rev C (2021) 104:024329. 10.1103/PhysRevC.104.024329
124.
YükselEPaarN. Implications of parity-violating electron scattering experiments on 48Ca (CREX) and 208Pb (PREX-II) for nuclear energy density functionals. Phys Lett B (2023) 836:137622. 10.1016/j.physletb.2022.137622
125.
ReedBTFattoyevFJHorowitzCJPiekarewiczJ. Density dependence of the symmetry energy in the post–PREX-CREX era. Phys Rev C (2024) 109:035803. 10.1103/PhysRevC.109.035803
126.
OrceJN. Polarizability effects in atomic nuclei. Int J Mod Phys E (2020) 29:2030002. 10.1142/S0218301320300027
127.
OrceJN. New formulas for the (−2) moment of the photoabsorption cross section, σ−2. Phys Rev C (2015) 91:064602. 10.1103/PhysRevC.91.064602
128.
von Neumann-CoselP. Comment on New formulas for the (−2) moment of the photoabsorption cross section, σ−2. Phys Rev C (2016) 93:049801. 10.1103/PhysRevC.93.049801
129.
SteinerAPrakashMLattimerJEllisP. Isospin asymmetry in nuclei and neutron stars. Phys Rep (2005) 411:325–75. 10.1016/j.physrep.2005.02.004
130.
TianJCuiHZhengKWangN. Effect of Coulomb energy on the symmetry energy coefficients of finite nuclei. Phys Rev C (2014) 90:024313. 10.1103/PhysRevC.90.024313
131.
CentellesMRoca-MazaXViñasXWardaM. Nuclear symmetry energy probed by neutron skin thickness of nuclei. Phys Rev Lett (2009) 102:122502. 10.1103/PhysRevLett.102.122502
132.
SchlimmeSAulenbacherKBaunackSBergerNDenigADoriaLet alThe MESA physics program. EPJ Web Conf (2024) 303:06002. 10.1051/epjconf/202430306002
133
GalesSTanakaKBalabanskiDNegoitaFStutmanDTesileanuOet alThe extreme light infrastructure-nuclear physics (ELI-NP) facility: new horizons in physics with 10 PW ultra-intense lasers and 20 MeV brilliant gamma beams. Rep Prog Phys (2018) 81:094301. 10.1088/1361-6633/aacfe8
134.
TanakaKASpohrKMBalabanskiDLBalascutaSCapponiLCernaianuMOet alCurrent status and highlights of the ELI-NP research program. Matter Radiat Extremes (2020) 5:024402. 10.1063/1.5093535
135.
WangH-WFanG-TLiuL-XXuH-HShenW-QMaY-Get alCommissioning of laser electron gamma beamline SLEGS at SSRF. Nucl Sci Tech (2022) 33:87. 10.1007/s41365-022-01076-0
136.
GheorgheIUtsunomiyaHStopaniKFilipescuDAri-izumiTBelyshevSet alUpdated neutron-multiplicity sorting method for producing photoneutron average energies and resolving multiple firing events. Nucl Instrum Methods A (2021) 1019:165867. 10.1016/j.nima.2021.165867
137.
ToganoYNakamuraTKobayashiTAumannTBabaHBoretzkyKet alElectric dipole strength of 52Ca. Il Nuovo Cim C (2024) 47:30. 10.1393/ncc/i2024-24030-3
138.
AumannTBertulaniCADuerMGalatyukTObertelliAPaninVet alNuclear structure opportunities with GeV radioactive beams at FAIR. Phil Trans R Soc A (2024) 382:2275. 10.1098/rsta.2023.0121
Summary
Keywords
dipole polarizability, neutron skin thickness, symmetry energy, density functional theory, ab initio calculations
Citation
von Neumann-Cosel P and Tamii A (2025) Electric dipole polarizability constraints on neutron skin and symmetry energy. Front. Phys. 13:1629987. doi: 10.3389/fphy.2025.1629987
Received
16 May 2025
Accepted
16 June 2025
Published
22 August 2025
Volume
13 - 2025
Edited by
Masayuki Matsuzaki, Fukuoka University of Education, Japan
Reviewed by
Shuichiro Ebata, Saitama University, Japan
Praveen C Srivastava, Indian Institute of Technology Roorkee, India
Updates
Copyright
© 2025 von Neumann-Cosel and Tamii.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peter von Neumann-Cosel, vnc@ikp.tu-darmstadt.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.