Simulation of Displacement Damage in CsPbBr 3 Induced by Neutron Irradiation Based on Monte Carlo Method

Abstract

Cesium lead bromide (CsPbBr3) is a fully inorganic halide perovskite material known for its excellent photoelectric properties, offering significant advantages for applications in aerospace and nuclear fields. To evaluate its radiation hardness under neutron exposure, the transport process of 1-14 MeV neutrons in CsPbBr3 was simulated using the Geant4 Monte Carlo toolkit. This study focuses on the primary damage characteristics, systematically analyzing the primary knock-on atom (PKA) spectrum and non-ionizing energy loss (NIEL). The simulation results indicate that most PKAs are distributed in the low-energy range. As the incident neutron energy increases, PKA types become more diverse, introducing transmutation products such as [78][79][80][81] Se, 133 Xe, 76 As, and 205 Hg. Furthermore, a distinct anomaly is observed at lower neutron energies (~3 MeV), where the displacement of Pb atoms exhibits a localized peak directly attributed to its prominent (n, n) elastic scattering resonance, although the overall macroscopic damage remains heavily dominated by the Br sublattice. Crucially, the calculated NIEL, number of displaced atoms (Nd), and displacements per atom (dpa) exhibit a non-monotonic dependence on incident neutron energy, initially increasing and then decreasing beyond ~10 MeV. This trend is primarily driven by the transition from elastic to inelastic scattering dominance, coupled with increased ionizing energy partitioning at higher PKA energies. This paper provides fundamental data on the primary damage state of CsPbBr3, establishing an essential source term basis for subsequent multiscale simulations of defect evolution.