Phononparticle coupling effects in oddeven mass differences of semimagic nuclei.
Abstract
A method to evaluate the particlephonon coupling (PC) corrections to the singleparticle energies in semimagic nuclei, based on a direct solving the Dyson equation with PC corrected mass operator, is used for finding the oddeven mass difference between 18 even Pb isotopes and their oddproton neighbors. The Fayans energy density functional (EDF) DF3a is used which gives rather high accuracy of the predictions for these mass differences already on the meanfield level, with the average deviation from the existing experimental data equal to 0.389 MeV. It is only a bit worse than the corresponding value of 0.333 MeV for the Skyrme EDF HFB17 which belongs to a family of Skyrme EDFs with the highest overall accuracy in describing the nuclear masses. Account for the PC corrections induced by the lowlaying phonons and significantly diminishes the deviation of the theory from the data till 0.218 MeV.
Phononparticle coupling effects in oddeven mass differences of semimagic nuclei.\sodtitlePhononparticle coupling effects in oddeven mass differences of semimagic nuclei. \rauthorE. E. Saperstein, M. Baldo, S. S. Pankratov, S. V. Tolokonnikov \sodauthorE. E. Saperstein, M. Baldo, S. S. Pankratov, S. V. Tolokonnikov \datesAugust 14, 2019* \PACS21.60.n; 21.65.+f; 26.60.+c; 97.60.Jd
The singleparticle (SP) spectrum of a nucleus essentially influences different nuclear properties. Therefore ability to describe nuclear SP spectra correctly is very important for any selfconsistent nuclear theory. Till now, experimental SP levels are known in detail only for magic nuclei [1]. Their analysis in [2] on the base of the energy density functional (EDF) DF3a [3], which is a small modification of the original Fayans EDF DF3 [4, 5], led to rather good description of the data, significantly better than the predictions of the popular Skyrme EDF HFB17 [6] which belongs to the family of the EDFs HFB17 – HFB27 [7] possessing the highest accuracy among all the selfconsistent calculations in reproducing nuclear masses.
Inclusion of the particlephonon coupling (PC) corrections, with account for the nonpole diagrams [8, 9], made the agreement even better. Socalled approximation for the PC corrections was applied, being the creation vertex of the phonon. The corresponding diagrams for the PC correction to the mass operator , in the representation of the SP states , are displayed in Fig. 1. The first one is the usual pole diagram, with obvious notation, whereas the second one represents the sum of all nonpole diagrams of the order . The latter is often named “the phonon tadpole” [10], as an analogue of the tadpolelike diagrams in the field theory [11].
In magic nuclei, the perturbation theory in is valid [2] for solving the Dyson equation with the mass operator . Another situation is often occurs in semimagic nuclei [12], due to a strong mixture of some SP states with those possessing the structure of a SP state + phonon. For such cases a method is developed in [12] which is based on a direct solving the Dyson equation with the mass operator . Each SP state splits to a set of solutions with the SP strength distribution factors . In [12] a method is proposed how to express the average SP energy and the average factor in terms of and . It is similar to the one used usually for finding the corresponding experimental values [1]. The experimental data in heavy nonmagic nuclei considered in [12] are practically absent, therefore no comparison with experiment was made in that work.
Fortunately, there is a massive of data which has a direct relevance to the set of the solutions under discussion. We mean the oddeven mass differences, that is the “chemical potentials” in the notation of the theory of finite Fermi systems (TFFS) [13]:
(1) 
(2) 
(3) 
(4) 
where is the binding energy of the corresponding nucleus. Evidently, they are equal to one nucleon separation energies [14] taken with the opposite sign. For example, we have or .
Indeed, let us write down the Lehmann spectral expansion for the Green function in the representation of the functions which diagonalize [13]:
(5) 
with obvious notation. The isotopic index in (5) is for brevity omitted. In both the sums, the summation is carried out for the exact states of nuclei with one added or removed nucleon. Explicitly, if is the ground state of the eveneven () nucleus, the states in the first sum correspond to the () one for and () for . Correspondingly, in the second sum they are () for and (Z1,) for . If is a ground state of the corresponding odd nucleus, the corresponding pole in (5) coincides with one the chemical potentials (1) – (4). At the mean field level, they can be attributed to the SP energies with zero excitation energy, whereas with account for the PC corrections they should coincide with the corresponding energies .
The calculation scheme used in [12] and in this work contains some approximations which are typical for the TFFS [13] and which are valid in heavy nuclei with accuracy of , . At the mean field level, we relate the poles of the Green function of the eveneven nucleus () to the SP levels of its odd neighbors, () or (). Namely, the effect of the core deformation by the odd particle or hole to the mass operator is neglected. Note that, according to the TFFS scheme, this odd particle induced deformation and the corresponding, say, quadrupole moment may be found explicitly by solving the equation for the effective field [13], which is similar to that of the quasiparticle random phase approximation (QRPA), see e.g. [15, 16]. It should be stressed that the main odd particle effect is the PC one, and we consider it explicitly. Another approximation we use at the PC stage and it concerns the characteristics of the phonons we consider. Namely, we use the QRPA solution for the phonon in the () nucleus for finding the PC corrections to the SP characteristics of of these odd nuclei. Thereby, we neglect the “blocking effect” of the odd particle (hole) in the QRPA equation for . Accuracy of this approximation can be estimated as , where is the number of the particlehole states which contribute effectively to the vertex . We deal with strongly collective and states in the even lead isotopes, for which we may estimate this number as . A method to take into account the effect of the odd particle in the problem under consideration is developed in [17]. It is mainly important for lighter nuclei. Note also that there is a direct, in general more accurate, method to find the mass differences (1) – (4) in terms of the binding energies of each nuclei entering these relations. However, the calculation of the PC corrections to binding energies is rather cumbersome [9] and up to now there is no systematic corresponding calculations.
Let us describe briefly the method [12] to solve the PC corrected Dyson equation for the quasiparticle Green function. We should solve the following equation:
(6) 
where is the quasiparticle Hamiltonian with the spectrum and wave functions .
In the case when several phonons are taken into account, the total PC variation of the mass operator in Eq. (6) is the sum over all phonons:
(7) 
We deal with the normal subsystem of the semimagic nucleus under consideration, correspondingly, is the mass operator of a normal Fermi system. In this case, the explicit expression for the pole term is well known and can be found in [2, 12]. As to the nonlocal term, we follow to the method developed by Khodel [8], who first considered such diagrams in the problem of PC corrections in nuclei, see also [9].
All lowlying phonons we deal with are of surface nature, the surface peak dominating in their creation amplitude:
(8) 
The first term in this expression is surface peaked, whereas the involume addendum is rather small. It is illustrated in Fig. 2 for the and states in Pb. If one neglects this involume term , very simple expression for the nonpole term can be obtained [9]:
(9) 
Just as in [2, 12], we will below neglect the involume term in (8) and use Eq. (9) for the nonpole term of .
In this work, we consider the chain of even lead isotopes, Pb, with account for two lowlying phonons, and . Their excitation energies and the coefficients in Eq. (8) are presented in Table 1. Comparison with existing experimental data [18] is given. We present only 3 decimal signs only of the latter to avoid a cumbersomeness of the table. On the whole, the values agree with the data sufficiently well. In more detail, for the interval of Pb, the theoretical excitation energies of the states are visibly less than the experimental ones. This is a signal of the fact that our calculations overestimate the collectivity of these states and, correspondingly, the PC effect in these nuclei. The opposite situation where is for the lightest Pb isotopes, , where we, evidently, underestimate the PC effect. The value defines the amplitude, directly in fm, of the surface vibration in the nucleus under consideration. We see that in the most cases both the phonons we consider are strongly collective, with fm. At small values of , the vibration amplitude behaves as [9, 19]. Both the PC corrections to the SP energy, pole and nonpole, are proportional to . The ghost state is also taken into account, although the corresponding correction for nuclei under consideration is very small, see [12], because it depends on the mass number as , .
A  

