Including particlevibration coupling in the Fayans functional. 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 the direct solution of the Dyson equation with PC corrected mass operator, is presented. It is used for finding the oddeven mass difference between even Pb and Sn isotopes and their oddproton neighbors. The Fayans energy density functional (EDF) DF3a is used which gives rather highly accurate predictions for these mass differences already at the meanfield level. In the case of the lead chain, account for the PC corrections induced by the lowlaying phonons and makes agreement of the theory with the experimental data significantly better. For the tin chain, the situation is not so definite. In this case, the PC corrections make agreement better in the case of the addition mode but they spoil the agreement for the removal mode. We discuss the reason of such a discrepancy.
pacs:
21.60.Jz, 21.10.Ky, 21.10.Ft, 21.10.ReI Introduction
Ability to describe nuclear singleparticle (SP) spectra correctly is an important feature of any selfconsistent nuclear theory. Last decade is characterized with a revival of the interest to study the particlephonon coupling (PC) corrections to the SP spectra within different selfconsistent nuclear approaches. Up to now, they were limited only to magic nuclei, see a study of this problem within the relativistic meanfield theory LitvRing (), and several ones within the Skyrme–Hartree–Fock method Bort (); Dobaczewski (); BaldoPC (). The consideration of this problem levels () on the basis of the selfconsistent theory of finite Fermi systems (TFFS) scTFFS () differs from those cited above by the account for the nonpole terms of the PC correction to the mass operator . In the nuclear PC problem, the nonpole diagrams of the operator were considered firstly by Khodel Khod76 (), see also scTFFS (). These diagrams are also sometimes named “the phonon tadpole” phontad (), by the analogy with the tadpolelike diagrams in field theory tad ().
There are several reasons, why the magic nuclei were chosen as a polygon for studying the PC corrections to SP levels. First, the experimental SP levels are known in detail only for magic nuclei levexp (). Second, these nuclei are nonsuperfluid, which makes the formulas for the PC correction to the mass operator much simpler than their analogue in the superfluid case phontad (). At last, in these nuclei, the PC strength is rather week and the perturbation theory in is valid, being the creation vertex of the phonon. Moreover, the perturbation theory in can be used for solving the Dyson equation with the mass operator .
Another situation occurs often in semimagic nuclei levsemi (), due to a strong mixture of some SP states with those possessing the structure of ( + phonon). For such cases a method was used in levsemi () based on the direct solution of the Dyson equation with the mass operator . Such a method was developed by Ring and Werner in 1973 RingWerner (). As result of such solution each SP state splits to a set of solutions with the SP strength distribution factors . An approximate expression for the nonpole term of was used in levsemi (), just as in levels (), which is valid for heavy nuclei.
A method was 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 levexp (). Unfortunately, the experimental data in heavy nonmagic nuclei considered in levsemi () are practically absent. However, there is a massive set of data which has a direct relevance to the solutions under discussion, namely, the oddeven mass differences, that is the “chemical potentials” in terms of the TFFS AB (). This quantity was analyzed in in a brief letter chimpots () for the lead chain. In this article, we present a detailed analysis of these results and carry out the analogous calculations for the tin chain.
Ii Direct solution of the Dyson equation with PC corrected mass operator
In this article, just as in Refs. levels (); levsemi (); chimpots (), we use the selfconsistent basis generated with the energy density functional (EDF) by Fayans Fay1 (); Fay4 (); Fay ():
(1) 
which depends on the equal footing on the normal density and the anomalous one . Note, that from the very beginning Fayans with the coauthors supposed that the parameters of the EDF should be chosen in such a way that they take into account various PC corrections on average.
To find the SP energies with account for the PC effects, we solve the following equation:
(2) 
where is the quasiparticle Hamiltonian with the spectrum and is the PC correction to the quasiparticle mass operator. This is equivalent to the Dyson equation for the oneparticle Green function with the mass operator .
For brevity, we present below explicit formulas for one phonon. In the case when several phonons are taken into account, the total PC variation of the mass operator in Eq. (2) is the sum over all phonons:
(3) 
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 .
In magic nuclei, the vertex in Fig. 1 obeys the RPAlike TFFS equation AB ()
(4) 
where is the particlehole propagator and is the LandauMigdal (LM) interaction amplitude.
The expression for the pole term in magic nuclei is well known LitvRing (); Bort (); levels (), but we present it here for completeness. In obvious symbolic notation, the pole diagram corresponds to , where is the phonon function:
(5) 
where is the phonon excitation energy.
Explicitly, one obtains
(6)  
where stands for the occupation numbers.
As to the nonpole term, we follow to the method developed by Khodel Khod76 (), see also scTFFS (). The second, nonpole, term in Fig. 1 is
(7) 
where can be found by variation of Eq. (4) in the field of the phonon:
(8)  
All the lowlying phonons we deal with are of surface nature, the surface peak dominating in their creation amplitude:
(9) 
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 scTFFS ():
(10) 
Just as in levels (); levsemi (), below we will neglect the involume term in (9) and use Eq. (10) for the nonpole term of .
All the lowlying phonons we consider have natural parity. In this case, the vertex possesses even parity. It is a sum of two components with spins and , respectively,
(11) 
where stand for the usual spinangular tensor operators BM1 (). The operators and have opposite parities, hence the spin component should be the odd function of the excitation energy, . In all the cases we consider the component is negligible and the component only will be retained, the index ’0’ being for brevity omitted: .
As it was discussed in the introduction, any semimagic nucleus consists of two subsystems with absolutely different properties, the normal and superfluid ones. Following to levsemi (), we consider the SP spectra only of the normal subsystem, e.g. the proton spectra for the lead isotopes. In this case, all the above formulas remain valid with one exception. Now the vertex obeys the QRPAlike TFFS equation AB (),
(12) 
where all the terms are matrices. The angular momentum projection , which is written down in Eq. (6) explicitly, is here and below for brevity omitted. In the standard TFFS notation, we have:
(13) 
(14) 
In (13), is the normal component of the vertex , whereas are two anomalous ones. In Eq. (14), is the usual LM interaction amplitude which is the second variation derivative of the EDF over the normal density . The effective pairing interaction is the second derivative of the EDF over the anomalous density . At last, the amplitude stands for the mixed derivative of over and .
The matrix consists of integrals over of the products of different combinations of the Green function and two Gor’kov functions and AB ().
As we need the proton vertex and the proton subsystem is normal, only the normal vertex is nonzero in this case. This is the meaning of the abbreviated notation in (6) and below.
For solving the above equations, we use the selfconsistent basis generated by the version DF3a DF3a () of the Fayans EDF DF3 Fay4 (); Fay (). The nuclear meanfield potential is the first derivative of in (1) over .
In magic nuclei levels (), the perturbation theory in with respect to can be used to solve this equation:
(15) 
where
(16) 
Another situation is typical for semimagic nuclei, where small denominators can appear in Eq. (6). We demonstrate in Table I the poor applicability of the perturbation theory on the example of the Pb nucleus with the use of two phonons, MeV and MeV. The value of in the column 6 is found as a sum of the values in three previous columns multiplied, in accordance with Eq. (15), by the total factor in the last column.
The most dangerous denominator is now MeV. The state 2 is an obvious victim of the use of a nonadequate perturbation theory in this case. Indeed, the inequality is a necessary condition for applicability of the perturbation theory. For the 2 state, we have the value of 0.95 for the left hand side of this inequality, which is almost equal to unit. In the table, there are other examples of poor applicability of the perturbation theory, but they are not so demonstrative.
’ghost’  

