Deformed Shell Model Applications to Weak Interaction Processes

: The deformed shell model (DSM), based on Hartree–Fock intrinsic states with angular momentum projection and band mixing, has been found to be quite successful in describing many spectroscopic properties of nuclei in the A = 60–100 region. More importantly, DSM has been used recently with good success in calculating nuclear structure factors, which are needed for a variety of weak interaction processes. In this article, in addition to giving an overview of this, we discuss the applications of DSM to obtain cross-sections for coherent and incoherent neutrino nucleus scattering on 96,98,100 Mo targets and also for obtaining two neutrino double beta decay nuclear transition matrix elements for 100 Mo.


Introduction
Over the last many years, the deformed shell model (DSM), based on Hartree-Fock intrinsic states with angular momentum projection and band mixing, has been established to be a good model to describe the properties of medium-mass nuclei-see the book [1] for details of DSM and for a summary of some of the applications of the model.DSM calculations are performed in the spirit of a spherical shell model, where one takes a model space and a suitable effective interaction (single-particle orbitals, single-particle energies, and a two-body effective interaction).This procedure has been found to be quite successful in the mass region A = 60-100 in describing spectroscopic properties like band crossing, shape coexistence, electro-magnetic transition probabilities, T = 0 and T = 1 bands in N = Z odd-odd nuclei with isospin projection, and so on [1].Going beyond spectroscopy, DSM has also been employed quite successfully for obtaining nuclear structure factors, which are needed for a variety of weak processes.These include (i) both two-neutrino and neutrinoless double beta decay, including positron double beta decay (all the examples studied so far are reviewed in [1]); (ii) µ → e conversion in the field of the nucleus in 72 Ge [2]; (iii) event rates for WIMP nucleus scattering off 23 Na, 40 Ar, 71,73 Ge, 75 As, 127 I, 133 Cs, and 133 Xe nuclei [3][4][5][6]; (iv) coherent and incoherent event rates for the neutrino-nucleus scattering from 28 Si, 64,70 Zn, 70,71,73,76 Ge, 75 As, 127 I, 133 Cs, and 133 Xe [4, 5,7,8].In the present paper, results are presented for neutrino-nucleus scattering cross-sections and two-neutrino double beta decay nuclear matrix elements obtained using DSM and involving Mo isotopes.Now, we give a preview.
In Section 2, we present some of the important steps in the formulation of DSM.In Section 3, we discuss its application to neutrino-nucleus scattering in 96,98,100 Mo isotopes.In Section 4, the application of the model to 2ν double beta decay in 100 Mo is presented.Finally, Section 5 gives conclusions.

