Effective extraction of photoneutron cross-section distribution using gamma activation and reaction yield ratio method

Zhi-Cai Li · Yue Yang · Zong-Wei Cao · Xin-Xiang Li · Yun Yuan · Zong-Qing Zhao · Gong-Tao Fan ·Hong-Wei Wang · Wen Luo

Abstract Photoneutron cross-section (PNCS) data are important in various current and emerging applications.Although a few sophisticated methods have been developed, there is still an urgent need to study the PNCS data.In this study, we propose the extraction of PNCS distributions using a combination of gamma activation and reaction yield ratio methods.To verify the validity of the proposed extraction method, experiments for generating 62,64Cu and 85m,87mSr isotopes via laser-induced photoneutron reactions were performed, and the reaction yields of these isotopes were obtained.Using the proposed extraction method, the PNCS distributions of 63Cu and 86Sr isotopes (leading to 85mSr isotope production) were successfully extracted.These extracted PNCS distributions were benchmarked against available PNCS data or TALYS calculations, demonstrating the validity of the proposed extraction method.Potential applications for predicting the PNCS distributions of the 30 isotopes are further introduced.We conclude that the proposed extraction method is an effective complement to the available sophisticated methods for measuring and evaluating PNCS data.

Keywords Photonuclear data · Photoneutron cross section · Gamma activation · Reaction yield · TALYS calculation

1 Introduction

Photonuclear data describe the interactions between photons and atomic nuclei.They are primarily composed of photoneutron cross-section (PNCS) data and photochargedparticle cross-section data.They play a vital role in a wide range of current and emerging applications [1].These applications include radiation shielding design, radiation transport analysis [2], absorbed dose calculations in the human body during radiotherapy, fission and fusion reactor physics and technology [3, 4], activation analyses, safeguards and inspection technologies, medical isotope production [5,6], nuclear waste transmutation, and astrophysical nucleosynthesis [7, 8].In the field of medical isotope production,the PNCS data can guide the production of radioisotopes for diagnostic and therapeutic purposes [9, 10].With the advent of new facilities that produce brilliant and intense photon beams, the photoproduction of neutron-deficient isotopes could become a competitive alternative to traditional methods that use neutrons generated by nuclear reactors or charged particle beams [11—13].

PNCS data have been experimentally obtained from various types of measurements, with bremsstrahlung and quasimonoenergetic photons from positron annihilation in flight,and more recently, from laser Compton scattering (LCS)[14—16].Most existing experimental data correlate with studies on PNCSs.Although PNCSs are obtained by counting the total number of emitted neutrons, the determination of partial PNCSs requires neutron multiplicity sorting [17].One advantage of bremsstrahlung measurements [18] is the high photon beam intensity.However, the technique has several disadvantages.First, it is necessary to know the bremsstrahlung spectrum sufficient for all electron energies.Next, measuring the reaction yield curve in a small energy step requires a stable accelerator and large counting statistics.Third, the process of subtracting the yield curves in the unfolding procedure may introduce correlations between the experimental data points,which can lead to fluctuations in the unfolded cross sections.In contrast to bremsstrahlung measurements employing a mathematical approach, positron annihilation in flight and LCS offered an apparatus method for producing them with variable energies [19, 20].However, the effective spectrum of the incident photon beam should be known in advance.Consequently,all these measurements require a precise scan of the incident photon energy and accurate knowledge of the incident photon spectrum and beam intensity.

Regarding the evaluation of the PNCS, the International Atomic Energy Agency has organized two coordinated research projects and photonuclear data of 219 isotopes were evaluated successfully because of the growing need for photonuclear data [21, 22].Despite the available experimental and evaluation efforts on photonuclear data, there is still an urgent need to study PNCS data for the following three reasons: (a)there is a lack of data for a number of cases; (b) it is not possible to produce a complete PNCS data file based on measured PNCSs alone; and (c) the experimental PNCS data often suffer from systematic discrepancies that are not easy to resolve [1].