1  1.21  1.01  0.32  0.04  0.37  1.58  0.64  0.82  0.56 
2  2.01  0.81  0.20  0.03  0.36  2.37  0.66  0.90  0.62 
1  3.36  0.43  0.27  0.04  0.09  3.45  0.74  0.98  0.73 
3  6.72  0.07  0.17  0.03  0.22  6.50  0.84  0.96  0.82 
2  7.51  2.02  0.17  0.03  0.10  7.40  0.05  0.97  0.05 
1  7.98  0.28  0.30  0.04  0.42  7.56  0.70  0.96  0.68 
2  8.80  0.33  0.18  0.03  0.26  8.541  0.57  0.80  0.50 
1  10.93  0.23  0.24  0.03  0.01  10.92  0.79  0.35  0.32 
For such cases, a method based on the direct solution of the Dyson equation (2) was developed in levsemi (). We describe it with the example of the same nucleus Pb which was considered above to demonstrate the nonapplicability of the perturbation theory. The model was suggested based on the fact that the PC corrections are important only for the SP levels nearby the Fermi surface. A model space was considered including two shells close to it, i.e. one hole and one particle shells. To avoid any misunderstanding, we stress that this restriction concerns only Eq. (2). In Eq. (6) for , we use rather wide SP space with energies MeV. The space involves 5 hole states (1g, 2d, 1h, 2d, 3s) and four particle ones (1g, 2f, 1i, 2f). We see that there is here only one state for each value. Therefore, we need only diagonal elements in (6), which simplifies very much the solution of the Dyson equation. As a consequence, Eq. (2) reduces as follows:
(17) 
The nonpole term does not depend on the energy, therefore only poles of Eq. (6) are the singular points of Eq. (17). They can be readily found from (6) in terms of and . It can be easily seen that the l.h.s of Eq. (17) always changes sign between any couple of neighboring poles, and the corresponding solution can be found with usual methods. In this notation, is just the index for the initial singleparticle state from which the state originated. The latter is a mixture of a singleparticle state with several particlephonon states. The corresponding singleparticle strength distribution factors (factors) can be found similar to (16):
(18) 
Evidently, they should obey the normalization rule:
(19) 
The accuracy of fulfillment of this relation is a measure of validity of the model space we use to solve the problem under consideration.
1  1  13.736  0.220 
2  11.596  0.777  
3  10.339  0.674  
4  8.862  0.484  
5  3.447  0.760  
6  2.217  0.199  
7  1.122  0.288  
3  1  9.877  0.608 
2  8.536  0.604  
3  6.493  0.839  
1  1  13.717  0.692 
2  11.690  0.256  
3  9.084  0.134  
4  7.509  0.814  
5  2.471  0.427  
6  1.095  0.473  
2  1  11.817  0.314 
2  11.150  0.139  
3  9.799  0.516  
4  8.580  0.312  
5  8.195  0.295  
6  7.404  0.171  
7  0.564  0.791  
A set of solutions for four states in the Pb nucleus is presented in Table II. The corresponding factors are displayed in Fig. 3. In three upper cases for a given there is a state with dominating value (). They are examples of “good” singleparticle states. In such cases, the following prescription looks natural for the PC corrected singleparticle characteristics:
(20) 
These are the analogous to Eqs. (15) and (16) in the perturbative solution.
The lowest panel in Fig. 3 represents a case of a strong spread where there are two or more numbers with comparable values of the spectroscopic factors . For such cases, the following generalization of Eq. (20) was suggested in levsemi ():
(21) 
where
(22) 
The average SP energies and the factors found according to the above prescriptions are presented in Table III.
1  1.21  0.14  1.07  0.97 

