First determination of \beta-delayed multiple neutron emission beyond A=100 through direct neutron measurement: The P{}_{2n} value of {}^{136}Sb.

First determination of -delayed multiple neutron emission beyond through direct neutron measurement: The P value of Sb.

R. Caballero-Folch TRIUMF, Vancouver, British Columbia V6T 2A3, Canada    I. Dillmann TRIUMF, Vancouver, British Columbia V6T 2A3, Canada GSI Helmholtz Center for Heavy Ion Research, D-64291 Darmstadt, Germany    J. Agramunt IFIC, CSIC - Universitat de València, E-46071 València, Spain    J.L. Taín IFIC, CSIC - Universitat de València, E-46071 València, Spain    A. Algora IFIC, CSIC - Universitat de València, E-46071 València, Spain Institute of Nuclear Research of the Hungarian Academy of Sciences, Debrecen H-4001, Hungary    J. Äysto University of Helsinki, Helsinki, Finland    F. Calviño Universitat Politècnica de Catalunya, E-08028 Barcelona, Spain    L. Canete University of Jyväskylä, FI-40014 Jyväskylä, Finland    G. Cortès Universitat Politècnica de Catalunya, E-08028 Barcelona, Spain    C. Domingo-Pardo IFIC, CSIC - Universitat de València, E-46071 València, Spain    T. Eronen University of Jyväskylä, FI-40014 Jyväskylä, Finland    E. Ganioglu University of Istanbul, 34134 Vezneciler, Turkey    W. Gelletly University of Surrey, Guildford GU2 7XH, United Kingdom    D. Gorelov University of Jyväskylä, FI-40014 Jyväskylä, Finland    V. Guadilla IFIC, CSIC - Universitat de València, E-46071 València, Spain    J. Hakala University of Jyväskylä, FI-40014 Jyväskylä, Finland    A. Jokinen University of Jyväskylä, FI-40014 Jyväskylä, Finland    A. Kankainen University of Jyväskylä, FI-40014 Jyväskylä, Finland    V. Kolhinen University of Jyväskylä, FI-40014 Jyväskylä, Finland    J. Koponen University of Jyväskylä, FI-40014 Jyväskylä, Finland    M. Marta GSI Helmholtz Center for Heavy Ion Research, D-64291 Darmstadt, Germany    E. Mendoza CIEMAT, E-28040 Madrid, Spain    A. Montaner-Pizá IFIC, CSIC - Universitat de València, E-46071 València, Spain    I. Moore University of Jyväskylä, FI-40014 Jyväskylä, Finland    C.R. Nobs University of Brighton, BN2 4AT Brighton, United Kingdom    S.E.A. Orrigo IFIC, CSIC - Universitat de València, E-46071 València, Spain    H. Penttilä University of Jyväskylä, FI-40014 Jyväskylä, Finland    I. Pohjalainen University of Jyväskylä, FI-40014 Jyväskylä, Finland    J. Reinikainen University of Jyväskylä, FI-40014 Jyväskylä, Finland    A. Riego Universitat Politècnica de Catalunya, E-08028 Barcelona, Spain    S. Rinta-Antila University of Jyväskylä, FI-40014 Jyväskylä, Finland    B. Rubio IFIC, CSIC - Universitat de València, E-46071 València, Spain    P. Salvador-Castiñeira Universitat Politècnica de Catalunya, E-08028 Barcelona, Spain    V. Simutkin University of Jyväskylä, FI-40014 Jyväskylä, Finland    A. Tarifeño-Saldivia IFIC, CSIC - Universitat de València, E-46071 València, Spain Universitat Politècnica de Catalunya, E-08028 Barcelona, Spain    A. Tolosa IFIC, CSIC - Universitat de València, E-46071 València, Spain    A. Voss University of Jyväskylä, FI-40014 Jyväskylä, Finland
July 16, 2019

Background: -delayed multiple neutron emission has been observed for some nuclei with A100, with Rb being the heaviest 2n emitter measured to date. So far, only 25 P values have been determined for the 300 nuclei that may decay in this way. Accordingly, it is of interest to measure P values for the other possible multiple neutron emitters throughout the chart of the nuclides. It is of particular interest to make such measurement for nuclei with A100 to test the predictions of theoretical models and simulation tools for the decays of heavy nuclei in the region of very neutron-rich nuclei. In addition, the decay properties of these nuclei are fundamental for the understanding of astrophysical nucleosynthesis processes such as the -process, and safety inputs for nuclear reactors.
Purpose: To determine for the first time the two neutron branching ratio, P value, for Sb through a direct neutron measurement, and to provide precise P values for Sb and Te.
Method: A pure beam of each isotope of interest was provided by the JYFLTRAP Penning trap at the IGISOL facility of the University of Jyväskylä, Finland. The purified ions were implanted into a moving tape at the end of the beam line. The detection setup consisted of a plastic scintillator placed right behind the implantation point after the tape to register the -decays, and the BELEN detector, based on neutron counters embedded in a polyethylene matrix. The analysis was based on the study of the - and neutron- growth-and-decay curves and the -one-neutron and -two-neutron time correlations, which allowed us the determination of the neutron-branching ratios.
Results: The P value of Sb was found to be 0.31(5)% and the measured P values for Sb and Te were found to be 31.9(15)% and 1.47(6)%, respectively. The measured P value is a factor 20 smaller than predicted by the FRDM+QRPA model used for -process calculations.

