Modeling of deuteron-induced reactions on molybdenum at low energies

.


I. INTRODUCTION
Accurate cross-section data of deuteron induced reactions are a pre-requisite also for reliable design of the accelerator-based neutron source facility IFMIF-DONES (International Fusion Materials Irradiation Facility -DEMO Oriented NEutron Source), e.g.[1,2] and Refs.therein, and material qualification of the European DEMO fusion reactor [3] beyond the ITER fusiondevice materials [4,5].Actually, the suitable account of deuteron-nucleus interactions is an important test for reaction mechanism models and nuclear-data evaluation within on-going strategic research programmes at largescale facilities using deuteron beams [6][7][8][9][10].The present work concerns in this respect the deuteron-induced reactions on Mo, which is used in stainless steels to provide a greater corrosion resistance [11,12].SS-316L steel, being a reference material for Li loop piping of DONES contains 2-3% of Mo [13].Since the natural Mo consists of seven stable isotopes, with abundances between 9-24%, the activation of number of enriched foils in separate experiments would be preferred but for obvious reasons less feasible.Thus, a new experiment performed with nat Mo targets in Center of Accelerators and Nuclear Analytical Methods (CANAM) infrastructure [14] of the Nuclear Physics Institute of the Czech Academy of Sciences (NPI CAS) and Neutrons for Science Facility at SPIRAL-2 [8] was less resource-demanding while the analysis (both experimental [10,15] and theoretical) becomes a quite complex task.Apart from the applications in fusion technology, knowledge of activation cross- * marilena.avrigeanu@nipne.rosections of deuteron-induced reactions on molybdenum is very important also for production of the medically relevant radionuclides (e.g.[16][17][18]).
The deuteron sub-library of the TALYS-based evaluated nuclear data library (TENDL) [19], related to the output of the TALYS nuclear model code system [20], is the current reference within the EUROfusion Consortium [21] project of the Early Neutron Source (ENS) [22].However, the need of its further improvement has been pointed out [6].The latest version of the Fusion Evaluated Nuclear Data Library (FENDL) [23,24] also includes for structural materials [25] the data of the first version of FENDL, which are from TENDL-2011.Meanwhile, recent advancements in deuteron reaction modeling in the TALYS code [26,27] are taken into account to provide more reliable data for TENDL evaluated files, to be reviewed and integrated step by step in the future [28].
The recent studies of deuteron activation of elements within the priority list of the structural materials [29,30] known to be problematic (e.g.[31]), namely Al, V, Cr, Mn, Fe, Co, Ni, Cu, Zr, and Nb [32][33][34][35][36][37][38][39][40], are continued by the present work on Mo.These nuclear-model analyses of deuteron-induced reactions have been completed with reference to the short reaction-times direct interactions (DI), i.e. breakup (BU) and direct reactions (DR), in addition to the statistical processes of pre-equilibrium emission (PE) and evaporation from compound nucleus (CN) at continuously increasing reaction times.A careful account of BU and DR contributions seems essential in this respect due to the still existing disagreement between recent measurements and calculated data [16,17], while due consideration has also been given to the weak points of an eventual evaluation.Moreover, the consistent analysis of all available data for competitive reac-tion channels as well as full stable isotopic chains and even neighboring elements, beyond a particular reaction analysis for one or only a couple of isotopes (e.g.[41]), strengthens the assessment of the model approach.
A consistent energy-dependent optical model potential (OMP) for deuterons on Mo isotopes is concerned in Sec.II A using the computer code SCAT2 [42].Then, theoretical framework of the deuteron BU mechanism (Sec.II B) and the DR analysis (Sec.II C) involving the computer codes TALYS-1.96[20] and FRESCO [43], respectively, is briefly mentioned.The PE and CN mechanism contributions to the population of various residual nuclei (Sec.II D), on the grounds of also TALYS-1.96, are completing the nuclear model discussion of the present work.The measured and calculated deuteron activation cross sections of natural Mo and its stable isotopes, as well as the corresponding evaluated data within the TENDL-2021 library [19], are then compared in Sec.III. Conclusions of this work are given in Sec.IV.

II. NUCLEAR MODEL FRAMEWORK
The weak binding energy of the deuteron is responsible for the high complexity of its interaction with nuclei, involving also a variety of reactions initiated by the nucleons following the deuteron breakup (breakup nucleons) [27,44,45].Its importance increases with the targetnucleus mass and charge, becoming dominant for heavy target nuclei at incident energies particularly around the Coulomb barrier [46,47].Overall, the accompanying breakup-nucleons induced reactions make substantially different the deuteron projectile among the other incident particles, with proved significant impact even for applications in nuclear transmutation of radioactive waste [48,49].Additional comments on systematic uncertainties in DR theories used to analyze particularly (d, p) reactions are referred to in Sec.II C.

A. Deuteron optical potential assessment
A simultaneous analysis of the deuteron elastic scattering as well as induced activation cross sections is essential as the deuteron OMP parameters are obtained by the former data fit, and then used within the analysis of the latter ones.The demand of a consistent input of nuclear model calculations of deuteron activation of Al, V, Cr, Fe, Co, Ni, Zr, Nb, Pa, U [32-40, 46, 47] was satisfied by the deuteron OMP of Daehnick et al. [57].It has been a first option of the present work, being obtained by use of a large experimental basis including the angular distributions of elastic scattered deuteron on the above-mentioned elements as well as on 100 Mo isotope.
Unfortunately, the systematics of elastic-scattering angular distributions of deuterons on the stable isotopes of Mo is scarce, despite the number of these isotopes.There is also no related measurement of the deuteron total-reaction cross section σ R .The comparison of the measured angular distributions of elastic-scattered deuterons on 92,94,96,98,100 Mo [50][51][52][53][54][55][56] at the incident energies of 14.5, 17.2, 21.14, 21.5, and 22 MeV, and the calculated values using the computer code SCAT2 [42] and Daehnick et al. [57] OMP is shown in Fig. 1.The good description of the measured data provided thus confidence in the further use of this potential within a consistent analysis of the deuteron activation cross sections of Mo isotopes through the BU, DR, PE, and CN reaction mechanisms.

