Photoneutron cross section measurements on ${}^{208}$ Pb
in the Giant Dipole Resonance region

I. Gheorghe ioana.gheorghe@nipne.ro National Institute for Physics and Nuclear Engineering, Horia Hulubei (IFIN-HH), 30 Reactorului, 077125 Bucharest-Magurele, Romania S. Goriely Institut d’Astronomie et d’Astrophysique, Université Libre de Bruxelles, Campus de la Plaine, CP-226, 1050 Brussels, Belgium N. Wagner Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany T. Aumann Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany GSI Helmholtzzentrum für Schwerionenforschung, 64291 Darmstadt,Germany Helmholtz Forschungsakademie Hessen für FAIR (HFHF), GSI Helmholtzzentrum für Schwerionenforschung, 64291 Darmstadt, Germany M. Baumann Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany P. van Beek Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany P. Kuchenbrod Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany H. Scheit Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany D. Symochko Institut für Kernphysik, Technische Universität Darmstadt, Darmstadt, 64289, Germany T. Ari-izumi Konan University, Department of Physics, 8-9-1 Okamoto, Higashinada, Kobe 658-8501,Japan F.L. Bello Garrote Department of Physics, University of Oslo, N-0316 Oslo, Norway T. Eriksen Department of Physics, University of Oslo, N-0316 Oslo, Norway W. Paulsen Department of Physics, University of Oslo, N-0316 Oslo, Norway L.G. Pedersen Department of Physics, University of Oslo, N-0316 Oslo, Norway F. Reaz Department of Physics, University of Oslo, N-0316 Oslo, Norway V. W. Ingeberg Department of Physics, University of Oslo, N-0316 Oslo, Norway S. Belyshev Lomonosov Moscow State University, Faculty of Physics, 119991 Moscow, Russia S. Miyamoto Laboratory of Advanced Science and Technology for Industry, University of Hyogo, 3-1-2 Kouto, Kamigori, Ako-gun, Hyogo 678-1205, Japan H. Utsunomiya hiro@konan-u.ac.jp Konan University, Department of Physics, 8-9-1 Okamoto, Higashinada, Kobe 658-8501, Japan

(March 18, 2024)

Abstract

Photoneutron reactions on ${}^{208}$ Pb in the Giant Dipole Resonance energy region have been investigated at the $\gamma$ -ray beam line of the NewSUBARU facility in Japan. The measurements made use of quasi-monochromatic laser Compton backscattering $\gamma$ -ray beams in a broad energy range, from the neutron threshold up to 38 MeV, and of a flat-efficiency moderated ${}^{3}$ He neutron detection system along with associated neutron-multiplicity sorting methods. We report absolute cross sections and mean photoneutron energies for the ${}^{208}$ Pb $(\gamma,\,inX)$ reactions with $i$ = 1 to 4. The fine structure present in the ${}^{208}$ Pb $(\gamma,\,n)$ cross sections at incident energies lower than 13 MeV has been observed. The photoabsorption cross section has been obtained as the sum of the $(\gamma,\,inX)$ reaction cross sections. By reproducing the measured ring-ratio values at excitation energies below the two neutron separation energy, we were able to extract estimations on the ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron energy spectra and on the partial photoneutron cross sections for leaving the residual ${}^{207}$ Pb in its ground and first two excited states. The present results are compared with data from the literature and statistical model calculations.

^†^†preprint: APS/123-QED

I Introduction

Photonuclear data, describing the response of atomic nuclei to photons, find use in fundamental nuclear physics and in a wide range of applications [1]. The photoabsorption cross section is used to directly determine the $\gamma$ -ray strength function ( $\gamma$ SF) [2], the key ingredient for computing $\gamma$ -ray cascades in nuclear reactions [3, 4]. Nuclear reactions induced by high energy photons of 10–20 MeV are dominated by $E1$ excitations known as the Giant Dipole Resonance (GDR) and understood as collective oscillations of the protons against the neutrons [5]. The photoabsorption cross section in the GDR energy range reveals information on important nuclear quantities [6], such as the symmetry energy, which is the restoring force against the separation of protons and neutrons.

The doubly magic ${}^{208}$ Pb is a benchmark case for theoretical modeling of the electric dipole response in nuclei [6]. It has been extensively investigated through photonuclear reactions in pioneering experiments at the Livermore [7, 8] and Saclay [9] positron in flight annihilation facilities. Recently, fine structures superimposed on the broad GDR resonance have been observed in great detail in high energy resolution inelastic hadron scattering experiments at RCNP GrandRaiden [10, 11] and iThemba [12] facilities, revealing single particle manifestations. Based on the $(p,\,p^{\prime})$ RCNP data and on ${}^{208}$ Pb $(^{3}$ He, ${}^{3}$ He ${}^{\prime}\gamma)$ Oslo-type data, the total ${}^{208}$ Pb $\gamma$ SF and its $E1$ , $M1$ and $E2$ components have been extracted in the 2.7 – 20 MeV excitation energy range [13].

However, there are unresolved systematic discrepancies [16] between photoabsorption and photoneutron cross sections obtained at the Saclay and Livermore facilities, with Saclay [9] $(\gamma,\,abs)$ and $(\gamma,\,n)$ cross sections systematically higher than the Livermore ones [7], and Livermore $(\gamma,\,2n)$ cross sections overestimating the Saclay ones. The different neutron multiplicity sorting procedures employed at the two laboratories have been suggested in Refs. [14, 15] as a discrepancy source. The need to remeasure the ${}^{208}$ Pb photoabsorption excitation function in the GDR energy range has also been pointed out in Ref. [6] in order to shed light on a discontinuity in the $E1$ moments between neighboring ${}^{208}$ Pb and ${}^{209}$ Bi nuclei.

New measurements of photoneutron reactions on ${}^{208}$ Pb in and above the GDR energy range have been performed at the laser Compton scattering (LCS) $\gamma$ -ray beam line of the NewSUBARU synchrotron radiation facility in Japan [17, 18]. Photoneutron cross sections have been measured at incident photon energies between 7.5 and 38 MeV. The investigations made use of a novel high-and-flat efficiency moderated neutron detection array (FED) and the associated neutron-multiplicity sorting methods [19, 20].

Decay experiments which measure the energy and angular distribution of $\gamma$ -rays and particles emitted from GDR states are in high demand for bringing information on the microscopic nature of the GDR [21]. However, these are very scarce, because of low detection efficiencies of neutron time of flight experiments. In lack of direct spectroscopic measurements, the present experiment provides indirect determinations of the average energy of neutrons emitted in the photoneutron reactions on ${}^{208}$ Pb. In the present work, we extract estimations for the neutron emission spectra in the low energy ${}^{208}$ Pb $(\gamma,\,n)$ reaction.

In Sect. II, we present the experimental technique and methodology, with focus on diagnostics of the incident LCS $\gamma$ -ray beams and neutron detection. The data analysis methods concerning neutron multiplicity sorting, extraction of information on the photoneutrons mean energy and spectra, as well as the energy unfolding are discussed in Sect. III. Results are discussed and compared with preceding data in Sect. IV and with theoretical calculations in Sect. V. These include the photoabsorption cross section, the photoneutron cross sections and average energies and also estimations on the ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron spectra and partial cross sections. A summary and conclusions are given in Sect. VI.

II Experimental method

Refer to caption — Figure 1: Diagram showing the BL01 LCS $\gamma$ -ray beam line and the experimental hutch GACKO at the NewSUBARU synchrotron radiation facility (not to scale). P ${}_{1}$ and the P ${}_{2}$ mark the focus point positions of the 1064 nm and, respectively the 532 nm laser beams at 1.8 and, respectively 3.8 m downstream of the eletron beam focus in the middle of the BL01.

The apparatus used for the present experiment is shown schematically in Fig. 1. Laser beams were guided into the BL01 straight beam-line section of the NewSUBARU synchrotron and were scattered by electrons circulating in the ring. The back-scattered $\gamma$ -ray photons that passed through a system of two lead collimators (C ${}_{1}$ and C ${}_{2}$ in Fig. 1) irradiated the target which was located in the GACKO (Gamma collaboration hutch of Konan University) experimental hutch, 7.5 m downstream of the second collimator.

The energy resolution and flux of the $\gamma$ -ray beam were monitored with a lanthanum bromide (LaBr ${}_{3}$ :Ce) and a NaI:Tl detector, respectively. The reaction neutrons were moderated and recorded with an array of ${}^{3}$ He counters placed in concentric rings around the target. The data were written in a triggerless list mode, using an eight-parameter 25 MHz digital data acquisition system which recorded the time and energy of the signals provided by the LaBr ${}_{3}$ :Ce and NaI detectors, the neutron detection time in each ring of ${}^{3}$ He counters and the external laser triggering signal.

We note that the present experiment provides the $\sigma_{\gamma,\,inX}$ sum cross sections of all photoneutron reactions with $i$ neutrons in the final state, accompanied or not by charged particle emission:

	$\displaystyle\sigma_{\gamma,\,inX}$	$\displaystyle\equiv\sigma(\gamma,\,inX)$
		$\displaystyle=\sigma(\gamma,\,in)+\sigma(\gamma,\,inp)+\sigma(\gamma,\,in% \alpha)+\dots$		(1)

However, the large Coulomb barrier in ${}^{208}$ Pb ( $\sim$ 26 MeV for $\alpha$ particles) hinders the emission of charged particles. Thus, for incident photon energies covering the GDR range up to $\sim$ 20 MeV, both the cross sections and the average photoneutron energies for the $(\gamma,\,inX)$ reactions generally coincide to or are a good approximation for the $(\gamma,\,in)$ neutron emission only ones [1].

II.1 LCS $\gamma$ -ray beams

Table 1 gives the main properties of the electron, laser and $\gamma$ -ray beams employed in the present study. The $\gamma$ beam energy was varied in 85 steps between 7.5 and 38.12 MeV by using two lasers with wavelengths of 1064 and 532 nm and tuning the electron energy between approximately 650 and 1050 MeV. The maximum energy of the LCS $\gamma$ -ray beam is directly determined by the known laser wavelength and by the electron beam energy, which has been calibrated with an accuracy of (5.5 – 9.4) $\times$ 10 ${}^{-5}$ [22].

