Beta-decay properties of neutron-rich Ge, Se, Kr, Sr, Ru, and Pd isotopes from deformed quasiparticle random-phase approximation

Beta-decay properties of neutron-rich Ge, Se, Kr, Sr, Ru, and Pd isotopes from deformed quasiparticle random-phase approximation


-decay properties of even and odd- neutron-rich Ge, Se, Kr, Sr, Ru, and Pd isotopes involved in the astrophysical rapid neutron capture process are studied within a deformed proton-neutron quasiparticle random-phase approximation. The underlying mean field is described self-consistently from deformed Skyrme Hartree-Fock calculations with pairing correlations. Residual interactions in the particle-hole and particle-particle channels are also included in the formalism. The isotopic evolution of the various nuclear equilibrium shapes and the corresponding charge radii are investigated in all the isotopic chains. The energy distributions of the Gamow-Teller strength as well as the -decay half-lives are discussed and compared with the available experimental information. It is shown that nuclear deformation plays a significant role in the description of the decay properties in this mass region. Reliable predictions of the strength distributions are essential to evaluate decay rates in astrophysical scenarios.

21.60.Jz, 23.40.Hc, 27.60.+j, 26.30.-k

I Introduction

The rapid structural changes occurring in the ground state and low-lying collective excited states of neutron-rich nuclei in the mass region have been extensively studied both theoretical and experimentally (see e.g. (1); (2) and references therein). From the theoretical side, the equilibrium nuclear shapes in this mass region have been shown to suffer rapid changes as a function of the number of nucleons with competing spherical, axially symmetric prolate and oblate, and triaxial shapes at close energies. Both relativistic (3); (4) and nonrelativistic (5); (6); (7); (8); (9) approaches agree in the general description of the nuclear structural evolution in this mass region, which is supported experimentally by spectroscopic studies (11); (12); (10), lifetime measurements (13); (14); (15) and quadrupole moments for rotational bands (15), as well as by laser spectroscopy measurements (16).

However, the nuclear structure richness is not the only attractive feature characterizing these nuclei. Another remarkable property of nuclei in this mass region is that they are involved in the astrophysical rapid neutron capture process (r process), which is considered to be one of the main nucleosynthesis mechanisms leading to the production of heavy neutron-rich nuclei in the universe (17); (18). The r-process nucleosynthesis involves many neutron-rich unstable isotopes, whose neutron capture rates, masses, and -decay half-lives () are crucial quantities to understand the possible r-process paths, the isotopic abundances, and the time scales of the process (18); (19); (20). Although much progress has been done measuring masses (see for example the Jyväskylä mass database (21)) and half-lives (22); (23); (24), most of the nuclear properties of relevance for the r process are experimentally unknown due to their extremely low production yields in the laboratory. Therefore, reliable nuclear physics models are required to simulate properly the r process.

The quasiparticle random-phase approximation (QRPA) is considered a well suited model to describe medium-mass open-shell nuclear properties and specifically -decay properties. QRPA calculations for neutron-rich nuclei have been carried out within different spherical formalisms, such as Hartree-Fock-Bogoliubov (HFB) (25), continuum QRPA with density functionals (26), and relativistic mean field approaches (27). However, the mass region we are dealing with requires nuclear deformation as a relevant degree of freedom to characterize the nuclear structure involved in the calculation of the -strength functions. The deformed QRPA formalism was developed in Refs. (28); (29); (30); (31), using phenomenological mean fields. A Tamm-Dancoff approximation with Sk3 interaction was also implemented in Ref. (32). More recently, deformed QRPA calculations using deformed Woods-Saxon potentials and realistic CD-Bonn residual forces have been performed in (33); (34). First-forbidden transitions were also considered in those references, showing that their effect in this mass region can be neglected. Various self-consistent deformed QRPA calculations to describe the -decay properties, either with Skyrme (35) or Gogny (36) interactions are also available in the literature.

In Refs. (37); (38) the decay properties of neutron-rich Zr and Mo isotopes were studied within a deformed proton-neutron QRPA based on a self-consistent Hartree-Fock (HF) mean field formalism with Skyrme interactions and pairing correlations in BCS approximation. Residual spin-isospin interactions were also included in the particle-hole and particle-particle channels (39); (40). In this work this study is extended to the neighboring regions including even and odd- neutron-rich Ge, Se, Kr, Sr, Ru, and Pd isotopes. These calculations are timely because they address a mass region which is at the borderline of present experimental capabilities for measuring half-lives at MSU and RIKEN (22); (23); (24). In addition, theoretical calculations can be tested with the available experimental information on half-lives providing simultaneously predictions for the underlying Gamow-Teller strength distributions and for the half-lives of more exotic nuclei not yet measured. Finally, this more comprehensive study allows one to judge better the extent to which the method is able to describe the decay properties of nuclei in a wider mass region that includes spherical, well deformed, and weakly deformed transitional isotopes, as well as isotopes exhibiting shape coexistence. Therefore, the theoretical method will be tested over a rich set of different nuclear structures that will reveal the limitations of the model.

The paper is organized as follows. In Sec. II a review of the theoretical formalism used is introduced. Section III contains the results obtained for the potential energy curves (PEC), Gamow-Teller (GT) strength distributions, and -decay half-lives, which are compared with the experimental data. Section IV summarizes the main conclusions.

Ii Theoretical Formalism

