TALYS calculations of prompt ﬁssion observables and independent ﬁssion product yields for the neutron-induced ﬁssion of 235 U

The TALYS nuclear reaction code’s Hauser-Feshbach statistical decay model has been adapted in order to calculate prompt ﬁssion neutron and γ -ray observables by iterating over deexciting ﬁssion fragments. Several ﬁs-sion fragment generators such as GEF, HF 3 D, and SPY were employed to provide TALYS with databases. These databases contain standardized tables with ﬁssion fragment yields, mean excitation energies and their widths, and average total kinetic energy, as a function of charge and mass number of primary ﬁssion fragments. The resulting calculations, including prompt particle multiplicities, spectra, average energies, and independent ﬁssion product yields, were compared with experimental and evaluated data. This work ﬁrst outlines the new methodology implemented in TALYS and examines the effects of three important parameters on the ﬁnal evaporation data. Furthermore, the neutron-induced ﬁssion of 235 U is investigated in detail as a function of incident energy. The results from TALYS, with input from GEF and HF 3 D, were compared with available experimental data and the results of the stand-alone GEF code. The proposed methodology contributes to an improved capability to model the ﬁssion process.


Introduction
Prompt particle emission from primary fission fragments during the nuclear fission process generates various correlated prompt fission quantities. These include independent a e-mail: fujio.k.aa@m.titech.ac.jp (corresponding author) fission product yields, prompt fission neutron multiplicities ν n , prompt fission neutron spectra (PFNS), prompt fission γ -ray multiplicitiesν γ , prompt fission γ -ray spectra (PFGS) and isomeric yield ratios. The characteristics of these prompt fission quantities contain valuable information for a detailed understanding of fission physics and are beneficial for many nuclear applications; e.g., the evaluation of delayed neutron yields for reactor safety [1], design of a transmutation systems for minor actinides in nuclear waste management [2], producing radioisotopes for diagnosis and therapy [3], and studying nuclear astrophysics, especially during the fission cycle in the r-process.
However, comprehensive modeling of these observables for different fissioning systems at various excitation energies remains challenging. Moreover, the estimation of accurate fission observables required for the aforementioned applications still relies on empirical models. For instance, Wahl systematics for fission yields [4][5][6], and Madland-England model [7] for isomeric ratios have been used in ENDF/B-VIII [8] and JENDL-4.0 [9] evaluations. The Los Alamos (Madland-Nix) model for PFNS [10] is based on physical considerations; however, it takes into account only one fragmentation (the so-called most probable fragmentation).
During the fission process, the compound nucleus undergoes deformation, therefore, may approach the scission point. In binary fission, the nucleus splits into two complementary fission fragments. Typically, this process may lead to several hundred different fission fragment pairs. After the mass split, the Coulomb repulsion between the nascent fragments accelerates them, back-to-back. The summation of the excitation energy of the compound nucleus and the Q-value of the reac-tion is divided between the total kinetic energy (TKE) and the total excitation energy (TXE). When the fragments reach their full acceleration, the highly excited fragments deexcite by emitting prompt fission neutrons and prompt fission γrays until they reach their metastables or ground states. Then the fission products undergo β − decay. β − decay takes place in neutron-rich independent fission products. The delayed neutrons and delayed γ -rays are emitted from the independent fission products to become the final stable or long-lived cumulative fission products.
Since the fission process consists of several different physical phenomena, the stages are described by different physics. The current modeling for the prompt fission observables often does not consider correlations among the independent fission product yield and the other prompt fission observables, which partly causes inconsistency in the model calculations [11]. Therefore, more recently, studies are also being conducted to evaluate prompt fission observables consistently.
In this study, we propose the use of a common, consistent, well-documented, and transparent nuclear deexcitation code to probe biases and systematic errors in the initial assumptions in fission physics. These systematic errors can be compensated by fine-tuning parameters during the deexcitation stage. We employ existing fission fragment generators to yield pre-neutron emission fission fragment mass and energy yield distributions, including both microscopic models [25] and phenomenological models [19,22,23]. The nuclear reaction model code TALYS [26] calculates the evaporation of neutrons and γ -rays from the excited fission fragment to obtain prompt fission observables. In the Hauser-Feshbach statistical decay calculation, TALYS adopts deterministic techniques similar to the HF 3 D model [22,23], while so far the calculation is limited to the first-chance fission only, in which no neutron is evaporated prior to fission.
We examine the detailed calculation of the neutroninduced fission on 235 U and show the calculated neutron and γ -ray multiplicities, PFNS, PFGS, and independent fission product yields, which are compared with the available experimental and evaluated data.