-delayed neutron emission, -decay, -process, nucleosynthesis, nuclear structure, neutron detector.
27.80.+w, 23.40.−s, 26.30.-k,21.10.-k

I Introduction

The Q value increases when going towards neutron-rich nuclei in the chart of nuclides. This makes the decay the dominant decay mode for this region. In very neutron-rich nuclei, the emission of one or more neutrons may also occur after a -decay. This process of -delayed neutron () emission was discovered in 1939 by Roberts et al. Roberts et al. (1939) and is energetically allowed when the Q value is positive, i.e. the Q value of the decay exceeds the neutron separation energy (S) of the daughter nucleus. This phenomenon becomes dominant when the populated state in the daughter nucleus, following the -decay, is higher in excitation energy than S, which gets lower going towards neutron-rich in the chart of nuclides. When the populated states lie even higher than the two-neutron separation energy, S, i.e. Q, two neutron emission may also occur. This phenomenon was predicted in 1960 by Goldansky Goldansky (1960). The first studies of multiple-neutron emission were carried out in the 1980s in which detection of two-neutron emission was observed for nuclei up to mass A=33 and predictions for the emission of three neutrons were made for masses around A=100 Lyutostansky et al. (1985). The -delayed two-neutron () emission probability, P value, has only been experimentally determined for 25 isotopes (two of them only approximations) Audi et al. (2017) of out of 300 potential multiple neutron emitter candidates Wang et al. (2017). The importance of providing more precise data for neutron emitting isotopes has been highlighted by the IAEA Coordinated Research Projects (CRP) on the Development of a Reference Database for Beta-Delayed Neutron Emission Dillmann et al. (2017) and on -delayed neutron emission evaluation Abriola et al. (2011). The latter emphasizes the importance of these data for safety improvement purposes for emerging nuclear power reactors, as well as for astrophysical studies. Indeed, these experimental data, involving isotopes in the neutron-rich region, are needed to achieve a better knowledge of the “freeze-out” of the rapid neutron capture (-process) Burbidge et al. (1957) when theoretical calculations of nucleosynthesis Arcones and Martínez-Pinedo (2011); Mumpower et al. (2016a); Surman et al. (2014) are applied. The delayed neutron emission is an important input for these models as it shifts the final isobaric solar system abundances of some species to lower masses, and introduces neutrons in the stellar environment that can be re-captured by other nuclei at later stages. This is of special interest in the regions of the -abundance peaks such as A130 near the doubly-magic Sn isotope (N=82 and Z=50). Up to now, the P data available for the heavy mass region is scarce, especially above  Caballero-Folch et al. (2016), and non-existent for multiple neutron emitters above .

This study presents a measurement of the P neutron branching ratio for Sb. With a Q window of 1884(6) keV Wang et al. (2017), Sb has been suggested to be a multiple-neutron emitter in several theoretical models Möller et al. (2003); Borzov (2016); Mumpower et al. (2016b) and by some experiments: these include estimates based on the Finite Fermi-system theory Lyutostansky and Panov (1983) and more recently an experiment involving several isotopes of mass in which the isobars could not be separated well enough to determine their neutron branching ratios precisely Testov et al. (2011). In order to get a successful measurement of multiple -delayed neutron emission, we used a high efficiency neutron detector and a beam free from contamination of other isotopes. The pure Sb beam was obtained with the Penning trap mass spectrometer, JYFLTRAP, at the IGISOL facility Äystö (2001); Hakala et al. (2012), located in the Accelerator Laboratory of the University of Jyväskylä, Finland. This manuscript describes the experimental setup in Section II, the details of the data analysis following the methodology reported in Agramunt et al. (2016) and the results in Section III. A summary and discussion is presented in Section IV, and the conclusions are given in Section V.

Ii Experimental setup

The experiment presented in this work was carried out at the Ion Guide Isotope Separator On-Line (IGISOL) facility of the University of Jyväskylä, Finland (JYFL). The isotopes of interest were produced in nuclear fission induced with a 25 MeV proton beam, with an intensity of around 7-10 , impinging on a natural uranium target. The fission fragments produced were extracted out from the helium gas cell using a sextupole ion guide (SPIG) Karvonen et al. (2008) and differential pumping, accelerated to 30 keV and mass-separated with a dipole magnet. The continuous beam was cooled and bunched using an RFQ cooler-buncher Nieminen et al. (2001) device before injecting the ions into the JYFLTRAP double Penning trap. The purification took place in the first trap, where a mass-selective buffer gas cooling technique Savard et al. (1991) was employed to resolve different isobars based on their cyclotron resonance frequencies , where q and m are the charge and the mass of the ion of interest and B is the magnetic field strength inside the trap. The description and the principles of the JYFLTRAP at IGISOL are reported in Ref. Eronen et al. (2012), and Figs. 1 and 2 of Ref. Kolhinen et al. (2013) detail the layout of the IGISOL facility.