B. Deuteron breakup
A detailed overview of the specific deuteron breakup in the Coulomb and nuclear fields of the target nucleus, including its complexity given by the addition to the primary deuteron-nucleus interaction of various nuclear reactions initiated by the breakup nucleons, was given recently [27,39].Therefore, only particular points are mentioned hereafter concerning the two distinct BU processes, i.e. the elastic breakup (EB) in which the target nucleus remains in its ground state and none of the deuteron constituents interacts with it, and the inelastic breakup or breakup fusion (BF), where one of these constituents interacts nonelastically with the target nucleus.
Empirical parametrization [27,44,45,58] has concerned the fractions f n/p BU and f EB of cross sections for the total breakup (EB+BF) nucleon-emission σ n/p BU and elastic breakup σ EB , respectively, to σ R .The experimental systematics of deuteron-induced reactions on target nuclei from 27 Al to 232 Th [59][60][61][62][63] was involved in this respect.Moreover, since equal BF nucleonemission cross sections σ n BF and σ p BF were assumed [63], the total breakup (BF+EB) cross section becomes σ BU =2σ n/p BF +σ EB , while the total breakup (BF+EB) nucleon-emission cross section is σ n/p BU =σ n/p BF +σ EB .At the same time, the BF fractions for each of the breakup nucleons are given by the difference An overall view of the incident-energy dependence of the BU, BF, and EB cross sections for deuterons on 92,100 Mo target nuclei is shown in Fig. 2(a,d).It is thus apparent the change of these excitation functions with the target-nucleus mass number, too.There are shown also both the BU and BF nucleon-emission excitation functions.These results as well as the earlier ones [32-40, 46, 47] point out the dominant role of the BF component, to be concerned in connection with the two opposite effects of the deuteron breakup on the deuteron activation cross sections [27,44,45].At first, the total-reaction cross section shared among various outgoing channels, is reduced by the value of the total breakup cross section σ BU .Next, this component brings contributions to different reaction channels through breakup-neutrons/protons interactions with target nucleus [27, 32-40, 46, 47]   tions, respectively.
It should be noted that the BF cross-section partition among various residual-nuclei population is triggered by the energy spectra of the breakup nucleons as well as the excitation functions of the reactions induced by these nucleons on the target nucleus [44,45], while the atomic mass and maybe also the atomic number of the related compound nuclei differ by one unit from that in deuteroninduced reactions.The BF related formalism involved in the present work is described in Sec.3.4 of Ref. [27], the corresponding reaction cross sections being calculated with the code TALYS-1.96[20] and its second option for the breakup model, i.e. including the BF enhancement ( [20], p. 40).
The enhancements due to (p, x) and (n, x) reactions induced by the breakup nucleons on 92,94,96,98,100,nat Mo are discussed in Sec.III and distinctly showed in Figs.9-23.It has been thus proved that BF enhancements are particularly important for the suitable description of the maximum as well as the high-energy side of the excitation functions for second and third chance emitted-particle channels [27].

C. Direct reactions
Apart from the breakup contributions to deuteron interactions, an increased attention has to be paid to the direct reactions so far poorly considered within deuteron activation analysis.For low-and medium-mass target nuclei and deuteron energies below and around the Coulomb barrier, the interaction process proceeds largely through DR mechanism, while PE+CN component become important with the incident-energy increase.Nev-ertheless, the assessment of DR cross sections is subject to available information on spectroscopic factors for populated states, outgoing particle angular distributions, or at least differential cross-section maximum values.
The appropriate assessment of stripping -(d, p), (d, n), and ( 3 He, d) -and pick-up -(d, t), (d, α), and (d, 3 He) -DR mechanism contributions has been performed through the distorted-wave Born approximation (DWBA) method, with the code FRESCO [43].The post/prior form distorted-wave transition amplitudes for stripping and pick-up reactions, respectively, and the finite-range interaction have been considered.A note may concern the finite range (FR) effects in deuteron stripping process which were formerly shown to be very small around the incident energy of 10 MeV, so that a zero-range (ZR) approximation was expected to be reasonable [64] and eventually completed (e.g.[65]) by an usual FR correction [66,67].More recently, ZR calculation has led to results very close to the full FR-DWBA calculation for deuteron reactions while FR effects become important for incident 6 Li ions [68].Nevertheless, it was also found these effects are important in (d, p) reactions at intermediate energies (20 MeV/u) [69] while they are not however among the actual major sources of uncertainties in a variety of reaction theories used to analyze (d, p) nuclear reactions in the incident energy range 10-20 MeV [70].
At the same time, a simple Gaussian form was considered for the n-p interaction in the deuteron [71], which reproduces well the binding energy of deuteron and the low-energy n-p scattering phase shifts.It is yet widely used (e.g.[68,72]), DWBA as well as BU cross sections being modified by a few percent when including a realistic NN interaction for the reaction at low en- n/p BF (dashed curves), and elastic breakup σEB (dash-dot-dotted curves) cross sections [45]; (b,e) DI excitation functions (thick solid curves) and its components: total BU (thin curves), DR (dashed curves), stripping (d, p) (dash-dotted curves) and (d, n) (dash-dot-dotted curves), and pick-up (d, α) (short-dotted curve) and (d, t) (dotted curves); (c,f) σR fractions of BU (thin solid curves), DR (dashed curves), DI (solid thick curves), and PE+CN (thick dotted curves) cross sections of deuteron interactions with (a,b,c) 92 Mo and (d,e,f) 100 Mo (see text) ergy, while at higher energy the cross section is not affect by details of the NN interaction [73][74][75][76][77].Following the original data analysis, similar form has been used for the d-n interaction in the triton [78], at the same time with a Woods-Saxon shape of the d-p interaction [71] in 3 He, as well as the d-d interaction [79] in the α particle.The transferred-nucleon and -deuteron bound states were generated in a Woods-Saxon real potential [32,36].
The number of N nodes corresponding to an L transferred angular momentum in the radial wave function was determined by the harmonic-oscillator energy conserva-tion rule: where n n/p and l n/p are the single-particle (neutron/proton) shell-model state quantum numbers.Particularly for the FRESCO code, the number of nodes includes the origin, so that N >0.The energy, spin, and parity of the residual-nuclei discrete levels within the ENSDF [90] and RIPL [91] libraries were used as starting input of the DWBA calculations for the stripping and pick-up reactions [92][93][94][95][96][97][98].
The analysis of the neutron angular distributions from 92,94 Mo(d, n) 93,95 Tc reactions [80] is shown in Fig. 3 for transitions to discrete levels of residual nuclei.There were thus provided the spectroscopic factors (given in the Appendix A) and then the corresponding calculated stripping cross sections, e.g.Fig. 2(b), Fig. 13(e,j,o), and Fig. 15(d,h,l) in Sec.III.
( 3 He, d) reactions.There are so scarce or even lacking measured neutron angular distributions for 95,96,98 Mo(d, n) 96,97,99 Tc stripping reactions, in order to obtain the stripped-proton spectroscopic factors for the assessment of the DR contributions Because of that, these spectroscopic factors have been got through deuteron angular-distribution analysis of ( 3 He, d) stripping reaction populating the same residual nucleus as the (d, n) process.Thus, description of the deuteron angular distributions from 95,96,98 Mo( 3 He, d) 96,97,99 Tc stripping reactions [81][82][83] (Figs. 4 and 5) validates the extracted stripped-proton spectroscopic factors.Then, they were involved in 95,96,98 Mo(d, n) 96,97,99  by the more recent analyzes of 92 Mo(d, p) 93 Mo [100] (Table III), and 98,100 Mo(d, p) 99,101 Mo reactions [99] (Tables XXIV and XXIII, respectively).The related Even worse than the scarce systematics of neutron angular distributions corresponding to (d, n) stripping processes, there exist measured α-particle angular distributions only for the 92 M o(d, α) 90 Nb pick-up reaction [84].Their analysis shown in Fig. 6 has allowed to extract the spectroscopic factors of the picked deuteron.The corresponding (d, α) excitation function is included also in Fig. 2(b), but particularly in Fig. 22(l) of Sec.III.On the other hand, spectroscopic factors corresponding to 94,97,98 Mo(d, α) 92,95,96 Nb pick-up reactions have been obtained following due consideration of the measured αparticle angular distributions [101]  (d, t) reactions.The analysis of measured triton angular distributions [85][86][87][88] of the 92,94,98,100 Mo(d, t) 91,93,97,99 Mo pick-up reactions has provided the spectroscopic factors for the corresponding transitions to discrete levels of residual nuclei.The suitable account of the available (d, t) angular distributions shown in Fig. 7 provides the modeling support of the calculated pick-up excitation functions shown in Figs.2(b,e) and 17(c).
(d, 3 He) reactions.The analysis of measured outgoing 3 He angular distributions of the 96,98 Mo(d, 3 He) 95,97 Nb pick-up processes [89] made possible further calculations of the corresponding pick-up cross-section contributions to the total activation cross sections of 95,97 Nb residual nuclei.Thus, appropriate description of these angular distributions (Fig. 8) validated the extracted spectroscopic factors which have led to pick-up excitation functions in Figs.20(f) and 21(c,h,m) of Sec.III.
The DR contributions, which have proved essential for describing the measured excitation functions corresponding to the first-chance emitted particles [32][33][34][35][36][37][38][39][40], are furthermore discussed in Sec.III.Nevertheless, it should be pointed out the significant effect of the maximum of the stripping excitation functions around 8-12 MeV, shown in Fig. 2(b,e), on the increasing low-energy side and the maximum of the corresponding total activation excitation functions.
To sum up, the deuteron total-reaction cross section that remains available for the PE+CN mechanisms has   97 Tc (bottom) [82] stripping transitions, at 33.6 and 18 MeV, respectively followed the due consideration of the incident flux leakage through the breakup, stripping and pick-up DI processes.
Using the present work notation as well as its first assessment of the deuteron break-up and DR transitions, the composite-formation cross section, e.g.Eq. (164) [20], becomes σ P E +σ CN =σ R -σ BU -σ DR .In brief, similarly to the depletion factor to account for the direct and PE ef-fects within the default TALYS computation, Eq. (127) [20], we have used finally a normalization factor of the deuteron total-reaction cross section, of the form: for the final assessment of the PE+CN cross section.The (deg) for an isotopic chain is triggered by the increase of the BU mechanism importance with increase of the targetnucleus mass number.Overall, the DI and statistical processes fractions vary around half of σ R [34,36,37,39] along the actual incident-energy range, pointing out the quite important DI role.