Fission models in TALYS
The TALYS code (version: 1.96) includes several modules for calculating fission observables. The first implementation of mass yields in TALYS was based on a revised version of the multi-modal random neck-rupture model (MM-RNRM). The original model was developed by Brosa [27] to calculate the properties of fission fragments at zero temperature. The temperature is added to the model calculation of the potential energy landscape of the nucleus [28]. A search for the fission channels in deformation space yields the superlong (SL), standard I (ST I), and standard II (ST II) fission barriers and pre-scission shapes as a function of temperature. The obtained temperature-dependent fission barrier and pre-scission shape parameters serve as an input option for the fragment mass distribution computations in TALYS. However, this model fails to agree with most experimental fission yields [29]. Another approach to enhance the preneutron mass distributions was to combine the GEF code with TALYS through direct code translation to Fortran. However, this method locks the code version and does not facilitate a smooth transition to an updated fission code version. It can also result in large discrepancies compared to recent versions of GEF [30].
A third approach was recently implemented [31,32] and makes use of phenomenological and empirical codes, such as GEF [19,30], HF 3 D [22,23], and SPY [25,33], as fission fragment generators to create a fission fragment database for TALYS. Employing these databases one can facilitate comprehensive model comparisons that address systematic errors and model defects. Figure 1 shows the workflow of the intended approach. The starting point involves the initial conditions defined by the fission fragment yield Y ff (Z , A) for each fragment charge Z and mass number A, the mean excitation energyĒ x of its Gaussian distribution, the width σ E x of the excitation energy distribution, and average total kinetic energy (TKE) for fragment pair. TALYS uses these databases to calculate evaporation data, e.g., neutron and γray emissions sequentially until both fragments are relaxed and reach their ground states. In this methodology, one may iterate the calculations to tweak parameters and perform sensitivity studies to reveal compensating trends and other correlated effects [34,35].

Generating fission fragment distribution database
This section describes two models of fission fragment generators used to produce fission fragment databases as input to TALYS. The first model, HF 3 D (Hauser-Feshbach Fission Fragment Decay), which is a statistical decay code [22,23], has been developed to study prompt neutron and γ -ray observables from the neutron-induced fission of 235 U, 238 U, and 239 Pu by the Hauser-Feshbach statistical decay code, Fig. 1 The proposed workflow, where various fission models are used to create pre-neutron emission databases for TALYS. The necessary inputs to TALYS are fission fragment yields Y ff (Z , A), the mean excitation energyĒ x , its width σ Ex , and average total kinetic energy (TKE). The output from the Hauser-Feshbach formalism of TALYS yields deexcitation data that can be compared with experimental data and with results from other stand-alone model codes. The methodology facilitates comprehensive model comparison addressing systematic errors and model defects BeoH, with the deterministic approach so that it can calculate contributions including fission fragments having very low yield and also to ensure that isomeric production comes from a proper decay scheme including discrete level information. In the HF 3 D model, the fission fragment distribution is generated by fitting the experimentally available Y ff (A), TKE(A) data and distributing excitation energy between the complement fission fragments based on a ratio of nuclear temperatures, the anisothermal model parameter R T [36]. The generated (Y ff (Z , A),Ē x , σ E x , TKE) are stored in the TALYS fission fragment database in tabulated format.
The second model, GEF [37], is a Monte Carlo-based phenomenological fission model that generates both pre-neutron emission quantities as well as deexcitation data, albeit by employing a simple evaporation scheme. In this work, GEF (version: 2021/1.2) was employed to produce the mass and charge yields of primary fission fragments and corresponding energies [30]. The Monte Carlo sampling was based on 1 million fission events for every reaction. The list-mode data feature was activated to allow for an event-by-event output of each fission simulation. In total, data files from 737 fissioning nuclei were calculated, ranging from 76 Os to 115 Mc, at excitation energies ranging from 0 to 20 MeV. The generated (Y ff (Z , A),Ē x , σ E x , TKE) were stored in the TALYS fission fragment database in a tabulated format. The data were modified to include only first-chance fission. In the future, TALYS will be responsible for assigning the multi-chance fission probability, thus allowing for interpolation between first-chance fission files.

