MINI REVIEW article

Front. Astron. Space Sci., 28 May 2021

Sec. High-Energy and Astroparticle Physics

Volume 8 - 2021 | https://doi.org/10.3389/fspas.2021.666383

Ordinary Muon Capture for Double Beta Decay and Anti-Neutrino Nuclear Responses

  • 1. Department of Physics, Faculty of Science, Universiti Teknologi Malaysia, Johor Bahru, Malaysia

  • 2. National Centre for Particle Physics, Universiti Malaya, Kuala Lumpur, Malaysia

  • 3. Research Center for Nuclear Physics, Osaka University, Osaka, Japan

Abstract

This is a brief review on ordinary muon capture (OMC) experiments at Research Center for Nuclear Physics (RCNP) Osaka University relevant for the study of double beta decays (DBDs) and astro anti-neutrinos (neutrino) nuclear responses. OMC usually leaves the nucleus in highly excited unbound state. OMC is a charge exchange reaction via the charged weak boson as given by (μ,vμ) reactions with μ and vμ being the muon and muon neutrino. Subjects discussed include 1) unique features of OMC for studying DBDs and astro anti-neutrino (neutrino) nuclear responses, 2) experiments of OMCs on 100Mo and natMo to study neutrino nuclear responses for DBDs and astro anti-neutrinos, 3) impact of the OMC results on neutrino nuclear responses for DBDs and astro anti-neutrinos. Remarks and perspectives on OMC experiments for neutrino nuclear responses are briefly described.

1 Introduction

Double beta decays (DBDs) and astro (solar and supernova) neutrinos and anti-neutrinos are of current interest. Neutrino-less DBDs are used to explore the fundamental neutrino properties such as the Majorana nature, the absolute mass scale, the mass hierarchy and the CP phases beyond the standard electroweak model. The zero-neutrino (0v) DBD rate is given by , where G0v is the phase space volume, meff is the effective neutrino mass and M0v is the zero-neutrino DBD nuclear matrix element, whose squared value is the nuclear response, i.e: . Here we consider the ground state (0+) to ground state (0+) transition AZXX. The DBD nuclear matrix element (NME) M0v is given by the coherent sum of the individual matrix elements M0vi which connect the initial and final ground states via the ith intermediate nucleus AZ+1Xi, where these extend up to about 100 MeV.

The supernova nucleo-synthesis rates induced by neutrinos and anti-neutrinos are proportional to the respective nuclear responseswhere are the NMEs, and (2JA + 1) is the spin factor for the initial nucleus. Note that the NMEs are the coherent sums of individual NMEs for each intermediate state (i) and in the respective isospin directions τ for ν or τ+ for . Note that the ν and NMEs for the ith intermediate state are given by the τ and τ+ NMEs of Mi and Mi+, respectively. Accordingly one needs the B0v response to derive the effective neutrino mass meff from the neutrino-less DBD rate and the astro neutrino (anti-neutrino) response BV () to derive the synthesis rate from the neutrino (anti-neutrino) flux. The neutrino nuclear responses are evaluated theoretically. However, the theoretical calculations for the neutrino nuclear responses for DBD and astro neutrino (anti-neutrino) depend much on the nuclear models and the nuclear parameters used. The neutrino nuclear responses and DBDs are extensively discussed in review articles and references therein (; ).

Experimental studies for the neutrino nuclear responses are interesting to provide the nuclear parameters to be used for theoretical calculations. Charge exchange reactions (CERs) using light ions have been used to study the neutrino nuclear responses. In particular, high energy-resolution (3He, t) CERs have been used to study the neutrino nuclear responses for nuclei of DBD and astro neutrino interests in the broad energy and momentum regions. The responses studied by these light-ion CERs are the τ-side ones, and the NMEs Mi are derived from CERs. On the other hand, the τ+-side responses and the NMEs Mi+ are not well studied since the (t,3He) CERs to be used for studying the τ+-side responses require the radioactive t beam, and thus high precision measurements are difficult. Medium energy-resolution measurements are made by using (d,2He) reactions. Neutrino nuclear responses studied by these light-ion CERs are discussed in the recent review articles and references therein (; ; ; ; ).

Recently ordinary muon capture reactions (OMC), where a negative muon in an atomic orbit is captured into the nucleus has been shown useful for studying the τ+-side responses, and the Mi+ NMEs at Research Center for Nuclear Physics (RCNP) Osaka University (; ; ). The present report is a brief mini-review on the OMC studies for the DBD and astro anti-neutrinos responses and related subjects on nuclear isotope productions at RCNP.