In this experiment, pure ion beams of Rb, I, Te and Sb were extracted from the trap without isobaric contaminants and transported to the implantation system. Rb and I were used for calibration purposes, as their P values are well known from previous experiments Rudstam et al. (1993); Pfeiffer et al. (2002); Liang et al. (2018). The implantation system consisted of a moving magnetic tape placed inside an aluminum tube, with a thickness of 1 mm and diameter of 46 mm, under vacuum. This tube linked JYFLTRAP and the end of the beamline, shown in Fig. 1.

Figure 1: (Color online) Aluminum tube linking JYFLTRAP and the end of the beamline. The rear side of the moving tape at the implantation position can be seen. The plastic scintillator detector used as a decay counter, shown in the bottom-right, was placed 6 mm behind the moving implantation tape at the end of the aluminum tube.

The moving tape system allowed control of the ion-implantation (beam on) and -decay (beam off) times according to the half-life of the isotope measured, in order to be able to reproduce the growth-and-decay curves in the analysis (see Section III). The detection system in this experiment consisted of a 3-mm-thick plastic scintillator counter, shown in Fig. 1 (bottom-right), placed at the end of the vacuum tube, surrounded by the Beta dELayEd Neutron (BELEN) detector Torner et al. (2014); Gómez-Hornillos et al. (2011, 2014). The latter consisted of 48 He counter tubes of 2.54 cm diameter manufactured by LND Inc. LND (), distributed in three concentric rings, and embedded in a high-density polyethylene (HDPE) matrix to moderate the neutrons, see Fig. 2. The BELEN was surrounded by 20 cm of HDPE shielding in order to moderate and absorb neutrons scattered from the surroundings (see also Fig. 2).

Figure 2: (Color online) The BELEN detector embedded in the high-density polyethylene (HDPE) matrix and shielding, and the data acquisition system.

Table 1 summarizes its characteristics.

Central Ring
hole Inner Middle Outer
Position (diameter) 60 120 230 340
Number of counters (10 atm) 0 8 0
Number of counters (8 atm) 6 10 24
Table 1: BELEN detector characteristics (units in mm)

The version of the BELEN detector used in this experiment was specifically designed and optimized by means of Monte Carlo MCNPX Pelowitz et al. (2005); Fishman (1996); Rubinstein and Kroese (1981); Riego (2016) and Geant4 Agostinelli et al. (2003); Tarifeño-Saldivia, A. (2017) simulations in order to achieve a high and flat detection efficiency in the range from 0.12 MeV (see Fig. 3).

Figure 3: BELEN detector efficiency from the GEANT4 simulation: total efficiency and for each one of the 3 rings of He proportional counters (dashed lines) Tarifeño-Saldivia, A. (2017), and neutron spectra of the measured isotopes (from ENSDF nnd ()) as a function of the neutron energy (colored lines). See text and Table 1 for details.

The most important constraint was in terms of efficiency: for the detection of two correlated neutrons, the efficiency, , is roughly proportional to the square of the detection efficiency for a single neutron, :


To avoid a low detection probability the central hole of BELEN was kept at 6 cm diameter, restricted to the diameter size of the vacuum tube of the implantation system. With this geometry, BELEN reached a one-neutron efficiency of 60%, and 36% for a two-neutron event, in the Q energy range of the isotopes of interest. The disadvantatge of this optimization compared to other designs was the impossibility of placing a HPGe detector right behind the implantation point to verify the identity of the isotope implanted by -ray detection. This latter technique was used in the prior experiment during the same experimental campaign with another version of BELEN detector, measuring isotopes in the same region of the nuclear chart Agramunt et al. (2017). The one-neutron efficiency for that version of BELEN detector was 40% in the same energy range (0.12 MeV) due to a different geometry, while using the same number of He neutron counters. For the present work, the efficiency obtained in the simulations, shown in Fig. 3, was experimentally validated at E=2.13 MeV using a Cf source, and with beams of isotopes with well-known P values, such as Rb and I, as reported in Section III.

Previous experiments using the BELEN detector demonstrated that it works well with the self-triggered data acquisition system specifically implemented for BELEN, named GASIFIC Agramunt et al. (2016); Caballero-Folch et al. (2017). This system integrates all signals from the - and the He neutron counters recording their energy, and a time-stamp with a clock of 10 ns resolution to be able to build the -neutron time-correlations over a certain time-window. In this experimental campaign, differential to single-ended converter modules, designed at JYFL, were added in the electronic chain. This made possible to link the output signals from the MPR-16-HV Mesytec preamplifiers, directly connected to the He tubes, to the SIS3316 Str () sampling ADC modules in the data acquisition system Agramunt et al. (2017). This improvement enabled the recording of the data without requiring the signal shapers used in previous experiments.

Regarding the response of the neutron detection system, the energy spectrum obtained for the neutron events comprises the range from a low-energy threshold at 191 keV up to the 764 keV peak. This energy is attributed to the kinetic energy released in the reaction:


The 764 keV peak corresponds to the sum of the collection of the total energy released by the two reaction products, a triton and a proton. The lower energy detection threshold is related to the partial collection of the energy and the wall effect Ravazzani et al. (2006). Gain-matching with a calibrated Cf source was carried out for all 48 He counters before the experiment. The stability of the overall detector response was checked regularly during the experiment. The accumulated spectrum for all 48 tubes during the Sb measurement is shown Fig. 4. The latter includes the uncorrelated neutron events, which were removed in the data analysis process together with other light particles and the noise at lower energies.