Generating the initial distributions
TALYS calculates fission fragment deexcitation by integrating the Hauser-Feshbach statistical decay calculation over the distribution for all fission fragments, instead of Monte Carlo sampling. The main advantage of the deterministic treatment is that this can take account into the small probability of fission events. A technique for this is explained in Refs. [22][23][24]. Figure 2 shows a schematic view of the multiple neutron and γ -ray emission process from an excited fission fragment (Z , A). The initial fragment (Z , A) is in the continuum excited state and it decays to either the continuum or discrete states of (Z , A − 1) by emitting a neutron and either the continuum or discrete states of (Z , A) by emitting γ -rays. This sequential process continues until each state reaches the ground or isomeric state.
To calculate such a deexcitation process, it requires the primary fission fragments' quantities. TALYS reads such initial distributions from the fission fragment databases described above and builds an excitation energy distribution G(E x ) and a spin-parity distribution R(J, π, E x ). The initial population P 0,Z ,A of a given fission fragment is expressed as follows: TALYS assumes a Gaussian fragment excitation energy distribution with the mean excitation energyĒ x and the width σ E x , Eur. Phys. J. A (2023) 59:178 The angular momentum generation remains a challenge for contemporary fission modeling [38][39][40][41][42]. TALYS adopts the spin-parity distribution R(J, π, E x ) by following the functional dependency of the level density formula [43]. Since the parity distribution is 1/2, the spin-parity distribution R(J, π, E x ) is assumed to be expressed as A scaling factor f 2 is introduced in order to assure a reasonable agreement with experimental data. This new scaling parameter controls the distribution of angular momenta for primary fission fragments.
Another important parameter is the spin cut-off parameter σ 2 (E x ) entering a level density model [44,45], defined as where Δ is the pairing energy correction,ã is the asymptotic level density parameter, and a is the level density parameter. TALYS uses another scaling factor, f s , which represents a global adjustment factor of the nuclear spin cut-off parameter, applicable for the level density model. This parameter is multiplied to σ 2 (E x ) and affects all independent fission products which are populated by the emission of prompt neutrons. The distributions G(E x ) and R(J, π, E x ) should satisfy the normalization condition, i.e., The other ingredients in the Hauser-Feshbach model calculations include the optical model potentials for neutron and charged particles, the level density parameters, the γray strength function, and the discrete level properties for all residual nuclei. For this study, we used the Koning-Delaroche global optical potential [46] for neutrons, a composite level density formula (Gilbert and Cameron model [47]) using the level density parameters and systematics from Ref. [48], IAEA-CRP SMLO 2019 tables and IAEA GSF CRP 2018 for E1 and M1 γ -ray strength functions [49][50][51], and the discrete level data from RIPL-3 [52].