Interactions between the unsynchronized laser and electron beams could take place along the entire 20 m length of the electron beamline, according to the laser-electron beam spatial overlap (see Fig. 6 of Ref. [24]). The resulting LCS $\gamma$ -ray beam had a pulsed time structure given by the $slow$ laser time structure (1 – 20 kHz, 20 – 40 ns pulse width) and the $fast$ electron beam one (500 MHz, 60 ps pulse width). A 100 ms macro-time structure of alternating 80 – 90 ms beam-on and 20 – 10 ms beam-off intervals was used for background monitoring, where Table 1 gives the corresponding 80 or 90 $\%$ fill factors for each measurement range. The NaI:Tl detector (8” diameter $\times$ 12” length) used as a flux monitor was placed in beam, downstream of the neutron detection system. The number of incident $\gamma$ -rays was obtained by applying the pile-up or Poisson fitting method [23] on the NaI:Tl detector response functions, which represented summed spectra for the Poisson distributed $\gamma$ photons in the 20 – 40 ns wide pulses. The mean numbers of $\gamma$ rays per pulse and the incident photon flux values are given in Table 1.

Table 1: Parameters for the laser, electron and LCS

\gamma

-ray beams and measurement conditions.

Z_{R}

is the laser Rayleigh length [24],

\Delta z

is the longitudinal displacement between the focal positions of the laser and electron beams.

f_{\mathrm{e^{-}}}

is the RF frequency of the NewSUBARU storage ring. We give the pulse width of the 198 electron beam bunches circulating in the ring.

E_{m}

is the maximum energy of the

\gamma

-ray beam [22].

100% linear polarization perpendicular to accelerator plane
	$(\gamma,\,n)$	$(\gamma,\,in)$
neutron		multiplicity
counting		sorting
Laser beam:
Laser	Inazuma	Talon
Wavelength (nm)	1064	532
$Z_{R}$ (m)	0.57	6.1
$\Delta z$ (m)	1.8	2.8
Power (W)	$<$ 40	$<$ 20
$f_{\mathrm{laser}}$ (kHz)	20	1
Pulse width (ns)	20	40
Beam on/off fill factor (%)	80	90
Electron beam:
$E_{\mathrm{e^{-}}}$ (MeV)	651.30 – 887.65	649.29 – 1050.63
$f_{\mathrm{e^{-}}}$ & pulse width	500 MHz & 60 ps
Emittance ( $\varepsilon_{x},\varepsilon_{y}$ ) (nm-rad):
– nominal at injection	(38,1-3.8)
– simulated	(50,5) – (70,7)
LCS $\gamma$ -ray beam:
$E_{m}$ (MeV)	7.50 – 13.86	13.84 – 38.02
Mean nb. of $\gamma$ /pulse	5 – 14	6 – 20
Incident flux ( $\gamma$ /s)	(8 – 22) $\times$ $10^{4}$	(0.5 – 1.8) $\times$ $10^{4}$
$\Delta$ E ${}_{\mathrm{FWHM}}$ (MeV)	0.2 – 0.5	0.2 – 1.1
$\Delta$ E ${}_{\mathrm{FWHM}}$ (%)	2.7 – 3.9	1.6 – 3.1
Measurement conditions:
Irradiation time (min)	5 – 30	$\sim$ 120
Target areal
density (g/cm ${}^{2}$ )	4.36	10.97

The spectral distribution of the incident LCS $\gamma$ -ray beam was monitored in between target irradiations with a 3.5” diameter $\times$ 4.0” length LaBr ${}_{3}$ :Ce detector placed in beam. For the incident spectra measurements, the laser was operated in a continuous wave (CW) mode at reduced power. The experimental detector response functions were reproduced by Monte Carlo simulations performed with the LCS $\gamma$ -source simulation code eliLaBr [26, 24, 25]. The simulated energy spectra of the incident beams used in the present experiment are shown in Fig. 2 by the $L(E_{\gamma},E_{m})$ path length weighted energy distributions for each LCS $\gamma$ -ray beam of $E_{m}$ maximum energy [24]. The $L(E_{\gamma},E_{m})$ distributions account for the $\gamma$ -beam self-attenuation in the irradiated sample material and for the secondary radiation generated by electromagnetic interaction of the $\gamma$ -beam with the target, which, for $\sim$ 30–40 MeV photon beams, has sufficiently high energies to induce nuclear reactions in the target.

By investigating the range of suitable values for the electron beam emittance and for small laser-electron beam transverse offsets, we estimated the uncertainty in the incident spectra determination, which is shown in Fig. 2 by the red and green bands for the $\gamma$ beams produced with the 1064 and respectively 532 nm lasers, respectively. Table 1 gives values for the main input parameters used in the LCS $\gamma$ -ray beam simulations and for the energy spread of the best fit $\gamma$ -ray beams.

II.2 Targets

Enriched metal powder of ${}^{208}$ Pb (98.4 $\%$ ${}^{208}$ Pb, 1.2 $\%$ ${}^{207}$ Pb, 0.1 $\%$ ${}^{206}$ Pb, 0.3 $\%$ ${}^{204}$ Pb) was pressed and shaped into two targets of 8 mm diameter and lengths of 4 and 6 mm. The neutron counting measurements below S ${}_{2n}$ were performed with the thinner target of 4.36 g/cm ${}^{2}$ areal density, for which the photon transmission varied between 79.4 and 82.1 $\%$ . The neutron multiplicity sorting measurements above S ${}_{2n}$ were performed with the two targets stacked, resulting in a 10.97 g/cm ${}^{2}$ total areal density. The photon transmission through the stacked targets varied between 45.0 and 55.3 $\%$ . The target alignment to the $\gamma$ -beam axis has been performed using the visible synchrotron radiation passing through the collimators C1 and C2 and has been checked by measurements with an X-ray MiniPIX camera [27, 28], which also confirmed the beam spot estimations of $\sim$ 4 mm diameter [29].

II.3 Neutron detection

The ${}^{208}$ Pb targets were placed in the center of a flat-efficiency moderated neutron detection array of ${}^{3}$ He proportional counters (10 atm., 2.5 cm diameter $\times$ 45 cm length active volume) [19]. The array consisted of three concentric rings of 4, 9, and 18 ${}^{3}$ He counters placed 5.5, 13.0 and 16.0 cm from the $\gamma$ -ray beam axis, respectively (see Figs. 1 and 2 of Ref. [19]). The ${}^{3}$ He tubes were embedded in a high-density polyethylene moderator block having a central axial hole that allowed the $\gamma$ -ray beam to pass through.

Figure 3(a) shows mcnp [30] simulations for the total detection efficiency (thick solid lines) and for the detection efficiencies of the inner ring (thin solid lines) and of the summed two outer rings (dotted lines) compared to an experimental calibration with a ${}^{252}$ Cf source of known activity [19]. Calculations are shown for monochromatic neutrons (black) and for evaporation neutron spectra described by the Weisskopf-Ewing function [31] (red lines) and are represented at the corresponding average neutron energies. We notice the kink in the monochromatic efficiency curve given by the 3 MeV resonance in the cross section of neutron elastic scattering on ${}^{12}$ C. The ${}^{3}$ He array geometrical configuration was optimized so that the total detection efficiency is (i) insensitive to the shape of the neutron spectrum and (ii) constant in the 10 keV to 5 MeV average neutron energy range, with a $\sim$ 5 $\%$ variation between 38 $\%$ and 33 $\%$ for both monochromatic and evaporation neutron spectra.

However, the detection efficiencies of the inner and outer rings do depend on the neutron energy: the inner ring efficiency decreases and, the outer rings efficiency increases with increasing neutron energy. This feature is used through the ring ratio technique to determine the average energy of the neutron emission spectrum [20]. Figure 3(b) shows the ratio of the simulated detection efficiencies of the outer and inner rings for evaporation and monochromatic spectra versus the average neutron energy. We notice that, as the $f_{RR-evap}$ evaporation and $f_{RR-mono}$ monochromatic ring ratio curves diverge starting with average neutron energies of $\sim$ 1 MeV, the ring ratio method is in fact sensitive to the neutron emission spectrum, which has to be properly described in order to extract accurate average neutron energy values.

II.4 Neutron coincidence events

The neutron-multiplicity sorting measurements above the two neutron separation energy $S_{2n}$ involved the recording of neutron coincidence events. By $i$ -fold coincidence neutron events we refer to events where $i$ neutrons were detected during the 1 ms time interval between two consecutive $\gamma$ -beam pulses. Figure 4 shows the histogram of the arrival time for neutrons emitted in the photoneutron reactions on ${}^{208}$ Pb at $E_{\gamma}$ = 27.72 MeV $<$ S ${}_{4n}$ incident photon energy. The neutron counts are discriminated by the (left) inner and (right) outer firing ring and by the coincidence order: (a) single, (b) double and (c) triple fold events. Black histograms show beam-on data and red ones show beam-off data multiplied by 9, the ratio of the measurement time with beam-on and beam-off of 90 ms and 10 ms, respectively. The flat background is due to neutrons generated by the continuous bremsstrahlung emitted by relativistic electrons circulating in the NewSUBARU ring and a small cosmic component. The background subtraction procedure discussed in Refs. [19, 32, 20] relies on fitting the time distribution with a sum of exponentials plus a flat background shown by the green lines in Fig. 4.

III Data analysis

In the following, we discuss the methods for neutron multiplicity sorting, for extraction of information on the photoneutrons mean energy and spectra and on the unfolding of the measured quantities by considering the spectral distribution of the incident photon beams.

III.1 Primary experimental quantities: $i$ -fold neutron cross sections and energies

