arxiv:nucl-th/ v1 20 Dec 2000

Transcription

1 Shell model study of the isobaric chains A = 5, A = 5 and A = 5 A. Poves and J. Sánchez-Solano arxiv:nucl-th/77v Dec Departamento de Física Teórica C-XI, Universidad Autónoma de Madrid, E 89 Madrid, Spain E. Caurier and F. Nowacki Groupe de Physique Théorique, Centre de Recherches Nucléaires, Institut National de Physique Nucléaire et de Physique des Particles, Centre National de la Recherche Scientifique, Université Louis Pasteur, Boîte Postale 8, F 677 Strasbourg Cedex, France Abstract Shell model calculations in the full pf-shell are carried out for the A=5, 5 and 5 isobars. The most frequently used effective interactions for the pf-shell, KB and FPD6 are revisited and their behaviour at the N=8 and Z=8 closures examined. Cures to their -relatively minor- defaults are proposed, and a new mass dependent version called KBG is released. Energy spectra, electromagnetic transitions and moments as well as beta decay properties are computed and compared with the experiment and with the results of the earlier interactions. A high quality description is achieved. Other miscellaneous topics are addressed; the Coulomb energy differences of the yrast states of the mirror pair 5 Mn- 5 Fe and the systematics of the magnetic moments of the N=8 isotones. Key words: A=5, A=5, A=5, Shell Model, Effective interactions, Full pf-shell spectroscopy, Level schemes and transition probabilities, Gamow-Teller beta decays, Half-lives, Magnetic moments, Coulomb energy differences. PACS:.. k, 7..+z,.6.Cs,.. s Introduction The pf-shell has been the focus of a lot of activity in nuclear structure during the last years. Prompted in some cases by the large scale shell model results, that indicated the presence of a region of deformation around 8 Cr [ ], a lot of new experiments and calculations have been carried out, addressing many Preprint submitted to Elsevier Preprint June 8

2 different issues; deformed bands and band termination [], yrast traps [5], high K isomers [6], coexistence of deformed bands of natural and non-natural parity [7], effects of neutron proton pairing [8], etc. A very recent highlight has been the discovery of an excited deformed band in 56 Ni [9] coexisting with the spherical states based in the doubly magic ground state. In addition to the exact shell model diagonalizations, the new Monte Carlo techniques, SMMC and MCSM have been extensively applied to this region [ ]. Mean field descriptions of various kinds have also been used to explore different issues concerning this deformation region [,5]. In this article we extend the full pf-shell calculations up to A=5. Detailed results for 5 Cr, 5 Mn and 5 Fe using KB have been already published in refs [5,6,7] and we will not deal with them here because the new interaction KBG gives equivalent results. In the cases of 5 Cr, 5 Cr, 5 Mn and 5 Mn we have carried out the full pf-shell calculation for the yrast states only. To perform detailed spectroscopy in the full space would have demanded a huge amount of computer time, not justified by the improvement on the results, as we have checked. Hence, for the non-yrast states we shall present results in a truncated (t=5) space (no more than 5 particles are allowed to excite from the f 7/ subshell). At this truncation level, the most relevant states are sufficiently converged. As we have discussed in detail elsewhere [] our usual effective interaction KB [8,9] produces a quasiparticle gap in 56 Ni about MeV too large. Approaching the doubly magic closure, the effects of this default become more visible. That is the reason why, in a recent study of the deformed excited band of 56 Ni [9], we used a preliminary modified version of KB in order to have the correct gap. This modified version of KB was also used in the study of the M strength functions of the N=8 isotones in ref. []. We shall examine the interaction issues in section II. In sections III, IV and V we present the spectroscopic results for A=5, 5 and 5 respectively. In section VI we gather the beta decay results. In section VII we discuss the Coulomb Energy Differences (CED) between the yrast states of the mirror pair A=5. Finally, in section VIII we study the behaviour of the magnetic moments of the N=8 isotones. We close with the conclusions. Throughout the paper f stands for f 7/ (except, of course, when we speak of thepf shell)andr,generically,foranyoralloftheothersubshells(p / p / f 5/ ). Spaces of the type f n n r n +f n n r n + +f n n t r n +t () represent possibletruncations: n isdifferent fromzero ifmorethan8neutrons (or protons) are present and when t = n n we have the full space (pf) n for A = +n.

3 The interaction KB is a (mostly) monopole modification of the original Kuo- Brown one [8]. The modifications are described in detail in []. In what follows, and unless specified otherwise, we use harmonic oscillator wave functions with b =.A /6 fm; bare electromagnetic factors in M transitions; effective charges of.5 e for protons and.5 e for neutrons in the electric quadrupole transitions and moments; Gamow-Teller (GT) strength defined through B(GT) = ( ) ga g V eff στ, στ = f kσ k t k ± i, () Ji + where the matrix element is reduced with respect to the spin operator only (Racah convention []), ± refers to β ± decay, t ± = (τ x ± iτ y )/, with t + p = n and (g A /g V ) eff is the effective axial to vector ratio for GT decays, ( ) ga g V eff =.77 ( ) ga g V bare, () with (g A /g V ) bare =.599(5) []; for Fermi decays we have B(F) = τ, τ = f kt k ± i Ji + ; () half-lives, T /, are found through (f A +f ǫ )T / = 66±6 (f V /f A )B(F)+B(GT). (5) We follow ref. [] in the calculation of the f A and f V integrals and ref. [] for f ǫ. The experimental energies are used. The intrinsic quadrupole moments Q are extracted from the spectroscopic ones through Q = (J +)(J +) K J(J +) Q spec(j), for K (6) or from the BE s through the rotational model prescription B(E,J J ) = 5 6π e JK J,K Q, for K,. (7) The diagonalizations are performed in the m-scheme using a fast implementation of the Lanczos algorithm through the code antoine [5] or in J-coupled scheme using the code nathan [6]. Some details may be found in ref. [7].