Calculated prompt fission observables
The statistical Hauser-Feshbach calculation is performed from the initial condition (J, π, E x ) of a fission fragment (Z , A) while the results for the emitted neutrons and γ -rays are ν n (Z , A, J, π, E x ) and ν γ (Z , A, J, π, E x ). These results are then weighted using the initial fission yields to calculate the final observables.
The prompt fission neutron multiplicityν n (Z , A) of a given fission fragment is calculated by adding the fission neutron production resulting from the initial population is the neutron multiplicity and PFNS (CMS)(Z ,A,J,π ) ( ) is the neutron kinetic energy spectrum in center-of-mass system (CMS). The PFNS (CMS)(Z ,A) ( ) produced by the fragment (Z , A) is given by The fission neutron spectrum of a fragment in the laboratory frame, PFNS (LAB)(Z ,A) (E), is converted from the centerof-mass system by Feather's formula [53][54][55]: where E f is the kinetic energy per nucleon of the fission fragment. The PFNS is usually plotted as ratio to a Maxwellian spectrum: In this work, we adopted T M = 1.32 MeV for comparison with experimental and evaluated data. The γ -ray observables can be obtained in a similar way to the neutron observables. The prompt fission γ -ray multiplicityν γ (Z , A) is calculated by using the initial population and the γ -ray spectrum PFGS (Z ,A,J,π ) ( ) in the center-of-mass system: The PFGS (Z ,A) (E) produced by the fragment (Z , A) is given by Due to prompt neutron emission, the mirror symmetry seen in the primary fission fragment yield is broken and characteristic peaks appear. The independent fission product yield Y (Z , A) is deduced from the primary fission fragment yields Y ff (Z , A): where P Z ,A (ν) is the probability to emit ν neutrons by the primary fragment (Z , A).

Sensitivity of fission observables on input parameters
We identified several key input parameters essential for controlling the statistical decay calculation of fission fragments.
To evaluate the observables from the fission fragment deexcitation, a sensitivity study was performed on the spin-parity distribution and the number of continuum states. The scaling factor f 2 in Eq. (3) modulates the angular momentum distribution. The other scaling factor adjusts the spin cutoff parameter, σ 2 (E x ), by the multiplication factor f s . The energy width of discretized continuum state Δ bins is given by where N is the number of continuum states, E max x is the maximum excitation energy, and E level

Sensitivity on spin cut-off parameters
In our previous investigation [22], the HF 3 D model indicated that f = 2.5 was necessary to reproduce the neutron observables in the neutron-induced fission of 235 U up to 5 MeV incident energy. We have chosen to set f 2 to a range of 4-6, based on rough estimates from the same study [22]. Additionally, another recent study [56] has suggested that setting f s within the range of 0.4-0.5 is optimal for reproducing isomeric yield ratios observed in experiments involving neutron capture and other reactions.
We conducted a sensitivity analysis by varying the f 2 parameter between 3 and 6 (changing in unit steps), and the f s parameter between 0.4 and 1.0 (changing by 0.2 in step size). The input for the TALYS calculation was based on (Y ff (Z , A),Ē x , σ E x , TKE) obtained from HF 3 D. Table 1 presents a summary of the calculated values forν γ ,ν n , γ , and n . Results obtained with f 2 = 4 yield better ν n values, but smallerν γ compared to those obtained with f 2 = 5, 6. The results indicate a clear trend in which γ decreases with increasing values of f 2 and f s , respectively.
n also decreases as f s increases. Figure 3 shows the calculated PFNS (LAB) as ratio to a Maxwellian spectrum with f 2 = 4 that gives goodν n . f s value has a significant impact on both the peak position and shape of the PFNS (LAB) above 2 MeV. Our priority is to accurately reproduce the neutron observables, particularlyν n . Therefore, the optimal spin cutoff parameters are f 2 = 4 and f s = 0.4 (see Table 1).

