arxiv: v1 [nucl-th] 25 Jan 2015

Transcription

1 Recet Advaces i the Alicatio of the Shell Model Mote Carlo Aroach to Nuclei arxiv: v1 [ucl-th] 25 Ja 215 Y. Alhassid 1, M. Boett-Matiz 1, A. Mukherjee 2, H. Nakada 3 ad C. Öze 4 1 Ceter for Theoretical Physics, Sloae Physics Laboratory, Yale Uiversity, New Have, Coecticut 652, USA 2 ECT*, Villa Tambosi, I Villazzao, Treto, Italy 3 Deartmet of Physics, Graduate School of Sciece, Chiba Uiversity, Iage, Chiba , Jaa 4 Faculty of Egieerig ad Natural Scieces, Kadir Has Uiversity, Istabul 3483, Turkey yoram.alhassid@yale.edu Abstract. The shell model Mote Carlo (SMMC) method is a owerful techique for calculatig the statistical ad collective roerties of uclei i the resece of correlatios i model saces that are may orders of magitude larger tha those that ca be treated by covetioal diagoalizatio methods. We review recet advaces i the develomet ad alicatio of SMMC to mid-mass ad heavy uclei. 1. Itroductio Most microscoic studies of heavier uclei are based o mea-field methods such as desity fuctioal theory. However, imortat correlatio effects beyod the mea field ca be missed. The cofiguratio-iteractio (CI) shell model aroach accouts for correlatios ad shell effects but covetioal diagoalizatio methods are limited by the size of the may-article model sace. The auxiliary-field Mote Carlo (AFMC) method eables calculatio of thermal observables ad groud state roerties i model saces that are may orders of magitude larger tha model saces that ca be treated by covetioal diagoalizatio methods. I uclear hysics this method is also kow as the shell model Mote Carlo (SMMC) method [1, 2]. SMMC has rove to be a owerful method i calculatig statistical ad collective roerties of uclei [3, 4, 5, 6, 7]. Here we review recet advaces i the develomet ad alicatio of SMMC to uclei. I articular, we describe a method to circumvet the sig roblem that arises from the rojectio o a odd umber of articles ad its alicatio to calculate airig gas i mid-mass uclei ad level desities of odd-mass isotoes. We also discuss the alicatio of SMMC to describe the microscoic emergece of collectivity i heavy rare-earth uclei ad reset results for the collective ehacemet factors of level desities i such uclei. 2. Auxiliary-field Mote Carlo (AFMC) method The AFMC method utilizes the Hubbard-Stratoovich trasformatio [8], i which the imagiary-time roagator e βh of a ucleus that is described by a Hamiltoia H at iverse