2  2.01  0.23  2.24  0.87 
1  3.36  0.17  3.19  0.96 
3  6.72  0.23  6.49  0.84 
2  7.51  0.05  7.46  0.94 
1  7.98  0.25  7.73  0.95 
2  8.80  0.17  8.63  0.92 
1  10.93  0.13  10.80  0.95 
Comparison of Table III with the perturbation theory Table I shows limited common features. In particular, the state 2 is now absolutely “healed” with a big value of the factor and a small shift of the PC corrected energy from the initial value.
Iii PC corrections to the oddeven mass differences of semimagic nuclei
As it was discussed in Introduction, in heavy semimagic nuclei the experimental data on the SP energies are practically absent, but there is another kind of observables, the oddeven mass differences, which has a direct relevance to the set of the solutions of the Dyson equation in the form (17). In terms of the (TFFS) AB () they are named “chemical potentials”:
(23) 
(24) 
(25) 
(26) 
where is the binding energy of the corresponding nucleus. Evidently, they are equal to one nucleon separation energies BM1 () taken with the opposite sign. For example, we have or .
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 
Indeed, let us write down the Lehmann spectral expansion for the Green function in the representation of the functions which diagonalize AB ():
(27) 
with obvious notation. The isotopic index in (27) 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 (27) coincides with one the chemical potentials (23) – (26). 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 .
nucl.  DF3a  DF3a +  DF3a +  exp mass ()  

