Source author record

Pawel Danielewicz

Pawel Danielewicz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions

Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

preprint2021arXiv

Deblurring for Nuclei: 3D Characteristics of Heavy-Ion Collisions

Observables from nuclear and high-energy experiments can be degraded by detector performance and/or methodology in extracting the observables, such as of the final-state characteristics of heavy-ion collisions in relation to a coarsely estimated reaction-plane direction. We propose the use of deblurring methods, such as in optics, to correct for observable degradation. Our main focus is the restoration of triple-differential particle distributions in heavy-ion collisions. We demonstrate that these could be extracted from collision measurements following the Richardson-Lucy deblurring method from optics. We illustrate basic features of the restoration methodology in a schematic model assuming either ideal or more realistic particle detection. The inferred three-dimensional (3D) distributions for collisions may easier to interpret in terms of collision dynamics and sought properties of bulk matter than the currently employed Fourier coefficients, that combine information from different azimuthal angles relative to the reaction plane.

preprint2019arXiv

Dynamics of one-dimensional correlated nuclear systems within non-equilibrium Green's function theory

Theory of non-equilibrium Green's function (NGF) provides a practical framework for studying quantum many-body systems out of equilibrium. Extending the previous mean field approach developed for nuclear systems in one dimension with NGF, we introduce isospin degrees of freedom to the Green's functions and incorporate short-range two-body interactions in the second-order self-consistent approximation to correlations, which represents the scattering of momentum orbitals in the Born approximation. We discuss the preparation of a finite nuclear system and examine the impact of correlations on the ground state. We also excite a finite symmetric nuclear system to oscillate in an isovector dipole mode and explore the dissipation effects in the oscillation. Finally, we demonstrate how to boost a slab to a constant and stable motion in a box, based on Galilean covariance of the theory. The studies in this paper lay the ground for the future exploration of collisions of correlated nuclear systems in one dimension.

preprint2016arXiv

Symmetry Energy III: Isovector Skins

Isoscalar density is a sum of neutron and proton densities and isovector is a normalized difference. Here, we report on the experimental evidence for the displacement of the isovector and isoscalar surfaces in nuclei, by $\sim$$0.9 \, \text{fm}$ from each other. We analyze data on quasielastic (QE) charge exchange (p,n) reactions, concurrently with proton and neutron elastic scattering data for the same target nuclei, following the concepts of the isoscalar and isovector potentials combined into Lane optical potential. The elastic data largely probe the geometry of the isoscalar potential and the (p,n) data largely probe a relation between the geometries of the isovector and isoscalar potentials. The targets include $^{48}$Ca, $^{90}$Zr, $^{120}$Sn and $^{208}$Pb and projectile incident energy values span the range of (10-50)$\,\text{MeV}$. In our fit to elastic and QE charge-exchange data, we allow the values of isoscalar and isovector radii, diffusivities and overall potential normalizations to float away from those in the popular Koning and Delaroche parametrization. We find that the best-fit isovector radii are consistently larger than isoscalar and the best-fit isovector surfaces are steeper. Upon identifying the displacement of the potential surfaces with the displacement of the surfaces for the densities in the Skyrme-Hartree-Fock calculations, and by supplementing the results with those from analysing excitation energies to isobaric analog states in the past, we arrive at the slope and value of the symmetry energy at normal density of $70<L<101\, \text{MeV}$ and $33.5<a_a^V<36.4\, \text{MeV}$, respectively.

preprint2016arXiv

Towards Quantum Transport for Central Nuclear Reactions

Nonequilibrium Green's functions represent a promising tool for describing central nuclear reactions. Even at the single-particle level, though, the Green's functions contain more information that computers may handle in the foreseeable future. In this study, we explore slab collisions in one dimension, first in the mean field approximation and demonstrate that only function elements close to the diagonal in arguments are relevant, in practice, for the reaction calculations. This bodes well for the application of the Green's functions to the reactions. Moreover we demonstrate that an initial state for a reaction calculation may be generated through adiabatic transformation of interactions. Finally, we report on our progress in incorporating correlations into the dynamic calculations.

preprint2016arXiv

Understanding transport simulations of heavy-ion collisions at 100 and 400 AMeV: Comparison of heavy ion transport codes under controlled conditions

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-Molecular-Dynamics-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.

preprint2016arXiv

White Paper on Nuclear Astrophysics

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

preprint2013arXiv

Symmetry Energy II: Isobaric Analog States

Using excitation energies to isobaric analog states (IAS) and charge invariance, we extract nuclear symmetry coefficients, from a mass formula, on a nucleus-by-nucleus basis. Consistently with charge invariance, the coefficients vary weakly across an isobaric chain. However, they change strongly with nuclear mass and range from a_a~10 MeV at mass A~10 to a_a~22 MeV at A~240. Following the considerations of a Hohenberg-Kohn functional for nuclear systems, we determine how to find in practice the symmetry coefficient using neutron and proton densities, even when those densities are simultaneously affected by significant symmetry-energy and Coulomb effects. These results facilitate extracting the symmetry coefficients from Skyrme-Hartree-Fock (SHF) calculations, that we carry out using a variety of Skyrme parametrizations in the literature. For the parametrizations, we catalog novel short-wavelength instabilities. In comparing the SHF and IAS results for the symmetry coefficients, we arrive at narrow (+-2.4 MeV) constraints on the symmetry energy values S(rho) at 0.04<rho<0.13 fm^-3. Towards normal density the constraints significantly widen, but the normal value of energy a_a^V and the slope parameter L are found to be strongly correlated. To narrow the constraints, we reach for the measurements of asymmetry skins and arrive at a_a^V=(30.2-33.7) MeV and L=(35-70) MeV, with those values being again strongly positively correlated along the diagonal of their combined region. Inclusion of the skin constraints allows to narrow the constraints on S(rho), at 0.04<rho<0.13 fm^-3, down to +-1.1 MeV. Several microscopic calculations, including variational, Bruckner-Hartree-Fock and Dirac-Bruckner-Hartree-Fock, are consistent with our constraint region on S(rho).

preprint2011arXiv

Symmetry Energy from Systematic of Isobaric Analog States

Excitation energies to isobaric states, that are analogs of ground states, are dominated by contributions from the symmetry energy. This opens up a possibility of investigating the symmetry energy on nucleus-by-nucleus basis. Upon correcting energies of measured nuclear levels for shell and pairing effects, we find that the lowest energies for a given isospin rise in proportion to the square of isospin, allowing for an interpretation of the coefficient of proportionality in terms of a symmetry coefficient for a given nucleus. In the (A,Z) regions where there are enough data, we demonstrate a Z-independence of that coefficient. We further concentrate on the A-dependence of the coefficient, in order to learn about the density dependence of symmetry energy in uniform matter, given the changes of the density in the surface region. In parallel to the analysis of data, we carry out an analysis of the coefficient for nuclei calculated within the Skyrme-Hartree-Fock (SHF) approach, with known symmetry energy for uniform matter. While the data from isobaric analog states suggest a simple interpretation for the A-dependent symmetry coefficient, in terms of the surface and volume symmetry coefficients, the SHF results point to a more complicated situation within the isovector sector than in the isoscalar, with much stronger curvature effects in the first. We exploit the SHF results in estimating the curvature contributions to the symmetry coefficient. That assessment is hampered by instabilities of common Skyrme parameterizations of nuclear interactions.

preprint2010arXiv

Towards a nonequilibrium Green's function description of nuclear reactions: one-dimensional mean-field dynamics

Nonequilibrium Green's function methods allow for an intrinsically consistent description of the evolution of quantal many-body body systems, with inclusion of different types of correlations. In this paper, we focus on the practical developments needed to build a Green's function methodology for nuclear reactions. We start out by considering symmetric collisions of slabs in one dimension within the mean-field approximation. We concentrate on two issues of importance for actual reaction simulations. First, the preparation of the initial state within the same methodology as for the reaction dynamics is demonstrated by an adiabatic switching on of the mean-field interaction, which leads to the mean-field ground state. Second, the importance of the Green's function matrix-elements far away from the spatial diagonal is analyzed by a suitable suppression process that does not significantly affect the evolution of the elements close to the diagonal. The relative lack of importance of the far-away elements is tied to system expansion. We also examine the evolution of the Wigner function and verify quantitatively that erasing of the off-diagonal elements corresponds to averaging out of the momentum-space details in the Wigner function.

preprint2009arXiv

Towards Quantum Transport for Nuclear Reactions

Nonequilibrium Green's functions represent a promising tool for describing central nuclear reactions. Even at the single-particle level, though, the Green's functions contain more information that computers may handle in the foreseeable future. In this study, we investigate whether all the information contained in the Green's functions is necessarily relevant when describing the time evolution of nuclear reactions. For this, we carry out mean-field calculations of slab collisions in one dimension.