2 Unique Features of OMCs for Neutrino Nuclear Responses

In OMC some of the unique feature are as follows:

  • (1) OMC can be used to study the τ+-responses of proton (p)neutron(n) in the nucleus, and the NME Mi+ for DBD and astro anti-neutrino.

  • (2) OMC transfers energy between 0–50 MeV and a momentum up to 100 MeV/c to the nucleus that are similar to those involved in neutrino-less DBD and supernova anti-neutrino.

  • (3) A negative muon is finally captured into the nucleus via the weak interaction. The capture probability in the medium and heavy nuclei is around 95% after a mean lifetime of about 100 ns in the atomic orbit. Low momentum (a few 10 MeV/c) beam muons at intensities of 103–104 muons/s are used for OMC studies.

  • (4) The OMC on AZX, where A and Z are the mass and the atomic numbers, excites the nucleus AZ−1X up to around 50 MeV, which decays by emitting a number (x) of mostly neutrons and gamma rays to the ground state of the residual nucleus AxZ−1X. The number of neutrons reflects the excitation energy. Accordingly, the relative strength as a function of the excitation energy is evaluated by measuring the number (x) of the neutrons, i.e. the mass distribution Ax of the residual nuclei.

  • (5) Absolute OMC rate is obtained from the measured lifetime of the trapped muons, and thus the neutrino nuclear responses are derived from the OMC rate.

3 Neutrino Nuclear Responses for Mo Isotopes by OMC

The intense 400 MeV proton beam with an intensity around 1 μA from the RCNP cyclotron is used to produce pions. The momentum ( 30–50 MeV/c) negative muons produced by the π-decays are guided by the MuSIC beamline to the target port. They are stopped in the target and are trapped in the inner orbit of the atom. Then, after around 100 ns, the muon is mainly (around 95%) captured into the nucleus via the OMC. The μ-beam spot and the muon beam intensity at the target port are around 6 cm × 6 cm and 5 × 103 per second (). The μ capture rate is derived by measuring the electrons’ time signal from the weak decay of the μ trapped in the atomic orbit by plastic scintillation detectors. The residual isotopes are identified by measuring characteristic prompt γ-rays of the residual nucleus by HPGe detectors online and delayed γ-rays from β-decays of residual isotopes by HPGe detectors offline.

Recently OMCs were studied on natMo and 100Mo, which are interesting for DBD and astro neutrino studies (). The OMC on A42Mo isotopes with the mass-number A = 92–100 produces A41Nb with the excitation energy 0–50 MeV. The excited Nb isotope decays mostly by emitting a number 10(x) of neutrons until the final state gets particle-bound. Then it decays by emitting prompt γ rays to the ground state of AxNb, which is followed by beta decay with comparatively long half-lives.

The OMC strength distribution as a function of the excitation energy E is derived from the mass-number (Ax) distribution (x distribution) through the particle cascade model (; ; ). The excited states in Nb isotopes decay by emitting mostly neutrons since the Coulomb barrier much suppresses proton emission. Thus, the neutron cascade emission model (NEM) (; ) was developed to deduce the excitation energy in the initial isotope of ANb from the number 10(x) of the emitted neutrons, i.e., the mass-number Ax of the residual isotope AxNb after the x neutron emission. The neutron energy spectrum for the first neutron consists of the pre-equilibrium (PEQ) and equilibrium (EQ) stages (). It is given aswhere TEQ(E) and TPEQ(E) are the EQ and PEQ nuclear temperatures, respectively. TEQ(E) is given as a function of excitation energy E (). The ratio of TPEQ(E)/TEQ(E) = 3 for the medium excitation 1040 MeV. After one neutron emission, the emission takes place only via the EQ stage. The OMC strength distribution for 100Nb is derived from the observed mass-number (Ax) distribution for the OMC on 100Mo.

The NEM analysis on the observed mass-number distribution of 100−xNb shows preferential excitation (muon giant resonance, μ-GR) at 10–14 MeV region and a broad bump at the higher excitation region of 25–40 MeV. The NEM analyses on other medium-heavy nuclei by () using experimental data from (; ) show similar features of the preferential excitation of the 10–15 MeV. This preferential excitation reflects the large branch of the one neutron (x = 1) emission. The NME analyses on the mass-number distributions for OMCs on light nuclei show a preferential excitation around 4–8 MeV region by () using experimental data from (; ). This observation is consistent with the calculation in ().

