1 Introduction
Proton radioactivity, experimentally observed as a decay from the ground state, at GSI in 1981, has provided very important information on the structure of nuclei beyond the proton drip-line. The more complicated decaying mode, two-proton radioactivity, proposed 50 years ago in a classical article[1] opened a new window to investigate nucleon-nucleon correlations and the structure of atomic nuclei. In 2002, the simultaneous emission of two-protons was for the first time observed in the decay of 45Fe by Pfutzner, Giovinazzoin experiments at GSI and GANIL[2,3]. Research in the field flourished after this breakthrough, and to date 54Zn[4], 48Ni[5], 19Mg[6], 16Ne[7], 17Ne[8], 18Ne[9], 10C[10], 14O[11] and 29S[12] have been found to exhibit two-proton emission. Several theoretical approaches such as Diproton model[13,14], R-matrix approach[15], continuum shell model[16], adiabatic hyperspherical approach[17], and the quantum three body cluster approach[18], where the tunneling through the barrier is treated in a dynamical way, were applied to the problem.
There are two different decay modes for simultaneous two-proton emission: (1) three-body direct breakup involving an uncorrelated emission of the two protons, usually referred to as democratic emission. (2) 2He cluster emission where a pair of protons, correlated in a quasi-bound 1S configuration, breakup, when emitted into two protons (diproton emission). The two protons have strong angular and energy correlations. The 2He appears as a resonance at 20 MeV/c in the two-proton relative momentum distribution[19]. The microscopic calculations for the one- and two-proton decays of the 6.15 MeV 1- state of 18Ne had been presented in the Ref.[20]. It was found that for the two-proton the sequential decay through a ghost of the 1/2+ state is within a factor of three of the observed width obtained with the assumption of democratic decay. The calculated width for diproton emission is only about a factor of two smaller than that for sequential decay indicating that the observed decay may be a combination of the two processes. In the excitation-energy spectrum of 18Ne in the Ref.[9], it’s strange that some states can be seen in the two-proton emission 18Ne→16O+2p channel and not in the one-proton emission 18Ne→17F+p channel. That means that in these states we cannot find the 17F in ground state. So the sequential decay for two protons is most likely to occur in these states.
In this paper we present the microscopic shell-model calculations for sequential two-proton decay from excited states in 18Ne by some different Hamiltonians.
2 Calculation and discussion
The spectroscopic factor is the most important quantity needed to obtain the decay width. In order to calculate it, we perform a shell-model calculation to get the wave functions for 18Ne. The model space was used, including the 0s, 0p, 1s0d and 1p0f orbits. 16O is treated as a s4p12 closed shell, and the low-lying positive parity states of 17F and 18Ne are taken as s4p12(sd)1 and s4p12(sd)2. The low-lying negative parity states of 17F and 18Ne are treated as 1ħω excitations of the form s4p11(sd)2, s4p12(pf)1 and s4p11(sd)3, s4p12(sd)1(pf)1. So the emitted protons in the 18Ne and 17F are coming from (sd)(pf) shells. Two Hamiltonians designed for those types of model space are chosen for calculating the wave functions, namely the WBP and WBT interactions[21]. We use a simply shell-model code by our group, in this code the spurious states are removed by the usual method[22] by adding a center-of-mass Hamiltonian to the interaction.
The calculated excited energies of these low-lying states are shown in Fig.1. Some states are in reasonable agreement with the energies found in 18Ne. The low-lying negative states are dominated by the s4p11(sd)3 configuration, but the smaller s4p12 (sd)1(pf)1 component is the one responsible for one- and two-proton decay. The shell-model spectroscopic factors are obtained by the wave functions of 18Ne and17F. The decays from the positive states of 18Ne to the positive states of 17F and from the negative states of 18Ne to the negative states of 17F can go by 0d-shell wave emission or 1s-shell wave emission. The decays from the positive states of 18Ne to the negative states of 17F and from the negative states of 18Ne to the positive states of 17F can go by 0f-shell wave emission or 1p-shell wave emission. Because the s4p11(sd)3 component in 18Ne is quite larger than that of s4p12(sd)1(pf)1, the spectroscopic factors are larger in the channel of positive states in 18Ne.
According to the scattering theory, the half-life for decay from initial state i to a final state f by one particle emission is given by:
where the decay width can be found from the relation[25,26]:
The total width for decay is a sum of partial widths:
The branching ratios are simply the ratio between a partial decay width and the total one:
For the fourth 1- state at 7.94 MeV in 18Ne, we find that the spectroscopic factors decaying to the 1/2- third excited state (Q1p = 0.914 1 MeV) is quite larger than that decaying to the 5/2+ ground state (Q1p = 4.018 4 MeV) of 17F. The spectroscopic factors and the widths for each of the channels are shown in Table 1. In this table, we can find the branch ratio that decays to 1/2- state is larger than those decays to the ground state and the first ground state because it has larger spectroscopic factor even though it has smaller single-particle width. We can conclude that the 7.94 MeV 1- state is most likely to be the best candidate for two-proton sequential decay. That is why it can be seen in the two-proton emission channel and not in the one-proton one. There are other states in this situation, like the 3- state around 9–10 MeV (9.809 MeV for WBP and 10.099 MeV for WBT) and the 5- state near 13 MeV (13.412 MeV for WBP and 13.200 MeV for WBT).
WBP | WBT | Expt. | ||
Ex / MeV | 7.648 0 | 7.698 6 | 7.94 | |
Channel | Spectroscopic factor | |||
WBP | WBT | Γsp / keV | ||
5/2+⊗0f7/2 | 0.010 61 | 0.003 52 | 129 | |
5/2+⊗0f5/2 | 0.010 11 | 0.007 65 | 101 | |
5/2+⊗1p3/2 | 0.001 01 | 0.004 28 | 2 818 | |
1/2+⊗1p3/2 | 0.003 34 | 0.002 23 | 2 239 | |
1/2+⊗1p1/2 | 0.000 13 | 0.000 13 | 2 188 | |
1/2-⊗0d3/2 | 0.001 06 | 0.005 25 | 2 | |
1/2-⊗1s1/2 | 0.091 86 | 0.240 19 | 122 | |
ΓWBP | ΓWBT | BrWBP | BrWBT | |
5/2+⊗(0f+1p) | 5 | 13 | 0.216 | 0.277 |
1/2+⊗1p | 8 | 5 | 0.32 | 0.11 |
1/2-⊗(0d+1s) | 11 | 29 | 0.464 | 0.613 |
Total Γ | 24 | 47 | Expt. ≤ 50 keV |
3 Conclusion
We have presented some preliminary results of the proton decay branch ratios and decay width in 18Ne using a shell-model calculation. The results obtained for the branch ratio from the 1- state at 7.94 MeV in 18Ne show that this state is most likely to be a candidate for sequential two-proton decay. The 3- and the 5- states can also be candidates for the same process.