D. Statistical emission
The statistical PE+CN mechanisms, which complete the deuteron interaction analysis along an enlarged nuclear-interaction time scale, become important with the increase of the incident energy above the Coulomb barrier.The corresponding reaction cross sections have been calculated using TALYS-1.96code [20], taking into account also the above-discussed BU approach as well as the stripping and pick-up mechanisms through the nor-     [57], Becchetti-Greenlees [103], and Avrigeanu et al. [104] for neutrons and protons, deuterons, tritons, and α-particles, respectively, (b) the model for breakup reactions including the BF enhancement [27,45], (c) the back-shifted Fermi gas formula for the nuclear level density (NLD), and (d) the PE exciton model with analytical transition rates with energy-dependent matrix element, and spin distribution based on particle-hole state densities.Furthermore, likewise to the analysis of 87 Y m,g activation by deuterons incident on nat Zr [39], better account of the measured ground and isomeric state excitation functions of 91 Nb, 93 Mo, and 93,94,96,97 Tc residual nuclei has been obtained within actual state-of-art TALYS calculations by amending the spin distribution cut-off parameter by a factor of 0.25.
Next to the above-mentioned sharing of around half of the deuteron total-reaction cross section by the PE+CN mechanisms (Fig. 2), the weight of every reaction mechanism is concerned in the following analysis of the available data for deuteron-induced reactions on 92,94,95,96,97,98,100,nat Mo at energies up to 50 MeV.