In this study, we propose an effective method for the extraction of PNCS distributions using gamma activation and the reaction yield ratio method.Note that gamma activation can be accomplished through either laser plasmabased bremsstrahlung radiation or LCS photon irradiation of a solid target.This extraction method was verified by reasonably reproducing the experimental/evaluated and calculated PNCS data and can potentially be used to predict the PNCS distributions of 30 stable isotopes.The remainder of this study is organized as follows: in Sect.2, the extraction method is introduced to obtain the peak cross sections and further predict PNCS distributions.The benchmark results obtained using the extraction method are presented in Sect.3.The merits of the extraction method and the possibility of predicting PNCS distributions for 30 isotopes are discussed in Sect.4.The conclusions and outlook are presented in Sect.5.

2 Extraction of PNCS distribution

Relativistic laser-plasma interactions are bright, incoherent sources of gamma radiation with energies up to several tens of megaelectronvolts.They can be produced by compact tabletop laser systems that deliver laser pulses with energies of less than 100 fs.The spectra of the laser-accelerated electrons and bremsstrahlung photons can be described by the Boltzmann distribution [23]

whereTe,γis the electron or photon temperature,kBis the Boltzmann’s constant, andne,γ(0) is the number of electrons or photons within the energy bin atEe,γ=0 MeV.It can be seen that the intensity of laser-accelerated electrons or bremsstrahlung photons decreases exponentially with energy.The threshold energiesEth are generally a few MeV for photoneutron reactions.To efficiently trigger such reactions, the photon temperatureTγmust reach the MeV/kB.

For nuclides with a single giant dipole resonance (GDR)peak, the PNCS distribution can be described by a Lorentzian-like shape [24].

where the peak energyEm, peak cross sectionσmand widthΓare the three Lorentzian distribution parameters.These parameters can be acquired from TALYS calculations,TENDL evaluations, and experimental data available in the EXFOR database.Figure 1 shows that the Lorentzian curves accurately reproduce the experimental PNCSs for the four representative nuclides:63Cu ,88Sr ,127I, and197Au.The TALYS calculation underestimated the measuredσmvalues for63Cu and127I, whereas the TENDL evaluations aligned with the experimental data.In addition, visible discrepancies were found between the experimental data and the TENDL evaluations.

The reaction yieldYdepends on the temperatureTγand cross sectionσ(E).The expression for the reaction yieldYcan be written as:

wherenis the number of target nuclei per unit area, anddis the thickness of the target.For a specific target comprising isotopesAandBirradiated by a laser-generated bremsstrahlung photon beam, the residual nucleiA′andB′are produced through photoneutron reactions.When isotopesAandBhave very similar mass numbers, their PNCS distributions could be very close in terms of peak energyEmand widthΓbut vary in terms of the peak cross sectionσm.Figure 2 shows the PNCS distributions of Cu, Sr, I, and Au nuclei, whose isotopes have three adjacent mass numbers.The three Lorentzian parameters of their PNCS distributions are listed in Table 1.It can be seen that, although the peak cross sectionsσmfor these nuclides with adjacent mass numbers have large relative differences of 10—30%, the peak energiesEmare very close to each other because their relative differences are lower than 3%, and the relative deviation of the widthΓis approximately 10%.Substituting Eq.(2) into Eq.(3), the isotope yield ratio can be written as

Fig.1 (Color online) PNCS distributions of four nuclides 63Cu ,88Sr , 127 I and 197Au.Experimental data in the figure are taken from Refs.[24—27].The hollow black circle and hollow pink triangle indicate the ( γ , n) cross sections obtained with TALYS(version 1.8) calculations and TENDL (2021) evaluations[28], respectively.The red line shows the Lorentzian curve fitted with Eq.(2)

Fig.2 (Color online) PNCS distributions of 62-64Cu isotopes(a), 87-89Sr isotopes (b), 126-128 I isotopes (c), and 196-198 Au isotopes (d) computed with TALYS software (version 1.8).Note that experimental data for these PNCS distributions are not available, and the TALYS calculations are used for showing the difference of the PNCS distributions of adjacent isotopes

Table 1 Lorentzian parameters for photoneutron reactions on Cu, Sr,I and Au nuclei