A summary of the theoretical framework used in this paper to describe the -decay properties in neutron-rich isotopes is shown in this section. More details of the formalism can be found elsewhere (39); (40). The method starts from a self-consistent calculation of the mean field by means of a deformed Skyrme Hartree-Fock procedure with pairing correlations in BCS approximation. Single-particle energies, wave functions, and occupation amplitudes are generated from this mean field. The Skyrme interaction SLy4 (41) is used as a representative of modern Skyrme forces. It has been very successful describing nuclear properties all along the nuclear chart and has been extensively studied (6); (7); (42).

The solution of the HF equation, assuming time reversal and axial symmetry, is found by using the formalism developed in Ref. (43). The single-particle wave functions are expanded in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates, using twelve major shells. The pairing gap parameters for protons and neutrons in the BCS approximation are determined phenomenologically from the odd-even mass differences (44). In a further step, constrained HF calculations with a quadratic constraint are performed to construct the PECs, analyzing the nuclear binding energies in terms of the quadrupole deformation parameter . Calculations for GT strengths are performed subsequently for the various minima in the energy curves indicating the equilibrium shapes of each nucleus. Since decays connecting different shapes are disfavored, similar shapes are assumed for the ground state of the parent nucleus and for all populated states in the daughter nucleus. The validity of this assumption was discussed for example in Refs. (28); (30).

To describe GT transitions, a separable spin-isospin residual interaction in the particle-hole (ph) and particle-particle (pp) channels is added to the mean field and treated in a deformed proton-neutron QRPA (28); (29); (30); (31); (32); (39); (40); (45). An optimum set of coupling strengths could be chosen following a case by case fitting procedure and one will finally get different answers depending on the nucleus, shape, and Skyrme force. However, since the purpose here is to test the ability of QRPA to account for the GT strength distributions in this mass region with as few free parameters as possible, the same coupling strengths are used for all the nuclei considered in this paper, which are taken from previous works (37); (38). We use MeV and MeV for the residual interaction in the ph and pp channels, respectively.

The sensitivity of the GT strength distributions to the various ingredients contributing to the deformed QRPA calculations, namely to the nucleon-nucleon effective force, to deformation, to pairing correlations, and to residual interactions, have been investigated in the past (39); (40); (46); (47); (48). In this work the most reasonable choices found in those references are used. Summarizing the various sensitivities, the conclusion is that the main features of the GT strength distributions are in general very robust against the Skyrme force used, showing some more sensitivity in the spherical cases, where the location of the single-particle energies is more critical to determine the excitation energies of the GT transitions. Deformation has been shown to be an important issue to describe the profiles of the GT strength distributions. First, because the degeneracy of the spherical shells is broken making the GT strength distributions more fragmented than the corresponding spherical ones. Secondly, because the energy levels of deformed orbitals cross each other in a way that depends on the magnitude of the quadrupole deformation as well as on the oblate or prolate character. This level crossing may lead in some instances to sizable differences in the GT profiles, a fact that has been exploited to learn about the nuclear shape from the measured -decay pattern (49); (50); (51). Pairing correlations are also important to describe nuclei out of closed shells. Their influence on the GT profiles was studied in Ref. (40), concluding that the main effect is to decrease slightly the strength at low energies and to create new transitions, mainly at high energies, that are forbidden in the absence of such correlations. The effect of the ph and pp residual interactions is also well known. The repulsive ph interaction redistributes the GT strength by shifting it to higher excitation energies causing a displacement of the GT resonance. It also reduces somewhat the total strength. The attractive pp interaction moves the strength to lower energies. Its effect on the GT resonance is in general negligible, but nevertheless, the changes induced in the low-energy region are of great relevance in the calculation of the -decay half-lives, which are only sensitive to the strength contained in the energy region below the -energy window.

The GT transition amplitudes in the intrinsic frame connecting the ground state of an even-even nucleus to one phonon states in the daughter nucleus are found to be






in terms of the occupation amplitudes for neutrons and protons () and the matrix elements of the spin operator, , connecting proton and neutron single-particle states, as they come out from the HF+BCS calculation. and are the forward and backward amplitudes of the QRPA phonon operator, respectively.

Once the intrinsic amplitudes in Eq. (1) are calculated, the GT strength in the laboratory system for a transition can be obtained as


in units. To obtain this expression, the initial and final states in the laboratory frame have been expressed in terms of the intrinsic states using the Bohr-Mottelson factorization (52).

The specific treatment of odd- systems has been described (31); (47) by considering two types of GT contributions. One type is due to phonon excitations in which the odd nucleon acts only as a spectator. The transition amplitudes in the intrinsic frame are in this case basically the same as in the even-even case, but with the blocked spectator excluded from the calculation. The other type of transitions involves the odd nucleon and is treated perturbatively by taking into account phonon correlations to first order in the quasiparticle transitions. The excitation energies of the GT states with respect to the ground state in the daughter nuclei have been discussed in Ref. (47) for both types of transitions in terms of the QRPA phonon energy and the quasiparticle energies.