Sensitivity on the number of continuum states
A sensitivity analysis was performed to study the role of the number of continuum states, N . The HF 3 D model employs Table 1 Sensitivity of prompt neutron and γ -ray multiplicities (ν n andν γ , respectively) and average energies ( n and γ , respectively) to changes in the spin cut-off parameters for the 235 U(n th ,f) reaction. n is given in the laboratory frame (LAB), while γ is given in the center-of-mass system (CMS)  Table 2 Sensitivity of prompt neutron and γ -ray multiplicities (ν n andν γ , respectively) and average energies ( n and γ , respectively) to changes in N for the 235 U(n th ,f) reaction. n is given in the laboratory frame (LAB), while γ is given in the center-of-mass system (CMS) a constant Δ bins size of 100 keV for both the primary fission fragment and the decaying nucleus, while TALYS partitioned the excitation energy into the specified number of N . To ensure a maximum Δ bins size of 500 keV, we set the N range to be between 50 and 300. Regarding γ -ray observables, bothν γ and the prompt fission γ -ray spectrum (PFGS) multiplied byν γ exhibit sensitivity to N . As seen in Table 2, the higher N value leads to more emitted γ -rays. A broader spin distribution also increases the average spin of the fission fragment. Fission fragments with high angular momenta are unable to reach low-excitation levels directly since the γ -ray transition is generally governed by dipole transitions. Therefore, larger values of initial spins lead to an increase in the number of emitted γ -rays. This stands in contrast to neutron emissions, highlighting a significant difference between the two types of emissions. Theν γ also becomes large as N increases because the number of continuum states increases as well. TALYS sequentially gives the same number of bins to the residual nucleus as the input, therefore, the size of bins becomes smaller as neutrons are emitted. Figure 4a shows ν γ (A), and (b) shows PFGS with different number of N . It is suggested that the prominent peak below 0.2 MeV in PFGS is mainly caused by the number of continuum states effect. We have verified that the impact of the bins parameter on neutron observables is negligible, and thus we select N = 300 for improved consistency with experimental and evaluated data. Regarding independent fission product yield, we have confirmed that the effects of these parameters are sufficiently insignificant. Consequently, we have identified the input parameter set, consisting of f 2 = 4, f s = 0.4, and N = 300, as the optimal selection for obtaining accurate prompt neutron multiplicity and better γ -ray observables at thermal energy. The results of the stand-alone GEF code are also plotted, they being denoted as GEF. The comparison between TALYS(GEF) and stand-alone GEF highlights important differences in the deexcitation procedure, and might shed light on valuable model discrepancies, given that the very same input pre-neutron data was employed.  [58]. Table 3 shows a summary of the calculatedν n ,ν γ , average energy of neutrons and γ -rays, denoted as n and γ , respectively, obtained from TALYS(GEF), TALYS(HF 3 D), and GEF, along with experimental and evaluated data for comparison. Theν n from TALYS(GEF) is lower than that from GEF becauseν n (A) is lower over a wide mass range. Comparing TALYS(GEF) and GEF with experimental data, theν n (A) underestimates in the light fragment and overestimates in the heavy fragment. These deviations cancel each other, therefore, both ofν n are in good agreement with the evaluated data. For n , both TALYS results agree with the evaluated value. ν n (A) of 235 U(n,f) at the incident energies of 0.5 and 5.55 MeV are plotted in Fig. 5b. The results from TALYS(GEF) are in better agreement with the experimental data from Ref. [59] at both 0.5 and 5.55 MeV, compared to the thermal data. This could probably be due to the important role these fission data play in fine-tuning the GEF excitation energy sharing. The energy-sorting mechanism in GEF [37] is responsible for the higher average neutron multiplicity from the heavy fragments. As measured by Müller et al. [59], the added excitation energy is mainly shared among the heavy fragments, leading to an increase inν n (A) from the heavy fragments. Since TALYS(GEF) uses the excitation energies of GEF directly, it reproduces this effect where a significant increase inν n (A) is observed for the heavy fragment mass region. In contrast, the results obtained from TALYS(HF 3 D) do not exhibit this observed trend. This discrepancy is attributed to the fact that the energy-sorting in the HF 3 D model is adjusted to reproduce the energy dependence ofν n [22,23]. TALYS reflects the difference in the energy-sorting mechanism from GEF and HF 3 D model.  Figure 6 represents theν n as a function of incident neutron energy ranging from thermal up to 5 MeV. The calculations show an enhanced number of emitted neutrons as a function of incident energy, which is consistent with known data. TALYS(HF 3 D) successfully reproduces the evaluated value at thermal energy. On the other hand, TALYS(GEF) underestimates the evaluated data at thermal energy by about 0.1, but it agrees with the data as the incident energy increases.

