Photonuclear production of nuclear isomers using bremsstrahlung induced by laser‑wakefield electrons

Hao‑Yang Lan ·Di Wu ·Jia‑Xin Liu ·Jian‑Yao Zhang ·Huan‑Gang Lu ·Jian‑Feng Lv ·Xue‑Zhi Wu ·Wen Luo ·Xue‑Qing Yan

Abstract In this study, we theoretically investigate the feasibility of using laser-wakefield accelerated (LWFA) electrons for the photonuclear measurement of nuclear isomers according to the characteristics of the electrons obtained from LWFA experiments conducted at the Compact Laser–Plasma Accelerator (CLAPA) laboratory.The experiments at the CLAPA show that a stable electron beam with an energy of 78–135 MeV and a charge of 300–600 pC can be obtained.The bremsstrahlung spectra were simulated using Geant4, which suggests that a bremsstrahlung source with a peak intensity of 1019 photons/s can be generated.Theoretical calculations of isomer production cross sections from the photonuclear reactions on six target nuclei,197Au, 180Hf, 159Tb, 115In, 103Rh, and 90Zr, were performed and compared with the available experimental data in EXFOR,which suggest that further experiments are required for a series of photonuclear reaction channels.Flux-averaged cross sections and isomer ratios (IR) resulting from such bremsstrahlung sources are theoretically deduced.The results suggest that IR measurements can be used to constrain nuclear components, such as γ strength function and optical model potential.In addition, the detection of the decay characteristics was evaluated with Geant4 simulations.The use of the LWFA electron beam and its bremsstrahlung for photonuclear studies involving nuclear isomers is anticipated.

Keywords Photonuclear reactions ·Laser–plasma acceleration ·Flux-averaged isomer ratio

1 Introduction

Since the theoretical prediction by Soddy [1] in 1917 and the first experimental observation by Hahn [2] in 1921, nuclear isomers have received intensive attention from the perspective of fundamental research and applicational development[3, 4].In particular, the effective population and manipulation of nuclear isomers has been an active topic in recent years, driven by the prospect of utilizing their intrinsic or depleted decays for the development of nuclear clocks,nuclear batteries, and gamma lasers [5–9].Moreover, the exploration of nuclei with isomeric states in nuclear astrophysics allows us to determine how they affect the creation of elements in the universe and how they eventually contribute to the makeup of life in our cosmos [10].

In recent decades, the worldwide development of highintensity and high-repetition-rate laser systems has opened up new research opportunities for particle acceleration[11–17] and nuclear photonics [18].In particular, the ultrashort pulse radiation source and plasma environment generated by the interaction between the laser and target provide unique physical environments for the study of nuclear isomers.Using a table-top hundred TW laser, Feng et al.[19]demonstrated that the femtosecond pumping of83mKr with a peak efficiency of 2.34 × 1015s−1can be achieved through the Coulomb excitation of ions with quivering electrons during the interaction between laser pulses and clusters at nearly solid densities.Wang et al.[20] proposed a novel approach to populate the nuclear clock isomer229mTh by combining laser-driven electron recollision and electronic nuclear excitation, for example, Coulomb excitation, which does not require precise knowledge of the isomeric energy.The common idea of these schemes is to make use of the dense low-energy electrons produced in laser–plasma interactions and the relatively high electronic excitation cross sections of ionized atoms in plasma.A comprehensive comparison of the possible physical processes for isomer population and manipulation, including nuclear excitation by electron capture, excitation by electron transition, electron bridge, and photoactivation, was also given in Ref.[21],which indicates that nuclear excitations assisted by electrons would be far more efficient than straightforward photon excitation.This is reasonable for isomer manipulation or triggered isomer depletion circumstances, where the target nucleus must remain unchanged, and the excitation energy required is as low as the difference between the atomic electron shells.However, this is somehow not the case when it comes to the isomer population, where the isomer can be populated without extra scruples on the disintegration of the target nucleus.Specifically, the cross section of photonuclear reactions can be significantly enhanced by giant dipole resonances (GDR) to hundreds of millibarns above the neutron separation threshold, whereas only the narrow discrete resonance from the isolated bound states is generally mentioned in most of the discussions for photoactivation of nuclear isomers.Utilizing (γ,γ′) scattering of the collective excitation mode, Pan et al.[22] have depicted the prospect of accumulating isomers of medical interest to a clinically applicable activity with an intense monochromaticγ-ray beam.Meanwhile, quasi-monochromatic high-intensity electron pulses generated by laser-wakefield acceleration(LWFA) in gas jets are commonly available in high-power laser facilities, which can be converted into bremsstrahlung radiation, thus providing vast research opportunities for photonuclear studies.Although previous studies have utilized laser-driven photonuclear reactions for electron diagnosis[23], nuclear waste transmutation [24], and medical isotope production [25], little attention has been paid to nuclear isomers.Consequently, a systematic investigation is required to evaluate the potential of the isomer population by laserdriven photonuclear excitation.

