# First-principles results for nuclear – decay properties of shell nuclei

###### Abstract

We evaluate the allowed - decay properties of nuclei with systematically under the framework of the nuclear shell model with the use of the valence space Hamiltonians derived from modern methods, such as in-medium similarity renormalization group and coupled-cluster theory. For comparison we also show results obtained with the newly developed shell-model interactions from chiral effective field theory and phenomenological USDB interaction. In this work, we have tested predictive power of effective interactions by comparing calculated results with the experimental data for – decay properties of shell nuclei. We have performed calculations for O F, F Ne, Ne Na, Na Mg, Mg Al, Al Si, Si P and P S transitions. Theoretical results of , log values and half-lives, are discussed and compared with the available experimental data. calculations of -decay properties are very limited for the shell nuclei, thus present comprehensive study will add more information to it.

###### pacs:

21.60.Cs - shell model, 23.40.-s --decay## I Introduction

The study of -decay properties of unstable-nuclei have been extensively investigated in the past brownwild (). There are several experimental data available for half-lives, log values, Gamow-Teller () strengths, values and branching fractions nndc (); ame2012 (); ame2016 (). On the other hand with the recent development in approaches it is possible to predict these properties. The -decay half-lives are very important for understanding -process nucleosynthesis. The strength is one of the important tools to study the structure of atomic nuclei. The experimental strengths can be obtained from -decay studies and charge-exchange (CE) reactions. The reliable estimates of strength distributions in neutron-rich nuclei can be of great interest for description of -decay properties.

Study of -decays based on methods is very limited, and used to be available only for few-body systems and several light nuclei. Study of transitions with the inclusion of the effects of three-nucleon forces and two-body currents from chiral effective field theory are reported for N and O in Ref. andres (). The origin of anomalous long life time for the -decay of C to N was investigated using no-core shell model with the Hamiltonian from the chiral effective field theory including three-nucleon force terms moris (). The no-core shell model results for strengths for -shell nuclei with , using two-nucleon () and three-nucleon () interactions derived from chiral effective field theory, were reported in Ref. PRL990425012007 (). The chiral low-energy constants and are constrained by means of accurate calculations of the binding energies of the system and half-life of triton gazit (). The study of the tritium -decay with the chiral effective field theory was reported in Ref. baroni (). The uncertainties in constraining low-energy constants from H decay were reported in Ref. klos ().

In very early years, features of the beta-decay of neutron-rich shell nuclei with five or more excess neutrons were predicted by Wildenthal et al wcb (). The comprehensive study of -decay properties of -shell nuclei for = 17–39 were reported by Brown and Wildenthal in Ref. brownwild (). The Gamow-Teller beta-decay rates for A18 nuclei were reported in Ref. chou (). For the -shell nuclei, the calculated Gamow-Teller matrix elements of 64 decays of nuclei in the mass range A=41–50 were reported in Ref. pinedo (). In these studies, a substantial quenching of the axial-vector coupling in the strength was found in -shell and -shell nuclei about by 20 and also in -shell nuclei by 25-26. Theoretical studies on the quenching of the strengths were carried out for nuclei with closed-core 1 nucleon Arima (); Towner (), and the effects of coupling to 2p-2h configurations and roles of two-body meson-exchange currents were found to be important.

In later years, the shell model calculations for -decay properties of neutron-rich nuclei with were reported by Li and Ren in Ref. li (). The nuclear -decay half-lives for and shell nuclei were reported in Ref. kumar_jpg (). The shell model description of strengths in -shell nuclei were available in Ref. vik_epja (). The systematic shell-model study of -decay properties and strength distributions in neutron-rich nuclei were reported in Ref. otsuka (). Theoretical calculations for half-lives of medium-mass and heavy mass neutron-rich nuclei from QRPA based on the Hartree-Fock Bogoliubov theory or other global models were available in the literature suhonen (); q1 (); q2 (); q3 (); ser ().

More recently, results of the study on and double-beta decays of heavy nuclei within a framework of an effective theory were presented in Ref. perez (). The calculations of strengths in shell nuclei for 13 different nuclear transitions including electron-capture reaction rates for NaNe and MgNa were reported in Ref. archana_gt (). Here, a need for the quenching of the strength was reported also for the shell-model interactions.

In the present work we have reported – decay properties of nuclei using nuclear shell model based on interactions and newly developed shell-model interactions from chiral effective field theory. The purpose of the present work is to study how well the recent and newly developed shell-model interactions based on chiral interactions can describe the -decay properties in -shell, and also to find how much quenching is necessary for these interactions by comparing with many more experimental data than in Ref. archana_gt (). This work will add more information to earlier works wcb (); brownwild (); li (), where shell model results with phenomenological effective interactions were reported.

This paper is organized as follows. In Sec. II, we present details of interactions. The formalism of the calculations for -decay properties are presented in Sec. III. In Sec. IV, we present theoretical results along with the experimental data. Finally, a summary and conclusions are drawn in Sec. V.

## Ii Ab initio Hamiltonians

To calculate , log values and half-lives for the shell nuclei, we have performed shell-model calculations using two interactions : IM-SRG stroberg () and CCEI jansen (); jan1 (). Also we have performed calculations with newly developed shell-model interaction derived from the chiral effective field theory Huth (). For comparison, we have also performed calculations with the phenomenological USDB effective interaction usdb () in addition to the above three interactions. For the diagonalization of matrices we used J-scheme shell-model code NuShellXnushellx ().