Pb  1  3.513  3.185  3.321  — 
Pb  1  2.942  2.564  2.695  — 
Pb  2.360  1.906  2.093  1.527(0.094)  
Pb  1.767  1.293  1.441  1.010(0.021)  
Pb  1.172  0.683  0.806  0.461(0.031)  
Pb  0.577  0.027  0.141  0.112(0.020)  
Pb  0.017  0.528  0.420  0.596(0.022)  
Pb  0.608  1.167  1.058  1.107(0.023)  
Pb  1.193  1.760  1.658  1.615(0.023)  
Pb  1.769  2.316  2.212  2.036(0.025)  
Pb  2.327  2.757  2.648  2.453(0.026)  
Pb  5.806  5.177  5.317  5.480(0.039)  
Pb  3.356  3.567  3.447  3.244(0.006)  
Pb  3.818  3.911  3.771  3.558(0.004)  
Pb  4.232  4.064  3.959  3.799(0.003)  
Pb  4.670  4.653  4.566  4.419(0.007)  
Pb  5.111  5.152  4.980  4.972(0.007)  
Pb  5.555  5.686  5.523  5.460(0.017)  
0.465  0.261  0.252 
In the sets of solutions for four states in Pb presented in Table II, 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 (27), is as follows:
(28) 
(29) 
iii.1 The lead chain
In this subsection, we consider the chain of even lead isotopes, Pb, with account of the two lowlying phonons, and . Their excitation energies and the coefficients in Eq. (9) are presented in Table IV. Comparison with existing experimental data expomega () is given. We present only 3 decimal signs 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 smaller than the experimental ones.
This is a signal of the fact that our calculations may overestimate the collectivity of these states and, correspondingly, the PC effect in these nuclei. The opposite situation is present 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 scTFFS (); BM2 (). 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 levsemi (), because it depends on the mass number as , .
nucl.  DF3a  DF3a +  DF3a +  exp mass ()  

