Recent developments in heavyion fusion reactions around the Coulomb barrier
Résumé
The nuclear fusion is a reaction to form a compound nucleus. It plays an important role in several circumstances in nuclear physics as well as in nuclear astrophysics, such as synthesis of superheavy elements and nucleosynthesis in stars. Here we discuss two recent theoretical developments in heavyion fusion reactions at energies around the Coulomb barrier. The first topic is a generalization of the Wong formula for fusion cross sections in a singlechannel problem. By introducing an energy dependence to the barrier parameters, we show that the generalized formula leads to results practically indistinguishable from a full quantal calculation, even for light symmetric systems such as C+C, for which fusion cross sections show an oscillatory behavior. We then discuss a semimicroscopic modeling of heavyion fusion reactions, which combine the coupledchannels approach to the stateoftheart nuclear structure calculations for lowlying collective motions. We apply this method to subbarrier fusion reactions of Ni+Ni and Ca+Ni systems, and discuss the role of anharmonicity of the lowlying vibrational motions.
Recent developments in heavyion fusion reactions around the Coulomb barrier
K. Hagino, N. Rowley, et J.M. Yao
1Department of Physics, Tohoku University, Sendai 9808578, Japan
\@textsuperscript2Research Center for Electron Photon Science, Tohoku University, 121 Mikamine, Sendai 9820826, Japan
\@textsuperscript3National Astronomical Observatory of Japan, 2211 Osawa,
Mitaka, Tokyo 1818588, Japan
\@textsuperscript4Institut de Physique Nucléaire, UMR 8608, CNRSIN2P3 et Université
de Paris Sud, 91406 Orsay Cedex, France
\@textsuperscript5School of Physical Science and Technology, Southwest University, Chongqing 400715, China
\@textsuperscript6Department of Physics and Astronomy, University of North Carolina,
Chapel Hill, North Carolina 275163255, USA

Résumé.
1 Introduction
Fusion is defined as a reaction in which two separate nuclei combine together to form a (hot) compound nucleus. It is indispensable to understand its reaction dynamics in order to understand synthesis of superheavy elements and nucleosynthesis in stars. Moreover, fusion reactions inherently contain rich physics, as there is a strong interplay in it, especially at energies around the Coulomb barrier, between the nuclear reaction and the nuclear structure. That is, it has been known well that fusion cross sections at subbarrier energies are largely enhanced relative to a prediction of a onedimensional potential model, due to the couplings of the relative motion between the colliding nuclei to several nuclear intrinsic motions [1, 2, 3, 4, 5]. This is in marked contrast to high energy nuclear reactions, in which the reaction dynamics is much simpler and the distorted wave Born approximation (DWBA) often suffices its treatment. Fusion also offers a unique opportunity to study quantum tunneling with many degrees of freedom. That is, fusion takes place only by quantum tunneling at energies below the Coulomb barrier, and the subbarrier enhancement of fusion cross sections can be viewed as a consequence of a coupling assisted tunneling. Heavyion fusion reactions are unique in this respect because a variety of intrinsic degrees of freedom are involved, such as a surface vibration with several multipolarities, various sort of nuclear deformations and the associated rotational motion, and several types of particle transfer processes. This is in contrast to atomic and molecular collisions, in which only a limited types of intrinsic motion are involved. Also, in heavyion fusion reactions, the incident energy can be varied in order to study the energy dependence of the tunneling probability, whereas the energy is fixed in many other tunneling phenomena in nuclear physics, such as alpha decays.
In this contribution, we shall discuss two recent theoretical developments in heavyion fusion reactions. The first topic is a generalization of the Wong formula for fusion cross sections [6]. Although this celebrated formula yields almost exact results for singlechannel calculations for relatively heavy systems such as O+Sm, it tends to overestimate the cross section for light systems such as C+C. We generalize the formula to take account of the energy dependence of the barrier parameters and show that the energydependent version reproduces almost perfectly results of a full quantal calculation. We also discuss the deviations arising from the discrete nature of the intervening angular momenta, whose effect can lead to an oscillatory contribution to the excitation function.
The second topic is a semimicroscopic modeling of heavyion fusion reactions [7]. For this purpose, we first describe microscopically lowlying collective excitations of atomic nuclei with the multireference covariant density functional theory, and combine them with coupledchannels calculations. We use the calculated transition strengths among several collective states as inputs to the coupledchannels calculations. This approach provides a natural way to describe anharmonic multiphonon excitations as well as a deviation of rotational excitations from a simple rigid rotor. We apply this method to subbarrier fusion reactions of Ni+Ni and Ca+Ni systems, and discuss the role of anharmonicity in subbarrier fusion of these systems.
2 Generalization of the Wong formula for fusion cross sections
Let us first discuss the wellknown Wong formula for fusion cross sections. The simplest approach to heavyion fusion reactions is to use the potential model, assuming that both the projectile and the target nuclei are inert. Fusion cross sections are then obtained by calculating the matrix, , for each partial wave. If all the flux crossing the Coulomb barrier fuses, the fusion cross sections read
(1) 
where is the wave number associated with the energy ( being the reduced mass) and is the penetrability of the Coulomb barrier for the th partial wave.
Wong has derived a compact approximate expression for Eq. (1) [8], which has by now been known as the Wong formula. To this end, Wong introduced the following three approximations:

the parabolic approximation. The total (that is, the nuclear + the Coulomb) internucleus potential for the wave is approximated by a parabolic function, that is,
(2) where , , and are the barrier height, the barrier position, and the “barrier curvature”, respectively.

the independent barrier position and curvature. The barrier position, , and the barrier curvature, , are assumed to be independent of . In this case, the effective potential for the th partial wave reads,
(3) The penetrability can then be calculated analytically with the HillWheeler formula [9] as,
(4) 
the continuum approximation for . The angular momentum is treated as a continuous variable and the sum in Eq. (1) is replaced by the integral,
(5)
With these approximations, the fusion cross sections are obtained as [8, 1, 3],
(6) 
Notice that this formula yields the classical fusion cross sections,
(7) 
at energies well above the Coulomb barrier, .
The left panel of Fig. 1 shows the performance of the Wong formula for the O+Sm system. The solid line is obtained quantum mechanically by solving the Schrödinger equation for a given internucleus potential, while the dashed line is obtained with the Wong formula, Eq. (6). Even though the Wong formula slightly overestimates the fusion cross sections at energies well above the Coulomb barrier, reflecting the fact that the barrier position gradually decreases as a function of , the overall agreement is satisfactory. In contrast, the situation drastically changes for light systems. The right panel of Fig. 1 shows a comparison between the quantal and the approximate fusion cross sections for the C+C system. Here, in order to simplify the discussion, for the moment we ignore the Bose symmetry of this identical spin0 system and sum over all even and odd partial waves. One can see that the Wong formula largely overestimates fusion cross sections, and the agreement is much worse than the heavier system, O+Sm. This is because the centrifugal potential plays a more important role in light systems, and the Coulomb barrier is not rigid against a variation of angular momentum [6].
One can improve the Wong formula by introducing the energy dependence to the barrier parameters . To this end, we first introduce the grazing angular momentum , at which the barrier height of the effective potential becomes identical to the incident energy, . If we denote the barrier position for as , and are related to each other as,
(8) 
where is the sum of the Coulomb and the nuclear potentials at . An idea of the generalized Wong formula is to replace the barrier parameters with the energy dependent ones, , where is the barrier curvature for the grazing angular momentum . That is, the generalized version of the Wong formula now reads [6],
(9) 
Notice that we still assume an independent barrier position and curvature in integrating Eq. (5). However, the generalization is such that these parameters are evaluated at the grazing angular momentum for each energy , rather than using the values for the wave. The filled circles in Fig. 1 show the fusion cross sections obtained with the generalized Wong formula, Eq. (9). One can see that the quantal results are well reproduced by introducing the energy dependence to the barrier parameters, both for the O+Sm and the C+C systems.
One can also go beyond the continuum approximation for the angular momentum sum. Using the Poisson sum formula [10], the next order correction to the Wong formula reads [6, 11],
(10) 
where is the smooth part of fusion cross sections given by the Wong formula, Eq. (9), while the oscillatory contribution is given by [6, 11],
(11) 
with
(12) 
Usually, the oscillatory contribution, , is much smaller than the smooth part given by . There are, however, certain circumstances where the oscillatory contribution becomes significantly large. This is for a light symmetric system of two identical spin0 bosons, in which all odd partial waves disappear due to the symmetrization of the wave function. In this case, the fusion cross section is given by twice the sum over the even partial waves. The resultant formula now becomes [6, 11],
(13) 
whereas the smooth part is still given by the Wong formula, Eq. (9). The oscillatory part is now significantly larger than before, since the negative exponent has been reduced by a factor of 2, that is .
Figure 2 shows fusion cross sections for the C+C system, obtained by taking only the even partial waves. The meaning of each line is the same as in Fig. 1. One can see that the fusion cross sections significantly oscillate as a function of . The generalized Wong formula, supplemented with the oscillatory correction, gives an excellent approximation to the exact results.
3 Semimicroscopic modeling of heavyion fusion reactions with a beyondmean field method
The potential model presented in the previous section works for light systems, such as N+C, but it largely underestimates fusion cross sections for heavier systems, such as O+Sm, at energies below the Coulomb barrier [1]. It has been well recognized by now that this large enhancement of subbarrier fusion cross sections is caused by the couplings of the relative motion between the colliding nuclei to several nuclear intrinsic degrees of freedom, such as lowlying collective excitations in the colliding nuclei as well as several nucleon transfer processes [1, 2, 3, 4]. A natural framework for heavyion subbarrier fusion has thus been the coupledchannels method with relevant degrees of freedom [1, 12]. This approach has not only successfully accounted for the subbarrier enhancement of fusion cross sections for many systems but has also provided a natural interpretation of the so called fusion barrier distributions [1, 2, 13, 14].
In the coupledchannels approach, if the projectile nucleus is assumed to be inert, one expands the total wave function of the system in terms of the eigenfunctions of the collective states in the target nucleus, , as
(14) 
where and are the angular momenta for the target state and for the relative motion, respectively. Here, we have introduced the isocentrifugal approximation [1], and have assumed that and (thus also ) do not change by the excitation of the target nucleus. Substituting Eq. (14) to the projected Schrödinger equation, , the coupledchannels equations for the radial wave functions read [1],
(15) 
where is the energy of the target state . are the coupling matrix elements given by
(16) 
where is the coupling potential and is the excitation operator with a multipolarity .
In heavyion fusion reactions at energies around the Coulomb barrier, multiple excitations to higher members of collective modes, such as multiphonon states and highspin states in the ground state rotational band, often play an important role [2, 14, 16]. As an example, Fig. 3 shows fusion cross sections for the Ni+Ni system obtained with the coupledchannels calculations which take into account the vibrational coupling to quadrupole phonon states. One can see that the coupling to the first excited state is insufficient in order to account for the subbarrier enhancement of fusion cross sections for this system. In order to take into account the multiple excitations in coupledchannels calculations, one usually uses the simple harmonic oscillator model for vibrational nuclei [1]. In this approximation, the phonon multiplet states are combined together to form a single effective channel as [1]
(17) 
where is the phonon creation operator and is the vacuum state. The energy of this channel is given by , where is the oscillator quantum, and the matrix elements of are related to the phonon numbers as . The solid line in Fig. 3 shows the result of the two phonon couplings. One can see that the subbarrier fusion enhancement is now well reproduced.
Even though the harmonic oscillator approximation appears to work well for the Ni+Ni system, most nuclei, including Ni, do not have a pure harmonic oscillator spectrum in reality. Concerning the Ni nucleus, the degeneracy of the experimentally known twophonon triplet is considerably broken (see Fig. 4), and also a recent theoretical calculation indicates that the strengths among the collective levels in Ni deviate largely from the harmonic oscillator limit [17]. Such anharmonic vibrations can be described, e.g., a multireference densityfunctional theory (MRDFT), which has been rapidly developed for the past decade [18, 19]. This method is based on the so called beyondmeanfield approximation, which incorporates the angular momentum and particle number projections for a meanfield wave function. The quantum fluctuation of the meanfield wave function is also taken into account with the generator coordinate method (GCM).
Figure 4 shows the spectrum of the Ni nucleus [7] constructed with the MRDFT calculation with the covariant density functional with the PCPK1 interaction [22]. One can see that the excitation energies of the 0, 2, and 4 states are about twice the energy of the state. On the other hand, it is interesting to notice that the overall pattern of values is quite different from what would be expected for a harmonic vibrator, in which the value from any of the twophonon triplet states to the 2 state is exactly twice the value from the 2 state to the ground state. In particular, the transition from the to the states is much smaller than that for the and the transitions. Instead, the state has a strong transition from the state, which clearly indicates that the state is not a member of the twophonon triplet. Compared to the state, the transition strength for the transition is much larger and is comparable to that for the and the transitions. This suggests that the state is a better candidate for a member of the twophonon triplets than the state, even though the excitation energy is a little large.
In order to take into account such anharmonic nature of the vibrational excitations of Ni in the coupledchannels calculations, we first replace the operator in Eq. (16) with the corresponding microscopic multipole operator. That is, we replace the matrix elements of as,
(18) 
where is the microscopic multipole operator, and and are the radius and the charge number of the target nucleus, respectively. on the right hand side is a manybody wave function evaluated with the MRCDFT method. We still use a phenomenological potential for and , and we call our method a semimicroscopic approach.
In heavyion fusion reactions, the coupling of the ground state to the lowestlying collective state is most important. The strength for such coupling can often be estimated from an experimental transition probability. Since it is too much to expect that a MRCDFT calculation agrees perfectly with experimental data, we introduce an overall scaling factor to the matrix elements, Eq. (18), so that the transition from the lowestlying collective state to the ground state is consistent with experimental data. The MRCDFT calculation then provides the relative strengths among collective levels, which are often not available experimentally. The excitation energy, on the other hand, is often known well for many levels, and we simply take them in the calculations whenever they are available experimentally.
Figures 5 shows the fusion cross section and the fusion barrier distribution [2, 13] for the Ni+Ni reaction so obtained. The dashed line shows the result of the coupledchannels calculations including up to the double phonon states in the harmonic oscillator limit. All the mutual excitations between the projectile and the target nuclei are included. On the other hand, the solid line in the figure is obtained with the coupling strengths evaluated with the MRCDFT method. To this end, we include the 0, 2, 0, 2, and 4 states in Ni in the coupledchannels calculations. Again, all the mutual excitation channels are taken into account. For comparison, the figure also shows the result of nocoupling limit by the dotted line. One can see that the calculations in the harmonic limit overestimate fusion cross sections at the two lowest energies, while the MRCDFT calculations underpredict fusion cross sections around 95 MeV. The energy dependence of fusion cross sections can be investigated with the fusion barrier distribution, which is plotted in the right panel of the figure. Even though the effect of anharmonicity is not large, one observes that the MRCDFT calculation leads to a minor improvement by considerably smearing each peak.
We carry out similar calculations also for the Ca + Ni system [7] (see Fig. 6). The effect of channel coupling is much smaller than in the Ni + Ni system due to a smaller charge product, . Nevertheless, we again observe that the anharmonicity effect in Ni smears the fusion barrier distribution, leading to a somewhat better agreement with the experimental fusion barrier distribution as compared to the results in the harmonic oscillator limit.
4 Summary
The research field of heavyion subbarrier fusion reactions started in the late ’70s, when a large enhancement of fusion cross sections was experimentally discovered with respect to the prediction of a simple potential model. The Wong formula has been widely used to estimate fusion cross sections for a given singlechannel potential, providing reference cross sections in the absence of channel couplings to discuss the subbarrier enhancement of cross sections. In this contribution, we have first extended the Wong formula by including the energy dependence of the parameters entering the formula, that is, the barrier height, the barrier position, and the barrier curvature. Evaluating these parameters for the grazing angular momentum at each energy, rather than at , we have shown that the energydependent version of Wong’s formula reproduces the exact result well, even for light systems. We have also derived a compact analytic formula for the oscillatory part of fusion cross sections, which originate from the discrete nature of angular momentum. These oscillatory parts become important in light symmetric systems, for which the symmetrization of the system amplifies the oscillations.
We have then proposed a semimicroscopic approach to heavyion subbarrier fusion reactions. The basic idea of this approach is to combine the stateoftheart nuclear structure calculations with coupledchannels calculations. For this purpose, we have used a multireference density functional theory (MRDFT) based on the beyondmeanfield approach. The MRDFT provides transition strengths among collective states without resorting to the harmonic oscillator model. We have applied this approach to the Ni+Ni and Ca+Ni fusion reactions, and have found that the anharmonicities smear the fusion barrier distributions, somewhat improving the agreement with the experimental data. It would be an interesting future problem to extend this treatment to heavyion elastic and inelastic scattering, for which the doublefolding approach is applicable and thus a fully microscopic approach can be developed using the multireference density functional theory. That approach would also be useful in applications to nuclear data as well as the problem of nuclear transmutation [24].
Acknowledgements
This work was partially supported by the National Natural Science Foundation of China under Grant Nos. 11305134, 11105111, and the Fundamental Research Funds for the Central University (XDJK2013C028).
Références

[1]
K. Hagino and N. Takigawa, Prog. Theor. Phys.
128, 1001 (2012);
http://www.nucl.phys.tohoku.ac.jp/ hagino/fusioncorr.pdf  [2] M. Dasgupta, D.J. Hinde, N. Rowley and A.M. Stefanini, Annu. Rev. Nucl. Part. Sci. 48, 401 (1998).
 [3] A.B. Balantekin and N. Takigawa, Rev. Mod. Phys. 70, 77(1998).
 [4] B.B. Back, H. Esbensen, C.L. Jiang and K.E. Rehm, Rev. Mod. Phys. 86, 317 (2014).
 [5] L.F. Canto, P.R.S. Gomes, R. Donangelo, J. Lubian and M.S. Hussein, Phys. Rep. 596, 1 (2015).
 [6] N. Rowley and K. Hagino, Phys. Rev. C91, 044617 (2015).
 [7] K. Hagino and J.M. Yao, Phys. Rev. C91, 064606 (2015).
 [8] C.Y. Wong, Phys. Rev. Lett. 31, 766 (1973).
 [9] D.L. Hill and J.A. Wheeler, Phys. Rev. 89, 1102 (1953).
 [10] D.M. Brink, SemiClassical Methods for NucleusNucleus Scattering, (University Press, Cambridge, 1985).
 [11] N. Poffé, N. Rowley and R. Lindsay, Nucl. Phys. A410, 498 (1983).
 [12] K. Hagino, N. Rowley, and A.T. Kruppa, Comp. Phys. Comm. 123, 143 (1999).
 [13] N. Rowley, G.R. Satchler, and P.H. Stelson, Phys. Lett. B254, 25 (1991).
 [14] J.R. Leigh et al., Phys. Rev. C52, 3151 (1995).
 [15] M. Beckerman et al., Phys. Rev. C23, 1581 (1981).
 [16] A.M. Stefanini et al., Phys. Rev. Lett. 74, 864 (1995).
 [17] J.M. Yao, M. Bender, and P.H. Heenen, Phys. Rev. C91, 024301 (2015).
 [18] M. Bender, P.H. Heenen, and P.G. Reinhard, Rev. Mod. Phys. 75, 121 (2003).
 [19] J.M. Yao, K. Hagino, Z.P. Li, J. Meng, and P. Ring, Phys. Rev. C89, 054306 (2014).
 [20] National Nuclear Data Center, http://www.nndc.bnl.gov/.
 [21] J.M. Allmond et al., Phys. Rev. C90, 034309 (2014).
 [22] P.W. Zhao, Z.P. Li, J.M. Yao, and J. Meng, Phys. Rev. C82, 054319 (2010).
 [23] D. Bourgin et al., Phys. Rev. C90, 044610 (2014).
 [24] E.Sh. Sukhovitskii, Y.O. Lee, J. Chang, S. Chiba, and O. Iwamoto, Phys. Rev. C62, 044605 (2000).