# Random Matrix Analysis of the Monopole Strength Distribution in Pb

###### Abstract

We study statistical properties of the spectrum of Pb in the energy region MeV. We use the Skyrme interaction SLy4 as our model Hamiltonian to create a single-particle spectrum and to analyze excited states. The finite-rank separable approximation for the particle-hole interaction enables us to perform the calculations in large configuration spaces. We show that while the position of the monopole resonance centroid is determined by one phonon excitations of , the phonon-phonon coupling is crucial for the description of a strength distribution of the spectrum. In fact, this coupling has an impact on the spectral rigidity which is shifted towards the random matrix limit of the Gaussian orthogonal ensembles.

###### pacs:

24.60.Lz, 21.60.Jz, 27.80.+w## I Introduction

Nuclear Giant Resonances (GR) are the subject of numerous investigations over several decades dstadta (). Some of the basic features such as centroids and collectivity (in terms of the sum rules) are reasonably well understood within microscopic theories BM (); M (); S (). As yet we have no answer on the question how a collective mode like the GR dissipates its energy.

According to accepted wisdom, GRs are essentially excited by an external field through a one-body interaction. It is, therefore, natural to describe these states as collective 1p-1h states. Once excited, the GR progresses to a fully equilibrated system via direct particle emission and by coupling to more complicated configurations (2p-2h, 3p-3h, etc). The former mechanism gives rise to an escape width. It is expected that the decay evolution along the hierarchy of more complex configurations till compound states determines spreading widths. A full description of this decay represents a fundamental problem which is, however, difficult to solve (if even is possible at all ?) due to existence of many degrees of freedom for many-body quantum system such as a nucleus. Therefore, to gain an insight into the nature of GR spreading we have to introduce various approximations or a model, which validity depends on a primal feasibility.

In general, the description of GR properties is based on the analysis of the electromagnetic strength distribution in an energy interval which is large enough to catch hold of basic GR features that is under investigation. An obvious requirement to the model consideration is to use configurations with various degree of complexity. Evidently, the more complex configuration is considered, the more cumbersome task should be solved. Therefore, a natural question arises: what type of a particular complex configuration should be enough in order to understand a gross structure of a particular GR which data are available in modern experiments? In addition, once this complex configuration is defined one can further ask about statistical properties of states that compose a GR strength distribution. As a result, one could illuminate the role of various correlations that diminish an importance of a specific excitation that determines a centroid position of a specific GR.

To answer these questions we will employ two approaches. On the one hand, the random matrix theory (RMT) B81 (); M91 (); GMW (); HW (); Gom () provides necessary tools to shed light on the spectral properties and the distribution of transition-strength properties, when specific features become not of a primary importance. The RMT assumes only that the nuclear Hamiltonian belongs to an ensemble of random matrices that are consistent with the fundamental symmetries of the system. In particular, in the case of the time-reversal symmetry, the relevant ensemble is the Gaussian orthogonal ensemble (GOE). On the other hand, to understand the fragmentation of high-lying states it is necessary also to exploit nuclear structure models that are based on the microscopic many-body theory, where the effects of the residual interaction on the statistics must be studied in large model spaces. Introducing a residual interaction in general implies a transition to GOE-properties above some excitation energy A ().

The quasiparticle-phonon model (QPM) S () offers an attractive framework for such studies. The separable form of the residual interaction of a model Hamiltonian allows to diagonalize it in a space spanned by states composed of one, two and three phonons considered in the random phase approximation (RPA). We would like to mention here the RMT analysis of statistical properties of a pygmy dipole resonance within the QPM, based on the Woods-Saxon potential VP (). It is desirable, however, to use a unified approach in which a mean field and a residual interaction are treated at the same footing in order to avoid any artifacts sfkh78 (); BNY (). For our purposes we choose the modern development of the QPM, a finite rank separable approximation (FRSA) gsv98 (); ssvg02 (); svg08 (). The FRSA follows the basic QPM ideas, but the single-particle (sp) spectrum and the residual interaction are calculated with the Skyrme forces. This approach enables us to consider a coupling between the one- and two-phonon components of the wave functions svg04 (). It was successfully used to study the properties of the low-lying states and giant resonances within the RPA and beyond gsv98 (); ssvg02 (); svg08 (); svg04 (); sap12 (); sapw14 ().