Figure 1: Potential energy curves for even-even neutron-rich Ge, Se, Kr, Sr, Zr, Mo, Ru, and Pd isotopes obtained from constrained HF+BCS calculations with the Skyrme force SLy4.
Figure 2: (Color online) Isotopic evolution of the quadrupole deformation parameter (a) and charge radius (b) corresponding to the energy minima obtained from the Skyrme interaction SLy4 for Ge isotopes. Ground state results are encircled.
Figure 3: (Color online) Same as in Fig. 2, but for Se isotopes.
Figure 4: (Color online) Same as in Fig. 2, but for Kr isotopes. Experimental charge radii are from (59).
Figure 5: (Color online) Same as in Fig. 2, but for Sr isotopes. Experimental charge radii are from (59).
Figure 6: (Color online) Same as in Fig. 2, but for Ru isotopes.
Figure 7: (Color online) Same as in Fig. 2, but for Pd isotopes.
Figure 8: (Color online) QRPA-SLy4 Gamow-Teller strength distributions for Ge isotopes as a function of the excitation energy in the daughter nucleus. The calculations correspond to the various equilibrium deformations found in the PECs.
Figure 9: (Color online) Same as in Fig. 8, but for Se isotopes.
Figure 10: (Color online) Same as in Fig. 8, but for Kr isotopes.
Figure 11: (Color online) Same as in Fig. 8, but for Sr isotopes.
Figure 12: (Color online) Same as in Fig. 8, but for Ru isotopes.
Figure 13: (Color online) Same as in Fig. 8, but for Pd isotopes.

The -decay half-life is obtained by summing all the allowed transition strengths to states in the daughter nucleus with excitation energies lying below the corresponding -energy, , written in terms of the nuclear masses and the electron mass (), and weighted with the phase space factors ,


with  s and , where 0.77 is a standard quenching factor. The same quenching factor is included in all the figures shown later for the GT strength distributions. The bare results can be recovered by scaling the results in this paper for and with the square of this quenching factor.

The Fermi integral is computed numerically for each value of the energy including screening and finite size effects, as explained in Ref. (53),




where ; ; is the fine structure constant and the nuclear radius. is the total energy of the particle, is the total energy available in units, and is the momentum in units.

Iii Results and discussion

In this section I present first the PECs in the isotopic chains studied. Quadrupole deformation parameters as well as charge r.m.s. radii () are analyzed as a function of the mass number. Then, energy distributions of the GT strength corresponding to the local minima in the PECs are calculated. Finally, half-lives are evaluated and compared with the experiment.

iii.1 Structural isotopic evolution

In Fig. 1 the PECs, i.e., the energies relative to that of the ground state, are plotted as a function of the quadrupole deformation for the neutron-rich Ge, Se, Kr, Sr, Zr, Mo, Ru, and Pd isotopes. The results correspond to the SLy4 interaction. The isotopes covered in this study include middle-shell nuclei with proton numbers between shell closures and , namely and neutron numbers between shell closures (as in Ge) and (as in the heaviest Pd).

In most of the isotopic chains one can see the appearance of several equilibrium nuclear shapes, whose relative energies change with the number of neutrons. In Ge isotopes, prolate shapes that are ground states in the lighter isotopes are found with the only exception of Ge, where a spherical shape is found in accordance with its semi-magic character. At (Ge) oblate and prolate shapes are practically degenerate in energy and oblate shapes become ground states for heavier isotopes. The case of Se isotopes is similar with oblate and prolate minima all along the isotopic chain. The lighter (heavier) isotopes have prolate (oblate) ground states with transitional isotopes around (Se). In this case the energy barriers are more pronounced than in the case of Ge isotopes. Kr isotopes show competing shapes in the lighter isotopes that become oblate at (Kr) and then turn into prolate shapes beyond Kr. In the heavier isotopes, as in the case of Se isotopes, shape coexistence is found with very well developed oblate and prolate minima separated with high energy barriers. Sr isotopes show a transition from oblate at (Sr) to prolate at (Sr) with a two minima structure for heavier isotopes. The cases of Zr and Mo isotopes were discussed in Refs. (37); (38). Both oblate and prolate minima are observed in the lighter isotopes of Zr and Mo with prolate ground states. Whereas the prolate shape remains ground state in most of the heavier Zr isotopes, oblate shapes are lower in energy for the heavier Mo isotopes. Finally, Ru and Pd isotopes show oblate and prolate minima in the lighter isotopes and a gradual transition into spherical shapes as one approaches the shell closure at .

In summary, a large diversity of nuclear structures are found in this mass region, from spherical to well deformed shapes, passing through soft transitional nuclei and even possible shape-coexistence structures. This rich variety of shapes represents a challenge to any theoretical model trying to describe them in a unified manner. In the next subsections the results are compared with the available experimental data, which are restricted at present to -decay half-lives. Then, the theoretical approach will be tested against this information and the limitations of the model will be established.

These results are in qualitative agreement with similar calculations obtained in this mass region from different theoretical approaches, including macroscopic-microscopic methods based on liquid drop models with shell corrections (54); (55), relativistic mean fields (56), as well as nonrelativistic calculations with Skyrme (5) and Gogny (8); (9); (57) interactions. Thus, a consistent theoretical description emerges, which is supported by the still scarce experimental information available (2); (11); (12); (10); (13); (14); (15); (16); (58).

The isotopic evolution can be better appreciated in Figs. 27, where quadrupole deformations (a) and r.m.s. charge radii (b) of the various energy minima are plotted as a function of the mass number . The deformation corresponding to the ground state for each isotope is encircled. Also shown in these figures for odd- isotopes, are the spin and parity of the different shapes and the experimental assignments (44). The experimental assignments based on systematics estimated from trends in neighboring nuclides have not been included.

In Fig. 2 for Ge isotopes one can see clearly the shape transition at from prolate shapes with to oblate shapes with . Charge r.m.s. radii have not been measured in these isotopes, but it is expected from these calculations a very smooth behavior given that the magnitude of in the prolate and oblate sectors are very similar. Fig. 3 shows the analogous results for Se isotopes. In this case one can see the transition from prolate () to oblate () at . The prolate shape grows in the heavier isotopes (), but they are never ground states and then, the expected increasing in the charge radii is smooth. Kr isotopes in Fig. 4 show first a shape transition from prolate () to oblate () and a subsequent transition from oblate to prolate () shapes. The radii are sensitive to this transitions, although the measured radii (59) seem to favored prolate shapes in the lighter isotopes. Sr isotopes in Fig. 5 show a clear transition from oblate to strong prolate () deformations at . This shape transition is well correlated with the trend change observed in the charge radii that shows a sizable jump between Sr and Sr both theoretical and experimentally (59). In the case of Ru (Pd) isotopes shown in Fig. 6 (7), one can see a smooth transition from deformed oblate (prolate) solutions in the lighter isotopes to spherical shapes in the heavier ones. This change is felt in the trends of the radii, but no experiments are yet available to compare with.

Spins and parities in odd- isotopes can be compared with their experimental assignments. In the Ge isotopes the calculations agree reasonably well with the assignments taking into account that oblate and prolate shapes are very close in energy and that a isomer is observed experimentally in Ge at 248 keV. In the lighter Se isotopes, and states are obtained, whereas experimental assignments are . In both isotopes, states very close in energy to the ground states are also obtained, although somewhat above. Similarly, in the lighter Kr isotopes the experimental assignments are obtained very close in energy to the ground states, although slightly above. On the other hand, a isomer is experimentally observed in Kr at 355 keV that corresponds to the ground state here. Sr isotopes exhibit a nice agreement. The measured spin and parities of ground states in Sr correspond to the prolate calculations. A state is also observed experimentally in Sr at 56 keV. In Sr the observed ground state appears as an excited state. It is also worth noting that the prolate ground state for this isotope is observed experimentally at 645 keV. In the case of Ru isotopes the measured are difficult to reproduce. They are found in the calculations, but not as ground states. On the other hand, the negative parity states found in the calculations are also seen experimentally at low energies. In particular, an isomeric state at an undetermined energy has been seen in Ru. Finally, in Pd isotopes the negative parity isomers, which are oblate in this description, are reproduced in the calculations, but not the ground states.

Figure 14: (Color online) QRPA-SLy4 accumulated GT strengths in Ge isotopes calculated for the various equilibrium shapes. and energies are shown by solid and dashed vertical arrows, respectively.
Figure 15: (Color online) Same as in Fig. 14, but for Se isotopes.
Figure 16: (Color online) Same as in Fig. 14, but for Kr isotopes.
Figure 17: (Color online) Same as in Fig. 14, but for Sr isotopes.
Figure 18: (Color online) Same as in Fig. 14, but for Ru isotopes.
Figure 19: (Color online) Same as in Fig. 14, but for Pd isotopes.
Figure 20: (Color online) Measured -decay half-lives for Ge isotopes compared to theoretical QRPA-SLy4 results calculated from different shapes. Circles are experimental values (open circles are experimental values from systematics) (44).
Figure 21: (Color online) Same as in Fig. 20, but for Se isotopes.
Figure 22: (Color online) Same as in Fig. 20, but for Kr isotopes. Experimental half-lives are from (44); (23).
Figure 23: (Color online) Same as in Fig. 20, but for Sr isotopes. Experimental half-lives are from (44); (23).
Figure 24: (Color online) Same as in Fig. 20, but for Ru isotopes.
Figure 25: (Color online) Same as in Fig. 20, but for Pd isotopes.
Figure 26: (Color online) Ratio of calculated to experimental -decay half-lives for two sets of calculations, with the spherical approximation (open dots) and with the deformation that corresponds to the minimum of the PECs (solid dots). The ratios are plotted as a function of the experimental half-lives (a) and as a function of the quadrupole deformation at the minimum of the PECs (b).

iii.2 Gamow-Teller strength distributions

In the next figures, the energy distributions of the GT strength corresponding to the various deformed equilibrium shapes are shown for each isotopic chain. The results are obtained from QRPA with the force SLy4 with pairing correlations and with residual interactions with the parameters written in Sec. II. The GT strength is plotted versus the excitation energy of the daughter nucleus with a quenching factor 0.77. Zr and Mo isotopes were already studied in Refs. (37); (38) and are not repeated here.

Figures 813 contain the results for Ge, Se, Kr, Sr, Ru, and Pd isotopes. The energy distributions of the individual GT strengths corresponding to the ground state shapes are shown, together with continuous distributions for the ground state shapes as well as for the other possible shapes, obtained by folding the strength with 1 MeV width Breit-Wigner functions. values are shown with vertical arrows. In both cases, even and odd isotopes, the values increase with the number of neutrons in each isotopic chain and the values in the odd- isotopes are about 2-3 MeV larger than the values in the neighbor even-even isotopes . The general structure of the GT distributions is characterized by the existence of a GT resonance, which is placed at increasing excitation energy as the number of neutrons increases in a given isotopic chain. The total GT strength also increases with , as it is expected to fulfill the Ikeda sum rule. The various shapes produce quite similar GT strength distributions on a global scale. Nevertheless, the small differences among the various shapes at the low energy tails (below the ) of the GT strength distributions that can be appreciated because of the logarithmic scale, lead to sizable effects in the -decay half-lives.