In this work, we theoretically investigate the potential to accumulate nuclear isomers by laser-driven photonuclear excitation through Monte Carlo simulations based on the electron characteristics obtained from LWFA experiments.LWFA experiments at the Compact Laser–Plasma Accelerator (CLAPA) laboratory of Peking University are introduced,and the conversion from electrons to bremsstrahlung is simulated with the Geant4 [26].Then, theoretical calculations on the photonuclear cross sections are performed.A series of nuclear isomers were selected as the calculation candidates,considering both the nuclear flow of photonuclear reactions and the availability of stable natural target nuclei.The corresponding photonuclear cross sections for isomer generation were evaluated using TALYS [27, 28].To verify the reliability of the TALYS inputs, a benchmark of the theoretical predictions against the available experimental data extracted from the experimental nuclear reaction database,EXFOR [29] is performed.The flux-averaged cross sections and flux-averaged isomer ratios were deduced according to the theoretical photonuclear cross sections and simulated bremsstrahlung spectra.Furthermore, the photonuclear cross sections of different reaction channels were incorporated into the simulation toolkit Geant4-GENBOD [30–32]to estimate the feasibility of detecting the nuclear isomers of interest.In the simulations, the typical laser-driven electron beam parameters extracted from the experiments are considered, as well as the physical processes of the radioactive decay of various reaction residuals and energy deposition in commercially available LaBr3scintillation detectors.A comprehensive discussion is provided in the context of the production and detection of the nuclear isomers of interest.

2 LWFA electrons and bremsstrahlung conversion at CLAPA

We conducted experiments to accelerate electrons for bremsstrahlung generation with the 200 TW laser facility at the CLAPA of Peking University, which can deliver laser pulses with a duration of 30 fs, energy of 4 J, and repetition rate of 5 Hz.The laser beam was focused on a spot of 21 ×23 μm full width at maximum onto the leading edge of a 4 ×1 mm supersonic gas jet.The effective energy on the target was approximately 1.8 J with a laser intensity of approximately 1.5 × 10−19W/cm−2.In the interaction between the laser pulses and the mixed gas made of 5% N2and 95%He2, monochromatic energetic electrons can be produced by LWFA.The electrons propagate through a thin Be window,heading out of the vacuum chamber to the magnetic spectrometer placed 100 cm downstream of the gas jet, where the deflected electron signals can be diagnosed using fluorescent screens and CCD cameras.The electron charge was measured by a Turbo Integrating Current Transformer (Turbo-ICT) assembled right out of the Be window.

The back pressure of the gas jet was varied from 33 bar to 39 bar, obtaining LWFA electrons with a peak energy of 78–135 MeV.One hundred continuous shots were recorded for each gas pressure to test the stability of the accelerated electrons.The repetition rate remained at 0.25 Hz to ensure that the vacuum level of the laser–gas interaction chamber could be recovered for the LWFA.Typical electron distributions recorded by the magnetic spectrometer are shown in Fig.1a–c.Electron energies of 78 ± 10, 103 ±14, and 135 ± 20 MeV were obtained at 39, 36, and 33 bar,respectively.The divergence of the electrons was within 4.3 mrad, and the directional stability was better than 5.2 mrad.The electron charge was measured to be 300–600 pC.

