Next Article in Journal
Nothing into Something and Vice Versa: A Cosmological Scenario
Next Article in Special Issue
Integral Fluxes of Neutrinos and Gamma-Rays Emitted from Neighboring X-ray Binaries
Previous Article in Journal
Analogue Gravitational Lensing in Bose-Einstein Condensates
Previous Article in Special Issue
Silicon Photomultipliers for Neutrino Telescopes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Spin-Dependent Short-Range Correlations on Nuclear Matrix Elements for Neutrinoless Double Beta Decay of 48Ca

1
Department of Physics, Indian Institute of Technology Roorkee, Roorkee 247 667, India
2
Faculty of International Studies, Osaka University of Economics and Law, Yao 581-0853, Japan
*
Author to whom correspondence should be addressed.
Universe 2023, 9(10), 444; https://doi.org/10.3390/universe9100444
Submission received: 16 August 2023 / Revised: 23 September 2023 / Accepted: 27 September 2023 / Published: 3 October 2023

Abstract

:
Neutrinoless double beta decay is a pivotal weak nuclear process that holds the potential to unveil the Majorana nature of neutrinos and predict their absolute masses. In this study, we delve into examining the impact of spin-dependent short-range correlations (SRCs) on the nuclear matrix elements (NMEs) for the light neutrino-exchange mechanism in neutrinoless double beta ( 0 ν β β ) decay of 48 Ca, employing an extensive interacting nuclear shell model. All computations are performed employing the effective shell model Hamiltonian GXPF1A, encompassing the entire f p model space through the closure approximation. Our investigation examines the NMEs’ dependencies on factors such as the number of intermediate states, coupled spin-parity attributes of neutrons and protons, neutrino momentum, inter-nucleon separation, and closure energy. This scrutiny is performed with respect to both the conventional Jastrow-type approach of SRCs, employing various parameterizations, and the spin-dependent SRC paradigm. Our findings illuminate a discernible distinction in NMEs induced by spin-dependent SRCs, differing by approximately 10–20% from those computed through the conventional Jastrow-type SRCs, incorporating distinct parameterizations.

1. Introduction

Neutrinoless double beta decay ( 0 ν β β ) is a rare and crucial weak nuclear decay that occurs in certain even–even nuclei such as 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 124 Sn, 130 Te, 126 Xe, etc. In this process, two neutrons inside the decaying nucleus are converted into two protons and two electrons without emitting any neutrinos that violate the lepton number conservation. The neutrino appears as a virtual intermediate Majorana particle [1,2,3,4,5]. Although predicted long back in 1939 by Wolfgang Furry [6,7] based on E. Majorana’s symmetric theory for fermion and anti-fermion [8] followed by G. Racah’s chain reactions [9] in 1937, this rare process has yet to be observed by experiment. However, observing this rare process would prove the existence of the Majorana neutrino, which is indeed widely favored in many beyond the standard model (BSM) physics theories for the explanation of the smallness of neutrino mass [10,11,12]. Furthermore, there is no definite value of absolute neutrino mass measured yet by experiment. An upper limit for the absolute neutrino mass of 0.8 eV has been derived through the tritium beta decay experiment-KATRIN [13] and the precise upper limit on the effective Majorana neutrino mass of 36–156 meV is predicted in 0 ν β β decay experiment KamLAND-Zen on 136 Xe [14]. In this context, the 0 ν β β process can also provide information on the absolute mass scale of neutrinos [5,15], which makes this process more interesting.
A number of particle physics mechanisms have been proposed for the 0 ν β β decay process. One of them is the widely studied standard light neutrino-exchange mechanism [16,17], which will also be the focus of this study. Other mechanisms include the heavy neutrino-exchange mechanism [7], the left-right symmetric mechanisms [18,19], the supersymmetric particle exchange mechanisms [20,21], etc.
The half-life for 0 ν β β decay is connected with absolute neutrino mass through a quantity called the nuclear matrix elements (NMEs), which are theoretically calculated in different nuclear many-body models [4] such as quasiparticle random phase approximation (QRPA) [7], the interacting shell-model (ISM) [22,23,24,25,26], the interacting boson model (IBM) [27,28], the generator coordinate method (GCM) [29,30,31], energy density functional (EDF) theory [29,30], relativistic energy density functional (REDF) theory [30,31], and the projected Hartree–Fock–Bogolibov model (PHFB) [32], the ab initio variational Monte Carlo (VMC) technique [33,34,35], etc. The accuracy of NME calculations plays an important role in predicting the correct value for the half-life of 0 ν β β decay and absolute neutrino mass. Hence, calculating the NMEs accurately is one of the central themes of research for studying the 0 ν β β decay.
Around the world, several major experimental investigations are being carried out for different 0 ν β β decaying candidates such as GERDA ( 76 Ge), MAJORANA DEMONSTRATOR ( 76 Ge), LEGEND ( 76 Ge), CUORE ( 130 Te), CUPID ( 82 Se/ 100 Mo/ 130 Te), AMoRE ( 100 Mo), EXO-200 ( 136 Xe), nEXO ( 136 Xe), NEXT ( 136 Xe), PandaX-III ( 136 Xe), KamLAND-Zen and KamLAND2-Zen ( 136 Xe), CANDLES ( 48 Ca), and Super-NEMO ( 82 Se) [1]. Our particular interest in this study is the 0 ν β β decay of 48 Ca, which is one of the simplest decay candidates with significant experimental interest. In particular, we aim to examine the effect of spin-dependent SRCs on the NMEs for the 0 ν β β decay of 48 Ca using the nuclear shell model.
0 ν β β decay in 48 Ca proceeds as:
48 Ca 48 Ti + e + e .
The computation of NMEs corresponding to the light neutrino-exchange mechanism of the 0 ν β β decay in 48 Ca was carried out in several studies ([23,24,25,36,37,38,39,40,41,42]). In these investigations, the influences of SRCs were integrated through a Jastrow-type approach [43,44], characterizing the radial-dependence effects of SRCs. An alternative method, known as the Unitary Correlation Operator Method (UCOM) [45,46,47], also accounted for SRCs via radial dependence. More recently, in a notable advancement [48], the interplay of spin in conjunction with SRCs was incorporated to compute NMEs for the light neutrino-exchange 0 ν β β decay in 48 Ca. This study, however, operated within a confined valence space, focusing exclusively on the f 7 / 2 orbital within the p f -shell. The main findings in Ref. [48] emphasize that accounting for spin-dependent SRCs leads to an approximate 20% reduction in the NMEs when compared to the no-SRC scenario within the pure shell model framework.
Motivated by these advances, our current endeavor aims to expand the computational scope. Specifically, we engage the complete f p model space encompassing the orbitals 0 f 7 / 2 , 0 f 5 / 2 , 1 p 3 / 2 , and 1 p 1 / 2 for 48 Ca. This comprehensive approach enables us to thoroughly assess the cumulative effect of spin on SRCs for the 0 ν β β decay in 48 Ca. Furthermore, we undertake an in-depth exploration of the variances between the spin-dependent SRC scenario and the conventional Jastrow-type SRC approach. We systematically investigate the dependencies of NMEs on factors such as the coupled spin-parity of protons and neutrons, the number of intermediate states, neutrino momentum, inter-nucleon separation, and closure energy. Our goal is to uncover how these dependencies diverge under the influence of spin-dependent SRCs and to contrast them with the outcomes of the traditional Jastrow-type SRC approach. As our investigation progresses, we envision extending this approach to encompass other decay candidates within the 0 ν β β realm, thereby offering a more comprehensive understanding of this intriguing nuclear process.
Driven by this motivation, the rest of this paper is structured as follows. Section 2 outlines the formulation of the decay rate, NMEs, and transition operators pertinent to 0 ν β β decay. In Section 3, we delve into the integration of SRC effects within both the Jastrow-type approach and the spin-dependent approach. The methodology for calculating NMEs using the shell model within the closure method is detailed in Section 4. Subsequently, Section 5 presents the comprehensive results and initiates discussions around the observed outcomes. Lastly, Section 6 encapsulates the synthesis of our findings and delineates potential avenues for future exploration.

