# Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.

###### Abstract

A method is presented to evaluate the particle-phonon coupling (PC) corrections to the single-particle energies (SPEs) in semi-magic nuclei. In such nuclei always there is a collective low-lying phonon, and a strong mixture of single-particle and particle-phonon states often occurs. As in magic nuclei, the so-called approximation, where is the vertex of the -phonon creation, can be used for finding the PC correction to the initial mass operator . In addition to the usual pole diagram, the phonon “tadpole” diagram is also taken into account. In semi-magic nuclei, the perturbation theory in with respect to is often invalid for finding the PC corrected SPEs. Instead, the Dyson equation with the mass operator is solved directly, without any use of the perturbation theory. Results for a chain of semi-magic Pb isotopes are presented.

Last decade, there was a revival of the interest within different self-consistent nuclear approaches to study the particle-phonon coupling (PC) effects in the single-particle energies (SPEs) of magic nuclei. We cite here such studies within the relativistic mean-field theory [1], within the Skyrme–Hartree–Fock method [2, 3, 4] and on the basis of the self-consistent theory of finite Fermi systems (TFFS) [5]. The Fayans energy density functional (EDF) was used in the last case to find the self-consistent basis. In all the references cited above double-magic nuclei were considered. There are several reasons for such choice. First, these nuclei are non-superfluid which simplifies the theoretical analysis. Second, the so-called approximation is, as a rule, valid in magic nuclei, being the vertex of the -phonon creation. Moreover, the perturbation theory in terms of the PC correction to the mass operator ia also applicable, which makes evaluation of PC corrected SPEs rather simple. At last, there is a lot of experimental data on SPEs in these nuclei [6].

In this work, we extend the field of this problem to semi-magic nuclei. Unfortunately, the experimental data on the SPEs in semi-magic nuclei are rather scarce. Indeed, the single-particle spectroscopic factor of an excited state under consideration should be known. In addition, its value should be rather large, in order that one can interpret this state as a single-particle one. Extraction of the spectroscopic factors from the reaction data is a complicated theoretical problem, therefore the list of known spectroscopic factors is rather limited. However, the PC corrections to the SPEs are necessary not only by themselves, they are also usually important ingredients of the procedure of finding PC corrections to other nuclear characteristics, e.g., magnetic moments and transitions in odd nuclei [7, 8, 9]. PC corrections to the double odd-even mass differences found in the approach starting from the free potential [10, 11] is another example where the PC corrections to SPEs are of primary importance.

A semi-magic nucleus consists of two sub-systems with different properties. One of them, magic, is normal, whereas the non-magic counterpart is superfluid. We will consider the SPEs of the normal sub-system only. Therefore, the main part of the formalism for description of the PC corrections developed for magic nuclei [5] remains valid. One difference with double-magic nuclei is that the vertex obeys the QRPA-like TFFS equation in superfluid nuclei, in contrast to simple RPA-like equation in magic nuclei. This complication is not very serious as we developed the necessary method with the use of the Fayans EDF in a previous paper [12]. A real difficulty arises in non-magic nuclei due to appearance of low-lying states which is a characteristic feature of such nuclei. As a result, small denominators appear regularly in the expressions for the PC corrections which makes unapplicable a plane perturbation theory.

To find the SPEs with account for the PC effects, we solve the following equation:

(0) |

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 one-particle Green function with the mass operator .

In magic nuclei [5, 10, 11] the perturbation theory in with respect to was used to solve this equation:

(0) |

where

(0) |

In this article, we will solve Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) directly, without any additional approximations. As to the -approximation for the PC correction to the mass operator, it remains valid in semi-magic nuclei, and only the next step from (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) to (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) becomes invalid. In the case when several -phonons are taken into account, the total PC variation of the mass operator in Eqs. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.)–(Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) is just the sum:

(0) |

The diagrams for the operator within the -approximation are displayed in Fig. 1. The first one is the usual pole diagram, with obvious notation, whereas the second, “tadpole” diagram represents the sum of all non-pole diagrams of order .

Explicit expression for the pole term is well known [1, 5, 10], and we present it just for completeness:

(0) | |||||

(0) |

where is the excitation energy of the -phonon and are the particle occupation number (remind that we deal with the normal subsystem of a semi-magic nucleus).

The vertex obeys the QRPA-like TFFS equation [13],

(0) |

where all the terms are matrices. The angular momentum projection , which is written down in Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) explicitly, is here and below for brevity omitted. In the standard TFFS notation, we have:

(0) |

(0) |

In (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.), is the normal component of the vertex , whereas are two anomalous ones. In Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.), is the usual Landau–Migdal 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 functios and [13].

As we need the proton vertex and the proton subsystem is normal, only the normal vertex is non-zero in this case. This is explicit meaning of the short notation in (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) and below.

For solving the above equations, we use the self-consistent basis generated by the version DF3-a [14] of the Fayans EDF [15, 16]. The nuclear mean-field potential is the first derivative of over .

All low-lying phonons we deal with are of surface nature, the surface peak dominating in their creation amplitude:

(0) |

The first term in this expression is surface peaked, whereas the in-volume term is rather small. It is illustrated in Fig. 2 for the and states in Pb. In this work, just as in [5, 10, 11], we neglect the in-volume term in (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) when considering the tadpole PC term of . In the result, it is reduced, see Ref. [5], to rather simple form:

(0) |

Note that the above scheme for the ghost () phonon results in an explicit expression for the “recoil effect.” Details can be found in [5].

In this work, we limit ourselves with four even lead isotopes, Pb. In all cases we consider two low-lying phonons, and . Their excitation energies are presented in Table 1. As one can see, they agree with existing experimental data sufficiently well. The ghost is also taken into account, although the corresponding correction for nuclei under consideration is very small, because it depends on the mass number as .

nucleus | ||||
---|---|---|---|---|

Pb | 0.789 | 1.026 | 2.620 | - |

Pb | 0.823 | 0.960 | 2.704 | 2.517 |

Pb | 0.882 | 0.899 | 2.785 | 2.621 |

Pb | 0.945 | 0.803 | 2.839 | 2.648 |

Let us go to the description of the method we use to solve Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.). As the PC corrections are important only for the SPEs 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. To avoid any misunderstanding, we stress that this restriction concerns only Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.). In Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) for , we use rather wide single-particle space with energies MeV. We take as an example to illustrate the method the nucleus Pb which is sufficiently distant from the double-magic Pb to be considered as a typical semi-magic. 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 (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.), which simplifies very much the solution of the Dyson equation. In the result, Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) reduces as follows:

(0) |

The tadpole term does not depend on the energy, therefore only poles of Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) are the singular points of Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.). They can be readily found from (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) in terms of and . It can be easily seen that the lhs of Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) 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 single-particle state from which the state originated. The latter is a mixture of a single-particle state with several particle-phonon states. The corresponding single-particle strength distribution factors (-factors) can be found similar to (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.):

(0) |

Evidently, they should obey the normalization rule:

(0) |

Accuracy of fulfillment of this relation is a measure of validity of the model space we use to solve the problem under consideration. In the major part of the cases we will consider it is fulfilled with 1–2% accuracy.

A set of solutions for four states in Pb is presented in Table 2. 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” single-particle states. In such cases, the following prescription looks natural for the PC corrected single-particle characteristics:

(0) |

This is an analog of Eqs. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) and (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.) in the perturbative solution.

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 | |

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 | |

3 | 1 | -9.877 | 0.608 |

2 | -8.536 | 0.604 | |

3 | -6.493 | 0.839 | |

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 | |

Nucleus | ||||||
---|---|---|---|---|---|---|

Pb | 1i | -0.26 | 0.39 | 0.13 | -0.13 | 0.96 |

2f | -1.05 | 0.24 | -0.30 | -1.35 | 0.83 | |

1h | -2.33 | 0.33 | 0.12 | -2.21 | 0.93 | |

3s | -5.81 | 0.20 | 0.01 | -5.80 | 0.89 | |

2d | -6.67 | 0.21 | 0.17 | -6.50 | 0.89 | |

1h | -7.06 | 0.37 | 0.25 | -6.81 | 0.93 | |

2d | -7.88 | 0.21 | 0.28 | -7.60 | 0.88 | |

1g | -9.97 | 0.29 | 0.08 | -9.89 | 0.91 | |

Pb | 1i | -0.74 | 0.41 | 0.13 | -0.61 | 0.95 |

2f | -1.52 | 0.25 | -0.29 | -1.81 | 0.83 | |

1h | -2.86 | 0.34 | 0.13 | -2.73 | 0.93 | |

3s | -6.26 | 0.21 | 0.01 | -6.25 | 0.89 | |

2d | -7.09 | 0.22 | 0.17 | -6.92 | 0.89 | |

1h | -7.52 | 0.38 | 0.25 | -7.27 | 0.93 | |

