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

New look at Geiger-Nuttall law and α clustering of heavy nuclei

  • The Geiger-Nuttall (GN) law of α decay is commonly explained in terms of the quantum tunneling phenomenon. In this study, we show that such an explanation is actually not enough regarding the α particle clustering. Such a conclusion is drawn after exploring the involved coefficients of the GN law based on the conventional description of α decay, namely the formation of an α cluster and its subsequent penetration. The specific roles of the two former processes, in the GN law, manifest themselves via the systematical analysis of the calculated and experimental α decay half-lives versus the decay energies across the Z=82 and N=126 shell closures. The α-cluster preformation probability is then found to behave in a GN-like pattern. This previously ignored point is explicitly demonstrated as the product of an interplay between the mean-field and pairing effect, which in turn reveals the structural influence on the formation of the α cluster in a simple and clear manner. In addition to providing an effective method to evaluate the amount of surface α clustering in heavy nuclei, the present conjecture supports other theoretical treatments of the α preformation probability.
  • Dated back to the early stages of nuclear physics, Geiger and Nuttall discovered an empirical relationship between the traveling range of an emitted α particle from the parent nucleus and the period of transformation [1]. Accordingly, the partial half-life T1/2 of the α decay is related to the corresponding decay energy,

    log10T1/2=A(Z)Q1/2α+B(Z),

    (1)

    where the coefficients A and B, interpreted as functions of the atomic number Z, are determined by fitting the experimental data for each isotopic chain. This simple law has been confirmed since then and still holds generally for various isotopes, despite the increasing amount of α decay data from heavy and superheavy nuclei [2-8]. Consequently, a natural question arises: what is the physical mechanism behind this simple and reliable relationship? Leveraging the foundations of quantum mechanics, the α decay process was described in terms of the quantum tunneling phenomenon by Gamow [9] and independently by Condon and Gurney [10, 11] in 1928, properly explaining the linear dependence of Q1/2α in the above formula. This was a remarkable milestone in modern physics and the first successful application of the quantum theory to nuclear physics.

    Currently, the original GN law has been developed and extended to make it universal or adjust it to different external environments [12-19]. In an earlier study, Buck et al. pointed out that the emitted α particle should orbit around the core nucleus with different global quantum numbers via the GN plots of isotopic chains [12]. The effects of these quantum numbers were embedded in a new version of the GN law, successfully reproducing the α decay half-lives of various isotopes across the N=126 shell closure [14]. Based on the Thomas expression of α emission, Qi et al. explored the validity of the GN law and identified its microscopic basis [15]. Very recently, the debate over the influence of a strong electromagnetic field on the α decay process was also initiated from a proposed shifted GN law [16, 17]. It is well known that the α decay is usually considered as a two-simultaneous-step process of preformation of α cluster at the nuclear surface and its penetration through the Coulomb barrier. Given the success of the GN law, it is expected that these two physical steps respectively play an important role in this simple relationship (1). However, the GN plot has been always attributed to the penetration part since the Gamow scenario, while the α preformation probability is supposed to vary smoothly in the off-shell region. Although the experimental α decay half-lives can be well reproduced through this strategy, the systematic deviation between the calculated and measured half-lives still exists, not to mention the situation across the shell closure. Therefore, there must be something missing, rooted in the preformation of the α cluster. The present study precisely addresses this subject in view of the GN law, which can not only lead to a novel understanding of an old rule but also provide a new perspective of α clustering in heavy nuclei.

    The present challenge is how to designate the role of the two abovementioned processes that act in the GN law. Meanwhile, the fully microscopic description of the formation of the α cluster, in particular for a large range of nuclei, is an open problem [20]. Hence, a practical alternative is to present the systematics of calculated and experimental α decay half-lives versus Q1/2α in parallel, thereby leading to the indirect description of the behavior of the α preformation factor. To this end, the available α-decaying nuclei are divided into three regions:

    (i) N126, Z82;

    (ii) N126, Z>82;

    (iii) N>126, Z>82,

    guided by present experimental facts. Note that the N>126, Z82 region is not considered here, because there is only one α emitter identified up to now, namely 210Pb. Let us now recall the relationship between the decay half-life and the decay width in logarithm scale,

    log10T1/2=log10ln2PαΓ=log10Γlog10Pα+log10(ln2).

    (2)

    Evidently, the first two terms are responsible for the Q1/2α dependence in the GN law of one isotopic chain. As mentioned above, the Γ term or related penetration probability, is indeed crucial for such an idea, based on the quantum tunneling interpretation. The α preformation factor, especially in the logarithm scale of the GN law, is usually believed to vary quite limitedly as compared to the decay width or the penetrability. The GN law is therefore well understood within the tunneling explanation, while its involved coefficients are generally constant for a prescribed isotopic chain. Of physical interest is the deviation from this judgement, especially when it comes to the extreme proton-neutron ratios or various isotopes across the shell closure [8, 21, 22]. The implementation of the present objective starts from the calculation of partial decay width based on the Thomas expression [23],

    Γ=2kμR2ϕnj(R)2G(R)2+F(R)2,

    (3)

    where the wave number k = 2μQα/, and μ is the reduced mass AdAα/(Ad+Aα) in the unit of nucleon mass. The angular momentum , carried by the emitted α particle, equals zero for the present case as only ground-state (0+) to ground-state transitions of even-even α emitters are considered. R is the α-core relative distance, while ϕnj and G (F) are respectively the radial wave function of the α-core relative motion and the regular (irregular) Coulomb function. For a large enough R beyond the nuclear surface, the final Γ value is actually not dependent at all on the choice of R considering the asymptotic behavior of the former functions. In such an α cluster model, the fundamental ingredient is the α-core relative motion, which is constrained by the large global quantum number G. The internal nodes n in the above wave function are connected to this quantum number by the Wildermuth condition [12],

    G=2n+=4i=1gi,

    (4)

    presenting the effect of the Pauli principle; gi denotes here the oscillator quantum number of nucleons to form the α particle, and its value is set as gi= 4 for nuclei with 50Z,N82, gi= 5 for nuclei with 82<Z,N126, and gi= 6 for nuclei above the N= 126 neutron shell closure. In this sense, all the nucleons in the α cluster should occupy the states above the Fermi surface of the residual daughter nucleus, and the cluster-daughter scheme is then somewhat constructed on the building blocks of the shell model. The α decay width Γ can then be obtained by following the above procedure once the α-core potential is constructed via the double-folding integral (see details in Refs. [24-26]). The log10ln2Γ in the right side of Eq. (2) is here written as log10Tcalc1/2. Once the left side term in Eq. (2) is chosen as the experimental value, namely log10Texpt1/2, the Pα term can be extracted as the deviation between the log10Texpt1/2 and log10Tcalc1/2, namely

    log10Pα=log10Tcalc1/2log10Texpt1/2.

    (5)

    To analyze the contribution of the penetration part in the GN law, we calculated the α decay width of various even-even isotopes with N126. These calculated α decay half-lives follow the GN law with different coefficients A and B fitted for each isotopic chain, as shown in Fig. 1, including elements from hafnium to thorium. By contrast, for the next major shell above N=126, the global quantum number G changes from 20 to 22. This can result in a deviation from the initial linear dependence of log10Tcalc1/2 on Q1/2α to a certain extent, which is crucial in explaining the different GN plots for N<126 and N>126 [12]. To present this key factor of the α-cluster model, the calculated α decay half-lives are plotted as a function of Q1/2α for the whole Rn and Ra isotopic chains. As expected, two linear relations can be clearly distinguished for regions (ii) and (iii) in Fig. 2; note that the deviation between them is quite limited. However, the GN plots of experimental α decay data are believed to change in an evident manner especially when crossing the N=126 shell closure. Hence, one may doubt whether it is enough that the well-known GN law is entirely described in terms of the tunneling phenomenon.

    Figure 1

    Figure 1.  (color online) The logarithm of calculated α decay half-lives versus Q1/2α for the even-even Hf-Th isotopes, in which the preformation factor is excluded, namely Pα=1. Note that the straight line guides the GN law for each isotopic chain.

    Figure 2

    Figure 2.  Same as Fig. 1 but only for the Rn and Ra nuclei including the regions with N>126 (square) and N126 (circle). Note that the two straight lines correspond to these two regions with different G values.

    To clarify the above ambiguity, we systematically analyzed the GN relationship of log10Texpt1/2 values with the decay energy term Q1/2α for available even-even α emitters in the three regions mentioned. Accordingly, the coefficients A and B were determined for one isotopic chain in a prescribed region. The same procedure is also valid for the Tcalc1/2 case. All these coefficients A and B are presented as functions of the atomic number Z of parent nuclei in Fig. 3. For convenience, the comparative analyses of experimental data Texpt1/2 and calculated results Tcalc1/2 are respectively denoted as case I and case II. A clear point is that different sets of coefficients (divided by the vertical dotted line) in the GN plot are requested for regions (i) (with Z82) and (ii) (with Z>82), regardless of considering case I or II. This is expected owing to the shell effect on the α decay, which is not a new phenomenon. However, the unexpected and remarkable issue is that there is an unambiguous discrepancy between the coefficients A (and B) of case I and those of case II, regardless of the region considered. We may therefore conclude that tunneling, corresponding to case II, is definitely not the whole explanation for the α decay process. The only way to fill this gap should be a dependency on the α preformation probability. According to Eq. (5), the α preformation factor Pα is supposed to follow an exponential law with Q1/2α, i.e.,

    Figure 3

    Figure 3.  (color online) Coefficients A and B of even-even nuclei involved in the GN law versus the atomic number Z for regions (i) plus (ii) at the left panel and region (iii) at the right panel. The square indicates the computational case II, whereas the experimental case I is denoted by the circle. The dashed lines guide the eye for presenting the various dependences in different regions.

    log10Pα=A(Z)Q1/2α+B(Z).

    (6)

    Furthermore, by a careful analysis on the specific values of A and B for both cases I and II, one can combine two panels of this figure to illustrate the following two aspects about the preformation probability of the emitted α particle:

    i) The two terms on the right side of the above Pα expression appear to compete with each other considering their signs and absolute values. What are their microscopic basements?

    ii) The slope of linear relationship A or B appears to be different for the three focused regions. Moreover, the values of A and B in region (ii) are obviously larger than those in the two other regions.

    To answer these questions and gain a deep understanding of clustering at the surface of heavy nuclei, we recall the formation amplitude in the conventional formalism, namely [20, 27]

    F(R)=dˆRdξαdξd[Ψd(ξd)ϕα(ξα)Yl(ˆR)]Ψp,

    (7)

    where ξd and ξα are respectively the internal degrees of freedom for the daughter nucleus and the α particle, and Ψd(ξd) and Ψp(ξd,ξα,ˆR) denote the wave functions of the daughter and parent nuclei, respectively. The intrinsic α particle wave function behaves in an n=l=0 harmonic oscillator eigenstate [28]

    ϕα(ξα)=18(ναπ)9/4exp[να(r2nn+r2pp+r2np)/4]Sα,

    (8)

    with the relative distances rnn (rpp) of neutron-neutron (proton-proton) pairs plus the distance rnp between the mass centers of the nn and pp pairs. The quantity να is the αparticle harmonic oscillation parameter, and Sα is its spin-isospin function. Very recently, given the extreme case induced by four non-interacting protons and neutrons, the particle decay unit (p.d.u.) was defined to set a simple averaged single-particle limit for the αparticle formation amplitude [29],

    Fα;pdu=8ν3/4απ7/4R3,

    (9)

    when the ground-state to ground-state (g.s.) transitions of even-even nuclei are considered. The collectivity involved in the α decay process is then revealed via the ratio between the extracted α particle formation amplitude and the corresponding p.d.u. value. Inspired by this, we assume that the α preformation probability can be written as

    Pα=PcαP0α,

    (10)

    where P0α denotes the contribution of the "mean-field"-like part in terms of uncorrelated nucleons, corresponding to the |Fα;pdu|2 value. The former Pcα part corresponds to the enhancement factor due to the collectivity's degree of freedom. By contrast, as shown above, Fα;pdu actually relates to the nuclear volume, i.e., the mass number of the parent nuclei. With this in mind, the B(Z) term in Eq. (6) is expected to present the effect of single particle pattern, while A(Z)Q1/2α should be derived from the collective behavior of corresponding parent nuclei. The latter Q1/2α dependence of Pα can be unravelled in the shell-model context. It is believed that the clustering of α particles at the nuclear surface is governed by the pairing force acting among the involved neutrons and protons. For two particles in a non-degenerate system with a constant pairing strength G, the corresponding two-particle wave function amplitudes Xi should be related to the G value. By contrast, the α (four-particle) formation amplitude is connected to the Xi values. This can also be confirmed through the multistep shell-model method (MSM), where the MSM basis consists of the tensorial product of the two-particle correlated states [30]. The g.s. energy of even-even nuclei can be presented as a function of the pairing strength G, which determines the α decay energy via the difference between g.s. energies of parent and daughter nuclei. As a result, a connection between Qα and the final α preformation probability is expected, namely PcαXiGQα. In this sense, the pairing correlations can be understood as the origin of the term A(Z)Q1/2α in the Pα expression, i.e., Eq. (6).

    In a nutshell, the interplay between the mean-field and pairing effects leads to the final clustering of two neutrons and two protons at the nuclear surface, expressing a simple but meaningful relationship similar to the GN law. After the above discussion on issue i), the second issue can be explained as follows. The valence neutrons and protons are supposed to be located in the same major shell in region (ii) with 82<Z,N<126, leading to larger overlaps of single particle radial wave functions and relative stronger pairing correlation. Hence, the absolute values of coefficients A and B in region (ii) are relatively larger than those in the two other regions. Moreover, despite this discrepancy, the comparable competition between these two terms in Eq. (6) causes a smooth varying trend of Pα off the shell closure. This is the deep reason why the Geiger-Nuttal law can be successfully explained in terms of the quantum tunneling phenomenon while excluding the contribution of the α preformation process. Let us mention here that the four-nucleon correlations (quartetting) may play an important role in the collective dynamics of clustering [31-33]. Considering this, we plot the coefficients A and B in Eq. (6), generated by the difference between the systematics of experimental data and that of the calculated results based on Eq. (5), versus the number of quartets Nq in Fig. 4. The Nq value is determined by the number of proton pairs above the closest proton magic number Zmag, namely Nq=(ZZmag)/2, as suggested in Ref. [33]. This figure shows that the involved coefficients in the Pα systematics vary more evidently with the quartet numbers in region (ii) as compared to those in region (i). These α decaying nuclei in region (ii) are closer to the doubly-magic nuclei, implying that the large quartetting effect is present in emitters near the shell closure. This point is consistent with other studies [32, 33], and the quartetting effect is supposed to diminish greatly due to the Pauli principle when more valence nucleons are added.

    Figure 4

    Figure 4.  Coefficients A and B extracted from the Pα formula, i.e., Eq. (6), based on Eq. (5), versus the quartet number for α emitters in regions (ii) Z>82 and (i) Z82, respectively.

    Last but not least, we would like to discuss the consistency of the present conjecture about Pα with other theoretical studies. Recently, the cluster formation model (CFM) was proposed to evaluate the α preformation probability in an effective way [34, 35]. Within the CFM, the α formation amplitude is connected with the formation and decay energies, thereby supporting the above conclusion. Given the two-particle formation treatment in the BCS approach, the similarity between the neutron/proton pairing gaps and the α preformation amplitude is also presented [5], supporting the present GN-like law of Pα to some extent. Dated back to the earlier microscopic calculation on the absolute α decay width in the shell model plus the cluster component [36], the g.s. energy of 212Po, with respect to the four-nucleon threshold, was found to be very sensitive to the shell-model basis states, including the cluster states. This g.s. energy, after dividing the α particle energy, equals the α decay energy, while the α formation is regulated by the components of the state vector in the standard cluster basis. Consequently, it is expected that the final α preformation probability is related to the decay energy accompanied by the Z dependence. Moreover, the penetration process is well addressed regardless of whether the microscopic or semi-classical (such as the WKB approximation) approach is used, which is actually the basis for various phenomenological formulas of α decay half-lives after considering some extra factors for different types of parent nuclei [13, 18, 19]. However, through the present analysis, onecan find that the missing point is the similar behavior of the α preformation factor, as demonstrated here, which is another integral aspect for the success of these empirical evaluations.

    In conclusion, surface α-clustering is an important and profound phenomenon for both the cluster structure and decay properties of heavy neutron-deficient nuclei. We derived an effective relationship for the α-cluster preformation probability by carefully analyzing the famous Geiger-Nuttall law of α decay, thereby shedding new light on this topic. By comparably analyzing the involved coefficients of the GN law from the systematics of experimental and calculated α decay half-lives, it was found that the amount of surface α clustering should follow a GN-like pattern. This conjecture is requested by the full description of the well-known GN law rather than its simple description in terms of quantum tunneling. Within the shell-model context, the behavior of the α preformation factor is determined by the mean-field and pairing effects, which explicitly manifest themselves via the two Z dependent terms in the GN-like formula of Pα. Furthermore, the balance between these two components produces a relative smooth value of α preformation probability, leading to the success of most of the effective approaches on α decay. Moreover, the results of the present study are consistent with other available theoretical results on α clustering of heavy nuclei. Encouraged by this, we hope that the present derivation can offer a novel window to understand the α correlations when probing fundamental nuclear interactions.

    [1] H. Geiger and J. M. Nuttall, Philos. Mag. 22, 613 (1911) doi: 10.1080/14786441008637156
    [2] G. Audi, F. G. Kondev, M. Wang et al., Chin. Phys. C 41, 030001 (2017) doi: 10.1088/1674-1137/41/3/030001
    [3] M. Wang, G. Audi, F. G. Kondev et al., Chin. Phys. C 41, 030003 (2017) doi: 10.1088/1674-1137/41/3/030003
    [4] NNDC at the Brookhaven National Laboratory, ENSDF database http://www.nndc.bnl.gov, Last updated 3rd, September 2020
    [5] A. N. Andreyev, M. Huyse, P. Van Duppen et al., Phys. Rev. Lett. 110, 242502 (2013) doi: 10.1103/PhysRevLett.110.242502
    [6] Yu. Ts. Oganessian and V. K. Utyonkov, Rep. Prog. Phys. 78, 036301 (2015) doi: 10.1088/0034-4885/78/3/036301
    [7] S. Hofmann, J. Phys. G: Nucl. Part. Phys. 42, 114001 (2015) doi: 10.1088/0954-3899/42/11/114001
    [8] Z. Y. Zhang, Z. G. Gan, H. B. Yang et al., Phys. Rev, Lett. 122, 192503 (2019) doi: 10.1103/PhysRevLett.122.192503
    [9] G. Gamow, Z. Phys. 52, 510 (1929) doi: 10.1007/BF01339451
    [10] E. U. Condon and R. W. Gurney, Nature (London) 122, 438 (1928) doi: 10.1038/122438a0
    [11] R. W. Gurney and E. U. Condon, Phys. Rev. 33, 127 (1929) doi: 10.1103/PhysRev.33.127
    [12] B. Buck, A. C. Merchant, and S. M. Perez, Phys. Rev. Lett 65, 2975 (1990) doi: 10.1103/PhysRevLett.65.2975
    [13] D. N. Poenaru, R. A. Gherghescu, and W. Greiner, Phys. Rev. C 83, 014601 (2011) doi: 10.1103/PhysRevC.83.014601
    [14] Y. Ren and Z. Ren, Phys. Rev. C 85, 044608 (2012) doi: 10.1103/PhysRevC.85.044608
    [15] C. Qi, A. N. Andreyev, M. Huysed et al., Phys. Lett. B 734, 203 (2014) doi: 10.1016/j.physletb.2014.05.066
    [16] D. S. Delion and S. A. Ghinescu, Phys. Rev. Lett. 119, 202501 (2017) doi: 10.1103/PhysRevLett.119.202501
    [17] Adriana Pálffy and Sergey V. Popruzhenko, Phys. Rev. Lett. 124, 212505 (2020) doi: 10.1103/PhysRevLett.124.212505
    [18] Z. Z. Ren, C. Xu, and Z. J. Wang, Phys. Rev. C 70, 034304 (2004)
    [19] D. D. Ni, Z. Z. Ren, T. K. Dong et al., Phys. Rev. C 78, 044310 (2008) doi: 10.1103/PhysRevC.78.044310
    [20] C. Qi, R. Liotta, and R. Wyss, Prog. Part. Nucl. Phys. 105, 214 (2019) doi: 10.1016/j.ppnp.2018.11.003
    [21] K. Auranen, D. Seweryniak, M. Albers et al., Phys. Rev. Lett. 121, 182501 (2018) doi: 10.1103/PhysRevLett.121.182501
    [22] J. Khuyagbaatar, A. Yakushev, Ch. E. Düllmann et al., Phys. Rev. Lett. 115, 242502 (2015) doi: 10.1103/PhysRevLett.115.242502
    [23] R. G. Thomas, Prog. Theor. Phys. 12, 253 (1954) doi: 10.1143/PTP.12.253
    [24] A. M. Kobos, B. A. Brown, P. E. Hodgson et al., Nucl. Phys. A 384, 65 (1982) doi: 10.1016/0375-9474(82)90305-0
    [25] Y. B. Qian and Z. Z. Ren, Phys. Lett. B 777, 298 (2018) doi: 10.1016/j.physletb.2017.12.046
    [26] D. S. Delion, Z. Z. Ren, A. Dumitrescu et al., J. Phys. G: Nucl. Part. Phys. 45, 053001 (2018) doi: 10.1088/1361-6471/aaac52
    [27] R. G. Lovas, R. J. Liotta, A. Insolia et al., Phys. Rep. 294, 265 (1998) doi: 10.1016/S0370-1573(97)00049-5
    [28] I. Tonozuka and A. Arima, Nucl. Phys. A 323, 45 (1979) doi: 10.1016/0375-9474(79)90415-9
    [29] C. Qi, R. Liotta, and R. Wyss, arXiv: 2003.00859v1
    [30] G. Dodig-crnković, F. A. Janouch, and R. J. Liotta, Nucl. Phys. A 501, 533 (1989) doi: 10.1016/0375-9474(89)90146-2
    [31] G. Röpke, P. Schuck, Y. Funaki et al., Phys. Rev. C 90, 034304 (2014) doi: 10.1103/PhysRevC.90.034304
    [32] S. Yang, C. Xu, G. Röpke et al., Phys. Rev. C 101, 024316 (2020)
    [33] D. S. Delion and A. Dumitrescu, Phys. Rev. C 102, 014327 (2020) doi: 10.1103/PhysRevC.102.014327
    [34] S. M S. Ahmed, R. Yahaya, S. Radiman et al., J. Phys. G: Nucl. Part. Phys. 40, 065105 (2013) doi: 10.1088/0954-3899/40/6/065105
    [35] D. M. Deng, Z. Z. Ren, D. D. Ni et al., J. Phys. G: Nucl. Part. Phys. 42, 075106 (2015) doi: 10.1088/0954-3899/42/7/075106
    [36] K. Varga, R. G. Lovas, and R. J. Liotta, Phys. Rev. Lett. 69, 37 (1992) doi: 10.1103/PhysRevLett.69.37
  • [1] H. Geiger and J. M. Nuttall, Philos. Mag. 22, 613 (1911) doi: 10.1080/14786441008637156
    [2] G. Audi, F. G. Kondev, M. Wang et al., Chin. Phys. C 41, 030001 (2017) doi: 10.1088/1674-1137/41/3/030001
    [3] M. Wang, G. Audi, F. G. Kondev et al., Chin. Phys. C 41, 030003 (2017) doi: 10.1088/1674-1137/41/3/030003
    [4] NNDC at the Brookhaven National Laboratory, ENSDF database http://www.nndc.bnl.gov, Last updated 3rd, September 2020
    [5] A. N. Andreyev, M. Huyse, P. Van Duppen et al., Phys. Rev. Lett. 110, 242502 (2013) doi: 10.1103/PhysRevLett.110.242502
    [6] Yu. Ts. Oganessian and V. K. Utyonkov, Rep. Prog. Phys. 78, 036301 (2015) doi: 10.1088/0034-4885/78/3/036301
    [7] S. Hofmann, J. Phys. G: Nucl. Part. Phys. 42, 114001 (2015) doi: 10.1088/0954-3899/42/11/114001
    [8] Z. Y. Zhang, Z. G. Gan, H. B. Yang et al., Phys. Rev, Lett. 122, 192503 (2019) doi: 10.1103/PhysRevLett.122.192503
    [9] G. Gamow, Z. Phys. 52, 510 (1929) doi: 10.1007/BF01339451
    [10] E. U. Condon and R. W. Gurney, Nature (London) 122, 438 (1928) doi: 10.1038/122438a0
    [11] R. W. Gurney and E. U. Condon, Phys. Rev. 33, 127 (1929) doi: 10.1103/PhysRev.33.127
    [12] B. Buck, A. C. Merchant, and S. M. Perez, Phys. Rev. Lett 65, 2975 (1990) doi: 10.1103/PhysRevLett.65.2975
    [13] D. N. Poenaru, R. A. Gherghescu, and W. Greiner, Phys. Rev. C 83, 014601 (2011) doi: 10.1103/PhysRevC.83.014601
    [14] Y. Ren and Z. Ren, Phys. Rev. C 85, 044608 (2012) doi: 10.1103/PhysRevC.85.044608
    [15] C. Qi, A. N. Andreyev, M. Huysed et al., Phys. Lett. B 734, 203 (2014) doi: 10.1016/j.physletb.2014.05.066
    [16] D. S. Delion and S. A. Ghinescu, Phys. Rev. Lett. 119, 202501 (2017) doi: 10.1103/PhysRevLett.119.202501
    [17] Adriana Pálffy and Sergey V. Popruzhenko, Phys. Rev. Lett. 124, 212505 (2020) doi: 10.1103/PhysRevLett.124.212505
    [18] Z. Z. Ren, C. Xu, and Z. J. Wang, Phys. Rev. C 70, 034304 (2004)
    [19] D. D. Ni, Z. Z. Ren, T. K. Dong et al., Phys. Rev. C 78, 044310 (2008) doi: 10.1103/PhysRevC.78.044310
    [20] C. Qi, R. Liotta, and R. Wyss, Prog. Part. Nucl. Phys. 105, 214 (2019) doi: 10.1016/j.ppnp.2018.11.003
    [21] K. Auranen, D. Seweryniak, M. Albers et al., Phys. Rev. Lett. 121, 182501 (2018) doi: 10.1103/PhysRevLett.121.182501
    [22] J. Khuyagbaatar, A. Yakushev, Ch. E. Düllmann et al., Phys. Rev. Lett. 115, 242502 (2015) doi: 10.1103/PhysRevLett.115.242502
    [23] R. G. Thomas, Prog. Theor. Phys. 12, 253 (1954) doi: 10.1143/PTP.12.253
    [24] A. M. Kobos, B. A. Brown, P. E. Hodgson et al., Nucl. Phys. A 384, 65 (1982) doi: 10.1016/0375-9474(82)90305-0
    [25] Y. B. Qian and Z. Z. Ren, Phys. Lett. B 777, 298 (2018) doi: 10.1016/j.physletb.2017.12.046
    [26] D. S. Delion, Z. Z. Ren, A. Dumitrescu et al., J. Phys. G: Nucl. Part. Phys. 45, 053001 (2018) doi: 10.1088/1361-6471/aaac52
    [27] R. G. Lovas, R. J. Liotta, A. Insolia et al., Phys. Rep. 294, 265 (1998) doi: 10.1016/S0370-1573(97)00049-5
    [28] I. Tonozuka and A. Arima, Nucl. Phys. A 323, 45 (1979) doi: 10.1016/0375-9474(79)90415-9
    [29] C. Qi, R. Liotta, and R. Wyss, arXiv: 2003.00859v1
    [30] G. Dodig-crnković, F. A. Janouch, and R. J. Liotta, Nucl. Phys. A 501, 533 (1989) doi: 10.1016/0375-9474(89)90146-2
    [31] G. Röpke, P. Schuck, Y. Funaki et al., Phys. Rev. C 90, 034304 (2014) doi: 10.1103/PhysRevC.90.034304
    [32] S. Yang, C. Xu, G. Röpke et al., Phys. Rev. C 101, 024316 (2020)
    [33] D. S. Delion and A. Dumitrescu, Phys. Rev. C 102, 014327 (2020) doi: 10.1103/PhysRevC.102.014327
    [34] S. M S. Ahmed, R. Yahaya, S. Radiman et al., J. Phys. G: Nucl. Part. Phys. 40, 065105 (2013) doi: 10.1088/0954-3899/40/6/065105
    [35] D. M. Deng, Z. Z. Ren, D. D. Ni et al., J. Phys. G: Nucl. Part. Phys. 42, 075106 (2015) doi: 10.1088/0954-3899/42/7/075106
    [36] K. Varga, R. G. Lovas, and R. J. Liotta, Phys. Rev. Lett. 69, 37 (1992) doi: 10.1103/PhysRevLett.69.37
  • 加载中

Cited by

1. Chandran, M., Anjali, V.K., Santhosh, K.P. Systematic study on heavy-particle radioactivity of superheavy nuclei 297–300119[J]. Nuclear Science and Techniques, 2024, 35(11): 196. doi: 10.1007/s41365-024-01558-3
2. Khuyagbaatar, J.. Ternary fission with the emission of long-range α particles in fission of the heaviest nuclei[J]. Physical Review C, 2024, 110(1): 014311. doi: 10.1103/PhysRevC.110.014311
3. Majekodunmi, J.T., Bhuyan, M., Anwar, K. et al. Cluster decay dynamics of actinides yielding non-Pb-daughter within relativistic mean field formalism[J]. European Physical Journal A, 2024, 60(5): 101. doi: 10.1140/epja/s10050-024-01324-4
4. Tang, S., Qian, Y., Ren, Z. Model-independent analysis on the regular behavior of α preformation probability in heavy nuclei[J]. Physical Review C, 2023, 108(6): 064303. doi: 10.1103/PhysRevC.108.064303
5. Blank, B., Page, R.D. Charged-Particle Radioactive Decays[J]. Handbook of Nuclear Physics: Volume 1-5, 2023. doi: 10.1007/978-981-19-6345-2_44
6. 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
7. Manjunatha, H.C., Vidya, Y.S., Gupta, P.S.D. et al. An Attempt to Unify Half-Lives of Proton Decay, Beta Decay and Electron Capture Through GN Law[J]. Iranian Journal of Science and Technology, Transaction A: Science, 2022, 46(4): 1301-1308. doi: 10.1007/s40995-022-01319-0
8. Deng, J.-G., Zhang, H.-F., Sun, X.-D. New behaviors of α-particle preformation factors near doubly magic 100Sn[J]. Chinese Physics C, 2022, 46(6): 061001. doi: 10.1088/1674-1137/ac5a9f
9. Fan, J., Xu, C. Exploring the half-lives of extremely long-lived α emitters[J]. Chinese Physics C, 2022, 46(5): 054105. doi: 10.1088/1674-1137/ac500d
10. Ghorbani, F., Alavi, S.A., Dehghani, V. Alpha decay half-lives of even-even nuclei using alpha-folding interaction[J]. European Physical Journal A, 2022, 58(2): 12. doi: 10.1140/epja/s10050-022-00670-5
11. Zou, Y.-T., Pan, X., Li, X.-H. et al. Systematic study of two-proton radioactivity with a screened electrostatic barrier[J]. Chinese Physics C, 2021, 45(10) doi: 10.1088/1674-1137/ac1b96
12. Soylu, A., Qi, C. Extended universal decay law formula for the α and cluster decays[J]. Nuclear Physics A, 2021. doi: 10.1016/j.nuclphysa.2021.122221
13. Zou, Y.-T., Pan, X., Liu, H.-M. et al. Systematic studies on a decay half-lives of neptunium isotopes[J]. Physica Scripta, 2021, 96(7): 075301. doi: 10.1088/1402-4896/abf795