We will here discuss the numbers of $i$ -fold neutron coincidence events $n_{i}$ and the average energies of neutrons recorded in $i$ -fold coincidences $E_{i}$ , which are the key experimental observables necessary for extracting the photoneutron cross sections $\sigma_{inX}$ and the average neutron energies $E_{inX}$ through neutron multiplicity sorting methods.

III.1.1 i-fold neutron cross sections: $N_{i}$

The number of $i$ -fold neutron coincidence events $n_{i}$ is obtained as

n_{i}(E_{m})=\cfrac{r_{i}^{out}(E_{m})+r_{i}^{in}(E_{m})}{i},

(2)

where $E_{m}$ is the maximum energy for each of the 85 incident LCS $\gamma$ -ray beams used in the present experiment. $r_{i}^{out}$ and $r_{i}^{in}$ are the numbers of neutrons recorded in $i$ -fold coincidence events by the outer two rings and by the inner ring of counters, respectively. They have been obtained by integrating the corresponding background subtracted arrival time histograms.

We further introduce the $i$ -fold neutron cross sections ( $N_{i}$ ) for conveniently representing the $i$ -fold neutron counts on the entire excitation energy range here studied, without depending on the particular target thickness and photon flux values for each experimental point. Thus, $N_{i}$ is defined as:

N_{i}(E_{m})=\cfrac{n_{i}(E_{m})}{N_{\gamma}(E_{m})n_{T}\xi(E_{m})}

(3)

where $n_{T}$ is the concentration of target nuclei and $N_{\gamma}$ is the incident photon number. $\xi=[1-\mathrm{exp}(-\mu L)]/\mu$ is a thick target correction factor given by the target thickness $L$ and attenuation coefficient $\mu$ . Figs. 5(a–d) show the experimental $i$ -fold cross sections (empty black dots) for the ${}^{208}$ Pb $(\gamma,\,inX)$ reactions with $i$ ranging from (a) 1 to (d) 4. The non-zero $N_{3}$ and $N_{4}$ values below the three and four neutron separation energies, respectively, show the presence of pile-ups, or multiple-firing neutron events which will be discussed in Sect. III.2.

III.1.2 $i$ -fold average neutron energies above S ${}_{2n}$ : $E_{i}$

Information on the neutron spectra emitted in photoneutron reactions on ${}^{208}$ Pb have been obtained based on the experimental $RR_{i}$ ring ratios for neutrons recorded in $i$ -fold coincidence events, which are defined as:

RR_{i}(E_{m})=\cfrac{r_{i}^{out}(E_{m})}{r_{i}^{in}(E_{m})}.

(4)

For incident energies above $S_{2n}$ , the average energy of neutrons recorded in $i$ -fold coincidence events has been directly determined by evaluating the $f_{RR-evap}^{-1}$ neutron evaporation ring ratio function at the experimental $RR_{i}(E_{m})$ ring ratio value (see Fig. 3)

E_{i}^{exp}=f_{RR-evap}^{-1}(RR_{i})\mathrm{\quad for\quad}E_{m}>S_{2n}.

(5)

Figure 5 shows the ring-ratio extracted average energies $E_{i}$ of neutrons recorded in $i$ -fold events for i = 1 (e) to i = 4 (h).

III.1.3 $(\gamma,\,n)$ photoneutron energies below S ${}_{2n}$

Figure 6 shows a partial decay scheme of ${}^{208}$ Pb illustrating the photoexcitation and decay of ${}^{208}$ Pb to states in ${}^{207}$ Pb by 1n-emission. E1 excitations populate $1^{-}$ states in ${}^{208}$ Pb with probabilities given by the photoabsorption cross section. The ${}^{207}$ Pb levels are well separated, with the level spacing much larger than the level width, up to excitation energies higher than 4 MeV. Based on the analysis of the cumulative number of levels, the RIPL3 recommendation [33] is to consider a complete level scheme of 66 discrete excited levels and start the continuum at 4.25 MeV excitation energy.

For incident photon energies below $S_{n}$ + 0.57 MeV = 7.94 MeV, only monochromatic neutrons can be emitted, populating the ground state of ${}^{207}$ Pb. Such is the case for the three lowest incident energy settings at 7.50, 7.66 and 7.80 MeV, for which the experimental neutron energies have been directly determined by evaluating the $f_{RR-mono}^{-1}$ function at the measured ring ratio values $RR_{1}$ :

E_{n}^{exp}=f_{RR-mono}^{-1}(RR_{1})\mathrm{\quad for\quad}E_{m}<\mathrm{7.94~% {}MeV}.

(6)

Figures 7(a) and (b) compare the experimental ring ratios and neutron energies (red points), respectively, with the values given by the kinematics of the two-body breakup ${}^{208}$ Pb $\rightarrow$ n + ${}^{207}$ Pb, $E_{n}=(E_{\gamma}-S_{n})\cdot 207/208$ . We notice that the experimental results are consistent, within error bars, with the two-body breakup kinematics.

Starting at photon beam energies larger than 7.94 MeV, neutrons can be emitted also to excited states in ${}^{207}$ Pb. With the incident photon energy increase, the photoneutron decay transitions from emission of discrete energy neutrons to statistical emission of neutrons with continuous spectrum. In order to (i) correct the ring ratio extracted average neutron energies by accounting for the differences between the actual neutron emission spectra and ideal monochromatic or evaporation ones and (ii) obtain indirect information on the neutron spectra, we performed an iterative search for ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron spectra for which the $RR_{1}^{calc}$ calculated ring ratio values reproduce the experimental ones.

Assuming the $Y(E,E_{m})$ spectrum of photoneutrons emitted at an incident energy $E_{m}$ , the calculated ring ratio is obtained as:

RR_{1}^{calc}(E_{m})=\cfrac{\sum_{E}Y(E,E_{m})\cdot\varepsilon^{out}_{mono}(E)% }{\sum_{E}Y(E,E_{m})\cdot\varepsilon^{in}_{mono}(E)},

(7)

where $\varepsilon^{in}_{mono}(E)$ and $\varepsilon^{out}_{mono}(E)$ are the inner ring and respectively, the outer rings detection efficiencies for monochromatic emission of $E$ energy neutrons. We used the CERN minuit package to perform a $\chi^{2}$ minimization procedure and determine the $Y(E,E_{m})$ neutron spectrum for which the calculated ring ratio reproduces the experimental one. Each bin content was considered as a free parameter, where the spectra binning varied between 33 keV and 100 keV, depending on the incident photon energy. The minimized quantity $\chi_{E_{m}}^{2}$ was defined as:

\chi_{E_{m}}^{2}=\cfrac{(RR_{1}^{calc}(E_{m})-RR_{1}(E_{m}))^{2}}{\sigma_{RR_{% 1}}^{2}(E_{m})}.

(8)

and has been calculated for each incident photon energy $E_{m}$ between 7.96 and 13.86 MeV. Here $\sigma_{RR_{1}}(E_{m})$ is the statistical uncertainty for the experimental ring ratio.

The minimization procedure has been applied considering 32 sets of starting values for the neutron spectra calculated using the EMPIRE reaction code [4] with combinations of different nuclear inputs:

•

$\gamma$ -ray strength function: Standard Lorentzian, Modified Lorentzian (1, 2, and 3), Enhanced Generalized Lorentzian and Simplified Modified Lorentzian models with RIPL-3 parameters [33];
•

spherical optical model potentials (OMP) for neutrons of Koning-Delaroche (RIPL IDs 2405 and 1467) and of Weisel (RIPL ID 121) and the coupled channels one of Vonach (RIPL ID 2) [33];
•

level densities: enhanced generalized superfluid model [4], Gilbert-Cameron and microscopic HFB level densities [33] with the $\tilde{a}$ parameter for the ${}^{207}$ Pb level density varied in an interval of 80 $\%$ to 100 $\%$ of the RIPL-3 recommended value,

and discrete levels retrieved from the RIPL-3 database. Sect. V gives more details on the EMPIRE calculations.

Figure 8(a) shows the EMPIRE calculations for ${}^{208}$ Pb $(\gamma,\,n)$ neutron spectra at 12.30 MeV incident photon energy. For all sets of input parameters in the statistical model calculations, we notice the typical neutron spectrum, with a continuous component due to statistical evaporation of low energy neutrons, and discrete high energy neutron transitions due to population of low lying discrete states in the residual. The best fit neutron spectra shown in Fig. 8(b) show that the minimization procedure systematically decreased the continuous component of low energy neutrons, thus increasing the contributions of intermediate energy neutrons of $\sim$ 2 MeV and of the high energy neutrons populating the low discrete states in ${}^{207}$ Pb.

Figure 7(a) shows the experimental ring ratios (black dots) and the ones calculated for the EMPIRE neutron spectra used as starting values in the minimization (solid lines), where each line corresponds to the one of the 32 sets of nuclear inputs. We notice significant differences between the experimental ratios and the ones computed for the starting values spectra. However, the minimization procedure identified $Y_{k}(E,E_{m})$ best-fit photoneutron spectra for which the calculated $RR_{1,k}^{calc}$ ring ratios (empty dots) reproduce well the experimental values, where $k$ indicates the nuclear input set and cycles up to 32.

The $E_{1,k}(E_{m})$ average neutron energies were obtained as the first moments of the $Y_{k}(E,E_{m})$ best-fit photoneutron spectra:

E_{1,k}(E_{m})=\cfrac{\sum_{E}E\cdot Y_{k}(E,E_{m})}{\sum_{E}Y_{k}(E,E_{m})}.

(9)

For each $k$ set of starting values, Fig. 7(b) shows the $E_{1,k}(E_{m})$ average neutron energies (solid lines) as function of the incident photon energy $E_{m}$ , with each line corresponding to one of the 32 sets of minimization results. Good agreement was obtained between the average neutron energy results for each incident photon energy value.

