arXiv:1507.07442v1 [nucl-th] 27 Jul 2015 ae lc [ place beta-decays and when takes captures phases neutron late between competition at a important duration particularly r-process becomes whole the but during occurs emission tron etn h iesaefrter-process. the for scale numbers, time atomic higher the towards flow setting matter of speed the aflvso h atcptn uli[ nuclei participating the of half-lives ihtmeaue u o1G)adnurndensities neutron and GK) 1 relatively to ( involving (up environments temperatures that explosive accepted high in commonly is site occurs it astrophysical identified, it its been though yet Even not than has r-process. heavier the elements by of synthesis iron the with deals physics h itiuino lmna bnacsi the is abundances elemental on of impact direct distribution a One having the properties chart. nuclear nuclear crucial sev- whole the in the of observables across re- of nuclei number thousand interactions, a eral of weak knowledge between dy- and the interplay complex, quiring electromagnetic a delicate strong, is a the It involving process large required. the input namical to nuclear due of rep- challenge amount r-process difficult the particularly a modeling, resents astrophysical the in volved n n n ftecretyms ciefilso ula astro- nuclear of fields active most currently the of One ag cl vlaino eadcyrtso -rcs nuc r-process of rates beta-decay of evaluation scale Large > 10 2 ntttfu enhsk(hoiznrm,Tcnsh Un Technische (Theoriezentrum), f¨ur Kernphysik Institut 20 ASnmes 11.g 16.z 34.s 26.30.Hj 23.40.-s, 21.60.Jz, 21.10.Tg, numbers: PACS consisten account models. into nucleosynthesis emission taken element neutron heavy be beta-delayed must and and half-lives, beta-decay rate decay total the well general nuclei. Conclusions: short-lived in for are neutron particular half-lives in the nuclei, experimental in transition The nuclei first-forbidden stability. the 5409 of for contribution calculated significant been have neutrons ihntepoo-eto eaiitcqaiatcepha Hartree-B Results: quasiparticle relativistic relativistic the proton-neutron with the calculated within is nuclei all of Method: t asses to rates. calculation reaction nucleosynthesis interaction r-process self-con an fully in a used within calculated transitions, forbidden Background: mato rtfridntastoso ea ae a not indic has are rates Purpose: there decay on though transitions even first-forbidden Additionally, i of making impact input. nuclei this relevant to all rich for els neutron values provides of that half-lives exists beta-decay are important ularly cm g 3 3 S emotznrmfu cwroeeoshn,Planck f¨ur Schwerioneneforschung, Helmholtzzentrum GSI .Priual motn r h beta- the are important Particularly ]. .INTRODUCTION I. 1 − hsc eatet aut fSine nvriyo Zagre of University Science, of Faculty Department, Physics 3 .Aatfo h opeiisin- complexities the from Apart ). The eueaflysl-ossetcvratdniyfunctiona density covariant self-consistent fully a use We opoieatbeof table a provide To β ncranrgoso h ula hr,fis-obde tra first-forbidden chart, nuclear the of regions certain In -rcs uloytei oesrl,b eest,o nu on necessity, by rely, models nucleosynthesis R-process dcyhalf-lives, -decay .Marketin, T. 1 , 2 .I determines It ]. rtfridntransitions first-forbidden β dlydneu- -delayed β 1 β dcyhl-ie and half-lives -decay .Huther, L. dlydnurneiso rbblte n h vrg nu average the and probabilities emission neutron -delayed Dtd uut8 2016) 8, August (Dated: β -decay 2 n .Mart´ınez-Pinedo G. and aflvsapocigtetidrpoespa [ peak r-process for third information the valuable approaching half-lives provided have GSI at periments o -rcs uliaebsdo h uspril random quasiparticle calculations the half-lives on closed based the are near of nuclei most r-process nuclei for reason to this For restricted shells. are shell- calculations cost, computational model increasing the to r-process due for However, half-lives of has around calculation and nuclei the transitions to applied first-forbidden been include [ to model shell extended interacting half- The theoretical lives. on mainly the rely at calculations cleosynthesis nuclei r-process for N half-lives experimental no far numbers neutron magic around N nuclei of half-lives decay N ie.Tu atracmltsaon ai numbers magic around half- accumulates stability N beta-decay matter the Thus larger to substantially lives. closer have moves nuclei flow where consequence matter a As r-process rates. the capture neutron low result- rather numbers in magic ing these around discontinuities show h iiiyo h ai numbers magic the of vicinity the distribution. r-process system solar tpeetol e aflvso -rcs uliin nuclei r-process of half-lives few a only present At vri¨tDrsat 48 amtd,Germany Darmstadt, 64289 iversit¨at Darmstadt, itn irsoi hoeia rmwr.Terslsar results The framework. theoretical microscopic sistent 2 hl lsr r vial.Hne -rcs nu- r-process Hence, available. are closure shell 126 = 2[ 82 = 50 = 50 = eapoiain(pn-RQRPA). approximation se otettldcyrt nnce a rmtevle of valley the from far nuclei in rate decay total the to s β gluo RB oe,adectdsae r obtained are states excited and model, (RHB) ogoliubov icl ots h estvt fncesnhssmod- nucleosynthesis of sensitivity the test to difficult t dlydnurneiso rbblte,icuigfirst- including probabilities, emission neutron -delayed rbblte aeantcal mato h eut of results the on impact noticeable a have probabilities esniiiyo ev lmn uloytei oweak to nucleosynthesis element heavy of sensitivity he uli tpeetol igesseai calculation systematic single a only present At nuclei. ensseaial tde ihnacnitn model. consistent a within studied systematically been erdcd ohfree-vn d- n odd-odd and odd-A even-even, for both reproduced, tosta hi otiuinmyb infiat the significant, be may contribution their that ations , tae1 49 amtd,Germany Darmstadt, 64291 1, straße , 2ad16pouigteosre ek nthe in peaks observed the producing 126 and 82 8 2ad16 h eto eaainenergies separation neutron The 126. and 82 – l nmdr vlain fhl-ie.Bt the Both half-lives. of evaluations modern in tly rc eino h ula hr.W bev a observe We chart. nuclear the of region -rich 11 aebe esrd vntog h ex- the though Even measured. been have ] hoy(DT rmwr.Tegon state ground The framework. (CDFT) theory l N ,100Zge,Croatia Zagreb, 10000 b, 0 2ad16selcoue [ closures shell 126 and 82 50, = la tutr oesfript Partic- input. for models structure clear stoscnttt ag rcinof fraction large a constitute nsitions e ihteicuinof inclusion the with lei ,3 2, 17 N a eetybeen recently has ] 0[ 50 = bro emitted of mber 12 4 – – 7 16 and ] ,so ], 18 ]. e 2 phase approximation on top of semiempirical global mod- and discussion of obtained results. Finally, in Section IV els [19–21], the Hartree-Fock-Bogoliubov [22], or the rel- we provide the concluding remarks. ativistic Hartree-Bogoliubov [23–25]. The quasiparticle random phase approximation (QRPA) approach based on a global effective interaction presents a valid alternative to the interacting shell model. It has been successfully employed in numerous applications [26] and can provide a systematic description of β-decay properties of arbitrar- II. THEORETICAL FORMALISM ily heavy nuclei. The importance of performing calcula- tions based on the self-consistent mean-field models, as A. QRPA calculations opposed to the empirical mean-field potentials has been already emphasized [20]. Currently, however, the only available large-scale calculation of β-decay properties of The calculation of beta-decay rates requires the cal- nuclei used in heavy element nucleosynthesis simulations culation of both the nuclear ground state and excited is based on an QRPA calculation based on a schematic in- states of the daughter nucleus and the transitions be- teraction on top of the Finite Range Droplet Model [19]. tween them, together with the evaluation of the lep- Additionally, in β-decay calculations the impact of first- ton phase space involved in the transition. Because forbidden transitions on the total decay rate is mostly the calculation requires a good description of a wide unexplored. While these transitions and their higher or- range of physical quantities, and because the goal is der counterparts have been included in the investigations to obtain the decay rates for a very large range of nu- of electron and muon capture rates and neutrino-nucleus clei, we employ a fully microscopic theoretical frame- interactions, their use in β-decay studies has been either work based on the relativistic nuclear energy density limited to particular isotopic chains [27] or have been de- functional (RNEDF). The RNEDF based framework em- scribed using a different theoretical foundation such as ploys the self-consistent mean field for nucleons and a the gross theory [19]. Furthermore, due to the existence minimal set of meson fields; the isoscalar scalar σ me- of a single data table of β-decay half-lives, it is not pos- son (J π = 0+,T = 0), the isoscalar vector ω meson sible to perform meaningful sensitivity studies on the in- (J π = 1−,T = 0) and the isovector vector ρ meson fluence of beta-decay rates on r-process nucleosynthesis. (J π = 1−,T = 1), supplemented with the electromag- netic field. The meson-nucleon interaction is included In an effort to avoid effects coming from combining with a minimal set of the interaction terms, where the models built on differing theoretical bases, and with the vertex functionals include explicit dependence on the aim of providing a new, high precision data table of β- nucleon vector density [28]. The nuclear ground state decay properties, in this study we utilize a fully self- properties are described using the relativistic Hartree- consistent theoretical framework based on the relativistic Bogoliubov model (RHB), which properly describes the Hartree-Bogoliubov (RHB) model [28] for the descrip- pairing effects in open-shell nuclei. For the model param- tion of the ground state of open- and closed-shell nuclei eters that determine the density-dependent couplings and with the proton-neutron relativistic QRPA (pn-RQRPA) ∗ the meson masses, in this work the D3C parametrization where the residual interaction is derived from the same is used [24], as it was previously shown that it provides density functional as was used for the ground state cal- a good description of β-decay half-lives in medium and culations. In this way, we ensure the consistency among heavy nuclei. The pairing correlations in open shell nu- the various required properties, from the masses and Q β clei are described by the finite range Gogny interaction, values to the nuclear response strength functions. The with the D1S parametrization [33]. RHB + pn-RQRPA was already successfully employed in the study of Gamow-Teller and higher order reso- The excited states are obtained using the proton- nances [29, 30], β-decay half-lives [23, 24] and stellar neutron relativistic quasiparticle random phase approxi- weak-interaction processes [31, 32]. This framework also mation (pn-RQRPA), formulated in the canonical single- enables the treatment of first-forbidden transitions on an nucleon basis of the RHB model [34] and extended equal footing with the Gamow-Teller transitions, allow- to the description of charge-exchange excitations (pn- ing for a meaningful investigation of the contribution of RQRPA) [29]. The RHB + RQRPA model is fully self- the J > 0 modes of nuclear excitation and for studying consistent in both the ph and pp channels. The same the regions of interest where this contribution is signif- interactions are used in the RHB equations that deter- icant. We use this model to calculate the β-decay half- mine the canonical quasiparticle basis, and in the matrix lives and the β-delayed neutron emission probabilities of equations of the RQRPA. Transitions between the 0+ 5409 neutron-rich nuclei with 8 Z 124, including ground state of a spherical parent nucleus and the J π even-nuclei but also nuclei with an≤ odd≤ number of parti- excited state of the corresponding daughter nucleus are cles. In Section II we present the theoretical formalism induced by a charge-exchange operator T JM . Assuming of the model, while Section III contains all the neces- spherical symmetry of the nuclear system, the quasipar- sary expressions for the evaluation of matrix elements ticle pairs can be coupled to good angular momentum and nuclear properties, together with the presentation and the matrix equations of the pn-RQRPA read [35]: 3

with values VL = 160.0 MeV, VD = 15.0 MeV, a = 7.2 and b = 0.3 adjusted to obtain the best possible de- AJ BJ XλJ 1 0 XλJ − ∗ ∗ = E , scription of available half-life data [36]. B J A J Y λJ λ 0 1 Y λJ      −    25 (1) where the matrices A and B are defined in the canonical 20 Ag λJ λJ basis [35]. For each energy Eλ, X and Y in Eq. (1) 15 [MeV] denote the corresponding forward- and backward-going - D3C*

β 10 FRDM QRPA amplitudes, respectively. The transition matrix Q HFB-21 5 element between the ground state of the parent nucleus exp. and the excited state of daughter nucleus, induced by the 0 operator T JM , reads 20 Cd 15 2 [MeV] - ± J λJ J λJ β 10 BλJ = p T n Xpn upvn + ( 1) Ypn vpun . Q h || || i − 5 pn X
(2) 0 64 68 72 76 80 84 88 92 96 100 104 108 where the X and Y amplitudes are obtained from solving neutron number the pn-RQPA equation (1). In the T = 1 channel of the pn-RQRPA we use the FIG. 1. (color online) β-decay Q-values for the silver and cad- pairing part of the Gogny force: mium isotopic chains. Values calculated in this work, marked with red squares, are compared with the values obtained from − 2 2 the FRDM mass table [37], denoted with blue diamonds, the V pp(1, 2) = e r12/µi (W + B P σ H P τ i i − i HFB-21 model [38] and with existing experimental data [36] i=1,2 denoted with black circles. X σ τ MiP P ), (3) − The model employed in this study is well suited to σ τ with r12 = r1 r2 , P , P the spin and isospin ex- describe the properties of the ground state and of the | − | change operators and the D1S set of parameters µi, Wi, excited states of even-even nuclei. However, as the r- Bi, Hi and Mi (i = 1, 2) taken from [33] . This force process naturally includes odd-A and odd-odd nuclei as has been carefully adjusted to pairing properties of fi- well, it becomes necessary to provide at least an approx- nite nuclei all over the periodic table. In particular, the imate description of the decay properties of these nuclei. basic advantage of the Gogny force is the finite range, With this aim, we compute the ground state of odd nuclei which automatically guarantees a proper cut-off in the by employing the same model as for even-even nuclei, but momentum space. For the T = 0 proton-neutron pairing constraining the expectation value of the particle number interaction in open shell nuclei we use a form consist- operator to an odd number of protons and/or neutrons. ing of a short-range repulsive Gaussian combined with a In this way an even RHB state is obtained, whose energy weaker longer-range attractive Gaussian. is different from true odd nucleus ground state energy by energy of the odd quasi-particle [39, 40]. To quantify 2 2 − r12 µ2 the validity of this approximation we show the computed V = V g e j Πˆ . (4) 12 − 0 j S=1,T =0 Qβ values for the silver and cadmium (Z = 45) isotopic j=1 X chains in Fig. 1, calculated by approximating the binding energy of the daughter nucleus [22]: where Πˆ S=1,T =0 projects onto states with S = 1 and T = 0. The ranges µ1 = 1.2 fm and µ2 = 0.7 fm of the two Gaussians are taken from the Gogny interaction. Q = ∆MnH + λn λp E2qp − − The relative strengths g1 = 1 and g2 = 2 are chosen = Mex(Z,N) Mex(Z +1,N 1), (6) so that the force is repulsive at small distances.− The − − only remaining free parameter is V0, the overall strength. where λn and λp are the neutron and proton Fermi ener- This interaction, with a constant value of V0 was used gies in the ground state of the parent nucleus, E2qp are in the non-relativistic QRPA calculation [22] of β-decay the lowest quasiparticle energies. ∆MnH = 0.782 MeV rates for spherical neutron-rich r-process waiting-point is the mass difference between a neutron and the hydro- nuclei. Two relativistic calculations of β-decay half-lives gen atom, and Mex(Z,N) is the mass excess. In silver of neutron-rich nuclei [23, 24] have shown that a single isotopes (odd Z), the Qβ values are well reproduced for value of the overall interaction strength cannot be suc- nuclei with an even number of neutrons, and only slightly cessfully used in different regions of the nuclear chart. underestimated for odd-odd nuclei. The predicted values Thus, we take the ansatz proposed in [25] closely follow the jump above the closed neutron shell, and begin to deviate from the FRDM values after the V shell closure, while the HFB-21 results are found be- V = V + D , (5) 0 L 1+ ea+b(N−Z) tween the other two models. For the cadmium isotopes 4 the model overestimates the Qβ values in nuclei with an where the factors k, ka, kb and kc are defined as [45]: even number of neutrons in lighter isotopes, but above 1 1 4 N = 78 the results reproduce the experimental values ac- k = ζ2 + w2 + ζ2 + (x + u)2 µ γ u(x + u) 0 9 1 9 − 9 1 1 curately. Again, the model very successfully reproduces  (0)  the changes around the closed neutron shell. For heavier 1 2 2 1 2 isotopes there are no available data, but all the models + W0 (2x + u) λ2(2x u) (11a) 18 − 18 − (1) coincide in their predictions of the Qβ values. We can  1 conclude from these results that this approach to esti- + z2 W 2 λ , mating the Q-values provides a good agreement with the 12 0 − 2  (2) available data and does not hinder our study of weak- 4 1
1 interaction processes. ka = uY W 4x2 +5u2 W z2 (11b), −3 − 9 0 − 6 0  (1)  (2) 2
kb = µ γ [ζ w] + [ζ (x + u)] , (11c) 3 1 1 − 0 (0) 1 (1) 1 n o B. Beta-decay half-lives kc = 8u2 + (2x + u)2 + λ (2x u)2 (11d) 18 2 − (1) 1  + (1 + λ ) z2 , In the present study of beta-decay half-lives we include 12 2 (2) both allowed (L = 0) and first-forbidden (L =0, 1) tran- with   sitions. The beta-decay rate for a transition between an ′ ′ 1 initial and final nuclear state is equal to [41] V = ξ v + ξw , ζ0 = V + wW0, 3 (12) ′ ′ ′ 1 Y = ξ y ξ(u + x ), ζ1 = Y + 3 (u x)W0. ln 2 p0 − − λ = p2 (W W )2 F (Z, W )C(W )dp , (7) K e 0 − e The numbers in parenthesis after the closing brackets Z0 denote the rank of the operators inside. The parameter 2 γ1 is given by 1 (αZ) . For the coulomb wave func- 2 − where W is the electron energy in units of mec with tions we use the approximations µ 1 and λ 1 [46]. p 1 ≈ 2 ≈ W0 being the maximum electron energy that is equal The quantity ξ = αZ/(2R) where we choose R = r 2. to the difference in nuclear masses between initial and In the Condon-Shortley phase convention [43] theh ma-i final nuclear states, W = (M M )/m , and p is the − p 0 i − f e e trix elements for β transitions are: electron momentum in units of mec. We approximate the Ck σk 0 tk maximum electron energy with [22]: f k rk 1 − i w = g √3 ⊗ , (13a) A D E − P √2Ji +1  Mi Mf λn λp + ∆mnp EQRPA. (8) k k f rkC t− i − ≈ − − x = k 1 , (13b) − √ P 2Ji +1 The constant K is measured in superallowed beta-decay Ck σk 1 tk to be, K = 6144 2 s [42]. f k rk 1 − i ± u = g √2 ⊗ , (13c) A D E The shape factor C(W ) differs for various decays. For − P √2Ji +1  − allowed decays it is energy independent. In the case of β Ck σk 2 tk f k rk 1 − i decay of neutron-rich nuclei the shape factor is simply the z =2g ⊗ , (13d) A D E Gamow-Teller reduced transition probability, C(W ) = P √2Ji +1  B(GT ), with k k 0 k f r I(1, 1, 1, 1, r ) C σ t− i ′ 2 k k k 1 ⊗ w = gA (13e), 2 − √3 D P √2Ji +1  E f σktk i 2 k − Cktk B(GT )= gA . (9) ′ 2 f k rkI(1, 1, 1, 1, rk) 1 − i (2Ji + 1) x = , (13f) P −3 √ P 2Ji +1 k k 1 k f r I(1, 1, 1, 1, r ) C σ t− i The matrix element is reduced with respect to the spin ′ 2√2 k k k 1 u = g ⊗ (13g). operator σ only using the Condon-Shortley phase convec- A D E − 3 P √2Ji +1  tion [43]. The sum runs over all nucleons. For the isospin lowering operator we use the convention t− n = p . Fi- The matrix elements connected with the relativistic cor- | i | i nally, gA = 1.2701(25) [44] is the weak axial coupling rections are − k k constant. ′ f γ t− i ξ v = g k 5 , (14a) For first-forbidden transitions the shape factor reads: − A √ P2Ji +1 k k f α t− i ξ′y = , (14b) C(W )= k + kaW + kb/W + kcW 2. (10) − √ 2Ji +1

5

4 α 10 with γ5 and being the Dirac matrices. Due to the fact 3 10 exp. (a) (b) that our formalism is relativistic we do not perform a 2 D3C* 10 FRDM non-relativistic reduction in the evaluation of the matrix 1

[s] 10 0

C 1/2 elements in equations (14). The quantity LM equals 10 T -1 10 Sr Ag -2 4π 10 C 4 LM = YLM , (15) 10 2L +1 3 r 10 (c) (d) 2 10 where YLM are the spherical harmonics 1

[s] 10 As the matrix elements are independent of the elec- 0 1/2 10 T tron and neutrino energies, the integrals over the elec- -1 10 Y Cd -2 tron phase space are evaluated independently and appear 10 -3 10 only as multiplicative factors. The function I(1, 1, 1, 1, r) 52 56 60 64 68 72 76 64 68 72 76 80 84 88 92 takes into account the nuclear charge distribution, and in neutron number neutron number the approximation of the uniform spherical distribution it has the form [45] FIG. 2. (color online) Comparison of the calculated half- lives with the FRDM results and the experimental data for 1 r 2 the strontium (upper left), yttrium (lower left), silver (upper 3 1 5 R 0 r R I(1, 1, 1, 1, r)= − 3 ≤ ≤ . (16) right) and cadmium (lower right) isotopic chains. In cases 2 R 1 R r > R where they are not visible, the error bars are smaller than the ( r 5 r
− data marker. Systematic calculations of beta-decay
half-lives deter- mined by Gamow-Teller transitions have shown that the 100 Sr D3C* (a)Ag (b) theoretical matrix elements need to be quenched by a 80 FRDM factor q that has been found to be independent of the 60 (%) tot

particular transition considered and approximately con- λ

/ 40 ff stant over the nuclear chart [47, 48]. There is also ev- λ idence that the total first-forbidden transition strength 20 is overestimated by theoretical approaches. In addi- 0 tion, the quenching factor seems to depend on the model 80 Y (c)Cd (d) space used and on the particular first-forbidden opera- 60 (%) tot

tor [18, 49–53]. So far no study has addressed the quench- λ

/ 40 ff ing of first-forbidden transitions of global calculations of λ beta-decay half-lives. In our calculations, we account 20 0 for the quenching of both the Gamow-Teller and first- 52 56 60 64 68 72 76 80 84 88 72 76 80 84 88 92 96 forbidden transitions by using the same effective value of neutron number neutron number gA = 1.0. − FIG. 3. (color online) Percentage of first-forbidden transitions contributing to the decay of Sr, Y, Ag and Cd isotopes. In the cases where it was not possible to extract the value for a III. RESULTS particular model, the data point was not plotted.

A. β-decay half-lives Q-values, however for isotopes more neutron-rich than In Fig. 2 we compare the calculated half-lives for four N = 58 we obtain very good agreement with the data, isotopic chains, including both the even and odd atomic with a slight tendency to underestimate the half-lives of numbers, where the half-lives are obtained from the decay yttrium isotopes. The FRDM+QRPA approach displays rate in Eq. (7): erratic jumps in the half-lives and odd-even staggering that is not present in the data. For very neutron-rich ln 2 T1/2 = . (17) isotopes the two models provide very similar results due λ to the increase of Q-beta values. The right panels compares the present calculation and In the four panels we show a comparison of the present the the FRDM+QRPA model with data for nuclei close results with the data and the FRDM results for the stron- to the Z = 50 closed proton shell, including the silver and tium (Sr, Z = 38), yttrium (Y, Z = 39), silver (Ag, cadmium isotopic chains. Similar to the two lighter iso- Z = 47), and cadmium (Cd, Z = 48) isotopic chains. topic chains, the D3C* interaction provides an excellent The two lightest isotopic chains, strontium and yttrium, description of the decay properties of neutron-rich nu- are found between two closed proton shells Z = 28 and clei, with the results correctly following the trend along Z = 50. Our calculation has difficulties reproducing the the chain. In the case of silver, the present work re- half-lives of the longest lived nuclei due to the rather low produces the trend quite well in comparison with the 6

FRDM. Most likely, because the FRDM reproduces the simulations. We will discuss this aspect in more detail in half-life of 114Ag it cannot provide a good trend for most the next section. of the isotopic chain. The half-life of this nucleus is un- It is also possible to compare the present results with characteristically short compared with the surrounding the newest data available to assess the quality of the isotopes due to the J π = 1+ angular momentum and extrapolation of half-lives into currently inaccessible re- parity of the ground state which quickly decays into the gions. Fig. 4 compares our calculations with the re- 0+ ground state of 114Cd. For cadmium nuclei both mod- els do quite well, although the FRDM shows unexpected 10 changes in the half-lives along the chain, as in the cases D3C* of other elements. FRDM In the following figure, Fig. 3, we plot a comparison of the contribution of the first-forbidden transitions in exp. the total decay rate, in percent, between the two models. 1/2 For the FRDM, we have extracted the values by com- /T 1 calc. paring the results of Ref. [54] and Ref. [19]. Due to the 1/2 changes in the details of the calculation between the two T publications, the values presented here are not precise. We expect them to be a good qualitative measure of the actual values. The FRDM calculation predicts, in two 0.1 out of four isotopic chains, a large contribution (more 96 98 100 102 104 106 108 110 112 114 116 118 than 50%) of first-forbidden transitions in nuclei close to number of nucleons stability, and a decrease in more neutron-rich nuclei. For FIG. 4. (color online) The ratio of the calculated and mea- nuclei in the isotopic chains of yttrium and silver, FRDM sured half-lives for nuclei in the 36 ≤ Z ≤ 43 isotopic chains predicts a steady, small contribution below 20% for most from Ref. [9], for D3C* interaction (full line) and the FRDM of the nuclei, with odd-even staggering. (dashed line). Different colors denote the isotopic chains, The present results predict a rather small contribution while the thick black lines denote deviations by a factor 2. of the forbidden transitions in nuclei close to the stability, with a smooth increase with additional neutrons. This cently measured beta-decay half-lives [9] for a number of increase is particularly visible at the N = 82 neutron isotopic chains with atomic number in the range 36 shell closures for Ag and Cd isotopes, where the neutrons Z 43. The results obtained with the D3C* interaction≤ ≤ begin to occupy the 2f7/2 and the 1h9/2 orbits, enabling are in agreement with the experiment within a factor 2. the strong ν2f7/2 π1g7/2 and ν1h9/2 π1g9/2 transi- The FRDM+QRPA approach, on the other hand, tends → → tions, which significantly increase the contribution of the to predict too long half-lives and displays jumps in the first-forbidden transitions as these transitions appear at ratio for specific nuclei. low excitation energies in the daughter nucleus. In the two lighter isotopic chains, the contribution is negligi- 1000 ble for the lightest isotopes and increases smoothly up exp. D3C* to 30% where it reaches a plateau. The main differ- FRDM ≈ 94 114 ence in the transitions between e.g. Sr and Sr is the 100 occupied neutron 2d3/2 state, which contributes to the total rate by decaying to the proton 2p states. Above the N = 82 closed neutron shell, the contribution of the [s] 10 1/2 first-forbidden transitions again rises with the occupation T of the 2f7/2 and the 1h9/2 states, just as in the silver and cadmium isotopes. 1 As discussed above, our calculations are only strictly applicable to even-even nuclei. However, we have ex- 0.1 tended them to odd-A and odd-odd nuclei in order to 194 Re 196 Re 200 Os 199 Ir 202 Ir 204 Pt 211 Tl 213 Tl 219 Bi provide a complete coverage for r-process simulations. 195 Re 199 Os 198 Ir 201 Ir 203 Pt 204 Au 212 Tl 218 Bi Our approach cannot account for mismatches between the angular momenta of the ground states of parent and FIG. 5. (color online) Comparison of the data (denoted with daughter nucleus in the decay. However, this aspect be- black circles) with the results of the present work (red squares) comes less and less important as the Q-value for the decay and the FRDM (blue diamonds). increases. This is one of the reasons why our approach becomes better as we move far from stability. Our predic- As we move to heavier nuclei both models have signifi- tions for half-lives are certainly not worse that those com- cant difficulties while attempting to reproduce the recent puted by the FRDM+QRPA approach [19] that is cur- data. Fig. 5 compares the two models with measurements rently the standard model for half-lives used in r-process [12, 13]. FRDM significantly overestimates the half-lives 7 of most of the nuclei, by more than an order of mag- with the very recent measurement of half-lives of 110 nitude. It is only for nuclei around the Z = 82 closed neutron-rich isotopes, 40 of which have not been previ- proton shell that provides a reasonable description of the ously measured. In the upper panel we plot the results of data. An interesting effect occurs in the case of rhenium the present study, from the isotopes of rubidium to tin, isotopes, where the half-life of 194Re is overestimated, and from N = 65 to N = 89. Almost all of the values but half-lives of the two heavier isotopes are reproduced. can be found within a factor of 2 from the data, with The discrepancy seems to be related to the change of pro- very little scatter for different isotopic chains, except for late to oblate deformation in moving from 194Re to the the case of tin where we underestimate the half-lives by heavier isotopes [37]. With the exception of 194Re the a factor of 5. Within a particular chain the data is re- D3C* calculations reproduce all beta-decay half-lives of produced smoothly, with very weak odd-even staggering. the Z 79 nuclei shown in figure 5 within an order of In the case of FRDM the scatter is much larger, and we magnitude.≤ Larger differences appear for heavier nuclei. also observe some regularities. The half-lives of nuclei in Taking into account that the nuclei shown in figure 5 have isotopic chains with an even atomic number are system- half-lives in the range 1 s to 100 s one should conclude atically overestimated by a factor of two or more - these that both FRDM and D3C* approaches provide a fair are the Zr, Mo, Ru, Pd, Cd and Sn isotopes. reproduction of data. Certainly, the agreement is worse that the one obtained in figure 4 but this should not be consider as a drawback of the two models discussed to re- produce half-lives of heavier nuclei. It simply reflects the fact that for heavy nuclei the half-lives are only known for relatively long lived nuclei close to the stability that normally have small Q-beta decay windows. The theo- retical description of these decays is a challenge for any theoretical model as the decay rate depends on very few selected transitions.

10

Rb D3C* Sr Y Zr Nb Mo Tc Ru exp. Rh

/T 1 Pd Ag calc.

T Cd In Sn

0.1 64 68 72 76 80 84 88 neutron number 10

Rb FRDM Sr Y Zr Nb Mo Tc Ru exp. Rh

/T 1 Pd Ag calc.

T Cd In FIG. 7. (color online) The logarithm of the ratio of the calcu- Sn lated and experimental half-lives is plotted for all neutron-rich nuclei with available data, with half-lives shorter than 1 hour. Black squares indicate that the model did not provide a finite

0.1 half-life. In the top panel we plot the results of the present 64 68 72 76 80 84 88 neutron number study, while in the bottom panel we plot the ratios obtained with the FRDM. FIG. 6. (color online) The ratio of the calculated and experi- mental half-lives for the recently measured 110 isotopes [55]. Fig. 7 shows the ratio of the calculated and measured In the top panel we present the results obtained in the present half-lives for neutron-rich nuclei with half-lives shorter study, while in the bottom panel we present the FRDM re- than 1 hour. We show the results for the FRDM+QRPA sults [19]. approach and our own calculations (bottom and top pan- els, respectively). The present calculation has difficulties Finally, in Fig. 6 we compare the theoretical results in reproducing the half-lives for nuclei close to simultane- 8 ous neutron and proton shell closure. In particular, the greatest discrepancies appear in the regions “north-west” of the doubly closed nuclei 78Ni and 132Sn. The amount of low energy GT strength for such nuclei is very sensitive to correlations [18] that are not captured by the QRPA approach. For an improved description of the low-energy strength distributions of these nuclei, complex configura- tions need to be included to shift the transition strength from the GT resonance to low energies [56, 57]. Cou- pling to collective phonon excitations is known to have a strong impact on the single particle structure around shell closures [58], and will enhance the decay rate of these nuclei [59]. There are also three regions where the D3C* calcula- tions predicts shorter decay half-lives, i.e. the regions around the 90Kr (Z = 36, N = 54) and 140Ba (Z = 56, N = 84) nuclei, and the Z . 82 N > 126 region. The first two cases correspond to regions of low defor- mation [60]. However, there are indications that nuclei FIG. 8. (color online) The ratio of the calculated and the on that region are gamma soft [61, 62]. So far the im- experimental half-lives [36] versus the experimental half-life. pact of gamma-softness on beta-decay rates has not been Different types of nuclei are denoted with different symbols. explored. It is interesting that the FRDM+QRPA ap- Even-even nuclei are denoted with blue squares, odd-N with proach also predicts shorter half-lives for these two re- orange triangles, odd-Z with black diamonds, and odd-odd gions. In the case of the region Z . 82 N > 126, we nuclei are denoted with red circles. are dealing with nuclei with half-lives longer than 100 s. The decay rate is determined by few transitions that depend very sensitively on the predicted single-particle calculated and experimental half-lives structure. In particular, for Hg and Tl isotopes, above calc. the N = 126 closed neutron shell the decay is domi- T1/2 ri = log10 exp. , (18) nated by the Gamow-Teller ν1i11/2 π1i13/2 and the T1/2 − → 1 ν1i11/2 π1h9/2 configurations. These transitions are not possible→ for isotopes with N 126, but above with the average value of the ratio providing a measure ≤ the neutron shell closure the νi11/2 is becoming signifi- of the global deviation from experimental values cantly occupied (occupation probability is already 0.2 210 211 ≈ N for Hg and Tl) at an energy of Ei 2.5 MeV. 1 11/2 ≈− r¯ = r , (19) N i Fig. 8 shows the ratio of the calculated and the exper- i=1 imental half-lives versus the experimental half-life. Dif- X ferent symbols are used for even-even, the odd-N, odd-Z while the standard deviation provides a measure of the and the odd-odd nuclei. The data is taken from Ref. [36]. spread [63] For long lived nuclei we observe a very large spread of the 1/2 ratios due to the high sensitivity of the decay rate to the N 1 2 details of the low energy transition strength. Short-lived σ = (ri r¯) . (20) N − " i=1 # nuclei have large Q-values and generally several transi- X tions contribute to the decay rate. This makes these nu- The average deviation and the spread obtained using the clei less sensitive to the detailed position and strength of ∗ D3C and the FRDM model [19] are shown. The top particular transitions. Our results show that the ratio of part of the table provides data for nuclei with the ex- half-lives converges towards one with decreasing half-live, perimental half-life shorter than a particular value. The for even-even, odd-A and odd-odd nuclei equally. This average value of deviation from experimentr ¯ is compara- behaviour points to the reliability of the model for all ble between the two models. There are small differences types of nuclei taking part in the r-process, as almost all for particular sets of nuclei, but they are not significant. of the predictions lie within an order of magnitude from ∗ The resulting standard deviation is larger with the D3C the data, with a large majority within a factor of 2. for long-lived nuclei, confirming that the description of In Table I we provide a quantitative analysis of the nuclei with small Qβ values is challenging for the model. results, by using several measures of the model precision. For nuclei with half-lives shorter than 1 s, the models are Because both the observable (half-life) and the ratio of comparable both in the average and the standard devia- the calculated and experimental values span orders of tion, and for even less stable nuclei D3C∗ even manages magnitude we define ri as the logarithm of the ratio of to provide a smaller spread of ratios. 9

TABLE I. Comparison of the average deviations (Eq. (19)) and their standard deviations (Eq. (20)) between the model used in the present work and the FRDM. In the upper part of the table, we show the results structured with respect to the experimental half-life of the nuclei included, while in the lower part of the table we show the results for nuclei with half-lives shorter than 1 s, but separated into groups of even-even, odd-Z, odd-N and odd-odd nuclei. ∗ D3C FRDM r¯ σ r¯ σ Texp. < 1000 s 0.011 0.889 0.021 0.660 < 100 s 0.057 0.791 0.040 0.580 < 10 s 0.061 0.645 0.046 0.515 < 1 s 0.011 0.436 0.019 0.409 < 0.1 s 0.041 0.195 0.021 0.354 nucleus type even-even −0.037 0.331 0.333 0.226 odd-Z 0.054 0.328 −0.128 0.288 odd-N −0.086 0.387 0.124 0.436 odd-odd 0.089 0.582 −0.179 0.409 total 0.011 0.436 0.019 0.409

In the bottom part of the table we provide the results, 50%, respectively. The standard deviations are compa- including nuclei with the experimental half-life lower or rable between the two models for odd-A nuclei, but are equal to 1 s, but shown for each type of nuclei, i.e. even- smaller for even-even and especially for odd-odd nuclei. even, odd-Z, odd-N and odd-odd, separately. From this Still, for both models the odd-N and odd-odd nuclei are point of view, the results of the two models are sig- the most difficult to describe and result with the largest nificantly different. In the present work, the half-lives spreads. So, even though the averager ¯ of the FRDM of even-even and odd-N nuclei are somewhat underesti- model is close to zero, it is actually a result of the can- mated by roughly 10% - 20%, and the half-lives of odd-Z cellation of overestimated and underestimated half-lives and odd-odd nuclei are, on average, overestimated by in different types of nuclei. This behaviour has also ap- 15% and 20%, respectively. These results agree with the peared in previous evaluations of the β-decay rates (see assumption that the model is capable of describing all Ref. [54], Table B.), where the half-lives of short-lived types of nuclei equally well, as the half-lives of odd-A even-even nuclei were, on average, overestimated by al- and odd-odd nuclei are reproduced as well as half-lives most a factor of 4, while the half-lives of odd-A and odd- of even-even nuclei. The difference between various types odd nuclei were underestimated by a factor of 2. In com- of nuclei becomes more evident by looking at the stan- parison, the results of the present work are more con- dard deviations. Here, the spread of values for even-even sistently close to reproducing the data, even though the and odd-A nuclei is comparable, but the standard devi- larger standard deviations highlight some unresolved is- ation of the results for odd-odd nuclei is almost twice as sues. large, pointing to potential angular momentum issues in the ground state of an odd-odd nucleus. It is interesting that our model does particularly well for even-N nuclei that are the most relevant nuclei for r-process nucleosyn- thesis as these nuclei are favored by neutron-captures. The FRDM+QRPA half-lives show a distinctly dif- ferent behaviour. As opposed to the D3C∗ model, the FRDM overestimates the half-lives of even-even and odd- N nuclei, by a factor of 2.15 and 1.33, respectively. That the predicted half-lives of even-even nuclei are on average twice the experimental value is surprising as, similarly to our calculation, the model should be more suitable for these nuclei. This fact may have important consequences for r-process nucleosynthesis and is might be responsible FIG. 9. (color online) The ratio of halflives obtained by the for the differences in r-process nucleosynthesis discussed two models. on section IIID. However, waiting point nuclei also in- clude odd-Z, even-N nuclei, and in both models there is Finally, in Fig. 9 we plot the ratio of half-lives obtained a cancellation between even-even and odd-Z nuclei. in the present study with the FRDM half-lives [19]. For On the other hand, the half-lives of odd-Z and odd- light and medium-heavy nuclei, the results are quite simi- odd nuclei are underestimated by approximately 30% and lar. Close to the valley of stability, the D3C∗ model tends 10

2 to provide longer half-lives than the FRDM for particu- 10 lar nuclei, especially so in the regions “north-west” of 1 the doubly closed nuclei 78Ni and 132Sn. Further from 10 0 the valley of stability the two models provide compara- 10 ble results, with the difference that the present calcula- [s]

1/2 -1

tion predicts smoother increase in the decay rates. For T 10 very heavy nuclei with N > 126 however, the present -2 (a) study predicts half-lives to be shorted by more than an 10 order of magnitude than the FRDM predictions. This is -3 10 especially clear in nuclei with high atomic numbers, i.e. 24 26 28 30 32 1 Z ' 95 where the difference is as large as three orders of 10 magnitude. This may have significant consequences on Shell model FRDM (GT + ff) the dynamics of the r-process nucleosynthesis in neutron 0 star merger conditions [64] 10 D3C* + RQRPA

[s] DF3 + QRPA exp.

1/2 -1

T 10 B. Impact of First-forbidden transitions (b) -2 10 The results presented in this manuscript are the first 42 44 46 48 50 global calculations of beta-decay half-lives for neutron- 1 rich nuclei that treat Gamow-Teller and first-forbidden 10 transitions on an equal footing. In this section, we pro- 0 vide a detailed analysis on the impact of first-forbidden 10 transitions with particular focus on nuclei with magic

[s] -1 neutron numbers N = 50, 82 and 126 that have been the 10 1/2

subject of many theoretical calculations. Additionally, T -2 for the N = 50 and 82 experimental data is available 10 that has contributed to constrain the calculations. (c) -3 Fig. 10 compares the beta-decay half-lives for r-process 10 nuclei with neutron magic number N = 50, 82 and 126 66 68 70 72 74 with the FRDM+QRPA model [19], the shell-model cal- proton number culations of ref. [18] and data [36]. For the N = 50 iso- tones (upper panel) with Z 28 the present calculations FIG. 10. (color online) Beta-decay half-lives for N = 50 fail to reproduce the measured≥ half-lives. In these nuclei (a), N = 82 (b), and N = 126 (c) isotopic chains. The present results are compared with those of FRDM [19], the the f proton shell is fully occupied suppressing low en- 7/2 shell-model [18], DF3 [27] and data [36]. ergy Gamow-Teller transitions. The decay consequently proceeds mainly by first-forbidden transitions resulting in long decay half-lives. As a consequence the model predicts a large contribution of forbidden transitions to tions for Z = 46, 47, 48 and similar values for the other the decay rate (see fig. 11). This is probably a limitation isotones. The odd-even staggering present in the beta- of the QRPA type calculations that are not able to pro- decay half-lives is also manifested on the contribution duce enough correlations around the proton magic num- of forbidden transitions. For N = 126 all models pre- ber Z = 28. This problem does not seem to be present dict an enhanced contribution of forbidden transitions. in the FRDM+QRPA approach. However, one should Ref. [18] argues that the contribution of forbidden tran- kept in mind that in this approach the Gamow-Teller sition should decrease with decreasing proton number. and first-forbidden contributions are not derived based This reduces the number of protons in the h11/2 increas- on the same microscopic model. As discussed above the ing the role of νh9/2 πh11/2 Gamow-Teller transitions. → FRDM+QRPA approach tends to produce strong odd- While all models predict such a reduction, the magnitude even effects that are not present in the data nor in the varies from model to model. The present calculations other approaches. predict a relatively minor reduction from 90% at Z = 73 The situation is different for the N = 82 and N = 126 to 70% at Z = 66, while the shell-model predicts a much isotonic chains. In these nuclei there is no proton shell larger reduction reaching below 20% at Z = 66. Exper- closure that hinders low energy Gamow-Teller transi- imental data will be very important to further constrain tions and in principle both Gamow-Teller and forbidden the theoretical calculations. transitions are possible. The present study predicts a Fig 12 shows the contribution of beta-decay half-lives very similar contribution of first-forbidden transitions to for all nuclei with predicted beta-decay half-lives shorter the shell-model, 20%, for all the N = 82 isotones. than 1 s. In general, first-forbidden transitions do not FRDM+QRPA predicts∼ rather large forbidden contribu- noticeably impact the total decay rate in nuclei close to 11

100 120 110 100 100 80 80 (a) 90 80 60 60 (%) 70 40

tot 60

λ 20 / 40 50 ff

λ 40 0 λff 20 proton number 30 100 λ +λ 20 GT ff 0 10 24 26 28 30 32 0 0 20 40 60 80 100 120 140 160 180 200 220 240 100 neutron number 80 Shell model (b) DF3 + QRPA FIG. 12. (color online) The contribution of the first-forbidden FRDM + QRPA 60 transitions in the total decay rate for nuclei with predicted (%) D3C* + RQRPA half-lives shorter than 1 s. Black lines indicate the position tot

λ of (predicted) closed proton and neutron shells. / 40 ff λ 20 ity and consequently Gamow-Teller and forbidden tran- 0 42 44 46 48 50 sitions contribute in similar amounts to the decay. The 100 situation is different in the region between N = 50 and N = 82. Here, once protons completely fill the pf shell 80 (Z > 40) forbidden transitions are suppressed and the decay is dominated by Gamow-Teller transitions. 60 (%) (c) An interesting phenomena occurs once we move to nu-

tot clei with N > 126 with a clearer manifestation once we λ

/ 40 ff reach N = 184. For these nuclei, neutrons and protons λ 20 occupy orbits that differ in two units of the main oscil- lator quantum number, =2n + l. Under the assump- 0 N 66 68 70 72 74 tion of isospin symmetry the single particle orbits for proton number neutrons and protons will be identical and consequently Gamow-Teller transitions will be exactly zero. However, FIG. 11. (color online) The contribution of the first-forbidden due to isospin breaking mainly due to the coulomb inter- transitions in the total decay rate for the N = 50, N = 82 and action the proton and neutron single particle states are N = 126 isotonic chains, in the top, central and bottom pan- not identical allowing for Gamow-Teller transitions, even els, respectively. Present results, denoted with black squares, if somewhat suppressed. This is the reason why forbid- are compared with the FRDM calculation [19] (denoted with den transitions dominate in this whole region with GT green triangles), DF3 calculation [20] (denoted with blue di- transitions representing around 20% of the decay rate. amonds) and a recent shell model calculation [18] (denoted with red circles). C. β-delayed neutron emission the valley of stability with Z 50. Here, protons and ≤ β-delayed neutron emission is another important com- neutrons fill orbits with the same parity and the decay is ponent in the late stages of r-process nucleosynthesis. dominated by Gamow-Teller transitions. It is interesting Here we approximate the probability of emission of x neu- to notice the similar behaviour for nuclei located around trons as the ratio of the rates between S and S 78 132 xn (x+1)n Ni and Sn. Moving in the “south-west” direction separation energies to the total decay rate, i.e. (Z < 28, N < 50 and Z < 50, N < 82 respectively) forbidden transitions are initially suppressed. However, min(Qβ ,S(x+1)n) once protons start to occupy orbits with different parity λi

(Z < 20 or Z . 40) the contribution of forbidden tran- i,Ei=Sxn sitions grows to values of 40% to 60% depending on the Pxn = X . (21) λi detailed nuclear structure. Moving in the “east” direction i (N > 50 and N > 82 respectively), we see that the con- Neutron emission probabilitiesP are, therefore, very sensi- tribution of forbidden transitions is rather independent tive to the neutron separation energies and display strong of the charge number and does not vary when crossing odd-even effects. As mean-field models cannot repro- the proton shell closure at Z = 28 and Z = 50, respec- duce this staggering we employed the global nuclear mass tively. In these regions, for neutrons and protons there model [60] to obtain the neutron separation energies. Ad- are valence orbits with both positive and negative par- ditionally, only for the purpose of calculating the emis- 12 sion probabilities, we smeared the beta strength distribu- tions with a Lorentzian in order to reduce the sensitiv- ity of the neutron emission probabilities on the positions of particular transitions. The width of the Lorentzian was determined by calculating the average number of de- layed neutrons emitted by the fission fragments of 325U, n =0.0158 0.0005 [65]. In this way the width was ad- hjustedi to Γ =± 65 keV. From the probabilities, we obtain the average number of emitted neutrons after the decay of a nucleus as

n = iP . (22) h i in i X

In Fig. 13 we plot the average number of emitted neu- trons for neutron-rich nuclei included in this study. Close to the valley of stability the Q-values are smaller than the one neutron separation energies, making it impos- sible for the daughter nucleus to de-excite via neutron FIG. 13. (color online) Average number of emitted neutrons emission. With additional neutrons both the Q-value in- ∗ after β-decay calculated with the FRDM and the D3C mod- creases and the separation energies decrease enabling the els. Note that because the available FRDM results only con- emission of one, two or more neutrons after decay. This tain the emission probabilities for up to 3 neutrons emitted, is particularly evident above neutron shell closures where and in the present work the probabilities are obtained for up the average number of emitted neutrons increases due to to 5 emitted neutrons, the colour scales are different. a drop in the separation energies. Another characteristic feature of n is the odd-even staggering as a function h i of atomic number. Isotopic chains with an even num- D. Consequences for heavy-element ber of protons typically have a smaller average number nucleosynthesis of emitted neutrons, i.e. smaller probabilities of emit- ting more neutrons, than their odd-proton neighbouring chains which mimics the behaviour of the Q-values. A To explore the impact of the obtained results on heavy smaller odd-even staggering is also present within any element nucleosynthesis, we have performed r-process nu- isotopic chain as a consequence of the staggering of the cleosynthesis calculations based on the “hot” and “cold” neutron separation energies. r-process conditions from ref. [3]. In the “hot” r-process scenario the r-process evolves under an (n,γ) ⇄ (γ,n) While both the D3C∗ and the FRDM models repro- equilibrium that breaks down once the r-process freeze- duce these general features, in other aspects there are out is reached. For the “cold” r-process scenario the significant differences. Firstly, the FRDM evaluation evolutions proceeds by a competition between neutron- was limited to a maximum of 3 emitted neutrons per captures and beta-decays. The neutron capture and pho- decay, while the present study takes into account the todissociation rates used on the network calculations are possibility for up to 5 neutrons to be emitted in a de- based on the statistical model approach [66] using the cay. The consequence of this change is evident in very FRDM mass model [37]. The only difference between the neutron-rich nuclei where the neutron separation ener- calculations is the use of the FRDM+QRPA beta-decay gies are very small, and where the D3C* model predicts half-lives and beta-delayed neutron emission probabili- a larger average number of emitted neutrons which con- ties [19] or those computed in the present study. tinuously increases up to the drip-line, while the FRDM In Fig. 14 we compare the abundances obtained with model predicts a saturation of n . On the other hand, the new decay half-lives and delayed neutron emission for isotopes closer to stability theh i D3C∗ model predicts probabilities with the results obtained using the FRDM smaller vales of n and a more gradual increase towards values for both the hot and cold r-process trajectories. the neutron drip-line.h i Finally, in very neutron-rich nuclei One can notice important differences particularly in the the FRDM displays decidedly unphysical, sudden oscil- region of the third r-process peak, A 195 between the lations between the values of n = 1 and n = 2 while calculations. There are also differences∼ in the abundances the present study predicts ah smoothi increaseh i with the around the second peak but the astrophysical conditions neutron number. explored by the present trajectories are expected to con- 13

-3 2 10 FRDM 1 D3C* 10 FRDM, hot SN solar 0 D3C*, hot SN -4 10 -1 10 -2 -5 10 seed log Y(A) -3 10 /Y n -4 -6 Y 10 -5 10 -6 -7 10 120 140 160 180 200 220 -7 A 10 -3 -8 FRDM 10 D3C* 1 solar 10 FRDM, cold SN -4 0 D3C*, cold SN 10 -1 10 -5 -2 10 log Y(A)

seed -3 10 /Y -6 n -4 Y 10 -5 10

-7 -6 120 140 160 180 200 220 10 A -7 10 -8 FIG. 14. (color online) The abundances of heavy nuclei re- 10 sulting from a hot r-process (upper panel) or a cold r-process 0.1 1 10 (lower panel) calculation using the half-lives obtained from t [s] the FRDM and the current model. The solid circles denote the solar r-process abundances. FIG. 15. (color online) The evolution of the neutron-to-seed ratio during the r-process. The upper panel shows the results for the hot r-process while the lower panel for the cold r- process. Dashed lines denote the results obtained with the tribute mainly to the third r-process peak region, see FRDM decay data, and full lines denote results obtained in ref. [3]. To get further insight in the origin of the differ- the present study. ences in the abundances we show in figure 15 the evolu- tion of the neutron-to-seed ratio in the present calcula- tions. We clearly see that for those calculations that use beta-flow equilibrium by the freeze-out time. Under these the D3C* the r-process freeze-out, defined as the moment conditions the product of the average beta-decay rate when the neutron-to-seed ratio reaches a value of 1, oc- times the abundance for each isotopic chain is constant, curs at an earlier time. An earlier freeze-out has impor- i.e. λβ(Z)Y (Z) = constant, with Y (Z) = A Y (Z, A) tant consequences for the shaping of the final abundance and pattern during the decay to stability. This is because it P occurs at a higher value of the density that speeds up the freeze-out producing a faster and larger drop in the 1 1 λβ (Z)= = λβ(Z, A)Y (Z, A). (23) neutron-to-seed ratio. As a consequence, there is less τβ(Z) Y (Z) A time for neutron-captures during the freeze-out. These X neutron captures are in fact responsible for shifting the This means that the isotopic abundances, Y (Z), are pro- peaks to higher A values as discussed in ref. [3]. portional to the lifetime, τ(Z), and the material accumu- The evolution proceeds faster with the new half-lives lates in regions with the longest lifetimes that correspond because they tend to be shorter than those computed to the magic shell closures. Figure 16, shows the isotopic by the FRDM+QRPA approach. This is particularly lifetimes at the freeze-out for the different calculations. the case near to the magic shell closures N = 82 and One can clearly see that the lifetimes are substantially N = 126 where the r-process moves closer to the stabil- longer for the FRDM+QRPA approach particularly close ity and reaches the longest living nuclei. Both for the hot to the magic neutron numbers N = 82, corresponding to and cold r-process conditions, the r-process has reached Z 50, and N = 126, corresponding to Z 70. ∼ ∼ 14

15 FRDM, hot SN D3C*, hot SN 0.1

(Z) [s] 10 β

τ 0.01 [MeV] 0.001 e

_ 5 E

FRDM, cold SN 0.1 D3C*, cold SN 0 (Z) [s]

β 0 5 10 15 20 25 30 35

τ 0.01 15

0.001 40 50 60 70 80 90 proton number 10 FIG. 16. (color online) The isotopic lifetimes defined in Eq. [MeV] _ (23) as a function of the atomic number, at the moment of ν

_ 5 r-process freze-out. In the top panel we show the results for E the hot r-process, and in the bottom panel the results for the cold r-process. 0 0 5 10 15 20 25 30 35 15

We expect that the variations discussed above will re- 10 main for a broad range of r-process conditions and in particular also for neutron-star merger conditions. This [MeV] aspect will be explored in more detail in a future pub- γ _

E 5 lication. An important aspect for r-process calculations under r-process conditions is the energy generation dur- ing the r-process. This depends mainly on the beta-decay 0 rates and determines the evolution of temperature dur- 0 5 10 15 20 25 30 35 ing the r-process [67]. The nuclear energy production Q [MeV] is also important to determine the evolution of material ejected in highly eccentric loosely bound orbits. Depend- FIG. 17. (color online) Average energy of the electrons, neu- ing on the energy produced by the r-process this mate- trinos and photons produced after beta-decay rial can become unbound affecting the late time evolu- tion of short gamma-ray burst afterglows assuming they values. are powered by matter fallback [68]. At later times the radioactive heating from the beta-decay of r-process ma- 1 1 E¯1 (Q)= Q, E¯2 (Q)= Q. (24) terial will also affect the neutron-star merger remnant e,ν,γ¯ 2 e,ν,γ¯ 3 evolution [69] and the electromagnetic transient, kilo- nova, lightcurve [70–72]. For all these applications it is One can see that the ratio of the average electron en- important to know how much of the energy produced ergy and the nuclear Q-value is concentrated around one by beta-decay is lost in the form of neutrinos and how third, with the same artio for the antineutrinos is in the the remaining energy is distributed between the electrons range 0.3-0.5. r-process nucleosynthesis is not expected and gamma-rays emitted from nuclear deexcitation. We to be sensitive to this range of variations in the heating have determined for each decay these quantities based on rate [67]. However, kilonova light curves are more sensi- the computed beta-decay strengths. A small part of the tive to this range of variation [70] and our calculations are energy will be carried away by the neutrons emitted af- the first to provide this information for future modeling. ter the beta-decay. Our gamma channel does in fact also include this energy as it is assumed that the emitted neu- trons will thermalize in a relatively short timescale. All IV. CONCLUSION AND OUTLOOK this information is included in the supplemental material that also includes the relevant beta-decay rates. In this study we have performed a large scale calcu- lation of β-decay rates and β-delayed neutron emission Figure 17 shows the average energy of electrons, neu- probabilities for neutron-rich nuclei with 8 Z 110. A trinos and photons produced after beta-decay as a func- fully self-consistent theoretical framework≤ was employed≤ tion of the decay Q-value. The red lines denote the one- based on the relativistic nuclear energy density func- half and one-third ratios between the energy and the Q- tional. The advantages of this approach are: (i) self- 15 consistent modeling of all relevant transition matrix el- ments is a consequence of a cancelation of a significant ements without any adjustment of the model parame- overestimate of the half-lives of even-Z (especially even- ters to the nucleus under consideration, and (ii) equal even) nuclei with a smaller underestimate of the half-lives treatment of allowed and first-forbidden transitions. The of odd-Z nuclei. In both cases, the description improves r-process involves even-even nuclei, but also odd-A and for shorter lived nuclei, due to the importance of the odd-nuclei. It is, therefore, very important that all types lepton phase space that increases sharply with larger Q- of nuclei are treated within a single theoretical frame- values. work. We have tested the quality of the description of The consistent inclusion of the first-forbidden transi- odd nuclei by comparing the β-decay Q-values with those tions in this study has allowed for a systematic study of of the FRDM and the data for the silver and cadmium the contribution of parity changing transitions to the de- isotopic chains, where the results show very good agree- cay rate for a large number of nuclei. In specific cases, ment with the available data. Additionally, in agreement namely the N = 50, N = 82 and N = 126 isotonic chains, with data there is no noticeable odd-even staggering in the results of the present study agree well with date and the obtained half-lives, and the statistical analysis shows results based on other models, except in particular cases only a small tendency to overestimate the half-lives of such as the nuclei above the Z = 28 proton shell where odd-Z and odd-odd nuclei. Additionally, even though we predict a strong contribution of the forbidden tran- the model used assumes a spherical shape of the nucleus sitions. For N = 126 isotones, the model predicts an under consideration, it has been shown that it can ac- increase of the contribution of first forbidden transitions curately describe the reaction rates both in the regions to the total decay rate, in agreement with the results of the nuclear chart close to closed shells and in open- obtained by other models. For nuclei with N < 126, we shell nuclei known to be well-deformed. Thus, the model predict a smooth increase of the contribution of the parity employed is shown to reproduce β-decay data across the changing transitions from the valley of stability towards whole neutron-rich side of the nuclear chart with good the neutron drip-line, and a saturation at approximately precision, regardless of the type of nucleus. 40%-60% of the total rate. Nuclei with N > 126 show In general, the key features of the decay rates across a suppression of Gamow-Teller transitions that became the nuclear chart have been reproduced. We obtain a only possible by isospin breaking terms mainly due to very good agreement with the data for isotopic chains coulomb. with both an even and odd atomic number, without Closely related to β-decay is the process of β-delayed any noticeable odd-even staggering along the chain and neutron emission, which also affects the late stages of the closely following the trend of the measurements. This r-process by contributing to the neutron flux and the re- characteristic ensures a smoother behaviour during the fi- distribution of mass. The characteristic quantity is the nal decay towards the line of stability in the late stages of average number of emitted neutrons which is a measure the heavy element nucleosynthesis. The FRDM results, of the decay rate that falls above the neutron separation however, often displays erratic jumps in the half-lives energy. This quantity depends on both the Q-value and without an underlying physical reason, as is evident in the one, two and more neutron separation energies and e.g. strontium and cadmium isotopic chains. In compar- displays odd-even staggering both between the isotopic ison to recent data, both models produce similar discrep- chains and along a particular chain. We have extracted ancies. In the light nuclei region, both models reproduce the probabilities for emission of one or more neutrons and the data within a factor of 2, for a wide selection of nuclei the average number of emitted neutrons from the previ- with atomic numbers ranging from Z = 36 to Z = 43. In ously calculated decay data. The results show a physical the region of very heavy nuclei, the FRDM overestimates saturation of the average number of emitted neutrons the half-lives of nuclei below the Z = 82 closed proton closer to the neutron drip-line which indicates that the shell, while the present work provides a reasonable de- decaying matter could provide a significant number of scription of the decay properties. However, around the neutrons. In contrast, the FRDM predicts a similar but N = 82 shell, the D3C∗ model significantly underesti- overall lower number of emitted neutrons per decay, with mates the half-lives. Extending the analysis to the whole oscillating behaviour close to the neutron drip line. neutron-rich part of the nuclear chart, we obtain very The present results have been applied in an r-process good agreement with the data, except in the regions just calculation in a hot and cold r-process scenarios corre- above a closed proton shell, but below the next closed sponding to matter ejected under moderately high en- neutron shell, i.e. the same regions where the model pre- tropy. We find substantial changes in the r-process abun- dicts high contribution of the first-forbidden transitions. dances that are mainly due to the fact that for r-process We also performed a statistical analysis, where the stud- nuclei the present half-lives are systematically shorter ied the quality of the half-life predictions in contrast to than those computed by the FRDM+QRPA approach. the measurements, and depending on the type of nuclei. This study has shown that the theoretical framework We obtain a very good average reproduction of the data, based on the relativistic nuclear energy density functional with a tendency to overestimate the half-lives of odd-Z is a mature one, capable of providing a good description nuclei. In contrast with this result, the very good overall of sensitive physical quantities. However, the description agreement of the FRDM predictions with the measure- of half-lives of a large range of nuclei remains challenging, 16 especially for nuclei close to the valley of stability. Fur- the detailed structure of the transition strength. ther advances in this field will require effort along three possible directions: (i) improving the relativistic energy density functional will serve to provide more accurate Qβ ACKNOWLEDGMENTS values and enhance the description of the energy gener- ation during the decay released in the form of electron, This work was supported in part by the Helmholtz antineutrinos and gamma-rays; (ii) a consistent treat- International Center for FAIR within the framework ment of deformations will be necessary to enhance the of the LOEWE program launched by the State of description of the transitions in particular open-shell nu- Hesse, by the BMBF-Verbundforschungsprojekt num- clei; and (iii) inclusion of complex-configurations in the ber 06DA7047I, the Helmholtz Association through the QRPA basis will enrich the transition spectrum, espe- Nuclear Astrophysics Virtual Institute (VH-VI-417), by cially so in the low-energy region, and provide a better the Ministry of Science, Education and Sport (MZOS) description of beta-decays, but also give information on Project No. 1191005-1010, by the IAEA Research Con- tract No. 18094/R0, and by the DAAD-MZOS bilateral collaboration program.

[1] K. Langanke and G. Martinez-Pinedo, [9] S. Nishimura, Z. Li, H. Watanabe, K. Yoshinaga, Rev. Mod. Phys. 75, 819 (2003). T. Sumikama, T. Tachibana, K. Yamaguchi, M. Kurata- [2] M. Arnould, S. Goriely, and K. Takahashi, Nishimura, G. Lorusso, Y. Miyashita, A. Odahara, Phys. Repts. 450, 97 (2007). H. Baba, J. S. Berryman, N. Blasi, A. Bracco, F. Cam- [3] A. Arcones and G. Mart´ınez-Pinedo, era, J. Chiba, P. Doornenbal, S. Go, T. Hashimoto, Phys. Rev. C 83, 045809 (2011). S. Hayakawa, C. Hinke, E. Ideguchi, T. Isobe, [4] P. T. Hosmer, H. Schatz, A. Aprahamian, O. Arndt, Y. Ito, D. G. Jenkins, Y. Kawada, N. Kobayashi, R. R. C. Clement, A. Estrade, K.-L. Kratz, S. N. Liddick, Y. Kondo, R. Kr¨ucken, S. Kubono, T. Nakano, H. J. P. F. Mantica, W. F. Mueller, F. Montes, A. C. Mor- Ong, S. Ota, Z. Podoly´ak, H. Sakurai, H. Scheit, ton, M. Ouellette, E. Pellegrini, B. Pfeiffer, P. Reeder, K. Steiger, D. Steppenbeck, K. Sugimoto, S. Takano, P. Santi, M. Steiner, A. Stolz, B. E. Tomlin, W. B. Wal- A. Takashima, K. Tajiri, T. Teranishi, Y. Wakabayashi, ters, and A. Wohr, Physical Review Letters 94, 112501 P. M. Walker, O. Wieland, and H. Yamaguchi, (2005). Physical Review Letters 106, 052502 (2011). [5] P. Hosmer, H. Schatz, A. Aprahamian, O. Arndt, [10] I. Dillmann, K.-L. Kratz, A. W¨ohr, O. Arndt, R. R. C. Clement, A. Estrade, K. Farouqi, K.-L. Kratz, B. A. Brown, P. Hoff, M. Hjorth-Jensen, S. N. Liddick, A. F. Lisetskiy, P. F. Mantica, P. M¨oller, U. K¨oster, A. N. Ostrowski, B. Pfeiffer, D. Sew- W. F. Mueller, F. Montes, A. C. Morton, M. Ouellette, eryniak, J. Shergur, and W. B. Walters, E. Pellegrini, J. Pereira, B. Pfeiffer, P. Reeder, P. Santi, Physical Review Letters 91, 162503 (2003). M. Steiner, A. Stolz, B. E. Tomlin, W. B. Walters, and [11] B. Pfeiffer, K. L. Kratz, F. K. Thielemann, and W. B. A. W¨ohr, Physical Review C 82, 025806 (2010). Walters, Nuclear Physics A 693, 282 (2001). [6] M. Quinn, A. Aprahamian, J. Pereira, R. Sur- [12] J. Benlliure, H. Alvarez-Pol, T. Kurtukian-Nieto, man, O. Arndt, T. Baumann, A. Becerril, T. El- A. I. Morales, L. Audouin, F. Becker, B. Blank, liot, A. Estrade, D. Galaviz, T. Ginter, M. Haus- I. Borzov, E. Casarejos, D. Cortina-Gil, B. Fer- mann, S. Hennrich, R. Kessler, K.-L. Kratz, G. Lorusso, nandez, J. Giovinazzo, D. Henzlova, B. Jurado, P. F. Mantica, M. Matos, F. Montes, B. Pfeiffer, L. Langanke, G. Martinez-Pinedo, P. Napoli- M. Portillo, H. Schatz, F. Schertz, L. Schnorrenberger, tani, J. Pereira, F. Rejmund, K. H. Schmidt, E. Smith, A. Stolz, W. B. Walters, and A. W¨ohr, O. Yordanov, and the RISING Collaboration, Physical Review C 85, 035807 (2012). Journal of Physics: Conference Series 337, 012070 (2012). [7] M. Madurga, R. Surman, I. N. Borzov, R. Grzywacz, [13] G. Benzoni, A. Morales, J. Valiente-Dobn, A. Gottardo, K. P. Rykaczewski, C. J. Gross, D. Miller, D. W. A. Bracco, F. Camera, F. Crespi, A. Corsi, S. Leoni, Stracener, J. C. Batchelder, N. T. Brewer, L. Cartegni, B. Million, R. Nicolini, O. Wieland, A. Gadea, S. Lu- J. H. Hamilton, J. K. Hwang, S. H. Liu, S. V. nardi, P. Boutachkov, A. Bruce, M. Grska, J. Grebosz, Ilyushkin, C. Jost, M. Karny, A. Korgul, W. Kr´olas, S. Pietri, Z. Podolyak, M. Pftzner, P. Regan, H. We- A. Ku´zniak, C. Mazzocchi, A. J. Mendez, K. Miernik, ick, J. A. Nez, A. Algora, N. Al-Dahan, G. de An- S. W. Padgett, S. V. Paulauskas, A. V. Ramayya, J. A. gelis, Y. Ayyad, N. Alkhomashi, P. Allegro, D. Baz- Winger, M. Woli ´nska Cichocka, and E. F. Zganjar, zacco, J. Benlliure, M. Bowry, M. Bunce, E. Casare- Physical Review Letters 109, 112501 (2012). jos, M. Cortes, A. D. Bacelar, A. Deo, C. Domingo- [8] F. Montes, A. Estrade, P. T. Hosmer, S. N. Liddick, P. F. Pardo, M. Doncel, Z. Dombradi, T. Engert, K. Ep- Mantica, A. C. Morton, W. F. Mueller, M. Ouellette, pinger, G. Farrelly, F. Farinon, E. Farnea, H. Geissel, E. Pellegrini, P. Santi, H. Schatz, A. Stolz, B. E. Tom- J. Gerl, N. Goel, E. Gregor, T. Habermann, R. Hoischen, lin, O. Arndt, K.-L. Kratz, B. Pfeiffer, P. Reeder, W. B. R. Janik, S. Klupp, I. Kojouharov, N. Kurz, S. Man- Walters, A. Aprahamian, and A. Wohr, Physical Review dal, R. Menegazzo, D. Mengoni, D. Napoli, F. Naqvi, C 73, 035801 (2006). C. Nociforo, A. Prochazka, W. Prokopowicz, F. Recchia, 17

R. Ribas, M. Reed, D. Rudolph, E. Sahin, H. Schaffner, [29] N. Paar, T. Nikˇsi´c, D. Vretenar, and P. Ring, Physical A. Sharma, B. Sitar, D. Siwal, K. Steiger, P. Strmen, Review C 69, 054303 (2004). T. Swan, I. Szarka, C. Ur, P. Walker, and H.-J. Woller- [30] T. Marketin, G. Mart´ınez-Pinedo, N. Paar, and D. Vrete- sheim, Physics Letters B 715, 293 (2012). nar, Physical Review C 85, 054313 (2012). [14] C. Domingo-Pardo, R. Caballero-Folch, J. Agramunt, [31] N. Paar, D. Vretenar, T. Marketin, and P. Ring, Physical A. Algora, A. Arcones, F. Ameil, Y. Ayyad, J. Benl- Review C 77, 024608 (2008). liure, M. Bowry, F. Calvi˜no, D. Cano-Ott, G. Cort´es, [32] Y. F. Niu, N. Paar, D. Vretenar, and J. Meng, T. Davinson, I. Dillmann, A. Estrade, A. Evdokimov, Physical Review C 83, 045807 (2011). T. Faestermann, F. Farinon, D. Galaviz, A. Garc´ıa-Rios, [33] J. Berger, M. Girod, and D. Gogny, Nuclear Physics A H. Geissel, W. Gelletly, R. Gernh¨auser, M. B. G´omez- 428, 23 (1984). Hornillos, C. Guerrero, M. Heil, C. Hinke, R. Kn¨obel, [34] N. Paar, P. Ring, T. Nikˇsi´c, and D. Vretenar, Physical I. Kojouharov, J. Kurcewicz, N. Kurz, Y. Litvinov, Review C 67, 034312 (2003). L. Maier, J. Marganiec, M. Marta, T. Mart´ınez, [35] P. Ring and P. Schuck, The Nuclear Many-Body Problem G. Mart´ınez-Pinedo, B. S. Meyer, F. Montes, I. Mukha, (Springer, Berlin, Heidelberg, New York, 1980). D. R. Napoli, C. Nociforo, C. Paradela, S. Pietri, [36] G. Audi, F. G. Kondev, M. Wang, B. Pfeiffer, X. Sun, Z. Podoly´ak, A. Prochazka, S. Rice, A. Riego, B. Ru- J. Blachot, and M. MacCormick, Chinese Physics C 36, bio, H. Schaffner, C. Scheidenberger, K. Smith, E. Sokol, 1157 (2012). K. Steiger, B. Sun, J. L. Ta´ın, M. Takechi, D. Testov, [37] P. M¨oller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, H. Weick, E. Wilson, J. S. Winfield, R. Wood, Atomic Data and Nuclear Data Tables 59, 185 (1995). P. Woods, and A. Yeremin, ArXiv e-prints (2013), [38] S. Goriely, N. Chamel, and J. Pearson, arXiv:1309.3047 [nucl-ex]. Physical Review C 82, 035804 (2010). [15] T. Kurtukian-Nieto, J. Benlliure, K.-H. Schmidt, L. Au- [39] T. Duguet, P. Bonche, P.-H. Heenen, and J. Meyer, douin, F. Becker, B. Blank, I. Borzov, E. Casare- Physical Review C 65, 014310 (2001). jos, F. Farget, M. Fern´andez-Ord´o˜nez, J. Giov- [40] T. Duguet, P. Bonche, P.-H. Heenen, and J. Meyer, inazzo, D. Henzlova, B. Jurado, K. Langanke, Physical Review C 65, 014311 (2001). G. Mart´ınez-Pinedo, J. Pereira, and O. Yordanov, [41] H. Behrens and J. J¨anecke, Landoldt-B¨ornstein Tables, Eur. Phys. J. A 50, 135 (2014). Gruppe I, Vol. 4 (Springer, Berlin, 1969). [16] A. I. Morales, J. Benlliure, T. Kurtuki´an-Nieto, K.- [42] J. C. Hardy and I. S. Towner, H. Schmidt, S. Verma, P. H. Regan, Z. Podoly´ak, Physical Review C 79, 055502 (2009). M. G´orska, S. Pietri, R. Kumar, E. Casarejos, N. Al- [43] A. R. Edmonds, Angular Momentum in Quantum Me- Dahan, A. Algora, N. Alkhomashi, H. Alvarez-Pol,´ chanics (Princeton University Press, Princeton, New Jer- G. Benzoni, A. Blazhev, P. Boutachkov, A. M. Bruce, sey, 1960). L. S. C´aceres, I. J. Cullen, A. M. Denis Bacelar, P. Door- [44] J. Beringer, J. F. Arguin, R. M. Barnett, K. Copic, nenbal, M. E. Est´evez-Aguado, G. Farrelly, Y. Fujita, O. Dahl, D. E. Groom, C. J. Lin, J. Lys, H. Mu- A. B. Garnsworthy, W. Gelletly, J. Gerl, J. Grebosz, rayama, C. G. Wohl, W. M. Yao, P. A. Zyla, C. Am- R. Hoischen, I. Kojouharov, N. Kurz, S. Lalkovski, sler, M. Antonelli, D. M. Asner, H. Baer, H. R. Band, Z. Liu, C. Mihai, F. Molina, D. M¨ucher, B. Rubio, T. Basaglia, C. W. Bauer, J. J. Beatty, V. I. Be- H. Shaffner, S. J. Steer, A. Tamii, S. Tashenov, J. J. lousov, E. Bergren, G. Bernardi, W. Bertl, S. Bethke, Valiente-Dob´on, P. M. Walker, H. J. Wollersheim, and H. Bichsel, O. Biebel, E. Blucher, S. Blusk, G. Brooij- P. J. Woods, Phys. Rev. Lett. 113, 022702 (2014). mans, O. Buchmueller, R. N. Cahn, M. Carena, A. Cec- [17] E. Caurier, G. Mart´ınez-Pinedo, cucci, D. Chakraborty, M. C. Chen, R. S. Chivukula, F. Nowacki, A. Poves, and A. P. Zuker, G. Cowan, G. D’Ambrosio, T. Damour, D. de Florian, Reviews of Modern Physics 77, 427 (2005). A. de Gouvˆea, T. DeGrand, P. de Jong, G. Dissertori, [18] Q. Zhi, E. Caurier, J. J. Cuenca-Garc´ıa, K. Lan- B. Dobrescu, M. Doser, M. Drees, D. A. Edwards, S. Ei- ganke, G. Mart´ınez-Pinedo, and K. Sieja, delman, J. Erler, V. V. Ezhela, W. Fetscher, B. D. Physical Review C 87, 025803 (2013). Fields, B. Foster, T. K. Gaisser, L. Garren, H. J. Gerber, [19] P. M¨oller, B. Pfeiffer, and K.-L. Kratz, Physical Review G. Gerbier, T. Gherghetta, S. Golwala, M. Goodman, C 67, 055802 (2003). C. Grab, A. V. Gritsan, J. F. Grivaz, M. Gr¨unewald, [20] I. N. Borzov, Nuclear Physics A 777, 645 (2006). A. Gurtu, T. Gutsche, H. E. Haber, K. Hagiwara, [21] D.-L. Fang, B. A. Brown, and T. Suzuki, C. Hagmann, C. Hanhart, S. Hashimoto, K. G. Hayes, Physical Review C 88, 034304 (2013). M. Heffner, B. Heltsley, J. J. Hern´andez-Rey, K. Hikasa, [22] J. Engel, M. Bender, J. Dobaczewski, W. Nazarewicz, A. H¨ocker, J. Holder, A. Holtkamp, J. Huston, J. D. and R. Surman, Physical Review C 60, 014302 (1999). Jackson, K. F. Johnson, T. Junk, D. Karlen, D. Kirkby, [23] T. Nikˇsi´c, T. Marketin, D. Vretenar, N. Paar, and S. R. Klein, E. Klempt, R. V. Kowalewski, F. Krauss, P. Ring, Physical Review C 71, 014308 (2005). M. Kreps, B. Krusche, Y. V. Kuyanov, Y. Kwon, O. La- [24] T. Marketin, D. Vretenar, and P. Ring, Physical Review hav, J. Laiho, P. Langacker, A. Liddle, Z. Ligeti, T. M. C 75, 024304 (2007). Liss, L. Littenberg, K. S. Lugovsky, S. B. Lugovsky, [25] Z. Niu, Y. Niu, H. Liang, W. Long, T. Nikˇsi´c, D. Vrete- T. Mannel, A. V. Manohar, W. J. Marciano, A. D. Mar- nar, and J. Meng, Physics Letters B 723, 172 (2013). tin, A. Masoni, J. Matthews, D. Milstead, R. Miquel, [26] N. Paar, D. Vretenar, E. Khan, and G. Colo, Reports K. M¨onig, F. Moortgat, K. Nakamura, M. Narain, P. Na- on Progress in Physics 70, 691 (2007). son, S. Navas, M. Neubert, P. Nevski, Y. Nir, K. A. [27] I. N. Borzov, Physical Review C 67, 025802 (2003). Olive, L. Pape, J. Parsons, C. Patrignani, J. A. Pea- [28] D. Vretenar, A. Afanasjev, G. Lalazissis, and P. Ring, cock, S. T. Petcov, A. Piepke, A. Pomarol, G. Punzi, Physics Reports 409, 101 (2005). A. Quadt, S. Raby, G. Raffelt, B. N. Ratcliff, P. Richard- 18

son, S. Roesler, S. Rolli, A. Romaniouk, L. J. Rosen- naga, Physical Review Letters 114, 192501 (2015). berg, J. L. Rosner, C. T. Sachrajda, Y. Sakai, G. P. [56] T. Marketin, E. Litvinova, D. Vretenar, and P. Ring, Salam, S. Sarkar, F. Sauli, O. Schneider, K. Scholberg, Physics Letters B 706, 477 (2012). D. Scott, W. G. Seligman, M. H. Shaevitz, S. R. Sharpe, [57] E. Litvinova, B. Brown, D.-L. Fang, T. Marketin, and M. Silari, T. Sj¨ostrand, P. Skands, J. G. Smith, G. F. R. Zegers, Physics Letters B 730, 307 (2014). Smoot, S. Spanier, H. Spieler, A. Stahl, T. Stanev, [58] E. Litvinova, Physical Review C 85, 021303 (2012). S. L. Stone, T. Sumiyoshi, M. J. Syphers, F. Takahashi, [59] Y. F. Niu, Z. Niu, G. Col`o, and E. Vigezzi, M. Tanabashi, J. Terning, M. Titov, N. P. Tkachenko, Physical Review Letters 114, 142501 (2015). N. A. T¨ornqvist, D. Tovey, G. Valencia, K. van Bib- [60] M. Liu, N. Wang, Y. Deng, and X. Wu, ber, G. Venanzoni, M. G. Vincter, P. Vogel, A. Vogt, Phys. Rev. C 84, 014333 (2011). W. Walkowiak, C. W. Walter, D. R. Ward, T. Watari, [61] C. Bauer, T. Behrens, V. Bildstein, A. Blazhev, G. Weiglein, E. J. Weinberg, L. R. Wiencke, L. Wolfen- B. Bruyneel, J. Butterworth, E. Cl´ement, L. Coquard, stein, J. Womersley, C. L. Woody, R. L. Workman, A. Ya- J. L. Egido, A. Ekstr¨om, C. R. Fitzpatrick, C. Fransen, mamoto, G. P. Zeller, O. V. Zenin, J. Zhang, R. Y. Zhu, R. Gernh¨auser, D. Habs, H. Hess, J. Leske, T. Kr¨oll, G. Harper, V. S. Lugovsky, and P. Schaffner (Particle R. Kr¨ucken, R. Lutter, P. Marley, T. M¨oller, T. Otsuka, Data Group), Phys. Rev. D 86, 010001 (2012). N. Patronis, A. Petts, N. Pietralla, T. R. Rodr´ıguez, [45] H. Behrens and W. B¨uhring, N. Shimizu, C. Stahl, I. Stefanescu, T. Stora, P. G. Thi- Nuclear Physics A 162, 111 (1971). rolf, D. Voulot, J. van de Walle, N. Warr, F. Wenander, [46] H. Behrens and W. B¨uhring, Electron radial wave func- and A. Wiens, Phys. Rev. C 86, 034310 (2012). tions and nuclear β-decay (Clarendon, Oxford, 1982). [62] T. R. Rodr´ıguez, Phys. Rev. C 90, 034306 (2014). [47] B. H. Wildenthal, M. S. Curtin, and B. A. Brown, [63] P. M¨oller and J. Randrup, Nuclear Physics A 514, 1 Phys. Rev. C 28, 1343 (1983). (1990). [48] G. Mart´ınez-Pinedo, A. Poves, E. Caurier, and A. P. [64] M. Eichler, A. Arcones, A. Kelic, O. Korobkin, K. Lan- Zuker, Phys. Rev. C 53, R2602 (1996). ganke, T. Marketin, G. Martinez-Pinedo, I. V. Panov, [49] E. K. Warburton, Phys. Rev. C 42, 2479 (1990). T. Rauscher, S. Rosswog, C. Winteler, N. T. Zin- [50] E. K. Warburton, Physical Review C 44, 233 (1991). ner, and F.-K. Thielemann, ArXiv e-prints (2014), [51] E. K. Warburton, J. A. Becker, B. A. Brown, and D. J. arXiv:1411.0974 [astro-ph.HE]. Millener, Annals of Physics 187, 471 (1988). [65] G. R. Keepin, T. F. Wimett, and R. K. Zeigler, [52] E. K. Warburton and I. S. Towner, Physical Review 107, 1044 (1957). Phys. Repts. 242, 103 (1994). [66] T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data [53] T. Suzuki, T. Yoshida, T. Kajino, and T. Otsuka, Tables 75, 1 (2000). Phys. Rev. C 85, 015802 (2012). [67] J. J. Mendoza-Temis, G. Mart´ınez-Pinedo, K. Langanke, [54] P. M¨oller, J. R. Nix, and K.-L. Kratz, Atomic Data and A. Bauswein, and H.-T. Janka, ArXiv e-prints (2014), Nuclear Data Tables 66, 131 (1997). arXiv:1409.6135 [astro-ph.HE]. [55] G. Lorusso, S. Nishimura, Z. Y. Xu, A. Jung- [68] B. D. Metzger, A. Arcones, claus, Y. Shimizu, G. S. Simpson, P.-A. S¨oderstr¨om, E. Quataert, and G. Mart´ınez-Pinedo, H. Watanabe, F. Browne, P. Doornenbal, G. Gey, Mon. Not. Roy. Ast. Soc. 402, 2771 (2010). H. S. Jung, B. Meyer, T. Sumikama, J. Taprogge, [69] S. Rosswog, O. Korobkin, A. Arcones, F.-K. Thielemann, Z. Vajta, J. Wu, H. Baba, G. Benzoni, K. Y. Chae, and T. Piran, Mon. Not. Roy. Ast. Soc. 439, 744 (2014). F. C. L. Crespi, N. Fukuda, R. Gernh¨auser, N. In- [70] B. D. Metzger, G. Mart´ınez-Pinedo, S. Darbha, abe, T. Isobe, T. Kajino, D. Kameda, G. D. Kim, E. Quataert, A. Arcones, D. Kasen, R. Thomas, Y.-K. Kim, I. Kojouharov, F. G. Kondev, T. Kubo, P. Nugent, I. V. Panov, and N. T. Zinner, N. Kurz, Y. K. Kwon, G. J. Lane, Z. Li, A. Montaner- Mon. Not. Roy. Ast. Soc. 406, 2650 (2010). Piz´a, K. Moschner, F. Naqvi, M. Niikura, H. Nishibata, [71] L. F. Roberts, D. Kasen, W. H. Lee, and E. Ramirez- A. Odahara, R. Orlandi, Z. Patel, Z. Podoly´ak, H. Saku- Ruiz, Astrophys. J. 736, L21 (2011). rai, H. Schaffner, P. Schury, S. Shibagaki, K. Steiger, [72] A. Bauswein, S. Goriely, and H.-T. Janka, H. Suzuki, H. Takeda, A. Wendt, A. Yagi, and K. Yoshi- Astrophys. J. 773, 78 (2013).