Product yields for the Photofission of ²³²Th, ^234,238U, ²³⁷Np and ^239,240,242Pu actinides at various incident photon energies

Product yields for the Photofission of 232 Th, 234,238 U, 237 Np and 239,240,242 Pu actinides at various incident photon energies 1 corresp first-author Corresponding author, m.pahlavani@umz.ac.ir Corresponding author, m.pahlavani@umz.ac.ir m.pahlavani@umz.ac.ir 1 Department of Physics, Faculty of Basic Science, University of Mazandaran, P.O.Box 47415-416, Babolsar, Iran Department of Physics, Faculty of Basic Science, University of Mazandaran, P.O.Box 47415-416, Babolsar, Iran Photofission fragments mass yield for 232 Th, 234,238 U, 237 Np, and 239,240,242 Pu isotopes are investigated. The calculations are done using a developed approach based on Gorodisskiy’s phenomenological formalism. The Gorodisskiy’s method is developed to be applied for neutron induced fission. Here we revised it for application to photofission. The effect of emitted neutron prior to fission on fission fragment mass yields has also been studied. The peak to valley ratio is extracted for the 240 Pu isotope as a function of energy. Obtained results of present formalism are compared with the available experimental data. Satisfactory agreement is achieved between results of present approach and experimental data. Fission fragments Fragment mass yield Photofission Neutron fission Heavy nuclei Peak to vally 1 Introduction 1 In spite of peaceful application in energy production, different features of various kinds of fission are not known yet. Photofission is one way to induce fission on actinides (U, Pu, Np, etc.). Although photofission was discovered more than fifty years ago [ 1 1 ], its theoretical modeling was produced with many challenges. Various theoretical approaches [ 2 2 - 7 7 2-7 ] have been presented to evaluate fragment mass yields for neutron/proton induced fission, but there is no comprehensive method to study photofission. The distribution of fragment mass yield plays a significant role in examining the validity of every theoretical method of fission. If the calculated fragment mass yields are well fitted with measured data, this means that the selected model is suitable to be applied for the fission process. In the statistical model of Wilkins [ 8 8 , 9 9 ] and its developed version by Moreau [ 10 10 ], the probability of fission for neutron-induced fission is presented as follows: <math display="block" id="M1"><mrow><mi>P</mi><mo>∼</mo><msup><mtext>e</mtext><mrow><mrow><mo>(</mo><mrow><mo>−</mo><msub><mi>E</mi><mrow><mtext>total</mtext></mrow></msub><mo>/</mo><mi>T</mi></mrow><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mtext>e</mtext><mrow><mrow><mo>(</mo><mrow><mo>−</mo><mrow><mo>(</mo><mrow><msub><mi>E</mi><mrow><mtext>mac</mtext></mrow></msub><mo>+</mo><msub><mi>E</mi><mrow><mtext>mic</mtext></mrow></msub></mrow><mo>)</mo></mrow><mo>/</mo><mi>T</mi></mrow><mo>)</mo></mrow></mrow></msup><mn>,</mn></mrow></math> where P is the probability of fission and E total includes all kinds of energies presented in the evaluation of compound nucleus from saddle to scission point. E mac and E mic are the macroscopic and microscopic energies, respectively. The relation between fission fragments mass yield and fragment mass number indicates that the probability appears in the Gaussian form. Therefore, the theoretical neutron induced fission models are developed in the form of a multi-Gaussian shape and are formulated as a sum of several Gaussian type function [ 11 11 - 13 13 ]. The old formalism presented by Wang Fu-cheng [ 5 5 ] was applicable for reproducing low-energy fission and contains five Gaussian terms. The formalism of Robert Mills[ 14 14 ], in which the dependence on energy was neglected, includes four Gaussian terms in the quadratic structure. The GEF computer code [ 15 15 ] is developed based on this method with various assumptions and many adjustable parameters. The complicated Langevin-Brownian method has been employed to study neutron induced fission and lunched some success in recent years [ 16 16 - 20 20 16-20 ]. The Gorodisskiy phenomenological method [ 21 21 ] has been recently successful in reproducing the experimental data for neutron induced fission of actinides. It is clear that the photofission process is similar to neutron induced fission except in multi-polarity absorptions and time scale of the fission process. These absorptions happen in the beginning of the process, therefore, after formation of the compound nucleus, the fission process of the neutron induced fission will become similar to photofission. According to the Bohr’s hypothesis, the compound nucleus rapidly loses its total formation memory except conserved degrees of freedom. Therefore the mode of decay of the compound nucleus does not depend on the way the compound nucleus is formed. Thus for the same excited configuration of the compound nucleus, the photofission fragment mass yield is expected to be approximately similar to the neutron induced fission. In the present study, we examined an extended version of the Gorodisskiy’s method to generate the photofission fragment mass yield. We modified the relative contributions of the symmetric and the asymmetric yields, ( Y s and Y a ) of the Gorodisskiy’s method and the variance of asymmetric fission ( σ a ) as presented in the next section. The paper is constructed as follows: Theoretical model used to calculate fission fragments mass distribution is presented in Sect. 2. In Sect. 3, theoretical results are compared with experimental data for different excitation energies and different numbers of emitted neutron prior to fission for 232 Th, 234,238 U, 237 Np, and 239,240,242 Pu actinides. Finally, a short conclusion is presented in Sect. 4. 2 Description of theoretical Model 1 Studies of fission fragments mass distribution of most fissioning nuclei indicate that the yield of fission fragments as a function of fragment mass number ( A ) is a Gaussian function with 4 [ 23 23 ] or 3 [ 22 22 ] modes. These modes are symmetric ( SL ) or asymmetric ( S 1, S 2, and S 3) related to the type of the compound nucleus and its excitation energy. In Gorodisskiy original approach that has developed for neutron induced fission, two main parts are considered to generate the mass distribution of fission fragments. The yield of fragment with mass number higher than the half of compound nucleus is indicated by Y H and its conjugate with the mass number lower than the half of compound nucleus is shown by Y L . These modes are presented as an exponential function with a small deviation from Gaussian shape. Therefore, it is convenient to use the Charlier’s distribution [ 21 21 , 24 24 , 25 25 ] to construct the fission fragment mass yield <math display="block" id="M2"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>u</mi><mo stretchy="false">)</mo><mo>=</mo><mn>1</mn><mo>−</mo><msub><mi>γ</mi><mn>1</mn></msub><mtext> </mtext><mrow><mo>(</mo><mrow><mn>1</mn><mo>/</mo><mn>2</mn><mtext> </mtext><mi>u</mi><mo>−</mo><mn>1</mn><mo>/</mo><mn>6</mn><mtext> </mtext><msup><mi>u</mi><mn>3</mn></msup></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>γ</mi><mn>2</mn></msub><mtext> </mtext><mrow><mo>(</mo><mrow><mn>1</mn><mo>/</mo><mn>24</mn><mtext> </mtext><msup><mi>u</mi><mn>4</mn></msup><mo>−</mo><mn>1</mn><mo>/</mo><mn>4</mn><mtext> </mtext><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo>/</mo><mn>8</mn></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where <math id="M3"><mrow><mi>u</mi><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><mi>M</mi><mo>−</mo><mo><</mo><mi>M</mi><mo>></mo><mo stretchy="false">)</mo></mrow><mi>σ</mi></mfrac></mrow></math> and coefficients γ 1 and γ 2 are respectively called dissymmetry and excess and are given by <math display="block" id="M4"><mtable columnalign="left"><mtr><mtd><msub><mi>γ</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mo><</mo><mo stretchy="false">(</mo><mi>M</mi><mo>−</mo><mo><</mo><mi>M</mi><msup><mrow><mo>></mo><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mo>></mo></mrow><mrow><msup><mi>σ</mi><mn>3</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>γ</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mo><</mo><mo stretchy="false">(</mo><mi>M</mi><mo>−</mo><mo><</mo><mi>M</mi><msup><mrow><mo>></mo><mo stretchy="false">)</mo></mrow><mn>4</mn></msup><mo>></mo></mrow><mrow><msup><mi>σ</mi><mn>4</mn></msup></mrow></mfrac><mo>−</mo><mn>3.</mn></mtd></mtr></mtable></math> When γ 1 = γ 2 =0, the value of Charlier’s distribution is equal to one. Values of these parameters for different actinides with their given atomic number are presented in Table 1 Table 1 [ 21 21 ]. 1 1 Z cn 1 1 90 1 1 91 1 1 92 1 1 93 1 1 94 1 1 95 1 1 γ 1 1 1 0 1 1 -0.08 1 1 0.36 1 1 0.23 1 1 0.3 1 1 0.38 1 1 γ 2 1 1 -0.36 1 1 -0.07 1 1 -0.34 1 1 -0.27 1 1 -0.3 1 1 -0.34 In the original approach of Gorodisskiy, the yield of heavy fragment is defined by the following relation <math display="block" id="M5"><mrow><msub><mi>Y</mi><mi>H</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>π</mi></mrow></msqrt></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>Y</mi><mi>s</mi></msub><mtext> </mtext><msup><mtext>e</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>2</mn><mtext> </mtext><msub><mi>u</mi><mi>s</mi></msub><msup><mrow/><mn>2</mn></msup></mrow></msup></mrow><mrow><msub><mi>σ</mi><mi>s</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mi>Y</mi><mi>a</mi></msub><mtext> </mtext><msup><mtext>e</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>2</mn><mtext> </mtext><msub><mi>u</mi><mi>a</mi></msub><msup><mrow/><mn>2</mn></msup></mrow></msup><mi>f</mi><mo stretchy="false">(</mo><mi>u</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>σ</mi><mi>a</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where Y s and Y a respectively correspond to the relative contributions of the symmetric and asymmetric modes of fission in fragment yield, respectively. Parameters γ 1 and γ 2 are adjusted with the experimental data for each compound nucleus formed after absorbtion of photon [ 21 21 ]. Also, u a , u s , σ a , and σ s are defined as follows: <math display="block" id="M6"><mrow><msub><mi>σ</mi><mtext>a</mtext></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mtext>pre</mtext></mrow></msub></mrow><mo>)</mo></mrow><mtext> </mtext><mrow><mo>(</mo><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><mn>73</mn></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mn>0.074</mn><mo>+</mo><mn>0.0296</mn><mroot><mi>E</mi><mn>4</mn></mroot></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac><mn>,</mn></mrow></math> <math display="block" id="M7"><mrow><msub><mi>σ</mi><mtext>s</mtext></msub><mo>=</mo><mfrac><mrow><mn>0.031</mn><mtext> </mtext><mrow><mo>(</mo><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mtext>pre</mtext></mrow></msub></mrow><mo>)</mo></mrow><mtext> </mtext><mroot><mi>E</mi><mn>4</mn></mroot></mrow><mrow><msqrt><mrow><mn>90.54</mn><mo>−</mo><mn>1.9</mn><mtext> </mtext><mfrac><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub><msup><mrow/><mn>2</mn></msup></mrow><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mtext>pre</mtext></mrow></msub></mrow></mfrac></mrow></msqrt></mrow></mfrac><mo>+</mo><mn>9.64,</mn></mrow></math> <math display="block" id="M8"><mrow><msub><mi>u</mi><mtext>s</mtext></msub><mo>=</mo><mfrac><mrow><mi>A</mi><mo>−</mo><mrow><mo>(</mo><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mtext>pre</mtext></mrow></msub></mrow><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow><mrow><msub><mi>σ</mi><mtext>s</mtext></msub></mrow></mfrac></mrow></math> and <math display="block" id="M9"><mrow><msub><mi>u</mi><mtext>a</mtext></msub><mo>=</mo><mfrac><mrow><mi>A</mi><mo>−</mo><mi>α</mi></mrow><mrow><msub><mi>σ</mi><mtext>a</mtext></msub></mrow></mfrac><mn>.</mn></mrow></math> where <math id="M10"><mrow><mi>α</mi><mo>=</mo><mn>54</mn><mfrac><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mtext>pre</mtext></mrow></msub></mrow><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac></mrow></math> for Z cn =90-91 and <math id="M11"><mrow><mi>α</mi><mo>=</mo><mn>28.6</mn><mfrac><mrow><msub><mi>A</mi><mrow><mi>c</mi><mi>n</mi></mrow></msub><mo>−</mo><msub><mi>υ</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi></mrow></msub></mrow><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac><mo>+</mo><mn>0.708</mn><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></math> for Z cn ≥ 92. A cn and Z cn are respectively the mass and charge numbers of compound nucleus. υ pre is the average pre-scission neutron multiplicity and E is the gamma-rays energy. Also, A is the mass number of the nascent fragment. In order to calculate the values of Y a and Y s in Eq. (3), the ratio <math id="M12"><mrow><mfrac><mrow><msub><mi>Y</mi><mtext>a</mtext></msub></mrow><mrow><msub><mi>Y</mi><mtext>s</mtext></msub></mrow></mfrac></mrow></math> is obtained through the fitting method with experimental data as follows: <math display="block" id="M13"><mrow><mfrac><mrow><msub><mi>Y</mi><mi>a</mi></msub></mrow><mrow><msub><mi>Y</mi><mi>s</mi></msub></mrow></mfrac><mo>=</mo><mn>1.244</mn><mtext> </mtext><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mtext>e</mtext><mrow><mo>−</mo><mn>0.0027</mn><mtext> </mtext><msup><mrow><mrow><mo>(</mo><mrow><mrow><mo>|</mo><mrow><mi>E</mi><mo>−</mo><mn>5.7</mn></mrow><mo>|</mo></mrow></mrow><mo>)</mo></mrow></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></msup></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mroot><mi>E</mi><mn>4</mn></mroot><mo>+</mo><mn>100</mn><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac><mo>−</mo><mn>0.4</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> and the values of Y a and Y s are obtained by replacing it in the above equations <math display="block" id="M14"><mrow><msub><mi>Y</mi><mtext>a</mtext></msub><mo>=</mo><mn>200</mn><mtext> </mtext><msup><mrow><mrow><mo>(</mo><mrow><mrow><mo>|</mo><mrow><mfrac><mrow><msub><mi>Y</mi><mtext>a</mtext></msub></mrow><mrow><msub><mi>Y</mi><mtext>s</mtext></msub></mrow></mfrac></mrow><mo>|</mo></mrow><mo>+</mo><mn>2</mn></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> and <math display="block" id="M15"><mrow><msub><mi>Y</mi><mtext>s</mtext></msub><mo>=</mo><mn>200</mn><mo>−</mo><mn>2</mn><msub><mi>Y</mi><mtext>a</mtext></msub><mn>.</mn></mrow></math> In a similar way, the yield of light conjugate fragment( Y L )is also obtained by replacing M L with A cn - A . Finally, the fission fragment mass yield is evaluated using the following equation <math display="block" id="M16"><mrow><mi>Y</mi><mo>=</mo><msub><mi>Y</mi><mtext>H</mtext></msub><mo>+</mo><msub><mi>Y</mi><mtext>L</mtext></msub><mn>.</mn></mrow></math> The Eq. (9) is revised by considering two adjustable parameters β and δ , instead of their fixed values in the original formalism of Gorodisskiy as, <math display="block" id="M17"><mrow><mfrac><mrow><msub><mi>Y</mi><mtext>a</mtext></msub></mrow><mrow><msub><mi>Y</mi><mtext>s</mtext></msub></mrow></mfrac><mo>=</mo><mn>1.244</mn><mtext> </mtext><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mtext>e</mtext><mrow><mo>−</mo><mn>0.0027</mn><mtext> </mtext><msup><mrow><mrow><mo>(</mo><mrow><mrow><mo>|</mo><mrow><mi>E</mi><mo>−</mo><mn>5.7</mn></mrow><mo>|</mo></mrow></mrow><mo>)</mo></mrow></mrow><mi>δ</mi></msup></mrow></msup></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mroot><mi>E</mi><mi>β</mi></mroot><mo>+</mo><mn>100</mn><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac><mo>−</mo><mn>0.4</mn></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> The parameters β and δ are obtained using the fitting method with experimental data. The developed approach is used to calculate the photofission fragment mass yield of various actinides at different energies. 3 Results and discussion 1 The revised formalism presented in this paper is applied to calculate the fragments mass distribution in photofission reactions of 232 Th at 16 MeV, 238 U at 17.2 MeV, 237 Np at 24 MeV, and 242 Pu at 20 MeV and compared with results of Gorodisskiy’s formalism as well as experimental data in Figs. 1 Figs. 1 , 2 2 , 3 3 , 4 4 . As it is clear from Figs. 1 Figs. 1 , 2 2 , 3 3 , 4 4 , the fission fragment mass yields of the present formalism are in better agreement with experimental data than Gorodisskiy’s original formalism. However there is some inconsistency between theoretical and experimental data especially in most probable fragmentation at maxima of fragment yields. In Figs. 5 Figs. 5 , 6 6 , 7 7 , and 8 8 the calculated results of fission fragment yields considering 1, 2, 3 and 4 pre-scission neutron multiplicity prior to fission are compared with experimental data for photofission reactions of 232 Th at 14 MeV, 238 U at 17.2 MeV, 237 Np at 24 MeV, and 242 Pu at 20 MeV. These figures indicate that the emission of neutrons prior to fission does not considerably affect the photofission fragment mass yield. This is one of the Gorodisskiy formalism drawbacks that is not able to properly indicate the role played by neutron emission in photofission. In order to calculate the fission fragment mass yield as a function of fragment atomic number, ( Z ), following the semi-empirical formula between fragments mass, ( A ) and its atomic number, ( Z ) is employed <math display="block" id="M18"><mrow><mi>A</mi><mo>=</mo><mfrac><mrow><msub><mi>A</mi><mrow><mtext>cn</mtext></mrow></msub></mrow><mrow><msub><mi>Z</mi><mrow><mtext>cn</mtext></mrow></msub></mrow></mfrac><mtext> </mtext><mrow><mo>(</mo><mrow><mi>Z</mi><mo>−</mo><mn>2.5</mn></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> The fragment mass yield of photofission for 238 U at 11.39 MeV and 13.39 MeV are calculated and indicated as a function of fragments atomic number in Figs. 9(a) and 9(b) Figs. 9(a) and 9(b) , respectively. This figure shows that the fragment mass yield of the present formalism for 238 U is better fitted with experimental data [ 30 30 ]. Generally, high-energy photons used in photofission studies, are generated through the Bremsstrahlung effect of accelerated electron beam crossing dense conversion targets. The maximum electron energy is called end-point energy and the average compound nucleus excitation energy corresponding to such end-point energy is an important factor to reproduce experimental data. Naik and his co-authors [ 29 29 ] determined the photofission fragment mass yield for the 240 Pu compound nucleus at 10 MeV photon end-point energy, which is equivalent to the average excitation energy equal to 7.61 MeV. In Fig. 10 Fig. 10 , the calculated results of fragment mass yield using the present formalism are compared with the experimental data. As it is clear from this figure, by considering the end-point energy instead of average excitation energy, the calculated fragment mass yield is better fitted with experimental data. The average mass of the light (〈 A L 〉) and heavy (〈 A H 〉) photofission fragments for 240 Pu and 232 Th actinides are also obtained using the following equation <math display="block" id="M19"><mrow><mo>〈</mo><msub><mi>A</mi><mtext>L</mtext></msub><mo>〉</mo><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mo stretchy="false">(</mo><mi>A</mi><mo>×</mo><msub><mi>Y</mi><mtext>L</mtext></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><msub><mi>Y</mi><mtext>L</mtext></msub></mrow></mfrac><mn>,</mn><mo>〈</mo><msub><mi>A</mi><mtext>H</mtext></msub><mo>〉</mo><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mo stretchy="false">(</mo><mi>A</mi><mo>×</mo><msub><mi>Y</mi><mtext>H</mtext></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><msub><mi>Y</mi><mtext>H</mtext></msub></mrow></mfrac><mn>,</mn></mrow></math> and the calculated results are compared respectively with experimental data of Thierens et al. [ 26 26 ] and Naik et al.[ 30 30 ] in Table 2 Table 2 . As it can be seen from Table 2 Table 2 , there is a fair agreement between results of our formalism and experimental data for 240 Pu isotope. A small inconsistency between results of the present formalism and experimental data [ 27 27 ] for 232 Th isotope is produced due to neutron emission through the fission process at high excitation energies. Photofission fragment mass yields for the 238 U isotope at 300 and 500 MeV are presented in Figs. 11(a) and 11(b) Figs. 11(a) and 11(b) . It is clear from these figures that, the calculated results of the present formalism are in good agreement with experimental data [ 28 28 ]. It should be noted that the original approach of Gorodisskiy was not successful in reproducing experimental data of high incident photons. 1 1   1 1 Our formalism 1 1 Thierens [ 26 26 ] 1 1 Our formalism 1 1 Naik [ 27 27 ] 1 1   1 1 240 Pu (30 MeV) 1 1   1 1 232 Th (80 MeV) 1 1 〈 A L 〉 1 1 100.57 1 1 100.29 1 1 92.93 1 1 91.74 1 1 〈 A H 〉 1 1 139.42 1 1 139.71 1 1 139.14 1 1 136.75 The peak to valley ratio of fission fragment yield for the 240 Pu compound nucleus is calculated and compared with experimental data of Naik [ 29 29 , 31 31 ] in Fig. 12 Fig. 12 . It is clearly indicated in this figure that the peak to valley ratio decreases as the Bremsstrahlung end-point photon energy grows. The photofission fragment mass yields for the 239 Pu isotope at 25 MeV incident photon energy and 234 U isotope at 5.77 and 6.11 MeV energies are compared with experimental data in Figs. 13 Figs. 13 and 14 14 . These Figures indicate satisfactory agreement between our results and experimental data except around maximum fragment mass yields where the results of the present formalism does not match with experimental data properly. 4 Conclusion 1 A phenomenological formula with some adjustable parameters is developed to calculate fission fragment mass yield for 232 Th, 234,238 U, 237 Np, and 239,240,242 Pu actinides at various energies. Calculated results are compared with results of Gorodisskiy’s formalism as well as experimental data. Satisfactory agreement has been achieved between results of the present formalism and experimental data especially at low and intermediate photon energies. In this research, the role played by neutron emission prior to fission has also been investigated. It has been stated earlier that, the present approach is not able to consider the effect of neutron multiplicity in photofission. The peak to valley ratio for photofission of the 240 Pu isotope is also obtained and indicated in Fig. 12 Fig. 12 . Figure 12 Figure 12 shows that, the values of peak to valley ratio is decreased with the increase of the Bremsstrahlung end-point energies. Generally, it has been shown that the calculated results of the present study are in good agreement with experimental data as compared with data reported in Ref. [ 29 29 ]. The average mass of light (〈 A L 〉) and heavy (〈 A H 〉) fragments for 232 Th and 240 Pu actinides are obtained and compared with experimental data [ 26 26 , 27 27 ]. Good agreement between the results of the present approach with experimental data has also been achieved. References W. Haxby , W. Shoupp , W. Stephens et al ., Photo-Fission of Uranium and Thorium , Phys. Rev . 59 57 ( 1941 ). doi: 10.1103/PhysRev.59.57 http://doi.org/10.1103/PhysRev.59.57 Photo-Fission of Uranium and Thorium C. Romano , Y. Danon , R. Block , J. Thompson , E. Blain , and E. Bond , Fission fragment mass and energy distributions as a function of incident neutron energy measured in a lead slowing-down spectrometer , Phys. Rev. C 81014607 ( 2010 ). doi: 10.1103/PhysRevC.81.014607 http://doi.org/10.1103/PhysRevC.81.014607 Fission fragment mass and energy distributions as a function of incident neutron energy measured in a lead slowing-down spectrometer D. M. Gorodisskiy , S. I. Mulgin , A. Ya. Rusanov and S. V. Zhdanov , Isotopic Invariance of Fission Fragment Charge Distributions for Actinide Nuclei at Excitation Energies above 10 MeV , Phys. At. Nucl . 66 1190 ( 2003 ). doi: 10.1134/1.1586436 http://doi.org/10.1134/1.1586436 Isotopic Invariance of Fission Fragment Charge Distributions for Actinide Nuclei at Excitation Energies above 10 MeV D. M. Gorodisskiy , S. I. Mulgin , V. N. Okolovich , A. Y. Rusanov , S. V. Zhdanov , Isotopic and isotonic effects in fission-fragment mass yields of actinide nuclei , Phys. Lett. B 548 , 45 ( 2002 ). doi: 10.1016/S0370-2693(02)02838-1 http://doi.org/10.1016/S0370-2693(02)02838-1 Isotopic and isotonic effects in fission-fragment mass yields of actinide nuclei W. Fu Cheng and H. Ji-Min , A study of the multimode fission model , Jour. of Phys. G: Nuclear and Particle physics 15 829 ( 1989 ). doi: 10.1088/0954-3899/15/6/013 http://doi.org/10.1088/0954-3899/15/6/013 A study of the multimode fission model P. Moller and J. Randrup , Calculated fission-fragment yield systematics in the region 74≤Z≤94 and 90≤N≤150 , Phys. Rev. C 91 044316 ( 2015 ). doi: 10.1103/PhysRevC.91.044316 http://doi.org/10.1103/PhysRevC.91.044316 Calculated fission-fragment yield systematics in the region 74≤Z≤94 and 90≤N≤150 M. R. Pahlavani and D. Naderi , Study of fusion cross-section in heavy-ion fusion-fission reactions at around fusion barrier energies using the Langevin dynamical approach , The Eur. Phy. J. A 48 129 ( 2012 ). doi: 10.1140/epja/i2012-12129-y http://doi.org/10.1140/epja/i2012-12129-y Study of fusion cross-section in heavy-ion fusion-fission reactions at around fusion barrier energies using the Langevin dynamical approach B. D. Wilkins and E. P. Steinberg , Semi-empirical interpretation of nuclear fission based on deformed-shell effects , Phys. Lett. B 42 141 ( 1972 ). doi: 10.1016/0370-2693(72)90046-9 http://doi.org/10.1016/0370-2693(72)90046-9 Semi-empirical interpretation of nuclear fission based on deformed-shell effects B. D. Wilkins , E. P. Steinberg and R. R. Chasman , Scission-point model of nuclear fission based on deformed-shell effects , Phys. Rev. C 14 1832 ( 1976 ). doi: 10.1103/PhysRevC.14.1832 http://doi.org/10.1103/PhysRevC.14.1832 Scission-point model of nuclear fission based on deformed-shell effects J. Moreau , K. Heyde , and M. Waroquier , Nuclear temperature effects in the scission-point model of nuclear fission , Phys. Rev. C 28 1640 ( 1983 ). doi: 10.1103/PhysRevC.28.1640 http://doi.org/10.1103/PhysRevC.28.1640 Nuclear temperature effects in the scission-point model of nuclear fission A. R. DeL. Musgrove , J. L. Cook and G. D. Trimble , in Proceedings Fifteenth Annual ACM , Bologna , V. II in IAEA-169 163 ( 1974 ). A. C. Wahl , Nuclear-charge distribution and delayed-neutron yields for thermal-neutron-induced fission of 235U, 233U, and 239Pu and for spontaneous fission of 252Cf , Atomic Data and Nuclear Data Tables 39 1 ( 1988 ). doi: 10.1016/0092-640X(88)90016-2 http://doi.org/10.1016/0092-640X(88)90016-2 Nuclear-charge distribution and delayed-neutron yields for thermal-neutron-induced fission of 235U, 233U, and 239Pu and for spontaneous fission of 252Cf J. Katakura , in Report JAERI-Research (2003-2004) . Report JAERI-Research (2003-2004) R. W. Mills , Fission product yield evaluation, Ph.D thessis , School of Physics and Space Research, University of Birmingham ( 1995 ). Fission product yield evaluation, Ph.D thessis K.-H. Schmidt , B. Jurado , C. Amouroux , C. Schmitt , General description of fission observables: GEF model code , Nucl. Data Sheets 131 107 ( 2016 ). doi: 10.1016/j.nds.2015.12.009 http://doi.org/10.1016/j.nds.2015.12.009 General description of fission observables: GEF model code M. R. Pahlavani and S. M. Mirfathi , Neutron-induced fission of even- and odd-mass plutonium isotopes within a four-dimensional Langevin framework , Phys. Rev. C 96 , 014606 ( 2017 ). doi: 10.1103/PhysRevC.96.014606 http://doi.org/10.1103/PhysRevC.96.014606 Neutron-induced fission of even- and odd-mass plutonium isotopes within a four-dimensional Langevin framework M. R. Pahlavani and S. M. Mirfathi , Dynamical simulation of neutron-induced fission of uranium isotopes using four-dimensional Langevin equations , Phys. Rev. C 93 , 044617 ( 2016 ). doi: 10.1103/PhysRevC.93.044617 http://doi.org/10.1103/PhysRevC.93.044617 Dynamical simulation of neutron-induced fission of uranium isotopes using four-dimensional Langevin equations M. R. Pahlavani and S. M. Mirfathi , Dynamics of neutron-induced fission of 235U using four-dimensional Langevin equations , Phys. Rev. C 92 , 024622 ( 2015 ). doi: 10.1103/PhysRevC.92.024622 http://doi.org/10.1103/PhysRevC.92.024622 Dynamics of neutron-induced fission of 235U using four-dimensional Langevin equations M. R. Pahlavani , S. M. Mirfathi , Dynamical simulation of neutron-induced fission of uranium isotopes using four-dimensional Langevin equations , Eur. Phys. J. A 52 95 ( 2016 ). doi: 10.1103/PhysRevC.93.044617 http://doi.org/10.1103/PhysRevC.93.044617 Dynamical simulation of neutron-induced fission of uranium isotopes using four-dimensional Langevin equations J. Randrup and P. Moller , Brownian Shape Motion on Five-Dimensional Potential-Energy Surfaces: Nuclear Fission-Fragment Mass Distributions . Phys. Rev. Lett . 106 132503 ( 2011 ). doi: 10.1103/PhysRevLett.106.132503; http://doi.org/10.1103/PhysRevLett.106.132503; Brownian Shape Motion on Five-Dimensional Potential-Energy Surfaces: Nuclear Fission-Fragment Mass Distributions J. Randrup , P. Moller and A. J. Sierk , Calculated fission yields of neutron-deficient mercury isotopes , Phys. Rev. C 85 , 024306 ,( 2012 ). doi: 10.1103/PhysRevC.85.024306 http://doi.org/10.1103/PhysRevC.85.024306 Calculated fission yields of neutron-deficient mercury isotopes D. M. Gorodisskiy , K. V. Kovalchuk , S. I. Mulgin et al ., Systematics of fragment mass yields from fission of actinide nuclei induced by the 5-200MeV protons and neutrons , Annals of Nuclear Energy 35 238 ( 2008 ). doi: 10.1016/j.anucene.2007.06.002 http://doi.org/10.1016/j.anucene.2007.06.002 Systematics of fragment mass yields from fission of actinide nuclei induced by the 5-200MeV protons and neutrons U. Brosa , S. Grossmann , A. Mller , Nuclear scission . Phys. Reports . 197 167 ( 1990 ). doi: 10.1016/0370-1573(90)90114-H http://doi.org/10.1016/0370-1573(90)90114-H . Nuclear scission D. M. Gorodisskiy , S. I. Mulgin , A. YA. Rusanov et al ., The new-revealed regularitieas of fragments mass yield from the proton induced fission of actinide nuclei ., Proceedings of the 5th International Conference on Dynamical Aspects of Nuclear Fission ( 2002 ): pp. 287 . doi: 10.1142/9789812776723-0023 http://doi.org/10.1142/9789812776723-0023 The new-revealed regularitieas of fragments mass yield from the proton induced fission of actinide nuclei S. I. Mulgin , V. N. Okolovich and S. V. Zhdanov , Observation of new channel in the proton-induced low-energy fission of nuclei from 233Pa to 245Bk , Phys. Lett. B 462 29 ( 1999 ). doi: 10.1016/S0370-2693(99)00859-X http://doi.org/10.1016/S0370-2693(99)00859-X Observation of new channel in the proton-induced low-energy fission of nuclei from 233Pa to 245Bk I. V. Pokrovsky , M. G. Itkis , J. M. Itkis et al ., Fission modes in the reaction 208Pb(18O,f) , Phys. Rev. C 62 014615 ( 2000 ). doi: 10.1103/PhysRevC.62.014615 http://doi.org/10.1103/PhysRevC.62.014615 Fission modes in the reaction 208Pb(18O,f) H. Thierens A. De. Clercq , E. Jacobs et al ., Kinetic energy and fragment mass distributions for 240Pu(s.f.), 239Pu(nth,f), and 240Pu(γ,f) , Phys. Rev. C 23 2104 ( 1981 ). doi: 10.1103/PhysRevC.23.2104 http://doi.org/10.1103/PhysRevC.23.2104 Kinetic energy and fragment mass distributions for 240Pu(s.f.), 239Pu(nth,f), and 240Pu(γ,f) H. Naik , G. N. Kim , R. Schwengner et al ., Fission product yield distribution in the 12, 14, and 16 MeV bremsstrahlung-induced fission of 232Th , The Eur. Phys. J. A 51 150 ( 2015 ). doi: 10.1140/epja/i2015-15150-8 http://doi.org/10.1140/epja/i2015-15150-8 Fission product yield distribution in the 12, 14, and 16 MeV bremsstrahlung-induced fission of 232Th B. Schroder , G. Nydahl and B. Forkman , High-energy photofission in 238U, 232Th and 209Bi , Nucl. Phys. A 143 449 ( 1970 ). doi: 10.1016/0375-9474(70)90541-5 http://doi.org/10.1016/0375-9474(70)90541-5 High-energy photofission in 238U, 232Th and 209Bi H. Naik , V. T. Nimje , D. Raj et al ., Mass distribution in the bremsstrahlung-induced fission of 232Th, 238U and 240Pu , Nucl. Phys. A 853 1 ( 2011 ). doi: 10.1016/j.nuclphysa.2011.01.009 http://doi.org/10.1016/j.nuclphysa.2011.01.009 Mass distribution in the bremsstrahlung-induced fission of 232Th, 238U and 240Pu H. Naik , F. Carrel , G. N. Kim et al ., Mass yield distributions of fission products from photo-fission of 238U induced by 11.5→17.3 MeV bremsstrahlung , Europ. Phys. Jour. A: Hadrons and Nuclei 49 94 ( 2013 ). doi: 10.1140/epja/i2013-13094-7 http://doi.org/10.1140/epja/i2013-13094-7 Mass yield distributions of fission products from photo-fission of 238U induced by 11.5→17.3 MeV bremsstrahlung S. S. Belyshev , B. S. Ishkhanov , A. A. Kuznetsov et al ., Mass yield distributions and fission modes in photofission of 238U below 20 MeV , Phys. Rev. C 91 034603 ( 2015 ). doi: 10.1103/PhysRevC.91.034603 http://doi.org/10.1103/PhysRevC.91.034603 Mass yield distributions and fission modes in photofission of 238U below 20 MeV E. Jacobs , H. Thierens , D. De. Frenne et al ., Product yields for the photofission of 238U with 12, 15, 20, 30, and 70 MeV bremsstrahlung , Phys. Rev. C 19 422 ( 1979 ). doi: 10.1103/PhysRevC.19.422 http://doi.org/10.1103/PhysRevC.19.422 Product yields for the photofission of 238U with 12, 15, 20, 30, and 70 MeV bremsstrahlung M. Ya. Kondrat’ko , V. N. Korinets and K. A. Petrzhak , The fragment yields in Np-237 photofission, Sov . Atom. Ene . 35 862 ( 1973 ). doi: 10.1007/BF01164117 http://doi.org/10.1007/BF01164117 The fragment yields in Np-237 photofission, Sov H. Thierens , E. Jacobs , P. D’hondt et al ., Fragment mass and kinetic energy distributions for 242Pu(sf), 241Pu(nth,f), and 242Pu(γ,f) , Phys. Rev. C 29 498 ( 1984 ). doi: 10.1103/PhysRevC.29.498 http://doi.org/10.1103/PhysRevC.29.498 Fragment mass and kinetic energy distributions for 242Pu(sf), 241Pu(nth,f), and 242Pu(γ,f) M. Ya. Kondratko A.V. Mosesov , K. A. Petrzhak et al ., Product yields of photofission of Pu-239, Sov . Atom. Ene . 50 , 41 ( 1981 ). doi: 10.1007/BF01141251 http://doi.org/10.1007/BF01141251 Product yields of photofission of Pu-239, Sov P. David , J. Debrus , U. Kim et al ., High-energy photofission of gold and uranium , Nucl. Phys. A 197 163 ( 1972 ). doi: 10.1016/0375-9474(72)90753-1 http://doi.org/10.1016/0375-9474(72)90753-1 High-energy photofission of gold and uranium Vo Dak Bang , Yu. S. Zamyatnin , Chan Dyk Tkhien et al ., Yield of Pu-242 photofission fragments, Sov . Atom. Ene . 58 , 320 ( 1985 ). doi: 10.1007/BF01207227 http://doi.org/10.1007/BF01207227 Yield of Pu-242 photofission fragments, Sov A. Gk , M. Chernykh , C. Eckardt et al ., Fragment characteristics from fission of 238U and 234U induced by 6.5 to 9.0 MeV bremsstrahlung , Nucl. Phys. A 851 , 1 ( 2011 ). doi: 10.1016/j.nuclphysa.2010.12.012 http://doi.org/10.1016/j.nuclphysa.2010.12.012 Fragment characteristics from fission of 238U and 234U induced by 6.5 to 9.0 MeV bremsstrahlung 10.1007/s41365-018-0482-1 2018-01-02 2018-05-15 2018-06-04 2018-08-31 2018-10-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0482-1 http://dx.doi.org/10.1007/s41365-018-0482-1 . 2018 Nuclear Science and Techniques 10 29