The averaged energy and divergence distributions of the electrons obtained at 39, 36, and 33 bar were incorporated into Geant4 to simulate the resulting bremsstrahlung characteristics.In the simulations, the electrons were directed to a 2-mm-thick Ta converter 1 m away from the electron generation point, inducing bremsstrahlung radiation.The electron charge was set to 300 pC, which is a conservative estimate for the number of LWFA electrons produced per shot.The bremsstrahlung photons were recorded in a plane perpendicular to the propagation direction assembled 20 cm downstream from the converter.The simulated bremsstrahlung energy distributions are shown in Fig.1d.The photon intensity shows a typical decaying trend with the increase in photon energy, which can fully cover the GDR region of the selected target nuclei.The conversion efficiency from electrons to photons increases with increasing electron energy.The total photon flux can reach the order of 4.35 × 108photons/shot, which corresponds to a peak photon flux of 1.24× 1019photons/s, considering a pulse duration spread of 35 ps inside the converter.The correlation between the photon energy and polar angle for the incident electrons at 33 bar is shown in Fig.1e.It is shown that the angular spread of the photons tends to narrow with an increase in photon energy.The angular spread of most of the photons is within 0.4 rad.With such a photon beam with a high peak flux and ultrashort pulse duration, one can potentially measure the nuclear isomers with short half-lives produced in photonuclear reactions.Accordingly, the reaction data covering the GDR and quasi-deuteron regions were probed.

3 Production route and photonuclear cross sections

3.1 Theory

The nuclear flow of photonuclear reactions and mechanism for isomer generation are depicted in Fig.2.When a target nucleus is excited above the particle separation energy the excited compound subsequently decays through the emission of particles (such as neutrons, protons, and alpha particles) or photons to different energy levels of the residual,including the isomer state.At energies of 10–40 MeV, the GDR resulting from the out-of-phase oscillation between the proton fluid and neutron fluid would clearly give rise to the odds of photon absorption and compound nucleus formation.

Fig.1 (Color online) Panels a, b, and c are the exemplary LWFA electron spectra obtained from the back pressure of 33, 36, and 39 bar, respectively.d The bremsstrahlung energy spectrum obtained from the Geant4 simulations considering the averaged energy spectra and angular spread of 100 LWFA electron shots.e The correlation of bremsstrahlung energy against its angular spread for electrons at 33 bars

Fig.2 (Color online) The nuclear flow of photonuclear reactions and the mechanism for isomer generation.Here, Sn denotes the neutron separation energy of the target nucleus 180 Hf

At 40–130 MeV, the photon interacts with the nucleus by the quasi-deuteron mechanism (QD), where the compound nucleus decays predominantly via the emission of several nucleons, but mostly neutrons.The decay of the compound nucleus, which would lead to the population of isomers, can be well described by the Hauser–Feshbach reaction model,which is based on the fundamental Bohr hypothesis that the reaction occurs by means of the intermediary formation of a compound nucleus that can reach a state of thermodynamic equilibrium.The Hauser–Feshbach model averages over many resonances in the formation of a compound nucleus and therefore can be applied if the nuclear level density(NLD) in the compound nucleus is sufficiently high at the compound formation energy, which is the excitation energy at which the compound nucleus is formed.

3.2 Verification of photonuclear cross sections

Generally, photonuclear reactions tend to push the involved nucleus toward the down-left side of the nuclear chart.This places a clear limit on the possible populated isomers when considering the stability of the target nucleus.According to the nuclear flow of photonuclear reactions and the stability of the target nucleus, we selected a series of stable nuclides, including197Au,180Hf,159Tb,115In,103Rh, and90Zr as candidates for theoretical calculations of reaction cross sections.The cross sections for typical photonuclear reactions, such as (γ,γ′), (γ,xn),(γ, p), and (γ,α), are calculated using the software package TALYS, which provides a complete description of all reaction channels and observables, and many state-of-theart nuclear models covering all the main reaction mechanisms encountered in light-particle-induced nuclear reactions are included.The nuclear structure components used in the TALYS calculations were as follows: The nuclear masses are taken from the Atomic Mass Evaluation 2016(AME2016) [67] whenever available, whereas the HFB-27 nuclear masses [68] are taken into account when the AME2016 mass data are not available.The discrete experimental levels compiled in the RIPL-3 library [69] and the continuum level spectrum represented by the nuclear level densities (NLDs) are both considered in the calculations.The NLDs were obtained from predictions in the Fermi gas model plus a constant temperature.The photon strength functions (SFs) obtained from the Brink–Axel Lorentzian model were used to calculate the electromagnetic transmission coefficients for the photon channel.The phenomenological optical model potentials (OMPs),i.e., the local and global parameterizations of Koning and Delaroche, are employed to determine the transmission coefficient for the particle channels.