2d | -8.34 | 0.22 | 0.30 | -8.04 | 0.87 | |

1g | -10.46 | 0.30 | 0.08 | -10.38 | 0.91 | |

Pb | 1i | -1.21 | 0.32 | 0.14 | -1.07 | 0.97 |

2f | -2.01 | 0.20 | -0.23 | -2.24 | 0.87 | |

1h | -3.36 | 0.27 | 0.17 | -3.19 | 0.96 | |

3s | -6.72 | 0.17 | 0.23 | -6.49 | 0.84 | |

2d | -7.51 | 0.17 | 0.05 | -7.46 | 0.94 | |

1h | -7.98 | 0.30 | 0.25 | -7.73 | 0.95 | |

2d | -8.80 | 0.18 | 0.17 | -8.63 | 0.92 | |

1g | -10.93 | 0.24 | 0.13 | -10.80 | 0.95 | |

Pb | 1i | -1.67 | 0.30 | 0.07 | -1.60 | 0.89 |

2f | -2.51 | 0.19 | -0.30 | -2.81 | 0.82 | |

1h | -3.82 | 0.25 | 0.19 | -3.63 | 0.97 | |

3s | -7.18 | 0.16 | 0.08 | -7.10 | 0.89 | |

2d | -7.91 | 0.16 | 0.11 | -7.80 | 0.86 | |

1h | -8.42 | 0.27 | 0.37 | -8.05 | 0.88 | |

2d | -9.28 | 0.16 | 0.12 | -9.16 | 0.95 | |

1g | -11.36 | 0.22 | 0.24 | -11.12 | 0.90 |

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 . In such cases, we suggest the following generalization of Eq. (Phonon-particle coupling effects in single-particle energies of semi-magic nuclei.):

(0) |

where

(0) |

In both the above sums, only the states with appreciable values of are included. In practice, we include in these sums the states with . The value of is just the centroid of the single particle energy distribution.

The results for SPEs and -factors are presented in Table 3. For each segment of the table which relates to the even core Pb, the lower part, till the state, describes the proton holes, i. e. the states of the Tl nucleus. On the contrary, the upper states are the proton particle states, i.e. they belong to the Bi nucleus. The tadpole term is given separately. The recoil correction is not presented explicitly as it is very small in the lead region, the maximum value is MeV. Thus, we may approximately include it to the pole correction . We see that the tadpole correction is always positive, whereas the pole one is, as a rule, negative. Their absolute values are of comparable magnitude. What is more, the tadpole term often dominates. In such cases, neglect with the tadpole term leads to a non-correct sign of the total SPE correction. In other cases, calculations of the PC corrections to SPEs without the tadpole term lead to correct signs of but overestimate their absolute values significantly. Thus, for the semi-magic nuclei we see the tendency noted before in magic nuclei [5].

In Fig. 4, we display the SPEs corresponding to the ground states of the corresponding nuclei. In this case, the experimental values are known as they can be found in terms of mass values of neighboring odd and even nuclei [17]. Explicitly, they are equal to , the proton separation energies taken with the opposite sign. The results are compared of two sets of calculations. In the first case, the single phonon is taken into account, whereas both the phonons, and , participate in the second set of calculations. We see, firstly, that the role of both phonons is comparable. Secondly, the total PC correction is rather small, and, finally, often it makes the agreement with the data worse. This is not strange as we deal with the EDF method containing phenomenological parameters fitted mainly to nuclear masses, which determine the experimental SPEs under discussion.

In general, the problem of explicit consideration of the PC corrections within the EDF method or any other self-consistent approach operating with phenomenological parameters is rather delicate. Indeed, these parameters include different PC effects implicitly. Therefore, a regular inclusion of the PC corrections inevitably should be accompanied with a readjustment of the initial parameters. In such a situation, it is a more promising method to separate the fluctuating part of the PC corrections, which changes in a non-regular way from a nucleus to another. Such a strategy was chosen, e.g., in [18] to explain an anomalous dependence of charge radii of heavy calcium isotopes found recently by the ISOLDE collaboration [19].

A different situation occurs typically for the approaches which start from the free interactions. In these cases, the PC induced corrections should be just added to the main terms found with a free potential, and the PC correction to the SPEs is one of the necessary ingredients for such calculations. Such approach was rather popular in the last decade for the nuclear pairing problem [20, 21, 22, 23]. An analogous method was developed in [24, 25] to find the double odd-even mass differences of magic nuclei or semi-magic ones, for the normal subsystems in the last case. In this problem, a method to find the PC corrections was developed for magic nuclei in [10, 11]. In that case, a plain perturbation theory was used to find the PC corrections to the SPEs.