Stroberg et al. stroberg () derived a mass-dependent valence-space Hamiltonians using the IM-SRG Tsukiyama (); bogner () including the chiral and nucleon interactions (see Refs. hergertpr (); hergertps (); lecture () for more details about IM-SRG). In the IM-SRG, a continuous unitary transformation, parametrized by the flow parameter s, is applied to the initial normal-ordered A-body Hamiltonian:

(1) |

where the and are the diagonal and off-diagonal part of the Hamiltonian, respectively. The off-diagonal Hamiltonian become zero as

(2) |

The resulting Hamiltonian is used in the shell-model calculations for shell nuclei. In the present work we use effective interactions with = 24 MeV (see details in Ref. stroberg () for other parameters).

We use another approach to study -decay properties, named as coupled-cluster effective interactions (CCEI) jansen (); jan1 ().

We can expand the Hamiltonian for the suitable model-space using the valence-cluster expansion lisetskiy () given as

(3) |

Here A is the mass of the nucleus for which we are doing calculations, is the core Hamiltonian, is the valence one-body Hamiltonian, and is the two-body Hamiltonian. The two-body term is derived from Eq. (3) by using the Okubo-Lee-Suzuki (OLS) similarity transformation okubo (); suzuki (). After using this unitary transformation the effective Hamiltonian becomes non-Hermitian. For changing the non-Hermitian to Hermitian effective Hamiltonian the matrix operator [] = is used, where S is a matrix that diagonalize the Hamiltonians (see Ref. nbarrett () for further details ). The Hermitian shell-model Hamiltonian is obtained as .

The and parts are taken from a next-to-next-to-next-to leading order (N3LO) chiral nucleon-nucleon interaction, and a next-to-next-to leading order (N2LO) chiral three-body interaction, respectively. For both IM-SRG and CCEI, we use = 500 MeV for chiral N3LO interaction entem (); machleidt (), and = 400 MeV for chiral N2LO interaction navratil (), respectively. The low-energy constants (LECs) of 3NF are given as = 0.098 and navratil ().

In the CCEI to achieve faster model-space convergence, the similarity renormalization group transformation has been used to evolve two-body and three-body forces to the lower momentum scale (see Ref. Jurgenson () for further details). Also, for the coupled-cluster calculations, a Hartree-Fock basis built from thirteen major harmonic-oscillator orbitals with frequency = 20 MeV have been used. Using IM-SRG targeted for a particular nucleus imsrg_ragnar () and CCEI interactions, the shell model results for spectroscopic factors and electromagnetic properties are reported in Refs. pcs_prc (); archana_prc ().

Recently, L. Huth et al. Huth () derived a shell-model interaction from chiral effective field theory. The valence-space Hamiltonian for shell is constructed as a general operators having two low energy constants (LECs) at leading order (LO) and seven new LECs at next-to-leading order (NLO) and fitted the LECs of CEFT operators directly to 441 ground- and excited-state energies. For the chiral EFT interaction they have taken the expansion in terms of power of based on Weinberg power counting wibnernpb (), where Q is a low-momentum scale or pion mass and 500 MeV is the chiral-symmetry-breaking scale.

## Iii Formalism

In the beta decay, the value corresponding to transition from the initial state of the parent nucleus to the final state in the daughter nucleus is expressed as Piechaczek ()

(4) |

where is the Gamow-Teller transition strength, and is the axial vector phase space factor that contains the lepton kinematics. In this work, we have calculated the phase space factor with parameters given by Wilkinson and Macefield wilkinson () together with the correction factors given in Refs. Sirlin (); wga (). The values are very large, so they are defined in term of “log” values. The log is expressed as log= log.

The total half-life is related to the partial half-life as

(5) |

where runs over all the possible daughter states that are populated through transitions.

The partial half-life is related to the total half-life of the allowed -decay as

(6) |

where, is called the branching ratio for the transition with partial half-life .

The Gamow-Teller strength is calculated using the following expression:

(7) |

where (=-1.260) is the axial-vector coupling constant of the weak interaction, and and are the initial and final state shell-model wave functions, respectively, here the refers to isospin operator for the decay, for the -decay we use the convention = , is the initial-state angular momentum, and is the quenching factor.

Following Refs. brownwild (); chou (); pinedo (), we define

(8) |

which is independent of the direction of the transitions.

values are defined as

(9) |

where the total strength is defined by

(10) |

In the -decay the endpoint energy of the electron (in units of MeV) is an essential quantity to calculate the phase space factor . is given by the expression brownwild ()

(11) |

where the Q is the -decay Q value, and and are excitation energies of the initial and final states. Here, we have taken Q values from the experimental data ame2016 ().

## Iv Results and Discussions

In Table 1 we compare calculated and experimental values of the matrix elements . Calculated values of presented here are those with =1. The -decay energies (E (decay)), branching ratios () and log values as well as the values of W are given in Table 1. The quenching factors are obtained by chi-squared fitting of the theoretical values to the corresponding experimental values. The quenching factors as well as the root-mean-square (RMS) deviations for the effective interactions considered here are given in Table 2. The quenching factors are slightly different for different effective interactions. Their values are in the range of =0.62-0.77. The value for USDB, =0.770.02, is consistent with the one, =0.764, in Ref. Richter (). The RMS deviations for the and CEFT interactions are enhanced compared with the USDB by 25-32 and 52, respectively.

We have plotted the experimental values with respect to the theoretical values for the shell nuclei in Fig.1. For further calculations of observables, we take the quenching factors from the Table 2.

M(GT) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

Z | Z | 2J,2T | E(decay) (MeV) | (%) | log(exp.) | EXPT. | USDB | IM-SRG | CCEI | CEFT | W |

O() | F | 7,1 | 0.442 | 0.0984(30) | 3.86(17) | 2.262 | 3.406 | 3.910 | 3.640 | 2.334 | 9.259 |

5,1 | 4.622 | 45.4(15) | 5.38(15) | 0.393 | 0.243 | 0.256 | 0.416 | 0.593 | |||

3,1 | 3.266 | 54.4(12) | 4.62(10) | 0.939 | 1.245 | 1.468 | 1.556 | 0.596 | |||

O() | F | 2,2 | 2.757 | 99.97(3) | 3.73(6) | 1.072 | 1.104 | 1.399 | 1.527 | 0.887 | 4.365 |

2,2 | 0.325 | 0.027(3) | 3.64(6) | 1.190 | 1.281 | 0.989 | 0.963 | 0.025 | |||

O() | F | 3,3 | 6.380 | 37.2(12) | 5.22(2) | 0.473 | 0.464 | 0.367 | 0.360 | 0.223 | 11.953 |

O() | F | 2,4 | 4.860 | 31(5) | 4.6(1) | 0.394 | 0.403 | 0.509 | 0.473 | 0.573 | 5.346 |

2,4 | 3.920 | 68(8) | 3.8(1) | 0.989 | 1.234 | 0.066 | 1.686 | 0.675 | |||

O() | F | 2,6 | 9.700 | 40(4) | 4.3(1) | 0.556 | 0.990 | 0.782 | 0.857 | 1.024 | 6.173 |

F() | Ne | 4,0 | 5.390 | 99.99(8) | 4.97(11) | 0.575 | 0.687 | 0.679 | 0.651 | 0.694 | 6.901 |

F() | Ne | 7,1 | 3.938 | 16.1(10) | 4.72(3) | 0.840 | 0.959 | 1.046 | 1.123 | 0.655 | 9.259 |

5,1 | 5.333 | 74.1(22) | 4.65(1) | 0.911 | 0.982 | 1.143 | 1.338 | 0.858 | |||

3,1 | 5.684 | 9.6(30) | 5.67(16) | 0.281 | 0.371 | 0.421 | 0.515 | 0.360 | |||

F() | Ne | 10,2 | 3.480 | 8.7(4) | 4.70(2) | 1.053 | 1.250 | 1.372 | 1.362 | 0.940 | 13.094 |

8,2 | 7.461 | 3.1(6) | 6.7(1) | 0.105 | 0.106 | 0.156 | 0.099 | 0.066 | |||

8,2 | 5.500 | 53.9(6) | 4.79(1) | 0.950 | 1.174 | 1.173 | 1.365 | 1.048 | |||

8,2 | 4.670 | 7.0(3) | 5.34(2) | 0.504 | 0.406 | 0.107 | 0.287 | 0.706 | |||

8,2 | 3.477 | 2.45(22) | 5.30(4) | 0.528 | 0.733 | 0.939 | 1.280 | 0.357 | |||

6,2 | 5.177 | 16.4(7) | 5.26(2) | 0.553 | 0.678 | 0.649 | 0.660 | 0.382 | |||

F() | Ne | 5,3 | 8.480 | 30(8) | 5.72(16) | 0.266 | 0.377 | 0.333 | 0.321 | 0.428 | 11.953 |

3,3 | 6.660 | 10.9(19) | 5.66(11) | 0.285 | 0.262 | 0.256 | 0.305 | 0.334 | |||

3,3 | 5.050 | 15.2(12) | 4.96(8) | 0.637 | 0.866 | 1.028 | 1.171 | 0.881 | |||

3,3 | 4.650 | 25(4) | 4.58(12) | 0.987 | 0.200 | 0.260 | 0.272 | 0.044 | |||

F() | Ne | 4,6 | 16.170 | 36(7) | 4.6(1) | 0.682 | 1.145 | 1.079 | 1.065 | 0.858 | 10.691 |

0,6 | 18.900 | 36.5(60) | 4.9(1) | 0.483 | 0.733 | 0.756 | 0.714 | 0.719 | |||

Ne() | Na | 5,1 | 3.950 | 32.0(13) | 5.38(2) | 0.393 | 0.372 | 0.392 | 0.753 | 0.656 | 9.259 |

3,1 | 4.383 | 66.9(13) | 5.27(1) | 0.446 | 0.355 | 0.585 | 0.888 | 0.516 | |||

Ne() | Na | 2,2 | 1.994 | 92.1(2) | 4.35(1) | 0.525 | 0.571 | 0.442 | 0.664 | 0.053 | 4.365 |

2,2 | 1.120 | 7.9(2) | 4.39(2) | 0.502 | 0.542 | 0.848 | 1.068 | 0.986 | |||

Ne() | Na | 3,3 | 7.160 | 76.6(20) | 4.41(2) | 0.693 | 0.751 | 0.751 | 0.697 | 0.604 | 6.901 |

1,3 | 6.180 | 19.5(20) | 4.70(6) | 0.496 | 0.597 | 0.419 | 0.547 | 0.724 | |||

1,3 | 2.960 | 0.53(15) | 4.82(16) | 0.432 | 0.595 | 0.713 | 1.037 | 0.657 | |||

Ne() | Na | 2,4 | 7.258 | 91.6(2) | 3.87(6) | 0.913 | 1.110 | 1.302 | 1.247 | 1.036 | 5.346 |

2,4 | 5.829 | 4.2(4) | 4.8(1) | 0.313 | 0.669 | 0.022 | 0.309 | 0.666 | |||

2,4 | 4.619 | 1.9(4) | 4.7(1) | 0.351 | 0.616 | 0.547 | 0.818 | 0.892 | |||