Pb  3s  1.119  0.571  0.793  0.938(0.054) 
Pb  3s  1.610  1.023  1.268  1.316(0.021) 
Pb  3s  2.104  1.450  1.727  1.753(0.022) 
Pb  3s  2.592  1.906  2.152  2.213(0.032) 
Pb  3s  3.072  2.356  2.561  2.661(0.019) 
Pb  3s  3.543  2.750  2.945  3.103(0.023) 
Pb  3s  4.005  3.265  3.440  3.572(0.020) 
Pb  3s  4.461  3.673  3.838  4.019(0.024) 
Pb  3s  4.911  4.111  4.268  4.494(0.025) 
Pb  3s  5.358  4.569  4.716  4.999(0.031) 
Pb  3s  5.806  5.177  5.317  5.480(0.039) 
Pb  3s  6.258  5.612  5.753  6.050(0.018) 
Pb  3s  6.717  6.357  6.493  6.637(0.003) 
Pb  3s  7.179  6.976  7.105  7.254(0.003) 
Pb  3s  7.611  7.778  7.633  8.004(0.007) 
Pb  3s  8.030  7.971  8.055  8.379(0.010) 
Pb  3s  8.446  8.276  8.481  8.758(0.044) 
Pb  3s  8.857  8.620  8.865  9.254(0.029) 
0.344  0.370  0.220 
The values, similar to that of (28) for all the chain under consideration, are given in Table V. In the last line, the average deviation is given of the theoretical predictions from existing experimental data:
(30) 
where . For comparison, we calculated the corresponding value for the “champion” Skyrme EDF HFB17 HFB17 () using the table HFBsite () of the nuclear binding energies. It is equal to MeV, that it is a bit worse than the DF3a value even without PC corrections. It agrees with the original Fayans’s idea Fay1 (); Fay () to 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. In this case, the phonon plays the main role, and additional improvement of the agreement due to the phonon is not essential.
Let us go to the values similar to that of (29). For all the lead chain, from Pb till Pb, they are given in Table VI. Again the last line contains the average deviation of the theoretical predictions from the experimental data found from the relation similar to (30) with the change of , and . In this case, the average of the HFB17 EDF predictions (HFB17)=0.253 MeV is better than that of the DF3a one, but again the account of the PC corrections permits to surpass the accuracy of HFB17. A relative role of the two phonons under consideration is now different from that in the case of , Table V. In the case of , addition of the corrections due to the phonon alone spoils the agreement a bit, but inclusion of both the phonons makes the agreement perfect. This confirms the statement made in the Introduction about a delicate status of the problem under consideration.
In conclusion of this subsection, we present the overall average deviation of the theoretical predictions from the experimental data, summed over both the and values, . In this case, on the meanfield level, the average accuracy of the predictions for the mass differences of the Fayans EDF MeV is worse than that (0.333 MeV) for the Skyrme EDF HFB17, but only a bit. Account for the PC corrections due to the lowlaying phonons and makes the agreement significantly better, MeV.
iii.2 The tin chain
Let us consider now the Sn chain, from Sn till Sn. To begin with, we consider validity of the perturbation theory recipe (15) for the nucleus Sn, which is in the middle of the chain. Table VII is completely similar to Table I for the nucleus Pb. As it is clear from the discussion of Table I, the factors contain the main information on the applicability of the perturbation theory for solving Eq. (2). Table VII contains all necessary information for this nucleus. The first seven columns contains a detailed information on the PC corrections on the SP energies, whereas the last three ones, on the factors. The values in columns 8, 9 and 10 are found from Eq. (16) by substituting there the values of , and their sum , correspondingly, in accordance with (3).
’ghost’  

1  1.86  1.99  0.54  0.06  0.53  2.39  0.57  0.54  0.38 
3  2.50  0.28  0.34  0.05  0.13  2.38  0.19  0.91  0.18 
2  2.64  6.20  0.36  0.08  0.19  2.83  0.03  0.92  0.03 
1  3.67  1.27  0.47  0.07  0.35  4.02  0.50  0.90  0.47 
2  4.49  2.09  0.37  0.05  0.90  5.39  0.60  0.84  0.54 
1  9.86  0.13  0.51  0.05  0.43  9.43  0.63  0.97  0.62 
We see that, just as for the nucleus Pb, Table I, there is a case of a catastrophic behavior with the factor close to zero, and again this is the state 2. From the above consideration for the lead chain, it is clear that for finding the and values we are interested mainly on the first particle and first hole states 2 and 1, correspondingly. Six states in Table VII create the model space in this case. In contrast to the lead case, where particle and hole states have mainly opposite parities, now all the members of are of positive parity, with the only exception of the state 1 which is rather distant and plays no important role in the PC corrections. This occurs because of a unique feature of the 1 level which alone is, in fact, a separate shell being an “intruder state” from the upper shell. Usually, an intruder state occurs among the states of the previous shell of states with opposite parity, the shell in this case. However, the latter is rather distant from the 1 level, and the inclusion of the shell into does not practically changes the PC effects we analyze. Such a peculiar position of the hole shell consisted from the single 1 state leads to another peculiarity of the PC corrections in this state. Usually, the pole correction , see other states in this table and the Table I for Pb, is rather big and negative, whereas the nonpole one is always positive but smaller than the pole correction in absolute value, so that the total PC correction remains to be negative. In this case the sign of the pole term is positive! Why it happens can be easily seen from Eq. (6). Indeed, all close states of the positive parity which are coupled with the phonon, are now on the other side of the Fermi level, the second term of this expression, which explains the positive sign of . As a consequence, the total PC correction to this SP level becomes positive.
nucl.  DF3a  DF3a +  exp mass ()  