The strength distribution is fitted by the sum of the μ-GR strengths of B1(μ,E1) and B2(μ,E2) given bywhere EGi and Гi with i = 1,2 are the resonance energy and the width for the ith GR, and the constant Bi(μ) is given by σiГi/(2π) with σi being the total strength integrated over the excitation energy. The obtained GR energies are 12 MeV for the first GR and 32.5 MeV for the second GR, and the widths for both GRs are 8 MeV, as shown in Figure 1B. Here we note the GR2 is not entirely clear beyond the statistical error and maybe a broad bump. The μ-GR energy peaks at around 12 MeV, which is lower than the giant dipole resonance (GDR) energy of 14 MeV deduced from a photon-induced reaction (). It also features a broader width of about 8 MeV compared to the 5 MeV width of the latter. This is, of course due to the additional multipole components of higher spin states excited by τ+, i. e.,: (n, p)-type, muon capture.

FIGURE 1

). (A) When OMC on 100Mo occurs, 100Nb is excited to some energy up to around 50 MeV, with Bn(Bp) = 5.68 MeV (9.46 MeV), 6.87 MeV (8.34 MeV), 5.99 MeV (7.87 MeV) and 8.07 MeV (7.46 MeV) for NbA with A = 100, 99, 98, and 97, respectively. (B) The excitation energy region of 8–16 MeV, which is strongly excited, decays by emitting mainly 1 neutron, and partly by 2 neutrons, depending on the neutron energy. (C) The excitation region of 26–34 MeV day by emitting 3–4 neutrons and also partly by emitting 3 neutrons and 1 proton. (D) The muon strength distribution, B(μ,E) for OMC on 100Mo from the excitation energy distribution in the daughter nucleus of Nb and Zr isotopes following the muon capture process. To reproduce the 100Mo experimental RI mass distribution, the EG1 and EG2 are observed at 12 MeV, and 30 MeV with an intensity ratio of EG1/EG2 = 1/6 ().

Recently a proton-neutron emission model (PNEM) has been developed by considering the probability of proton emission as well (). The proton energy spectrum is given by S(Ep) = m × exp(−n×Ep), where m and n are the density parameter constant with the value of 0.832 MeV and 0.163 MeV−1. Here, the proton is assumed to be emitted in the case of with BP, E and Bn being the proton binding energy, the excitation energy and the neutron binding energy, because the neutron decays are forbidden and gamma decays are much smaller than the proton emission even below the Coulomb barrier. The PNEM is shown schematically in Figure 1A. The analysis of the mass-number distribution of 100−xNb shows the same GR1 and GR2 as in case of the NEM analysis. The missing Zr isotopes produced by a proton emission in experimental data are predicted by PNEM for both natMo () and 100Mo (). The OMC strength distribution shown in Figure 1A reproduces the observed radioactive isotope (RI) mass distribution.

Since then, various calculations have been made for understanding the formation of the GR populated by muon capture reactions in 98, 209 (). These calculations have been compared with previous experimental works by OMC on the nuclei reported in references (; ; ). The parameters of EG1 and EG2 as a function of A are given as EG1 = 30A−1/5 and EG2 = 75A−1/5 for OMC on natMo, 100Mo, 107Pd, 108Pd, 127I and 209Bi. They are obtained from a comparison of PNEM calculations with experimental data. The observed mass distributions agree with the predictions as shown in Figure 2.

FIGURE 2

). The black bar represents the experimental observation while the gray bar represents calculated value by Proton neutron emission model (PNEM). Errors obtain for the experiment and calculation are included based on statistical and systematical study for each parameters assumption. The value of model error are are assumed to be small since the assumption made are suitable for particle evaporation from medium-heavy nuclei (A)natMo (isotopic enrichment of 9.7% of 100Mo, 24.3% of 98Mo, 9.6% of 97Mo, 16.7% of 96Mo, 15.9% of 96Mo), 9.2% of 95Mo and 14.6% of 94Mo) (), (B)100Mo (isotopic enrichment of 96% of 100Mo, 2.5% of 98Mo, 0.5% of 97Mo and 1% of 96Mo) (), (C)127I (isotopic enrichment of 100% of 127I) () and (D)209Bi (isotopic enrichment of 100% of 219Bi) ().

4 Implications for Double Beta Decays

