Loading [MathJax]/jax/output/HTML-CSS/jax.js

Determination of isospin asymmetry effects on α-decay

  • In this study, we compared the effect of the isospin asymmetry of proton and neutron density distributions in the neutron skin-type (NST) case and in the Hartree-Fock formalism (HF) on the half-life of alpha emitters with the atomic number in the range of 82Z92. The NST case and HF formalism based on the Skyrme-SLy4 effective interaction reveal different isospin asymmetries for selected alpha emitters. Furthermore, the obtained results reveal an increase in the α-decay widths of about 30% for the NST case in comparison with the equivalent values obtained by HF formalism. The standard deviations for calculated half-lives within the NST case and HF formalism are about 0.438 and 0.391, respectively.
  • The α-decay process is essentially a quantum tunneling effect, which was first described by Gamow and by Condon and Gurney [1, 2]. Several theoretical calculations have been performed to understand the microscopic nature of this phenomenon, obtain the α-decay width and gain information about the structure of the nuclei. These studies are based on various theoretical models such as the shell model, fission-like model, and cluster model [310]. In general, the α-decay constant λ is the product of the preformation factor Pα, alpha particle penetration through the potential barrier P0, and the assault frequency ν0. The nuclear values of these parameters are generally dependent on nuclear structure and internuclear potential.

    Because the formation of the alpha particle, including valence nucleons, is closely related to the nuclear structure and alpha preformation factor containing most of the nuclear structure information, the preformation factor is an important component of the α-decay study. Consequently, the parent nucleus can be approximated as a simple two-particle system before emission, such that the alpha emitter is assumed as a system consisting of a core nucleus and a single alpha cluster [11, 12]. The surface nucleons and isospin asymmetry of these is expected to play a more important role in the preformation of the alpha nucleus [1116].

    The neutron skin thickness is defined as the difference between neutron and proton root-mean-squares; ΔrnpRnRp, which can be related to the difference in nuclear asymmetric equation states inside and on the surface of the nucleus. Various studies have been conducted on the determination of proton and neutron density distributions, the experimental root-mean-square (rms) measurements of nuclei and the effect of neutron skin thickness on α-decay [17-19]. In some studies, to determine the Δrnp, it is assumed that the proton diffuseness ap and neutron diffuseness an are the same, while in fact, at the nucleus surface this is not exactly the case. Although this assumption is true for light nuclei with N=Z, it is not validated for the nuclei with N>Z, as the proton and neutron density distributions of heavy nuclei are not exactly the same on the nuclear surface. Moreover, the shape of the neutron density distribution is more diffuse than the charge density distribution [18, 20, 21].

    The HF formalism and Skyrme forces parameterized to describe various ground state properties of nuclei are successful theories for the determination of more realistic nuclei densities [16, 17, 2224]. The effects of the isospin asymmetry are embedded so well in this formalism that it can be applied to surface studies of the nuclei. The specifics of the effects of these different isospin asymmetries on the half-lives of the α-decay can be expected within the different Δrnp originated by HF formalism and the NST case.

    For this purpose, using isospin asymmetry data obtained from the HF formalism and the NST case, we aim to determine their effects on the α-decay width, half-life, and on the preformation factor of α-decay for even-even isotope groups with 82Z92.

    This paper is organized as follows: In Sec. 2, the theoretical framework of the double folding (DF) model is presented in the calculation of nuclear potential, half-life, and the preformation factor. Sec. 3 presents the results and discussion. A summary is provided in Sec. 4.

    The total alpha core potential is constructed by nuclear and Coulomb terms, and it is given by

    VT(r)=VC(r)+VN(r).

    (1)

    Nuclear and Coulomb potentials are obtained from the DF model by folding the densities of the alpha particle and the daughter nuclei with effective M3Y interactions [2529],

    VN(r)=λVF(r)=λρ1(r1)ν(s=|r+r2r1|)×ρ2(r2)d3r1d3r2,

    (2)

    where s is the distance between the two interacting nucleons, and ρ1(r1) and ρ2(r2) are the density distributions of the alpha particle and daughter nuclei, respectively. The ν(s=|r+r2r1|) is the effective nucleon-nucleon interaction, where we employ M3Y Paris-NN forces with zero range approximation for the nucleon-nucleon interaction [28].

    The parameter λ in Eq. (2) is a re-normalization factor. To reproduce the experimental α-decay energies, the λ-value is determined under the Bohr-Sommerfeld quantization condition [30],

    r2r12μ2|V(r)Q|dr=(2n+1)π2=(GL+1)π2,

    (3)

    where r1 and r2 (and later r3) are the classical return points, determined from the relation V(R)=Q (decay energy of alpha). The global quantum number G is estimated by the Wildermuth rule [31],

    G=2N+L=4i=1(2ni+li)=4i=1gi,

    (4)

    where N is the number of nodes in the alpha core wave function, L is the angular momentum of the cluster motion, and gi is the oscillation quantum number of the cluster nucleon. For N and Z, the parent nucleus is for (Z,N)82, gi=4, for 82<(Z,N)126, gi=5, and for (Z,N)>126, gi=6, corresponding to the order on the oscillator layers of 4ω, 5ω, and 6ω, respectively. In the quasiclassical approximation, the α-decay width Γα is obtained by [30, 32]

    Γα=PαF24μexp(2r3r2k(r)dr),

    (5)

    where Pα is the alpha preformation factor, which can be extracted from the ratio of the experimental decay width to the calculated one [27].

    The normalization factor is,

    Fr2r1dr2k(r)=1,

    (6)

    where k(r)=2μ2|V(r)Q| is the wave number. Then, the α-decay half-life is related to the width by T1/2=ln2/Γα.

    In Eq. (2), the mass density distribution of the spherical alpha particle ρ1(r1) has a Gaussian form given by

    ρ1(r1)=0.4229exp(0.7024r21),

    (7)

    and the mass density distribution of the daughter nucleus is written as

    ρτ2(r2)=ρτ01+exp[(r2Rτ1/2)/aτ],

    (8)

    where τ=n or p and the value of ρτ0 is fixed by integrating the density distribution equivalent to the mass number of the daughter nucleus. The half-density radius Rτ1/2 is related to the mass number of the nucleus. The neutron and proton rms radii of the daughter nuclei are

    R(S)τ≡<r2>1/2=[ρτ2(r)r4drρτ2(r)r2dr]12.

    (9)

    Considering the importance of deformation of the nuclei on the calculation of this parameter and α-decay study [6, 26, 33], we have employed the axially symmetric deformed distribution of the daughter nucleus in the DF model calculations using the following profile [33],

    ρτ2(r2,θ)=ρτ01+exp{[r2Rτ1/2(θ)]/aτ},

    (10)

    where the half-density radius Rτ1/2(θ) is parameterized as

    Rτ1/2(θ)=rτA1/3[1+β2Y20(θ)+β4Y40(θ)],

    (11)

    where β2 and β4 are the quadrupole and hexadecapole deformation parameters of daughter nuclei, respectively, with their values taken from Ref. [34]. The rms neutron and proton radii using the deformed nuclear distributions are calculated by

    R(D)τ[ρτ2(r,θ)r4drdΩρτ2(r,θ)r2drdΩ]12.

    (12)

    The multipole expansion of the deformed distribution ρτ2(r2) is used and then the orientation-dependent DF potential can be evaluated as the sum of different multipole components [18, 33]. The effect of nuclear deformations on the decay width could also be evaluated by integrating the partial width along the direction,

    Γ=π/20Γ(β)sinβdβ,

    (13)

    where β is the angle between the symmetry axis of the daughter nucleus and the direction of the α-particle emission.

    We studied two cases. In the first case, the radius and diffuseness parameters were obtained by HF calculations based on the Skyrme-SLy4 effective interaction that successfully explains the ground state properties of nuclei [3537]. The other case is the NST with the assumption rn>rp and an=ap=0.54 that the radius parameters rn and rp have determined in terms of the neutron and protons rms radii. The proton root-mean-square radius Rp is extracted by empirical data, and the neutron rms radius Rn is determined by the relation Rn=Rp+Δrnp, in which the Δrnp is obtained by

    Δrnp=(0.04±0.03)+(1.01±0.15)δ,

    (14)

    which presents the linear relationship between the isospin asymmetry δ=(NZ)/(N+Z) and Δrnp [38].

    For the nuclei with N>Z, the rms of the neutron density distribution is larger than the rms of the proton density distribution, hence the neutron density distribution can be different from the proton density distribution. Accordingly, to calculate the nucleus density distribution, it is not generally correct to assume that ρn=NZρp when N>Z. Instead, it is better to separately employ the neutron-to-proton densities.

    The ratio of the neutron-to-proton density versus distance from the nucleus center was plotted for some neutron-rich nuclei 210Pb, 218Po, 214Rn, 226Ra, 230Th, and 236U in Fig. 1. The ratio of densities is equal to N/Z in the range R=04 fm, followed by a sharp rise. This rising trend is more dramatic in HF compared to NST. To determine the effect of the neutron skin thickness on the α-decay, the neutron skin thickness was calculated by the NST and HF formalism versus the isospin asymmetry parameter δ compared in Fig. 2. This figure illustrates the variation between the calculated values of Δrnp within these two different approaches.

    Figure 1

    Figure 1.  (color online) Ratio of neutron-to-proton density versus distance from nucleus center.

    Figure 2

    Figure 2.  (color online) Neutron skin thickness Δrnp versus asymmetry parameter δ for isotope groups of elements Pb,Po,Rn,Ra,Th and U for two cases NST and HF.

    Because of the important role of the neutron and proton density distributions in calculating the nuclear potential, typically, the total potential for the 216Po, whose deformation parameters are zero, in the NST and HF cases were plotted in Fig. 3(a). In the range of R=9.912.5(fm), the potential value in the HF case is greater by 0.15% to 2.5% in comparison with the NST case. Because of the effect of the neutron skin thickness in NST and HF cases on the total potential, it is expected that the potential variations will affect the penetration probabilities and consequently affect the α-decay width. Therefore, the α-decay width was investigated in NST and HF cases for isotope groups of elements Pb,Po,Rn,Ra,Th, and U. Typically, the treatment of the calculated α-decay width for the Po isotope group in terms of the neutron number was presented in Fig. 3(b). In this figure, the α-decay widths of the NST case are greater by about 30% than their equivalent values of HF case. Therefore, this can affect the α-decay half-life and the extracted preformation factor. Within the calculated potentials, α-decay half-lives were calculated for the even-even isotope groups of elements Pb,Po,Rn,Ra,Th, and U by employing relations 6 to 8. In the NST case, the diffuseness parameters were fixed, and the radius parameters rn and rp were determined in terms of the neutron and proton rms radii. The rms proton radius Rp was taken from the experimental rms charge radii, and the rms neutron radius Rn was obtained from the relationship Rn=Rp+Δrnp. Since the experimental rms charge radii for some isotopes was not available, we have assumed the values of rms charge radii for these isotopes from a simple relationship, given by Rp=(r0+r1A2/3+r2A4/3)A1/3 with r0=0.891(2) fm, r1=1.52(3) fm, and r2=2.8(1) fm. This relationship was adapted from a systematic study on available data of the experimental rms charge radii [39].

    Figure 3

    Figure 3.  (color online) (a) Total potential in terms of distance from center of nucleus for 216Po, assuming that HF and NST and (b) is the α-decay width versus neutron number for even-even Po isotopes.

    In Fig. 4, the logarithm of the preformation factor, calculated from the relation Pα=Γexpα/Γcalα [27], is plotted according to the parent neutron number for even-even isotopes with the atomic number in the range 82Z92. In general, the variation treatments in NST and HF cases are the same. Results indicate that the curve of HF has a value closer to the experimental data. The standard deviation from the experimental data obtained by σ={Ni=1[log10(Tical/Tiexp)]2/N}1/2 is equal to 0.438 and 0.391 for NST and HF cases, respectively. This confirms that the calculated half-lives in the HF case are closer to experimental data.

    Figure 4

    Figure 4.  (color online) Logarithms of preformation factor of α-decay half-life versus parent nucleus neutron number, for even-even isotope groups of Pb,Po,Rn,Ra,Th, and U in HF and NST cases.

    In this study, we investigated the effects of isospin asymmetries resulting from different proton and neutron density distributions on the total potentials, α-decay widths, and half-lives of even-even alpha emission nuclei with the atomic number in the range 82Z92. The neutron and proton density distributions of the daughter nuclei were calculated in the two-parameter Fermi form of nucleon density distributions. The radius and diffuseness parameter values of the nuclei were obtained from the NST and HF formalism, leading to different isospin asymmetries. The results indicated that the values of Δrnp for the two NST and HF cases affect the α-decay half-life and the preformation factor.

    The standard deviation of the calculated half-lives from experimental data reached about 0.438 for the NST case. However, the results illustrated that the isospin asymmetry within HF formalism lowered the standard deviation in comparison with its equivalent value obtained from NST. The standard deviation for the HF case is about 0.391, which indicates that the calculated half-life, in this case, is closer to the value obtained by experimental data.

    [1] G. Z. Gamow, Phys., 51: 204 (1928) doi: 10.1007/BF01343196
    [2] R. W. Gurney and E. U. Condon, Nature, 122: 439 (1928) doi: 10.1038/122439a0
    [3] K. Varga, R. Lovas, and R. Liotta, Phys. Rev. Lett., 69: 37 (1992) doi: 10.1103/PhysRevLett.69.37
    [4] B. A. Brown, Phys. Rev. C, 46: 811 (1992) doi: 10.1103/PhysRevC.46.811
    [5] B. Buck et al, At. Data Nucl. Data Tables, 54: 53 (1993) doi: 10.1006/adnd.1993.1009
    [6] C. Xu and Z. Ren, Phys. Rev. C, 74: 014304 (2006)
    [7] G. Royer and H. F. Zhang, Phys. Rev. C, 77: 037602 (2008) doi: 10.1103/PhysRevC.77.037602
    [8] H. F. Zhang, Y. Gao, N. Wang et al, Phys. Rev. C, 85: 14325 (2012) doi: 10.1103/PhysRevC.85.014325
    [9] Y. Ren and Z.Ren, Phys. Rev. C, 85: 044608 (2012) doi: 10.1103/PhysRevC.85.044608
    [10] Y. Qian and Z. Ren, Phys. lett. B, 738: 87 (2014) doi: 10.1016/j.physletb.2014.09.024
    [11] M. Ismail and A. Adel, J. Phys. G: Nucl. Phys., 44: 125106 (2017) doi: 10.1088/1361-6471/aa957f
    [12] D. Ni and Z. Ren, Phys. Rev. C, 80: 014314 (2009) doi: 10.1103/PhysRevC.80.014314
    [13] D. Ni and Z. Ren, Nucl. Phys. A, 825: 145 (2009) doi: 10.1016/j.nuclphysa.2009.04.010
    [14] E. Shin, Y. Lim, C. H. Hyun et al, Phys. Rev. C, 94: 024320 (2016) doi: 10.1103/PhysRevC.94.024320
    [15] Y. Qian and Z. Ren, J. Phys. G: Nucl. Phys., 43: 065102 (2016) doi: 10.1088/0954-3899/43/6/065102
    [16] D. T. Akrawy, H. Hassanabadi, S. S. Hosseini et al, Nucl. Phys. A, 975: 19 (2018) doi: 10.1016/j.nuclphysa.2018.04.001
    [17] W. M. Seif and A. Abdurrahman, Chinese Phys. C, 42: 014106 (2018) doi: 10.1088/1674-1137/42/1/014106
    [18] D. Ni and Z. Ren, Phys. Rev. C, 92: 054322 (2015) doi: 10.1103/PhysRevC.92.054322
    [19] S. Typel, Phys. Rev. C, 89: 064321 (2014) doi: 10.1103/PhysRevC.89.064321
    [20] G. R. Satchler and W. G. Love, Phys. Rep., 55: 183 (1979) doi: 10.1016/0370-1573(79)90081-4
    [21] D. Ni and Z. Ren, Phys. Rev. C, 93: 054318 (2016) doi: 10.1103/PhysRevC.93.054318
    [22] W. M. Seif, M. Shalaby, and M. F. Alrakshy, Phys. Rev. C, 84: 064608 (2011) doi: 10.1103/PhysRevC.84.064608
    [23] M. Ismail and A. Adel, Phys. Rev. C, 86: 014616 (2012) doi: 10.1103/PhysRevC.86.014616
    [24] J. C. Pei, F. R. Xu, Z. J. Lin et al, Phys. Rev. C, 76: 044326 (2007) doi: 10.1103/PhysRevC.76.044326
    [25] C. Xu and Z. Z. Ren, Nucl. Phys. A, 753: 174 (2005) doi: 10.1016/j.nuclphysa.2005.02.125
    [26] M. Ismail, A. Y. Ellithi, M. M. Botros et al, Phys. Rev. C, 81: 024602 (2010) doi: 10.1103/PhysRevC.81.024602
    [27] Y. Qian and Z. Ren, Phys. lett. B, 777: 298 (2018) doi: 10.1016/j.physletb.2017.12.046
    [28] D. T. Khoa, G. R. Satchler, and W. von Oertzen, Phys. Rev. C, 56: 954 (1997) doi: 10.1103/PhysRevC.56.954
    [29] D. N. Basu et al, J. Phys. G: Nucl. Part. Phys., 29: 2079-2085 (2003)
    [30] B. Buck, A. C. Merchant, and S. M. Perez, Phys. Rev. C, 45: 2247 (1992) doi: 10.1103/PhysRevC.45.2247
    [31] K. Wildermuth and Y. C. Tang, A Unified Theory of the Nucleus(Academic Press, New York, 1977)
    [32] S. A. Gurvitz and G.\begin{document}${\small{\rm{\ddot K}}}$\end{document}albermann, Phys. Rev. Lett., 59: 262 (1987)
    [33] D. Ni, Z. Ren, T. Dong, and Y. Qian, Phys. Rev. C, 87: 024310 (2013) doi: 10.1103/PhysRevC.87.024310
    [34] P. Möller, A. J. Sierk, T. Ichikawa et al, At. Data Nucl. Data Tables, 109-110: 1 (2016) doi: 10.1016/j.adt.2015.10.002
    [35] D. E. Ward, B. G. Carlsson, and S.Åberg, Phys. Rev. C, 88: 064316 (2013) doi: 10.1103/PhysRevC.88.064316
    [36] W. M. Seif, Phys. Rev. C, 91: 014322 (2015) doi: 10.1103/PhysRevC.91.014322
    [37] S. Goriely, M. Samyn, M. Bender et al, Phys. Rev. C, 68: 054325 (2003) doi: 10.1103/PhysRevC.68.054325
    [38] A. Trzcińska, J. Jastrzȩbski, P. Lubiński et al, Phys. Rev. Lett., 87: 082501 (2001) doi: 10.1103/PhysRevLett.87.082501
    [39] I. Angeli, At. Data Nucl. Data Tables, 87: 185 (2004) doi: 10.1016/j.adt.2004.04.002
  • [1] G. Z. Gamow, Phys., 51: 204 (1928) doi: 10.1007/BF01343196
    [2] R. W. Gurney and E. U. Condon, Nature, 122: 439 (1928) doi: 10.1038/122439a0
    [3] K. Varga, R. Lovas, and R. Liotta, Phys. Rev. Lett., 69: 37 (1992) doi: 10.1103/PhysRevLett.69.37
    [4] B. A. Brown, Phys. Rev. C, 46: 811 (1992) doi: 10.1103/PhysRevC.46.811
    [5] B. Buck et al, At. Data Nucl. Data Tables, 54: 53 (1993) doi: 10.1006/adnd.1993.1009
    [6] C. Xu and Z. Ren, Phys. Rev. C, 74: 014304 (2006)
    [7] G. Royer and H. F. Zhang, Phys. Rev. C, 77: 037602 (2008) doi: 10.1103/PhysRevC.77.037602
    [8] H. F. Zhang, Y. Gao, N. Wang et al, Phys. Rev. C, 85: 14325 (2012) doi: 10.1103/PhysRevC.85.014325
    [9] Y. Ren and Z.Ren, Phys. Rev. C, 85: 044608 (2012) doi: 10.1103/PhysRevC.85.044608
    [10] Y. Qian and Z. Ren, Phys. lett. B, 738: 87 (2014) doi: 10.1016/j.physletb.2014.09.024
    [11] M. Ismail and A. Adel, J. Phys. G: Nucl. Phys., 44: 125106 (2017) doi: 10.1088/1361-6471/aa957f
    [12] D. Ni and Z. Ren, Phys. Rev. C, 80: 014314 (2009) doi: 10.1103/PhysRevC.80.014314
    [13] D. Ni and Z. Ren, Nucl. Phys. A, 825: 145 (2009) doi: 10.1016/j.nuclphysa.2009.04.010
    [14] E. Shin, Y. Lim, C. H. Hyun et al, Phys. Rev. C, 94: 024320 (2016) doi: 10.1103/PhysRevC.94.024320
    [15] Y. Qian and Z. Ren, J. Phys. G: Nucl. Phys., 43: 065102 (2016) doi: 10.1088/0954-3899/43/6/065102
    [16] D. T. Akrawy, H. Hassanabadi, S. S. Hosseini et al, Nucl. Phys. A, 975: 19 (2018) doi: 10.1016/j.nuclphysa.2018.04.001
    [17] W. M. Seif and A. Abdurrahman, Chinese Phys. C, 42: 014106 (2018) doi: 10.1088/1674-1137/42/1/014106
    [18] D. Ni and Z. Ren, Phys. Rev. C, 92: 054322 (2015) doi: 10.1103/PhysRevC.92.054322
    [19] S. Typel, Phys. Rev. C, 89: 064321 (2014) doi: 10.1103/PhysRevC.89.064321
    [20] G. R. Satchler and W. G. Love, Phys. Rep., 55: 183 (1979) doi: 10.1016/0370-1573(79)90081-4
    [21] D. Ni and Z. Ren, Phys. Rev. C, 93: 054318 (2016) doi: 10.1103/PhysRevC.93.054318
    [22] W. M. Seif, M. Shalaby, and M. F. Alrakshy, Phys. Rev. C, 84: 064608 (2011) doi: 10.1103/PhysRevC.84.064608
    [23] M. Ismail and A. Adel, Phys. Rev. C, 86: 014616 (2012) doi: 10.1103/PhysRevC.86.014616
    [24] J. C. Pei, F. R. Xu, Z. J. Lin et al, Phys. Rev. C, 76: 044326 (2007) doi: 10.1103/PhysRevC.76.044326
    [25] C. Xu and Z. Z. Ren, Nucl. Phys. A, 753: 174 (2005) doi: 10.1016/j.nuclphysa.2005.02.125
    [26] M. Ismail, A. Y. Ellithi, M. M. Botros et al, Phys. Rev. C, 81: 024602 (2010) doi: 10.1103/PhysRevC.81.024602
    [27] Y. Qian and Z. Ren, Phys. lett. B, 777: 298 (2018) doi: 10.1016/j.physletb.2017.12.046
    [28] D. T. Khoa, G. R. Satchler, and W. von Oertzen, Phys. Rev. C, 56: 954 (1997) doi: 10.1103/PhysRevC.56.954
    [29] D. N. Basu et al, J. Phys. G: Nucl. Part. Phys., 29: 2079-2085 (2003)
    [30] B. Buck, A. C. Merchant, and S. M. Perez, Phys. Rev. C, 45: 2247 (1992) doi: 10.1103/PhysRevC.45.2247
    [31] K. Wildermuth and Y. C. Tang, A Unified Theory of the Nucleus(Academic Press, New York, 1977)
    [32] S. A. Gurvitz and G.\begin{document}${\small{\rm{\ddot K}}}$\end{document}albermann, Phys. Rev. Lett., 59: 262 (1987)
    [33] D. Ni, Z. Ren, T. Dong, and Y. Qian, Phys. Rev. C, 87: 024310 (2013) doi: 10.1103/PhysRevC.87.024310
    [34] P. Möller, A. J. Sierk, T. Ichikawa et al, At. Data Nucl. Data Tables, 109-110: 1 (2016) doi: 10.1016/j.adt.2015.10.002
    [35] D. E. Ward, B. G. Carlsson, and S.Åberg, Phys. Rev. C, 88: 064316 (2013) doi: 10.1103/PhysRevC.88.064316
    [36] W. M. Seif, Phys. Rev. C, 91: 014322 (2015) doi: 10.1103/PhysRevC.91.014322
    [37] S. Goriely, M. Samyn, M. Bender et al, Phys. Rev. C, 68: 054325 (2003) doi: 10.1103/PhysRevC.68.054325
    [38] A. Trzcińska, J. Jastrzȩbski, P. Lubiński et al, Phys. Rev. Lett., 87: 082501 (2001) doi: 10.1103/PhysRevLett.87.082501
    [39] I. Angeli, At. Data Nucl. Data Tables, 87: 185 (2004) doi: 10.1016/j.adt.2004.04.002
  • 加载中

