G. Colò

1. Effect of particle-vibration coupling on the single-particle states: a consistent study within the Skyrme framework JAPAN-ITALY EFES Workshop Torino, Sep. 6-8 2010 G. Colò

2. Co-workers • P.F. Bortignon (Università degli Studi and INFN, Milano, Italy) • H. Sagawa (The University of Aizu, Japan) Topic of this talk • What is the status of modern particle-vibration coupling (PVC) calculations ? • How well can we reproduce s.p. spectra ?

3. Z N The problem of the single-particle states The description in terms of indipendent nucleons lies at the basis of our understanding of the nucleus, but in many models the s.p. states are not considered (e.g., liquid drop, geometrical, or collective models). There is increasing effort to try to describe the s.p. spectroscopy. The shell model can describe well the np-nh couplings, less well the coupling with shape fluctuations. There is an open debate to which extent density-functional methods can describe the s.p. spectroscopy (without dynamical effects).

4. 1-body density matrix Slater determinant Energy density functionals (EDFs) • The minimization of E can be performed either within the nonrelativistic or relativistic framework → Hartree-Fock or Hartree equations • In the former case one often uses a two-body effective force and defines a starting Hamiltonian; in the latter case a Lagrangian is written, including nucleons as Dirac spinors and effective mesons as exchanged particles. • 8-10 free parameters (typically). Skyrme/Gogny vs. RMF/RHF. • The linear response theory describes the small oscillations, i.e. the Giant Resonances (GRs) or other multipole strength → (Quasiparticle) Random Phase Approximation or (Q)RPA • Self-consistency !

5. EDFs vs. many-body approaches In the standard theory of quantum many-particle systems one has a more general picture in which dynamical effects play a significant role. A set of closed equations for G, Π(0), W, Σ, Γ can be written (v12 given). They can be found e.g. in the famous paper(s) by L. Hedin in the case of the Coulomb force – they hold more generally (except…). The Dyson equation reads in terms of the one-body Green’s function EDF = static limit

6. NPA 553, 297c (1993) A. Gade et al., PRC 77 (2008) 044306 Experiment: (e,e’p), as well as (hadronic) transfer or knock-out reactions, show the fragmentation of the s.p. peaks. • Problems: • Ambiguities in the definition: use of DWBA ? Theoretical cross section have ≈ 30% error. • Consistency among exp.’s. • Dependence on sep. energy ? S ≡ Spectroscopic factor

7. + + … = Particle-vibration coupling (PVC) for nuclei 2nd order PT: ε + <Σ(ε)> In the Dyson equation we assume the self-energy is given by the coupling with RPA vibrations In a diagrammatic way Particle-vibration coupling

8. + … + = G P. Papakonstantinou et al., Phys. Rev. C 75, 014310 (2006) W • For electron systems it is possible to start from the bare Coulomb force: • In the nuclear case, the bare VNN does not describe well vibrations ! Phys. Stat. Sol. 10, 3365 (2006) TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: AAAAAAAAA

9. Second-order perturbation theory In most of the cases the coupling is treated phenomenologically. In, e.g., the original Bohr-Mottelson model, the phonons are treated as fluctuations of the mean field δU and their properties are taken from experiment. No treatment of spin and isospin. • THE MAIN PROBLEM (in the nuclear case): • A LOT OF UNCONTROLLED APPROXIMATIONS HAVE BEEN MADE WHEN IMPLEMENTING THE THEORY IN THE PAST !

10. C. Mahaux et al., Phys. Rep. 120, 1 (1985) Despite quantitative differences old calculations agree that one needs to introduce dynamical effects to explain the density of s.p. levels. This is associated to the effective mass. One calculates the expressions corresponding to the diagrams using standard rules. The signs of the denominators are such that “as a rule” particle states (hole states)close to the Fermi energy are shifted downwards (upwards).

11. RPA microscopic Vph • RMF + PVC calculations by P. Ring et al.: they also approximate the phonon part. • The most “consistent” calculations which are feasible at present start from Hartree or Hartree-Fock with Veff, by assuming this includes short-range correlations, and add PVC on top of it. • Very few ! • Pioneering Skyrme calculation by V. Bernard and N. Van Giai in the 80s (neglect of the velocity-dependent part of Veff in the PVC vertex, approximations on the vibrational w.f.)

12. A consistent study within the Skyrme framework We have implemented a version of PVC in which the treatment of the coupling is exact, namely we do not wish to make any approximation in the vertex. The whole phonon wavefunction is considered, and all the terms of the Skyrme force enter the p-h matrix elements Our main result: the (t1,t2) part of Skyrme tend to cancel quite significantly the (t0,t3) part. We have also compared with the Landau-Migdal approximation.

13. 40Ca (neutron states) Δεi is the expectation value of Re Σi

14. The (t0,t3) part of the interaction, reads, for IS phonons The velocity-dependent part is, using the Landau-Migdal approximation Use of the Landau-Migdal approximation Exact = -0.95 MeV

15. 40Ca (neutron states) • The tensor contribution is in this case negligible, whereas the PVC provides energy shifts of the order of MeV. • The r.m.s. difference between experiment and theory is: • σ(HF+tensor) = 0.95 MeV • σ(including PVC) = 0.62 MeV

16. 208Pb (neutron states) • The tensor contribution plays a role in this case. In principle all parameters should be refitted after PVC. • The r.m.s. difference between experiment and theory is: • σ(HF+tensor) = 1.51 MeV • σ(including PVC) = 1.21 MeV

17. Still to be done… • Do we learn in this case by looking at isotopic trends ? We need pairing and proper continuum treatment if we wish to go towards weakly bound systems (cf. K. Mizuyama talk). • We have a quite large model space of density vibrations. Do we miss important states which couple to the particles ? • Do we need to go beyond perturbation theory ? • Is the Skyrme force not appropriate ? TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: AAAAAAAAA

18. A few conclusions • The aim of this contribution consists in showing the feasibility of fully MICROSCOPIC calculations including the particle-vibration coupling. • Our approach takes into account Skyrme forces consistently. Results are (slightly) better than simple EDFs. • We can also discuss s.p. spectroscopic factor. • Technical progress in our calculations is still underway.

19. Backup slides

20. Difference between self-consistent mean field (SCMF) and energy density functionals In the self-consistent mean field (SCMF) one starts really from an effective Hamiltonian Heff = T + Veff, and THEN builds < Φ | Heff | Φ > and defines this as E. In DFT, one builds directly E[ρ]. → More general !

21. Hedin’s equations (natural units)

22. A reminder on effective mass(es) E-mass: m/mE k-mass: m/mk TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: AAAAAAAAA

23. The (time-consuming) calculations is much simplified if the interaction is simply a density-dependent delta-force: And becomes

24. Contribution of phonons with different multipolarity Upper panel: particle states. Lower panel: hole states. Signs fixed by energy denominators ω-Eint+iη

25. How to compare EDF and PVC ? ωn Since the phonon wavefunction is associated to variations (i.e., derivatives) of the denisity, one could make a STATIC approximation of the PVC by inserting terms with higher densities in the EDF.