By means of this approach and by the RMT tools we attempt in this paper to understand the complex structure observed in the spectrum of the doubly magic nucleus Pb in the region of the isoscalar giant monopole resonance (ISGMR). This strength distribution is extensively studied in many experiments Y99 (); Y04 (); P13 (); P14 (). The experimental properties have been described within the RPA with the Skyrme interactions (for a review see, for example, Ref.rev2007 ()). In this system the effect of the anharmonicity is expected to be small. Contrary to the expectations, we will show the importance of the phonon-phonon coupling (PPC) effects for the statistical properties of the spectrum calculated with the aid of the Skyrme SLy4 interaction, taken as an example.

## Ii The Model

For the analysis of the doubly magic nucleus we impose a spherical symmetry on the sp wave functions in our HF calculations. The continuous part of the sp spectrum is discretized by diagonalizing the HF Hamiltonian on a harmonic oscillator basis. The cut-off of the continuous part is at the energy of 100 MeV. As the parameter set, we use the Skyrme force SLy4 sly () which was adjusted to reproduce the nuclear matter properties, as well as nuclear charge radii, binding energies of doubly magic nuclei. The residual particle-hole interaction is obtained as the second derivative of the energy density functional with respect to the particle density. By means of the standard procedure tebdns05 () we obtain the familiar RPA equations in the 1p-1h configuration space. The eigenvalues of the RPA equations are found numerically as the roots of a relatively simple secular equation within the FRSA gsv98 (). Since the FRSA enables to us to use the large 1p-1h space, there is not need in effective charges.

Using the basic QPM ideas in the simplest case of the configuration mixing S (), we construct the wave functions from a linear combination of one-phonon and two-phonon configurations states as

(1) |

where is the RPA excitation having energy ; denotes the total angular momentum and is its z-projection in the laboratory system. The ground state is the RPA phonon vacuum . The normalization condition for the wave functions (1) yields

(2) |

The variational principle leads to a set of linear equations for the unknown amplitudes and svg04 ():

(3) |

(4) |

The rank of the set of linear equations is equal to the number of one- and two-phonon configurations included in the wave function (1). To resolve this set it is required to compute the coupling matrix elements

(5) |

between one- and two-phonon configurations (see details in Ref. svg04 ()). Evidently, the nonzero matrix elements result in the inclusion of the PPC effects. Eqs.(3) and (4) have the same form as the QPM equations S (). It is important to stress, however, that the sp spectrum and the parameters of the residual p-h interaction are calculated with the chosen Skyrme forces, without any further adjustments.

The excitation operator of the ISGMR is defined as

(6) |

The wave functions (1) allow us to determine the transition probabilities . The matrix elements for direct excitation of two-phonon components from the ground state are about two orders of magnitude smaller as compared to the excitation of one-phonon components S (). Therefore, they are neglected in our calculation of the transition probabilities. The RPA analysis of the ISGMR shows that 96 % of the non energy-weighted sum rules (NEWSR) is located in the energy region MeV. To build the wave functions (1) of the excited states up to 20 MeV we take into account all one- and two-phonon configurations that are constructed from the , , , and phonons with energies below 25 MeV for computational convenience. The high-energy configurations plays a minor role in our calculations. It is noteworthy that the pair-transfer mode (see, e.g., Refs. J87 (); H15 ()) is outside the present work.

Properties of the low-energy two-phonon states are reflected in the deviation from the harmonic picture for the multiphonon excitations L97 (); P99 (). It is interesting to study the energies, reduced transition probabilities of the , and states which are the important ingredients of our calculations of the two-phonon states of Pb. The results obtained within the one-phonon approximation are compared with the experimental data H82 (); S83 () in Table 1. There is a satisfactory description of the reduced transition probabilities. The overestimate of the experimental energies indicates on some missing mechanisms. In our consideration we consider the PPC as the one that might improve the description.

