Sub-barrier fusion of Si+Si systems : does the deformation of 28 Si play a role?

Similar documents
Recent experiments in inverse kinematics with the magnetic spectrometer PRISMA

Quasi-elastic reactions : an interplay of reaction dynamics and nuclear structure

Fusion Reactions with Carbon Isotopes

(Multi-)nucleon transfer in the reactions 16 O, 3 32 S Pb

Probing Fusion Dynamics with Scattering Experiments

The heavy-ion magnetic spectrometer PRISMA

Nuclear Reactions. Shape, interaction, and excitation structures of nuclei. scattered particles. detector. solid angle. target. transmitted particles

Heavy-ion sub-barrier fusion reactions: a sensitive tool to probe nuclear structure

Response function of the magnetic spectrometer PRISMA

Annax-I. Investigation of multi-nucleon transfer reactions in

Study of multi-nucleon transfer reactions with light nuclei

Sub-barrier fusion enhancement due to neutron transfer

Testing the shell closure at N=82 via multinucleon transfer reactions at energies around the Coulomb barrier

SUB-BARRIER FUSION of HEAVY IONS

Inclusive breakup measurements of the 7 Li+ 119 Sn system.

PRE-EQUILIBRIUM α-particle EMISSION AS A PROBE TO STUDY α-clustering IN NUCLEI

Isospin influence on Fragments production in. G. Politi for NEWCHIM/ISODEC collaboration

arxiv:nucl-th/ v1 4 Nov 2003

Isospin influence on the decay modes of compound nuclei produced in the 78,86 Kr + 40,48 Ca at 10 MeV/nucleon

Accreting neutron stars provide a unique environment for nuclear reactions

Preliminary results of the indirect study of the 12 C( 12 C,α) 20 Ne reaction via the THM applied to the 16 O( 12 C,α 20 Ne )α reaction

X-ray superburst ~10 42 ergs Annual solar output ~10 41 ergs. Cumming et al., Astrophys. J. Lett. 559, L127 (2001) (2)

PRE-EQUILIBIRUM α-particle EMISSION AS A PROBE TO STUDY α-clustering IN NUCLEI

Measuring fusion excitation functions with RIBs using the stacked target technique: Problems and possible solutions

Study and comparison of the decay modes of the systems formed in the reactions

Probing surface diffuseness of nucleus-nucleus potential with quasielastic scattering at deep sub-barrier energies

Microscopic Fusion Dynamics Based on TDHF

arxiv: v1 [nucl-th] 21 Apr 2007

Extracting spectroscopic information from ( 18 O, 16 O) two-neutron transfer reactions. Manuela Cavallaro INFN -LNS (Italy)

Nuclear contribution into single-event upset in 3D on-board electronics at moderate energy cosmic proton impact

Strange Hadron Production from STAR Fixed-Target Program

Fusion of light halo nuclei

Study of multinucleon transfer (MNT) reactions of 136 Xe Pt for production of exotic nuclei

The 22 Ne(α,n) 25 Mg reaction at astrophysical energies studied via the Trojan Horse Method applied to the 2 H( 25 Mg, α 22 Ne) 1 H reaction

How to do C.C. calculations if there is only limited experimental information on intrinsic degrees of freedom?

DIFFUSENESS OF WOODS SAXON POTENTIAL AND SUB-BARRIER FUSION

Open problems in formation and decay of composite systems in heavy ion reactions

Status of the magnetic spectrometer PRISMA

Chapter 6. Summary and Conclusions

Heavy-ion fusion reactions for superheavy elements Kouichi Hagino

This is a repository copy of How well do we understand the reaction rate of C burning?.

Nonstatistical fluctuations for deep inelastic processes

Reaction dynamics for fusion of weakly-bound nuclei

Subbarrier fusion of carbon isotopes ~ from resonance structure to fusion oscillations ~

Direct reactions at low energies: Part II Interactions and couplings

Static versus energy-dependent nucleus-nucleus potential for description of sub-barrier fusion

Pre-equilibrium -particle emission as a probe to study -clustering in nuclei

Physics with Exotic Nuclei

Dipartimento di Fisica e Astronomia, Università degli Studi di Catania, Catania, Italy 3

Fusion Barrier Distribution: Concept and Experiments

Role of the Nucleus-Nucleus Potential in Sub-barrier Fusion