Ne() | Na | 5,5 | 12.590 | 59.5(30) | 4.40(4) | 0.992 | 1.229 | 1.042 | 1.059 | 1.258 | 11.548 |

Ne() | Na | 2,6 | 12.280 | 55(5) | 4.2(1) | 0.624 | 1.185 | 1.104 | 1.039 | 1.147 | 6.173 |

2,6 | 10.350 | 1.7(4) | 5.3(1) | 0.176 | 0.816 | 0.734 | 0.352 | 0.679 | |||

2,6 | 10.160 | 20.1(12) | 4.2(1) | 0.624 | 0.341 | 0.484 | 0.677 | 0.333 | |||

2,6 | 9.570 | 8.5(6) | 4.5(1) | 0.442 | 0.471 | 0.394 | 1.293 | 0.189 | |||

Na() | Mg | 8,0 | 1.392 | 99.855(5) | 6.11(1) | 0.208 | 0.338 | 0.125 | 0.065 | 0.355 | 9.259 |

Na() | Mg | 7,1 | 2.223 | 9.48(14) | 5.03 | 0.588 | 0.642 | 0.629 | 0.577 | 0.471 | 9.259 |

5,1 | 3.835 | 62.5(20) | 5.26 | 0.451 | 0.558 | 0.636 | 1.169 | 0.685 | |||

3,1 | 2.860 | 27.46(22) | 5.04 | 0.581 | 0.708 | 0.842 | 0.748 | 0.633 | |||

3,1 | 1.033 | 0.247(4) | 5.25 | 0.457 | 0.683 | 0.712 | 1.010 | 0.910 | |||

Na() | Mg | 6,2 | 5.413 | 1.31(4) | 5.87(1) | 0.242 | 0.226 | 0.377 | 0.316 | 0.580 | 11.548 |

6,2 | 5.004 | 3.17(7) | 5.33(1) | 0.450 | 0.714 | 0.536 | 1.142 | 0.504 | |||

6,2 | 3.229 | 1.72(4) | 4.74(1) | 0.887 | 1.051 | 0.801 | 1.713 | 1.549 | |||

4,2 | 7.545 | 87.80(7) | 4.71(1) | 0.918 | 1.073 | 0.492 | 0.865 | 0.836 | |||

4,2 | 6.416 | 0.05(4) | 7.60(4) | 0.033 | 0.129 | 0.404 | 0.115 | 0.047 | |||

4,2 | 5.022 | 1.65(3) | 5.62(1) | 0.322 | 0.411 | 1.037 | 0.971 | 0.552 | |||

4,2 | 4.519 | 2.738(19) | 5.25(1) | 0.493 | 0.558 | 0.149 | 0.650 | 0.648 |

M(GT) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

Z | Z | 2J,2T | E(decay) (MeV) | (%) | log(exp.) | EXPT. | USDB | IM-SRG | CCEI | CEFT | W |

Na() | Mg | 5,3 | 7.310 | 11.3(7) | 4.99(3) | 0.616 | 0.602 | 0.406 | 0.478 | 0.650 | 11.953 |

3,3 | 8.030 | 85.8(11) | 4.30(15) | 1.363 | 1.747 | 1.683 | 1.562 | 1.509 | |||

Na() | Mg | 4,4 | 12.556 | 11(6) | 5.1(2) | 0.384 | 0.294 | 0.438 | 0.520 | 0.245 | 9.259 |

2,4 | 9.469 | 3.2(4) | 5.1(1) | 0.384 | 0.536 | 0.583 | 0.592 | 0.925 | |||

0,4 | 14.030 | 60(5) | 4.6(1) | 0.682 | 0.840 | 0.702 | 0.773 | 0.882 | |||

0,4 | 10.168 | 20.1(19) | 4.42(1) | 0.839 | 1.116 | 1.174 | 0.961 | 0.787 | |||

Na() | Mg | 3,5 | 13.272 | 24(8) | 5.06(15) | 0.464 | 0.786 | 0.737 | 0.832 | 0.739 | 11.548 |

Na() | Mg | 4,6 | 15.790 | 9.5(11) | 5.86(6) | 0.206 | 0.408 | 0.451 | 0.509 | 0.411 | 13.803 |

Mg() | Al | 3,1 | 1.596 | 29.06(9) | 4.934(16) | 0.381 | 0.450 | 0.373 | 0.754 | 0.468 | 5.346 |

1,1 | 1.766 | 70.94(9) | 4.73(10) | 0.480 | 0.597 | 0.178 | 0.473 | 0.766 | |||

Mg() | Al | 2,2 | 0.459 | 94.8(10) | 4.45(9) | 0.468 | 0.624 | 0.379 | 0.819 | 0.581 | 4.365 |

2,2 | 0.211 | 4.9(10) | 4.57(9) | 0.408 | 0.495 | 0.454 | 0.291 | 1.029 | |||

Mg() | Al | 5,3 | 7.613 | 27(8) | 5.32(14) | 0.344 | 0.579 | 0.403 | 0.409 | 0.829 | 9.760 |

5,3 | 4.551 | 6.0(16) | 4.93(13) | 0.539 | 1.064 | 1.194 | 1.769 | 1.193 | |||

5,3 | 4.428 | 28(5) | 4.21(9) | 1.234 | 1.098 | 0.367 | 0.518 | 0.954 | |||

3,3 | 5.389 | 21(6) | 4.73(13) | 0.678 | 0.658 | 0.341 | 0.918 | 0.187 | |||