4 The strength functions are obtained using Whitehead s prescription [8], explained and illustrated in refs. [9 ]. All the experimental results for which no explicit credit is given come from the electronic version of Nuclear Data Sheets compiled by Burrows []. The Interactions Following ref. [], we shall treat the effective interaction as a sum of a monopole and a multipole part. The monopole part is responsible for the energies of the closed shells (CS) and the closed shells plus or minus one particle (CS±). In our valence space, the starting point (or vacuum) is the Ca core, and the single particle energies are provided by the levels of Ca. The harmonic oscillator closure should be 8 Zr and its corresponding hole states in 79 Y. Nevertheless, and taking care of the effect of correlations, 8 Ca and 56 Ni can be also taken as reference closed shells. In this we are lucky, because the information given by 8 Zr is rather useless. For N,Z> the influence of the g 9/ orbit is very strong and the pf valence space is no longer valid. Around 8 Zr, the occupation of the orbits g 9/ and d 5/ drives the nuclei into deformed shapes [,5]. There are no experimental results available for the single hole states in 79 Y. Even if some excited levels were accessible, a particle plus rotor spectrum should be expected, coming from the coupling of the holes to the 8 Zr deformed core. With this information missing, we have to rely on indirect indications, as those coming from spherical Hartree Fock calculations (not very accurate) or from the new monopole formulae fitted to the real shell closures and their particle and hole partners [6]. Hence, our main guidance to adjust the monopole part of our interactions comes from the gaps around 8 Ca and 56 Ni and from the single particle states of 9 Ca and 57 Ni. The multipole part of the interactions will be also analysed in terms of coherent multipoles, as proposed in ref. []. This part of the interaction is left unchanged. The quasiparticle gap for 8 Ca is defined as = BE( 8 Ca) BE( 9 Ca) BE( 7 Ca) (8) and similarly for 56 Ni. In table we compare these quantities for the most popular effective interactions used in the pf-shell; KB [8,9] and FPD6 [7]. Both interactions use very similar single particle energies; ǫ 7/ =.MeV, ǫ / =.MeV, ǫ / =.MeV and ǫ 5/ =6.5MeV for KB and ǫ 7/ =.MeV, ǫ / =.89MeV, ǫ / =.9MeV and ǫ 5/ =6.9MeV for FPD6. The two body matrix elements

5 Table 8 Ca and 56 Ni gaps and single particle energies, with the interactions KB and FPD6 (energies in MeV). ǫ5 ǫ ǫ ǫ A=8 KB FPD6 KB FPD6 KB FPD6 t= Full A=56 KB FPD6 KB FPD6 KB FPD6 t= t= defining FPD6 are scaled with the mass number A as (/A).5, while KB does not incorporate any mass dependence. Notice that KB does definitely better than FPD6 for the single particle spectra of 9 Ca and 57 Ni, while for the 56 Ni gap, FPD6 is better. In addition the modification of KB s gap, we want to to make it mass dependent, in order to be able to use it safely around and beyond 56 Ni, therefore we have to adjust the monopoles anew. As the gaps are subjet to variation when correlations are allowed, some trial and error fitting of the monopole changes is needed. The final modifications of KB, defining KBG for A= (we stick to the (/A) / mass dependence) are the following: Vfp T= (KBG) = Vfp T= (KB) 5keV, Vfp T= (KBG) = Vfp T= (KB) kev, Vff T= 5/ (KBG) = Vff T= 5/ (KB) kev, V T= ff 5/ (KBG) = V T= ff 5/ (KB) 5keV, V T pp (KBG) = V T pp (KB)+keV, where p denotes any of the orbits p and p. The T= and T= modifications are different in order to recover simultaneouslythegoodgapsaround 8 Caand 56 Ni.Themodificationoftheppcentroids bears no relationship with the gaps, it is aimed to give a single hole spectrum of 8 Zrinaccordwiththepredictions ofref[6].itisimportanttoremarkthat, for the nuclei already studied with KB, KBG produces equivalent results. 5

6 Table Same as in table with the modified interactions KBG and FPD6*. ǫ5 ǫ ǫ ǫ A=8 KBG FPD6* KBG FPD6* KBG FPD6* t= Full A=56 KBG FPD6* KBG FPD6* KBG FPD6* t= t= In the case of FPD6, the gaps are nearly correct and the monopole defects amount to having the f 5/ orbit too low and the p / orbit too high, both in 9 Ca and 57 Ni. The modifications that repair that are: V T=, ff 5/ (FPD6*) = V T=, ff 5/ (FPD6)+5keV, V T=, fp / (FPD6*) = V T=, fp / (FPD6) 5keV. In table we present the resulting values around 8 Ca and 56 Ni. Notice that we have not aimed to an exact agreement. In particular, we know [9] that a t= calculations is not fully converged in 57 Ni, therefore we have let room for the f 5/ to move up and for the p / to move down a little. We also show the evolution of the quasiparticle neutron gaps through the N=8 isotones for both modified interactions in table, The gaps produced by FPD6 do not differ appreciably from those of FPD6*. The agreement with the experimental results is the best that can be reasonably aimed to. Once the monopole part under control, we can analyse the multipole hamiltonian using the method of ref. []. It amounts to diagonalize the monopolefree two body matrix elements in the particle-particle representation or in the particle-hole representation (i.e. after a Racah transformation). In the former case, the physically relevant terms are the isovector and isoscalar pairing terms, while in the latter all the multipoles will show up with different coherent strengths given by the lowest eigenvalues of the respective matrices. In table, we have listed these numbers for the two interactions and for the most important multipoles. The differences are small; basically consist in FPD6 being 5-% more intense than KB in all the channels. It is however satisfying the closeness of the values for two effective interactions derived by very different procedures. Remember that for its multipole part, KB is just equivalent to the original Kuo and Brown G-matrix [8] including the bubble correction, 6

