Search NASA⌕ Search

Engineering topics

Danielewicz, P.

Publications and source records attributed to Danielewicz, P..

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of fragment formation on shear viscosity in the nuclear liquid-gas phase transition region

Within the improved quantum molecular dynamics (ImQMD) model we follow the evolution of nuclear matter for planar Couette flow in a periodic box. We focus on the region of liquid-gas phase transition and extract the shear viscosity coefficient from the local stress tensor, directly following the viscosity definition. Here, by switching on and off the mean field and thus inducing the phase transition, we are able to observe the impact of clumping in the phase transition region on the viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microscopic optical potentials for calcium isotopes

We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for 40 Ca and 48 Ca. For the computation of the ground states of 40 Ca and 48 Ca, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = ±1 nuclei we use particle-attached and particle-removed equation-of-motion method truncated at two-particle–one-hole and one-particle–two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction NNLO sat , which reproduces the charge radii of 40 Ca and 48 Ca, and the chiral nucleon-nucleon interaction NNLO opt . In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. Furthermore, the latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter η in the Green's function, practical results can be further improved.

74 ATOMIC AND MOLECULAR PHYSICS↗

Pion production in rare-isotope collisions

Pion energy spectra are presented for central collisions of neutron-rich 132 Sn + 124 Sn and neutron-deficient 108 Sn + 112 Sn systems using simulations with the Boltzmann-Uehling-Uhlenbeck transport model. These calculations, which incorporate isospin-dependent mean-field potentials for relevant baryons and mesons, display a sensitivity to the pion spectra that could allow significant constraints on the density dependence of the symmetry energy and its mean-field potential at suprasaturation densities. As a result, the predicted sensitivity increases with the isospin asymmetry of the total system and decreases with incident energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optical potential from first principles

Here, we develop a method to construct a microscopic optical potential from chiral interactions for nucleon-nucleus scattering. The optical potential is constructed by combining the Green’s function approach with the coupled-cluster method. To deal with the poles of the Green’s function along the real energy axis we employ a Berggren basis in the complex energy plane combined with the Lanczos method. Using this approach, we perform a proof-of-principle calculation of the optical potential for the elastic neutron scattering on 16 O. For the computation of the ground-state of 16 O, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A ±1 nuclei we use particle-attached/removed equation-of-motion method truncated at two-particle-one-hole and one-particle-two-hole excitations, respectively. We verify the convergence of the optical potential and scattering phase shifts with respect to the model-space size and the number of discretized complex continuum states. We also investigate the absorptive component of the optical potential (which reflects the opening of inelastic channels) by computing its imaginary volume integral and find an almost negligible absorptive component at low-energies. To shed light on this result, we computed excited states of 16 O using equation-of-motion coupled-cluster method with singles-and- doubles excitations and we found no low-lying excited states below 10 MeV. Furthermore, most excited states have a dominant two-particle-two-hole component, making higher-order particle-hole excitations necessary to achieve a precise description of these core-excited states. We conclude that the reduced absorption at low-energies can be attributed to the lack of correlations coming from the low-order cluster truncation in the employed coupled-cluster method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