The final $E_{1}$ average energies of ${}^{208}$ Pb $(\gamma,\,n)$ photoneutrons for excitation energies below $S_{2n}$ are shown by the full black dots in Fig. 7(b) and have been obtained as the average of the $E_{1,k}(E_{m})$ results for each incident photon energy. The error bars account for the mean standard deviation of the set of minimization results for each incident energy, the statistical uncertainty in the neutron counts for each ring of counter and a 5 $\%$ uncertainty in the ring ratio average neutron energy determination.

Figure 7(b) also compares the present $E_{1}$ average neutron energies to the ones directly extracted through the ring ratio method assuming either evaporation (empty black dots) or monochromatic (full red dots) neutron emission spectra. For incident photon energies below 8.5 MeV there is good agreement between all three curves, which follows from the overlap between the $f_{RR-evap}$ and $f_{RR-mono}$ ring ratio curves for low neutron energies below $\sim$ 1 MeV. At incident photon energies between 8.5 and 10 MeV, the best-fit $E_{1}$ values reproduce well the monochromatic ones and are well below the evaporation assumptions. For photon energies above 10 MeV, the best-fit $E_{1}$ results follow closely the evaporation curve, suggesting a transition to a statistical neutron emission. We notice that the ring ratio minimization procedure applied only small corrections of up to $\sim$ 10 $\%$ to the ${}^{208}$ Pb $(\gamma,\,n)$ average neutron energies in the transition excitation energy range between 7.94 MeV and $S_{2n}$ .

The $Y(E,E_{m})$ final estimations for the ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron spectra at $E_{m}$ incident photon energy have been obtained by averaging each neutron energy bin content of the set of 32 best-fit $Y_{k}(E,E_{m})$ results. This can be observed in Fig. 8(b) for 12.30 MeV incident photon energy. The resulting $Y(E,E_{m})$ neutron spectra are shown and discussed in Sect. IV along with the partial cross sections for the ${}^{208}$ Pb $(\gamma,\,n)$ reactions that populate the residual nucleus ${}^{207}$ Pb in its ground state and in the first and second excited states.

III.2 Neutron multiplicity sorting

Let us consider a photon beam of energy $E_{\gamma}$ incident on the ${}^{208}$ Pb target, where $S_{\texttt{N}n}<E_{\gamma}<S_{(\texttt{N}+1)n}$ and $S_{\texttt{N}n}$ is the separation energy of N neutrons, which are given in Table 2 for neutron multiplicities up to 5. Photoneutron $(\gamma,\,inX)$ reactions with $i$ cycling up to N will be induced in the target, each characterized by the $\sigma_{inX}$ cross section and by the $E_{inX}$ average energy of the photoneutron spectrum. For (i) $single$ -firing conditions in which no more than one reaction is induced in the target by the same photon pulse, and (ii) constant neutron detection efficiency, the Direct Neutron Multiplicity (DNM) sorting method relates the $\sigma_{inX}$ cross sections to the raw $N_{i}$ neutron coincidence experimental observables as [19]:

N_{i}=\sum_{x=i}^{\texttt{N}}\sigma_{xnX}\cdot{}_{x}C_{i}\varepsilon^{i}(1-% \varepsilon)^{x-i},

(10)

and the $E_{inX}$ photoneutron energies to the $E_{i}$ average energies of neutrons detected in $i$ -fold coincidences as [20]:

E_{i}=\sum_{x=i}^{\texttt{N}}E_{xnX}\sigma_{xnX}\cdot{}_{x}C_{i}\varepsilon^{i% }(1-\varepsilon)^{x-i}/N_{i},

(11)

where $\varepsilon$ is the constant neutron detection efficiency and ${}_{x}C_{i}$ is the binomial coefficient.

Table 2:

{}^{208}

Pb separation energies for

i

neutrons

S_{in}

in MeV.

S ${}_{n}$	S ${}_{2n}$	S ${}_{3n}$	S ${}_{4n}$	S ${}_{5n}$
7.368	14.107	22.194	28.927	37.323

For $single$ -firing conditions to be met, low reaction rates were maintained by conveniently limiting the target thickness and incident photon flux. Still, as the number of photons in the LCS $\gamma$ -ray pulses followed Poisson distributions with mean values between 6 and 20 photons per pulse, the incidence of $multiple$ -firings has only been minimized, but not completely eliminated. Multiple firings are visible in Figs. 5(c) and (d) from the non-zero $N_{3}$ and $N_{4}$ cross sections, respectively, below the three neutrons separation energy $S_{3n}$ . Explicitly, the 3-neutron coincidence events below $S_{3n}$ originate from $(\gamma,\,n)$ and $(\gamma,\,2n)$ reactions induced by the same $\gamma$ -ray pulse.

Thus, we applied the statistical treatment of neutron coincidence events described in Ref. [20], which models the $multiple$ -firing of all available combinations of photoneutron $(\gamma,\,inX)$ reactions with corresponding contributions given by the reaction cross sections, areal density of the irradiated samples and the multiplicity of incident photons per pulse. Using the CERN minuit package, a $\chi^{2}$ minimization procedure has been performed to determine the $multiple$ -firing corrected ${}^{208}$ Pb cross sections and average photoneutron energies from the measured $i$ -fold cross sections $N_{i}$ and average energies $E_{i}$ . The ${}^{208}$ Pb $(\gamma,\,inX)$ cross sections and average neutron energies are the free parameters in the minimization procedure, while the target ( $n_{T}$ , $\xi$ ) and LCS $\gamma$ -ray beam ( $N_{\gamma}$ , mean number of photons per pulse) characteristics have been fixed to the experimentally determined values. The best-fit $i$ -fold cross sections $N_{i}^{M\!F}$ and average neutron energies $E_{i}^{M\!F}$ are shown in Fig. 5 in comparison with the raw measured values. We notice that the raw measured data are well reproduced by the $multiple$ -firing results on the entire excitation energy region and for all observed neutron multiplicities.

Figure 9 shows the present neutron multiplicity sorting results for the ${}^{208}$ Pb $(\gamma,\,inX)$ reactions with (a) $i$ = 1, (b) $i$ = 2, (c) $i$ = 3 and (d) $i$ = 4. The DNM cross sections $\sigma_{inX}^{D\!N\!M}$ and average photoneutron energies $E_{inX}^{D\!N\!M}$ obtained by solving the systems of equations (10), and (11), respectively, are shown by the blue crosses. The $multiple$ -firing (MF) corrected ones $\sigma_{inX}^{M\!F}$ and $E_{inX}^{M\!F}$ are shown by the empty red dots and correspond to the best-fit $N_{i}^{M\!F}$ and $E_{i}^{M\!F}$ shown in Fig. 5.

$\sigma_{inX}$ results.

The reasonable fulfillment of the $single$ firing conditions is indicated by the increasing behavior of the $\sigma^{D\!N\!M}_{inX}$ with i = 2, 3, 4 starting from the reactions threshold and their drop with the opening of the channel with the immediately higher neutron emission multiplicity. Still, $multiple$ -firing corrections have produced visible changes to the cross sections, especially in the excitation energy regions just above $S_{in}$ reaction thresholds, bringing a general improvement in the scatter of the cross sections for all channels. The $(\gamma,\,2n)$ cross sections in the vicinity of $S_{2n}$ have been decreased by taking into account the multiple firings of the strong $(\gamma,\,1nX)$ channel. The non-zero $\sigma^{D\!N\!M}_{3nX}$ and $\sigma^{D\!N\!M}_{4nX}$ values below $S_{3n}$ were reallocated to the $(\gamma,\,1nX)$ and $(\gamma,\,2nX)$ channels.

$E_{inX}$ results.

The $multiple$ -firing corrections applied to the $E_{inX}^{D\!N\!M}$ energies are generally small and within error bars. Exceptions are seen in the incident energy regions just above the $S_{2n}$ and $S_{3n}$ thresholds, where the weak newly opened $(\gamma,\,2n)$ and $(\gamma,\,3n)$ channels have small photoneutron energies compared with the stronger channels of lower neutron emission multiplicities. Explicitly, the $multiple$ firings of $(\gamma,\,1nX)$ reactions with high $E_{1nX}$ values artificially increase the average energy of double neutron events to values higher than the maximum energies given by kinematics $207/208\cdot(E_{\gamma}-S_{2n})$ . We notice that the average neutron energies show strong statistical fluctuations at high excitation energies, which confirms the need of a flat efficiency neutron detector for precise reaction cross section measurements.

III.3 Energy unfolding

The $(\gamma,\,inX)$ cross sections and average photoneutron energies obtained through the neutron multiplicity sorting method discussed above are referred to as monochromatic approximations representing in fact the folding of the true, energy-dependent quantities with the spectral distribution of the incident photon beams. The measured cross sections are given by the folding of the excitation function $\sigma_{inX}(E_{\gamma})$ with the beam spectral distribution $L(E_{\gamma},E_{m})$ :

\sigma_{inX}^{M\!F}(E_{m})=\cfrac{1}{\xi}\int_{0}^{E_{m}}L(E_{\gamma},E_{m})% \sigma_{inX}(E_{\gamma})\,dE_{\gamma},

(12)

where $\xi$ is the thick target correction factor introduced in Sec. III.1.1. The measured average energies are obtained by folding the incident energy-dependent function $E_{inX}(E_{\gamma})$ with the beam spectral distribution and the excitation function:

E_{inX}^{M\!F}(E_{m})=\cfrac{\int_{0}^{E_{m}}\!\!E_{inX}(E_{\gamma})L(E_{% \gamma},\!E_{m})\sigma_{inX}(E_{\gamma})dE_{\gamma}}{\sigma_{inX}^{M\!F}(E_{m}% )\xi}.

(13)

The $\sigma_{inX}^{M\!F}(E_{m})$ and $E_{inX}^{M\!F}(E_{m})$ monochromatic approximations are not connected to a specific incident photon energy $E_{\gamma}$ , but are functions of the LCS $\gamma$ -beams maximum energy $E_{m}$ .

An iterative unfolding procedure described in Refs. [34, 35] has been applied to the quantities defined in Eqs. (12) and (13) in order to extract the energy-dependent photoneutron cross sections and average energies. The method was applied independently for each quantity $\sigma_{inX}(E_{\gamma})$ and $E_{inX}(E_{\gamma})$ with $i$ = 1 to 4, first for the cross sections and then for the average neutron energies.