2. The Decay Rate and NMEs for 0 ν β β Decay

The decay rate for the light neutrino-exchange mechanism of 0 ν β β decay can be written as [17]
[ T 1 2 0 ν ] 1 = G 0 ν | M 0 ν | 2 m β β m e 2 ,
where G 0 ν is the phase-space factor that can be calculated accurately [49], M 0 ν is the total nuclear matrix element for the light neutrino-exchange mechanism, and m β β is the effective Majorana neutrino mass defined by the neutrino mass eigenvalues m k and the neutrino mixing matrix elements U e k
m β β = k m k U e k 2 .
The total nuclear matrix element M 0 ν is the sum of Gamow–Teller ( M G T ), Fermi ( M F ), and tensor ( M T ) matrix elements, as given by [5]:
M 0 ν = M G T g V g A 2 M F + M T ,
where g V and g A are the vector and axial-vector constants, respectively. In this study, g V = 1 and the bare value of g A = 1.27 is used. The matrix elements M G T , M F , and M T of transition operator O 12 α of 0 ν β β decay are expressed as [25]:
M α = f | O 12 α | i
where α F , G T , T , | i corresponds to the 0 + ground state of the parent nucleus 48 Ca in the present study, and | f corresponds to the 0 + ground state of the granddaughter nucleus 48 Ti.
The calculation of two-body matrix elements (TBMEs) for 0 ν β β decay scalar two-particle transition operators O 12 α have both spin and radial neutrino potential parts. These operators are given by [24]:
O 12 G T = τ 1 τ 2 ( σ 1 · σ 2 ) H G T ( r , E k ) , O 12 F = τ 1 τ 2 H F ( r , E k ) , O 12 T = τ 1 τ 2 S 12 H T ( r , E k ) ,
where, τ is the isospin annihilation operator, r = r 1 r 2 is the inter-nucleon distance of the decaying nucleons, and r = | r | . The operator S 12 is defined as S 12 = 3 ( σ 1 · r ^ ) ( σ 2 · r ^ ) ( σ 1 · σ 2 ) . For the light neutrino-exchange mechanism of 0 ν β β decay, the radial neutrino potential with explicit dependence on the energy of the intermediate states is given by [24]:
H α ( r , E k ) = 2 R π 0 f α ( q , r ) q d q q + E k ( E i + E f ) / 2 ,
where R is the radius of the parent nucleus, q is the momentum of the virtual Majorana neutrino, E i , E k , and E f are the energies of the initial, intermediate, and final nuclei, f α ( q , r ) = j p ( q , r ) h α ( q , r ) with j p ( q , r ) is a spherical Bessel function with p = 0 for Fermi and Gamow–Teller type NMEs and h α ( q , r ) is the term that accounts for the effects of higher-order currents (HOC) and finite nucleon size (FNS) [17] given by [17,43]:
h F ( q 2 ) = g V 2 ( q 2 ) , h G T ( q 2 ) = g A 2 ( q 2 ) g A 2 1 2 3 q 2 q 2 + m π 2 + 1 3 q 2 q 2 + m π 2 2
+ 2 3 g M 2 ( q 2 ) g A 2 q 2 4 m p 2 , h T ( q 2 ) = g A ( q 2 ) g A 2 2 3 q 2 q 2 + m π 2 1 3 q 2 q 2 + m π 2 2
+ 1 3 g M 2 ( q 2 ) g A 2 q 2 4 m p 2 .
In this regard, the g V ( q 2 ) , g A ( q 2 ) , and g M ( q 2 ) form factors, which account for FNS effects, are used. In the dipole approximation, the form factors are given by the following Equations [17,23]:
g V ( q 2 ) = g V 1 + q 2 M V 2 2 ,
g A ( q 2 ) = g A 1 + q 2 M A 2 2 ,
g M ( q 2 ) = ( μ p μ n ) g V ( q 2 ) .
where g V , g A , μ p , and μ n are the vector, axial-vector, and magnetic moment coupling constants for the nucleon, and M V and M A are the vector and axial-vector meson masses, respectively. The values of M V and M A are 850 MeV and 1086 MeV, respectively, while a μ p μ n value of 4.7 is used in the calculation [23]. The masses of the proton and pion are denoted by m p and m π , respectively.

3. The Effects of Short-Range Correlations on 0 ν β β Decay

In the computation of the NMEs for 0 ν β β decay, it becomes imperative to factor in the influences of SRCs, as discussed previously. Over time, various techniques have been employed to integrate SRC effects into NME calculations for 0 ν β β decay. One of the extensively utilized and classical methodologies is known as the Jastrow approach, as discussed in references [43,44]. More recently, a novel spin-dependent approach has emerged, as outlined in reference [48]. Both of these methodologies are delineated comprehensively below.

3.1. The Jastrow Approach

A standard method to include SRCs is via a phenomenological Jastrow-like function [43,44]. By including the SRC effect in the Jastrow approach, one can write the NMEs of 0 ν β β defined in Equation (5) as [43]
M α = f | f J a s t r o w ( r ) O 12 α f J a s t r o w ( r ) | i ,
where the Jastrow-type SRC function is defined as
f J a s t r o w ( r ) = 1 c e a r 2 ( 1 b r 2 ) .
In the literature, three different SRC parameterizations are used: Miller–Spencer, Charge-Dependent-Bonn (CD-Bonn), and Argonne V18 (AV18) to parameterize a , b , and c[23]. The parameters a, b, and c in different SRC parameterizations are given in Table 1. This approach of using a Jastrow-like function to include the effects of SRCs is extensively used in Refs. [23,41,42].