III. RESULTS AND DISCUSSION
The excitation functions of residual nuclei from deuterons interaction with nat Mo are compared in Figs.9-23 with the available data corresponding also to irradiation of all Mo stable isotopes [50], the similar TENDL-2021 evaluation [19], and the results of calculations using FRESCO [43] and TALYS [20] codes.A particular attention is given, beyond the more recent data [16-18, 106, 107], to new measured excitation functions at NPI CAS in Řež [15] for the activation of 97 Nb, 98 Nb m , 97 Tc m , and 101 Tc residual nuclei.
The higher number of Mo stable isotopes, i.e. 92 Mo (14.7%), 94 Mo (9.2%), 95 Mo (15.9%), 96 Mo (16.7%), 97 Mo (9,6%), 98 Mo (24.2%) and 100 Mo (9.7%), made necessary the analysis of their detailed contributions to the activation of a certain residual nucleus.The involved reaction mechanisms are particularly illustrated in order to compare their strength.While some contributions are not shown for all Mo stable isotopes in Figs.9-23, they were all considered within the cross-section calculation for the natural Mo target.
The activation of a certain residual nucleus is presented following four reaction sequences: nat Mo(d, xn) 93−101   88,89 Zr (Fig. 23).The following discussion concerns firstly the heavier residual nuclei since fewer stable isotopes of Mo contribute to their population, while more isotopes, processes, and emitted particles are involved going towards the lighter products.The analysis of 101 Tc residual nucleus activation by deuterons incident on nat Mo concerns only the deuteron interaction with the neutron-richest stable isotope 100 Mo.Moreover, 101 Tc excitation function has been for the first time reported [15] for a natural molybdenum target [Fig.9(a)], in addition to the previously measured 100 Mo(d, n) 101 Tc activation cross sections [105] [Fig.9(b)].In fact, there are two ways of 101 Tc residual nucleus population, i.e. via (d, n) reaction and β − decay of 101 Mo (T 1/2 =14.61 min) which is activated through nat Mo(d, p) reaction (Sec.III B 1).The NPI-reported 101 Tc activation cross sections [15] have been corrected for the 101 Mo decay contribution, while the same correction was noticed also by Randa et al. within their reported data for the activation of 100 Mo isotope.Moreover, amending Randa et al. data by 100 Mo natural abundance, it results a good agreement of the two data sets below the incident energy of 10 MeV.However, the maximum of the wholly resulting excitation function in Fig. 9(a) is maybe followed by a too sharp decrease at the higher incident energies of the NPI measurement [15].
The calculated cross sections for 100 Mo(d, n) 101 Tc reaction have proved a significant DR stripping (d, n) contribution only for incident energies above 20 MeV, becoming larger than the statistical PE+CN from ∼30 MeV [Fig.9(b)].At the same time the inelastic breakup enhancement through (p, γ) reaction remains much lower in the whole energy range.Overall, it results a notable underestimation of the measured cross-section maximum [105], at the energies where PE+CN component is prevailing.A rather similar trend has the TENDL evaluation [19], particularly at the energies where the DR weight is larger.
On the other hand, the addition of the 101 Mo decay has led for the cumulative calculated cross sections to an improved account of the reported 100 Mo(d, n) 101 Tc excitation-function maximum [105].Its further decrease is however less quick, as it is that of the cumulative TENDL-2021 evaluation in Fig. 9(b), too.
Nevertheless, the present approach has provided the rather good agreement of the newly measured nat Mo(d, n) 101 Tc cross sections shown in Fig. 9(a).This suitable description stands also as a proof of the appropriate correction for the removal of the 101   The nat Mo(d, xn) 99 Tc m activation cross sections measured most recently at NPI [15] are in the good agreement with the recent data [16,17,106] shown in Fig. 10(a).It should be noted again that the population of 99 Tc m is also cumulative, including the β − decay of 99 Mo (T 1/2 =65.9 h) activated through 98 Mo(d, p) reaction (Sec.III B).The appropriate correction has been carried out within the NPI-measured values [15], as well as reported by Lebeda et al. [106], Elbinawi et al. [17], and Tárkányi et al. [16] on both nat Mo and 100 Mo [18].The same correction was noticed by Randa et al. [105] for 98 Mo(d, n) 99 Tc m reaction, but not by Zarubin et al. [108] for a similar measurement [Fig.10(b)].
The calculated results of this work for the nat Mo(d, xn) 99 Tc m activation are in a rather good agreement with all data available until an incident energy of ∼20 MeV [Fig.10(a)].More important within this energy range is firstly the (d, n) reaction on 98 Mo nucleus which, however, becomes half of order of magnitude lower than that of (d, 3n) reaction at 20 MeV.The statistical PE+CN component is the major one for the former reaction, below ∼30 MeV, while then this role is taken by the stripping DR [Fig.10(b)].It is not the same case of the latter reaction, where BF contribution increases significantly at higher energies but yet below the PE+CN one shown in Fig. 10(c).
Thus, it has been of interest to understand a sizable underestimation of 98 Mo(d, n) 99 Tc m data of Zarubin et al. [108] at the same time with a suitable account of the similar data of Randa et al. [105] in Fig. 10(b).Addition   96 Mo(d, n) 97 Tc m , (c) 97 Mo(d, 2n) 97 Tc m , (d) 98 Mo(d, 3n) 97 Tc m , and (e) 100 Mo(d, 5n) 97 Tc m reaction cross sections, with (a) contributions of (d, n), (d, 2n), (d, 3n), and (d, 5n) reactions (dash-dot-dotted, dash-dotted, dashed, and thin solid curves, respectively) of the above-mentioned decay of 99 Mo residual nucleus, in a similar way to 100 Mo(d, p) 101 Tc reaction analysis, has led anyhow to a cumulative excitation function in agreement with the former data, too.
However, a suitable account of the measured 100 Mo(d, 3n) 99 Tc m excitation function [18] is found only until ∼30 MeV, the arrow in Fig. 10(c) pointing out the threshold of this reaction.Moreover, an overestimation is then apparent by even the statistical PE+CN component alone.The BF enhancement being already comparable with the experimental data, its strong importance becomes obvious for the third-chance particle emission channel, in comparison to (p, γ) BF added to (d, n) oneparticle emission.Overall, the model calculations lead to better data account versus TENDL evaluation while a particular extension of the measurements to incident energies over 30 MeV is highly requested.The measured nat Mo(d, xn) 97 Tc m excitation function has recently been reported for the first time [15] while no measurement exists for 96,97,98,100 Mo isotopes involved in this respect.Its account by the present calculations is re-lated firstly to results for the (d, n) and (d, 2n) reactions on 96 Mo and 97 Mo isotopes, respectively [Fig.11(a)].The (d, 3n) reaction channel for 98 Mo target nucleus becomes also significant close to the incident energy of 20 MeV, while the (d, 5n) reaction on 100 Mo is relevant for 97 Tc m activation only above 50 MeV.
A routine comment may concern the PE+CN major contribution to the population of 97 Tc m .However, the (d, n) stripping component is shown in Fig. 11(b) to be even more important at incident energies over ∼25 MeV.These DR results are, at their turn, exceeded at the low energies ≤9 MeV by the BF contribution through the (p, γ) reaction even though followed by a fast decrease.
Concerning the involved reaction mechanisms, PE+CN plays an important role except the (d, n) reaction [Fig.12(b,g,l)].In a similar way to the previous analysis of 97 Tc m activation, there is a strong competition of the DR stripping and PE+CN mechanisms beyond, e.g. 10 MeV.The former becomes dominant at higher energies.The BF enhancement, due to the (p, γ) reaction induced by breakup protons, also exceeds the stripping contribution at low incident energies, although it decreases faster.
Moreover, there is a sound increase, by about two order of magnitudes, of the BF contributions through (p, n) and (p, 2n) reactions [Fig.12(c,d)].They become comparable with the statistical PE+CN at energies around 50 MeV.The (p, 3n) BF contribution to the (d, 4n) forthchance particle emission is strong, too [Fig.12(e,j,o)], but yet lower than (p, n) and (p, 2n) reactions.Nevertheless, the lack of data related to specific Mo isotopes beyond the only one set corresponding to 95 Mo(d, n) 96 Tc reaction [105] restrains further comments on the theoretical approach.
5. nat Mo(d, xn) 95 Tc g+m,g,m The 95 Tc g+m,m activation cross sections include significant contributions of 94,95,96,97 Mo isotopes, through reactions from (d, n) to (d, 4n) respectively.The weight of the BU, DR, PE, and CN mechanisms, shown for each of these reactions in Fig. 13(b-e,g-j,l-o), may prove an appropriate account and validation of the present analysis for these reactions.
However, similarly to the 96 Tc g+m,g,m activation, very scarce related data of specific Mo isotopes exist only for 95 Tc g,m populated through (d, n) reaction [105,111].The present model analysis describes the measured excitation function for the latter state but overestimate the former [Fig.13(g,l)].At the same time, amending 94 Mo(d, n) 95 Tc g activation cross sections [111] by natural abundance of 94 Mo isotope, a discrepancy between them and nat Mo(d, xn) 95 Tc g data systematics becomes apparent in Fig. 13(f).
Again, the PE+CN dominant role within 95 Tc g+m,g,m population at lower energies is visible at incident energies above 20-30 MeV [Fig.13(c-e,h-j,m-o  elastic breakup contribution, through the (p, γ) reaction, and the DR stripping into the (d, n) outgoing channel is also obvious in Fig. 13(b,g,l).The former exceeds the latter one at low incident energies, with an excitation function maximum higher by an order of magnitude at ∼8 MeV.Then, a faster decrease is leading to the opposite case around 50 MeV.So, at the same energy, the DR contribution is however well below the PE+CN one.
Once more, the need to enlarge the data systematics related to Mo stable isotopes activation, particularly at higher energies, should be pointed out.The actual model analysis shows a very good agreement with the measured excitation function for 94 Tc m activation [Fig.14(i)].However, the 94 Tc g activation data are well described just above the threshold as well as at incident energies above 20 MeV, while there is an overestimation between 10-20 MeV in Fig. 14(e).
As already noticed above, scarce measurements for activation of specific Mo isotopes consist of two earlier data sets corresponding in Fig. 14(f,j) to 94 Mo(d,2n) 94 Tc g,m reactions [105,111].Unfortunately, they were carried out only at the lower incident energies ≤12 MeV.The present model calculations describe the latter [105] but slightly overestimate the former [111].
It is thus proved that similar overestimation of the 94 Tc g activation is due to the (d, 2n) reaction on 94 Mo, despite the lowest natural abundance of this isotope.One may further note the statistical nature of the related difference between experimental and calculated results, the only additional BF contribution through the (p, n) reaction being lower by roughly one order of magnitude.Certainly, any further measurements of 94 Tc activation on specific Mo isotopes as well as measured excitation functions of nat Mo(d, xn) 94 Tc g+m,g over 20 MeV incident energy would be most helpful for detailed analysis of these discrepancies.The model calculations describe the excitation functions of 93 Tc g+m,g,m activation by deuterons on natural Mo, as well as the data corresponding to 94 Mo(d, n) 93 Tc m reaction [105] in Fig. 15

(j).
There is however an overestimation of the measured 92 Mo(d, n) 93 Tc g+m,g excitation functions [105] in Fig. 15(b,f) below an incident energy of 12 MeV.Moreover, amending the activation cross sections of the 92 Mo(d, n) reaction [105] by the natural abundance of 92 Mo, the corresponding derived data in Fig. 15(a,e) are lower, too, than existing systematics for natural Mo target.
It should be noted also in this case the BF enhancement of the (d, n), (d, 3n), and (d, 4n) cross sections through the (p, γ), (p, 2n), and (p, 3n) reactions, respectively, induced by the breakup protons.However, 12 but for (a-d) 94 Tc, (e-h) 94 Tc g , and (i-l) 94 Tc m activation [15,16,50,105,106,111], with (a,e,i) contributions of (d, 2n), (d, 3n), and (d, 4n) reactions (dash-dotted, dashed, and thin solid curves, respectively).Similarly to the case of 101 Tc residual nucleus, the activation of 101 Mo follows the deuteron interaction with only the neutron-richest 100 Mo stable isotope.The model calculations, taking into account the BU, DR, and statistical PE+CN mechanisms, describe rather well the two experimental excitation functions [15,106] shown in Fig. 16(a).There are also measured activation cross sections corresponding to 100 Mo(d, p) 101 Mo reaction [112], slightly overestimated in Fig. 16(b).Amending these cross sections by the natural abundance of 100 Mo isotope, it was obtained a derived data set also lower than both newer measurements in the former figure.
The stripping (d, p) process provides now a much strong contribution to 101 Mo activation comparing with (d, n) one shown in Fig. 9(b) for 101 Tc activation.Actually, the description of nat Mo(d, p) 101 Mo excitation function, with contributions only from a single Mo isotope, is relevant for the importance of the stripping (d, p) mechanism.This proof is just in line with the previous analyses of the experimental (d, p) excitation functions which can not be described as long as the substantial DR stripping contribution is neglected [33-37, 39, 40].Moreover, the apparent underestimation by TENDL evaluation of both nat,100 Mo(d, p) 101 Mo reactions (Fig. 16) follows entirely the lack of taking into consideration the DR stripping process.

nat Mo(d, xnp) 99 Mo
The analysis of 99 Mo activation by deuterons interaction with natural Mo and its stable isotopes 98,100 Mo could take advantage of the existing systematics for Mo target [15-17, 50, 106, 107, 109] shown in Fig. 17(a) as well as the measured excitation functions for 98 Mo [108,112] and 100 Mo [18,113] target nuclei, respectively [Fig.17(b,c)].In fact, it involves two direct reactions of stripping 98 Mo(d, p) 99 Mo and pick-up 100 Mo(d, t) 99 Mo, along with the breakup and PE+CN statistical emission.Their weight is however quite distinct for each of the two isotopes and even within different energy ranges, so that their specific discussion is required.
Thus, the DR stripping mechanism prevails over PE+CN in the case of 98 Mo(d, p) 99 Mo excitation func-tion above the incident energy of 15 MeV, while at lower energies they are close to each other.The measured 98 Mo(d, p) 99 Mo data [108,112] are thus described rather well [Fig.17(b)], while the BF enhancement through breakup-neutron induced reaction 98 Mo(n, γ) 99 Mo is essentially negligible.Unfortunately, lack of data at energies above 15 MeV makes not possible a more detailed analysis.
Comparatively, the BF contribution through (n, 2n) and (p, d) reactions is shown in Fig. 17(c) to be strong for 100 Mo(d, 2np) 99 Mo reaction.It is crossing the pick-up (d, t) excitation function at ∼20 MeV, becoming higher by one order of magnitude at ∼30 MeV.However, it should be pointed out this pick-up reaction as the only contribution to 99 Mo activation by deuterons on 100 Mo at the incident energies ≤ 14 MeV.Nevertheless, it is by an order of magnitude lower than the (d, p) stripping reaction on 98 Mo, and thus not visible in the case of 99 Mo activation by deuterons on the natural Mo target.
Finally, the sum of BF, DR, and PE+CN contributions overestimates the dispersed 100 Mo(d, 2np) 99 Mo activation data [16,113].Consequently, a similar overestimation exists above 30 MeV for the nat Mo(d, p) 99 Mo excitation function [Fig.17(a)].On the other hand, the suitable account of this excitation function up to ∼20 MeV does validate the (d, p) stripping modeling.New measurement at higher energies may show if there is a major question of the actual analysis in this respect.The analysis for activation of the high-spin 21 + /2 isomeric state of 93 Mo, at 2425.9 keV excitation energy, takes advantage of a close agreement [Fig.18(a)] between the NPI recent measurement [15] and the data obtained by Tárkányi et al. [16] within their latest irradiation in 2007, during fifteen years (1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007) studies at different accelerators.In fact, a large spread between these results and the earlier ones (1995) [16] has been just at the energies of the latest measurement [15], while more approaching data were found at higher energies.Otherwise, a model analysis is restrained by the lack of any measurement of this activation within deuteron-induced reactions on either 92 Mo, 94 Mo, or 95 Mo isotopes that may have major related yields.
Furthermore, a complementary modeling should take into consideration the DR mechanism for the striping 92 Mo(d, p) as well as pick-up 94 Mo(d, t) reactions, in addition to also more complex BF enhancements and, finally, the statistical PE+CN emission.However, there is no measured angular distribution corresponding to this 93 Mo isomeric state either in (d, p) stripping, or in (d, t) pick-up reactions, in order to extract the corresponding spectroscopic factors of the stripped/picked neutron needed for the calculations of DR corresponding yield.Moreover, the calculation results for the total stripping and pick-up cross sections, shown in Fig. 2(b), do not in-  clude excited states of 93 Mo which may populate 93 Mo m through their decay.Therefore, no DR component has been possible to be added to the BF+PE+CN calculated cross sections in Fig. 18(b,c).Consequently, an underestimated excitation function could be expected at low incident energies in the absence of both stripping and pick-up contributions.
At the same time it should be noted that the model results shown in Fig. 18 have been obtained using an additional factor of 0.25 for the NLD spin cut-off parameter of the residual nucleus 93 Mo.Comparing with the analysis of 87 Y g,m activated by deuteron interaction with nat Zr target [39], there is no measurement of 93 Mo g excitation function to confirm this spin cut-off adjustment on the basis of the corresponding g/m ratio.This adjustment could be the reason of the apparent differences between the present and TENDL-2021 predictions in Fig. 18, too.
Nevertheless, one may note that the large amount of measured isomeric cross sections for the 93 Nb(d, 2n) 93 Mo m [40] and 94 Mo(n, 2n) 93 Mo m [114] reactions was already well described using the same breakup approach and no spin cut-off adjustment but the particle-hole state density (PSD) [115,116] within the geometry-dependent hybrid PE model [117].This PSD formula has included the PE spin distribution that was early discussed by Feshbach et al. [118] and further detailed by Fu [119].These results are in line with the more recent conclusion that reduced values of the spin cut-off parameter, obtained from isomeric cross-section analysis, were artificial and resulted from the improper use of the CN spin distribution also for the PE spin distribution [120].The present analysis describes the data of Tárkányi et al. [16] in Fig. 19(a).While only the 92 Mo target nucleus [Fig.19(b)] plays a role in this respect, versus 94 Mo isotope [Fig.19(c)], a contribution due to the decay of 90 Tc residual nucleus (T 1/2 =50.7 s) may also be noted.It follows the (d, 4n) reaction on 92 Mo, being however lower by at least two orders of magnitude than 90 Mo activation data and the corresponding cumulative excitation  90 Mo activation [16], with contributions of the (d, 3np) and (d, 5np) reactions (dashed and thin solid curves, respectively), as well as (a,b) 90 Tc decay (dotted) to cumulative population (thin curves + cross) of 90 Mo function.
Moreover, there are no measurements of the 90 Mo activation by deuterons incident on the distinct 92,94 Mo isotopes.Thus, the analysis corresponding to the natural Mo is the only way for validation of the contributions from BF and PE+CN mechanisms, the latter being shown in Fig. 19(b,c) to be obviously the dominant one.
The TENDL-2021 predictions in Fig. 19(a,b) are also cumulative, including the 90 Tc decay.However, the underestimation by the TENDL-evaluation at higher energies may be just due to the yet neglected contribution of the inelastic breakup enhancement.A rather equally-well account of the excitation functions measured either recently [15] for nat Mo(d, xn2p) 98 Nb m reaction [Fig.20(a)] or formerly for 100 Mo(d, 2n2p) 98 Nb m reaction by Tárkányi et al. [18] [Fig.20(c)] has been shown by the actual analysis.Actually, only the two heaviest Mo isotopes may contribute to the activation of this isomeric state, through 98 Mo(d, 2p) 98 Nb m and 100 Mo(d, 2n2p) 98 Nb m reactions.The latter reaction channel corresponds, especially at the energies lower than the (d, 2n2p) reaction threshold, also to the (d, α) reaction with a maximum around the incident energy of 16 MeV.
Because of the significantly higher effective threshold of the (d, 2p) reaction on 98 Mo [Fig.20(b)], the 98 Nb m activation for the natural Mo target comes fully from this (d, α) reaction.The obvious correspondence between the two above-mentioned data sets is emphasized within the actual BF+PE+CN analysis except the less usual steep increase of the latter excitation function at incident energies over ∼45 MeV.However, this model analysis has been proved difficult due to the lack of any spectroscopic information concerning the pick-up 100 Mo(d, α) 98 Nb m reaction.Therefore, beyond the dominant PE+CN statistical mechanisms, only the inelastic breakup has been involved in this particular case.Its correspondence to a breakup-neutron induced reaction (n, 2np) led to crosssection contribution at incident energies only above 35 MeV and too low to account for the above-mentioned increase in Fig. 20(c).
action by breakup neutrons through 98 Mo(n, p) 98 Nb m reaction [Fig.20(b)].It exceeds by three times the PE+CN contribution for the (d, 2p) reaction on same target nucleus.This proton-emission enhancement is essential for the analysis of deuteron radiation-damage estimation in the structural materials, due to the additional hydrogen accumulation known as the "gas bubble accumulation".
2. nat Mo(d, xn2p) 97 Nb Similarly to the above, there is a suitable account of 97 Nb activation data first-time recently reported [15] for deuterons incident on natural Mo [Fig.20(d)] as well as of the earlier measurement [18] of the (d, 3n2p) reaction on 100 Mo [Fig.20(g)].There is also a correspondence in the last case with a first excitation-function maximum around the incident energy of 23 MeV due to a (d, nα) reaction component with an obviously related low-energy side of the data for the natural Mo.
There is also a straightforward strong increase of the BF enhancement, being by far the dominant mechanism for 97 Mo(d, 2p) 97 Nb reaction [Fig.20(e)].It overcomes the PE+CN contribution in 98 Mo(d, n2p) 97 Nb process too, for the incident energies ≥45 MeV shown in Fig. 20(f).
In fact, the availability of spectroscopic data related to the 98 Mo(d, 3 He) 97 Nb pick-up reaction has been followed by a DR cross-section assessment.It has been thus proved the dominant role of this mechanism at incident energies ≤25 MeV, while it is still higher than the BF component up to ∼37 MeV [Fig.20(f)].The analysis of the 96 Nb activation is of prime interest for the present analysis due to measured excitation function by Tárkányi et al. [18], with the lower incident energies ≤20 MeV corresponding fully to the (d, α) reaction on the most abundant isotope 98 Mo [Fig.20(h,k)].There are also concerned contributions from 96,97,98,100 Mo isotopes [108,121], through various BF, DR, PE, and CN mechanisms [Fig.20(h-l)].
The model calculations describe much better the more recently measured nat Mo(d, xn2p) 96 Nb data [18] than the earlier 98 Mo(d, 2n2p) 96 Nb excitation function [108,121] at the energies of the only (d, α) reaction.There are thus validated the major CN+PE statistical emission as well as the smaller but yet notable DR component in Fig. 20(k).
Nevertheless, a main point is still the BF mechanism, with the (n, p) enhancement bringing the dominant contribution to the 96 Mo(d, 2p) 96 Nb activation.It completes thus the 98,97 Nb activation in Fig. 20(b,e,i).The basic role of inelastic breakup and pick-up reaction contributions is confirmed also by the data underestimation of TENDL-2021 predictions [19] which are yet not including the account of these processes.
The same is true for the data reported more recently by Tárkányi et al. [16], including an additional excitation function of nat Mo(d, xn2p) 95 Nb g activation [Fig.21(f)].A weak contribution of 95 Zr decay (T 1/2 =65.9 h), visible only in Fig. 21(a), has been considered too, with no real difference for the 95 Nb g+m,g cumulative activation.
There are three points of this analysis which demands additional comments.
First, the most abundant isotope 98 Mo brings the dominant contribution to 95 Nb g+m,g,m activation by (d, 3n2p) reaction on natural Mo above the deuteron energy of ∼16 MeV.A first maximum of the (d, 2n2p) reaction, at lower energies, corresponds in fact to the one-particle emission (d, α) reaction on the nearly least abundant isotope 97 Mo [Fig.21(d,i,n)].Then, the two-particle emission (d, nα) sort of (d, 3n2p) reaction on 98 Mo has an excitationfunction maximum at 22-25 MeV [Fig.21(e,j,o)].At the same time there are earlier data sets [121,122] of (d, α) reaction on 97 Mo which are largely overestimated.However, it should be noted that the closely related recent data for natural Mo are well accounted for.
Second, it could be again pointed out the important contributions of BF enhancement to (d, 2p) reaction cross sections, for the 95 Mo target nucleus [Fig.21(b,j,l)], similarly to the above-discussed 96,97,98 Nb activation.Thus, the BF component exceeds the PE+CN statistical emission at the deuteron energies above 35 MeV.The latter becomes however dominant with the increase of additional neutron emission in (d, xn2p) reactions with x≥2.This is true also for deuterons on 100 Mo, with results shown only in Fig. 21(a,f,k) but similar mechanism contributions as for 98 Mo isotope.Third, cross-section increases due to the (d, 3 He) and (d, α) pick-up DRs are noteworthy particularly for the former.This one is even larger than CN+PE component, at the lower energies ≤13 MeV, and BF below 30-34 MeV [Fig.21(c,h,m)].However, the TENDL-2021 evaluation [19] has followed the low-energy side of the earlier measured excitation function [106,121,122] in Fig. 21(a,d,n).Nevertheless, its underestimation of the nat Mo(d, xn2p) 95 Nb g,m reaction cross sections could be related to the missing of BF and DR pick-up contributions.FIG.22: (Color online) As Fig. 21 but for (a-e) 92 Nb m , (f-j) 91 Nb m , and (k-n) 90 Nb activation [15,16,50,106,121], with (a,f,k) contributions of (d, 2p), (d, n2p), (d, 2n2p), (d, 3n2p), (d, 4n2p), and (d, 5n2p) reactions (short-dotted, short-dotted, dash-dotdotted, dash-dotted, dashed, and thin-solid curves, respectively), and (k) 90 Mo decay (short-dash-dotted) The nat Mo(d, xn2p) 92 Nb m activation cross sections, measured recently at NPI with higher accuracy [15], as well as previous similar data [16,106] are rather well described by the actual model analysis except the lowest deuteron energies below 4 MeV [Fig.22(a)].A similar agreement exists for an earlier measurement [121] at lower energies of the (d, α) reaction in Fig. 22(c), as the first maximum of the 94 Mo(d, 2n2p) excitation function, apart from data at ∼12 MeV.Nevertheless, there is a valuable account of the whole nat Mo systematics at the same energy [Fig.22(a)], with the largest contribution of 94 Mo isotope in spite of its smallest natural abundance.The next important and main contribution up to the deuteron energy of 20 MeV, belongs to 95 Mo isotope.It corresponds also to a first maximum of the two-particle (d, nα) side of the (d, 3n2p) reaction channel [Fig.22(d)].
At the same time, the isotope 92 Mo has the main role within 92 Nb m activation above 20 MeV, through the (d, 2p) reaction.The corresponding cross sections are again much increased by the BF enhancement due to the 92 Mo(n, p) reaction induced by breakup neutrons [Fig.22(b)].This BF component becomes even higher than the CN+PE statistical processes at incident energies above 35 MeV.It is however overridden, up to the same energy, by the pick-up (d, α) reaction on 94 Mo nucleus [Fig.22(c)].Despite an increased uncertainty above its effective threshold, the available data of 91 Nb m activation [15,16] are rather well described taking into account the contributions of four isotopes and three reaction mechanisms [Fig.22(f-j)].The dominant one in the whole energy range comes from the semi-magic nucleus 92 Mo due to a BF enhancement close to main CN+PE statistical emission in Fig. 22(g).It is exceeded only at energies below 20 MeV by the first maximum of the (d, 3n2p) excitation function for the target nucleus 94 Mo, corresponding to the sequential emission of only a neutron and an αparticle through actually a (d, nα) reaction (at variance with the sequential five-nucleon emission).
Unfortunately, there is no measurement of this activation on any of the Mo stable isotopes.Overall, the results of this analysis for 90 Nb activation are in rather good agreement with the data measured for natural Mo and 92 Mo targets [15,16,106,121].A note should concern the main contribution which is by far given by the two maxima of 92 Mo(d, 2n2p) 90 Nb excitation-function [Fig.22(l)].The former, corresponding to the α emission channel, is also additionally increased by the (d, α) pick-up DR, while the BF has the  22(k-m) but for (a-c) 89 Zr, (d-f) 89 Zr m , and (g-i) 88 Zr activation [15,16,50,106], with contributions of (a,d,g) the 94,92 Mo target isotopes (dashed and thin-solid curves, respectively), as well as (a,d) the 89 Mo and 89 Nb decay (dotted and short-dotted curves respectively), and (g) the 88 Nb decay (short-dotted curve) same effect for the latter one, following the sequential nucleons evaporation.Actually, the CN+PE statistical emission mainly contributes to 90 Nb activation while the breakup and pick-up processes have also a significant addition of ∼15%.
The calculated results for 90 Nb activation above 28 MeV had to include a steady contribution due to the decay of 90 Tc (T 1/2 =50.7 s) and consequently 90 Mo (T 1/2 =5.56 h), leading finally to a calculated cumulative 90 Nb excitation function.This addition, involved similarly for the comparison with the TENDL-2021 evaluation, has been mostly useful to account for the data at deuteron energies above 20 MeV.The measured cross sections [15,16,106] for the activation of the residual nucleus 89 Zr [Fig.23(a)] and 89 Zr m state [Fig.23(d)] are cumulative, i.e. they include the decays of 89 Tc (T 1/2 =12.9 s) and 89 Mo (T 1/2 =2.11 min), and 89 Nb (T 1/2 =2.03 h) [90].Consequently, the present model analysis has taken into account, in addition to the 92,94 Mo isotopes contributions, these decays to 89 Zr residual nucleus of 89 Mo, including the 89 Tc decay, as well as the decay of 89 Nb.Thus, the experimental data are rather well described in Fig. 23(a,d) while it should be pointed out the strongest contributions of 89 Nb residualnucleus decay versus the deuteron-induced reactions on 92,94 Mo isotopes.The same addition has been considered for the comparison with the TENDL-2021 evaluations.
Concerning the reaction mechanisms contributing to 89 Zr residual-nucleus activation, the dominant PE+CN statistical emission is additionally increased at higher energies by the BF enhancement [Fig.23(b,c)].Particularly, this BF contribution is relevant to 92 Mo(d, 2n3p) 89 Zr reaction.Actually, the first maximum of the 89 Zr excitation function for both 92,nat Mo targets corresponds to α-particle emission within (d, pα) reaction (Q=1.491MeV), so important for the gas accumulation assessment.

nat Mo(d, xn3p) 88 Zr
The model analysis for the activation of 88 Zr residual nucleus involves fewer issues and thus may have a more reliable outcome.Thus, it concerns contributions due to deuterons interaction mainly with 92,94 Mo isotopes, as well as following 88 Nb decay (T 1/2 =7.78 min) [90] taken into account for the present analysis results as well as well as in the case of the TENDL-2021 evaluations.A further contribution of 88 Mo decay [90] to 88 Zr activation cross sections is not significant in the energy range of this work.
The modeling results describe satisfactory the measured cumulative 88 Zr excitation function [16] shown in Fig. 23(g).The largest contribution is related to the PE+CN mechanisms [Fig.23(h,i)], and also to the 92 Mo isotope.An yet visible addition is due to the inelastic breakup enhancement in the case of 92 Mo(d, 3n3p) reaction.The so weak contribution of 88 Nb decay to 88 Zr activation in Fig. 23(g), in addition to 92,94 Mo iso-topes components, has pointed out the statistical emission model validation on this experimental basis.

IV. CONCLUSIONS
A consistent analysis has concerned all reaction mechanisms involved in the complex deuteron interaction process, e.g.elastic scattering, breakup, stripping, pickup, pre-equilibrium and compound nucleus, including the contributions of specific Mo isotopes too.The brief overview of the activation data for the deuteron-induced reactions on nat Mo, the related TENDL-2021 evaluation [19], and model calculations carried out within present work is shown in Figs.24 and 25.It is thus underlined comparatively the neutron emission nat Mo(d, xn) and the neutron emission possibly accompanied by H and He isotopes emission nat Mo(d, xnyp), if emerging nucleon clusters are given due consideration.They are related to the high requests of the radioactivity risks and radiation damages estimation, critical benchmark for selecting and validating the best structural materials and a number of key technologies.Thus, accurate estimation of (d, xn) activation is essential for shielding design of nuclear installations while H, T, and He gas emission leading to "gas bubbles accumulation", are of high interest for damage studies of the structural materials.Overall, the due account of most experimental data has thoroughly validated the present theoretical framework and highlighted some prevalent features.  9Mo, 95 Nb g+m,g , 90 Nb, 88 Zr, and 89 Zr cumulative activation including the decay of 90 Tc, 95 Zr, 90 Mo, 88 Nb, 89 Nb (dashed) and 89 Mo (dash-dotted), respectively (see text) Hence strong neutron emission goes through (d, 2n) and (d, 3n) reaction channels (Figs.[10][11][12][13][14][15], while the charged-particle emission (protons, tritons, 3,4 He) is almost an order of magnitude less (Fig. 25).There is only one different case in this respect, of the strong protons emission within 100 Mo(d, p) 101 Mo reaction.Then, higher protons and even tritons emission leads to 99 Mo residual nucleus, with a maximum of ∼65 mb for the corresponding excitation function below 10 MeV [Fig.25(b)]. Inceased emission of neutrons, protons, 3 He and αparticles concern the activation of Nb and Zr isotopes, whose estimation is important for both theory validation and structural materials comparative analysis.

3 .
nat Mo(d, xn) 97 Tc m )].Nevertheless, the neglected additions of BF and especially DR processes might contribute to the apparent underestimation of nat Mo(d, xn) 97 Tc m excitation function by TENDL-2021 in Fig. 11(a).
)].It is yet close to the inelastic breakup enhancement for (d, xn) reactions, with x=2-4.The competition between the in-