Sn  2  0.295  0.172  0.58851(0.00924) 
Sn  2  0.527  1.110  1.47008(0.01065) 
Sn  2  1.338  1.901  2.28372(0.02263) 
Sn  2  2.142  2.838  3.05006(0.01777) 
Sn  2  2.936  3.484  3.73504(0.017) 
Sn  2  3.720  4.230  4.83153(0.00853) 
Sn  2  4.490  4.969  5.1103(0.0082) 
Sn  2  5.245  5.748  5.78895(0.00371) 
Sn  2  5.982  6.474  6.57225(0.00452) 
Sn  2  6.700  7.066  7.31099(0.00361) 
Sn  2  7.400  7.747  7.97316(0.01557) 
Sn  2  8.084  8.342  8.55629(0.03889) 
Sn  2  8.963  9.067  9.13802(0.00424) 
Sn  2  9.892  9.838  9.66757(0.00608) 
0.722  0.286 
Another conclusion, which can be made from the analysis of Table VII, is the main role of the phonon for PC corrections of the tin nuclei. Indeed, the influence of a phonon can be estimated from the value of the difference of . We see, that for all the states except 1, the role of the phonon is rather moderate. Therefore, we shall find the oddeven mass differences and for all the tin chain with account only of the phonon, sometimes checking the role of the phonon.
Table VIII contains the characteristics of the phonon which are used in the above calculation scheme. The excitation energies are in reasonable agreement with experimental data. As in the lead chain, we begin from the values.
nucl.  DF3a  DF3a +  exp mass ()  

Sn  1  4.786  4.264  5.00208(0.01534) 
Sn  1  5.652  5.014  5.79467(0.01651) 
Sn  1  6.537  5.918  6.64302(0.0178) 
Sn  1  7.424  6.669  7.55406(0.00408) 
Sn  1  8.288  7.682  8.48049(0.00182) 
Sn  1  9.102  8.544  9.27862(1.0809E4) 
Sn  1  9.859  9.320  9.99878(0.00538) 
Sn  1  10.585  10.024  10.68806(0.00821) 
Sn  1  11.298  10.743  11.39433(0.02982) 
Sn  1  12.007  11.578  12.09275(0.02084) 
Sn  1  12.713  12.291  12.82742(0.03747) 
Sn  1  13.420  13.083  13.75269(0.03885) 
Sn  1  14.128  13.905  14.58394(0.00482) 
Sn  1  14.842  14.696  15.8073(0.00557) 
0.324  0.740 
In the last line, just as in Tables V and VI, the average deviation is shown of the theoretical predictions from the experimental data, found from Eq. (30). As for the Pb chain, we compare the average accuracy of our calculations with those of the Skyrme EDF HFB17. In this case we have MeV, which is significantly better than the mean field prediction of the Fayans EDF DF3a. However, the account for the PC corrections makes the agreement for the latter significantly better, even better than that of the HFB17 EDF.
106  1.316  1.207  0.33 

