Semi Empirical Formula For Neutrinoless Double Beta Decay

A Semi empirical formula for both phase space factor and Nuclear Matrix Element (NME) is developed for neutrinoless double beta decay, and the formula is used to compute the neutrinoless double beta decay half lives. The computed half lives for neutrinoless double beta decay are compared with the corresponding experimental values and with those predicted by QRPA model. The semi empirical formula predictions are found to be in good agreement with experimental data. The semi empirical formula is used to predict neutrinoless double beta decay of various isotopes Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm that exhibiting single beta decay. As our semi empirical formula predictions agree with the experimental data we hope that the present work will be useful for the future experiments. Keyword: Neutrinoless double beta decay, Nuclear matrix element


INTRODUCTION
Double beta decay is a radioactive decay process where a nucleus releases two beta rays as a single process. Here two neutrons in the nucleus are converted in to two protons and in the process two electrons and two electron antineutrinos are emitted. In order for beta decay to be possible the final nucleus must have larger binding energy than the original nucleus. Double beta decay is difficult to study in most practically interesting cases, because both beta decay and double beta decay are possible, with probability favouring beta decay. The double beta decay is usually studied only for beta stable nuclei. Like single beta decay, double beta decay does not change the mass number A. More than 60 naturally occurring isotopes are capable of undergoing double beta decay. violates lepton number conservation and is therefore forbidden in the standard electroweak theory. According to this theory neutrinos are massless. The observation of neutrino mass and oscillation is a clear example of a phenomenon at variance with the standard model.
There are different models for explaining the double beta decay process. Among them, two methods are mainly used to calculate Nuclear Matrix Elements (NME) for 2β(0ν) decays. One is the family of Quasi Particle Random Phase Approximation (QRPA) [21]. This method has been used by different groups and varieties of techniques are employed with results for most of the possible emitters [24]. The other method concerned to double beta decay process is the Interacting Shell Model (ISM) [8]. It has been shown that as the difference in deformation between parent and daughter grows, the NME's of both the neutrinoless and two neutrino mode decreases rapidly.
The interest in double beta decay spans more than six decades. In 1937 Racah [16] following the fundamental suggestion of Majorana [11], discussed the possibility of a neutrinoless transformation of two neutrons into two protons plus two electrons. Even earlier Geoppert-Mayer [28] evaluated the decay rate of 2β(2ν) mode and realized that the corresponding half lives could exceed 10 20 years. Furry [39] shortly afterwards estimated that 2β(0ν) should be much faster than 2β(2ν) decay. Thus the stage was set for the realization that observation of the 2β(0ν) decay would establish that the neutrino is a massive Majorana particle. In 1982 J. Schechter-Valle [23], while considering 2β(0ν) decay, suggested the existence of Majorana mass of the neutrino in the frame work of Gauge theories. In 1984 Fiorini et al [9] introduced a program to develop low temperature detectors for 2β decay search. Next year Doi et al [27] made a fundamental theoretical analysis of 2β decay to obtain the main formulae for probability of decay, energy and angular electron spectra. In 1986, using QRPA model Vogel et al [31] got satisfactory agreement between theoretical and experimental 2β(2ν) half life values. In the review work of Gmez-Cadenas et al [22], the authors provide an answer to the fundamental question of whether neutrinos are Dirac or Majorana. Neutrinoless double beta decay is of great interest for studying the fundamental properties of neutrino beyond the standard electro-weak theory. High sensitivity 2β(0ν) studies are the unique and practical ways for studying  [4,5,7,20,32].
For the double beta decay processes, two crucial ingredients are the phase space factors and the Nuclear Matrix Elements (NME). A general theory of phase space factors was developed by Doi et al. [25,26] following the previous work of Primakoff and Rosen [17], and Konopinski [10]. It was reformulated by Tomoda [37] by approximating the electron wave functions at the nuclear radius and without inclusion of electron screening. Accurate values of the phase space factors are necessary ingredients for theorists, to improve the double beta decay lifetime predictions and for experimentalists to plan their set-ups [30]. The Nuclear Matrix Element depends on the nuclear structure of the nuclei involved in the decay. Frank T Avignone et al [14] have done a detailed study on nuclear matrix elements and also pointed out the increasing sensitivity of experiments and improvements in nuclear theory make the future exciting for this field at the interface of nuclear and particle physics. The expression for Nuclear Matrix Element can be written in general as the sum of three components [29] as , are the Gamow-Teller, Fermi and tensor components respectively. g A is the axial vector coupling constant and g ν is the vector coupling constant. The present work aims to develop a semi empirical formula for both phase space factor and Nuclear Matrix Element for computing the neutrinoless double beta decay half life. By using this formula we would like to predict the possibility of 2β(0ν) decay from various isotopes exhibiting single beta decay. The details of the semi empirical formula are given in Section 2 and results, discussion and conclusion are given in Section 3.

THE SEMI EMPIRICAL FORMULA
In the standard scenario, when 2β(0ν) decay process occurs by exchange of light Majorana neutrinos between two nucleons inside the nucleus, and in the presence of left handed weak interactions, the life time expression can be written as a product of three factors and is given as [3] T G M m m e where G 0ν is the phase space factor for this decay mode, 〈 〉 m ν is the effective neutrino mass parameter, m e is the electron mass and M 0ν are the Nuclear Matrix Elements depending on the nuclear structure of the nuclei involved in the decay.
The phase space factor depends on the decay energy Q and nuclear charge Z and in the present investigation we studied the dependence of phase space factor with ZQ 3 and Z 2 Q 6 for various isotopes undergoing neutrinoless double beta decay. From the observed dependence of phase space factor taken from ref [34], on ZQ 3 and Z 2 Q 6 , we have developed a semi empirical formula for the phase space factor. Using ZQ 3 , Z 2 Q 6 and Z 3 Q 9 as variables, a new formula is obtained and is given as, The constants are, a E = − 2 48904 26 Due to the two-body nature of the transition operator, the NMEs can also be expressed as a sum of product of two-body transition densities (TBTDs) and matrix elements of the two-body transition operators for two-particle states.
We have studied the dependence of nuclear matrix element values taken from [38] with Z -1/3 for various isotopes undergoing neutrinoless double beta decay and a new formula is obtained by making least-squares fit to the nuclear matrix elements data and is given as, The constants are, The comparison of the computed nuclear matrix elements using the present formula with the values of Ref [38] and comparison of computed phase space factor with the values of Ref [34] Ref. [

RESULTS, DISCUSSION AND CONCLUSION
The Q value for double beta decay of mother nuclide with mass m m to the daughter nuclide with mass m d is given by the mass difference Q m m m d = − In the present work Q values are computed using the experimental binding energies of Audi and Wapstra [15]. The present empirical formula is applied for all the observed neutrinoless double beta decay isotopes. Column 7 of Table 1 represents the computed half-lives for neutrinoless double beta decay of various isotopes and is compared with the experimental values given in column 8 and QRPA values [38] in column 9. It is found from the table that our formula predictions are in good agreement with the experimental values and the QRPA values. The value of <m ν > is taken as 50MeV and is obtained from Rodin et al [38]. We have applied the present formula for computing the phase space factor, Nuclear Matrix Element and half lives for various isotopes that exhibiting single beta decay. Tables 2 represents the computed Q values, Phase space factors, Nuclear Matrix Elements and half lives for neutrinoless double beta decay of various Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm isotopes. As our semi empirical formula prediction agree with the experimental data we hope that our prediction on neutrinoless double beta decay of various Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm isotopes will be a guide for future experiments.