Deformed Shell Model
The details of this model have been described in many of our earlier publications (see, for example, ref. [1]).In this model, for a given nucleus, starting with a model space consisting of a given set of single-particle (sp) orbitals and effective two-body Hamiltonian (TBME + spe), the lowest-energy intrinsic states are obtained by solving the axially symmetric Hartree-Fock (HF) single-particle equation self-consistently.Excited intrinsic configurations are obtained by making particle-hole excitations over the lowest intrinsic state.These intrinsic states χ K (η) do not have definite angular momenta.Hence, states of good angular momentum projected from an intrinsic state χ K (η) can be written in the form where N JK is the normalization constant, given by In Equation (1), Ω represents the Euler angles (α, β, γ) and R(Ω), which is equal to exp(−iαJ z ) exp(−iβJ y ) exp( −iγJ z ), represents the general rotation operator.The good angular momentum states projected from different intrinsic states are not, in general, orthogonal to each other.Hence, they are orthonormalized and then band-mixing calculations are performed.The resulting eigenfunctions are of the form 3. Cross-Sections of Neutral-Current Neutrino Scattering on 96,98,100 Mo Isotopes Some years ago, the neutral-current neutrino-nucleus channel was detected for the first time in the COHERENT experiment [9], more than four decades after its theoretical prediction by Freedman [10].The measurement of coherent elastic event rates of neutrino scattering on a sodium-doped CsI detector may help in answering some key questions related to neutrino properties and provide an understanding of some theories beyond the standard model of the electroweak interactions [11].Several promising nuclear detectors are on their way to being employed in designed and ongoing experiments [12][13][14][15][16][17][18][19][20][21][22].Following all these, very recently, cross-sections for neutrino scatterings off 98,100 Mo were calculated by us, using DSM for nuclear structure factors [23].In addition to briefly describing these results, here, we also report the first results for 96 Mo.
The details of the relevant formulation that provides the neutral-current ν-nucleus scattering differential cross-sections have been discussed earlier in [24][25][26].Applying a multipole analysis on the weak Hadronic current [25,26], the neutrino-nucleus scattering cross-section becomes The δ-function on the right-hand side of the above equation ensures the energy conservation, so that the excitation energy of the nucleus ω is given by E i and E f represent the energy of the initial and final states of the studied nucleus, and ϵ i and ϵ f are the incoming and outgoing energies of the neutrino.In Equation (4), the (−) sign corresponds to the scattering of the neutrinos and the (+) to the scattering of the antineutrinos.The terms C V and C A include a summation over the contributions coming from the polar-vector and axial-vector multipole operators.Similarly, the interference term C VA in Equation (4) contains the product of transverse polar-vector and transverse axial-vector matrix elements.All these involve square of the four-momentum transfer q 2 µ and the scattering angle θ.Note that the square of the four-momentum transfer q 2 µ and the magnitude of the three-momentum transfer q ≡| q | can be written in terms of the kinematic parameters, i.e., laboratory scattering angle θ and neutrino energies ϵ i = ϵ ν and ϵ f = (ϵ i − ω), as The reduced matrix elements that enter in the formulas for C V(A) and C VA are evaluated using DSM band-mixed wave function Φ J M (η), defined in Equation (3).Full details of the various operators and the method employed for the calculation of various reduced matrix elements are given in [23].
Using the above effective interaction and the single-particle energies, we first carried out an axially symmetric HF calculation by solving the HF equation self-consistently.The detailed spectroscopic results obtained using DSM for 98 Mo and 100 Mo are given in [23].The observed low-lying energy levels for both the nuclei are well-reproduced in our calculation, as are the lowest few B(E2) values.Beyond these two isotopes, here, we report briefly the results obtained for 96 Mo.For 96 Mo, the axially symmetric HF calculation is carried out using the above-described model space and effective interaction.It is found that the lowest prolate and oblate HF solutions are almost degenerate, with the oblate solution being lower by ∼0.25 MeV.We obtained excited HF configurations by making particlehole excitations over these lowest HF configurations.For positive parity, we consider 33 intrinsic configurations (17 oblate and 16 prolate) and, for negative parity, 24 intrinsic configurations (12 oblate and 12 prolate).As described before, good angular momentum states are projected from each of these intrinsic states.These states are not orthogonal to each other.Hence, they are orthonormalized and a band-mixing calculation is performed.States having similar structure and, thus, similar electromagnetic properties are classified in one band.For 96 Mo, not many collective bands have been observed.In Figure 1, we compare the calculated ground band with the experimentally observed band and we see that the agreement is reasonable.Further details regarding other spectroscopic properties of 96 Mo will be reported elsewhere.
One of our goals is to calculate the incoherent differential cross-sections for each J π state of Mo isotopes with A = 92-100.Here, we report some results for 96,98,100 Mo isotopes and they are obtained by evaluating the transition matrix elements, defining C V , C A and C VA , and connecting the ground 0 + state to various excited J π states through the use of the DSM wave function Φ J M (η) defined in Equation (3).The incoming neutrino energy is assumed to be ϵ ν = 15 MeV.We use a quenching factor of 0.35.Results for 98,100 Mo are given in our recent publication [23] and they are summarized in Figure 2. In addition, recently, we have also considered neutrino scattering off 96 Mo and details of this study will be presented in a separate publication.Here, we summarize, in Table 1, the calculated ratio of the coherent to the total (coherent + incoherent) cross-section for the three Mo isotopes.These lie, as seen from Table 1, between 85% and 91%.We note that our cross-sections correspond only to the Standard Model (SM) ν-nucleus interactions.It is worth remarking that, in other BSM electroweak processes where the coherent channel is also possible, e.g., the muon-to-electron conversion in nuclei (see e.g., refs.[29,30]), the portion of the coherent branching ratio into the total one represents also about 85-90 % of the total rate.The data are taken from ref. [23].
The contribution of the excitation to the J = 1 + state is represented by blue, to J = 2 + by red, to J = 1 − by cyan, and to J = 2 − by black.
Table 1.Coherent, incoherent, and total differential cross-sections (in units of 10 −42 MeV −1 cm 2 ) of neutrino scattering off the 96,98,100 Mo detector isotopes.The portion of the coherent into the total cross-section is also listed.Results for 98,100 Mo are taken from [23].