Thermodynamical coordinates of excited nuclear systems:

Search for 24 Mg resonances inside the Gamowwindow for 12 C+ 12 Cfusion. Neven Soić, Ruđer Bošković Institute, Zagreb, Croatia

Effect of Barrier Height on Nuclear Fusion

SOME ASPECTS OF TRANSFER REACTIONS IN LIGHT AND HEAVY ION COLLISIONS

A Comparison between Channel Selections in Heavy Ion Reactions

FIssion of Excited Systems at Tandem-Alpi

Presence of Barrier Distributions in Heavy Ion Fusion

A program for coupled-channels calculations with all order couplings for heavy-ion fusion reactions arxiv:nucl-th/ v1 30 Mar 1999

Development of Secondary Electron Time Detector for Ion Beams

Production of superheavy elements. Seminar: Key experiments in particle physics Supervisor: Kai Schweda Thorsten Heußer

Capture barrier distributions and superheavy elements

Accreting Neutron Stars

Fusion08 : Summary talk. L.Corradi

arxiv:nucl-th/ v1 15 Feb 2007

PHL424: Nuclear fusion

High-spin studies and nuclear structure in three semi-magic regions of the nuclide chart High-seniority states in Sn isotopes

Microscopic (TDHF and DC-TDHF) study of heavy-ion fusion and capture reactions with neutron-rich nuclei

Use of fragmentation beams at LNS with CHIMERA detector

The rotational γ -continuum in the mass region A 110

PROTON-PROTON FEMTOSCOPY AND ACCESS TO DYNAMICAL SOURCES AT INTERMEDIATE ENERGIES

Reaction dynamics and nuclear structure of moderately neutron-rich Ne isotopes by heavy-ion reactions

N/Z influence on the level density parameter

New application of the quasi-free reaction mechanism to study neutron induced reactions at low energy

Transfer reaction studies in inverse kinematics with the magnetic spectrometer

Multinucleon transfer processes in 64 Ni 238 U

Fission fragment mass distributions via prompt γ -ray spectroscopy

A new analysis procedure to extract fusion excitation function with large beam energy dispersions: application to the 6 Li+ 120 Sn and 7 Li+ 119 Sn

Neutrino Physics and Nuclear Astrophysics: the LUNA-MV project at Gran Sasso

NN th INTERNATIONAL CONFERENCE ON NUCLEUS-NUCLEUS COLLISSIONS June 21-26, 2015, Catania, Italy

Reaction dynamics and exotic systems: a focus on fast processes

At low excitation, the compound nucleus de-excites by statistical emission of light particles (n,p,α)

Sunday Monday Thursday. Friday

Reactions of neutron-rich Sn isotopes investigated at relativistic energies at R 3 B

He+ 6,7 Li and 6 He+ 12 C reactions. Matko Milin Ruđer Bošković Institute Zagreb, Croatia

CENBG, Bordeaux, France 2. CEA/DAM-DIF, Arpajon,France 3. CEA/DEN, Saint Paul Lez Durance, France 4. IPN Orsay, Orsay France 5

Charged-particle evaporation from hot 164 Yb compound nuclei and the role of 5 He emission

Citation EPJ Web of Conferences (2014), provided the original work is prope

O WILEY- MODERN NUCLEAR CHEMISTRY. WALTER D. LOVELAND Oregon State University. DAVID J. MORRISSEY Michigan State University

Sub-barrier fusion of neutron-rich nuclei: 132 Sn+ 64 Ni

Introduc7on: heavy- ion poten7al model for sub- barrier fusion calcula7ons

Method of active correlations in the experiment 249 Cf+ 48 Ca n

Geiger-Marsden experiments: 100 years on

Measurement of the plasma astrophysical S factor for the 3 He(D, p) 4 He reaction in exploding molecular clusters

Decay studies of 170,171 Au, Hg, and 176 Tl

The Origin of the Elements between Iron and the Actinides Probes for Red Giants and Supernovae

Isospin symmetry breaking in mirror nuclei. Experimental and theoretical methods

PoS(NIC XII)184. Neutron capture on the s-process branch point nucleus 63 Ni

Performance of the ALICE Muon Trigger system in Pb Pb collisions

Role of projectile breakup effects and intrinsic degrees of freedom on fusion dynamics

Transcription:

: does the deformation of 28 Si play a role?, G. Montagnoli, M. Faggian, A. Goasduff, M. Mazzocco, F. Scarlassara, C. Stefanini, E. Strano, M. Urbani Dipartimento di Fisica e Astronomia, Universitá di Padova, and INFN, I-353 Padova, Italy. E-mail: giulia.colucci@studenti.unipd.it A. M. Stefanini, L. Corradi, E. Fioretto, F. Galtarossa, G. L. Zhang INFN, Laboratori Nazionali di Legnaro, I-35020 Legnaro, Italy. D. Bourgin, S. Courtin, F. Haas IPHC, CNRS-IN2P3, Université de Strasbourg, F-67037 Strasbourg Cedex 2, France. P. Colovic, S. Szilner Ruder Bošković Institute, HR-0002 Zagreb, Croatia. This contribution reports on the results of measurements of near- and sub-barrier fusion cross sections in the system performed at the Laboratori Nazionali di Legnaro of INFN. The 30 Si beam of the XTU Tandem accelerator in the energy range of 47-90 MeV, was delivered on a metallic 30 Si target (50 µg/cm 2 ) enriched to 99.64 % in mass 30. A beam electrostatic deflector was used for the detection of evaporation residues. The excitation function obtained for has been compared with the previous data on and Coupled Channels (CC) calculations using the M3Y+repulsion potential, taking into account the low lying 2 + and 3 excitations. Reproducing the low-energy excitation function was only possible by using a weak imaginary potential, probably simulating the oblate deformation of this nucleus. On the contrary, the data on are nicely reproduced by the CC calculations without any imaginary potential ( 30 Si has a spherical shape). The astrophysical S-factor does not show a maximum, so that there is no evidence for hindrance in the measured energy range. The logarithmic derivative of the two excitation functions highlights the difference between the two systems. Even above the barrier the two systems behave differently. This is best seen by comparing the two barrier distributions where the high energy peak observed for is not found for. This is presumably due to the stronger couplings present in 28 Si and further theoretical analyses are in progress. The 26th International Nuclear Physics Conference -6 September, 206 Adelaide, Australia Speaker. Also Dipartimento di Fisica e Scienze della Terra, Univ. di Ferrara, Ferrara, Italy. c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). http://pos.sissa.it/

. Introduction Several experiments have shown that nuclear structure is very important in heavy ion fusion dynamics near the barrier. In order to clearly indicate this sensitivity, a valuable and often employed tool is the comparison of fusion data for neighboring isotopes. Well-known results were obtained by comparing the fusion data for different nickel [, 2] and calcium [3] isotopes, where the influence of transfer channels in fusion reactions was shown, or in the fusion of 6 O with different samarium isotopes [4, 5], Figure : Measured fusion excitation functions of where the change of the structure from spherical vibrational to strongly deformed, 28 Si+ 30 Si [6] and previous data for [7]. The lines are the results of CC calculations (see text). nuclei was clearly observed. A comparative study of Si+Si systems [6] has revealed the influence of transfer on the fusion of the asymmetric case 28 Si+ 30 Si, well explained by including one- and successive two-neutron transfer channels in the coupling scheme, and highlighted a different behaviour of the two symmetric systems and (see Fig. ). These systems are very interesting because of the different shape of the silicon isotopes: indeed 30 Si is nearly spherical, whereas 28 Si is strongly deformed with an oblate shape. The excitation function of has an unusual behavior, because the cross sections are hindered just below the barrier, but this effect disappears at lower energies. It was further surprising that the low-energy data were well reproduced only by artificially applying a weak, short-ranged imaginary potential, probably simulating the effect of the oblate deformation. This feature has to be further investigated and the present study of has been very significant in this respect. No positive Q-value transfer channels are available in. This allows us to consider only couplings to vibrational states, neglecting the rotational ones which instead are present in the deformed nucleus 28 Si. This difference has an impact, as we shall see, on the energy dependence of fusion cross sections at sub-barrier energies for the two symmetric systems. In view of all this, it is very tempting to compare the two reactions. Unfortunately the lack of experimental fusion data below 4 mb for prevented to perform a meaningful comparison and also to obtain information on the presence of hindrance phenomenon. A fusion experiment was then performed at the Laboratori Nazionali di Legnaro (LNL) of INFN with the purpose to extend the data on. Here we present the results of this experiment on from well below to well above the Coulomb barrier and we compare the results with the previous data on.

2. Experimental set-up Fusion cross sections can be determined by direct detection of the evaporation residues (ER). The ER, which are forward peaked, must be physically separated from the direct beam and the intense flux of elastically scattered beam particles. In order to achieve this separation at 0 and at small angles, an electrostatic field perpendicular to the direction of the particles is applied between the target and the detection system. Such an electrostatic beam deflector has been used in several sub- and near-barrier experiments at LNL and it was employed also in the study of fusion reaction. This set-up is shown in Fig. 2 and allows fast and reliable measurements of relative and absolute cross sections. The electrostatic deflector assures the separation of Figure 2: The experimental set-up (up) and its schematic ER from the transmitted beam, while a design (down). The reaction chamber, the electrostatic deflector and the telescope E-TOF- E are drawn from the left. subsequent E-TOF- E telescope allows their identification. This telescope is made up of two micro-channel plates (MCP) detectors, a transverse field ionization chamber and a silicon detector placed in the same gas (CH 4 ) volume of the chamber. The ionization chamber provides a loss of energy signal E, while the two MCPs yield the time of flights together with the silicon detector. In this configuration the silicon detector measures the residual energy of the ER and provides both the trigger for data acquisition and the start signal for the time of flight. In Fig. 3 typical graphs reporting the time of flight (TOF) as a function of the residual energy E are shown. A good separation of ER events from the residual beam is achieved at energies both above and below the Coulomb barrier. Analogous two-dimensional spectra reporting the time of flight as a function of the loss of energy ( E) and the E vs. residual energy were used in the data analysis. The sliding seal reaction chamber allows to measure angular distributions. In the reaction chamber four silicon detectors are placed symmetrically around the beam direction at the same scattering angle θ lab = 6. They are used to monitor the beam and normalize to the Mott scattering cross section. During the experiment three angular distributions were measured at the energies of 58, 72 and 80 MeV in the range from -6 to +9. The total fusion cross sections were derived by integrating those distributions, and by simple interpolations or extrapolations for all other energies where ER measurements were taken only at 2 (or 3 for low energies). The XTU Tandem accelerator in LNL provided 30 Si beams in the energy range of 47-90 MeV, with intensities of 5-30 pna. The target consisted of 50 µg/cm 2 30 Si evaporated on 30 µg/cm 2 carbon 2

Figure 3: The graphs report the time of flight TOF (ordinate) as a function of the residual energy E (abscissa). The spectra were obtained above the Coulomb barrier (80 MeV, σ f us =527 mb, left panel) and at a sub-barrier energy of 50 MeV (σ f us =5 µb, right panel). backings facing the beam. The isotopic enrichment of 30 Si was 99.64 %, where the small amounts of 29 Si and 28 Si did not affect the fusion cross section because of their higher Coulomb barriers with respect to 30 Si. The carbon backing and the silicon target introduced an average beam energy loss of around 750-850 kev, which was taken into account in the analysis. The measured energy range allowed to extend the excitation function down to 4 µb, as shown Fig. 4 (left panel). The reported errors are statistical uncertainties, which reach a minimum value of - 2 % at the higher energies and increase to 0-20 % at energies below the barrier. The systematic errors on the absolute cross section are of ± 7-8 %, due to the geometrical solid angle uncertainties, angular distribution integrations and mainly to the deflector transmission ( 4 %). The fusion cross sections of the two symmetric sistems are reported in Fig. 4 (right panel). The two excitation functions behave differently at low energies with a flatter slope for. A coupled channels analysis is necessary to understand the origin of this difference. σ fus (mb) 000 00 0 0. 0.0 LNL Bozek σ fus (mb) 000 00 0 0. 0.0 0.00 0.00 24 28 32 36 40 44 0.000 22 24 26 28 30 32 34 36 38 Figure 4: Fusion excitation function of compared to previous experimental data [7] (left panel) and to system (right panel). 3. Coupled Channel calculation A preliminary coupled-channels (CC) analysis of the sub-barrier excitation function was per- 3

000 00 000 00 σ fus (mb) 0 0. 0.0 0.00 Ch W=0 Ch0 W=0 Ch0 W=5 25 30 35 40 45 σ fus (mb) 0 0. 0.0 0.00 0.000 Ch0 W=5 Ch0 W=0 Ch W=0 0.8 0.85 0.9 0.95.05..5.2 E /V (a.u.) B Figure 5: Fusion cross sections for compared to the no-coupling Ch calculation and to coupledchannels calculations at ten channels that are based on the predicted M3Y+rep potential with (Ch0 w5) and without (Ch0 w0) the addition of an imaginary potential (left panel). In the right panel the comparison with is reported including the coupled channel calculation adopted for this system, which employs a M3Y+rep potential with an imaginary part. formed, using the M3Y+repulsion potential [8]. The CC calculations made use of the so-called isocentrifugal approximation, which simplifies the calculations by reducing the number of channels. The density of the 30 Si nucleus that determines the potential was adjusted to optimize the fit to the fusion data. The complete CC calculation includes all one- and two-phonon excitations as well as mutual excitations of the low-lying 2 + and 3 states in both projectile and target, reaching a total of 5 channels. In order to simplify the calculations, the mutual excitation of the 2 + and 3 states in the same nucleus were ignored and so also the excitations of states above 0 MeV. These corrections allowed to neglegt the two-phonon and mutual excitations of the 3 states and reduced the number of channels to 0 (refereed to as Ch0). The 30 Si nucleus is nearly spherical, indeed the measured quadrupole moment of the 2 + state is Q 2 = -0.05(6) barn, which is consistent with zero. This feature makes the calculation Ch0 more robust than in the other case because the location of minimum of channels potential is essentially the same in all channels, as it will be explained later. The adopted coupling strengths are shown in Table I of [6]. As shown in the previous section, the comparison of the excitation functions of and indicates a different behavior at energies below the Coulomb barrier. The reason for this was found [6] in the CC analysis for the two systems and originates in the shape of the two silicon isotopes. Indeed, the deformation makes the channel potentials for the excited states of the two nuclei very different. In the case of 28 Si the minimum of the potential for the 2 + state is shifted to a smaller separation distance than the one observed in the entrance potential. In the case of 30 Si the minimum of all channel potentials are located essentially at the same separation distance. Since the incoming-wave boundary condition (IWBC) is imposed at the minimum of the entrance channel potential, this implies that for the system the fusion in the 2 + channel is cut off at an energy that is higher than the minimum of the 2 + channel potential, causing a suppression of the 4

0.3 0.25 0.2 MeV.5 MeV 0.35 0.3 0.25 MeV.5 MeV BD(E) (MeV - ) 0.5 0. 0.05 0 BD(E) (MeV - ) 0.2 0.5 0. 0.05-0.05 0-0. 24 26 28 30 32 34 36 38-0.05 22 24 26 28 30 32 34 36 38 Figure 6: Barrier distribution extracted for (left panel) and (right panel) system, with energy steps of.0 MeV (red points) and.5 MeV (blue points). calculated cross sections at low energies. By applying a weak imaginary potential it was possible to reproduce the low-energy data. On the contrary, beacuse of the spherical shape of 30 Si, the IWBC is imposed at the minimum of each channel potential, providing a more consistent treatment of the fusion in the different reaction channels. As a consequence, there is no need for an imaginary potential at low energy. This is shown in Fig. 5 where the results of CC calculation without and with the imaginary potential are reported. One sees that there is essentially no difference between the two calculations. 4. Barrier distribution The evaluation of the barrier distribution [9] for the system and its comparison with the one of is very useful in order to obtain more information on the role of the deformation of 28 Si. The barrier distributions were obtained using the three-point difference formula [0], with two energy steps of.0 MeV and.5 MeV, and a meaningful difference between the two cases has been observed. The system shows a single well-defined peak at low energies (Fig. 6, left panel), while a double-peak structure appears for with a clearly defined peak at higher energies (Fig. 6, right panel). These dissimilar behaviours arise from the different structures of the two nuclei and from the stronger couplings present in due to the oblate shape of 28 Si. 5. Hindrance The hindrance phenomenon is a suppressions of fusion cross sections with respect to standard CC calculations, occurring far below the Coulomb barrier [, 2]. The investigation of this effect usually makes use of two representations of the excitation functions, i.e. the astrophysical S-factor [3] and the logarithmic derivative L(E), which allow to reveal the presence of hindrance without using model calculations. Indeed, the S factor shows a maximum vs. energy in the presence of hindrance at the energy where the logarithmic slope L(E) reaches the value L CS expected for a constant astrophysical S factor. 5