To resume, a method is developed to find the PC corrections to SPEs for semi-magic nuclei beyond the perturbation theory in the PC correction to the mass operator with respect to . Instead, the Dyson equation with the mass operator is solved directly, without any use of the perturbation theory. The method is checked for a chain of even Pb isotopes. This makes it possible to extend to semi-magic nuclei the field of consistent consideration of the PC corrections to the double odd-even mass differences and some another problems. For semi-magic nuclei under consideration, the tadpole correction to the SPEs turned out to be of primary importance.

We acknowledge for support the Russian Science Foundation, Grants Nos. 16-12-10155 and 16-12-10161. The work was also partly supported by the Russian Foundation for Basic Research Grants 14-02-00107-a, 14-22-03040-ofi_m and 16-02-00228-a. Calculations were partially carried out on the Computer Center of Kurchatov Institute.

- 1. E. Litvinova and P. Ring, Phys. Rev. C 73, 044328 (2006).
- 2. Li-Gang 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. H. Grawe, K. Langanke, and G. Martínez-Pinedo, Rep. Prog. Phys. 70, 1525 (2007).
- 7. E.E. Saperstein, S. Kamerdzhiev, S. Krewald, J. Speth, and S.V. Tolokonnikov, EPL, 103, 42001 (2013).
- 8. E.E. Saperstein, S.P. Kamerdzhiev, S. Krewald, J. Speth, and S.V. Tolokonnikov, JETP Lett 98, 562 (2013).
- 9. E.E. Saperstein, O.I. Achakovskiy, S. Kamerdzhiev, S. Krewald, J. Speth, and S.V. Tolokonnikov. Phys. At. Nucl. 77, 1033 (2014).
- 10. E.E. Saperstein, M. Baldo, N.V. Gnezdilov, S.V. Tolokonnikov, JETP Lett., 103, 3 (2016).
- 11. E.E. Saperstein, M. Baldo, N.V. Gnezdilov, S.V. Tolokonnikov, Phys. Rev. C 93, 034302 (2016).
- 12. S.V. Tolokonnikov, S. Kamerdzhiev, D. Voytenkov, S. Krewald, and E. E. Saperstein, Phys. Rev. C 84, 064324 (2011).
- 13. A.B. Migdal Theory of finite Fermi systems and applications to atomic nuclei (Wiley, New York, 1967).
- 14. S.V. Tolokonnikov and E.E. Saperstein, Phys. At. Nucl. 73, 1684 (2010).
- 15. A.V. Smirnov, S.V. Tolokonnikov, S.A. Fayans, Sov. J. Nucl. Phys. 48, 995 (1988).
- 16. S.A. Fayans, S.V. Tolokonnikov, E.L. Trykov, and D. Zawischa, Nucl. Phys. A 676, 49 (2000).
- 17. M. Wang, G. Audi, A.H. Wapstra, F.G. Kondev, M. MacCormick, X. Xu, and B. Pfeiffer, Chinese Physics C, 36, 1603 (2012).
- 18. E.E. Saperstein, i.n. Borzov, and S.V. Tolokonnikov, JETP Lett., 104, 216 (2016).
- 19. R.F. Garcia Ruiz, M.L. Bissell, K. Blaum et al. (Collaboration), Nature Physics 12, 594 (2016).
- 20. F. Barranco, R.A. Broglia, G. Colo, G. Gori, E. Vigezzi, and P.F. Bortignon, Eur. Phys. J. A 21, 57 (2004).
- 21. A. Pastore, F. Barranco, R.A. Broglia, and E. Vigezzi, Phys. Rev. C 78, 024315 (2008).
- 22. S.S. Pankratov, M. Baldo, M.V. Zverev, U. Lombardo, E.E. Saperstein, and S.V. Tolokonnikov, JETP Lett. 90, 612 (2009).
- 23. S.S. Pankratov, M.V. Zverev, M. Baldo, U. Lombardo, and E.E. Saperstein, Phys. Rev. C 84, 014321 (2011).
- 24. N.V. Gnezdilov and E.E. Saperstein, JETP Lett. 95, 603 (2012).
- 25. E.E. Saperstein, M. Baldo, N.V. Gnezdilov, U. Lombardo, and S.S. Pankratov, EPJ Web. Conf. 38, 05002 (2012).