Search

Source functions are obtained from $p$-$p$ and $d$--$\alpha$ correlation functions by applying the Richardson-Lucy (RL) deblurring to the Koonin-Pratt (KP) equation. To prevent fitting of noise in the correlation function, total-variation (TV) regularization is employed that has been effective in ordinary image restoration. TV alone cannot ensure normalization of the source functions. To ensure the latter, we propose a maximum-entropy regularized RL algorithm (MEM-RL). We outline the MEM-RL formalism and optimization strategy for the KP equation, demonstrating its effectiveness on both simulated and experimental data, including the $p$-$p$ and $d$-$\alpha$ correlation functions., Comment: 22 pages, 6 figures

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., Comment: 17 pages, 9 figures

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., Comment: 5 pages, 3 figures

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 $\Delta$ 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 $\pi^-/\pi^+$ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced $\pi^-/\pi^+$ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total $\pi^-/\pi^+$ 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\sigma$ 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., Comment: 39 pages, 17 figures

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.

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., Comment: 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

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, Comment: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

$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 $\rho_{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 $\rho_{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 $\Delta R_{n p}$. $Conclusions:$ Results illustrate that the charge-exchange ($p, n$) reaction is a sensitive probe of $\rho_{n}$. The results in this paper can offer useful guides for future experiments of neutron characteristics.

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., Comment: 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

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.

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., Comment: 10 pages, 9 figures, 2 tables; Physical Review C, in press

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., Comment: 114 pages, 14 figures, 479 references, accepted for publication in Progress of Particle and Nuclear Phsics

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., Comment: 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

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., Comment: 17 pages,11 figures, 26 references

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.

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$\pi$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., Comment: 8 pages, 4 figures, 1 table (accepted for publication in PLB)

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 ($\Lambda$) 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 $\Lambda$. The $\Lambda$ 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.

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.

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., Comment: 9 pages, 5 figures

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$\rho_{0}$ to $4.5\rho_{0}$. While the uncertainties in the symmetry pressure are large, they can be reduced with new experimental and astrophysical results., Comment: arXiv admin note: substantial text overlap with arXiv:1901.07673

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, $\Lambda$, 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 $\Lambda$, 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 $\Lambda$ 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 $\Lambda$ and total pressure providing guidance to laboratory nucleus-nucleus collision experiments.

We compare ten transport codes for a system confined in a box, aiming at improved handling of the production of $\Delta$ resonances and pions, which is indispensable for constraining high-density symmetry energy from observables such as the $\pi^-/\pi^+$ 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\Delta$ and $\Delta\leftrightarrow N\pi$ 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 $\Delta$ and $\pi$, 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 $\Delta$ and $\pi$. 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 $\pi^-/\pi^+$ ratio for the system initialized at n/p = 1.5, after letting the existing $\Delta$ 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., Comment: 36 pages, 27 figures; a new Fig. 21 and revised results from some codes, achieving improved and consistent understanding

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., Comment: 10 pages, 3 figures. arXiv admin note: text overlap with arXiv:1807.06571

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., Comment: arXiv admin note: substantial text overlap with arXiv:1807.06571

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 $\eta$ in the Green's function, practical results can be further improved., Comment: 11 pages, 10 figures, typos corrected, accepted for publication in Phys. Rev. C

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., Comment: 11 pages, 3 figures

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., Comment: 11 pages, 7 figures

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., Comment: Contribution to Proc. Frontiers of Quantum and Mesoscopic Thermodynamics, 9-15 July 2017, Prague, Czech Republic, ed. Vaclav Spicka; 10 pages, 7 figures

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., Comment: 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

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., Comment: Contribution to Proc. 5th Conference on Nuclei and Mesoscopic Physics, E Lansing, 6-10 March 2017; 9 pages, 8 figures

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 $\delta$ 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., Comment: 36 pages, 8 figures, submitted to Phys. Rev. C

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.

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., Comment: 31 pages, 12 figures

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., Comment: 12 pages, 12 figures, a few points clarified, typos corrected, accepted for publication in Phys. Rev. C

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

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., Comment: 7 pages, 7 figures, accepted for publication in Eur. Phys. J. A

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., Comment: 26 pages, 13 figures, a manuscript on the heavy-ion collision part of transport code comparison project