The strength distribution of ISGMR is displayed in Fig.1. Both, experimental P13 () and theoretical results show the fragmentation and splitting of the ISGMR strength. The coupling between the one- and two-phonon states yields a noticeable redistribution of the ISGMR strength in comparision with the RPA results. In particular, the coupling decreases the NEWSR till 78 % in the ISGMR region ( MeV). Also, the PPC induces the 1 MeV downward shift of the main peak. There are the low-energy part, the main peak and the high-energy tail. The coupling produces a shift of order 11 % (7 %) of the NEWSR from the ISGMR region to the high (lower) energy region. The strength distribution of the ISGMR obtained within the PPC is rather close to the experimental distribution P13 (). Our analysis shows that the major contribution to the strength distribution is brought about by the coupling between the and components. We recall that the importance of the complex configurations for the interpretation of basic peculiarities of the ISGMR strength distribution of Pb was already qualitatively discussed in the framework of simple model GB79 (); BB81 (). Our calculations give the same tendency.

We turn now to the mechanism that dominates in the low-energy part of the spectrum. There is no state below 10.2 MeV. The extension of the variational space from the standard RPA to two-phonon configurations result in a formation of the low-lying states. The state is the lowest excitation which leads to the minimal two-phonon energies and the maximal matrix elements coupling between one- and two-phonon configurations. Since the PPC induces a downward shift of the energy, the energy state at 6.5 MeV is very close to the value with (see Table 1). Our analysis suggests the dominance (85%) of the octupole , and configurations in the wave functions of the excited , and states, respectively. The collective character of the state is mainly due to their coupling to the ISGMR, produced by the state. In particular, the wave function normalization of the state contains 4% of the . This small change in structure has a large impact on the value, see Fig.1. The lowest two-phonon state was first observed as the lowest-spin member of the multiplet in Ref. Y96 (). This fact was confirmed by the QPM analysis P99 ().

## Iii Spectral statistics

Let us study statistical properties of the spectrum up to 20 MeV. We examine the spectra calculated with and without the PPC effects, i.e., the cases of and , respectively. Fig. 2 displays the PPC impact on the energies. As was mentioned above, the coupling shifts down the part of states and modifies the level density. To elucidate the role of the residual interaction we also consider the spectrum of unperturbed 1p-1h and 2p-2h states (see Fig. 2a). Note that the level density of the unperturbed 3p-3h states is much smaller than the 2p-2h ones. As can be seen from Fig. 2, the difference between the unperturbed p-h and the spectra is remarkable. The downward shift of the spectrum is due to the residual interaction in the RPA framework. The coupling does not lead to visible spectrum changes. However, it brings important correlations that affect the spectral statistics.

The three spectra are analyzed within the RMT that enables to us to study the statistical laws governing fluctuations that, in general, can have very different origins. Starting from the spectrum , one can construct the staircase function which is defined as the state number below the energy . The function can be separated in a smooth part and the fluctuating part , where the integral of is zero. The function can be determined either from semiclassical arguments or using a polynomial for . To get a constant mean spacing of levels, we employ the unfolded spectrum defined by the mapping . Note that the values are introduced as the spacings. We use two typical measures to analyze the fluctuation properties of unfolded spectrum: the nearest-neighbor spacing distribution (NNSD) and the spectral rigidity of Dyson and Metha, the statistics DM63 ().

Due to the unfolding we have

(7) |

If the unfolded energies are in a regular system then the NNSD is known as the Poisson distribution,

(8) |

In the GOE, i.e. the energies are in a chaotic system, the NNSD is approximately given as the Wigner distribution M91 (),

(9) |

As can be seen from Fig. 3, for the unperturbed p-h spectrum we obtain a behavior close to the Poisson distribution, expected for uncorrelated energies. For the case of , the statistics is changing to the GOE limit. This fact indicates the onset of correlations that redistribute a strength over two-phonon components constructed by phonons with the other multipolarities. Indeed, the comparison of the NNSD without and with the coupling illuminates this fact evidently. At the spectrum is characterized by the Poisson (uncorrelated) statistics. The coupling between the one- and two-phonon components modifies the spectrum, and the NNSD becomes close to the Wigner surmise.

Another measure of correlations is the statistics defined as

(10) |