2 temerature β = 1/T is rereseted by a fuctioal itegral of oe-body roagators of oiteractig ucleos over exteral auxiliary fields e βh = D[σ]G σ U σ. (1) Here G σ is a Gaussia weight ad U σ is the roagator of o-iteractig ucleos movig i time-deedet auxiliary fields σ(τ). The thermal exectatio value of a observable O is O T = Tr(e βh O) Tre βh = D[σ]Wσ Φ σ O σ, (2) D[σ]Wσ Φ σ where W σ = G σ Tr U σ is a ositive-defiite weight fuctio ad Φ σ = Tr U σ / Tr U σ is the Mote Carlo sig fuctio. The quatity O σ = Tr (OU σ )/Tr U σ is the thermal exectatio value of the observable for a give cofiguratio σ of the auxiliary fields. I SMMC, we choose samles σ k distributed accordig to W σ ad estimate the observables from O T k O σ k Φ σk / k Φ σ k. AFMC ofte suffers from a sig roblem at low temeratures. However, there are classes of good-sig iteractios that are free of the sig roblem ad are comosed of the domiat collective comoets of effective uclear iteractios [9]. Small bad-sig comoets ca be treated by the extraolatio method of Ref. [2]. 3. Circumvetig the odd-article sig roblem Alicatios of SMMC to odd-eve ad odd-odd uclei have bee hamered by a sig roblem that origiates from the rojectio o a odd umber of articles. Such a sig roblem occurs at low temeratures eve for good-sig iteractios ad makes it difficult to estimate accurately the groud-state eergy of a odd-article system. A breakthrough was a method we itroduced to calculate accurately the groud-state eergy of a odd-article system that circumvets this sig roblem [1]. It is based o the scalar imagiary-time sigle-article Gree s fuctios G ν (τ) = T a νm (τ)a νm(). (3) m Here T deotes time orderig, ν = (lj) deotes the radial quatum umber, orbital agular mometum l ad total si j of the sigle-article levels, ad m is the magetic quatum umber of j. I a asymtotic regime of sufficietly large τ but τ β, the Gree s fuctios of a eve-eve ucleus A have the form G ν (τ) e β[e j(a±1) E g.s.(a)] τ, where E j (A ± 1) is the lowest eergy of si J = j i the eighborig eve-odd uclei A±1 ad E g.s. is the groud-state eergy of the ucleus A. Thus E j (A ± 1) ca be determied from the sloe of the calculated l G ν ad the groud-state eergy of the eve-eve ucleus. The groud-state eergy of the eve-odd ucleus is the foud by miimizig E j (A ± 1) over j. We used this method to calculate eutro airig gas i mid-mass uclei from the secodorder eergy differece i eutro umber = ( ) N [E(Z, N +1)+E(Z, N 1) 2E(Z, N)]/2, where E(Z, N) is the groud-state eergy of the ucleus with Z rotos ad N eutros. Fig. 1 shows eutro airig gas for mid-mass isotoes. The calculatios were carried out i the comlete fg 9/2 shell with the iteractio of Ref. [3]. The SMMC results (solid circles) are i overall good agreemet with the exerimetal values (oe circles). 4. Level desities SMMC has bee successful i calculatig uclear state desities ρ(e x ) versus excitatio eergy E x i the resece of correlatios. The desity calculated i SMMC is the state desity, i which each level with si J is couted 2J + 1 times. However, ofte the desity measured i the exerimets is the level desity, i which each level with si J is couted oce.

3 (MeV) Ti Cr Fe Z Ge A Figure 1. Neutro airig gas as a fuctio of mass umber A i families of fg 9/2 -shell isotoes. The gas calculated i SMMC (solid circles) are comared with the exerimetal gas (oe circles). Adated from Ref. [1]. We recetly itroduced a method to calculate the level desity directly i SMMC [11]. Deotig by ρ M the level desity at give value M of the agular mometum comoet J z, the level desity ρ is give by ρ = ρ M= for eve-mass uclei ad by ρ = ρ M=1/2 for odd-mass uclei. J z rojectio ca be carried out exactly i SMMC as discussed i Ref. [5]. Level desities were recetly measured by roto evaoratio sectra i a family of ickel isotoes [12]. We calculated these level desities [13] withi the fg 9/2 shell usig the Hamiltoia of Ref. [3]. The calculatio of the groud-state eergies of the odd-mass ickel isotoes, ad hece their excitatio eergies, was made ossible by the method discussed i Sec. 3. Fig. 2 shows the level desities of The SMMC results are i close agreemet with the level desities determied from various exerimetal data sets ~ (MeV _ 1) ρ E x (MeV) Figure 2. Level desities ρ i a family of ickel isotoes. The SMMC level desities (solid circles) are comared with level desities determied from various exerimetal data sets: roto evaoratio sectra (oe squares ad quasi-cotiuous lies), level coutig data (histograms) ad eutro resoace data whe available (triagles). Take from Ref. [13]. 5. Microscoic emergece of collectivity i heavy uclei Heavy uclei are kow to exhibit various tyes of collectivity. Tyically, uclei ear shell closure are sherical ad dislay vibratioal collectivity, while mid-shell uclei are strogly