7 Table Theoretical quasiparticle neutron gaps of the N=8 isotones (in MeV) compared with experiment. For Z> we list t= results, for the rest, full pf-shell results (energies in MeV). δ = (exp) (th) Nucleus FPD6* EXP. KBG FPD6* KBG 8 Ca Sc Ti V Cr Mn Fe Co Ni Table Strengths (in MeV) of the coherent terms of the multipole Hamiltonian. Interaction particle-particle particle-hole JT= JT= λτ = λτ = λτ = KB FPD GOGNY while FPD6 was obtained via a potential fit to selected energy levels in the pf-shell. This supports of the conclusions in [] about the universality of the multipole shell model Hamiltonian. We have purposely left to the end the line labeled GOGNY. It comes from the same analysis applied to the pf-shell two body matrix elements obtained using the density dependent interaction of Gogny []. The calculation was carried out using the single particle wave function obtained in a spherical Hartree- Fock calculation for 8 Cr in the uniform filling approximation in ref. []. In spite of the rather hybrid approach, the most important terms are again very similar to those arising from a G-matrix or from a shell model fit. In particular the agreement for the quadrupole and the spin-isospin terms is excellent. When it comes to pairing, this way of looking to the interaction 7

8 can help to overcome some language barriers between mean field and shell model practitioners. As it is evident from the table, the Gogny interaction has essentially the same amount of isovector and isoscalar pairing than the realistic interactions. Therefore, it contains the right proton-neutron pairing. Thus, if there is something to blame in the mean field calculations for N=Z nuclei, it should be rather the mean field approximations and not the interaction itself. In what follows we shall use mainly the interaction KBG. We will compare in some cases with KB in order to evaluate the importance of the changes made in the monopole behaviour. The comparison with FPD6 will show to what extent the residual differences between good behaved interactions can have spectroscopic consequences. We have not attempted to compute all the states not even to draw those we have computed. Those not shown here can be obtained on request to the authors. The isobars A=5 We have carried out the full pf-shell calculations for all the isobars. The results for 5 Cr and 5 Mn have been already published in refs. [6,,7]. The experimental data for which no specific credit is given are taken from ref. [].. Spectroscopy of 5 Ca The experimental data are compared with the calculation in fig.. All the experimental statesareplotted. Wealsoplotall thecalculatedonesupto5mev. Beyond, only the yrast states are drawn. The scarcity of spin assignment precludes a more detailed analysis. The + excitation energy is well reproduced as well as the gap in the spectrum between and MeV of excitation energy. KB gives equivalent results. No experimental information on transitions is available.. Spectroscopy of 5 Sc Figure reflects the experimental situation and our calculated spectrum. We have plotted all the experimental states even if only the lowest ones have assigned spin. All the calculated states up to the first 7 + state are shown. Aboveit,onlytheyrastband.Thegroundstatemultiplet,(f 7/ ) 9,p / iswell reproduced. The first excited state must be a + because the multiplet cannot contain two + s. The triplet at MeV, belonging to the configurations 8

9 (8 +,( )) (,(8 + )) ( & 5 + ), +, , , KBG Fig.. Energy levels of 5 Ca ,, + 5 +, +, + ( + ) , KBG Fig.. Energy levels 5 Sc. (f 7/ ) 9,(f 5/,p / ), is also well given. Beyond, the level density increases rapidly and no spin assignments are available. Notice that the experimental density of states is nicely reproduced up to MeV. 9

10 9 () ( + ) ( ) 5+ + ( + ) + ( + ) +, FPD6 KBG. Spectroscopy of 5 Ti Fig.. Energy levels of 5 Ti. 5 Ti is the stable member of the isobaric multiplet and the best known experimentally. It is semi-magic and therefore any defects in the neutron gap would be evident in the spectrum. In fig the experimental level scheme is compared with the full pf-shell results using the interactions KBG and FPD6. Up to MeV all the states are plotted. Between MeV and 5.6MeV only those with unambiguous experimental spin assignment and beyond, only the yrast. The agreement is impressive for the two interactions, perhaps with a bonus for KBG. Notice the excellent reproduction of the second and third + states as well as the perfect location of the high spin states. The second + corresponds to the excitation of two neutrons to the p orbit, while the third has a more complex structure. Table 5 Transitions in 5 Ti. B(E) (e fm ) Th (9) () 86 () The quadrupole and the magnetic moments of the yrast J = + and J =

11 states are known []. Their values, µ exp ( + ) =.89(5)µ N [], Q exp ( + ) = +8(6)efm and µ exp ( ) = +9.()µ N agree quite well with our predictions; µ( + ) = +.5µ N, Q( + ) = +6efm and µ( ) = +8.µ N. A more detailed discussion on the magnetic moments is given in section 8. Similarly the E transitions of the yrast J = +, y states are well reproduced by the calculation (see table 5).. Energy levels of 5 V The experimental information is also very rich for this nucleus. Given the high level density, we have represented in fig. all the states up to.6 MeV only. Above, just the high spins and a couple of + s and s that may have experimental counterparts. The calculation gives a good reproduction of the ground state quintuplet and locates correctly the high spin states 9 +, and. The only discrepancy affects to the bunch of, + and + states at around.5 MeV that are placed 5keV too low. 5 ( ) ( ) ( + ) +, (9 + ) +, (+),5 (+),6 (+) (8) , (7) +, 5 + +, , , + KBG Fig.. Energy levels of 5 V. The magnetic and the quadrupole moment of the J = ground state have been measured []. Their values (µ exp ( ) = ()µ N, Q exp ( ) = +.9()efm )areverywellreproducedbythecalculations(µ( ) = +.8µ N, Q( ) = +9.6efm ). The electromagnetic transitions of the yrast states are also known (table 6). Notice the excellent agreement found for the -dominant- M transitions. For the J=, E transitions the huge error bars make the

12 comparison meaningless. On the contrary, for the only J= transition measured, the agreement is very good. Table 6 Transitions in 5 V. Th. B(M) (µ N ) (µ N ) 7 +.() ()..9(). B(E) (e fm ) (e fm ) () The isobars A=5. Spectroscopy of 5 Ca The experimental data become rarer as we go far from stability. In 5 Ca only the ground state J = is assigned in the experimental scheme. It is what is expected in any reasonable calculation. This nucleus was also studied in ref. [], using the interactions KB and FPD6. As no comparison with the data is possible, we can examine the predictions of the different interactions. In fig. 5 we compare KB and KBG. As the gap is very similar for the two interactions in this nucleus, the differences must be due to the changes in the pp interaction. This is clearly seen in the figure. The four states /, /, 5/ and / are dominantly p states and remain unchanged between KB and KBG. However, the states 7/ and 5/ are swapped with a relative change of about MeV. This is due to their different structures; the 5/ is f 8 7p f 5, thus, with KBG it gains 8keV relative to the f 8 7p states, the 7/ has a structure f 7 7 p and it looses about.mev. Similar arguments explain the lowering of the bunch of states at about 5MeV excitation energy. Therefore a better experimental spectrum will be of much help in refining these monopolar changes.