For each unfolding process, the method starts with a constant trial function of 10 keV binning spanning from the corresponding $S_{in}$ threshold up to the maximum investigated energy of 38.02 MeV. The trial function is iteratively adjusted based on the difference between its folding with the beam spectral distribution and the measured quantity. The iterations stop when the folded trial function reproduces the measured cross sections (average neutron energies). The adjustment of the trial function has been performed independently for each energy bin using linear interpolations between the measured/folded values as functions of $E_{m}$ . To prevent spurious fluctuations being introduced by the unfolding method, after each iteration we applied an energy-dependent (0.1 – 1.4) MeV smoothing factor related to the FWHM energy distribution of the LCS $\gamma$ -ray beams listed in Table 1.

Particular attention has been paid to the unfolding of the resonant $\sigma_{\gamma,\,n}$ cross sections in the energy excitation region below $S_{2n}$ . Here, we found that smoothing factors smaller than the experimental FHWM spectral distributions were needed in order to reproduce the rapid changes in the monochromatic approximation of the $(\gamma,\,n)$ cross section. To validate the results, we performed a separate unfolding process in which we described the photoabsorption cross section as a sum of analytical functions: two Lorentz functions for the entire GDR region and eight Gauss functions for describing the resonant structure below $E_{\gamma}$ = 13 MeV. Using the CERN minuit package, we searched for an excitation function that, folded with the $\gamma$ beam spectral distributions, optimally reproduced the monochromatic $(\gamma,\,n)$ cross sections below $S_{2n}$ . The strength, centroid, and width of the Lorentz and Gauss functions were considered as free parameters in the minimization. The resulting $\sigma_{\gamma,\,n}$ , shown by the empty green triangles in Fig. 9(a1), reproduced within error bars the cross sections obtained through the iterative method of Ref. [34]. Finally, we adopted the results obtained by the iterative method of Ref. [34], which makes no assumptions on the shape of the unfolded cross sections.

Figure 9 shows the energy unfolded cross sections and average neutron energies (full black dots) along with the monochromatic approximation values. The spectral distribution of the incident $\gamma$ -ray beams has already been considered in the procedure applied for extracting the $E_{1}$ average energy of $(\gamma,\,n)$ neutrons, and thus the unfolded values (black dots) coincidence with the ones obtained in Sect. III.1.3 (red empty dots). The error bars for the energy unfolded results account for the statistical uncertainties in the neutron detection and for uncertainties of 3 $\%$ for the neutron detection efficiency [19, 32], 3 $\%$ for the photon flux determination, 1 $\%$ for the target thickness and the incident photon spectra uncertainty. The uncertainty in the incident photon spectra, shown by the red and green bands in Fig. 2, has been propagated by applying the unfolding procedure separately for the upper and lower limit of the incident spectra.

IV Experimental results

IV.1 Photoabsorption cross sections

Based on the negligible contributions of $(\gamma,\,p)$ , $(\gamma,\,\alpha)$ , etc. charged particle only emission reactions, which are highly suppressed by the large Coulomb barrier, the total photoabsorption cross section in ${}^{208}$ Pb is well approximated by the sum cross section for the neutron emission reactions:

	$\displaystyle\sigma_{abs}$	$\displaystyle\equiv\sigma_{\gamma,\,abs}$
		$\displaystyle\approx\sigma_{1nX}+\sigma_{2nX}+\sigma_{3nX}+\sigma_{4nX}$		(14)

The so obtained ${}^{208}$ Pb photoabsorption cross sections are shown in Fig. 10 in comparison with preceding data. The present results are in good agreement with the Saclay data of Veyssiere $et$ $al.$ [9] and with the indirectly determined photoabsorption cross sections extracted from RCNP proton inelastic scattering experiments by Tamii $et$ $al.$ [10]. The Livermore cross sections obtained by Harvey $et$ $al.$ [7] for the entire GDR energy range strongly underestimate all four other sets of results shown in Fig. 10. The cross sections remeasured by Berman $et$ $al.$ [8] at Livermore using a natural lead sample are only slightly lower than the present results, by $\sim$ 5 $\%$ for the GDR peak value. We note instead that the centroid and width of the giant dipole resonance is in good agreement for all experimental data sets.

IV.2 Photoneutron cross sections

Figure 11(a) shows the present $\sigma_{1nX}$ photoneutron cross sections for ${}^{208}$ Pb compared with previous measurements. As for the photoabsorption, we notice a good agreement with the Saclay data of Veyssiere $et$ $al.$ [9], a strong underestimation of the present results by the Livermore data of Harvey $et$ $al.$ [7] and $\sim$ 5 $\%$ higher present GDR peak cross section values than the Berman $et$ $al.$ [8] data. The present $\sigma_{1nX}$ are also in good agreement the recent measurements performed by Kondo $et$ $al.$ [38] (green empty triangles) using LCS $\gamma$ -ray beams with energies up to 13 MeV and a moderated neutron detection array of ${}^{3}$ He counters. Good agreement with bremmstrahlung monochromator data is observed for the unpublished results of Young [39] (full green dots) in the GDR peak region and with those of Calarco [40] only on the increasing slope of the $\sigma_{1nX}$ . The $\sigma_{1nX}$ cross sections of Alarcon $et$ $al.$ [41] (magenta empty diamonds), which were obtained by integrating differential cross sections over neutron energy and angles, show GDR values lower than the present data.

The present ${}^{208}$ Pb $(\gamma,\,2nX)$ cross sections shown in Fig. 11(b) are systematically higher than the Veyssiere $et$ $al.$ [9] and Berman $et$ $al.$ [8] ones. A better agreement is found with the Harvey $et$ $al.$ [7] results below the $\sim$ 18 MeV. Above $\sim$ 18 MeV, it is difficult to make a comparison because of the large statistical fluctuations of the Livermore results.

Figure 11(c) shows the present ${}^{208}$ Pb $(\gamma,\,3nX)$ and ${}^{208}$ Pb $(\gamma,\,4nX)$ cross sections, which are in good agreement with the Saclay results in the excitation energy region up to 38 MeV.

The present photoneutron cross sections are available in numerical format in Ref. [37]. The EMPIRE and TALYS calculations shown in Fig. 11 by the blue and green lines, respectively, are discussed in Sect. V.

IV.3 Resonant structures in low-energy ${}^{208}$ Pb $(\gamma,\,n)$ cross sections

Figure 12 shows the ${}^{208}$ Pb $(\gamma,\,n)$ cross sections in the low excitation region below the $S_{2n}$ threshold. The present cross sections are plotted along with existing high resolution measurements in order to compare the results on the resonant structures below $E_{\gamma}$ = 13 MeV. The peaks at 7.65, 8.10 and 8.45 MeV reproduce well the structures observed in the Saclay experiment [9]. The broad shoulder at 9.45 MeV has also been observed both in the Saclay and RCNP [10] experiments and in the tagged bremmstrahlung measurements of Young [39] and Alarcon $et$ $al.$ [41]. For the prominent $\sim$ 10 MeV peak we obtained a 10.06 MeV centroid value, which is in agreement with the tagged bremmstrahlung experiments of Bell $et$ $al.$ [42] and Alarcon, but about $\sim$ 150 keV higher than the peak position obtained in the Saclay experiment. The 10.70 MeV peak has been previously observed only in the Bell measurements, while the Saclay and RNCP experiments described it as a low shoulder. The present 11.31 MeV peak position confirms the tagged bremmstrahlung results of Alarcon and Calarco [40], but is again $\sim$ 100 keV higher than the Saclay and RCNP centroid values. The 12.20 MeV shoulder is also observed in all the other measurements.

IV.4 Average photoneutron energies

Figure 13 shows the average energies of photoneutron energies emitted in the (a) $(\gamma,\,1nX)$ , (b) $(\gamma,\,2nX)$ , (c) $(\gamma,\,3nX)$ and $(\gamma,\,4nX)$ reactions in ${}^{208}$ Pb.

The energy of the $(\gamma,\,1nX)$ photoneutrons shows a steep increase starting from $S_{n}$ up to 11.5 MeV excitation energy followed by a fast drop until the opening of the $(\gamma,\,2n)$ channel. Above $S_{2n}$ , the average neutron energy resumes its increase, which confirms the unpublished results of Mutchler [43] obtained by the photon difference method using bremsstrahlung photon beams of 13, 14 and 15 MeV and time of flight measurements at 3 angles relative to the photon beam. We notice that the average neutron energies steeply increase around 22 MeV excitation energy.

The $(\gamma,\,2nX)$ reaction is characterized by a long rise in the average photoneutron energy $E_{2nX}$ , which increases quickly above $S_{2n}$ . The growth slows down at $\sim$ 18 MeV excitation energy and continues until the $S_{4n}$ threshold, above which the experimental $E_{2nX}$ shows a slow decrease with strong fluctuations which follow from the statistical scatter in the unfolded average neutron energy values shown in Fig. 9(b2). The $E_{3nX}$ average neutron energies, starting from $S_{3n}$ , increase up to $\sim$ 5 MeV above $S_{4n}$ , where they start decreasing.The $(\gamma,\,4nX)$ reaction shows a rather abrupt increase in the average neutron energy values, starting $\sim$ 4 MeV above the reaction threshold.

The present photoneutron average energies are available in numerical format in Ref. [37]. The comparison with reaction model codes is discussed in Sect. V.

IV.5 Estimations for $(\gamma,\,n)$ photoneutron spectra and partial cross sections

IV.5.1 ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron energy spectra

Figure 14 shows the present estimations for the ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron spectra $Y(E,E_{m})$ obtained by reproducing the experimental ring ratio values. The spectra represent branching ratios normalized so that their integral is equal to 1000. The error bar for each energy bin content is the quadratic sum of the ring ratio statistical error and of the standard deviation of the mean over $Y_{k}(E,E_{m})$ , with $k$ indicating the set of nuclear inputs used for the minimization starting values.

