Activation cross section data of deuteron induced nuclear reactions on rubidium up to 50 MeV

Activation cross sections of the natRb(d,xn)87m,85m,85g,83,82Sr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\mathrm {nat}}\hbox {Rb (d,xn)}^{87\mathrm {m,85m,85g,83,82}}\hbox {Sr}$$\end{document}, natRb(d,x)86,84,83,82mRb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\mathrm {nat}}\hbox {Rb(d,x)}^{\mathrm {86,84,83,82m}}\hbox {Rb}$$\end{document} and natRb(d,x)85mKr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\mathrm {nat}}\hbox {Rb(d,x)}^{85\mathrm {m}}\hbox {Kr}$$\end{document} nuclear reactions have been measured for the first time through an activation method combining the stacked foil irradiation technique and gamma-ray spectrometry. The provided cross sections from the present investigation are all new, in such a way contribute to the completeness of the experimental database. The experimental cross sections were compared with the theoretical prediction in the TENDL-2019 TALYS based library and with our calculation using ALICE-D and EMPIRE-D model codes in order the improve their predictivity. Thick target production yields were calculated form the new cross sections for all investigated radioisotopes. Practical applications of the results are shortly discussed.


Introduction
Integral excitation functions for the production of residual nuclides through light charged particle activation constitute basic data for various applications. Nowadays, deuteron induced reactions start to play an important role as the stripping process generates high production yields and the stopping power for deuterons is relatively low. Some years ago, to meet requirements of practical applications, we started to establish an experimental activation database by performing new experiments and a systematic survey of existing data of deuteron induced cross-sections up to 50 MeV. This database is essential for accelerator and target technology to produce high-energy, high-intensity neutron fluxes, wear and a e-mail: ditroi@atomki.hu (corresponding author) material studies through TLA (Thin Layer Activation), medical radioisotope production, space applications (resistance of electronics, shielding, etc.), monitoring of deuteron beam intensities and energies, etc.
Unlike for proton induced reactions, the status (reliability and completeness) of the experimental data for deuteron induced reactions was rather poor, especially above [15][16][17][18][19][20] MeV. We embarked hence on a systematic program of extending the study to more targets and up to 50 MeV incident deuteron energy. In this work we present the excitation functions for reactions on rubidium (stable isotopes 85 Rb and 87 Rb) leading to radionuclides of Sr, Rb and Kr.
The measured excitation functions are compared with the results of three nuclear reaction model codes. These comparisons can show the present status of the predictivity of these codes and contribute to their development. The reliability of the presently used theoretical codes for deuteron induced reactions is low, compared to proton and alpha particle induced reactions, due to the modeling problems of the deuteron stripping and pickup.
Concerning activation cross sections of deuteron induced reactions on Rb, no experimental data were found in the literature. Only integral physical yields at 22 MeV deuteron energy were reported by Dmitriev et al. [1] for production of 84,86 Rb and 85 Sr.

