Improved initial colliding nuclei density profile method for QMD-type transport models

Accurate modeling of the initial density profile is essential for studying heavy-ion collisions (HICs) with a transport model. Within the framework of the quantum molecular dynamics (QMD)-type model, a novel method for generating nuclear density distributions based on a Fourier series expansion (FSE) is proposed. In this approach, the objective density distribution is expanded into a Fourier series to construct a filter function, which is then applied to select the randomly sampled nucleon coordinates in phase space to generate a three-dimensional nuclear density distribution that matches the desired profile. This new initialization method is further incorporated into the ultrarelativistic quantum molecular dynamics (UrQMD) model, and the bubble-like density distribution of $^{96}$Ru is constructed, showing good stability. Then, by simulating $^{96}$Ru+$^{96}$Ru collisions at $E_\mathrm{lab}=1500$ MeV/nucleon with different equations of state (EoS) and initialization methods, the effects of the initialization method on the final state observables and the constrained information of EoS are analyzed. It is found that the maximum system density increases when the new FSE initialization method is adopted, and results in an enhanced collective flow. Moreover, a relatively stiff EoS with $K_0>280$ MeV is favored when adopting the Woods-Saxon initialization method, whereas an EoS of $K_0$ = 200-280 MeV is supported when using the FSE initialization. These results indicate that, within QMD-like transport models, the FSE filtering method provides a reliable means to sample nuclei with exotic density profiles, offering new insight to investigate nuclear structure and dense nuclear matter EoS through HICs.