The theoretical photonuclear reaction cross sections based on TALYS for the selected targets and their corresponding experimental data available in the EXFOR database are shown in Fig.3.It can be seen that the (γ, n) reaction cross section is generally dominant in the GDR region, with a maximum of hundreds of mb.As the excitation energy continues to increase, multiple emission reactions take over,whose maximum cross section decreases with an increase in the number of emitted particles.The theoretical calculations agree with the EXFOR data, which is obtained from previous measurements using photon sources, such as (tagged)bremsstrahlung [33–39], laser-Comptonγrays [40–44], and positron annihilation in flight [45–53].This indicates the reliability of the TALYS inputs and calculations.In addition,it is noticeable that the previous measurements have not fully covered all the photonuclear reaction channels, although the(γ,xn) reaction data for the most studied nucleus197Au are rather abundant.Future experiments are expected to provide new data for these reaction channels.The isomer production cross sections resulting from photonuclear reactions on the six selected target nuclei are shown in Fig.4.Isomer production from (γ, n) reactions is generally greater than that from the other reaction channels, with a maximum cross section of hundreds of ∼100 mb.It is also noted that experimental isomer production cross sections are very rare,which can be potentially measured using state-of-the-art LCS facilities [70–72] in the future.

Fig.3 (Color online) The calculated ( γ , xn) cross sections for 197Au, 180Hf, 159Tb, 115In, 103Rh, 90Zr.The available experimental data from the EXFOR database are also shown for comparison [33–53]

Fig.4 (Color online) The isomer production cross sections for the six selected target nuclei.The available experimental data from the EXFOR database are also shown for comparison [54–60]

4 Production yields and flux‑averaged isomer ratios

4.1 Isomer production yields

In bremsstrahlung-driven photonuclear measurements,isomers and their decay have long been interesting topics because they can provide constraints for nuclear structures.Under bremsstrahlung radiation, the flux average cross sections (FACS) of a given photonuclear reaction can be given by:

whereσ(E) is the reaction cross section at photon energyE,Φ(E) is the photon flux,Eeis the electron energy, andEthis the photonuclear reaction threshold.

Using the bremsstrahlung driven by laser–plasmaaccelerated electron beams described in Sect.2, we conducted experimental measurements on the (γ,xn) cross sections of197Au and compared them against the theoretical estimations based on Eq.1.The experimental details and data analysis can be found in Ref.[73].The agreement between experimental results and theory was found, which demonstrates that the procedure of measuring FACS with laser-driven bremsstrahlung is feasible.With this premise, we further evaluated the expected FACS of possible isomer residuals when laser-driven bremsstrahlung beams were used for the photoactivation of197Au,180Hf,159Tb,115In,103Rh, and90Zr, as shown in Fig.5.Noticeably, for the FACS to produce the nuclear isomers196m1Au,195mAu,179m1Hf,158mTb,114mIn,101mRh, and89mZr were relatively larger than the others, exceeding the value of 10 mb.Moreover, the available experimental FACS to produce nuclear isomers for these target nuclei is very limited.The evaluated FACS for the production of196m2Au from the(γ, n) reaction on197Au agrees well with the experimental measurements made by Naik et al.[61], while those for the production of102m,101m,100m,99mRh produced from(γ,xn) reactions on103Rh are all lower than the experimental values obtained by Rahman et al.[62].The FACS for the isomers of the other isotopes produced from reactions, such as (γ, p) and (γ,α) reactions, are also shown in Supplementary Material, which are less abundant in the experimental data.These results show that experimental efforts to measure FACS involving isomer residuals are expected to resolve discrepancies and deficiencies in photonuclear data in the QD region.

