Quantum molecular dynamics model based on relativistic mean field theory for light nucleus fragmentation in hadron therapy

This study evaluates the accuracy of nuclear fragmentation simulations using a quantum molecular dynamics (QMD) model based on relativistic mean field (RMF) theory for an energy range of 50-400 MeV/u, relevant to hadron therapy. A total of 16 parameter sets within the RMF framework are assessed based on their ability to reproduce ground-state properties such as the mean squared radius and binding energy, as obtained in QMD simulations. Among these, the NS2 parameter set is identified as the most suitable for describing stable nuclei over a wide mass range, with the use of an adaptive Gaussian wave packet width. Fragmentation cross sections of carbon ion projectiles on light nuclei targets (H, C, O, Al, Ti, and Cu) are simulated at incident energies of 50, 95, 290, and 400 MeV/u and compared with experimental data. The results indicate that the RQMD.RMF model provides superior reproductions for fragmentation at lower energies (50 and 95 MeV/u) compared to the Light Ion QMD (LIQMD) model implemented in Geant4 version 11.2. At higher energies (290 and 400 MeV/u), the RQMD.RMF model performs comparably to the LIQMD. This study demonstrates that the RQMD.RMF model provides a reliable framework for analyzing nuclear fragmentation and holds potential for applications in the planning and quality assurance of hadron therapy.