In nuclear collisions induced by stable or radioactive neutron-rich nuclei a transient state of nuclear matter with an appreciable isospin asymmetry as well as thermal and compressional excitation can be created. This offers the possibility to study the properties of nuclear matter in the region between symmetric nuclear matter and pure neutron matter. In this review, we discuss recent theoretical studies of the equation of state of isospin-asymmetric nuclear matter and its relations to the properties of neutron stars and radioactive nuclei. Chemical and mechanical instabilities as well as the liquid-gas phase transition in asymmetric nuclear matter are investigated. The in-medium nucleon-nucleon cross sections at different isospin states are reviewed as they affect significantly the dynamics of heavy ion collisions induced by radioactive beams. We then discuss an isospin-dependent transport model, which includes different mean-field potentials and cross sections for the proton and neutron, and its application to these reactions. Furthermore, we review the comparisons between theoretical predictions and available experimental data. In particular, we discuss the study of nuclear stopping in terms of isospin equilibration, the dependence of nuclear collective flow and balance energy on the isospin-dependent nuclear equation of state and cross sections, the isospin dependence of total nuclear reaction cross sections, and the role of isospin in preequilibrium nucleon emissions and subthreshold pion production.

Simulations by transport codes are indispensable to extract valuable physical information from heavy-ion collisions. In order to understand the origins of discrepancies among different widely used transport codes, we compare 15 such codes under controlled conditions of a system confined to a box with periodic boundary, initialized with Fermi-Dirac distributions at saturation density and temperatures of either 0 or 5 MeV. In such calculations, one is able to check separately the different ingredients of a transport code. In this second publication of the code evaluation project, we only consider the two-body collision term; i.e., we perform cascade calculations. When the Pauli blocking is artificially suppressed, the collision rates are found to be consistent for most codes (to within 1% or better) with analytical results, or completely controlled results of a basic cascade code. PHYSICAL REVIEW C 97, 034625 (2018) to reach that goal, it was necessary to eliminate correlations within the same pair of colliding particles that can be present depending on the adopted collision prescription. In calculations with active Pauli blocking, the blocking probability was found to deviate from the expected reference values. The reason is found in substantial phase-space fluctuations and smearing tied to numerical algorithms and model assumptions in the representation of phase space. This results in the reduction of the blocking probability in most transport codes, so that the simulated system gradually evolves away from the Fermi-Dirac toward a Boltzmann distribution. Since the numerical fluctuations are weaker in the Boltzmann-Uehling-Uhlenbeck codes, the Fermi-Dirac statistics is maintained there for a longer time than in the quantum molecular dynamics codes. As a result of this investigation, we are able to make judgements about the most effective strategies in transport simulations for determining the collision probabilities and the Pauli blocking. Investigation in a similar vein of other ingredients in transport calculations, like the mean-field propagation or the production of nucleon resonances and mesons, will be discussed in the future publications.

Transport simulations are very valuable for extracting physics information from heavy-ion collision experiments. With the emergence of many different transport codes in recent years, it becomes important to estimate their robustness in extracting physics information from experiments. We report on the results of a transport code comparison project. 18 commonly used transport codes were included in this comparison: 9 Boltzmann-Uehling-Uhlenbeck-type codes and 9 Quantum-MolecularDynamics-type codes. These codes have been required to simulate Au+Au collisions using the same physics input for mean fields and for in-medium nucleon-nucleon cross sections, as well as the same initialization set-up, the impact parameter, and other calculational parameters at 100 and 400 AMeV incident energy. Among the codes we compare one-body observables such as rapidity and transverse flow distributions. We also monitor non-observables such as the initialization of the internal states of colliding nuclei and their stability, the collision rates and the Pauli blocking. We find that not completely identical initializations constitute partly for different evolutions. Different strategies to determine the collision probabilities, and to enforce the Pauli blocking, also produce considerably different results. There is a substantial spread in the predictions for the observables, which is much smaller at the higher incident energy. We quantify the uncertainties in the collective flow resulting from the simulation alone as about 30% at 100 AMeV and 13% at 400 AMeV, respectively. We propose further steps within the code comparison project to test the different aspects of transport simulations in a box calculation of infinite nuclear matter. This should, in particular, improve the robustness of transport model predictions at lower incident energies where abundant amounts of data are available.

No abstract