Cited by

1. Prathapan, K., Deneshan, P., Rajan, M.K.P. et al. Nearly Symmetric, Asymmetric, Very Asymmetric, and Super Asymmetric Fission Half-Lives of Heavy and Superheavy Elements Using a Modified Geiger-Nuttall Law[J]. Brazilian Journal of Physics, 2024, 54(2): 39. doi: 10.1007/s13538-023-01412-8
2. Cheng, S., Wu, W., Cao, L. et al. Isospin effects on α decay and cluster radioactivity[J]. European Physical Journal A, 2022, 58(9): 168. doi: 10.1140/epja/s10050-022-00825-4
3. Dehghani, V., Alavi, S.A., Razavi, R. et al. Cluster decay half-lives using asymmetry dependent densities[J]. Chinese Physics C, 2022, 46(4): 044104. doi: 10.1088/1674-1137/ac4035
4. Ghodsi, O.N., Amiri, M.M. Exploring decay properties in the superheavy region through the double-folding formalism and Skyrme interactions[J]. Physical Review C, 2021, 104(4): 044618. doi: 10.1103/PhysRevC.104.044618
5. Gheshlagh, B.A.. Studying the α -Decay Half-Lives of Nuclei in the Range 68 ≤ Z ≤ 118 by Double-Folding and Cluster Formation Formalisms[J]. Physics of Atomic Nuclei, 2021, 84(5): 675-682. doi: 10.1134/S1063778821050070
6. Seif, W.M., Adamian, G.G., Antonenko, N.V. et al. Correlations of -decay properties and isospin-asymmetry[J]. Physical Review C, 2021, 104(1): 014317. doi: 10.1103/PhysRevC.104.014317
7. He, Y., Yu, X., Zhang, H.-F. Improved empirical formula for α particle preformation factor[J]. Chinese Physics C, 2021, 45(1): 014110. doi: 10.1088/1674-1137/abc684
8. Nejati, S., Ghodsi, O.N. Study of the dependence of alpha decay half-life on the surface symmetry energy[J]. International Journal of Modern Physics E, 2020, 29(9): 2050070. doi: 10.1142/S0218301320500706