3.2. Spin-Dependent Approach

Recently [48], for the first time, the short nature of nucleon–nucleon correlations was accounted for based on the dependence of the SRC operator in spin and coordinate space. The new spin-dependent SRC correlation function is defined as
f S D ( r ) = f ( r ) + g ( r ) σ 1 . σ 2 .
The transition operators for 0 ν β β decay incorporating spin-dependent SRCs are written as
O ˜ 12 α = f S D ( r ) O 12 α f S D ( r ) .
Using the relation ( σ 1 . σ 2 ) 2 = ( 3 2 ( σ 1 . σ 2 ) ) , the new Fermi operator of 0 ν β β decay is written as
O ˜ 12 F = [ f ( r ) 2 + 3 g ( r ) 2 ] H F ( r ) + ( σ 1 . σ 2 ) 2 g ( r ) [ f ( r ) g ( r ) ] H F ( r )
and the Gamow–Teller operator of 0 ν β β decay is written as
O ˜ 12 G T = [ f ( r ) 2 4 g ( r ) f ( r ) + 7 g ( r ) 2 ] H G T ( r ) ( σ 1 . σ 2 ) + 6 g ( r ) [ f ( r ) g ( r ) ] H G T ( r )
The functions f ( r ) and g ( r ) are defined as [48,50,51]
f ( r ) = a 1 b 1 e c 1 r 2 + d 1 e e 1 ( r f 1 ) 2 ,
g ( r ) = a 2 e b 2 r 2 ( 1 + c 2 r + d 2 r 2 ) .
The different parameters of f ( r ) and g ( r ) are presented in Table 2. These parameters are good for full f p and other configuration spaces. The authors of [48] may have considered only the 0 f 7 / 2 orbital as they wanted only to acquire initial insights through the pure shell model. Explicit details of using this SRC approach can be found in Ref. [50]. The derivation of the above formalism can be found in Refs. [50,51].
In addition to the above, the authors of [52,53] have recently proposed another method, namely, the Unitary Correlation Operator Method [45,46,47], to estimate the effects of SRCs, which leads to a much smoother correction [44]. In the present study, we only focus on using a standard Jastrow-type approach and a spin-dependent approach to estimate the effects of SRCs. Detailed descriptions of incorporating the SRC effects in different approaches can be found in Refs. [43,44].
It is worthwhile to mention that in a recent work [54], the short-range (SR) matrix element, a leading-order factor in light neutrino-exchange 0 ν β β decay, as established in a previous notable study [55], was meticulously computed. This calculation employed the generalized contact formalism (GCF) in conjunction with the shell model. The resulting total NME is composed of a long-range term, denoted as M L 0 ν = M G T + M F + M T , and a short-range term, M S 0 ν . The investigation revealed a substantial contribution of the SR term to the total NME. Specifically, for nuclei with A = 48, M S 0 ν accounted for approximately 35–60% of M L 0 ν . Moreover, it was observed that the SR term amplifies the total NME by 25–40% in heavier nuclei, indicating a significant impact on the half-life of 0 ν β β decay, while an intriguing avenue of inquiry lies in assessing how spin-dependent SRC effects might influence the SR matrix element, it should be noted that this current manuscript does not delve into that particular investigation, as it falls outside the scope of the present study. For those interested in a more comprehensive understanding of SR matrix elements, detailed information can be found in Refs. [54,55,56].

4. The Closure Method of NME Calculation for 0 ν β β Decay through the ( n 2 )  Channel

Now we discuss, including all the effects, how NMEs are calculated in the nuclear shell model. First, there are two approaches, nonclosure and closure. The nonclosure approach [24,36,37] considers the effects of excitation energies of a large number of intermediate states of the true virtual nucleus ( 48 Sc in the present case) explicitly through the denominator of Equation (7). In closure approximation [23], one approximates the term E k ( E i + E f ) / 2 in the denominator of the neutrino potential of Equation (7) with a constant closure energy ( E ) value such that
[ E k ( E i + E f ) / 2 ] E .
Closure approximation avoids the complexity of calculating a large number of intermediate states, which can be computationally challenging for nuclear shell models for higher mass isotopes. The difficult part of closure approximation is picking the correct closure energy, which has no definite method yet and can greatly influence the accuracy of the calculated NMEs. Recently, there have been many studies of the nonclosure approach for 0 ν β β decay of 48 Ca [24,36,37], which predicted that E = 0.5 MeV is an optimal value of closure energy that gives values of the NMEs in the closure method very close to those in the nonclosure method. Hence, we use the closure method of NME calculation with a closure energy value E = 0.5 MeV.
Based on the closure approach, the M G T , M F , and M T matrix elements of the scalar two-body transition operator O 12 α of 0 ν β β can be expressed as the sum over the product of the two-body transition density (TBTD) and anti-symmetric two-body matrix elements ( k 1 , k 2 , J T | O 12 α | k 1 , k 2 , J T A ) . For the first, we write the partial NME as a function of the coupled spin-parity ( J π ) of protons and neutrons and a number of intermediate states ( N m ) of the TBTD calculations as [25]
M α ( J π , N m ) = k 1 k 2 , k 1 k 2 TBTD ( f , m , i , J π ) k 1 , k 2 , J π T | O 12 α | k 1 , k 2 , J π T A ,
such that the final required Fermi, GT, and tensor-type NME of Equation (5) can be written as
M α = J , N m N c M α ( J π , N m ) ,
where, α = ( F , G T , T ) , J π is the coupled spin-parity of two decaying neutrons or two final created protons, A denotes that the two-body matrix elements are obtained using anti-symmetric two-nucleon wave functions, k stands for the set of spherical quantum numbers ( n ; l ; j ) , N c is cutoff on number of states of intermediate nucleus ( 46 Ca in the present case) for TBTD calculations. In our case, | i is 0 + g.s. of the parent nucleus 48 Ca, | f is the 0 + g.s. of the granddaughter nucleus 48 Ti, and k has the spherical quantum numbers for 0 f 7 / 2 , 0 f 5 / 2 , 1 p 3 / 2 , and 1 p 1 / 2 orbitals. The TBTD can be expressed as [25]
TBTD ( f , m , i , J ) = f | | [ A + ( k 1 , k 2 , J ) A ˜ ( k 1 , k 2 , J ) ] ( 0 ) | | i ,
where,
A + ( k 1 , k 2 , J ) = [ a + ( k 1 ) a + ( k 2 ) ] M J 1 + δ k 1 k 2 ,
and
A ˜ ( k 1 , k 2 , J ) = ( 1 ) J M A + ( k 1 , k 2 , J , M )
are the two-particle creation and annihilation operators of rank J, respectively.
To evaluate TBTD, one needs a large number of two-nucleon transfer amplitudes (TNAs). The TNAs are calculated with a large set of intermediate states | m of the (n-2) nucleon system ( 46 Ca in the present study), where n is the number of nucleons for the parent nucleus. The TBTD in terms of TNAs is expressed as [25]
TBTD ( f , m , i , J ) = m TNA ( f , m , k 1 , k 2 , J m ) TNA ( i , m , k 1 , k 2 , J m ) ,
where the TNAs are given by
TNA ( f , m , k 1 , k 2 , J m ) = f | | A + ( k 1 , k 2 , J ) | | m 2 J 0 + 1 .
Here, J m is the spin of the allowed states of 46 Ca, J 0 is spin of | i and | f , and J m = J when J 0 = 0 [25].

