Advanced breakup-nucleon enhancement of deuteron-induced reaction cross sections

Following the EUROfusion PPPT-programme action for an advanced modeling approach of deuteron-induced reaction cross sections, as well as specific data evaluations in addition of the TENDL files, an assessment of the details and corresponding outcome for the latter option of TALYS for the breakup model has been carried out. The breakup enhancement obtained in the meantime within computer code TALYS, by using the evaluated nucleon-induced reaction data of TENDL-2019, is particularly concerned. Discussion of the corresponding results, for deuteron-induced reactions on 58\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{58}$$\end{document}Ni, 96\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{96}$$\end{document}Zr, and 231\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{231}$$\end{document}Pa target nuclei up to 200 MeV incident energy, includes limitations still existing with reference to the direct-reaction account.


Introduction
The design analyses of the high intensity accelerator-based neutron source International Fusion Materials Irradiation Facility (IFMIF) [1], for ITER fusion device [2] and European DEMO fusion reactor [3] material qualification, require also high-quality cross-section data of deuteron induced reactions. Related activities are currently conducted within the Power Plant Physics and Technology (PPPT) programme of the EUROfusion Consortium [4], following a previous partnership agreement of Fusion for Energy (F4E) [5] with the Consortium on Nuclear Data Development and Analysis [6]. A sub-library of the latest version 3.2 [7] of Fusion Evaluated Nuclear Data Library (FENDL) [8], developed at IAEA/NDS, is going to include the most recent corresponda e-mail: marilena.avrigeanu@nipne.ro (corresponding author) b e-mail: dimitri-alexandre.rochman@psi.ch c e-mail: a.koning@iaea.org ing results, while there is a world-wide use of the deuteron sub-library of the current TENDL-2019 evaluated nuclear data library [9] based on the output of the TALYS nuclear model code system [10] for direct use in both basic research and applications.
The TENDL deuteron data library was in fact the current reference used in the PPPT project of the Early Neutron Source (ENS), the need of its further improvement being recently pointed out [11]. The dedicated effort already allocated in the PPPT programme including also contributions from all involved deuteron reaction mechanisms [12] has been concerned in this respect. It followed specific data evaluations which have complemented and updated the TENDL files which are based to a large extend on automated nuclear model calculations with the TALYS code using default models and parameters [13].
Actually, the mechanisms of deuteron-induced reactions are more complex than those of nucleon-induced ones due to the deuteron weakly bound nature. Thus, owing to its binding energy B d =2.224 MeV, the deuteron easily breaks up in the Coulomb and nuclear fields of a target nucleus. Then, this already intricate process involves also a variety of reactions initiated by the neutron and proton that follow the deuteron breakup (BU) [14,15], i.e. the elastic breakup (EB) in which the target nucleus remains in its ground state and none of the deuteron constituents interacts inelastically with it, and the inelastic breakup or breakup fusion (BF) [16], where one of these deuteron constituents interacts nonelastically with the target nucleus throughout breakup states in the continuum. Therefore, the usual low-energy reaction mechanisms of direct reaction (DR), pre-equilibrium emission (PE) and compound-nucleus (CN) de-excitation should be considered not only for the incident deuterons on a target nucleus but also for reactions induced by breakup-nucleons on the same nucleus.
On the other hand, these specific issues made possible projects of highly efficient fast-neutron sources by using deuterons accelerated on light nuclei (isotopes of Li, Be, C, etc), including IFMIF, NFS [17], and SARAF [18]. At the same time, deuteron activation of various structural materials (Al, Fe, Ni, etc.) within deuteron-accelerator components needs also a suitable assessment, despite of a scarce experimental data basis, too. Therefore, accurate and comprehensive deuteron nuclear data libraries over wide ranges of incident energies and target mass from light elements to structural materials, are mandatory for the engineering design of a deuteron-accelerator neutron source facility as IFMIF.
Nevertheless, TENDL-2017 deuteron sub-library improvement is based on the use in TALYS of the latest phenomenological BU model by Kalbach [19], following the analysis of an extensive database of double-differential cross sections for the breakup of deuterons, 3 He ions, and αparticles. Thus, global predictions for the absorptive breakup cross section have been involved, at once with Kalbach's pick-up and stripping models [20] for all combinations of projectiles and ejectiles. However, while it was decided to adopt TENDL-2017 for FENDL-3.2 charged-particle induced activation files, data from the incident chargedparticle medical isotope cross section data [21], if available, are taken in preference to TENDL-2017 [22]. This decision followed the recommendation of full parametrization of the available data, by a least-squares fit with Padé approximations of variable order, for deuteron monitor reactions and therapeutic radionuclides-production cross sections [23,24].
In fact, deficiencies in describing the elastic and especially the inelastic BU that dominates the deuteron reaction cross section at beam energies comparable with the Coulomb barrier motivated this recommendation [24], however with no predictive power. Earlier, a simple phenomenological description of the available stripping (d, p) experimental data was also proposed and included in the codes ALICE-D and EMPIRE-D [25]. It includes six parameters obtained by analysis of the excitation functions of four target nuclei, from 59 Co to 197 Au, and estimated accuracy of about 25-30%. However, its further use has not been reported yet while another energy-dependent enhancement factor K has most recently been mentioned without either definite values or model analysis [26,27].
Following the above-mentioned PPPT-programme action for an advanced modeling approach and specific data evaluations in addition of the TENDL files [11,12], an assessment of the details and corresponding outcome related to a latter option of TALYS for the breakup model makes the object of this work. An overview of the model approaches related to deuteron nuclear data for the last decade and a half, including the activities conducted recently in the PPPT programme as a first step of the flow chart linking modern theoretical nuclear physics to large-scale nuclear energy systems [28], is given in Sect. 2 of this paper. The consequent assumptions and add-ons proved useful in the meantime within TALYS, concerning particularly the breakup enhancement of the deuteron-induced reaction cross sections, are subject of Sect. 3. The corresponding results and their discussion including limitations still existing with reference to the direct-reaction account, are shown in Sect. 4, as well as conclusions in Sect. 5.

Deuterons on light nuclei
In order to improve the calculation of the d-Li neutron source term of IFMIF, an analysis of the elastic-scattering angular distributions of deuterons on 6,7 Li, for energies from 3 up to 50 MeV, using a phenomenological optical model potential (OMP) [29] was involved within an updated d+ 6,7 Li data evaluation [30,31]. The OMP analysis was completed by the coupled reaction channels (CRC) method within the code FRESCO [32] for consideration of the elastic transfer additional mechanism [33]. The reliability of the deduced particular and average OMP parameters was proved at variance with the extrapolation of the deuteron global potentials over a wider range of mass and energy domains. The key output of this analysis has been the deuteron total-reaction cross sections that are involved subsequently within all related data evaluation. A clear improvement of the prediction accuracy for the IFMIF neutron source term simulation was proved through benchmark calculations against experimental thick lithium target neutron yields [30,31].
Also in connection with the IFMIF project, Ye et al. analyzed in Ref. [34] the deuteron elastic scattering from 6,7 Li in the energy range from 10 to 50 MeV by the continuum discretized coupled channels (CDCC) method ( [35] and Refs. therein). They paid thus due attention to the deuteron breakup by a quantum-mechanical approach, taking into account that intermediate-energies semiclassical models were used previously [30] at a relatively low energy and neglecting the EB contribution. The deuteron breakup process is indeed involved within CDCC method by explicitly using a phenomenological three-body Hamiltonian in which the nucleon-nucleus interaction is represented by the OMP at half the deuteron incident energy and an effective nucleonnucleon potential is used for the p-n interaction. A good agreement was obtained at forward angles to the same extent as the optical model calculation with the phenomenological deuteron OMP [29], while a reduction of the elasticscattering angular distributions at large angles was found due to the BU account.
The analysis of the deuteron EB in terms of the CDCC method was then completed in Ref. [36] by that of the nonelastic breakup (with the notation STR for nucleonstripping in the continuum, at that time) using the semiclassical approach of the Glauber model [37] to account for angular distribution and double differential cross sections (DDX) of (d, xn) and (d, x p) reactions on 7 Li at deuteron energies of 40 and 100 MeV, respectively. This investigation of the Glauber model applicability at incident energies below 40 MeV/nucleon proved successful at forward angles, showing also that the nonelastic breakup component is more important than EB in deuteron reactions on 7 Li. Shapes and magnitudes of the experimental (d, x p) spectra at 100 MeV have been reproduced at forward angles (<20 • ) also for the target nuclei 9 Be and 12 C [38]. Results of a CDCC extension [39] for the DDX data of (d, xn) reaction on 7 Li at 40 MeV seems to be only slightly better than the Glauber-model ones.
Moreover, a modification of the Glauber model that was found necessary at low energies, to take into account the much stronger effects of the nuclear and Coulomb potentials than those at high energies, corresponding to a projectile curved trajectory [40], improved the total-reaction and EB cross sections for six target nuclei, from 7 Li to 208 Pb, at energies up to 100 MeV/nucleon. However, despite of remarkable results, the theoretical foundation of this semi-classical model was considered still questionable in the low energy region where the wavelength of the projectile becomes large and the classical trajectory approximation becomes worse.
In addition to the EB and nonelastic breakup modeled by CDCC method and Glauber model [39], respectively, the account of the stripping reactions to bound states of the residual nuclei through zero-range distorted wave Born approximation (DWBA) was considered within angular-distribution analysis for the (d, p) reaction on 12 C as well as heavier nuclei 27 Al, 40 Ca, and 58 Ni, at deuteron energies from 7 to 56 MeV [41]. This BU + DR analysis was completed by PE + CN contributions provided by the exciton and Hauser-Feshbach statistical models. DDX data integrated over the outgoing energy, for incident energies from 56 to 100 MeV, were similarly analyzed as well as angular-distributions for the 12 C(d, n) 13 N reaction at incident energies from 7 to 18 MeV [42]. Of particular interest has been a similar, even if less clear, energy dependence of the spectroscopic factors extracted from DWBA analyses for the (d, p) reaction on all the four target nuclei, under the assumption that the incidentenergy dependence is same for all final states. An empirical expression of this energy dependence has been obtained, with various scaling factors for each target nucleus and residual state, while finally was considered that it might not be possible to deduce such a simple empirical formula.
The above-mentioned CDCC-Glauber-DWBA analysis completed by PE + CN contributions within deuteroninduced reaction analysis code system DEURACS, for DDX data of (d, xn) reactions on 9 Be and 12 C, at incident energies up to 50 MeV, took into account in Ref. [43] also the PE + CN decay following the absorption by the target nucleus of either a BU-neutron or a BU-proton in addition to a deuteron itself. An incoherent summation of the EB and nonelastic-breakup components was thus explicitly enlarged by three statistical-decay components, corresponding to the three distinct compound nuclei and weighted by nonelasticbreakup formation fractions R d , R n , and R p calculated with the Glauber model [39]. The additional (n, xn) and ( p, xn) components were calculated for nucleon-incident energies at half the deuteron incident energy for saving computational time, while the DDX data were well reproduced quantitatively over a wide range of neutron emission energy except the energies below a few MeV.
Recently, energy-dependent formation factors R n (E n ) and R p (E p ) have been involved in Ref. [44] within the same DEURACS framework, as given by Glauber model [39] for the nonelastic-breakup (NBU, replacing the authors' former notation STR). The related analysis of DDX data for (d, xn) reactions on 9 Be and 12 C, at incident energy of 102 MeV, has proved the applicability of DEURACS to the (d, xn) reactions on light nuclei at incident energies around 100 MeV, as well as the essential NBU role of for accurate predictions. On the other hand, it has been pointed out that the related NBU model may not be valid for low incident energy close to the Coulomb barrier, so that it needs to be additionally tested for such data. Nevertheless, a deuteron nuclear data library up to 200 MeV for 7 Li, 9 Be, and 12,13 C has been developed on this basis, with particular attention to neutron production in Ref. [45]. The above-mentioned underestimation below a few MeV, which was attributed to an absent particle-decay from discrete levels in residual light nuclei [43,44], has also been amended.

Breakup-parametrization based model analysis
Former analysis of deuteron-induced activation data for IFMIF calculations showed that the deuterons are much more important than the neutrons (about a factor of 70) in activating elements other than lithium, so that refinement of the model calculations and further experimental measurements are needed if the deuteron libraries are to approach the standard of the established neutron ones [46], for a suitable input into IFMIF engineering validation and engineering design activities (EVEDA) phase.
Along this line, an increased attention was firstly provided to BU effects as well as the simultaneous analysis of the deuteron elastic-scattering and induced activation cross sections as for 27 Al in Ref. [47]. Thus, a parametrization of both EB cross section σ E B and total BU proton-emission cross section σ p BU =σ E B +σ p B F was obtained by analysis of experimental systematics of deuteron-induced reactions on target nuclei from 27 Al to 232 Th [14,16,[48][49][50]. Equal BF nucleon-emission cross sections σ n B F and σ p B F were assumed [16], so that the total BU cross section is σ BU = σ E B +2σ n/ p B F . It should be underlined that the experimental data [14,16,[48][49][50] were obtained through a distinct analysis of the BU, DR, PE, and CN reaction mechanisms.
Consequently σ BU was subtracted from deuteron totalreaction cross section σ R that should be shared among different DR + PE + CN outgoing channels, while the BF component brings additional contributions to different (d, x) reaction channels through (n, x) and ( p, x) reactions on the same target nucleus 27 Al. This breakup enhancement of various (d, x) reaction cross sections by breakup-neutrons through (n, x) reactions was taken as the product between σ n B F and the ratio σ (n,x) /σ non , where σ non is the neutron non-elastic cross section. The similar effect of the ( p, x) reactions was given by the equivalent σ p B F multiplied by the ratio σ ( p,x) /σ R , where σ R is the proton total-reaction cross section. The two ratios were expressed as a function of the deuteron incident energy using the Kalbach [51] formula for the center-of-mass system centroid energy of the deuteron-breakup peak energies of the emitted constituents.
At the same time there have been used: (i) a semimicroscopic OMP based on a double-folding model (DFM) real part, and phenomenological imaginary and spin-orbit terms together with dispersive corrections to the DFM results, and (ii) a related phenomenological OP real part replacing the DFM one, both of them being proved by analysis of deuteron elastic-scattering data at energies from 5 to 59 MeV, and σ R at ∼13 and 25 MeV, as well as (iii) the DR contributions obtained with code FRESCO [32] and DWBA spectroscopic factors from data analysis of (d, p) angular distributions for 35 final states of the residual nucleus 28 Al up to 5.135 MeV, at incident energies of 6, 12, and 23 MeV, and also (d, n) angular distributions for 18 final states of 28 Si up to 10.38 MeV, at incident energy of 6 MeV [47], and (iv) usual PE + CN contributions of the code STAPRE-H [52]. This simultaneous analysis of the deuteron elastic-scattering and induced activation cross sections including new data measured at incident energies from 3 up to 20 MeV was finally providing suitable account of all available excitation functions, up to the incident energy of 60 MeV, not only for the (d, p) but (d, α), (d, 2n), (d, 2 p), and (d, pα) reactions, too [53].
A further point concerned an apparent decrease of the fraction f E B = σ E B /σ R at energies beyond the range E <30 MeV of the EB data [14,16], unlike the trend of both the fraction f p BU = σ p BU /σ R and σ R . Consequently, the correctness of this EB extended parametrization has been checked in Ref. [54] through the comparison of its predictions and results of the CDCC method for the 63 Cu and 93 Nb tar-get nuclei, at incident energies up to 50 MeV. Following a CDCC suitable account of experimental elastic-scattering angular distributions of deuterons on 63 Cu and 93 Nb at energies between 12 and 52 MeV, a good agreement was found between the EB data [14] and the corresponding CDCC results as well as empirical parametrization [47] at the energies of the available data. However, f E B becomes significantly lower at higher energies. Therefore, extrapolation of f E B empirical parametrization beyond the energies of the corresponding data have to be considered with caution, at the same time with the challenging CDCC calculations for each target/energy of interest. In the meantime a normalization of the EB fraction for the energies beyond the maximum of the former parametrization [47] has been adopted in Ref. [55] by taking unchanged the ratio of the EB and BU fractions at the incident energies above the f E B maximum.
However, the former breakup-enhancement broad approximation [47] using only the centroid energies of the breakup nucleons [51], was at variance with quite large widths of the assumed Gaussian distribution of the corresponding deuteron-breakup nucleon energies, for a given deuteron incident energy. Therefore it has been replaced by a convolution of the each one of energy-dependent ratio σ (n,x) /σ non or σ ( p,x) /σ R , and the related Gaussian line shape. The corresponding distributions and convolution results, at various deuteron energies, shown for deuterons incident on Cu [56], 231 Pa [57], and 93 Nb [58] proved the need for this improvement, with the detailed formula given elsewhere [59,60]. Additionally, the key role of BU and DR processes in deuteron-induced reactions has been pointed out in Ref. [12] by the differences between larger CN contributions obtained from measured neutron angular distributions, with no BU consideration, and BU + DR + PE + CN analyses.
Actually, the same approach beyond elastic-scattering analysis for OMP validation, using BU parametrization [47,55] for the related deuteron σ R decrease as well as breakup-enhancement of various (d, x) reaction cross sections, DWBA spectroscopic factors from data analysis with FRESCO of available particle-emission angular distributions, and PE + CN framework of either STAPRE-H or then TALYS codes has been consistently involved within suitable account of all available excitation-function data of deuterons up to 60 MeV on 27 [65]. Due consideration of all EB, BF, DR, PE, and CN reaction mechanisms [60,61] is indeed found crucial in Ref. [66] for a consistent analysis of the deuteronreaction measured data and even high production of protonrich nuclei [67], while insufficient treatment and separation between different reaction mechanisms such as DR and BU components [68] may be related to deviations between measurements and advanced surrogate reaction studies [69]. On the other hand, the concerned choice for the BU fractions has thus been validated.

Breakup-semi-classical account
On the other hand, at once with the work on deuteron-induced reactions on Li, Be, and C isotopes, the above-mentioned CDCC results [39] reproduced also experimental elasticscattering angular distribution of deuterons on 58 Ni at 200 MeV/nucleon, while the Glauber model was found good for light targets, as expected, but not for EB and nonelastic breakup for heavy targets. At the same time, the particular attention paid to the deuteron BU over a wide mass range of target nuclei from 9 Be to 238 U [38] has proved shapes and magnitudes of the experimental (d, x p) spectra at 100 MeV which are reproduced at forward angles (<20 • ) for target nuclei less than 58 Ni but underestimated the specific bump around half the incident energy as the target mass increased. It was thus concluded that more quantummechanical approaches beyond the Glauber model would be necessary to account for (d, x p) reactions for heavier-mass nuclei at relatively low incident energies below 100 MeV.
Similar results have been obtained for elastic-scattering angular distribution of deuterons on 27 Al and 58 Ni at incident energies between ∼ 11-79 MeV, as well as for corresponding DDX data of (d, x p) reactions at 56 and 100 MeV [70]. An additional note may concern in this case the inclusion of PE + CN contributions within the exciton and Hauser-Feshbach models, corresponding to the absorption by the target nucleus of either a BU-neutron or a BU-proton in addition to a deuteron itself. The related NBU formation fractions R n , R p , and R d were calculated with the Glauber model [39,40] and the assumption that the incident energies of proton and neutron are half the deuteron incident energy. Actually, the additional two statistical-decay components correspond to the former breakup-enhancement approach [53] that involved BF cross sections [47] and BU-nucleon centroid energies of Kalbach [51].
The above-mentioned CDCC-Glauber-DWBA analysis, including the PE + CN components as above, of the (d, p) reaction also on 27 Al, 40 Ca, and 58 Ni [41] has reproduced also the experimental angular-distributions of the (d, p) reaction on 27 Al at 6, 12, and 23 MeV [71] as well as 27 Al(d, p) 28 Al excitation function from the threshold to 20 MeV. More recent DEURACS analysis of (d, x p) reactions on 27 Al at incident energies of 56 and 100 MeV [72] has reproduced fairly well the DDX data including both the NBU broad peak at half the incident energy, and the peak structure due to the stripping DR to bound states. The same approach has been proved successful for the account of 45 Sc(d, p) 46 Sc excitation function for incident energies up to 50 MeV.
At the same time, DDX data of (d, xd) reactions on 27 Al and 58 Ni at incident energies of 80 and 100 MeV, respectively [44], have been better reproduced than the original PE + CN models due particularly to inclusion of the Kalbach's semi-empirical inelastic scattering model [73]. Energy-dependent formation factors R n (E n ) and R p (E p ), given by the Glauber model [39,40] for the NBU processes, have been also involved in a way rather similar to the above-mentioned BU enhancement based on the convolution of breakup-nucleons energy distributions and the related nucleon-induced reaction excitation functions [56][57][58][59][60].
A conclusive assessment is that, 'taking into account that the CDCC method and the Glauber model have been widely used in analyses at incident energies above 200 MeV/nucleon, it is expected that the framework of DEURACS for the breakup processes is applicable to the incident energy range of 100-200 MeV/nucleon' [74]. This approach has thus been considered in Refs. [74,75] as being complementary to the method using the deuteron-breakup empirical formulas based on data up to 80 MeV [47] and proved successful for reactions at relatively low incident energies below 60 MeV [53,[56][57][58]60,62,63].

Empirical parametrization
Present analysis of the deuteron breakup mechanism is based on parametrization [12,47,54,55] of both the total breakup (EB + BF) and EB data, under the assumption that the inelastic-breakup cross section for neutron emission σ n B F is the same as that for the proton emission σ p B F [16]. Thus, the total breakup cross sections σ BU is given by the sum σ E B + 2σ n/ p B F . The parametrization has concerned the total BU nucleon-emission and EB fractions, i.e. f n/ p BU = σ n/ p BU /σ R and f E B =σ E B /σ R , respectively, where σ R is the deuteron totalreaction cross section. The dependence of these fractions on the deuteron incident energy E and the target-nucleus atomic Z and mass A numbers was obtained [47,55] through analysis of the experimental systematics of deuteron-induced reactions on target nuclei from 27 Al to 232 Th and incident energies up to 80 MeV for the former [14,[48][49][50], but within a more restricted energy range up to 30 MeV [14,16,[48][49][50] for the latter: The comparison of experimental data and parametrization results shown in Fig. 1 for deuterons incident on nuclei from 27 Al to 232 Th, at energies up to 80 MeV, has proved a suitable  62 Ni (e) 20  agreement [55]. On the other hand, it was found an apparent decrease of the fraction f E B at energies beyond the range E<30 MeV of the EB data [14,16], unlike the trend of both the fraction f p BU and the total-reaction cross section. Because of that, the correctness of the extrapolation of EB parametrization has been checked [54] by comparison with results of the microscopic CDCC method [35]. Thus, a normalization factor has been introduced [55] for extrapolation of f E B at energies beyond the available data, in agreement with the behavior of f p BU [47] and the CDCC calculation results [54]. Hence, we have chosen to keep unchanged the ratio of the EB and BU fractions at incident energies above an energy E max of the f E B maximum [55], by means of the relation: Therefore, the normalized EB fraction (dotted curves) follows the behavior of the total BU nucleon-emission fraction shown in Fig. 1 to be in agreement with the CDCC results [54]. Despite the EB component is less than 10% of total BU cross section, this f E B normalization is of particular interest at deuteron energies above ∼ 50 MeV and especially for heavier target nuclei, for the inelastic breakup fraction f n/ p as well as the total breakup fraction under the above-mentioned assumption of equal neutron-and proton-emission BU cross sections [16,47,55]. On the other hand, since the corresponding data errors amount to 10-15% [47], a similar accuracy has been considered also for the present parametrization.

Additional BU constraint for A > 200
The total BU proton-emission data systematics for heaviest nuclei (A>200) at incident energies around the Coulomb barrier, of great interest for deuteron interaction with actinides nuclei , includes only one single point for 232 Th at E = 15 MeV [14]. The parametrization [55] describes it properly as well as the value at the incident energy of 70 MeV reported by Wu et al. [49] (Fig. 1). However, following the EB fraction normalization by Eqs. (3) and (5), the total BU fraction corresponding to this target nucleus may exceed unity. To avoid this unphysical overrun of the total-reaction cross section, and taking into account the above-mentioned systematics accuracy, an additional constraint was adopted for A > 200, namely 0.9 for the f BU maximum value. Nevertheless, the empirical parametrization should be confirmed by further data measurements and also advanced theoretical modeling.

Additive empirical parametrization and microscopic studies
The recent detailed analyzes of deuteron both EB and BF components by the DWBA method [76,77] performed a successful description of proton spectra and angular distributions for the (d, p) reaction on 27 Al, 58 Ni, 93 Nb, and 118 Sn at incident energies from 15 MeV to 100 MeV. The corresponding calculated EB and BF cross sections would be, however, also quite useful for the comparison with experimental data [14,[48][49][50]. More recently, Neoh et al. [76] applied a CDCC extension, using also microscopic optical potentials, to the analysis of the EB and neutron removal cross sections at 28 MeV/nucleon on various target nuclei from 12 C to 209 Bi, of interest for actinides nuclei [57,[78][79][80] and further studies of unstable-nuclei structure. Consequently, the measured total BU proton-emission cross sections σ p BU at 15, 25.5, 56, 70 and 80 MeV deuteron energies and for target nuclei from 12 C to 232 Th [14,[48][49][50] are compared in Fig. 2 with the above-described parametrization and results of CDCC extension [76] and DWBA method [77]. Because the absolute cross sections may depend Thus it results that increase of σ p BU with the target-nucleus mass is well described by the empirical parametrization for all deuteron energies from 15 to 80 MeV. A similar trend of the microscopic results for medium-mass nuclei with 40 <A<120, but an overestimation of the measured data are apparent for light nuclei (A<40) at both 25 and 56 MeV incident energies, as well as an underestimation for heavier ones (A>120). Moreover, the importance of the BU mechanism, shown by f p BU , is increasing with the target-nucleus mass, from 27 Al up to 232 Th, at the lower incident energies of 15 and 25.5 MeV. This increase is less significant at the energy of 56 MeV, and even reversed at 70-80 MeV so that it seems that the fraction f p BU has reached its maximum at 56 MeV for the target nuclei with A>120. Then, this maximum moves at energies over 56 MeV for 40<A<120, and is still increasing with the incident energy even at 80 MeV for light target nuclei (A<40). The energy dependence of the measured f p BU are satisfactorily described by the empirical parametrization, while the microscopic results are almost constant for all nuclei at 25 MeV [77], show a steep decrease from A = 12 up to A ∼ 120 at 56 MeV [76,77], apart from the data, and underestimate but describe however the A-dependence for A>120.
Nevertheless, the recently increased interest on the theoretical analysis of the breakup components, e.g. Refs. [81][82][83][84][85], may lead eventually to improvement of technologically relevant deuteron-breakup empirical parametrization, already involved within successful analysis of all available data for various nuclei [18,53,56,58,[60][61][62][63][64][65], as well as additional validation and progress of microscopic calculations. It is noticeable in this respect the challenging task found by Potel et al. [84] for the future of the latter topic to be discrimination of the different non-elastic breakup contributions at least in energy regions where they are competing.

Deuteron inelastic-breakup enhancement
While the former effect of BU consideration should be the decrease of σ R that is shared among different outgoing channels, by the σ BU value, the latter one concerns the BF contributions to different reaction channels. Thus, the absorbed proton or neutron following the deuteron breakup contributes to enhancement of the corresponding (d, xn) or (d, x p) reaction cross sections, respectively. The compound nuclei in reactions induced by BF nucleons differ by one unit of the atomic mass and maybe of also the atomic number than that one in deuteron-induced reactions on the same target nucleus. Moreover, the partition of the BF cross section among various residual-nuclei population is triggered by the energy spectra of the BF nucleons and the excitation functions of the reactions induced by these nucleons on the same target nucleus.
Hence, in order to calculate the BF enhancement of, e.g., the (d, xn) reaction cross sections, the BF proton-emission cross section σ p B F should be (i) multiplied by the protonenhancing ratios σ ( p,x) /σ p R , (ii) convoluted with the Gaussian line shape distribution of the BF-proton energy E p for a given deuteron incident energy E [56][57][58], and (iii) integrated over the BF-proton energy. Consequently, the BFenhancement cross section has the form [59][60][61]64]: where σ p R is the proton total-reaction cross section, x stands for various γ , n, d, or α outgoing channels, while the Gaussian distribution parameters w and E 0 p given by Kalbach [51] were used.
Similarly, the BF enhancement of the (d, x p) reaction cross sections brought by breakup neutrons is given by the following expression: where σ n T is the neutron total-reaction cross section, and σ (n,x) /σ n T are the neutron-enhancing ratios. Interpolation of experimental nucleon-induced reaction cross sections from the EXFOR library [86] or from newest TENDL library has been involved within estimation of the BU enhancement [53,[56][57][58][60][61][62][63][64][65]80,87], in order to reduce as much as possible the supplementary uncertainties brought by additional theoretical calculations. In the present work the evaluated neutrons and protons cross sections have been taken from the TENDL-2019 library [9].
On the whole, the enhancing effect of the breakup mechanism is important mainly for describing the excitation functions of second and third chance emitted-particle channels [53,[56][57][58][60][61][62][63][64][65]80,87]. This is shown in Fig. 3 for the BF enhancements brought by the breakup-nucleon interactions with various target nuclei from 55 Mn up to 231 Pa. Particular comments may concern the BF contribution that is even larger than PE + CN ones for the (d, 2 p) reaction at all incident energies (Fig. 3d) as well as for (d, 2n) reaction above 30 MeV (Fig. 3b) and (d, 3n) reaction just above the effective threshold (Fig. 3c). Moreover, the entire excitation functions of also (d, nα), (d, pα), and (d, t) would be not described without BF consideration despite of the dominant nature of the corresponding PE + CN components shown in Fig. 3e-h.

Discussion
Selection of the breakup formalism in TALYS code is available by means of the breakupmodel keyword, with value 1 for Kalbach's [19] and 2 for Avrigeanu et al. models. The latter uses Avrigeanu et al. parametrization [12,47,54,55] to calculate the total nucleon-emission, elastic, and inelastic breakup fractions by Eqs. 1, 2, 4, respectively, the corresponding cross sections σ as well as the BU neutron-and proton-emission spectra on the basis of the former Gaussian parametrization of Kalbach [51]. A distinct note should concern the maximum energy of the BU nucleons spectra which is the difference between the deuteron incident energy and its binding energy.
Verification and validation (V&V) [7] of the advanced breakup account in addition to TALYS-1.95 has been carried out for deuteron-induced reactions on 58 Ni, 96 Zr, and 231 Pa target nuclei, up to 200 MeV incident energy. Of particular interest for the present work is the comparative analysis of the results of TALYS-1.95 using both the option breakupmodel 1, and option breakupmodel 2 with the additional breakup enhancement detailed within present work.
First, an outline of the additional breakup enhancement is made possible by comparison of the TALYS-1.95 default output and results of the implicit BU-enhancement addition in the cases of activation excitation functions of residual nuclei 64 Cu and 96 Nb by (d, xn) reactions on natural Ni and Zr, respectively (Fig. 4a,c). The latter values are increased by 57% and 86%, respectively, with respect to the former ones at a reference deuteron energy of 40 MeV. Nevertheless, results identical with the TALYS-1.95 default output are obtained by using the option breakupmodel 1 also within latest code version, which includes the BU enhancement only for the breakupmodel 2 option.
The BF explicit contribution by means of the ( p, n) reactions induced by the BU-protons on the same target nuclei may be shown for distinct isotopes, e.g. the heaviest stable ones, in Fig. 4b,d. Actually they were selected due to the corresponding larger (x, 2n) reaction cross sections. These BF contributions amount to ∼ 36% and 46%, respectively, of the total activation cross sections of the two residual nuclei at the same reference deuteron energy of 40 MeV.
The comparative analysis of these functions reveals both complementary effects of the breakup mechanism: overall decreasing of the excitation functions due the leakage of initial deuteron flux to the breakup process, and the inelastic breakup enhancement. Due to the restrictive energy range of the experimental BU fractions [14,[48][49][50] involved for settling our parametrization [55], the corresponding excitation functions for BF enhancement shown in Fig. 5 are not going over 80 MeV. Nevertheless, the 80 MeV deuteron incidentenergy limit is still appropriate for the EUROfusion objectives, while an eventual switch to the default TALYS-1.95 results seems practical at higher energies.
On the other hand, the default TALYS-1.95 results and the statistical PE + CN contributions corresponding to latter breakup-model option, including the additional breakup enhancement, are rather close (dash-dotted and dash-dotdotted curves, respectively, in Fig. 5). The effects of breakup enhancement, that are also is evidenced by the difference between the total excitation functions and the PE + CN contributions, are larger for the (d, 2n) reactions and lighter target nuclei. However, the breakup enhancement becomes significant also for the heavier nuclei at deuteron energies above ∼ 40 MeV.

Summary and conclusions
Subsequently to the EUROfusion PPPT-programme action for an advanced modeling approach of deuteron reaction cross sections [11], as well as specific data evaluations [53,[56][57][58][60][61][62][63][64][65] in addition of the TENDL files, an assessment of the details and corresponding outcome related to a latter option of TALYS for the breakup model [12,47,55] has been carried out. It is based on an empirical parametrization of both the total breakup (EB + BF) and EB data, under the assumption that the inelastic-breakup cross section for neutron emission σ n B F is the same as that for the proton emission σ p B F [16]. Actually, the parametrization has concerned the total BU nucleon-emission as well as EB fractions, i.e. f n/ p BU = σ n/ p BU /σ R and f E B = σ E B /σ R , respectively. The dependence of these fractions on the deuteron incident energy E and the target-nucleus atomic Z and mass A numbers was obtained [47,55] through analysis of the experimental systematics of deuteron-induced reactions on target nuclei from 27 Al to 232 Th and incident energies up to 80 MeV for the former [14,[48][49][50]. Moreover, a normalization factor has been introduced [55] for extrapolation of f E B at energies beyond the available data, in agreement with the behavior of f p BU [47] and the CDCC calculation results [54]. At the same time, an additional constraint was adopted for A > 200, namely 0.9 for the f BU maximum value.
On the other hand, while the σ R value is decreased by the total-breakup cross section, the nucleons following the deuteron breakup and absorbed by the same target nucleus contributes to the enhancement of the (d, x) reaction cross sections through corresponding (n, x) and ( p, x) reactions. The partition of the BF cross section among various residualnuclei population is triggered by the energy spectra of the BF nucleons and the excitation functions of the reactions induced by these nucleons on the same target nucleus. Consequently, in order to calculate the BF enhancement of, e.g., the (d, xn) reaction cross sections, the BF proton-emission cross section σ p B F should be (i) multiplied by the proton-enhancing ratios σ ( p,x) /σ p R , with the evaluated proton-induced reaction data of TENDL-2019 concerned, (ii) convoluted with the Gaussian line shape distribution [51] of the BF-proton energy E p for a given deuteron incident energy E [56][57][58], and (iii) integrated over the BF-proton energy. Actually, the former breakupenhancement account [53] corresponds to the 2nd and 3rd statistical-decay components [70] within former DEURACS framework, weighted by nonelastic-breakup formation fractions R n and R p calculated with the Glauber model at half of the incident-deuteron energy. At the same time, the final one [59,60] corresponds to the DEURACS use of energydependent fractions [44] R n (E n ) and R p (E p ) also within Glauber model.
On the whole, the enhancing effect of the breakup mechanism is important mainly for describing the excitation functions for second and third chance emitted-particle channels [53,[56][57][58][60][61][62][63][64][65]80,87]. Particular comments concerned the BF contribution even larger than the PE + CN ones for the (d, 2 p) reaction for heavier nuclei at all incident energies, as well as (d, 2n) reaction above 30 MeV, and even (d, 3n) reaction just above the effective threshold (Fig. 5). Moreover, the entire excitation functions of also (d, nα), (d, pα), and (d, t) would be not described without BF consideration despite of the dominant nature of the corresponding PE + CN components (Fig. 3).
Verification and validation of the advanced breakup account in addition to TALYS-1.95 has been carried out for deuteron-induced reactions on 58 Ni, 96 Zr, and 231 Pa target nuclei. Thus it has been found that the BU-enhancement addition increases by 57% and 86% the activation cross sections of residual nuclei 64 Cu and 96 Nb in (d, xn) reactions on natural Ni and Zr, respectively, at a reference deuteron energy of 40 MeV. The BF contribution through ( p, n) reaction induced by the BU-protons on the heaviest stable isotopes 64 Ni and 96 Zr, corresponding to larger (d, 2n) reaction cross sections, amounts to ∼ 36% and 46%, respectively, of their total activation cross sections at the same reference deuteron energy of 40 MeV. Moreover, the effects of breakup enhancement are larger for the (d, 2n) reactions and lighter target nuclei. However, they become significant also for the heavier nuclei at deuteron energies above ∼ 40 MeV.
Actually, the same approach beyond elastic-scattering analysis for OMP validation, using BU parametrization [47,55] for the related deuteron σ R decrease as well as breakup-enhancement of various (d, x) reaction cross sections, DWBA spectroscopic factors and related uncertainties [65] from data analysis with FRESCO of available particleemission angular distributions, and PE + CN framework of either STAPRE-H or then TALYS codes has been consistently involved within suitable account of all available excitation functions of deuterons up to 60 MeV on 27 [65]. Due consideration of all EB, BF, DR, PE, and CN reaction mechanisms [60,61] is indeed found crucial [66] for a consistent analysis of the deuteron-reaction measured data and even high production of proton-rich nuclei [67], while insufficient treatment and separation between different reaction mechanisms such as DR and BU components [68] may be related to deviations between measurements and advanced surrogate reaction studies [69].
Nevertheless, the updated approach detailed in this work can be used to provide more accurate deuteron-reaction data that are also eventually required during the deuteron crosssection V&V activities within EUROfusion/ENS, with addi-tion of only the particular (d, p) and (d, t) reaction contribution to be obtained by using the direct-reaction DWBA method as also it is performed within the code system DEU-RACS for deuterons incident on light as well as increasingmass nuclei.
Acknowledgements This work has been partly supported by Autoritatea Nationala pentru Cercetare Stiintifica (Project PN-19060102) and carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014-2018 and 2019-2020 under Grant agreement No 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a model analysis of the published data, and results of all calculations are displayed in the figures and can be reproduced by using the parameter values given or properly referenced. Additional data related to this work are available from the authors on request.] 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/.