3,3 | 4.747 | 7.8(15) | 4.90(10) | 0.558 | 0.669 | 0.806 | 0.036 | 0.833 | |||

1,3 | 6.215 | 7(3) | 5.49(19) | 0.283 | 0.283 | 0.198 | 0.036 | 0.067 | |||

1,3 | 4.180 | 3.0(9) | 5.06(14) | 0.464 | 0.710 | 0.383 | 0.571 | 0.380 | |||

Mg() | Al | 2,4 | 6.274 | 68(20) | 3.96(13) | 0.823 | 1.203 | 1.167 | 0.994 | 1.337 | 5.346 |

2,4 | 4.549 | 7(1) | 4.30(7) | 0.556 | 0.870 | 0.783 | 1.333 | 0.800 | |||

Mg() | Al | 2,6 | 10.150 | 55 | 4.4 | 0.496 | 1.596 | 1.531 | 1.550 | 1.450 | 6.173 |

2,6 | 7.380 | 24.6(8) | 4.1 | 0.700 | 0.303 | 0.385 | 0.063 | 0.190 | |||

2,6 | 6.950 | 10.7(10) | 4.4 | 0.496 | 0.013 | 0.072 | 0.101 | 0.025 | |||

Al() | Si | 4,0 | 2.863 | 99.99(1) | 4.87(4) | 0.764 | 0.945 | 0.353 | 0.983 | 0.920 | 8.166 |

Al() | Si | 3,1 | 2.406 | 89.9(3) | 5.05(5) | 0.575 | 0.924 | 0.237 | 0.388 | 0.821 | 9.259 |

3,1 | 1.253 | 6.3(2) | 5.03(15) | 0.591 | 0.589 | 0.591 | 1.657 | 1.059 | |||

Al() | Si | 6,2 | 3.730 | 6.6(2) | 4.985(17) | 0.669 | 1.001 | 0.673 | 0.917 | 0.921 | 11.548 |

6,2 | 3.329 | 2.6(2) | 5.17(4) | 0.541 | 0.657 | 0.229 | 0.869 | 0.167 | |||

4,2 | 6.326 | 17.1(9) | 5.619(25) | 0.322 | 0.362 | 0.379 | 0.531 | 0.157 | |||

4,2 | 5.063 | 67.3(11) | 4.578(12) | 1.069 | 1.417 | 1.192 | 0.933 | 1.189 | |||

4,2 | 3.751 | 5.7(2) | 5.06(19) | 0.614 | 0.762 | 0.067 | 0.476 | 0.561 | |||