5. Results and Discussion

The TNAs are evaluated utilizing the nuclear shell model. For this purpose, we employ the KSHELL code [57], which facilitates the computation of essential wave functions and energies concerning the initial, intermediate, and final nuclei involved in the 0 ν β β decay process of 48 Ca. Subsequently, these calculated wave functions are harnessed to determine the TNAs, a pivotal component within the expression of NMEs for the 0 ν β β decay.
For these calculations, we utilize the shell model Hamiltonian GXPF1A [58], tailored for the f p model space, serving as a foundational input for our computations.
The f p model space comprises the f 7 / 2 , p 3 / 2 , f 5 / 2 , and p 1 / 2 orbitals. The computation of the two-body matrix element (TBME) components within the NMEs is executed using a proprietary code developed by us. The outcomes of distinct NME types, computed through the shell model approach with the Jastrow methodology employing various SRC parameterizations, and spin-dependent SRCs are presented in Table 3.
It is important to note that the presentation of tensor-type NMEs is omitted from Table 3 due to the lack of a solvable expression for the TBMEs inclusive of the spin-dependent SRC function specified in Equation (16). This is not a critical concern as the Gamow–Teller (GT) and Fermi-type NMEs overwhelmingly dominate over tensor-type NMEs. Therefore, in Table 3, the calculation of total NMEs solely considers the Fermi and GT-type contributions.
The NMEs corresponding to the application of spin-dependent SRCs are approximately 10% diminished in comparison to those calculated with AV18-type NMEs, as well as SRC-none scenarios, across Fermi, Gamow–Teller (GT), and total NMEs. In contrast, when gauged against CD-Bonn-type, the NMEs with spin-dependent SRCs experience a more pronounced reduction, approximately 20%, across Fermi, GT, and total NMEs.
Remarkably, among the various SRC parameterizations, the Miller–Spencer-type NMEs manifest the most pronounced impact of SRC. Notably, NMEs with spin-dependent SRC showcase an expansion of approximately 20% relative to Miller–Spencer type NMEs, encompassing the realms of Fermi, GT, and total NMEs.
As mentioned above, the main findings of Ref. [48] is that accounting for spin-dependent SRsC leads to an approximate 20% reduction in NME values when compared to the no-SRC scenario within the pure shell model framework. In our analysis, we observe a comparable 10% decrease in NME values under the same conditions relative to the no-SRC case, but in the interacting shell model. These differences may be attributed to our consideration of the entire f p model space, encompassing 0 f 7 / 2 , 0 f 5 / 2 , 1 p 3 / 2 , and 1 p 1 / 2 orbitals in contrast to the approach in Ref. [48], which employed a simplified shell model with a single active orbital, 0 f 7 / 2 .

5.1. Dependence of NMEs on Coupled Spin-Parity ( J π ) of Two Decaying Neutrons and Two Created Protons

Next, we examine the comparison of the contributions from different coupled spin-parities in the NMEs for CD-Bonn-type NMEs and spin-dependent SRCs. To do this, we decompose NMEs in terms of partial nuclear matrix elements of coupled spin-parity ( J π ) of two decaying neutrons or two created protons as
M α ( J π ) = N m N c M α ( J π , N m ) ,
where one can define M α 0 ν ( J π ) using Equation (23). It should be noted that j π in the above equation also represents the spin-parity of the states of the intermediate nucleus 46 Ca.
The contributions of NMEs through different J π are shown in Figure 1. We observe that for all types of NMEs, the most dominant contributions come from the 0 + and 2 + states. Additionally, the contribution from the 0 + and 2 + states have opposite signs, leading to a reduction in the total NMEs. There are also small contributions from the 4 + and 6 + states, with almost negligible contributions from odd- J π states. This is due to the pairing effect, which is responsible for the dominance of even- J π contributions [25]. Note that spin-dependent SRCs are smaller in magnitude than CD-Bonn-type SRC parameterizations. However, the overall pattern of dependence on J π is very similar for both SRCs.

5.2. Variation in NMEs for 0 ν β β with the Cutoff Number of States ( N c ) of 46 Ca

To assess the impact of the number of states on the calculated NMEs, we examine the dependence of the NMEs on the cutoff number of states ( N c ) for each allowed J m π of 46 Ca, which acts as an intermediate nucleus for TNA calculations. Shell model code packages like KSHELL do not provide options for setting a cutoff for the excitation energy, which may be more relevant, but they do provide options to set a cutoff for the number of states for each spin-parity of the intermediate states of 46 Ca. Hence, we consider a certain number of intermediate states for each spin-parity of 46 Ca with a cutoff according to the computational power available for the calculations. We express the NMEs as a function of N c in the closure method as
M α = J , N m N c M α ( J π , N m )
where M α ( J π , N m ) is the same as that defined in Equation (23).
The dependence of the different types of NMEs on N c is shown in Figure 2. We compare the results for CD-Bonn-type SRC parameterization with spin-dependent SRCs. We find that the first few low-lying states contribute constructively and destructively, but after N c = 30 , the different types of NMEs reach a stable value. For larger N c , the NMEs become mostly constant. To obtain NMEs with negligible uncertainty, we were able to consider N c = 100 for each allowed J m π of 46 Ca. We note that a similar dependence of NMEs on N c is seen for other SRC parameterizations. The overall pattern of dependence of NMEs on N c for spin-dependent SRCs and CD-Bonn-type SRCs is similar, but spin-dependent SRC NMEs are much smaller in magnitude compared to NMEs for CD-Bonn-type SRCs.