The neutron emission to the ground (0 keV, 1/2 ${}^{-}$ ), first (569.7 keV, 5/2 ${}^{-}$ ) and second (897.7 keV, 3/2 ${}^{-}$ ) excited states, which are recognized by the three highest energy discrete neutron emissions, show significant yields up to incident photon energies of $\sim$ 13 MeV. Starting with incident photon energies above 10 MeV, we notice significant population yields for the higher excited states. This is in agreement with the behavior of the average energy of the total ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron spectra, which, as shown in Fig. 7(b) follow closely the evaporation ring ratio average energy curve for excitation energies higher than 10 MeV. We also notice that the present estimations obtained by reproducing the experimental ring ratio values confirm the transition to a statistical neutron emission towards the $(\gamma,\,2n)$ reaction threshold.

IV.5.2 Partial photoneutron cross sections for $E_{\gamma}<S_{2n}$

Figure 15 shows the partial cross sections for ${}^{208}$ Pb $(\gamma,\,n)$ photoneutron reactions populating the residual nucleus ${}^{207}$ Pb in its (a) ground state, (b) first and (c) second excited states and (d) in any excited state higher than the second. These were obtained by multiplying the branching ratio corresponding to the respective state with the total photoneutron cross section.

The present estimations for the partial cross section show reasonably good agreement with the results of Bell $et$ $al.$ [42] (full red dots) obtained at the Illinois bremmstrahlung monochromator using time of flight measurements at seven angles relative to the photon beam. The present $\sigma(\gamma,\,n_{0})$ and $\sigma(\gamma,\,n_{1})$ cross sections for leaving the residual ${}^{207}$ Pb in its ground state and in its first excited state, respectively, reproduce the narrow peak at 10 MeV excitation energy and the structure at $\sim$ 10.5 MeV. However, the present $\sigma(\gamma,\,n_{2})$ cross section for populating ${}^{207}$ Pb in its second excited state doesn’t reproduce the strong structure observed by Bell at 10.5 – 11 MeV. The present $\sigma(\gamma,\,n_{>})$ for leaving ${}^{207}$ Pb in an excited state higher then the second is in agreement with the increasing behavior of the Bell data, however with systematically higher absolute values. This follows from the low energy neutron detection threshold of $\sim$ 600 keV from the Bell experiment, in which low energy neutrons emitted to high excited states were not recorded.

The present $\sigma(\gamma,\,n_{0})$ estimations are further compared with the high resolution measurements of Sherman et al. [44], which were performed using bremsstrahlung beams and the time of flight technique with a single neutron detector placed at 90 ${}^{\circ}$ to the photon beam. The approximate angle integrated cross sections for populating the ground state of ${}^{207}$ Pb shown by the orange histogram in Fig. 15(a) were obtained by assuming $W(\theta)=2-P_{2}(\cos\theta)$ as the neutron angular distribution. Overall good agreement is found between the present $\sigma(\gamma,\,n_{0})$ estimations and the data of Sherman with the exception of the 10.5 – 11.5 MeV excitation energy region, where the present results are systematically higher than the ones of Sherman.

The partial cross sections of Calarco [40] shown by the green crosses in Fig. 15 (a), (b) and (c) are also approximations. They have been obtained by multiplying by 4 $\pi$ the differential cross section measured at 115 ${}^{\circ}$ , in the assumption of isotropic neutron emission. Considering only dipole photoexcitations, s-wave neutron emission is possible from the 1 ${}^{-}$ states in ${}^{208}$ Pb to the 1/2 ${}^{-}$ ground state and to the 3/2 ${}^{-}$ second excited state in ${}^{207}$ Pb. Indeed, there is remarkably good agreement between the present $\sigma(\gamma,\,n_{2})$ cross sections and the ones of Calarco. However, the present $\sigma(\gamma,\,n_{0})$ cross sections are systematically higher than the ones of Calarco. For the $\sigma(\gamma,\,n_{1})$ cross section, the Calarco results confirm the $\sim$ 11.25 MeV peak position also observed in the present estimations but do not show the structures at 10.5 – 11 MeV present in the Bell measurements and in the present results.

V Statistical model calculations

The new experimental data are now compared with statistical model calculations obtained by the EMPIRE [4] and TALYS [3] codes. Since the entrance channel plays a fundamental role for an accurate description of the various reaction channels, the fit to experimental photoabsorption cross sections has been tested with several Lorentzian-type closed-forms (SLO, MLO1, SMLO) plus the quasi-deuteron contribution for the E1 $\gamma$ -ray strength functions [33, 36]. We found that the SMLO model, properly tuned, reproduces well, though not perfectly, the experimental data. The resulting photon strength function is obtained by the sum of three SMLO-type Lorentzians for which the parameters are given in Table 3. In addition, the quasideuteron contribution to the photo-absorption cross section is included with a normalisation factor of 60% to ensure a proper description of the high-energy tail of the cross section at the highest energies considered here. The resulting SMLO fit is shown in Fig. 10 to be rather accurate though some fine structures (e.g. around 20-24 MeV) could not be described.

The specific optical model potential (OMP) for ${}^{206-208}$ Pb (RIPL ID - 102) [33, 45] is used in the TALYS code to obtain the transmission coefficients for neutron emission. Since the EMPIRE code can only consider a single OMP common to all nuclei in the reaction chain, the Koning-Delaroche general spherical OMP (RIPL ID - 2405) is used in the EMPIRE calculations. The latter was also used for light Pb isotopes within the TALYS code. The level densities at the equilibrium deformationare described with the enhanced generalized superfluid model [4] in the EMPIRE code and with the temperature-dependent Hartree-Fock-Bogolyubov (HFB) plus combinatorial model [46] in the TALYS code. In the latter code, the first 40 discrete excited levels coming from experiment [33] are adopted for all target and residual nuclei.

In both approaches, the level densities have been adjusted to reproduce the low-lying discrete level scheme and the $s$ -wave resonance spacing at the neutron separation energy. In order to reproduce the present experimental results, we found it necessary in the EMPIRE calculations to reduce the so determined ${}^{207}$ Pb level density by lowering the $\tilde{a}$ parameter by 20%. This increased the $(\gamma,\,1nX)$ cross section tail at the expense of the newly opened $(\gamma,\,2n)$ channel and lowered the $E_{1nX}$ average neutron energies. In the EMPIRE calculations, we also lowered the starting energy for the continuum in ${}^{207}$ Pb from the 4.25 MeV RIPL-3 recommendation to 3 MeV. This improved the description of the $(\gamma,\,n)$ average neutron energies at 10 to 12 MeV excitation energy.

The high photon energies considered in the present experiment also induce an important contribution from the pre-equilibrium cross section. In the TALYS code, it is given by the default two-component exciton model [3] where only the single-particle state densities of the Pb isotopes have been adjusted to improve the description of the present experimental cross sections. We found that a reduction of about 30% is needed, except for ${}^{207}$ Pb for which the reduction factor is about 80%. Such a high reduction factor can be attributed to the strong shell effect in this doubly magic region which is known to strongly affect the level density parameter.

Table 3:

{}^{208}

Pb GDR parameters adopted within the SMLO model for the three Lorentzians.

\sigma

corresponds to the peak cross section,

E

to the centroid energy and

\Gamma

to the full width at half maximum.

i

is the Lorentzian index.

$i$	$\sigma_{\rm GDR}^{(i)}$ (mb)	$E_{\rm GDR}^{(i)}$ [MeV]	$\Gamma_{\rm GDR}^{(i)}$ [MeV]
1	260.65	12.20	3.251
2	526.13	13.93	3.06
3	8.44	25.56	1.97

The new experimental photoneutron cross sections are compared in Fig. 11 with EMPIRE and TALYS calculations. The ${}^{208}$ Pb $(\gamma,1nX)$ cross section is rather well described by both calculations, though, as for the photoabsorption cross section, some detailed pattern could not be reproduced with the SMLO input photon strength function. Both the EMPIRE and TALYS calculations tend to overestimate the ${}^{208}$ Pb $(\gamma,2nX)$ cross section in the 15 to 20 MeV incident energy region, but a proper description of the energy dependence in the vicinity of 25 MeV is ensured by the inclusion of the third Lorentzian in the SMLO photon strength function (see Table 3). While the ${}^{208}$ Pb $(\gamma,3nX)$ cross section can be correctly reproduced by both reaction codes, the $(\gamma,4nX)$ channel is underestimated by TALYS but correctly described by EMPIRE.

In a similar way, the newly measured average photoneutron energies are compared in Fig. 13 with EMPIRE and TALYS calculations. Here, we plot the average energies of the exclusive neutron spectra obtained in the EMPIRE and TALYS calculations, which contain the summed contributions of reactions with emission of a given $i$ neutron multiplicity. The EMPIRE and TALYS model calculations describe qualitatively well the $E_{inX}$ behavior at the lowest incident photon energies. However, the calculations could not reproduce the detailed patterns observed in the energy dependence of $E_{inX}$ , e.g. the position of the inflection point at $\sim$ 11 MeV in the $E_{1nX}$ evolution, suggesting that a better knowledge of the discrete levels in ${}^{207}$ Pb is required to be able to describe the high average neutron energies in the 10 to 12 MeV excitation energy region. Above 25 MeV, the $E_{1nX}$ average neutron energy shows a long decrease, a tendency qualitatively reproduced by the EMPIRE calculations, which however overestimate the $E_{1nX}$ values between 17 and 35 MeV excitation energy and cannot describe the structures present at 16, 20 and 29 MeV. At these high energies, above their corresponding thresholds, EMPIRE also predicts much larger photoneutron energies compared with TALYS, for each of these $E_{inX}$ channels. EMPIRE calculations follow closely the increasing $E_{2nX}$ and $E_{3nX}$ experimental values, but fail to reproduce their average energy drop at high incident photon energies above $S_{4n}$ . A decrease in the $E_{2nX}$ is predicted in the EMPIRE calculations at incident energies of 35 MeV, but this follows a continuous increase in average neutron energies up to 6 MeV, well above experimental values. For the $(\gamma,1nX)$ and $(\gamma,2nX)$ reactions, the drop in mean neutron energies predicted by EMPIRE calculations is due to the contribution of charged particle emission channels. This can be seen from the continuous increasing trend of the average energies of neutrons emitted in the reactions without emission of charged particles $(\gamma,\,n)$ and $(\gamma,\,2n)$ , represented by thin solid lines in Figs. 13(a, b). The sudden increase of $E_{4nX}$ some 4 MeV above the $S_{4n}$ threshold could not be reproduced by the model calculations. The difference between both EMPIRE and TALYS predictions stems mainly from the different nuclear level density prescriptions considered and the impact of their associated energy-dependent shell effect. In particular, a smaller shell effect in the level density tends to increase $E_{1nX}$ and decrease $E_{2nX}$ predictions.

