Semi Empirical Formula For Neutrinoless Double Beta Decay

Semi Empirical Formula For Neutrinoless Double Beta Decay Abstract 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. Double beta decay is of two types; the two neutrino and neutrinoless double beta decay. The two neutrino double beta decay [2ÃŽ ²(2ÃŽ ½)] which involves the transformation of two neutrons into two protons conserves not only the electric charge but also the lepton number. On the other hand neutrinoless double beta decay [2ÃŽ ²(0ÃŽ ½)] 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) [1]. This method has been used by different groups and varieties of techniques are employed with results for most of the possible emitters [2]. The other method concerned to double beta decay process is the interacting shell model (ISM) [3]. 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 [4] following the fundamental suggestion of Majorana [5], discussed the possibility of a neutrinoless transformation of two neutrons into two protons plus two electrons. Even earlier Geoppert-Mayer [6] evaluated the decay rate of 2ÃŽ ²(2ÃŽ ½) mode and realized that the corresponding half lives could exceed 1020years. Furry [7] 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 while regarding 2ÃŽ ²(0ÃŽ ½) decay suggested the existence of Majorana mass of the neutrino in the frame work of Gauge theories [8]. In 1984 Fiorini et al [9] introduced a program to develop low temperature detectors for 2ÃŽ ² decay search. Next year Doi et al [10] made a fundamental theoretical analysis of 2ÃŽ ² deca y to obtain the main formulae for probability of decay, energy and angular electron spectra. In 1986 using QRPA model Vogel et al [11] gave satisfactory agreement between theoretical and experimental 2ÃŽ ²(2ÃŽ ½) half life values. 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 the Majorana nature of neutrinos, the neutrino mass spectrum, the absolute neutrino mass scale, the majorana CP phases and other fundamental properties of neutrinos in the foreseeable future. The first experiment [12] to claim 2ÃŽ ²(0ÃŽ ½) is the Klapdor, HM experiment done in the year 2001. Numerous experiments like COBRA, GERDA etc have been carried out to search neutrinoless double beta decay and 48Ca, 76Ge, 82Se, 96Zr, 100Mo, 116Cd, 128Te, 150Nd, 238U are some of the isotopes exhibiting neutrinoless double beta decay [13-17]. 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. [18, 19] following the previous work of Primakoff and Rosen [20], and Konopinski [21]. It was reformulated by Tomoda [22] by approximating the electron wave functions at the nuclear radius and without inclusion of electron screening. The Nuclear Matrix Element depends on the nuclear structure of the nuclei involved in the decay. The expression for Nuclear Matrix Element can be written in general as the sum of three components [23] as (1) Where, , , are the Gamow-Teller, Fermi and tensor components respectively. is the axial vector coupling constant and 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 [24] (2) Where G0ÃŽ ½ is the phase space factor for this decay mode, is the effective neutrino mass parameter, me is the electron mass and M0ÃŽ ½ are the Nuclear Matrix Elements depending on the nuclear structure of the nuclei involved in the decay. The phase space factor depends on the energy decay QÃŽ ²ÃŽ ² and nuclear charge Z and studied the dependence of phase space factor with ZQ3 and Z2Q6 for various isotopes undergoing neutrinoless double beta decay. From the observed dependence of phase space factor taken from ref [25], with ZQ3 and Z2Q6 we have developed a semi empirical formula for the phase space factor. Using ZQ3, Z2Q6 and Z3Q9 as variables, a new formula is obtained and is given as, (3) The constants are, , , 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 [26] 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, (4) The constants are, , ,, , , The comparison of the computed nuclear matrix elements using the present formula with the values of Ref [26] and comparison of computed phase space factor with the values of Ref [25] are shown in Table 1. Results, discussion and conclusion The Q value for double beta decay of mother nuclide with mass mm to the daughter nuclide with mass md is given by the mass difference [27] which in turn can be written as a function of frequency ratio and the electron mass me; (5) In the present work Q values are computed using the experimental binding energies of Audi and Wapstra [28]. 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 [26] 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 ÃŽ ½> is taken as 50meV and is obtained from Rodin et al [25]. 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.  ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­Ã‚ ­ Table 1. The computed, and for neutrino less double beta decay of various isotopes and their comparison with the experimental, QRPA and Ref [26] values _____________________________________________________________________________________ _____________________________________________________________________________________ _____________________________________________________________________________________ Table 2. The computed Q values, Phase space factors, nuclear matrix elements and the predicted half lives for neutrino less double beta decay of various Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm, Gd and Pt isotopes ____________________________________________________________ ____________________________________________________________ _____________________________________________________________ Table 2. Continued†¦.. ____________________________________________________________ ____________________________________________________________ 103Mo 6408.8 1.10330E-11 3.47678 7.83E+23 104Mo 7759.0 6.71620E-11 3.47678 1.29E+23 105Mo 8588.0 1.72446E-10 3.47678 5.01E+22 106Mo 10067 7.43636E-10 3.47678 1.16E+22 107Mo 11430 2.37009E-09 3.47678 3.65E+21 109Pd 901.0 1.92669E-15 2.56879 8.22E+27 110Pd 2004.0 7.81588E-15 2.56879 2.03E+27 111Pd 3253.5 2.48567E-14 2.56879 6.37E+26 112Pd 4244.5 2.49383E-13 2.56879 6.35E+25 113Pd 5359.3 2.60657E-12 2.56879 6.07E+24 114Pd 6523.9 1.74643E-11 2.56879 9.06E+23 115Pd 7690.5 8.20373E-11 2.56879 1.93E+23 116Pd 8759.0 2.73683E-10 2.56879 5.78E+22 117Pd 9895.0 8.38279E-10 2.56879 1.89E+22 118Pd 11239 2.67934E-09 2.56879 5.91E+21 114Cd 540.1 1.30076E-15 2.37189 1.43E+28 115Cd 1945.5 7.60027E-15 2.37189 2.44E+27 116Cd 2809.1 1.46594E-14 2.37189 1.27E+27 117Cd 3975.0 1.46594E-14 2.37189 1.25E+26 118Cd 4947.1 1.48733E-13 2.37189 1.37E+25 119Cd 6158.4 1.35555E-12 2.37189 1.61E+24 120Cd 7131.1 1.15360E-11 2.37189 4.01E+23 121Cd 8144.1 4.62494E-11 2.37189 1.17E+23 122Cd 9215.9 1.59314E-10 2.37189 3.73E+22 123Cd 10510.5 4.97662E-10 2.37189 1.12E+22 124Cd 11526.8 1.65584E-09 2.37189 4.84E+2