5.3. Neutrino Momentum (q) and Radial (r) Distribution of NMEs

We have also examined the neutrino momentum transfer (q) distribution of NMEs. One can define q dependent distribution ( C α ( q ) ) such that the NME defined in Equation (5) can be written as [59]
M α = 0 C α ( q ) d q .
The distribution of different C α ( q ) is shown in Figure 3. Here, NMEs are compared for CD-Bonn-type SRC parameterization and spin-dependent SRCs. It is found that most of the contributions of the NMEs come from q below 500 MeV. The peak contributions come from q around 10 MeV and 180 MeV with a peak value around 0.003 MeV 1 for GT-type NMEs with CD-Bonn-type SRCs and a peak value around 0.0027 MeV 1 for GT-type NMEs with spin-dependent SRCs. The peak value for Fermi-type NMEs is around −0.001 MeV 1 for both CD-Bonn-type and spin-dependent SRCs. The overall pattern of the dependence of NMEs for both CD-Bonn-type and spin-dependent SRCs is similar.
The radial distribution of NMEs is also explored. One can write the radial-dependent NME distribution ( C α ( r ) ) [42,59,60] such that
M α = 0 C α ( r ) d r .
The distribution of different C α ( r ) is shown in Figure 4. Here, NMEs are calculated using the closure method for CD-Bonn-type and spin-dependent SRCs with E = 0.5 MeV. It is found that most of the contribution comes from r less than 4 fm. The NMEs peak around 1 fm for both CD-Bonn-type and spin-dependent SRCs. The maximum value for C α ( r ) is near 1.2 fm 1 For GT-type and around −0.5 fm 1 for Fermi-type NMEs for both CD-Bonn-type and spin-dependent SRCs. Here also, the overall pattern of variation for CD-Bonn-type and spin-dependent SRCs is similar, but they differ slightly in magnitude.

5.4. Variation in NMEs with Closure Energy

Finally, we also show the variation in different types of total NMEs with closure energy E in Figure 5. The NMEs shown here are the total NMEs in different parameterizations of the Jastrow-type approach and spin-dependent approach. In total NMEs, the Fermi- and GT-type NMEs are only considered, and tensor-type NME is not included for the reason mentioned above. In changing E = 0 to 15 MeV, there are about 10–15% decrements of total NMEs in the closure method. Other than varying in magnitude, the spin-dependent SRCs vary very similarly to other SRC parameterizations of the Jastrow-type approach.

6. Summary and Conclusions

In summary, our investigation has centered on comprehensively exploring the effects of spin-dependent SRCs on the NMEs associated with the light neutrino-exchange 0 ν β β decay of 48 Ca. This exploration was juxtaposed with the traditional Jastrow-type approach to SRCs, employing diverse parameterizations. All computations were performed within the framework of a nuclear shell model, incorporating the closure approximation and utilizing the GXPF1A Hamiltonian to encapsulate the complete f p model space. The nuclear wave functions were derived through shell model diagonalization, leveraging the KSHELL code. These derived wave functions formed the bedrock for evaluating a multitude of TNAs, integral to the formulation of the NMEs. The segment concerning TBMEs within the NMEs was computed through code developed by us.
Beyond the comparison of the NMEs of different types (Fermi, Gamow, and total) originating from spin-dependent SRCs and Jastrow-type SRCs, our analysis delved into the comparative dependencies of NMEs linked to attributes such as the coupled spin-parity of protons and neutrons, the number of intermediate states of TNA calculation, neutrino momentum, radial separation, and closure energy.
Our findings prominently highlight a discernible reduction in NMEs resulting from the application of spin-dependent SRCs, exhibiting, for instance, a 10–20% decrease in NMEs when contrasted with Jastrow-type SRCs with diverse parameterizations. This substantial variation may bear noticeable implications for half-life and neutrino mass considerations. At the same time, while exhibiting a different magnitude, the general pattern of variation in NMEs under the influence of spin-dependent SRCs still bears a resemblance to that observed with the Jastrow-type approach.
In our forthcoming research, our aim is to extend the influence of the spin-dependent SRC approach to encompass other candidates for 0 ν β β decay, exploring a spectrum of mechanisms and scenarios.

Author Contributions

Validation: S.S. and Y.I. Original Draft Preparation: S.S. and Y.I. Writing—Review and Editing: S.S. and Y.I. Both authors have made equal contributions throughout multiple phases of calculations and the formulation of the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

S.S. is thankful to the Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Government of India, for the award of the SERB-National Postdoctoral Fellowship under Grant No. PDF/2022/003729.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