4 deformed ad dislay rotatioal collectivity. Differet tyes of collectivity are well described by emirical models but much less is uderstood microscoically. Here we use SMMC to describe the emergece of collectivity i heavy uclei withi the framework of the CI shell model aroach Thermal sigatures of collectivity A give tye of collectivity is usually idetified by its associated sectrum. While SMMC is a owerful method for calculatig thermal ad groud-state roerties i very large model saces, it does ot rovide detailed sectroscoic iformatio. To overcome this difficulty, we have idetified a thermal observable whose low-temerature behavior is sesitive to the tye of collectivity. Such a observable is J 2 T, where J is the total agular mometum. Assumig a eve-eve ucleus with a rotatioal or a vibratioal groud-state bad ad a excitatio eergy E 2 + of the first 2 + state, we have [6, 7] J 2 T 3 e E 2 + /T (1 e E 2 + /T ) 2 vibratioal bad 6 E T 2 + rotatioal bad. (4) <J 2 > T Dy <J 2 > T Sm.1 T (MeV) T (MeV) Figure 3. Left: J 2 T i 162 Dy. The SMMC results (solid circles) are comared with fits to the rotatioal model (solid lie) ad vibratioal model (dashed-dotted lie). Adated from Ref. [6]. Right: J 2 T i 148 Sm. The SMMC results (oe circles) are comared with fits to the vibratioal model (solid lie) ad rotatioal model (dashed lie). We carried out calculatios for rare-earth uclei usig the 5 82 shell lus 1f 7/2 for rotos, ad the shell lus h 11/2, 1g 9/2 for eutros. The Hamiltoias we used are give i Refs. [6, 7]. The left ael of Fig. 3 shows J 2 T for the deformed rare-earth ucleus 162 Dy. The SMMC results are comared with fits to both the rotatioal ad vibratioal models. The agreemet of the SMMC results with the rotatioal model cofirms that we are able to describe the rotatioal character of 162 Dy i the framework of a trucated sherical shell model aroach. The fit to the rotatioal formula gives E 2 + = 84.5 ± 8.9 kev, i agreemet with the exerimetal result of E 2 + = 8.6 kev. The right ael of Fig. 3 shows J 2 T for the sherical rare-earth ucleus 148 Sm. Here the SMMC results follow closely the vibratioal formula fit with E 2 + = 538 ± 31 kev, i agreemet with the exerimetal value of E 2 + = 55 kev The crossover from vibratioal to rotatioal collectivity Heavy uclei are kow to exhibit a crossover from vibratioal to rotatioal character as the umber of eutros icreases from shell closure towards mid-shell. I a mea-field aroximatio, this is described by a hase trasitio from sherical to deformed uclei. We have studied two such families of rare-earth isotoes, Sm ad ad [7].

5 <J 2 > T T (MeV) Figure 4. J 2 T as a fuctio of temerature T for a family of eve eodymium isotoes. The solid lies describe the exerimetal values extracted usig Eq. (5). For we oly use the discrete levels i Eq. (5) (dashed lies). Adated from Ref. [7]. Fig. 4 shows J 2 T for a family of eve eodymium isotoes. The oe circles are the SMMC results. We observe a crossover from a soft resose to temerature i , tyical of a vibratioal ucleus, to a rigid resose i 152, tyical of a rotatioal ucleus. The solid lies are extracted from exerimetal results usig J 2 T = ( 1 N ) J i (J i + 1)(2J i + 1)e Ei/T + de x ρ(e x ) J 2 Ex e Ex/T Z(T ) i E N, (5) where Z(T ) = N i (2J i + 1)e Ei/T + E N de x ρ(e x )e Ex/T is the corresodig exerimetal artitio fuctio. Here E i are the measured low-lyig eergy levels with si J i, which we assume to be comlete u to to certai threshold eergy E N. Above the threshold eergy, we use a backshifted Bethe formula for the level desity ρ(e x ) whose arameters are determied from a fit to level coutig data at low eergy ad eutro resoace data at the eutro resoace threshold. We use i Eq. (5) the si cutoff model estimate J 2 Ex = 3σ 2, where σ 2 = IT/ h 2 is calculated assumig a rigid-body momet of iertia I. We observe a overall good agreemet betwee the SMMC values of J 2 T ad the exerimetally determied values. 6. Collective ehacemet factors of level desities i heavy uclei SMMC state desities i rare-earth uclei were foud to be i good agreemet with exerimetal data [14]. Collective states i these uclei lead to ehacemet of their state desity, which is described by a collective ehacemet factor K. The decay of K with excitatio eergy is oe of the least uderstood toics i the modelig of level desities [15], ad it is ofte arameterized by emirical formulas. We defie a collective ehacemet factor by K = ρ SMMC /ρ HFB, where ρ SMMC is the total SMMC state desity ad ρ HFB is the level desity calculated i the fiitetemerature Hartree-Fock-Bogoliubov (HFB) aroximatio. The latter describes the desity of itrisic states ad thus the above ratio is a measure of the ehacemet i the umber of states due to collective states that are built o to of the itrisic states. Fig. 5 shows the calculated ehacemet factor K as a fuctio of excitatio eergy E x for the eve eodymium isotoes We also show by arrows the eergies of the various thermal hase trasitios that occur i the HFB aroximatio. The thi arrows corresod to the airig trasitios (for both rotos ad eutros), while the thick arrows describe the shae trasitios. The sherical uclei do ot suort rotatioal bads, so the ehacemet factor K is due to vibratioal states aloe. We observe that this vibratioal ehacemet factor decays to 1 at excitatio eergies above the roto ad eutro airig trasitios. The heavier eodymium isotoes are deformed i their groud state. We observe that K has a local miimum above the airig trasitios, which we attribute to the iterlay betwee