It characterizes the deviation of the staircase function from a straight line, and a rigid unfolded spectrum corresponds to smaller values of , while a soft spectrum has a larger . In fact, for a given L, smaller values of imply stronger long-range correlations between the levels.

For the sake of convenience, the function , averaged over intervals

(11) |

can be easily calculated from the number statistics, , which is the number of levels in an energy interval of length

(12) |

(13) |

For an uncorrelated spectrum one has

(14) |

while for the GOE it is

(15) |

in the limit. Fig.4 demonstrates the evolution of the measure from the uncorrelated states to the GOE limit, when the PPC effects are only responsible for the statistical correlations. These correlations dissolve the collective ISGMR in the sea of the fragmented two-phonon components created by the other multipolarities.

## Iv Summary

With the aid of a finite rank separable approximation we have analyzed the strength distribution of states ( MeV) of Pb. To simulate the mean-field we have used the SLy4 Skyrme interaction. To analyse the excitations we take into account all RPA states with = , , , , . Within the RPA approach the centroid location of the ISGMR is found at MeV. On the other hand, we have demonstrated that the coupling between one- and two-phonon terms in the wave functions of excited states is crucially important for the interpretation of the strength distribution of the ISGMR in the energy interval MeV. The results of the calculated transition-strength distribution are generally in a reasonable agreement with the experimental data.

The RMT measures such as the NNSD and the function indicate a transition towards GOE as soon as the coupling is switched on. It appears that the presence of two-phonon components in our wave function, in addition to the one-phonon ones, already enables to us to describe the gross strength distribution of the ISGMR in the experimentally available energy interval. We observed that the major contribution that evolves the system under consideration to the GOE limit is brought about by two-pnonon components of the octupole nature. In fact, their number exceeds essentially the numbers of two-phonon components that are constructed from phonons with the other multipolarities. A further systematic statistical studies of the impact of the phonon-phonon coupling on the vibrational spectra and the -transition strengths is clearly necessary and is in progress.

## Acknowledgments

A.P.S. thanks for the hospitality the Division of Mathematical Physics, Lund University, where a part of this work has been done. This work was partly supported by Russian Foundation for Basic Research under Grant nos. 16-52-150003 and 16-02-00228.

## References