Unfortunately, comparison with experiment is still not possible for the GT strength distributions, the measured half-lives will be compared to the calculations in the next subsection. Comparison with calculated GT distributions from other theoretical approaches is also restricted to the few cases where these results have been published (33); (35). In Refs. (33); (60) the authors performed QRPA calculations with deformed Woods-Saxon potentials and realistic CD-Bonn residual forces using the -matrix formalism and compared these results with the results obtained from separable forces. While in Ref. (33) the comparison between the results obtained from realistic or separable residual interactions is restricted to the half-lives, in Ref. (60) the authors compared those results in the context of two-neutrino double-beta decay, concluding that both approaches, realistic and separable, lead to similar results. On the other hand, in Ref. (35) the Skyrme force SLy4 was used to generate the mean field as it is done in this work. The residual interaction in the ph channel was self-consistently introduced and not reduced to a separable form. Finally the pp residual interaction was written as a contact force with a coupling strength fitted to reproduce the half-life in Zr. The GT strength distributions in neutron-rich Zr isotopes obtained from this approach were compared with the corresponding distributions obtained with separable forces in Figs. 5-6 in Ref. (35). From this comparison one can conclude that in many aspects the main characteristics of the consistent force are maintained by a separable force with a much lower computational cost. The comparison of the half-lives shows also a remarkable agreement between both approaches.

In the next figures, Figs. 1419, one can see in more detail the accumulated GT strength in the energy region below the corresponding energy of each isotope, which is the relevant energy range for the calculation of the half-lives. The vertical solid (dashed) arrows show the () energies, taken from experiment (44). In these figures the sensitivity of these distributions to deformation can be appreciated and one can understand that measurements of the GT strength distribution from -decay can be, in particular cases, an additional source of information about the nuclear deformation, as it was shown in Refs. (49); (50); (51). The GT strength distribution in odd- isotopes is found to be displaced to higher energies (typically about 2-3 MeV) with respect to the even-even case. The shift corresponds roughly to the breaking of a neutron pair and therefore it amounts to about twice the neutron pairing gap. Below this energy only transitions involving the odd nucleon are possible.

The energy distribution of the GT strength is fundamental to constrain the underlying nuclear structure. For a theoretical model, it represents a more demanding test than just reproducing half-lives or total GT strengths that are integral quantities obtained from these strength distributions properly weighted with phase factors (see Eq. 6). These quantities might be reproduced even with wrong strength distributions. This is of especial importance in astrophysical scenarios of high densities and temperatures that cannot be reproduced in the laboratory. Given that the phase factors in the stellar medium are different from those in the laboratory, the stellar half-lives become dependent on the electron distribution in the stellar plasma that eventually may block the -particle emission (61). Therefore, to describe properly the decay rates under extreme conditions of density and temperature, it is not sufficient to reproduce the half-lives in the laboratory. One needs, in addition, to have a reliable description of the GT strength distributions (62); (63).

iii.3 Beta-decay Half-lives

The calculation of the half-lives in Eq. (6) involves knowledge of the GT strength distribution and of the energies (), which are evaluated by using values obtained from the mass differences between parent and daughter nuclei obtained from SLy4 with a zero-range pairing force and Lipkin-Nogami obtained from the code HFBTHO (64).

In Figs. 2025 the measured -decay half-lives (solid dots, open dots stand for experimental values from systematics) (44); (23) are compared with the theoretical results obtained with the prolate, oblate, and spherical equilibrium shapes, for the various isotopic chains. In Fig. 20 one can see the half-lives for Ge isotopes. The lighter isotopes are not well reproduced, being largely overestimated. This point will be discussed later. The half-lives obtained from oblate shapes are larger than the corresponding prolate ones. This feature is correlated with the GT strength contained below the energy in Figs. 8 and 14. Prolate shapes, which are closer to experiment, are also the ground states in this range of masses according to the calculations (see Fig. 2). For heavier isotopes, the half-lives for oblate and prolate shapes are very similar. In the case of Se isotopes in Fig. 21, the calculations also overestimate the half-lives of the lighter isotopes, but the agreement with experiment is in this case much better. In the middle region the experimental half-lives, which are taken from systematics, are reasonably well reproduced. The half-lives of heavier isotopes exhibit a rather flat behavior. Half-lives of Kr isotopes are shown in Fig. 22. As in the previous figures, the half-lives from the oblate shapes are larger than the prolate ones in the lighter Kr isotopes, but the situation is reversed at Kr. This is again nicely correlated with the GT strength at low excitation energies shown in Fig. 16. In general, the half-lives in the middle region are well described. This is also true for Sr and Ru isotopes in Figs. 23 and 24, respectively, where the trends observed experimentally are well reproduced, except for the lighter Sr isotopes that are clearly underestimated and the heavier Ru isotopes, where the data from systematics fall down faster than the calculations. Finally, in the case of Pd isotopes, shown in Fig. 25, the calculations underestimate (overestimate) the measured half-lives in the lighter (heavier) isotopes.

All in all, the agreement with experiment is reasonable, especially in the middle regions. These regions contain in general well deformed nuclei, where the present approach is more suitable. On the other hand, weakly deformed transitional isotopes, such as light Ge and Se isotopes and heavy Ru and Pd isotopes are not so well described. Furthermore, in the light isotopes of all the isotopic chains, which are closer to the valley of stability, the half-lives are larger because of the small energies involved. In these cases the half-lives are determined exclusively by the very low energy tail of the GT strength distribution contained in the narrow window below . Therefore, tiny variations in the description of the GT strength distribution in the low-lying energy region can drive sizable effects in the half-lives. Of course it is also important to describe the half-lives of the long-lived isotopes, but their significance to constrain the GT strength distribution is minor since the half-lives are insensitive to most of this distribution.

