Available online at www.sciencedirect.com ScienceDirect Physics Procedia 47 (2013 ) 53 59 Scientific Worshop on Nuclear Fission Dynamics and the Emission of Prompt Neutrons and Gamma Rays, Biarritz, France, 28-30 November 2012 Energy partition in fission within a microscopic approach M. Mirea Horia Hulubei National Institute for Physics and Nuclear Engineering, P.O. Box MG-6, Bucharest, Romania Abstract The intrinsic excitation energy of fission fragments is dynamically evaluated in terms of the time-dependent pairing equations. These equations are corroborated with two conditions. One of them fixes the number of particles and the other separates the pairing active spaces associated to the two fragments in the vicinity of the scission configuration. The excitation energy in a wide distribution of fission fragments is calculated for the 234 U parent nucleus. c 2013 The TheAuthors. Published by by Elsevier Elsevier B.V. B.V. Open access under CC BY-NC-ND license. Selection and peer-review under responsibility of of Joint Joint Research Centre Centre - Institute - Institute for Reference for Reference Materials Materials and Measurements and Measurements. Keywords: Fission, Excitation energy, 234 U PACS: 24.75.+i, 21.60.Cs 1. Introduction Under the action of a mutual Coulomb repulsion, at scission the fission fragments are accelerated in opposite directions. These fragments are highly excited. The maximal inetic energy issued in the process amounts to the Q-value in the case of cold fission. The fragments decay on their ground states mainly by evaporation of neutrons and by radiation emission. It is nown that the motion of any physical system is governed by conservative forces and by frictional ones that give rise to dissipation. Consequently, the excitation energy of the fragments must depend on the dynamics of the nuclear system in its path to scission. In this contribution a microscopic model for the intrinsic energy partition between two complementary fission fragments in terms of the time-dependent pairing equations (TDPE) (Mirea, 2011, 2012) is presented. These equations are corroborated with two conditions. One of them fixes the number of particles and the other separates the pairing active spaces associated to the two fragments in the vicinity of the scission configuration. It is shown that the available intrinsic dissipated energy is not shared proportionally to the masses of the two fission fragments (Gönnenwein, 2003). If the heavy fragment possesses nucleon numbers close to the magic ones, the accumulated intrinsic excitation energy is lower than that of the light fragment. Recently, this problem was intensively investigated within a wide range of methods, including Hartree-Foc theories (Younes amd Gogny, 2011), the statistical mechanics (Schmidt and Jurado, 2010, 2011), the phenomenological point by point model (Tudora et al., 2012; Manailescu et al., 2011), empirical evaluations (Young-Jing et al,, 2012), or the shell model in the sudden approximation (Carjan et al., 2012). Corresponding author. Tel.: +4-021-404-2300 ; fax: +4-021-457-4440. E-mail address: mirea@nipne.ro 1875-3892 2013 The Authors. Published by Elsevier B.V. Open access under CC BY-NC-ND license. Selection and peer-review under responsibility of Joint Research Centre - Institute for Reference Materials and Measurements doi:10.1016/j.phpro.2013.06.009
54 M. Mirea / Physics Procedia 47 ( 2013 ) 53 59 2. Model First of all, the main steps in deriving the microscopic equations of motion constrained by the projections of number of particles are presented. The starting point of the model is the energy functional L = ϕ H i t λ(n 1 ˆN 2 N 2 ˆN 1 ) ϕ (1) where H(t)= ε (t)(a + a + a + a ) G a + a+ a iaī (2) >0,i>0 is the Hamiltonian with pairing residual interactions, ε being single particle energies, ( ) ϕ = u + v a + a+ 0 (3) is the Bogoliubov wave function with BCS occupation v and vacancy u amplitudes. λ represents a Lagrange multiplier and Nˆ 1 = (a 1 a + 1 + a 1 a + ), 1 1 Nˆ 2 = (a 2 a + 2 + a 2 a + ) (4) 2 2 are operators for the number of particles in the pairing active level space pertaining to the nuclei A 1 and A 2, respectively. The summation index, 1 and 2 run over all states considered in the pairing active space, the states of the fragment one, and the those of the fragments two in the same space, respectively. N 1 and N 2 are the correspondent number of particles belonging to the two fission fragments. Obviously, the identity ϕ (N 1 ˆN 2 N 2 ˆN 1 ) ϕ = 0 can be used to fix the number of particle N 1 and N 2 in the two fission fragments. The condition is added to the energy functional (1) by mean of the Lagrange multiplier as an implicit condition. Therefore, the equation (1) is subject to a supplementary condition that project the number of particles onto the two fragments. From Eq. (1), the next system is obtained for the TDPE, eventually: i ρ 1 = κ 1 Δ 1 κ 1 Δ 1, i ρ 2 = κ 2 Δ 2 κ 2 Δ 2, i κ 1 =(2ρ 1 1)Δ 1 2κ 1 (ε 1 + λn 2 ), (5) i κ 2 =(2ρ 2 1)Δ 2 2κ 2 (ε 2 λn 1 ) where ρ = v 2 are occupation probabilities, κ = u v are pairing moment components, and Δ = G κ is the pairing gap. Δ 1 = G 1 1 κ 1 + G 12 2 κ 2 and Δ 2 = G 12 1 κ 1 + G 2 2 κ 2 are the gaps for the two fragments. It must be noticed that if G 12 =0, then we obtain two sets of non coupled equations, one set for each fission fragment (Mirea, 2011). This last condition separates the pairing active space into two sub-spaces. A particular form of these equations was deduced for the first time in Refs. (Koonin and Nix, 1976; Bloci and Flocard, 1976). They were developed in the last years (Mirea, 2012b; Mirea et al., 2004) to tae into account the Landau-Zener effect (Mirea, 2008, 2003), the pair breaing mechanism (Mirea, 2009), and the dissipation in the effective mass (Mirea, 2011b, 2010). The time dependent pairing equations are solved simply in terms of a monopole pairing force with strengths determined by a renormalization of the BCS equations as function of the number of levels. Recently, a realistic treatment was made with a density functional for the pairing interaction (Avez et al., 2008; Ebata et al., 2010). These TDPE offer a measure of the average dissipated energy at a given deformation during the disintegration process provided that the velocity of the deformations are nown. The difference between the total energy value E obtained within the TDPE and E 0 given by the static BCS-equations represents an approximate measure for the dissipation E : E = E E 0. (6)
M. Mirea / Physics Procedia 47 (2013 ) 53 59 55 E is expressed simply in terms of ρ and κ E = 2 ε ρ G κ 2 G ρ 2, (7) and E 0 corresponds to the values ρ 0 and κ0 associated to the lower energy state, that is, obtained from BCS equations. This definition was introduced in Ref. (Koonin and Nix, 1976). So, the TDPE provide the values of ρ and κ in each fragment. By calculating the BCS solutions of the same fragment in the average deformation, in the corresponding pairing active space, the intrinsic excitation is obtained with formula (6). 3. Results The low energy fission of 234 U will be investigated. In the Macroscopic-Microscopic Method (MMM), the whole system is characterized by some collective coordinates that vary in time leading to a split of the nucleus into two fragments of masses A 1 and A 2. The macroscopic-microscopic model provides a single particle potential in which the nucleons move independently. The basic ingredient of the MMM is the nuclear shape parametrization. In our wor, an axial symmetric nuclear shape is obtained by smoothly joining two spheroids of semi-axis a i and b i (i=1,2) with a nec surface generated by the rotation of an arc of circle. By imposing the condition of volume conservation we are left with five independent generalized coordinates that are associated to five degrees of freedom: the elongation R denoting the distance between the centers of the fragments; the necing parameter C = S/R 3 related to the curvature of the nec; the eccentricities ε i =[1 (b i /a i ) 2 ] 1/2 if a i b i (or ε i = [1 (a i /b i ) 2 ] 1/2 if a i < b i ) associated to deformations of the two fragments (a i and b i denoting the semi-axis) and the mass asymmetry parameter η = a 1 /a 2 given by the ratio of the major semi-axis. This parametrization was widely used previously in the calculation concerning cluster decay (Mirea et al., 2011, 2012, 2011b; Mirea, 2012c), barriers for fission (Mirea et al., 2010; Mirea and Tassan-Got, 2011; Companis et al., 2011) or superheavy elements synthesis (Mirea et al., 2009). A Woods-Saxon two center shell model (Mirea, 2008) compatible with this nuclear shape parametrization is used to compute the single particle levels and to disentangle the wave functions before that the scission is reached (Mirea, 2011). The Woods-Saxon potential, the Coulomb interaction and the spin-orbit term are diagonalized in the orthogonal wave function basis of the semi-symmetric two center oscillator (Mirea, 1996, 2006). Using the minimal action principle (Hill et al., 1972) a trajectory in the configuration space beginning with the parent ground state and ending in the exit point of the barrier could be determined. In order to minimize the action integral, a numerical procedure able to determine the dynamical fission path is used. This procedure was widely used in fission calculations in the past (Mirea and Tassan-Got, 2011). For the effective mass, the semi-adiabatic craning approach is used (Mirea, 2011b, 2010). This model gives values of the inertia close to that obtained within the Gaussian overlap method (Gozdz et al., 1985). The fission barrier along the minimal action path is presented in Fig. 1. The heights are in good agreement with experimental data. As remared in Ref. (Randrup and Möller, 2011), it is not yet understood how the compound nucleus is transformed in a variety of different fragmentations. It is also believed that the models for mass distributions have limited predictive power. To overcome these difficulties, some simple assumptions are made in our wor. As in Ref. (Mirea et al., 2010), the mass distribution of the fragments was relatively well reproduced by considering that the variation of the mass asymmetry is linear from the saddle configuration of the outer barrier up to the exit point. These partitions are selected from the maximal yields given in Ref. (Wahl, 1988). The action integral gives the probabilities to obtain these partitions and therefore the relative yields can be evaluated. The relative yields are plotted in Fig. 2. The most important behavior of the mass distribution is reproduced: the distribution exhibits two clearly separated symmetric mountains and the heavy pea is centered at A 1 134, while the experimental value is approximately 136 (Nishio et al., 1988). Another prominent feature of the mass distribution is an oscillation in yield. A such phenomenon was evidenced experimentally (Grant, 1976). The ingredients needed to obtain the dissipated energy are the single particle energies and the internuclear velocity. These energies are calculated within the Woods-Saxon two center model and the internuclear velocity is considered to be 10 6 m/s. This velocity gives a tunneling time of 10 21 s. A similar value resorts in calculations devoted to alpha (Mirea, 2001) and cluster decays (Mirea, 1998, 1999, 2003). This internuclear velocity is in agreement with the results obtained in Ref. (Serot et al., 1994). In the previous reference a parent metastable state is generated and a tunneling time of the order of (10 21-10 22 ) s is computed.
56 M. Mirea / Physics Procedia 47 ( 2013 ) 53 59 Fig. 1. 234 U fission barrier V for a final partition 102 Zr+ 132 Te determined along the minimal action trajectory. Some particular shapes related to the ground state, the extremes of the barrier, the exit point and the scission point are inserted in the plot. The distances for the elongation R that characterizes the shapes are 4.17, 7.7, 10.43, 12.64, 15.53, 17.53 and 20.2 fm. Fig. 2. Relative yields Y in a logarithmic scale as function of the masses of the fragments A 1 and A 2.
M. Mirea / Physics Procedia 47 (2013 ) 53 59 57 Fig. 3. The total excitation energies E of the fission fragments are plotted versus their masses with a thic full line. A thin full line is used for the deformation energy, a thin dashed line gives the neutron dissipated energy while the thin dot dashed line gives the dissipated energy for the proton space. The sum of excitation energies of the fission partners is plotted as function of the heavy fragment mass with a thic dashed line. Having the dynamical occupation probabilities together with the BCS solutions and the level scheme for each nucleus, we evaluated the excitation energy for all individual fragment. It is not possible to extract directly the energy partition from experimental data. Nevertheless, the main de-excitation process is the neutron evaporation. Therefore, indirect information can be obtained from neutron multiplicities, for which accurate results are available (Nishio et al., 1988). The excitation energy of each fragment is computed within the relation (6). The results are plotted in Fig. 3. Several experimental features are reproduced by the theoretical data. The deeply minimum in the neutron multiplicity occurs close to the mass of the doubly magic nucleus 132 Sn. A maximal value of the neutron multiplicity is obtained for the mass 116, that is in the symmetric fission region. In general, the excitation energy of the light fragment is larger than that of the heavy one. We can compare also the behavior of the average neutron multiplicity with the sum of total excitation energies of partners. The average neutron multiplicity is defined as the sum of neutrons from both fragments (Nishio et al., 1988). The minimal value of the average neutron multiplicity appears around 130 u. This fact is reproduced by theoretical data. Our results plotted with a thic dashed line in Fig. 3 give a minimum of the excitation. An enhanced neutron emission is obtained in the symmetric fission and the very asymmetric region. Between these extreme, the total neutron emission shows a structure with a maximal value around A 1 =140. These behavior are also exhibited by the total excitation energy (Nishio et al., 1988). From another respect, it must be mentioned that the excitation energy exhibits a structure with strong fluctuations, especially for A 2 < 82. A possible explanation for this peculiar aspect could be related to the existence of another fission path associated to large mass asymmetries. Therefore, the deformations and the dynamics assumed in this region could be not realistic. In this description, the deformations of the complementary fragments is intimately related to their excitations. Just before scission, the single particle level scheme of the compound nucleus could be considered as the overlap of the level schemes of the two fragments. If the pairing residual interaction is neglected, the levels are filled with nucleons up to the Fermi energy in the compound nucleus. For a zero excitation energy, the levels schemes of the two fragments must be also filled up to the Fermi their energies. But the distributions of the levels in the two fragments are strongly dependent of deformations. Only for particular deformations the levels of both fragments could be filled up to the Fermi energies. If some deviations from these optimal deformations are produced, an inherent single particle excitation must be produced. In the treatment that involves the dynamical projection of the number of particles in term of the time dependent pairing equations, the excitations are obtained by calculating the occupation probabilities
58 M. Mirea / Physics Procedia 47 ( 2013 ) 53 59 of the single particle states in the common active pairing space, that is for a finite number of levels above and below the Fermi energy. When the scission is reached, the levels of the two separated fragments are superimposed. The dynamical occupation probabilities of the levels depend on the way in which the levels of complementary fragments are interlaced and bounced in shells. On another hand, the ground state BCS probabilities in each nucleus depends on the succession and the energies of the levels. These considerations emphasize the role of the deformations. If the deformation of one of the nucleus is changed, the levels change their energies leading to a modification of the intermission with the single particle states of the complementary fragment, and even the number of levels in the common pairing space of the two fragments could not be the same. Due to these properties, for some deformations, it is even impossible to find solutions of the equations or the dissipated energy found in one fragment is very large. Acnowledgement This wor was supported by CNCS-UEFISCDI, project number PN-II-ID-PCE-2011-3-0068. References Avez, B., Simenel, C., Chomaz, Ph., 2008. Physical Review C 78, 044318. Bloci, J., Flocard, H., 1976. Nuclear Physics A 273, 45. Brac, M., Damgaard, J., Jensen, A.S., Pauli, H.C., Strutinsy, V.M., Wong, C.Y., 1972. Review of Modern Physics 44, 320. Carjan, N., Hambsch, F.-J., Rizea M., Serot, O., 2012. Physical Review C 85, 044601. Companis, I., Mirea, M., Isbasescu, A., 2011. Romanian Journal of Physics 56, 63. Ebata, S., Naatsuasa, T., Inaura, T., Yoshida, K., Hashimoto, Y., Yabana, K., 2010. Physical Review C 82, 034306. Gönnenwein, F., 2004. Neutrons from Fission, Seminar on Fission: Pont d Oye, Habay la Neuve Belgium, September 16-19 2003, World Scientific page 113. Gozdz, A., Pomorsi, K., Brac, M., Werner, E., 1985. Nuclear Physics A 442, 26. Grant, I.S., 1976. Report Progress in Physics 39, 955. Koonin, S.E., Nix J.R., 1976. Physical Review C 13, 209. Manailescu, C., Tudora, A., Hambsch, F.-J., Morariu, C., Oberstedt, S., 2011. Nuclear Physics A 867, 12. Mirea, M., 1998. Physical Review C 57, 2484. Mirea, M., 1999. European Physical Journal A 4, 335. Mirea, M., 2001. Physical Review C 63, 034603. Mirea, M., 2003. Modern Physics Letters A 18, 1809. Mirea, M., Tassan-Got, L., Stephan, C., Bacri, C.O., 2004. Nuclear. Physics. A 735, 21. Mirea, M., 2008. Physical Review C 78, 044618. Mirea, M., 2009. Physics Letters B 680, 316. Mirea, M., Bobulescu, R.C., 2010. Journal of Physics G 37, 055106. Mirea, M., 2011. Physical Review C 83, 054608. Mirea, M., 2011b. Romanian Reports in Physics 63, 676. Mirea, M., 2012. Physics Letters B 717, 252. Mirea, M., 2012b. International Journal of Modern Physics E 21, 1250035. Mirea, M., Delion, D.S., Sandulescu, A., 2009. EPL 85, 12001. Mirea, M., Sandulescu, A., Delion, D.S., 2011. Nuclear Physics A 870-871, 23. Mirea, M., Sandulescu, A., Delion, D.S., 2011b. Proceedings of the Romanian Academy Series A 12, 203. Mirea, M., Sandulescu, A., Delion, D.S., 2012. European Physical Journal A 48, 86. Mirea, M., 2012c, Romanian Journal of Physics 57, 372. Mirea, M., Delion, D.S., Sandulescu, A., 2010. Physical Review C 81, 044317. Mirea, M., Tassan-Got, L., 2011. Central European Journal of Physics 9, 116. Mirea, M., 1996. Physical Review C 54, 302. Mirea, M., 2006. Nuclear Physics A 780, 13. Nishio, K., Naashima, M., Kimura, I., Naagome, Y., 1988. Journal of Nuclear Science and Technologies 35, 631. Morariu, C., Tudora, A., Hambsch, F.J., Oberstedt, S., Manailescu, C., 2012. Journal of Physics G 39, 055103. Randrup, J., Möller, P., 2011. Physical Review Letters 106, 132503. Schmidt, K.-H., Jurado, B, 2010. Physical Review Letters 104, 212501. Schmidt, K.-H., Jurado, B., 2011. Physical Review C 83, 014607. Serot, O., Carjan N., Strottman, D., 1994. Nuclear Physics A 569, 562. Wahl, A.C., 1988. Atomic Data and Nuclear Data Tables 39, 1.
M. Mirea / Physics Procedia 47 (2013 ) 53 59 59 Younes, W., Gogny, D., 2011. Physical Review Letters 107, 132501. Yong-Jing, C., Jing, Q., Ting-Jin, L., Zhu-Xia, L., Xi-Zhen, W., Neng-Chuan, S., 2012. International Journal of Modern Physics E 21, 1250073.