- (1) M. N. Harakeh and A. van der Woude, Giant Resonances: Fundamental High-Frequency Modes of Nuclear Excitation (Clarendon Press, Oxford, 2001).
- (2) A. Bohr and B. Mottelson, Nuclear Structure, Vol. 2 (World Scientific, Singapore, 1998).
- (3) A. B. Migdal, Teoriya Konechnykh Fermi-Sistem i Svoistva Atomnih Yader (2nd edn., Nauka, Moscow, 1983).
- (4) V. G. Soloviev, Theory of Atomic Nuclei: Quasiparticles and Phonons (Institute of Physics, Bristol and Philadelphia, 1992).
- (5) T. A. Brody, J. Flores, J. B. French, P. A. Mello, A. Pandey, and S. S. M. Wong, Rev. Mod. Phys. 53, 385 (1981).
- (6) M. L. Metha, Random Matrices (2nd edn., Academic, New York, 1991).
- (7) T. Guhr, A. Müller-Groeling, and H. A. Weidenmüller, Phys. Rep. 299, 189 (1998).
- (8) H. A. Weidenmüller and G. E. Mitchell, Rev. Mod. Phys. 81, 539 (2009).
- (9) J.M.G. Gómez, K. Kar, V.K.B. Kota, R.A. Molina, A. Relaño, and J. Retamosa, Phys. Rep. 499, 103 (2011).
- (10) S. Åberg, Phys. Rev. Lett., 64, 3119 (1990).
- (11) J. Enders, T. Guhr, A. Heine, P. von Neumann-Cosel, V.Yu. Ponomarev, A. Richter, and J. Wambach, Nucl. Phys. A 741, 3 (2004).
- (12) E. E. Saperstein, S. A. Fayans, and V. A. Khodel, Fiz. Elem. Chastits At. Yadra 9, 221 (1978) [Sov. J. Part. Nucl. 9, 91 (1978)].
- (13) see discussion on self-consistent effects in various mesoscopic systems, for example, in J.L. Birman, R.G. Nazmitdinov, V.I. Yukalov, Phys. Rep. 526, 1 (2013).
- (14) Nguyen Van Giai, Ch. Stoyanov, and V. V. Voronov, Phys. Rev. C 57,1204 (1998).
- (15) A. P. Severyukhin, Ch. Stoyanov, V. V. Voronov, and Nguyen Van Giai, Phys. Rev. C 66, 034304 (2002).
- (16) A. P. Severyukhin, V. V. Voronov, and Nguyen Van Giai, Phys. Rev. C 77, 024322 (2008).
- (17) A. P. Severyukhin, V. V. Voronov, and Nguyen Van Giai, Eur. Phys. J. A22, 397 (2004).
- (18) A. P. Severyukhin, N. N. Arsenyev, and N. Pietralla, Phys. Rev. C 86, 024311 (2012).
- (19) A. P. Severyukhin, N. N. Arsenyev, N. Pietralla, and V. Werner, Phys. Rev. C 90, 011306(R) (2014).
- (20) D. H. Youngblood, H. L. Clark, and Y.-W. Lui, Phys. Rev. Lett. 82, 691 (1999).
- (21) D. H. Youngblood, Y.-W. Lui, H. L. Clark, B. John, Y. Tokimoto, and X. Chen, Phys. Rev. C 69, 034315 (2004).
- (22) D. Patel, U. Garg, M. Fujiwara, T. Adachi, H. Akimune, G. P. A. Berg, M. N. Harakeh, M. Itoh, C. Iwamoto, A. Long, J. T. Matta, T. Murakami, A. Okamoto, K. Sault, R. Talwar, M. Uchida, and M. Yosoi, Phys. Lett. B 726, 178 (2013).
- (23) D. Patel, U. Garg, M. Itoh, H. Akimune, G. P. A. Berg, M. Fujiwara, M. N. Harakeh, C. Iwamoto, T. Kawabata, K. Kawase, J. T. Matta, T. Murakami, A. Okamoto, T. Sako, K. W. Schlax, F. Takahashi, M. White, and M. Yosoi, Phys. Lett. B 735, 387 (2014).
- (24) N. Paar, D. Vretenar, E. Khan, and G. Colò, Rep. Prog. Phys. 70, 691 (2007).
- (25) E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A 635, 231 (1998).
- (26) J. Terasaki, J. Engel, M. Bender, J. Dobaczewski, W. Nazarewicz, and M. Stoitsov, Phys. Rev. C 71, 034310 (2005).
- (27) R. Julin, J. Kantele, J. Kumpulainen, M. Luontama, A. Passoja, W. Trzaska, E. Verho, and J. Blomqvist, Phys. Rev. C 36, 1129 (1987).
- (28) A. Heusler, T. Faestermann, R. Hertenberger, H.-F. Wirth, and P. von Brentano, Phys. Rev. C 91, 044325 (2015).
- (29) E. G. Lanza, M. V. Andrés, F. Catara, Ph. Chomaz, and C. Volpe, Nucl. Phys. A 613, 445 (1997).
- (30) V. Yu. Ponomarev and P. von Neumann-Cosel, Phys. Rev. Lett. 82, 501 (1999).
- (31) J. Heisenberg, J. Lichtenstadt, C. N. Papanicolas, and J. S. McCarthy, Phys. Rev. C 25, 2292 (1982).
- (32) R. H. Spear, W. J. Vermeer, M. T. Esat, J. A. Kuehner, A. M. Baxter, and S. Hinds, Phys. Lett. B 128, 29 (1983).
- (33) G. F. Bertsch, P. F. Bortignon, R. A. Broglia, and C. H. Dasso, Phys. Lett. B 80, 161 (1979).
- (34) P. F. Bortignon and R. A. Broglia, Nucl. Phys. A 371, 405 (1981).
- (35) M. Yeh, P. E. Garrett, C. A. McGrath, S. W. Yates, and T. Belgya, Phys. Rev. Lett. 76, 1208 (1996).
- (36) F.J. Dyson and M.L. Mehta, J. Math. Phys. 4, 701 (1963).