A general framework for describing photofission observables of actinides at an average excitation energy below 30 MeV

A reasonable prediction of photofission observables plays a paramount role in understanding the photofission process and guiding various photofission-induced applications, such as short-lived isotope production, nuclear waste disposal, and nuclear safeguards. However, the available experimental data for photofission observables are limited, and the existing models and programs have mainly been developed for neutron-induced fission processes. In this study, a general framework is proposed for characterizing the photofission observables of actinides, including the mass yield distributions (MYD) and isobaric charge distributions (ICD) of fission fragments and the multiplicity and energy distributions of prompt neutrons (n_p) and prompt γ rays (γ_p). The framework encompasses various systematic neutron models and empirical models considering the Bohr hypothesis and does not rely on the experimental data as input. These models are then validated individually against experimental data at an average excitation energy below 30 MeV, which shows the reliability and robustness of the general framework. Finally, we employ this framework to predict the characteristics of photofission fragments and the emissions of prompt particles for typical actinides including ²³²Th, ^{235, 238}U and ²⁴⁰Pu. It is found that the ²³⁸U(γ, f) reaction is more suitable for producing neutron-rich nuclei compared to the ²³²Th(γ, f) reaction. In addition, the average multiplicity number of both n_p and γ_p increases with the average excitation energy.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.