Figure 4: BELEN energy spectrum acquired with the 48 He tubes for the Sb measurement.

Iii Determination of the neutron branching ratios

The equation that describes the total -delayed neutron branching probability, P, of a nucleus is defined as the sum of all, single and multiple, neutron emission contributions present in the decay:


For nuclei with only one-neutron emission energetically allowed, i.e. Q0 and Q0, the P value is equal to P. This is the case for the Rb, I, and Te isotopes measured in this experiment. In order to determine their P values the methodology successfully implemented in a previous experiment with a similar setup at the IGISOL facility was followed Agramunt et al. (2016), where the P values were deduced as:


where N is the number of the net and neutron time-correlated events, N the number of decays registered, is the mean efficiency, the averaged efficiency above the S weighted according to the neutron energy spectrum (see Fig. 5), and the neutron efficiency (see Fig. 3).

Figure 5: The MCNPX simulation of the counter efficiency (dashed line), and neutron energy spectra of the measured isotopes (colored lines). The lower efficiency at low energies affects the detection of correlated events for those isotopes with a low Q value (see text for details).

The analysis of the acquired data directly provides three parameters that are needed to evaluate and determine the neutron branching ratios. These are the number of particles and the number of neutrons detected, together with the time-correlation events. Figure 6 shows the analysis to determine the and neutron integrals for each measured isotope, Rb, I, Te, and Sb, from the growth-and-decay curves, using the Bateman equations Bateman (1910).

(a) Rb analysis
(b) I analysis
(c) Te analysis
(d) Sb analysis
Figure 6: Analysis of the and neutron events accumulated during the implantation-and-decay cycles in the growth-and-decay curves for each isotope measured. In the analysis of I (b) the fit function also includes a parameter to reflect the loss of the xenon isotopes in the chain. As xenon is a noble gas, a fraction of the nuclei rapidly escape from the implantation tape. Details of this effect are presented in Ref. Agramunt et al. (2016). In the case of Rb (a) the trap purification cycle can be identified. In all cases the ions were extracted out from the trap in bunches, once in every 181 ms according to the length of this cycle. This was repeated as many times as needed during the “beam on” time. So the beam was not continuous but rather a quasi-continuous.

The particular growth (implantation) and decay times in the curves were precisely set according to the half-lives, taken form the literature nnd (), of the isotopes present in the decay chain. The flat area at the beginning of the timing cycle (see fig. 6) corresponds to the characterization of the and neutron background. Table 2 provides the Q values of interest together with half-lives and the details of their implantation and decay cycle times. The numerical results from Fig. 6 for the precursor isotopes are shown in later tables.

Nuclei Rb I Te Sb
t(s) 0.3777(8) 24.5(2) 17.63(8) 0.923(14)
Q(keV) 9228(20) 6027(8) 5120(14) 9916(7)
Q(keV) 4881(20) 2002(8) 1292(6) 5147(7)
Q(keV) 1951(20) 6082(8) 6509(15) 1884(6)
Impl. time (s) 1.27 81.3 79.5 3.80
Decay time (s) 2.54 171.5 68.5 3.71
Table 2: Half-lives nnd () and Q values Audi et al. (2017) of the investigated isotopes, and characteristics of the implantation-and-decay cycle, corresponding to the growth-and-decay curves.

In the present work I and Rb beams were used as calibration isotopes because of their well-known P values Rudstam et al. (1993); Pfeiffer et al. (2002); Liang et al. (2018). In addition, these two isotopes have very different neutron energy spectra (see Fig. 3). In the case of Rb, its neutron energy spectrum is mainly centered at lower energies, with a large Q value of 4.8 MeV. This indicates that the lower detection probability at lower energies (see Fig. 5) does not miss correlated events in the analysis. Therefore, assuming in Eq. 4, it is possible to determine its P value to a first approximation with the equation:


The resulting P value for Rb, considering a constant value of 62.98%, taken from the simulation and weighted with the neutron spectrum, is 8.6(4)%, in agreement with the literature value of 8.94(37)% Liang et al. (2018). On the other hand, I has a neutron spectrum spread over the whole energy range up to its endpoint at 2 MeV (see Fig. 3). In this case its analysis is affected due to the non-constant efficiency at low energies (see Fig 5), and the simplified expression detailed in Eq. 5 is not accurate enough. Thus, we defined an expression to avoid the and dependence in Eq. 4 as in Ref. Agramunt et al. (2016), using the total number of detected neutrons, N, independently of the coincidence or not with the :


from which, it is then possible to determine the ratio with the well-known values of the Rb and I calibration isotopes. The ratio obtained can be used to determine the remaining P values in this experiment. Table 3 summarizes the integral values of the correlated events, the number of events, N, and the number of neutrons, N, for the Rb and I measurements, together with the calculated / ratio.

Nuclei P(%) N N N /
Rb 8.94(37) Liang et al. (2018) 33011 610229 90445 1.66(7)
I 7.66(14) Liang et al. (2018) 21888111Biased value due to the lower at lower energies. 592009 72031 1.59(3)
Table 3: The experimental values measured for Rb and I and their respective / ratio.

