Front. Phys. Frontiers in Physics Front. Phys. 2296-424X Frontiers Media S.A. 1487954 10.3389/fphy.2024.1487954 Physics Original Research Mixing of isoscalar and isovector characteristics in the low-energy dipole mode Inakura and Ebata 10.3389/fphy.2024.1487954 Inakura Tsunenori 1 2 * Ebata Shuichiro 3 * 1 Office of Institutional Research and Decision Support, Tokyo Institute of Technology, Tokyo, Japan 2 Laboratory for Zero-Carbon energy, Institute of Innovative Research, Tokyo Institute of Technology, Tokyo, Japan 3 Graduate School of Science and Engineering, Saitama University, Saitama, Japan

Edited by: Oscar Moreno, Complutense University of Madrid, Spain

Reviewed by: Chen Ji, Central China Normal University, China

Francesco Giovanni Celiberto, University of Alcalá, Spain

 *Correspondence: Tsunenori Inakura, inakura@gmail.com; Shuichiro Ebata, ebata@mail.saitama-u.ac.jp 06 12 2024 2024 12 1487954 29 08 2024 05 11 2024 Copyright © 2024 Inakura and Ebata. 2024 Inakura and Ebata

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.

We investigated isospin splitting in low-energy dipole (LED) states of spherical nuclei such as 40Ca, 90Zr, 132Sn, 208Pb, and several N = 50 isotones using self-consistent Hartree–Fock plus random phase approximation calculations. Our analysis of isovector dipole (IVD) and isoscalar dipole (ISD) strengths, along with transition densities, reveals a clear energy-dependent relationship between IS and IV modes in 40Ca and 90Zr. For 208Pb and 132Sn, LED states show mixed IS + IV characteristics due to different neutron and proton shell structures. In N = 50 isotones, E 1 modes exhibit varying IS and IV properties with smooth transitions as neutron excess increases. Our results suggest that compressional ISD strengths could provide valuable insights into the slope parameter of the nuclear equation of state. The observed dependence on nuclear shell structures and neutron–proton correlations highlights the need for precise measurements and further research in nuclear physics.

 electric dipole mode PDR isoscalar and isovector characteristics RPA energy density functional Hartree–Fock method section-at-acceptance Nuclear Physics​ 1 Introduction

Low-energy dipole (LED) excitation is one of the key ingredients to investigate nuclear properties. The LED states appear in low excitation energy compared with the excitation energy of giant dipole resonance (GDR), having sizable strengths up to several percentage of the Thomas–Reiche–Kuhn sum rule. They have been observed in wide nuclear mass regions, e.g., 48Ca [1], 68Ni [2], 90Zr [3, 4], 132Sn [5], and 208Pb [6]. The systematic calculations for the dipole mode ( E 1 ) [7–9] show that the LED state appears in many isotopes. The LED state on a neutron-rich nucleus has been called the pygmy dipole resonance, which is often interpreted as a collective oscillation of the neutron skin and the remnant core [10, 11] (skin-core oscillation). There have been many studies for the LED to probe the properties of neutron matter, nuclear symmetry energy, and the slope parameter in the equation of state of neutron matter. However, its nature is still under debate.

It is found from recent experiments [12–19] that the LED state has two components in its lower- and higher-energy regions. The lower LED has an isoscalar (IS) + isovector (IV) characteristic, whereas the higher LED has IV dominant characteristic. These characteristics are indicated by comparing the E 1 strengths induced by ( α , α ′ ) scattering for the IS characteristic and ( γ , γ ′ ) scattering for the IV characteristic. This phenomenon is called LED isospin splitting. Several theoretical calculations [20–31] using the (quasi-particle) random-phase approximation (RPA) show that the calculated LED state in magic nuclei has both IS and IV characteristics in its strengths and transition densities. However, the underlying structure of LED isospin splitting is not well-understood yet.

To investigate LED isospin splitting, we systematically study the LED states of several spherical nuclei and N = 50 isotones, using the RPA calculation, which describes the IV dipole (IVD) and the compressional IS dipole (ISD) modes. We obtain the pure IS characteristic (and pure IV character) in the E 1 modes of 40Ca and 90Zr, while the isospin-mixed LED strengths appear in N ≫ Z nuclei. Isospin-mixing is induced by proton and neutron shell structures near the Fermi levels. Our calculation for N = 50 isotone shows that the E 1 mode with the pure IS characteristic varies smoothly to the lower component in LED isospin splitting. The LED is expected to be a substitute for observing the slope parameter L of the nuclear equation of state (EOS). There have been many studies to investigate which observable LED has the strongest correlation with L , but we have not yet pinned down the value of L . We show that low-energy compressional ISD strengths might be a suitable probe of L .

The manuscript is organized as follows: Section 2 reviews briefly the Hartree–Fock (HF) and the RPA calculation. In Section 3, we calculate the LED states in 40Ca, 90Zr, 208Pb, and 132Sn and analyze the properties of LED isospin splitting. Changes in the LED properties in N = 50 isotones are shown. We mention the relation between the ISD strength and L . Conclusions are given in Section 4.

 2 Methods

We employ the HF + RPA approach to describe LED states. The numerical code used in the present calculation is a revised version of the code developed in [32].

