Showing 1–50 of 133 results for author: Danielewicz, P

Analysis of Iterative Deblurring: No Explicit Noise

Authors: Sinethemba Neliswa Mamba, Pawel Danielewicz

Abstract: Iterative deblurring, notably the Richardson-Lucy algorithm with and without regularization, is analyzed in the context of nuclear and high-energy physics applications. In these applications, probability distributions may be discretized into a few bins, measurement statistics can be high, and instrument performance can be well understood. In such circumstances, it is essential to understand the de… ▽ More Iterative deblurring, notably the Richardson-Lucy algorithm with and without regularization, is analyzed in the context of nuclear and high-energy physics applications. In these applications, probability distributions may be discretized into a few bins, measurement statistics can be high, and instrument performance can be well understood. In such circumstances, it is essential to understand the deblurring first without any explicit noise considerations. We employ singular value decomposition for the blurring matrix in a low-count pixel system. A strong blurring may yield a null space for the blurring matrix. Yet, a nonnegativity constraint for images built into the deblurring may help restore null-space content in a high-contrast image with zero or low intensity for a sufficient number of pixels. For low-contrast images, the control over null-space content may be gained through regularization. When the regularization is applied, the blurred image is, in practice, restored to an image that is still blurred but less than the starting one. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 17 pages, 9 figures

MSC Class: Primary: 65R30; 65R32; Secondary: 65Z05

Side Ridge in Ar + KCl Collisions at 1.8 GeV/nucleon with Reaction-Plane Deblurring

Authors: Pawel Danielewicz, Herbert Stroebele, Pierre Nzabahimana

Abstract: Reaction-plane deblurring is applied to the triple-differential distributions of protons, deuterons and negatively charged pions from central Ar + KCl collisions at 1.8 GeV/nucleon, measured with the LBL-GSI Streamer Chamber, to yield distributions relative to the true reaction plane of the collisions. Within the reaction plane and in the forward cm rapidity region a side ridge away from the beam… ▽ More Reaction-plane deblurring is applied to the triple-differential distributions of protons, deuterons and negatively charged pions from central Ar + KCl collisions at 1.8 GeV/nucleon, measured with the LBL-GSI Streamer Chamber, to yield distributions relative to the true reaction plane of the collisions. Within the reaction plane and in the forward cm rapidity region a side ridge away from the beam axis is observed: the distributions peak away from the beam direction for protons and deuterons but not for pions. These findings are consistent with a source of particles moving at an in-plane transverse velocity of ~0.1 c. Transport pBUU calculations yield results in semi-quantitative agreement with those from the data. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 5 pages, 3 figures

Comparing pion production in transport simulations of heavy-ion collisions at $270A$ MeV under controlled conditions

Authors: Jun Xu, Hermann Wolter, Maria Colonna, Mircea Dan Cozma, Pawel Danielewicz, Che Ming Ko, Akira Ono, ManYee Betty Tsang, Ying-Xun Zhang, Hui-Gan Cheng, Natsumi Ikeno, Rohit Kumar, Jun Su, Hua Zheng, Zhen Zhang, Lie-Wen Chen, Zhao-Qing Feng, Christoph Hartnack, Arnaud Le Fèvre, Bao-An Li, Yasushi Nara, Akira Ohnishi, Feng-Shou Zhang

Abstract: Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at $270A$ MeV, and put particular emphasis on the production of pions and $Δ$ resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleon… ▽ More Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at $270A$ MeV, and put particular emphasis on the production of pions and $Δ$ resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleonic evolution of the collision and the nucleonic observables in these codes do not completely converge, but the differences among the codes can be understood as being due to several reasons: the basic differences between BUU and QMD models in the representation of the phase-space distributions, computational differences in the mean-field evaluation, and differences in the adopted strategies for the Pauli blocking in the collision integrals. For pionic observables, we find that a higher maximum density leads to an enhanced pion yield and a reduced $π^-/π^+$ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced $π^-/π^+$ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total $π^-/π^+$ yield ratio but reduces the ratio at high pion energies, although differences in its implementations do not have a dominating role in the differences among the codes. Taking into account only the results of codes that strictly follow the homework specifications, we find a convergence of the codes in the final charged pion yield ratio to a $1σ$ deviation of about $5\%$. However, the uncertainty is expected to be reduced to about $1.6\%$ if the same or similar strategies and ingredients, i.e., an improved Pauli blocking and calculation of the non-linear term in the mean-field potential, are similarly used in all codes. △ Less

Submitted 14 March, 2024; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 39 pages, 17 figures

Journal ref: Physical Review C 109, 044609 (2024)

Source Function from Two-Particle Correlation Through Deblurring

Authors: Pierre Nzabahimana, Pawel Danielewicz

Abstract: In heavy-ion collisions, low relative-velocity two-particle correlations have been a tool for assessing space-time characteristics of particle emission. Those characteristics may be cast in the form of a relative emission source related to the correlation function through the Koonin-Pratt (KP) convolution formula that involves the relative wave-function for the particles in its kernel. In the lite… ▽ More In heavy-ion collisions, low relative-velocity two-particle correlations have been a tool for assessing space-time characteristics of particle emission. Those characteristics may be cast in the form of a relative emission source related to the correlation function through the Koonin-Pratt (KP) convolution formula that involves the relative wave-function for the particles in its kernel. In the literature, the source has been most commonly sought by parametrizing it in a Gaussian form and fitting to the correlation function. At times the source was more broadly imaged from the function, still employing a fitting. Here, we propose the use of the Richardson-Lucy (RL) optical deblurring algorithm for deducing the source from a correlation function. The RL algorithm originally follows from probabilistic Bayesian considerations and relies on the intensity distributions for the optical object and its image, as well as the convolution kernel, being positive definite, which is the case for the corresponding quantities of interest within the KP formula. △ Less

Submitted 14 August, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of mu… ▽ More The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Long Range Plan: Dense matter theory for heavy-ion collisions and neutron stars

Authors: Alessandro Lovato, Travis Dore, Robert D. Pisarski, Bjoern Schenke, Katerina Chatziioannou, Jocelyn S. Read, Philippe Landry, Pawel Danielewicz, Dean Lee, Scott Pratt, Fabian Rennecke, Hannah Elfner, Veronica Dexheimer, Rajesh Kumar, Michael Strickland, Johannes Jahan, Claudia Ratti, Volodymyr Vovchenko, Mikhail Stephanov, Dekrayat Almaalol, Gordon Baym, Mauricio Hippert, Jacquelyn Noronha-Hostler, Jorge Noronha, Enrico Speranza , et al. (39 additional authors not shown)

Abstract: Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theo… ▽ More Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science △ Less

Submitted 7 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

Report number: LA-UR-22-31648

Exploring sensitivity of charge-exchange ($p, n$) reactions to the neutron density distribution

Authors: Jian Liu, Yunsheng Wang, Yonghao Gao, Pawel Danielewicz, Chang Xu, Zhongzhou Ren

Abstract: $Background:$ The determination of the nuclear neutron properties suffers from uncontrolled uncertainties, which attracted considerable attention recently, such as in the context of the PREX experiment. $Purpose:$ Our aim is to analyze the sensitivity of charge-exchange ($p, n$) reactions to the neutron density distribution $ρ_{n}… ▽ More $Background:$ The determination of the nuclear neutron properties suffers from uncontrolled uncertainties, which attracted considerable attention recently, such as in the context of the PREX experiment. $Purpose:$ Our aim is to analyze the sensitivity of charge-exchange ($p, n$) reactions to the neutron density distribution $ρ_{n}$ and constrain the neutron characteristics in the nuclear structure models. $Method:$ By combing the folding and the mean-field models, the nucleon-nucleus ($NA$) potential can be obtained from the nuclear density distribution. Further, the ($p, p$) and ($p, n$) cross sections for $^{48}$Ca and $^{208}$Pb are calculated following the distorted-wave Born approximation (DWBA) method. $Results:$ Compared with the ($p, p$) cross section, the effects of $ρ_{n}$ variation on the ($p, n$) cross section are significant, which is due to the impact of isovector properties. Based on the global folding model analyses of data, it is found that $^{48}$Ca and $^{208}$Pb have relatively large neutron skin thickness $ΔR_{n p}$. $Conclusions:$ Results illustrate that the charge-exchange ($p, n$) reaction is a sensitive probe of $ρ_{n}$. The results in this paper can offer useful guides for future experiments of neutron characteristics. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Optical potentials for the rare-isotope beam era

Authors: C. Hebborn, F. M. Nunes, G. Potel, W. H. Dickhoff, J. W. Holt, M. C. Atkinson, R. B. Baker, C. Barbieri, G. Blanchon, M. Burrows, R. Capote, P. Danielewicz, M. Dupuis, Ch. Elster, J. E. Escher, L. Hlophe, A. Idini, H. Jayatissa, B. P. Kay, K. Kravvaris, J. J. Manfredi, A. Mercenne, B. Morillon, G. Perdikakis, C. D. Pruitt , et al. (4 additional authors not shown)

Abstract: We review recent progress and motivate the need for further developments in nuclear optical potentials that are widely used in the theoretical analysis of nucleon elastic scattering and reaction cross sections. In regions of the nuclear chart away from stability, which represent a frontier in nuclear science over the coming decade and which will be probed at new rare-isotope beam facilities worldw… ▽ More We review recent progress and motivate the need for further developments in nuclear optical potentials that are widely used in the theoretical analysis of nucleon elastic scattering and reaction cross sections. In regions of the nuclear chart away from stability, which represent a frontier in nuclear science over the coming decade and which will be probed at new rare-isotope beam facilities worldwide, there is a targeted need to quantify and reduce theoretical reaction model uncertainties, especially with respect to nuclear optical potentials. We first describe the primary physics motivations for an improved description of nuclear reactions involving short-lived isotopes, focusing on its benefits for fundamental science discoveries and applications to medicine, energy, and security. We then outline the various methods in use today to build optical potentials starting from phenomenological, microscopic, and ab initio methods, highlighting in particular the strengths and weaknesses of each approach. We then discuss publicly-available tools and resources facilitating the propagation of recent progresses in the field to practitioners. Finally, we provide a set of open challenges and recommendations for the field to advance the fundamental science goals of nuclear reaction studies in the rare-isotope beam era. △ Less

Submitted 10 March, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: This paper is the outcome of the Facility for Rare Isotope Beams Theory Alliance (FRIB - TA) topical program "Optical Potentials in Nuclear Physics" held in March 2022 at FRIB. Its content is non-exhaustive, was chosen by the participants and reflects their efforts related to optical potentials

Report number: LLNL-JRNL-840487

Deconstructing experimental decay energy spectra: the $^{26}$O case

Authors: Pierre Nzabahimana, Thomas Redpath, Thomas Baumann, Pawel Danielewicz, Pablo Giuliani, Paul Guèye

Abstract: In nuclear reaction experiments, the measured decay energy spectra can give insights into the shell structure of decaying systems. However, extracting the underlying physics from the measurements is challenging due to detector resolution and acceptance effects. The Richardson-Lucy (RL) algorithm, a deblurring method that is commonly used in optics and has proven to be a successful technique for re… ▽ More In nuclear reaction experiments, the measured decay energy spectra can give insights into the shell structure of decaying systems. However, extracting the underlying physics from the measurements is challenging due to detector resolution and acceptance effects. The Richardson-Lucy (RL) algorithm, a deblurring method that is commonly used in optics and has proven to be a successful technique for restoring images, was applied to our experimental nuclear physics data. The only inputs to the method are the observed energy spectrum and the detector's response matrix also known as the transfer matrix. We demonstrate that the technique can help access information about the shell structure of particle-unbound systems from the measured decay energy spectrum that isn't immediately accessible via traditional approaches such as chi-square fitting. For a similar purpose, we developed a machine learning model that uses a deep neural network (DNN) classifier to identify resonance states from the measured decay energy spectrum. We tested the performance of both methods on simulated data and experimental measurements. Then, we applied both algorithms to the decay energy spectrum of $^{26}\mathrm{O} \rightarrow ^{24}\mathrm{O}$ + n + n measured via invariant mass spectroscopy. The resonance states restored using the RL algorithm to deblur the measured decay energy spectrum agree with those found by the DNN classifier. Both deblurring and DNN approaches suggest that the raw decay energy spectrum of $^{26}\mathrm{O}$ exhibits three peaks at approximately 0.15~MeV, 1.50~MeV, and 5.00~MeV, with half-widths of 0.29~MeV, 0.80~MeV, and 1.85~MeV, respectively. △ Less

Submitted 15 February, 2023; v1 submitted 30 September, 2022; originally announced October 2022.

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

Authors: X. G. Den, P. Danielewicz, Y. G. Ma, H. Lin, Y. X. Zhang

Abstract: Within the improved quantum molecular dynamic (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 viscosity definition. By switching on and off the mean field and thus inducing the phase transition, we are… ▽ More Within the improved quantum molecular dynamic (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 viscosity definition. 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 onto the viscosity. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 10 pages, 9 figures, 2 tables; Physical Review C, in press

Journal ref: Physical Review C 105, 064613 (2022); Editors' Suggestion

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

Authors: Hermann Wolter, Maria Colonna, Dan Cozma, Pawel Danielewicz, Che Ming Ko, Rohit Kumar, Akira Ono, ManYee Betty Tsang, Jun Xu, Ying-Xun Zhang, Elena Bratkovskaya, Zhao-Qing Feng, Theodoros Gaitanos, Arnaud Le Fèvre, Natsumi Ikeno, Youngman Kim, Swagata Mallik, Paolo Napolitani, Dmytro Oliinychenko, Tatsuhiko Ogawa, Massimo Papa, Jun Su, Rui Wang, Yong-Jia Wang, Janus Weil , et al. (27 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 4 May, 2022; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 114 pages, 14 figures, 479 references, accepted for publication in Progress of Particle and Nuclear Phsics

Journal ref: Prog. Part. Nucl. Phys. 125 (2022) 103962

Transverse momentum analysis of collective motion in relativistic nuclear collisions

Authors: P. Danielewicz, G. Odyniec

Abstract: Novel transverse-momentum technique is used to analyse charged-particle exclusive data for collective motion in the Ar+KCl reaction at 1.8 GeV/nucl. Previous analysis of this reaction, employing the standard sphericity tensor, revealed no significant effect. In the present analysis, collective effects are observed, and they are substantially stronger than in the Cugnon cascade model, but weaker th… ▽ More Novel transverse-momentum technique is used to analyse charged-particle exclusive data for collective motion in the Ar+KCl reaction at 1.8 GeV/nucl. Previous analysis of this reaction, employing the standard sphericity tensor, revealed no significant effect. In the present analysis, collective effects are observed, and they are substantially stronger than in the Cugnon cascade model, but weaker than in the hydrodynamical model. △ Less

Submitted 11 September, 2021; originally announced September 2021.

Comments: Preprint version, serving as a reference for another work (14 pages, 3 figures). This version contains, in particular, calculation of second-order azimuthal moments of particle distributions, in the elements of the sphericity tensor, not in the journal. The editor insisted on shortening of the paper, while the referee insisted on additional explanations, so something had to go

Report number: LBL-18600

Journal ref: Phys. Lett. B 157, 146 (1985)

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

Authors: Pawel Danielewicz, Mizuki Kurata-Nishimura

Abstract: 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 restor… ▽ More 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. △ Less

Submitted 29 December, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 17 pages,11 figures, 26 references

Journal ref: Phys. Rev. C 105, 034608 (2022)

Comparison of Heavy-Ion Transport Simulations: Mean-field Dynamics in a Box

Authors: Maria Colonna, Ying-Xun Zhang, Yong-Jia Wang, Dan Cozma, Pawel Danielewicz, Che Ming Ko, Akira Ono, Manyee Betty Tsang, Rui Wang, Hermann Wolter, Jun Xu, Zhen Zhang, Lie-Wen Chen, Hui-Gan Cheng, Hannah Elfner, Zhao-Qing Feng, Myungkuk Kim, Youngman Kim, Sangyong Jeon, Chang-Hwan Lee, Bao-An Li, Qing-Feng Li, Zhu-Xia Li, Swagata Mallik, Dmytro Oliinychenko , et al. (4 additional authors not shown)

Abstract: Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared… ▽ More Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. The significance of these results for the description of real heavy-ion collisions is discussed. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Journal ref: Phys. Rev. C 104, 024603 (2021)

Symmetry energy investigation with pion production from Sn+Sn systems

Authors: G. Jhang, J. Estee, J. Barney, G. Cerizza, M. Kaneko, J. W. Lee, W. G. Lynch, T. Isobe, M. Kurata-Nishimura, T. Murakami, C. Y . Tsang, M. B. Tsang, R. Wang, D. S. Ahn, L. Atar, T. Aumann, H. Baba, K. Boretzky, J. Brzychczyk, N. Chiga, N. Fukuda, I. Gasparic, B. Hong, A. Horvat, K. Ieki , et al. (55 additional authors not shown)

Abstract: In the past two decades, pions created in the high density regions of heavy ion collisions have been predicted to be sensitive at high densities to the symmetry energy term in the nuclear equation of state, a property that is key to our understanding of neutron stars. In a new experiment designed to study the symmetry energy, the multiplicities of negatively and positively charged pions have been… ▽ More In the past two decades, pions created in the high density regions of heavy ion collisions have been predicted to be sensitive at high densities to the symmetry energy term in the nuclear equation of state, a property that is key to our understanding of neutron stars. In a new experiment designed to study the symmetry energy, the multiplicities of negatively and positively charged pions have been measured with high accuracy for central $^{132}$Sn+$^{124}$Sn, $^{112}$Sn+$^{124}$Sn, and $^{108}$Sn+$^{112}$Sn collisions at $E/A=270~\mathrm{MeV}$ with the S$π$RIT Time Projection Chamber. While the uncertainties of individual pion multiplicities are measured to 4\%, those of the charged pion multiplicity ratios are measured to 2\%. We compare these data to predictions from seven major transport models. The calculations reproduce qualitatively the dependence of the multiplicities and their ratios on the total neutron to proton number in the colliding systems. However, the predictions of the transport models from different codes differ too much to allow extraction of reliable constraints on the symmetry energy from the data. This finding may explain previous contradictory conclusions on symmetry energy constraints obtained from pion data in Au+Au system. These new results call for better understanding of the differences among transport codes, and new observables that are more sensitive to the density dependence of the symmetry energy. △ Less

Submitted 13 December, 2020; originally announced December 2020.

Comments: 8 pages, 4 figures, 1 table (accepted for publication in PLB)

Journal ref: Physics Letters B 813 (2021) 136016

Impact of the neutron-star deformability on equation of state parameters

Authors: C. Y. Tsang, M. B. Tsang, Pawel Danielewicz, W. G. Lynch, F. J. Fattoyev

Abstract: We use a Bayesian inference analysis to explore the sensitivity of Taylor expansion parameters of the nuclear equation of state (EOS) to the neutron star dimensionless tidal deformability ($Λ$) on 1 to 2 solar masses neutron stars. A global power law dependence between tidal deformability and compactness parameter (M/R) is verified over this mass region. To avoid superfluous correlations between t… ▽ More We use a Bayesian inference analysis to explore the sensitivity of Taylor expansion parameters of the nuclear equation of state (EOS) to the neutron star dimensionless tidal deformability ($Λ$) on 1 to 2 solar masses neutron stars. A global power law dependence between tidal deformability and compactness parameter (M/R) is verified over this mass region. To avoid superfluous correlations between the expansion parameters, we use a correlation-free EOS model based on a recently published meta-modeling approach. We find that assumptions in the prior distribution strongly influence the constraints on $Λ$. The $Λ$ constraints obtained from the neutron star merger event GW170817 prefer low values of $L_\text{sym}$ and $K_\text{sym}$, for a canonical neutron star with 1.4 solar mass. For neutron star with mass $<1.6$ solar mass, $L_\text{sym}$ and $K_\text{sym}$ are highly correlated with the tidal deformability. For more massive neutron stars, the tidal deformability is more strongly correlated with higher order Taylor expansion parameters. △ Less

Submitted 11 September, 2020; originally announced September 2020.

Journal ref: Phys. Rev. C 102, 045808 (2020)

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

Authors: Hao Lin, Hossein Mahzoon, Arnau Rios, Pawel Danielewicz

Abstract: 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 approximat… ▽ More 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. △ Less

Submitted 24 December, 2019; v1 submitted 23 December, 2019; originally announced December 2019.

Spinodal Instability at the Onset of Collective Expansion in Nuclear Collisions

Authors: Pawel Danielewicz, Hao Lin, Jirina R. Stone, Yoritaka Iwata

Abstract: Using transport theory to model central Au + Au collisions in the energy region of 20 - 110 MeV/u, at impact parameters b <= 5 fm, we predict a measurable impact of spinoidal instability as the collective expansion sets in with energy. Two transport models are employed, the pBUU model, solving a Boltzmann-Uehling-Uhlenbeck equation, and the Brownian Motion (BM) model, solving a set of Langevin equ… ▽ More Using transport theory to model central Au + Au collisions in the energy region of 20 - 110 MeV/u, at impact parameters b <= 5 fm, we predict a measurable impact of spinoidal instability as the collective expansion sets in with energy. Two transport models are employed, the pBUU model, solving a Boltzmann-Uehling-Uhlenbeck equation, and the Brownian Motion (BM) model, solving a set of Langevin equations to describe the motion of individual nucleons in a noisy nuclear medium. We find without ambiguity, for the first time, that a combination of delayed equilibration, onset of collective expansion and the spinodal instability produces a pair of transient ring structures, made of the projectile and target remnants, with spectator nucleons predicted to end in the entities reminiscent of stones in jewelry, on the rings. The ring structures, calculated in the configuration space and mapped onto the velocity space, could be detected in experimental collective flow data. △ Less

Submitted 17 December, 2020; v1 submitted 23 October, 2019; originally announced October 2019.

Comments: 9 pages, 5 figures

Journal ref: Phys. Lett. B 811, 135957 (2020)

Symmetry energy constraints from GW170817 and laboratory experiments

Authors: M. B. Tsang, W. G. Lynch, P. Danielewicz, C. Y. Tsang

Abstract: The LIGO-Virgo collaboration detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not elucidate its origins, by not distinguishing the contribution to the pressure from symmetry energy which governs much of the internal structure of a… ▽ More The LIGO-Virgo collaboration detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not elucidate its origins, by not distinguishing the contribution to the pressure from symmetry energy which governs much of the internal structure of a neutron star. By combining the neutron star EoS extracted from the GW170817 event and the EoS of symmetric matter from nucleus-nucleus collision experiments, we extract the symmetry pressure, which is the difference in pressure between neutron and nuclear matter over the density region from 1.2$ρ_{0}$ to $4.5ρ_{0}$. While the uncertainties in the symmetry pressure are large, they can be reduced with new experimental and astrophysical results. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1901.07673

Insights on Skyrme parameters from GW170817

Authors: C. Y. Tsang, M. B. Tsang, Pawel Danielewicz, W. G. Lynch, F. J. Fattoyev

Abstract: The binary neutron-star merger event, GW170817, has cast a new light on nuclear physics research. Using a neutron-star model that includes a crust equation of state (EoS), we calculate the properties of a 1.4 solar-mass neutron star. The model incorporates more than 200 Skyrme energy density functionals, which describe nuclear matter properties, in the outer liquid core region of the neutron star.… ▽ More The binary neutron-star merger event, GW170817, has cast a new light on nuclear physics research. Using a neutron-star model that includes a crust equation of state (EoS), we calculate the properties of a 1.4 solar-mass neutron star. The model incorporates more than 200 Skyrme energy density functionals, which describe nuclear matter properties, in the outer liquid core region of the neutron star. We find a power-law relation between the neutron-star tidal deformability, $Λ$, and the neutron-star radius, R. Without an explicit crust EoS, the model predicts smaller R and the difference becomes significant for stars with large radii. To connect the neutron star properties with nuclear matter properties, we confront the predicted values for $Λ$, against the Taylor expansion coefficients of the Skyrme interactions. There is no pronounced correlation between Skyrme parameters in symmetric nuclear matter and neutron star properties. However, we find the strongest correlation between $Λ$ and $K_{sym}$, the curvature of the density dependence of the symmetry energy at saturation density. At twice the saturation density, our calculations show a strong correlation between $Λ$ and total pressure providing guidance to laboratory nucleus-nucleus collision experiments. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comparison of heavy-ion transport simulations: Collision integral with pions and $Δ$ resonances in a box

Authors: Akira Ono, Jun Xu, Maria Colonna, Pawel Danielewicz, Che Ming Ko, Manyee Betty Tsang, Yong-Jia Wang, Hermann Wolter, Ying-Xun Zhang, Lie-Wen Chen, Dan Cozma, Hannah Elfner, Zhao-Qing Feng, Natsumi Ikeno, Bao-An Li, Swagata Mallik, Yasushi Nara, Tatsuhiko Ogawa, Akira Ohnishi, Dmytro Oliinychenko, Jun Su, Taesoo Song, Feng-Shou Zhang, Zhen Zhang

Abstract: We compare ten transport codes for a system confined in a box, aiming at improved handling of the production of $Δ$ resonances and pions, which is indispensable for constraining high-density symmetry energy from observables such as the $π^-/π^+$ yield ratio in heavy-ion collisions. The system in a box is initialized with nucleons at saturation density and at 60 MeV temperature. The reactions… ▽ More We compare ten transport codes for a system confined in a box, aiming at improved handling of the production of $Δ$ resonances and pions, which is indispensable for constraining high-density symmetry energy from observables such as the $π^-/π^+$ yield ratio in heavy-ion collisions. The system in a box is initialized with nucleons at saturation density and at 60 MeV temperature. The reactions $NN\leftrightarrow NΔ$ and $Δ\leftrightarrow Nπ$ are implemented, but the Pauli blocking and the mean-field potential are deactivated in the present comparison. Results are compared to those from the two reference cases of a chemically equilibrated ideal gas mixture and of the rate equation. In the results of the numbers of $Δ$ and $π$, deviations from the reference values are observed in many codes, and they depend significantly on the size of the time step. These deviations are tied to different ways in ordering the sequence of collisions and decays, that take place in the same time step. Better agreements are seen in the reaction rates and the number ratios among the isospin species of $Δ$ and $π$. These are, however, affected by the correlations, which are absent in the Boltzmann equation, but are induced by the way particle scatterings are treated in transport calculations. The uncertainty in the transport-code predictions of the $π^-/π^+$ ratio for the system initialized at n/p = 1.5, after letting the existing $Δ$ resonances decay, is found to be within a few percent, which is sufficiently small so that it does not strongly impact constraining the high-density symmetry energy from heavy-ion collisions. Most of the sources of uncertainties have been understood, and individual codes may be further improved. This investigation will be extended in the future to heavy-ion collisions to ensure the problems identified here remain under control. △ Less

Submitted 27 September, 2019; v1 submitted 5 April, 2019; originally announced April 2019.

Comments: 36 pages, 27 figures; a new Fig. 21 and revised results from some codes, achieving improved and consistent understanding

Journal ref: Phys. Rev. C 100, 044617 (2019)

Towards a better understanding of the symmetry energy within neutron stars

Authors: C. Y. Tsang, M. B. Tsang, P. Danielewicz, W. G. Lynch, F. J. Fattoyev

Abstract: The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not, by itself, elucidate its microscopic origins, including the pressure arising from the symmetry energy, that governs much of the inte… ▽ More The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has expanded efforts to understand the Equation of State (EoS) of nuclear matter. These measurements provide new constraints on the overall pressure, but do not, by itself, elucidate its microscopic origins, including the pressure arising from the symmetry energy, that governs much of the internal structure of a neutron star. To correlate microscopic constraints from nuclear measurements to the GW170817 constraints, we calculate neutron star properties with more than 200 Skyrme energy density functionals that describe properties of nuclei. Calculated neutron-star radii (R) and the tidal deformabilities which show a strong correlation with pressure at twice saturation density. By combining the neutron star EoS extracted from the GW170817 event and the EoS of symmetric matter from nucleus-nucleus collision experiments, we extract the density dependence of the symmetry pressure from 1.2 to 4.5 times saturation density. While the uncertainties in the symmetry pressure are large, they can be reduced with new experimental and astrophysical results. △ Less

Submitted 18 January, 2019; originally announced January 2019.

Comments: 10 pages, 3 figures. arXiv admin note: text overlap with arXiv:1807.06571

Laboratory Probes of the Neutron-Matter Equation of State

Authors: M. B. Tsang, C. Y. Tsang, P. Danielewicz, W. G. Lynch, F. J. Fattoyev

Abstract: To relate constraints from nuclear physics to the tidal deformabilities of neutron stars, we construct a neutron star model that accepts input from a large collection of Skyrme density functions to calculate properties of 1.4 solar-mass neutron stars. We find that restricting this set of Skyrme to density functions that describe nuclear masses, isobaric analog states, and low energy nuclear reacti… ▽ More To relate constraints from nuclear physics to the tidal deformabilities of neutron stars, we construct a neutron star model that accepts input from a large collection of Skyrme density functions to calculate properties of 1.4 solar-mass neutron stars. We find that restricting this set of Skyrme to density functions that describe nuclear masses, isobaric analog states, and low energy nuclear reactions does not sufficiently restrict the predicted neutron-star radii and the tidal deformabilities. However, pressure constraints on the EoS around twice saturation density ($2\times2.74\times10^{14}g/cm^3$), obtained from high energy nucleus-nucleus collisions, does constrain predicted tidal deformabilities with uncertainties smaller than those obtained from the analysis of GW170817. We also found that the density-pressure constraint on the EoS obtained from a recent analysis of the neutron-star merger event agree very well with the density pressure constraints obtained from nuclear physics experiments published in 2002. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: arXiv admin note: substantial text overlap with arXiv:1807.06571

Microscopic optical potentials for calcium isotopes

Authors: J. Rotureau, P. Danielewicz, G. Hagen, G. R. Jansen, F. M. Nunes

Abstract: We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$. For the computation of the ground-state of $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = $\pm 1$ nuclei we use particle-attached/removed equation-of-… ▽ More We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$. For the computation of the ground-state of $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = $\pm 1$ nuclei we use particle-attached/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 $\rm{NNLO_{sat}}$, which reproduces the charge radii of $^{40}$Ca and $^{48}$Ca, and the chiral nucleon-nucleon interaction $\rm{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. 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. △ Less

Submitted 29 October, 2018; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: 11 pages, 10 figures, typos corrected, accepted for publication in Phys. Rev. C

Journal ref: Phys. Rev. C 98, 044625 (2018)

Constraining neutron-star equation of state using heavy-ion collisions

Authors: C. Y. Tsang, M. B. Tsang, P. Danielewicz, W. G. Lynch, F. J. Fattoyev

Abstract: The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has intensified efforts towards the understanding of the equation of state (EoS) of nuclear matter. In this letter, we compare directly the density-pressure constraint on the EoS obtained from a recent analysis of the neutron-star merger event to density-pressure constraints obtained from nucl… ▽ More The LIGO-Virgo collaboration ground-breaking detection of the binary neutron-star merger event, GW170817, has intensified efforts towards the understanding of the equation of state (EoS) of nuclear matter. In this letter, we compare directly the density-pressure constraint on the EoS obtained from a recent analysis of the neutron-star merger event to density-pressure constraints obtained from nuclear physics experiments. To relate constraints from nuclear physics to the radii and the tidal deformabilities of neutron stars, we use a large collection of Skyrme density functionals that describe properties of nuclei to calculate properties of 1.4 solar maass neutron stars. We find that restricting this set of Skyrme equations of state to density functionals that describe nuclear masses, isobaric analog states, and low energy nuclear reactions does not sufficiently restrict the predicted neutron-star radii and the tidal deformabilities. Including pressure constraints on the EoS around twice saturation density, obtained from high energy nucleus-nucleus collisions, does constrain predicted radii and tidal deformabilities to be consistent with the results obtained from the analysis of GW170817. We discuss how new measurements of nucleus-nucleus collisions can improve these constraints on the EoS to be more restrictive than the current constraints from the GW170817 merger event. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: 11 pages, 3 figures

One-body Langevin dynamics in heavy-ion collisions at intermediate energies

Authors: Hao Lin, Pawel Danielewicz

Abstract: We present a new framework to treat the dissipation and fluctuation dynamics associated with nucleon-nucleon scattering in heavy-ion collisions. The two-body collision processes are effectively described in terms of the diffusion of nucleons in viscous nuclear media, governed by a set of Langevin equations in momentum space. The new framework combined with the usual mean field dynamics can be used… ▽ More We present a new framework to treat the dissipation and fluctuation dynamics associated with nucleon-nucleon scattering in heavy-ion collisions. The two-body collision processes are effectively described in terms of the diffusion of nucleons in viscous nuclear media, governed by a set of Langevin equations in momentum space. The new framework combined with the usual mean field dynamics can be used to simulate heavy-ion collisions at intermediate energies. As a proof of principle, we simulate Au + Au reactions and obtain results consistent with other existing codes under the same constrained conditions. We also study the formation of fragments in Sn + Sn reactions at 50 MeV/nucleon, and results are discussed and compared with two other models commonly employed for collisions. △ Less

Submitted 16 July, 2018; originally announced July 2018.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. C 99, 024612 (2019)

Quantum Dynamics of Nuclear Slabs: Mean Field and Short-Range Correlations

Authors: Hossein Mahzoon, Pawel Danielewicz, Arnau Rios

Abstract: Computational difficulties aside, nonequilibrium Green's functions appear ideally suited for investigating the dynamics of central nuclear reactions. Many particles actively participate in those reactions. At the two energy extremes for the collisions, the limiting cases of the Green's function approach have been successful: the time-dependent Hartree-Fock theory at low energy and Boltzmann equati… ▽ More Computational difficulties aside, nonequilibrium Green's functions appear ideally suited for investigating the dynamics of central nuclear reactions. Many particles actively participate in those reactions. At the two energy extremes for the collisions, the limiting cases of the Green's function approach have been successful: the time-dependent Hartree-Fock theory at low energy and Boltzmann equation at high. The strategy for computational adaptation of the Green's function to central reactions is discussed. The strategy involves, in particular, incremental progression from one to three dimensions to develop and assess approximations, discarding of far-away function elements, use of effective interactions and preparation of initial states for the reactions through adiabatic switching. At this stage we concentrate on inclusion of correlations in one dimension, where relatively few approximations are needed, and we carry out reference calculations that can benchmark approximations needed for more dimensions. We switch on short-range interactions generating the correlations adiabatically in the Kadanoff-Baym equations to arrive at correlated ground states for uniform matter. As the energy of the correlated matter does not quite match the expectations for nuclear matter we add mean field to arrive at the match in energy. From there on, we move to finite systems. In switching on the correlations we observe emergence of extended tails in momentum distributions and evolution of single particle occupations away from 1 and 0. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: Contribution to Proc. Frontiers of Quantum and Mesoscopic Thermodynamics, 9-15 July 2017, Prague, Czech Republic, ed. Vaclav Spicka; 10 pages, 7 figures

Comparison of heavy-ion transport simulations: Collision integral in a box

Authors: Ying-Xun Zhang, Yong-Jia Wang, Maria Colonna, Pawel Danielewicz, Akira Ono, Betty Tsang, Hermann Wolter, Jun Xu, Lie-Wen Chen, Dan Cozma, Zhao-Qing Feng, Subal Das Gupta, Natsumi Ikeno, Che-Ming Ko, Bao-An Li, Qing-Feng Li, Zhu-Xia Li, Swagata Mallik, Yasushi Nara, Tatsuhiko Ogawa, Akira Ohnishi, Dmytro Oliinychenko, Massimo Papa, Hannah Petersen, Jun Su , et al. (5 additional authors not shown)

Abstract: Simulations by transport codes are indispensable to extract valuable physics information from heavy ion collisions. In order to understand the origins of discrepancies between 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 temperatu… ▽ More Simulations by transport codes are indispensable to extract valuable physics information from heavy ion collisions. In order to understand the origins of discrepancies between 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 after eliminating the correlations within the same pair of colliding particles. 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 towards a Boltzmann distribution. 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. △ Less

Submitted 17 April, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 23 pages, 10 figures. The abstract is changed a little bit due to the limits character number. The complete abstract is in the paper. Accepted by Phys.Rev.C

Journal ref: Phys. Rev. C 97, 034625 (2018)

Correlations within the Non-Equilibrium Green's Function Method

Authors: M. H. Mahzoon, P. Danielewicz, A. Rios

Abstract: Non-equilibrium Green's Function (NGF) method is a powerful tool for studying the evolution of quantum many-body systems. Different types of correlations can be systematically incorporated within the formalism. The time evolution of the single-particle Green's functions is described in terms of the Kadanoff-Baym equations. The current work initially focuses on introducing the correlations within i… ▽ More Non-equilibrium Green's Function (NGF) method is a powerful tool for studying the evolution of quantum many-body systems. Different types of correlations can be systematically incorporated within the formalism. The time evolution of the single-particle Green's functions is described in terms of the Kadanoff-Baym equations. The current work initially focuses on introducing the correlations within infinite nuclear matter in one dimension and then in a finite system in the NGF approach. Starting from the harmonic oscillator Hamiltonian, by switching on adiabatically the mean-field and correlations simultaneously, a correlated state with ground-state characteristics is arrived at within the NGF method. Furthermore the use of cooling to for improving the adiabatic switching is explored. △ Less

Submitted 23 July, 2017; originally announced July 2017.

Comments: Contribution to Proc. 5th Conference on Nuclei and Mesoscopic Physics, E Lansing, 6-10 March 2017; 9 pages, 8 figures

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions

Authors: J. R. Stone, P. Danielewicz, Y. Iwata

Abstract: 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_{\rm beam} \leq 800 \, \text{MeV/nucl}$ using pBUU and TDHF models. The symmetric $^\text{40}$Ca +$^\text{40}$Ca, $^\text{48}$Ca +$^\text{48}$Ca, $^\text{100}$Sn +$^\text{100}$Sn and $^\text{120}$Sn + $^\text{120}$Sn and asymmetric… ▽ More 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_{\rm beam} \leq 800 \, \text{MeV/nucl}$ using pBUU and TDHF models. The symmetric $^\text{40}$Ca +$^\text{40}$Ca, $^\text{48}$Ca +$^\text{48}$Ca, $^\text{100}$Sn +$^\text{100}$Sn and $^\text{120}$Sn + $^\text{120}$Sn and asymmetric $^\text{40}$Ca +$^\text{48}$Ca and $^\text{100}$Sn +$^\text{120}$Sn systems were chosen for the simulations. 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 proton-neutron 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, but cannot be consistently incorporated into semiclassical transport models. △ Less

Submitted 5 June, 2017; originally announced June 2017.

Comments: 36 pages, 8 figures, submitted to Phys. Rev. C

Journal ref: Phys. Rev. C 96, 014612 (2017)

Pion Production in Rare Isotope Collisions

Authors: M. B. Tsang, J. Estee, H. Setiawan, W. G. Lynch, J. Barney, M. B. Chen, G. Cerizza, P. Danielewicz, J. Hong, P. Morfouace, R. Shane, S. Tangwancharoen, K. Zhu, T. Isobe, M. Kurata-Nishimura, J. Lukasik, T. Murakami, the SπRIT collaboration

Abstract: Pion energy spectra are presented for central collisions of neutron-rich 132Sn+124Sn and neutron-deficient 108Sn+112Sn systems using simulations with 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 dens… ▽ More Pion energy spectra are presented for central collisions of neutron-rich 132Sn+124Sn and neutron-deficient 108Sn+112Sn systems using simulations with 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 supra-saturation densities. The predicted sensitivity increases with the isospin asymmetry of the total system and decreases with incident energy. △ Less

Submitted 12 March, 2017; v1 submitted 20 December, 2016; originally announced December 2016.

Journal ref: Phys. Rev. C 95, 044614 (2017)

Shear viscosity from nuclear stopping

Authors: Brent Barker, Pawel Danielewicz

Abstract: Within a Boltzmann transport model, we demonstrate correlation between stopping observables and shear viscosity in central nuclear collisions at intermediate energies (on the order of 10 to 1000 MeV/nucleon). The correlation allows us to assess the viscosity of nuclear matter, by tuning the in-medium nucleon-nucleon cross section in our transport model to agree with nuclear stopping data. We also… ▽ More Within a Boltzmann transport model, we demonstrate correlation between stopping observables and shear viscosity in central nuclear collisions at intermediate energies (on the order of 10 to 1000 MeV/nucleon). The correlation allows us to assess the viscosity of nuclear matter, by tuning the in-medium nucleon-nucleon cross section in our transport model to agree with nuclear stopping data. We also calculate the ratio of shear viscosity to entropy density to determine how close the system is to the universal quantum lower limit proposed in the context of ultrarelativistic heavy ion collisions. △ Less

Submitted 12 March, 2019; v1 submitted 14 December, 2016; originally announced December 2016.

Comments: 31 pages, 12 figures

Journal ref: Phys. Rev. C 99, 034607 (2019)

Optical potential from first principles

Authors: J. Rotureau, P. Danielewicz, G. Hagen, F. Nunes, T. Papenbrock

Abstract: 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. Us… ▽ More 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}{\rm O}$. For the computation of the ground-state of $^{16}{\rm O}$, we use the coupled-cluster method in the singles-and-doubles approximation, while for the $A=\pm 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}{\rm 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. △ Less

Submitted 29 January, 2017; v1 submitted 14 November, 2016; originally announced November 2016.

Comments: 12 pages, 12 figures, a few points clarified, typos corrected, accepted for publication in Phys. Rev. C

Journal ref: Phys. Rev. C 95, 024315 (2017)

Symmetry Energy III: Isovector Skins

Authors: Pawel Danielewicz, Pardeep Singh, Jenny Lee

Abstract: 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 th… ▽ More 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. △ Less

Submitted 23 November, 2016; v1 submitted 6 November, 2016; originally announced November 2016.

Comments: 58 pages, 21 figures, 3 tables; version accepted for Nuclear Physics A

Report number: INT-PUB-16-036

Journal ref: Nucl. Phys. 958 (2017) 147

Light cluster production at NICA

Authors: N. -U. Bastian, P. Batyuk, D. Blaschke, P. Danielewicz, Yu. B. Ivanov, Iu. Karpenko, G. Röpke, O. Rogachevsky, H. H. Wolter

Abstract: Light cluster production at the NICA accelerator complex offers unique possibilities to use these states as "rare probes" of in-medium characteristics such as phase space occupation and early flow. In order to explain this statement, in this contribution theoretical considerations from the nuclear statistical equilibrium model and from a quantum statistical model of cluster production are suppleme… ▽ More Light cluster production at the NICA accelerator complex offers unique possibilities to use these states as "rare probes" of in-medium characteristics such as phase space occupation and early flow. In order to explain this statement, in this contribution theoretical considerations from the nuclear statistical equilibrium model and from a quantum statistical model of cluster production are supplemented with a discussion of a transport model for light cluster formation and with results from hydrodynamic simulations combined with the coalescence model. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Comments: 7 pages, 7 figures, accepted for publication in Eur. Phys. J. A

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

Authors: Jun Xu, Lie-Wen Chen, ManYee Betty Tsang, Hermann Wolter, Ying-Xun Zhang, Joerg Aichelin, Maria Colonna, Dan Cozma, Pawel Danielewicz, Zhao-Qing Feng, Arnaud Le Fevre, Theodoros Gaitanos, Christoph Hartnack, Kyungil Kim, Youngman Kim, Che-Ming Ko, Bao-An Li, Qing-Feng Li, Zhu-Xia Li, Paolo Napolitani, Akira Ono, Massimo Papa, Taesoo Song, Jun Su, Jun-Long Tian , et al. (6 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 26 March, 2016; originally announced March 2016.

Comments: 26 pages, 13 figures, a manuscript on the heavy-ion collision part of transport code comparison project

Journal ref: Phys.Rev.C93:044609,2016

White Paper on Nuclear Astrophysics

Authors: Almudena Arcones, Dan W. Bardayan, Timothy C. Beers, Lee A. Berstein, Jeffrey C. Blackmon, Bronson Messer, B. Alex Brown, Edward F. Brown, Carl R. Brune, Art E. Champagne, Alessandro Chieffi, Aaron J. Couture, Pawel Danielewicz, Roland Diehl, Mounib El-Eid, Jutta Escher, Brian D. Fields, Carla Fröhlich, Falk Herwig, William Raphael Hix, Christian Iliadis, William G. Lynch, Gail C. McLaughlin, Bradley S. Meyer, Anthony Mezzacappa , et al. (18 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2016; v1 submitted 4 March, 2016; originally announced March 2016.

Comments: 131 pages, community white paper based on 2012 JINA Town Meeting in Detroit, MI, and 2014 APS Town Meeting in College Station, TX, corrected author names

Towards Quantum Transport for Central Nuclear Reactions

Authors: Pawel Danielewicz, Arnau Rios, Brent Barker

Abstract: 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… ▽ More 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. △ Less

Submitted 25 February, 2016; originally announced February 2016.

Comments: Contribution to Proc Workshop on Progress in Nonequilibrium Green's Functions VI, Lund U, Sweden, Aug 17-21, 2015, ed. C. Verdozzi; 13 pages, 13 figures. arXiv admin note: text overlap with arXiv:0906.4841

Scaling properties of light-cluster production

Authors: Z. Chajecki, M. Youngs, D. D. S. Coupland, W. G. Lynch, M. B. Tsang, D. Brown, A. Chbihi, P. Danielewicz, R. T. deSouza, M. A. Famiano, T. K. Ghosh, B. Giacherio, V. Henzl, D. Henzlova, C. Herlitzius, S. Hudan, M. A. Kilburn, Jenny Lee, F. Lu, S. Lukyanov, A. M. Rogers, P. Russotto, A. Sanetullaev, R. H. Showalter, L. G. Sobotka , et al. (5 additional authors not shown)

Abstract: We show that ratios of light-particle energy spectra display scaling properties that can be accu- rately described by effective local chemical potentials. This demonstrates the equivalence of t/3He and n/p spectral ratios and provides an essential test of theoretical predictions of isotopically resolved light-particle spectra. In addition, this approach allows direct comparisons of many theoretica… ▽ More We show that ratios of light-particle energy spectra display scaling properties that can be accu- rately described by effective local chemical potentials. This demonstrates the equivalence of t/3He and n/p spectral ratios and provides an essential test of theoretical predictions of isotopically resolved light-particle spectra. In addition, this approach allows direct comparisons of many theoretical n/p spectral ratios to experiments where charged-particle spectra but not neutron spectra are accurately measured. Such experiments may provide much more quantitative constraints on the density and momentum dependence of the symmetry energy. △ Less

Submitted 21 February, 2014; originally announced February 2014.

Subthreshold pion production within a transport description of central Au+Au collisions

Authors: Jun Hong, P. Danielewicz

Abstract: Mapping out the equation of state (EOS) of nuclear matter is a long standing problem in nuclear physics. Recent emphasis is onto the density dependence of the symmetry energy, with experiments needing dedicated symmetry-energy observables. Towards the latter goal, we employ pBUU transport model to simulate pion production in heavy ion collision (HIC). We find that the net pion yield tests the mome… ▽ More Mapping out the equation of state (EOS) of nuclear matter is a long standing problem in nuclear physics. Recent emphasis is onto the density dependence of the symmetry energy, with experiments needing dedicated symmetry-energy observables. Towards the latter goal, we employ pBUU transport model to simulate pion production in heavy ion collision (HIC). We find that the net pion yield tests the momentum dependence of nuclear mean field (MF). In exploring the sensitivity of pion observables to the symmetry energy at higher than normal densities, we find that our calculations of pion ratios contradict, at some level, predictions from both IBUU and ImIQMD models. We propose to employ the pion ratio in the high-energy tail of spectra in future experiments, to distinguish between different variants of high-density symmetry energy. △ Less

Submitted 1 June, 2014; v1 submitted 29 July, 2013; originally announced July 2013.

Comments: 12 pages, 13 figures

Journal ref: Phys. Rev. C 90, 024605 (2014)

Symmetry Energy II: Isobaric Analog States

Authors: Pawel Danielewicz, Jenny Lee

Abstract: 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 consider… ▽ More 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). △ Less

Submitted 15 July, 2013; originally announced July 2013.

Comments: 101 pages, 27 figures, 2 tables; submitted to Nuclear Physics A

Journal ref: Nucl. Phys. A 922, 1 (2014)

Interplay of anisotropies of momentum distribution and mean field in heavy-ion collisions

Authors: C. H. Simon, P. Danielewicz

Abstract: Two important parametrizations of momentum-dependent nucleonic fields, proposed for the simulations of central heavy-ion collisions, one by Gale et al. and the other by Welke et al., suffer from practical limitations. The first gives rise to mean fields isotropic in momentum, even when underlying momentum distributions are anisotropic, making descriptions of early nonequilibrium stages of collisio… ▽ More Two important parametrizations of momentum-dependent nucleonic fields, proposed for the simulations of central heavy-ion collisions, one by Gale et al. and the other by Welke et al., suffer from practical limitations. The first gives rise to mean fields isotropic in momentum, even when underlying momentum distributions are anisotropic, making descriptions of early nonequilibrium stages of collisions unrealistic. The second parametrization gives rise to anisotropic mean fields, but is computationally expensive, because the mean field has to be computed separately for every location of a nucleon in phase space, through folding. Here we construct a parametrization of the nucleonic mean field that yields an anisotropic mean field for an anisotropic momentum distribution and is inexpensive computationally. To demonstrate the versatility of our parametrization, we take the case of results from the parametrization by Welke et al. and attempt to approximate them. In arriving at a suitable anisotropic mean-field potential, we draw, on one hand, from the idea behind the parametrization of Gale et al., of a separable expansion of the potential energy, and, on the other, from the idea of a parallel expansion of the energy and mean field in anisotropy. We show that using our novel parametrization we can qualitatively and partially quantitatively reproduce the features of the mean-field parametrization of Welke et al.. This opens up the possibility of exploring the effects of mean-field anisotropy in collisions, without the penalty of computational cost behind the folding parametrization. △ Less

Submitted 22 August, 2013; v1 submitted 18 December, 2012; originally announced December 2012.

Comments: 18 pages, 14 figures, expanded justification for the effort discussed in the paper, introduced other minor corrections

Journal ref: Phys. Rev. C 87, 054619 (2013)

Probing the symmetry energy with isospin ratio from nucleons to fragments

Authors: Yingxun Zhang, Zhuxia Li, Chengshuang Zhou, Jixian Chen, M. Colonna, P. Danielewicz, M. B. Tsang

Abstract: Within the framework of ImQMD05, we study several isospin sensitive observables, such as DR(n/p) ratios, isospin transport ratio (isospin diffusion), yield ratios for LCPs between the projectile region and mid-rapidity region for the reaction systems Ni+Ni, Zn+Zn, Sn+Sn at low-intermediate energies. Our results show that those observables are sensitive to the density dependence of symmetry energy,… ▽ More Within the framework of ImQMD05, we study several isospin sensitive observables, such as DR(n/p) ratios, isospin transport ratio (isospin diffusion), yield ratios for LCPs between the projectile region and mid-rapidity region for the reaction systems Ni+Ni, Zn+Zn, Sn+Sn at low-intermediate energies. Our results show that those observables are sensitive to the density dependence of symmetry energy, and also depend on the cluster formation mechanism. By comparing these calculations to the data, the information of the symmetry energy and reaction mechanism is obtained. △ Less

Submitted 31 August, 2012; originally announced August 2012.

Comments: Talk given by Yingxun Zhang at the 11th International Conference on Nucleus-Nucleus Collisions (NN2012), San Antonio, Texas, USA, May 27-June 1, 2012. To appear in the NN2012 Proceedings in Journal of Physics: Conference Series (JPCS)

Constraints on the symmetry energy and neutron skins from experiments and theory

Authors: M. B. Tsang, J. R. Stone, F. Camera, P. Danielewicz, S. Gandolfi, K. Hebeler, C. J. Horowitz, Jenny Lee, W. G. Lynch, Z. Kohley, R. Lemmon, P. Moller, T. Murakami, S. Riordan, X. Roca-Maza, F. Sammarruca, A. W. Steiner, I. Vidaña, S. J. Yennello

Abstract: The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints… ▽ More The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined. △ Less

Submitted 2 April, 2012; originally announced April 2012.

Comments: 28 pages, 4 figures

Symmetry Energy from Systematic of Isobaric Analog States

Authors: Pawel Danielewicz, Jenny Lee

Abstract: 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 is… ▽ More 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. △ Less

Submitted 1 November, 2011; originally announced November 2011.

Comments: 6 pages, 3 figures; talk given at IX Latin American Symposium on Nuclear Physics and Applications, July 18-22, 2011, Quito, Ecuador

Angular Dependence in Proton-Proton Correlation Functions in Central $^{40}Ca+^{40}Ca$ and $^{48}Ca+^{48}Ca$ Reactions

Authors: V. Henzl, M. A. Kilburn, Z. Chajecki, D. Henzlova, W. G. Lynch, D. Brown, A. Chbihi, D. Coupland, P. Danielewicz, R. deSouza, M. Famiano, C. Herlitzius, S. Hudan, Jenny Lee, S. Lukyanov, A. M. Rogers, A. Sanetullaev, L. Sobotka, Z. Y. Sun, M. B. Tsang, A. Vander Molen, G. Verde, M. Wallace, M. Youngs

Abstract: The angular dependence of proton-proton correlation functions is studied in central $^{40}Ca+^{40}Ca$ and $^{48}Ca+^{48}Ca$ nuclear reactions at E=80 MeV/A. Measurements were performed with the HiRA detector complemented by the 4$π$ Array at NSCL. A striking angular dependence in the laboratory frame is found within p-p correlation functions for both systems that greatly exceeds the measured and e… ▽ More The angular dependence of proton-proton correlation functions is studied in central $^{40}Ca+^{40}Ca$ and $^{48}Ca+^{48}Ca$ nuclear reactions at E=80 MeV/A. Measurements were performed with the HiRA detector complemented by the 4$π$ Array at NSCL. A striking angular dependence in the laboratory frame is found within p-p correlation functions for both systems that greatly exceeds the measured and expected isospin dependent difference between the neutron-rich and neutron-deficient systems. Sources measured at backward angles reflect the participant zone of the reaction, while much larger sources observed at forward angles reflect the expanding, fragmenting and evaporating projectile remnants. The decrease of the size of the source with increasing momentum is observed at backward angles while a weaker trend in the opposite direction is observed at forward angles. The results are compared to the theoretical calculations using the BUU transport model. △ Less

Submitted 11 August, 2011; originally announced August 2011.

Comments: 8 pages, 3 figures, submitted to PRC

Journal ref: Phys.Rev. C85 (2012) 014606

Influence of Transport Variables on Isospin Transport Ratios

Authors: D. D. S. Coupland, W. G. Lynch, M. B. Tsang, P. Danielewicz, Yingxun Zhang

Abstract: The symmetry energy in the nuclear equation of state affects many aspects of nuclear astrophysics, nuclear structure, and nuclear reactions. Recent constraints from heavy ion collisions, including isospin diffusion observables, have started to put constraints on the symmetry energy below nuclear saturation density, but these constraints depend on the employed transport model and input physics othe… ▽ More The symmetry energy in the nuclear equation of state affects many aspects of nuclear astrophysics, nuclear structure, and nuclear reactions. Recent constraints from heavy ion collisions, including isospin diffusion observables, have started to put constraints on the symmetry energy below nuclear saturation density, but these constraints depend on the employed transport model and input physics other than the symmetry energy. To understand these dependencies, we study the influence of the symmetry energy, isoscaler mean field compressibility and momentum dependence, in-medium nucleon-nucleon cross sections, and light cluster production on isospin diffusion within the pBUU transport code. In addition to the symmetry energy, several uncertain issues strongly affect isospin diffusion, most notably the cross sections and cluster production. In addition, there is a difference in the calculated isospin transport ratios, depending upon whether they are computed using the isospin asymmetry of either the residue or of all forward moving fragments. Measurements that compare the isospin transport ratios of these two quantities would help place constraints on the input physics, such as the density dependence of the symmetry energy. △ Less

Submitted 19 July, 2011; originally announced July 2011.

Comments: 31 pages, 15 figures, submitted to PRC

Constraints on the density dependence of the symmetry energy from heavy ion collisions

Authors: M. B. Tsang, Z. Chajecki, D. Coupland, P. Danielewicz, F. Famiano, R. Hodges, M. Kilburn, F. Lu, W. G. Lynch, J. Winkelbauer, M. Youngs, YingXun Zhang

Abstract: Constraints on the Equation of State for symmetric matter (equal neutron and proton numbers) have been extracted from energetic collisions of heavy ions over a range of energies. Collisions of neutron-deficient and neutron-rich heavy ions now provide initial constraints on the EoS of neutron-rich matter at sub-saturation densities from isospin diffusions and neutron proton ratios. This article rev… ▽ More Constraints on the Equation of State for symmetric matter (equal neutron and proton numbers) have been extracted from energetic collisions of heavy ions over a range of energies. Collisions of neutron-deficient and neutron-rich heavy ions now provide initial constraints on the EoS of neutron-rich matter at sub-saturation densities from isospin diffusions and neutron proton ratios. This article reviews the experimental constraints on the density dependence of Symmetry Energy at sub-saturation density. △ Less

Submitted 19 January, 2011; originally announced January 2011.

Comments: 7 pages, 3 figures, to be published in: Progress in Particle and Nuclear Physics

Journal ref: Prog.Part.Nucl.Phys.66:400-404,2011

Deducing spectroscopic factors from wave-function asymptotics

Authors: P. Capel, P. Danielewicz, F. M. Nunes

Abstract: In a coupled-channel model, we explore the effects of coupling between configurations on the radial behavior of the wave function and, in particular, on the spectroscopic factor (SF) and the asymptotic normalization coefficient (ANC). We evaluate the extraction of a SF from the ratio of the ANC of the coupled-channel model to that of a single-particle approximation of the wave function. We perform… ▽ More In a coupled-channel model, we explore the effects of coupling between configurations on the radial behavior of the wave function and, in particular, on the spectroscopic factor (SF) and the asymptotic normalization coefficient (ANC). We evaluate the extraction of a SF from the ratio of the ANC of the coupled-channel model to that of a single-particle approximation of the wave function. We perform this study within a core + n collective model, which includes two states of the core that connect by a rotational coupling. To get additional insights, we also use a simplified model that takes a delta function for the coupling potential. Calculations are performed for 11Be. Fair agreement is obtained between the SF inferred from the single-particle approximation and the one obtained within the coupled-channel models. Significant discrepancies are observed only for large coupling strength and/or large admixture, that is, a small SF. This suggests that reliable SFs can be deduced from the wave-function asymptotics when the structure is dominated by one configuration, that is, for a large SF. △ Less

Submitted 8 December, 2010; v1 submitted 15 September, 2010; originally announced September 2010.

Comments: Title corrected

Journal ref: Phys.Rev.C82:054612,2010

The Influence of in-medium NN cross-sections, symmetry potential and impact parameter on the isospin observables

Authors: Yingxun Zhang, D. D. S. Coupland, P. Danielewicz, Zhuxia Li, Hang Liu, Fei Lu, W. G. Lynch, M. B. Tsang

Abstract: We explore the influence of in-medium nucleon-nucleon cross section, symmetry potential and impact parameter on isospin sensitive observables in intermediate-energy heavy-ion collisions with the ImQMD05 code, a modified version of Quantum Molecular Dynamics model. At incident velocities above the Fermi velocity, we find that the density dependence of symmetry potential plays a more important role… ▽ More We explore the influence of in-medium nucleon-nucleon cross section, symmetry potential and impact parameter on isospin sensitive observables in intermediate-energy heavy-ion collisions with the ImQMD05 code, a modified version of Quantum Molecular Dynamics model. At incident velocities above the Fermi velocity, we find that the density dependence of symmetry potential plays a more important role on the double neutron to proton ratio $DR(n/p)$ and the isospin transport ratio $R_i$ than the in-medium nucleon-nucleon cross sections, provided that the latter are constrained to a fixed total NN collision rate. We also explore both $DR(n/p)$ and $R_i$ as a function of the impact parameter. Since the copious production of intermediate mass fragments is a distinguishing feature of intermediate-energy heavy-ion collisions, we examine the isospin transport ratios constructed from different groups of fragments. We find that the values of the isospin transport ratios for projectile rapidity fragments with $Z\ge20$ are greater than those constructed from the entire projectile rapidity source. We believe experimental investigations of this phenomenon can be performed. These may provide significant tests of fragmentation time scales predicted by ImQMD calculations. △ Less

Submitted 20 January, 2012; v1 submitted 10 September, 2010; originally announced September 2010.

Comments: 24 pages, 9 figures, to be published in Phys. Rev. C