Half-lives for neutron-rich Kr, Sr, Zr, and Mo isotopes calculated from self-consistent deformed QRPA calculations with the Gogny D1M interaction and experimental values of (36) agree with the results in this work within the uncertainties of the calculations. The agreement is also very reasonable between the calculated half-lives and those obtained from deformed QRPA calculations using deformed Woods-Saxon potentials to generate the mean field and complemented with realistic CD-Bonn residual forces (33); (34). The agreement is also good with the results in Ref. (35) using the Skyrme force SLy4 with consistent residual interactions in the ph channel as mentioned earlier. Fig. 7 in that reference displays this comparison.

It is also worth noticing that the worst agreement with experiment occurs in the light Ge isotopes, as well as in heavy Pd isotopes. In these cases the calculations overestimate the experiment leaving room for contributions coming from first forbidden (FF) transitions. One can understand from simple qualitative arguments that the role of FF transitions is expected to be more important in lighter Ge and in Pd isotopes. Thus, for Ge, Se, Kr, and Sr isotopes, the last occupied proton orbitals come basically from the and negative-parity spherical shells. On the other hand the neutrons in Ge isotopes occupy orbitals belonging to the and positive-parity spherical shells. Therefore, in the -decay one neutron in a positive-parity state is transformed into a proton that would sit in a negative-parity state, thus suppressing GT and favoring FF transitions in the low-lying transitions. This is particularly true for the lighter Ge isotopes. In the heavier ones, other neutron states with negative parity () have to be considered because of deformation effects. The same argument can be applied to the lighter Se, Kr, and Sr isotopes, but in these cases proton states from positive parity () are closer in energy and would participate in the decay favoring GT transitions. The situation is different in the case of Ru and Pd isotopes. Now the available proton states for the decay are of positive parity (), whereas most of the last occupied neutrons belong to negative-parity states (), thus favoring FF transitions. According to calculations (26); (33) of the FF transitions in this mass region, minor effects are expected from them. Nevertheless, it would be very interesting in the future to study systematically the FF contributions in all the isotopes in this mass region.

Another feature observed in the present calculations is the existence of some odd-even staggering effect in the calculated half-lives, which is not observed experimentally. This effect is particularly evident in Ru and Pd isotopes. There are not many calculations involving simultaneously even-even and odd- isotopes, but some of them exhibit some sort of staggering effect as well (33). The appearance of this effect in the half-lives suggests some deficiency in the model that might be related to the determination of ground-state energies in the odd- systems (65). Unfortunately, there are more sources of uncertainty related to the odd- systems that should be considered as well (66), such as the spin and parity assignments, the blocking procedures or the treatment of the 1qp excitations involving the odd nucleon. This issue will be the subject of a future investigation in this direction.

It is also interesting to look for the simultaneous appearance of structural effects that eventually can appear in different observables. One example can be seen in the evolution of the experimental half-lives with the number of neutrons in the isotopic chains. At some points one observes discontinuities in the general trends of behavior, such as in the mass regions Se, Kr, Sr, and Ru. These experimental findings on the half-lives are correlated with the shape transitions in Figs. 36 predicted in the model. One cannot state firmly that these sharp changes in the behavior of the half-lives are signatures of shape transitions, but certainly this correlation cannot be discarded given that a change of the deformation in the nuclear system involves a structural change to whom the half-lives are also sensitive.

Finally, the impact of deformation on the decay properties can be better appreciated in a systematic comparison of the half-lives calculated with both the spherical approximation and the deformation that corresponds to the minimum of the PEC for each isotope. Then, Fig. 26 shows the ratios of the calculated and experimental half-lives for two sets of data corresponding to a spherical calculation (open dots) and to a deformed calculation (solid dots) at the self-consistent deformation that gives the minimum of the PECs. These ratios are plotted as a function of the experimental half-lives (a) and as a function of the quadrupole deformation at the minimum of the PECs (b). To increase the size of the sample, besides the isotopes considered in this work with measured half-lives, I have also included the set of Zr and Mo neutron-rich isotopes studied in Ref. (38) with measured half-lives. In the upper panel of Fig. 26 (a) one can see how deformation improves the description of the half-lives. Practically all the full black dots are contained within the horizontal lines defining the region of one order of magnitude agreement. On the other hand, the results from the spherical calculation are more spread out with larger discrepancy with experiment. One can also see that the results are better in both spherical and deformed calculations for shorter half-lives, whereas the results for larger half-lives show sizable deviations. The latter correspond to isotopes close to the valley of stability with small -values, where the half-lives are only sensitive to the small portion of the GT strength distribution at low excitation energies below . In the lower panel (b) one can see the results from a different point of view and it can be studied whether deformation improves the results evenly in the whole range of deformations or whether its effect is stronger at large deformations. Three regions of accumulation of results can be distinguished. Two of them correspond to well deformed nuclei located at and . In these regions the deformed calculations clearly improve the results from the spherical ones that show a tendency to underestimate the experiment. The other region corresponds to values, where nuclei are softer or transitional and the deformed formalism should be improved. In this region the results are more scattered than in the well deformed regions, but the deformed calculations show deviations that rarely exceed one order of magnitude, still representing an improvement over the spherical results.

In order to have a quantitative estimation of the quality of the various calculations, following the analysis made in Ref. (29), the logarithms of the ratios of the calculated and experimental half-lives are introduced through the quantities