Figures(4)

Get Citation
Yi-Bin Qian and Zhong-Zhou Ren. New look at Geiger-Nuttall Law and α clustering of heavy nuclei[J]. Chinese Physics C. doi: 10.1088/1674-1137/abce14
Yi-Bin Qian and Zhong-Zhou Ren. New look at Geiger-Nuttall Law and α clustering of heavy nuclei[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abce14 shu
Milestone
Received: 2020-10-09
Article Metric

Article Views(2962)
PDF Downloads(84)
Cited by(13)
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:

New look at Geiger-Nuttall law and α clustering of heavy nuclei

    Corresponding author: Yi-Bin Qian, qyibin@njust.edu.cn
    Corresponding author: Zhong-Zhou Ren, zren@nju.edu.cn
  • 1. Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, China
  • 2. School of Physics Science and Engineering, Tongji University, Shanghai 200092, China

Abstract: The Geiger-Nuttall (GN) law of α decay is commonly explained in terms of the quantum tunneling phenomenon. In this study, we show that such an explanation is actually not enough regarding the α particle clustering. Such a conclusion is drawn after exploring the involved coefficients of the GN law based on the conventional description of α decay, namely the formation of an α cluster and its subsequent penetration. The specific roles of the two former processes, in the GN law, manifest themselves via the systematical analysis of the calculated and experimental α decay half-lives versus the decay energies across the Z=82 and N=126 shell closures. The α-cluster preformation probability is then found to behave in a GN-like pattern. This previously ignored point is explicitly demonstrated as the product of an interplay between the mean-field and pairing effect, which in turn reveals the structural influence on the formation of the α cluster in a simple and clear manner. In addition to providing an effective method to evaluate the amount of surface α clustering in heavy nuclei, the present conjecture supports other theoretical treatments of the α preformation probability.

    HTML

    I.   INTRODUCTION
    • Dated back to the early stages of nuclear physics, Geiger and Nuttall discovered an empirical relationship between the traveling range of an emitted α particle from the parent nucleus and the period of transformation [1]. Accordingly, the partial half-life T1/2 of the α decay is related to the corresponding decay energy,

      log10T1/2=A(Z)Q1/2α+B(Z),

      (1)

      where the coefficients A and B, interpreted as functions of the atomic number Z, are determined by fitting the experimental data for each isotopic chain. This simple law has been confirmed since then and still holds generally for various isotopes, despite the increasing amount of α decay data from heavy and superheavy nuclei [2-8]. Consequently, a natural question arises: what is the physical mechanism behind this simple and reliable relationship? Leveraging the foundations of quantum mechanics, the α decay process was described in terms of the quantum tunneling phenomenon by Gamow [9] and independently by Condon and Gurney [10, 11] in 1928, properly explaining the linear dependence of Q1/2α in the above formula. This was a remarkable milestone in modern physics and the first successful application of the quantum theory to nuclear physics.

      Currently, the original GN law has been developed and extended to make it universal or adjust it to different external environments [12-19]. In an earlier study, Buck et al. pointed out that the emitted α particle should orbit around the core nucleus with different global quantum numbers via the GN plots of isotopic chains [12]. The effects of these quantum numbers were embedded in a new version of the GN law, successfully reproducing the α decay half-lives of various isotopes across the N=126 shell closure [14]. Based on the Thomas expression of α emission, Qi et al. explored the validity of the GN law and identified its microscopic basis [15]. Very recently, the debate over the influence of a strong electromagnetic field on the α decay process was also initiated from a proposed shifted GN law [16, 17]. It is well known that the α decay is usually considered as a two-simultaneous-step process of preformation of α cluster at the nuclear surface and its penetration through the Coulomb barrier. Given the success of the GN law, it is expected that these two physical steps respectively play an important role in this simple relationship (1). However, the GN plot has been always attributed to the penetration part since the Gamow scenario, while the α preformation probability is supposed to vary smoothly in the off-shell region. Although the experimental α decay half-lives can be well reproduced through this strategy, the systematic deviation between the calculated and measured half-lives still exists, not to mention the situation across the shell closure. Therefore, there must be something missing, rooted in the preformation of the α cluster. The present study precisely addresses this subject in view of the GN law, which can not only lead to a novel understanding of an old rule but also provide a new perspective of α clustering in heavy nuclei.

    II.   THEORETICAL FRAMEWORK FOR α DECAY WIDTH
    • The present challenge is how to designate the role of the two abovementioned processes that act in the GN law. Meanwhile, the fully microscopic description of the formation of the α cluster, in particular for a large range of nuclei, is an open problem [20]. Hence, a practical alternative is to present the systematics of calculated and experimental α decay half-lives versus Q1/2α in parallel, thereby leading to the indirect description of the behavior of the α preformation factor. To this end, the available α-decaying nuclei are divided into three regions:

      (i) N126, Z82;

      (ii) N126, Z>82;

      (iii) N>126, Z>82,

      guided by present experimental facts. Note that the N>126, Z82 region is not considered here, because there is only one α emitter identified up to now, namely 210Pb. Let us now recall the relationship between the decay half-life and the decay width in logarithm scale,

      log10T1/2=log10ln2PαΓ=log10Γlog10Pα+log10(ln2).

      (2)

      Evidently, the first two terms are responsible for the Q1/2α dependence in the GN law of one isotopic chain. As mentioned above, the Γ term or related penetration probability, is indeed crucial for such an idea, based on the quantum tunneling interpretation. The α preformation factor, especially in the logarithm scale of the GN law, is usually believed to vary quite limitedly as compared to the decay width or the penetrability. The GN law is therefore well understood within the tunneling explanation, while its involved coefficients are generally constant for a prescribed isotopic chain. Of physical interest is the deviation from this judgement, especially when it comes to the extreme proton-neutron ratios or various isotopes across the shell closure [8, 21, 22]. The implementation of the present objective starts from the calculation of partial decay width based on the Thomas expression [23],

      Γ=2kμR2ϕnj(R)2G(R)2+F(R)2,

      (3)

      where the wave number k = 2μQα/, and μ is the reduced mass AdAα/(Ad+Aα) in the unit of nucleon mass. The angular momentum , carried by the emitted α particle, equals zero for the present case as only ground-state (0+) to ground-state transitions of even-even α emitters are considered. R is the α-core relative distance, while ϕnj and G (F) are respectively the radial wave function of the α-core relative motion and the regular (irregular) Coulomb function. For a large enough R beyond the nuclear surface, the final Γ value is actually not dependent at all on the choice of R considering the asymptotic behavior of the former functions. In such an α cluster model, the fundamental ingredient is the α-core relative motion, which is constrained by the large global quantum number G. The internal nodes n in the above wave function are connected to this quantum number by the Wildermuth condition [12],

      G=2n+=4i=1gi,

      (4)

      presenting the effect of the Pauli principle; gi denotes here the oscillator quantum number of nucleons to form the α particle, and its value is set as gi= 4 for nuclei with 50Z,N82, gi= 5 for nuclei with 82<Z,N126, and gi= 6 for nuclei above the N= 126 neutron shell closure. In this sense, all the nucleons in the α cluster should occupy the states above the Fermi surface of the residual daughter nucleus, and the cluster-daughter scheme is then somewhat constructed on the building blocks of the shell model. The α decay width Γ can then be obtained by following the above procedure once the α-core potential is constructed via the double-folding integral (see details in Refs. [24-26]). The log10ln2Γ in the right side of Eq. (2) is here written as log10Tcalc1/2. Once the left side term in Eq. (2) is chosen as the experimental value, namely log10Texpt1/2, the Pα term can be extracted as the deviation between the log10Texpt1/2 and log10Tcalc1/2, namely

      log10Pα=log10Tcalc1/2log10Texpt1/2.

      (5)
    III.   NEW PERSPECTIVE FROM THE GEIGER-NUTTALL PLOT AND ORIGIN OF THE CLUSTER FORMATION
    • To analyze the contribution of the penetration part in the GN law, we calculated the α decay width of various even-even isotopes with N126. These calculated α decay half-lives follow the GN law with different coefficients A and B fitted for each isotopic chain, as shown in Fig. 1, including elements from hafnium to thorium. By contrast, for the next major shell above N=126, the global quantum number G changes from 20 to 22. This can result in a deviation from the initial linear dependence of log10Tcalc1/2 on Q1/2α to a certain extent, which is crucial in explaining the different GN plots for N<126 and N>126 [12]. To present this key factor of the α-cluster model, the calculated α decay half-lives are plotted as a function of Q1/2α for the whole Rn and Ra isotopic chains. As expected, two linear relations can be clearly distinguished for regions (ii) and (iii) in Fig. 2; note that the deviation between them is quite limited. However, the GN plots of experimental α decay data are believed to change in an evident manner especially when crossing the N=126 shell closure. Hence, one may doubt whether it is enough that the well-known GN law is entirely described in terms of the tunneling phenomenon.

      Figure 1.  (color online) The logarithm of calculated α decay half-lives versus Q1/2α for the even-even Hf-Th isotopes, in which the preformation factor is excluded, namely Pα=1. Note that the straight line guides the GN law for each isotopic chain.

      Figure 2.  Same as Fig. 1 but only for the Rn and Ra nuclei including the regions with N>126 (square) and N126 (circle). Note that the two straight lines correspond to these two regions with different G values.

      To clarify the above ambiguity, we systematically analyzed the GN relationship of log10Texpt1/2 values with the decay energy term Q1/2α for available even-even α emitters in the three regions mentioned. Accordingly, the coefficients A and B were determined for one isotopic chain in a prescribed region. The same procedure is also valid for the Tcalc1/2 case. All these coefficients A and B are presented as functions of the atomic number Z of parent nuclei in Fig. 3. For convenience, the comparative analyses of experimental data Texpt1/2 and calculated results Tcalc1/2 are respectively denoted as case I and case II. A clear point is that different sets of coefficients (divided by the vertical dotted line) in the GN plot are requested for regions (i) (with Z82) and (ii) (with Z>82), regardless of considering case I or II. This is expected owing to the shell effect on the α decay, which is not a new phenomenon. However, the unexpected and remarkable issue is that there is an unambiguous discrepancy between the coefficients A (and B) of case I and those of case II, regardless of the region considered. We may therefore conclude that tunneling, corresponding to case II, is definitely not the whole explanation for the α decay process. The only way to fill this gap should be a dependency on the α preformation probability. According to Eq. (5), the α preformation factor Pα is supposed to follow an exponential law with Q1/2α, i.e.,

      Figure 3.  (color online) Coefficients A and B of even-even nuclei involved in the GN law versus the atomic number Z for regions (i) plus (ii) at the left panel and region (iii) at the right panel. The square indicates the computational case II, whereas the experimental case I is denoted by the circle. The dashed lines guide the eye for presenting the various dependences in different regions.

      log10Pα=A(Z)Q1/2α+B(Z).

      (6)

      Furthermore, by a careful analysis on the specific values of A and B for both cases I and II, one can combine two panels of this figure to illustrate the following two aspects about the preformation probability of the emitted α particle:

      i) The two terms on the right side of the above Pα expression appear to compete with each other considering their signs and absolute values. What are their microscopic basements?

      ii) The slope of linear relationship A or B appears to be different for the three focused regions. Moreover, the values of A and B in region (ii) are obviously larger than those in the two other regions.

      To answer these questions and gain a deep understanding of clustering at the surface of heavy nuclei, we recall the formation amplitude in the conventional formalism, namely [20, 27]

      F(R)=dˆRdξαdξd[Ψd(ξd)ϕα(ξα)Yl(ˆR)]Ψp,

      (7)

      where ξd and ξα are respectively the internal degrees of freedom for the daughter nucleus and the α particle, and Ψd(ξd) and Ψp(ξd,ξα,ˆR) denote the wave functions of the daughter and parent nuclei, respectively. The intrinsic α particle wave function behaves in an n=l=0 harmonic oscillator eigenstate [28]

      ϕα(ξα)=18(ναπ)9/4exp[να(r2nn+r2pp+r2np)/4]Sα,

      (8)

      with the relative distances rnn (rpp) of neutron-neutron (proton-proton) pairs plus the distance rnp between the mass centers of the nn and pp pairs. The quantity να is the αparticle harmonic oscillation parameter, and Sα is its spin-isospin function. Very recently, given the extreme case induced by four non-interacting protons and neutrons, the particle decay unit (p.d.u.) was defined to set a simple averaged single-particle limit for the αparticle formation amplitude [29],

      Fα;pdu=8ν3/4απ7/4R3,

      (9)

      when the ground-state to ground-state (g.s.) transitions of even-even nuclei are considered. The collectivity involved in the α decay process is then revealed via the ratio between the extracted α particle formation amplitude and the corresponding p.d.u. value. Inspired by this, we assume that the α preformation probability can be written as

      Pα=PcαP0α,

      (10)

      where P0α denotes the contribution of the "mean-field"-like part in terms of uncorrelated nucleons, corresponding to the |Fα;pdu|2 value. The former Pcα part corresponds to the enhancement factor due to the collectivity's degree of freedom. By contrast, as shown above, Fα;pdu actually relates to the nuclear volume, i.e., the mass number of the parent nuclei. With this in mind, the B(Z) term in Eq. (6) is expected to present the effect of single particle pattern, while A(Z)Q1/2α should be derived from the collective behavior of corresponding parent nuclei. The latter Q1/2α dependence of Pα can be unravelled in the shell-model context. It is believed that the clustering of α particles at the nuclear surface is governed by the pairing force acting among the involved neutrons and protons. For two particles in a non-degenerate system with a constant pairing strength G, the corresponding two-particle wave function amplitudes Xi should be related to the G value. By contrast, the α (four-particle) formation amplitude is connected to the Xi values. This can also be confirmed through the multistep shell-model method (MSM), where the MSM basis consists of the tensorial product of the two-particle correlated states [30]. The g.s. energy of even-even nuclei can be presented as a function of the pairing strength G, which determines the α decay energy via the difference between g.s. energies of parent and daughter nuclei. As a result, a connection between Qα and the final α preformation probability is expected, namely PcαXiGQα. In this sense, the pairing correlations can be understood as the origin of the term A(Z)Q1/2α in the Pα expression, i.e., Eq. (6).

      In a nutshell, the interplay between the mean-field and pairing effects leads to the final clustering of two neutrons and two protons at the nuclear surface, expressing a simple but meaningful relationship similar to the GN law. After the above discussion on issue i), the second issue can be explained as follows. The valence neutrons and protons are supposed to be located in the same major shell in region (ii) with 82<Z,N<126, leading to larger overlaps of single particle radial wave functions and relative stronger pairing correlation. Hence, the absolute values of coefficients A and B in region (ii) are relatively larger than those in the two other regions. Moreover, despite this discrepancy, the comparable competition between these two terms in Eq. (6) causes a smooth varying trend of Pα off the shell closure. This is the deep reason why the Geiger-Nuttal law can be successfully explained in terms of the quantum tunneling phenomenon while excluding the contribution of the α preformation process. Let us mention here that the four-nucleon correlations (quartetting) may play an important role in the collective dynamics of clustering [31-33]. Considering this, we plot the coefficients A and B in Eq. (6), generated by the difference between the systematics of experimental data and that of the calculated results based on Eq. (5), versus the number of quartets Nq in Fig. 4. The Nq value is determined by the number of proton pairs above the closest proton magic number Zmag, namely Nq=(ZZmag)/2, as suggested in Ref. [33]. This figure shows that the involved coefficients in the Pα systematics vary more evidently with the quartet numbers in region (ii) as compared to those in region (i). These α decaying nuclei in region (ii) are closer to the doubly-magic nuclei, implying that the large quartetting effect is present in emitters near the shell closure. This point is consistent with other studies [32, 33], and the quartetting effect is supposed to diminish greatly due to the Pauli principle when more valence nucleons are added.

      Figure 4.  Coefficients A and B extracted from the Pα formula, i.e., Eq. (6), based on Eq. (5), versus the quartet number for α emitters in regions (ii) Z>82 and (i) Z82, respectively.

      Last but not least, we would like to discuss the consistency of the present conjecture about Pα with other theoretical studies. Recently, the cluster formation model (CFM) was proposed to evaluate the α preformation probability in an effective way [34, 35]. Within the CFM, the α formation amplitude is connected with the formation and decay energies, thereby supporting the above conclusion. Given the two-particle formation treatment in the BCS approach, the similarity between the neutron/proton pairing gaps and the α preformation amplitude is also presented [5], supporting the present GN-like law of Pα to some extent. Dated back to the earlier microscopic calculation on the absolute α decay width in the shell model plus the cluster component [36], the g.s. energy of 212Po, with respect to the four-nucleon threshold, was found to be very sensitive to the shell-model basis states, including the cluster states. This g.s. energy, after dividing the α particle energy, equals the α decay energy, while the α formation is regulated by the components of the state vector in the standard cluster basis. Consequently, it is expected that the final α preformation probability is related to the decay energy accompanied by the Z dependence. Moreover, the penetration process is well addressed regardless of whether the microscopic or semi-classical (such as the WKB approximation) approach is used, which is actually the basis for various phenomenological formulas of α decay half-lives after considering some extra factors for different types of parent nuclei [13, 18, 19]. However, through the present analysis, onecan find that the missing point is the similar behavior of the α preformation factor, as demonstrated here, which is another integral aspect for the success of these empirical evaluations.

    IV.   SUMMARY
    • In conclusion, surface α-clustering is an important and profound phenomenon for both the cluster structure and decay properties of heavy neutron-deficient nuclei. We derived an effective relationship for the α-cluster preformation probability by carefully analyzing the famous Geiger-Nuttall law of α decay, thereby shedding new light on this topic. By comparably analyzing the involved coefficients of the GN law from the systematics of experimental and calculated α decay half-lives, it was found that the amount of surface α clustering should follow a GN-like pattern. This conjecture is requested by the full description of the well-known GN law rather than its simple description in terms of quantum tunneling. Within the shell-model context, the behavior of the α preformation factor is determined by the mean-field and pairing effects, which explicitly manifest themselves via the two Z dependent terms in the GN-like formula of Pα. Furthermore, the balance between these two components produces a relative smooth value of α preformation probability, leading to the success of most of the effective approaches on α decay. Moreover, the results of the present study are consistent with other available theoretical results on α clustering of heavy nuclei. Encouraged by this, we hope that the present derivation can offer a novel window to understand the α correlations when probing fundamental nuclear interactions.

Reference (36)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return