Al() | Si | 4,4 | 11.080 | 4.7((13) | 5.29(13) | 0.308 | 0.181 | 0.122 | 0.317 | 0.414 | 9.259 |

4,4 | 8.790 | 3.0(8) | 5.00(12) | 0.430 | 0.590 | 0.754 | 0.449 | 0.625 | |||

0,4 | 13.020 | 85(5) | 4.36(3) | 0.899 | 1.093 | 0.846 | 0.952 | 1.157 | |||

0,4 | 8.040 | 4.3(11) | 4.66(12) | 0.637 | 1.229 | 1.336 | 1.206 | 0.938 | |||

Al() | Si | 3,5 | 11.960 | 88(2) | 4.3 | 1.363 | 2.225 | 1.966 | 1.491 | 2.080 | 14.143 |

Si() | P | 1,1 | 1.491 | 99.94(7) | 5.525(8) | 0.272 | 0.318 | 0.076 | 0.600 | 0.389 | 7.560 |

Si() | P | 2,2 | 0.227 | 100 | 8.21(6) | 0.006 | 0.069 | 0.136 | 0.908 | 0.384 | 4.365 |

Si() | P | 1,3 | 5.845 | 93.7(7) | 4.96(17) | 0.520 | 0.264 | 0.394 | 0.313 | 0.329 | 9.760 |

Si() | P | 2,4 | 2.984 | 100 | 3.3 | 1.759 | 0.639 | 0.372 | 0.845 | 0.232 | 5.346 |

P() | S | 0,0 | 1.710 | 100 | 7.90(2) | 0.015 | 0.038 | 0.140 | 0.543 | 0.352 | 5.346 |

P() | S | 3,1 | 0.248 | 100 | 5.022(7) | 0.343 | 0.295 | 0.025 | 0.336 | 0.337 | 5.346 |

P() | S | 4,2 | 3.255 | 14.8(20) | 4.93(6) | 0.467 | 0.742 | 0.533 | 0.593 | 0.837 | 7.560 |

4,2 | 1.268 | 0.31(6) | 4.88(9) | 0.494 | 0.082 | 0.313 | 1.573 | 0.229 | |||

0,2 | 5.383 | 84.8(21) | 5.159(12) | 0.358 | 0.480 | 0.180 | 1.545 | 0.269 |

Interaction | RMS deviations | |
---|---|---|

USDB | 0.770.02 | 0.0356 |

IM-SRG | 0.750.03 | 0.0469 |

CCEI | 0.620.03 | 0.0440 |

CEFT | 0.730.04 | 0.0541 |

F(5/2) Ne | EXPT | USDB | IM-SRG | CCEI | CEFT |
---|---|---|---|---|---|

log (5/2) | 4.650.01 | 4.81 | 4.70 | 4.73 | 4.97 |

half-life(s) | 4.1580.0205 | 5.30 | 5.19 | 5.41 | 6.93 |

Ne(0) Na | |||||

log (1) | 4.20.01 | 3.87 | 3.95 | 4.17 | 3.94 |

half-life(ms) | 201 | 13.91 | 18.65 | 25.64 | 17.50 |

Na(1) Mg | |||||

log (0) | 4.60.01 | 4.65 | 4.82 | 4.91 | 4.65 |

half-life(ms) | 30.50.4 | 36.65 | 55.16 | 69.89 | 42.07 |

Mg(0) Al | |||||

log (1) | 4.4 | 3.61 | 3.67 | 3.83 | 3.74 |

half-life(ms) | 865 | 22.84 | 26.12 | 37.53 | 30.91 |

Mg(0) Al() | |||||

log (1) | 4.450.09 | 4.43 | 4.88 | 4.38 | 4.54 |

half-life(h) | 20.90.009 | 5.92 | 1.49 | 3.59 | 22.96 |

In Fig. 2 we show the distribution of calculated log values with the experimental data for some –decays nuclei for which experimental log values are available for excited states also. In case of F, although results of the interactions for excitation energy for the excited and states slightly differ from the experimental data, all the interactions give calculated log values close to the experimental data. For Ne, the calculated log values with the CCEI are better in comparison to other interactions. The calculated value for excitation energy for state is in good agreement for all the interactions for Na. For this nucleus all the four interactions give reasonable results for log values.

In Table 3, we compare the theoretical and the experimental log and -decay half-lives of several selected shell nuclei. Calculated and experimental half-lives are in fairly good agreement in the cases shown in Fig. 2. For F(5/2) Ne, Ne(0) Na and Na(1) Mg transitions, the difference in the half-lives is almost determined by the difference of the B(GT) value of the transition to the state with the largest branching ratio, and also difference in phase space factors. For F(5/2) Ne, the transition to the g.s. (3/2) gives some contributions. For Mg(0) Al transition, the agreement between calculations and experiment is not so good for log and half-lives. This is due to large values of the calculated B(GT) strength compared with the experimental one. For Mg(0) Al() transition, the phase space factor, which is estimated to be roughly proportional to (decay energy), depends very much on the interactions as the value for this transition is as small as 1.832 MeV. The excitation energies for the 1 state of Al obtained for the interactions are smaller than the experimental one, =1.373 MeV, by 0.175, 0.571, 0.251 and 0.018 MeV for USDB, IM-SRG, CCEI and CEFT, respectively, which leads to an enhancement of the phase space factor by nearly 10 times for IM-SRG. Though the difference of the B(GT) values is within a range of a factor of about 3, large difference in the phase space factors leads to larger difference in the half-lives.

Here, we make some general comments on the half-lives. (1) For O and F isotopes, calculated half-lives are in fair agreement with the experimental values within a factor of 2.1-2.2, except for O obtained with IM-SRG. (2) The discrepancy between calculated and experimental half-life becomes large (a) when the discrepancy between the calculated and experimental is large, or (b) when the transition with the dominant branching ratio is different between the calculation and the experiment, or (c) when the value for the transition is small and the difference between the calculated and experimental excitation energies is large enough to lead to a substantial change of the phase space factor for the transition. In case of O with IM-SRG, a large discrepancy comes from combined effects of (a) and (b). Nuclei in the island of inversion such as Mg can not be well described for both the and phenomenological interactions due to the reason (a). Mg discussed above corresponds to the case (c). (3) For isotopes with =10-13 (=14-15), there are one or two cases (or more cases) for each in which the calculated half-lives differ from the experimental ones by a factor more than 3 due to the reasons (a), (b) or (c) in case of IM-SRG and CCEI. More extensive study of the half-lives as well as branching ratios and values will be presented elsewhere.

## V Summary and conclusions

In the present work we have performed shell model calculations using approaches along with interaction based on chiral effective field theory and phenomenological USDB interaction, and evaluated , log values and half-lives for the shell nuclei. We also obtained quenching factors corresponding to different interactions for the calculation of the strengths.

All the interactions as well as the interaction based on chiral effective theory considered here need certain quenching of the strengths, as large as by 44-62, as for the phenomenological USDB interaction.
It would be an interesting problem to study how much of them can be explained by the contributions from two-body currents.
The effects of the two-body currents have been studied for transitions in tritium gazit (); baroni () as well as in C and O andres ().
The effects have been studied also for electromagnetic moments and transitions in few-body and light nuclei with 9 Pia (); Pastore ().
The two-body currents have been taken into account in the study of electromagnetic moments and transitions in -, - and -shell nuclei Parzu ()
though they are limited to the induced part.
It is desirable to extend the study of the two-body currents systematically to nuclei with medium and heavy masses.

## Acknowledgments

We would like to thank P. Navrátil, S. R. Stroberg, Jason Holt and G. R. Jansen for useful discussions on effective interactions. We would also like to thank Vikas Kumar and Archana Saxena for useful discussions during this work. AK acknowledges financial support from MHRD for his Ph.D. thesis work. TS would like to thank a support from JSPS KAKENHI under Grant No. 15K05090.

## References

- (1) B. A. Brown and B. H. Wildenthal, “Experimental and theoretical Gamow-Teller beta-decay observables for the -shell nuclei,” At. Data Nucl. Data Tables 33, 347 (1985).
- (2) ENSDF database, http://www.nndc.bnl.gov/ensdf/
- (3) G. Audi, F.G. Kondev, M. Wang, B. Pfeiffer, X. Sun, J. Blachot and M. MacCormick, “The NUBASE2012 evaluation of nuclear properties,” Chinese Phys. C 36, 1157 (2012).
- (4) M. Wang, G. Audi, F.G. Kondev, W.J. Huang, S. Naimi and X. Xu, “The AME2016 atomic mass evaluation,” Chinese Phys. C 41, 030003 (2017).
- (5) A. Ekström, G.R. Jansen, K.A. Wendt, G. Hagen, T. Papenbrock, S. Bacca, B. Carlsson, and D. Gazit, “Effects of Three-Nucleon Forces and Two-Body Currents on Gamow-Teller Strengths,” Phys. Rev. Lett. 113, 262504 (2014).
- (6) P. Maris, J. P. Vary, P. Navrátil, W.E. Ormand, H. Nam, and D. J. Dean, “Origin of the Anomalous Long Lifetime of ,” Phys. Rev. Lett. 106, 202502 (2011).
- (7) P. Navrátil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand, and A. Nogga, “Structure of Nuclei with Two- Plus Three-Nucleon Interactions from Chiral Effective Field Theory,” Phys. Rev. Lett. 99, 042501 (2007).
- (8) D. Gazit, S. Quaglioni, and P. Navrátil, “Three-nucleon low-energy constants from the consistency of interactions and currents in chiral effective field theory,” Phys. Rev. Lett. 103, 102502 (2009).
- (9) A. Baroni, L. Girlanda, A. Kievsky, L. E. Marcucci, R. Schiavilla, and M. Viviani, “Tritium decay in chiral effective field theory,” Phys. Rev. C 94, 024003 (2016).
- (10) P. Klos, A. Carbone, K. Hebeler, J. Menéndez, and A. Schwenk, “Uncertainties in constraining low-energy constants from H decay,” Eur. Phys. J. A 53, 168 (2017).
- (11) B. H. Wildenthal, M. S. Curtin and B. A. Brown, “Predicted features of the beta decay of neutron-rich -shell nuclei,” Phys. Rev. C 28, 1343 (1983).
- (12) W-T Chou, E. K. Warburton and B. A. Brown, “Gamow-Teller beta-decay rates for A18 nuclei,” Phys. Rev. C 47, 163 (1993).
- (13) G. Martínez-Pinedo, A. Poves, E. Caurier and A.P. Zuker, “Effective in the shell,” Phys. Rev. C 53, R2602 (1996).
- (14) A. Arima, K. Shimizu, W. Bentz and H. Hyuga, “Nuclear magnetic properties and Gamow-Teller transitions,” Adv. Nucl. Phys. 18, 1 (1987).
- (15) I. S. Towner, “Quenching of spin matrix elements in nuclei,” Phys. Rep. 155, 263 (1987).
- (16) H. Li and Z. Ren, “Shell model calculations for the -decays of nuclei,” J. Phys. G: Nucl. Part. Phys. 40, 105110 (2013).
- (17) V. Kumar, P.C. Srivastava and H. Li, “Nuclear -decay half-lives for and shell nuclei,” J. Phys. G: Nucl. Part. Phys. 43, 105104 (2016).
- (18) V. Kumar and P.C. Srivastava, “Shell model description of Gamow-Teller strengths in -shell nuclei”, Eur. Phys. J. A : Hadrons and Nuclei 52, 181 (2016).
- (19) S. Yoshida, Y. Utsuno, N. Shimizu and T. Otsuka, “Systematic shell-model study of -decay properties and Gamow-Teller strength distributions in neutron-rich nuclei”, Phys. Rev. C 97, 054321 (2018) .
- (20) J. Suhonen, From Nucleons to Nucleus: Concepts of Microscopic Nuclear Theory (Springer, Berlin, 2007).
- (21) J. Engel, M. Bender, J. Dobaczewski, W. Nazarewicz, and R. Surman, “ decay rates of r-process waiting-point nuclei in a self-consistent approach,” Phys. Rev. C 60, 014302 (1999).
- (22) P Möller, J.R. Nix and K.-L. Kratz, “Nuclear Properties for Astrophysical and Radioactive-ion-beam,” At. Data Nucl. Data Tables 66, 131 (1997).
- (23) P Möller, B. Pfeiffer, and K.-L. Kratz, “New calculations of gross -decay properties for astrophysical applications: Speeding-up the classical process” Phys. Rev. C 67, 055802 (2003) .
- (24) P. Sarriguren, A. Algora, and G. Kiss, ”-decay properties of neutron-rich Ca, Ti, and Cr isotopes”, Phys. Rev. C 98, 024311 (2018) .
- (25) E.A. Coello Pérez, J. Menéndez, and A. Schwenk, “Gamow-Teller and double-beta decays of heavy nuclei within an effective theory,” Phys. Rev. C 98, 045501 (2018).
- (26) A. Saxena, P.C. Srivastava and T. Suzuki, “ calculations of Gamow-Teller strengths in the shell,” Phys. Rev. C 97, 024310 (2018).
- (27) S. R. Stroberg, H. Hergert, J. D. Holt, S. K. Bogner, and A. Schwenk, “Ground and excited states of doubly open-shell nuclei from valence-space Hamiltonians,” Phys. Rev. C 93, 051301(R) (2016).
- (28) G. R. Jansen, J. Engel, G. Hagen, P. Navratil, and A. Signoracci, “ Coupled-Cluster Effective interactions for the Shell Model: Application to Neutron-Rich Oxygen and Carbon Isotopes,” Phys. Rev. Lett. 113, 142502 (2014).
- (29) G. R. Jansen, M. D. Schuster, A. Signoracci, G. Hagen, and P. Navrátil, “Open -shell nuclei from first principles,” Phys. Rev. C 94, 011301(R) (2016).
- (30) L. Huth, V. Durant, J. Simonis, ans A. Schwenk, “Shell-model interactions from chiral effective field theory,” Phys. Rev. C 98, 044301 (2018).
- (31) B. A. Brown and W. A. Richter, “New USD Hamiltonians for the shell,” Phys. Rev. C 74, 034315 (2006).
- (32) B. A. Brown and W. D. M. Rae, “The shell-model code NuShellX@MSU,” Nucl. Data Sheets 120, 115 (2014).
- (33) K. Tsukiyama, S. K. Bogner, and A. Schwenk, “In-Medium Similarity Renormalization Group for Nuclei,” Phys. Rev. Lett. 106, 222502 (2011).
- (34) S. K. Bogner, H. Hergert, J. D. Holt, A. Schwenk, S. Binder, A. Calci, J. Langhammer, and R. Roth, “Nonperturbative Shell-Model Interactions from the In-Medium Similarity Renormalization Group,” Phys. Rev. Lett. 113, 142501 (2014).
- (35) H. Hergert, S. K. Bogner, T.D. Morris, A. Schwenk, and K. Tsukiyama, “The In-Medium Similarity Renormalization Group: A novel ab initio method for nuclei,” Phys. Rep. 621, 165-222 (2016).
- (36) H. Hergert, “In-medium similarity renormalization group for closed and open-shell nuclei,” Phys. Scr. 92, 023002 (2017).
- (37) H. Hergert, S. K. Bogner, J. G. Lietz, T. D. Morris, S. J. Novario, N. M. Parzuchowski, and F. Yuan, Chapter 10, Lecture Notes in Physics (Springer, 2017), ”An advanced course in computational nuclear physics: Bridging the scales from quarks to neutron stars”, edited by M. Hjorth-Jensen, M. P. Lombardo, and U. van Kolck.
- (38) A. F. Lisetskiy, B. R. Barrett, M. K. G. Kruse, P. NavrÃ¡til, I. Stetcu, and J. P. Vary, “ shell model with a core,” Phys. Rev. C 78, 044302 (2008).
- (39) S. Okubo, “Diagonalization of Hamiltonian and Tamm-Dancoff Equation,” Prog. Theor. Phys. 12, 603 (1954).
- (40) K. Suzuki, “Construction of Hermitian Effective Interaction in Nuclei,” Prog. Theor. Phys. 68, 246 (1982).
- (41) P. NavrÃ¡til and B. R. Barrett, “No-core shell-model calculations with starting-energy-independent multivalued effective interactions,” Phys. Rev. C 54, 2986 (1996).
- (42) D. R. Entem and R. Machleidt, “Accurate charge-dependent nucleon-nucleon potential at fourth order of chiral perturbation theory,” Phys. Rev. C 68, 041001 (2003).
- (43) R. Machleidt and D. Entem, “Chiral effective field theory and nuclear forces,” Phys. Rep. 503, 1 (2011).
- (44) P. Navratil, “Local three-nucleon interaction from chiral effective field theory,” Few-Body Syst. 41, 117 (2007).
- (45) E. D. Jurgenson, P. NavrÃ¡til, and R. J. Furnstahl, “Evolution of Nuclear Many-Body Forces with the Similarity Renormalization Group,” Phys. Rev. Lett. 103, 082501 (2009).
- (46) S.âR. Stroberg, A. Calci, H. Hergert, J.âD. Holt, S.âK. Bogner, R. Roth, and A. Schwenk, “Nucleus-Dependent Valence-Space Approach to Nuclear Structure,” Phys. Rev. Lett. 118, 032502 (2017).
- (47) P.C. Srivastava and V. Kumar, “Spectroscopic factor strengths using approaches,” Phys. Rev. C 94, 064306 (2016).
- (48) A. Saxena and P.C. Srivastava, “First-principles results for electromagnetic properties of shell nuclei,” Phys. Rev. C 96, 024316 (2017).
- (49) S. Weinberg, “Effective chiral lagrangians for nucleon-pion interactions and nuclear forces,” Nucl. Physics B 363, 3 (1991).
- (50) A. Piechaczek et al, “ Beta-decay of Mg,” Nucl. Phys. A 584, 509-531 (1995).
- (51) D. H. Wilkinson and B. E. F. Macefield, “A Parametrization of the Phase Space Factor for Allowed -decay,” Nucl. Phys. A 232, 58 (1974).
- (52) A. Sirlin and R. Zucchini, “Accurate Verification of the Conserved-Vector-Current and Standard-Model Predictions,” Phys. Rev. Lett. 57, 1994 (1986).
- (53) D. H. Wilkinson, A. Gallmann and D. E. Alburger, “Super-allowed Fermi decay:Half-lives of O and ,” Phys. Rev. C 18, 401 (1978).
- (54) W. A. Richter, S. Mkhize, and B. A. Brown, “-shell observables for USDA and USDB Hamiltoninas,“ Phys. Rev. C 78, 064302 (2008).
- (55) M. Piarulli, L. Girlanda, L. E. Marcucci, S. Pastore, R. Schiavilla, and M. Viviani, “Electromagnetic structure of = 2 and 3 nuclei in chiral effective field theory,” Phys. Rev. C 87, 014006 (2013).
- (56) S. Pastore, S. C. Pieper, R. Schiavilla, and R. B. Wiringa, “Quantum Monte Carlo calculations of electromagnetic moments and transitions in 9 nuclei with meson-exchange currents derived from chiral effective field theory,” Phys. Rev. C 87, 035503 (2013).
- (57) N. M. Parzuchowski, S. R. Stroberg, P. Navrátil, H. Hergert, and S. K. Bogner, ” electromagnetic observables with the in-medium similarity renormalization group,” Phys. Rev. C 96, 034324 (2017).