3.5 3 0 7 0 6 L(E) (MeV - ) 2.5 2.5 L CS S (a.u.) 0 5 0 4 000 00 0.5 0 0.8 0.85 0.9 0.95.05..5.2 E /V (a.u.) B 0.8 0.9..2 E /V (a.u.) B Figure 7: In the left panel the logarithmic derivative of the excitation function of the two symmetric systems is reported, where L CS (black line) is the slope expected for a constant astrophysical S-factor. The S-factor extracted from experimental data for both cases is reported in the right panel. The investigation of hindrance in the fusion reaction is interesting also because the Q- value for the fusion of this system is positive. It has been pointed out that in such cases the S factor does not necessarily develop a maximum [2]. In the case of the S factor does not show any maximum vs. energy (see Fig. 7, right panel), even if it appears to saturate at the lowest measured energies. The logarithmic derivative does not cross the L CS curve (left panel). This allows us to conclude that there is no evidence of hindrance down to around 4 µb. The system shows a different behavior of the S factor and the slope. Indeed, there is evidence of hindrance at energies slightly below the Coulomb barrier [6] but the effect disappears at lower energies. On the other hand, the logarithmic derivative saturates far below the barrier. 6. Summary The fusion excitation function of has been measured in a wide energy range down to 4 µb. In this range there is no evidence for hindrance, when observing the logarithmic slope and S-factor. The comparison with shows that the two systems behave differently both near the barrier (the high energy peak observed in the barrier distribution of is not found for ) and farther below, where the regular trend of the heavier system can contrasted with the unusual one of. The weaker couplings in and its spherical shape make the difference. This is confirmed by coupled-channels calculations using the M3Y potential that are in good agreement with experimental data. At variance with the case of, the addition of a weak short-ranged imaginary potential is not necessary to reproduce the sub-barrier energy fusion cross sections. 7. Acknowledgments We are very grateful to the XTU Tandem staff and to M.Loriggiola for preparing targets of excellent quality. This work has been supported in the part by the Croatian Science Foundation under the project 794. 6

References [] M. Beckerman, M. Salomaa, A. Sperduto, J. D. Molitoris, and A. Di Rienzo, Sub-barrier fusion of 58,64 Ni with 64 Ni and 74 Ge, Phys. Rev. C 25 (982) 837. [2] M. Beckerman et al., Dynamic Influence of Valence Neutrons upon the Complete Fusion of Massive Nuclei, Phys. Rev. Lett. 45 (980) 472. [3] G. Montagnoli et al., Fusion of 40 Ca+ 40 Ca and other Ca+Ca systems near and below the barrier, Phys. Rev. C 85 (202) 024607. [4] R. Stokstad et al., The sub-barrier fusion of 40 Ar with 44,48,54 Sm, Z. Phys. A 295 (980) 269. [5] H. Esbensen, Fusion and zero-point motion, Nucl. Phys. A 352 (982) 47. [6] G. Montagnoli et al., Fusion of 28 Si+ 28,30 Si: Different trends at sub-barrier energies, Phys. Rev. C 90 (204) 044608. [7] E. Bozek et al., De-excitation of 58,60,62 Ni compound nuclei formed via symmetric and asymmetric entrance channels, Nucl. Phys. A 45 (986) 7. [8] Ş. Mişicu and H. Esbensen, Signature of shallow potentials in deep sub-barrier fusion reactions, Phys. Rev. C 75 (2007) 034606. [9] N. Rowley, G. R. Satchler, and P. H. Stelson, On the distribution of barrier interpretation of heavy-ion fusion, Phys. Lett. B 254 (99) 25. [0] H. Timmers et al., A case study of collectivity, transfer and fusion enhancement, Nucl. Phys. A 633, (998) 42. [] C. L. Jiang, et al., Unexpected Behavior of Heavy-Ion Fusion Cross Sections at Extreme Sub-Barrier Energies, Phys. Rev. Lett. 89 (2002) 05270. [2] C. L. Jiang, H. Esbensen, B. B. Back, R. V. F. Janssens, and K. E.Rehm, Analysis of heavy-ion fusion reactions at extreme sub-barrier energies, Phys. Rev. C 69 (2004) 04604. [3] E.Margaret Burbidge, G.R. Burbidge, William A. Fowler and F. Hoyle, Synthesis of the Elements in Stars, Rev. Mod. Phys. 29 (957) 547. 7