108  1.231  1.206  0.37 
110  1.162  1.212  0.35 
112  1.130  1.257  0.40 
114  1.156  1.300  0.34 
116  1.186  1.294  0.32 
118  1.217  1.230  0.32 
120  1.240  1.171  0.33 
122  1.290  1.141  0.33 
124  1.350  1.132  0.27 
126  1.405  1.141  0.27 
128  1.485  1.169  0.23 
130  1.610    0.18 
132  4.329  4.041  0.22 
Let us go to the values which correspond to removing of a particle from an even Sn nucleus, i.e. to adding a hole to the state 1, see Table X. The above discussion of Table VII indicates that one could expect something wrong in this case. And a bad event does occur: the PC corrections spoil the agreement with the data, and significantly. This is especially disturbing as in this case the meanfield predictions of the DF3a EDF are rather good, on average better than those of the HFB17 EDF. Indeed, in this case the corresponding value of (HFB17)=0.484 is significantly worse than the DF3a value.
Our failure in the case of 1 holes confirms the above remarks that inclusion of the PC corrections from the lowlying phonons to the mean field observable values found for the EDFs which contain the phenomenological parameters involving such contributions on average is rather delicate procedure. In this case, there is a special reason why it occurs, as far as this intruder state alone generates a shell separated with the Fermi level from the particle shell of states of the same parity. Such a situation is essentially different from the “normal” one when two shells adjacent to the Fermi level, the particle and hole ones, contain mainly the states of opposite parities.
Iv Conclusion
In this article, a method, developed recently levsemi () to find the PC corrections to the SP energies of semimagic nuclei is used for finding the oddeven mass difference between the even members of a semimagic chain and their oddproton neighbors. In terms of TFFS AB (), these differences are the chemical potentials or for the addition or removal mode, correspondingly. The method of levsemi () is based on the direct solution of the Dyson equation with the PC corrected mass operator .
The Fayans EDF DF3a DF3a () is used for finding the initial values and for calculating the corresponding PC corrections. This is a small modification of the original Fayans EDF DF3 Fay4 (); Fay (), the spinorbit and effective tensor parameters being changed only. It should be stressed that, from the very beginning Fay1 (); Fay4 (); Fay (), the Fayans EDF was constructed in such a way that it should take into account all the PC effects on average. And, indeed, the DF3a EDF turned out to be rather accurate in predicting the values on the mean field level. We compare our results with predictions of the Skyrme EDF HFB17 HFB17 () fitted to nuclear masses with highest accuracy among the selfconsistent calculations. We analyzed four sets of data, and values in the Pb and Sn chains.
mode, chain  HFB17  DF3a  DF3a + PC 