The $(\gamma,\,n)$ partial cross sections to the first excited states in ${}^{207}$ Pb are compared with EMPIRE and TALYS calculations in Fig. 15. As statistical codes, the fine structures observed experimentally cannot be reproduced. However, both TALYS and EMPIRE satisfactorily reproduce the global trend and order of magnitude of the partial cross sections of populating the residual ${}^{207}$ Pb in its ground and first two excited states. In particular, the TALYS calculations describe well these three partial sections in the region of high incident energies, above 13 MeV, where the statistical neutron emission is dominant. The experimental cross sections $\sigma(\gamma,n_{>})$ for populating ${}^{207}$ Pb in any state higher than the second state are relatively well described by both codes.

Finally, based on the newly measured $E1$ photoneutron cross sections, it is possible to estimate the three main moments of the $E1$ distributions, namely

•

the integrated cross section $\Sigma_{TRK}$ defined in terms of the Thomas-Reiche-Kuhn (TRK) sum rule $\sigma_{TRK}=60NZ/A$ mb MeV, i.e.

\Sigma_{TRK}=\frac{1}{\sigma_{TRK}}\times\int_{0}^{\infty}\sigma_{abs}(\omega)d\omega

(15)

•

the centroid energy

E_{c}=\frac{\int_{0}^{\infty}\sigma_{abs}(\omega)d\omega}{\int_{0}^{\infty}% \sigma_{abs}(\omega)/\omega~{}d\omega}

(16)

•

and the polarizability

\alpha_{D}=\frac{\hbar c}{2\pi^{2}e^{2}}\int_{0}^{\infty}\frac{\sigma_{abs}(% \omega)}{\omega^{2}}d\omega~{}.

(17)

As detailed in Ref. [6], to estimate the moments, the measured photoabsorption cross sections above the neutron separation energy has been supplemented with the SMLO values below the neutron separation energy and the quasideuteron contribution has been excluded. The results are given in Table 4 where the statistical uncertainties have been simply calculated considering the maximum and minium values of the measured cross section. Table 4 also compares the present $E1$ moments with the values obtained in Ref. [6] on the basis of experimental data from Ref. [10]. The resulting moments are seen to be in rather good agreement and confirm the previous mysterious kink found between ${}^{208}$ Pb and ${}^{209}$ Bi polarizabilities in Ref. [6].

Table 4:

{}^{208}

Pb experimental integrated cross section

\Sigma_{\rm TRK}

, centroid energy

E_{c}

, and polarizability with their estimated uncertainties (Err) on the basis of the present measurements and comparison with the values obtained in Ref. [6]

	Present	Ref.[6]
$\Sigma_{\rm TRK}$	$1.27\pm 0.10$	$1.29\pm 0.03$
$E_{c}$ [MeV]	$14.20\pm 0.12$	$14.53\pm 0.29$
$\alpha_{D}$ [fm ${}^{3}$ / $e^{2}$ ]	$20.00\pm 1.30$	$19.82\pm 0.49$

VI Summary and conclusions

New measurements of photoneutron reactions on ${}^{208}$ Pb have been performed in the GDR energy region using 7.50 MeV – 38.02 MeV quasi-monochromatic laser Compton scattering $\gamma$ -ray beams of the NewSUBARU synchrotron radiation facility. A high-and-flat efficiency moderated neutron detection array of ${}^{3}$ He counters together with an associated neutron-multiplicity sorting method have been employed.

We obtained the cross sections for the $(\gamma,\,1nX)$ , $(\gamma,\,2nX)$ , $(\gamma,\,3nX)$ and $(\gamma,\,4nX)$ photoneutron reaction channels. The photoabsorption cross sections, extracted as the sum of the photoneutron cross sections, confirm the Saclay results of Ref. [9]. The present experiment extended the area of investigation of the photoabsorption excitation function above the maximum limit of 21 MeV reached in the Saclay experiment. Fine structures have been observed in the ${}^{208}$ Pb $(\gamma,\,n)$ reaction at incident energies lower than 13 MeV and compared with previous high resolution measurements. Average energies of neutrons emitted in each reaction have also been extracted based on the ring ratio data. Low average energies for the $(\gamma,\,1nX)$ and $(\gamma,\,2nX)$ neutrons at excitation energies above $S_{3n}$ indicate a significant contribution of charged particle emission reactions. The experimental photoneutron cross sections and average energies have been satisfactorily reproduced with both EMPIRE and TALYS calculations by slight adjustments of model parameters.

Based on the measured ring ratio values, we have extracted estimations on the total neutron emission spectra in the ${}^{208}$ Pb $(\gamma,\,n)$ reaction at incident energies lower than $S_{2n}$ , which confirm a gradual transition from discrete energy neutron emission to a statistical neutron emission towards $S_{2n}$ . Present estimations for the partial photoneutron cross sections for populating the ${}^{207}$ Pb residual in its ground and first two excited states reproduce the resonant structures observed in the previous experiments and significantly extend the previously investigated energy range. Until detailed time of flight measurements are available, such systematic estimations on a wide energy range are useful for microscopic descriptions of the GDR structure.

Finally, the $E1$ moments extracted from the present data are in rather good agreement with those obtained from the previous measurement of Ref. [10] and confirm the mysterious kink previously found between ${}^{208}$ Pb and ${}^{209}$ Bi polarizabilities in Ref. [6].

All the experimental results obtained in the present paper are available in numerical format in the Supplemental Material [37].

VII Acknowledgments

The authors are grateful to H. Ohgaki of the Institute of Advanced Energy, Kyoto University for making a large volume LaBr ${}_{3}$ :Ce detector available for the experiment. I.G. expresses her gratitude to Prof. M. Sin of the University of Bucharest, Romania, for her interest in this paper and valuable discussions. This work was supported by a grant of the Ministry of Research, Innovation and Digitization, CNCS - UEFISCDI, project number PN-III-P1-1.1-PD-2021-0468, within PNCDI III. This work was also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) Project IDs 279384907 (SFB 1245) and 499256822 (GRK 2891) and by the State of Hesse under grant ”Nuclear Photonics” within the LOEWE program (LOEWE/2/11/519/03/04.001(0008)/62). S.G. acknowledges financial support from the Fonds de la Recherche Scientifique (F.R.S.-FNRS) and the Fonds Wetenschappelijk Onderzoek - Vlaanderen (FWO) under the EOS Project nr O000422F. F.L.B.G., T.E., V.W.I. and L.G.P. acknowledges funding from the Research Council of Norway through its grant to the Norwegian Nuclear Research Centre (Project No. 341985).

References