Figures(4)

Get Citation
O. N. Ghodsi and B. A. Gheshlagh. Isospin asymmetry effect determinations on α-decay[J]. Chinese Physics C. doi: 10.1088/1674-1137/43/12/124105
O. N. Ghodsi and B. A. Gheshlagh. Isospin asymmetry effect determinations on α-decay[J]. Chinese Physics C.  doi: 10.1088/1674-1137/43/12/124105 shu
Milestone
Received: 2019-05-06
Revised: 2019-09-16
Article Metric

Article Views(3030)
PDF Downloads(34)
Cited by(8)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Determination of isospin asymmetry effects on α-decay

    Corresponding author: O. N. Ghodsi, o.nghodsi@umz.ac.ir
  • Department of Physics, Faculty of Science, University of Mazandaran, P.O. Box 47415-416, Babolsar, Iran

Abstract: In this study, we compared the effect of the isospin asymmetry of proton and neutron density distributions in the neutron skin-type (NST) case and in the Hartree-Fock formalism (HF) on the half-life of alpha emitters with the atomic number in the range of 82Z92. The NST case and HF formalism based on the Skyrme-SLy4 effective interaction reveal different isospin asymmetries for selected alpha emitters. Furthermore, the obtained results reveal an increase in the α-decay widths of about 30% for the NST case in comparison with the equivalent values obtained by HF formalism. The standard deviations for calculated half-lives within the NST case and HF formalism are about 0.438 and 0.391, respectively.

    HTML

    1.   Introduction
    • The α-decay process is essentially a quantum tunneling effect, which was first described by Gamow and by Condon and Gurney [1, 2]. Several theoretical calculations have been performed to understand the microscopic nature of this phenomenon, obtain the α-decay width and gain information about the structure of the nuclei. These studies are based on various theoretical models such as the shell model, fission-like model, and cluster model [310]. In general, the α-decay constant λ is the product of the preformation factor Pα, alpha particle penetration through the potential barrier P0, and the assault frequency ν0. The nuclear values of these parameters are generally dependent on nuclear structure and internuclear potential.

      Because the formation of the alpha particle, including valence nucleons, is closely related to the nuclear structure and alpha preformation factor containing most of the nuclear structure information, the preformation factor is an important component of the α-decay study. Consequently, the parent nucleus can be approximated as a simple two-particle system before emission, such that the alpha emitter is assumed as a system consisting of a core nucleus and a single alpha cluster [11, 12]. The surface nucleons and isospin asymmetry of these is expected to play a more important role in the preformation of the alpha nucleus [1116].

      The neutron skin thickness is defined as the difference between neutron and proton root-mean-squares; ΔrnpRnRp, which can be related to the difference in nuclear asymmetric equation states inside and on the surface of the nucleus. Various studies have been conducted on the determination of proton and neutron density distributions, the experimental root-mean-square (rms) measurements of nuclei and the effect of neutron skin thickness on α-decay [17-19]. In some studies, to determine the Δrnp, it is assumed that the proton diffuseness ap and neutron diffuseness an are the same, while in fact, at the nucleus surface this is not exactly the case. Although this assumption is true for light nuclei with N=Z, it is not validated for the nuclei with N>Z, as the proton and neutron density distributions of heavy nuclei are not exactly the same on the nuclear surface. Moreover, the shape of the neutron density distribution is more diffuse than the charge density distribution [18, 20, 21].

      The HF formalism and Skyrme forces parameterized to describe various ground state properties of nuclei are successful theories for the determination of more realistic nuclei densities [16, 17, 2224]. The effects of the isospin asymmetry are embedded so well in this formalism that it can be applied to surface studies of the nuclei. The specifics of the effects of these different isospin asymmetries on the half-lives of the α-decay can be expected within the different Δrnp originated by HF formalism and the NST case.

      For this purpose, using isospin asymmetry data obtained from the HF formalism and the NST case, we aim to determine their effects on the α-decay width, half-life, and on the preformation factor of α-decay for even-even isotope groups with 82Z92.

      This paper is organized as follows: In Sec. 2, the theoretical framework of the double folding (DF) model is presented in the calculation of nuclear potential, half-life, and the preformation factor. Sec. 3 presents the results and discussion. A summary is provided in Sec. 4.

    2.   Theoretical formalism
    • The total alpha core potential is constructed by nuclear and Coulomb terms, and it is given by

      VT(r)=VC(r)+VN(r).

      (1)

      Nuclear and Coulomb potentials are obtained from the DF model by folding the densities of the alpha particle and the daughter nuclei with effective M3Y interactions [2529],

      VN(r)=λVF(r)=λρ1(r1)ν(s=|r+r2r1|)×ρ2(r2)d3r1d3r2,

      (2)

      where s is the distance between the two interacting nucleons, and ρ1(r1) and ρ2(r2) are the density distributions of the alpha particle and daughter nuclei, respectively. The ν(s=|r+r2r1|) is the effective nucleon-nucleon interaction, where we employ M3Y Paris-NN forces with zero range approximation for the nucleon-nucleon interaction [28].

      The parameter λ in Eq. (2) is a re-normalization factor. To reproduce the experimental α-decay energies, the λ-value is determined under the Bohr-Sommerfeld quantization condition [30],

      r2r12μ2|V(r)Q|dr=(2n+1)π2=(GL+1)π2,

      (3)

      where r1 and r2 (and later r3) are the classical return points, determined from the relation V(R)=Q (decay energy of alpha). The global quantum number G is estimated by the Wildermuth rule [31],

      G=2N+L=4i=1(2ni+li)=4i=1gi,

      (4)

      where N is the number of nodes in the alpha core wave function, L is the angular momentum of the cluster motion, and gi is the oscillation quantum number of the cluster nucleon. For N and Z, the parent nucleus is for (Z,N)82, gi=4, for 82<(Z,N)126, gi=5, and for (Z,N)>126, gi=6, corresponding to the order on the oscillator layers of 4ω, 5ω, and 6ω, respectively. In the quasiclassical approximation, the α-decay width Γα is obtained by [30, 32]

      Γα=PαF24μexp(2r3r2k(r)dr),

      (5)

      where Pα is the alpha preformation factor, which can be extracted from the ratio of the experimental decay width to the calculated one [27].

      The normalization factor is,

      Fr2r1dr2k(r)=1,

      (6)

      where k(r)=2μ2|V(r)Q| is the wave number. Then, the α-decay half-life is related to the width by T1/2=ln2/Γα.

      In Eq. (2), the mass density distribution of the spherical alpha particle ρ1(r1) has a Gaussian form given by

      ρ1(r1)=0.4229exp(0.7024r21),

      (7)

      and the mass density distribution of the daughter nucleus is written as

      ρτ2(r2)=ρτ01+exp[(r2Rτ1/2)/aτ],

      (8)

      where τ=n or p and the value of ρτ0 is fixed by integrating the density distribution equivalent to the mass number of the daughter nucleus. The half-density radius Rτ1/2 is related to the mass number of the nucleus. The neutron and proton rms radii of the daughter nuclei are

      R(S)τ≡<r2>1/2=[ρτ2(r)r4drρτ2(r)r2dr]12.

      (9)

      Considering the importance of deformation of the nuclei on the calculation of this parameter and α-decay study [6, 26, 33], we have employed the axially symmetric deformed distribution of the daughter nucleus in the DF model calculations using the following profile [33],

      ρτ2(r2,θ)=ρτ01+exp{[r2Rτ1/2(θ)]/aτ},

      (10)

      where the half-density radius Rτ1/2(θ) is parameterized as

      Rτ1/2(θ)=rτA1/3[1+β2Y20(θ)+β4Y40(θ)],

      (11)

      where β2 and β4 are the quadrupole and hexadecapole deformation parameters of daughter nuclei, respectively, with their values taken from Ref. [34]. The rms neutron and proton radii using the deformed nuclear distributions are calculated by

      R(D)τ[ρτ2(r,θ)r4drdΩρτ2(r,θ)r2drdΩ]12.

      (12)

      The multipole expansion of the deformed distribution ρτ2(r2) is used and then the orientation-dependent DF potential can be evaluated as the sum of different multipole components [18, 33]. The effect of nuclear deformations on the decay width could also be evaluated by integrating the partial width along the direction,

      Γ=π/20Γ(β)sinβdβ,

      (13)

      where β is the angle between the symmetry axis of the daughter nucleus and the direction of the α-particle emission.

      We studied two cases. In the first case, the radius and diffuseness parameters were obtained by HF calculations based on the Skyrme-SLy4 effective interaction that successfully explains the ground state properties of nuclei [3537]. The other case is the NST with the assumption rn>rp and an=ap=0.54 that the radius parameters rn and rp have determined in terms of the neutron and protons rms radii. The proton root-mean-square radius Rp is extracted by empirical data, and the neutron rms radius Rn is determined by the relation Rn=Rp+Δrnp, in which the Δrnp is obtained by

      Δrnp=(0.04±0.03)+(1.01±0.15)δ,

      (14)

      which presents the linear relationship between the isospin asymmetry δ=(NZ)/(N+Z) and Δrnp [38].

    3.   Results and discussion
    • For the nuclei with N>Z, the rms of the neutron density distribution is larger than the rms of the proton density distribution, hence the neutron density distribution can be different from the proton density distribution. Accordingly, to calculate the nucleus density distribution, it is not generally correct to assume that ρn=NZρp when N>Z. Instead, it is better to separately employ the neutron-to-proton densities.

      The ratio of the neutron-to-proton density versus distance from the nucleus center was plotted for some neutron-rich nuclei 210Pb, 218Po, 214Rn, 226Ra, 230Th, and 236U in Fig. 1. The ratio of densities is equal to N/Z in the range R=04 fm, followed by a sharp rise. This rising trend is more dramatic in HF compared to NST. To determine the effect of the neutron skin thickness on the α-decay, the neutron skin thickness was calculated by the NST and HF formalism versus the isospin asymmetry parameter δ compared in Fig. 2. This figure illustrates the variation between the calculated values of Δrnp within these two different approaches.

      Figure 1.  (color online) Ratio of neutron-to-proton density versus distance from nucleus center.

      Figure 2.  (color online) Neutron skin thickness Δrnp versus asymmetry parameter δ for isotope groups of elements Pb,Po,Rn,Ra,Th and U for two cases NST and HF.

      Because of the important role of the neutron and proton density distributions in calculating the nuclear potential, typically, the total potential for the 216Po, whose deformation parameters are zero, in the NST and HF cases were plotted in Fig. 3(a). In the range of R=9.912.5(fm), the potential value in the HF case is greater by 0.15% to 2.5% in comparison with the NST case. Because of the effect of the neutron skin thickness in NST and HF cases on the total potential, it is expected that the potential variations will affect the penetration probabilities and consequently affect the α-decay width. Therefore, the α-decay width was investigated in NST and HF cases for isotope groups of elements Pb,Po,Rn,Ra,Th, and U. Typically, the treatment of the calculated α-decay width for the Po isotope group in terms of the neutron number was presented in Fig. 3(b). In this figure, the α-decay widths of the NST case are greater by about 30% than their equivalent values of HF case. Therefore, this can affect the α-decay half-life and the extracted preformation factor. Within the calculated potentials, α-decay half-lives were calculated for the even-even isotope groups of elements Pb,Po,Rn,Ra,Th, and U by employing relations 6 to 8. In the NST case, the diffuseness parameters were fixed, and the radius parameters rn and rp were determined in terms of the neutron and proton rms radii. The rms proton radius Rp was taken from the experimental rms charge radii, and the rms neutron radius Rn was obtained from the relationship Rn=Rp+Δrnp. Since the experimental rms charge radii for some isotopes was not available, we have assumed the values of rms charge radii for these isotopes from a simple relationship, given by Rp=(r0+r1A2/3+r2A4/3)A1/3 with r0=0.891(2) fm, r1=1.52(3) fm, and r2=2.8(1) fm. This relationship was adapted from a systematic study on available data of the experimental rms charge radii [39].

      Figure 3.  (color online) (a) Total potential in terms of distance from center of nucleus for 216Po, assuming that HF and NST and (b) is the α-decay width versus neutron number for even-even Po isotopes.

      In Fig. 4, the logarithm of the preformation factor, calculated from the relation Pα=Γexpα/Γcalα [27], is plotted according to the parent neutron number for even-even isotopes with the atomic number in the range 82Z92. In general, the variation treatments in NST and HF cases are the same. Results indicate that the curve of HF has a value closer to the experimental data. The standard deviation from the experimental data obtained by σ={Ni=1[log10(Tical/Tiexp)]2/N}1/2 is equal to 0.438 and 0.391 for NST and HF cases, respectively. This confirms that the calculated half-lives in the HF case are closer to experimental data.

      Figure 4.  (color online) Logarithms of preformation factor of α-decay half-life versus parent nucleus neutron number, for even-even isotope groups of Pb,Po,Rn,Ra,Th, and U in HF and NST cases.

    4.   Summary
    • In this study, we investigated the effects of isospin asymmetries resulting from different proton and neutron density distributions on the total potentials, α-decay widths, and half-lives of even-even alpha emission nuclei with the atomic number in the range 82Z92. The neutron and proton density distributions of the daughter nuclei were calculated in the two-parameter Fermi form of nucleon density distributions. The radius and diffuseness parameter values of the nuclei were obtained from the NST and HF formalism, leading to different isospin asymmetries. The results indicated that the values of Δrnp for the two NST and HF cases affect the α-decay half-life and the preformation factor.

      The standard deviation of the calculated half-lives from experimental data reached about 0.438 for the NST case. However, the results illustrated that the isospin asymmetry within HF formalism lowered the standard deviation in comparison with its equivalent value obtained from NST. The standard deviation for the HF case is about 0.391, which indicates that the calculated half-life, in this case, is closer to the value obtained by experimental data.

Reference (39)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return