180  1.415  1.168(1)  0.31  2.008  –  0.35 
182  1.284  0.888  0.31  1.836  –  0.35 
184  1.231  0.702  0.32  1.839  –  0.36 
186  1.133  0.662  0.33  1.881  –  0.34 
188  1.028  0.724  0.34  1.968  –  0.34 
190  0.930  0.774  0.36  2.052  –  0.33 
192  0.849  0.854  0.35  2.160  –  0.32 
194  0.792  0.965  0.35  2.272  –  0.32 
196  0.764  1.049  0.35  2.390  2.471(?)  0.31 
198  0.762  1.064  0.35  2.506  –  0.31 
200  0.789  1.027  0.30  2.620  –  0.31 
202  0.823  0.961  0.31  2.704  2.517  0.31 
204  0.882  0.899  0.22  2.785  2.621  0.31 
206  0.945  0.803  0.16  2.839  2.648  0.32 
208  4.747  4.086  0.33  2.684  2.615  0.09 
210  1.346  0.800  0.07  2.183  1.870(10)  0.19 
2.587  2.828(10)  0.17  
212  1.444  0.805  0.17  1.788  1.820(10)  0.36 
214  1.125  0.835(1)  0.19  1.469  –  0.37 
, MeV  
3  1  8.701  0.144 
2  7.270  0.667  
3  4.716  0.741  
4  2.078  0.194  
1  1  11.974  0.250 
2  9.952  0.895  
3  8.795  0.749  
4  7.357  0.127  
5  2.212  0.636  
6  0.466  0.272  
7  0.199  0.139  
8  2.750  0.481  
nucl.  DF3a  DF3a +  DF3a +  exp [20]  

Pb  , 1h  3.513  3.185  3.321  — 
, 3s  1.119  0.571  0.793  0.938(0.054)  
Pb  , 1h  2.942  2.564  2.695  — 
, 3s  1.610  1.023  1.268  1.316(0.021)  
Pb  , 1h  2.360  1.906  2.093  1.527(0.094) 
, 3s  2.104  1.450  1.727  1.753(0.022)  
Pb  , 1h  1.767  1.293  1.441  1.010(0.021) 
, 3s  2.592  1.906  2.152  2.213(0.032)  
Pb  , 1h  1.172  0.683  0.806  0.461(0.031) 
, 3s  3.072  2.356  2.561  2.661(0.019)  
Pb  , 1h  0.577  0.027  0.141  0.112(0.020) 
, 3s  3.543  2.750  2.945  3.103(0.023)  
Pb  , 1h  0.017  0.528  0.420  0.596(0.022) 
, 3s  4.005  3.265  3.440  3.572(0.020)  
Pb  , 1h  0.608  1.167  1.058  1.107(0.023) 
, 3s  4.461  3.673  3.838  4.019(0.024)  
Pb  , 1h  1.193  1.760  1.658  1.615(0.023) 
, 3s  4.911  4.111  4.268  4.494(0.025)  
Pb  , 1h  1.769  2.316  2.212  2.036(0.025) 
, 3s  5.358  4.569  4.716  4.999(0.031)  
Pb  , 1h  2.327  2.757  2.648  2.453(0.026) 
, 3s  5.806  5.177  5.317  5.480(0.039)  
Pb  , 1h  5.806  5.177  5.317  5.480(0.039) 
, 3s  6.258  5.612  5.753  6.050(0.018)  
Pb  , 1h  3.356  3.567  3.447  3.244(0.006) 
, 3s  6.717  6.357  6.493  6.637(0.003)  
Pb  , 1h  3.818  3.911  3.771  3.558(0.004) 
, 3s  7.179  6.976  7.105  7.254(0.003)  
Pb  , 1h  4.232  4.064  3.959  3.799(0.003) 
, 3s  7.611  7.778  7.633  8.004(0.007)  
Pb  , 1h  4.670  4.653  4.566  4.419(0.007) 
, 3s  8.030  7.971  8.055  8.379(0.010)  
Pb  , 1h  5.111  5.152  4.980  4.972(0.007) 
, 3s  8.446  8.276  8.481  8.758(0.044)  
Pb  , 1h  5.555  5.686  5.523  5.460(0.017) 
, 3s  8.857  8.620  8.865  9.254(0.029)  
0.385  0.321  0.218  
As the nonregular PC corrections to the SP energies we examine are important only for the states nearby the Fermi surface, we limit ourselves with a model space including two shells close to it, i.e., one hole and one particle shells, and besides we retain only the negative energy states. Note that for finding the pole term we use essentially wider SP space with energies MeV. To illustrate the method, we take for example the nucleus Pb . The space involves 5 hole states (, , , , ) and three particle ones (, , ). We see that there is here only one state for each value. Therefore, we need only diagonal elements of in Eq. (6). In the result, it reduces as follows:
(9) 
Details of finding the solutions of Eq. (9) can be found in [12]. In this notation, is just the index for the initial SP state from which the state originated. The corresponding SP strength distribution factors (factors) are:
(10) 
They should obey the normalization rule:
(11) 
Accuracy of fulfillment of this relation is a measure of the completeness of the model space we use to solve the problem under consideration. Two examples of the sets of solutions for four states in Pb are presented in Table 2. They originate from the first hole and the first particle states in the model space . In this case, our prescription for the oddeven mass differences, in accordance with the Lehmann expansion (5), is as follows:
(12) 
(13) 
As an illustration, we displayed in Fig. 3 the SP strength distributions (factors) of the similar two states nearby the Fermi level in the neighboring nucleus Pb.
The similar values for all the chain under consideration are given in Table 3. In the last line, the average deviation is given of the theoretical predictions from existing experimental data:
(14) 
where . For comparison, we calculated the corresponding value for the “champion” Skyrme EDF HFB17 [6] using the table [7] of the nuclear binding energies. It is equal to MeV. We see that accuracy of the Fayans EDF DF3a without PC in predicting the oddeven mass differences is only a bit worse. It agrees with the original Fayans’s idea [4, 5] develop an EDF without PC corrections. However, account for the PC corrections due to two lowlaying collective phonons makes agreement with the data significantly better.
To resume, a method, developed recently [12] to find the PC corrections to SP energies of semimagic nuclei based on the direct solution of the Dyson equation with the PC corrected mass operator, is used for finding the oddeven mass difference between the even Pb isotopes and their oddproton neighbors. The Fayans EDF DF3a is used for generating the mean field basis. On the meanfield level, the average accuracy of the predictions for the mass differences MeV is only a bit worse than that (0.333 MeV) for the Skyrme EDF HFB17 fitted to nuclear masses with highest accuracy among the selfconsistent calculations. Account for the PC corrections due to the lowlaying phonons and makes the agreement significantly better, MeV.
We acknowledge for support the Russian Science Foundation, Grants Nos. 161210155 and 161210161. The work was also partly supported by the RFBR Grant 160200228a. This work was carried out using computing resources of federal center for collective usage at NRC “Kurchatov Institute”, http://ckp.nrcki.ru. EES thanks the Academic Excellence Project of the NRNU MEPhI under contract by the Ministry of Education and Science of the Russian Federation No. 02. A03.21.0005.
References
 H. Grawe, K. Langanke, and G. MartínezPinedo, Rep. Prog. Phys. 70, 1525 (2007).
 N.V. Gnezdilov, I.N. Borzov, E.E. Saperstein, and S.V. Tolokonnikov, Phys. Rev. C 89, 034304 (2014).
 S.V. Tolokonnikov and E.E. Saperstein, Phys. At. Nucl. 73, 1684 (2010).
 A.V. Smirnov, S.V. Tolokonnikov, S.A. Fayans, Sov. J. Nucl. Phys. 48, 995 (1988).
 S.A. Fayans, S.V. Tolokonnikov, E.L. Trykov, and D. Zawischa, Nucl. Phys. A 676, 49 (2000).
 S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. Lett. 102, 152503 (2009).