S.S. expresses profound gratitude to Rajdeep Chatterjee for the invaluable mentorship received throughout the tenure as a SERB-National Postdoctoral Fellow at the Indian Institute of Technology Roorkee. S.S. also acknowledges the National Supercomputing Mission (NSM) for providing computing resources of ‘PARAM Ganga’ at the Indian Institute of Technology Roorkee, which is implemented by C-DAC and supported by the Ministry of Electronics and Information Technology (MeitY) and Department of Science and Technology (DST), Government of India.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agostini, M.; Benato, G.; Detwiler, J.A.; Menéndez, J.; Vissani, F. Toward the discovery of matter creation with neutrinoless ββ decay. Rev. Mod. Phys. 2023, 95, 025002. [Google Scholar] [CrossRef]
  2. Dolinski, M.J.; Poon, A.W.P.; Rodejohann, W. Neutrinoless Double-Beta Decay: Status and Prospects. Ann. Rev. Nucl. Part. Sci. 2019, 69, 219–251. [Google Scholar] [CrossRef]
  3. Vergados, J.D.; Ejiri, H.; Šimkovic, F. Neutrinoless double beta decay and neutrino mass. Int. J. Mod. Phys. E 2016, 25, 1630007. [Google Scholar] [CrossRef]
  4. Engel, J.; Menéndez, J. Status and future of nuclear matrix elements for neutrinoless double-beta decay: A review. Rep. Prog. Phys. 2017, 80, 046301. [Google Scholar] [CrossRef]
  5. Avignone, F.T.; Elliott, S.R.; Engel, J. Double beta decay, Majorana neutrinos, and neutrino mass. Rev. Mod. Phys. 2008, 80, 481–516. [Google Scholar] [CrossRef]
  6. Furry, W.H. On Transition Probabilities in Double Beta-Disintegration. Phys. Rev. 1939, 56, 1184–1193. [Google Scholar] [CrossRef]
  7. Vergados, J.D.; Ejiri, H.; Šimkovic, F. Theory of neutrinoless double-beta decay. Rep. Prog. Phys. 2012, 75, 106301. [Google Scholar] [CrossRef] [PubMed]
  8. Majorana, E. Symmetric theory of electron and positron. Il Nuovo Cimento (1924–1942) 1937, 14, 171. [Google Scholar] [CrossRef]
  9. Racah, G. Sulla simmetria tra particelle e antiparticelle. Il Nuovo Cimento 1937, 14, 322. [Google Scholar] [CrossRef]
  10. Deppisch, F.F.; Hirsch, M.; Päs, H. Neutrinoless double-beta decay and physics beyond the standard model. J. Phys. G Nucl. Part. Phys. 2012, 39, 124007. [Google Scholar] [CrossRef]
  11. Schechter, J.; Valle, J.W. Neutrinoless double-β decay in SU (2)× U (1) theories. Phys. Rev. D 1982, 25, 2951. [Google Scholar] [CrossRef]
  12. Rodejohann, W. Neutrino-less double beta decay and particle physics. Int. J. Mod. Phys. E 2011, 20, 1833–1930. [Google Scholar] [CrossRef]
  13. Aker, M.; Beglarian, A.; Behrens, J.; Berlev, A.; Besserer, U.; Bieringer, B.; Block, F.; Bobien, S.; Bottcher, M.; Bornschein, B.; et al. Direct neutrino-mass measurement with sub-electronvolt sensitivity. Nat. Phys. 2022, 18, 160–166. [Google Scholar] [CrossRef]
  14. KamLAND-Zen Collaboration. Search for the Majorana Nature of Neutrinos in the Inverted Mass Ordering Region with KamLAND-Zen. Phys. Rev. Lett. 2023, 130, 051801. [Google Scholar] [CrossRef]
  15. Tomoda, T. Double beta decay. Rep. Prog. Phys. 1991, 54, 53. [Google Scholar] [CrossRef]
  16. Rodin, V.; Faessler, A.; Šimkovic, F.; Vogel, P. Assessment of uncertainties in QRPA 0νββ-decay nuclear matrix elements. Nucl. Phys. A 2006, 766, 107–131. [Google Scholar] [CrossRef]
  17. Šimkovic, F.; Pantis, G.; Vergados, J.D.; Faessler, A. Additional nucleon current contributions to neutrinoless double β decay. Phys. Rev. C 1999, 60, 055502. [Google Scholar] [CrossRef]
  18. Mohapatra, R.N.; Senjanović, G. Neutrino Mass and Spontaneous Parity Nonconservation. Phys. Rev. Lett. 1980, 44, 912–915. [Google Scholar] [CrossRef]
  19. Mohapatra, R.N.; Vergados, J.D. New Contribution to Neutrinoless Double Beta Decay in Gauge Models. Phys. Rev. Lett. 1981, 47, 1713–1716. [Google Scholar] [CrossRef]
  20. Mohapatra, R.N. New contributions to neutrinoless double-beta decay in supersymmetric theories. Phys. Rev. D 1986, 34, 3457–3461. [Google Scholar] [CrossRef]
  21. Vergados, J. Neutrinoless double β-decay without Majorana neutrinos in supersymmetric theories. Phys. Lett. B 1987, 184, 55–62. [Google Scholar] [CrossRef]
  22. Caurier, E.; Menéndez, J.; Nowacki, F.; Poves, A. Influence of Pairing on the Nuclear Matrix Elements of the Neutrinoless ββ Decays. Phys. Rev. Lett. 2008, 100, 052503. [Google Scholar] [CrossRef]
  23. Horoi, M.; Stoica, S. Shell model analysis of the neutrinoless double-β decay of 48Ca. Phys. Rev. C 2010, 81, 024321. [Google Scholar] [CrossRef]
  24. Sen’kov, R.A.; Horoi, M. Neutrinoless double-β decay of 48Ca in the shell model: Closure versus nonclosure approximation. Phys. Rev. C 2013, 88, 064312. [Google Scholar] [CrossRef]
  25. Brown, B.A.; Horoi, M.; Sen’kov, R.A. Nuclear Structure Aspects of Neutrinoless Double-β Decay. Phys. Rev. Lett. 2014, 113, 262501. [Google Scholar] [CrossRef]
  26. Iwata, Y.; Shimizu, N.; Otsuka, T.; Utsuno, Y.; Menéndez, J.; Honma, M.; Abe, T. Large-Scale Shell-Model Analysis of the Neutrinoless ββ Decay of 48Ca. Phys. Rev. Lett. 2016, 116, 112502. [Google Scholar] [CrossRef] [PubMed]
  27. Barea, J.; Iachello, F. Neutrinoless double-β decay in the microscopic interacting boson model. Phys. Rev. C 2009, 79, 044301. [Google Scholar] [CrossRef]
  28. Barea, J.; Kotila, J.; Iachello, F. Limits on Neutrino Masses from Neutrinoless Double-β Decay. Phys. Rev. Lett. 2012, 109, 042501. [Google Scholar] [CrossRef]
  29. Rodríguez, T.R.; Martínez-Pinedo, G. Energy Density Functional Study of Nuclear Matrix Elements for Neutrinoless ββ Decay. Phys. Rev. Lett. 2010, 105, 252503. [Google Scholar] [CrossRef]
  30. Song, L.S.; Yao, J.M.; Ring, P.; Meng, J. Relativistic description of nuclear matrix elements in neutrinoless double-β decay. Phys. Rev. C 2014, 90, 054309. [Google Scholar] [CrossRef]
  31. Yao, J.M.; Song, L.S.; Hagino, K.; Ring, P.; Meng, J. Systematic study of nuclear matrix elements in neutrinoless double-β decay with a beyond-mean-field covariant density functional theory. Phys. Rev. C 2015, 91, 024316. [Google Scholar] [CrossRef]
  32. Rath, P.K.; Chandra, R.; Chaturvedi, K.; Raina, P.K.; Hirsch, J.G. Uncertainties in nuclear transition matrix elements for neutrinoless ββ decay within the projected-Hartree-Fock-Bogoliubov model. Phys. Rev. C 2010, 82, 064310. [Google Scholar] [CrossRef]
  33. Pastore, S.; Carlson, J.; Cirigliano, V.; Dekens, W.; Mereghetti, E.; Wiringa, R.B. Neutrinoless double-β decay matrix elements in light nuclei. Phys. Rev. C 2018, 97, 014606. [Google Scholar] [CrossRef]
  34. Wang, X.B.; Hayes, A.; Carlson, J.; Dong, G.; Mereghetti, E.; Pastore, S.; Wiringa, R.B. Comparison between variational Monte Carlo and shell model calculations of neutrinoless double beta decay matrix elements in light nuclei. Phys. Lett. B 2019, 798, 134974. [Google Scholar] [CrossRef]
  35. Cirigliano, V.; Dekens, W.; de Vries, J.; Graesser, M.L.; Mereghetti, E.; Pastore, S.; Piarulli, M.; van Kolck, U.; Wiringa, R.B. Renormalized approach to neutrinoless double-β decay. Phys. Rev. C 2019, 100, 055504. [Google Scholar] [CrossRef]
  36. Sarkar, S.; Iwata, Y.; Raina, P.K. Nuclear matrix elements for the λ mechanism of 0νββ decay of 48Ca in the nuclear shell-model: Closure versus nonclosure approach. Phys. Rev. C 2020, 102, 034317. [Google Scholar] [CrossRef]
  37. Sarkar, S.; Kumar, P.; Jha, K.; Raina, P.K. Sensitivity of nuclear matrix elements of 0νββ of 48Ca to different components of the two-nucleon interaction. Phys. Rev. C 2020, 101, 014307. [Google Scholar] [CrossRef]
  38. Ahmed, F.; Horoi, M. Interference effects for 0νββ decay in the left-right symmetric model. Phys. Rev. C 2020, 101, 035504. [Google Scholar] [CrossRef]
  39. Horoi, M.; Neacsu, A. Shell model study of using an effective field theory for disentangling several contributions to neutrinoless double-β decay. Phys. Rev. C 2018, 98, 035502. [Google Scholar] [CrossRef]
  40. Horoi, M. Shell model analysis of competing contributions to the double-β decay of 48Ca. Phys. Rev. C 2013, 87, 014320. [Google Scholar] [CrossRef]
  41. Neacsu, A.; Stoica, S.; Horoi, M. Fast, efficient calculations of the two-body matrix elements of the transition operators for neutrinoless double-β decay. Phys. Rev. C 2012, 86, 067304. [Google Scholar] [CrossRef]
  42. Menendez, J.; Poves, A.; Caurier, E.; Nowacki, F. Disassembling the Nuclear Matrix Elements of the Neutrinoless beta beta Decay. Nucl. Phys. A 2009, 818, 139–151. [Google Scholar] [CrossRef]
  43. Šimkovic, F.; Faessler, A.; Müther, H.; Rodin, V.; Stauf, M. 0νββ-decay nuclear matrix elements with self-consistent short-range correlations. Phys. Rev. C 2009, 79, 055501. [Google Scholar] [CrossRef]
  44. Vogel, P. Nuclear structure and double beta decay. J. Phys. G Nucl. Part. Phys. 2012, 39, 124002. [Google Scholar] [CrossRef]
  45. Feldmeier, H.; Neff, T.; Roth, R.; Schnack, J. A unitary correlation operator method. Nucl. Phys. A 1998, 632, 61–95. [Google Scholar] [CrossRef]
  46. Neff, T.; Feldmeier, H. Tensor correlations in the unitary correlation operator method. Nucl. Phys. A 2003, 713, 311–371. [Google Scholar] [CrossRef]
  47. Roth, R.; Neff, T.; Hergert, H.; Feldmeier, H. Nuclear structure based on correlated realistic nucleon–nucleon potentials. Nucl. Phys. A 2004, 745, 3–33. [Google Scholar] [CrossRef]
  48. Benhar, O.; Biondi, R.; Speranza, E. Short-range correlation effects on the nuclear matrix element of neutrinoless double-β decay. Phys. Rev. C 2014, 90, 065504. [Google Scholar] [CrossRef]
  49. Kotila, J.; Iachello, F. Phase-space factors for double-β decay. Phys. Rev. C 2012, 85, 034316. [Google Scholar] [CrossRef]
  50. Speranza, E. Correlation Effects on the Nuclear Matrix Element of Neutrinoless Double β-Decay. Master’s Thesis, Sapienza University of Rome, Roma, Italy, December 2011. [Google Scholar]
  51. Valli, M. Shear Viscosity of Neutron Matter from Realistic Nucleon-Nucleon Interactions. Ph.D. Thesis, Sapienza University of Rome, Roma, Italy, 2007. [Google Scholar]
  52. Kortelainen, M.; Civitarese, O.; Suhonen, J.; Toivanen, J. Short-range correlations and neutrinoless double beta decay. Phys. Lett. B 2007, 647, 128–132. [Google Scholar] [CrossRef]
  53. Kortelainen, M.; Suhonen, J. Improved short-range correlations and 0νββ nuclear matrix elements of 76Ge and 82Se. Phys. Rev. C 2007, 75, 051303. [Google Scholar] [CrossRef]
  54. Weiss, R.; Soriano, P.; Lovato, A.; Menendez, J.; Wiringa, R.B. Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism. Phys. Rev. C 2022, 106, 065501. [Google Scholar] [CrossRef]
  55. Cirigliano, V.; Dekens, W.; de Vries, J.; Graesser, M.L.; Mereghetti, E.; Pastore, S.; van Kolck, U. New Leading Contribution to Neutrinoless Double-β Decay. Phys. Rev. Lett. 2018, 120, 202001. [Google Scholar] [CrossRef] [PubMed]
  56. Jokiniemi, L.; Soriano, P.; Menéndez, J. Impact of the leading-order short-range nuclear matrix element on the neutrinoless double-beta decay of medium-mass and heavy nuclei. Phys. Lett. B 2021, 823, 136720. [Google Scholar] [CrossRef]
  57. Shimizu, N.; Mizusaki, T.; Utsuno, Y.; Tsunoda, Y. Thick-Restart Block Lanczos Method for Large-Scale Shell-Model Calculations. Comput. Phys. Commun. 2019, 244, 372–384. [Google Scholar] [CrossRef]
  58. Caurier, E.; Nowacki, F.; Poves, A.; Sieja, K. Collectivity in the light xenon isotopes: A shell model study. Phys. Rev. C 2010, 82, 064304. [Google Scholar] [CrossRef]
  59. Menéndez, J. Neutrinoless ββ decay mediated by the exchange of light and heavy neutrinos: The role of nuclear structure correlations. J. Phys. G Nucl. Part. Phys. 2017, 45, 014003. [Google Scholar] [CrossRef]
  60. Šimkovic, F.; Faessler, A.; Rodin, V.; Vogel, P.; Engel, J. Anatomy of the 0νββ nuclear matrix elements. Phys. Rev. C 2008, 77, 045503. [Google Scholar] [CrossRef]