Fig.5 (Color online) The FACS for ( γ , xn) reactions for the six selected target nuclei.The available experimental data from the EXFOR database are also shown for comparison [61, 62]

4.2 Dependence of isomer ratio on nuclear structure ingredients

The ratio of the probabilities of forming the isomer and the ground state is called the isomeric ratio, which can be given by:

The results are predicted using six sets of NLDs (1–Constant temperature Fermi gas model [74], 2–Back-shifted Fermi gas model [75], 3–Generalized superfluid model[76], 4–HFB-Skyrme model [77], 5–HFB-Skyrme model with combinatorial method [78], and 6–Temperature(T)-dependent HFB-Gogny model [79]), two sets of OMPs(0–phenomenological Wood–Saxon potential [80], 1–semi-microscopic JLMB potential with HFB-Skyrme matter density [68, 81, 82, 83]), and eight sets of PSFs (1–Generalized Lorentzian [84], 2–Brink–Axel Lorentzian[85, 86], 3–HFBCS-QRPA model [87], 4–HFB-Skyrme-QRPA model [88], 5–Hybrid model [89], 6–T-dependent HFB-Skyrme-QRPA model [88, 90], 7–T-dependent RMF model [91], and 8–HFB-D1M-QRPA model [92]) which are shown in Fig.6.It is shown that the IRs grow rapidly with the increase in incident electron energy at low energies and tend to saturate at great energies of tens to hundreds MeV, where discrepancies resulting from different models occur.Specifically, IRm1/gis not sensitive to the three nuclear ingredients, whereas both IRm2/(m1+g)and IRm2/gshow relatively strong sensitivity to the SFs and NLDs.This suggests that the existing experimental data and further measurements aimed at IRs, especially in high-energy regions, can be used to constrain these nuclear structure ingredients.