Prompt neutron observables
The calculated PFNS (LAB) as ratio to a Maxwellian spectrum at T M = 1.32 MeV is shown in Fig. 7. Both PFNS (LAB) results of TALYS(GEF) and TALYS(HF 3 D) do not describe the shape of experimental data above 2 MeV outgoing energy, TALYS(HF 3 D) has a harder tail and TALYS(GEF) has a softer tail of PFNS (LAB) . Moreover, the results differ much between TALYS(GEF) and GEF, which reflects the difference in the evaporation scheme in these codes.

Prompt γ -ray observables
The experimental data on the number of emitted γ -rays, ν γ (A), is unfortunately still scarce. In recent years, however, the accuracy of theν γ measurements has improved due to the development in instrumentation and novel detector systems [60]. As a consequence, recent investigations [64] shed more light on the saw-tooth shape inν γ (A).
The mass-dependent γ -ray multiplicity,ν γ (A), calculated in this work, shows some significant differences between TALYS(GEF) and GEF as seen in Fig. 8a. While GEF has a flatterν γ as a function of the fragment mass, TALYS results exhibit a more pronounced saw-tooth shape. The results agree fairly well with the experimental data, albeit smaller γ -ray multiplicities are observed around fragment mass numbers A = 100 to 110. Other contemporary fission codes, such as FIFRELIN [65], PbP, and DSE [66], also verify the saw-tooth shape in several fissioning systems.
Another important fission observable is the evolution of the total number of released γ -rays as a function of the incident neutron energy, which can be seen in Fig. 8b. A higher excitation energy leads to an increasedν γ similar to the trend observed in the prompt fission neutrons. TALYS agrees generally with GEF and with the evaluated data files. The current evaluation data indicate higherν γ as it is based on the current experimental and theoretical underlying data. Taking into account these factors, TALYS tends to underestimate theν γ . This underestimation is not fully understood yet but could arise from a non-optimal spin cut-off parameter, which is known to strongly influence the number of emitted γ -rays. Another possible explanation is that older experimental data pointed to a lower total number of emitted γ -rays, which is around 6.5 γ -rays/fission as can be seen in Table 3. In contrast, this value has been increasing as a function of time, where recent experimental studies indicated a number as high as 8.2 γ -rays/fission. The older data might still be used in fine-tuning of excitation energies and spin distributions, which could affect the current calculations.
Another aspect is the energy spectrum of the γ -rays. Figure 9 shows calculated PFGS multiplied byν γ compared with experimental and evaluated data. A pronounced peak is observed around 0.2 MeV in both TALYS(HF 3 D) and TALYS(GEF) results. This peak is due to the treatment of bins in the TALYS as mentioned in Sect. 4.1.2. As a consequence, the average energy of emitted γ -rays tends to be lower for TALYS results compared to other data, as shown in Table 3. Figure 10 shows the calculated independent fission product yield as a function of mass Y (A) for 235 U(n th ,f). The results show an agreement between TALYS(GEF) and stand-alone Table 3 Multiplicities of prompt neutrons and γ -rays,ν n andν γ , and the average energies of neutrons and γ -rays, n and γ for the 235 U(n th ,f) reaction. n is given in the laboratory frame (LAB), while γ is given in the center-of-mass system (CMS)    reproduce the peak while the original HF 3 D model reproduces [22]. The calculated Y (Z , A) for several fission products were compared with the experimental data from Rudstam et al. [67], as shown in Fig. 11. Both the TALYS(GEF) and TALYS(HF 3 D) results align well with the experimental data, with the calculated/experimental (C/E) ratio approaching 1.0 for yields above 0.001. However, there is greater dispersion in the C/E for Y (Z , A) below approximately 0.001. Both TALYS results, TALYS(GEF) and TALYS(HF 3 D), are roughly consistent with the tendency shown in the experimental data where experimental uncertainties are also large.