S. Goriely, http://wwwastro.ulb.ac.be/bruslib/
nucdata/  V. A. Khodel’, Sov. J. Nucl. Phys. 24, 282 (1976).
 V. A. Khodel, E. E. Saperstein, Phys. Rep. 92, 183 (1982).
 S. Kamerdzhiev and E. E. Saperstein, EPJA 37, 333 (2008).
 S. Weinberg, Phys. Rev. Lett. 31, 494 (1973).
 E.E. Saperstein, M. Baldo, S.S. Pankratov, and S.V. Tolokonnikov, JETP Lett. 104, 609 (2016).
 A.B. Migdal Theory of finite Fermi systems and applications to atomic nuclei (Wiley, New York, 1967).
 A. Bohr and B.R. Mottelson, Nuclear Structure (Benjamin, New York, Amsterdam, 1969.), Vol. 1.
 S.V. Tolokonnikov, S. Kamerdzhiev, D. Voytenkov, S. Krewald, and E.E. Saperstein, Phys. Rev. C 84, 064324 (2011).
 S.V. Tolokonnikov, S. Kamerdzhiev, S. Krewald, E.E. Saperstein, and D. Voitenkov, Eur. Phys. J. A 48, 70 (2012).
 M. Baldo, P.F. Bortignon, G. Coló, D. Rizzo, and L. Sciacchitano, J. Phys. G: Nucl. Phys. 42, 085109 (2015).
 S. Raman, C.W. Nestor Jr., and P. Tikkanen, At. Data Nucl. Data Tables 78, 1 (2001).
 A. Bohr and B.R. Mottelson, Nuclear Structure (Benjamin, New York, 1974.), Vol. 2.
 M. Wang, G. Audi, A.H. Wapstra, F.G. Kondev, M. MacCormick, X. Xu, and B. Pfeiffer, Chinese Physics C, 36, 1603 (2012).