Background: The distribution of protons and neutrons in the matter created in heavy-ion collisions is one of the main points of interest for the collision physics, especially at supranormal densities. We report results of the first systematic simulation of proton and neutron density distributions in central heavy-ion collisions within the beam energy range of E beam ≤ 800 MeV/nucl. Method: The Boltzmann-Uhlenbeck-Uehling (pBUU) transport model with four empirical models for the density dependence of the symmetry energy was employed. Results of simulations using pure Vlasov dynamics were added for completeness. In addition, the Time Dependent Hartree Fock (TDHF) model, with the SV-bas Skyrme interaction, was used to model the heavy-ion collisions at E beam ≤ 40 MeV/nucl. Maximum proton and neutron densities ρ max p and ρ max n , reached in the course of a collision, were determined from the time evolution of ρ p and ρ n .Results: The highest total densities predicted at E beam = 800 MeV/nucl were of the order of ∼ 2.5 ρ 0 (ρ 0 = 0.16 fm −3 ) for both Sn and Ca systems. They were found to be only weakly dependent on the initial conditions, beam energy, system size and a model of the symmetry energy.The proton-neutron asymmetry δ = (ρ max n −ρ max p )/(ρ max n +ρ max p ) at maximum density does depend, though, on these parameters. The highest value of δ found in all systems and at all investigated beam energies was ∼ 0.17.
Summary:We report the first systematic simulation of proton and neutron densities in matter produced in heavy-ion collisions of Ca and Sn systems at beam energies below 800 MeV/nucl using pBUU and TDHF models. We find limits on the maximum proton and neutron densities and the related proton-neutron asymmetry δ as a function of the initial state, beam energy, system size and a symmetry energy model. While the maximum densities are almost independent of these parameters, our simulation reveals, for the first time, their subtle impact on the protonneutron asymmetry. Most importantly, we find that variations in the proton-neutron asymmetry at maximum densities are related at most at 50% level to the details in the symmetry energy at supranormal density. The reminder is due to the details in the symmetry energy at subnormal densities and its impact on proton and neutron distributions in the initial state. This result puts to forefront the need of a proper initialization of the nuclei in the simulation, but also brings up the question of microscopy, such as shell effects, that affect initial proton and neutron densities, 2 but cannot be consistently incorporated into semiclassical transport models.