Fig.6 (Color online) The isomer ratios of IRm2/g (a–c, IRm2/g d–f, and IRm2/(m1+g) g–i for the 197Au(γ , n) reactions calculated with six sets of eight sets of SFs, two sets of OMPs, and NLDs available in TALYS 1.9.The nuclear structure ingredients used for the calculations are indexed by numbers and described in the main text.The experimental data from EXFOR are also shown for comparison [61, 63, 64, 65, 66]

5 Detection of characteristic decays

In photonuclear experiments, there are two ways to measure the reaction cross sections: online detection of photoinduced particle emissions (neutrons, protons, andαparticles) and offline detection of the characteristicγdecays of radioactive residuals.The former is generally used in experiments conducted with energy-tunable, quasi-monochromatic probes, such as laser-Comptonγ-ray sources and positron annihilation in flight because the energy spectra of the emitted particles would exhibit easily identifiable reaction kinematics.In bremsstrahlung-driven photonuclear experiments, offline detection is typically performed.The detection of residuals depends not only on their decay characteristics and production yields but also on the energy resolution of the detector.

We performed Geant4 simulations to evaluate the feasibility of detecting the considered isomers produced in photonuclear reactions.The simulation setup is arranged as follows.The energy distribution and angular spread of the bremsstrahlung produced by the LWFA electrons are set according to the experimental results obtained at a gas pressure of 36 bar.The electron charge was set as 100 shots at 300 pC.The bremsstrahlung photons irradiated the activation targets with a thickness of 5 mm, placed 20 cm downstream of the converter, inducing photonuclear reactions and producing radioactive residuals.Theγdecays of the radioactive residuals, including the isomers, were recorded by a LaBr3detector with dimensions ofφ8 cm × 10 cm assembled at 90◦with respect to the beam propagation.The energy resolution is set as 3%, which is achievable for commercially available LaBr3detectors.To simulate isomer production and decay radiation detection in a one-stop manner, we adopted the Geant4-GENBOD [30–32] toolkit and further implemented the radioactive decay process of the radioactive residuals (including the isomers).The photonuclear reaction cross sections of the related isotopes were calculated with a routine similar to Sect.3.2.

Fig.7 (Color online) a The recorded γ-ray spectra for the photonuclear products of 90 Zr when the cooling time tc is 1 ms, 10 s, 10 m,1 h, and 1 d.The data acquisition time is set to 14 days.Here, the offline detection of delayed radiation is imitated by applying a time window for the energy deposition in the detector.b The dependence of recorded counts on the energy deposition and detection time

Table 1 The expected counts,half-lives, γ-ray energies, and intensities of the principal decays from the photonuclear products of 197Au, 180Hf, 159Tb,115In, 103Rh, and 90Zr

Table 1 (continued)

An exemplary correlation of energy deposition versus detection time for the radioactive residuals induced by photonuclear reactions on90Zr is given in Fig.7.In activation experiments based on laser–plasma acceleration, there is a period (cooling timetc) for opening the reaction chamber before the activated samples are handed to the detector for data acquisition.Iftcis larger than 1 h and 10 s, the decay information of the products89mZr (T1∕2= 4.16 m) and90mZr(T1∕2= 809 ms) will be lost.To detect isomers with short half-lives, such as90mZr the detector is placed near the activation target in an “online”manner.The backgrounds induced by the backscattered bremsstrahlung, electrons, and electromagnetic pulses can be cut off by applying a time window with a delay signal generator.The correlation of energy deposition versus detection time for the radioactive residuals induced by the photonuclear reactions on all six selected targets is presented in Supplementary Material.The principal decay characteristics and expected counts within three half-lives for the isomers produced from the photonuclear reactions of the six selected isotopes are shown in Table 1.The photonuclear products196m1Au,194m1Au,157Tb,156m1Tb, and103mRh are not measurable for the LaBr3detector because their decay energy is too little and the attenuation in the target is strong.The measurement of195m2Au is unlikely because its decay energy and intensity has not yet been determined.For180Hf, the decay of180mHf is difficult to measure because the (γ,γ′) reaction cross section is very low.In addition, most of the isomer residuals can be measured with a net count of 102−6within 14 days of data acquisition.The decay energies of some products, such as193Au and193mAu, are very close and cannot be disentangled from the recorded energy spectrum of a detector with poor energy resolution.However, this can be remedied by applying a fit to the recorded counts in a few data acquisition cycles.

6 Conclusion

In this study, we calculated the production cross sections of nuclear isomers from the photonuclear reaction on197Au,180Hf,159Tb,115In,103Rh, and90Zr with TALYS.Deficiencies in the corresponding experimental data were found in the comparison between theoretical estimations and the available data from EXFOR.Laser–plasma electron acceleration experiments carried out at CLAPA are introduced.Stable LWFA electrons with energies of 78 ± 10, 103 ±14, and 135 ± 20 MeV can be obtained with a charge of 300−600 pC per shot.The resulting peak bremsstrahlung flux is expected to be 1.24 × 1019photons/s.To evaluate the prospect of using such a photon beam for photonuclear studies regarding nuclear isomers, the FACSs and IRs for isomer production were deduced and compared with the data found in EXFOR.The sensitivity of IRs against the nuclear ingredients SF and NLD was determined, which highlights the importance of measuring photonuclear reactions involving isomers.Moreover, the prospect of photonuclear measurement of nuclear isomers using LWFA electrons at CLAPA is foreseen with Geant4 simulations.

AcknowledgementsThe authors thank the staff of the 200 TW laser of CLAPA laboratory for the smooth operation of the machine.

Author contributionsAll authors contributed to the study conception and design.Material preparation, data collection, and analysis were performed by Hao-Yang Lan, Di Wu, Jia-Xin Liu, Jian-Yao Zhang,Huan-Gang Lu, Jian-Feng Lv, and Xue-Zhi Wu.The first draft of the manuscript was written by Hao-Yang Lan and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.

Data availabilityThe data that support the findings of this study are openly available in Science Data Bank at https:// doi.org/ 10.57760/ scien cedb.j00186.00066 and https:// cstr.cn/ 31253.11.scien cedb.j00186.00066.

Conflict of interestXue-Qing Yan 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.