The average of the / ratio, calculated from Rb and I measurements in Table 3, is 1.62(7). As these two isotopes have a the large Q values this ratio value is expected to be alike. Thus, we assume that this value also applies for all the other isotopes in the present experiment.

iii.1 The P value of Te

With a Q value of 1292(6) keV, Te has an energy window that allows -delayed one-neutron emission. In this case the expected energy spectrum of the emitted neutrons is also affected by the non-constant at low energies described above for I. Thus, the equation to determine its neutron branching ratio can be defined by Eq. 6, using the / ratio determined with the calibration isotopes. With the integral values obtained in the analysis of the growth-and-decay curves in Fig. 6 and reported in Table 4, the analysis yields a P of 1.47(6)%. This value is slightly higher but in fair agreement with those reported in the literature and the IAEA evaluation: 1.31(5)% Rudstam et al. (1993); Liang et al. (2018) and 1.26(20)% Pfeiffer et al. (2002). Using the simplified Eq. 5, with the correlated detected events derived from the analysis of Fig. 7, the P value would be 0.90(5)%. Comparing this value with the one obtained using the efficiency ratio correction we observe an increase of about 63%. This gives an idea of the fraction of correlated events missed because of the lower at low energies when measuring an isotope with a neutron energy spectrum like Te (see Fig. 5). The data for this isotope in this measurement were recorded for 165 cycles which included a background measurement of 10 s, 79.5 s of implantation time, and 68.5 s of decay for each one, being 158.5 seconds per cycle (see Fig. 6 and Table 2), i.e. of 7.3 hours of beamtime. A total of 2082 net events were registered for this isotope.

Figure 7: Analysis of the correlation events for Te. According to the background level on the left side, the neutron moderation time in the polyethylene for this experiment was determined to be 400 s.

iii.2 P and P values of Sb

As described in Section I, when the state populated after the decay is higher than the two-neutron separation energy, S, in the daughter nucleus (Q), double neutron emissions are allowed. This is the case for Sb with Q=1884(6) keV. Its measurement in this experiment comprised 26465 recorded implantation-and-decay cycles of 9.51 seconds (3 days). As can be observed in Fig. 6, in each cycle the first 2 s were used to characterize the and neutron background. Table 2 details its implantation-and-decay cycle times. The neutron branching ratios of this isotope have two contributions to the aforementioned general P equation (Eq. 3): the P and the P values.

In order to determine its neutron branching ratios, several expressions based on the analysis were used. Eq. 7 describes the total number of neutrons, N, measured with the neutron background subtracted:


It is worth mentioning that the total number of neutrons has contributions from the one- and two-neutron emission branching ratios of Sb plus a small contribution from its daughter Te. The latter is present in the decay chain and the relevant number of neutrons can be determined with the information extracted from the Sb -decay analysis (Fig. 6) and its P value. In order to calculate the contributions of one-neutron and two-neutron events from the Sb decay in Eq. 7, we define the following expressions:




These expressions relate the number of events, N, and neutron events, N, and their detection efficiencies, to the unknown P and P values. However, from the available information, the analysis only provides the net number of N events for each specie in the decay chain, and the net number of the total neutron events, N, after the subtraction of the background. It does not distinguish between single- and double- neutron emission, and it is not possible to solve this using these equations alone. In order to address this issue we tried to evaluate the two-neutron correlated events from a background measurement with the Sb data. As shown in Fig. 8, this method is not useful because it is not possible to subtract the background and get the number of double neutron correlations, N, coming only from the decays of interest. Thus the P value had to be determined directly from the correlation events recorded.

Figure 8: (Color online) Background level evaluation of double neutron events, N, for the Sb setting. As can be observed, the number of events in a background measurement is of the same order as in the setting of interest. See text for details.

The large background can be associated with different kinds of interactions between particles, either cosmic or reactions of the beam with materials, and the BELEN detector.

Therefore, the strategy to determine the neutron branching ratios was based on the evaluation of the direct correlation events obtained in the data analysis, presented in Fig. 9.

Figure 9: (Color online) Sb correlation events. In blue are shown the first correlated neutrons after a decay, in red the second neutron in the correlation. A total of 49 net events were recorded.

The total number of net correlation events observed was 49. If we were to assume that the mean energy of the neutron spectrum is not affected by the lower at low energies, the P value could be calculated with a simplified expression such as Eq. 10, analogously to Eq. 5 used for the case of Rb for its P value:


With this assumption, the P value would be 0.22(4)%. However, due to the expected low in the energy region of interest for the two-neutron energy spectrum (Q value 1884(6) keV, see Fig. 5), it is necessary to apply the correction based on the (Eq. 11) in the same way as Eq. 4 for the case of one-neutron emission (see Section III):


The unknown neutron energy spectrum for the Sb two-neutron emission suggests the use of the calculated value for I, as its Q value, 2002 keV, is very close to the Sb Q value. So, using the / coefficient, 1.62(7), previously determined, the P value is calculated to be 0.31(5)%. Having obtained the P, the P value, has been determined to be 31.9(15)%.

Nuclei N N N N P(%) P(%)
Te 2082 - 360655 8645 1.47(6) -
Sb 9326222Includes twice the events. 49 57590 30455333Includes 2n events, and 1n events from Te 31.9(15) 0.31(5)
Table 4: Te and Sb neutron emission probabilitiy results.

Iv Summary and discussion

The P and P values obtained in this work are summarized in Table 5 together with the previously available experimental values and theoretical predictions from several models, including FRDM+QRPA Möller et al. (2003), DF3+cQRPA+RQRPA Borzov (2006), KTUY Koura et al. (2005), RHB+RQRPA Marketin et al. (2016), and QRPA-HF Mumpower et al. (2016b), and the semi-empirical effective density model Miernik (2014). The P value obtained for Te, 1.47(6)%, is higher but in fair agreement with those reported in Refs. Rudstam et al. (1993); Pfeiffer et al. (2002) and Liang et al. (2018). Concerning Sb, the analysis yields a P value of 31.9(15)%. This is higher than previous experimental results also taken from Refs. Rudstam et al. (1993); Pfeiffer et al. (2002) and Liang et al. (2018). In the case of the P we have obtained a value of 0.31(5)%. This is lower than the 1.4% estimated in a measurement with isobaric contamination Testov et al. (2011), and one order of magnitude lower than the predictions of the FRDM+QRPA and the DF3+cQRPA+RQRPA models.

Isotope t (s) nnd () P(%) P Ref. t(s) P(%) P(%) Ref.
(Experimental) (Theory)
Te 17.63(8) 1.47(6) - (This work) 10.166 2.43 0.0  Möller et al. (2003)
1.31(5) -  Liang et al. (2018) 21.189 1.80 -  Borzov (2013)
1.31(5) -  Rudstam et al. (1993) 0.600 0.41 -  Koura et al. (2005)
1.26(20) -  Pfeiffer et al. (2002) 0.548 0.4 0.1  Marketin et al. (2016)
1.7(8)444Updated values from the evaluation performed in Rudstam et al. (1993) -  Crancon et al. (1978) - 9 0.0  Mumpower et al. (2016b)
0.7(4) -  Rudolph et al. (1977) - 2.8 0.0  Miernik (2014)
Sb 0.923(14) 31.9(15) 0.31(5) (This work) 1.998 33.5 6.19  Möller et al. (2003)
- 1.4(2)555Measurement with isobaric contamination  Testov et al. (2011) 0.8 10.5(51) 4.15(1.05)  Borzov (2016)
18.7(18) 1  Liang et al. (2018) 0.760 33.20 0.0  Koura et al. (2005)
16.3(32) -  Rudstam et al. (1993) 0.175 3.8 0.2  Marketin et al. (2016)
23.2(68) -  Pfeiffer et al. (2002) - 40.0 1.0  Mumpower et al. (2016b)
44(57)4 -  Rudolph et al. (1977) - 37.3 0.0  Miernik (2014)
33(40)4 -  Crancon et al. (1978) 0.46 17.1 0.28  Lyutostansky and Panov (1983)
Table 5: Results

As shown in Table 5, the values of most of the theoretical predictions are far from the experimental values. Some of them are compatible for Sb but not for Te, and vice versa. With the current results it is not easy to decide which of the models reproduces better the experimental results. The DF3+cQRPA Borzov (2013, 2016) approaches better the magnitudes of the P and half-life for Te, but in the case of Sb, although the predicted half-life is close to the experimental value, it underestimates P value and overestimates P value. The model which better reproduces the Sb data is the KTUY Koura et al. (2005) but the half-life predicted for Te remains far from the experimental value.

To date, several calculations and estimates for multiple-neutron emission beyond A=100 have also been reported. Table 5 reports the available ones for the isotopes measured in this study. Some of them follow the main theoretical models describing the strength functions and some others are estimates and extrapolations like Miernik (2014). Concerning experimental measurements, there is a recent P value reported for Sb Moon et al. (2017). This value was determined through an indirect measurement based on relative intensities of -rays observed from transitions that were identified as belonging to Te. The value reported is “about” 8% and no uncertainty was given. This suggests again the need for direct neutron measurements in order to obtain more precise data on neutron emission branching ratios.

V Conclusions and outlook

We have determined the neutron branching ratios for Te and Sb through a direct neutron measurement. This represents the first experimental multiple-neutron emission value above A=100 and an improvement for the values available so far for these isotopes. As discussed above, some of the theoretical predictions agree well with the values obtained, but none shows agreement for all of the parameters present in the decay. This, together with the discrepancies with the experimental data available, indicates that more measurements with pure beams and direct neutron detection are needed in order to provide further input for the models in this region, and to study the nuclear properties above the neutron separation energies. Consequently, this experimental campaign in which six more isotopes included in the IAEA priority list for reactor physics Dillmann et al. (2014) were also measured Agramunt et al. (2017), and present and future campaigns like the BRIKEN Tarifeño-Saldivia et al. (2017) are an opportunity to increase the amount of data available in this field.

It is also of interest to determine whether the two neutrons are emitted simultaneously or sequentially in the -delayed two neutron-emission process. In the first case an angular correlation between both neutrons and the -ray emitted by the final nucleus would be expected. In the second case it would be necessary to correlate these events with very narrow time-windows to be able to confirm the phenomenon. In this work we were not able to address this question because of the moderation of the neutrons in the polyethylene matrix and the statistics available.

This work was supported by the National Research Council of Canada (NSERC) Discovery Grants SAPIN-2014-00028 and RGPAS 462257-2014 at TRIUMF, and by the Spanish Ministerio de Economia y Competitividad under grants: FPA2010-17142, AIC-D-2011-0705, FPA2011-28770-C03-03, FPA2011-24553, FPA2014-52823-C2-1-P, FPA2014-52823-C2-2-P and the program Severo Ochoa (SEV-2014-0398). It is also supported by the European Commission under the FP7/EURATOM contract 605203 and by the Academy of Finland under the Finnish Centre of Excellence Programme 2012–2017 (Project No. 213503, Nuclear and Accelerator-Based Physics Research at JYFL). A.K. also acknowledges the Academy of Finland grants No. 275389 and 284516. I.D. and M.M. acknowledge the support of the German Helmholtz Association via the Young Investigators Grant No. VH-NG 627 (LISA- Lifetime Spectroscopy for Astrophysics). W.G. acknowledges the support of the UK Science & Technology Faculties Council (STFC) under grant No. ST/F012012/1 and the University of Valencia.


  • Roberts et al. (1939) R. Roberts, R. Meyer,  and P. Wang, Physical Review 55, 510 (1939).
  • Goldansky (1960) V. Goldansky, Nuclear Physics 19, 482 (1960).
  • Lyutostansky et al. (1985) Y. S. Lyutostansky, V. Sirotkin,  and I. Panov, Physics Letters B 161, 9 (1985).
  • Audi et al. (2017) G. Audi, F. Kondev, M. Wang, W. Huang,  and S. Naimi, Chinese physics C 41, 030001 (2017).
  • Wang et al. (2017) M. Wang, G. Audi, F. Kondev, W. Huang, S. Naimi,  and X. Xu, Chinese physics C 41, 030003 (2017).
  • Dillmann et al. (2017) I. Dillmann, P. Dimitriou,  and B. Singh, Development of a Reference Database for Beta-Delayed Neutron Emission, Tech. Rep. (INDC(NDS)-0735 - IAEA, 2017).
  • Abriola et al. (2011) D. Abriola, B. Singh,  and I. Dillmann, Beta-Delayed Neutron Emission Evaluation, Tech. Rep. (INDC(NDS)-0599 - IAEA, 2011).
  • Burbidge et al. (1957) E. Burbidge, G. Burbidge, W. Fowler,  and F. Hoyle, Rev. Mod. Phys 29, 547 (1957).
  • Arcones and Martínez-Pinedo (2011) A. Arcones and G. Martínez-Pinedo, Physical Review C 83, 045809 (2011).
  • Mumpower et al. (2016a) M. Mumpower, R. Surman, G. McLaughlin,  and A. Aprahamian, Progress in Particle and Nuclear Physics 86, 86 (2016a).
  • Surman et al. (2014) R. Surman, M. Mumpower, J. Cass, I. Bentley, A. Aprahamian,  and G. C. McLaughlin, in EPJ Web of Conferences, Vol. 66 (EDP Sciences, 2014) p. 07024.
  • Caballero-Folch et al. (2016) R. Caballero-Folch et al., Physical Review Letters 117, 012501 (2016).
  • Möller et al. (2003) P. Möller, B. Pfeiffer,  and K.-L. Kratz, Physical Review C 67, 055802 (2003).
  • Borzov (2016) I. Borzov, Physics of Atomic Nuclei 79, 910 (2016).
  • Mumpower et al. (2016b) M. Mumpower, T. Kawano,  and P. Möller, Physical Review C 94, 064317 (2016b).
  • Lyutostansky and Panov (1983) Y. S. Lyutostansky and I. Panov, Zeitschrift für Physik A Hadrons and Nuclei 313, 235 (1983).
  • Testov et al. (2011) D. Testov et al., “Esp-rus congress,”  (2011), Talks_Presentations/Parallel_Sessions/Seminar_S214/09-11-11/Testov.pdf.
  • Äystö (2001) J. Äystö, Nuclear Physics A 693, 477 (2001).
  • Hakala et al. (2012) J. Hakala et al., Physical Review Letters 109, 032501 (2012).
  • Agramunt et al. (2016) J. Agramunt et al., Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 807, 69 (2016).
  • Karvonen et al. (2008) P. Karvonen, I. Moore, T. Sonoda, T. Kessler, H. Penttilä, K. Peräjärvi, P. Ronkanen,  and J. Äystö, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 266, 4794 (2008).
  • Nieminen et al. (2001) A. Nieminen, J. Huikari, A. Jokinen, J. Äystö, E. Collaboration, et al., in Atomic Physics at Accelerators: Mass Spectrometry (Springer, 2001) pp. 527–530.
  • Savard et al. (1991) G. Savard, G. Bollen, H.-J. Kluge, R. Moore, T. Otto, L. Schweikhard, H. Stolzenberg, U. Wiess, et al., Physics Letters A 158, 247 (1991).
  • Eronen et al. (2012) T. Eronen et al., in Three decades of research using IGISOL technique at the University of Jyväskylä (Springer, 2012) pp. 61–81.
  • Kolhinen et al. (2013) V. Kolhinen et al., Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 317, 506 (2013).
  • Rudstam et al. (1993) G. Rudstam, K. Aleklett,  and L. Sihver, Atomic Data and Nuclear Data Tables 53, 1 (1993).
  • Pfeiffer et al. (2002) B. Pfeiffer, K.-L. Kratz,  and P. Möller, Progress in Nuclear Energy 41, 39 (2002).
  • Liang et al. (2018) J. Liang et al., Submitted to Nuclear Data Sheets  (2018), Evaluation of Beta-Delayed Neutron Emission Probabilities and Half-Lives for Z 28.
  • Torner et al. (2014) A. Torner, J. Agramunt, A. Algora, L. Batet, R. Caballero-Folch, F. Calviño, D. Cano-Ott, A. García, G. Cortés, I. Dillmann, C. Domingo-Pardo, M. Gómez-Hornillos, V. Gorlychev, M. Marta, T. Martínez, A. Poch, C. Pretel, A. Riego,  and J. Taín, Technical Design Report of the Beta-Delayed Neutron Detector (BELEN) for NUSTAR (DESPEC), Tech. Rep. (Universitat Politècnica de Catalunya, Barcelona (Spain)., 2014).
  • Gómez-Hornillos et al. (2011) M. Gómez-Hornillos et al.Journal of Physics: Conference Series 312, 052008 (2011).
  • Gómez-Hornillos et al. (2014) M. Gómez-Hornillos et al.Hyperfine Interactions 223, 185 (2014).
  • (32) LND Inc,
  • Pelowitz et al. (2005) D. B. Pelowitz et al.MCNPX user’s manual version 2.5.0, Tech. Rep. (Los Alamos National Laboratory, 2005).
  • Fishman (1996) G. Fishman, Monte-Carlo concepts, Algorithms and Applications (1996).
  • Rubinstein and Kroese (1981) R. Y. Rubinstein and D. P. Kroese, Simulation and the Monte Carlo Method (1981).
  • Riego (2016) A. Riego, Design of the BELEN detector for wide energy range with flat and high detection efficiency, Ph.D. thesis, Universitat Politècnica de Catalunya, Barcelona (Spain). (2016).
  • Agostinelli et al. (2003) S. Agostinelli et al., Nuclear instruments and methods in physics research section A: Accelerators, Spectrometers, Detectors and Associated Equipment 506, 250 (2003).
  • Tarifeño-Saldivia, A. (2017) Tarifeño-Saldivia, A.,  (2017), (private communication).
  • (39) Evaluated Nuclear Structure Data File ENSDF.,,
  • Agramunt et al. (2017) J. Agramunt et al., in EPJ Web of Conferences, Vol. 146 (EDP Sciences, 2017) p. 01004.
  • Caballero-Folch et al. (2017) R. Caballero-Folch et al., Phys. Rev. C 95, 064322 (2017).
  • (42) Struck Innovative Systems,
  • Ravazzani et al. (2006) A. Ravazzani, A. F. Para, R. Jaime, M. Looman, M. M. Ferrer, P. Peerani,  and P. Schillebeeckx, Radiation measurements 41, 582 (2006).
  • Bateman (1910) H. Bateman, Proc. Camb. Philos. Soc. 15, 423 (1910).
  • Borzov (2006) I. Borzov, Nuclear Physics A 777, 645 (2006).
  • Koura et al. (2005) H. Koura, T. Tachibana, M. Uno,  and M. Yamada, Progress of theoretical physics 113, 305 (2005).
  • Marketin et al. (2016) T. Marketin, L. Huther,  and G. Martínez-Pinedo, Physical Review C 93, 025805 (2016).
  • Miernik (2014) K. Miernik, Physical Review C 90, 054306 (2014).
  • Borzov (2013) I. Borzov, “-decay and delayed neutron emission of very neutron-rich nuclei,” in Fission and Properties of Neutron-Rich Nuclei (World scientific, 5International Conference on Fission and properties of neutron-rich nuclei(ICFN5), 2012, 2013) pp. 530–537, + (private communication).
  • Crancon et al. (1978) J. Crancon, C. Ristori, H. Ohm, W. Rudolph, K.-L. Kratz,  and M. Asghar, Zeitschrift für Physik A Hadrons and Nuclei 287, 45 (1978).
  • Rudolph et al. (1977) W. Rudolph, K.-L. Kratz,  and G. Herrmann, Journal of Inorganic and Nuclear Chemistry 39, 753 (1977).
  • Moon et al. (2017) B. Moon et al., Physical Review C 95, 044322 (2017).
  • Dillmann et al. (2014) I. Dillmann, P. Dimitriou,  and B. Singh, Development of a Reference Database for Beta-delayed Neutron Emission Evaluation, Tech. Rep. (INDC(NDS)-0643 - IAEA, 2014).
  • Tarifeño-Saldivia et al. (2017) A. Tarifeño-Saldivia et al., Journal of Instrumentation 12, P04006 (2017).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description