Because of the complexity of the integral ratio in Eq.(4), we cannot deduce the integral ratio directly; thus, an approximate expression is provided.To verify the validity of the approximation, the integral ratio was calculated using the summation method (i.e., the integral was replaced by the summation).Considering the slope temperature of the photon spectrumTγ=5 MeV/kB(which can be readily achieved by state-of-the-art laser plasma-based bremsstrahlung sources [29, 30]), the integral ratio is approximated to be 1.0 and the root mean square deviation is 0.067, which is calculated over a wide range ofAandZ(note that 26

After irradiating a target composed of isotopesAandBwith laser-generated bremsstrahlung photons,γ-ray counting was performed for residual nucleiA′andB′.Theγlines unique to the residual nuclei are clearly identified in theγ-ray spectra, which are generally measured using an energyand efficiency-calibrated high-purity germanium (HPGe)detector.The reaction yieldsYA′andYB′were obtained using theγ-ray spectrum.For example, the expressions forYA′orYB′are as follows [31]:

whereNdetis the peak count of characteristicγrays used to identify the residual nuclei,tris the time for real measurement,tmis the measurement time,tcis the cooling time,Iγis the branching intensity,εis the source-peak detection efficiency, andis the correction factor for the selfabsorption of the characteristicγrays in the sample thicknessdwith absorption coefficientμ.According to Eq.(4),the peak cross sectionsσm,Acan be obtained as follows:

This indicates that when the distribution of the PNCS is known for isotopeB(orA), the distribution of the PNCS for isotopeA(orB) can be readily extracted if their peak energiesEmand widthsΓare very close to each other.The above method used to extract the PNCS distribution is valid only for nuclei with a single GDR peak.In the case of deformed nuclei, the GDR peak may split into two, where each corresponds to the major or minor axis of the ellipsoid [19],and the resultingσ(E) is given by the sum of the two GDR peaks, which complicates the proposed extraction method.

3 Benchmarking results

To verify the aforementioned extraction method using gamma activation and the isotope yield ratio, experiments for generating62,64Cu and85m,87mSr isotopes via laserinduced photonuclear reactions were performed at the Xing-Guang III laser facility of the Laser Fusion Research Center(LFRC) at Mianyang.Large-charge megaelectronvolt electron (e-) beams were generated using 100 TW picosecond(ps) laser pulses.Thee-beams then impinge on metal stacks composed of a Ta foil and activation plates (Cu or SrCl target) of interest, producing high-energy bremsstrahlung radiation and isotopes62,64Cu or85m,87mSr.The bremsstrahlung radiations produced have a Boltzmann-like distribution withTγ≈ 4.6 MeV/kB.The obtained reaction yields ofY62,64Cuare obtained 4.86 × 106and 2.65 × 106per laser shot, respectively [32].The reaction yieldsY85m,87mSrwere determined to be 1.8 × 104and 2.1 × 105per laser shot, respectively.In the following section, we introduce how to obtain the peak cross sections of the isotopes of interest and then predict the PNCS distributions for two exemplary isotopes,63Cu and86Sr.

3.1 PNCS distribution of 63Cu

The half-lives of62,64Cu isotopes are 9.7 min and 12.7 h,respectively.They have principle characteristic emissionsEγ=511.0 keV with branching intensityIγ= 195.7% for62Cu and ofIγ=35.2% for64Cu.Although the62,64Cu isotopes have the same characteristic emissions at 511 keV,their reaction yields can be reasonably extracted because of the significant differences in their half-lives [32].According to the experimental reaction yieldsY62,64Cuand the natural abundance of Cu isotopes, the ratio of the peak cross sectionσm,63Cu/σm,65Cuwas obtained as 0.82.Experimental crosssectional data [17] and evaluated data (IAEA-2019) [33] are available for65Cu(γ, n) reaction, as shown in Fig.3a.The Lorentzian-like shape of Eq.(2) was used to fit the data.The Lorentzian parameters (Em,65,63Cu,σm,65,63Cu, andΓ65,63Cu)were then obtained (see Table 2).Using experimental and evaluated PNCS data, theσm,65Cuare fitted to be 96.63 mb and 93.03 mb, respectively.According to Eq.(6), the peak cross sectionsσm,63Cuwere obtained as 79.24±5.44 mb and 76.28±5.24 mb, respectively.Finally, the PNCS distribution for63Cu can be predicted by the three Lorentzian parametersEm,65Cu,σm,63Cu, andΓ65Cu(orΓ63Cuwhen available).Figure 3b shows two PNCS distributions for63Cu (i.e., PNCS #1 and PNCS #2), which were extracted from the Lorentzian parameters fitted from the experimental and evaluated data of63Cu ,respectively.The experimental and evaluated PNCS data for63Cu are shown in Fig.3b.The uncertainty of the predicted PNCS distribution for63Cu was obtained using the error propagation formula, considering a statistical uncertainty of 1.5% of the measuredY62,64Cuand a systematic uncertainty of 6.7% induced by the proposed extraction method.PNCS #1 and PNCS #2 have an uncertainty of 6.9%.The benchmark results show that the extracted PNCS distribution is consistent with the evaluated data (IAEA-2019).Consequently, the validity of the proposed extraction method was confirmed.

Fig.3 (Color online) a Experimental and evaluated PNCS data for 65Cu [25, 33] and the resultant fitted Lorentzian-like distributions and b predicted PNCS distributions for 63Cu and the experimental and evaluated PNCS data for comparison.PNCS #1 and PNCS #2 represent the PNCS distributions extracted from the Lorentzian parameters fitted from experimental and evaluated data, respectively.The shadow area indicates the uncertainty of the extracted PNCS distributions for63Cu

3.2 PNCS distribution of 86Sr (leading to 85mSr isotope production)

The half-lives of85m,87mSr isotopes are 67.6 min and 2.8 h,respectively.They have principle characteristic emissionsEγ=231.9 keV with branching intensityIγ=83.9% for85mSr andEγ=388.5 keV withIγ=82.2% for87mSr.Similarly, the ratio of the peak cross sectionσm,86Sr/σm,88Srwas 0.72.Because experimental PNCS data were not available for88Sr(γ, n)87mSr and86Sr(γ, n)85mSr reactions, TALYS calculations were performed, and the calculated curves are shown in Fig.4a.These curves were then fitted to a Lorentzian-like shape, and the obtained Lorentzian parameters (Em,86,88Sr,σm,86,88Sr, andΓ86,88Sr) are listed in Table 3.The fittedσm,88Sris 163.13 mb.According to Eq.(6), the peak cross sectionσm,86Srwas obtained to be 117.11 ±7.96 mb, respectively.Finally, the PNCS distribution of86Sr can be predicted using the following three Lorentzian parameters:Em,88Sr,σm,86SrandΓ88Sr.The results are presented in Fig.4b.The cross-sectional data for the86Sr(γ, n)85mSr reaction were calculated using the TALYS software and are also shown in Fig.4b.Similarly, the uncertainty in the predicted PNCS distribution for86Sr is 6.8%.The extracted PNCS distribution was consistent with the TALYS calculations, validating the proposed extraction method.

Table 2 Lorentzian parameters fitted using experimental and evaluated PNCS data

Fig.4 (Color online) a TALYS-calculated PNCS data for 88Sr [34]and the resultant fitted Lorentzian-like distributions, and b predicted PNCS distributions for 86Sr and the TALYS-calculated PNCS data for comparison.PNCS represents the PNCS distribution extracted from the Lorentzian parameters fitted from TALYS data.The shadow area indicates the uncertainty of the predicted PNCS distribution for 86Sr

Table 3 Lorentzian parameters fitted using TALYS-calculated PNCS data

Fig.5 (Color online) 30 isotopes with corresponding PNCS distributions that are suitable for extraction using our proposed method.The transversal and longitudinal coordinates indicate the name of these isotopes and their natural abundances,respectively.The color bar presents the half-lives of the residual nuclei

Table 4 Detailed information for 30 isotopes with corresponding PNCS distributions that are suitable for extraction using our proposed method.A large uncertainty means the experimental data have an uncertainty larger than 10%

4 Discussions

The proposed extraction method has the following advantages: First, an accurate knowledge of the bremsstrahlung spectrum induced by intense lasers is not required because the reaction yield ratio method is used.When employing intense laser-accelerated electrons to produce bremsstrahlung radiation, they can still be subtracted using our extraction method although shot-to-shot variations exist in the bremsstrahlung spectrum and intensity.Second, it is not necessary to scan the incident photon energy, which generally happens in both the apparatus method using a quasi-monoenergetic photon beam and the mathematical method using classical bremsstrahlung measurements.Third, PNCS distributions can be successfully extracted within a very short time, within seconds, when using a stateof-the-art laser plasma-based bremsstrahlung source [29,30].Finally, it should be noted that any gamma source with a continuous spectrum covering the central GDR region is suitable for gamma activation.In general, gamma rays can be produced using either bremsstrahlung or LCS.When the spectral distribution of the gamma source isnγ(E)=nγ(0)f(E) , wheref(E) is an energy-dependent function, Eq.(4) can then be written asThe approximations given in Eq.(4) remain valid because of the similar values ofEmandΓfor the same nuclides as adjacent A, as discussed above.

However, the PNCS distributions of the studied isotopes must have Lorentzian-like shapes.This was a prerequisite for validating the proposed extraction method.Here, it is interesting to discuss how precisely the PNCS distribution can be described by a Lorentzian shape.Generally, a Lorentzianlike shape accurately describes the PNCS distribution in the central GDR energy range but fails in the high-energy tail[16].In our study, the relative difference between the predicted PNCS (integrated over the 10—25 MeV energy range)and the benchmarked PNCS was determined to be 3.2% for63Cu and 3.9% for86Sr.However, it increases to more than 50% in the high-energy range of 20—25 MeV.Meanwhile,successful extraction of the PNCS distribution relies on the fact that the residual nuclide should have a suitable halflife (note that a relatively long half-life is also acceptable),which is helpful for detecting characteristic emissions after target activation.To obtain a high reaction yield and reduce statistical uncertainty, the offline-detected characteristic emissions should have relatively large branching intensities and detection efficiencies.In addition, isotopes undergoing photonuclear reactions should have acceptable natural abundances because natural metal targets are commonly used for irradiation.Based on these considerations, the PNCS distributions of 30 isotopes were predicted using the proposed extraction method.Figure 5 shows the natural abundances of these isotopes and half-lives of the residual nuclei.Approximately half of them had PNCS distributions, leading to isomeric states of the residual nuclei.More detailed information regarding the target and residual nuclides is presented in Table 4.Their experimental PNCS data for 18 isotopes exhibited large uncertainty (>10%).The PNCS data for85Rb ,107Ag ,136Xe , and153Eu have acceptable uncertainties because these isotopes are categorized as either medical or other types of materials [1], which requires an uncertainty of less than 10%.

5 Conclusion

We proposed an effective method for obtaining the peak cross section and extracting the PNCS distribution using the gamma activation and reaction yield ratio method.The PNCS distributions of63Cu and86Sr isotopes (leading to85mSr isotope production) were successfully extracted using the proposed extraction method.The uncertainty of this method was maintained within 7%.The extracted PNCS distributions and the available PNCS data (or TALYS calculations) were compared to validate the proposed extraction method.The prerequisites and merits of the extraction method and the possibility of predicting PNCS distributions for 30 isotopes were also discussed.The proposed extraction method could be complementary to the available sophisticated methods, including the mathematical method using bremsstrahlung measurements and the apparatus method with positron annihilation in flight and LCS photon spectrum.In the near future, we plan to perform gamma activation experiments using laser plasma-driven bremsstrahlung sources and extract the PNCS distributions for a list of isotopes with relatively large abundances and significant characteristic emissions after short-term irradiation.

Acknowledgements We express our gratitude to the XingGuang III operation team for their support.

Author contributions All authors contributed to the study conception and design.Material preparation, data collection and analysis were performed by Zhi-Cai Li, Yue Yang, Zong-Wei Cao, Xin-Xiang Li,Yun Yuan, Zong-Qing Zhao, Gong-Tao Fan, Hong-Wei Wang and Wen Luo.The first draft of the manuscript was written by Zhi-Cai Li,and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.

Data availability The data that support the findings of this study are openly available in Science Data Bank at https://doi.org/10.57760/sciencedb.j00186.00302 and https://cstr.cn/31253.11.sciencedb.j00186.00302.

Declarations

Conflict of interest Hong-Wei Wang is an editorial board member for Nuclear Science and Techniques and was not involved in the editorial review, or the decision to publish this article.All authors declare that there are no competing interests.