Coulomb and nuclear potentials between deformed nuclei applied to the fusion process

NUCLEAR REACTIONS

Coulomb and nuclear potentials between deformed nuclei applied to the fusion process

L. C. Chamon^I; G. P. A. Nobre^I; D. Pereira^I; C. P. Silva^I; E. S. Rossi Jr.^I; L. R. Gasques^II; B. V. Carlson^III

^IDepartamento de Física Nuclear, Instituto de Física da Universidade de São Paulo, Caixa Postal 66318, 05315-970, São Paulo, SP, Brazil

^IIDepartment of Physics & The Joint Institute for Nuclear Astrophysics, University of Notre Dame, Notre Dame, Indiana, 46556

^IIIDepartamento de Física, Instituto Tecnológico de Aeronáutica, Centro Técnico Aeroespacial, São José dos Campos, SP, Brazil

ABSTRACT

We obtained theoretical fusion cross sections based on a new and accurate procedure to calculate the Coulomb and nuclear interactions between deformed nuclei. The calculations were performed within the context of the barrier penetration model with zero point motion for the ⁶⁴Ni + ⁶⁴Ni system.

I. INTRODUCTION

It is well known that fusion cross section data for heavy-ion systems show large enhancements at sub-barrier energies in comparison with theoretical predictions from the unidimensional barrier penetration (BP) model [1]. Several attempts to improve the BP model have been performed by considering permanent and/or vibrational deformations of the nuclei (e.g. [2]). Most works in this direction have used approximations for the corresponding deformed potentials. In recent papers [3, 4], we proposed a new procedure to calculate accurately the Coulomb and nuclear interactions between deformed nuclei. We demonstrated that, in certain cases, the usual approximations for the deformed potentials significantly differ from the more accurate results. In the present paper, we compare the results for the fusion cross section when calculating it using the approximate or the more accurate deformed potentials. The calculations were performed for the ⁶⁴Ni + ⁶⁴Ni system (data from [5–7]), because this system present very large enhancements in comparison with the predictions of undeformed (unidimensional) BP calculations.

II. THE COULOMB AND NUCLEAR INTERACTIONS

Hereafter, we assume the two parameter Fermi (2pF) distribution to describe the deformed density of a given nucleus i

where Q_i is the angle between and the symmetry axis of the deformation. The corresponding deformation parameter is defined from: d_li = b_li R_0i, and the densities are normalized to

where X is the number of protons Z or nucleons A for the charge and matter distributions, respectively. We limit the present study to l = 2 and 3. In the vibrational model, the b_l parameter vary according a Gaussian distribution of probability with standard deviation s_l. Considering a previous systematics [3], we have adopted the following deformation parameters for the ⁶⁴Ni: s₂ = 0.214 and s₃ = 0.262 .

The Coulomb interaction between two charge distributions is given by:

where is the position vector of the center of mass of nucleus 2 measured from that of the nucleus 1, and describes their relative motion. The direction of the symmetry axis of each deformation relative to is defined by two angles: q_i and f_i.

The computational resolution of Eq. (4) is very time-consuming due to the six-dimensional integral involved in the calculation. This becomes a problem in studies where it should be solved several times. Therefore, faster methods to obtain good approximations for the deformed potential are quite convenient. According [3, 4], the Coulomb potential can be obtained with good accuracy from:

where is the non-deformed Coulomb potential, obtained with the corresponding non-deformed densities. The (extensive) expressions for V_Cor(R) can be found in Refs. [3, 4].

Within the non local model [8–10], the nuclear interaction V_N is connected with the folding potential V_F through

where c is the speed of light and v is the local relative velocity between the two nuclei,

The folding potential depends on the densities of the nuclei involved in the collision:

Usually, in Eq. (8) r_i are the nucleon densities of the nuclei, and v_NN is the effective nucleon-nucleon interaction. In many works, the Paris and Reid versions of the M3Y interaction [11] have been assumed for v_NN. In [10], we demonstrated that the folding type interaction, Eq. (9), produces very similar results in comparison with those from the M3Y:

where V₀ = –456 MeV fm³ and r_m is the matter density of the nucleon. Based on the intrinsic charge distribution of the proton in free space, which has been determined by electron scattering experiments, an exponential shape has been assumed for the matter density of the nucleon [10]. Due to the delta function, the folding in Eq. (9) is named as the zero-range approach. We have also defined [10] the matter density of the nucleus (r_M) by folding the corresponding nucleon density (r_N) with the matter density of the nucleon

Thus, we distinguish the matter density of the nucleus from the nucleon distribution of the nucleus by taking into account the finite size of the nucleon. By inserting Eqs. (9) and (10) in Eq.(8), the folding potential can be recast in the following form

Due to the six-dimensional integral, the numerical resolution of Eq. (8) implies a similar problem as that for the Coulomb potential. Eq. (11) is much faster to calculate because the zero-range approach reduces it to a three-dimensional integral.

In many works (e.g. [2]), the Coulomb and nuclear deformed potentials have been calculated through the following approximate expressions:

where (R) is the folding potential obtained with the undeformed densities and a_P is the diffuseness of the potential. Taking into account the approximate exponential shape of the potential in the surface region, the diffuseness involved in Eq. (13) can be estimated by

However, we have demonstrated that such approximations are not much precise in certain cases [3]. In the next section, we show that the use of these approximations result in fusion cross sections quite different from the values obtained with the more precise deformed potentials.

III. BARRIER PENETRATION CALCULATIONS

In the context of the barrier penetration model, the effective potential is a sum of the Coulomb, nuclear and centrifugal parts:

The BP cross section is associated with the transmitted flux through

For deformed densities, the effective potential and the corresponding BP cross sections depend on the deformation parameters of the nuclei: b_li, q_i, and f_i. In this work, we assume the sudden approximation where such parameters remain constant during the collision. However, from one collision to another, these parameters can vary according their distributions of probability:

We have calculated BP cross sections by obtaining the corresponding average over the possible configurations of the deformation parameters. We have used a Monte Carlo method, where the values of b_li, q_i, and f_i are obtained randomly according their distributions of probability. Fig. 1 presents the fusion data for the ⁶⁴Ni + ⁶⁴Ni system. There are different data sets with cross sections that differ by a factor of about three in the sub-barrier energy region. These different data sets have been represented by different symbols. The dotted line in Fig. 1 represents the cross sections obtained with the undeformed potential, while the dashed and solid lines represent the results obtained with the approximate and with the more accurate expressions for the deformed potentials, respectively. Clearly, the deformation largely enhances the cross section in the sub-barrier region and almost does not affect it in the region above the barrier. The use of the approximations results in fusion cross sections quite different from the values obtained with the more precise deformed potentials.

Our calculations with the accurate potentials show that the vibrational model is in a reasonable agreement with the fusion data in a level so down as about 1 mb. Below this level, for extreme sub-barrier energies, the data present a decrease much steeper than the theoretical predictions. Similar hindrance of the data compared to predictions of coupled-channel calculations has already been observed for several systems [7, 12, 13].

IV. SUMMARY AND CONCLUSION

We have used an accurate method for calculating deformed potentials to study the fusion process for the ⁶⁴Ni + ⁶⁴Ni systems within the context of the barrier penetration model with zero point motion. We demonstrated that the accurate method provides fusion cross sections that significantly differ from those obtained with the usual approximations for the deformed potentials. The accurate method gives cross section predictions that reproduce the order of magnitude of the data to a level about 1 mb. The hindrance of the fusion data below this level, observed in extreme sub-barrier energies, is still an open problem. This phenomenon is quite important for reactions of astrophysical interest, as these usually take place in such very low energy region.

Acknowledgments

This work was partially supported by Financiadora de Estudos e Projetos (FINEP), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and The Joint Institute for Nuclear Astrophysics (JINA) NSF PHY 0216783.

[4] B. V. Carlson, L. C. Chamon, and L. R. Gasques, to appear in Phys. Rev. C.

Received on 22 December, 2004

• [1] L. C. Vaz and J. M. Alexander, Phys. Rev. C 18, 2152 (1978).
• [2] H. Esbensen, Nucl. Phys. A 352, 147 (1981).
• [3] L. C. Chamon, G. P. A. Nobre, D. Pereira, E. S. Rossi Jr., C. P. Silva, L. R. Gasques, and B. V. Carlson, Phys. Rev. C 70, 014604 (2004).
• [5] M. Beckerman, M. Salomaa, A. Sperduto, J. D. Molitoris, and A. DiRienzo, Phys. Rev. C 25, 837 (1982).
• [6] M. Stefanini, K. M. Varier, H. Zhang, F. Scarlassara, S. Beghini, G. Montagnoli, L. Miller, G. F. Segato, F. Soramel, and C. Signorini, Nucl. Phys. A 609, 91 (1996).
• [7] C. L. Jiang, K. E. Rehm, R. V. F. Janssens, H. Esbensen, I. Ahmad, B. B. Back, P. Collon, C. N. Davids, J. P. Greene, D. J. Enderson, G. Mukherjee, R. C. Pardo, M. Paul, T. O. Pennington, D. Seweryniak, S. Sinha, and Z. Zhou, Phys. Rev. Lett. 93, 012701 (2004).
• [8] M. A. Cândido Ribeiro, L. C. Chamon, D. Pereira, M. S. Hussein and D. Galetti, Phys. Rev. Lett. 78, 3270 (1997).
• [9] L. C. Chamon, D. Pereira, M. S. Hussein, M. A. Cândido Ribeiro and D. Galetti, Phys. Rev. Lett. 79, 5218 (1997).
• [10] L. C. Chamon, B. V. Carlson, L. R. Gasques, D. Pereira, C. De Conti, M. A. G. Alvarez, M. S. Hussein, M. A. Candido Ribeiro, E. S. Rossi Jr., and C. P. Silva, Phys. Rev. C 66, 014610 (2002).
• [11] G. R. Satchler and W. G. Love, Phys. Rep. 55, 183 (1979).
• [12] C. L. Jiang, H. Esbensen, K. E. Rehm, B. B. Back, R. V. F. Janssens, J. A. Caggiano, P. Collon, J. Greene, A. M. Heinz, D. J. Henderson, I. Nishinaka, T. O. Pennington, and . Seweryniak, Phys. Rev. Lett. 89, 052701 (2002).
• [13] C. L. Jiang, H. Esbensen, B. B. Back, R. V. F. Janssens, and K. E. Rehm, Phys. Rev. C 69, 014604 (2004).

Publication Dates

• Publication in this collection
  07 Nov 2005
• Date of issue
  Sept 2005

History

• Received
  22 Dec 2004