Then, the average position of the points, , the standard deviation, , and the total error, , are defined as


and their corresponding factors , , and . The analysis of the results shown in Fig. 26 involving nuclei leads to the values , , and in the spherical case and , , and in the deformed one, showing clearly the improvement achieved with the deformed formalism.

Iv Conclusions

A microscopic approach based on a deformed QRPA calculation on top of a self-consistent mean field obtained with the SLy4 Skyrme interaction has been used to study the nuclear structure and the decay properties of even and odd neutron-rich isotopes in the mass region . The nuclear model and interaction have been successfully tested in the past providing good agreement with the available experimental information on bulk properties all along the nuclear chart. Decay properties in different mass regions have been well reproduced as well. The structural isotopic evolution has been studied from their PECs. Depending on the isotopic chain, a large variety of nuclear shapes is found, including spherical shapes, well developed deformed shapes, and transitional soft shapes. Charge radii have been also investigated, showing the connection between a discontinuous behavior in the isotopic trend with a shape transition and comparing the results with the available measurements from laser spectroscopy. Then, Gamow-Teller strength distributions and -decay half-lives have been computed for the equilibrium shapes.

The isotopic evolution of the GT strength distributions exhibits some typical features, such as GT resonances increasing in energy and strength as the number of neutrons increases. Effects of deformation are hard to see on a global scale, but they become apparent in the low excitation energy below energies, a region that determines the half-lives. Half-lives have been calculated using energies calculated with the force SLy4. In general, a reasonable agreement with experiment is obtained, especially in the short-lived nuclei of Ge, Se, Kr, Sr, and Ru isotopes. The results are comparable to other calculations using different approaches for the mean field and/or residual interactions. Special difficulties are found to describe properly the half-lives of the lighter Ge isotopes and the Pd isotopes. These are examples of transitional nuclei where the nuclear structure is more involved and the concept of a well defined shape might not be meaningful.

A systematic comparison of the ratios of the calculated and experimental half-lives has been done using both spherical and deformed calculations, showing that the inclusion of deformation improves significantly the description of the decay properties.

Experimental information on the energy distribution of the GT strength is a valuable piece of knowledge about nuclear structure in this mass region. The study of these distributions is within the current experimental capabilities in the case of the lighter isotopes considered in this work. Here, I have presented theoretical predictions for them based on microscopic calculations. Similarly, measuring the half-lives of the heavier isotopes will be highly beneficial to model the r process and to constrain theoretical nuclear models. This possibility is also open within present capabilities at RIKEN.