We use the Skyrme effective interaction and the minimal nuclear energy density functional (EDF) [33]. The Skyrme interaction is an EDF often used in mean-field model calculations and is characterized by being expanded in terms of delta functions in space. The parameter set is determined according to each protocol. The adopted Skyrme interactions are SkM * [34], SLy4 [35], UNEDF1 [36], SkI2, and SkI3 [37] parameter sets. The minimal EDF has been recently presented as an EDF with the smallest number of possible phenomenological parameters. The SeaLL1 parameter set is used in this work. The pairing correlation is ignored in this paper because its contribution to LEDs is small, as shown in [9]. We only restrict our consideration to spherical nuclei to avoid complicating the discussion with nuclear deformation and pairing effects. Therefore, the HF + RPA approach is a sufficient model for our purpose, when we discuss the LED in closed-shell and spherical nuclei.

Here, two operators for IVD and compressional ISD are introduced to investigate isovector and isoscalar characteristics in LED. The IVD operator D ̂ is expressed as D ̂ = N A ∑ i ∈ p r i Y 1 Ω i − Z A ∑ i ∈ n r i Y 1 Ω i , where Y ( 1 ) denotes the spherical harmonics for the dipole mode ℓ = 1 . The index i ∈ p ( i ∈ n ) indicates that the sum runs over protons (neutrons). For the compressional ISD mode, the operator, D ̂ IS = ∑ i ∈ n , p r i 2 − 5 3 ⟨ r 2 ⟩ r i Y 1 Ω i , is adopted where 5 ⟨ r 2 ⟩ / 3 is introduced to eliminate the center-of-mass contribution.

The E 1 transition strength from the ground state | 0 〉 to the excited state | α 〉 by D ̂ or D ̂ IS in an even-even nucleus is B ( E 1 ; α ) = 〈 α | D ̂ | 0 〉 2 or B ( IS E 1 ; α ) = 〈 α | D ̂ IS | 0 〉 2 . The transition amplitude 〈 α | D ̂ | 0 〉 can be calculated within HF + RPA for the one-particle-one-hole excitations. Using them, the strength function is given by smearing out with the Lorentzian and is written by S E 1 , E = γ π ∑ α 1 E − E α 2 + γ 2 − 1 E + E α 2 + γ 2   ⟨ α | D ̂ | 0 ⟩ 2 , where E α is the energy of the excited state | α 〉 and the strength function is denoted as S ( I S E 1 ; α ) for D ̂ IS . We adopt 2 γ = 1 MeV throughout this paper.

The excited states are analyzed by using particle–hole (ph) contribution C p h and transition density δ ρ , which are obtained in the HF + RPA calculations. The ph contribution of an excitation from the i orbit to the m orbit is evaluated, with the forward and backward amplitudes X i ( r ) and Y i ( r ) , as C m i = X i m 2 − Y i m 2 , in which X i m is extracted from X i ( r ) as X i m = ∫ d r   ϕ m * r X i r and similar for Y i m . Here, ϕ m is the m th wave function and the spin index is omitted for simplicity. The transition density δ ρ is expressed as δ ρ r = ∑ i ∈ n , p ϕ i * r X i r + Y i * r ϕ i r .

The radial dipole transition density δ ρ L = 1 ( r ) is calculated from δ ρ as δ ρ L = 1 r = ∫ d Ω   r Y 1 Ω   δ ρ r .

The size of the RPA matrix is reduced by assuming the reflection symmetry of the ground state with respect to x = 0 , y = 0 , and z = 0 planes. We adopt the three-dimensional coordinate (3D) representation within a sphere of radius R box = 25 fm with a uniform mesh spacing of 0.6 fm.

Furthermore, we introduce an index for the isospin structure ( I iss ) of the excited state to quantify the isospin characteristic of the dipole mode. I iss is obtained from the transition densities of the neutron and proton at the state with an excitation energy E as I iss E = ∫ d r   Z δ ρ ν r ; E − N δ ρ π r ; E ∫ d r   Z δ ρ ν r ; E + N δ ρ π r ; E , where δ ρ π ( ν ) denotes the transition density on protons (neutrons). If I iss ( E ) ≃ 0 , it means that the transition density has almost the same distributions in both the proton and neutron, and the state with energy E has a pure IS characteristic. On the contrary, if I iss ( E ) ≃ 1 , the transition densities are out of phase over the space r , and the state has a pure IV characteristic. Then, I iss ( E ) with less than or over 0.5 means the state has an IS or IV dominant characteristic, and the I iss ( E ) ∼ 0.5 state has an IS + IV mixed characteristic.

 3 Results and discussion

We calculate the E1 strengths and transition densities of the IVD and ISD modes, focusing mainly on their LED states. First, we show the LED modes in 40Ca and 90Zr, which hold isospin symmetry in their excited states. Next, we investigate the LED states in both stable (208Pb) and unstable (132Sn) nuclei, where N ≫ Z , demonstrating the emergence of isospin mixing in the excited states. We also investigate the interaction dependence of the LED properties for 90Zr and 132Sn. Then, we show how the LED properties change in N = 50 isotones: 100Sn, 90Zr, 82Ge, 78Ni, and 70Ca. Finally, we evaluate the relations between the slope parameter L and some quantities of dipole strengths. We can guess that the strength of compressional ISD modes is a more suitable value to probe the slope parameter than the IVD strength, dipole polarization, and partial cross sections of the LED.

3.1 <sup>40</sup>Ca and <sup>90</sup>Zr

For the N = Z nucleus, 40Ca, Figures 1A, B show the IVD and ISD transition strengths, B ( E 1 ; E ) and B ( IS E 1 ; E ) , and the strength functions S ( E 1 ; E ) and S ( IS E 1 ; E ) . They are calculated with the SkM * interaction. The E 1 mode with a large ISD strength appears at the excitation energy 7.2 MeV, which is the lowest in the LED states, while it has a negligibly small IVD strength. The radial transition densities r 2 δ ρ L = 1 of the lowest state are plotted in Figure 1C. Proton ( π ) and neutron ( ν ) transition densities, δ ρ L = 1 π and δ ρ L = 1 ν , are in phase and almost the same, indicating pure IS characteristic. Its index for isospin structure I iss = 0.044 also signifies the pure IS characteristic. A coherent superposition of proton and neutron excitations π , ν : 1 d 3 / 2 → 2 p 3 / 2 and 1 d 5 / 2 → 1 f 7 / 2 generates the state at 7.2 MeV. This isospin-symmetric mechanism is a characteristic transition on the E 1 excitation in N = Z nuclei, due to the similar shell structures near the Fermi levels in both the proton and neutron. It should be noted that the low-energy ISD mode is also obtained in another isospin-symmetric N = Z stable nucleus, 16O [21, 31].

IVD and ISD strength and strength functions of 40Ca are shown in (A, B), respectively. The arrow indicates the neutron threshold energy. (C) Neutron (solid) and proton (dashed) transition densities r 2 δ ρ L = 1 of the E 1 mode at 7.2 MeV.

Figures 2A, B show E 1 strengths for 90Zr, same as Figure 1. Similar to 40Ca, the E 1 state at 7.7 MeV has a large ISD strength and a negligibly small IVD strength. Its transition density is displayed in Figure 2C, which seems to have a pure IS characteristic and indicates I iss = 0.030 . The state is generated by coherent superposition of proton and neutron excitations, π , ν : 1 f 5 / 2 → 2 d 5 / 2 and 2 p 1 / 2 → 3 s 1 / 2 , as listed in Table 1 where we list large contributions C ph of the decomposed ph excitations. The highest occupied neutron orbit is 1 g 9 / 2 , which only contributes to the LED state at 7.7 MeV. The E 1 transition between 1 g 9 / 2 and 1 h 11 / 2 is the one major shell excitation and a node-less excitation, which needs nodes to be orthogonal to the center-of-mass mode. Then, its contribution C ph gets small, C ph = 0.052 (Table 1). Therefore, although 90Zr has 10 more neutrons than protons, its LED properties are almost the same as those of the N = Z nuclide (40Ca), which has a large ISD strength.

Same as Figure 1 but for 90Zr. (A) IVD strength, (B) ISD strength, (C, D) transition densities of the 7.7, and 11.9-MeV states, respectively.

Calculated ph contributions C ph of LED states in 90Zr.

7.7-MeV state ( I iss = 0.030)
ν : 1 f 5 / 2 → 2 d 5 / 2 0.377 π : 2 p 1 / 2 → 3 s 1 / 2 0.117
ν : 2 p 1 / 2 → 3 s 1 / 2 0.220 π : 1 f 5 / 2 → 2 d 5 / 2 0.056
ν : 1 g 9 / 2 → 1 h 11 / 2 0.052 π : 1 f 7 / 2 → 1 g 9 / 2 0.038
ν : 1 f 5 / 2 → 2 d 3 / 2 0.034 π : 1 f 7 / 2 → 1 g 7 / 2 0.011
11.9-MeV state ( I iss = 0.856)
ν : 2 p 3 / 2 → 3 s 1 / 2 0.189 π : 2 p 3 / 2 → 3 s 1 / 2 0.096
ν : 2 p 3 / 2 → 2 d 3 / 2 0.082 π : 2 p 3 / 2 → 2 d 3 / 2 0.082
ν : 2 p 3 / 2 → 2 d 5 / 2 0.044 π : 2 p 3 / 2 → 2 d 5 / 2 0.057
ν : 2 p 1 / 2 → 2 d 3 / 2 0.014 π : 2 p 1 / 2 → 2 d 3 / 2 0.048

The LED states with large IVD and small ISD strengths appear at approximately 12 MeV. Three peaks appear in the IVD strengths, but their transition densities and C ph s are similar to each other. The transition densities of the E 1 mode at 11.9 MeV are shown in Figure 2D. The aspect of the out-of-phase δ ρ ν and δ ρ π indicates a result of the fragmentation of the GDR due to their two bumps. The E 1 mode has the transition densities of a typical pure IV characteristic, also shown in its I iss = 0.856 .

The E 1 modes with pure IV and IS characteristics with energetic differences appear in the LED strengths of 90Zr. The LED modes with large ISD and small IVD strengths appear in lower energy (at 7.2 MeV), and the LED modes with the opposite E 1 strength distribution in higher energy (at approximately 12 MeV). This property of strength is consistent with that of LED isospin splitting. Furthermore, similar results can be obtained in other interactions shown in Figure 3. For SLy4, UNEDF1, and SeaLL1 interactions, the lower-LED modes with large ISD and negligibly small IVD strengths appear at ≈ 8 − 10 MeV, and the higher-LED modes with large IVD and small ISD strengths appear around the neutron threshold energy.

IVD and ISD strengths in 90Zr are displayed in upper and lower panels, respectively. The results were obtained using SLy4, UNEDF1, and SeaLL1 parameter sets.

The mechanism to induce LED is the same for all interactions in this work. The interaction dependence appears in the strengths and energy position of the LED, which is a result of the chemical potential and nuclear shell structure.

 3.2 <sup>208</sup>Pb and <sup>132</sup>Sn

Here, we show E 1 strengths and transition densities of LED states in N ≫ Z nuclei: 208Pb and 132Sn. Figure 4 shows the E 1 strengths and transition densities in 208Pb, similar to Figure 1, calculated with the SkM * interaction. The lowest E 1 mode at 6.9 MeV has both large ISD and IVD strengths, which is clearly different from 40Ca and 90Zr. The second-lowest E 1 mode appears at 7.6 MeV, with a large IVD strength and small but non-negligible ISD strength. Their transition densities of 6.9 MeV and 7.6 MeV states, displayed in Figures 4C, D, can be interpreted as the skin-core oscillation. This E 1 mode of 208Pb differs from that of N = Z nuclei, in which protons and neutrons act coherently. The differences come from the nuclear structures near the Fermi levels of nucleons. The E 1 excitation mode in 208Pb is described by the superposition of neutron (proton) single ph excitations over the N = 126 ( Z = 82 ) gap. This mechanism makes it difficult to induce purely the in-phase transition densities of neutrons and protons, unlike 40Ca and 90Zr, and causes mixing of the ISD and IVD modes, namely, the IS + IV characteristic in the LED states. The ph contributions of the 6.9 and the 7.6 MeV states, listed in Table 2, are not much different. In both states, significant contributions from four neutron ph excitations, namely, 2 f 5 / 2 → 3 d 5 / 2 , 1 h 9 / 2 → 2 g 9 / 2 , 3 p 3 / 2 → 3 d 5 / 2 , and 2 f 5 / 2 → 3 d 3 / 2 , exhaust 50% of the total. I iss s of the 6.9 MeV and 7.6 MeV states are 0.237 and 0.417, respectively. Both modes should be interpreted as IS + IV characteristic modes, although the transition densities of the 7.6 MeV state seem to have a dominant IV characteristic trend.

Same as Figure 1 but for 208Pb. (A) IVD strength, (B) ISD strength, (C) transition density of the 6.9-MeV state, and (D) transition density of the 7.6-MeV state.

Calculated ph contribution C ph of LED states in 208Pb.

6.9-MeV state ( I iss = 0.237)
ν : 2 f 5 / 2 → 3 d 5 / 2 0.250 π : 1 h 11 / 2 → 1 i 13 / 2 0.052
ν : 1 h 9 / 2 → 2 g 9 / 2 0.187 π : 1 h 11 / 2 → 1 i 11 / 2 0.034
ν : 1 i 13 / 2 → 1 j 15 / 2 0.103 π : 2 d 3 / 2 → 3 p 3 / 2 0.020
ν : 3 p 3 / 2 → 3 d 5 / 2 0.048 π : 1 g 7 / 2 → 2 f 7 / 2 0.014
ν : 2 f 5 / 2 → 3 d 3 / 2 0.031 π : 3 s 1 / 2 → 3 p 3 / 2 0.011
7.6-MeV state ( I iss = 0.417)
ν : 1 h 9 / 2 → 2 g 9 / 2 0.199 π : 2 d 5 / 2 → 1 h 9 / 2 0.038
ν : 3 p 3 / 2 → 3 d 5 / 2 0.173 π : 1 h 11 / 2 → 1 i 13 / 2 0.033
ν : 2 f 5 / 2 → 3 d 5 / 2 0.167 π : 1 h 11 / 2 → 1 i 11 / 2 0.030
ν : 3 p 3 / 2 → 3 d 3 / 2 0.088 π : 3 s 1 / 2 → 3 p 3 / 2 0.016
ν : 2 f 5 / 2 → 3 d 3 / 2 0.075 π : 2 d 5 / 2 → 2 f 7 / 2 0.013

Figure 5 shows the LED states of a neutron-rich unstable nucleus 132Sn which are calculated with the SkM * interaction. Similar to 208Pb, some E 1 modes appear around the neutron threshold energy with large IVD and ISD strengths. We pick up three E 1 states with large strengths at 7.2 MeV, 8.0 MeV, and 8.6 MeV. The 7.2-MeV state has the largest ISD strength lower than 10 MeV, and its I iss is 0.237. The 8.0-MeV state has large ISD and IVD strengths and I iss = 0.448 , indicating the IS + IV characteristic. The 8.6-MeV state has large IVD and small ISD strengths with I iss = 0.388 , which also indicates the IS + IV characteristic. At a glance, three states have the characteristics of IS, IS + IV, and IV dominances in their E 1 strengths, respectively, but their I iss signify IS + IV characteristics of three modes. To investigate the isospin characteristic of the E 1 state, not only the strength but also the transition density are necessary. It is not easy to separate pure IS and IV characteristics of the LED modes in N ≫ Z nuclei because different proton and neutron shell structures near the Fermi levels hardly induce in-phase transition densities of protons and neutrons, irrespective of the interactions used.

Same as Figure 1 but for 132Sn. (A) IVD strength, (B) ISD strength, (C–E) transition densities of the 7.2, 8.0, and 8.6-MeV states, respectively.

Figure 6 shows the LED strengths calculated with other interactions to compare with them. The LED modes in 132Sn occupy ∼ 2 % of the ISD energy-weighted sum rule value, regardless of which interaction is used. We can see similar LED strengths of the SkM * in other interactions, although some detail differences exist, similar to the case of 90Zr. From the E 1 strength distribution viewpoint, LED isospin splitting is less clear than SkM* in the results using SLy4, UNEDF1, SkI2, SkI3, and SeaLL1 due to large ISD and IVD strengths in the states. These interaction differences in LED isospin splitting are attributed to different structures of single-particle states near the Fermi levels and to the correlation between protons and neutrons, which induces the IS + IV characteristic of the LED modes.

Same as Figure 3 but for 132Sn. The results were obtained using SLy4, UNEDF1, SeaLL1, SkI2, and SkI3 parameter sets.

 3.3 Transition of LED isospin splitting in <inline-formula id="inf186"> <mml:math id="m194"> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>=</mml:mo> <mml:mn>50</mml:mn> </mml:mrow> </mml:math> </inline-formula> isotones

The LED strengths and its transition densities for 40Ca, 90Zr, and N ≫ Z nuclei are investigated in the above sections. For the LED states of 90Zr, we obtained LED isospin splitting, where the lower LED mode exhibits large ISD and small IVD strengths, while the higher LED mode shows the opposite. Furthermore, the transition densities of the lower and higher LED modes can take pure IS and pure IV characteristics due to the similar proton and neutron shell structures. On the other hand, in the N ≫ Z nuclei, the isospin characteristics of the LED modes mix due to the difference between proton and neutron shell structures near the Fermi levels. This is a general property, irrespective of the used interactions.

Here, we show the transition of LED isospin splitting in N = 50 isotones from N = Z to N ≫ Z . Figure 7 shows the E 1 strengths calculated with the SkM * interaction in N = 50 isotones, 100Sn, 90Zr, 82Ge, 78Ni, and 70Ca. Interestingly, the lowest ISD mode on the proton-drip-line nucleus 100Sn is essentially the same as that of the stable nucleus 90Zr, as well as with its higher LED mode at 12 MeV. The transition densities of the lower LED modes for N = 50 isotones are shown in Figure 8, and the higher LED modes, in Figure 9. Table 3 shows the excitation energies of lower and higher LED modes and I iss s. From the viewpoints of the strengths and I iss s, the E 1 states of 100Sn, 90Zr, and 82Ge have the characteristic of LED isospin splitting, which is also shown in their transition densities from Figures 8, 9. The higher LED has out-of-phase two-bump transition densities regarding the fragmentation of the GDR for 100Sn, 90Zr, and 82Ge in Figure 9. However, a small mixing of IS and IV characteristics is seen in the transition densities of the LED states in 82Ge [Figures 8C, 9C]. δ ρ π and δ ρ ν have different behaviors inside the nucleus ( r ≲ 5 fm). Considering 78Ni, the E 1 transition strengths at 8.3 MeV seem to have the IS dominant characteristic and its transition density and I iss indicate the IS + IV characteristic, like the case of 132Sn. The higher LED mode at 12.4 MeV has both large ISD and IVD strengths, but this is the neutron emission mode from the g 9 / 2 orbit to the continuum, which has a long tail of neutron transition density. For the LED modes of 70Ca, the ISD and IVD strength distributions are similar. The lower LED at 5.8 MeV, which is composed mainly of the neutron excitation from 1 f 5 / 2 to the resonant 2 d 3 / 2 orbit, has I iss ∼ 0.5 , indicating the IS + IV characteristic, while the higher LED at 10.9 MeV keeps the IV characteristic ( I iss = 0.86 ) .

E 1 strengths in N = 50 isotones, 100Sn, 90Zr, 82Ge, 78Ni, and 70Ca, with the SkM* interaction. For 90Zr, the figure shows the same results in Figures 2A, B.

Transition densities of lower LED modes, (A) at 8.2 MeV in 100Sn, (B) 7.7 MeV in 90Zr, (C) 8.6 MeV in 82Ge, (D) 8.3 MeV in 78Ni, and (E) 5.8 MeV in 70Ca.

Transition densities of higher LED modes, (A) at 11.9 MeV in 100Sn, (B) 11.9 MeV in 90Zr, (C) 12.4 MeV in 82Ge, (D) 12.1 MeV in 78Ni, and (E) 10.9 MeV in 70Ca.

I iss of lower and higher LED modes for N = 50 isotones: Excitation energy (MeV) and I iss of the lower LED mode and of the higher LED mode, for each isotone.

Lower LED E I iss Higher LED E I iss
100Sn 8.2 0.074 11.9 0.900
90Zr 7.7 0.030 11.9 0.856
82Ge 8.6 0.166 12.4 0.918
78Ni 8.3 0.330 12.1 0.517
70Ca 5.8 0.289 10.9 0.432

As one approaches from N = Z nuclei to N ≫ Z nuclei, the characteristic of the LED modes changes from the pure IS or pure IV characteristic to the IS + IV characteristic, although many of the higher LED modes in N ≫ Z drip-line nuclei are neutron emission modes.

 3.4 ISD strength and the slope parameter <italic>L</italic>

We investigate the relation between the low-energy E 1 strengths in 132Sn and the slope parameter L of EOS. The LED has attracted much attention because the LED observables strongly correlate with L [38, 39, 40]. However, the value of L has not yet been completely determined. The candidates observable have been suggested to access and determine the value of L . They are the summation of LED strength ∑ B ( E 1 ) , the LED cross section σ LED , and the LED dipole polarizability α LED . σ LED = 16 π 3 e 2 9 ℏ c ∫ 0 E dip d E   E   S E 1 , E , α LED = 8 π e 2 9 ∫ 0 E dip d E   S E 1 , E E , where E dip is defined as the energy corresponding to the minimum of the strengths S ( E 1 ; E ) and S ( I S E 1 ; E ) that exist between the LED peak and the GDR on each nucleus. The upper limit of the summation of LED strength is also E dip . The E dip values for IVD and ISD of 132Sn were all within 10 ± 1 MeV, as shown in Figure 6. Although the values depend on the smearing parameter γ , they hardly change in the range γ = 0.4 − 0.6 MeV. Table 4 lists the L values and the calculated LED observables, σ LED , α LED , and ∑ B ( E 1 ) , sorted by the value of L . In this work, we add the summation of low-energy compressional ISD strength ∑ B ( I S E 1 ) to the candidates.

Slope parameter L and LED observables for 132Sn; σ LED , α LED , and the summations of IVD and ISD strength ∑ E dip B ( E 1 ) , ∑ E dip B ( I S E 1 ) . The interactions are sorted by the order of L value.

Interaction L [MeV] σ LED [mb] α LED [ f m 2 / M e V ] ∑ B ( E 1 ) [ e 2 fm 2 ] ∑ B ( I S E 1 ) [ e 2 f m 6 ]
SeaLL1 32.4 33.8 0.555 1.076 1.25 × 1 0 4
UNEDF1 40.0 29.0 0.574 1.012 1.03 × 1 0 4
SkM * 45.8 26.4 0.402 0.808 1.18 × 1 0 4
SLy4 45.9 17.0 0.230 0.492 1.18 × 1 0 4
SkI3 100 94.5 1.078 2.508 1.34 × 1 0 4
SkI2 104 104 1.398 2.992 1.40 × 1 0 4

Table 4 shows linear relations between E 1 strength observables and values of L in the Skyrme parameter sets, excluding minimal EDF SeaLL1. The linear correlations between σ LED , α LED , and ∑ B ( E 1 ) and L can be confirmed, although the linearity deviates locally. However, the sequence of L values and ∑ B ( I S E 1 ) calculated with Skyrme interactions aligns in Figure 10. Its Pearson correlation coefficient is r = 0.93 and the standard deviation is σ = 624 [ e 2 f m 6 ] in this work. More detailed data are required, but this suggests that ∑ B ( I S E 1 ) might also be a suitable quantity for investigating L . D ̂ IS contains r 2 , so it can be expected to be an operator that probes the nuclear surface more than the IVD operator D ̂ . Therefore, the strength of the ISD of neutron-rich nuclei may be an important quantity for extracting L information.

∑ E dip B ( I S E 1 ) of 132Sn and the parameter L calculated with Skyrme interactions. The Pearson correlation coefficient is denoted as r . The shaded area means the standard deviation σ of the linear chi-squared fit.

 4 Conclusion

We investigated LED isospin splitting in spherical nuclei (40Ca, 90Zr, 132Sn, 208Pb, and several N = 50 isotones). Using fully self-consistent HF + RPA calculations with Skyrme interactions and minimal EDF in 3D Cartesian coordinates, we analyzed the isospin splitting mechanism through IVD and ISD strength functions and transition densities, introducing the isospin structure index I iss for analysis.

Their LED modes vary with excitation energy for 40Ca and 90Zr. Below 10 MeV, E 1 modes show large ISD and minimal IVD strengths, with transition densities in-phase for neutrons and protons (IS characteristic, I iss ∼ 0 ). In contrast, higher-energy E 1 modes exhibit out-of-phase densities (IV characteristic, I iss ∼ 1 ). This indicates a fundamental energy-position relation between IS and IV modes in 40Ca and 90Zr, as confirmed by HF + RPA calculations with various interactions.

We investigated the LED states of nuclei with more neutrons than protons, stable 208Pb, and unstable 132Sn. The states differ from those in nuclei with N = Z . In both cases, the LED modes below 9 MeV show both IS and IV characteristics, with an IS + IV mixed characteristic ( I iss ∼ 0.5 ) . This mixed characteristic is expected due to different neutron and proton shell structures near the Fermi levels. In 208Pb and 132Sn, the transition densities with the IS characteristic are less favored due to the distinct shell structures. We also examined how different interactions affect the LED states of 132Sn. Our HF + RPA results with SLy4 and SeaLL1 show unclear isospin splitting, whereas SkM * and UNEDF1 better reproduce the LED isospin splitting.

We calculated the LED states for spherical N = 50 isotones: 100Sn, 90Zr, 82Ge, 78Ni, and 70Ca. We found a smooth evolution in the transition densities of E1 modes with the largest strengths below 10 MeV. In 100Sn and 90Zr, the E 1 modes show pure IS or IV characteristic. In 82Ge, there is a slight mixing of IS and IV due to neutron excess. For 78Ni, the LED mode exhibits skin-core oscillation, while at 70Ca, it shows neutron emission. As we approach the neutron drip line, the IS characteristic of the LED modes gradually shifts to a mixed IS + IV characteristic.

