Author
https://www.sciencedirect.com/science/article/pii/S0370269310003734
Received 22 May 2009, Revised 12 March 2010, Accepted 16 March 2010, Available online 18 March 2010.
Editor: W. Haxton
Keywords
Superdeformed decay
Two-state model of SD decay
Quantum-mechanical path integrals
1. Introduction
Since their first experimental observation [1], superdeformed (SD) nuclear states, with their strong ellipsoidal deformation and special set of shell closures, have offered a tantalizing and unique window into subatomic physics. Their rapid decay-out, in particular, has been the subject of great interest (e.g., Refs. [2–21]). In the standard theoretical approach [2,3], this process is modeled by a two-well potential function of deformation: Here, the nucleus is a single quantum mechanical particle, which tunnels between the two wells, and can escape the system via electromagnetically induced decay from either. Because the barrier between the SD and normal-deformed (ND) wells is a direct consequence of nucleon–nucleon interactions, an understanding of its shape for various nuclei and angular momenta would be of considerable importance to the study of nuclear structure. Thus, a common objective of theoretical studies is to bridge the gap between measured experimental data, such as lifetimes and nuclear spins, and the shape of this barrier. In this Letter, we show that the rapid decrease in barrier height with decreasing nuclear spin explains the SD decay mechanism.
It was previously shown [15,21] that an elegant, two-state model [11] of SD decay-out is sufficient to give an excellent picture of the system's time-evolution. One of the principal advantages of such a straightforward technique was the extraction from experiment of such important quantities as the tunneling matrix element V and the net rate of tunneling through the barrier . The purpose of the present Letter is to move beyond these phenomenological quantities, and extract the height of the barrier itself as a function of nucleus and nuclear spin, as shown in Fig. 1. Previous approaches [8] used a semiclassical treatment that did not allow for a systematic computation of the tunneling rate prefactor. As we will see, this can introduce significant deviation into physically meaningful quantities, such as the rate of decay-out. In this Letter, we go beyond previous treatments by computing the tunneling rate in a systematic and controlled fashion using a functional integral approach [22].
2. Path-integral approach to tunneling
Absent additional information about the nuclear shell potential, the smoothest potential describing both the tunnel barrier and the SD well (which lies far above the ND yrast line at decay out) is a cubic polynomial:
(1)
where x is a coordinate describing the quadrupole deformation of the nucleus ( corresponds to the bottom of the SD well, while the ND well occurs for ), ω is the oscillator frequency of the SD well, and is the inertia of the quadrupole vibrational mode.
In the absence of coupling to the electromagnetic field, the two-state model's Hamiltonian possesses only two real parameters [11]. Thus, when connecting to this model, it is appropriate that also has only two effective degrees of freedom, and λ. This phenomenological approach is to be contrasted with a microscopic approach [5–7,12], which may yield a more complex potential and mass tensor. In other words, a suitable choice of the two effective parameters and λ provides a maximum entropy (i.e., containing the least possible externally assigned information) approximation to the unknown nuclear potential barrier. Note that for this simple potential, the barrier frequency .
The use of Euclidean complex-time path integrals over the tunneling coordinate provides a systematic approach for determining quantum tunnelling rates at arbitrary temperature and dissipation [22], as an asymptotic expansion in ℏ. This method allows calculation not only of the leading-order exponential dependence of the tunneling rate on potential parameters, but also the more computationally difficult subdominant asymptotics (i.e., the prefactor term). For incoherent tunneling through the potential barrier (1) at zero temperature and friction [23–26], the mean-square tunneling matrix element out of the SD state is found to be:
(2)
where is the barrier height (as measured from the bottom of the SD well), and is the mean level spacing in the ND band at the energy of the decaying SD state. The action S to tunnel out of the SD state through the barrier is proportional to the argument of the exponential function in Eq. (2):
(3)
where the factor of two is due to the power of V.
To make contact with experiment, the tunneling matrix element may be estimated using the two-state model of SD decay [11,15,20], which assumes the decay-out process is dominated by coupling of each SD state with its nearest-lying energy level in the ND band. The branching ratios and for decay out and intraband decay, respectively, are determined by three rates [15]:
(4)
where and are the electromagnetic decay rates of the SD and ND states, respectively, and is the nucleus' net tunneling rate through the barrier, with , and Δ the energy difference between the unperturbed (i.e., ) SD and ND states. Given the experimentally determined branching ratios and the electromagnetic widths, is known [20]:
(5)
The tunneling matrix element V may then be determined statistically [15,20], assuming the SD and ND levels are uncorrelated, and that the ND levels obey the Wigner surmise. The mean-square tunneling matrix element is found to be [20]
(6)
where a numerically negligible correction whose relative size is has been omitted.
Eqs. (2) and (6) may be combined to yield an expression for the tunneling width in terms of the properties of the nuclear potential barrier:
(7)
Note that this result for the net tunneling width, which includes tunneling and electromagnetic decay on an equal footing, differs by a factor of from the bare tunneling width into an infinitely broadened, fully continuum ND spectrum. From the values of , , and listed in Table 1, one sees that usage of such a bare ND-continuum result [7,27] could result in a significant difference.nucleus(I) | (meV) | (meV) | (eV) | (meV) | ℏω (MeV) | Refs. | |||
---|---|---|---|---|---|---|---|---|---|
192Hg(12) | 0.26 | 0.128 | 0.613 | 135. | 0.049 | 0.24 | 1.8 | 6.5 | [18,28] |
192Hg(10) | 0.92 | 0.050 | 0.733 | 89. | 0.37 | 0.24 | 1.3 | 4.7 | [18,28] |
192Pb(16) | <0.01 | 0.487 | 0.192 | 1,362. | <0.0050 | 0.24 | >1.8 | >6.5 | [13,17] |
192Pb(14) | 0.02 | 0.266 | 0.201 | 1,258. | 0.0056 | 0.24 | 1.7 | 6.2 | [13,17] |
192Pb(12) | 0.34 | 0.132 | 0.200 | 1,272. | 0.10 | 0.24 | 1.3 | 4.6 | [13,17] |
192Pb(10) | 0.88 | 0.048 | 0.188 | 1,410. | 1.9 | 0.24 | 0.76 | 2.7 | [13,17] |
192Pb(8) | >0.75 | 0.016 | 0.169 | 1,681. | >0.067 | 0.24 | <1.2 | <4.3 | [13,17] |
194Hg(12) | 0.42 | 0.097 | 4.8 | 16.3 | 0.071 | 0.24 | 2.3 | 8.4 | [29–32] |
194Hg(10) | >0.91 | 0.039 | 4.1 | 26.2 | >0.44 | 0.24 | <2.0 | <7.1 | [29–32] |
194Hg(12) | 0.40 | 0.108 | 21. | 344. | 0.072 | 0.24 | 2.1 | 7.6 | [33] |
194Hg(10) | 0.97 | 0.046 | 20. | 493. | 1.6 | 0.24 | 1.6 | 5.8 | [33] |
194Hg(12) | 0.40 | 0.086 | 1.345 | 19. | 0.060 | 0.24 | 2.2 | 7.8 | [18,31] |
194Hg(10) | ⩾0.95 | 0.033 | 1.487 | 14. | ⩾1.1 | 0.24 | ⩽1.8 | ⩽6.5 | [18,31] |
194Hg(15) | 0.10 | 0.230 | 4. | 26.5 | 0.026 | 0.24 | 2.4 | 8.6 | [31,32] |
194Hg(13) | 0.16 | 0.110 | 4.5 | 19.9 | 0.021 | 0.24 | 2.5 | 8.9 | [31,32] |
194Hg(11) | >0.93 | 0.048 | 6.4 | 7.2 | >0.71 | 0.24 | <2.2 | <7.8 | [31,32] |
194Pb(10) | 0.10 | 0.045 | 0.08 | 21,700. | 0.0053 | 0.24 | 1.1 | 4.1 | [32,34–36] |
194Pb(8) | 0.38 | 0.014 | 0.50 | 2,200. | 0.0087 | 0.24 | 1.6 | 5.8 | [32,34–36] |
194Pb(6) | >0.91 | 0.003 | 0.65 | 1,400. | >0.032 | 0.24 | <1.5 | <5.5 | [32,34–36] |
194Pb(12) | <0.01 | 0.125 | 0.476 | 236. | <0.0013 | 0.24 | >2.2 | >8.0 | [17,32] |
194Pb(10) | 0.10 | 0.045 | 0.470 | 244. | 0.0051 | 0.24 | 2.0 | 7.2 | [17,32] |
194Pb(8) | 0.35 | 0.014 | 0.445 | 273. | 0.0077 | 0.24 | 1.9 | 6.9 | [17,32] |
194Pb(6) | >0.96 | 0.003 | 0.405 | 333. | >0.088 | 0.24 | <1.5 | <5.4 | [17,32] |
152Dy(28) | 0.40 | 10. | 17. | 220. | 11. | 0.56 | 1.6 | 5.8 | [33] |
152Dy(26) | 0.81 | 7. | 17. | 194. | 140. | 0.56 | 1.2 | 4.4 | [33] |
The energy barriers (in units of the barrier frequency) obtained by solving Eq. (7) for all SD decays for which the four parameters, , , , and , are known are listed in Table 1 (see also Fig. 1). Also listed is the tunneling action , which has been chosen as a reasonable value for this characteristic measure of barrier strength [7,27]. Note that depends only weakly (logarithmically) on the barrier frequency ω. In the literature, the value has been used [3,27], but we shall determine ω self-consistently below.
For almost all decay-out sequences, we find that the barrier height decreases with decreasing angular momentum. The sole exception is the odd-spin 194Hg sequence, for which the two highest-spin calculated barriers are so close that statistical fluctuations about the mean-square matrix elements of Eqs. (2) and (6) are sufficient to reverse the trend. This could occur, for example, due to an accidental near-degeneracy of the SD and ND states in 194Hg(15), which would lead to a larger than expected branching ratio .
3. Centrifugal tunnel barrier
Finally, we address whether the decrease in the tunnel barrier with decreasing spin (cf. Fig. 1 and Table 1) is consistent with the centrifugal barrier of a spinning nucleus. If the superdeformed nucleus and the saddle configuration at the top of the energy barrier are described as rigid rotors with moments of inertia and , respectively, then the barrier height as a function of the angular momentum quantum number I is simply the sum of the barrier and its rotational increase, minus the rotational increase of the bottom of the SD well (from which W is measured):
(8)
Although the rigid-rotor model is a simplification, nevertheless the decrease in the barrier height between successive SD states can be rigorously expressed in terms of the kinetic moments of inertia:
(9)
The kinetic moments of inertia of several SD yrast states in the 150 and 190 mass regions have been measured. For 152Dy, and the aspect ratio [33], with a and b the smaller and larger radii of the nucleus, respectively. For 192Hg, and the aspect ratio [37].
The moment of inertia of the barrier configuration is not measured, but must be determined theoretically. This could be done by applying the Strutinsky shell correction method to the cranking model [38]. However, we employ a phenomenological two-fluid model [39] in which only the region outside the largest possible central sphere contributes to the moment of inertia. Within this two-fluid model, we find that the moment of inertia is [40]
(10)
where the nucleus has been taken as a prolate ellipsoid of revolution with aspect ratio η and atomic mass number A, is the mass of a nucleon, and . With these parameters, the measured kinetic moments of inertia of 152Dy and 192Hg at decay out are reproduced to within 1% accuracy. To leading order in the quadrupole deformation parameter ε (see Ref. [39]) Eq. (10), gives
(11)
For 152Dy, the barrier occurs at an aspect ratio of [41], so that and . Assuming a constant barrier frequency, and comparing to the results from Table 1 (), implies a barrier frequency .
For 192Hg, the barrier is estimated to occur at an aspect ratio of [37], so that and . Assuming a constant barrier frequency, and comparing to the results from Table 1 (), implies a barrier frequency . Because the logarithmic dependence of W on ω almost completely cancels out in such a calculation, the differences and ΔS are nearly independent of the choice of that parameter.
4. Conclusions
In conclusion, we have determined the barrier height W and tunneling action S for decay-out of a superdeformed band by combining an asymptotically exact quantum tunneling calculation with a two-state dynamical model. The Table presents our numerical results for all superdeformed decays for which sufficient experimental data are known. We find that the tunnel barrier decreases significantly with decreasing spin during the decay-out process. The spin-dependence of the barrier is explained quantitatively in terms of the variation of the centrifugal barrier within a two-fluid model of nuclear rotation, which in turn allows us to self-consistently predict the tunnel barrier's curvature.
The approach we suggest here makes use of a systematic, controlled expansion for the tunneling dynamics of the decaying nucleus. The advantage of such an approach is that the expansion has a well defined accuracy, which can be improved by keeping higher-order sub-dominant terms in Eq. (2). As such, the model can be used to judge the appropriateness of approximations used in other models of SD decay-out. For example, it is to be noted that, while some of the physical results presented here differ strongly from those of the bare ND continuum model [7,27], similarities can be found in the results of the statistical mixing model of Ref. [9].
The results presented in this Letter complete the chain of reasoning needed to connect the intriguing phenomenology of the decay-out process in superdeformed nuclei with an understanding of the underlying nuclear structure. Our results indicate that the rapidity and universality of the decay-out profiles can be explained straightforwardly within our two-state dynamical model by the decrease of the centrifugal barrier between the super-deformed and normal-deformed energy wells with decreasing spin.
Acknowledgements
The authors thank Teng Lek Khoo and Anna Wilson for useful discussions, and TRIUMF and the Institute for Nuclear Theory at the University of Washington for their hospitality and the Department of Energy for partial support during the completion of this work. B.R.B. was supported in part by NSF grant PHY-0555396. D.L.S. was supported in part by NSF grant PHY-06501077.
References
- [1]
- P.J. Twin, B.M. Nyakó, A.H. Nelson, J. Simpson, M.A. Bentley, H.W. Cranmer-Gordon, P.D. Forsyth, D.Howe, A.R. Mokhtar, J.D. Morrison, J.F. Sharpey-Schafer, G. SlettenPhys. Rev. Lett., 57 (1986), p. 811
- [2]
- E. Vigezzi, R.A. Broglia, T. DøssingNucl. Phys. A, 520 (1990), p. 179c
- [3]
- E. Vigezzi, R.A. Broglia, T. DøssingPhys. Lett. B, 249 (1990), p. 163
- [4]
- P.J. TwinNucl. Phys. A, 520 (1990), p. 17c
- [5]
- W. Satuła, S. Ćwiok, W. Nazarewicz, R. Wyss, A. JohnsonNucl. Phys. A, 529 (1991), p. 289
- [6]
- Y.R. Shimizu, F. Barranco, R.A. Broglia, T. Døssing, E. VigezziPhys. Lett. B, 274 (1992), p. 253
- [7]
- Y.R. Shimizu, E. Vigezzi, T. Døssing, R.A. BrogliaNucl. Phys. A, 557 (1993), p. 99c
- [8]
- R. Krücken, A. Dewald, P. von Brentano, D. Brazzacco, C. Rossi-AlvarezPhys. Rev. C, 54 (1996), p. 1182
- [9]
- J. Gu, H.A. WeidenmüllerNucl. Phys. A, 660 (1999), p. 197
- [10]
- S. ÅbergPhys. Rev. Lett., 82 (1999), p. 299
- [11]
- C.A. Stafford, B.R. BarrettPhys. Rev. C, 60 (1999), p. 051305
- [12]
- Y.R. Shimizu, M. Matsuo, K. YoshidaNucl. Phys. A, 682 (2001), p. 464c
- [13]
- A.N. Wilson, G.D. Dracoulis, A.P. Byrne, P.M. Davidson, G.J. Lane, R.M. Clark, P. Fallon, A. Görgen, A.O. Macchiavelli, D. WardPhys. Rev. Lett., 90 (2003), p. 142501
- [14]
- A.Y. Dzyublik, V.V. UtyuzhPhys. Rev. C, 68 (2003), p. 024311
- [15]
- D.M. Cardamone, C.A. Stafford, B.R. BarrettPhys. Rev. Lett., 91 (2003), p. 102502
- [16]
- C. Andreoiu, T. Dssing, C. Fahlander, I. Ragnarsson, D. Rudolph, S. Berg, R.A. Austin, M.P. Carpenter, R.M. Clark, R.V. Janssens, T.L. Khoo, F.G. Kondev, T. Lauritsen, T. Rodinger, D.G. Sarantites, D.Seweryniak, T. Steinhardt, C.E. Svensson, O. Thelen, J.C. WaddingtonPhys. Rev. Lett., 91 (2003), p. 232502
- [17]
- A.N. Wilson, P.M. DavidsonPhys. Rev. C, 69 (2004), p. 041303(R)
- [18]
- A.N. Wilson, A.J. Sargeant, P.M. Davidson, M.S. HusseinPhys. Rev. C, 71 (2005), p. 034319
- [19]
- E. Caurier, F. Nowacki, A. PovesPhys. Rev. Lett., 95 (2005), p. 042502
- [20]
- D.M. Cardamone, B.R. Barrett, C.A. StaffordPhys. Lett. B, 661 (2008), p. 233
- [21]
- A.N. Wilson, S.S. Szigeti, P.M. Davidson, J.I. Rogers, D.M. CardamonePhys. Rev. C, 79 (2009), p. 014312
- [22]
- P. Hänggi, P. Talkner, M. BorkovecRev. Mod. Phys., 62 (1990), p. 251
- [23]
- A. Caldeira, A. LeggettAnn. Phys. (NY), 149 (1983), p. 374
- [24]
- E. Freidkin, P. Riseborough, P. HänggiZ. Phys. B, 64 (1986), p. 237
- [25]
- E. Freidkin, P. Riseborough, P. HänggiPhys. Rev. B, 34 (1986), p. 1952
- [26]
- D.L. SteinBraz. J. Phys., 35 (2005), p. 242
- [27]
- T.L. Khoo, T. Lauritsen, I. Ahmad, M.P. Carpenter, P.B. Fernandez, R.V.F. Janssens, E.F. Moore, F.L.H.Wolfs, P. Benet, P.J. Daly, K.B. Beard, U. Garg, D. Ye, M.W. DrigertNucl. Phys. A, 557 (1993), p. 83c
- [28]
- T. Lauritsen, T.L. Khoo, I. Ahmad, M.P. Carpenter, R.V.F. Janssens, A. Korichi, A. Lopez-Martens, H.Amro, S. Berger, L. Calderin, T. Døssing, S.M. Fischer, G. Hackman, F. Hannachi, C.J. Lister, E.F.Moore, D.T. Nisius, C. Schück, S. SiemPhys. Rev. C, 62 (2000), p. 044316
- [29]
- T.L. Khoo, M.P. Carpenter, T. Lauritsen, D. Ackermann, I. Ahmad, D.J. Blumenthal, S.M. Fischer, R.V.F.Janssens, D. Nisius, E.F. Moore, A. Lopez-Martens, T. Døssing, R. Kruecken, S.J. Asztalos, J.A.Becker, L. Bernstein, R.M. Clark, M.A. Deleplanque, R.M. Diamons, P. Fallon, L.P. Farris, F. Hannachi, E.A. Henry, A. Korichi, I.Y. Lee, A.O. Macchiavelli, F.S. StephensPhys. Rev. Lett., 76 (1996), p. 1583
- [30]
- R. Kühn, A. Dewald, R. Krücken, C. Meier, R. Peusquens, H. Tiessler, O. Vogel, S. Kasemann, P. von Brentano, D. Bazzacco, C. Rossi-Alvarez, S. Lunardi, J. de BoerPhys. Rev. C, 55 (1997), p. R1002
- [31]
- E.F. Moore, T. Lauritsen, R.V.F. Janssens, T.L. Khoo, D. Ackermann, I. Ahmad, H. Amro, D. Blumenthal, M.P. Carpenter, S.M. Fischer, G. Hackman, D. Nisius, F. Hannachi, A. Lopez-Martens, A. Korichi, S.Asztalos, R.M. Clark, M.A. Deleplanque, R.M. Diamond, P. Fallon, I.Y. Lee, A.O. Macchiavelli, F.S.Stephens, J.A. Becker, L. Bernstein, L.P. Farris, E.A. HenryPhys. Rev. C, 55 (1997), p. R2150
- [32]
- R. Krücken, A. Dewald, P. von Brentano, H.A. WeidenmüllerPhys. Rev. C, 64 (2001), p. 064316
- [33]
- T. Lauritsen, M.P. Carpenter, T. Døssing, P. Fallon, B. Herskind, R.V.F. Janssens, D.G. Jenkins, T.L.Khoo, F.G. Kondev, A. Lopez-Martens, A.O. Macchiavelli, D. Ward, K.S.A. Saleem, I. Ahmad, R. Clark, M. Cromaz, J.P. Greene, F. Hannachi, A.M. Heinz, A. Korichi, G. Lane, C.J. Lister, P. Reiter, D.Seweryniak, S. Siem, R.C. Vondrasek, I. WiedenhöverPhys. Rev. Lett., 88 (2002), p. 042501
- [34]
- P. Willsau, H. Hubel, W. Korten, F. Azaiez, M.A. Deleplanque, R.M. Diamond, A.O. Macchiavelli, F.S.Stephens, H. Kluge, F. Hannachi, J.C. Bacelar, J.A. Becker, M.J. Brinkman, E.A. Henry, A. Kuhnert, T.F. Wang, J.A. Draper, E. RubelZ. Phys. A, 344 (1993), p. 351
- [35]
- A. Lopez-Martens, F. Hannachi, A. Korichi, C. Schück, E. Geuorguieva, C. Vieu, B. Haas, R. Lucas, A.Astier, G. Baldsiefen, M. Carpenter, G. de France, R. Duffait, L. Ducroux, Y.L. Coz, C. Finck, A. Gorgen, H. Hübel, T.L. Khoo, T. Lauritsen, M. Meyer, D. Prévost, N. Redon, C. Rigollet, H. Savajols, J.F.Sharpey-Schafer, O. Stezowski, C. Theisen, U.V. Severen, J.P. Vivien, A.N. WilsonPhys. Lett. B, 380 (1996), p. 18
- [36]
- K. Hauschild, L.A. Bernstein, J.A. Becker, D.E. Archer, R.W. Bauer, D.P. McNabb, J.A. Cizewski, K.-Y.Ding, W. Younes, R. Krücken, R.M. Diamond, R.M. Clark, P. Fallon, I.-Y. Lee, A.O. Macchiavelli, R.MacLeod, G.J. Schmid, M.A. Deleplanque, F.S. Stephens, W.H. KellyPhys. Rev. C, 55 (1997), p. 2819
- [37]
- R.V.F. Janssens, T.L. KhooAnnu. Rev. Nucl. Part. Sci., 41 (1991), p. 321
- [38]
- V.M. StrutinskyNucl. Phys. A, 122 (1968), p. 1
- [39]
- S.G. Nilsson, I. RagnarssonShapes and Shells in Nuclear StructureCambridge University Press, Cambridge (1995)
- [40]
- B.R. Barrett, C.A. Stafford, unpublished (2009)
- [41]
- P.J. Nolan, P.J. TwinAnnu. Rev. Nucl. Part. Sci., 38 (1988), p. 533
- 1
- Present address: Department of Physics and Astronomy, California State University, Sacramento, CA 95819, United States.
- 2
- Present address: Department of Physics and Astronomy, University of California, Irvine, CA 92697, United States.
Copyright © 2010 Elsevier B.V.
For further detauils log on website :https://www.sciencedirect.com/science/article/pii/S0370269310003734
No comments:
Post a Comment