Conclusion
A new approach for calculating nuclear fission data is proposed. This procedure consists of the use of fission codes to provide input for the Hauser-Feshbach formalism of TALYS. Fig. 11 Comparison of the calculated independent fission product yield Y (Z , A) with the experimental data reported by Rudstam et al. [67] in 235 U(n,f) at thermal energy To illustrate this method, two fission codes, GEF and HF 3 D, were utilized to create databases that were subsequently used as input for TALYS. The TALYS results are compared with both experimental/evaluated data and the results of the standalone GEF code. The development of such a database has streamlined the inclusion of fission data from other code developers due to the implementation of a standardized format that contains information on fission fragment yields and their excitation energy distributions. The adoption of this new standard format encourages fission model developers to prepare and supply TALYS with input data.
The neutron-induced fission of 235 U at incident energies from thermal up to 5 MeV was chosen for the validation of this procedure. We investigated the sensitivity of fission observables on mainly three parameters, namely the number of continuum states (N ), the scaling factor for the spin cut-off parameter in the level density formula ( f s ), and the scaling factor for the primary fission fragment angular momentum population ( f 2 ). We chose f 2 = 4 and f s = 0.4 because they showed a good agreement with neutron observables. Moreover, we selected N = 300 particularly fine-tuned for the average total γ -ray multiplicity.
During the investigation, TALYS showed a decent agreement with the experimental data of prompt neutron observables, e.g., the saw-tooth shapes ofν n (A) and the total average number of prompt neutrons as a function of incident energy. Significant differences were however observed between TALYS(GEF) and TALYS(HF 3 D) results, highlighting again the importance of different physical considerations and assumptions on the energy partition in fission. The prompt neutron spectrum results of both TALYS(GEF) and TALYS(HF 3 D) do not describe the shape exhibited by the experimental data, pointing to the need for further investigations and model adjustments.
Regarding the prompt γ -ray observables, the TALYS results tend to show a saw-tooth shape ofν γ (A), which agrees with recent experimental findings. Despite this, theν γ is underestimated in absolute scale. The total prompt fission γ -ray spectrum (PFGS) exhibits a strong peak around 0.2 MeV, which is believed to be a binning effect in the TALYS code, and it underestimates the experimental data in some regions e.g., below 2 MeV and after 6 MeV.
The independent fission product yield Y (A) of TALYS(GEF) is in overall good agreement with the experimental data and the recent evaluations, reproducing the pronounced peaks (at A = 134, 138, 100, 94) and dips in the independent fission product yield structure. For Y (Z , A), both TALYS results are generally consistent with the literature data.
In the future, more correlated physical quantities will be studied, such as cumulative fission product yield, γ -ray observables, and isomeric yield ratios. TALYS has been earlier employed, to calculate isomeric yield ratios based on different assumed angular momentum generation in the nascent fission fragments [68][69][70][71]. However, these calculations focus on individual isomeric yield ratios and utilize an external spin-energy matrix implementation. Moreover, no particular emphasis was put on the other fission observables. The outlook from this work is to exploit the new implementation of internal looping to generate systematic global isomeric yield ratio investigations as a function of mass and excitation energy taking into account the correlation between different fission observables. Technology, Tokyo, Japan by the Nuclear Regulation Agency of Japan, for supporting this work. A. A. would like to acknowledge Liljewalch travel scholarships and Ingegerd Berghs stiftelse for their research grants.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The datasets generated and/or analyzed during the current study are available and reproducible in/by the open-source nuclear reaction model code, TALYS.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.