Finally, we showed the behaviors of the low-energy E 1 strengths in 132Sn and the slope parameter L of EOS. It was shown that the sum of the compressional ISD strengths in the low-energy region might be an important quantity for exploring the slope parameter L . This is also conjectured from the form of the ISD operator, which includes a term that probes the nuclear surface.

LED isospin splitting is characterized by the isospin (IS, IV, and IS + IV) on the dipole transition, which strongly depends on the nuclear shell structures and the neutron–proton correlation. Therefore, to measure and analyze more precisely the isospin splitting will give the fundamental knowledge of the neutron–proton correlation, and pure IS or neutron characteristic might also be extracted. In the future work, we should proceed to study isospin splitting, considering the nuclear deformation and pairing correlation theoretically.

 Data availability statement

The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding authors.

 Author contributions

TI: writing–original draft and writing–review and editing. SE: writing–original draft and writing–review and editing.

 Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the MEXT Leading Initiative for Excellent Young Researchers Grant in Japan.

The authors thank K. Hagino for the fruitful discussions.

 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.

 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 Hartmann T Babilon M Kamerdzhiev S Litvinova E Savran D Volz S Microscopic nature of the pygmy dipole resonance: the stable Ca isotopes. Phys Rev Lett (2004) 93:192501. 10.1103/PhysRevLett.93.192501 Wieland O Bracco A Camera F Benzoni G Blasi N Brambilla S Search for the Pygmy Dipole Resonance in 68Ni at 600 MeV/nucleon. Phys Rev Lett (2009) 102:092502. 10.1103/PhysRevLett.102.092502 Schwengner R Rusev G Tsoneva N Benouaret N Beyer R Erhard M Pygmy dipole strength in 90Zr. Phys Rev C (2008) 78:064314. 10.1103/PhysRevC.78.064314 Iwamoto C Utsunomiya H Tamii A Akimune H Nakada H Shima T Separation of pygmy dipole and M1 resonance in 90Zr by a high-resolution inelastic proton scattering near 0°. Phys Rev Lett (2012) 108:262501. 10.1103/PhysRevLett.108.262501 Adrich P Klimkiewicz A Fallot M Boretzky K Aumann T Cortina-Gil D Evidence for pygmy and giant dipole resonances in 130Sn and 132Sn. Phys Rev Lett (2005) 95:132501. 10.1103/PhysRevLett.95.132501 Tamii A Poltoratska I von Neumann-Cosel P Fujita Y Adachi T Bertulani CA Complete electric dipole response and the neutron skin in 208Pb. Phys Rev Lett (2011) 107:062502. 10.1103/PhysRevLett.107.062502 Ebata S Nakatsukasa T Inakura T Yoshida K Hashimoto Y Yabana K . Canonical-basis time-dependent Hartree-Fock-Bogoliubov theory and linear-response calculations. Phys Rev C (2010) 82:034306. 10.1103/physrevc.82.034306 Inakura T Nakatsukasa T Yabana K . Emergence of pygmy dipole resonances: magic numbers and neutron skins. Phys Rev C (2011) 84:021302(R). 10.1103/physrevc.84.021302 Ebata S Nakatsukasa T Inakura T . Systematic investigation of low-lying dipole modes using the canonical-basis time-dependent Hartree-Fock-Bogoliubov theory. Phys Rev.C (2014) 90:024303. 10.1103/physrevc.90.024303 Vretenar D Paar N Ring P Lalazissis GA . Collectivity of the low-lying dipole strength in relativistic random phase approximation. Nucl Phys A (2001) 692:496–517. 10.1016/s0375-9474(01)00653-4 Paar N Vretenar D Khan E Colò G . Exotic modes of excitation in atomic nuclei far from stability. Rep Prog Phys (2007) 70:691. 10.1088/0034-4885/70/5/R02 Savran D Babilon M van den Berg AM Harakeh MN Hasper J Matic A Nature of the pygmy dipole resonance in 140Ce studied in (alpha, alpha' gamma) experiments. Phys Rev Lett (2006) 97:172502. 10.1103/PhysRevLett.97.172502 Endres J Savran D van den Berg AM Dendooven P Fritzsche M Harakeh MN Splitting of the pygmy dipole resonance in 138Ba and 140Ce observed in the (alpha, alpha’ gamma) reaction. Phys Rev C (2009) 80:034302. 10.1103/PhysRevC.80.034302 Endres J Litvinova E Savran D Butler PA Harakeh MN Harissopulos S Isospin character of the pygmy dipole resonance in 124Sn. Phys Rev Lett (2010) 105:212503. 10.1103/PhysRevLett.105.212503 Roming C Beller J Glorius J Isaak J Kelley JH Kwan E Low-lying dipole strength of the open-shell nucleus 94Mo. Phys Rev C (2013) 88:044331. 10.1103/PhysRevC.88.044331 Derya V Savran D Endres J Harakeh M Hergert H Kelley J Isospin properties of electric dipole excitations in 48Ca. Phys Lett B (2014) 730:288–92. 10.1016/j.physletb.2014.01.050 Negi D Wiedeking M Lanza EG Litvinova E Vitturi A Bark RA Nature of low-lying electric dipole resonance excitations in 74Ge. Phys Rev C (2016) 94:024332. 10.1103/PhysRevC.94.024332 Nakatsuka N Baba H Aumann T Avigo R Banerjee SR Bracco A Observation of isoscalar and isovector dipole excitations in neutron-rich 20O. Phys Lett B (2017) 768:387. 10.1016/j.physletb.2017.03.017 Savran D Derya V Bagchi S Endres J Harakeh MN Isaak J Multi-messenger investigation of the pygmy dipole resonance in 140Ce. Phys Lett B (2018) 786:16. 10.1016/j.physletb.2018.09.025 Paar N Niu YF Vretenar D Meng J . Isoscalar and isovector splitting of pygmy dipole structures. Phys Rev Lett (2009) 103:032502. 10.1103/physrevlett.103.032502 Papakonstantinou P Ponomarev VY Roth R Wambach J . Isoscalar dipole coherence at low energies and forbidden E1 strength. Eur Phys J A (2011) 47:14. 10.1140/epja/i2011-11014-7 Yüksel E Khan E Bozkurt K . Analysis of the neutron and proton contributions to the pygmy dipole mode in doubly magic nuclei. Nucl Phys A (2012) 877:35–50. 10.1016/j.nuclphysa.2012.01.006 Papakonstantinou P Hergert H Ponomarev VY Roth R . Low-energy dipole strength and the critical case of 48Ca. Phys Lett B (2014) 709:270–5. 10.1016/j.physletb.2012.02.024 Vretenar D Niu YF Paar N Meng J . Low-energy isovector and isoscalar dipole response in neutron-rich nuclei. Phys Rev C (2012) 85:044317. 10.1103/physrevc.85.044317 Nakada H Inakura T Sawai H . Crossover from skin mode to proton-neutron mode inE1excitations of neutron-rich nuclei. Phys Rev C (2013) 87:034302. 10.1103/physrevc.87.034302 Papakonstantinou P Hergert H Ponomarev VY Roth R . Low-energy electric dipole response of Sn isotopes. Phys Rev C (2014) 89:034306. 10.1103/physrevc.89.034306 Lanza EG Vitturi A Litvinova E Savran D . Dipole excitations via isoscalar probes: the splitting of the pygmy dipole resonance in124Sn. Phys Rev C (2014) 89:041601. 10.1103/physrevc.89.041601 Knapp F Lo Iudice N Veselý P Andreozzi F De Gregorio G Porrino A . Dipole response inSn132within a self-consistent multiphonon approach. Phys Rev C (2014) 90:014310. 10.1103/physrevc.90.014310 Papakonstantinou P Hergert H Roth R . Isoscalar and neutron modes in theE1spectra of Ni isotopes and the relevance of shell effects and the continuum. Phys Rev C (2015) 92:034311. 10.1103/physrevc.92.034311 Hamamoto I Sagawa H . Interplay between isoscalar and isovector correlations in neutron-rich nuclei. Phys Rev C (2017) 96:064312. 10.1103/physrevc.96.064312 Inakura T Togano Y . Skyrme random-phase approximation analysis of low-energy dipole states in oxygen isotopes. Phys Rev C (2018) 97:054330. 10.1103/physrevc.97.054330 Inakura T Imagawa H Hashimoto Y Mizutori S Yamagami M Matsuyanagi K . Mixed representation RPA calculation for octupole excitations on superdeformed states in the 40Ca and neutron-rich sulfur regions. Nucl Phys A (2006) 768:61–79. 10.1016/j.nuclphysa.2006.01.008 Bulgac A Forbes MM Jin S Perez RN Schunck N . Minimal nuclear energy density functional. Phys Rev C (2018) 97:044313. 10.1103/physrevc.97.044313 Bartel J Quentin P Brack M Guet C Håkansson HB . Towards a better parametrisation of Skyrme-like effective forces: a critical study of the SkM force. Nucl Phys A (1982) 386:79–100. 10.1016/0375-9474(82)90403-1 Chabanat E Bonche P Heenen P Meyer J Schaeffer R . A skyrme parametrization from subnuclear to neutron star densities part II. Nuclei far from stabilities Nucl Phys A (1998) 635:231. 10.1016/S0375-9474(98)00180-8 Kortelainen M McDonnell J Nazarewicz W Reinhard P-G Sarich J Schunck N Nuclear energy density optimization: large deformations. Phys Rev C (2012) 85:024304. 10.1103/physrevc.85.024304 Reinhard P-G Flocard H . Nuclear effective forces and isotope shifts. Nucl Phys A (1995) 584:467–88. 10.1016/0375-9474(94)00770-n Carbone A Col ò G Bracco A Cao L-G Bortignon PF Camera F Constraints on the symmetry energy and neutron skins from pygmy resonances in 68Ni and 132Sn. Phys Rev C (2010) 81:041301(R). 10.1103/PhysRevC.81.041301 Reinhard P-G Nazarewicz W . Information content of a new observable: the case of the nuclear neutron skin. Phys Rev C (2010) 81:051303. 10.1103/physrevc.81.051303 Rossi DM Adrich P Aksouh F Alvarez-Pol H Aumann T Benlliure J Measurement of the dipole polarizability of the unstable neutron-rich nucleus 68Ni. Phys Rev Lett (2013) 111:242503. 10.1103/PhysRevLett.111.242503