6 K E x (MeV) Figure 5. Collective ehacemet factor K (oe circles) of the eve eodymium isotoes versus excitatio eergy E x (see text for the defiitio of K). The thi arrows corresod to the roto ad eutro airig trasitios, while the thick arrows describe the shae trasitios. Adated from Ref. [14]. the decay of vibratioal collectivity ad the ehacemet due to rotatioal collectivity. The lateau that follows the local miimum describes the ehacemet from rotatioal states. This rotatioal ehacemet decays i the viciity of the shae hase trasitio. 7. Coclusio The SMMC is a owerful method for calculatig statistical ad collective roerties of uclei i the framework of the CI shell model i very large model saces. We reviewed a recet method to circumvet the odd-article sig roblem ad its alicatio to the calculatio of airig gas ad level desities of odd-mass uclei. We discussed recet alicatios of SMMC to heavy rare-earth uclei ad i articular the microscoic emergece of collectivity i such uclei. Ackowledgmets It is a leasure to dedicate this article to Aldo Covello i recogitio of his imortat cotributios to uclear structure hysics. This work was suorted i art by the DOE grat No. DE-FG-291-ER-468, by the JSPS Grat-i-Aid for Scietific Research (C) No , ad by the Turkish Sciece ad Research Coucil (TÜBİTAK) grat No. ARDEB T973. Comutatioal cycles were rovided by the NERSC high erformace comutig facility ad by the High Performace Comutig Ceter at Yale Uiversity. Refereces [1] Lag G H, Johso C W, Kooi S E ad Ormad W E 1993 Phys. Rev. C [2] Alhassid Y, Dea D J, Kooi S E, Lag G H, ad Ormad W E 1994 Phys. Rev. Lett [3] Nakada H ad Alhassid Y 1997 Phys. Rev. Lett [4] Alhassid Y, Liu S ad Nakada H 1999 Phys. Rev. Lett [5] Alhassid Y, Liu S ad Nakada H 27 Phys. Rev. Lett [6] Alhassid Y, Fag L ad Nakada H 28 Phys. Rev. Lett [7] Öze C, Alhassid Y ad Nakada H 213 Phys. Rev. Lett [8] Hubbard J 1959 Phys. Rev. Lett. 3 77; Stratoovich R L 1957 Dokl. Akad. Nauk. S.S.S.R [9] Dufour M ad Zuker A P 1996 Phys. Rev. C [1] Mukherjee A ad Alhassid Y 212 Phys. Rev. Lett [11] Alhassid Y, Boett-Matiz M, Liu S ad Nakada H arxiv: [12] Voiov A V, Grimes S M, Brue C R, Massey T ad Schiller A 212 EPJ Web of Cofereces 21 52; Voiov A V (rivate commuicatio) [13] Boett-Matiz M, Mukherjee A ad Alhassid Y 213 Phys. Rev. C Raid Comm [14] Alhassid Y, Öze C ad Nakada H 214 Nuclear Data Sheets [15] Caote R et al. 29 Nuclear Data Sheets , ad refereces therei