1 Introduction
Besides the well-known α, β, and γ radioactivity decays, exotic radioactivity modes also exist in very proton-rich nuclei [1-4]. Two-proton emission is one of the most interesting phenomena in nuclei beyond or close to the proton drip line [5, 6]. Generally speaking, there are three different mechanisms for proton-rich nuclei to emit two protons: (i) two-body sequential emission, (ii) three-body simultaneous emission, and (iii) diproton emission (also called 2He cluster emission). The third mode is an extreme case with the emission of two strongly correlated protons. The 2He cluster can only exist for a short while and then separates after penetrating the Coulomb barrier.
The two-particle momentum correlation is influenced by the nuclear force between two particles [7]; consequently, the proton-proton momentum correlation plays an important role in the emission mechanism and causes the two-proton relative momentum (qpp) and opening angle (θpp) to be quite different compared with other emission mechanisms [8, 9]. Generally, the two-proton correlation of diproton emission is much stronger than that of the other two mechanisms. Additionally, the emission time difference between two protons for sequential emission is long, compared to the other two modes. In [10], the three-body decay of two excited proton-rich nuclei, namely 23Al p + p + 21 Na and 22Mg p + p + 20 Ne, has been measured at the RIKEN RI Beam Factory. It has been noted that the emission mechanism of the Isospin Analogue State (IAS) of 22Mg has strong diproton emission probability, based on the analysis of (qpp) and (θpp). However, it is difficult to determine the emission mechanism of the excited 23Al. For three-body simultaneous emission, the two protons are emitted almost at the same time, while the emission time of two-body sequential emission is quite different. Comparing the experimental data with the theoretical simulations, the source size and proton emission time can be extracted [11-13]. In [14], the p-p momentum correlation function (Cpp(q)) was studied for these two decay channels, and the source size and emission time information were extracted, as well as for the emission mechanisms of two protons from 23Al and 22Mg.
There are two main factors that affect the p-p momentum correlation function. One is the source size, the other is the emission time difference between the two protons. An increase of the source size will decrease the strength of Cpp(q), and the emission time difference will also have a similar effect on Cpp(q). However, the source size and emission time determining Cpp(q) are unique, or a different combination of source size and emission time could give the same Cpp(q). In this study, the p-p momentum correlation function was investigated systematically by using the code Correlation After Burner (CRAB), which is a widely used method for calculating the momentum correlation function in nuclear Physics [15]. Assuming the first proton being emitted at time t = 0 and the second proton being emitted at time t, a Gaussian source with the form
2 Calculation results
2.1 Spherical Gaussian Source
We first calculated the p-p momentum correlation function Cpp(q) by assuming a spherical Gaussian source. In Fig. 1, the results of Cpp(q) with source size of r0=0.5 fm and r0=2.5 fm at different values of τ are presented. The figure shows that Cpp(q) decreases as τ increases for fixed r0. For larger r0, the difference of Cpp(q) between different τ becomes larger.
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F001.jpg)
In Fig. 2, we can see that Cpp(q) decreases as r0 increases for a fixed value of τ. For larger τ, the difference of Cpp(q) between different r0 becomes smaller. In these figures, Cpp(q) has the maximum value at around qpp = 20 MeV/c.
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F002.jpg)
The proton-proton momentum correlation function Cpp(q) increases with qpp and saturates at around 1. Two parameters were used to describe Cpp(q) for studying the effect of source size and emission time systematically. One is the maximum value of Cpp(q) at approximately qpp = 20 MeV/c (Cmax(q)), the other is the full width at half maximum (FWHM) of Cpp(q) determined by the difference of the two qpp values of Cpp(q)=[Cmax(q)-1]/2 located at the left and right side of the Cpp(q) maximum. The r0 dependence of Cmax(q) and FWHM for the p-p momentum correlation function with different τ are given in Fig. 3. As shown in Fig. 3(a), Cmax(q) decreases gradually with increasing r0. For a large τ, the change of Cmax(q) is very small. As shown in Fig. 3(b), the FWHM is inversely proportional to r0. For different τ, the behavior of the FWHM with r0 is very similar.
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F003.jpg)
Similarly, the τ dependence of Cmax(q) and the FWHM for the p-p momentum correlation function at different r0 are shown in Fig. 4. The dependence of Cmax(q) on τ is quite similar with r0, except for the FWHM results. The difference of the FWHM values is much larger for different source sizes.
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F004.jpg)
Since both r0 and τ affect the correlation function, it is interesting to see whether or not different r0 and τ combinations result in the same proton-proton momentum correlation function Cpp(q). To find the answer, contour plots of Cmax(q) and the FWHM extracted from Cpp(q) at different r0 and τ values are given in Fig. 5, in which Cmax(q) and the FWHM are the Z axis. From this figure, we can see clearly that each isoline of Cmax(q) or the FWHM has only one intersection point with each other. This indicates that a set of Cmax(q) and the FWHM values or one proton-proton momentum correlation function has only a uniquely determined r0 and τ combination. This is due mainly to the monotonic dependence of Cmax(q) or the FWHM on r0 and τ, respectively. Based on the above results, it is shown that the source size r0 and proton emission time difference τ could be determined uniquely by fitting the experimental Cpp(q) with the CRAB calculation, as demonstrated in [14].
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F005.jpg)
2.2 Ellipsoidal Gaussian Source
Considering the possible deformation of the nucleus, we calculated the proton-proton correlation function Cpp(q) for a deformed nucleus through the CRAB code. The deformation is described by the parameter ϵ=ΔR/R, where R is the nuclear radius with no deformation and ΔR is the difference in radius before and after the deformation.
After considering the deformation of the source using CRAB, the calculated Cpp(q) results are shown in Fig. 6. The Cmax(q) and the FWHM of the ellipsoidal Gaussian source and the spherical Gaussian source are almost same in the range of ϵ from -0.10 to 0.10. In fact, the Cpp(q) of the ellipsoidal Gaussian source and the spherical Gaussian source are almost identical, with the same effective source radius r0. This indicates that the deformation (not very large) of the nucleus has little effect on the p-p momentum correlation function.
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F006.jpg)
2.3 Double Gaussian Source
The α cluster structure is one of the most common aspects of a nucleus [18]. If two protons are emitted from this kind of nucleus, the protons may come from the same or different α cluster within it. It would be interesting to see the effect of the α cluster on the p-p momentum correlation function. To study the cluster structure in the nucleus, a double Gaussian source was used in CRAB to simulate two clusters inside the nucleus.
We assume that the source of two protons is not distinguished. Thus, the two emitted protons may come from the same cluster or from two different clusters. Define ΔCmax(q) and ΔFWHM as the difference of Cmax(q) and FWHM between the spherical Gaussian source and the double Gaussian source with the same effective source size r0. The results of ΔCmax(q) and ΔFWHM are presented in Fig. 7. In Fig. 7(a), ΔCmax(q) first increases and then decreases with r0, and the maximum value appears near r0=1.5 fm, i.e., the p-p momentum correlation function has the largest difference for two protons emitted from ordinary nuclei and from nuclei with clusters. ΔCmax(q) has the largest value when the source size is near 1.5 fm. As shown in Fig. 7(b), ΔCmax(q) decreases as τ increases. For different r0, the ΔCmax(q) values are very large when τ is small, but they are quite close when τ is large enough. We can also see that ΔFWHM gradually increases with r0, but there is almost no change with the increase of τ, as shown in Fig. 7(c) and Fig. 7(d).
-202005/1001-8042-31-05-010/alternativeImage/1001-8042-31-05-010-F007.jpg)
ΔCmax(q) decreases with increasing τ, while ΔFWHM does not change with τ, which indicates that the difference of Cpp(q) between the double Gaussian source and the spherical Gaussian source decreases and nears the same value with increasing τ. ΔCmax(q) decreases with increasing r0, while ΔFWHM gradually increases with r0, which indicates that the difference of Cpp(q) between the double Gaussian source and the spherical Gaussian source becomes significantly larger as r0 increases. These results indicate that the p-p momentum correlation function for a source with or without cluster structure will have large systematical differences with the variance of r0 and τ. This may provide a possible method for experimentally observing the cluster structure in proton-rich nuclei. For practical applications, further investigations are necessary.
3 Summary and Outlook
In summary, the proton-proton momentum correlation functions (Cpp(q)) for a sphere, ellipsoid, and double Gaussian source were investigated using the CRAB code. From systematical studies of the varied effects of the source size r0 and emission time τ on Cpp(q), it was shown that one Cpp(q) distribution corresponds to unique values of r0 and τ. There is almost no difference in Cpp(q) between a spherical Gaussian source and an ellipsoidal Gaussian source, i.e., a small nuclear deformation has very little effect on Cpp(q). The proton-proton momentum correlation function has the largest difference between ordinary nuclei and clustered nuclei when the source size is near 1.5 fm; this may give a possible method for experimentally observing the cluster structure of proton-rich nuclei. Recently, artificial neural networks have been widely used in the research of many practical problems that are difficult for modern computers to solve [19-22]. Extracting the source size and emission time of two particles from experimental data is relatively difficult. It may be interesting to study systematically the p-p momentum correlation functions by using artificial neural networks in future studies.
Radioactive decays at limits of nuclear stability
. Rev. Mod. Phys. 84, 567 (2012). https://dx.doi.org/10.1103/RevModPhys.84.567Two-proton radioactivity
. Rep. Prog. Phys. 71, 046301 (2008). https://dx.doi.org/10.1088/0034-4885/71/4/046301Landscape of two-proton radioactivity
. Phys. Rev. Lett. 110, 222501 (2013). https://dx.doi.org/10.1103/PhysRevLett.110.222501Observation of long-range three-body coulomb effects in the decay of 16Ne
. Phys. Rev. Lett. 113, 232501 (2014). https://dx.doi.org/10.1103/PhysRevLett.113.232501On neutron-deficient isotopes of light nuclei and the phenomena of proton and two-proton radioactivity
. Nucl. Phys. 19, 482 (1960). https://dx.doi.org/10.1016/0029-5582(60)90258-3Implantation-decay method to study the β-delayed charged particle decay
. Nucl. Sci. Tech. 29, 98 (2018). https://dx.doi.org/10.1007/s41365-018-0438-5Measurements of momentum correlation and interaction parameters between antiprotons
. Nucl. Sci. Tech. 27, 152 (2016). https://dx.doi.org/10.1007/s41365-016-0147-xTwo-proton emission from the ground state of 12O
. Phys. Rev. Lett. 74, 860 (1995). https://dx.doi.org/10.1103/PhysRevLett.74.860Experimental evidence of 2He decay from 18Ne excited states
. Phys. Rev. Lett. 100, 192503 (2008). https://dx.doi.org/10.1103/PhysRevLett.100.192503Different mechanism of two-proton emission from proton-rich nuclei 23Al and 22Mg
. Phys. Lett. B 743, 306 (2015). https://dx.doi.org/10.1016/j.physletb.2015.02.066Observation of lifetime effects in two-proton correlations for well-characterized sources
. Phys. Rev. Lett. 71, 2863 (1993). https://dx.doi.org/10.1103/PhysRevLett.71.2863Correlations and characterization of emitting sources
. Eur. Phys. J.A 30, 81 (2006). https://dx.doi.org/10.1140/epja/i2006-10109-6Two-pion correlations in heavy ion collisions
. Phys. Rev. C 29, 2173 (1984). https://dx.doi.org/10.1103/PhysRevC.29.2173Proton-proton correlations in distinguishing the two-proton emission mechanism of 23Al and 22Mg
. Phys. Rev. C 94, 044621 (2016). https://dx.doi.org/10.1103/PhysRevC.94.044621Testing transport theories with correlation measurements
. Nucl. Phys. A 566, 103c (1994). https://dx.doi.org/10.1016/0375-9474(94)90614-9A theoretical study on different cluster configurations of the 9Be nucleus by using a simple cluster model
. Nucl. Sci. Tech. 28, 86 (2017). https://dx.doi.org/10.1007/s41365-017-0239-2Natural orbital description of the halo nucleus 6He
. Nucl. Sci. Tech. 28, 179 (2017). https://dx.doi.org/10.1007/s41365-017-0332-6Nuclear clusters and nuclear molecules
. Phys. Rep. 432, 43 (2006). https://dx.doi.org/10.1016/j.physrep.2006.07.001Recovery of saturated signal waveform acquired from high-energy particles with artificial neural networks
. Nucl. Sci. Tech. 30, 148 (2019). https://dx.doi.org/10.1007/s41365-019-0677-0Slice-wise reconstruction for low-dose cone-beam CT using a deep residual convolutional neural network
. Nucl. Sci. Tech. 30, 59 (2019). https://dx.doi.org/10.1007/s41365-019-0581-7Determination of Gamma point source efficiency based on a back-propagation neural network
. Nucl. Sci. Tech. 29, 61 (2018). https://dx.doi.org/10.1007/s41365-018-0410-4Particle dispersion modeling in ventilated room using artificial neural network
. Nucl. Sci. Tech. 28, 5 (2017). https://dx.doi.org/10.1007/s41365-016-0159-6