Experiments and data evaluation
The main experimental parameters and data evaluation methods used in this work are shown in Table 1. The cross section data were measured by using the activation method, stacked foil irradiation technique and high resolution gamma-ray Fig. 1 Excitation function of the simultaneously measured 27 Al(d,x) 24,22 Na monitor reactions spectrometry. Beam current and energy scale are based on monitor reactions, re-measured simultaneously over the whole covered energy range (Fig. 1).
The Rb-targets were obtained by deposition of Rb 2 SO 4 (99.8%, Sigma Aldrich) using a sedimentation method [2] on a 50 μm thick, high purity Al backing in an 11 mm diameter spot. Actual thickness of the layer was determined by measuring the area and differential weighing of the material deposited on the Al backing. The pellets were covered by a 10 μm Al foil for protection. The metal foils were high purity products from Goodfellow.
The stack was irradiated with deuterons of 50 MeV incident energy at an external beam line of the Cyclone 110 cyclotron of the Université Catholique in Louvain la Neuve (LLN): (protons 20-80 MeV; d, α 2.3-27 MeV/nucleon; heavy ions, radioactive ions).
The present stack contained 19 blocks of Al-Rb 2 SO 4 -Al targets alternated with Al-RuCl 2 -Al targets and further 6 blocks of Al-Rb 2 SO 4 -Al targets alternated with Al-RuCl 2 -Al targets combined with 20 μm Ti monitors. The nominal current during the 40 min irradiation was 27 nA. The targets were mounted in a target holder (Faraday cup like) provided with a long collimator defining a 5 mm beam diameter.
Gamma-spectra were measured with Canberra GX1520 HPGe detectors of 15% efficiency and 1.9 keV resolution at the 1332 keV 60 Co peak at the VUB (Free University of Brussels) laboratory. The detectors were vertically arranged with adjustable distance sample holders from the detector surface to several meters. The detector head was put into 10 cm lead shielding with 2 mm copper internal lining. The individual foils or a group of target foils were measured in such a distance that the dead time was always under 5%. To make the evaluation easier not more than 5 distances were selected. The samples were measured first after a short cooling time of several hours for the short-lived components up to after several months for the longer-lived radioisotopes. The measuring time was chosen to acquire good statistics for the expected gamma-peaks, but of course it cannot be optimized for every radioisotope in question. The used decay data, taken from the online version of NUDAT 2.6 [3], and the energy threshold values in MeV, obtained from the Q-value calculator [4], are presented in Table 1.
The median beam energies in the individual targets were preliminarily determined by a degradation calculation (see [5]) and were corrected on the basis of the fitted monitor reaction [6]. The activities of the individual radioisotopes were determined from the measurement of the corresponding gamma-peaks (Eq. 1): where A E O B is the activity at the End of Bombardment, T is the net area of the gamma-peak, λ is the decay constant, t real , t live and t c are the real-time, live-time of the measurement and the cooling time, I γ is the gamma-line abundance (intensity) and ε is the detector efficiency. From the measured activities the corresponding cross section was calculated according to Eq. 2.
where z is the projectile charge, I (A) is beam current, N A is Avogadro's number (6.02214 × 10 23 mol −1 ), M (g/mol) should be molar mass of the chemical compound used as the target material, e is the electron charge (1.6x10 −19 C), ρ (g/cm 2 ) is the mass density of the target material, s (cm) is the thickness of the measured target foil, ν is the number of the particular element atoms in the compound molecule and f is the abundance of the particular isotope (entering the particular nuclear reaction) in the element. σ (E) (cm 2 ) is the cross section at bombarding energy E, where E is the median energy in the measured foil of the stack. By using both Eqs. 1 and 2 one must take care of the measuring units. Equation 2 is valid for isotopic cross sections. In our case mainly elemental cross sections were determined, where f = 1.
Uncertainties of the median energies were obtained taking into account cumulative effects of possible uncertainties (primary energy, target thickness, energy straggling, correction to monitor reaction). Uncertainty of cross-sections was determined by considering the sum in quadrature of all independent contributions following the ISO recommendations (cf. [7]): beam current (7 %), target thickness (6 %), detector efficiency (5 %), nuclear decay data (3 %), full energy peak area determination and counting statistics (1-20 %).
Natural Rb has two stable isotopes: 85 Rb (72.17 %) and 87 Rb (27.83 %). The experimental parameters and the data evaluation methods are summarized in Table 1.

Nuclear model calculation
The cross sections of the investigated reactions were calculated using the pre-compound model codes ALICE-IPPE [13] and EMPIRE-II [14]. The experimental data are also compared with the cross section data in the TENDL-2019 on-line library [15]. The TENDL-2019 library contains the results of the blind and adjusted calculations on the basis of the latest version (1.6) of the TALYS nuclear model code system [16]. For showing the evolution of this code and database also the results of the TENDL-2017 version are indicated and nearly no or only marginal changes in the results can be noticed.
In our previous works ALICE-IPPE and EMPIRE-II were used successfully for the description of a large amount of reaction cross sections induced by light charged particles. However, during the recent analyses of the (d,p) reactions on the isotopes, 114 Cd [17]; 169 Tm [18]; 192 Os [19] and some others, we were confronted with a large underestimation of the measured cross sections. We come to the conclusion that the experimentally observed cross sections of the (d,p) reaction cannot be reproduced below 20-30 MeV with the available statistical model codes. To achieve a better description of available data for (d,p) reactions with the ALICE and EMPIRE a phenomenological simulation of direct (d,p) and (d,t) transitions were introduced to the above-mentioned codes. A phenomenological enhancement factor K in these relations was taken as energy dependent and estimated to describe the whole set of the observed (d,p) cross sections for medium and heavy nuclei. By this improvement, in the ALICE IPPE-D and EMPIRE-D code versions for the deuteron induced reactions, the direct (d,p) channel is increased strongly and this is reflected in changes for all other reaction channels in both codes. As ALICE-IPPE calculates only the total cross section, for estimation of isomeric states from the ALICE code the isomeric ratios calculated by EMPIRE-D were applied.

Cross sections for production of strontium radioisotopes
The radioisotopes of strontium are produced via direct (d,xn) reactions on one or both of the stable Rb isotopes.

nat Rb(d,xn) 87m Sr
The isotope 87 Sr has a stable ground state and a rather shortlived isomeric state 87m Sr (T 1/2 = 2.815 h) that was measured with poor statistics because of the long cooling time. It is produced in our experimental circumstances only through the 87 Rb(d,2n) reaction. The agreement of our indicative results (large uncertainties, huge difference in decay corrections between different target foils) with predictions of the theoretical codes is acceptable (Fig. 2). There is no visible difference between the two versions of TENDL.

nat Rb(d,xn) 85m Sr
The radionuclide 85 Sr has a longer-lived ground state (T 1/2 = 64.849 d) and a shorter-lived isomeric state (T 1/2 = 67.63 min) decaying partly to the ground state. The excitation function for production of the metastable state is shown in Fig. 3. There is a good agreement with the TENDL-2019 (TENDL-2017) predictions at lower energies, representing the 85 Rb(d,2n) contribution. But around the weak maximum of the 87 Rb(d,4n) reaction the predictions of the three codes  To determine the activity of 85g Sr (T 1/2 = 64.849 d) with gamma-spectrometry, it is necessary to separate the 514 keV gamma-line from the strong annihilation peak. The measured cross sections include the complete decay contribution from the short-lived metastable state (T 1/2 = 67.63 min, IT 83.6 %) (see above). Comparing with code results the magnitudes are similar, but the shapes differ significantly (Fig. 4).

nat Rb(d,xn) 83g Sr(m+)
As it is obvious from Table 2, the very short-lived isomeric state of 83 Sr (T 1/2 = 4.95 s, IT: 100 %) decayed completely to the ground state (T 1/2 = 32.41 h) before the spectra measurements started. The experimental and theoretical data for the excitation curve of 83g Sr(m+) are systematically shifted in energy and differ in maximum cross section value (Fig. 5).

nat Rb(d,xn) 82 Sr
The longer-lived radionuclide 82 Sr (T 1/2 = 25.35 d) has no gamma-lines. The activity was measured through the 776 keV gamma-line of its daughter 82g Rb (T 1/2 = 1.2375 min). The comparison of the experimental and the theoretical data (in our experimental conditions only contribution of 85 Rb(d,5n) 82 Sr) are presented in Fig. 6, showing an excellent agreement with EMPIRE-D while an energy shift with the other codes appear.

Calculations for strontium isotopes
The statistical codes EMPIRE-D and TALYS describe rather well the near-threshold parts of the neutron emission reactions. The ALICE-D results are less satisfactory for this region through overly simplified modeling of the low-lying discrete levels of the considered nuclei. For calculations of the cross sections in the region of maximums, a ratio of the level densities of competing reaction channels is crucial. The existing differences in calculations reflect differences of the default level-density parameters of the corresponding codes. Discrepancies between calculations with experiment can be

Excitation functions for production of rubidium radioisotopes
The radioisotopes of rubidium are produced via direct (d,pxn) reactions and through the ε and/or β + decay of Sr parent radioisotopes.

nat Rb(d,x) 86g Rb (m+)
The longer-lived ground state 86 Rb (T 1/2 = 18.642 d) is produced via direct (d,pxn) reaction on 85,87 Rb. The measured data are so called (m+), including the complete decay of the isomeric state (T 1/2 = 1.017 min, IT 100 %) (Fig. 7). The description of 87 Rb(d,p) part is not satisfactory in all model codes. The experimental data are in good agreement with the systematics for a (d,p) reaction in this mass region [20].

nat Rb(d,x) 84g Rb(m+)
The 84g Rb (T 1/2 = 32.82 d) and 84m Rb (T 1/2 = 20.26 min, IT 100 %) isomers are produced only directly on nat Rb targets via 85 Rb(d,p2n) and 87 Rb(d,p4n) reactions. The cross sections for production of 84g Rb were deduced from spectra measured after complete isomeric decay of the shortlived metastable state. Among the model predictions only the EMPIRE-D descriptions are acceptable over the energy region studied (Fig. 8).

nat Rb(d,x) 83 Rb(cum)
The measured production cross section of 83 Rb (T 1/2 = 86.2 d) are cumulative. It contains the complete decay of its 83 Sr (T 1/2 = 32.41 h, ε : 100 %) parent isotope (Fig. 9). The model predictions are very different and cannot even describe rather well the cumulative maximum. EMPIRE-D is acceptable up to 38 MeV, and the cumulative TENDL version give the better results above this energy.

Calculations for rubidium isotopes
The general comments on the (d, pxn) reactions are the same as those discussed in Sect. 4.1.6 above. However, a special consideration is required for the (d, p) reaction. It is well known that this reaction is dominated by the mechanism of a direct breakup without the formation of any intermediate preequilibrium state of the nucleus. A satisfactory description of the experimental data on this reaction can be achieved only on the basis of the phenomenological systematics proposed  first(early) series of gamma spectra measurements and with low statistics. We had possibility to deduce the cross section data for 85m Kr via subtracting the contribution of the 85m Sr. The contribution was significant and taking into account that the cross section of the 85m Kr is low, the statistics of the separated 151 of the 85m Kr was poor. The experimental data, slowly rising as a result of the combination of the two contributing reactions without pronounced maximum in the studied energy region, and their curve is described acceptable well by the used model codes (Fig. 11).

The thick target and integral yields of the investigated radioisotopes
The thick target yields as a function of the bombarding energy, calculated from integration of a fit to our experimental excitation functions (Eq. 3), are shown in Figs. 12 and 13 compared to values obtained in direct thick target measurements [1]. All thick target and integral yields are deduced for metallic Rb targets of natural composition. The integral yields represent so called physical yields (see [11] and [12]). The only experimental values are for 85g Sr and 86g,84g Rb available [1] and shows acceptable agreement with our curve at 22 MeV deuteron energy.
Where E th is the threshold energy of the nuclear reaction in question, and F is the conversion constant, which contains all the constants as density, molar weight, decay constant, electron charge, Avogadro's number, etc. and the conversion factors necessary to get the integral yield (Y (E)) in proper units. Instead of integral along a numerical integration was performed along depth with Δx(E,dE), where dE is equidistant, and in such a way Δx is not.

Applications of deuteron induced reactions on rubidium for nuclear medicine
Many investigated activation products are relevant for nuclear medicine ( 87m Sr, 85g Sr, 83g Sr, 82 Sr ( 82 Rb), 82m Rb, 85m Kr). Through the presently investigated reactions, the production of 82m Rb is carrier added and the production yields of 85m Kr and 82 Sr (in the investigated energy range) are very low. Only production of the 87m Sr, 85 Sr, 83 Sr and 82 Sr radioisotopes are hence discussed shortly in more detail. More detailed discussion, including production parameters (yield, impurities, targetry, etc.) will be prepared and published in a dedicated journal.  [2] or proton induced reaction on enriched 87 Sr target [29,[34][35][36][37], proton induced reaction on 88 Sr target [29,[38][39][40] or proton induced reaction on natural strontium target [41][42][43] or deuteron induced reactions on nat Sr [44,45]. In practical medical isotope production, the role of the 3 He induced reaction is minimal due to lack of accelerator and high price of 3 He and the high stopping power. The situation is practically same in case of alpha particle induced reactions, except the price of 4 He gas. Alpha particles are used only when no low energy possibility with protons and deuterons is available ( 211 At). In practice proton induced reactions are the favorites due to the large cross sections, availability of high intensity beams, and lower stopping power. But during our systematic study of deuteron induced reactions it was shown that some isotopes cannot be produced by protons, but they are available with deuterons via (d,n) and (d,p) reactions, and at medium and heavy mass targets the production yields of (d,2n) are significantly higher comparing to (p,n). The 85g Sr (T 1/2 = 64.849 d) had importance in nuclear medicine as diagnostic radioisotope, later it was proposed to use it for the endotherapy and as important disturbing impurity in the application of the 82 Sr and 83g Sr. It can be produced directly or indirectly through 85 Y parent. In the literature the following direct and indirect production routes have been found. Direct production routes: 85 Rb(p,n) 85 Sr [29,[46][47][48], nat Rb(p,xn) 85 Sr [22,24,[49][50][51][52][53], nat Rb(d,xn) 85 Sr (this work), nat Kr(α, n) 85m Sr [26,54], nat Kr(α, n) 85g Sr [26,54], nat Kr( 3 He, n) 85m Sr [26,54], nat Kr( 3 He, n) 85g Sr [26,54]. Indirect production routes: 86 Sr(p,2n) 85 Y [2,29], nat Sr(d,xn) 85 Y [45].

Production of 83 Sr
The 83g Sr is a positron emitter analog of the β − emitting 89 Sr (T 1/2 = 50.5 d) used for endotherapy.
The following production routes have been investigated in the literature: nat Rb(p,xn) 83 Sr [22,24,49], 85 Rb(p,3n) 83 Sr [29,48,55]   ature data for production of 84,86 Rb and 85 Sr. The status of cross section and yield data for production of some medically important radioisotopes is also discussed. The obtained experimental data also provide a basis for improved model calculations and for applications in radioisotope production.
Funding Open Access funding provided by ELKH Institute for Nuclear Research Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: Cross section data is given in the manuscript in table form, both cross section and yield data will be deposited in the IAEA EXFOR database after publication, but is a question of time.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.