, Pb  0.486  0.465  0.252 
, Pb  0.253  0.344  0.220 
, Sn  0.346  0.722  0.286 
, Sn  0.484  0.324  0.740 
Table XI contains the values of the average deviation of the theoretical predictions from the data for two EDFs, HFB17 and DF3a at the meanfield level and for the DF3a EDF with phonon corrections. We see that in two cases, the removal mode in Pb and the addition one in Sn, the EDF HFB17 exceeds DF3a in accuracy, but PC corrections to DF3a make the agreement better compared to the HFB17 EDF. In the two other cases, the addition mode in Pb and the removal one in Sn, the mean field predictions of the EDF DF3a turned out to be more accurate than those of HFB17. For the addition mode in Pb chain PC corrections make the accuracy of DF3a even better, but for the one in Sn chain a failure occurs: PC corrections spoil the agreement, and significantly.The reason is discussed in detail in the previous Section, but, in short, it is connected with a peculiar feature of the state which alone represents a separate shell. We think that such cases of “bad” behavior of the PC corrections will be seldom, but this case demonstrates, that it is not clear a priory if the PC corrections will improve agreement with experiment or not. This agrees qualitatively with the analysis of Dobaczewski () where it was found that the possibility of improvement of the meanfield description of the SP energies with account for the PC corrections depends significantly on the kind of the EDF used in these calculations.
Acknowledgements.
We acknowledge for support the Russian Science Foundation, Grants Nos. 161210155 and 161210161. The work was also partly supported by the RFBR Grant 160200228a. E.S. and S.P. thank the INFN, Sezione di Catania, for hospitality. This work has been carried out using computing resources of the federal collective usage center Complex for Simulation and Data Processing for Megascience Facilities 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
 (1) E. Litvinova and P. Ring, Phys. Rev. C 73, 044328 (2006).
 (2) LiGang Cao, G. Colò, H. Sagawa, and P.F. Bortignon, Phys. Rev. C 89, 044314 (2014).
 (3) D. Tarpanov, J. Dobaczewski, J. Toivanen, and B.G. Carlsson, Phys. Rev. Lett. 113, 252501 (2014).
 (4) M. Baldo, P.F. Bortignon, G. Coló, D. Rizzo, and L. Sciacchitano, J. Phys. G: Nucl. Phys. 42, 085109 (2015).
 (5) N.V. Gnezdilov, I.N. Borzov, E.E. Saperstein, and S.V. Tolokonnikov, Phys. Rev. C 89, 034304 (2014).
 (6) V. A. Khodel, E. E. Saperstein, Phys. Rep. 92, 183 (1982).
 (7) V. A. Khodel’, Sov. J. Nucl. Phys. 24, 282 (1976).
 (8) S. Kamerdzhiev and E. E. Saperstein, EPJA 37, 333 (2008).
 (9) S. Weinberg, Phys. Rev. Lett. 31, 494 (1973).
 (10) H. Grawe, K. Langanke, and G. MartínezPinedo, Rep. Prog. Phys. 70, 1525 (2007).
 (11) E.E. Saperstein, M. Baldo, S.S. Pankratov, and S.V. Tolokonnikov, JETP Lett. 104, 609 (2016).
 (12) P. Ring and E. Werner, Nucl. Phys. A 211, 198 (1973).
 (13) A.B. Migdal Theory of finite Fermi systems and applications to atomic nuclei (Wiley, New York, 1967).
 (14) E.E. Saperstein, M. Baldo, S.S. Pankratov, and S.V. Tolokonnikov, JETP Lett., to be published, 106, (2017).
 (15) A.V. Smirnov, S.V. Tolokonnikov, S.A. Fayans, Sov. J. Nucl. Phys. 48, 995 (1988).
 (16) I. N. Borzov, S. A. Fayans, E. Kromer, and D. Zawischa, Z. Phys. A 355, 117 (1996).
 (17) S.A. Fayans, S.V. Tolokonnikov, E.L. Trykov, and D. Zawischa, Nucl. Phys. A 676, 49 (2000).
 (18) S.V. Tolokonnikov and E.E. Saperstein, Phys. At. Nucl. 73, 1684 (2010).
 (19) A. Bohr and B.R. Mottelson, Nuclear Structure (Benjamin, New York, Amsterdam, 1969.), Vol. 1.
 (20) S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. Lett. 102, 152503 (2009).

(21)
S. Goriely, http://wwwastro.ulb.ac.be/bruslib/
nucdata/  (22) S.V. Tolokonnikov, S. Kamerdzhiev, D. Voytenkov, S. Krewald, and E.E. Saperstein, Phys. Rev. C 84, 064324 (2011).
 (23) S.V. Tolokonnikov, S. Kamerdzhiev, S. Krewald, E.E. Saperstein, and D. Voitenkov, Eur. Phys. J. A 48, 70 (2012).
 (24) M. Baldo, P.F. Bortignon, G. Coló, D. Rizzo, and L. Sciacchitano, J. Phys. G: Nucl. Phys. 42, 085109 (2015).
 (25) S. Raman, C.W. Nestor Jr., and P. Tikkanen, At. Data Nucl. Data Tables 78, 1 (2001).
 (26) A. Bohr and B.R. Mottelson, Nuclear Structure (Benjamin, New York, 1974.), Vol. 2.
 (27) M. Wang, G. Audi, A.H. Wapstra, F.G. Kondev, M. MacCormick, X. Xu, and B. Pfeiffer, Chinese Physics C, 36, 1603 (2012).