Figure 1. (Color online) Contribution of various coupled spin-parity states ( J π ) of two initial neutrons or two final created protons for (a) Fermi-type and (b) Gamow–Teller-type NMEs. The illustrated NMEs pertain to the light neutrino-exchange mechanism in the context of 0 ν β β decay of 48 Ca. Calculations have been carried out utilizing the nuclear shell model, considering both CD-Bonn and spin-dependent SRCs.
Figure 1. (Color online) Contribution of various coupled spin-parity states ( J π ) of two initial neutrons or two final created protons for (a) Fermi-type and (b) Gamow–Teller-type NMEs. The illustrated NMEs pertain to the light neutrino-exchange mechanism in the context of 0 ν β β decay of 48 Ca. Calculations have been carried out utilizing the nuclear shell model, considering both CD-Bonn and spin-dependent SRCs.
Universe 09 00444 g001
Figure 2. (Color online) The graph illustrates the variation in (a) Fermi and (b) Gamow-Teller NMEs for the 0 ν β β (the light neutrino-exchange mechanism) of 48 Ca with the cutoff number of states ( N c ) of the intermediate nucleus 46 Ca (for TNA calculation). The NMEs are calculated using the GXPF1A Hamiltonian for CD-Bonn-type and spin-dependent SRCs.
Figure 2. (Color online) The graph illustrates the variation in (a) Fermi and (b) Gamow-Teller NMEs for the 0 ν β β (the light neutrino-exchange mechanism) of 48 Ca with the cutoff number of states ( N c ) of the intermediate nucleus 46 Ca (for TNA calculation). The NMEs are calculated using the GXPF1A Hamiltonian for CD-Bonn-type and spin-dependent SRCs.
Universe 09 00444 g002
Figure 3. (Color online) Distribution of NMEs ( C α ( q ) ) with neutrino momentum (q) transfer. The NMEs are computed employing the closure method, considering both CD-Bonn-type SRC parameterization and spin-dependent SRCs. These calculations are executed utilizing the GXPF1A Hamiltonian with closure energy E = 0.5 MeV.
Figure 3. (Color online) Distribution of NMEs ( C α ( q ) ) with neutrino momentum (q) transfer. The NMEs are computed employing the closure method, considering both CD-Bonn-type SRC parameterization and spin-dependent SRCs. These calculations are executed utilizing the GXPF1A Hamiltonian with closure energy E = 0.5 MeV.
Universe 09 00444 g003
Figure 4. (Color online) Radial distribution of NMEs ( C α ( r ) ) with inter-nucleon distance (r). The NMEs have been computed utilizing the closure method, considering both CD-Bonn-type SRC parameterization and spin-dependent SRCs. These calculations are performed using the GXPF1A Hamiltonian with a closure energy E = 0.5 MeV.
Figure 4. (Color online) Radial distribution of NMEs ( C α ( r ) ) with inter-nucleon distance (r). The NMEs have been computed utilizing the closure method, considering both CD-Bonn-type SRC parameterization and spin-dependent SRCs. These calculations are performed using the GXPF1A Hamiltonian with a closure energy E = 0.5 MeV.
Universe 09 00444 g004
Figure 5. (Color online) Variation in total NMEs pertaining to the 0 ν β β decay via the light neutrino-exchange mechanism in 48 Ca, in relation to the closure energy E . These computations employ the GXPF1A Hamiltonian and adopt different SRC parameterizations through the closure method. It is worth noting that the calculation of total NMEs does not encompass the tensor-type NMEs.
Figure 5. (Color online) Variation in total NMEs pertaining to the 0 ν β β decay via the light neutrino-exchange mechanism in 48 Ca, in relation to the closure energy E . These computations employ the GXPF1A Hamiltonian and adopt different SRC parameterizations through the closure method. It is worth noting that the calculation of total NMEs does not encompass the tensor-type NMEs.
Universe 09 00444 g005
Table 1. Parameters for the Jastrow-type function of Equation (15).
Table 1. Parameters for the Jastrow-type function of Equation (15).
SRC Typeabc
Miller–Spencer1.100.681.00
CD-Bonn1.521.880.46
AV181.591.450.92
Table 2. Parameters of f ( r ) and g ( r ) for the spin-dependent SRCs.
Table 2. Parameters of f ( r ) and g ( r ) for the spin-dependent SRCs.
Parameters for f ( r ) ValueParameters for g ( r ) Value
a 1 1.00 a 2 0.04
b 1 0.92 b 2 1.39
c 1 2.56 c 2 2.92
d 1 0.33 d 2 −5.97
e 1 0.57
f 1 −0.94
Table 3. NMEs for 0 ν β β decay (the light neutrino-exchange mechanism) of 48 Ca calculated with nuclear shell model for spin-dependent SRCs and Jastrow-type SRCs with different parameterizations. Calculations are conducted within the closure approximation in the ( n 2 ) channel, employing the GXPF1A effective shell model Hamiltonian.
Table 3. NMEs for 0 ν β β decay (the light neutrino-exchange mechanism) of 48 Ca calculated with nuclear shell model for spin-dependent SRCs and Jastrow-type SRCs with different parameterizations. Calculations are conducted within the closure approximation in the ( n 2 ) channel, employing the GXPF1A effective shell model Hamiltonian.
NME TypeSRC TypeNME Value
M F None−0.215
M F Miller–Spencer−0.144
M F CD-Bonn−0.232
M F AV18−0.213
M F Spin-Dependent SRCs−0.190
M G T None0.774
M G T Miller–Spencer0.540
M G T CD-Bonn0.806
M G T AV180.740
M G T Spin-Dependent SRCs0.674
M 0 ν None0.873
M 0 ν Miller–Spencer0.629
M 0 ν CD-Bonn0.950
M 0 ν AV180.872
M 0 ν Spin-Dependent SRCs0.792
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sarkar, S.; Iwata, Y. Effect of Spin-Dependent Short-Range Correlations on Nuclear Matrix Elements for Neutrinoless Double Beta Decay of 48Ca. Universe 2023, 9, 444. https://doi.org/10.3390/universe9100444

AMA Style

Sarkar S, Iwata Y. Effect of Spin-Dependent Short-Range Correlations on Nuclear Matrix Elements for Neutrinoless Double Beta Decay of 48Ca. Universe. 2023; 9(10):444. https://doi.org/10.3390/universe9100444

Chicago/Turabian Style

Sarkar, Shahariar, and Yoritaka Iwata. 2023. "Effect of Spin-Dependent Short-Range Correlations on Nuclear Matrix Elements for Neutrinoless Double Beta Decay of 48Ca" Universe 9, no. 10: 444. https://doi.org/10.3390/universe9100444

APA Style

Sarkar, S., & Iwata, Y. (2023). Effect of Spin-Dependent Short-Range Correlations on Nuclear Matrix Elements for Neutrinoless Double Beta Decay of 48Ca. Universe, 9(10), 444. https://doi.org/10.3390/universe9100444

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop