Towards a universal nuclear structure model Xavier Roca-Maza Congresso del Dipartimento di Fisica Milano, June 28 29, 217 1
Table of contents: Brief presentation of the group Motivation Model and selected results Conclusions 2
Brief presentation of the group: The group interests focus on theoretical nuclear structure studies: that is, on the study of i) the strong interaction that effectively acts between nucleons in the medium; and ii) suitable many-body techniques to understand the very diverse nuclear phenomenology. P. F. Bortignon (PO) X. Roca-Maza (RTD) G. Colò (PO) E. Vigezzi (Dir. Ric.) 3
Some open lines of research in which the group is involved: Density functional theory: successful approach to the study of the nuclear matter EoS or of the mass, size and collective excitations in nuclei Nuclear field theory: successful approach to the study of fragmentation of sp states, resonance widths or half-lives Parity violating and conserving e-n scattering: characterize the electric and weak charge distribution in nuclei Superfluidity in nuclei: driven by the pairing interaction Nucleon transfer reactions: probe spectroscopic properties and superfluidity Astrophysical appliactions: neutron stars, electron capture, β decay, Gamow-Teller resonances 4
Today... I will concentrate on some of the basic complexities of the nuclear many-body problem and briefly explain the way along our group is working 5
Motivation: The Nuclear Many-Body Problem Nucleus: from few to more than 2 strongly interacting and self-bound fermions (neutrons and protons). Complex systems: spin, isospin, pairing, deformation,... 3 of the 4 fundamental forces in nature are contributing to the nuclear phenomena (as a whole driven by the strong interaction). β decay: weak process Nuclei: self-bound system by the strong interaction [In the binding energy B coul B strong /(3 to 1)] α decay: interplay between the strong and electromagnetic interaction 6
Motivation: The Nuclear Many-Body Problem Underlying interaction: the so called residual strong interaction=nuclear force, the one acting effectively between nucleons, has not been derived yet from first principles as QCD is non-perturbative at the low-energies relevant for the description of nuclei. The nuclear force in practice: effective potential fitted to nucleon-nucleon scattering data in the vacuum. 3 Body force are needed. 4 Body? 7
Motivation: The Nuclear Many-Body Problem State-of-the-art many-body calculations based on these potentials are not conclusive yet: Which parametrization of the residual strong interaction should we use? Which many-body technique is the most suitable? Exp. B/A(.16 fm 3 ) = 16 MeV 8
So, two important points to remember: Working with the exact Hamiltonian and most general wave function is not possible at the moment Best attempts so far show too large discrepancies and can be applied only to nuclei up to mass number 1-2 approx. 9
Motivation: So what to do? Effective interactions solved at the Hartree-Fock or Mean-Field but fitted to experimental data in many-body system have been shown to be successful in the description of bulk properties of all nuclei (masses, nuclear radii, deformations, Giant Resonances...) These effective models can be understood as an approximate realization of a nuclear energy density functional E[ρ]. Density Functional Theory rooted on the Hohenberg-Kohn theorem exact functional exists. The advantage of these approximate E[ρ] is that they are quite versatile: nowadays our unique tool to selfconsistently access the ground state and some excited state properties of ALL atomic nuclei 1
Applicability of nuclear E[ρ] as compared to other methods 11
Nuclear mean-field models (or EDFs) Ψ H Ψ Φ H eff Φ = E[ρ] where Φ is a Slater determinant and ρ is a one-body density matrix we expect reliable description for the expectation value of one-body operators. Main types of models: Relativistic based on Lagrangians where effective mesons carry the interaction (π,σ,ω...). Non-relativistic based on effective Hamiltonians (Yukawa, Gaussian or zero-range two body forces) Both give similar results phenomenological models difficult to connect to the residual strong interaction 12
Examples: Binding energies Relativistic model by Milano and Barcelona groups 6 4 DD-MEδ (E th -E exp )/ E exp *1 [%] 2-2 +.5% -.5% -4 (c) 4 8 12 16 2 24 28 Mass number A [PHYSICAL REVIEW C 89, 5432 (214)] Remarkable accuracy on thousands of measured binding energies 13
Examples: Charge radii Theory-lines/Experiment - circles 5 4.8 Charge radii r c [fm] 4.6 4.4 4.2 4 3.8 3.6 Sr (Z=38) Kr (Z=36) Se (Z=34) B (b) 3.5 2 4 6 8 1 Neutron number N 6 Cm (exp) Pu (exp) Charge radii r c [fm] 3 2.5 2 1.5 1 2 3 4 Neutron number N A (c) Charge radii r c [fm] 5.8 5.6 5.4 5.2 5 Pb (Z=82) Hg (Z=8) Pt (Z=78) C Cm (Z=96) Pu (Z=94) U (Z=92) 8 1 12 14 Neutron number N (d) [PHYSICAL REVIEW C 89, 5432 (214)] 14
Examples: Giant Monopole and Dipole Resonances Non-relativistic model by Milano and Aizu (Japan) groups 1 2 R ISGMR [fm 4 MeV 1 ] 12 9 6 3 SLy5 SAMi Exp. 28 Pb SLy5 SAMi Exp. 28 Pb 25 2 15 1 5 R IVGDR [fm 2 MeV 1 ] 1 12 14 16 18 E x [MeV] 8 12 16 E x [MeV] R=nuclear response function (in dipole resonance is related with the probability of a photon absorption by the nucleus or σ γ ) Good description excitation energy and integrated R but not the width of the resonance. 15
Examples: Gamow Teller Resonance Gamow-Teller Resonance driven by strong nuclear force (analogous transitions to β decay). R GT [MeV 1 ] 6 4 2 28 Pb SAMi Exp. PRC 85, 6466 (212) Exp. PRC 52, 64 (1995) 5 1 15 2 25 E x [MeV] [Phys.Rev. C86 (212) 3136] 16
Mean-field approach: overall good description of ground state and excited state properties in nuclei It is beyond the MF approach: accurate description of the fragmentation of the single-particle and collective states 17
Examples: Observables beyond the MF approach Single particle (sp) states: Density of the system is well described within the MF approach (E[ρ]) while sp are not satisfactorily reproduced. [J. Phys. G. 44 (217) 4516] 18
Examples: Observables beyond the MF approach Spectroscopic factor S is associated to n,j,l: Removal probability for valence protons * From theory: For a bound A 1 final state S theo = d p Φ A 1 a p Φ A 2 * Transfer reaction: A A 1 dσ = dω exp dσ S exp dω theo [Nuclear Physics A 553 (1993) 297-38] 19
Possible solution to these problems: PVC There is NO explicit interplay of collective motion (e.g. giant resonances) and single-particle motion in the mean-field approach Solution: take into account their interplay For example, in spherical nuclei: mainly sp+surface vibrations Particle Vibration Coupling model (PVC) while in deformed nuclei: mainly sp+rotations PVC model in a diagramtic way (e.g. Σ or sp potential): Selection of most relevant diagrams same V eff at all vertices (fitted at Mean-Field level) 2
Some selected examples 21
Examples: Observables beyond the MF approach Single particle states: E [MeV] 8 4-4 -8-12 -16-2 -24 2f 5/2 1h 9/2 3p 1/2 3p 3/2 2f 7/2 RMF QVC EXP 4 2-2 -4-6 3s 1/2 1h 11/2 2d 3/2 2d 5/2 1g 7/2 1h 11/2 3s -8 1/2 2d 3/2-1 2d 1g 5/2 9/2-12 1g 7/2 2p -14 1/2 1g 9/2 2p 3/2-16 Neutrons 12-18 Sn [Phys. Rev. C 85, 2133(R)] RMF QVC EXP Protons 22
Examples: Observables beyond the MF approach Total single particle strength associated to Collective strength n, j, l: transfer reaction photoabsorption cross section in 12 Sn - ρ lj [MeV -1 ].8.6.4.2.4.3.2.1 1.8.6.4.2 1.8.6.4.2 12 Exp.data HF+PVC(2 +,3 -,4 +,5 - ) s 1/2 1 8 4 Ca SLy5 d 3/2 [g.s. of 39 Ca] Excitation energy in 39 Ca [MeV] 6 4 d 5/2 p 3/2 2.8.6.4.2.3.25.2.15.1.5.6.4.2.8.6.4.2.8.6.4.2.8.6.4.2 f 5/2 g 9/2 d 5/2 p 1/2 p 3/2 f 7/2 [g.s. of 41 Ca] 2 4 6 8 1 12 Excitation energy in 41 Ca [MeV] [Phys. Rev. C 86, 34318] (Remember spectroscopic factor) 23
Examples: Observables beyond the MF approach Optical potential model based on PVC. Neutron elastic cross section by 16 O at 28 MeV dσ/dω [barn/sr] 1 1 1 1-1 1-2 1-3 E = 28 MeV Exp. cpvc HF 1-4 3 6 9 12 15 18 θ [deg] [Phys. Rev. C 86, 4163(R) (212)] 24
PVC: How to renormalize the effective interacrion? PVC includes many-body correlations not explicitely included at the Mean-Field level. While V eff fitted at the Mean-Field to experimental data. This implies double counting since parameters contain correlations beyond Mean-Field Solution: refit the interaction, that is, renormalize the theory. Divergences may also appear in beyond Mean-Field calculations if zero-range V eff are used. The group is now studiyng different strategies for the renormalization of V eff. Essentially: In a simplified case (2nd order term is kept, no summation up to infinity). Cutoff and dimensional renormalization have been performed In the full PVC approach by appliyng the subtraction method 25
Conclusions Effective interactions solved at Hartree-Fock or Mean-Field level have been shown to be successful in the description of all nuclei (masses, nuclear sizes, deformations, Giant Resonances... These effective models can be understood as an approximate realization of a nuclear energy density functional E[ρ] exact functional exist. An accurate description of the fragmentation of single-particle and collective states is reached beyond the MF approach Particle Vibration model improve the description of these observables, although renormalization of the interaction needs to be investigated. 26
Thank you! 27
Simplified PVC: Cutoff renormalization on EoS E potential = V + ν V 2 E E ν ν where is the GS and ν an excited state [Phys. Rev. C 94, 34311 (216)] E/A [MeV] E (2) /A [MeV] -5-1 -15-2 -25-3 (a) -35,5,1,15,2,25,3-5 -1-15 -2-25 SLy5 mean field Sec. Order, Λ=.5 fm -1 Sec. Order, Λ=1 fm -1 Sec. Order, Λ=1.5 fm -1 (b) -3 Sec. Order, Λ=2 fm -1.5.1.15.2.25.3 ρ [fm -3 ] E/A [MeV] -3-6 -9-12 -15 SLy5 mean field Sec. Order, Λ=.5 fm -1 Sec. Order, Λ=1 fm -1 Sec. Order, Λ=1.5 fm -1 Sec. Order, Λ=2 fm -1-18.5.1.15.2.25.3 ρ [fm -3 ] Note: since we do not know the TRUE EoS we chose one calculation of EoS as a benchmark 28
Full PVC: Subtraction method Subtraction method is based on a simple idea: modify the theory so that the expectation value of any one-body operator (idealy accurate at the MF) do not change with respect the MF prediction. * Its realization is simple ( Σ sub (E) = Σ Nosub (E) Σ HF induced eff. interaction) * Unfortunatelly corrects only approximately some of the studied observables (such as different moments of the response function) * One (different) subtraction recipe should be applied for each observable that needs to be renormalized. 1 3 xs(quadrupole) (fm 4 / MeV) 4 3 2 1 RPA PVC - Sub. PVC - No Sub. SLy5 28 Pb 5 1 15 2 E (MeV) [Phys. Rev. C 94, 34311 (216)] 29