[1] T. Kawano, Y.S. Cho, P. Dimitriou, D. Filipescu, N. Iwamoto, et al., IAEA Photonuclear Data Library 2019, Nuclear Data Sheets 163, 109 (2020).
[2] S. Goriely, P. Dimitriou, M. Wiedeking, T. Belgya, R. Firestone , J. Kopecky, M. Krtička, V. Plujko, R. Schwengner, S. Siem, H. Utsunomiya, S. Hilaire, S. Péru, Y. S. Cho, D. M. Filipescu, N. Iwamoto, T. Kawano, V. Varlamov and R. Xu, Reference Database for Photon Strength Functions, Eur. Phys. J. A 172, 55 (2019).
[3] A. Koning, S. Hilaire, S. Goriely, TALYS: modeling of nuclear reactions, Eur. Phys. J. A 59, 131 (2023).
[4] M. Herman, R. Capote, B.V. Carlson, P. Obložinský, M. Sin, A. Trkov, H. Wienke, V. Zerkin, EMPIRE: Nuclear Reaction Model Code System for Data Evaluation, Nuclear Data Sheets 108, 2655-2715 (2007).
[5] A. Bracco, E.G. Lanza, A. Tamii, Isoscalar and isovector dipole excitations: Nuclear properties from low-lying states and from the isovector giant dipole resonance, Progress in Particle and Nuclear Physics 106, 360-433 (2019).
[6] S. Goriely, S. Péru, G. Colò, X. Roca-Maza, I. Gheorghe, D. Filipescu, and H. Utsunomiya, E1 moments from a coherent set of measured photoneutron cross sections, Phys. Rev. C 102, 064309 (2020).
[7] R.R. Harvey, J.T. Caldwell, R.L. Bramblett, S.C. Fultz, Photoneutron Cross Sections of Pb ${}^{206}$ , Pb ${}^{207}$ , Pb ${}^{208}$ , and Bi ${}^{209}$ , Phys. Rev. 136, B126–B131 (1964).
[8] B.L. Berman, R.E. Pywell, S.S. Dietrich, M.N. Thompson, K.G. McNeill, J.W. Jury, Absolute photoneutron cross sections for Zr, I, Pr, Au, and Pb, Phys. Rev. C 36, 1286–1292 (1987).
[9] A. Veyssiere, H. Beil, R. Bergere, P. Carlos, A. Lepretre, Photoneutron cross sections of ${}^{208}$ Pb and ${}^{197}$ Au, Nuclear Physics A 159, 561–576 (1970).
[10] A. Tamii, I. Poltoratska, P. von Neumann-Cosel, Y. Fujita, T. Adachi, et al., Complete Electric Dipole Response and the Neutron Skin in ${}^{208}\mathrm{Pb}$ , Phys. Rev. Lett. 107, 062502 (2011).
[11] I. Poltoratska, R.W. Fearick, A.M. Krumbholz, E. Litvinova, H. Matsubara, P. von Neumann-Cosel, V.Yu. Ponomarev, A. Richter, A. Tamii, Fine structure of the isovector giant dipole resonance in ${}^{208}\mathrm{Pb}$ : Characteristic scales and level densities, Phys. Rev. C 89, 054322 (2014).
[12] M. Jingo, E.Z. Buthelezi, J. Carter, G.R.J. Cooper, R.W. Fearick, Studies of the Giant Dipole Resonance in ${}^{27}$ Al , ${}^{40}$ Ca , ${}^{56}$ Fe , ${}^{58}$ Ni and ${}^{208}$ Pb with high energy-resolution inelastic proton scattering under 0 ${}^{\circ}$ , Eur. Phys. J. A 54, 234 (2018).
[13] S. Bassauer, P. von Neumann-Cosel, A. Tamii, $\gamma$ strength function and level density of ${}^{208}\text{Pb}$ from forward-angle proton scattering at 295 MeV, Phys. Rev. C 94, 054313 (2016).
[14] E. Wolynec, A.R.V. Martinez, P. Gouffon, Y. Miyao, V.A. Serrão, M.N. Martins, Comment on photoneutron cross sections, Phys. Rev. C 29, 1137–1139 (1984).
[15] E. Wolynec, M. N. Martins, Discrepancies between Saclay and Livermore photoneutron cross sections, Revista Brasileira de Fisica 17, 56–87 (1987).
[16] V.V. Varlamov, A.I. Davydov, V.N. Orlin,New evaluated data on ${}^{206,207,208}$ Pb photodisintegration, Eur. Phys. J. A 57, 287 (2021).
[17] S. Amano, K. Horikawa, K. Ishihara, S. Miyamoto, T. Hayakawa, T. Shizuma, T. Mochizuki, Several-MeV $\gamma$ -ray generation at NewSUBARU by laser Compton backscattering, Nucl. Instrum. Methods Phys. Res. A 602, 337 (2009).
[18] K. Horikawa, S. Miyamoto, S. Amano, T. Mochizuki, Measurements for the energy and flux of laser Compton scattering $\gamma$ -ray photons generated in an electron storage ring: NewSUBARU, Nucl. Instrum. Methods Phys. Res. A 618, 209 – 215 (2010).
[19] H. Utsunomiya, I. Gheorghe, D.M. Filipescu, T. Glodariu, S. Belyshev, et al., Direct neutron-multiplicity sorting with a flat-efficiency detector, Nucl. Instrum. Methods Phys. Res. A 871, 135 (2017).
[20] I. Gheorghe, H. Utsunomiya, K. Stopani, D. Filipescu, T. Ari-izumi, et al., Updated neutron-multiplicity sorting method for producing photoneutron average energies and resolving multiple firing events, Nucl. Instrum. Methods Phys. Res. A 1019, 165867 (2021).
[21] W.L. Lv, Y.F. Niu and G. Colò, Learning about the structure of giant resonances from their $\gamma$ decay, Phys. Rev. C 103, 064321 (2021).
[22] H. Utsunomiya, T. Shima, K. Takahisa, D.M. Filipescu, O. Tesileanu, et al., Energy Calibration of the NewSUBARU Storage Ring for Laser Compton-Scattering Gamma Rays and Applications, IEEE Tran. Nucl. Sci. 61, 1252 (2014).
[23] H. Utsunomiya, T. Watanabe, T. Ari-izumi, D. Takenaka, T. Araki, et al., Photon-flux determination by the Poisson-fitting technique with quenching corrections, Nucl. Instrum. Meth. Phys. Res. A 896, 103 (2018).
[24] D. Filipescu, I. Gheorghe, K. Stopani, S. Belyshev, S. Hashimoto, S. Miyamoto, H. Utsunomiya, Spectral distribution and flux of $\gamma$ -ray beams produced through Compton scattering of unsynchronized laser and electron beams, Nucl. Instrum. Methods Phys. Res. A 1047, 167885 (2023), arXiv:2211.14650.
[25] D. Filipescu, Monte Carlo simulation method of polarization effects in Laser Compton Scattering on relativistic electrons, JINST 17 P11006 (2022), arXiv:2210.14669.
[26] D. Filipescu and I. Gheorghe, eliLaBr code, https://github.com/dan-mihai-filipescu/eliLaBr (2022).
[27] ADVACAM, MiniPIX (2022) https://advacam.com/minipix
[28] C. Granja, J. Jakubek, P. Soukup, M. Jakubek, D. Turecek, et al., Spectral and directional sensitive composition characterization of mixed-radiation fields with the miniaturized radiation camera MiniPIX Timepix2, JINST 17 C11014 (2022).
[29] T. Ari-Izumi, I. Gheorghe, D. Filipescu, S. Hashimoto, S. Miyamoto, H. Utsunomiya, Spatial profiles of collimated laser Compton-scattering $\gamma$ -ray beams, JINST 18 T06005 (2023).
[30] X-6 Monte Carlo Team, MCNP–A general N-particle transport code, version 6, Volume I: Overview and theory, Los Alamos National Laboratory, Los Alamos, NM, http://laws.lanl.gov/vhosts/mcnp.lanl.gov
[31] V. Weisskopf, Statistics and Nuclear Reactions, Phys. Rev. 52, 295 (1937).
[32] I. Gheorghe, H. Utsunomiya, S. Katayama, D. Filipescu, S. Belyshev, et al., Photoneutron cross-section measurements in the ${}^{209}$ Bi( $\gamma,\,xn$ ) reaction with a new method of direct neutron-multiplicity sorting, Phys. Rev. C 96, 044604 (2017) and Erratum Phys. Rev. C 99, 059901(E) (2019).
[33] R. Capote, M. Herman, P. Obložinský, P.G. Young, S. Goriely, et al., RIPL – Reference Input Parameter Library for Calculation of Nuclear Reactions and Nuclear Data Evaluations, Nuclear Data Sheets 110, 3107–3214 (2009).
[34] T. Renstrøm, H. Utsunomiya, H.T. Nyhus, A.C. Larsen, M. Guttormsen, et al., Verification of the detailed balance for $\gamma$ absorption and emission in Dy isotopes, Phys. Rev. C 98, 054310 (2018).
[35] A.C. Larsen, G.M. Tveten, T. Renstrøm, H. Utsunomiya, E. Algin, et al., New experimental constraint on the ${}^{185}$ W $(n,\,\gamma)^{186}$ W, Phys. Rev. C 108, 025804 (2023), arXiv:2301.13301.
[36] V.A. Plujko, O.M. Gorbachenko, R. Capote, P. Dimitriou, Giant dipole resonance parameters of ground-state photoabsorption: Experimental values with uncertainties, At. Data Nucl. Data Tables 123-124, 1–85 (2018).
[37] See Supplemental Material for the ${}^{208}$ Pb experimental data points: $(\gamma,\,1nX)$ , $(\gamma,\,2nX)$ , $(\gamma,\,3nX)$ and $(\gamma,\,4nX)$ cross sections and average photoneutron energies; as well as for estimations for the cross sections of the ${}^{208}$ Pb $(\gamma,\,n)$ reactions that populate the ${}^{207}$ Pb residual nucleus in its ground, first and second excited states and in any state higher than the second.
[38] T. Kondo, H. Utsunomiya, S. Goriely, I. Daoutidis, C. Iwamoto, Total and partial photoneutron cross sections for Pb isotopes, Phys. Rev. C 86, 014316 (2012).
[39] Lloyd Martin Young, Photoneutron cross sections and spectra from monoenergetic photons on Y, Pr, Pb and Bi in the Giant Resonance, PhD thesis, University of Illinois at Urbana-Champaign, 1972. https://core.ac.uk/download/pdf/29157805.pdf
[40] John Richard Calarco, Photoneutron spectra from monoenergetic photons on ${}^{208}$ Pb in the Giant Resonance, PhD thesis, University of Illinois, 1969.
[41] R. Alarcon, P.L. Cole, D.S. Dale, P.T. Debevec, L.J. Morford, Statistical and nonstatistical neutron decay of the giant electric dipole resonance of ${}^{208}\mathrm{Pb}$ , Phys. Rev. C 43, R2470–R2473 (1991).
[42] Z.W. Bell, L.S. Cardman, P. Axel, Fine structure in the ${}^{208}$ Pb photoneutron cross section between 9.9 and 11.2 MeV, Phys. Rev. C 25, 791–803 (1982).
[43] Gordon S. Mutchler, The Angular Distributions and Energy Spectra of Photoneutrons from Heavy Elements, PhD thesis, Massachusetts Institute of Technology, 1966. https://www.osti.gov/servlets/purl/4566721
[44] N.K. Sherman, H.M. Ferdinande, K.H. Lokan, C.K. Ross, Photoneutron Fine Structure in ${}^{208}\mathrm{Pb}$ and Multipole Resonance, Phys. Rev. Lett. 35, 1215–1219 (1975).
[45] R.W. Finlay, J. R. M. Annand, T. S. Cheema, J. Rapaport, F. S. Dietrich, Energy dependence of neutron scattering from ${}^{208}$ Pb in the energy range 7–50 MeV, Phys. Rev. C 30, 796 (1984).
[46] S. Hilaire, M. Girod, S. Goriely, A.J. Koning, Temperature-dependent combinatorial level densities with the D1M Gogny force, Phys. Rev. C 86, 064317 (2012).

Photoneutron cross section measurements on 208208{}^{208}start_FLOATSUPERSCRIPT 208 end_FLOATSUPERSCRIPTPb in the Giant Dipole Resonance region