2ν Double Beta Decay of 100 Mo
The two-neutrino double beta decay (2νDBD) is a second-order process of weak interaction.This decay has been observed in more than 10 nuclei and is consistent with the standard model.On the other hand, neutrinoless double beta decay violates lepton number conservation and has not been observed.The successful description of the halflives (or nuclear matrix elements) for 2νββ decay provides a stringent test of the goodness of the nuclear model used, and, with confidence, we can use it to predict the neutrinoless double beta decay matrix elements.A number of nuclear models has been used to study 2νDBD, namely, spherical shell model [28,[32][33][34] , QRPA and its variants [35], interacting boson model [36], deformed shell model [1], and so on.Recently, Patel et al. [34] have employed large-scale shell model calculations to study the two-neutrino double beta decay in 82 Se, 94 Zr, 108 Cd, 124 Sn, 128 Te, 130 Te, 136 Xe, and 150 Nd using different effective interactions.Deformed shell model calculations have been performed to study the double beta decay for different nuclei in the A = 60-90 region and the results have been summarized in ref. [1].The nucleus 100 Mo is a famous double beta decay emitter and has been employed in various experiments like NEMO [37] and MOON [38] .Different variants of QRPA have been used to study ββ decay for this nucleus [39].But there are very few shell model calculations for this nucleus because of the large matrices to be diagonalized.Recently, Coraggio et al [33] have performed shell model calculation for this nucleus using effective 2νββ operators.Here, we discuss the application of DSM for two-neutrino double beta decay of 100 Mo (0ν ββ decay will be considered elsewhere).
Half-life for the 0 + i → 0 + f 2νDBD is given by g A corresponds to the axial-vector coupling constant.The kinematical factor G 2ν is independent of nuclear structure and its value G 2ν = 310.6 × 10 −20 yr −1 [40].On the other hand, the nuclear transition matrix elements (NTME) M 2ν are nuclear-model-dependent and they are given by where 0 + i , 0 + f , and |1 + m ⟩ are the initial, final, and virtual intermediate states, respec- tively, and E m are the energies of the intermediate nucleus 100 Tc.Similarly, E i and E f are the ground-state energies of the parent and daughter nuclei.We have, from ref. [41], the atomic masses of 100 Mo, 100 Tc, and 100 Ru as −86.193, −86.0215, and −89.227MeV, respectively (for nuclear mass, we need to subtract the mass of the electrons from the atomic mass).Experimentally, the lowest 1 + state in 100 Tc is the ground state.We finally obtain MeV. DSM is used to calculate the energies E 1 + of the 1 + states of 100 Tc and the corresponding wave functions.Similarly, we obtained ground-state (with 0 + ) wave functions of 100 Mo and the daughter nucleus 100 Ru.With these, the reduced matrix elements in Equation ( 8) are obtained.We have included 29 1 + states (up to 5.7 MeV) in the calculation.The calculated half life is compared with experiments in Table 1 and the calculated half life is about 20 times larger compared to the experimental value.Employing the renormalized GT operator or, equivalently, using an effective g A value, as mentioned in Table 2, will improve the calculated half life, just as in the shell model [33].

Conclusions
We have described the application of DSM to evaluate incoherent and coherent crosssections for 96,98,100 Mo.As a test of DSM, we have also reported the half life for two neutrino double beta decays of 100 Mo and compared them with experimental values.The recently observed neutral-current coherent elastic neutrino-nucleus scattering data, as well as those expected to be measured at the designed CEνNS experiments with other nuclear detectors, offer the possibility to constrain the parameter space of several models, progressing beyond the standard model, which assumes vector and scalar mediators.Our present predictions on the 96,98,100 Mo may be of help towards this direction since, as is known, these isotopes are prominent detection media.They have been employed previously in appreciably sensitive detectors (for the detection of neutrinoless double beta decay and dark matter) in the MOON and NEMO experiments.
Finally, it may be useful to perform DSM calculations for nuclear structure factors with an improved effective interactions and also by including many more deformed intrinsic states.This will be attempted in the future.

Figure 1 .
Figure 1.The ground band observed for 96 Mo is compared with the DSM predicted values.The experimental data are taken from [31].

Figure 2 .
Figure 2.The differential cross-section as a function of the excitation energy ω for 98,100 Mo at incoming neutrino energy ϵ ν = 15 MeV for different excited states.The data are taken from ref.[23].The contribution of the excitation to the J = 1 + state is represented by blue, to J = 2 + by red, to J = 1 − by cyan, and to J = 2 − by black.