13 8 / / (/,/,5/) (/,/,5/) / 9/ 7/ / / 9/ 5/ / 5/ / 7/ / 9/ / 5/ 7/ / / 9/ 7/ 5/ / 5/ / 5/ / (/ ) KB KBG / /. Spectroscopy of 5 Sc Fig. 5. Energy levels of 5 Ca. As in 5 Ca, this nucleus was studied in ref. []. Some numerical errors present in this reference were discussed and corrected in refs. [5 7]. Figure 6 shows all theknown experimental states. Thecalculated onesupto MeVenergy are also shown. Our results using KBG are slightly better than those obtained with KB and substantially better than those produced by FPD6 [5,7], specially up to MeV of excitation energy where the spacing of the calculated levels follows closely that of the experiment.. Spectroscopy of 5 Ti In figure 7 we have plotted all the levels up to MeV and a few more up to MeV which have spin and parity assignments. No high spin states are known beyond this energy. The agreement is quite good, exception made of a couple of doublets that come out inverted (J = 5 and J = 7 at.5mev and the second J = 7 with the first J = at around.5mev). The experimental spin assignment for the state around.7mev is either J = or J = 7. Our calculation predicts a J = close to this energy, while the first J = 7 state is predicted MeV higher. If we compute the 5 Ti spectrum with FPD6 there are not significant differences, because the dominant configurations are (f 8 7 p ) ν (f 7) π, therefore the f5 does not play an

14 (/,5/ ) (/,/,5/) (/,5/ ) (/,5/ ) (/ ) (/ ) (/,5/) (/,5/,7/ + ) (/) (/) / 5/ / 7/ 5/ / 7/ 5/ 5/ 9/ /,/ / 7/ 5/ 9/ / / 9/ (7/) 7/ KBG Fig. 6. Energy levels of 5 Sc. (/,7/) (5/,7/,9/) / (7/,9/) (5/,7/) / (5/ ) (7/,9/) 7/ (/ ) / 5/ (5/) 7/ / / / 7/ 9/ / 5/ 5/ 9/ / 7/ / 5/ 7/ 5/ / / / KBG Fig. 7. Energy levels of 5 Ti. important role. We can interpret most of the 5 Ti level scheme as the result of the coupling 5 Ti p. Some electromagnetic transitions are experimentaly known that are well reproduced by the calculation(see table 7). Again in this case, the E transitions

15 with J= are poorly determined experimentally, due to the dominance of the M transitions. The calculated J= B(E) s agree extremely well with the measured values. Table 7 Transitions in 5 Ti. B(M) (µ N ) (µ N ) 5 Th..6().7 B(E) (e fm ) (e fm ) () 88 5() 88 95(7) 95 6() 5. Spectroscopy of 5 V 5 V is the stable isotope in the isobar chain A = 5 and corresponds to the N=8 neutron shell closure. The experimental information is very rich, extending up to MeV. Figure 8 shows the experimental yrast band with all the known high spin assignments compared with the calculated one including three more predicted spins. The agreement is extremely good. For a more complete analysis below MeV energy, where more experimental information is available, we have plotted in figure 9 all the levels up to that energy and a few more above that can be put in correspondence with the experiment. We have calculated with KBG and FPD6, both reproduce nicely the level scheme but definitely the quality of the agreement with KBG is better. The magnetic moments of the ground state J = 7 and first excited state J = 5 as well as the quadrupole moment of the ground state are known []. Their values are: µ exp ( 7 ) = (8)µ N Q exp ( 7 ) = 5.()efm µ exp ( 5 ) = +.86()µ N while the calculation gives: µ( 7 ) = +.99µ N Q( 7 ) = 6.5efm µ( 5 ) = +.6µ N 5

16 / 9/ 7/ (5/ ) (/ ) (/ ) 9/ 7/ / 5/ 9/ / 5/ / / 9/ 7/ / 5/ 9/ / 7/ 7/ KBG Fig. 8. Yrast band of 5 V. 9/,/ 9/,/ 9/ 9/ 9/ (9/,/ ),/ (5/ ) (5/ ) 5/ / / 9/ / 9/ 5/ / 5/ / / 9/ 5/ 9/,/ / 5/ 5/ / 9/ / 9/ / 9/ / / / / 5/ 5/ 5/ 7/ 7/ 7/ FPD6 KBG Fig. 9. Energy levels of 5 V. in good agreement with the data. Furthermore, some electromagnetic transitions, mainly between yrast states are known experimentally []. The calculation does also quite well, as it is shown in table 8. 6

17 Table 8 Transitions in 5 V. B(M) (µ N ) (µ N ) 9 7 Th..6(). 6 <.77. B(E) (e fm ) (e fm ) (6) 95(8) <.8. 66(6) 78.5 Spectroscopy of 5 Cr We present in fig. the results for the yrast band of 5 Cr in a t = 5 truncation and in the full space using KBG. We have also included the low-lying J =, and 5 states, and a second state of each spin from J = 9 on. The figure shows that the full calculation do not modify substantially the results obtained at t = 5. The agreement with the experimental data is very (5/, 7/ ) (/, 5/ ) 7/ 7/ 5/ 5/ / / 7/ 7/ 5/ 5/ / / 6 5 (/, / ) (/ ) (/ ) / / 9/ / / 9/ (9/ ) 9/ 9/ (7/) / 5/ / 9/ 5/ / / 7/ 7/ / 5/ / 5/ 9/ / / 7/ 7/ / 5/ / 9/ 5/ / / 7/ KBG (t=5) KBG (full) Fig.. Yrast band of 5 Cr; experiment, t = 5 and full calculation. 7

18 satisfactory. Notice that our prediction for the spins of the triplet of states around 6 MeV is shifted down by one unit relative to the preliminary experimental assignment, that should be revised in the light of our results. The detailed comparison with all the levels up to.8mev, and with the high spin states above this energy up to J = 9 is carried out with a t = 5 calculation. The results for both KB and KBG are plotted in figure. The agreement is very good in spite of some local inversions. We can also notice that KBG improves systematically the KB results. The sequence of states J = 7, 9,J = belong to the configuration (f7), while the doublets J =, and J = 5, 7 located at 8keV and.5mev, correspond to a neutron jump to the orbit p, a sort of coupling 5 Cr p. As a consequence of the reduction of the gap, KBG put these doublets at lower energy, while leaving the yrast ones unchanged. This provides a good illustration of a direct monopole effect. A similar though more pronounced effect is seen in the lowering of the doublet J =, 5 experimentally at MeV predicted by KB 5keV too high, again a gap effect solved by KBG. Besides, the high spin states were also too high in KB and come now very close to their experimental position. 5/,7/ (9/ ) / (7/) / 9/ (/ ) / 9/ 7/ 5/ 5/ / 7/ / 5/ 9/ / / 9/ / 7/ / 9/ / 9/ / 5/,7/ 5/,/ / 9/ 7/ / 5/ 9/ / / 9/ / 7/ / 9/ / / / 5/ 7/ 9/ / 5/ / 7/ 9/ 5/ / / 7/ 7/ 7/ KB KBG Fig.. Energy levels of 5 Cr, t = 5. The magnetic moments of the ground state and the first excited state are known, (µ exp ( 7 ) =.9(5)µ N and µ exp ( ) =.86()µ N ). Our prediction for the ground state is quite accurate µ( 7 ) =.886µ N while for the J = state our result µ( ) =.µ N is off by a factor two. 8

19 In table 9 we compare the experimental and calculated transition probabilities. The B(M) s and the J= B(E) s are well reproduced. However this is not the case for the J= B(E) s, in particular the discrepancies are huge for the transitions 9 and. We have carried out a consistency test using the experimental M/E branchings, δ, and we have found that the experimental δ values are inconsistent with the experimental B(E) s but consistent with the calculated values. This is not surprising, because of the complete dominance of the M transition, that may cause large errors in the extraction of the B(E) value. In order to discard any other origin of the discrepancy we have repeated the calculation with FPD6 and the situation is the same or even worse. Table 9 Transitions in 5 Cr. B(M) (µ N ) (µ N ) Th..68().9.5(57) (969).56 B(E) (e fm ) (e fm ) (56) 67() 7 8() 8 6() 5 () 5.6 Spectroscopy of 5 Mn As for 5 Cr, we present in fig. the comparison between the experimental level scheme and our calculations in the full space and with a t = 5 truncation for the yrast band and the, doublet around MeV. Again no relevant difference is appreciated between both calculated schemes and the agreement with the experiment is very good. Therefore, a t = 5 truncation is used to get a more detailed spectroscopy for all the states below MeV and for the high spins up to J = 9 (fig. ). New data are available from a Gammasphere experiment aimed to the study of Coulomb energy differences in mirror nuclei [8]. The agreement is very good in both regions, although the high spin states are slightly shifted up. 9

20 / / / 5/ 7/ / / 5/ 7/ / / 5 9/ 7/ 5/ / 9/ 7/ 5/ / 9/ 7/ 5/ / / / / /,5/,7/ / / / / / 9/ 9/ 9/ 7/ 5/ 7/ 5/ 7/ 5/ KBG (t=5) KBG (full) Fig.. Yrast band of 5 Mn; experiment, t = 5 and full calculation. (9/ ) (7/ ) (5/ ) (/ ) ( ) / 5/ / / ( ) /,/ (7/) 5/ ( ) (5/ ) / / /,5/,7/ / (/) (9/) 9/ 7/ 5/ /,9/ 7/ 5/,/ / 5/ 7/ / / / / 9/ 7/ 5/ 7/ 5/ KBG Fig.. Energy levels of 5 Mn, t = 5. Good agreement is also obtained for the electromagnetic moments of the ground state. The experimental values: µ exp ( 5 ) =.568()µ N Q exp ( 5 ) = (7)efm

21 are very well reproduced by the KBG calculation: µ( 5 ) =.µ N ) = 5efm Q( 5 In addition, there is experimental information on electromagnetic transitions, that we compare with the calculation in table. The accord is excellent for themtransitions. FortheE sthesituationissimilar totheotherisobars;for the J= transitions the agreement is quite good, while some discrepancies are present in the J= cases, that can be very large, for instance in the transition from J = to J = 9. As in the 5 Cr case, we have found an inconsistency between the experimental values for the B(M) and B(E) and the experimental mixing parameter δ. On the contrary the experimental value is fully compatible with the computed transition probabilities. Notice the abrupt change (more than two orders of magnitude decrease) in the transition probabilities either M or E of the J = 7 state, that are spectacularly reproduced by the calculation. A very intuitive physical explanation of this isomerism will be given later, when discussing the Coulomb energy differences between the yrast states of the mirror pair 5 Fe- 5 Mn. Table Transitions in 5 Mn. B(M) (µ N ) (µ N ) Th..7().77.6(5)..66(8)..(). > B(E) (e fm ) (e fm ) (6) 5 69(67) 8 () 6(67) 5.6(5) 9.6(7).5

22 5 The isobars A=5 As in the A = 5 isobaric multiplet we have performed full hω calculations for all the isotopes in 5 Cr and 5 Mn the full calculation is limited to the yrast levels. More detailed spectroscopy is carried out with a t=5 truncation. The results of the full pf-shell calculation for 5 Fe have been already published in [5]. The experimental values are taken from ref. []. 5. Spectroscopy of 5 Ca The experimental information is scarce, as can be seen in figure. Both KB and KBG give the same excitation energy for the doublet J =, + at.5 MeV, and reproduce correctly the energy of the first excited state J = +. The states whose leading configuration has four p-particles behave like the ground state under the KBG changes, therefore we expect that their excitation energies do not change between KB and KBG. This is the case for the above mentioned doublet, whose leading configuration is f 8 7 for the second + whose configuration is f 8 7 p locate them at the same place. p p p and, therefore, both interactions ( + ) KB KBG Fig.. Energy levels of 5 Ca. The multiplet J =, +, +, predicted at 5MeV by KB has a configuration f 8 7 p f 5, consequently it is less affected by the monopole corrections

23 ( V pp ) made in KBG than the ground state, and moves down MeV in excitation energy. The states J = 5 +, J =, J = 7 +, J = 8 + and J = 9 + are also lowered by KBG, due to the occupation of the f 5 orbit. 5. Spectroscopy of 5 Sc The dominant configuration in the ground state of 5 Sc is (f 7) π (f 8 7 p ) ν. It produces the multiplet of states with J going from + to 5 +, below MeV. Experimentally little is known, the ground state +, one state without spin assignedatabout.6mevandfour statesseeninthebetadecayof 5 Ca[9] (see fig. 5). In the figure we include the calculated ground state multiplet and the high spins up to J = 8 +, as well as all the J = states present below.5mev. According to the KBG results, the state at.6mev should be the + member ofthegroundstatemultiplet. Inref.[5]it isarguedthatthisstate could be a, based on the FPD6 results. However this has to be attributed to the fact that FPD6 puts the f5 orbit too low at N=8. If this state were a, it would have been seen in the decay of 5 Ca, and this is not the case. The calculation produces six states inthe regionwhere experimentally only four have been found. We will discuss the decay pattern of 5 Ca in section , ,, + KB KBG Fig. 5. Energy levels of 5 Sc.

24 5. Spectroscopy of 5 Ti The experimental level scheme is compared in figure 6 with the full pf-shell results using the interactions KBG and FPD6. The location of the yrast states from J = + to J = corresponding to the dominant configuration (f 8 7 p ) ν (f 7 ) π is correct for both interactions. The lowest states of 5 Ti can be described essentially in terms of the coupling 5 Ti p. Notice the right position of the triplet J = +,, + at.mev. However, the doublet J =, + at.5mev is clearly better placed by KBG than by FPD6. 5 ( + ), +, + +, ( ) + + ( ) , , , , , , FPD6 KBG Fig. 6. Energy levels of 5 Ti. These states are dominantly f 8 7f 5 and they are shifted down by FPD6. This is even more evident for the second J = + state that is reasonably positioned by KBG but strongly shifted down by FPD6. The richness of the calculated level scheme calls for an improvement of the experimental spectrum. There are six experimental transitions known [], though with large uncertainties. The calculated values are in accord with the data, as can be seen in table.

25 Table Transitions in 5 Ti. B(M) (µ N ) (µ N ) >.6.5 B(E) (e fm ) (e fm ) >7 66 () 8 Th. 5. Spectroscopy of 5 V In figure 7 we have plotted all the states with parity and spin assignments up to MeV. Above this energy only the levels with spins equal or greater than are shown. The situation is similar to 5 Sc. The dominant configuration ) π (f 8 7 p (f 7 going from to 5 +. ) ν in the ground state produces the multiplet of states with J (9) + (6) () (5)+ +, KB KBG Fig. 7. Energy levels of 5 V. 5

26 The calculated ground state multiplet is more dilated that its experimental counterpart. Although the more relevant features are well accounted for we do not predict the correct ground state spin. The J = +, doublet around 85keV comes out inverted. These states can be described as the result of the coupling 5 V p. The high spin states J = 7 + and J = 9 + are shifted up, making the agreement with experiment worse than in the other nuclei studied in this work. Only two electromagnetic transitions are experimentally known []: the E connecting the states J = 9 + and J = 7 + (B(E) = 8(6)e fm ) and the M conecting the first excited J = + state and the J = + ground state (B(M) =.7()µ N ). For the first one the calculation produces a value B(E) = 87e fm and for the second B(M, + + ) =.7µ N, a rather impresive agreement. 5.5 Spectroscopy of 5 Cr Figure8showstheyraststatesof 5 Cr,uptothebandtermination,calculated in the full pf-shell and in a t = 5 truncation, compared with the experimental values. It is quite evident that the states obtained with a t = 5 truncation are nearly converged. The accord of the KBG results with the experiment is excellent KBG (t=5) KBG (full) Fig. 8. Yrast band of 5 Cr; experiment, t = 5 and full calculation. More detailed spectroscopic results are shown in figures 9 and in a t = 5 6

27 truncated calculation. In fig. 9 all the states experimentally known up to MeV [] have been included as well as the corresponding ones obtained with KBG and KB. The levels above this energy are plotted in figure, where the experimental data come from a recent experiment [5]. The agreement for the non yrast states is also very satisfactory. Notice the good positioning of the doublets J = +, at.7mev and J =, + at.5mev and the one to one correspondence between the three J = + and the three J = 5 + states measured in the excitation energy range MeV MeV and the calculated ones. Beyond MeV the yrare states and the odd yrast are also fairly well accounted for. It is clear from fig. that KBG produces a more satisfactory level scheme than KB, that shifts up all the levels with spins greater than J = 8 +, the maximum value attainable within the (f 7/ ) configuration. This is again a manifestation of the too large N=8 gap produced by this interaction. (5) KB KBG Fig. 9. Energy levels of 5 Cr up to MeV, t = 5. The experimental magnetic and quadrupole moments for the first excited state J = + are: and the calculated ones: µ exp ( + ) = +.(5)µ N Q exp ( + ) = 8.(6)efm µ( + ) = +.5µ N Q( + ) = 9.efm 7

28 + + (+) 9 (+) 8 7 (+) (+) (+) (9 + ) 8 (+) (8) + 7 (+) 6 (+) 8 (+) ( ) 5 (+) KB KBG + Fig.. Energy levels of 5 Cr above MeV, t = 5. In table the known experimental values for the electromagnetic transitions between the states of the yrast sequence are given. The correspondence with the predicted ones is good with the exception of the E transition J = to J = +, whose measured value is far outside the range expected in this mass region. Given its very large error bar, we would rather disregard this measure. Table Transitions in 5 Cr. B(M) (µ N ) (µ N ) (759). B(E) (e fm ) (e fm ) + + (6) (65) 7 59() () ().6 Th. 8

29 5.6 Spectroscopy of 5 Mn The full space calculation for the yrast states of the odd-odd nucleus 5 Mn produces only minor differences with the t = 5 results (fig. ). Notice the full correspondence between the states of the multiplet below MeV. The experimental data for the spins beyond the band termination ( to ) come from a recent experiment [5]. The agreement between experiment and theory is spectacular. There is only three experimental states with spin assignment not represented in the figure: a second J = (5 + ) at.mev, a third J = (5 + ) at.68mev and a second J = ( ) at.96mev. The t = 5 calculation places the second J = 5 + at.7mev, the third J = 5 + at.5mev and the second J = at.9mev ( ) ( ) + + (9 + ) (8 + ) KBG (t=5) KBG (full) Fig.. Yrast band of 5 Mn; experiment, t = 5 and full calculation. The experimental magnetic and quadrupole moments for the ground state are known: µ exp ( ) = +.6()µ N Q exp ( ) = +5(7)efm the calculated values reproduce nicely these values: µ( ) = +.958µ N Q( ) = +5efm 9

30 Some electromagnetic transitions along the yrast sequence have been measured. The calculations are in fair agreement with them (see table ). Table Transitions in 5 Mn. B(M) (µ N ) (µ N ) 7 +.5(56) > B(E) (e fm ) (e fm ) > > (6) 5 Th. 6 Half-lives and other β-decay properties Once the level schemes have been analysed, we study the β decays of these nuclei, to complete the description of this mass region. We present the results with the same ordering used for the level schemes. Because of the increasing sizes of the calculations we have limited this study to the isotopes of Calcium, Scandium, Titanium and Vanadium. The calculated values for KB and KBG gathered in the tables correspond to full hw calculation unless otherwise indicated. The Q β values are also included. The errors attached to the calculated values proceed from the errors in the experimental Q β values. We compute the half-lives by making the convolution of the strength function produced by the Lanczos method with the Fermi function, increasing the number of iterations until convergence is achieved. 6. β decays in the isobar chain A = 5 The experimental and calculated half-lives in the isobar chain A = 5 are given in table. The agreement is quite satisfactory.

31 Table Half-lives of A = 5 isobars. Q β values from []. T A J π KB KBG Q exp β (MeV) 5 Ca +... s.+.s.9(6)s.966(7).±.s 5 Sc s + s.5(5)s 6.888(6) s 5 Ti Stable 5 V th forbidden C.E. The Gamow-Teller strength function for 5 Ca is shown in fig.. The amount of strength below the Q β value is very small and it is concentrated at MeV. Notice that there is a large amount of strength close above the Q β threshold. Thus, small variations in the excitation energies of the daughter states could influence (although not drastically) the final value of the half-life. In 5 Sc (see figure ) the situation is similar. Only a tiny fraction of the total strength lies below the Q β value <στ> Fig.. Gamow-Teller strength after 5 iterations for the 5 Ca β decay. The < στ > values are summed up in MeV bins. The dashed line indicates the experimental Q β value. We have computed the intensity (in percentage) of the ground state decay to the levels within the Q β window of the daughter nucleus. 5 Ca decays 99.7% to theexcited statej = at.78mev in 5 Sc, what agreesperfectly withthe experimental situation (99.()% decay to the J = state at.85mev).

32 6 5 <στ> Fig.. Gamow-Teller strength after 5 iterations for the 5 Sc β decay. The < στ > values are summed up in MeV bins. The dashed line indicates the experimental Q β value. Analogously, figure shows the experimental and calculated intensities for the β decay of 5 Sc. Both experimental peaks are reproduced by the calculation although they are slightly shifted. Since the 5 Ca decay feeds an unique state and 5 Sc decays to just two, it is possible to use the experimental energies of the states to calculate the half-lives and to gauge the effect of the phase space. Proceeding in such a way, the new calculated half-lives for 5 Ca and 5 Sc are, respectively, T =. +.s and = s which improve T the agreement with the experimental half-lives. These results are also shown in table. 9 8 Teor. 7 6 I β (%) Fig.. Percentages of the β decay of 5 Sc. Experimental data from [].

33 6. β decays in the isobar chain A = 5 Table 5 gathers the experimental and calculated half-lives for the isobar chain A = 5. The calculated values are quite close to the experimental ones. Table 5 Half-lives of A = 5 isobars. Q β values from []. T A J π KB KBG Q exp β (MeV) 5 Ca 5 7 Sc 5 Ti s s.(8)s 7.(9)..s..s.()s 6.58().±.s 6.9.min 8..min 5.76()min.76(5) 6.97±.min 5 V Stable Although within the Q β energy window there is a large amount of Gamow- Tellerstrength(seefigures5,6and7),thesituationlookslikeintheA = 5 chain and much strenght is located around the Q β values <στ> Fig. 5. Gamow-Teller strength after 5 iterations for the 5 Ca β decay. The < στ > values are summed up in MeV bins. The dashed line indicates the experimental Q β value. In figures 8 and 9 we compare the experimental percentages for the decays

34 <στ> Fig. 6. Gamow-Teller strength after 5 iterations for the 5 Sc β decay. The < στ > values are summed up in MeV bins. The dashed line indicates the experimental Q β value <στ> Fig. 7. Gamow-Teller strength after 5 iterations for the 5 Ti β decay. The < στ > values are summed up in MeV bins. The dashed line indicates the experimental Q β value. with the calculation for both 5 Ca and 5 Sc. In the latter, the agreement is very good. 5 Ca demands some comments. To compare with the experimental data we must keep in mind that only states up to.8mev in 5 Sc have been observed. Inthefirst MeV thereisasingle experimental statefedbythedecay. This fact is well reproduced by the calculation. Up to.5mev there are seven states experimentally observed in the decay, the same number is given by the calculation. However, the distribution of the intensity is not the same. In the calculation all the intensity is concentrated between.5 MeV and. MeV, while experimentally there are two states strongly fed at lower energy.

35 Teor. I β (%) Fig. 8. Percentages of the β decay of 5 Ca. Experimental data from []. Teor. I β (%) 5 6 Fig. 9. Percentages of the β decay of 5 Sc. Experimental data from []. In the case of 5 Ti the decay only feeds two states in the daughter nucleus, a J = at.98mev (8.()%) and a J = 5 at.mev (9.9()%). The calculation agrees with it since it predicts the decay to a J = state at.55mev (6.5%) and to another J = 5 state at.57mev (9.5%). In view of the comparison of the intensities for the decays of 5 Ca, 5 Sc and 5 Ti, it seems plausible to make a new analysis in the same terms we did in the previous section, i.e. to take the experimental energies for the four states strongly fed in 5 Sc and the corresponding two in 5 Ti. The new half-lives we obtain are T =. +.s for 5 Sc and T = min for 5 Ti, in much better agreement with the experimental values shown in table 5. For 5 Ca this analysis is not so easy since the number of states fed is larger and it is 5

36 not simple to make a one to one correspondence between experimental and theoretical states. 6. β decays in the isobar chain A = 5 In table 6 the experimental data and the calculated results for the half-lives of the nuclei in the A = 5 isobar chain are shown. For 5 V we have only made a t = calculation, because, to be consistent, the calculation in the daugther has to be t = 5. Going beyond that would have demanded an enormous computational effort. The theoretical half-lives compare reasonably well with the experimental numbers, although the discrepancy in 5 Ca is severe. We shall now analyse the intensity distributions and try to recompute the halflives using as much as possible the experimental excitation energies. Table 6 Half-lives of A = 5 isobars. Q β values from []. T A J π KB KBG Q exp β (MeV) 5 Ca Sc s s s s.6()s 7.9(5) s 8.()s 9.(6) 5 Ti min.8 +.7min.7()min.97(8).±.min 5 V min min.7(5)min.9756().87±.min Figure shows the percentages for the 5 Ca β decay and figure the same for 5 Sc. The 5 Ca case demands a special attention. When we analyzed the 5 Sc spectrum in section 5. we mentioned that there was one state without any experimental spin assignment at 675 kev (fig. 5). In [5] the calculation was made with the FPD6 interaction, producing a first J = at 7keV. It was suggested that it could correspond to the 675 kev state. Our interpretation of this state as the + member of the ground state multiplet is experimentally supported by the fact that if this state were a J =, it would have been observed in the β decay of 5 Ca. This is not the case, as can been seen in figure. On the contrary, most of the decay (86.8%) goes to the first experimental state with J = assignment, corresponding to the one we obtain in the calculation. Thus, it is discarded that this unassigned experimental state 6

37 9 8 Teor. 7 6 I β (%) Fig.. Percentages of the β decay of 5 Ca. Experimental data from []. Teor. I β (%) 5 6 Fig.. Percentages of the β decay of 5 Sc. Experimental data from []. could be a J =. The first J = corresponds to a configuration with a neutron in the f5 orbit. As we said in the introduction, FPD6 locates the f5 orbit too low, producing a very low J = state. The other three experimental J = states are also fed through the decay of 5 Ca although with smaller intensity. The percentage of the second at.7mev is.% while for the third at.5mev it is.6% and.% for the last one at.mev. These data help us to make a one to one correspondence between the states above MeV in figure 5. Our third state J = at.6mev corresponds to the second experimental one, our fourth J = at.mev to the third and the fourth experimental J = to our sixth at.mev. The situation for 5 Sc in figure is very satisfactory. The decay to the J = + excited state at MeV and to the triplet around.mev is well reproduced, 7

38 although in the calculation the triplet is more expanded than the experimental one. There is one state at.mev without assigned spin, that seems to correspond to the + state near MeV (the calculation predicts a.7% branch to it). The next branch goes to the second at.mev. Beyond, it becomes difficult to establish a detailed correspondence. The decay of 5 Ti feeds (%) a single state in the daughter nucleus 5 V, the J = at kev. This fact is reproduced properly by the calculation which predicts a % feeding of the J = state at kev in the daughter nucleus for the β decay of 5 Ti. On the other hand, 5 V decays 99.(5)% to the J = + excited state of 5 Cr located at.mev, what agrees with the calculation (a 98.9% feeds the J = + state at.5mev). Finally, we recalculate the half-lives using the experimental energies of the single state fed in the decay of 5 Ti and 5 V, and of the two states fed in the decay of 5 Ca. For 5 Sc the analysis in these terms is not so easy. The new calculated half-lives are T =. +.min for 5 Ti and T =.87 +.min for 5 V. Both results improve the previous ones shown in table 6. For 5 Ca we get a half-life T = s, reducing substantially the initial discrepancy. 7 Coulomb energy differences (CED) in the mirror pair 5 Mn- 5 Fe It is well known that if the nucleon-nucleon interaction were charge symmetric the mirror nuclei would have the same level scheme. Thus, the small differences (normally a few tens of kev) in their excitation energies, are due to the isospin symmetry breaking Coulomb interaction. The analysis of this effect has been recentlypushedupinmass, withthestudyofthea=7anda=9mirrorpairs [5,5]. It was found that the CED s are extremely sensitive to the structure of the nuclear wave functions. The next step has been to measure the 5 Mn- 5 Fe mirror pair, the heaviest one in which high spin states up to the band termination have been observed [8]. In figure the calculated spectra for both nuclei, including the Coulomb interaction are compared with the experimental data. The Coulomb matrix elements are those denoted A in ref [5]. At this scale, very slight differences (hardly any) are appreciable between them. However, looking more carefully at the CED s an abrupt change is observed at J = 7 (fig. ). The effect is also present in the calculation, that shows exactly the same trends as the experiment, although with enhanced values. The large increase in the CED can be interpreted as due to the alignment of one proton pair in 5 Fe which does not occur in 5 Mn. As a consequence, the Coulomb energy is sharply reduced in 5 Fe but not in 5 Mn, therefore the CDE increases dramatically (in absolute value). The breaking of a proton pair occurs in 5 Fe because it has 8