This work was supported by Ministerio de Economía y Competitividad (Spain) under Contract No. FIS2011–23565 and the Consolider-Ingenio 2010 Program CPAN CSD2007-00042.


  1. J. L. Wood, K. Heyde, W. Nazarewicz, M. Huyse, and P. Van Duppen, Phys. Rep. 215, 101 (1992).
  2. K. Heyde and J. L. Wood, Rev. Mod. Phys. 83, 1467 (2011).
  3. J. Xiang, Z. P. Li, Z. X. Li, J. M. Yao, and J. Meng, Nucl. Phys. A873, 1 (2012).
  4. H. Mei, J. Xiang, J. M. Yao, Z. P. Li, and J. Meng, Phys. Rev. C 85, 034321 (2012).
  5. P. Bonche, H. Flocard, P.-H. Heenen, S. J. Krieger, and M. S. Weiss, Nucl. Phys. A443, 39 (1985).
  6. M. Bender, G. F. Bertsch, and P.-H. Heenen, Phys. Rev. C 78, 054312 (2008).
  7. M. Bender, K. Bennaceur, T. Duguet, P.-H. Heenen, T. Lesinski, and J. Meyer, Phys. Rev. C 80, 064302 (2009).
  8. R. Rodriguez-Guzman, P. Sarriguren, L. M. Robledo, and S. Perez-Martin, Phys. Lett. B 691, 202 (2010).
  9. R. Rodriguez-Guzman, P. Sarriguren, and L. M. Robledo, Phys. Rev. C 82, 044318 (2010); Phys. Rev. C 82, 061302(R) (2010); Phys. Rev. C 83, 044307 (2011).
  10. T. Sumikama et al., Phys. Rev. Lett. 106, 202501 (2011).
  11. C. Y. Wu et al., Phys. Rev. C 70, 064312 (2004).
  12. W. Urban et al., Eur. Phys. J. A 22, 241 (2004).
  13. H. Mach et al., Phys. Lett. B 230, 21 (1989); Phys. Rev. C 41, 350 (1990).
  14. C. Goodin et al., Nucl. Phys. A787, 231c (2007).
  15. W. Urban et al., Nucl. Phys. A689, 605 (2001).
  16. B. Cheal and K. T. Flanagan, J. Phys. G: Nucl. Part. Phys. 37, 113101 (2010).
  17. E. M. Burbidge, G. M. Burbidge, W. A. Fowler, and F. Hoyle, Rev. Mod. Phys. 29, 547 (1959).
  18. J. J. Cowan, F.-K. Thielemann, and J. W. Truran, Phys. Rep. 208, 267 (1991).
  19. K.-L. Kratz, J.-P. Bitouzet, F.-K. Thielemann, P. Möller, and B. Pfeiffer, Ap. J. 403, 216 (1993).
  20. N. Nishimura, T. Kajino, G. J. Mathews, S. Nishimura, and T. Suzuki, Phys. Rev. C 85, 048801 (2012).
  22. J. Pereira et al., Phys. Rev. C 79, 035806 (2009).
  23. S. Nishimura et al., Phys. Rev. Lett. 106, 052502 (2011).
  24. M. Quinn et al., Phys. Rev. C 85, 035807 (2012).
  25. J. Engel, M. Bender, J. Dobaczewski, W. Nazarewicz, and R. Surman, Phys. Rev. C 60, 014302 (1999).
  26. I. N. Borzov, J. J. Cuenca-García, K. Langanke, G. Martínez-Pinedo, and F. Montes, Nucl. Phys. A814, 159 (2008).
  27. T. Niksic, T. Marketin, D. Vretenar, N. Paar, and P. Ring, Phys. Rev. C 71, 014308 (2005).
  28. J. Krumlinde and P. Möller, Nucl. Phys. A417, 419 (1984).
  29. P. Möller, B. Pfeiffer, and K.-L. Kratz, Phys. Rev. C 67, 055802 (2003).
  30. H. Homma, E. Bender, M. Hirsch, K. Muto, H. V. Klapdor-Kleingrothaus, and T. Oda, Phys. Rev. C 54, 2972 (1996).
  31. K. Muto, E. Bender, T. Oda, and H. V. Klapdor-Kleingrothaus, Z. Phys. A 341, 407 (1992).
  32. F. Frisk, I. Hamamoto, and X.Z. Zhang, Phys. Rev. C 52, 2468 (1995).
  33. D. L. Fang, B. A. Brown, and T. Suzuki, Phys. Rev. C 88, 024314 (2013).
  34. D. Ni and Z. Ren, Phys. Rev. C 89, 064320 (2014); J. Phys. G: Nucl. Part. Phys. 39, 125105 (2012).
  35. K. Yoshida, Prog. Theor. Exp. Phys. 113D02 (2013).
  36. M. Martini, S. Péru, and S. Goriely, Phys. Rev. C 89, 044306 (2014).
  37. P. Sarriguren and J. Pereira, Phys. Rev. C 81, 064314 (2010).
  38. P. Sarriguren, A. Algora, and J. Pereira, Phys. Rev. C 89, 034311 (2014).
  39. P. Sarriguren, E. Moya de Guerra, A. Escuderos, and A. C. Carrizo, Nucl. Phys. A635, 55 (1998).
  40. P. Sarriguren, E. Moya de Guerra, and A. Escuderos, Nucl. Phys. A691, 631 (2001).
  41. E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A635, 231 (1998).
  42. M. V. Stoitsov, J. Dobaczewski, W. Nazarewicz, S. Pittel, and D. J. Dean, Phys. Rev. C 68, 054312 (2003).
  43. D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626 (1972); D. Vautherin, Phys. Rev. C 7, 296 (1973).
  44. G. Audi et al., Chinese Physics C 36, 1157 (2012); 1603 (2012).
  45. J.-U. Nabi and H. V. Klapdor-Kleingrothaus, At. Data Nucl. Data Tables 88, 237 (2004).
  46. P. Sarriguren, E. Moya de Guerra, and A. Escuderos, Nucl. Phys. A658, 13 (1999).
  47. P. Sarriguren, E. Moya de Guerra, and A. Escuderos, Phys. Rev. C 64, 064306 (2001).
  48. P. Sarriguren, R. Alvarez-Rodríguez, and E. Moya de Guerra, Eur. Phys. J. A 24, 193 (2005).
  49. E. Poirier et al., Phys. Rev. C 69, 034307 (2004).
  50. E. Nácher et al., Phys. Rev. Lett. 92, 232501 (2004).
  51. A.B. Pérez-Cerdán et al., Phys. Rev. C 88, 014324 (2013).
  52. A. Bohr and B. Mottelson, Nuclear Structure, Vols. I and II, (Benjamin, New York 1975).
  53. N. B. Gove and M. J. Martin, Nucl. Data Tables 10, 205 (1971).
  54. J. Skalski, S. Mizutori, and W. Nazarewicz, Nucl. Phys. A617, 281 (1997).
  55. P. Möller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, At. Data Nucl. Data Tables 59, 185 (1995).
  56. G. A. Lalazissis, S. Raman, and P. Ring, At. Data Nucl. Data Tables 71, 1 (1999).
  57. S. Hilaire and M. Girod, Eur. Phys. J. A 33, 237 (2007);
  58. M. Albers et al., Phys. Rev. Lett. 108, 062701 (2012).
  59. I. Angeli, At. Data Nucl. Data Tables 87, 185 (2004).
  60. D. Fang, A. Faessler, V. Rodin, M. S. Yousef, and F. Simkovic, Phys. Rev. C 81, 037303 (2010).
  61. K. Langanke and G. Martinez-Pinedo, Nucl. Phys. A673, 481 (2000).
  62. P. Sarriguren, Phys. Rev. C 79, 044315 (2009).
  63. P. Sarriguren, Phys. Lett. B 680, 438 (2009).
  64. M. V. Stoitsov, J. Dobaczewski, W. Nazarewicz, and P. Ring, Comp. Phys. Comm. 167, 43 (2005).
  65. T. Duguet, P. Bonche, P.-H. Heenen, and J. Meyer, Phys. Rev. C 65, 014310 (2001); 014311 (2001).
  66. N. Schunck, J. Dobaczewski, J. McDonnell, J. Moré, W. Nazarewicz, J. Sarich, and M. V. Stoitsov, Phys. Rev. C 81, 024316 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description