The experimental OMC rate for 100Mo as a function of the excitation energy is compared with the theoretical calculation by using proton-neutron quasi-particle random-phase approximation (pn-QRPA) (). The pn-QRPA well reproduces the GR at around 10–14 MeV with mainly low multipoles of Jπ = 0±,1± and 2±. The observed OMC rate is the same order of magnitude as the empirical value suggested by Primakoff () but is lower by a factor 5 than the pn-QRPA (; ). This suggests a quenching coefficient for the axial-vector weak coupling of 0.5 (; ) for μ-NMEs, being consistent with the quenching factors for Gamow-Teller (GT) and spin dipole (SD) NMEs (; ). Here the quenching is common in the wide momentum region (; ). The quenching effects are also discussed theoretically (; ; ).

It is interesting to note that the OMC results together with the (3He,t) CERs and data suggest severe quenching of experimental axial-vector NMEs of Mi and Mi+ for multipoles of Jπ = 1+, 2, 3+ (; ). The common quenching coefficients of kNM = gAeff/gA for them are around 0.4–0.6 with respect to the pn-QRPA. In this case, the axial-vector DBD NMEs are quenched by the coefficient k2NM and the axial-vector anti-neutrino NMEs by kNM. On the other hand, the recent pn-QRPA calculations for OMC rates show no quenching () or only weak quenching (). The different quenching coefficients are partly due to the different model space and the different muon-wave functions used in their QRPA calculations (; ; ; ). Further studies are needed to see how pn-QRPA calculations reproduce both the relative and absolute OMC rates.

5 Remarks and Perspectives

The relative strength distribution of OMC on Mo isotopes show the μ-GR around 12 MeV consistent with the pn-QRPA calculation (). However, the absolute strength derived from the lifetime is much smaller than the model value, suggesting a similar quenching of the Mi+ NMEs as the Mi NMEs derived from light-ion CERs. Another pn-QRPA () shows no severe quenching. The quenching for μ NMEs remains to be studied.

The extensive experimental programs on OMC for other nuclei of DBD and supernova anti-neutrino interests are under progress at RCNP and Paul Scherrer Institute (PSI), Switzerland, by the join group of Joint Institute for Nuclear Research (JINR), Dubna, RCNP, Osaka and Universiti Teknologi Malaysia (UTM) (). The pn-QRPA theoretical calculation is also in progress at Jyvaskyla ().

Finally, it is remarked that muon capture isotope production (MuCIP) is used for producing efficiently nuclear isotopes with the atomic number Z−1, less by one than the atomic number Z of the target isotope. By using AZX target isotopes, isotopes of A−1Z−1X are preferentially produced, and several isotopes with A,A−2,A−3,A−4 are also produced (; ). There is potential for using them for basic and applied science. One RI used for medical checks is the 99Mo, which is well produced by OMC on 100Mo, as studied at RCNP (). MuCIP is complementary to photon capture reactions where isotopes of A−1ZX are well produced (; ). Note that ppb-level nuclei (impurities) are identified by measuring gamma rays from OMC, which are characteristic of the nuclei, as explained in ().

Statements

Author contributions

IH and HE are equally accounted for the content of the work.

Funding

This work was financially supported by the Universiti Teknologi Malaysia research grants (Q.J130000.3026.01M14 and R.J130000.7854.5F227) and the Ministry of Higher Education Malaysia Fundamental Research Grant Scheme (FRGS/1/2019/STG02/UTM/02/6).

Acknowledgments

We are grateful to RCNP, UTM, and JINR colleagues for the collaboration and valuable discussions.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

ordinary muon capture, muon charge exchange reaction, neutrino nuclear response, double beta decay, supernova neutrino, nuclear matrix element, neutrino mass

Citation

Hashim IH and Ejiri H (2021) Ordinary Muon Capture for Double Beta Decay and Anti-Neutrino Nuclear Responses. Front. Astron. Space Sci. 8:666383. doi: 10.3389/fspas.2021.666383

Received

10 February 2021

Accepted

05 May 2021

Published

28 May 2021

Volume

8 - 2021

Edited by

Theocharis S. Kosmas, University of Ioannina, Greece

Reviewed by

Dieter Frekers, University Muenster, Germany

Frank Franz Deppisch, University College London, United Kingdom

Updates

Copyright

*Correspondence: Izyan Hazwani Hashim,

This article was submitted to High-Energy and Astroparticle Physics, a section of the journal Frontiers in Astronomy and Space Sciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics